New insights from cross-correlation studies between solar activity indices and cosmic-ray flux during Forbush decrease events

This study utilizes cross-correlation analysis of data from the SOHO/ERNE instrument and neutron monitor networks across two solar cycles to demonstrate that power exponents modeling solar energetic particle fluence spectra are superior predictors of Forbush decrease magnitude compared to coronal mass ejection velocities.

The Gist

The Big Picture: The Sun's "Space Weather"

Imagine the Sun isn't just a ball of fire, but a giant, active sprinkler system. Sometimes, it sneezes. When it sneezes, it shoots out huge clouds of gas (Coronal Mass Ejections or CMEs) and blasts of fast particles (Solar Energetic Particles or SEPs).

When these solar sneezes hit Earth, two main things happen:

1. The "Ground Level Enhancement" (GLE): The fast particles hit our atmosphere and cause a brief, dangerous spike in radiation. Think of this as a sudden, blinding flash of light.
2. The "Forbush Decrease" (FD): After the flash, the solar wind acts like a giant shield or a vacuum cleaner, sweeping away the normal background radiation (Galactic Cosmic Rays) that usually drifts in from deep space. This causes a sudden drop in radiation levels. Think of it as a sudden fog rolling in that blocks out the stars.

What Did These Scientists Do?

The researchers wanted to understand the relationship between the sneezes (the CMEs), the flashes (the particle bursts), and the fog (the Forbush Decrease).

Specifically, they asked: Can we look at the shape of the "flash" (the particle data) to predict how thick the "fog" will be?

Usually, scientists try to predict how bad the "fog" (Forbush Decrease) will be by measuring how fast the "sneeze" (the CME) was moving. It's like trying to predict how big a wave will be by looking at how fast a boat is moving.

The New Discovery: The "Recipe" of the Particles

The team looked at data from two solar cycles (about 22 years of history). They didn't just look at speed; they looked at the energy spectrum of the particles.

The Analogy:
Imagine the solar particles are a bag of mixed candies.

• Old way: You just weigh the bag to see how heavy the sneeze was.
• New way: You look at the ratio of small candies to big candies. Is it mostly tiny sprinkles? Or mostly giant gumballs?

The scientists found that the "recipe" (the mathematical shape of the energy spectrum) of the particles tells a better story than just the speed of the CME. They used two numbers (called power exponents, a and b) to describe this recipe.

The "Magic Number" (6%)

Here is the most interesting part of their findings. They found that the relationship between the particle "recipe" and the "fog" (Forbush Decrease) changes depending on how big the event is.

They discovered a tipping point at 6%:

• Big Events (>6% drop): When the "fog" is very thick (a massive drop in radiation), the particle "recipe" is a perfect predictor. If you know the recipe, you can accurately guess how thick the fog will be in space, even better than knowing the speed of the CME.
• Small Events (<6% drop): When the "fog" is light, the recipe doesn't help much. In fact, the relationship gets messy or even backwards.

Why Does This Matter?

Think of Earth's magnetic field as a bouncer at a club.

• The CME is the storm outside.
• The Forbush Decrease is how many people get kicked out of the club.
• The Magnetic Field is the bouncer who decides who gets in or out based on the storm.

The scientists realized that the "particle recipe" is actually a better way to predict how many people get kicked out outside the club (in deep space) than just measuring the storm's speed. It helps them separate the effect of the storm itself from the effect of Earth's bouncer (our magnetic field).

The Takeaway

This paper suggests that instead of just asking "How fast was the solar wind?", we should ask "What did the particles look like?"

By analyzing the shape of the particle energy (the recipe), scientists can now predict the strength of space weather events more accurately, especially for the big, dangerous ones. It's like realizing that to predict a storm's damage, you shouldn't just look at the wind speed, but also at the type of rain falling.

In short: They found a new "secret code" in the solar particles that helps us predict space weather storms better than the old methods, especially when the storms are really big.

Technical Summary

1. Problem Statement

Forbush decreases (FDs) are transient reductions in galactic cosmic ray (CR) intensity caused by interplanetary coronal mass ejections (ICMEs) and their associated shockwaves. While the relationship between ICME parameters (specifically CME velocity) and FD magnitude is well-established, predicting the exact magnitude of CR modulation in the interplanetary magnetic field (IMF) remains challenging.

Current predictive models often rely on CME velocities, but these do not fully account for the complex interaction between the IMF, the geomagnetic field, and the specific spectral characteristics of Solar Energetic Particles (SEPs) that accompany these events. Furthermore, distinguishing the effects of the IMF from geomagnetic shielding at Earth is difficult. The authors aim to investigate whether the spectral shape of SEP fluence (specifically power-law exponents) offers a better predictive capability for FD magnitudes in the interplanetary medium than traditional CME velocity metrics.

2. Methodology

Data Sources:

• SEP Data: Proton flux data from the SOHO/ERNE instrument (Low-Energy and High-Energy Detectors) covering Solar Cycles 23 and 24.
• FD & Solar Parameters: Data from the IZMIRAN database, including FD magnitudes ($M$ and $M_M$), CME velocities ($V_{mean}$, $V_{meanC}$, $V_{max}$), and solar wind speeds.
• Geomagnetic Indices: $Kp_{max}$, $Ap_{max}$, and $Dst_{min}$.
• Ground Observations: Neutron monitor (NM) data to validate FD events.

Event Selection:

• The study focused on non-recurrent, ICME-induced FDs with magnitudes $>4\%$ for 10 GV rigidity particles.
• A strict temporal correlation was enforced to ensure that the SEP flux increase, the FD, and the associated CME were linked to the same solar event.
• Final Dataset: 21 high-quality events were selected after filtering out events with overlapping flux time series or high uncertainty in fluence integration.

Parametrization of SEP Spectra:

1. Fluence Calculation: Proton flux time series were integrated over specific time windows defined by the onset and recovery of the FD event.
2. Spectral Modeling: The energy spectra were fitted using a double power-law function with an exponential cutoff (Band function).
   • The model includes a "knee" energy ($E_b$) where the spectral slope changes.
   • Two power-law exponents were derived: $\alpha$ (slope below the knee) and $\beta$ (slope above the knee).
3. Knee Energy Determination: Instead of visual estimation, $E_b$ was calculated quantitatively by fitting a power function to the relationship between integral fluence and knee energy ($E_b = aJ^b$).

Correlative Analysis:

• Pearson correlation coefficients ($r$) were calculated between the spectral exponents ($\alpha, \beta$) and various parameters: FD magnitudes ($M$, $M_M$), CME velocities, and geomagnetic indices.
• $M$ represents the raw FD magnitude, while $M_M$ is the magnitude corrected for magnetospheric effects (using the $Dst$ index) to approximate the FD magnitude in the interplanetary magnetic field.

3. Key Contributions

• Novel Parameterization: The study introduces the power-law exponents ($\alpha$ and $\beta$) of SEP fluence spectra as new, quantitative parameters for analyzing solar-terrestrial connections.
• Decoupling IMF and Geomagnetic Effects: By utilizing the magnetospheric-corrected FD magnitude ($M_M$), the authors isolate the CR disturbance occurring in the interplanetary medium from local geomagnetic shielding effects.
• Discovery of Event Classes: The analysis reveals a bifurcation in event behavior based on the magnitude of the corrected FD ($M_M$), suggesting two distinct classes of solar events with different physical drivers or propagation characteristics.

4. Key Results

Correlation Findings:

• SEP Exponents vs. CME Velocity: Strong correlations were found between the spectral exponents ($\alpha, \beta$) and mean CME velocities ($r \approx 0.75 - 0.77$).
• SEP Exponents vs. FD Magnitude ($M$): The correlation with raw FD magnitude ($M$) was significant but lower than with CME velocities ($r \approx 0.67$).
• SEP Exponents vs. Corrected FD Magnitude ($M_M$): Crucially, the correlation between the spectral exponent $\alpha$ and the magnetospheric-corrected magnitude ($M_M$) was found to be stronger ($r \approx 0.64$) than the correlation between $M_M$ and CME velocities ($r \approx 0.57$).

The "Two-Class" Phenomenon:
When the dataset was split based on the corrected FD magnitude ($M_M$):

1. High Magnitude Group ($M_M \geq 6\%$): Exhibited very strong positive correlations between spectral exponents and FD magnitudes ($r \approx 0.77$ for $\alpha$ vs. $M_M$). The correlation with CME velocity was also high.
2. Low Magnitude Group ($M_M < 6\%$): Exhibited negative or near-zero correlations. For this group, spectral exponents and CME velocities showed little to no predictive power regarding the FD magnitude.

Predictive Capability:
The study concludes that SEP fluence power exponents are potentially better predictor variables for the magnitude of Forbush decreases in the interplanetary magnetic field than CME velocities, particularly for events where $M_M \geq 6\%$.

5. Significance

• Improved Space Weather Modeling: The findings suggest that the spectral shape of SEPs contains critical information about the efficiency of particle acceleration and the structure of the interplanetary shock, which directly influences how much CR flux is modulated.
• Refined Prediction Models: Traditional models relying solely on CME speed may be insufficient. Incorporating SEP spectral parameters could improve the accuracy of forecasting CR intensity drops, which is vital for radiation safety in aviation and space missions.
• Physical Insight: The distinction between high and low-magnitude events implies that different physical mechanisms or propagation conditions dominate depending on the severity of the solar eruption. The spectral exponents serve as a diagnostic tool to differentiate these regimes.

In summary, Savic et al. demonstrate that analyzing the spectral hardness (power exponents) of Solar Energetic Particles provides a more robust link to the interplanetary modulation of cosmic rays than traditional kinematic parameters like CME velocity, offering a new avenue for understanding and predicting Forbush decreases.