We performed 2D PIC simulations of Kelvin Helmholtz instability (KHI) with symmetric and asymmetric density and temperature profiles along the flow shear with a northward interplanetary magnetic field. The Magnetic Flux Transport method, field topology and magnetic field minimums are used to identify the reconnection X-lines. We start to observe the reconnection signatures such as magnetic field and flow reversals at the vortex edges in the nonlinear phase of the KHI when the vortices are rolling up. The number of reconnection regions increases at the turbulence phase. The signatures eventually decrease and finally disappear at very turbulent stages of KHI developments. Our results qualitatively agree with MMS observations of reconnection signatures at KHI along the magnetospheric flanks.

We carried out a statistical study of equatorial dipolarization fronts (DFs) detected by the Magnetospheric Multiscale mission (MMS) during the full 2017 Earth’s magnetotail season. We found that two DF classes are distinguished: class I (74.4%) corresponds to the standard DF properties and energy dissipation and a new class II (25.6%). This new class includes the six DF discussed in Alqeeq et al. (2022) and corresponds to a bump of the magnetic field associated with a minimum in the ion and electron pressures and a reversal of the energy conversion process. The possible origin of this second class is discussed. Both DF classes show that the energy conversion process in the spacecraft frame is driven by the diamagnetic current dominated by the ion pressure gradient. In the fluid frame, it is driven by the electron pressure gradient. In addition, we have shown that the energy conversion processes are not homogeneous at the electron scale mostly due to the variations of the electric fields for both DF classes.

In-situ measurements from the Magnetospheric Multiscale (MMS) mission are used to estimate electron density from spacecraft potential and investigate compressive turbulence in the Earth’s magnetosheath. During the MMS Solar Wind Turbulence Campaign in February 2019, the four MMS spacecraft were arranged in a logarithmic line constellation enabling the study of measurements from multiple spacecraft at varying distances. We estimate the electron density from spacecraft potential for a time interval in which the ion emitters actively control the potential. The derived electron density data product has a higher temporal resolution than the plasma instruments, enabling the examination of fluctuation for scales down to the sub-ion range. The inter-spacecraft separations range from 132 km to 916 km; this corresponds to scales of 3.5 to 24.1 ion inertial lengths. The derived density and magnetic field data are used to study fluctuations in the magnetosheath through time lags on a single spacecraft and spatial lags between pairs of spacecraft over almost one decade in scale. The results show an increase in anisotropy as the scale decreases.

The Kelvin-Helmholtz instability (KHI) at Earth’s magnetopause and associated turbulence are suggested to play a role in the transport of mass and momentum from the solar wind into Earth’s magnetosphere. We investigate electromagnetic turbulence observed in KH vortices encountered at the dusk flank magnetopause by the Magnetospheric Multiscale (MMS) spacecraft under northward interplanetary magnetic field (IMF) conditions in order to reveal its generation process, mode properties, and role. A comparison with another MMS event at the dayside magnetopause with reconnection but no KHI signatures under a similar IMF condition indicates that while high-latitude magnetopause reconnection excites a modest level of turbulence in the dayside low-latitude boundary layer, the KHI further amplifies the turbulence, leading to magnetic energy spectra with a power-law index –5/3 at magnetohydrodynamic scales even in its early nonlinear phase. The mode of the electromagnetic turbulence is analyzed with a single-spacecraft method based on Ampère’s law, developed by Bellan (2016), for estimating wave vectors as a function of spacecraft-frame frequency. The results suggest that the turbulence does not consist of propagating normal-mode waves, but is due to interlaced magnetic flux tubes advected by plasma flows in the vortices. The turbulence at sub-ion scales in the early nonlinear phase of the KHI may not be the cause of the plasma transport across the magnetopause, but rather a consequence of three-dimensional vortex induced reconnection, the process that can cause an efficient transport by producing tangled reconnected field lines.

Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth's bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.

Shock waves are common in the heliosphere and beyond. The collisionless nature of most astrophysical plasmas allows for the energy processed by shocks to be partitioned amongst particle sub-populations and electromagnetic fields via physical mechanisms that are not well understood. The electrostatic potential across such shocks is frame dependent. In a frame where the incident bulk velocity is parallel to the magnetic field, the deHoffmann-Teller frame, the potential is linked directly to the ambipolar electric field established by the electron pressure gradient. Thus measuring and understanding this potential solves the electron partition problem, and gives insight into other competing shock processes. Integrating measured electric fields is space is problematic since the measurements can have offsets that change with plasma conditions. The offsets, once integrated, can be as large or larger than the shock potential. Here we exploit the high-quality field and plasma measurements from NASA's Magnetospheric Multiscale mission to attempt this calculation. We investigate recent adaptations of the deHoffmann-Teller frame transformation to include time variability, and conclude that in practise these face difficulties inherent in the 3D time-dependent nature of real shocks by comparison to 1D simulations. Potential estimates based on electron fluid and kinetic analyses provide the most robust measures of the deHoffmann-Teller potential, but with some care direct integration of the electric fields can be made to agree. These results suggest that it will be difficult to independently assess the role of other processes, such as scattering by shock turbulence, in accounting for the electron heating.

Lower hybrid waves are investigated at the magnetosheath separatrix region in asymmetric guide-field reconnection by using the Magnetospheric Multiscale (MMS) mission. Three of the four MMS spacecraft observe clear wave activities around the lower hybrid frequency across the magnetosheath separatrix, where a density gradient is present. The observed waves are consistent with generation by the lower hybrid drift instability. The characteristic properties of these waves include: (1) the waves propagate toward the x-line in the spacecraft frame due to the large out-of-plane magnetic field, which is in the same direction of the diamagnetic drift of the x-line; (2) the wave potential is about 20\% of the electron temperature. These drift waves effectively produce cross-field particle diffusion, enabling the transport of magnetosheath electrons into the exhaust region. The observations presented in this study indicate unique features of asymmetric guide-field reconnection.

We present observations in Earth’s magnetotail by the Magnetospheric Multiscale spacecraft that are consistent with magnetic field annihilation, rather than magnetic topology change, causing fast magnetic-to-electron energy conversion in an electron-scale current sheet. Multi-spacecraft analysis for the magnetic field reconstruction shows that an electron-scale magnetic island was embedded in the observed electron diffusion region (EDR), suggesting an elongated shape of the EDR. Evidence for the annihilation was revealed in the form of the island growing at a rate much lower than expected for the standard collisionless reconnection, which indicates that magnetic flux injected into the EDR was not ejected from the X-point or accumulated in the island, but was dissipated in the EDR. This energy conversion process is in contrast to that in the standard EDR of a reconnecting current sheet where the energy of antiparallel magnetic fields is mostly converted to electron bulk-flow energy. Fully kinetic simulation also demonstrates that an elongated EDR is subject to the formation of electron-scale magnetic islands in which fast but transient annihilation can occur. Consistent with the observations and simulation, theoretical analysis shows that fast magnetic diffusion can occur in an elongated EDR in the presence of nongyrotropic electron effects. We suggest that the annihilation in elongated EDRs may contribute to the dissipation of magnetic energy in a turbulent collisionless plasma.

We report observations of a Bursty Bulk Flow (BBF) penetrating to the outer edge of the radiation belt. The turbulent BBF braking region is characterized by ion velocity fluctuations, magnetic field (B) variations, and intense electric fields (E). In this event, energetic (>100 keV) electron and ion fluxes are appreciably enhanced. Importantly, fluctuations in energetic electrons and ions suggest that they are locally energized. Using correlation distances and other observed characteristics of turbulent E, test-particle simulations support that local energization by E favors higher-energy electrons and leads to an enhanced energetic shoulder and tail in the electron distributions. The energetic shoulder and tail can be amplified to MeV energies by adiabatic transport into the radiation belt where |B| is higher. This analysis suggests that turbulence generated by BBFs can, in part, supply energetic particles to the outer radiation belt and that turbulence can be a significant contributor to particle acceleration.

Decomposing the electric field (E) into the contributions from generalized Ohm’s law provides key insight into both nonlinear and dissipative dynamics across the full range of scales within a plasma. Using high-resolution, multi-spacecraft measurements of three intervals in Earth’s magnetosheath from the Magnetospheric Multiscale mission, the influence of the magnetohydrodynamic, Hall, electron pressure, and electron inertia terms from Ohm’s law, as well as the impact of a finite electron mass, on the turbulent spectrum are examined observationally for the first time. The magnetohydrodynamic, Hall, and electron pressure terms are the dominant contributions to over the accessible length scales, which extend to scales smaller than the electron inertial length at the greatest extent, with the Hall and electron pressure terms dominating at sub-ion scales. The strength of the non-ideal electron pressure contribution is stronger than expected from linear kinetic Alfvén waves and a partial anti-alignment with the Hall electric field is present, linked to the relative importance of electron diamagnetic currents in the turbulence. The relative contribution of linear and nonlinear electric fields scale with the turbulent fluctuation amplitude, with nonlinear contributions playing the dominant role in shaping for the intervals examined in this study. Overall, the sum of the Ohm’s law terms and measured agree to within ~20% across the observable scales. These results both confirm general expectations about the behavior of in turbulent plasmas and highlight features that should be explored further theoretically.

Electrostatic solitary waves (ESWs) are a type of nonlinear time-domain plasma structure (TDS) generally defined by bipolar electric fields and propagation parallel to the local magnetic field. Formation mechanisms for TDSs in the magnetosphere have been studied extensively and are associated with plasma boundary layers and the braking of bursty bulk flows (BBFs). However, the rapid timescales over which these TDSs occur (< 2 ms) make them infeasible to count by eye over large time periods. Furthermore, high-cadence data are not always available. The Solitary Wave Detector (SWD) on NASA’s Magnetospheric Multiscale (MMS) mission quantifies the occurrence and amplitude of TDS throughout the constellation’s orbit; analysis of burst (65 kS/s) parallel electric field data indicates that the SWD captures appx. 60% of all bipolar TDS encountered in the tail region, enabling large-scale examination of their occurrence. Maps of TDS occurrence rates during several years of the MMS mission were generated from SWD data, showing enhanced TDS density in the tail region between 6-9 Re; enhance occurrence in or near shocks; and an unexpected enhancement in the dawn side of the tail and in the radiation belt.

Martian sub-solar electron temperatures obtained below 250 km are examined using data obtained by instruments on the Mars Atmosphere Evolution Mission (MAVEN) during the three sub-solar deep dip campaigns and a one-dimensional fluid model. This analysis was done because of the uncertainty in MAVEN low electron temperature observations at low altitudes and the fact that the Level 2 temperatures reported from the MAVEN Langmuir Probe and Waves (LPW) instrument are more than 400 Kelvin above the neutral temperatures at the lowest altitudes sampled (~120 km). These electron temperatures are well above those expected before MAVEN was launched. We find that an empirical normalization parameter, neutral pressure divided by local electron heating rate, organized the electron temperature data and identified a similar altitude (~160 km) and time scale (~2,000 s) for all three deep dips. We show that MAVEN data are not consistent with a plasma characterized by electrons in thermal equilibrium with the neutral population at 100 km. Because of the lack data below 120 km and the uncertainties of the data and the cross sections used in the one dimensional fluid model above 120 km, we cannot use MAVEN observations to prove that the electron temperature converges to the neutral temperature below 100 km. However, the lack of our understanding the electron temperature altitude profile below 120 km does not impact our understanding of the role of electron temperature in determining ion escape rates because ion escape is determined by electron temperatures above 180 km.