Radial Dependency of ICME-associated Particle Acceleration Processes: Statistical Multipoint Observations from 2016-2023

This study statistically analyzes 39 multipoint ICME events observed between 2016 and 2023 to reveal that shock acceleration efficiency for energetic particles consistently increases with heliocentric distance up to 0.7 au before declining at greater distances.

The Gist

Imagine the Sun as a giant, chaotic lighthouse that occasionally blasts out massive clouds of charged gas and magnetic fields. These are called Coronal Mass Ejections (CMEs). As these clouds travel through space, they act like a snowplow, pushing the solar wind ahead of them and creating a massive, invisible shockwave at the front.

When this shockwave hits, it acts like a cosmic particle accelerator, slamming into tiny particles (like protons and helium nuclei) and boosting them to incredibly high speeds. These high-speed particles are called Energetic Storm Particles (ESPs).

This paper is a statistical detective story. The authors wanted to answer a simple question: Does the "speed" of this particle accelerator change as the shockwave travels further away from the Sun?

The Setup: A Cosmic Relay Race

To solve this, the researchers didn't just look at one spot. They used a "distributed array" of spacecraft, which is like having a relay team of observers stationed at different distances from the Sun:

• Parker Solar Probe: The sprinter, closest to the Sun (as close as 0.045 AU).
• Solar Orbiter: The middle-distance runner (around 0.3 AU).
• STEREO-A, Wind, and ACE: The long-distance runners, sitting near Earth's orbit (1 AU).

Between 2016 and 2023, they tracked 39 specific events where these different spacecraft all saw the same shockwave pass by. They filtered this down to 23 events where the spacecraft were lined up well enough to compare notes.

The Investigation: Measuring the "Break"

When these particles are accelerated, their energy levels don't just go up in a straight line. If you graph their energy, the line usually goes up, hits a specific point, and then changes slope. The authors call this the "spectral break."

Think of the spectral break like a speed limit sign on a highway.

• Below the sign, cars (particles) are accelerating easily.
• At the sign, the rules change, and it becomes much harder to go faster.
• The higher the "speed limit" (the energy of the break), the more efficient the accelerator is at pushing particles to extreme speeds.

The researchers used complex math to find the exact location of this "speed limit" for different types of particles (mostly Helium-4) at different distances from the Sun.

The Surprise Discovery: The "Sweet Spot"

The team expected to see a simple story: as the shockwave moves away from the Sun, it gets weaker (like a sound fading as you walk away from a speaker). They expected the "speed limit" to drop steadily the further out you go.

But the data told a different story.

1. The Inner Loop (0 to 0.7 AU): As the shockwave traveled from the Sun out to about 70% of the distance to Earth, the "speed limit" actually went up. The accelerator became more efficient the further it traveled.

   • The Analogy: Imagine a runner starting a race. Instead of getting tired immediately, they find a "sweet spot" in the middle of the track where the wind is perfectly at their back, and they suddenly start running faster than they did at the starting line.
   • The Cause: The authors suggest this is due to particle trapping. As the shock moves, it creates a turbulent "foreshock" region (like a wake behind a boat). This region acts like a cage, trapping particles and giving them more time to bounce back and forth, gaining more energy before they escape.

2. The Outer Loop (Beyond 0.7 AU): Once the shockwave passed the 0.7 AU mark and headed toward Earth, the "speed limit" finally started to drop, just as the team originally expected.

   • The Analogy: The runner finally hits the headwind. The magnetic field weakens, the shock slows down, and the "cage" becomes less effective. Particles start to escape, and the maximum energy they can reach drops.

What They Didn't Find

The researchers also checked if the angle of the shockwave or the turbulence of the magnetic field was the main reason for these changes.

• They found that the angle of the shock (whether it was hitting head-on or glancing) didn't seem to be the main driver.
• They found that the "bounciness" of the magnetic field (turbulence) didn't have a simple, direct correlation with the energy changes in this specific dataset.

The Bottom Line

The paper concludes that the efficiency of the Sun's particle accelerator isn't a straight line. It has a peak performance zone between the Sun and about 70% of the way to Earth.

• Near the Sun: The accelerator is just getting warmed up.
• Middle Distance (0.2 – 0.7 AU): The accelerator hits its stride, trapping particles and boosting them to their highest energies.
• Far Distance (Near Earth): The accelerator starts to wind down as the shockwave weakens.

This finding is crucial because it changes how we predict space weather. If we want to know how dangerous a solar storm will be for satellites or astronauts near Earth, we can't just look at how strong the storm was when it left the Sun. We have to understand how the shockwave evolves and "traps" particles during its journey through the inner solar system.

Technical Summary

1. Problem Statement

Interplanetary Coronal Mass Ejections (ICMEs) drive shocks that accelerate Solar Energetic Particles (SEPs) and Energetic Storm Particles (ESPs). While the Diffusive Shock Acceleration (DSA) mechanism is well-established, the radial evolution of these acceleration processes remains poorly understood.

• The Gap: Theoretical models and numerical simulations suggest that as an ICME-driven shock propagates away from the Sun, the weakening magnetic field and shock speed should lead to a decrease in the maximum energy particles can reach (a decrease in spectral break energy, $E_0$).
• The Challenge: Definitively observing this radial evolution requires simultaneous, multipoint in-situ measurements at varying heliocentric distances. Historically, such data was scarce. With the advent of missions like Parker Solar Probe (PSP) and Solar Orbiter, combined with near-Earth assets (Wind, ACE, STEREO-A), a statistical analysis of radial dependencies is now possible.
• Objective: To statistically analyze how local shock conditions and particle acceleration efficiency evolve as an ICME propagates from the inner heliosphere ($<0.7$ au) to 1 au, specifically focusing on the spectral shape of upstream ESPs.

2. Methodology

Data Collection and Event Selection

• Database: The authors compiled a catalog of 39 multipoint ICME events occurring between 2016 and 2023.
• Spacecraft: Observations were drawn from Parker Solar Probe (PSP), Solar Orbiter, ACE, Wind, and STEREO-A.
• Selection Criteria:
  • Events must be observed by at least two spacecraft.
  • Radial Separation ($\Delta r$): Must be $> 0.2$ au to ensure significant radial evolution.
  • Longitudinal Separation ($\Delta \text{Long}$): Must be $< 45^\circ$ to ensure the spacecraft are observing the same shock structure rather than different parts of a wide CME.
  • Final Sample: After filtering for data availability and separation constraints, 23 events were selected for detailed analysis (yielding 52 individual sightings).

Spectral Analysis

• Target Population: The study focuses on upstream ESPs (particles trapped ahead of the shock), as they are less affected by downstream turbulence and better reflect the acceleration mechanism.
• Ion Species: Analysis included H, $^4$He, C, O, and Fe. However, due to data limitations and the need for consistent multipoint comparisons, $^4$He was the primary species used for multipoint statistical analysis.
• Fitting Models: Particle spectra were fitted using three models to determine the best representation:
  1. Pan-Spectrum Fitting Formula: A flexible model incorporating a spectral break region.
  2. Classic Double Power Law (CDPL): A sharp break at energy $E_0$.
  3. Single Power Law: Used when no break was observed within the instrument's energy range.
• Parameter Extraction: The study extracted the spectral break energy ($E_0$) and power-law indices ($\beta_1, \beta_2$). Local shock parameters (shock normal angle $\theta_{bn}$, magnetic turbulence $\delta B/B_0$, density compression $n_{comp}$, shock speed $V_{sh}$) were also calculated.

Statistical Approach

• Correlation Analysis: Used Kendall's $\tau$ coefficient to assess correlations between spectral parameters and shock conditions/distance.
• Relative Change Analysis: For multipoint events, the authors calculated the relative change in break energy ($dE_0$) between two observers to isolate radial evolution effects from event-to-event variability.

3. Key Results

A. Correlation with Local Shock Parameters

• Local Parameters: No statistically significant correlation was found between the spectral break energy ($E_0$) and local shock parameters (shock angle, magnetic turbulence, or density compression) across all species. This suggests that local shock conditions alone cannot predict the spectral break energy.
• Radial Distance: A significant positive correlation was found between $E_0$ and heliocentric distance for observations within 0.8 au.
  • When filtering for $r < 0.8$ au, the correlation became statistically significant for all ion species ($\tau = 0.338, p = 0.005$) and specifically for $^4$He ($\tau = 0.556, p = 0.045$).

B. Radial Evolution of Spectral Break Energy

The study revealed a non-monotonic evolution of the spectral break energy:

1. Inner Heliosphere (0.1 – 0.7 au): Contrary to simple theoretical expectations, the spectral break energy increases with distance. This implies that particle acceleration efficiency is increasing as the shock moves away from the Sun in this region.
2. Outer Region (> 0.7 au): Beyond approximately 0.7 au, the trend reverses, and the break energy decreases as the shock approaches 1 au. This aligns with traditional expectations of weakening shock strength and magnetic fields.

C. Spectral Indices

• Low-energy index ($\beta_1$): Showed no clear statistical variation with radial distance.
• High-energy index ($\beta_2$): Systematically decreased with distance (spectra became harder at high energies) up to ~1 au, though with high variability.

4. Discussion and Physical Interpretation

The authors propose that the initial increase in acceleration efficiency (and thus $E_0$) in the inner heliosphere is driven by particle trapping in the foreshock region.

• Mechanism: As the ICME propagates, it induces upstream turbulence. This creates a foreshock region that effectively traps high-energy particles, allowing them to undergo continuous acceleration over a longer timescale as they transit through the inner heliosphere.
• Transition Point: This trapping mechanism remains efficient until the shock strength and magnetic field intensity decline sufficiently (around 0.7–0.8 au). Beyond this point, the efficiency of trapping drops, and particle escape dominates, leading to the expected decrease in maximum energy.
• Relation to DME/IVA Events: This trend may explain "Delayed Maximum Energy" (DME) or "Inverse Velocity Arrival" (IVA) events, where high-energy particles arrive later than expected because they are being accelerated further out in the heliosphere rather than near the Sun.

5. Significance and Contributions

• Statistical Validation: This is one of the first statistical studies to use multipoint in-situ data to validate the radial evolution of ICME-driven shock acceleration, moving beyond single-event case studies.
• Challenging Simple Models: The results challenge the simple assumption that shock acceleration efficiency strictly decreases with distance. The finding of an increase in efficiency within 0.7 au highlights the complexity of particle trapping and foreshock dynamics.
• Space Weather Implications: Understanding that high-energy particle production can peak at 0.6–0.8 au rather than near the Sun has critical implications for space weather forecasting. It suggests that predictions based solely on near-Sun shock parameters may underestimate the intensity and energy of ESPs at Earth.
• Future Directions: The study highlights the need for more multipoint events with robust plasma data to further constrain the relationship between local turbulence and acceleration efficiency.

Conclusion

The paper concludes that ICME-driven shock acceleration efficiency is not a simple function of radial distance. Instead, it exhibits a "hump" profile: efficiency increases from the Sun to ~0.7 au due to foreshock particle trapping, followed by a decrease as the shock weakens further out. This finding necessitates a revision of current SEP/ESP prediction models to account for radial evolution and trapping effects.