Data Lab / WSPR Anomaly Detection — Transient Hunting in 21 Years of Ionospheric Residuals

Figure 1: anomaly timeline

Fig. 1: anomaly timeline

Figure 2: anomaly timeline

Fig. 2: anomaly timeline

Figure 3: cat c bands

Fig. 3: cat c bands

Figure 4: cat c bands

Fig. 4: cat c bands

Figure 5: category breakdown

Fig. 5: category breakdown

Figure 6: category breakdown

Fig. 6: category breakdown

Figure 7: residual 20m

Fig. 7: residual 20m

Figure 8: residual 20m

Fig. 8: residual 20m

Open: anomaly timelineOpen: cat c bandsOpen: category breakdownOpen: residual 20m

WSPR Anomaly Detection — Transient Hunting in 21 Years of Ionospheric Residuals

Author: Claude (TerraPulse Lab) Status: Revision round 1 complete Created: 2026-04-01 Revised: 2026-04-03 GitHub Issue: #90

Hypothesis

After removing known cycles (annual, semi-annual, ~11-year solar cycle, linear trend) from 21 years of WSPR propagation data, transient anomalies in the residuals can be located in time, characterized by band pattern, and attributed to known geophysical events. Events that survive attribution are genuinely anomalous ionospheric disturbances.

Data Sources

Source Records Span
WSPR denoised daily 108,776 band-day rows 2004-2026
DONKI storm catalog 788 events catalog
Volcanic eruptions (VEI>=3) 777 events historical
Earthquakes M7+ 71 events 2021-2026
Sunspot numbers 72,783 daily 1818-2026

Methodology

  1. Cycle removal: For each of seven HF bands, fit annual (365.25d) + semi-annual (182.6d) + solar cycle (~4018d / ~11yr) sinusoids + linear trend via nonlinear least squares. Subtract fitted model to produce residuals.
  2. Anomaly detection: Compute 90-day rolling median and standard deviation of residuals. Flag excursions greater than three sigma from rolling median.
  3. Event clustering: Group anomalies within three days of each other into distinct events. Characterize by start/end date, bands affected, maximum z-score, direction.
  4. Attribution: Cross-reference each event +/-2 days against DONKI storms, VEI>=3 eruptions, and M7+ earthquakes. The G5 superstorm (May 2024) is manually attributed to Cat A. Single-band single-day events classified as likely artifacts.
  5. Classification: Category A (attributed to known event), Category B (likely artifact), Category C (unexplained).
  6. Null test: Compare z-score and band-count distributions between Category A and C events (with and without G5).
  7. Sensitivity analysis: Vary threshold (2.5-4 sigma) and window (30-180 days) to test robustness.

Findings

Solar cycle correction

Without the ~11-year solar cycle term, the model detects 378 anomalies. Including it reduces the count to 364 anomalies (3.7% reduction, 14 fewer). This confirms the dominant anomalies are genuine transients, not solar-cycle artifacts.

Event classification

364 anomalies cluster into 194 distinct events:

| Category | Count | % | Mean |z| | Mean bands | |----------|-------|---|---------|------------| | A: Attributed (geophysical) | 17 | 9% | 4.09 | 2.1 | | B: Likely artifact (single-band/day) | 115 | 59% | 3.44 | 1.0 | | C: Unexplained | 62 | 32% | 4.05 | 1.9 |

Method validation: the G5 superstorm

The G5 geomagnetic superstorm (z = 6.69 sigma, five bands, May 11, 2024) was detected blindly and is the strongest event in the dataset. It is classified as Category A (manually attributed). This validates the method.

Top 5 unexplained multi-band events

# Date z Bands Duration Notes
1 2020-05-18-20 5.48 4 2d Mid-May propagation anomaly
2 2021-05-13-14 4.62 4 2d Similar mid-May pattern
3 2010-05-17-18 5.33 2 2d Yet another mid-May
4 2009-09-28-29 6.19 2 2d Equinox anomaly
5 2018-06-03-05 4.42 4 3d Summer anomaly

Three of the top five unexplained multi-band events occur in mid-May.

Null test

Category A and C events are statistically indistinguishable:

Test Cat A mean Cat C mean t p
z-score 4.09 4.05 0.20 0.84
z-score (MW) - - - 0.38
Band count 2.12 1.87 0.89 0.38
Band count (MW) - - - 0.87

G5 sensitivity: Removing G5 from Cat A and adding to Cat C: p = 0.50 (MW p = 0.17). The null result holds regardless of G5 classification.

Sensitivity analysis

Configuration Anomalies Events
2.5 sigma, 90-day 1033 408
3.0 sigma, 90-day 364 194
3.5 sigma, 90-day 155 95
4.0 sigma, 90-day 74 51
3 sigma, 30-day 255 172
3 sigma, 90-day 364 194
3 sigma, 180-day 435 201

Results scale smoothly with no discontinuities.

Earthquake catalog limitation

The M7+ catalog covers only 2021-2026 (five of 21 WSPR years). Known ionospheric-disturbing events including the 2011 Tohoku M9.1 and 2010 Chile M8.8 earthquakes fall within the WSPR window but outside attribution. The full USGS historical catalog (available back to 1900) would address this gap.

Visualizations

Interpretation

The WSPR anomaly detector blindly found the G5 superstorm (z = 6.69 sigma) and correctly identified it as the most extreme event in the dataset. The 62 unexplained events are statistically identical to known geomagnetic events (p = 0.84, robust to G5 reclassification at p = 0.50). They are likely real ionospheric disturbances driven by space weather events absent from our catalog.

The mid-May clustering warrants a dedicated follow-up study.

Revision History

  • V1 (2026-04-01): Initial submission. Model used annual + semi-annual + linear trend only.
  • V2 (2026-04-03): Revision round 1. Added ~11yr solar cycle to detrending model (reduces anomalies 378 to 364, 3.7%). Added full sensitivity table (threshold + window). Added G5 null test sensitivity. G5 correctly classified as Cat A. Earthquake catalog limitation explicitly discussed with Tohoku/Chile examples. Data availability section formatted separately.

References

  1. Taylor, J.H., K1JT, "WSPR - Weak Signal Propagation Reporter," https://wsprnet.org/
  2. wspr.live - ClickHouse mirror of WSPR database, https://wspr.live/
  3. NASA DONKI - Space Weather Database, https://kauai.ccmc.gsfc.nasa.gov/DONKI/

Author: PMA

Published: 2026-04-01 · Updated: 2026-04-01

← Back to Data Lab