AUTHOR=Blandin Matthew , Connor Hyunju K. , Öztürk Doğacan S. , Keesee Amy M. , Pinto Victor , Mahmud Md Shaad , Ngwira Chigomezyo , Priyadarshi Shishir 

TITLE=Multi-Variate LSTM Prediction of Alaska Magnetometer Chain Utilizing a Coupled Model Approach

JOURNAL=Frontiers in Astronomy and Space Sciences

VOLUME=9

YEAR=2022

URL=https://www.frontiersin.org/journals/astronomy-and-space-sciences/articles/10.3389/fspas.2022.846291

DOI=10.3389/fspas.2022.846291

ISSN=2296-987X

ABSTRACT=<p>During periods of rapidly changing geomagnetic conditions electric fields form within the Earth’s surface and induce currents known as geomagnetically induced currents (GICs), which interact with unprotected electrical systems our society relies on. In this study, we train multi-variate Long-Short Term Memory neural networks to predict magnitude of north-south component of the geomagnetic field (|<italic>B</italic><sub><italic>N</italic></sub>|) at multiple ground magnetometer stations across Alaska provided by the SuperMAG database with a future goal of predicting geomagnetic field disturbances. Each neural network is driven by solar wind and interplanetary magnetic field inputs from the NASA OMNI database spanning from 2000–2015 and is fine tuned for each station to maximize the effectiveness in predicting |<italic>B</italic><sub><italic>N</italic></sub>|. The neural networks are then compared against multivariate linear regression models driven with the same inputs at each station using Heidke skill scores with thresholds at the 50, 75, 85, and 99 percentiles for |<italic>B</italic><sub><italic>N</italic></sub>|. The neural network models show significant increases over the linear regression models for |<italic>B</italic><sub><italic>N</italic></sub>| thresholds. We also calculate the Heidke skill scores for d|<italic>B</italic><sub><italic>N</italic></sub>|/dt by deriving d|<italic>B</italic><sub><italic>N</italic></sub>|/dt from |<italic>B</italic><sub><italic>N</italic></sub>| predictions. However, neural network models do not show clear outperformance compared to the linear regression models. To retain the sign information and thus predict <italic>B</italic><sub><italic>N</italic></sub> instead of |<italic>B</italic><sub><italic>N</italic></sub>|, a secondary so-called polarity model is utilized. The polarity model is run in tandem with the neural networks predicting geomagnetic field in a coupled model approach and results in a high correlation between predicted and observed values for all stations. We find this model a promising starting point for a machine learned geomagnetic field model to be expanded upon through increased output time history and fast turnaround times.</p>