A discussion of the advantages and limitations of the concept of critical balance, as employed in turbulence phenomenologies, is presented. The incompressible magnetohydrodynamic (MHD) case is a particular focus. The discussion emphasizes the status of the original Goldreich & Sridhar (1995) critical balance conjecture relative to related theoretical issues and models in an MHD description of plasma turbulence. Issues examined include variance and spectral anisotropy, influence of a mean magnetic field, local and nonlocal effects, and the potential for effects of external driving. Related models such as Reduced MHD provide a valuable context in the considerations. Some new results concerning spectral features and timescales are presented in the course of the discussion. Also mentioned briefly are some adaptations and variations of critical balance.

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.

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.

We study the evolution of kinetic and magnetic energy spectra in magnetohydrodynamic flows in the presence of strong cross helicity. For forced turbulence, we find weak inverse transfer of kinetic energy toward the smallest wavenumber. This is plausibly explained by the finiteness of scale separation between the injection wavenumber and the smallest wavenumber of the domain, which here is a factor of 15. In the decaying case, there is a slight increase at the smallest wavenumber, which is probably explained by the dominance of kinetic energy over magnetic energy at the smallest wavenumbers. Within a range of wavenumbers covering almost an order of magnitude the decay is purely exponential, which is argued to be a consequence of a suppression of nonlinearity due to the presence of strong cross helicity.

We extend a two-component model for the evolution of fluctuations in the solar wind plasma so that it is fully three-dimensional (3D) and also coupled self-consistently to the large-scale magnetohydrodynamic (MHD) equations describing the background solar wind. The two classes of fluctuations considered are a high-frequency parallel-propagating wave-like piece and a low-frequency quasi-two-dimensional component. For both components, the nonlinear dynamics is dominanted by quasi-perpendicular spectral cascades of energy. Driving of the fluctuations, by, for example, velocity shear and pickup ions, is included. Numerical solutions to the new model are obtained using the Cronos framework, and validated against previous simpler models. Comparing results from the new model with spacecraft measurements, we find improved agreement relative to earlier models that employ prescribed background solar wind fields. Finally, the new results for the wave-like and quasi-two-dimensional fluctuations are used to calculate ab initio diffusion mean free paths and drift lengthscales for the transport of cosmic rays in the turbulent solar wind.

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus establish that the local structure function is not a measure of the energy spectrum. Evidently the local enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon must be understood in terms of higher order correlations, fourth-order and above.

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data.

Previous formulations of heating and transport associated with strong magnetohydrodynamic (MHD) turbulence are generalized to incorporate separate internal energy equations for electrons and protons. Electron heat conduction is included. Energy is supplied by turbulent heating that affects both electrons and protons, and is exchanged between them via collisions. Comparison to available Ulysses data shows that a reasonable accounting for the data is provided when (i) the energy exchange timescale is very long and (ii) the deposition of heat due to turbulence is divided, with 60% going to proton heating and 40% into electron heating. Heat conduction, determined here by an empirical fit, plays a major role in describing the electron data.

The heating of the lower solar corona is examined using numerical simulations and theoretical models of magnetohydrodynamic turbulence in open magnetic regions. A turbulent energy cascade to small length scales perpendicular to the mean magnetic field can be sustained by driving with low-frequency Alfven waves reflected from mean density and magnetic field gradients. This mechanism deposits energy efficiently in the lower corona, and we show that the spatial distribution of the heating is determined by the mean density through the Alfven speed profile. This provides a robust heating mechanism that can explain observed high coronal temperatures and accounts for the significant heating (per unit volume) distribution below two solar radius needed in models of the origin of the solar wind. The obtained heating per unit mass on the other hand is much more extended indicating that the heating on a per particle basis persists throughout all the lower coronal region considered here.

We determine that the sizes of bursts in mean-square current density in a reduced magnetohydrodynamic (RMHD)simulation follow power-law probability density function (PDF). The PDFs for burst durations and waiting time between bursts are clearly not exponential and could also be power-law. This suffices to distinguish their behaviour from the original Bak et al. sandpile model which had exponential waiting time PDFs. However, it is not sufficient to distinguish between turbulence, some other SOC-like models, and other red noise sources.