亚暴事件中磁尾多重等离子体团的数值研究（Ⅱ）

以宁静磁尾平衡解为初态,本文考虑介质的可压缩性,对GEOTAIL卫星于亚暴事件中,观测到多重等离子体团与高能离子爆相对应之特征作模拟研究．数值结果表明：重复形成的等离子体团是高温、高密度区,揭示了等离子体团与高能离子爆──对应的动力学原因．本文的模拟结果还表明：大尺度等离子体团的地向运动,与X中性点的尾向迁移及地向流动增强相对应,而上述现象与驱动入流沿着边界的分布形态（即电场E的分布）有关．     

亚暴事件中磁尾多重等离子体团的数值研究(Ⅱ)

以宁静磁尾平衡解为初态,本文考虑介质的可压缩性,对GEOTAIL卫星于亚暴事件中,观测到多重等离子体团与高能离子爆相对应之特征作模拟研究．数值结果表明：重复形成的等离子体团是高温、高密度区,揭示了等离子体团与高能离子爆──对应的动力学原因．本文的模拟结果还表明：大尺度等离子体团的地向运动,与X中性点的尾向迁移及地向流动增强相对应,而上述现象与驱动入流沿着边界的分布形态（即电场E的分布）有关．     

亚暴事件中磁尾多重等离子体团的数值研究（Ⅰ）

GEOTAIL卫星于1994年1月15日亚暴期间,在深磁尾（x＝96RE）观测到多重等离子体团及与之相对应的高能离子爆,作者以宁静磁尾平衡位形为初态,考虑介质的可压缩性,数值研究亚暴期间磁尾动力学过程．计算结果展现了等离子体团间歇性形成及其运动发展过程．体现了强亚暴事件中储存于碰尾的能量,通过多重等离子体团的排放而逐渐释放的进程．数值结果还表明：持续施加于边界上的晨昏电场及由此引发的驱动重联是导致等离子体团准周期形成的主要因素．此外,作者还考察尾瓣内任一点磁场强度及其分量随时间的演化,它与行进压缩区（TCRs）的观测特征基本相符．     

亚暴事件中磁尾多重等离子体团的数值研究（Ⅰ）

GEOTAIL卫星于1994年1月15日亚暴期间，在深磁尾（x＝96R_E）观测到多重等离子体团及与之相对应的高能离子爆，作者以宁静磁尾平衡位形为初态，考虑介质的可压缩性，数值研究亚暴期间磁尾动力学过程．计算结果展现了等离子体团间歇性形成及其运动发展过程．体现了强亚暴事件中储存于碰尾的能量，通过多重等离子体团的排放而逐渐释放的进程．数值结果还表明：持续施加于边界上的晨昏电场及由此引发的驱动重联是导致等离子体团准周期形成的主要因素．此外，作者还考察尾瓣内任一点磁场强度及其分量随时间的演化，它与行进压缩区（TCRs）的观测特征基本相符．     

亚暴膨胀相近磁尾位形不稳定性模型Ⅱ.中磁尾等离子体流与近磁尾亚暴活动的关系

在单流体MHD近似下,研究了有背景等离子体流存在时的近磁尾位形不稳定性．分析表明,等离子体地向流使两支漂移气球模DBM1和DBM2的不稳定增长更快．与卫星观测资料对比显示,诸多的观测观象与本模型的预言一致．这个结果对中磁尾磁场重联及高速等离子体流与亚暴膨胀相在近地触发之间的关系给出了新的解释．     

亚暴膨胀相近磁尾位形不稳定性模型Ⅱ.中磁尾等离子体流与近磁尾亚暴活动的关系

在单流体MHD近似下，研究了有背景等离子体流存在时的近磁尾位形不稳定性．分析表明，等离子体地向流使两支漂移气球模DBM1和DBM2的不稳定增长更快．与卫星观测资料对比显示，诸多的观测观象与本模型的预言一致．这个结果对中磁尾磁场重联及高速等离子体流与亚暴膨胀相在近地触发之间的关系给出了新的解释．     

亚暴膨胀相近磁尾位形不稳定性模型Ⅱ.中磁尾等离子体流与近磁尾亚暴活动的关系

在单流体MHD近似下，研究了有背景等离子体流存在时的近磁尾位形不稳定性．分析表明，等离子体地向流使两支漂移气球模DBM1和DBM2的不稳定增长更快．与卫星观测资料对比显示，诸多的观测观象与本模型的预言一致．这个结果对中磁尾磁场重联及高速等离子体流与亚暴膨胀相在近地触发之间的关系给出了新的解释．     

亚暴期间磁层和电离层相互耦合的计算

本文假设亚暴时跨尾电流被等离子体团中断,转向成沿磁力线流动的1区场向电流。根据文献[1]的模拟结果,算出场向电流在电离层中的分布。并以此场向电流为输入,计算电离层电势和电流。结果表明,场向电流强度可达1.5×106 A以上,与观测结果一致;它主要分布于70°N以北的区域。电离层电流主要以夜间西向电射流形式出现,主要成份是Hall电流;极区电位降最大可达190kV。本文还简要讨论了等离子体团的晨-昏向尺度对场向电流分布和电离层电导率对电离层电流的影响。     

亚暴期间磁层和电离层相互耦合的计算

本文假设亚暴时跨尾电流被等离子体团中断,转向成沿磁力线流动的1区场向电流。根据文献[1]的模拟结果,算出场向电流在电离层中的分布。并以此场向电流为输入,计算电离层电势和电流。结果表明,场向电流强度可达1.5×10⁶ A以上,与观测结果一致;它主要分布于70°N以北的区域。电离层电流主要以夜间西向电射流形式出现,主要成份是Hall电流;极区电位降最大可达190kV。本文还简要讨论了等离子体团的晨-昏向尺度对场向电流分布和电离层电导率对电离层电流的影响。     

2008年1月5日磁层亚暴期间磁尾近地重联事件的THEMIS观测与分析

本文主要应用THEMIS卫星的磁场和等离子体流观测数据,分析了2008年1月5日08∶51~08∶57 UT亚暴膨胀相期间磁尾的一个近地重联事件.在亚暴膨胀相期间,地面的全天空成像仪清楚地记录到了极光的极向扩展,THEMIS的P5卫星在地球同步轨道附近观测到了磁场的偶极化现象.在亚暴膨胀相末期的08∶51~08∶57 UT期间,P3(X_GSM~-9.12R_E) 和P4 (X_GSM ~-9.40R_E) 同时观测到了一对方向相反的高速等离子体流.这对方向相反的高速等离子体流是由磁尾的重联现象所引起.重联的位置被估计位于X_GSM ~-9.12R_E 和X_GSM~-9.40R_E之间较小的空间范围内.并且,在重联位置的两侧,重联的Hall效应被P3和P4两颗卫星观测到.因此,这一磁尾重联事件发生在距离地球非常近的空间范围内.     

引力场中磁场重联的数值研究（Ⅲ）等离子体团的形成和运动

应用解析和数值相结合的方法求得了含有闭场区、中性片和开场区的二维局域磁静力平衡日冕基态，进而在此基态下数值研究了由电阻撕裂模不稳定性引起的磁场重联过程．结果表明，在电流片长度远大于其半宽度的情况下，仍将发生具有两个Ｘ线的磁场重联，形成磁岛和高温高密度的等离子体团．等离子体团运动的方向取决于闭合区底部等离子体压力和磁压的比值β０．当β０较大时，等离子体团向下运动，引起结合不稳定性和磁岛的合并，等离子体团中的等离子体不断落入闭合磁场区两侧．当β０较小时，等离子体团向上运动，导致电流片中等离子体的喷发．结果还表明，磁张力的分布是决定等离子体团运动方向的主要因子．     

引力场中磁场重联的数值研究（Ⅲ）等离子体团的形成和运动

应用解析和数值相结合的方法求得了含有闭场区、中性片和开场区的二维局域磁静力平衡日冕基态,进而在此基态下数值研究了由电阻撕裂模不稳定性引起的磁场重联过程．结果表明,在电流片长度远大于其半宽度的情况下,仍将发生具有两个Ｘ线的磁场重联,形成磁岛和高温高密度的等离子体团．等离子体团运动的方向取决于闭合区底部等离子体压力和磁压的比值β０．当β０较大时,等离子体团向下运动,引起结合不稳定性和磁岛的合并,等离子体团中的等离子体不断落入闭合磁场区两侧．当β０较小时,等离子体团向上运动,导致电流片中等离子体的喷发．结果还表明,磁张力的分布是决定等离子体团运动方向的主要因子．     

引力场中磁场重联的数值研究（Ⅲ）等离子体团的形成和运动

应用解析和数值相结合的方法求得了含有闭场区、中性片和开场区的二维局域磁静力平衡日冕基态，进而在此基态下数值研究了由电阻撕裂模不稳定性引起的磁场重联过程．结果表明，在电流片长度远大于其半宽度的情况下，仍将发生具有两个Ｘ线的磁场重联，形成磁岛和高温高密度的等离子体团．等离子体团运动的方向取决于闭合区底部等离子体压力和磁压的比值β_０．当β_０较大时，等离子体团向下运动，引起结合不稳定性和磁岛的合并，等离子体团中的等离子体不断落入闭合磁场区两侧．当β_０较小时，等离子体团向上运动，导致电流片中等离子体的喷发．结果还表明，磁张力的分布是决定等离子体团运动方向的主要因子．     

中磁尾重联触发亚暴的单事例分析

磁层亚暴的发生与近磁尾（约6～8 R_E）电流片中断和中磁尾（约20～30 R_E）磁场重联密切相关,而极光的极向扩展、电流片中断和磁尾重联的时序过程对于认识亚暴的触发机制至关重要. 本文利用位于中磁尾的CLUSTER卫星,同步轨道附近LANL-01、LANL-97卫星,近磁尾POLAR和 极区IMAGE卫星的观测,分析了单个亚暴事例.结果表明,在此事件中,中磁尾磁场重联起始比近尾电流片中断早3 min发生,电流片中断发生4 min后,IMAGE卫星观测到极光增亮,同时AE指数突然增大,亚暴膨胀相起始. 观测结果与亚暴中性线模型较为吻合.     

亚暴期间等离子体对流速度及热层纬向风速度变化的统计学分析

本文利用2001年至2005年间CHAMP卫星及DMSP(F13,F15)卫星观测数据,对亚暴发生前后,亚暴初始位置所在磁地方时(Magnetic Local Time, MLT)东侧2 h(+2 h,MLT)至西侧4 h(-4 h,MLT)范围内等离子体对流速度(Vy)及热层纬向风速度(Uy)的变化进行了统计学分析.研究发现在亚暴发生后的1.5 h内,所有MLT区间的Vy均明显增大,且峰值位置的地磁纬度向赤道侧移动,1.5 h后,Vy减小,峰值的纬度向极区移动,表明亚暴的发生能显著增强等离子体对流速度;在亚暴发生位置的西侧0~2 h内,Vy增幅最大,这表明亚暴对热层的影响主要在西侧,影响最大的区域是西侧0~2 h MLT区间;Uy在亚暴发生后3 h持续增大,其对亚暴的响应相较于Vy有1.5 h的延迟.     

磁场偶极化过程中磁尾离子加速:Ⅱ上行离子速度和能量的变化

根据采用动力学方程对亚暴期间磁尾磁场向偶极形弛豫过程中离子分布函数的模拟结果 ,研究了磁尾来自电离层的O+,H+和He+离子的速度及能量随时间的变化 .主要结果为 :（1 ）离子的加速及能量变化主要发生在磁场偶极化过程的中期 ,对应的地心距离位于- 1 2RE到 - 8RE 之间 ;（2 ）垂直于磁场方向上离子加速及能量变化较快 ,平行方向上较慢 ;（3）轻离子较重离子加速及能量变化快 ,磁场偶极化终结 ,3种离子的能量均可增加 2 0 0倍左右 ;（4）初始能量较高时 ,离子加速及能量变化较快 ,离子最终获得的能量较大 .理论计算的磁尾离子能量在磁场偶极化过程终了可达 1 0 2 keV的量级 ,这与观测结果一致 .     

磁尾地向传播等离子体团的模拟研究

卫星观测证实了磁尾等离子体团与亚暴活动的相关性,除了具有北-南双极特征的尾向传播等离子体团外,还发现地向传播等离子体团,它们表现为南-北双极中性片事件和南-北双极瓣区讯号. 资料分析表明:南-北双极讯号的出现几率远低于北-南双极讯号,并且南-北双极事件主要发生于行星际磁场北向和地磁宁静条件,它们往往与小的孤立的地磁亚暴相关. 本文根据地磁宁静时期(IMF Bz北向且By≥Bz)越尾电场Ey分量的分布特点,对地向传播等离子体团作模拟研究. 两类算例的数值结果展示了通量绳磁结构及具有复杂闭合磁力线位形的等离子体团的基本特征,上述特征与尾向传播的等离子体团类似,与IMP 8卫星关于地向传播南-北中性片事件的观测特征大致相符. 数值结果还展示了与Schindler示意图相类似的磁力线拓扑位形,在一定程度上为南-北事件出现几率低作出了解释;并且揭示了磁尾中性片内越尾磁场分量By对磁重联发展的抑制作用. 本文的模拟研究说明:无论磁尾处于活动时期(IMF Bz为南向),还是宁静时期(IMF Bz为北向且By≥Bz),磁场重联均是磁尾等离子体加速和加热的通用机制.     

磁尾地向传播等离子体团的模拟研究

卫星观测证实了磁尾等离子体团与亚暴活动的相关性,除了具有北-南双极特征的尾向传播等离子体团外,还发现地向传播等离子体团,它们表现为南-北双极中性片事件和南-北双极瓣区讯号. 资料分析表明:南-北双极讯号的出现几率远低于北-南双极讯号,并且南-北双极事件主要发生于行星际磁场北向和地磁宁静条件,它们往往与小的孤立的地磁亚暴相关. 本文根据地磁宁静时期(IMF Bz北向且By≥Bz)越尾电场Ey分量的分布特点,对地向传播等离子体团作模拟研究. 两类算例的数值结果展示了通量绳磁结构及具有复杂闭合磁力线位形的等离子体团的基本特征,上述特征与尾向传播的等离子体团类似,与IMP 8卫星关于地向传播南-北中性片事件的观测特征大致相符. 数值结果还展示了与Schindler示意图相类似的磁力线拓扑位形,在一定程度上为南-北事件出现几率低作出了解释;并且揭示了磁尾中性片内越尾磁场分量By对磁重联发展的抑制作用. 本文的模拟研究说明:无论磁尾处于活动时期(IMF Bz为南向),还是宁静时期(IMF Bz为北向且By≥Bz),磁场重联均是磁尾等离子体加速和加热的通用机制.     

磁场偶极化过程中磁尾离子加速:Ⅱ上行离子速度和能量的变化

根据采用动力学方程对亚暴期间磁尾磁场向偶极形弛豫过程中离子分布函数的模拟结果 ,研究了磁尾来自电离层的O+,H+和He+离子的速度及能量随时间的变化 .主要结果为 :（1 ）离子的加速及能量变化主要发生在磁场偶极化过程的中期 ,对应的地心距离位于- 1 2RE到 - 8RE 之间 ;（2 ）垂直于磁场方向上离子加速及能量变化较快 ,平行方向上较慢 ;（3）轻离子较重离子加速及能量变化快 ,磁场偶极化终结 ,3种离子的能量均可增加 2 0 0倍左右 ;（4）初始能量较高时 ,离子加速及能量变化较快 ,离子最终获得的能量较大 .理论计算的磁尾离子能量在磁场偶极化过程终了可达 1 0 2 keV的量级 ,这与观测结果一致 .     

驱动重联过程中等离子体团的聚合

本文对驱动重联过程作了二维可压缩磁流体动力学模拟．对于长的计算域（１：４）,在持续入流的作用下形成了双重磁岛结构,这些磁岛在向外运动的过程中不断靠近,逐渐聚合成为一个大的等离子体团,并继续向外运动．这说明在地球磁层顶和磁尾所观测到的大尺度磁结构,很可能是在驱动重联过程中,由较小的等离子体团聚合后形成的．并初步揭示了在空间等离子体中,磁岛的聚合过程难以直接观测的动力学原因．