Estimating ionospheric currents by inversion from ground-based geomagnetic data and calculating geoelectric fields for studies of geomagnetically induced currents

Abstract

This research focuses on the inversion of geomagnetic variation field measurements to obtain the source currents in the ionosphere and magnetosphere, and to determine the geoelectric fields at the Earth’s surface. During geomagnetic storms, the geoelectric fields create geomagnetically induced currents (GIC) in power networks. These GIC may disturb the operation of power systems, cause damage to power transformers, and even result in power blackouts. In this model, line currents running east–west along given latitudes are postulated to exist at a certain height above the Earth’s surface. This physical arrangement results in the fields on the ground being composed of a zero magnetic east component and a nonzero electric east component. The line current parameters are estimated by inverting Fourier integrals (over wavenumber) of elementary geomagnetic fields using the Levenberg–Marquardt technique. The output parameters of the model are the ionospheric current strength and the geoelectric east component at the Earth’s surface. A conductivity profile of the Earth is adapted from a shallow layered-Earth model for one observatory, together with a deep-layer model derived from satellite observations. This profile is used to obtain the ground surface impedance and therefore the reflection coefficient in the integrals. The inputs for the model are a spectrum of the geomagnetic data for 31 May 2013. The output parameters of the model are spectrums of the ionospheric current strength and of the surface geoelectric field. The inverse Fourier transforms of these spectra provide the time variations on the same day. The geoelectric field data can be used as a proxy for GIC in the prediction of GIC for power utilities. The current strength data can assist in the interpretation of upstream solar wind behaviour.