When magnetohydrodynamic turbulence evolves in the presence of a large-scale mean magnetic field, an anisotropy develops relative to that preferred direction. The well-known tendency is to develop stronger gradients perpendicular to the magnetic field, relative to the direction along the field. This anisotropy of the spectrum is deeply connected with anisotropy of estimated timescales for dynamical processes, and requires reconsideration of basic issues such as scale locality and spectral transfer. Here analysis of high-resolution three-dimensional simulations of unforced magnetohydrodynamic turbulence permits quantitative assessment of the behavior of theoretically relevant timescales in Fourier wavevector space. We discuss the distribution of nonlinear times, Alfv\'en times, and estimated spectral transfer rates. Attention is called to the potential significance of special regions of the spectrum, such as the two-dimensional limit and the "critical balance" region. A formulation of estimated spectral transfer in terms of a suppression factor supports a conclusion that the quasi two-dimensional fluctuations (characterized by strong nonlinearities) are not a singular limit, but may be in general expected to make important contributions.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity.

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

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.

The Parker Solar Probe (PSP) primary mission extends seven years and consists of 24 orbits of the Sun with descending perihelia culminating in a closest approach of ($\sim 9.8~R_\odot$). In the course of these orbits PSP will pass through widely varying conditions, including anticipated large variations of turbulence properties such as energy density, correlation scales and cross helicities. Here we employ global magnetohydrodynamics simulations with self-consistent turbulence transport and heating \citep{usmanov2018} to preview likely conditions that will be encountered by PSP, by assuming suitable boundary conditions at the coronal base. The code evolves large-scale parameters -- such as velocity, magnetic field, and temperature -- as well as turbulent energy density, cross helicity, and correlation scale. These computed quantities provide the basis for evaluating additional useful parameters that are derivable from the primary model outputs. Here we illustrate one such possibility in which computed turbulence and large-scale parameters are used to evaluate the accuracy of the Taylor "frozen-in" hypothesis along the PSP trajectory. Apart from the immediate purpose of anticipating turbulence conditions that PSP will encounter, as experience is gained in comparisons of observations with simulated data, this approach will be increasingly useful for planning and interpretation of subsequent observations.

The solar corona and young solar wind may be characterized by critical surfaces -- the sonic, Alfv\'en, and first plasma-$\beta$ unity surfaces -- that demarcate regions where the solar wind flow undergoes certain crucial transformations. Global numerical simulations and remote sensing observations offer a natural mode for the study of these surfaces at large scales, thus providing valuable context for the high-resolution in-situ measurements expected from the soon-to-be-launched Parker Solar Probe (PSP). The present study utilizes global three-dimensional magnetohydrodynamic simulations of the solar wind to characterize the critical surfaces and investigate the flow in propinquitous regions. Effects of solar activity are incorporated by varying source magnetic dipole tilts and employing magnetogram-based boundary conditions. A magnetohydrodynamic turbulence model is self-consistently coupled to the bulk flow equations, enabling investigation of turbulence properties of the flow in the vicinity of critical regions. The simulation results are compared with a variety of remote sensing observations. A simulated PSP trajectory is used to provide contextual predictions for the spacecraft in terms of the computed critical surfaces. Broad agreement is seen in the interpretation of the present results in comparison with existing remote sensing results, both from heliospheric imaging and from radio scintillation studies. The trajectory analyses show that the period of time that PSP is likely to spend inside the $\beta=1$, sonic and Alfv\'en surfaces depends sensitively on the degree of solar activity and the tilt of the solar dipole and location of the heliospheric current sheet.

We study spectral features of ion velocity and magnetic field correlations in the solar wind and in the magnetosheath using data from the Magnetospheric Multi-Scale (MMS) spacecraft. High resolution MMS observations enable the study of transition of these correlations between their magnetofluid character at larger scales into the sub-proton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment and energy partitioning is examined over a suit- able range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.

Protons (ionized hydrogen) in the solar wind frequently exhibit distinct temperatures ($T_{\perp p}$ and $T_{\parallel p}$) perpendicular and parallel to the plasma's background magnetic-field. Numerous prior studies of the interplanetary solar-wind have shown that, as plasma beta ($\beta_{\parallel p}$) increases, a narrower range of temperature-anisotropy ($R_p\equiv T_{\perp p}\,/\,T_{\parallel p}$) values is observed. Conventionally, this effect has been ascribed to the actions of kinetic microinstabilities. This study is the first to use data from the Magnetospheric Multiscale Mission (MMS) to explore such $\beta_{\parallel p}$-dependent limits on $R_p$ in Earth's magnetosheath. The distribution of these data across the $(\beta_{\parallel p},R_p)$-plane reveals limits on both $R_p>1$ and $R_p<1$. Linear Vlasov theory is used to compute contours of constant growth-rate for the ion-cyclotron, mirror, parallel-firehose, and oblique-firehose instabilities. These instability thresholds closely align with the contours of the data distribution, which suggests a strong association of instabilities with extremes of ion temperature anisotropy in the magnetosheath. The potential for instabilities to regulate temperature anisotropy is discussed.

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.