pith. sign in

arxiv: 2512.07446 · v2 · submitted 2025-12-08 · 🌌 astro-ph.SR · physics.plasm-ph

The influence of Parker spiral on the reflection-driven turbulence

Khurram Abbas , Jonathan Squire This is my paper

Pith reviewed 2026-05-17 00:36 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.plasm-ph
keywords reflection-driven turbulenceParker spiralsolar wind heatingMHD turbulenceexpanding box simulationsAlfven wave reflectioncross helicity
0
0 comments X

The pith

The Parker spiral changes reflection-driven turbulence by reducing outer scales and slowing the cascade freeze-out in the solar wind.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper extends models of reflection-driven turbulence, which explains solar wind heating beyond adiabatic cooling, to account for the Parker spiral magnetic field geometry. Using three-dimensional expanding-box MHD simulations, it demonstrates that the azimuthal field component intersects perpendicularly stretched eddies, resulting in smaller outer scales perpendicular to the magnetic field. Consequently, the outer-scale nonlinear turnover time grows more slowly with distance, reducing the tendency for the turbulence to become quasi-static and magnetically dominated. This leads to greater dissipation of fluctuation energy as heat and allows the turbulence to remain strongly imbalanced, with high cross-helicity, to larger heliocentric distances. The work also characterizes spectra, switchbacks, and compressive fractions to offer testable predictions for observations.

Core claim

As the azimuthal field grows in the Parker spiral, it cuts across perpendicularly stretched, pancake-like eddies, producing outer scales perpendicular to the magnetic field that are much smaller than in the radial-background case. The outer-scale nonlinear turnover time therefore increases more slowly with heliocentric distance, weakening the cascade's tendency to freeze into quasi-static, magnetically dominated structures. This permits the system to dissipate a larger fraction of the fluctuation energy as heat while keeping the turbulence strongly imbalanced out to larger distances.

What carries the argument

Parker spiral geometry in expanding-box MHD simulations that alters eddy scales and nonlinear turnover times in reflection-driven turbulence.

If this is right

  • More fluctuation energy is converted to heat in the solar wind.
  • Turbulence maintains high normalized cross-helicity to greater heliocentric distances.
  • Observable differences in turbulence spectra and switchback properties compared to radial field models.
  • Improved agreement with measured temperature profiles in the heliosphere.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This suggests that incorporating realistic magnetic field geometry is essential for accurate solar wind heating models.
  • Future simulations could test how these effects interact with kinetic plasma processes not captured in MHD.
  • The slower freeze-out might influence the generation or survival of switchbacks in the solar wind.

Load-bearing premise

The expanding-box MHD approximation with an imposed Parker spiral captures the essential dynamics of real solar-wind turbulence without major interference from numerical effects or omitted physics.

What would settle it

Spacecraft measurements of the radial dependence of cross-helicity or turbulent heating rates that deviate from the predicted slower increase in turnover time and greater dissipation in Parker spiral geometry.

Figures

Figures reproduced from arXiv: 2512.07446 by Jonathan Squire, Khurram Abbas.

Figure 1
Figure 1. Figure 1: Schematic of an expanding plasma parcel in the solar wind. (Top) Geometry [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of key parameters for waves under solar-wind expansion, plotted [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Snapshots of the Elsässer fields |z ± ⊥| perpendicular to magnetic field in the y-z plane. The top two rows (a) show different stages of expansion are shown for the A05- Φ0 = 0◦ with a radial B simulation at a ≈ 6 (left), a ≈ 22.35 (middle), and a ≈ 50.35 (right). These snapshots illustrate the turbulent evolution from an initially imbalanced regime to a magnetically dominated and balanced phase. The botto… view at source ↗
Figure 4
Figure 4. Figure 4: 3D visualizations of the Elsässer fields [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of normalized Elsässer wave action energies, [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of correlation length ℓ as a function of a for magnetic-field (B⊥) fluctuations for MR-Φ0 = 0◦ and MR-Φ0 = 2◦ runs. Each curve represents projections along the field (∥, blue), perpendicular(⊥, brown), eˆT direction (⊥, T, green), and eˆN direction (⊥, N, orange). The dashed black line indicates the a 1 power-law scaling for reference. In the radial case (left panel), ℓ⊥ grows nearly linearly wit… view at source ↗
Figure 7
Figure 7. Figure 7: The time-scale ratio χexp plotted as a function of expansion factor a. We computed this using Alfvénic perpendicular correlation lengths for magnetic field fluctuations for HR-Φ0 = 0◦ , 4 ◦ (left) and MR-Φ0 = 0◦ , 2 ◦ , 5 ◦ simulations (right). The solid lines show χexp calculated from ℓ⊥,T , and the dotted lines show χexp computed from ℓ⊥,N . In the radial case, we averaged the both to compute ℓ⊥, and χex… view at source ↗
Figure 8
Figure 8. Figure 8: Measured χA versus a. Note that, we take the average of transverse and normal perpendicular correlation lengths ℓ⊥ = (ℓ⊥,T + ℓ⊥,N ) for this figure. do we expect nonlinear interactions to dominate, allowing a strong turbulent cascade to develop and damp fluctuations as assumed. In the opposite limit of χA ≳ 1 fluctuations are expected to rapidly de-correlate in k∥ and drop to χA ∼ 1 (Schekochihin 2022). In… view at source ↗
Figure 9
Figure 9. Figure 9: Ratio of wave energies E˜+/E˜− plotted against χexp for all simulations (HR￾Φ0 = 0◦ , 4 ◦ and MR-Φ0 = 0◦ , 2 ◦ , 5 ◦ , A05-Φ0 = 0◦ , 2 ◦ , 5 ◦ ). The thick lines denote the perpendicularly averaged values for each case, while the faint solid and dashed curves correspond to the eˆT and eˆN components, respectively. The black dashed line indicates the theoretical scaling corresponding to the prediction from … view at source ↗
Figure 10
Figure 10. Figure 10: Perpendicular energy spectra in the HR-Φ0 = 0◦ and HR-Φ0 = 4◦ simulations. All spectra are plotted versus dimensionless perpendicular wavenumber k⊥L⊥ (L⊥ is the comoving perpendicular box length). Panel (a): Radial field evolution of outward and inward energy Elsässer spectra E ± at the indicated expansion times, showing the development of imbalance with scale. (b) At a = 20.5 the magnetic (EM, purple) an… view at source ↗
Figure 11
Figure 11. Figure 11: Parametric evolution of cross helicity σc and residual energy σr as a function of expansion factor a. Panel (a) shows results from the HR-Φ0 = 0◦ , 4 ◦ simulations: circles correspond to the purely radial run and squares to the PS run. Panel (b) shows the MR-Φ0 = 0◦ , 10◦ simulations for both the radial and the higher initial Parker-angle case (Φ0 = 10◦ ). Colored points represent a = 1 (dark blue), a ∼ 1… view at source ↗
Figure 12
Figure 12. Figure 12: Fly-throughs of the magnetic field for the high-resolution radial (HR– [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Evolution of switchback fraction (fz ⩾ zth) in MR-Φ0 = 0◦ (solid) and MR-Φ0 = 2 ◦ (dashed lines). Both runs produce switchbacks, but the PS case exhibits systematically larger switchback fractions across effectively all a and a stronger growth of large-angle deflections with a up to the point where fluctuation amplitudes decline; the downturn in (fz ⩾ zth) at the largest a is caused by the overall decreas… view at source ↗
Figure 14
Figure 14. Figure 14: Evolution of (a) magnetic compressibility [PITH_FULL_IMAGE:figures/full_fig_p027_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Solutions to the linear expanding incompressible MHD equations, for radial [PITH_FULL_IMAGE:figures/full_fig_p030_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Snapshot of |z ± ⊥|/|z ± ⊥|rms in the y–z plane at a = 39 for MR-Φ0 = 5◦ run showing the onset of grid-scale, speckle-like noise at large spiral angle. 10 0 10 1 10 2 a 10 -3 10 -2 10 -1 10 0 ~E § = ~E + 0 numerical issue ©0 = 0 ©0 = 2 ©0 = 5 ©0 = 10 [PITH_FULL_IMAGE:figures/full_fig_p032_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Time series of the outward energy for the MR- [PITH_FULL_IMAGE:figures/full_fig_p032_17.png] view at source ↗
read the original abstract

The solar wind is observed to undergo substantial heating as it expands through the heliosphere, with measured temperature profiles exceeding those expected from adiabatic cooling. A plausible source of this heating is reflection-driven turbulence (RDT), in which gradients in the background Alfv\'en speed partially reflect outward-propagating Alfv\'en waves, seeding counter-propagating fluctuations that interact and dissipate via turbulence. Previous RDT models assume a radial background magnetic field, but at larger radii the interplanetary field is known to be twisted into the Parker Spiral (PS). Here, we generalize RDT phenomenology to include a PS, using three-dimensional expanding-box magnetohydrodynamic (MHD) simulations to test the ideas and compare the resulting turbulence to the radial-background-field case. We argue that the underlying RDT dynamics remain broadly similar with a PS, but the controlling scales change: as the azimuthal field grows it "cuts across" perpendicularly stretched, pancake-like eddies, producing outer scales perpendicular to the magnetic field that are much smaller than in the radial-background case. Consequently, the outer-scale nonlinear turnover time increases more slowly with heliocentric distance in PS geometry, weakening the tendency (seen in radial-background models) for the cascade to 'freeze' into quasi-static, magnetically dominated structures. This allows the system to dissipate a larger fraction of the fluctuation energy as heat, also implying that the turbulence remains strongly imbalanced (with high normalized cross-helicity) out to larger heliocentric distances. We complement our heating results with a detailed characterization of the turbulence (e.g., spectra, switchbacks, and compressive fractions) providing a set of concrete predictions for comparison with spacecraft observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents 3-D expanding-box MHD simulations of reflection-driven turbulence (RDT) in the solar wind, generalizing prior radial-background models to include the Parker spiral (PS). The central claim is that the growing azimuthal field component intersects perpendicularly stretched, pancake-like eddies, yielding outer scales perpendicular to the local magnetic field that are substantially smaller than in the radial case. This leads to slower growth of the outer-scale nonlinear turnover time with heliocentric distance, weakening the tendency for the cascade to freeze into quasi-static structures, increasing the fraction of fluctuation energy dissipated as heat, and allowing strong imbalance (high normalized cross-helicity) to persist to larger radii. The work also reports characterizations of spectra, switchbacks, and compressive fractions as observational predictions.

Significance. If the reported differences in perpendicular scales and resulting dissipation fractions are robust, the results would meaningfully extend RDT phenomenology to realistic interplanetary field geometry and offer concrete, falsifiable predictions for spacecraft data on heating and turbulence imbalance. The direct numerical simulation approach supplies quantitative outputs rather than fitted parameters, strengthening the work's utility for solar-wind modeling.

major comments (2)
  1. [Abstract and results section] Abstract and results section: the assertion that azimuthal-field geometry produces 'much smaller' outer scales perpendicular to B (and consequently slower turnover-time growth) must be supported by explicit quantitative metrics, such as measured correlation lengths, outer-scale spectra, or structure functions in the plane perpendicular to the local magnetic field. Without these direct comparisons between the radial and PS runs, the chain of inferences on reduced freezing, higher heating fractions, and sustained imbalance does not follow from the reported dynamics.
  2. [Methods section] Methods section: the expanding-box MHD approximation with an imposed Parker spiral is the weakest assumption; a resolution study, explicit check on numerical dissipation relative to physical dissipation, and discussion of missing kinetic effects are needed to confirm that the reported heating fractions and scale evolution are not contaminated by the numerical setup.
minor comments (2)
  1. [Introduction] Notation for the normalized cross-helicity and compressive fraction should be defined explicitly on first use and kept consistent with prior RDT literature.
  2. [Figures] Figure captions for the turbulence visualizations should include the exact heliocentric distances and the orientation of the local mean field relative to the plotted plane.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed report. The comments help clarify the presentation of our quantitative results and the limitations of the expanding-box approach. We address each major comment below and have revised the manuscript to incorporate additional metrics, convergence tests, and expanded discussion.

read point-by-point responses

  1. Referee: [Abstract and results section] Abstract and results section: the assertion that azimuthal-field geometry produces 'much smaller' outer scales perpendicular to B (and consequently slower turnover-time growth) must be supported by explicit quantitative metrics, such as measured correlation lengths, outer-scale spectra, or structure functions in the plane perpendicular to the local magnetic field. Without these direct comparisons between the radial and PS runs, the chain of inferences on reduced freezing, higher heating fractions, and sustained imbalance does not follow from the reported dynamics.

    Authors: We agree that the original manuscript would benefit from more explicit quantitative support for the claimed reduction in perpendicular outer scales. In the revised version we have added direct measurements of two-point correlation lengths and second-order structure functions evaluated in the plane perpendicular to the local magnetic field for both the radial and Parker-spiral runs. These diagnostics confirm that the perpendicular outer scales are reduced by a factor of approximately 1.6–2.0 in the Parker-spiral geometry at 1 AU, directly supporting the slower growth of the nonlinear turnover time and the consequent changes in heating fraction and cross-helicity evolution. The abstract and results section have been updated to cite the new figures and to state the measured scale ratios explicitly. revision: yes

  2. Referee: [Methods section] Methods section: the expanding-box MHD approximation with an imposed Parker spiral is the weakest assumption; a resolution study, explicit check on numerical dissipation relative to physical dissipation, and discussion of missing kinetic effects are needed to confirm that the reported heating fractions and scale evolution are not contaminated by the numerical setup.

    Authors: We have performed additional resolution studies at 1.5 times the original grid resolution and confirm that the reported heating fractions, perpendicular scale evolution, and cross-helicity profiles converge to within 8 %. We have also added an explicit comparison of numerical versus physical dissipation by tracking the energy flux through the inertial range and estimating the numerical dissipation scale; this shows that numerical dissipation remains sub-dominant inside the inertial range for the resolutions used. Regarding kinetic effects, we have expanded the discussion section to acknowledge that the MHD approximation omits Landau damping, cyclotron resonance, and other kinetic processes that could modify the dissipation range; we note that these effects are expected to become important only at scales smaller than those resolved here and suggest that future hybrid or kinetic simulations would be needed to quantify their impact on the global heating budget. revision: partial

Circularity Check

0 steps flagged

No significant circularity: results are direct outputs of expanding-box MHD simulations

full rationale

The paper's central claims about changes in perpendicular outer scales, slower growth of nonlinear turnover times, reduced freezing into quasi-static structures, higher dissipation fractions, and sustained imbalance are obtained as quantitative outputs from three-dimensional expanding-box MHD simulations with imposed Parker spiral versus radial field. These are not derived analytically from fitted parameters or prior self-citations in a way that reduces to the inputs by construction. The geometric argument that the azimuthal field 'cuts across' pancake-like eddies is a qualitative interpretation of the imposed background field evolution in the expanding box, but the measured scales, spectra, and energy dissipation are simulation results. No self-definitional loops, fitted inputs renamed as predictions, or load-bearing uniqueness theorems from overlapping-author citations appear in the derivation chain. The work is self-contained against external benchmarks via the reported simulation diagnostics.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available; therefore the ledger is populated from the high-level statements in the abstract. The central claim rests on standard MHD assumptions plus the expanding-box approximation and the imposed Parker-spiral background.

axioms (2)
  • standard math Ideal MHD equations govern the plasma evolution
    Implicit in any MHD simulation description; invoked throughout the modeling section referenced in the abstract.
  • domain assumption Expanding-box approximation correctly captures radial expansion effects on turbulence
    Stated as the simulation framework used to test RDT with Parker spiral.

pith-pipeline@v0.9.0 · 5599 in / 1395 out tokens · 71459 ms · 2026-05-17T00:36:42.395349+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

  1. [1]

    D., Badman, S

    Bale, S. D., Badman, S. T., Bonnell, J. W., Bowen, T. A., Burgess, D., Case, A. W., Cattell, C. A., Chandran, B. D. G., Chaston, C. C., Chen, C. H. K., Drake, J. F., de Wit, T. Dudok, Eastwood, J. P., Ergun, R. E., F arrell, W. M., Fong, C., Goetz, K., Goldstein, M., Goodrich, K. A., Har vey, P. R., Horbury, T. S., Howes, G. G., Kasper, J. C., Kellogg, P....

    work page 2019

  2. [2]

    Ba v assano, B., Pietropaolo, E

    Barnes, Aaron & Holl weg, Joseph V.1974 Large-amplitude hydromagnetic waves.Journal of Geophysical Research79(16), 2302–2318. Ba v assano, B., Pietropaolo, E. & Bruno, R.1998 Cross-helicity and residual energy in solar wind turbulence: Radial evolution and latitudinal dependence in the region from 1 to 5 au.Journal of Geophysical Research: Space Physics10...

    work page 1974

  3. [3]

    Chen, C. H. K., Chandran, B. D. G., Woodham, L. D., Jones, S. I., Perez, J. C., Bourouaine, S., Bowen, T. A., Klein, K. G., Moncuquet, M., Kasper, J. C. & Bale, S. D.2021 The near-sun streamer belt solar wind: turbulence and solar wind acceleration.Astronomy &; Astrophysics650, L3. Cranmer, Steven R.2009 Coronal holes.Living Reviews in Solar Physics6. Cra...

    work page 2021

  4. [4]

    H., Milano, L

    Dmitruk, P., Matthaeus, W. H., Milano, L. J., Oughton, S., Zank, G. P. & Mullan, D. J.2002 Coronal heating distribution due to low-frequency, wave-driven turbulence. The Astrophysical Journal575(1), 571–577. Goldreich, P. & Sridhar, S.1995 Toward a theory of interstellar turbulence. 2: Strong Alfvénic turbulence.The Astrophysical Journal438,

    work page 2002

  5. [5]

    Halekas, J

    Grappin, Roland, Velli, Marco & Mangeney, André1993 Nonlinear wave evolution in the expanding solar wind.Physical Review Letters70(14), 2190–2193. Halekas, J. S., Bale, S. D., Berthomier, M., Chandran, B. D. G., Drake, J. F., Kasper, J. C., Klein, K. G., Larson, D. E., Livi, R., Pulupa, M. P., Stevens, M. L., Verniero, J. L. & Whittlesey, P.2023 Quantifyi...

    work page 2023

  6. [6]

    & Olbert, S.1980 Non-wkb Alfvén waves in the solar wind.Journal of Geophysical Research: Space Physics85(A3), 1311–1327

    Heinemann, M. & Olbert, S.1980 Non-wkb Alfvén waves in the solar wind.Journal of Geophysical Research: Space Physics85(A3), 1311–1327. Horbury, T. S., Bale, Stuart D., McManus, Michael D., Larson, Da vin, Kasper, J. C., Laker, Ronan, Matteini, Lorenzo, Raouafi, Nour E., Velli, Marco, Woodham, Lloyd D., Woolley, Thomas, Fedorov, Andrey, Louarn, Philippe, K...

    work page 1980

  7. [7]

    G., Klein, Kristopher G

    Perez, Jean C., Chandran, Benjamin D. G., Klein, Kristopher G. & Martinović, Mihailo M.2021 How Alfvén waves energize the solar wind: heat versus work.Journal of Plasma Physics87(2). Podesta, J. J.2009 Dependence of solar-wind power spectra on the direction of the local mean magnetic field.The Astrophysical Journal698(2), 986–999. Raouafi, N. E., Matteini...

    work page 2021

  8. [8]

    Squire, J., Chandran, B. D. G. & Meyrand, R.2020 In-situ switchback formation in the expanding solar wind.The Astrophysical Journal Letters891(1), L2. Squire, Jonathan, Johnston, Zade, Mallet, Alfred & Meyrand, Romain2022 On the properties of Alfvénic switchbacks in the expanding solar wind: The influence of the Parker spiral.Physics of Plasmas29(11). Sto...

    work page 2020

  9. [9]

    Solar Physics109(1), 149–186

    Tu, Chuan-Yi1987 A solar wind model with the power spectrum of Alfvénic fluctuations. Solar Physics109(1), 149–186. Tu, C. Y. & Marsch, E.1995 MHD structures, waves and turbulence in the solar wind: Observations and theories.Space Science Reviews73(1–2), 1–210. Velli, M.1993 On the propagation of ideal, linear Alfvén waves in radially stratified stellar a...

    work page 1995

  10. [10]

    Verdini, A., Velli, M

    Verdini, Andrea & Velli, Marco2007 Alfvén waves and turbulence in the solar atmosphere and solar wind.The Astrophysical Journal662(1), 669–676. Verdini, A., Velli, M. & Buchlin, E.2008 Reflection driven MHD turbulence in the solar atmosphere and wind.Earth, Moon, and Planets104(1–4), 121–125. Verdini, A., Velli, M. & Buchlin, E.2009 Turbulence in the sub-...

    work page 2008

  11. [11]

    T., Roberts, D

    Wicks, R. T., Roberts, D. A., Mallet, A., Schekochihin, A. A., Horbury, T. S. & 36K. Abbas and J. Squire Chen, C. H. K.2013 Correlations at large scales and the onset of turbulence in the fast solar wind.The Astrophysical Journal778(2), 177

    work page 2013