Search for: All records

Abstract It is well known that the nonlinear evolution of magnetohydrodynamic (MHD) turbulence generates current sheets. In the solar wind turbulence, current sheets are frequently observed, and they are believed to be an important pathway for the turbulence energy to dissipate and heat the plasma. In this study, we perform a comprehensive analysis of current sheets in a high-resolution two-dimensional simulation of balanced, incompressible MHD turbulence. The simulation parameters are selected such that tearing mode instability is triggered and plasmoids are generated throughout the simulation domain. We develop an automated method to identify current sheets and accurately quantify their key parameters, including thickness (a), length (L), and Lundquist number (S). Before the triggering of tearing instability, the current sheet lengths are mostly comparable to the energy injection scale. After the tearing mode onsets, smaller current sheets with lower Lundquist numbers are generated. While power-law scaling relations betweenLandaand betweena/LandSare observed, no clear correlation is found between the upstream magnetic field strength and thicknessa. Finally, although the turbulence energy shows anisotropy between the directions parallel and perpendicular to the local magnetic field increment, we do not observe a direct correspondence between the shape of the current sheets and that of the turbulence “eddies.” These results suggest that one needs to be cautious when applying the scale-dependent dynamic alignment model to the analysis of current sheets in MHD turbulence. 

Abstract We perform a comprehensive superposed epoch analysis of more than 200 corotating interaction regions (CIRs) using WIND spacecraft observations at 1 au. The stream interfaces are identified by minimum variance analysis, and turbulence properties are evaluated using wavelet transforms over a wide range of temporal scales. The analysis of normalized cross helicity (σ_c) and normalized residual energy (σ_r) reveals distinct turbulence behaviors across frequencies. The spectral indices of both magnetic and velocity fluctuations smoothly transition from steeper in the slow wind to shallower in the fast wind, while a localized steepening of the velocity spectra near the interface indicates enhanced dissipation due to compression. Across broad frequency bands,σ_cshows a clear dip at the stream interface—signifying increased inward Alfvén wave energy—whereasσ_rdisplays a peak–valley–peak structure mainly driven by large-scale velocity shear. In lower-frequency ranges, velocity shear artificially enhances velocity fluctuation energy, producing strong peaks inσ_r, while higher-frequency ranges show a smooth increase ofσ_rfrom slow wind to fast wind. Nearly half of the analyzed CIRs are accompanied by a heliospheric current sheet (HCS), with many HCSs closely aligned with the stream interface, suggesting an intrinsic link between the two structures. Our findings offer valuable clues for reconciling discrepancies among earlier observational and simulation studies, and provide new insight into how compression and velocity shear modulate solar wind turbulence near CIRs. 

We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicityσ_cand compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, withS₂(b) − S₂(u) ∝ l_⊥or equivalently $-E_r\propto k_{\perp }^{-2}$. HereS₂(b) andS₂(u) are second-order structure functions of magnetic field and velocity respectively. In most runs,S₂(b) develops a scaling relation $S_2\left(\mathit{b}\right)\propto l_{\perp }^{1/2}$( $E_b\propto k_{\perp }^{-3/2}$). In contrast,S₂(u) is consistently shallower thanS₂(b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initialσ_cand are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales. 

Numerical simulations have been an increasingly important tool in space physics. Here, we introduce an open-source three-dimensional compressible Hall-Magnetohydrodynamic (MHD) simulation codeLAPS(UCLA-Pseudo-Spectral,https://github.com/chenshihelio/LAPS). The code adopts a pseudo-spectral method based on Fourier Transform to evaluate spatial derivatives, and third-order explicit Runge-Kutta method for time advancement. It is parallelized using Message-Passing-Interface (MPI) with a “pencil” parallelization strategy and has very high scalability. The Expanding-Box-Model is implemented to incorporate spherical expansion effects of the solar wind. We carry out test simulations based on four classic (Hall)-MHD processes, namely, 1) incompressible Hall-MHD waves, 2) incompressible tearing mode instability, 3) Orszag-Tang vortex, and 4) parametric decay instability. The test results agree perfectly with theory predictions and results of previous studies. Given all its features,LAPSis a powerful tool for large-scale simulations of solar wind turbulence as well as other MHD and Hall-MHD processes happening in space. 

Abstract Switchbacks, defined as Alfvénic reversals in magnetic field polarity, can dissipate their magnetic energy with heliocentric distance. To further investigate this, two distinct solar wind parcels tracing back to a similar solar source region were examined during a radial alignment between Parker Solar Probe (@25.8RS) and Solar Orbiter (@152RS). The one caveat was that the two probes were located on opposite sides of the heliospheric current sheet during the alignment. The two parcels contained a multitude of switchbacks—the parcel closer to the Sun was characterized as a switchback patch (SBP), where background proton velocity (v_p) is comparable to the pristine solar wind (v_sw), while the parcel farther from the Sun showed characteristics attributable to a microstream (MS;v_p > v_sw). It was found that (1) MS contains, on average, 30% fewer switchbacks than SBP, and (2) dynamic and thermal pressures decreased by up to 20% across switchback boundaries in SBP and relatively unchanged in MS. Magnetic relaxation can explain the lower number of switchbacks in MS compared to SBP. Switchback relaxation inside SBP can, in turn, accelerate plasma inside SBP over time and heliocentric distance, thus resulting inv_p>v_swin MS. Therefore, it is hypothesized that magnetic relaxation of switchbacks may cause SBPs to evolve into MSs over time and heliocentric distance. 

Abstract Magnetic switchbacks are rapid high-amplitude reversals of the radial magnetic field in the solar wind that do not involve a heliospheric current sheet crossing. First seen sporadically in the 1970s in Mariner and Helios data, switchbacks were later observed by the Ulysses spacecraft beyond 1 au and have been recently discovered to be a typical component of solar wind fluctuations in the inner heliosphere by the Parker Solar Probe spacecraft. While switchbacks are now well understood to be spherically polarized Alfvén waves thanks to Parker Solar Probe observations, their formation has been an intriguing and unsolved puzzle. Here we provide a simple yet predictive theory for the formation of these magnetic reversals: the switchbacks are produced by the distortion and twisting of circularly polarized Alfvén waves by a transversely varying radial wave propagation velocity. We provide an analytic expression for the magnetic field variation, establish the necessary and sufficient conditions for the formation of switchbacks, and show that the proposed mechanism works in a realistic solar wind scenario. We also show that the theoretical predictions are in excellent agreement with observations, and the high-amplitude radial oscillations are strongly correlated with the shear of the wave propagation speed. The correlation coefficient is around 0.3–0.5 for both encounter 1 and encounter 12. The probability of this being a lucky coincidence is essentially zero with p-values below 0.1%. 

Abstract The heliosphere is permeated with highly structured solar wind originating from the Sun. One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between in situ measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed in situ, and the PSP plasma measurements suffer from a limited field of view. Therefore, it lacks a universal tool to self-contextualize the in situ measurements. Here we develop a novel time series visualization method named Gaussianity Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify in situ structures that are possible remnants of solar atmospheric and magnetic structures spanning more than 7 orders of magnitude, from years to seconds, including polar and midlatitude coronal holes, as well as structures compatible with supergranulation, “jetlets” and “picoflares.” Furthermore, computer simulations of Alfvénic turbulence successfully reproduce the Gaussianization of magnetic magnitude, supporting the observed distribution. Building upon these discoveries, the Gaussianity Scalogram can help future studies to reveal the fractal-like fine structures in the solar wind time series from both PSP and a decades-old data archive. 

Abstract Parker Solar Probe observations reveal that the near-Sun space is almost filled with magnetic switchbacks (“switchbacks” hereinafter), which may be a major contributor to the heating and acceleration of solar wind. Here, for the first time, we develop an analytic model of an axisymmetric switchback with uniform magnetic field strength. In this model, three parameters control the geometry of the switchback: height (length along the background magnetic field), width (thickness along radial direction perpendicular to the background field), and the radial distance from the center of switchback to the central axis, which is a proxy of the size of the switchback along the third dimension. We carry out 3D magnetohydrodynamic simulations to investigate the dynamic evolution of the switchback. Comparing simulations conducted with compressible and incompressible codes, we verify that compressibility, i.e., parametric decay instability, is necessary for destabilizing the switchback. Our simulations also reveal that the geometry of the switchback significantly affects how fast the switchback destabilizes. The most stable switchbacks are 2D-like (planar) structures with large aspect ratios (length to width), consistent with the observations. We show that when plasma beta (β) is smaller than one, the switchback is more stable asβincreases. However, whenβis greater than 1, the switchback becomes very unstable as the pattern of the growing compressive fluctuations changes. Our results may explain some of the observational features of switchbacks, including the large aspect ratios and nearly constant occurrence rates in the inner heliosphere. 

Abstract Alfvén waves, considered one of the primary candidates for heating and accelerating the fast solar wind, are ubiquitous in spacecraft observations, yet their origin remains elusive. In this study, we analyze data from the first 19 encounters of the Parker Solar Probe and report the dominance of 2 minute oscillations near the Alfvén surface. The frequency-rectified trace magnetic power spectral density (PSD) of these oscillations indicates that the fluctuation energy is concentrated around 2 minutes for the “youngest” solar wind. Further analysis using wavelet spectrograms reveals that these oscillations primarily consist of outward-propagating, spherically polarized Alfvén wave bursts. Through Doppler analysis, we show that the wave frequency observed in the spacecraft frame can be mapped directly to the launch frequency at the base of the corona, where previous studies have identified a distinct peak around 2 minutes (~8 mHz) in the spectrum of swaying motions of coronal structures observed by the Solar Dynamics Observatory Atmospheric Imaging Assembly. These findings strongly suggest that the Alfvén waves originate from the solar atmosphere. Furthermore, statistical analysis of the PSD deformation beyond the Alfvén surface supports the idea of dynamic formation of the otherwise absent 1/frange in the solar wind turbulence spectrum. 

Abstract In this Letter, we report observations of magnetic switchback (SB) features near 1 au using data from the Wind spacecraft. These features appear to be strikingly similar to the ones observed by the Parker Solar Probe mission closer to the Sun: namely, one-sided spikes (or enhancements) in the solar-wind bulk speed V that correlate/anticorrelate with the spikes seen in the radial-field component B R . In the solar-wind streams that we analyzed, these specific SB features near 1 au are associated with large-amplitude Alfvénic oscillations that propagate outward from the Sun along a local background (prevalent) magnetic field B 0 that is nearly radial. We also show that, when B 0 is nearly perpendicular to the radial direction, the large-amplitude Alfvénic oscillations display variations in V that are two sided (i.e., V alternately increases and decreases depending on the vector Δ B = B − B 0 ). As a consequence, SBs may not always appear as one-sided spikes in V , especially at larger heliocentric distances where the local background field statistically departs from the radial direction. We suggest that SBs can be well described by large-amplitude Alfvénic fluctuations if the field rotation is computed with respect to a well-determined local background field that, in some cases, may deviate from the large-scale Parker field.