Cosmic-Ray-Constrained LSTM Model for Geomagnetic Storm Prediction
Zongyuan Ge, Chenwaner Zhang, Wei Zhou, Hongyu Zeng, Guiping Zhou
TL;DR
This work develops a physics-informed LSTM for geomagnetic storm prediction by explicitly incorporating cosmic-ray modulation from neutron monitors as a predictive precursor. Leveraging 1995–2020 multi-source data at 1-hour cadence, the model fuses solar-wind/IMF inputs with cosmic-ray background, transient anomalies, and inter-station coherence to extend lead times and improve 48-hour event-detection performance, achieving RMSEs from $5.106$ to $14.788$ nT and a $25.84\%$ relative gain in F1-score over a purely solar-wind-based model. A targeted optimization reduces input dimensionality and tunes the history window, with a single-layer, 50-unit LSTM yielding the best generalization and a runtime suitable for near-real-time forecasting. The findings support the value of physics-informed ML for space-weather, and point to future directions in multi-task learning, uncertainty quantification, and physics-constraint integration to further enhance reliability and interpretability.
Abstract
Geomagnetic storms (GSTs) driven by solar wind-magnetosphere coupling can severely disrupt technological systems, motivating the need for improved prediction accuracy and longer warning times. In this study, we develop a physics-informed Long Short-Term Memory (LSTM) model that incorporates cosmic-ray flux modulation as an additional precursor signal. As coronal mass ejection (CME)-driven disturbances propagate through the heliosphere, enhanced turbulence and magnetic-field compression reduce galactic cosmic-ray (GCR) flux measured by ground-based neutron monitors (Forbush decreases), providing early information that can precede near-Earth solar-wind signatures by 1--3 days. We integrate multi-source space-weather data, spanning 1995-2020, including cosmic-ray observations, solar wind plasma parameters, interplanetary magnetic-field data, and geomagnetic indices. Based on these data, we construct a 19-dimensional feature vector that includes flux background levels, decrease-related indicators, and inter-station correlation measures as model inputs. Employing a 50-unit LSTM architecture, the proposed model achieves root-mean-square errors (RMSE) of 5.106 nT, 8.315 nT, 10.854 nT, 12.883 nT, and 14.788 nT for 2-, 6-, 12-, 24-, and 48-hour predictions, respectively. Incorporating cosmic-ray information further improves 48-hour forecast skill by up to 25.84% (from 0.178 to 0.224). These results demonstrate the value of physics-informed deep learning and cosmic-ray precursors for advancing space-weather forecasting.
