The first complete estimation of the compressible energy cascade rate $|\varepsilon_\text{C}|$ at magnetohydrodynamic (MHD) and sub-ion scales is obtained in the Earth's magnetosheath using Magnetospheric MultiScale (MMS) spacecraft data and an exact law derived recently for {\it compressible} Hall MHD turbulence. A multi-spacecraft technique is used to compute the velocity and magnetic gradients, and then all the correlation functions involved in the exact relation. It is shown that when the density fluctuations are relatively small, $|\varepsilon_\text{C}|$ identifies well with its incompressible analogue $|\varepsilon_\text{I}|$ at MHD scales but becomes much larger than $|\varepsilon_\text{I}|$ at sub-ion scales. For larger density fluctuations, $|\varepsilon_\text{C}|$ is larger than $|\varepsilon_\text{I}|$ at every scale with a value significantly higher than for smaller density fluctuations. Our study reveals also that for both small and large density fluctuations, the non-flux terms remain always negligible with respect to the flux terms and that the major contribution to $|\varepsilon_\text{C}|$ at sub-ion scales comes from the compressible Hall flux.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

Three-dimensional direct numerical simulations are used to study the energy cascade rate in isothermal compressible magnetohydrodynamic turbulence. Our analysis is guided by a two-point exact law derived recently for this problem in which flux, source, hybrid, and mixed terms are present. The relative importance of each term is studied for different initial subsonic Mach numbers $M_S$ and different magnetic guide fields ${\bf B}_0$. The dominant contribution to the energy cascade rate comes from the compressible flux, which depends weakly on the magnetic guide field ${\bf B}_0$, unlike the other terms whose modulus increase significantly with $M_S$ and ${\bf B}_0$. In particular, for strong ${\bf B}_0$ the source and hybrid terms are dominant at small scales with almost the same amplitude but with a different sign. A statistical analysis made with an isotropic decomposition based on the SO(3) rotation group is shown to generate spurious results in presence of ${\bf B}_0$, when compared with an axisymmetric decomposition better suited to the geometry of the problem. Our numerical results are compared with previous analyses made with in-situ measurements in the solar wind and the terrestrial magnetosheath.

The first estimation of the energy cascade rate ${|\epsilon_C|}$ of magnetosheath turbulence is obtained using the CLUSTER and THEMIS spacecraft data and an exact law of compressible isothermal magnetohydrodynamics turbulence. ${|\epsilon_C|}$ is found to be of the order of ${10{^{-13}} J.m^{3}.s^{-1}}$, at least two orders of magnitude larger than its value in the solar wind (order of ${10{^{-16}} J.m^{3}.s^{-1}}$ in the fast wind). Two types of turbulence are evidenced and shown to be dominated either by incompressible Alfv\'enic or magnetosonic-like fluctuations. Density fluctuations are shown to amplify the cascade rate and its spatial anisotropy in comparison with incompressible Alfv\'enic turbulence. Furthermore, for compressible magnetosonic fluctuations, large cascade rates are found to lie mostly near the linear kinetic instability of the mirror mode. New empirical power-laws are evidenced and relate ${|\epsilon_C|}$ to the turbulent Mach number and the internal energy. These new finding have potential applications in distant astrophysical plasmas that are not accessible to in situ measurements.

We derive the exact law for three-dimensional (3D) homogeneous compressible isothermal Hall magnetohydrodynamics (CHMHD) turbulence, without the assumption of isotropy. The Hall current is shown to introduce new flux and sources terms that act at the small scales (comparable or smaller than the ion skin depth) to significantly impact the turbulence dynamics. The new law provides an accurate means to estimate for the first time the energy cascade rate over a broad range of scales covering the MHD inertial range and the sub-ion dispersive range in 3D numerical simulations and {\it in situ} spacecraft observations of compressible turbulence. This work is particularly relevant to astrophysical flows in which small scale density fluctuations cannot be ignored such as the solar wind, planetary magnetospheres and the interstellar medium (ISM).

We derive exact scaling laws for a three-dimensional incompressible helical two-fluid plasma, without the assumption of isotropy. For each ideal invariant of the two-fluid model, i.e. the total energy, the electron helicity and the proton helicity, we derive simple scaling laws in terms of two-point increments correlation functions expressed in terms of the velocity field of each species and the magnetic field. These variables are appropriate for comparison with \textit{in-situ} measurements in the solar wind at different spatial ranges and data from numerical simulations. Finally, with the exact scaling laws and dimensional analysis we predict the magnetic energy and electron helicity spectra for different ranges of scales.

In the solar wind, power spectral density (PSD) of the magnetic field fluctuations generally follow the so-called Kolmogorov spectrum f^-5/3 in the inertial range, where the dynamics is thought to be dominated by nonlinear interactions between counter-propagating incompressible Alfv\'en wave parquets. These features are thought to be ubiquitous in space plasmas. The present study gives a new and more complex picture of magnetohydrodynamics (MHD) turbulence as observed in the terrestrial magnetosheath. The study uses three years of in-situ data from the Cluster mission to explore the nature of the magnetic fluctuations at MHD scales in different locations within the magnetosheath, including flanks and subsolar regions. It is found that the magnetic field fluctuations at MHD scales generally have a PSD close to f^-1 (shallower than the Kolmogorov one f^-5/3) down to the ion characteristic scale, which recalls the energy containing scales of solar wind turbulence. The Kolmogorov spectrum is observed only away from the bow shock toward the flank and the magnetopause regions in 17% of the analyzed time intervals. Measuring the magnetic compressibility, it is shown that only a fraction (35%) of the observed Kolmogorov spectra were populated by shear Alfv\'enic fluctuations, whereas the majority of the events (65%) was found to be dominated by compressible magnetosonic-like fluctuations, which contrasts with well-known turbulence properties in the solar wind. This study gives a first comprehensive view of the origin of the f^-1 and the transition to the Kolmogorov inertial range; both questions remain controversial in solar wind turbulence.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

Estimation of the energy cascade rate in the inertial range of solar wind turbulence has been done so far mostly within the incompressible magnetohydrodynamics (MHD) theory. Here, we go beyond that approximation to include plasma compressibility using a reduced form of a recently derived exact law for compressible, isothermal MHD turbulence. Using in-situ data from the THEMIS/ARTEMIS spacecraft in the fast and slow solar wind, we investigate in detail the role of the compressible fluctuations in modifying the energy cascade rate with respect to the prediction of the incompressible MHD model. In particular, we found that the energy cascade rate: i) is amplified particularly in the slow solar wind; ii) exhibits weaker fluctuations in spatial scales, which leads to a broader inertial range than the previous reported ones; iii) has a power law scaling with the turbulent Mach number; iv) has a lower level of spatial anisotropy. Other features of solar wind turbulence are discussed along with their comparison with previous studies that used incompressible or heuristic (non exact) compressible MHD models.

Low frequency turbulence in Saturn's magnetosheath is investigated using in-situ measurements of the Cassini spacecraft. Focus is put on the magnetic energy spectra computed in the frequency range $\sim[10^{-4}, 1]$Hz. A set of 42 time intervals in the magnetosheath were analyzed and three main results that contrast with known features of solar wind turbulence are reported: 1) The magnetic energy spectra showed a $\sim f^{-1}$ scaling at MHD scales followed by an $\sim f^{-2.6}$ scaling at the sub-ion scales without forming the so-called inertial range; 2) The magnetic compressibility and the cross-correlation between the parallel component of the magnetic field and density fluctuations $ C(\delta n,\delta B_{||}) $ indicates the dominance of the compressible magnetosonic slow-like modes at MHD scales rather than the Alfv\'en mode; 3) Higher order statistics revealed a monofractal (resp. multifractal) behaviour of the turbulent flow downstream of a quasi-perpendicular (resp. quasi-parallel) shock at the sub-ion scales. Implications of these results on theoretical modeling of space plasma turbulence are discussed.

We present a first statistical study of subproton and electron scales turbulence in the terrestrial magnetosheath using the Cluster Search Coil Magnetometer (SCM) waveforms of the STAFF instrument measured in the frequency range [1,180] Hz. It is found that clear spectral breaks exist near the electron scale, which separate two power-law like frequency bands referred to as the dispersive and the electron dissipation ranges. The frequencies of the breaks f_b are shown to be well correlated with the electron gyroscale \rho_e rather than with the electron inertial length de. The distribution of the slopes below fb was found to be narrow and peaks near -2.9, while that of the slopes above fb was found broader, peaks near -5.2 and has values as low as -7.5. This is the first time that such steep power-law spectra are reported in space plasma turbulence. These observations provide strong constraints on theoretical modeling of kinetic turbulence and dissipation in collisionless magnetized plasmas.

Electron scale solar wind turbulence has attracted great interest in recent years. Clear evidences have been given from the Cluster data that turbulence is not fully dissipated near the proton scale but continues cascading down to the electron scales. However, the scaling of the energy spectra as well as the nature of the plasma modes involved at those small scales are still not fully determined. Here we survey 10 years of the Cluster search-coil magnetometer (SCM) waveforms measured in the solar wind and perform a statistical study of the magnetic energy spectra in the frequency range [$1, 180$]Hz. We show that a large fraction of the spectra exhibit clear breakpoints near the electon gyroscale $\rho_e$, followed by steeper power-law like spectra. We show that the scaling below the electron breakpoint cannot be determined unambiguously due to instrumental limitations that will be discussed in detail. We compare our results to recent ones reported in other studies and discuss their implication on the physical mechanisms and the theoretical modeling of energy dissipation in the SW.

The nature of solar wind (SW) turbulence below the proton gyroscale is a topic that is being investigated extensively nowadays. Although recent observations gave evidence of the dominance of Kinetic Alfv\'en Waves (KAW) at sub-ion scales with $\omega<{\omega_{ci}}$, other studies suggest that the KAW mode cannot carry the turbulence cascade down to electron scales and that the whistler mode (i.e., $\omega>\omega_{ci}$) is more relevant. Here, we propose to study key properties of the short wavelength plasma modes under realistic SW conditions, typically $\beta_i\gtrsim \beta_e\sim 1$ and for high oblique angles of propagation $80^\circ\leq \Theta_{\bf kB}<90^\circ$ as observed from the Cluster data. The linear properties of the plasma modes under these conditions are poorly known, which contrasts with the well-documented cold plasma limit and/or moderate oblique angles of propagation ($\Theta_{\bf kB} <80^\circ$). Based on linear solutions of the Vlasov kinetic theory, we discuss the relevance of each plasma mode (fast, Bernstein, KAW, whistler) in carrying the energy cascade down to electron scales. We show, in particular, that the shear Alfv\'en mode extends at scales $k\rho_i\gtrsim1$ following either a whistler mode ($\omega>\omega_{ci}$) or a KAW mode (with $\omega<\omega_{ci}$) depending on the anisotropy $k_\parallel/ k_\perp$. This contrasts with the well-accepted idea that the whistler branch develops as a continuation at high frequencies of the fast magnetosonic mode. We show, furthermore, that the whistler branch is more damped than the KAW one, which makes the latter a more relevant candidate to carry the energy cascade down to electron scales. We discuss how these new findings may facilitate resolution of the controversy concerning the nature of the small scale turbulence, and we discuss the implications for present and future spacecraft wave measurements in the SW.