AEGIS 将火星大气物理耦合到可微 Dinosaur 动力核心，十年积分稳定并复现季节 CO2 循环
相关研究与后续进展
核心概要
作者提出 AEGIS，一个把火星辐射、地表、CO2 霜与风化层物理耦合到 Dinosaur 谱动力核心的可微火星气候模型，在 T21/L12 配置下完成十个火星年的稳定积分，复现季节 CO2 循环并把大气—霜总 CO2 库存守恒到极小漂移，地表气压跟随 MOLA 地形，且自动微分梯度与有限差分一致，可用于物理标定与神经闭合训练。
Figure 1: AEGIS architecture. Planetary constants and boundary data configure a shared GCM. Dinosaur supplies resolved dynamics; Mars physical operators supply column tendencies and surface evolution; neural components enter as bounded residuals or component replacements. The coupled time step advances atmosphere and surface reservoirs. Trajectory objectives connect the simulated fields to physical-parameter calibration and neural training through automatic differentiation. Dashed arrows denote this gradient-based update path.
arXiv深度剖析
AEGIS 把行星单位重标定、耦合的大气—地表—霜—风化层状态、CO2 相变进入预报对数气压以及每步大气—霜库存投影加入 Dinosaur 核心，使质量交换与地表记忆成为可微轨迹的一部分。 此前的 Dinosaur 系模型（NeuralGCM、JCM）面向地球配置；AEGIS 在动力接口使用火星半径、自转与热力学尺度，并新增火星辐射、地表与 CO2 凝结循环算子。 方法部分给出耦合状态方程、CO2 凝结对预报对数气压的贡献式与逐步库存投影；附录 A 详述辐射、地表、边界层与 CO2 库存算子。
通过耦合轨迹的梯度经有限差分验证：54/54 物理敏感性检查通过，多诊断雅可比九个导数与中心差分一致，最大相对偏差很小。 既有火星 GCM 多为传统 Fortran 有限差分/有限体积实现，不暴露梯度；AEGIS 在同一实现中同时提供前向与反向模式微分。 表 1 列出四类梯度检查全部通过；多诊断实验在 0.25、1、5 sol 三个时间尺度上比较 AD 与有限差分，并指出气压与霜的导数等大反号，其和受库存投影约束。
梯度标定在十二窗口试点中降低各数据划分的大气温度 RMSE：训练 2.45→2.38 K、验证 2.49→2.44 K、测试 2.46→2.40 K，四个测试窗口全部改善。 这展示了可微轨迹损失能给出方向正确的优化信号，而传统参数扫描与集合标定通常计算昂贵。 表 2 给出物理基线、标定物理与耦合神经加热三种配置的分划分 RMSE；附录 B 说明十二个 ARCO-MACDA 窗口的划分与 50 次 Adam 更新协议。
长期积分稳定且守恒：一个火星年 86.6 分钟完成，十个火星年 11.4 小时完成，全部保存状态有限，大气—霜 CO2 库存最大相对漂移极小，季节循环在第十年与第九年几乎重合。 模型复现季节 CO2 循环（全球平均地表气压 485–654 Pa，约 30% 年振幅）与 MOLA 地形驱动的地表气压结构，十二次季节—参考对比的气压图样相关为 0.975–0.989。 性能与守恒数据来自单 GPU 运行统计；表 3 给出与 NASA Ames、MCD 6.1、ARCO-MACDA 三种参考产品的 RMSE 与相关；消融套件中十三种配置均以有限状态完成。
启示与展望
该工作面向火星大气这一具体物理域，在 T21 三角截断、12 层 sigma 与 300 s 步长的粗分辨率配置下评估，适用于物理敏感性研究、梯度标定与神经闭合的受控实验。作者指出框架已支持更广的观测对比、跨随机种子的重复神经实验与多年平衡研究，并可通过显式行星、辐射与储库接口引入改变的地表反照率、外加加热或指定挥发分库存变化，再用空间气候诊断评估；对大气成分与质量的更大改动，需要同时扩展行星热力学、不透明度数据、可凝结过程与适用的动力假设。
季节对比中近地表风速是各参考产品差异最大的场，因为产品在报告层次与平均方式上不同，作者据此把比较定位为模型与再分析对比；地表温度图样相关为 0.732–0.821，低于气压相关。神经辐射实验中，继续局部拟合在插值与高沙尘外推上均优于轨迹微调，提示局部与耦合轨迹目标之间的取舍仍是开放问题。消融在 44,400 s 长时程上自由轨迹达到混沌饱和，过程敏感性需从短时程读取。此外，正文中若干数值（如梯度偏差、库存漂移、闭合残差的具体量级）在文本中以符号或占位形式出现，读者若需精确数值应查阅原文表格与图件。