航天器轨道的精确计算在航天器测控、空间交会对接、空间目标碰撞预警等任务中具有重要作用。为了准确计算航天器在轨运行轨道，需要精确掌握航天器在轨运行中受到的摄动力，这些摄动力包括地球非球形引力、大气阻力、太阳光压、潮汐力等。针对这些摄动因素，研究者对其中绝大多数因素已经构建了精确度相当高的物理或经验模型，如重力场模型、潮汐模型等。但受到热层大气密度变化的不确定性影响，大气阻力模型的误差仍然比较高，成为影响低轨航天器精密定轨和轨道预报的最重要制约因素。
几十年来随着大气探测技术的发展和探测数据的日益增多，人们对热层大气的变化规律有了更加深入的认识，但由于热层大气的变化极其复杂，物理机制尚不清楚外，已有的大气经验模型自身也存在缺陷，地磁平静期各系列模型误差平均在15％，地磁暴期间误差甚至达到100%。因此，研究和掌握地磁暴期的大气密度变化特征和规律，是改进和建立高精度大气模型的基础。
为了研究地磁暴期间大气密度变化特征，中国科学院国家空间科学中心太阳活动与空间天气重点实验室空间天气预报中心苗娟研究员和博士生王昕，对2001-2008年期间的256个地磁暴事例，按照小、中、大3种磁暴强度对大气密度变化特征、大气密度与焦耳加热分布特征进行了详细分析和研究。
研究结果揭示了焦耳加热与大气密度的位置和时间变化关系：小磁暴和中等磁暴期间，两者峰值时间相差0-2hr，纬度相差0-10o；强磁暴期间时间相差3-5hr，纬度相差10-15o，纬度的差异可能是由中性风的影响造成的。
图1. (a)大气密度滞后焦耳加热峰值时间,（b）焦耳加热峰值持续时间
图2. 2003年11月20日地磁暴期间大气密度及焦耳加热位置变化：（a-d）焦耳加热变化，（e-h）大气密度变化量，(i-l) 大气密度误差
苗娟等在大气密度和焦耳加热变化关系研究基础上，进一步加入地磁Dst指数和AE指数与大气密度变化的特征分析，建立了综合的地磁暴大气密度修正函数，验证结果显示修正后的大气密度相对误差（MRE）由40% 下降至 10% 。
相关成果发表在国际知名期刊JGR和Space Weather上。JGR审稿评价认为“文章的结果与讨论，有利于提高对热层大气中加热和冷却机制之间复杂相互作用的理解，这是目前该领域研究非常感兴趣的一个话题，这些结论是当前研究中重要的新发现”。
图3. 大气密度修正值（蓝线），大气密度实测值（红线）和NRLMSISE-00模型计算值（黑线）的比较
Citation：
（1）Wang, X., Miao, J., Aa, E., Ren, T.,Wang, Y., & Liu, J., et al. (2020).Statistical analysis of Joule heating and thermosphere response during geomagnetic storms of different magnitudes. Journal of Geophysical Research: Space Physics, 125, e2020JA027966. . https://doi.org/10.1029/2020JA027966
（2）Wang, X., Miao, J. , Lu, X., Aa, E., Liu, J., Wang, Y., & Liu, S. (2021). Latitudinal impacts of Joule heating on the high-latitude thermospheric density enhancement during geomagnetic storms. Journal of Geophysical Research: Space Physics, 126, e2020JA028747. https://doi.org/10.1029/2020JA028747
（3）Wang, X., Miao, J., Lu, X., Aa, E., Luo, B., Liu, J., et al. (2022). Using temporal relationship of thermospheric density with geomagnetic activity indices and Joule heating as calibration for NRLMSISE-00 during geomagnetic storms. Space Weather, 20, e2021SW003017. https://doi.org/10.1029/2021SW003017
（供稿：天气室）
