🔭 astrophysics

Turbulence and its Potential Impact on Solar Chromospheric and Coronal Heating

Using particle-in-cell simulations and a transport model, this study demonstrates that low-frequency turbulence generated by mixed-polarity magnetic fields in the solar chromosphere can be efficiently transported and dissipated to provide sufficient heating for both the chromosphere and the corona, while also gradually heating spicules.

Original authors: Gary P. Zank, Xiaocan Li, Krishna Khanal, Alphonse C. Sterling, Masaru Nakanotani, Linging Zhao, Laxman Adhikari, Yalim Mehmet, Subramania Athiray Panchapakesan, Fan Guo, Ronald L. Moore

Published 2026-04-30

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Gary P. Zank, Xiaocan Li, Krishna Khanal, Alphonse C. Sterling, Masaru Nakanotani, Linging Zhao, Laxman Adhikari, Yalim Mehmet, Subramania Athiray Panchapakesan, Fan Guo, Ronald L. Moore

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Sun's surface not as a calm, glowing ball, but as a chaotic, bubbling pot of soup. Just below the visible surface (the photosphere), there is a layer called the chromosphere, and above that, the super-hot corona. For decades, scientists have been puzzled by a big question: How does the corona get so incredibly hot, and how does the chromosphere stay warm enough? It's like trying to figure out why the steam above a pot is hotter than the water inside.

This paper proposes that the answer lies in turbulence—specifically, a kind of magnetic chaos generated by the Sun's own "magnetic carpet."

Here is a simple breakdown of what the researchers did and what they found:

1. The Magnetic Carpet: A Tangled Mess

Think of the Sun's surface as being covered in a "carpet" made of tiny, invisible magnetic loops. These loops are constantly popping up, twisting, and disappearing. Some loops have their magnetic north pole pointing up, while others point down.

The Problem: When a "north" loop bumps into a "south" loop, they don't just sit there. They snap and reconnect, like two tangled rubber bands suddenly snapping apart.
The Result: This snapping (called magnetic reconnection) doesn't just release a little spark; it creates a massive, chaotic storm of magnetic turbulence.

2. The Simulation: A Digital Sandbox

To understand this, the scientists used a super-powerful computer simulation (called a "Particle-in-Cell" or PIC code). Imagine they built a tiny, digital box representing a piece of the Sun's atmosphere.

They filled this box with magnetic loops of mixed directions (some up, some down).
They watched what happened when these loops crashed into each other.

What they found:

The Chaos Takes Over: The original magnetic loops didn't just disappear; they were completely destroyed and replaced by a swirling, turbulent soup of tiny magnetic knots and sheets.
Two Types of Turbulence:
- In "Quiet" Areas (like the Quiet Sun): The magnetic field is a jumbled mess with no clear direction. The turbulence is isotropic, meaning it's chaotic in every direction, like a swarm of bees flying randomly.
- In "Open" Areas (like Coronal Holes): There is a stronger, dominant magnetic field pointing up. Here, the turbulence is anisotropic. It's like a crowd of people trying to walk through a narrow hallway; the movement is organized along the hallway but chaotic side-to-side.

3. The Transport Problem: How does the heat get up there?

Here is the tricky part: The turbulence is created low down in the chromosphere. But the corona is high up. How does the energy get from the bottom to the top?

No Big Wind: Unlike the solar wind that blows out into space, the chromosphere doesn't have a steady, strong upward wind to carry the heat.
The "Patchwork" Solution: The researchers realized the chromosphere is full of smaller, random "elevator rides."
- Spicules: These are like giant, fast-moving fountains of gas shooting up from the surface.
- Shockwaves: Like sound waves traveling through air, these are pressure waves moving up from the surface.
- Emerging Loops: As the magnetic carpet loops pop up, they drag gas with them.

The paper suggests these random, fast-moving flows act like trucks on a highway. They pick up the turbulent magnetic energy created at the bottom and carry it upward.

4. The Heating Mechanism: Friction in the Sky

Once this turbulence is carried up by the "trucks" (spicules and shocks), it doesn't just sit there. It starts to break down.

Imagine a large, swirling whirlpool. As it moves, it breaks into smaller and smaller whirlpools.
Eventually, these tiny whirlpools get so small that they create friction. In the Sun's atmosphere, this friction turns the magnetic energy into heat.
The Result: This process heats the chromosphere from the inside out. Any turbulence that isn't used up heating the lower layers gets carried all the way to the base of the corona, providing enough energy to keep the Sun's outer atmosphere scorching hot.

5. The "Spicule" Mystery

The paper also offers a new explanation for Type II spicules (those fast, hot fountains mentioned earlier).

Scientists have noticed that these spicules get hotter as they go higher, which is strange because usually things cool down as they rise.
The authors suggest that as these fast-moving spicules shoot up, they are "entraining" (dragging along) the magnetic turbulence. As the turbulence inside the spicule breaks down and dissipates, it heats the spicule itself.
The Analogy: Think of a car driving through a thick fog. The faster the car goes, the more the fog swirls around it, creating friction and heat. The spicule is the car, and the magnetic turbulence is the fog. The faster the spicule moves, the higher up it carries the turbulence before it finally burns off as heat.

The Bottom Line

The paper concludes that the Sun doesn't need a single, giant heater to warm its atmosphere. Instead, it relies on a constant, chaotic dance:

Magnetic loops constantly reconnect and create turbulence.
Random, fast-moving flows (spicules and shocks) act as elevators, carrying this turbulence upward.
As the turbulence breaks down into smaller and smaller pieces, it turns into heat.

This process generates enough energy to explain both the warmth of the chromosphere and the extreme heat of the corona, solving a puzzle that has stumped scientists for years. The authors emphasize that this is a model of energy transport and dissipation, not just a simple injection of hot plasma. The turbulence itself is the fuel, and the random flows are the delivery system.

1. Problem Statement

The solar chromospheric and coronal heating problems remain two of the most significant unsolved issues in solar physics. While energy requirements are well-constrained (approx. $4 \times 10^3$ J m $^{-2}$ s $^{-1}$ for the chromosphere and $1.4 \times 10^4$ J m $^{-2}$ s $^{-1}$ for the corona), the mechanisms responsible for dissipating this energy are debated.

The Gap: Traditional models often treat the chromosphere as a static atmosphere or rely solely on wave heating (Alfvén waves) or impulsive reconnection (nanoflares). However, the chromosphere is highly dynamic, filled with mixed-polarity magnetic fields ("magnetic carpet"), spicules, and shocks, yet lacks a steady, large-scale bulk flow to transport energy.
The Question: How is low-frequency turbulence generated in the chromosphere via magnetic reconnection, and how is this turbulence transported and dissipated throughout the chromosphere and into the low corona to provide the necessary heating?

2. Methodology

The authors employ a two-pronged approach combining kinetic simulations with statistical transport modeling.

A. Particle-in-Cell (PIC) Simulations

To accurately capture the physics of magnetic reconnection and the subsequent transition to turbulence (which Magnetohydrodynamics/MHD often fails to resolve at kinetic scales), the authors use the 3D VPIC code.

Setup: They simulate a force-free current sheet representing the emergence of mixed-polarity magnetic carpet loops from the photosphere into the chromosphere.
Cases: Two distinct magnetic environments are modeled:
1. Imbalanced Guide Field ( $B_g/B_0 = 1$ ): Surrogate for open coronal holes or large loops where a strong guide field exists.
2. Balanced Guide Field ( $B_g/B_0 = 0.2$ ): Surrogate for the quiet Sun where mixed polarity dominates and the guide field is weak.
Physics: The simulations track the evolution from initial reconnection (tearing mode instability) to fully developed turbulence, analyzing energy spectra, magnetic islands, flux ropes, and current sheets.

B. Statistical Transport Model

Using the simulation results, the authors develop a 1D transport model for turbulence advected through the chromosphere.

Turbulence Type: The model focuses on quasi-2D non-propagating structures (magnetic islands/flux ropes) and a minority slab component (Alfvén waves).
Advection Mechanism: Since no steady bulk wind exists in the chromosphere, turbulence is transported by a "patchwork" of random, energy-containing flows: magnetic carpet emergence, acoustic shocks, and Type I/II spicules.
Statistical Approach: Flow speeds ( $U$ ) are modeled using log-normal probability distribution functions (PDFs). The authors solve coupled transport equations for Elsässer energy density and correlation length, incorporating Kolmogorov-like dissipation laws.
Parameters: The model evaluates expectations for total energy density $\langle y \rangle$ , Elsässer specific energy $\langle \langle Z_\infty^2 \rangle \rangle$ , heating rates $\langle \dot{H} \rangle$ , and correlation lengths $\langle \lambda \rangle$ as functions of height ( $h$ ).

3. Key Contributions and Results

A. Nature of Reconnection-Driven Turbulence (PIC Results)

Rapid Transition: Mixed-polarity magnetic fields transition rapidly to a fully turbulent state dominated by small-scale nonlinear structures. The initial magnetic field is largely annihilated.
Anisotropy vs. Isotropy:
- Imbalanced ( $B_g/B_0=1$ ): Turbulence is anisotropic. It is dominated by quasi-2D magnetic fluctuations (flux ropes) perpendicular to the mean field, with a power spectrum $\propto k_\perp^{-3/2}$ . Slab turbulence (Alfvén waves) is a minority component ( $\propto k_\parallel^{-2}$ ).
- Balanced ( $B_g/B_0=0.2$ ): Turbulence is isotropic. Without a strong guide field, the fluctuations are randomly oriented non-propagating structures (flux ropes) with a Kolmogorov-like spectrum ( $\propto k^{-5/3}$ ).
Energy Dominance: In both cases, the fluctuating magnetic and kinetic energy densities significantly exceed the mean field energy. The quasi-2D component is energetically dominant over the slab component.

B. Transport and Heating Mechanisms

Advection by Patchwork Flows: The model demonstrates that turbulence generated at the photosphere (via magnetic carpet reconnection) and entrained photospheric turbulence can be transported upward by stochastic flows (spicules, shocks, emerging loops).
Heating Rates: The calculated turbulent heating rates $\langle \dot{H} \rangle$ and energy injection rates $\langle \dot{S} \rangle$ at the base of the corona exceed the required energy constraints ( $4 \times 10^3$ to $1.4 \times 10^4$ J m $^{-2}$ s $^{-1}$ ) for various values of the Kolmogorov parameter $\alpha_K$ .
Role of Spicules:
- Slow flows (e.g., magnetic carpet emergence, Type I spicules) dissipate turbulence efficiently within the lower chromosphere, heating the plasma there.
- Fast flows (e.g., Type II spicules, $U > 30-50$ km/s) advect turbulence to greater heights before significant dissipation occurs. This explains the observed gradual heating of Type II spicules as they rise through the transition region into the corona.
Alfvén Wave Reflection: The model suggests that Alfvén waves generated by interchange reconnection are strongly reflected (>95%) at the transition region. This reflection creates counter-propagating waves that interact to generate 2D zero-frequency modes, further enhancing the quasi-2D turbulence population in the chromosphere.

C. Correlation Length Evolution

The model predicts that the correlation length $\langle \lambda \rangle$ increases with height. For specific parameters (e.g., $\alpha_K = 0.1$ ), the predicted correlation length at the transition region ( $\sim 3$ Mm) aligns well with observational data from recent studies (e.g., Z. Bailey et al., 2025).

4. Significance

Unification of Heating Mechanisms: The paper bridges the gap between reconnection and turbulence. It posits that reconnection of the magnetic carpet generates the turbulence, which is then dissipated via a turbulent cascade to heat the plasma. This moves beyond the "nanoflare" (direct heating) or "wave" (pure wave dissipation) dichotomies.
Chromospheric Dynamics: It provides a physical mechanism for how a static-density atmosphere can be heated uniformly despite the lack of a steady bulk wind: the stochastic advection of turbulence by transient flows (spicules, shocks).
Spicule Heating: It offers a robust explanation for the heating of Type II spicules, attributing it to the delayed dissipation of entrained turbulence in high-speed flows.
Coronal Base Conditions: The study establishes that the turbulent energy injection from the chromosphere is sufficient to drive the solar wind and heat the corona, providing specific boundary conditions (energy density, correlation length) for future solar wind models.

In conclusion, the authors demonstrate that low-frequency, quasi-2D turbulence, generated by the reconnection of the ubiquitous magnetic carpet and transported by the dynamic, patchwork flows of the chromosphere, is a viable and sufficient mechanism to satisfy the energy requirements for heating both the solar chromosphere and the corona.