A Remarkably Accurate Predictor of Sunspot Cycle Amplitude

This paper demonstrates that the linear relationship between sunspot and white-light facular areas at the onset of solar cycles serves as a highly accurate, physically based predictor of cycle amplitude up to four years in advance, successfully confirming historical data and correctly forecasting a larger-than-expected Cycle 25.

The Gist

Imagine the Sun as a giant, churning pot of magnetic soup. Every 11 years or so, this soup boils over in a predictable pattern called a Sunspot Cycle. During these cycles, the Sun gets "angry" (more active) or "calm" (less active). Scientists have been trying to guess how "angry" the next cycle will be for over a century because it affects our satellites, power grids, and even our auroras.

This paper by Peter Foukal presents a remarkably accurate "crystal ball" for predicting the strength of these solar storms. Here is the story of how he did it, explained simply.

The Old Clue: The "Spot-to-Freckle" Ratio

Back in the 1980s, scientists noticed a strange pattern in the Sun's history. They found that if you look at the very beginning of a new solar cycle, there is a relationship between two things:

1. Sunspots: The big, dark, angry scars on the Sun's surface.
2. Faculae: The tiny, bright "freckles" or glowing patches that surround the spots.

The Analogy: Think of the Sun's magnetic field like a garden.

• Sunspots are the big, dark weeds.
• Faculae are the tiny, bright flowers growing around them.

The researchers discovered that the ratio of weeds to flowers at the start of the season tells you how big the garden will get by the end of the season.

• If you see a lot of weeds compared to the flowers (a steep slope), the garden is going to be huge and wild (a strong solar cycle).
• If you see few weeds compared to the flowers (a shallow slope), the garden will stay small and tame (a weak solar cycle).

The Problem: The Garden Photos Were Blurry

For a long time, scientists used old photos from the 1800s and 1900s to count these weeds and flowers. But those photos were like old, grainy snapshots. It was hard to tell the difference between a tiny flower and a smudge on the lens. Because of this, they stopped using this method to predict future cycles.

The Solution: High-Definition Space Cameras

The author, Peter Foukal, decided to try again using modern technology. He used two powerful space telescopes:

• SOHO (MDI): A camera that took pictures from 2001 to 2011.
• SDO (HMI): A super-sharp camera that started in 2011 and is still working today.

These cameras are like switching from a blurry 1990s webcam to a 4K Ultra-HD camera. Suddenly, the "freckles" (faculae) were crystal clear.

The Experiment: Testing the Crystal Ball

Foukal looked at the "start of the season" for two recent solar cycles: Cycle 24 (which started around 2009) and Cycle 25 (which started around 2020).

1. Cycle 24: He counted the weeds and flowers. The ratio predicted a medium-sized cycle. Result: The prediction was spot-on. It matched reality almost perfectly.
2. Cycle 25: He did the same thing. The ratio showed a steeper slope (more weeds relative to flowers) than Cycle 24. This meant Cycle 25 would be stronger.

The Big Prediction

In 2022, three years before the Sun was supposed to reach its peak, Foukal used this method to predict that Cycle 25 would be much stronger than other experts thought.

• The "Official" Guess: A panel of international experts (the IPP) predicted a moderate cycle.
• Foukal's Guess: He predicted a much larger, more energetic cycle.

The Result: As of mid-2025, the Sun has indeed been very active, confirming Foukal's prediction was closer to reality than the official panel's. While his specific number was a little high (he predicted a peak of 185, reality was around 160), he was the only one who correctly called the direction of the trend: It's going to be a big one.

Why Does This Work? (The Physics)

Why does the ratio of weeds to flowers matter?

• Sunspots are made of huge, thick magnetic ropes.
• Faculae are made of tiny, thin magnetic threads.

If the Sun's internal engine (the "dynamo") is pumping out a lot of energy, it tends to create more of those thick, powerful ropes (sunspots) relative to the thin threads. So, a high ratio of spots to faculae is a sign that the Sun's engine is revving up to high power.

The Bottom Line

This paper proves that by simply looking at the balance between dark spots and bright freckles at the very start of a solar cycle, we can predict how violent the rest of the cycle will be with incredible accuracy (within 4% error for most cycles).

It's like looking at the size of the first few waves at the beach to know if a tsunami is coming. It's a simple, physical rule that has worked for over 150 years, and thanks to modern space cameras, we can use it to protect our technology from the Sun's next big tantrum.

Technical Summary

Here is a detailed technical summary of the paper "A Remarkably Accurate Predictor of Sunspot Cycle Amplitude" by Peter Foukal.

1. Problem Statement

Predicting the amplitude of upcoming sunspot cycles is a critical challenge in solar physics with significant practical implications for space weather forecasting. Current dynamo models lack the precision to provide reliable physical predictions, forcing researchers to rely on empirical or semi-empirical models.

A specific, promising empirical relationship was identified by Brown and Evans (1980), which correlated the slope of the linear relation between sunspot areas ($A_s$) and white light (WL) facular areas ($A_f$) during the early stages of a cycle with the cycle's eventual peak amplitude. However, this relationship had not been utilized for cycles after Cycle 21 due to a lack of high-quality, consistent data. The Royal Greenwich Observatory (RGO) program, which provided the original data, ended in 1976. Subsequent ground-based data (Mt. Wilson, Kodaikanal) suffered from emulsion inhomogeneities and exposure variations that made accurate facular area measurements difficult.

2. Methodology

The study leverages modern space-based instrumentation to revive and validate the Brown and Evans (1980) method for Solar Cycles 24 and 25.

• Data Sources:
  • SOHO/MDI: Michelson Doppler Imager (narrow-band quasi-continuum at 6768 Å) for Cycle 24 (2009–2010).
  • SDO/HMI: Heliospheric and Magnetic Imager (continuum passband at 6173 Å) for Cycle 25 (2020–2021).
  • Historical Data: Re-calibrated International Sunspot Number ($R_{max}$) and slopes derived by Brown and Evans for Cycles 12–21.
• Image Processing:
  • Images were selected from the first two years of each cycle to capture the onset phase.
  • Projected areas of spots and faculae were measured manually (counting squares on printed images) and corrected for foreshortening, mimicking the RGO procedure.
  • Data was smoothed using a 13-month running mean to reduce short-term variability.
• Analysis:
  • Linear regression was performed to determine the slope ($dA_s/dA_f$) for each cycle.
  • A master correlation plot was constructed linking these slopes to the peak smoothed sunspot numbers ($R_{m-max}$) of Cycles 12–21.
  • The slopes for Cycles 24 and 25 were plugged into this correlation to predict their amplitudes.
• Cross-Validation: The authors addressed discrepancies between MDI and HMI resolutions by comparing co-temporal data during the 2010 overlap period to derive transformation relations, ensuring consistency across datasets.

3. Key Contributions

• Validation of Historical Method: The paper confirms that the linear relationship between early-cycle sunspot and facular areas, and its correlation with cycle amplitude, holds true for the modern space-based era (Cycles 24 and 25), extending the dataset by nearly 50 years.
• High-Precision Prediction: The study demonstrates that this method can predict cycle amplitudes with an accuracy of $\pm$4% RMS for Cycles 12–24.
• Early Prediction of Cycle 25: In 2022 (three years before the predicted maximum), the author used this method to predict a Cycle 25 amplitude ($R_{m-max} \approx 185$) significantly higher than estimates from the International Prediction Panel (IPP).
• Physical Interpretation: The paper provides a physical basis for the empirical correlation, linking the slope change to a shift in the photospheric magnetic flux tube size spectrum. A steeper slope indicates a shift toward larger-scale magnetic structures (larger spots) relative to small-scale faculae as cycle amplitude increases.

4. Results

• Cycle 24: The measured slope was 0.24, predicting a peak amplitude of 115. This matched the observed 13-month smoothed amplitude of 115 (achieved in 2014) almost perfectly.
• Cycle 25: The measured slope was 0.61, predicting a peak amplitude of 185.
  • Discrepancy: The actual observed peak (mid-2025) was approximately 160. The prediction was an overestimate of roughly 16%.
  • Analysis of Error: The authors attribute this overestimate to uncertainties in correcting for the angular resolution differences between HMI (high resolution) and the historical RGO/MDI data. While transformation equations were derived, the low activity levels during the overlap period prevented a definitive correction for facular areas.
• Overall Accuracy: For Cycles 12 through 24, the calculated amplitudes agreed with observed values within $\pm$4% RMS. Table I in the paper details these comparisons, showing deviations ranging from 0 to 15 units for most cycles.
• Comparison with IPP: Despite the overestimation of the absolute value for Cycle 25, the method correctly predicted that Cycle 25 would be larger than the IPP's estimates at the time of the prediction.

5. Significance

• Objective and Physically Based: Unlike purely statistical models, this method relies on a physically grounded relationship between magnetic flux tube sizes (spots vs. faculae) and global dynamo output.
• Lead Time: The technique offers a reliable forecast 3–4 years before a cycle's maximum, provided a stable source of full-disk continuum images is available.
• Limitations and Future Work:
  • The method does not predict the timing of the maximum (though it can be combined with Waldmeier's rule).
  • It cannot be easily applied to Cycles 22 and 23 using Ca K plage data because plages cover a broader range of flux tube sizes than WL faculae, reducing diagnostic sensitivity.
  • Future investigations could apply this method to Cycle 23 using enhanced contrast continuum images from the Kanzelhöhe Observatory.

In conclusion, Foukal's paper establishes a robust, empirically validated, and physically interpretable tool for forecasting sunspot cycle amplitudes, bridging the gap between historical ground-based observations and modern space-based solar monitoring.