行星际磁场北向时磁层顶区磁场重联的全球模式

在对背阳面磁层顶区局域磁场重联模拟的基础上提出了一个行星际磁场北向时磁层顶磁场重联的全球模式。行星际磁场北向时碰层顶磁场重联导致近地尾瓣的能量被输送到远磁尾,太阳风能量不在磁尾储存,向阳面磁层顶变厚,磁层受到一系列扰动。     

基于MHD模拟数据的正午午夜子午面磁层顶位形研究

通过分析太阳风-磁层-电离层系统的三维全球磁流体力学(MHD)模型的计算数据, 给出了正午-午夜子午面磁层顶位形的定量模型. 分析表明, 正午-午夜子午面磁层顶位形可以用文献[3]提出的基于卫星观测数据的、描述赤道面磁层顶位形的函数来描述. 与赤道面磁层顶不同, 正午-午夜子午面磁层顶位形更为复杂. 在忽略极尖区(cusp)的简化条件下, 磁层顶位形仍需利用两条曲线来拟合. 太阳风动压Dp与行星际磁场分量Bz是控制磁层顶位形的主要因素. 行星际磁场为北向时, 磁场增强, 日下点距离r0增大; 行星际磁场为南向时, 磁场增强, 磁层顶日下点距离r0减小. 整体而言, 行星际磁场分量Bz由南转北时, r0增大, 且Bz对r0的影响减弱. 太阳风动压Dp是控制磁层顶日下点的主要因素, Dp增大, r0减小. 磁层顶位形的另一个参数磁层顶磁尾张角α, 随着行星际磁场南向分量增强而增大, 即磁层顶张开程度更加显著, 更多的磁通量由向阳侧传输到夜侧; Dp增大, α略增大, 这意味着Dp对磁通量由日侧向夜侧的传输也有一定的贡献.     

太阳风向磁层输入电磁能量研究

利用全球磁流体力学（MHD）模拟结果,通过确立包含磁层顶的太阳风流线内边界来识别三维磁层顶位形,并以极尖区位置作为磁层顶日侧与夜侧的分界线,在此基础上定量研究了不同条件下穿过磁层顶向磁层内输入的电磁能量. 研究发现,磁层顶的能量传输与太阳风条件密切相关,磁重联是控制电磁能量传输的重要机制. 结果表明,当IMF（行星际磁场）南向时,极尖区后方的磁尾附近存在电磁能输入最大值,当IMF北向时,电磁能输入最大值发生在极尖区附近;南向IMF条件下,在IMF强度增大或太阳风密度增大时,磁层顶电磁能传输的电磁能量比北向IMF条件时增加更显著. 太阳风通过调节磁层顶面积间接影响到磁层顶能量传输大小. 研究还发现,北向IMF与南向IMF条件下穿过磁层顶的电磁能输入的比值范围约为10%～30%,此比值一定程度上反映了北、南方向IMF与地磁场磁重联效率的比值.     

磁层氧离子向磁鞘泄漏的机制和效应

１９９７年 １月 １０日磁暴期间, Ｇｅｏｔａｉｌ卫星在向阳侧的磁鞘中观测到了磁层氧离子突增事件．这些氧离子的出现和磁鞘中存在很强的南向行星际磁场有关．事件期间向阳面发生了准静态的磁重联,氧离子流存在由北向南的速度分量．通量突增过程具有逆向和正向能量色散现象,磁层内部只有氧离子有可能被梯度漂移输送到重联区,所以只有氧离子在磁鞘中持续地被观测到．估计氧离子的逃逸速率为 ０．６１× １０２３／ｓ,大约为环电流氧离子输入率的 ３３％．大量的环电流氧离子由磁层跑到了磁鞘,导致环电流指数 ＡＳＹ－Ｈ呈现明显的非对称性．     

MMS星座对磁层顶磁通量绳内离子惯性尺度结构的观测

利用MMS观测数据,对磁层顶通量绳内离子惯性尺度（d_i）的结构进行分析研究.结果发现,许多不同尺度（约1d_i至数十d_i）的通量绳内都存在具有d_i尺度的电流 j _m,其方向在磁层顶局地坐标系的-M方向,即与磁层顶查普曼-费拉罗电流同向,由电子在+M方向的运动（ v _em）携带.这些电流结构具有以下特征：磁鞘与磁层成分混合,磁场为开放形态;离子去磁化,电子与磁场冻结;N方向（即垂直于磁层顶电流片方向）的电场 E _n显著增大,幅度达到约20mV·m^-1,并伴有明显的尖峰状起伏,该增强和尖峰状起伏的电场对应于霍尔电场.分析表明,电流、电子与离子运动的偏离以及霍尔电场之间遵从广义欧姆定律,三者密切关联.进一步对磁层顶磁重联的探测数据进行分析发现,在很多重联区内也存在与通量绳内相似的结构,其尺度约为d_i量级,其中霍尔电场 E _N、电流 j _M和电子速度 v _eM均与通量绳内对应物理量的方向相同且幅度相近.基于上述观测事实,采用经典FTE通量绳模型,对通量绳内电流、电子运动和霍尔电场的起源进行了初步探讨,认为其来源于磁层顶无碰撞磁重联区内的相应结构,并且后者在离子尺度通量绳的形成过程中起到重要作用.     

磁层顶穿越事件自动识别算法

磁层顶是太阳风与磁层进行质量、动量、能量交换的关键区域.磁层顶穿越事件（MCEs）可通过对卫星探测到的粒子能谱和磁场数据图进行人工分析的方式来识别.因太阳风动压和行星际磁场的易变性,位于磁层顶附近的卫星经过长期观测可能会经历成千上万次的磁层顶穿越.人工分析的方法工作量巨大,而且识别速度慢.本文发展了一种新的日下点附近MCEs自动识别算法.此算法综合分析卫星探测到的粒子和磁场数据,能有效地减少误判的发生.为了验证算法的有效性,采用单CPU计算机对THEMIS卫星在2007—2018年靠近日下点附近观测到的数据进行MCEs自动识别,最终在约6h共识别出16758个MCEs.这些自动识别出来的MCEs样本可用于统计研究磁层顶相关的诸多物理问题,如凹陷磁层顶、太阳风与磁层相互作用,磁层顶磁场重联等.同时还分析了算法的精确性和局限性.     

火星磁层顶形成和变化的理论研究

考虑太阳风动压与行星电离层中的带电粒子热压及磁压之和平衡,建立了有大气（电离层）的行星磁层顶形成的理论模型,结合卫星对火星的观测数据,对子午面内向日侧火星磁层顶位形进行了数值计算和分析,研究了火星磁层顶位形及其与太阳风动压之间的变化关系．结果认为,火星磁层顶位形与地球磁层顶相似．太阳风动压越大,火星磁层顶越靠近火星;太阳风动压越弱,火星磁层顶越远离火星．根据火星内秉磁矩从古到今逐渐减小的观点,探索了大尺度磁场（内禀磁矩）对火星磁层顶的贡献作用,结果认为大尺度磁场越强,火星磁层顶越远离行星．这对于进一步研究火星磁层的长期演化以及其他行星磁层的位形变化都具有重要的意义．     

磁层顶附近流场剪切度的确定

磁层顶附近的流场剪切度与磁层顶附近能量转换的程度有关.很多磁层顶数值模拟用到流场剪切度这个输入参数,但一直是假设的.本文利用Cluster多卫星同时观测数据及独特的时空分辨功能,采用线性插值和重心坐标的方法确定了磁层顶附近流场的剪切度.通过对晨侧和昏侧磁层顶及附近磁层磁鞘流场剪切度的真实空间分布的研究结果表明,在平静的太阳风条件和地磁条件下磁层顶附近流场剪切度有时也很大,可达每百公里330 km/s的相对速度差.但在很多情况下流场是弱剪切的,在上千公里的距离上只有每秒几十公里的相对速度差.本文确定流场剪切度的方法可以推广用来确定任一位置的流场剪切度.     

磁层顶通量传输事件轴向统计分析

磁层顶通量传输事件（Flux Transfer Event,FTE）与磁重联相关,其典型特征为磁场法向分量的双极变化.在不同FTE模型里,FTE结构可能为重联的通量管、由多X线重联形成的闭合磁通量绳或者由单X线重联形成的开放磁场环,从而在磁层顶有不同的整体位形.使用一种新的轴向分析方法,对Cluster在一个向阳面磁层顶穿越季观测到的505个FTE进行统计研究.结果表明：在磁层顶中低纬度的侧翼,大多数FTE轴向均为沿磁层磁力线方向即南北方向,少数FTE轴向沿着不同于磁层磁力线方向的东西方向;在高纬磁层顶,大多数FTE轴向沿东西方向,少数FTE轴向沿着磁层磁力线方向即南北方向.这些统计特征有助于重新认识FTE的全球形态.     

Scientific visualization to study Flux Transfer Events at the Community Coordinated Modeling Center

In this paper we present recent additions to the visualization toolset offered by the Community Coordinated Modeling Center (CCMC). Two suites of visualization tools are available that can address different needs during the analysis of model simulations of the magnetosphere that are provided by the CCMC. The online, server-side visualization allows the user to quickly browse through simulation runs and now can create maps of magnetic field line topology in the magnetosphere. The second tool, SWX, can be used on the client computer after data have been downloaded. With this second tool the user can interact directly with the three-dimensional objects that are being rendered. We present results from a simulation of a Flux Transfer Event that was performed at the CCMC using a magnetohydrodynamic model of the Earth’s magnetosphere with a high resolution grid focused on the dayside magnetosheath and dayside magnetopause. The simulation shows that the FTE that results from localized magnetic reconnection is a complicated three-dimensional structure that requires modern visualization techniques. Visualization techniques that are presented here allow the researcher to fully appreciate the complexity contained in magnetospheric simulation results.     

非对称磁场重联电子扩散区附近的能量通量分析

磁场重联是等离子体中非常重要的能量转换过程,研究该过程中的能量转化形式及其分配十分重要。目前有很多关于对称重联中的能量通量的研究,而磁层顶的非对称重联的研究却很少,特别是多事件的统计研究。因此挑选了MMS卫星观测到的10个经过电子扩散区的磁层顶重联事件对其进行分析,结果表明在不同的事件中能量通量的分配情况不一样,但在大多数事件中,离子焓通量占主导地位,其次是坡印廷通量,离子热通量略小于坡印廷通量,离子动能通量、电子焓通量和电子热通量占比在10%以下。通过对事件进行归一化处理,得到了L方向和M方向的不同能量通量随L方向的磁场和离子速度的关系图,并分析了每种能量通量的分布特征。     

电流片的自发重联与无力场──磁岛的形成

采用二维三分量磁流体力学数值模拟方法,研究了线性无力场对电流片的自发重联的影响。结果表明:电流片的背景场(无力场)对重联的拓扑位形有重要的影响,当无力场因子达到一定值时,在模拟区域内出现大范围的磁岛,磁岛厚度随时间的发展而增加,而无力场因子较小时则没有磁岛形成。当无力场因子为零时,即回到普通的位场情形的单x线重联,位场情形磁场拓扑位形与Scholar的结果一致。     

火星磁场产生的原因及其分布

构建了一个可以得到火星赤道面上磁场分布的模型. 模型根据卫星观测数据, 提出了火星电离层、磁层顶和磁尾电流片上都各自通有电流的假设. 由电流的连续性条件可知, 这三种背景条件下的电流之间满足一定关系, 即火星磁层顶上的总电流是电离层上的总电流与磁尾电流片上的总电流之和. 这些电流产生的磁场与太阳风磁场共同构成了火星赤道面上磁场分布. 通过计算发现, 采用这种磁场模型得到的结果与目前卫星所观测的结果以及与采用其他方法得到的结果符合得较好.     

Hall MHD reconnection with an initial guide field By0

An uniform out-of-plane magnetic field component B_y0 is added to the equilibrium Harris sheet with plasma β = 0.5 and L_c = 0.5d_i (where L_c is the half-width of the equilibrium current layer and d_i is the ion inertial length). Driven by the continuous boundary inflows, the magnetic reconnections with the guide field B_y0/B₀ ranging from 0 to 4.0 are investigated using a 2.5D Hall magnetohydro-dynamic (MHD) code developed from a multi-step implicit scheme. The features of the reconnection field are substantially altered in the presence of the guide field. The openness of the magnetic separatrix angle is slightly reduced and the anti-symmetric quadrupolar structure of B_y field and the symmetric distribution of plasma pressure P are replaced by an asymmetric B_y four-wing structure and an asymmetric P plot as a non-zero B_y0 is added. The decoupling of electrons and ions also occurs near the X line in the case with a finite B_y0, but the effect of initial B_y0 on the electron flow is greater than that on the ion flow. The reconnection rates at the X-line drops from 0.151 to 0.06, namely, ∂A/∂t is reduced by a factor of 2.5 as B_y0/B₀ increases from 0 to 4.0. The reduction of reconnection rate might be related to the reducing openness of reconnection layer with the increasing B_y0.