Three-Dimensional Energy Transfer in Space Plasma Turbulence from Multipoint Measurement.

A novel multispacecraft technique applied to Magnetospheric Multiscale Mission data in the Earth's magnetosheath enables evaluation of the energy cascade rate from the full Yaglom's equation. The method differs from existing approaches in that it (i) is inherently three-dimensional, (ii) provides a statistically significant number of estimates from a single data stream, and (iii) allows visualization of energy flux in turbulent plasmas. This new "lag polyhedral derivative ensemble" technique exploits ensembles of tetrahedra in lag space and established curlometerlike algorithms.

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 (Jy0) in the CSFI (ß ~ 10). The magnitude of the Jy0 in the CSFI is about 5.5 nA/m2. The CSFI is highly turbulent, with the ratio of ∆J/J0 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/m2. 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.

Properties and Acceleration Mechanisms of Electrons Up To 200 keV Associated With a Flux Rope Pair and Reconnection X-Lines Around It in Earth's Plasma Sheet.

The properties and acceleration mechanisms of electrons (<200 keV) associated with a pair of tailward traveling flux ropes and accompanied reconnection X-lines in Earth's plasma sheet are investigated with MMS measurements. Energetic electrons are enhanced on both boundaries and core of the flux ropes. The power-law spectra of energetic electrons near the X-lines and in flux ropes are harder than those on flux rope boundaries. Theoretical calculations show that the highest energy of adiabatic electrons is a few keV around the X-lines, tens of keV immediately downstream of the X-lines, hundreds of keV on the flux rope boundaries, and a few MeV in the flux rope cores. The X-lines cause strong energy dissipation, which may generate the energetic electron beams around them. The enhanced electron parallel temperature can be caused by the curvature-driven Fermi acceleration and the parallel electric potential. Betatron acceleration due to the magnetic field compression is strong on flux rope boundaries, which enhances energetic electrons in the perpendicular direction. Electrons can be trapped between the flux rope pair due to mirror force and parallel electric potential. Electrostatic structures in the flux rope cores correspond to potential drops up to half of the electron temperature. The energetic electrons and the electron distribution functions in the flux rope cores are suggested to be transported from other dawn-dusk directions, which is a 3-dimensional effect. The acceleration and deceleration of the Betatron and Fermi processes appear alternately indicating that the magnetic field and plasma are turbulent around the flux ropes.

Collisionless relaxation of a disequilibrated current sheet and implications for bifurcated structures.

Current sheets are ubiquitous plasma structures that play the crucial role of being energy sources for various magnetic phenomena. Although a plethora of current sheet equilibrium solutions have been found, the collisionless process through which a disequilibrated current sheet relaxes or equilibrates remains largely unknown. Here we show, through analyses of phase-space distributions of single-particle orbit classes and particle-in-cell simulations, that collisionless transitions among the orbit classes are responsible for this process. Bifurcated current sheets, which are readily observed in geospace but whose origins remain controversial, are shown to naturally arise from the equilibration process and thus are likely to be the underlying structures in various phenomena; comparisons of spacecraft observations to particle-in-cell simulations support this fact. The bearing of this result on previous explanations of bifurcated structures is also discussed.

Self-consistent kinetic model of nested electron- and ion-scale magnetic cavities in space plasmas.

NASA's Magnetospheric Multi-Scale (MMS) mission is designed to explore the proton- and electron-gyroscale kinetics of plasma turbulence where the bulk of particle acceleration and heating takes place. Understanding the nature of cross-scale structures ubiquitous as magnetic cavities is important to assess the energy partition, cascade and conversion in the plasma universe. Here, we present theoretical insight into magnetic cavities by deriving a self-consistent, kinetic theory of these coherent structures. By taking advantage of the multipoint measurements from the MMS constellation, we demonstrate that our kinetic model can utilize magnetic cavity observations by one MMS spacecraft to predict measurements from a second/third spacecraft. The methodology of "observe and predict" validates the theory we have derived, and confirms that nested magnetic cavities are self-organized plasma structures supported by trapped proton and electron populations in analogous to the classical theta-pinches in laboratory plasmas.

The Effects of Upper-Hybrid Waves on Energy Dissipation in the Electron Diffusion Region.

Using a two-dimensional particle-in-cell simulation, we investigate the effects and roles of upper-hybrid waves (UHW) near the electron diffusion region (EDR). The energy dissipation via the wave-particle interaction in our simulation agrees with J · E ' measured by magnetospheric multiscale (MMS) spacecraft. It means that UHW contributes to the local energy dissipation. As a result of wave-particle interactions, plasma parameters which determine the larger-scale energy dissipation in the EDR are changed. The y-directional current decreases while the pressure tensor P y z increases/decreases when the agyrotropic beam density is low/high, where (x, y, z)-coordinates correspond the (L, M, N)-boundary coordinates. Because the reconnection electric field comes from -∂ P y z /∂ z, our result implies that UHW plays an additional role in affecting larger-scale energy dissipation in the EDR by changing plasma parameters. We provide a simple diagram that shows how the UHW activities change the profiles of plasma parameters near the EDR comparing cases with and without UHW.

High-Frequency Waves Driven by Agyrotropic Electrons Near the Electron Diffusion Region.

National Aeronautics and Space Administration's Magnetosphere Multiscale mission reveals that agyrotropic electrons and intense waves are prevalently present in the electron diffusion region. Prompted by two distinct Magnetosphere Multiscale observations, this letter investigates by theoretical means and the properties of agyrotropic electron beam-plasma instability and explains the origin of different structures in the wave spectra. The difference is owing to the fact that in one instance, a continuous beam mode is excited, while in the other, discrete Bernstein modes are excited, and the excitation of one mode versus the other depends on physical input parameters, which are consistent with observations. Analyses of dispersion relations show that the growing mode becomes discrete when the maximum growth rate is lower than the electron cyclotron frequency. Making use of particle-in-cell simulations, we found that the broadening angle Δ in the gyroangle space is also an important factor controlling the growth rate. Ramifications of the present finding are also discussed.

Statistics of Kinetic Dissipation in the Earth's Magnetosheath: MMS Observations.

A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy conversion becomes more complex than in the fluid case, and both pressure-dilatation and pressure strain interactions (Pi-D≡-Π_{ij}D_{ij}) become relevant and potentially important. Pi-D in plasma turbulence has been studied so far primarily using simulations. The present study provides a statistical analysis of Pi-D in the Earth's magnetosheath using the unique measurement capabilities of the Magnetospheric Multiscale (MMS) mission. We find that the statistics of Pi-D in this naturally occurring plasma environment exhibit strong resemblance to previously established fully kinetic simulations results. The conversion of energy is concentrated in space and occurs near intense current sheets, but not within them. This supports recent suggestions that the chain of energy transfer channels involves regional, rather than pointwise, correlations.

In Situ Observation of Hall Magnetohydrodynamic Cascade in Space Plasma.

We present estimates of the turbulent energy-cascade rate derived from a Hall-magnetohydrodynamic (MHD) third-order law. We compute the contribution from the Hall term and the MHD term to the energy flux. Magnetospheric Multiscale (MMS) data accumulated in the magnetosheath and the solar wind are compared with previously established simulation results. Consistent with the simulations, we find that at large (MHD) scales, the MMS observations exhibit a clear inertial range dominated by the MHD flux. In the subion range, the cascade continues at a diminished level via the Hall term, and the change becomes more pronounced as the plasma beta increases. Additionally, the MHD contribution to interscale energy transfer remains important at smaller scales than previously thought. Possible reasons are offered for this unanticipated result.

MMS SITL Ground Loop: Automating the Burst Data Selection Process.

Global-scale energy flow throughout Earth's magnetosphere is catalyzed by processes that occur at Earth's magnetopause (MP). Magnetic reconnection is one process responsible for solar wind entry into and global convection within the magnetosphere, and the MP location, orientation, and motion have an impact on the dynamics. Statistical studies that focus on these and other MP phenomena and characteristics inherently require MP identification in their event search criteria, a task that can be automated using machine learning so that more man hours can be spent on research and analysis. We introduce a Long-Short Term Memory (LSTM) Recurrent Neural Network model to detect MP crossings and assist studies of energy transfer into the magnetosphere. As its first application, the LSTM has been implemented into the operational data stream of the Magnetospheric Multiscale (MMS) mission. MMS focuses on the electron diffusion region of reconnection, where electron dynamics break magnetic field lines and plasma is energized. MMS employs automated burst triggers onboard the spacecraft and a Scientist-in-the-Loop (SITL) on the ground to select intervals likely to contain diffusion regions. Only low-resolution survey data is available to the SITL, which is insufficient to resolve electron dynamics. A strategy for the SITL, then, is to select all MP crossings. Of all 219 SITL selections classified as MP crossings during the first five months of model operations, the model predicted 166 (76%) of them, and of all 360 model predictions, 257 (71%) were selected by the SITL. Most predictions that were not classified as MP crossings by the SITL were still MP-like, in that the intervals contained mixed magnetosheath and magnetospheric plasmas. The LSTM model and its predictions are public to ease the burden of arduous event searches involving the MP, including those for EDRs. For MMS, this helps free up mission operation costs by consolidating manual classification processes into automated routines.

High-Frequency Wave Generation in Magnetotail Reconnection: Nonlinear Harmonics of Upper Hybrid Waves.

MMS3 spacecraft passed the vicinity of the electron diffusion region of magnetotail reconnection on 3 July 2017, observing discrepancies between perpendicular electron bulk velocities and E → × B → drift, and agyrotropic electron crescent distributions. Analyzing linear wave dispersions, Burch et al. (2019, https://doi.org/10.1029/2019GL082471) showed the electron crescent generates high-frequency waves. We investigate harmonics of upper-hybrid (UH) waves using both observation and particle-in-cell (PIC) simulation, and the generation of electromagnetic radiation from PIC simulation. Harmonics of UH are linearly polarized and propagate along the perpendicular direction to the ambient magnetic field. Compared with two-dimensional PIC simulation and nonlinear kinetic theory, we show that the nonlinear beam-plasma interaction between the agyrotropic electrons and the core electrons generates harmonics of UH. Moreover, PIC simulation shows that agyrotropic electron beam can lead to electromagnetic (EM) radiation at the plasma frequency and harmonics.

Structure of the Current Sheet in the 11 July 2017 Electron Diffusion Region Event.

The structure of the current sheet along the Magnetospheric Multiscale (MMS) orbit is examined during the 11 July 2017 Electron Diffusion Region (EDR) event. The location of MMS relative to the X-line is deduced and used to obtain the spatial changes in the electron parameters. The electron velocity gradient values are used to estimate the reconnection electric field sustained by nongyrotropic pressure. It is shown that the observations are consistent with theoretical expectations for an inner EDR in 2-D reconnection. That is, the magnetic field gradient scale, where the electric field due to electron nongyrotropic pressure dominates, is comparable to the gyroscale of the thermal electrons at the edge of the inner EDR. Our approximation of the MMS observations using a steady state, quasi-2-D, tailward retreating X-line was valid only for about 1.4 s. This suggests that the inner EDR is localized; that is, electron outflow jet braking takes place within an ion inertia scale from the X-line. The existence of multiple events or current sheet processes outside the EDR may play an important role in the geometry of reconnection in the near-Earth magnetotail.

Energy partitioning constraints at kinetic scales in low-ß turbulence.

Turbulence is a fundamental physical process through which energy injected into a system at large scales cascades to smaller scales. In collisionless plasmas, turbulence provides a critical mechanism for dissipating electromagnetic energy. Here we present observations of plasma fluctuations in low-ß turbulence using data from NASA's Magnetospheric Multiscale mission in Earth's magnetosheath. We provide constraints on the partitioning of turbulent energy density in the fluid, ion-kinetic, and electron-kinetic ranges. Magnetic field fluctuations dominated the energy density spectrum throughout the fluid and ion-kinetic ranges, consistent with previous observations of turbulence in similar plasma regimes. However, at scales shorter than the electron inertial length, fluctuation power in electron kinetic energy significantly exceeded that of the magnetic field, resulting in an electron-motion-regulated cascade at small scales. This dominance should be highly relevant for the study of turbulence in highly magnetized laboratory and astrophysical plasmas.

Multiscale Currents Observed by MMS in the Flow Braking Region.

We present characteristics of current layers in the off-equatorial near-Earth plasma sheet boundary observed with high time-resolution measurements from the Magnetospheric Multiscale mission during an intense substorm associated with multiple dipolarizations. The four Magnetospheric Multiscale spacecraft, separated by distances of about 50 km, were located in the southern hemisphere in the dusk portion of a substorm current wedge. They observed fast flow disturbances (up to about 500 km/s), most intense in the dawn-dusk direction. Field-aligned currents were observed initially within the expanding plasma sheet, where the flow and field disturbances showed the distinct pattern expected in the braking region of localized flows. Subsequently, intense thin field-aligned current layers were detected at the inner boundary of equatorward moving flux tubes together with Earthward streaming hot ions. Intense Hall current layers were found adjacent to the field-aligned currents. In particular, we found a Hall current structure in the vicinity of the Earthward streaming ion jet that consisted of mixed ion components, that is, hot unmagnetized ions, cold E × B drifting ions, and magnetized electrons. Our observations show that both the near-Earth plasma jet diversion and the thin Hall current layers formed around the reconnection jet boundary are the sites where diversion of the perpendicular currents take place that contribute to the observed field-aligned current pattern as predicted by simulations of reconnection jets. Hence, multiscale structure of flow braking is preserved in the field-aligned currents in the off-equatorial plasma sheet and is also translated to ionosphere to become a part of the substorm field-aligned current system.

Statistical Study of the Properties of Magnetosheath Lion Roars.

Lion roars are narrowband whistler wave emissions that have been observed in several environments, such as planetary magnetosheaths, the Earth's magnetosphere, the solar wind, downstream of interplanetary shocks, and the cusp region. We present measurements of more than 30,000 such emissions observed by the Magnetospheric Multiscale spacecraft with high-cadence (8,192 samples/s) search coil magnetometer data. A semiautomatic algorithm was used to identify the emissions, and an adaptive interval algorithm in conjunction with minimum variance analysis was used to determine their wave vector. The properties of the waves are determined in both the spacecraft and plasma rest frame. The mean wave normal angle, with respect to the background magnetic field (B 0), plasma bulk flow velocity (V b ), and the coplanarity plane (V b × B 0) are 23°, 56°, and 0°, respectively. The average peak frequencies were ~31% of the electron gyrofrequency (ω ce ) observed in the spacecraft frame and ~18% of ω ce in the plasma rest frame. In the spacecraft frame, ~99% of the emissions had a frequency < ω ce , while 98% had a peak frequency < 0.72 ω ce in the plasma rest frame. None of the waves had frequencies lower than the lower hybrid frequency, ω. From the probability density function of the electron plasma ß e , the ratio between the electron thermal and magnetic pressure, ~99.6% of the waves were observed with ß e < 4 with a large narrow peak at 0.07 and two smaller, but wider, peaks at 1.26 and 2.28, while the average value was ~1.25.

Wave-particle energy exchange directly observed in a kinetic Alfvén-branch wave.

Alfvén waves are fundamental plasma wave modes that permeate the universe. At small kinetic scales, they provide a critical mechanism for the transfer of energy between electromagnetic fields and charged particles. These waves are important not only in planetary magnetospheres, heliospheres and astrophysical systems but also in laboratory plasma experiments and fusion reactors. Through measurement of charged particles and electromagnetic fields with NASA's Magnetospheric Multiscale (MMS) mission, we utilize Earth's magnetosphere as a plasma physics laboratory. Here we confirm the conservative energy exchange between the electromagnetic field fluctuations and the charged particles that comprise an undamped kinetic Alfvén wave. Electrons confined between adjacent wave peaks may have contributed to saturation of damping effects via nonlinear particle trapping. The investigation of these detailed wave dynamics has been unexplored territory in experimental plasma physics and is only recently enabled by high-resolution MMS observations.

Cometary science. Birth of a comet magnetosphere: a spring of water ions.

The Rosetta mission shall accompany comet 67P/Churyumov-Gerasimenko from a heliocentric distance of >3.6 astronomical units through perihelion passage at 1.25 astronomical units, spanning low and maximum activity levels. Initially, the solar wind permeates the thin comet atmosphere formed from sublimation, until the size and plasma pressure of the ionized atmosphere define its boundaries: A magnetosphere is born. Using the Rosetta Plasma Consortium ion composition analyzer, we trace the evolution from the first detection of water ions to when the atmosphere begins repelling the solar wind (~3.3 astronomical units), and we report the spatial structure of this early interaction. The near-comet water population comprises accelerated ions (<800 electron volts), produced upstream of Rosetta, and lower energy locally produced ions; we estimate the fluxes of both ion species and energetic neutral atoms.