Statistical and case studies, as well as data‐mining reconstructions suggest that the magnetotail current in the substorm growth phase has a multiscale structure with a thin ion‐scale current sheet ...embedded into a much thicker sheet. This multiscale structure may be critically important for the tail stability and onset conditions for magnetospheric substorms. The observed thin current sheets are found to be too long to be explained by the models with isotropic plasmas. At the same time, plasma observations reveal only weak field‐aligned anisotropy of the ion species, whereas the anisotropic electron contribution is insufficient to explain the force balance discrepancy. Here we elaborate a self‐consistent equilibrium theory of multiscale current sheets, which differs from conventional isotropic models by weak ion anisotropy outside the sheet and agyrotropy caused by quasi‐adiabatic ion orbits inside the sheet. It is shown that, in spite of weak anisotropy, the current density perturbation may be quite strong and localized on the scale of the figure‐of‐eight ion orbits. The magnetic field, current and plasma density in the limit of weak field‐aligned ion anisotropy and strong current sheet embedding, when the ion scale thin current sheet is nested in a much thicker Harris‐like current sheet, are investigated and presented in an analytical form making it possible to describe the multiscale equilibrium in sharply stretched 2D magnetic field configurations and to use it in kinetic simulations and stability analysis.
Plain Language Summary
Conventional kinetic equilibria with isotropic pressures for ions and electrons aimed to describe the current sheet in Earth's magnetotail cannot reproduce its multiscale structure with the proton gyroradius‐scale current sheet being embedded into a much thicker sheet. They cannot explain either the formation of such thin current sheets sufficiently far from Earth. The embedding effect can be reproduced in case of anisotropic and agyrotropic plasmas because orbits of weakly magnetized ions near the current sheet deviate from the Larmor circle and become more like a figure of eight. However, the corresponding multiscale current sheet models have been studied so far for substantial and strong plasma anisotropy, while observations suggest that the tail plasmas are weakly anisotropic. Here we perform an analysis of a weakly anisotropic current sheet model, which transforms in the isotropic limit into a classical Harris sheet model, and show that the key observed embedding features can be reproduced.
Key Points
Kinetic equilibria of ion‐scale current sheets embedded into a thicker weakly anisotropic Harris‐like current sheet are investigated
The current density increase due to quasi‐adiabatic ion motions may be substantial in spite of weak plasma anisotropy
2D thin current sheets have aspect ratios consistent with observations and controlled by the embedding strength
Onset of reconnection in the tail requires the current sheet thickness to be of the order of the ion thermal gyroradius or smaller. However, existing isotropic plasma models cannot explain the ...formation of such thin sheets at distances where the X‐lines are typically observed. Here we reproduce such thin and long sheets in particle‐in‐cell simulations using a new model of their equilibria with weakly anisotropic ion species assuming quasi‐adiabatic ion dynamics, which substantially modifies the current density. It is found that anisotropy/agyrotropy contributions to the force balance in such equilibria are comparable to the pressure gradient in spite of weak ion anisotropy. New equilibria whose current distributions are substantially overstretched compared to the magnetic field lines are found to be stable in spite of the fact that they are substantially longer than isotropic sheets with similar thickness.
Plain Language Summary
Ion scale current sheets forming sufficiently far from Earth are necessary to explain its stretched magnetic field reconfiguration on the night side. However, isotropic plasmas form magnetic fields that inflate with distance from Earth and cannot reproduce the observed stretched geometry. We present kinetic simulations of current sheets that inflate more gradually with distance due to slight field‐aligned anisotropy of the ion species. Their formation is provided by a special population of suprathermal ions with figure‐of‐eight orbits. We find that the resulting current sheets are stable over a long time scale and have a thickness comparable to the size of these orbits.
Key Points
Two‐dimensional ion‐scale current sheets stretched way beyond the isotropic limit are reproduced in particle‐in‐cell simulations
Weak ion anisotropy and agyrotropy substantially modify the current density and the isotropic force balance
Ion‐scale current sheets are stable in spite of the fact that they are longer compared to isotropic sheets with similar thickness
Mining of substorm magnetic field data reveals the formation of two X‐lines preceded by the flux accumulation at the tailward end of a thin current sheet (TCS). Three‐dimensional particle‐in‐cell ...simulations guided by these pre‐onset reconnection features are performed, taking also into account weak external driving, negative charging of TCS and domination of electrons as current carriers. Simulations reveal an interesting multiscale picture. On the global scale, they show the formation of two X‐lines, with stronger magnetic field variations and inhomogeneous electric fields found closer to Earth. The X‐line appearance is preceded by the formation of two diverging electron outflow regions embedded into a single diverging ion outflow pattern and transforming into faster electron‐scale reconnection jets after the onset. Distributions of the agyrotropy parameters suggest that reconnection is provided by ion and then electron demagnetization. The bulk flow and agyrotropy distributions are consistent with MMS observations.
Plain Language Summary
The process of sudden changes of the magnetic field topology on the night side of the magnetosphere, a global magnetic bubble shielding our planet from the flow of solar wind particles, is simulated for the first time using several important observational constraints. They include the magnetic flux accumulation, formation of a thin current layer with the thickness comparable to the ion gyroradius, solar wind driving electric field, negative charging of the layer and domination of electrons as current carriers. Simulations resolve ion and electron motions beyond their fluid approximation and reveal interesting embedded structure of electron and ion watersheds preceding the topology change.
Key Points
Onset of reconnection is reproduced in three‐dimensional particle‐in‐cell simulations guided by magnetometer data mining and local plasma observations
It is preceded by the formation of two diverging electron outflow regions embedded into a single diverging ion outflow pattern
Reconnection is provided by ion and then electron demagnetization consistent with MMS data
This work builds on and extends our previous effort (Tsyganenko et al., 2003) to develop a dynamical model of the storm‐time geomagnetic field in the inner magnetosphere, using space magnetometer ...data taken during 37 major events in 1996–2000 and concurrent observations of the solar wind and interplanetary magnetic field (IMF). The essence of the approach is to derive from the data the temporal variation of all major current systems contributing to the distant geomagnetic field during the entire storm cycle, using a simple model of their growth and decay. Each principal source of the external magnetic field (magnetopause, cross‐tail current sheet, axisymmetric and partial ring currents, and Birkeland current systems) is driven by a separate variable, calculated as a time integral of a combination of geoeffective parameters NλVβBsγ, where N, V, and Bs are the solar wind density, speed, and the magnitude of the southward component of the IMF, respectively. In this approach we assume that each source has its individual relaxation timescale and residual quiet‐time strength, and its partial contribution to the total field depends on the entire history of the external driving of the magnetosphere during a storm. In addition, the magnitudes of the principal field sources were assumed to saturate during extremely large storms with abnormally strong external driving. All the parameters of the model field sources, including their magnitudes, geometrical characteristics, solar wind/IMF driving functions, decay timescales, and saturation thresholds, were treated as free variables, and their values were derived from the data. As an independent consistency test, we calculated the expected Dst variation on the basis of the model output at Earth's surface and compared it with the actual observed Dst. A good agreement (cumulative correlation coefficient R = 0.92) was found, in spite of the fact that ∼90% of the spacecraft data used in the fitting were taken at synchronous orbit and beyond, while only 3.7% of those data came from distances 2.5 ≤ R ≤ 4 RE. The obtained results demonstrate the possibility to develop a truly dynamical model of the magnetic field, based on magnetospheric and interplanetary data and allowing one to reproduce and forecast the entire process of a geomagnetic storm, as it unfolds in time and space.
Dipolarization fronts (DFs), characterized by a strong and steep increase of the tail magnetic field component Bz normal to the neutral plane and preceded by a much less negative dip of Bz, are ...reported in many observations of bursty bulk flows and substorm activations throughout the whole Earth's magnetotail. It is shown that similar structures appear in full‐particle simulations with open boundaries in a transient regime before the steady reconnection in the original Harris current sheet driven out of the equilibrium by the initial X‐line perturbation is established. Being secondary reconnection structures propagating with the Alfvén speed, DFs are different from the magnetic field pileup regions reported in earlier simulations with closed boundaries. They also differ from the secondary plasmoids with bipolar Bz changes reported in earlier fluid simulations and particle simulations with open boundaries. In spite of their transient nature, DFs are found to form when the force balance is already restored in the system, which justifies their interpretation as a nonlinear stage of the tearing instability developing in two magnetotail‐like structures on the left and on the right of the initial central X‐line. Both electrons and ions are magnetized at the front of the dipolarization wave. In contrast, in its trail, ions are unmagnetized and move slower compared to the E × B drift, whereas electrons either follow that drift being completely magnetized or move faster, forming super‐Alfvénic jets. In spite of the different motions of electrons and ions, the growth of the front is not accompanied by the corresponding growth of the electrostatic field and the energy dissipation in fronts is dominated by ions.
Magnetospheric substorms represent key explosive processes in the interaction of the Earth's magnetosphere with the solar wind, and their understanding and modeling are critical for space weather ...forecasting. During substorms, the magnetic field on the nightside is first stretched in the antisunward direction and then it rapidly contracts earthward bringing hot plasmas from the distant space regions into the inner magnetosphere, where they contribute to geomagnetic storms and Joule dissipation in the polar ionosphere, causing impressive splashes of aurora. Here we show for the first time that mining millions of spaceborne magnetometer data records from multiple missions allows one to reconstruct the global 3‐D picture of these stretching and dipolarization processes. Stretching results in the formation of a thin (less than the Earth's radius) and strong current sheet, which is diverted into the ionosphere during dipolarization. In the meantime, the dipolarization signal propagates further into the inner magnetosphere resulting in the accumulation of a longer lived current there, giving rise to a protogeomagnetic storm. The global 3‐D structure of the corresponding substorm currents including the substorm current wedge is reconstructed from data.
Plain Language Summary
Using several millions of historical magnetometer records and data mining techniques, we form virtual spacecraft constellations of tens of thousands of spacecraft to reconstruct the global shape of the terrestrial magnetosphere at the moments of its most dramatic reconfigurations responsible for major space weather disturbances.
Key Points
Substorm tail current sheet thinning and dipolarization are reproduced using novel data mining technique
Global 3‐D structure of substorm currents including the substorm current wedge is reconstructed from data
Substorms contribute to an accumulation of a longer‐lived thick current in the innermost part of the magnetosphere
The onset of reconnection in 2‐D current sheet equilibria that include an X line separating tail‐like regions with magnetized electrons is simulated with a full‐particle code. The onset is driven by ...a finite convection electric field applied outside the current sheet. In the case of tearing stable tails with no accumulated magnetic flux, the convection electric field penetrates the sheet near the X line. In contrast, in multiscale equilibria where the X line is framed by local areas of enhanced flux, the electric field avoids the X line, directly penetrates the areas of increased flux, and ejects them downstream. The ejecta form dipolarization fronts (DFs), sharp magnetic pileups with a thickness on the order of the ion inertial length, much smaller than the mesoscales of the initial flux increase regions. The DFs move with the reconnection outflows in the direction opposite the magnetic field stretching, while behind them new X lines, distinct from the original, form. Simulations with a reduced driving field suggest that DF formation shares properties with the ion tearing instability, which is consistent with its potential destabilization in multiscale equilibria. Weak driving of equilibria with tearing stable tails first forms flux accumulation regions, which then rapidly transform into DFs, making 2‐D equilibria inherently metastable. The results are compared with observations of DFs, the statistical visualization of Earth's magnetotail during substorm onset, and the bubble‐blob pair formation model.
Key Points
Formation of dipolarization fronts is a key feature of reconnection onset
Reconnection onset in 2‐D equilbria is consistent with tearing stability theory
The 2‐D reconnecting current sheet behaves as a metastable system
Explosive magnetotail activity has long been understood in the context of its auroral manifestations. While global models have been used to interpret and understand many magnetospheric processes, the ...temporal and spatial scales of some auroral forms have been inaccessible to global modeling creating a gulf between observational and theoretical studies of these phenomena. We present here an important step toward bridging this gulf using a newly developed global magnetosphere‐ionosphere model with resolution capturing
≲ 30 km azimuthal scales in the auroral zone. In a global magnetohydrodynamic (MHD) simulation of the growth phase of a synthetic substorm, we find the self‐consistent formation and destabilization of localized magnetic field minima in the near‐Earth magnetotail. We demonstrate that this destabilization is due to ballooning‐interchange instability which drives earthward entropy bubbles with embedded magnetic fronts. Finally, we show that these bubbles create localized field‐aligned current structures that manifest in the ionosphere with properties matching observed auroral beads.
Plain Language Summary
The aurora has long been used as a window onto the magnetosphere. However, auroral observations are inherently limited in trying to reconstruct global magnetospheric dynamics from the “magnetic shadow” they cast on Earth. For this reason modeling has been used in tandem with observations to better contextualize and understand the data. Substorms, the violent reconfiguration of the magnetotail and one of the most dynamic magnetospheric phenomena, have been known to be preceded by the formation of bead‐like structures in the aurora. The processes responsible for auroral beading and their causal versus correlative role with substorm onset have remained an enduring mystery. The vast disparity between the spatial scales of auroral beads and those of the global magnetosphere has greatly complicated the use of modeling in unraveling this mystery. We show here for the first time a demonstration of the self‐consistent formation of a magnetospheric configuration that becomes unstable during the period preceding the substorm onset and that this instability manifests in the ionosphere with similar morphology to auroral beads. The global context of the model shows that the magnetospheric processes responsible for beading are not necessarily causal to onset but a consequence of the slow magnetotail reconfiguration that precedes onset.
Key Points
We present the first global magnetosphere simulation to reveal ballooning‐interchange instability of a narrow Bz minimum in the near‐Earth magnetotail
The instability is prominent during the substorm growth phase and generates earthward entropy bubbles with embedded magnetic fronts
The bubbles drive mesoscale ionospheric field‐aligned currents and auroral structures (beads) with properties matching to those observed
The sufficient stability criterion of the collisionless ion tearing mode in the magnetotail current sheet, which was first obtained by Lembege and Pellat in 1982, is considered. For many conventional ...2D current sheet equilibria, this criterion is satisfied within the WKB approximation, which is commonly interpreted as stability of those equilibria with respect to tearing. However, this is not necessarily the case for equilibria with more than two characteristic spatial scales. An example for substantial tearing destabilization of an equilibrium with accumulation of the magnetic flux at the tailward end of a thin current sheet is presented. Similar equilibria are reported in Geotail and THEMIS observations prior to onsets of magnetospheric substorms and dipolarization fronts associated with bursty bulk flows.
Kink‐like flapping motions of current sheets are commonly observed in the magnetotail. Such oscillations have periods of a few minutes down to a few seconds and they propagate toward the flanks of ...the plasma sheet. Here, we report a short‐period (T≈25 s) flapping event of a thin current sheet observed by the Magnetospheric Multiscale spacecraft in the dusk‐side plasma sheet following a fast Earthward plasma flow. We characterize the flapping structure using the multi‐spacecraft spatiotemporal derivative and timing methods, and we find that the wave‐like structure is propagating along the average current direction with a phase velocity comparable to the ion velocity. We show that the wavelength of the oscillating current sheet scales with its thickness as expected for a drift‐kink mode. The decoupling of the ion bulk motion from the electron bulk motion suggests that the current sheet is thin. We discuss the presence of the lower hybrid waves associated with gradients of density as a broadening process of the thin current sheet.
Key Points
Kink‐like flapping ion‐scale current sheet (CS) propagating along the current direction is observed in the dusk‐side plasma sheet
The wavenumber k of the flapping oscillations scales with CS thickness h as kh ∼ 1.15, consistent with a drift‐kink instability
Lower hybrid drift waves are observed when the CS is the thinnest. Growth of the waves can be one of the factors limiting CS thinning