Search

Measuring the energy cascade rate in space plasmas is a challenging task for several reasons. This quantity is (i) inherently three-dimensional (ii) scale-dependent, (iii) anisotropic in the interplanetary plasma, and (iv) requires measurements of plasma parameters in at least four points. Here, we show how three of such problems have been addressed by applying the novel Lag Polyhedra Derivative Ensemble (LPDE) technique to the Magnetospheric Multiscale (MMS) mission in the Earth's magnetosheath.

It is common for researchers to record long, multiple time series from experiments or calculations. But sometimes there are no good models for the systems or no applicable mathematical theorems that can tell us when there are basic relationships between subsets of the time series data such as continuity, differentiability, embeddings, etc. The data is often higher dimensional and simple plotting will not guide us. At that point fitting the data to polynomials, Fourier series, etc. becomes uncertain. Even at the simplest level, having data that shows there is a function between the data subsets is useful and a negative answer means that more particular data fitting or analysis will be suspect and probably fail. We show here statistics that test time series subsets for basic mathematical properties and relations between them that not only indicate when more specific analyses are safe to do, but whether the systems are operating correctly. We apply these statistics to examples from reservoir computing where an important property of reservoir computers is that the reservoir system establishes an embedding of the drive system in order to make any other calculations with the reservoir computer successful., Comment: 15 pages, 10 figures

Using Parker Solar Probe data from orbits 8 through 17, we examine fluctuation amplitudes throughout the critical region where the solar wind flow speed approaches and then exceeds the Alfv\'en wave speed, taking account of various exigencies of the plasma data. In contrast to WKB theory for non-interacting Alfv\'en waves streaming away from the Sun, the magnetic and kinetic fluctuation energies per unit volume are not monotonically decreasing. Instead, there is clear violation of conservation of standard WKB wave action, which is consistent with previous indications of strong in-situ fluctuation energy input in the solar wind near the Alfv\'en critical region. This points to strong violations of WKB theory due to nonlinearity (turbulence) and major energy input near the critical region, which we interpret as likely due to driving by large-scale coronal shear flows.

We present a broad review of 1/f noise observations in the heliosphere, and discuss and complement the theoretical background of generic 1/f models as relevant to NASA's Polarimeter to Unify the Corona and Heliosphere (PUNCH) mission. First observed in the voltage fluctuations of vacuum tubes, the scale-invariant 1/f spectrum has since been identified across a wide array of natural and artificial systems, including heart rate fluctuations and loudness patterns in musical compositions. In the solar wind, the interplanetary magnetic field trace spectrum exhibits 1/f scaling within the frequency range from around 2e-6 Hz to 1e-4 Hz at 1 au. One compelling mechanism for the generation of 1/f noise is the superposition principle, where a composite 1/f spectrum arises from the superposition of a collection of individual power-law spectra characterized by a scale-invariant distribution of correlation times. In the context of the solar wind, such a superposition could originate from scale-invariant reconnection processes in the corona. Further observations have detected 1/f signatures in the photosphere and corona at frequency ranges compatible with those observed at 1 au, suggesting an even lower altitude origin of 1/f spectrum in the solar dynamo itself. This hypothesis is bolstered by dynamo experiments and simulations that indicate inverse cascade activities, which can be linked to successive flux tube reconnections beneath the corona, and are known to generate 1/f noise possibly through nonlocal interactions at the largest scales. Conversely, models positing in situ generation of 1/f signals face causality issues in explaining the low-frequency portion of the 1/f spectrum. Understanding 1/f noise in the solar wind may inform central problems in heliospheric physics, such as the solar dynamo, coronal heating, the origin of the solar wind, and the nature of interplanetary turbulence., Comment: 19 pages, 10 figures, submitted to Solar Physics

Plasma turbulence involves complex, nonlinear interactions of electromagnetic fields and charged particles across multiple scales. Studying these phenomena in space plasmas, like the solar wind, is facilitated by the intrinsic scale separations and the availability of in situ spacecraft observations. However, the single-point or single-scale configurations of current spacecraft limit our understanding of many properties of the turbulent solar wind. To overcome these limitations, multipoint measurements spanning a range of characteristic scales are essential. This paper prepares for the enhanced measurement capabilities of upcoming multispacecraft missions by demonstrating that higher-order statistics, specifically kurtosis, as a baseline for intermittency can be accurately measured. Using synthetic turbulent fields with adjustable intermittency levels, we achieve scale separations analogous to those in the solar wind and apply these techniques to the planned trajectories of the HelioSwarm mission. This approach promises significant advancements in our understanding of plasma turbulence.

We present an event observed by Parker Solar Probe at $\sim$0.2 au on March 2, 2022 in which imaging and \emph{in situ} measurements coincide. During this event, PSP passed through structures on the flank of a streamer blowout CME including an isolated flux tube in front of the CME, a turbulent sheath, and the CME itself. Imaging observations and \emph{in situ} helicity and principal variance signatures consistently show the presence of flux ropes internal to the CME. In both the sheath, and the CME interval, the distributions are more isotropic, the spectra are softer, and the abundance ratios of Fe/O and He/H are lower than those in the isolated flux tube, and yet elevated relative to typical plasma and SEP abundances. These signatures in the sheath and the CME indicate that both flare populations and those from the plasma are accelerated to form the observed energetic particle enhancements. In contrast, the isolated flux tube shows large streaming, hard spectra and large Fe/O and He/H ratios, indicating flare sources. Energetic particle fluxes are most enhanced within the CME interval from suprathermal through energetic particle energies ($\sim$ keV to $>10$ MeV), indicating particle acceleration, and confinement local to the closed magnetic structure. The flux-rope morphology of the CME helps to enable local modulation and trapping of energetic particles, particularly along helicity channels and other plasma boundaries. Thus, the CME acts to build-up energetic particle populations, allowing them to be fed into subsequent higher energy particle acceleration throughout the inner heliosphere where a compression or shock forms on the CME front., Comment: 41 pages, 19 figures, In Press

The transition from subAlfv\'enic to superAlfv\'enic flow in the solar atmosphere is examined by means of Parker Solar Probe (PSP) measurements during solar encounters 8 to 14. Around 220 subAlfv\'enic periods with a duration $\ge$ 10 minutes are identified. The distribution of their durations, heliocentric distances, and Alfv\'en Mach number are analyzed and compared with a global magnetohydrodynamic model of the solar corona and wind, which includes turbulence effects. The results are consistent with a patchy and fragmented morphology, and suggestive of a turbulent Alfv\'en zone within which the transition from subAlfv\'enic to superAlfv\'enic flow occurs over an extended range of helioradii. These results inform and establish context for detailed analyses of subAlfv\'enic coronal plasma that are expected to emerge from PSP's final mission phase, as well as for NASA's planned PUNCH mission.

A well-known property of solar wind plasma turbulence is the observed anisotropy of the autocorrelations, or equivalently the spectra, of velocity and magnetic field fluctuations. Here we explore the related but apparently not well-studied issue of the anisotropy of plasma density fluctuations in the energy-containing and inertial ranges of solar wind turbulence. Using 10 years (1998-2008) of in situ data from the Advanced Composition Explorer (ACE) mission, we find that for all but the fastest wind category, the density correlation scale is slightly larger in directions quasi-parallel to the large-scale mean magnetic field as compared to quasi-perpendicular directions. The correlation scale in fast wind is consistent with isotropic. The anisotropy as a function of the level of correlation is also explored. We find at small correlation levels, i.e., at energy-containing scales and larger, the density fluctuations are close to isotropy for fast wind, and slightly favor more rapid decorrelation in perpendicular directions for slow and medium winds. At relatively smaller (inertial range) scales where the correlation values are larger, the sense of anisotropy is reversed in all speed ranges, implying a more "slab-like" structure, especially prominent in the fast wind samples. We contrast this finding with published results on velocity and magnetic field correlations., Comment: 14 pages, 6 figures, accepted by the Astrophysical Journal

Solar energetic particle (SEP) events have been observed by the Parker Solar Probe (PSP) spacecraft since its launch in 2018. These events include sources from solar flares and coronal mass ejections (CMEs). Onboard PSP is the IS\(\odot\)IS instrument suite measuring ions over energies from ~ 20 keV/nucleon to 200 MeV/nucleon and electrons from ~ 20 keV to 6 MeV. Previous studies sought to group CME characteristics based on their plasma conditions and arrived at general descriptions with large statistical errors, leaving open questions on how to properly group CMEs based solely on their plasma conditions. To help resolve these open questions, plasma properties of CMEs have been examined in relation to SEPs. Here we reexamine one plasma property, the solar wind proton temperature, and compare it to the proton SEP intensity in a region immediately downstream of a CME-driven shock for seven CMEs observed at radial distances within 1 au. We find a statistically strong correlation between proton SEP intensity and bulk proton temperature, indicating a clear relationship between SEPs and the conditions in the solar wind. Furthermore, we propose that an indirect coupling of SEP intensity to the level of turbulence and the amount of energy dissipation that results is mainly responsible for the observed correlation between SEP intensity and proton temperature. These results are key to understanding the interaction of SEPs with the bulk solar wind in CME-driven shocks and will improve our ability to model the interplay of shock evolution and particle acceleration., Comment: 12 pages, 4 figures, and 2 tables

Small-scale interplanetary magnetic flux ropes (SMFRs) are similar to ICMEs in magnetic structure, but are smaller and do not exhibit ICME plasma signatures. We present a computationally efficient and GPU-powered version of the single-spacecraft automated SMFR detection algorithm based on the Grad-Shafranov (GS) technique. Our algorithm is capable of processing higher resolution data, eliminates selection bias caused by a fixed $\avg{B}$ threshold, has improved detection criteria demonstrated to have better results on an MHD simulation, and recovers full 2.5D cross sections using GS reconstruction. We used it to detect 512,152 SMFRs from 27 years (1996 to 2022) of 3-second cadence \emph{Wind} measurements. Our novel findings are: (1) the radial density of SMFRs at 1 au (${\sim}1$ per $\si{10^6\kilo\meter}$) and filling factor (${\sim}$35\%) are independent of solar activity, distance to the heliospheric current sheet (HCS), and solar wind plasma type, although the minority of SMFRs with diameters greater than ${\sim}$0.01 au have a strong solar activity dependence; (2) SMFR diameters follow a log-normal distribution that peaks below the resolved range ($\gtrsim 10^4$ km), although the filling factor is dominated by SMFRs between $10^5$ to $10^6$ km; (3) most SMFRs at 1 au have strong field-aligned flows like those from PSP measurements; (4) in terms of diameter $d$, SMFR poloidal flux $\propto d^{1.2}$, axial flux $\propto d^{2.0}$, average twist number $\propto d^{-0.8}$, current density $\propto d^{-0.8}$, and helicity $\propto d^{3.2}$. Implications for the origin of SMFRs and switchbacks are briefly discussed. The new algorithm and SMFR dataset are made freely available.

We examine dissipation and energy conversion in weakly collisional plasma turbulence, employing in situ observations from the Magnetospheric Multiscale (MMS) mission and kinetic Particle-in-Cell (PIC) simulations of proton-electron plasma. A previous result indicated the presence of viscous-like and resistive-like scaling of average energy conversion rates -- analogous to scalings characteristic of collisional systems. This allows for extraction of collisional-like coefficients of effective viscosity and resistivity, and thus also determination of effective Reynolds numbers based on these coefficients. The effective Reynolds number, as a measure of the available bandwidth for turbulence to populate various scales, links macro turbulence properties with kinetic plasma properties in a novel way., Comment: 9 pages, 2 figures, 5 tables

Magnetic flux ropes are a key component of coronal mass ejections, forming the core of these eruptive phenomena. However, determining whether a flux rope is present prior to eruption onset and, if so, the rope's handedness and the number of turns that any helical field lines make is difficult without magnetic field modelling or in-situ detection of the flux rope. We present two distinct observations of plasma flows along a filament channel on 4 and 5 September 2022 made using the \textit{Solar Orbiter} spacecraft. Each plasma flow exhibited helical motions in a right-handed sense as the plasma moved from the source active region across the solar disk to the quiet Sun, suggesting that the magnetic configuration of the filament channel contains a flux rope with positive chirality and at least one turn. The length and velocity of the plasma flow increased from the first to the second observation, suggesting evolution of the flux rope, with the flux rope subsequently erupting within $\sim$5~hours of the second plasma flow. The erupting flux rope then passed over the \textit{Parker Solar Probe} spacecraft during its Encounter 13, enabling \textit{in-situ} diagnostics of the structure. Although complex and consistent with the flux rope erupting from underneath the heliospheric current sheet, the \textit{in-situ} measurements support the inference of a right-handed flux rope from remote-sensing observations. These observations provide a unique insight into the eruption and evolution of a magnetic flux rope near the Sun., Comment: 19 pages, 11 figures, accepted for publication in The Astrophysical Journal

A novel multispacecraft technique applied to Magnetospheric Multiscale (MMS) mission data collected in the Earth's magnetosheath enables evaluation of the energy cascade rate solving the full Yaglom's equation in a turbulent space plasma. The method differs from existing approaches in that (i) it is inherently three-dimensional; (ii) it provides a statistically significant number of estimates from a single data stream; and (iii) it allows for a direct visualization of energy flux in turbulent plasmas. This new technique will ultimately provide a realistic, comprehensive picture of the turbulence process in plasmas.

Many questions remain about the compressibility of solar wind turbulence with respect to its origins and properties. Low plasma beta (ratio of thermal to magnetic pressure) environments allow for the easier generation of compressible turbulence, enabling study of the relationship between density fluctuations and turbulent Mach number. Utilizing Parker Solar Probe plasma data, we examine the normalized proton density fluctuations $\langle \delta n_p^2 \rangle ^{1/2}/\langle n_p\rangle = \delta {n_p}_{rms}/\langle n_p\rangle$ as a function of turbulent Mach number $M_t$ conditioned on plasma beta and cross helicity. With consideration of statistical error in the parameters computed from in-situ data, we find a general result that $\delta {n_p}_{rms}/\langle n_p\rangle \sim M_t^{1.18 \pm 0.04}$, consistent with both linear-wave theory, and nearly-incompressible turbulence in an inhomogeneous background field. We compare observational results conditioned on plasma beta and cross helicity with 3D magnetohydrodynamic simulations, and observe rather significant similarities with respect to how those parameters affect the proportionality between density fluctuations and turbulent Mach number. This study further investigates the complexity of compressible turbulence as viewed by the density scaling relationship, and may help better understand the compressible environment of the near-Sun solar wind., Comment: 8 pages, 3 figures, 1 table, submitted to ApJL

In turbulence, nonlinear terms drive energy transfer from large-scale eddies into small scales through the so-called energy cascade. Turbulence often relaxes toward states that minimize energy; typically these states are considered globally. However, turbulence can also relax toward local quasi-equilibrium states, creating patches or cells where the magnitude of nonlinearity is reduced and energy cascade is impaired. We show, for the first time, compelling observational evidence that this ``cellularization'' of turbulence can occur due to local relaxation in a strongly turbulent natural environment such as the Earth's magnetosheath.

Exploration of plasma dynamics in space, including turbulence, is entering a new era of multi-satellite constellation measurements that will determine fundamental properties with unprecedented precision. Familiar but imprecise approximations will need to be abandoned and replaced with more advanced approaches. We present a preparatory study of the evaluation of second- and third-order statistics, using simultaneous measurements at many points. Here, for specificity, the orbital configuration of the NASA Helioswarm mission is employed in conjunction with three-dimensional magnetohydrodynamics numerical simulations of turbulence. The Helioswarm 9-spacecraft constellation flies virtually through the turbulence to compare results with the exact numerical statistics. We demonstrate novel increment-based techniques for the computation of (1) the multidimensional spectra and (2) the turbulent energy flux. This latter increment-space estimate of the cascade rate, based on the third-order Yaglom-Politano-Pouquet theory, uses numerous increment-space tetrahedra. Our investigation reveals that Helioswarm will provide crucial information on the nature of astrophysical turbulence.

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 examine the radial evolution of correlation lengths perpendicular (\(\lambda_C^{\perp}\)) and parallel (\(\lambda_C^{\parallel}\)) to the magnetic-field direction, computed from solar wind magnetic-field data measured by Parker Solar Probe (PSP) during its first eight orbits, Helios 1, Advanced Composition Explorer (ACE), WIND, and Voyager 1 spacecraft. Correlation lengths are grouped by an interval's alignment angle; the angle between the magnetic-field and solar wind velocity vectors (\(\Theta_{\rm BV}\)). Parallel and perpendicular angular channels correspond to angles \(0^{\circ}~<~\Theta_{\rm BV}~<~40^{\circ}\) and \(50^{\circ}~<~\Theta_{\rm BV}~<~90^{\circ}\), respectively. We observe an anisotropy in the inner heliosphere within 0.40~au, with \(\lambda_C^{\parallel} / \lambda_C^{\perp} \approx 0.75\) at 0.10~au. This anisotropy reduces with increasing heliocentric distance and the correlation lengths roughly isotropize within 1~au. Results from ACE and WIND support a reversal of the anisotropy, such that \(\lambda_C^{\parallel} /\lambda_C^{\perp} \approx 1.29\) at 1~au. The ratio does not appear to change significantly beyond 1~au, although the small number of parallel intervals in the Voyager dataset precludes unambiguous conclusions from being drawn. This study provides insights regarding the radial evolution of the large, most energetic interacting turbulent fluctuations in the heliosphere. We also emphasize the importance of tracking the changes in sampling direction in PSP measurements as the spacecraft approaches the Sun, when using these data to study the radial evolution of turbulence. This can prove to be vital in understanding the more complex dynamics of the solar wind in the inner heliosphere and can assist in improving related simulations., Comment: 9 pages, 3 figures, 2 tables

Collisionless shocks and plasma turbulence are crucial ingredients for a broad range of astrophysical systems. The shock-turbulence interaction, and in particular the transmission of fully developed turbulence across the quasi-perpendicular Earth's bow shock, is here addressed using a combination of spacecraft observations and local numerical simulations. An alignment between the Wind (upstream) and MMS (downstream) spacecraft is used to study the transmission of turbulent structures across the shock, revealing an increase of their magnetic helicity content in its downstream. Local kinetic simulations, in which the dynamics of turbulent structures is followed through their transmission across a perpendicular shock, confirm this scenario, revealing that the observed magnetic helicity increase is associated with the compression of turbulent structures at the shock front., Comment: 9 Pages, 4 Figures

The subject of switchbacks, defined either as large angular deflections or polarity reversals of the magnetic field, has generated substantial interest in the space physics community since the launch of Parker Solar Probe (PSP) in 2018. Previous studies have characterized switchbacks in several different ways, and have been restricted to data available from the first few orbits. Here, we analyze the frequency of occurrence of switchbacks per unit distance for the first full eight orbits of PSP. In this work, are considered switchback only the events that reverse the sign of magnetic field relative to a regional average. A significant finding is that the rate of occurrence falls off sharply approaching the sun near 0.2 au (40 $R_\odot$), and rises gently from 0.2 au outward. The analysis is varied for different magnetic field cadences and for different local averages of the ambient field, confirming the robustness of the results. We discuss implications for the mechanisms of switchback generation. A publicly available database has been created with the identified reversals.