基于Fortran语言的地球外辐射带电子三维数据同化建模

开发地球电子辐射带的数据同化模型,对于理解辐射带电子的动态演化过程和辐射带空间天气预报具有重要意义.结合范阿伦卫星的辐射带电子观测数据和外辐射带三维扩散模型,采用卡尔曼滤波算法,本文开发了基于Fortran语言的外辐射带电子三维数据同化模型(Three-dimensional Data Assimilative Model of Outer Radiation belt Electrons,简称TDAMORE),实现对L^(*)=3~7、能量范围为0.1~5 MeV、投掷角范围为5°~90°的外辐射带电子时空变化过程的三维重构.通过对2018年8月期间外辐射带电子通量演化过程的重构,证实TDAMORE模型可以较好地重现不同能量和不同投掷角电子通量在磁暴前后的演化特征.通过分析电子通量的观测和同化结果之间的相关系数、平均误差、平均绝对误差和均方误差,发现对于能量低于4 MeV的电子,观测与同化结果之间的相关系数基本大于0.8且误差相对较低.而对于更高能量的电子,观测与同化结果之间的误差相对较高,这可能是同化模型忽略了电磁离子回旋波对电子的散射损失导致的.

地球电子外辐射带的数据同化建模与分析

地球电子外辐射带对太阳与地磁活动呈现高度动态变化的响应,了解外辐射带的全球动态变化过程对于近地空间粒子辐射环境的理解认知和预测预报具有重要意义.基于卡尔曼滤波数据同化方法,本文利用范阿伦A星、B星和GOES-13和GOES-15四颗卫星的辐射带电子观测数据,分别利用三种不同维度的辐射带物理模型,将观测结果与数值结果有机融合,对2013年3月地球外辐射带电子通量的径向分布与变化进行数据同化分析.结果表明,考虑了磁层波动与辐射带电子共振作用引起的径向扩散、投掷角扩散以及能量扩散过程的三维同化模型可有效、合理地重现外辐射带电子通量的径向分布.本文进一步利用该三维同化模型对2013年一整年外辐射带电子的相空间密度分布进行重构与分析,得到了不同绝热不变量和不同地磁活动条件下电子辐射带的时空演化过程,从而为深入理解外辐射带电子的变化过程和动力学机制提供了强有力信息.通过分析同化过程中的新息矢量以及度量同化过程中观测数据在多大程度上修改了物理模型结果,还有助于定量分析现有辐射带物理模型中的源项和损失项的相对贡献以及可能忽略的物理机制或过程.