关键词:
地球外辐射带
电子通量
电子相空间密度
多维数据同化
摘要:
基于Explorer 1飞船的观测数据,科学家们首次发现了地球辐射带,这是太空时代的重大发现,开启了磁层物理研究的大门。随着更多科学卫星的发射和更精密仪器的使用,人们对辐射带有了更深刻的认识。在多种物理机制的共同作用下,外辐射带中的高能电子呈现着非常复杂的剧烈动态变化特征。由于辐射带电子会对在轨航天器与宇航员的健康造成威胁,理解和预报辐射带电子的动态变化有着十分重要的意义。尽管已有大量的卫星观测数据可供使用,但卫星的局地观测受到时间和轨道的限制,测量结果较为稀疏;而基于物理的数值模型往往误差较大,且初边界条件很难确定。数据同化可以将卫星的局地观测与数值模型的全球建模有机结合,进而实现对地球辐射带电子时空分布与动态演化的重构与预测。因此,辐射带数据同化不仅是磁层物理领域的研究前沿,也是空间天气预报的重要支撑。
本文使用了范阿伦A星、B星和GOES-13、15四颗卫星的辐射带电子观测数据,以及基于Fokker-Planck方程的辐射带数值模型,采用基于卡尔曼滤波的数据同化方法,对地球外辐射带电子进行了多维数据同化建模与分析。主要结论如下:
(1)本文首先对2013年3月的辐射带电子通量进行了同化分析,通过在数值模拟阶段加入不同的扩散机制,来比较不同物理过程对辐射带电子演化过程的影响。本文将同化模型区分为:只考虑径向扩散的一维同化模型;考虑径向扩散和投掷角扩散的二维同化模型;考虑径向扩散、投掷角扩散和能量扩散的三维同化模型。其中,一维同化模型可以模拟出辐射带电子通量随地磁活动的基本变化趋势,但无法还原出电子动态演变的精细轮廓。二维同化模型加入了投掷角扩散机制,可以在等离子体层顶内有效地损失电子。三维同化模型加入了能量扩散机制,可以重现等离子体层顶外电子通量在磁暴恢复相的累积和峰值。通过比较不同地磁环境下的数值模拟结果和同化结果可以看出,数据同化过程对数值模拟结果进行了修正使其更加趋近于卫星观测。本文进一步对比了三组同化模型的数据同化结果,结果表明同时考虑了径向扩散、投掷角扩散和能量扩散的三维同化模型结果在所有L*和多组能量投掷角组合均与观测吻合较好,可以合理地重现辐射带电子通量的全球分布。
(2)本文进一步使用三维同化模型对2012-2015年辐射带电子通量的长期分布进行了同化分析。同化结果弥补了观测数据在时空分布上的缺失,完整地模拟出电子通量的径向分布和时空演化。为了区分绝热和非绝热效应以理解辐射带电子加速损失的平衡,本文进一步对辐射带电子相空间密度进行了三维同化分析。同化结果表明电子相空间密度的时空变化与磁暴事件具有很强的依赖性:磁暴主相期间电子相空间密度会在大多数径向范围上的快速降低,而在磁暴恢复相期间电子相空间密度一般会增强并达到峰值。通过对不同地磁活动条件下的辐射带电子相空间密度长期分布进行同化分析,本文所使用的三维同化模型可以准确地重构外辐射带电子相空间密度的时空演化过程,从而有效地应用于研究影响辐射带电子加速与损失过程的物理机制。
(3)为了量化同化过程对数值模拟的修改程度,本文计算了2012-2015年同化过程中的平均新息矢量(Innovation vector)。首先以2013年的平均新息矢量结果为例,比较了不同扩散机制下新息矢量的数值变化和变化范围。可以看出,新息矢量的正负会随加入的扩散机制不同改变,且变化趋势符合所添加扩散机制的原理,说明模拟过程中加入的扩散机制都是合理有效的。本文还计算了2012-2015年三维同化结果的新息矢量,结果表明数值模型在低L*略低于观测值,而在高L*略高于观测值,这可能是数值模拟阶段使用的波动模型不够准确导致的。新息矢量的数值与地磁活动有明显的相关性,说明本文使用数值模型的模拟效果在地磁活动剧烈时存在更大的不确定性。通过研究新息矢量,可以为后续模型的改进和优化提供有效的参考。