A statistical survey of spectral anisotropy of space plasma turbulence is performed using five years measurements from MMS in the magnetosheath. By measuring the five-point second-order structure functions of the magnetic field, we have for the first time quantified the three-dimensional anisotropies and scalings at sub-ion-scales ($<$ 100 km). In the local reference frame $(\hat L_{\perp}, \hat l_{\perp}, \hat l_{\parallel})$ defined with respect to local mean magnetic field $\boldsymbol {B}_0$ (Chen et al. 2012), the "statistical eddies" are found to be mostly elongated along $\boldsymbol {B}_0$ and shortened in the direction perpendicular to both $\boldsymbol {B}_0$ and local field fluctuations. From several $d_i$ (ion inertial length) toward $\sim$ 0.05 $d_i$, the ratio between eddies' parallel and perpendicular lengths features a trend of rise then fall, whereas the anisotropy in the perpendicular plane appears scale-invariant. Specifically, the anisotropy relations for the total magnetic field at 0.1-1.0 $d_i$ are obtained as $l_{\parallel} \simeq 2.44 \cdot l_{\perp}^{0.71}$, and $L_{\perp} \simeq 1.58 \cdot l_{\perp}^{1.08}$, respectively. Our results provide new observational evidence to compare with phenomenological models and numerical simulations, which may help to better understand the nature of kinetic scale turbulence.

Turbulence develops in any stressed flow when the scales of the forcing are much larger than those of the dissipation. In neutral fluids, it consists of chaotic motions in physical space but with a universal energy spectrum in Fourier space. Intermittency (non-Gaussian statistics of fluctuations) is another general property and it is related to the presence of coherent structures. Space plasmas are turbulent as well. Here, we focus on the kinetic plasma scales, which are not yet well understood. We address the following fundamental questions: (1) Do the turbulent fluctuations at kinetic scales form a universal spectrum? and (2) What is the nature of the fluctuations? Using measurements in the solar wind we show that the magnetic spectra of kinetic turbulence at 0.3, 0.6 and 0.9 AU from the Sun have the same shape as the ones close to the Earth orbit at 1 AU, indicating universality of the phenomenon. The fluctuations, which form this spectrum, are typically non-linearly interacting eddies that tend to generate magnetic filaments.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

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.

We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations ({\delta}B) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations ({\delta}V) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of {\delta}B in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes.

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.

We study the anisotropy with respect to the local magnetic field of turbulent magnetic fluctuations at magnetofluid scales in the solar wind. Previous measurements in the fast solar wind obtained axisymmetric anisotropy, despite that the analysis method allows nonaxisymmetric structures. These results are probably contaminated by the wind expansion that introduces another symmetry axis, namely, the radial direction, as indicated by recent numerical simulations. These simulations also show that while the expansion is strong, the principal fluctuations are in the plane perpendicular to the radial direction. Using this property, we separate 11 yr of Wind spacecraft data into two subsets characterized by strong and weak expansion and determine the corresponding turbulence anisotropy. Under strong expansion, the small-scale anisotropy is consistent with the Goldreich & Sridhar critical balance. As in previous works, when the radial symmetry axis is not eliminated, the turbulent structures are field-aligned tubes. Under weak expansion, we find 3D anisotropy predicted by the Boldyrev model, that is, turbulent structures are ribbons and not tubes. However, the very basis of the Boldyrev phenomenology, namely, a cross-helicity increasing at small scales, is not observed in the solar wind: the origin of the ribbon formation is unknown.

We perform the first statistical study of the effects of the interaction of suprathermal electrons with narrow-band whistler mode waves in the solar wind. We show that this interaction does occur and that it is associated with enhanced widths of the so called strahl component. The latter is directed along the inter- planetary magnetic field away from the Sun. We do the study by comparing the strahl pitch angle widths in the solar wind at 1AU in the absence of large scale discontinuities and transient structures, such as interplanetary shocks, interplanetary coronal mass ejections, stream interaction regions, etc. during times when the whistler mode waves were present and when they were absent. This is done by using the data from two Cluster instruments: STAFF data in frequency range between ~0.1 Hz and ~200 Hz were used for determining the wave properties and PEACE datasets at twelve central energies between ~57 eV (equivalent to ~10 typical electron thermal energies in the solar wind, E_T ) and ~676 eV (~113 E_T ) for pitch angle measurements. Statistical analysis shows that during the inter- vals with the whistler waves the strahl component on average exhibits pitch angle widths between 2 and 12 degrees larger than during the intervals when these waves are not present. The largest difference is obtained for the electron central energy of ~344 eV (~57 E_T ).

In the solar wind, magnetic field power spectra usually show several power-laws. In this paper, magnetic field data from the Cluster mission during an undisturbed interval of slow solar wind is analyzed at 0.28Hz, near the spectral break point between the ion inertial and dissipation/dispersion ranges. Assuming Taylor's frozen-in condition, it corresponds to a proton kinetic scale of $kv_A/\Omega_p \sim 0.38$, where $v_A$ and $\Omega_p$ are the Alfv\'en speed and proton angular gyrofrequency, respectively. Data show that the Cluster spacecraft passed through a series of wavepackets. A strong isolated wavepacket is found to be in accordance with the four Cluster satellites crossing an Alfv\'en vortex, a nonlinear solution to the incompressible MHD equations. A strong agreement is seen between the data from four satellites and a model vortex with a radius of the order of $40$ times the local proton gyro-radii. The polarization at different spacecraft is compared and is found to agree with the vortex model, whereas it cannot be explained solely by the linear plane wave approach.

In this paper we investigate spectral and phase coherence properties of magnetic fluctuations in the vicinity of the spectral transition from large, magnetohydrodynamic (MHD) to sub-ion scales using in-situ measurements of the Wind spacecraft in a fast stream. For the time interval investigated by Leamon et al. (1998) the phase-coherence analysis shows the presence of sporadic quasi-parallel Alfv\'en Ion Cyclotron (AIC) waves as well as coherent structures in the form of large-amplitude, quasi-perpendicular Alfv\'en vortex-like structures and current sheets. These waves and structures importantly contribute to the observed power spectrum of magnetic fluctuations around ion scales; AIC waves contribute to the spectrum in a narrow frequency range whereas the coherent structures contribute to the spectrum over a wide frequency band from the inertial range to the sub-ion frequency range. We conclude that a particular combination of waves and coherent structures determines the spectral shape of the magnetic field spectrum around ion scales. This phenomenon provides a possible explanation for a high variability of the magnetic power spectra around ion scales observed in the solar wind.

The nature of the magnetic field fluctuations in the solar wind between the ion and electron scales is still under debate. Using the Cluster/STAFF instrument, we make a survey of the power spectral density and of the polarization of these fluctuations at frequencies $f\in[1,400]$ Hz, during five years (2001-2005), when Cluster was in the free solar wind. In $\sim 10\%$ of the selected data, we observe narrow-band, right-handed, circularly polarized fluctuations, with wave vectors quasi-parallel to the mean magnetic field, superimposed on the spectrum of the permanent background turbulence. We interpret these coherent fluctuations as whistler mode waves. The life time of these waves varies between a few seconds and several hours. Here we present, for the first time, an analysis of long-lived whistler waves, i.e. lasting more than five minutes. We find several necessary (but not sufficient) conditions for the observation of whistler waves, mainly a low level of the background turbulence, a slow wind, a relatively large electron heat flux and a low electron collision frequency. When the electron parallel beta factor $\beta_{e\parallel}$ is larger than 3, the whistler waves are seen along the heat flux threshold of the whistler heat flux instability. The presence of such whistler waves confirms that the whistler heat flux instability contributes to the regulation of the solar wind heat flux, at least for $\beta_{e\parallel} \ge$ 3, in the slow wind, at 1 AU.

Solar wind is probably the best laboratory to study turbulence in astrophysical plasmas. In addition to the presence of magnetic field, the differences with neutral fluid isotropic turbulence are: weakness of collisional dissipation and presence of several characteristic space and time scales. In this paper we discuss observational properties of solar wind turbulence in a large range from the MHD to the electron scales. At MHD scales, within the inertial range, turbulence cascade of magnetic fluctuations develops mostly in the plane perpendicular to the mean field. Solar wind turbulence is compressible in nature. The spectrum of velocity fluctuations do not follow magnetic field one. Probability distribution functions of different plasma parameters are not Gaussian, indicating presence of intermittency. At the moment there is no global model taking into account all these observed properties of the inertial range. At ion scales, turbulent spectra have a break, compressibility increases and the density fluctuation spectrum has a local flattening. Around ion scales, magnetic spectra are variable and ion instabilities occur as a function of the local plasma parameters. Between ion and electron scales, a small scale turbulent cascade seems to be established. Approaching electron scales, the fluctuations are no more self-similar: an exponential cut-off is usually observed indicating an onset of dissipation. The nature of the small scale cascade and a possible dissipation mechanism are still under debate.

The description of the turbulent spectrum of magnetic fluctuations in the solar wind in the kinetic range of scales is not yet completely established. Here, we perform a statistical study of 100 spectra measured by the STAFF instrument on the Cluster mission, which allows to resolve turbulent fluctuations from ion scales down to a fraction of electron scales, i.e. from $\sim 10^2$ km to $\sim 300$ m. We show that for $k_{\perp}\rho_e \in[0.03,3]$ (that corresponds approximately to the frequency in the spacecraft frame $f\in [3,300]$ Hz), all the observed spectra can be described by a general law $E(k_\perp)\propto k_\perp^{-8/3}\exp{(-k_\perp \rho_e)}$, where $k_{\perp}$ is the wave-vector component normal to the background magnetic field and $\rho_e$ the electron Larmor radius. This exponential tail found in the solar wind seems compatible with the Landau damping of magnetic fluctuations onto electrons.

The turbulent spectrum of magnetic fluctuations in the solar wind displays a spectral break at ion characteristic scales. At electron scales the spectral shape is not yet completely established. Here, we perform a statistical study of 102 spectra at plasma kinetic scales, measured by the Cluster/STAFF instrument in the free solar wind. We show that the magnetic spectrum in the high frequency range, [1,400] Hz, has a form similar to what is found in hydrodynamics in the dissipation range ~Ak^(-\alpha)exp(-kl_d). The dissipation scale l_d is found to be correlated with the electron Larmor radius \rho_e. The spectral index \alpha varies in the range [2.2,2.9] and is anti-correlated with l_d, as expected in the case of the balance between the energy injection and the energy dissipation. The coefficient A is found to be proportional to the ion temperature anisotropy, suggesting that local ion instabilities may play some role for the solar wind turbulence at plasma kinetic scales. The exponential spectral shape found here indicates that the effective dissipation of magnetic fluctuations in the solar wind has a wave number dependence similar to that of the resistive term in collisional fluids ~k^2\delta B.

The anisotropy of turbulence in the fast solar wind, between the ion and electron gyroscales, is directly observed using a multispacecraft analysis technique. Second order structure functions are calculated at different angles to the local magnetic field, for magnetic fluctuations both perpendicular and parallel to the mean field. In both components, the structure function value at large angles to the field S_perp is greater than at small angles S_par: in the perpendicular component S_perp/S_par = 5 +- 1 and in the parallel component S_perp/S_par > 3, implying spatially anisotropic fluctuations, k_perp > k_par. The spectral index of the perpendicular component is -2.6 at large angles and -3 at small angles, in broad agreement with critically balanced whistler and kinetic Alfven wave predictions. For the parallel component, however, it is shallower than -1.9, which is considerably less steep than predicted for a kinetic Alfven wave cascade.

We investigate the spectral shape, the anisotropy of the wave vector distributions and the anisotropy of the amplitudes of the magnetic fluctuations in the Earth's magnetosheath within a broad range of frequencies. We present the first observations of a Kolmogorov-like inertial range of Alfvenic fluctuations in the magnetosheath flanks, below fci. In the vicinity of fci, a spectral break is observed, like in solar wind turbulence. Above the break, the energy of compressive and Alfvenic fluctuations generally follow a power law with a spectral index between -3 and -2. Concerning the anisotropy of the wave vector distribution, we observe a change in its nature in the vicinity of ion characteristic scales: if at MHD scales there is no evidence for a dominance of a slab (k||>kperp) or 2D (kperp>k||) turbulence, above the spectral break, (f>fci, kc/wpi>1) the 2D turbulence dominates. This 2D turbulence is observed in six selected one-hour intervals among which the average proton beta varies from 0.8 to 9. It is observed for both the transverse and compressive magnetic fluctuations, independently on the presence of linearly unstable modes at low frequencies or Alfven vortices at the spectral break. We then analyse the anisotropy of the magnetic fluctuations in a time dependent reference frame based on the field B and the flow velocity V directions. Within the range of the 2D turbulence, at scales [1,30]kc/wpi, and for any beta we find that the magnetic fluctuations at a given frequency in the plane perpendicular to B have more energy along the BxV direction. This non-gyrotropy of the fluctuations is consistent with gyrotropic fluctuations at a given wave vector, with kperp>k||, which suffer a different Doppler shift along and perpendicular to V in the plane perpendicular to B.

Magnetic fluctuations in the solar wind are distributed according to Kolmogorov's power law $f^{-5/3}$ below the ion cyclotron frequency $f_{ci}$. Above this frequency, the observed steeper power law is usually interpreted in two different ways: a dissipative range of the solar wind turbulence or another turbulent cascade, the nature of which is still an open question. Using the Cluster magnetic data we show that after the spectral break the intermittency increases toward higher frequencies, indicating the presence of non-linear interactions inherent to a new inertial range and not to the dissipative range. At the same time the level of compressible fluctuations raises. We show that the energy transfer rate and intermittency are sensitive to the level of compressibility of the magnetic fluctuations within the small scale inertial range. We conjecture that the time needed to establish this inertial range is shorter than the eddy-turnover time, and is related to dispersive effects. A simple phenomenological model, based on the compressible Hall MHD, predicts the magnetic spectrum $\sim k^{-7/3+2\alpha}$, which depends on the degree of plasma compression $\alpha$.