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

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

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

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

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

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

2004年3月13日TC-1和Cluster卫星联合观测的磁通量传输事件

分析了2004年3月13日12:15到12:25UT期间TC-1和Cluster卫星簇的磁通门磁力计(FGM)和电子/电流试验仪(PEACE)的联合观测数据.在此期间,TC-1卫星位于日下点以南的磁层顶附近的磁鞘中,并在12:19UT左右观测到了一个典型的先正后负的磁鞘磁通量传输事件(FTE);而Cluster卫星簇位于北半球日侧高纬磁层项附近,并于12:23UT左右穿出磁层顶进入磁鞘,且在12:21 UT左右也观测到了一个典型的先正后负的磁层FTE.比较分析发现此两个FTE具有类似的磁场结构和等离子体特征,可能是同一个北向运动的FTE先后被TC-1和Cluster卫星观测到.利用Cluster 4颗卫星的多点同时观测数据,采用最小方向微分法和时空微分方法,推断Cluster卫星观测的这个FTE是尺度大小约为1.21Re的准二维结构,其运动方向为东北方向,与Cooling模型预测方向基本一致.利用Cooling模型的预测,推算了TC-1卫星在12:19UT观测的FTE的运动速度和尺度,进而得出随着通量管的极向运动,其速度和尺度均有所增加.     

远磁尾磁层顶位形的偏转

利用WIND和ARTEMIS卫星观测数据，分析远磁尾磁层顶对行星际和太阳风变化的响应，尤其是偏离日地连线的太阳风速度改变对远磁尾磁层顶的影响.研究发现在2011年9月13日的事件中，P2卫星观测到高速且高密度的磁鞘流.利用最小变量法进行分析发现，磁层顶沿着偏离日地连线的太阳风速度方向发生偏转.根据相似三角形定理，推断出本次事件中磁层顶在y方向和z方向上的偏转幅度分别达到10R_e和6R_e.P1和P2卫星的相对位置也证实了这一观点.因此，偏离日地连线的太阳风速度对远磁尾磁层顶的位形影响很大.研究结果可为建立包含太阳风速度v_y和v_z效应的磁层顶模型提供观测证据.      

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

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

基于能量色散特征计算磁层顶重联位置的方法研究

磁层顶磁场重联是太阳风向磁层输入能量的主要方式.重联如何触发一直是空间物理研究的难点,其机制仍然有待深入研究.由于卫星穿越磁层顶时,很难恰好穿越重联发生的区域,因此难以观测到重联的触发条件.本文利用THEMIS卫星观测,确立了反演磁层顶重联点的方法.当重联刚开始发生时,卫星能够观测到离子的能量色散特征,可利用其计算卫星到重联发生位置的距离.沿着磁力线模型追踪该距离即可反演出磁层顶发生重联的位置.与其他方法进行了对比分析,结果显示本文方法比其他方法具有更高的精度.      

2.5维自适应磁场重联MHD模式

磁雷诺数(R_m)是影响磁场重联的重要因素. 真实的物理环境中R_m往往很高, 例如, 在行星际空间和太阳日冕中R_m通常大于104量级. 高R_m条件下的磁重联表现出很多异常特性, 然而高R_m条件下的磁场重联数值模拟需要很高的时空分辨率, 否则很难分辨出重联过程中形成的薄电流片. 本文基于自适应软件包PARAMESH将并行自适应网格技术引入磁场重联数值模拟, 建立了一个2.5维自适应磁场重联MHD模式, 研究高磁雷诺数条件下重联的动态演化过程, 进而将不同磁雷诺数的参数进行对比研究. 结果表明, 该模式可以自动捕捉到磁场重联产生的奇性电流片, 高磁雷诺数条件下产生的慢激波结构可提供一种快速磁能释放机制.      

TC-1和Cluster对向阳侧磁层顶通量传输事件的联合观测研究

2004年2至4月期间,探测一号(TC-1)卫星和Cluster卫星有25次同时处在向阳侧磁层顶附近的磁鞘内,TC-1卫星在低纬区,Cluster卫星在中高纬区.利用这一期间两卫星探测到的27个通量传输事件(FTE),分析行星际磁场(IMF)横向分量BT={By,Bz}对磁层顶重联发生位置的影响,以及分量重联的观测事实,得到如下主要结果.(1)当IMF南向分量Bz占优势(|Bz|＞|By|)时,FTE大多(约占87.5%)能在低纬观测到,而当IMF By分量占优势(|Bz|＜|By|)时,则FTE大部分能在中高纬观测到(占84.2%);(2)很少观测到相关联的事件(关联事件指在低纬生成的FTE,向高纬运动中先后被TC-1卫星和Cluster卫星探测到的事件),表明在低纬形成的FTE可能大多沿磁层顶两侧滑向磁尾,只有少数可能运动到高纬地区;(3)中纬地区探测到的FTE大多是以分量重联方式产生于该区,而非来自磁赤道附近成对形成的FTE.      

环电流区能量离子分布特性研究

为了考察环电流区离子的分布情况,采用环电流粒子理论模式,对环电流中10-100 keV的离子进行了模拟研究.这个模式能够根据近地注入区外边界处离子的分布函数得出磁暴主相期间环电流中的主要成分H+,O+,He+3种离子的通量分布.计算结果分析表明,在其他条件相同的情况下,不同种类离子的通量分布的形态结构十分相似.电场强度对环电流离子通量的空间分布具有决定性的作用;晨昏电场强度越强,离子的通量越高;晨昏电场越强,环电流离子的内边界越接近地球.10keV的离子在电场相当弱的时候还是存在着连续的通量分布,但他们的形态和结构随着电场的变化有明显的变化.电场很弱时,离子分布主要集中于内外两个环带,离子通量在晨侧的更多一些,离子通量的最大值基本上是在比较靠近地球的环带上;随着电场的增强,离子分布的内外两个环带逐步合并,离子的分布逐渐靠近地球,通量分布的最大值也移动到了昏侧.环电流离子投掷角分布具有各向异性,投掷角在90°左右的时候,离子通量能达到最大值.      

行星际磁场时钟角与地球磁层开放磁通的关系

地球磁层开放磁通F_pc是研究磁层动力学过程的重要参数之一,其与日侧和夜侧磁尾的磁场重联具有密切关系. 日侧重联率控制稳定状态下磁层开放磁通的大小,主要受各种太阳风条件的影响. 其中,行星际磁场（IMF）的时钟角是影响日侧重联率的一个重要因素. 通过全球MHD模拟,研究了行星际磁场时钟角θ_c与地球磁层开放磁通F_pc 之间的关系. 结果表明,开放磁通F_pc随着行星际磁场时钟角 θ_c逐渐接近180°（纯南向）而逐渐增加,两者之间的关系近似为F_pc∝sin3/2（θ_c/2）. 由于表征行星际磁场与地球磁场剪切程度的θ_c影响日侧重联率,从而控制F_pc,该关系反映了二者之间的物理联系.      

利用Cluster卫星磁尾单事例研究具有复杂核心磁场结构的磁通量绳

2002年8月28日09:50UT-10:50UT,Cluster卫星在地球磁尾观测到一次导向场磁场重联事件.卫星观测到磁场重联扩散区附近清晰的霍尔（Hall）四极型磁结构.由于导向场的存在,该四极型结构被扭曲变形.在该磁场重联事件中,卫星观测到多个磁通量绳,大部分磁通量绳的核心场极性与导向场极性一致.但是,其中一例磁通量绳的核心场结构极性较复杂.该例磁通量绳中心区域核心场强度出现峰值,核心场极性和导向场极性一致;中心以外区域的核心场极性和导向场极性相反.这种复杂核心场结构以前未见报道.通过最小方向导数法,发现该磁通量绳的轴向是弯曲的.C1和C3卫星穿越了磁通量绳弯曲部分,探测到核心场极性变化;C2和C4卫星位于C1和C3卫星的北侧,仅穿越了磁通量绳弯曲处的一部分,故核心场具有单极性.      

Theoretical model of steady-state magnetic reconnection in collisionless incompressible plasma based on the Grad–Shafranov equation solution

The problem of steady-state magnetic reconnection in an infinite current layer in collisionless, incompressible, nonresistive plasma, except of the electron diffusion region, is examined analytically using the electron Hall magnetohydrodynamics approach. It is found that this approach allows reducing the problem to the magnetic field potential finding, while last one has to satisfy the Grad–Shafranov equation. The obtained solution demonstrates all essential Hall reconnection features, namely proton acceleration up to Alfvén velocities, the forming of Hall current systems and the magnetic field structure expected. It turns out that the necessary condition of steady-state reconnection to exist is an electric field potential jump across the electron diffusion region and the separatrices. Besides, the powerful mechanism of electron acceleration in X-line direction is required. It must accelerate electrons up to the electron Alfvén velocity inside the diffusion region and on the separatrixes. This is a necessary condition for steady-state reconnection as well.     

Intrinsic time scale for reconnection on the dayside magnetopause

Recently much attention has been focused on the transient behavior of the magnetopause in response to pressure pulses and southward fluctuations of the interplanetary magnetic field. We examine the motion of the magnetopause behind the foreshock and conclude that this motion is affected by foreshock pressure variations but not by fluctuations in the direction of the magnetic field. Neither magnetopause erosion nor flux transfer event occurrence is controlled by the foreshock. On the contrary, flux transfer events occur at times of steady IMF and thier quasi-periodic behavior is controlled by the magnetopause or the magnetosphere and is not driven by the external boundary conditions. Since flux transfer events are clearly due to reconnection, this observation implies that the IMF must be southward some time perhaps as long as 7 minutes before flux transfer begins.     

Interconnection of high-latitude and low-latitude boundary layers when IMF BY is dominant

The antiparallel merging model places the location of the reconnection region for a dominant interplanetary magnetic field (IMF) B_Y at high latitudes at the dayside magnetopause and predicts that the low-latitude boundary layer (LLBL) is located on open field lines of the magnetospheric flanks. Interball-1 data obtained in the wide local time range near the low-latitude magnetopause makes it possible to analyze the LLBL plasma population and to find a link between possible reconnection at high latitudes and LLBL occurrence. We found that no boundary layer was observed in the regions which have no topological connection with the merging site. All cases of LLBL observations are located downstream from a specific boundary. This boundary coincides with the first magnetospheric field line touching the reconnection region and can be located in a wide local time region depending on the instant IMF direction. Even the LLBL on closed field lines shows the tendency to be concentrated in the vicinity of this boundary. Thus we show that all types of observed LLBLs are linked to reconnection sites predicted by the antiparallel merging model.     

Two-spacecraft observations of reconnection at the magnetopause: Model results and data comparison

We revisit an example of “quasi-steady” magnetic reconnection at the dayside magnetopause on February 11, 1998, observed by Equator-S and Geotail at the dawnside magnetopause. Phan et al. [Phan, T.D. et al., 2000. Extended magnetic reconnection at the Earth’s magnetopause from detection of bi-directional jets. Nature 404, 848–850.] reported oppositely directed jets at these spacecrafts and inferred a length of the reconnection line of about 38R_E. Pinnock et al. [Pinnock, M., Chisham, G., Coleman, I.J., Freeman, M.P., Hairston, M., Villain, J.-P., 2003. The location and rate of dayside reconnection during an interval of southward interplanetary magnetic field. Ann. Geophys. 21, 1467–1482.] used measurements from SuperDARN radars to show that the reconnection electric field was variable. Here we complement this work by obtaining snapshots of the reconnection electric field from the in situ observations. To do this, we apply a reconstruction method based on a model of compressible Petschek-type magnetic reconnection. This independent method uses magnetic field observations as input data to calculate the reconnection electric field. We obtain average values of E_rec in the range of 0.4–2.4 mV/m. Further we infer a distance perpendicular to the reconnection line of 0.4–0.6R_E. The model results are compared with the two studies mentioned above. It thus appears that while the transfer of momentum for this event is indeed large-scale, the actual rate depends on the time it is measured.     

Analytical MHD model of three-dimensional time-dependent reconnection in a compressible plasma

We describe and analyse a model of three-dimensional time-varying reconnection in which the effect of surface waves is neglected. Reconnection is assumed to be caused by a localized decrease of the plasma conductivity inside the diffusion region. The localized dissipation gives rise to an electric field E^*, which determines the reconnection rate. As a result, the current sheet decays into a system of large-amplitude MHD waves, which propagate along the current sheet and thereby induce perturbations in the surrounding medium. Our model is applied to the case of reconnection at the magnetopause, in particular to the signatures referred to as FTEs (flux transfer events).     

Global Hybrid Simulation of Magnetic Reconnection in the Magnetosheath

A three-dimensional (3-D) global hybrid simulation is carried out for the generation and structure of magnetic reconnection in the magnetosheath due to interaction of an interplanetary Tangential Discontinuity (TD) with the bow shock and magnetosphere. Runs are performed for solar wind TDs possessing different initial half-widths. As the TD propagates through the bow shock toward the magnetopause, it is greatly narrowed by a two-step compression processes, a "shock compression' followed by a subsequent ``convective compression'. In cases with a relatively thin solar wind TD, 3-D patchy reconnection is initiated in the transmitted TD, forming magnetosheath flux ropes. Multiple components of ion particles are present in the velocity distribution in the magnetosheath merging, accompanied by ion heating. For cases with a relatively wide initial TD, a dominant single X-line appears in the subsolar magnetosheath after the transmitted TD is narrowed. A shock analysis is performed for the detailed structure of magnetic reconnection in the magnetosheath. Rotational Discontinuity (RD)/Time-Dependent Intermediate Shock (TDIS) are found to dominate the reconnection layer, which and some weak slow shocks are responsible for the ion heating and acceleration.