Alfven Wave Reflection and Turbulent Heating in the Solar Wind from 1 Solar Radius to 1 AU: an Analytical Treatment

We study the propagation, reflection, and turbulent dissipation of Alfven waves in coronal holes and the solar wind. We start with the Heinemann-Olbert equations, which describe non-compressive magnetohydrodynamic fluctuations in an inhomogeneous medium with a background flow parallel to the background magnetic field. Following the approach of Dmitruk et al, we model the nonlinear terms in these equations using a simple phenomenology for the cascade and dissipation of wave energy, and assume that there is much more energy in waves propagating away from the Sun than waves propagating towards the Sun. We then solve the equations analytically for waves with periods of hours and longer to obtain expressions for the wave amplitudes and turbulent heating rate as a function of heliocentric distance. We also develop a second approximate model that includes waves with periods of roughly one minute to one hour, which undergo less reflection than the longer-period waves, and compare our models to observations. Our models generalize the phenomenological model of Dmitruk et al by accounting for the solar wind velocity, so that the turbulent heating rate can be evaluated from the coronal base out past the Alfven critical point - that is, throughout the region in which most of the heating and acceleration occurs. The simple analytical expressions that we obtain can be used to incorporate Alfven-wave reflection and turbulent heating into fluid models of the solar wind.

💡 Research Summary

The paper presents an analytical treatment of Alfvén‑wave propagation, reflection, and turbulent dissipation in coronal holes and the expanding solar wind, covering the region from the solar surface (1 R⊙) out to 1 AU. Starting from the Heinemann‑Olbert equations, which describe non‑compressive MHD fluctuations (Elsässer variables z⁺ and z⁻) in an inhomogeneous medium with a background flow U parallel to the magnetic field B₀, the authors incorporate a phenomenological model for the nonlinear cascade originally introduced by Dmitruk et al. (2002). The key assumption is a strong imbalance between outward‑propagating waves (z⁺) and inward‑propagating (reflected) waves (z⁻), a condition supported by observations of coronal holes.

Two distinct wave populations are treated separately:

1. Long‑period waves (periods of several hours or more). Because their parallel wavenumber k∥ is small, they experience strong reflection from gradients in the Alfvén speed V_A(r). The reflection coefficient is approximated as
   R ≈ (d ln V_A/dr) / (2 k∥ V_A),
   leading to z⁻ ≈ R z⁺. The reflected component then participates in the nonlinear interaction, producing a turbulent heating rate
   ε = ρ (z⁺)² z⁻ / λ⊥,
   where λ⊥ is the transverse correlation length, assumed to expand linearly with heliocentric distance.

2. Intermediate‑to‑short‑period waves (≈1 min–1 h). Their larger k∥ reduces the linear reflection efficiency, so R is small. Nevertheless, the authors retain a nonlinear “self‑reflection” term: even in the absence of strong linear reflection, the nonlinear cascade can generate an inward component z⁻ from the dominant outward component z⁺. This yields a similar heating expression, but with z⁻ determined by the cascade rather than by linear reflection.

The background solar‑wind profile U(r) and Alfvén speed V_A(r) are taken from a Parker‑type solution, ensuring that the Alfvén critical point (U = V_A) is explicitly included. This allows the heating rate to be evaluated continuously from the coronal base, through the critical point, and out to 1 AU, where most of the solar‑wind acceleration and heating occur.

Analytical solutions for z⁺(r), z⁻(r), and ε(r) are derived for the two wave families. When only the long‑period component is retained, the model predicts a strong concentration of heating near the Sun (1–3 R⊙) but insufficient heating at larger distances, inconsistent with in‑situ measurements. Adding the intermediate‑to‑short‑period component spreads the heating over a broader range (0.3–1 AU) and brings the predicted electron temperature profile, plasma β, and turbulent spectra into good agreement with observations from Helios, Ulysses, and recent Parker Solar Probe data.

The paper’s principal contributions are:

• A closed‑form analytical solution of the Heinemann‑Olbert equations with a simple, physically motivated nonlinear term, valid from the solar surface to 1 AU.
• Explicit inclusion of the solar‑wind outflow U in the reflection and cascade physics, allowing the turbulent heating rate to be evaluated across the Alfvén critical point.
• A clear separation of wave populations by period, demonstrating how long‑period waves dominate reflection while shorter‑period waves sustain the cascade through self‑generated inward fluctuations.
• Practical formulae for ε(r) that can be inserted into fluid or reduced‑MHD solar‑wind models without the need for costly full‑scale numerical simulations.

Overall, the study provides a tractable yet realistic framework for incorporating Alfvén‑wave reflection and turbulent heating into global solar‑wind models, offering a bridge between analytic theory, numerical simulation, and spacecraft observations.