基于范阿伦双星EMFISIS观测数据的NWC和NAA人工甚低频台站信号的内磁层全球统计分布

地球表面的人工甚低频台站信号可以穿透电离层泄漏进地球磁层导致内辐射带电子沉降到两极大气.因此研究人工甚低频台站信号的空间全球分布特性对于分析辐射带电子的损失具有重要科学意义.本文使用范阿伦双星从2013年到2018年共计6年的高质量的波动观测数据,统计了 NWC(19.8 kHz)、NAA(24.0 kHz)两个人工VLF台站信号的全球分布,分析了台站信号的电场功率谱密度对地理经纬度、磁壳值L、磁地方时MLT、地磁活动水平的依赖性.结果表明,在内磁层中,人工台站VLF信号主要沿着台站位置对应的磁力线传播,夜侧强度高于日侧,冬季高于夏季.这种日夜和夏冬差异的形成是因为夜侧和冬季的日照强度较弱,电离层电子密度较低,VLF信号较容易穿透电离层进入磁层.此外人工VLF台站信号的全球分布受地磁活动的影响很弱.这些统计观测结果给出了 NWC和NAA两个重要人工VLF台站信号强度的全球分布特征,为进一步分析人工VLF台站信号与地球辐射带电子的波粒相互作用提供了关键信息.     

基于Van Allen Probes EMFISIS波动仪器观测的内磁层下频带哨声合声波全球分布的统计分析

基于Van Allen Probes近三年的EMFISIS仪器波动观测数据,对内磁层下频带哨声模合声波幅度的全球分布特性对地磁活动水平的依赖性进行了详细的统计分析,着重研究下频带合声波平均场强幅度随磁壳值L、磁地方时、地磁纬度的分布特征及不同强度区间的合声波的发生概率.结果表明,下频带合声波的波动强度与地磁活动密切正相关,处于强磁扰期间的合声波具有更大的振幅,其发生率与地磁活动强度具有同样的正相关特性.下频带合声波主要发生于午夜至下午的磁地方时区间,其余的磁地方时时段下频带合声波较弱.赤道面附近的下频带合声波主要分布在夜侧至黎明这一时段内,随着磁纬度的增加逐步向日侧扩展.下频带合声波在午夜侧(21-03 MLT)主要出现在15°的磁纬范围内,在晨侧(03-09 MLT)可以到达15°磁纬甚至更高纬度.下频带合声波主要发生于L=~4.5的附近区域.随着地磁活动的增加,下频带合声波所覆盖的L-shell空间区域增大,趋势为向高、低L值区域同时扩展.建立的下频带哨声合声波的全球分布模型将有助于进一步深入理解该重要磁层波动对辐射带电子的波粒作用散射效应和对辐射带动力学过程的定量贡献.     

基于Van Allen Probes EMFISIS波动仪器观测的内磁层上频带哨声合声波全球分布的统计分析

基于Van Allen Probes近三年的EMFISIS仪器波动观测数据,针对内磁层上频带哨声模合声波幅度的全球分布特性对地磁活动水平的依赖性进行了详细的统计分析,着重研究上频带合声波平均场强幅度随磁壳值(L)、磁地方时(MLT)、地磁纬度(MLAT)的分布特征及不同强度区间的合声波的发生概率.结果表明,上频带合声波的平均场强幅度与地磁活动条件密切相关,在强磁扰期,平均幅度可达到40 pT以上.在外辐射带中心区域(L=4~6),上频带合声波的幅度最强;在L~3的区域,上频带磁层合声波没有分布.在夜侧至晨侧(22—09MLT),上频带合声波幅度最强;在下午侧至昏侧(15-19MLT),上频带合声波幅度最弱;日侧(10-14MLT)上频带合声波在不同地磁活动条件下都存在,幅度偏小.上频带合声波主要分布在|MLAT|10°,其中21-09MLT范围内、磁纬位于|MLAT丨5°的平均场强幅度最强,磁扰期间可达约100 pT.另外,统计而言,中等幅度(10~30 pT)的上频带合声波在夜侧至晨侧(23-09MLT)靠近磁赤道区域的发生率最高,可达15%左右.强幅度(30 pT)的上频带合声波普遍分布在夜侧(01-05MLT),发生率最小.本文建立的上频带哨声模合声波的全球分布模型结合已经建立的下频带合声波的全球分布模型,将有助于进一步深入理解该重要磁层等离子体波动对地球等离子体片、辐射带、环电流动力学过程的定量贡献.     

伴随等离子体密度下降的磁声波与辐射带电子的波粒相互作用及其散射效应

利用范阿伦卫星的高质量观测数据,我们报道了伴随等离子体密度下降的磁声波现象.通过选取分别发生于2013年7月26日(事件A)和2013年9月19日(事件B)的两个相应事件进行细致分析,我们开展试验粒子模拟计算了磁声波对辐射带电子的散射系数,并求解二维福克-普朗克扩散方程量化了磁声波散射导致的辐射带电子动态变化.结果表明,事件A中的磁声波的散射作用主要发生于投掷角范围为60°～80°、能量范围为20～200keV的辐射带电子,而事件B中的磁声波的散射作用主要发生于投掷角范围为50°～80°、能量范围为20～400keV的辐射带电子;两个事件中的磁声波均能导致辐射带电子的蝴蝶状投掷角分布,但是由于事件B的磁声波幅度更强,形成的电子蝴蝶状分布更明显.     

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

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

风云3C对不同类型辐射带电子分布变化的观测：2015年3月17日磁暴事件分析

风云3C卫星(简称FY-3C)测量的低地球轨道(Low Earth Orbit,简称LEO)高能粒子(电子、质子和重离子)通量,是空间环境监测预警业务的重要自主数据来源.深入分析风云3C卫星测量的不同种类高能电子的分布特征和变化规律,对辐射带研究和空间天气预报具有重要意义.本文使用FY-3C测量的空间粒子数据,系统分析了2015年3月17日强磁暴期间辐射带捕获、准捕获和沉降电子的演化特征及其主导性物理机制.结果表明,在磁暴主相,外辐射带电子通量(E2:0.35～0.65 MeV, E3:0.65～1.2 MeV)出现快速下降,达到两个数量级,E1能级(E1:0.15～0.35 MeV)电子通量在L<5区域因为亚暴注入先出现增长,之后在L>4的区域出现明显下降.在磁暴恢复相,由于合声波(chorus waves)先加速低能段电子,E1能段电子通量最先(小于五小时)开始增加,在磁暴结束时覆盖25 cm^-2·sr^-1·s^-1.电子能段越高通量增长得越慢...     

基于JUNO卫星WAVES仪器波动数据的木星磁层哨声波分析研究

等离子体波的空间分布在木星磁层高能电子动力学过程中起着重要的作用.现有大多数对木星磁层哨声波的观测仅限于|λ|≤15°的磁纬范围内,但是最新的JUNO卫星WAVES仪器提供的波动数据使得更高纬度、更广区域范围内的等离子体波动分布研究成为可能.本文通过对JUNO卫星WAVES仪器数据进行分析处理,详细研究了木星磁层哨声波的空间分布特性.观测表明,存在位于高L_J、高磁纬的木星磁层哨声波,它们广泛分布于距木星中心距离35～75个木星半径、磁纬为|λ|≤30°的空间区域.分析研究发现,WAVES仪器观测的木星磁层哨声波幅度一般为几个pT,远小于地球磁层哨声波的强度.木星磁层哨声波幅会随着L_J的增大缓慢增加,也会随着磁纬的减小趋向平缓变化.基于以上观测事实,本文利用指数幂函数分别拟合得到木星磁层哨声波幅随L_J和磁纬变化的经验模型.该模型将有助于精确理解哨声波对木星磁层高能电子动力学过程的重要影响.     

磁暴期间辐射带电子相空间密度的多卫星联合观测

地球外辐射带是一个高度动态变化的空间环境,辐射带电子通量的变化在磁暴期间尤为明显.要分析潜在的电子动态变化机制,需要排除绝热效应产生的影响.在以三个绝热不变量组成的相空间坐标中,利用相空间密度(PSD)可以反映电子的真实加速和损失情况.本文详细分析两颗范艾伦卫星和三颗GPS导航卫星在2013年3月的同步电子通量观测数据,发现在3月17日磁暴期间,当太阳风动压增大、行星际磁场南向时,辐射带电子通量会发生骤降.进一步将电子通量转换成电子相空间密度并利用不同第一、第二绝热不变量(μ,K)组合条件下PSD径向分布的差异性,深入探究磁暴期间辐射带电子的动态变化机制.结果表明:磁暴初期由于电子的局地加速导致PSD不断上升;磁暴主相期间,由于磁层顶阴影效应以及伴随的向外径向扩散损失导致PSD快速降低;位于不同空间位置的多颗卫星观测为明晰辐射带电子动态物理过程提供了重要的便利.     

第23太阳活动周期太阳风参数及地磁指数的统计分析

日冕物质抛射(Coronal Mass Ejection,简称CME)和共转相互作用区(Corotating Interaction Region,简称CIR)是造成日地空间行星际扰动和地磁扰动的两个主要原因,提供了地球磁暴的主要驱动力,进而显著影响地球空间环境.为深入研究太阳风活动及受其主导影响的地磁活动的时间分布特征,本文对大量太阳风参数及地磁活动指数的数据进行了详细分析.首先,采用由NASA OMNIWeb提供的太阳风参数及地磁活动指数的公开数据,通过自主编写matlab程序对第23太阳活动周期(1996-01-01—2008-12-31)的数据包括行星际磁场Bz分量、太阳风速度、太阳风质子密度、太阳风动压等重要太阳风参数及Dst指数、AE指数、Kp指数等主要的地磁指数进行统计分析,建立了包括269个CME事件和456个CIR事件列表的数据库.采用事例分析法和时间序列叠加法分别对两类太阳活动的四个重要太阳风参数(IMF Bz、太阳风速度、太阳风质子密度、太阳风动压)和三个主要地磁指数(Dst、AE、Kp)进行统计分析,并研究了其统计特征.其次,根据Dst指数最小值确定了第23太阳活动周期内的355个孤立地磁暴事件,并以Dst指数最小值为标准将这些磁暴进一步分类为145个弱磁暴、123个中等磁暴、70个强磁暴、12个剧烈磁暴和5个巨大磁暴.最后,采用时间序列叠加法对不同强度磁暴的太阳风参数和地磁指数进行统计分析.统计分析表明,对于CME事件,Nsw/Pdyn(Nsw表示太阳风质子密度,Pdyn表示太阳风动压)线性拟合斜率一般为正;对于CIR事件,Nsw/Pdyn线性拟合斜率一般为负,这可作为辨别CME和CIR事件的一种有效方法.从平均意义上讲,相较于CIR事件,CME事件有更大的南向IMF Bz分量、太阳风动压Pdyn、AE指数、Kp指数以及更小的Dstmin.一般情况下,CME事件有更大的可能性驱动极强地磁暴.总体而言,对于不同强度的地磁暴,Dst指数的变化呈现出一定的相似性,但随着地磁暴强度的增强,Dst指数衰减的速度变快.CME和CIR事件以及其各自驱动的地磁暴事件有着很多不同,因此,需要将CME事件驱动的磁暴及CIR事件驱动的磁暴分开研究.建立CME、CIR事件及地磁暴的数据库以及获取的统计分析结果,将为深入研究地球磁层等离子体片、辐射带及环电流对太阳活动的响应特征提供有利的帮助.