Search

Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere., Comment: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

We study the nature of pressure-strain interaction at reconnection sites, detected by NASA's Magnetospheric Multiscale (MMS) Mission. We employ data from a series of published case studies, including a large-scale reconnection event at the magnetopause, three small-scale reconnection events at the magnetosheath current sheets, and one example of the recently discovered electron-only reconnection. In all instances, we find that the pressure-strain shows signature of conversion into (or from) internal energy at the reconnection site. The electron heating rate is larger than the ion heating rate and the compressive heating is dominant over the incompressive heating rate in all cases considered. The magnitude of thermal energy conversion rate is close to the electromagnetic energy conversion rate in the reconnection region. Although in most cases the pressure-strain interaction indicates that the particle internal energy is increasing, in one case the internal energy is decreasing. These observations indicate that the pressure-strain interaction can be used as an independent measure of energy conversion and dynamics in reconnection regions, in particular independent of measures based on the electromagnetic work. Finally, we explore a selected reconnection site in a turbulent Particle-in-Cell (PIC) simulation which further supports the observational results., Comment: The following article has been accepted by Physics of Plasmas

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von K\'arm\'an Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady cascade rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. These results provide explicit evidence that reconnection itself is fundamentally an energy cascade process.

Motivated by prior remote observations of a transition from striated solar coronal structures to more isotropic ``flocculated'' fluctuations, we propose that the dynamics of the inner solar wind just outside the Alfv\'en critical zone, and in the vicinity of the first $\beta=1$ surface, is powered by the relative velocities of adjacent coronal magnetic flux tubes. We suggest that large amplitude flow contrasts are magnetically constrained at lower altitude but shear-driven dynamics are triggered as such constraints are released above the Alfv\'en critical zone, as suggested by global magnetohydrodynamic (MHD) simulations that include self-consistent turbulence transport. We argue that this dynamical evolution accounts for features observed by {\it Parker Solar Probe} ({\it PSP}) near initial perihelia, including magnetic ``switchbacks'', and large transverse velocities that are partially corotational and saturate near the local Alfv\'en speed. Large-scale magnetic increments are more longitudinal than latitudinal, a state unlikely to originate in or below the lower corona. We attribute this to preferentially longitudinal velocity shear from varying degrees of corotation. Supporting evidence includes comparison with a high Mach number three-dimensional compressible MHD simulation of nonlinear shear-driven turbulence, reproducing several observed diagnostics, including characteristic distributions of fluctuations that are qualitatively similar to {\it PSP} observations near the first perihelion. The concurrence of evidence from remote sensing observations, {\it in situ} measurements, and both global and local simulations supports the idea that the dynamics just above the Alfv\'en critical zone boost low-frequency plasma turbulence to the level routinely observed throughout the explored solar system.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins., Comment: submitted to ApJS, welcome comments

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., Comment: Submitted to special issue of ApJ for Parker Solar Probe

The solar wind proton temperature at 1-au has been found to be correlated with small-scale intermittent magnetic structures, i.e., regions with enhanced temperature are associated with coherent structures such as current sheets. Using Parker Solar Probe data from the first encounter, we study this association using measurements of radial proton temperature, employing the Partial Variance of Increments (PVI) technique to identify intermittent magnetic structures. We observe that the probability density functions of high-PVI events have higher median temperatures than those with lower PVI, The regions in space where PVI peaks were also locations that had enhanced temperatures when compared with similar regions suggesting a heating mechanism in the young solar wind that is associated with intermittency developed by a nonlinear turbulent cascade.n the immediate vicinity., Comment: 6 pages, 3 figures, part of ApJ special issue for PSP

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., Comment: Accepted for publication in The Astrophysical Journal Supplement, PSP special issue

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., Comment: Accepted for publication in The Astrophysical Journal Supplement, PSP special issue

The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell (PIC) simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power-law behavior, with a slope near -5/3 at wavenumbers smaller than the inverse ion inertial length, and a slope steeper than -8/3 for larger wavenumbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70 degrees to about 48 degrees (nearly isotropic) by the end of the simulation., Comment: 8 figures

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., Comment: White paper submitted to the PLASMA 2020 Decadal Survey Committee

Using kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion Larmor radius, the ions do not respond to the reconnection dynamics leading to ''electron-only'' reconnection with very large quasi-steady reconnection rates. The transition to more traditional ''ion-coupled'' reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfv\'en speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

In presence of an externally supported, mean magnetic field a turbulent, conducting medium, such as plasma, becomes anisotropic. This mean magnetic field, which is separate from the fluctuating, turbulent part of the magnetic field, has considerable effects on the dynamics of the system. In this paper, we examine the dissipation rates for decaying incompressible magnetohydrodynamic (MHD) turbulence with increasing Reynolds number, and in the presence of a mean magnetic field of varying strength. Proceeding numerically, we find that as the Reynolds number increases, the dissipation rate asymptotes to a finite value for each magnetic field strength, confirming the K\'arm\'an-Howarth hypothesis as applied to MHD. The asymptotic value of the dimensionless dissipation rate is initially suppressed from the zero-mean-field value by the mean magnetic field but then approaches a constant value for higher values of the mean field strength. Additionally, for comparison, we perform a set of two-dimensional (2DMHD) and a set of reduced MHD (RMHD) simulations. We find that the RMHD results lie very close to the values corresponding to the high mean-field limit of the three-dimensional runs while the 2DMHD results admit distinct values far from both the zero mean field cases and the high mean field limit of the three-dimensional cases. These findings provide firm underpinnings for numerous applications in space and astrophysics wherein von K\'arm\'an decay of turbulence is assumed., Comment: Accepted for publication in Physical Review X

Studies of solar wind turbulence traditionally employ high-resolution magnetic field data, but high-resolution measurements of ion and electron moments have been possible only recently. We report the first turbulence studies of ion and electron velocity moments accumulated in pristine solar wind by the Fast Particle Investigation instrument onboard the Magnetospheric Multiscale (MMS) Mission. Use of these data is made possible by a novel implementation of a frequency domain Hampel filter, described herein. After presenting procedures for processing of the data, we discuss statistical properties of solar wind turbulence extending into the kinetic range. Magnetic field fluctuations dominate electron and ion velocity fluctuation spectra throughout the energy-containing and inertial ranges. However, a multi-spacecraft analysis indicates that at scales shorter than the ion-inertial length, electron velocity fluctuations become larger than ion velocity and magnetic field fluctuations. The kurtosis of ion velocity peaks around few ion-inertial lengths and returns to near gaussian value at sub-ion scales., Comment: Accepted for publication in The Astrophysical Journal Supplement

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., Comment: Submitted to The Astrophysical Journal Letters

Using observational data from the \emph{Magnetospheric Multiscale} (MMS) Mission in the Earth's magnetosheath, we estimate the energy cascade rate using different techniques within the framework of incompressible magnetohydrodynamic (MHD) turbulence. At the energy containing scale, the energy budget is controlled by the von K\'arm\'an decay law. Inertial range cascade is estimated by fitting a linear scaling to the mixed third-order structure function. Finally, we use a multi-spacecraft technique to estimate the Kolmogorov-Yaglom-like cascade rate in the kinetic range, well below the ion inertial length scale. We find that the inertial range cascade rate is almost equal to the one predicted by the von K\'arm\'an law at the energy containing scale, while the cascade rate evaluated at the kinetic scale is somewhat lower, as anticipated in theory~\citep{Yang2017PoP}. Further, in agreement with a recent study~\citep{Hadid2018PRL}, we find that the incompressive cascade rate in the Earth's magnetosheath is about $1000$ times larger than the cascade rate in the pristine solar wind., Comment: Accepted for publication in Astrophysical Journal

Plasma turbulence is investigated using high-resolution ion velocity distributions measured by the Magnetospheric Multiscale Mission (MMS) in the Earth's magnetosheath. The particle distribution is highly structured, suggesting a cascade-like process in velocity space. This complex velocity space structure is investigated using a three-dimensional Hermite transform that reveals a power law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space cascade. The scaling theory is in agreement with observations, suggesting a new path for the study of plasma turbulence in weakly collisional space and astrophysical plasmas.

Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients., Comment: 28 pages, 10 figures

The analysis of the Parker-Moffatt problem, recently revisited in Pezzi et al. (2016), is here extended by including the Hall magnetohydrodynamics and two hybrid kinetic Vlasov-Maxwell numerical models. The presence of dispersive and kinetic features is studied in detail and a comparison between the two kinetic codes is also reported. Focus on the presence of non-Maxwellian signatures shows that - during the collision - regions characterized by strong temperature anisotropy are recovered and the proton distribution function displays a beam along the direction of the magnetic field, similar to some recent observations of the solar wind., Comment: Accepted for publication in J. Plasma Phys

The interaction of two colliding Alfv\'en wave packets is here described by means of magnetohydrodynamics (MHD) and hybrid kinetic numerical simulations. The MHD evolution revisits the theoretical insights described by Moffatt, Parker, Kraichnan, Chandrasekhar and Els\"asser in which the oppositely propagating large amplitude wave packets interact for a finite time, initiating turbulence. However, the extension to include compressive and kinetic effects, while maintaining the gross characteristics of the simpler classic formulation, also reveals intriguing features which go beyond the pure MHD treatment.

The kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations. In the magnetohydrodynamic regime this vortex leads rapidly to broadband turbulence. Significant differences from MHD arise at small scales, where the fluid scale energy dissipates into heat almost exclusively through the magnetic field because the protons are decoupled from the magnetic field. Although cyclotron resonance is absent, the protons heat preferentially in the plane perpendicular to the mean field, as in the corona and solar wind. Effective transport coefficients are calculated., Comment: 4 pages, 4 figures. Submitted to PRL

We present a statistical study of coherent structures at kinetic scales, using data from the Magnetospheric Multiscale mission in the Earths magnetosheath. We implemented the multi-spacecraft partial variance of increments (PVI) technique to detect these structures, which are associated with intermittency at kinetic scales. We examine the properties of the electron heating occurring within such structures. We find that, statistically, structures with a high PVI index are regions of significant electron heating. We also focus on one such structure, a current sheet, which shows some signatures consistent with magnetic reconnection. Strong parallel electron heating coincides with whistler emissions at the edges of the current sheet.