pub fn ice_tectonic_rate_mm_myr(tidal_stress_mpa: f64) -> f64 {
    0.8 + 2.5 * tidal_stress_mpa.max(0.0)
}

pub fn sublimation_rate_mm_yr(temperature_k: f64) -> f64 {
    (temperature_k - 50.0).max(0.0) * 0.01
}

pub fn space_weathering_index(exposure_myr: f64, solar_wind_factor: f64) -> f64 {
    (1.0 - (-exposure_myr / 60.0).exp()) * solar_wind_factor.max(0.0)
}
