A turbulence index independent framework for deriving solar wind speed and coronal electron density from radio spectral broadening

This paper presents a generalized, turbulence-index-independent framework that utilizes spectral broadening of spacecraft radio signals to simultaneously derive solar wind speeds and coronal electron densities, demonstrating its effectiveness through Akatsuki observations while highlighting the necessity of accounting for turbulence anisotropy in slow solar wind regions.

The Gist

Imagine the Sun isn't just a glowing ball of fire, but a giant, chaotic fountain. It constantly sprays a stream of invisible particles (plasma) into space. This stream is called the solar wind. Sometimes it flows gently like a slow river (the "slow wind"), and sometimes it rushes out like a high-speed jet (the "fast wind").

For a long time, scientists have wanted to know exactly how fast this wind is blowing and how dense it is right near the Sun. But you can't just stick a thermometer or a speedometer there; the heat would melt any instrument instantly.

This paper presents a clever new way to "feel" the wind without touching it, using a method that doesn't require guessing the rules of the game first.

The "Radio Flashlight" Analogy

Imagine you are standing on Earth holding a very powerful, steady flashlight (the Akatsuki spacecraft orbiting Venus). You shine this light across the solar system, aiming it so the beam just grazes the edge of the Sun before hitting a receiver back on Earth.

As the light beam passes through the Sun's atmosphere (the corona), it doesn't travel through a vacuum. It has to pass through a turbulent soup of invisible particles. Think of this soup like a foggy day where the air is churning with invisible whirlpools and eddies.

• The Problem: When your flashlight beam hits these invisible whirlpools, the light gets "jittery." It wobbles, spreads out, and gets a little fuzzy. In radio terms, this is called spectral broadening.
• The Old Way: Previously, scientists tried to calculate the wind speed and density by assuming the "whirlpools" in the solar wind were all the same size and shape (a specific mathematical rule called the "Kolmogorov spectrum"). It was like assuming every storm in the world has the exact same wind pattern. If the storm was actually different, their calculations would be slightly off.
• The New Way (This Paper): The authors in this paper said, "Let's stop guessing the shape of the whirlpools." Instead, they built a new mathematical framework that can measure the shape of the turbulence itself directly from the radio signal. It's like having a camera that doesn't just measure how blurry the photo is, but also figures out why it's blurry (is it wind? is it rain? is it a shaky hand?) and then corrects the calculation based on that specific reason.

How They Did It

The team looked at data from two different times when the Akatsuki spacecraft passed behind the Sun from Earth's perspective:

1. 2016 (The Slow Wind): They looked at a time when the solar wind was slow and coming from a region with giant magnetic loops (like a river flowing through a canyon).
2. 2022 (The Fast Wind): They looked at a time when the wind was fast and coming from a "hole" in the Sun's magnetic field (like a jet stream).

By analyzing the "fuzziness" of the radio signal in both cases, they could:

• Measure the Turbulence: They figured out the "texture" of the solar wind's turbulence (whether it was smooth, choppy, or chaotic).
• Calculate the Speed and Density: Once they knew the texture, they could accurately calculate how fast the wind was blowing and how many particles were in it, without needing to assume a fixed rule.

What They Found

• The Slow Wind (2016): Near the Sun, the turbulence was messy and "flat" (like a calm lake with big, slow ripples). As the wind moved further away, it became more chaotic and "steep" (like white-water rapids), eventually settling into a standard pattern. This showed that the slow wind takes a while to "get its act together" and become a smooth, fast stream.
• The Fast Wind (2022): This was surprising! Even very close to the Sun, the fast wind was already behaving like a mature, chaotic storm. It was "isotropic," meaning the turbulence was the same in all directions, like a perfectly mixed smoothie. This suggests the fast wind gets its speed boost very early, right near the Sun's surface.

Why This Matters

Think of this like weather forecasting. Before, we had to assume all storms were the same to predict the wind. Now, we have a tool that can look at a specific storm, understand its unique personality, and give us a much more accurate forecast.

This new "turbulence-independent" framework allows scientists to:

1. Get better data: They can now measure solar wind speed and density more accurately, even in weird or changing conditions.
2. Understand the Sun better: It helps explain how the Sun heats its outer atmosphere and accelerates the wind to supersonic speeds.
3. Protect Earth: Since the solar wind drives "space weather" (which can knock out satellites and power grids), understanding it better helps us predict and prepare for solar storms.

In short, the authors turned a radio signal into a turbulence microscope, allowing us to see the invisible churning of the solar wind with unprecedented clarity, without needing to guess the rules of the game first.

Technical Summary

1. Problem Statement

The solar wind originates from the solar corona, where plasma is heated and accelerated by turbulent energy transfer. Understanding the physical mechanisms driving this acceleration requires precise measurements of electron density ($N_e$) and solar wind velocity ($v$) in the near-Sun region ($1.4 - 10 R_\odot$).

Traditional radio occultation (RO) techniques rely on analyzing the spectral broadening of spacecraft signals caused by plasma turbulence. However, previous methods typically assumed a fixed turbulence spectral index ($p = 11/3$), corresponding to a fully developed, isotropic Kolmogorov spectrum.

• The Limitation: Observations suggest that the turbulence spectrum is not always isotropic or Kolmogorov-like. It varies with heliocentric distance, latitude, and solar wind type (fast vs. slow). Assuming a fixed $p$ introduces systematic errors in derived plasma parameters, particularly in regions where turbulence is anisotropic or in transition.
• The Goal: Develop a framework that derives $N_e$ and $v$ without assuming a fixed turbulence index, allowing for the direct use of measured spectral slopes to account for varying turbulence regimes.

2. Methodology

The authors utilized X-band (8.41 GHz) radio occultation data from the JAXA Akatsuki spacecraft during superior conjunctions with Earth in June 2016 (Solar Cycle 24, declining phase) and October 2022 (Solar Cycle 25, rising phase).

Data Processing

1. Signal Acquisition: Open-loop radio science data were recorded at the Usuda Deep Space Center.
2. Detrending: Raw Doppler frequency time series were detrended by removing deterministic effects (spacecraft motion, Earth rotation) and slow instrumental drifts using a second-order polynomial fit.
3. Power Spectral Density (PSD): The detrended residuals were analyzed using Welch's method with overlapping sliding windows to compute the PSD.
4. Spectral Index Extraction: The temporal spectral slope ($\alpha$) was determined by fitting a power law ($P(f) \propto f^{-\alpha}$) to the PSD in the inertial range ($10^{-3}$ to $10^{-1}$ Hz). The spatial turbulence index was then calculated as $p = \alpha + 3$.

Theoretical Framework

The study generalized existing empirical formulas (previously derived for S-band and fixed $p=11/3$) to accommodate arbitrary spectral indices ($p$) and observing wavelengths ($\lambda$).

• Scaling Relation: Spectral broadening ($B$) scales with wavelength as $B \propto \lambda^{2/(p-2)}$.
• Generalized Equations: By substituting the scaling relation into the standard empirical equations, the authors derived new expressions for $N_e$ and $v$ that explicitly include the measured spectral index $p$ (or $\alpha$):
  $$N_e \propto B^{5/6} \cdot \lambda^{-\frac{5}{3(p-2)}}$$
  $$v \propto B^{1/6} \cdot \lambda^{-\frac{1}{3(p-2)}}$$
  This allows the calculation of plasma parameters using the actual measured turbulence spectrum rather than an assumed one.

3. Key Contributions

• Turbulence-Independent Formulation: The primary contribution is a mathematical framework that decouples the derivation of solar wind parameters from the assumption of a fixed Kolmogorov spectrum. It enables the use of observed spectral slopes to correct for anisotropy and varying turbulence regimes.
• Unified Multi-Frequency Approach: The method provides a consistent way to interpret data across different radio frequencies (S-band, X-band, etc.) by incorporating wavelength-dependent scaling factors derived from the turbulence index.
• Comparative Analysis of Wind Regimes: The study applies this framework to two distinct solar wind environments:
  • 2016 (Slow Wind): Streamer belt regions.
  • 2022 (Fast Wind): Equatorial coronal hole regions.

4. Results

Turbulence Characteristics

• 2016 (Slow Wind): The spectral index $p$ showed significant radial evolution. Near the Sun ($1.45 R_\odot$), $p \approx 3.18$ (flattened spectrum, indicating energy injection from large-scale structures). Beyond $3 R_\odot$, $p$ steepened and stabilized around $3.6-3.7$, approaching the Kolmogorov value ($11/3 \approx 3.67$). This indicates a transition from energy-injection-dominated turbulence to a developed MHD cascade.
• 2022 (Fast Wind): The spectral index remained relatively constant and close to the Kolmogorov value ($p \approx 3.46 - 3.79$) across all sampled distances ($3.9 - 9.1 R_\odot$). This suggests that turbulence in coronal hole outflows is already well-developed and nearly isotropic close to the solar surface.

Derived Plasma Parameters

• Electron Density ($N_e$):
  • 2016: Decreased from $\sim 8.7 \times 10^{12} \, \text{m}^{-3}$ at $1.45 R_\odot$ to $\sim 7.7 \times 10^{10} \, \text{m}^{-3}$ at $8.85 R_\odot$.
  • 2022: Ranged from $\sim 2.0 \times 10^{11} \, \text{m}^{-3}$ at $3.9 R_\odot$ to $\sim 2.6 \times 10^{10} \, \text{m}^{-3}$ at $9.1 R_\odot$.
  • These trends align well with empirical coronal models (e.g., Guhathakurta & Holzer, 1994) and in-situ Parker Solar Probe data.
• Solar Wind Velocity ($v$):
  • 2016: Accelerated from $\sim 32 \, \text{km/s}$ to $\sim 155 \, \text{km/s}$.
  • 2022: Accelerated from $\sim 171 \, \text{km/s}$ to $\sim 363 \, \text{km/s}$.
  • The fast wind velocities in 2022 were significantly higher, consistent with the characteristics of coronal hole outflows.

Validation

The results derived using the variable turbulence index showed differences of less than 5% compared to previous studies that assumed a fixed Kolmogorov spectrum. However, the new framework revealed systematic deviations in the slow wind (2016) that were attributed to anisotropy, which the fixed-index model could not capture.

5. Significance

• Refined Diagnostics: The study demonstrates that explicitly considering near-coronal turbulence anisotropy is essential for accurate retrieval of solar wind parameters. Ignoring the variation in $p$ can lead to systematic biases, particularly in slow wind regions.
• Physical Insight: The contrasting behavior of $p$ in 2016 vs. 2022 provides observational evidence that fast solar wind turbulence is mature and isotropic near the Sun, whereas slow wind turbulence undergoes a distinct evolution from large-scale injection to a developed cascade.
• Future Applications: This methodology is scalable for future multi-frequency missions (e.g., Square Kilometre Array, Parker Solar Probe, Solar Orbiter). It bridges the gap between remote sensing (radio occultation) and in-situ measurements, offering a robust tool for mapping energy transport and wave-particle interactions in the inner heliosphere.

In conclusion, this paper establishes a robust, turbulence-index-independent framework that enhances the accuracy of remote sensing for solar wind parameters, revealing critical nuances in the nature of coronal turbulence across different solar wind regimes.