AEGIS: Differentiable Mars Climate Model with Neural Closures
Organizations: Department of Computer Science UC Santa Cruz Santa Cruz, CA 95064, USA
Abstract
General circulation models (GCMs) are the primary tool for simulating planetary atmospheres. They play a vital role in understanding Mars's atmosphere, as forecasting its unique weather is mission-critical for operations such as entry, descent, and landing. Mars poses unusual challenges for these models, as observations are sparse compared to Earth. In addition, a thin \co{} atmosphere alongside a radiatively active dust cycle creates a volatile atmosphere with large diurnal temperature swings and no true terrestrial analog for validation. Existing Mars GCMs, including the LMD PCM, the NASA Ames Mars GCM, and PlanetWRF, are mature and physically detailed but are implemented in legacy Fortran with finite-difference or finite-volume solvers, and they do not expose gradients for calibration or machine learning. Here we present AEGIS, a modular differentiable Mars climate model that couples Mars's unique atmospheric physics to the Dinosaur dynamical core, with interfaces for neural closures. We showcase stable ten-Mars-year simulations that reproduce the seasonal \co{} cycle while conserving the total \co{} inventory, capture realistic large-scale surface-temperature structure, and produce surface pressure that follows Mars Orbiter Laser Altimeter (MOLA) topography. Gradients through coupled trajectories agree with finite differences and support physical calibration and neural training. We compare with conventional GCMs, highlighting the framework's computational efficiency and differentiability.
Figures & tables
| Category | Objective | Horizon | Checks | Passed |
| Physical sensitivity | Final mean atmospheric temperature | 1–148 steps | 54 | 54/54 |
| Physical calibration | Scalar fitting objective | task-dependent | 9 | 9/9 |
| Neural fitting | Multi-field trajectory objective | 25 steps | 3 | 3/3 |
| Multi-diagnostic sensitivity | Mean surface , , and CO 2 frost | 74–1480 steps | 9 | 9/9 |
| Configuration | Training | Validation | Test |
|---|---|---|---|
| Physical baseline | 2.45 | 2.49 | 2.46 |
| Calibrated physics | 2.38 | 2.44 | 2.40 |
| Coupled neural heating | 2.42 | 2.49 | 2.46 |
| Pressure | Temperature | |||
| Reference | RMSE (Pa) | Correlation | RMSE (K) | Correlation |
| NASA Ames | 52.37 | 0.9864 | 20.64 | 0.7889 |
| MCD 6.1 | 57.18 | 0.9775 | 21.58 | 0.7872 |
| ARCO-MACDA | 40.72 | 0.9879 | 22.27 | 0.7715 |
| Training procedure | Interpolation | Higher dust |
|---|---|---|
| Initial local fit | 0.620 | 1.799 |
| Trajectory fine-tuning | 0.566 | 1.653 |
| Continued local fitting | 0.416 | 0.979 |
Appendix figures & tables3 assets
Supplementary material from the paper’s appendix.
Appendix
Climate Modeling
Climate modeling uses physical simulations and data-driven methods to understand, predict, and project Earth’s climate and its response to external forces.
See all 69 papers