First results from the Parker Solar Probe (PSP) mission have revealed ubiquitous coherent ion-scale waves in the inner heliosphere, which are signatures of kinetic wave-particle interactions and fluid-scale instabilities. However, initial studies of the circularly polarized ion-scale waves observed by PSP have only thoroughly analyzed magnetic field signatures, precluding a determination of solar-wind frame propagation direction and intrinsic wave-polarization. A comprehensive determination of wave-properties requires measurements of both electric and magnetic fields. Here, we use full capabilities of the PSP/FIELDS instrument suite to measure both the electric and magnetic components of circularly polarized waves. Comparing spacecraft frame magnetic field measurements with the Doppler-shifted cold-plasma dispersion relation for parallel transverse waves constrains allowable plasma frame polarizations and wave-vectors. We demonstrate that the Doppler-shifted cold-plasma dispersion has a maximum spacecraft frequency $f_{sc}^{*}$ for which intrinsically right-handed fast-magnetosonic waves (FMWs) propagating sunwards can appear left-handed in the spacecraft frame. Observations of left-handed waves with $|f|>f_{sc}^{*}$ are uniquely explained by intrinsically left-handed, ion-cyclotron, waves (ICWs). We demonstrate that electric field measurements for waves with $|f|>f_{sc}^{*}$ are consistent with ICWs propagating away from the sun, verifying the measured electric field. Applying the verified electric field measurements to the full distribution of waves suggests that, in the solar wind frame, the vast majority of waves propagate away from the sun, indicating that the observed population of coherent ion-scale waves contains both intrinsically left and right hand polarized modes.

We perform a statistical study of the turbulent power spectrum at inertial and kinetic scales observed during the first perihelion encounter of Parker Solar Probe. We find that often there is an extremely steep scaling range of the power spectrum just above the ion-kinetic scales, similar to prior observations at 1 AU, with a power-law index of around $-4$. Based on our measurements, we demonstrate that either a significant ($>50\%$) fraction of the total turbulent energy flux is dissipated in this range of scales, or the characteristic nonlinear interaction time of the turbulence decreases dramatically from the expectation based solely on the dispersive nature of nonlinearly interacting kinetic Alfv\'en waves.

We consider 2D joint distributions of normalised residual energy $\sigma_r(s,t)$ and cross helicity $\sigma_c(s,t)$ during one day of Parker Solar Probe's (PSP's) first encounter as a function of wavelet scale $s$. The broad features of the distributions are similar to previous observations made by HELIOS in slow solar wind, namely well correlated and fairly Alfv\'enic, except for a population with negative cross helicity which is seen at shorter wavelet scales. We show that this population is due to the presence of magnetic switchbacks, brief periods where the magnetic field polarity reverses. Such switchbacks have been observed before, both in HELIOS data and in Ulysses data in the polar solar wind. Their abundance and short timescales as seen by PSP in its first encounter is a new observation, and their precise origin is still unknown. By analysing these MHD invariants as a function of wavelet scale we show that MHD waves do indeed follow the local mean magnetic field through switchbacks, with net Elsasser flux propagating inward during the field reversal, and that they therefore must be local kinks in the magnetic field and not due to small regions of opposite polarity on the surface of the Sun. Such observations are important to keep in mind as computing cross helicity without taking into account the effect of switchbacks may result in spurious underestimation of $\sigma_c$ as PSP gets closer to the Sun in later orbits.

Understanding the physical processes in the solar wind and corona which actively contribute to heating, acceleration, and dissipation is a primary objective of NASA's Parker Solar Probe (PSP) mission. Observations of coherent electromagnetic waves at ion scales suggests that linear cyclotron resonance and non-linear processes are dynamically relevant in the inner heliosphere. A wavelet-based statistical study of coherent waves in the first perihelion encounter of PSP demonstrates the presence of transverse electromagnetic waves at ion resonant scales which are observed in 30-50\% of radial field intervals. Average wave amplitudes of approximately 4 nT are measured, while the mean duration of wave events is of order 20 seconds; however long duration wave events can exist without interruption on hour-long timescales. Though ion scale waves are preferentially observed during intervals with a radial mean magnetic field, we show that measurement constraints, associated with single spacecraft sampling of quasi-parallel waves superposed with anisotropic turbulence, render the measured quasi-parallel ion-wave spectrum unobservable when the mean magnetic field is oblique to the solar wind flow; these results imply that the occurrence of coherent ion-scale waves is not limited to a radial field configuration. The lack of strong radial scaling of characteristic wave amplitudes and duration suggests that the waves are generated {\em{in-situ}} through plasma instabilities. Additionally, observations of proton distribution functions indicate that temperature anisotropy may drive the observed ion-scale waves.

The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfv\'enic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfv\'en speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfv\'en point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

The question of the relative importance of coherent structures and waves has for a long time attracted a great deal of interest in astrophysical plasma turbulence research, with a more recent focus on kinetic scale dynamics. Here we utilize high-resolution observational and simulation data to investigate the nature of waves and structures emerging in a weakly collisional, turbulent kinetic plasma. Observational results are based on in situ solar wind measurements from the Cluster and MMS spacecraft, and the simulation results are obtained from an externally driven, three-dimensional fully kinetic simulation. Using a set of novel diagnostic measures we show that both the large-amplitude structures and the lower-amplitude background fluctuations preserve linear features of kinetic Alfven waves to order unity. This quantitative evidence suggests that the kinetic turbulence cannot be described as a mixture of mutually exclusive waves and structures but may instead be pictured as an ensemble of localized, anisotropic wave packets or "eddies" of varying amplitudes, which preserve certain linear wave properties during their nonlinear evolution.

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.

Stochastic heating refers to an increase in the average magnetic moment of charged particles interacting with electromagnetic fluctuations whose frequencies are smaller than the particles' cyclotron frequencies. This type of heating arises when the amplitude of the gyroscale fluctuations exceeds a certain threshold, causing particle orbits in the plane perpendicular to the magnetic field to become stochastic rather than nearly periodic. We consider the stochastic heating of protons by Alfv\'en-wave (AW) and kinetic-Alfv\'en-wave (KAW) turbulence, which may make an important contribution to the heating of the solar wind. Using phenomenological arguments, we derive the stochastic-proton-heating rate in plasmas in which $\beta_{\rm p} \sim 1-30$, where $\beta_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. (We do not consider the $\beta_{\rm p} \gtrsim 30$ regime, in which KAWs at the proton gyroscale become non-propagating.) We test our formula for the stochastic-heating rate by numerically tracking test-particle protons interacting with a spectrum of randomly phased AWs and KAWs. Previous studies have demonstrated that at $\beta_{\rm p} \lesssim 1$, particles are energized primarily by time variations in the electrostatic potential and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the electrostatic potential. In contrast, at $\beta_{\rm p} \gtrsim 1$, particles are energized primarily by the solenoidal component of the electric field and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the magnetic field.

It is widely reported that the power spectra of magnetic field and velocity fluctuations in the solar wind have power law scalings with inertial-range spectral indices of -5/3 and -3/2 respectively. Studies of solar wind turbulence have repeatedly demonstrated the impact of discontinuities and coherent structures on the measured spectral index. Whether or not such discontinuities are self-generated by the turbulence or simply observations of advected structures from the inner heliosphere has been a matter of considerable debate. This work presents a statistical study of magnetic field and velocity spectral indices over 10 years of solar-wind observations; we find that anomalously steep magnetic spectra occur in magnetically dominated intervals with negative residual energy. However, this increase in negative residual energy has no noticeable impact on the spectral index of the velocity fluctuations, suggesting that these intervals with negative residual energy correspond to intermittent magnetic structures. We show statistically that the difference between magnetic and velocity spectral indices is a monotonic function of residual energy, consistent with previous work which suggests that intermittency in fluctuations causes spectral steepening. Additionally, a statistical analysis of cross helicity demonstrates that when the turbulence is balanced (low cross-helicity), the magnetic and velocity spectral indices are not equal, which suggests that our observations of negative residual energy and intermittent structures are related to non-linear turbulent interactions rather than the presence of advected pre-existing flux-tube structures.

We present results from a three-dimensional particle-in-cell simulation of plasma turbulence, resembling the plasma conditions found at kinetic scales of the solar wind. The spectral properties of the turbulence in the subion range are consistent with theoretical expectations for kinetic Alfv\' en waves. Furthermore, we calculate the local anisotropy, defined by the relation $k_{\parallel}(k_{\perp})$, where $k_{\parallel}$ is a characteristic wave number along the local mean magnetic field at perpendicular scale $l_{\perp}\sim 1/k_{\perp}$. The subion range anisotropy is scale dependent with $k_{\parallel}<k_{\perp}$ and the ratio of linear to nonlinear time scales is of order unity, suggesting that the kinetic cascade is close to a state of critical balance. Our results compare favorably against a number of \emph{in situ} solar wind observations and demonstrate---from first principles---the feasibility of plasma turbulence models based on a critically balanced cascade of kinetic Alfv\' en waves.

We calculate the disruption scale $\lambda_{\rm D}$ at which sheet-like structures in dynamically aligned Alfv\'enic turbulence are destroyed by the onset of magnetic reconnection in a low-$\beta$ collisionless plasma. The scaling of $\lambda_{\rm D}$ depends on the order of the statistics being considered, with more intense structures being disrupted at larger scales. The disruption scale for the structures that dominate the energy spectrum is $\lambda_{\rm D}\sim L_\perp^{1/9}(d_e\rho_s)^{4/9}$, where $d_e$ is the electron inertial scale, $\rho_s$ is the ion sound scale, and $L_\perp$ is the outer scale of the turbulence. When $\beta_e$ and $\rho_s/L_\perp$ are sufficiently small, the scale $\lambda_{\rm D}$ is larger than $\rho_s$ and there is a break in the energy spectrum at $\lambda_{\rm D}$, rather than at $\rho_s$. We propose that the fluctuations produced by the disruption are circularised flux ropes, which may have already been observed in the solar wind. We predict the relationship between the amplitude and radius of these structures and quantify the importance of the disruption process to the cascade in terms of the filling fraction of undisrupted structures and the fractional reduction of the energy contained in them at the ion sound scale $\rho_s$. Both of these fractions depend strongly on $\beta_e$, with the disrupted structures becoming more important at lower $\beta_e$. Finally, we predict that the energy spectrum between $\lambda_{\rm D}$ and $\rho_s$ is steeper than $k_\perp^{-3}$, when this range exists. Such a steep "transition range" is sometimes observed in short intervals of solar-wind turbulence. The onset of collisionless magnetic reconnection may therefore significantly affect the nature of plasma turbulence around the ion gyroscale.

We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfv\'enic turbulence are destroyed by magnetic reconnection at a scale $\hat{\lambda}_{\rm D}$, larger than the dissipation scale predicted by models of intermittent, dynamically aligning turbulence. The reconnection process proceeds in several stages: first, a linear tearing mode with $N$ magnetic islands grows and saturates, and then the $X$-points between these islands collapse into secondary current sheets, which then reconnect until the original structure is destroyed. This effectively imposes an upper limit on the anisotropy of the structures within the perpendicular plane, which means that at scale $\hat{\lambda}_{\rm D}$ the turbulent dynamics change: at scales larger than $\hat{\lambda}_{\rm D}$, the turbulence exhibits scale-dependent dynamic alignment and a spectral index approximately equal to $-3/2$, while at scales smaller than $\hat{\lambda}_{\rm D}$, the turbulent structures undergo a succession of disruptions due to reconnection, limiting dynamic alignment, steepening the effective spectral index and changing the final dissipation scale. The scaling of $\hat{\lambda}_{\rm D}$ with the Lundquist (magnetic Reynolds) number $S_{L_\perp}$ depends on the order of the statistics being considered, and on the specific model of intermittency; the transition between the two regimes in the energy spectrum is predicted at approximately $\hat{\lambda}_{\rm D} \sim S_{L_\perp}^{-0.6}$. The spectral index below $\hat{\lambda}_{\rm D}$ is bounded between $-5/3$ and $-2.3$. The final dissipation scale is at $\hat{\lambda}_{\eta,\infty}\sim S_{L_\perp}^{-3/4}$, the same as the Kolmogorov scale arising in theories of turbulence that do not involve scale-dependent dynamic alignment.

We propose a simple statistical model of three-dimensionally anisotropic, intermittent, strong Alfv\'enic turbulence, incorporating both critical balance and dynamic alignment. Our model is based on log-Poisson statistics for Elsasser-field increments along the magnetic field. We predict the scalings of Elsasser-field conditional two-point structure functions with point separations in all three directions in a coordinate system locally aligned with the direction of the magnetic field and of the fluctuating fields and obtain good agreement with numerical simulations. We also derive a scaling of the parallel coherence scale of the fluctuations, $l_\parallel \propto \lambda^{1/2}$, where $\lambda$ is the perpendicular scale. This is indeed observed for the bulk of the fluctuations in numerical simulations.

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.

The interaction between Alfven-wave turbulence and the background solar wind affects the cross helicity in two ways. Non-WKB reflection converts outward-propagating Alfven waves into inward-propagating Alfven waves and vice versa, and the turbulence transfers momentum to the background flow. When both effects are accounted for, the total cross helicity is conserved. In the special case that the background density and flow speed are independent of time, the equations of cross-helicity conservation and total-energy conservation can be combined to recover a well-known equation derived by Heinemann and Olbert that has been interpreted as a non-WKB generalization of wave-action conservation. This latter equation (in contrast to cross-helicity and energy conservation) does not hold when the background varies in time.

We present numerical evidence that in strong Alfvenic turbulence, the critical balance principle---equality of the nonlinear decorrelation and linear propagation times---is scale invariant, in the sense that the probability distribution of the ratio of these times is independent of scale. This result only holds if the local alignment of the Elsasser fields is taken into account in calculating the nonlinear time. At any given scale, the degree of alignment is found to increase with fluctuation amplitude, supporting the idea that the cause of alignment is mutual dynamical shearing of Elsasser fields. The scale-invariance of critical balance (while all other quantities of interest are strongly intermittent, i.e., have scale-dependent distributions) suggests that it is the most robust of the scaling principles used to describe Alfvenic turbulence. The quality afforded by situ fluctuation measurements in the solar wind allows for direct verification of this fundamental principle.

We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.

The Solar Probe Plus (SPP) spacecraft will explore the near-Sun environment, reaching heliocentric distances less than $10 R_{\odot}$. Near Earth, spacecraft measurements of fluctuating velocities and magnetic fields taken in the time domain are translated into information about the spatial structure of the solar wind via Taylor's "frozen turbulence" hypothesis. Near the perihelion of SPP, however, the solar-wind speed is comparable to the Alfv\'en speed, and Taylor's hypothesis in its usual form does not apply. In this paper, we show that, under certain assumptions, a modified version of Taylor's hypothesis can be recovered in the near-Sun region. We consider only the transverse, non-compressive component of the fluctuations at length scales exceeding the proton gyroradius, and we describe these fluctuations using an approximate theoretical framework developed by Heinemann and Olbert. We show that fluctuations propagating away from the Sun in the plasma frame obey a relation analogous to Taylor's hypothesis when $V_{\rm sc,\perp} \gg z^-$ and $z^+ \gg z^-$, where $V_{\rm sc,\perp}$ is the component of the spacecraft velocity perpendicular to the mean magnetic field and $\bm{z}^+$ ($\bm{z}^-$) is the Elsasser variable corresponding to transverse, non-compressive fluctuations propagating away from (towards) the Sun in the plasma frame. Observations and simulations suggest that, in the near-Sun solar wind, the above inequalities are satisfied and $\bm{z}^+$ fluctuations account for most of the fluctuation energy. The modified form of Taylor's hypothesis that we derive may thus make it possible to characterize the spatial structure of the energetically dominant component of the turbulence encountered by SPP.

We develop an analytic model of intermittent, three-dimensional, strong, reduced magnetohydrodynamic turbulence with zero cross helicity. We take the fluctuation amplitudes to have a log-Poisson distribution and incorporate into the model a new phenomenology of scale-dependent dynamic alignment. The log-Poisson distribution in our model is characterized by two parameters. To calculate these parameters, we make use of three assumptions: that the energy cascade rate is independent of scale within the inertial range, that the most intense coherent structures at scale $\lambda$ are sheet-like with a volume filling factor proportional to $\lambda$, and that most of the cascade power arises from interactions between exceptionally intense fluctuations and much weaker fluctuations. We then compute the scalings of the power spectrum, the kurtosis, higher-order structure functions, and three different average alignment angles. These scalings appear to be consistent with numerous results from direct numerical simulations.

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 investigate the anisotropy of Alfv\'enic turbulence in the inertial range of slow solar wind and in both driven and decaying reduced magnetohydrodynamic simulations. A direct comparison is made by measuring the anisotropic second-order structure functions in both data sets. In the solar wind, the perpendicular spectral index of the magnetic field is close to -5/3. In the forced simulation, it is close to -5/3 for the velocity and -3/2 for the magnetic field. In the decaying simulation, it is -5/3 for both fields. The spectral index becomes steeper at small angles to the local magnetic field direction in all cases. We also show that when using the global rather than local mean field, the anisotropic scaling of the simulations cannot always be properly measured.