Study of three dimensional structure of the fast convection flow in the plasma sheet by MHD simulations on the basis of spontaneous fast reconnection model

Three dimensional structure of the fast convection flow in the plasma sheet is examined using magnetohydrodynamic (MHD) simulations on the basis of spontaneous fast reconnection model. The fast flow observed in the near-Earth magnetotail is one of the key phenomena in order to understand the causal relationship between magnetic substorm and magnetic reconnection. In this paper, we focus on this earthward fast flow in the near-Earth magnetotail. Our previous studies have shown that the fast reconnection produces the Alfvénic fast reconnection outflow and drastic magnetic field dipolarization in the finite extent. In this paper, the results of our simulations are compared with those of the in-situ observations in the geomagnetotail. They have consistent temporal profiles of the plasma quantities. It is suggested that the fast convection flows are caused by spontaneous fast reconnection.     

Generation mechanism of the whistler-mode waves in the plasma sheet prior to magnetic reconnection

The whistler-mode waves and electron temperature anisotropy play a key role prior to and during magnetic reconnection. On August 21, 2002, the Cluster spacecrafts encountered a quasi-collisionless magnetic reconnection event when they crossed the plasma sheet. Prior to the southward turning of magnetospheric magnetic field and high speed ion flow, the whistler-mode waves and positive electron temperature anisotropy are simultaneously observed. Theoretic analysis shows that the electrons with positive temperature anisotropy can excite the whistler-mode waves via cyclotron resonances. Using the data of particles and magnetic field, we estimated the whistler-mode wave growth rate and the ratio of whistler-mode growth rate to wave frequency. They are 0.0016f_ce (Electron cyclotron frequency) and 0.0086f_ce, respectively. Therefore the whistler-mode waves can grow quickly in the current sheet. The combined observations of energetic electron beams and waves show that after the southward turning of magnetic field, energetic electrons in the reconnection process are accelerated by the whistler-mode waves.     

Ion dynamics associated with Alfven wave in the near-Earth magnetotail: Two-dimensional global hybrid simulation

Ion dynamics in the near-Earth magnetotail region is examined during periods of fast Earthward flow with a two-dimensional (2-D) global-scale hybrid simulation. The simulation shows that shear Alfven waves are generated at x ∼ −10R_E, where the strong earthward flow is arrested by the dipole field, and propagate along field lines from the equator to both southern and northern polar ionosphere. Non-gyrotropic ion velocity distributions occur where the large-amplitude Alfven waves are dominant. The simulation indicates that the Alfven waves are generated by interaction of the fast earthward flow with the stationary near-Earth plasma. Beam ions are found to be pitch-angle scattered and trapped in the wave field, leading to the non-gyrotropic ion distributions in the high-latitude plasma sheet boundary. In addition, significant particle heating and acceleration are found to occur behind the dipolarization front due to the effect of wave turbulence.     

Study of waves in the magnetotail region with cluster and DSP

The study of the neutral sheet is of fundamental importance in understanding the dynamics of the Earth’s magnetosphere. From the earliest observation of the magnetotail, it has been found that the neutral sheet frequently appears to be in motion due to changing solar wind conditions and geomagnetic activity. Multiple crossings of the neutral sheet by spacecraft have been attributed to a flapping motion of the neutral sheet in the north–south direction, a wavy profile either along the magnetotail or the dawn–dusk direction. Cluster observations have revealed that the flapping motions of the Earth’s magnetotail are of internal origin and that kink-like waves are emitted from the central part of the tail and propagate toward the tail flanks. This flapping motion is shown here to propagate at an angle of ∼45° with x_GSM. A possible assumption that the flapping could be created by a wake travelling away from a fast flow in the current sheet is rejected. Other waves in the magnetotail are found in the ULF range. One conjunction event between Cluster and DoubleStar TC1 is presented where all spacecraft show ULF wave activity at a period of approximately 5 min during fast Earthward flow. These waves are shown to be Kelvin–Helmholtz waves on the boundaries of the flow channel. Calculations show that the conversion of flow energy into magnetic energy through the Kelvin–Helmholtz instability can contribute to a significant part of flow breaking between Cluster and DoubleStar TC1.     

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.     

Plasma flow in the near and distant geomagnetic tail

Measurements of the bulk flow of plasma in the outer magnetosphere were first made a little over a decade ago with Los Alamos instruments on the Vela satellites. During the intervening years, as flow measurements have been made with improved instruments and by other satellites they have come to play a crucial role in the development of our understanding of the structure and dynamics of the magnetosphere. For example, they were the means of discovery of the magnetosphere's boundary layer and of plasma vortices within the plasma sheet. They were the essential ingredient in the identification of signatures of magnetic reconnection at the magnetopause. And they were indispensible in clarifying the complex phenomena in the magnetotail accompanying substorms and in showing that these phenomena are consistent with a substorm model involving magnetic reconnection at a near-earth neutral line. Most recently, magnetotail plasma flow measurements by the ISEE-3 satellite at distances as great as 230 R_E have been instrumental in fixing the average location of the “distant” neutral line at ～ 60 to 120 R_E and in identifying plasmoids (i.e., severed sections of the plasma sheet), released during substorms and escaping down-tail. This paper reviews the features of magnetotail plasma flow, describes the most recent observations, and discusses their implications for magnetospheric physics.     

Evaluation of substorm models with Cluster observations of plasma flow reversal in the magnetotail

We evaluate two prevailing substorm models with an event of plasma flow reversal from tailward to Earthward detected by Cluster at the downstream distance of ∼19 R_E in the magnetotail during a substorm on August 22, 2001. We use the unique capability of Cluster measurements in determining gradients to examine the associated current density, Lorentz force, and current dissipation/dynamo term. In association with plasma flow reversal, it is found that (1) there was no clear quadrupole magnetic perturbation signature, (2) the x-component of the Lorentz force did not change sign, (3) the y-component of the product of the current density and the electric field was occasionally negative indicative of a dynamo effect, and (4) the timing sequence of flow reversal from the Cluster configuration did not match tailward motion of a single plasma flow source. These observations are consistent with the near-Earth initiation model for substorms with multiple current disruption sites moving progressively tailward near the late stage of substorm expansion.     

Numerical study of magnetic reconnection possibly occurring near heliospheric current sheet

The third-order accurate upwind compact difference scheme has been applied for the numerical study of the magnetic reconnection driven by a plasma blob impacting the heliospheric current sheet, under the framework of the two-dimensional compressible magnetohydrodynamics. The results show that the driven reconnection near the current sheet could occur in about 10–30 min for the interplanetary high magnetic Reynolds number, R_M = 2000–10,000, a stable magnetic reconnection structure can be formed in hour order of magnitude, and there appear some basic properties such as the multiple X-line reconnections, vortex structures, filament current systems, splitting and collapse of the high-density plasma blob. These results are helpful in understanding and identifying the magnetic reconnection phenomena possibly occurring near the heliospheric current sheets.     

Dynamics of the plasma sheet in the magnetotail: Interrelation of turbulent flows and thin current sheet structures

Dynamics of the magnetotail involves elementary processes of magnetic field merging (reconnection layer formation) occurring on medium spatial scales. Every such process features two different stages, a fast one and a subsequent slower one. The corresponding short time scale _T1$T_1$ is associated with disturbances propagating in the tail lobes. The longer time scale _T2$T_2$ is associated with plasma motions in the plasma sheet. A disturbance appearing in the magnetotail on the time scale _T1$T_1$ results in a loss of equilibrium in the plasma sheet. By means of theoretical argument and numerical simulation, it is shown that the relaxation process which follows on the time scale _T2$T_2$, produces extremely thin embedded current sheets, along with generation of fast plasma flows. The process provides an effective mechanism for transformation of magnetic energy accumulated in the magnetotail, into energy of plasma flows. The fast flows may drive turbulent motions on shorter spatial scales. In their turn, those motions can locally produce very thin current sheets; after that, nonlinear tearing process leads to generation of neutral lines, and reconnection. The latter produces new fast disturbances on the time scale _T1$T_1$ closing the feedback loop.     

Electron temperature anisotropy effects on tearing mode in ion-scale current sheets

Recent two-dimensional (2-D) particle-in-cell (PIC) simulations have shown that there is a critical thickness of a current sheet, above which no significant saturation amplitude of the 2-D tearing (TI) mode can be expected. Here, we have introduced the initial electron temperature anisotropy (α_e0 = T_e⊥/T_e|| > 1), which is known to raise significantly the linear growth rates, and inspected if α_e0 > 1 can change the saturation level of the TI in a super-critical current sheet. Varying α_e0 and D (D: the current sheet half-thickness) systematically, we have found that while α_e0 boosts up the linear growth rate in both sub- and super-critical current sheets, macroscopic effects are obtained only in sub-critical current sheets, that is, energy transfer from the fastest growing short wavelength modes to longer wavelength modes are available only in the sub-critical regime. Since the critical thickness is a fraction of the ion inertial length, the tearing mode assisted by the electron temperature anisotropy alone, despite its significant boost in the linear growth rate, cannot be the agent for reconnection triggering in a current sheet of ion-scale thickness.     

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.     

等离子体彗尾动力学的数值研究

本文基于可压缩磁流体动力学模型,数值研究了尾瓣巾具有超Alfven速流动的等离子体彗尾的动力学特征。结果表明,等离子体片和尾瓣之间的剪切等离子体流动将会激发流动撕裂模不稳定性,引起彗尾等离子体片中发生磁场重联,形成磁岛和高密度的等离子体团。进而模拟了太阳风引起的局部驱动力对等离子体彗尾中磁场重联的影响,其特征时间远大于流动撕裂模。我们认为一些观测到的等离子体彗尾中的四块和彗尾截断事件可能主要与彗尾中剪切等离子体流动所引起的流动撕裂模不稳定性有关。      

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

A fast tearing mode instability driven by agyrotropic electron pressure

The collisionless plasma environment at the current sheet of the Earth’s magnetotail is subjected to fast dynamic evolutions such as tearing instability. By considering agyrotropic pressure for electron and ion components of a collisionless plasma, we analytically investigate the dynamics of tearing mode instability, in which, breaking the frozen-in condition can either be provided by the electron inertia or by agyrotropic electron pressure. A set of linearized Hall-Magnetohydrodynamic (MHD) equations describes the evolution of tearing mode in a sheared force-free field. The presented scaling analysis shows that if the plasma-β$\beta$ exceeds a specified value, then the main mechanism of magnetic reconnection process is the nongyrotropic electron pressure. In this regime, the role played by agyrotropic ion pressure inside the reconnection layer is out of significance. Therefore, the electron-MHD framework, adequately, describes the dynamics of tearing instability with a growth rate which is much faster compared to the cases with a dominated bulk inertia or a gyrotropic plasma pressure.     

双带耀斑磁重联过程的数值研究

本文用多步隐格式求解包含电阻的磁流体力学方程组, 对双带耀斑的主相作数值模拟, 清晰地展示了中性片区由撕裂模线性重联向准稳态重联的过渡以及后随耀斑环的产生和等离子体团的喷发过程.对于在能量方程中计及和忽略焦耳加热两种情况, 分别作了计算.结果表明, 计及焦耳加热时, 电流片中等离子体的温度显著增加(是初始温度的2—3倍), 但等离子体的运动速度却变化不大.两种情况的计算结果均表明:等离子体的运动速度低于声速, 因此不会形成快激波.计及焦耳加热的计算结果显示了两个新的特征:其一是中性片高密度等离子体的受热膨胀, 增大了电流片的有效厚度, 它使重联速率降低, 并逐渐趋于饱和, 其二是同时形成上升和沉降等离子体团, 后者与耀斑环碰撞, 并合并于后随耀斑环内.      

The plasma sheet and boundary layers under northward IMF: A multi-point and multi-instrument perspective

During conditions of northward interplanetary magnetic field (IMF), the near-tail plasma sheet is known to become denser and cooler, and is described as the cold-dense plasma sheet (CDPS). While its source is likely the solar wind, the prominent penetration mechanisms are less clear. The two main candidates are solar wind direct capture via double high-latitude reconnection on the dayside and Kelvin–Helmholtz/diffusive processes at the flank magnetopause. This paper presents a case study on the formation of the CDPS utilizing a wide variety of space- and ground-based observations, but primarily from the Double Star and Polar spacecraft on December 5th, 2004. The pertinent observations can be summarized as follows: TC-1 observes quasi-periodic (∼2 min period) cold-dense boundary layer (compared to a hot-tenuous plasma sheet) signatures interspersed with magnetosheath plasma at the dusk flank magnetopause near the dawn-dusk terminator. Analysis of this region suggests the boundary to be Kelvin–Helmholtz unstable and that plasma transport is ongoing across the boundary. At the same time, IMAGE spacecraft and ground based SuperDARN measurements provide evidence of high-latitude reconnection in both hemispheres. The Polar spacecraft, located in the southern hemisphere afternoon sector, sunward of TC-1, observes a persistent boundary layer with no obvious signature of boundary waves. The plasma is of a similar appearance to that observed by TC-1 inside the boundary layer further down the dusk flank, and by TC-2 in the near-Earth magnetotail. We present comparisons of electron phase space distributions between the spacecraft. Although the dayside boundary layer at Polar is most likely formed via double high-altitude reconnection, and is somewhat comparable to the flank boundary layer at Double Star, some differences argue in favour of additional transport that augment solar wind plasma entry into the tail regions.     

Acceleration theory for 5–40 keV ions at interplanetary shocks

Power-law spectra f(E)∝E^?2.7 of < 40 keV suprathermal ions within ～10⁷ km of propagating interplanetary shocks are explained by diffusive scattering near a plane shock. The theory fits the 25 November 1977 event with a mean free path perpendicular to the shock with is 0.01 AU in front of the shock and less than .0003 AU behind it, for 1 keV ions. The theory predicts a steepening spectrum at higher energies, of the form $\text{f(v)∝v}{}^{\mathrm{?4}}\text{exp(??λdv/ur)}\text{where}\text{u = (ΔV)}{}^2\text{/2V}{}_{\mathrm{W}}$ depends on the plasma velocity jump ΔV and the plasma speed V_W and mean free path λ in front of the shock     

On the large-scale structure of the tail current as measured by THEMIS

The magnetic field structure and the spatial characteristics of the large-scale currents in the magnetospheric tail were studied during quiet and moderately disturbed geomagnetic conditions in 2009. The magnetic field of the currents other than the tail current was calculated in terms of a paraboloid model of the Earth’s magnetosphere, A2000, and was subtracted from measurements. It was found on the base of obtained tail current magnetic field radial distribution that the inner edge of the tail current sheet is located in the night side magnetosphere, at distances of about 10 R_E and of about 7 R_E during quiet and disturbed periods respectively. During the disturbance of February 14, 2009 (Dst_min ∼ −35 nT), the B_x and the B_z component of the tail current magnetic field near its inner edge were about 60 nT, and −60 nT that means that strong cross-tail current have been developed. The tail current parameters at different time moments during February 14, 2009 have been estimated. Solar wind conditions during this event were consistent with those during moderate magnetic storms with minimum Dst of about −100 nT. However, the magnetospheric current systems (magnetopause and cross-tail currents) were located at larger geocentric distances than typical during the 2009 extremely quiet epoch and did not provide the expected Dst magnitude. Very small disturbance on the Earth’s surface was detected consistent with an “inflated” magnetosphere.     

Determination of reconnected flux via remote sensing

Magnetic reconnection is one of the most fundamental processes in the magnetosphere. We present here a simple method to determine the essential parameters of reconnection such as reconnected flux and location of the reconnection site out of single spacecraft data via remote sensing. On the basis of a time-dependent reconnection model, the dependence of the reconnected flux on the magnetic field z-component B_z is shown. The integral of B_z over time is proportional to the reconnected flux and depends on the distance between the reconnection site and the actual position where B_z is measured. This distance can be estimated from analysis of magnetic field B_z data. We apply our method to Cluster measurements in the Earth’s magnetotail.