Local inversions are often observed in the heliospheric magnetic field (HMF), but their origins and evolution are not yet fully understood.Parker Solar Probe has recently observed rapid, Alfvenic, HMF inversions in the inner heliosphere, known as 'switchbacks', which have been interpreted as the possible remnants of coronal jets. It has also been suggested that inverted HMF may be produced by near-Sun interchange reconnection; a key process in mechanisms proposed for slow solar wind release. These cases suggest that the source of inverted HMF is near the Sun, and it follows that these inversions would gradually decay and straighten as they propagate out through the heliosphere. Alternatively, HMF inversions could form during solar wind transit, through phenomena such velocity shears, draping over ejecta, or waves and turbulence. Such processes are expected to lead to a qualitatively radial evolution of inverted HMF structures. Using Helios measurements spanning 0.3-1 AU, we examine the occurrence rate of inverted HMF, as well as other magnetic field morphologies, as a function of radial distance A, and find that it continually increases. This trend may be explained by inverted HMF observed between 0.3-1 AU being primarily driven by one or more of the above in-transit processes, rather than created at the Sun. We make suggestions as to the relative importance of these different processes based on the evolution of the magnetic field properties associated with inverted HMF. We also explore alternative explanations outside of our suggested driving processes which may lead to the observed trend.

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 paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

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.

Kinetic Alfv\'{e}n waves (KAWs) are the short-wavelength extension of the MHD Alfv\'{e}n-wave branch in the case of highly-oblique propagation with respect to the background magnetic field. Observations of space plasma show that small-scale turbulence is mainly KAW-like. We apply two theoretical approaches, collisional two-fluid theory and collisionless kinetic theory, to obtain predictions for the KAW polarizations depending on $\beta_\mathrm{p}$ (the ratio of the proton thermal pressure to the magnetic pressure) at the ion gyroscale in terms of fluctuations in density, bulk velocity, and pressure. We perform a wavelet analysis of MMS magnetosheath measurements and compare the observations with both theories. We find that the two-fluid theory predicts the observations better than kinetic theory, suggesting that the small-scale KAW-like fluctuations exhibit a fluid-like behavior in the magnetosheath although the plasma is weakly collisional. We also present predictions for the KAW polarizations in the inner heliosphere that are testable with Parker Solar Probe and Solar Orbiter.

The development of knowledge of how the coronal origin of the solar wind affects its in situ properties is one of the keys to understanding the relationship between the Sun and the heliosphere. In this paper, we analyse ACE/SWICS and WIND/3DP data spanning >12 years, and test properties of solar wind suprathermal electron distributions for the presence of signatures of the coronal temperature at their origin which may remain at 1AU. In particular we re-examine a previous suggestion that these properties correlate with the oxygen charge state ratio O7+/O6+; an established proxy for coronal electron temperature. We find only a very weak but variable correlation between measures of suprathermal electron energy content and O7+/O6+. The weak nature of the correlation leads us to conclude, in contrast to earlier results, that an initial relationship with core electron temperature has the possibility to exist in the corona, but that in most cases no strong signatures remain in the suprathermal electron distributions at 1AU. It can not yet be confirmed whether this is due to the effects of coronal conditions on the establishment of this relationship, or to the altering of the electron distributions by processing during transport in the solar wind en route to 1AU. Contrasting results for the halo and strahl population favours the latter interpretation. Confirmation of this will be possible using Solar Orbiter data (cruise and nominal mission phase) to test whether the weakness of the relationship persists over a range of heliocentric distances. If the correlation is found to strengthen when closer to the Sun, then this would indicate an initial relationship which is being degraded, perhaps by wave-particle interactions, en route to the observer

Observations in the solar wind suggest that the compressive component of inertial-range solar-wind turbulence is dominated by slow modes. The low collisionality of the solar wind allows for non-thermal features to survive, which suggests the requirement of a kinetic plasma description. The least-damped kinetic slow mode is associated with the ion-acoustic (IA) wave and a non-propagating (NP) mode. We derive analytical expressions for the IA-wave dispersion relation in an anisotropic plasma in the framework of gyrokinetics and then compare them to fully-kinetic numerical calculations, results from two-fluid theory, and MHD. This comparison shows major discrepancies in the predicted wave phase speeds from MHD and kinetic theory at moderate to high $\beta$. MHD and kinetic theory also dictate that all plasma normal modes exhibit a unique signature in terms of their polarization. We quantify the relative amplitude of fluctuations in the three lowest particle velocity moments associated with IA and NP modes in the gyrokinetic limit and compare these predictions with MHD results and in-situ observations of the solar-wind turbulence. The agreement between the observations of the wave polarization and our MHD predictions is better than the kinetic predictions, suggesting that the plasma behaves more like a fluid in the solar wind than expected.

We measure the local anisotropy of numerically simulated strong Alfv\'enic turbulence with respect to two local, physically relevant directions: along the local mean magnetic field and along the local direction of one of the fluctuating Elsasser fields. We find significant scaling anisotropy with respect to both these directions: the fluctuations are "ribbon-like" --- statistically, they are elongated along both the mean magnetic field and the fluctuating field. The latter form of anisotropy is due to scale-dependent alignment of the fluctuating fields. The intermittent scalings of the $n$th-order conditional structure functions in the direction perpendicular to both the local mean field and the fluctuations agree well with the theory of Chandran et al. 2015, while the parallel scalings are consistent with those implied by the critical-balance conjecture. We quantify the relationship between the perpendicular scalings and those in the fluctuation and parallel directions, and find that the scaling exponent of the perpendicular anisotropy (i.e., of the aspect ratio of the Alfv\'enic structures in the plane perpendicular to the mean magnetic field) depends on the amplitude of the fluctuations. This is shown to be equivalent to the anticorrelation of fluctuation amplitude and alignment at each scale. The dependence of the anisotropy on amplitude is shown to be more significant for the anisotropy between the perpendicular and fluctuation-direction scales than it is between the perpendicular and parallel scales.

Intensive studies have been conducted to understand the anisotropy of solar wind turbulence. However, the anisotropy of Els\"asser variables ($\textbf{Z}^\pm$) in 2D wave-vector space has yet to be investigated. Here we first verify the transformation based on the projection-slice theorem between the power spectral density PSD$_{2D}(k_\parallel,k_\perp )$ and the spatial correlation function CF$_{2D} (r_\parallel,r_\perp )$. Based on the application of the transformation to the magnetic field and the particle measurements from the WIND spacecraft, we investigate the spectral anisotropy of Els\"asser variables ($\textbf{Z}^\pm$), and the distribution of residual energy E$_{R}$, Alfv\'en ratio R$_{A}$ and Els\"asser ratio R$_{E}$ in the $(k_\parallel,k_\perp)$ space. The spectra PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$ (the larger of $\textbf{Z}^\pm$) show a similar pattern that PSD$_{2D}(k_\parallel,k_\perp )$ is mainly distributed along a ridge inclined toward the $k_\perp$ axis. This is probably the signature of the oblique Alfv\'enic fluctuations propagating outwardly. Unlike those of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$, the spectrum PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is distributed mainly along the $k_\perp$ axis. Close to the $k_\perp$ axis, $\left| {E}_{R}\right|$ becomes larger while R$_{A}$ becomes smaller, suggesting that the dominance of magnetic energy over kinetic energy becomes more significant at small $k_\parallel$. R$_{E}$ is larger at small $k_\parallel$, implying that PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is more concentrated along the $k_\perp$ direction as compared to that of $\textbf{Z}_{major}$. The residual energy condensate at small $k_\parallel$ is consistent with simulation results in which E$_{R}$ is spontaneously generated by Alfv\'en wave interaction.

We present a measurement of the scale-dependent, three-dimensional structure of the magnetic field fluctuations in inertial range solar wind turbulence with respect to a local, physically motivated coordinate system. The Alfvenic fluctuations are three-dimensionally anisotropic, with the sense of this anisotropy varying from large to small scales. At the outer scale, the magnetic field correlations are longest in the local fluctuation direction, consistent with Alfven waves. At the proton gyroscale, they are longest along the local mean field direction and shortest in the direction perpendicular to the local mean field and the local field fluctuation. The compressive fluctuations are highly elongated along the local mean field direction, although axially symmetric perpendicular to it. Their large anisotropy may explain why they are not heavily damped in the solar wind.

We present the first measurement of the scale-dependent power anisotropy of Elsasser variables in imbalanced fast solar wind turbulence. The dominant Elsasser mode is isotropic at lower spacecraft frequencies but becomes increasingly anisotropic at higher frequencies. The sub-dominant mode is anisotropic throughout, but in a scale-independent way (at higher frequencies). There are two distinct subranges exhibiting different scalings within what is normally considered the inertial range. The low Alfven ratio and shallow scaling of the sub-dominant Elsasser mode suggest an interpretation of the observed discrepancy between the velocity and magnetic field scalings. The total energy is dominated by the latter. These results do not appear to be fully explained by any of the current theories of incompressible imbalanced MHD turbulence.

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.

A relationship is derived between power anisotropy and wavevector anisotropy in turbulent fluctuations. This can be used to interpret plasma turbulence measurements, for example in the solar wind. If fluctuations are anisotropic in shape then the ion gyroscale break point in spectra in the directions parallel and perpendicular to the magnetic field would not occur at the same frequency, and similarly for the electron gyroscale break point. This is an important consideration when interpreting solar wind observations in terms of anisotropic turbulence theories. Model magnetic field power spectra are presented assuming a cascade of critically balanced Alfven waves in the inertial range and kinetic Alfven waves in the dissipation range. The variation of power anisotropy with scale is compared to existing solar wind measurements and the similarities and differences are discussed.

We measure the power and spectral index anisotropy of high speed solar wind turbulence from scales larger than the outer scale down to the ion gyroscale, thus covering the entire inertial range. We show that the power and spectral indices at the outer scale of turbulence are approximately isotropic. The turbulent cascade causes the power anisotropy at smaller scales manifested by anisotropic scalings of the spectrum: close to k^{-5/3} across and k^{-2} along the local magnetic field, consistent with a critically balanced Alfvenic turbulence. By using data at different radial distances from the Sun, we show that the width of the inertial range does not change with heliocentric distance and explain this by calculating the radial dependence of the ratio of the outer scale to the ion gyroscale. At the smallest scales of the inertial range, close to the ion gyroscale, we find an enhancement of power parallel to the magnetic field direction coincident with a decrease in the perpendicular power. This is most likely related to energy injection by ion kinetic modes such as the firehose instability and also marks the beginning of the dissipation range of solar wind turbulence.

We investigate the spatial correlation properties of the solar wind using simultaneous observations by the ACE and WIND spacecraft. We use mutual information as a nonlinear measure of correlation and compare this to linear correlation. We find that the correlation lengthscales of fluctuations in density and magnetic field magnitude vary strongly with the solar cycle, whereas correlation lengths of fluctuations in B field components do not. We find the correlation length of |B| ~ 120 Re at solar minimum and ~ 270 Re at maximum and the correlation length of density ~ 75 Re at minimum and ~ 170 Re at minimum. The components of the B field have correlation lengths ~ correlation length |B| at minimum.