The team analyzed measurements taken by Solar Orbiter on 2 March 2022, when the spacecraft was roughly 0.55 astronomical units from the Sun, about 120 solar radii. During a six-hour window, the spacecraft recorded five intervals in which the radial component of the magnetic field reversed direction. Crucially, the suprathermal electron beam, the population of energetic electrons that streams away from the Sun along magnetic field lines, continued to point sunward throughout these reversals. That detail matters enormously: if the field had genuinely flipped its polarity, the electron beam would have reversed too. Instead, the field was folded back on itself, the defining signature of a switchback rather than a true sector boundary crossing.

What elevates this event from another statistical entry in the switchback catalog is the heavy ion data. Solar Orbiter’s Solar Wind Analyser, equipped with a Heavy Ion Sensor, measured the charge states of oxygen and carbon ions within and around the switchback. The ratios of O7+ to O6+, C6+ to C4+, and C6+ to C5+ all rose sharply and simultaneously with the largest switchback, which spanned from roughly 20:25 to 21:55 UTC. Because ions freeze in their charge state once the coronal density drops too low for further ionization or recombination, these ratios act as a fossil thermometer, recording the electron temperature of the plasma’s source region millions of kilometers back along the solar wind stream.

The elevated charge-state ratios pointed to a source region far hotter than the typical open magnetic field of a coronal hole. Adding to the picture, the ratio of helium ions to protons dropped measurably inside the switchback. Helium depletion is a well-known signature of plasma that has spent time on closed coronal loops, where wave heating and collisional transport fractionate elements differently than on open field lines. Taken together, the hot ion charge states and the helium deficit indicated that the switchback plasma had originated in hot, closed magnetic loop structures rather than in the cooler, open flux of the coronal hole itself.

To pin down the source more precisely, the researchers combined magnetic field modeling with remote sensing. Using a potential field source surface model coupled to a Parker spiral, they back-traced the field lines connecting Solar Orbiter to the solar surface, testing source surface heights between 1.5 and 3 solar radii. The footpoints converged on a large equatorial coronal hole, consistent with earlier connectivity studies for this period. Meanwhile, the Solar Dynamics Observatory’s Atmospheric Imaging Assembly imaged the same portion of the Sun in six extreme ultraviolet passbands, allowing the team to apply a differential emission measure technique and extract the most probable electron temperature and density of the low corona at the connected footpoints.

Those remote measurements then fed a semi-empirical coronal model. The team set the measured temperature and density, around 1.1 to 1.2 million kelvin, as the lower boundary condition of a coronal hole flux tube model and solved the non-equilibrium ionization equation outward through the corona, using ionization and recombination rates from the CHIANTI atomic database. The predicted frozen-in charge-state ratios matched the values Solar Orbiter measured at the peak of the switchback. Mapping the model across every pixel of the solar image revealed two candidate source regions: the boundary of the coronal hole, or small magnetic loops embedded within the coronal hole itself. Applying a constant-pressure scaling law, the team estimated the loop lengths at roughly 8 to 12 megameters, about a third of a supergranule, consistent with the observed size of coronal bright points.

The physical interpretation that emerges is interchange reconnection. In this process, closed magnetic loops carrying hot plasma reconnect with open field lines of the coronal hole, releasing the loop plasma into the solar wind and exchanging the connectivity of the field. The hot, helium-depleted plasma measured inside the switchback is exactly what interchange reconnection would deliver, and the mechanism has previously been invoked to explain the fast solar wind emerging from coronal holes. The loop sizes found here, tens of megameters, align with models in which interchange reconnection injects both energy and mass into the corona to create the solar wind.

Yet the story does not end at the Sun. The team also examined how the switchback propagated through interplanetary space. Inside the major switchback, the radial proton speed increased by approximately 50 kilometers per second, nearly equal to the local Alfvén speed, while the magnetic field strength remained roughly constant. The relative drift speeds between heavy ions and protons changed sign across the switchback boundaries, with all differences remaining below the local Alfvén speed. These observations are broadly consistent with an outward-propagating, large-amplitude, spherically polarized Alfvén-like fluctuation, suggesting that once launched by reconnection, the switchback travels outward as a wave-like structure, riding on the solar wind at roughly the Alfvén speed above the bulk flow.

This dual character, a reconnection birth and an Alfvénic propagation, is perhaps the study’s most consequential message. The solar physics community has long been split between two paradigms for heating the corona and accelerating the solar wind: reconnection on one side, and waves and turbulence on the other. The new results suggest these mechanisms are not mutually exclusive competitors but partners in a single system. Interchange reconnection appears to govern where and how switchback plasma is released from the corona, while wave and turbulence physics governs how the resulting structures propagate and evolve. Switchbacks, in this view, are a living demonstration of the elaborate interplay between the two mechanisms.

The authors caution that this is a single case study, and that the analyzed switchback was unusually long-lived compared with the typical duration distribution. Larger statistical surveys, including intervals with multiple coronal holes on the solar disk and future missions carrying more capable heavy ion instrumentation, will be needed to establish how often switchbacks preserve their coronal fingerprints. Still, the implications are striking. If switchbacks reliably retain measurable signatures of their source, they become messengers from regions of the solar atmosphere that are extraordinarily difficult to observe directly, encoding information about heating in the chromosphere and corona, kinetic energization of ions, and the fundamental energy transfer processes that drive the solar wind. Tracing these folded field lines from the photosphere, where the energy input begins, to the point where spacecraft intercept them may finally allow scientists to quantify the balance between reconnection and turbulence that powers our star’s ever-flowing wind.

Subject of Research: The coronal origin of magnetic switchbacks in the solar wind traced by Solar Orbiter heavy ion observations

Article Title: On the coronal origin of magnetic switchbacks in the solar wind

Article References: Coburn, J. T., Yardley, S. L., Ngampoopun, N., Dewey, R., Rivera, Y. J., De Marco, R., Verscharen, D., Owen, C. J., Suen, G. H. H., Alterman, B., Baker, D., Fazakerley, A., Galvin, A. B., He, J., Horbury, T. S., Hou, C., Ioannou, C., Kistler, L. M., Lepri, S. T., … Rouillard, A. P. (2026). On the coronal origin of magnetic switchbacks in the solar wind. Nature Astronomy. https://doi.org/10.1038/s41550-026-02928-0

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02928-0

Keywords: solar wind, magnetic switchbacks, Solar Orbiter, corona, interchange reconnection, Alfvén waves, heavy ions, charge states, coronal holes, coronal loops, heliophysics, Nature Astronomy

The study, led by Etienne Gavazzi of UiT The Arctic University of Norway in Tromsø and published in Annales Geophysicae, mines measurements from the VISIONS-2 sounding rocket mission. Two rockets were launched from Ny-Ålesund on Svalbard on 7 December 2018, flying south-west through the dayside auroral region at around 13.5 magnetic local time. The lower of the two payloads, which reached an apogee of just over 600 kilometers, carried a top-hat electrostatic analyser capable of sweeping through 49 energy steps from 3 electronvolts to 28.6 kiloelectronvolts every 50 milliseconds, while simultaneously resolving the directions of incoming electrons across twenty pitch-angle channels. That combination of speed and angular resolution proved decisive.

The physical culprit behind these auroras is a family of electromagnetic disturbances known as dispersive Alfvén waves. Ordinary magnetohydrodynamic Alfvén waves cannot carry a parallel electric field, the ingredient needed to accelerate electrons along magnetic field lines. But when the waves acquire small perpendicular scales, comparable to the electron inertial length in the tenuous plasma below a few Earth radii, or to the ion acoustic gyroradius farther out, they become dispersive and develop exactly that field-aligned electric field. Through resonant processes, electrons surfing these waves can be driven to parallel velocities approaching, and in some regimes exceeding, the local Alfvén speed, producing the kiloelectronvolt electrons that light up the sky. Previous spacecraft statistics suggested that Alfvén waves power on average a quarter to nearly forty percent of all auroral electron energy, dominating under active conditions near noon and midnight.

The fingerprint of this wave-driven precipitation is distinctive. In time-energy spectrograms, the electrons arrive broadband in energy and dispersed in time: the fastest, most energetic electrons reach the spacecraft first, with progressively slower ones trailing behind. The VISIONS-2 low-flyer recorded dozens of such dispersed structures during its flight, with peak energies from a few hundred electronvolts to several kiloelectronvolts and durations typically under half a second. Crucially, the same intervals showed semicircular dispersions in pitch-angle, where field-aligned electrons arrived first and arrival time grew with angle, together with broadband electromagnetic fluctuations and bursts of downward Poynting flux. A cross-spectral analysis of the perpendicular electric and magnetic fields yielded a ratio of about 2500 to 3000 kilometers per second between 1 and 8 hertz, closely matching the local Alfvén speed of roughly 2730 kilometers per second, and cross-phases consistent with waves propagating downward into the ionosphere.

To convert these fleeting signatures into altitudes, the team applied the classical time-of-flight technique: assuming electrons of all energies were released simultaneously from a single height, a linear fit of arrival time against inverse velocity yields the source distance. They also tested a quadratic fit, previously shown to describe dayside dispersions better, and introduced a logarithmic fit that appears as a straight line on semi-logarithmic energy-time plots. Using the corrected Akaike Information Criterion to compare models across 29 well-resolved events, the non-linear fits won out overwhelmingly, with 14 events best described by the logarithmic curve, 12 by the quadratic, and only 3 by the linear model. That statistical preference itself carries meaning: it supports a picture in which acceleration is smeared across a range of altitudes as the wave travels downward, rather than confined to a single point.

The real methodological innovation, however, lies in a second approach that the authors developed to escape the simultaneous-release assumption altogether. Instead of fitting energy-time curves, they built a forward model of pitch-angle-time dispersions. Electrons of a single energy are launched from a candidate source altitude with a spread of pitch angles, their trajectories computed under magnetic mirroring using the IGRF-14 field model, and their arrival times convolved with a Gaussian injection profile and the detector’s sampling behaviour. Three parameters, source altitude, injection duration, and a time offset, are then adjusted for each energy channel independently by minimizing the squared error between modelled and observed dispersions. Because each energy bin receives its own release height, the method relaxes the core assumption of the energy-time techniques while exploiting the rocket’s unusually fine angular resolution.

The two families of methods, resting on different assumptions, converged. For most analysed structures, the inferred source altitudes fell between 1000 and 3000 kilometers, with a clear trend across events that higher-energy electrons were released higher up: a few hundred electronvolts mapped to roughly 1000 to 3000 kilometers, while kiloelectronvolt electrons traced back to about 2000 to 5000 kilometers. Where the methods disagreed, the differences rarely exceeded 1000 kilometers. The fitted time offsets also increased systematically for lower-energy electrons, indicating later release times, exactly what one expects if the acceleration region rides downward with a propagating wave. One event at 500 seconds flight time even yielded a source altitude near 6400 kilometers for 4300-electronvolt electrons, with modelled source pitch angles confined within about 25 degrees of field alignment.

Placed against theory, the results sit at the lower end of expectations but remain physically coherent. When the team compared their altitude-energy scatter with inertial Alfvén wave velocity profiles computed from published plasma density models, the observations aligned best with profiles featuring relatively small oxygen ion scale heights, meaning the plasma transitions from heavy O+ dominance to light H+ dominance at comparatively low altitudes. Almost none of the events matched the outside-cavity profile of Lysak and Song, while many were consistent with profiles associated with density cavities, the depleted plasma regions long linked to intense parallel electric fields in auroral zone observations. The event-to-event variability in source altitude, even within a single rocket flight, points to genuinely different plasma and wave conditions along neighbouring magnetic field lines, possibly structured by the ionospheric Alfvén resonator or by the aftermath of earlier precipitation.

Perhaps the most exciting prospect is what the technique makes measurable. By fitting the inferred release altitudes and energies to a parametrized inertial Alfvén speed profile, the authors showed that the oxygen and hydrogen density profiles and the perpendicular wavelength of the wave can be reconstructed from rocket data alone. Applied across their events, about two-thirds produced reasonably constrained fits, with best-fit parameters remarkably close to values reported in earlier statistical studies. Because the topside ionosphere near the O+-to-H+ transition is one of the most poorly sampled regions of near-Earth space, a method that turns a brief rocket traversal into a density profile and a wave scale is a genuinely valuable addition to the toolkit. Future instruments with higher time, energy, and angular resolution could sharpen these estimates further, and the approach could be extended to satellites, trading some precision for coverage across a far wider range of geomagnetic conditions.

For now, the study delivers something auroral physics has long wanted: a consistent, multi-method measurement of where the sky’s most restless lights get their power. The electrons that set the dayside aurora ablaze appear to be flung Earthward from effective release altitudes of one to three thousand kilometers most of the time, energized progressively by inertial Alfvén waves as those waves sweep down magnetic field lines. Each dispersed burst recorded over half a second of rocket flight is thus a snapshot of a wave-particle interaction unfolding thousands of kilometers overhead, and with the new pitch-angle technique, that snapshot can now be read as a probe of the plasma environment itself, turning the aurora from a spectacle into an instrument.

Subject of Research: Estimating the source altitudes of Alfvénic auroral electron acceleration in the dayside ionosphere using sounding rocket dispersion measurements

Article Title: Estimating the source altitude of auroral precipitation from dispersed Alfvén waves in the dayside ionosphere

Article References: Gavazzi, E., Spicher, A., Gustavsson, B., Vierinen, J., Clemmons, J., Pfaff, R., & Rowland, D. (2026). Estimating the source altitude of auroral precipitation from dispersed Alfvén waves in the dayside ionosphere. Annales Geophysicae, 44(2), 903-920. https://doi.org/10.5194/angeo-44-903-2026

Image Credits: AI Generated

DOI: 10.5194/angeo-44-903-2026

Keywords: aurora, Alfvén waves, sounding rockets, VISIONS-2, space plasma physics, electron acceleration, dayside ionosphere, pitch-angle dispersion, time-of-flight analysis, magnetosphere-ionosphere coupling, plasma density profiles, Annales Geophysicae

The scale separation in space is staggering. Reynolds numbers, which measure the ratio of inertial to viscous effects, routinely reach 10 billion or larger in astrophysical flows, meaning turbulence is essentially unavoidable. Unlike water or air, astrophysical plasmas are almost perfectly conducting, so their dynamics are governed by the MHD equations, coupling currents, magnetic fields and the Lorentz force. In our Galaxy, magnetic fields of around 5 microgauss reach rough equipartition with turbulent kinetic energy; in galaxy clusters, fields of 1 to 3 microgauss sit near a twentieth of equipartition. These fields are not decorative: they accelerate particles, fill the Universe with non-thermal radiation, and even enable the accretion of matter onto black holes through the magnetorotational instability, a process estimated to be the most potent energy source in the cosmos, exceeding thermonuclear burning in stars.

The mathematical backbone of the subject remains the Kolmogorov cascade picture. Energy injected at large scales passes, without dissipation, through an inertial range of scales until viscosity finally wins at the Kolmogorov scale. Dimensional analysis then yields the famous spectrum in which energy content scales as the wavenumber to the minus five-thirds power, with a dimensionless Kolmogorov constant that experiments and simulations place near 1.6. Beryesnyak emphasizes a rigorous tool called scaling convergence: because the underlying equations contain no preferred scale, small-scale statistics from simulations of different resolutions should collapse onto a single universal curve when properly normalized, a technique that drives statistical error virtually to zero and has even resolved the tiny intermittency correction of about 0.04 to the spectral slope.

MHD turbulence, however, refuses to behave like its hydrodynamic cousin. A large-scale magnetic field cannot be transformed away, and it remains dynamically important on every scale. Linearizing the MHD equations reveals four wave modes, of which the transverse Alfvén mode dominates the nonlinear cascade. Early theorists Iroshnikov and Kraichnan imagined weak, wave-like interactions, but later work showed that turbulence becomes stronger, not weaker, as it cascades. Energy transfer proceeds preferentially perpendicular to the field, producing extreme anisotropy. Goldreich and Sridhar proposed that this anisotropy is capped by critical balance, where the cascade time matches the wave period, yielding a perpendicular spectrum of minus five-thirds and the relation that the parallel wavenumber scales as the perpendicular wavenumber to the two-thirds power.

High-resolution direct numerical simulations have now put these ideas to demanding tests. Using the scaling convergence method on simulations up to 4096 cubed grid points, Beresnyak found that the perpendicular spectrum converges best near a slope of minus 1.7, consistent with the Kolmogorov picture rather than competing minus three-halves models. The Alfvénic Kolmogorov constant was measured as 3.3 with a total value near 4.2 when the passively advected slow mode is included. Remarkably, the residual energy, the difference between magnetic and kinetic energy, turns out to be a constant fraction, about 15 percent, of the total energy throughout the inertial range, with a corresponding Alfvén ratio of roughly 0.74, resolving earlier conceptual difficulties with theories that predicted scale-dependent behavior.

An elegant theoretical result connects the parallel spectrum to Lagrangian statistics. Because oppositely directed Alfvén wave packets propagate along magnetic field lines at a fixed speed, measuring fluctuations along the field is mathematically equivalent to following a fluid element in time. This argument yields a parallel spectrum proportional to the wavenumber to the minus second power, scaled by the inverse of the Alfvén speed, without ever invoking critical balance. Numerical measurements along the local magnetic field overwhelmingly confirm this minus-two law, matching observations from the solar wind, where spacecraft such as Helios 2 have recorded clean power-law spectra over decades of frequency.

Imbalanced turbulence, where waves traveling one direction dominate, presents a harder puzzle, and it is the norm in the solar wind and near astrophysical jets. Because critical balance cannot hold simultaneously for counter-propagating waves of unequal amplitude, several competing models were proposed. Simulations with systematically varied imbalance show that the Lithwick-Goldreich-Sridhar model captures the spectra at small imbalances, while the Beresnyak-Lazarian model, which relaxes locality for the dominant component, best matches the energy ratios and the diverging anisotropies of the two populations at strong imbalance.

Perhaps the most consequential result concerns the small-scale dynamo, the process by which turbulence amplifies weak magnetic fields. Once the kinematic, exponential phase ends, magnetic energy grows linearly in time as turbulence converts a fixed fraction of cascade power into magnetism. The measured efficiency constant is small, about 0.05, but the implications are enormous: in galaxy clusters, this keeps the ratio of magnetic to thermal energy constant at roughly 40 over the past 10 billion years, in agreement with Faraday rotation observations. The review also issues a caution to simulators. Because numerical Reynolds numbers are vastly smaller than astrophysical ones, starting a simulated young object such as a collapsing cloud with a vanishing field can artificially delay magnetization, and implicit large-eddy codes with zero initial field produce no field at all, in gross contradiction with nature, where the dynamo always jump-starts itself.

The frontier now extends to supersonic turbulence in molecular clouds, where Mach numbers near 10 produce densities varying by orders of magnitude and log-normal probability distributions sculpted by slow-mode shocks and sheared by Alfvénic motions, with direct consequences for star formation theory. It extends, too, to magnetic reconnection, where current sheets tear and spawn their own strong, critically balanced turbulence, producing reconnection rates of about 1.5 percent of the Alfvén speed that are independent of resistivity. From solar flares to black hole jets to the magnetization of the cosmic web, MHD turbulence has become the connective tissue of modern astrophysics, and the convergence of theory, simulation and spacecraft measurement suggests the field is finally converging on universal answers.

Subject of Research: Theory and numerical simulation of magnetohydrodynamic turbulence and its astrophysical applications

Article Title: MHD turbulence

Article References: Beresnyak, A. (2019). MHD turbulence. Living Reviews in Computational Astrophysics, 5(1), Article 2. https://doi.org/10.1007/s41115-019-0005-8

Image Credits: AI Generated

DOI: 10.1007/s41115-019-0005-8

Keywords: MHD turbulence, magnetohydrodynamics, astrophysical plasmas, Kolmogorov cascade, Alfvén waves, turbulent dynamo, solar wind, galaxy clusters, magnetic reconnection, interstellar medium, numerical simulations, star formation