First observation of fluid-like eddy-dominant bursty bulk flow turbulence in the Earth's tail plasma sheet

Turbulence is a ubiquitous phenomenon in neutral and conductive fluids. According to classical theory, turbulence is a rotating flow containing vortices of different scales. Eddies play a fundamental role in the nonlinear cascade of kinetic energy at different scales in turbulent flow. In conductive fluids, the Alfvénic/kinetic Alfvénic wave (AW/KAW) is the new "cell" of magnetohydrodynamic (MHD) turbulence (frozen-in condition). Wave energy, which has equal kinetic and magnetic energy, is redistributed among multiple-scale Fourier modes and transferred from the large MHD scale to the small kinetic scale through the collision of counter-propagating Alfvénic wave packages propagating along the magnetic field line. Fluid-like eddy-dominant plasma flow turbulence has never been found in space since the launch of the first satellite in 1957. In this paper, we report the first observation of eddy-dominant turbulence within magnetic reconnection-generated fast flow in the Earth's tail plasma sheet by the Magnetospheric Multiscale Spacecraft (MMS). In eddy-dominant turbulent reconnection jet, ions dominate the flow field while electrons dominate current and magnetic fluctuations. Our findings shed new light on the nonlinear kinetic and magnetic energy cascade in MHD turbulence.

Turbulent current sheet frozen in bursty bulk flow: observation and model

Utilizing four-point joint observations by Magnetospheric Multiscale Spacecraft (MMS), we investigate the main features of the current sheet frozen in (CSFI) the bursty bulk flow. Typical event on the steady long-lasting BBF on July 23, 2017 shows the enhanced dawn-dusk current (Jy₀) in the CSFI (β ~ 10). The magnitude of the Jy₀ in the CSFI is about 5.5 nA/m². The CSFI is highly turbulent, with the ratio of ∆J/J₀ of ~ 2 (where ∆J is perturbed J). The turbulent CSFI is characterized by intermittent current coherent structures. The magnitude of the spiky-J at coherent structures is typically above 30 nA/m². Spectrum analysis exhibits that BBF turbulence follows distinct dissipation laws inside and outside the CSFI. Based on MMS observations, we propose a new model of the BBF in the framework of magnetohydrodynamics. In this model, the BBF is depicted as a closed plasma system with the localized current sheet frozen at the center of the flow (Taylor’s hypothesis). In the light of principle of Helmholtz-decomposition, the BBF motion in the tail plasma sheet is explained. The model also predicts the thermal expansion of the BBF after leaving the reconnection source region.

Magnetotail Ion Structuring by Kinetic Ballooning-Interchange Instability

By combining three-probe THEMIS observations and 3-D Particle-in-Cell simulations, we identify key structures on the ion gyroradius scale that occur in connection with ballooning-interchange instability heads in the Earth's magnetotail. The mesoscale structures occur at sites of strong ion velocity shear and vorticity where the thermal ion Larmor radius is about half of the width of the head. Finer structures occur at the smaller scales characterizing the wavelength of the electromagnetic ion cyclotron waves generated at the heads. These two processes act to erode and thin the current sheet, thereby forming a local magnetotail configuration that is favorable for reconnection.

Characteristics of Plasma Dynamics in Current Sheets Formed in Helium Plasma

The characteristic features of plasma acceleration in the current sheets are discussed on the basis of an analysis of the structure of electrodynamic forces at successive stages of the evolution of the current sheets formed in the plasma with helium ions. Of particular interest is the generation of reverse currents at the side edges of the sheet and the appearance of forces, which are braking previously accelerated plasma flows.

Multifractal Characteristics of Geomagnetic Field Fluctuations for the Northern and Southern Hemispheres at Swarm Altitude

This paper explores the spatial variations of the statistical scaling features of low to high latitude geomagnetic field fluctuations at Swarm altitude. The data for this study comes from the vector field magnetometer onboard Swarm A satellite, measured at low resolution (1 Hz) for one year (from 9 March 2016, to 9 March 2017). We estimated the structure-function scaling exponents using the p-leaders discrete wavelet multifractal technique, from which we obtained the singularity spectrum related to the magnetic fluctuations in the North-East-Center (NEC) coordinate system. From this estimation, we retain just the maximal fractal subset, associated with the Hurst exponent H. Here we present thresholding for two levels of the Auroral Electrojet index and almost the whole northern and southern hemispheres, the Hurst exponent, the structure-function scaling exponent of order 2, and the multifractal p-exponent width for the geomagnetic fluctuations. The latter quantifies the relevance of the multifractal property. Sometimes, we found negative values of H, suggesting a behavior similar to wave breaking or shocklet-like propagating front. Furthermore, we found some asymmetries in the magnetic field turbulence between the northern and southern hemispheres. These estimations suggest that different turbulent regimes of the geomagnetic field fluctuations exist along the Swarm path.

Long-Term Variations of Plasmaspheric Total Electron Content from Topside GPS Observations on LEO Satellites

The plasmasphere is located above the ionosphere with low-energy plasma, which is an important component of the solar-terrestrial space environment. As the link between the ionosphere and the magnetosphere, the plasmasphere plays an important role in the coupling process. Therefore, it is of great significance to study the electron content variation of the plasmasphere for the solar-terrestrial space environment. Nowadays, the topside global positioning system (GPS) observations on Low Earth Orbit (LEO) satellites provide a unique opportunity to estimate and study variations in the plasmasphere. In this paper, the plasmaspheric total electron content (PTEC) is estimated, and its long-term variations are studied from topside GPS observations onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC). The PTEC in the daytime is higher than that in the nighttime, with the peak between 14:00 and 17:00 in the magnetic local time, while the minimum value of PTEC in the belt appears between 3:00 and 6:00 in the magnetic local time before sunrise. For seasonal variations, the PTEC is the highest in spring of the northern hemisphere and the lowest in summer of the northern hemisphere regardless of the state of the solar activity. The long-term variation in PTEC is further analyzed using 11-year COSMIC GPS observation data from 2007 to 2017. A high correlation between PTEC and the F10.7 indices is found. Particularly in the geomagnetic high-latitude region during the daytime, the correlation coefficient reaches 0.93. The worst case occurs during the nighttime in the geomagnetic middle-latitude region, but the correlation coefficient is still higher than 0.88. The long-term variations of plasmaspheric TEC are mainly related to the solar activity.

Contribution of Bursty Bulk Flows to the Global Dipolarization of the Magnetotail During an Isolated Substorm

This paper addresses the question of the contribution of azimuthally localized flow channels and magnetic field dipolarizations embedded in them in the global dipolarization of the inner magnetosphere during substorms. We employ the high-resolution Lyon-Fedder-Mobarry global magnetosphere magnetohydrodynamic model and simulate an isolated substorm event, which was observed by the geostationary satellites and by the Magnetospheric Multiscale spacecraft. The results of our simulations reveal that plasma sheet flow channels (bursty bulk flows, BBFs) and elementary dipolarizations (dipolarization fronts, DFs) occur in the growth phase of the substorm but are rare and do not penetrate to the geosynchronous orbit. The substorm onset is characterized by an abrupt increase in the occurrence and intensity of BBFs/DFs, which penetrate well earthward of the geosynchronous orbit during the expansion phase. These azimuthally localized structures are solely responsible for the global (in terms of the magnetic local time) dipolarization of the inner magnetosphere toward the end of the substorm expansion. Comparison with the geostationary satellites and Magnetospheric Multiscale data shows that the properties of the BBFs/DFs in the simulation are similar to those observed, which gives credence to the above results. Additionally, the simulation reveals many previously observed signatures of BBFs and DFs, including overshoots and oscillations around their equilibrium position, strong rebounds and vortical tailward flows, and the corresponding plasma sheet expansion and thinning.

Explosive Magnetotail Activity

Modes and manifestations of the explosive activity in the Earth's magnetotail, as well as its onset mechanisms and key pre-onset conditions are reviewed. Two mechanisms for the generation of the pre-onset current sheet are discussed, namely magnetic flux addition to the tail lobes, or other high-latitude perturbations, and magnetic flux evacuation from the near-Earth tail associated with dayside reconnection. Reconnection onset may require stretching and thinning of the sheet down to electron scales. It may also start in thicker sheets in regions with a tailward gradient of the equatorial magnetic field B z ; in this case it begins as an ideal-MHD instability followed by the generation of bursty bulk flows and dipolarization fronts. Indeed, remote sensing and global MHD modeling show the formation of tail regions with increased B z , prone to magnetic reconnection, ballooning/interchange and flapping instabilities. While interchange instability may also develop in such thicker sheets, it may grow more slowly compared to tearing and cause secondary reconnection locally in the dawn-dusk direction. Post-onset transients include bursty flows and dipolarization fronts, micro-instabilities of lower-hybrid-drift and whistler waves, as well as damped global flux tube oscillations in the near-Earth region. They convert the stretched tail magnetic field energy into bulk plasma acceleration and collisionless heating, excitation of a broad spectrum of plasma waves, and collisional dissipation in the ionosphere. Collisionless heating involves ion reflection from fronts, Fermi, betatron as well as other, non-adiabatic, mechanisms. Ionospheric manifestations of some of these magnetotail phenomena are discussed. Explosive plasma phenomena observed in the laboratory, the solar corona and solar wind are also discussed.

Energetic Ion Injections Inside Geosynchronous Orbit: Convection- and Drift-Dominated, Charge-Dependent Adiabatic Energization (W = qEd)

Particle injection, a major mode of plasma transport and energization throughout the magnetosphere, has been studied for decades. Nonetheless, the physical processes that lead to the acceleration and transport of very energetic ions in the inner magnetosphere during injection events are still under debate. In this paper, we analyze several injection events occurring near the Van Allen Probes apogee. Our analysis shows that the highest energy of an injected ion population depends on the charge state of that population. We show that most of the helium injected is doubly ionized (He++), while oxygen charge states are consistent with the presence of both ionospheric (O+) and solar wind (O6+) source populations. Based on the findings of our data analysis and with the use of a simple model, we demonstrate that the behavior of each injection of energetic ions near the Van Allen Probes apogee (5 < L < 7 R E) is well explained by simple adiabatic or nearly adiabatic transport within flow channels from higher L (≥10 R E) with velocities at 10 R E ranging between ~200 and 2,000 km/s and falling with inward transport consistent with fixed potential drops across the flow channels. Gradient/curvature drift during transport limits the highest energy/charge observed for each injection at the Van Allen Probes. Even at the highest measured ion energies where gyroradius and scattering effects might be expected to appear, energization depends on charge state but not on ion mass.

Magnetotail Fast Flow Occurrence Rate and Dawn-Dusk Asymmetry at XGSM ∼ -60 RE

As a direct result of magnetic reconnection, plasma sheet fast flows act as primary transporter of mass, flux, and energy in the Earth's magnetotail. During the last decades, these flows were mainly studied within X_GSM>-60R_E , as observations near or beyond lunar orbit were limited. By using 5 years (2011-2015) of ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moons Interaction with the Sun) data, we statistically investigate earthward and tailward flows at around 60 R_E downtail. A significant fraction of fast flows is directed earthward, comprising 43% (v_x >400 km/s) to 56% (v_x >100 km/s) of all observed flows. This suggests that near-Earth and midtail reconnection are equally probable of occurring on either side of the ARTEMIS downtail distance. For fast convective flows (v_⊥x >400 km/s), this fraction of earthward flows is reduced to about 29%, which is in line with reconnection as source of these flows and a downtail decreasing Alfvén velocity. More than 60% of tailward convective flows occur in the dusk sector (as opposed to the dawn sector), while earthward convective flows are nearly symmetrically distributed between the two sectors for low AL (>-400 nT) and asymmetrically distributed toward the dusk sector for high AL (<-400 nT). This indicates that the dawn-dusk asymmetry is more pronounced closer to Earth and moves farther downtail during high geomagnetic activity. This is consistent with similar observations pointing to the asymmetric nature of tail reconnection as the origin of the dawn-dusk asymmetry of flows and other related observables. We infer that near-Earth reconnection is preferentially located at dusk, whereas midtail reconnection (X >- 60R_E ) is likely symmetric across the tail during weak substorms and asymmetric toward the dusk sector for strong substorms, as the dawn-dusk asymmetric nature of reconnection onset in the near-Earth region progresses downtail.

Magnetotail energy dissipation during an auroral substorm

Violent releases of space plasma energy from the Earth's magnetotail during substorms produce strong electric currents and bright aurora. But what modulates these currents and aurora and controls dissipation of the energy released in the ionosphere? Using data from the THEMIS fleet of satellites and ground-based imagers and magnetometers, we show that plasma energy dissipation is controlled by field-aligned currents (FACs) produced and modulated during magnetotail topology change and oscillatory braking of fast plasma jets at 10-14 Earth radii in the nightside magnetosphere. FACs appear in regions where plasma sheet pressure and flux tube volume gradients are non-collinear. Faster tailward expansion of magnetotail dipolarization and subsequent slower inner plasma sheet restretching during substorm expansion and recovery phases cause faster poleward then slower equatorward movement of the substorm aurora. Anharmonic radial plasma oscillations build up displaced current filaments and are responsible for discrete longitudinal auroral arcs that move equatorward at a velocity of about 1km/s. This observed auroral activity appears sufficient to dissipate the released energy.

Problems with mapping the auroral oval and magnetospheric substorms

Accurate mapping of the auroral oval into the equatorial plane is critical for the analysis of aurora and substorm dynamics. Comparison of ion pressure values measured at low altitudes by Defense Meteorological Satellite Program (DMSP) satellites during their crossings of the auroral oval, with plasma pressure values obtained at the equatorial plane from Time History of Events and Macroscale Interactions during Substorms (THEMIS) satellite measurements, indicates that the main part of the auroral oval maps into the equatorial plane at distances between 6 and 12 Earth radii. On the nightside, this region is generally considered to be a part of the plasma sheet. However, our studies suggest that this region could form part of the plasma ring surrounding the Earth. We discuss the possibility of using the results found here to explain the ring-like shape of the auroral oval, the location of the injection boundary inside the magnetosphere near the geostationary orbit, presence of quiet auroral arcs in the auroral oval despite the constantly high level of turbulence observed in the plasma sheet, and some features of the onset of substorm expansion.

Rapid acceleration of protons upstream of earthward propagating dipolarization fronts

[1] Transport and acceleration of ions in the magnetotail largely occurs in the form of discrete impulsive events associated with a steep increase of the tail magnetic field normal to the neutral plane (Bz ), which are referred to as dipolarization fronts. The goal of this paper is to investigate how protons initially located upstream of earthward moving fronts are accelerated at their encounter. According to our analytical analysis and simplified two-dimensional test-particle simulations of equatorially mirroring particles, there are two regimes of proton acceleration: trapping and quasi-trapping, which are realized depending on whether the front is preceded by a negative depletion in Bz . We then use three-dimensional test-particle simulations to investigate how these acceleration processes operate in a realistic magnetotail geometry. For this purpose we construct an analytical model of the front which is superimposed onto the ambient field of the magnetotail. According to our numerical simulations, both trapping and quasi-trapping can produce rapid acceleration of protons by more than an order of magnitude. In the case of trapping, the acceleration levels depend on the amount of time particles stay in phase with the front which is controlled by the magnetic field curvature ahead of the front and the front width. Quasi-trapping does not cause particle scattering out of the equatorial plane. Energization levels in this case are limited by the number of encounters particles have with the front before they get magnetized behind it.

Energy partition in magnetic reconnection in Earth's magnetotail

The partition of energy flux in magnetic reconnection is examined experimentally using Cluster satellite observations of collisionless reconnection in Earth's magnetotail. In this plasma regime, the dominant component of the energy flux is ion enthalpy flux, with smaller contributions from the electron enthalpy and heat flux and the ion kinetic energy flux. However, the Poynting flux is not negligible, and in certain parts of the ion diffusion region the Poynting flux in fact dominates. Evidence for earthward-tailward asymmetry is ascribed to the presence of Earth's dipole fields.

New features of electron phase space holes observed by the THEMIS mission

Observations of electron phase-space holes (EHs) in Earth's plasma sheet by the THEMIS satellites include the first detection of a magnetic perturbation (deltaB_{ parallel}) parallel to the ambient magnetic field (B0). EHs with a detectable deltaB_{ parallel} have several distinguishing features including large electric field amplitudes, a magnetic perturbation perpendicular to B0, high speeds ( approximately 0.3c) along B0, and sizes along B0 of tens of Debye lengths. These EHs have a significant center potential (Phi approximately k_{B}T_{e}/e), suggesting strongly nonlinear behavior nearby such as double layers or magnetic reconnection.

Observations of double layers in earth's plasma sheet

We report the first direct observations of parallel electric fields (E_{ parallel}) carried by double layers (DLs) in the plasma sheet of Earth's magnetosphere. The DL observations, made by the THEMIS spacecraft, have E_{ parallel} signals that are analogous to those reported in the auroral region. DLs are observed during bursty bulk flow events, in the current sheet, and in plasma sheet boundary layer, all during periods of strong magnetic fluctuations. These observations imply that DLs are a universal process and that strongly nonlinear and kinetic behavior is intrinsic to Earth's plasma sheet.