IMF北向时太阳风粒子向磁层输运的试验粒子模拟研究.docxVIP

? ? IMF北向时太阳风粒子向磁层输运的试验粒子模拟研究* ? ? 赵明现 (中国气象局空间天气重点开放实验室/国家卫星气象中心(国家空间天气监测预警中心) 北京 100081) (许健民气象卫星创新中心 北京 100081) 0 引言 太阳风和行星际磁场（IMF）与地球内禀磁场相互作用，形成了向阳面被压缩、背阳面被拉伸的地球磁层。磁层中的等离子体主要有内外两个来源，内源为电离层，外源为太阳风[1]。一般来说，太阳风等离子体对磁层中的远地区域起支配作用，而来自电离层的粒子对磁层中的近地区域起主导作用[2,3]。 太阳风等离子体向磁层输运是太阳风-磁层相互作用的一个重要过程，对这个过程物理机制的探讨一直是一个重要的研究方向。Dungey[4]提出的磁场重联机制是目前普遍接受的一种理论，在向阳面磁层顶日下点发生磁场重联，从而剥离向阳面磁层中的磁力线，并将之拖向背阳面磁尾。日下点的磁场重联受IMF 控制，IMF 南向时的重联率要显著的高于IMF北向时[5]，此时太阳风粒子主要通过晨侧磁层顶直接进入磁层中的等离子体片和环电流区域[2]。在IMF北向期间，磁重联发生在高纬的南北极尖区背阳面一侧[6,7]，这个过程会在向阳面磁层中形成新的闭合磁力线，闭合磁力线携带太阳风粒子沉入磁层内部，并沿方位角方向往晨昏两侧形成低纬边界层（Low Latitude Boundary Layer，LLBL）。此 外，IMFBx和By在太阳风-磁层能量耦合过程中也起着重要的作用，同样条件下IMF 径向时进入磁层的能量比IMF北向时要大[8]。 Axford 和Hines[9]提出，太阳风与磁层顶侧面等离子体的黏性相互作用会在磁层中形成对流，从而使太阳风动量传入磁层内。这种机制不受IMF 方向的影响，不过对太阳风向磁层的能量传输贡献并不大，大概在10%～20%[10]。K-H 不稳定（Kelvin-Helmholtz Instability,KHI）也是一种可能的传输机制。观测表明，在低纬的晨昏两侧，太阳风等离子体相对于磁层顶内的等离子体具有很高的速度，K-H 波在边界层会形成局地的小尺度涡旋[11,12]。虽然具体的物理过程尚无定论，但是从卫星观测数据分析，在IMF 北向时，磁尾晨昏两侧的KHI 会导致磁鞘中的太阳风向磁层内的传输[13]。另外，磁鞘中的动力阿尔芬波也能让磁鞘等离子体进入磁层[14,15]，如果磁鞘中存在动量足够大的不均匀结构体，也可以直接穿过磁层顶进入磁层[16,17]。不同的传输机制之间并不互斥，在实际的物理过程中，可能同时起作用[18]。 试验粒子方法已经在空间物理和天文学中得到广泛的应用，常用于日地空间中粒子的起源、传播和损耗等过程的研究。试验粒子法需要给定背景场，包括电场和磁场等，通过跟踪单个或多个粒子的动力学过程，来研究复杂系统中的物理过程。由于不考虑试验粒子对背景场的反馈以及各试验粒子之间的相互作用，所以这个方法是不自洽的。但是在研究大尺度的复杂问题时，自洽模型会带来巨大的运算量，而作为系统一阶近似的试验粒子方法足够研究很多粒子的动力学问题，是大尺度磁流体模拟的有益补充[19]。也可以在计算试验粒子轨道时，使用时间回推方法，研究粒子起源[20]。 试验粒子法又可以分为粒子轨道法、蒙特卡罗法等[19]，粒子轨道法通过单个粒子或一群粒子的运动轨道，来研究系统中的粒子输运过程，本文研究使用的是试验粒子轨道法。Speiser[21]使用试验粒子轨道法基于Dungey 磁层模型研究了太阳风粒子向磁层电流片中的输运和极光粒子的加速，但其使用的磁层背景场设定太过简化，与真实物理过程相去甚远。Delcourt[22]通过单粒子的轨道变化研究了磁层中氧离子在地磁亚暴期间的加速过程。Li等[23,24]研究亚暴期间能量粒子注入磁层的过程，Tang等[25]对磁暴期间能量粒子进入磁层做了数值模拟研究，这些研究着重于能量粒子，而不是能量较低的太阳风粒子。Moore等[2]使用四阶龙格库塔方法求解全洛伦兹方程，以LFM 磁层模型的磁流体（MHD）模拟结果为背景场，研究了太阳风和极区电离层向磁层中的粒子输运特征，不过其所使用的MHD 背景场是选择了IMF 南向和北向两个时刻，是固定的背景场，此外该MHD模拟并未考虑地磁轴倾角的影响。Ding[3,26]使用试验粒子法模拟计算了太阳风粒子向磁层等离子体片的输运，背景场使用BATS-R-US 的MHD 模拟结果，背景场随时间变化且考虑了地磁倾角，不过MHD 模拟使用的是人工设定的理想太阳风条件，并非实际观测得到的太阳风数据，另外粒子轨道的计算采用了引导中心近似，而引导中心近似方法得到的粒子运动轨迹与全洛伦兹方程相比会存在一定的误差[3]。Guo等[27]基于自主开发的磁层MHD 模型和试验粒子法，研