Polar ionospheric currents and high temporal resolution geomagnetic field models

Estimating high resolution models of the Earth's core magnetic field and its time variation in the polar regions requires that one can adequately account for magnetic signals produced by polar ionospheric currents, which vary on a wide range of time and length scales. Limitations of existing ionospheric field models in the challenging polar regions can adversely affect core field models, which in turn has important implications for studies of the core flow dynamics in those regions. Here we implement a new approach to co-estimate a climatological model of the ionospheric field together with a model of the internal and magnetospheric fields within the CHAOS geomagnetic field modelling framework. The parametrization of the ionospheric field exploits non-orthogonal magnetic coordinates to efficiently account for the geometry of the Earth's magnetic field and scales linearly with external driving parameters related to the solar wind and the interplanetary magnetic field. Using this approach we derive a new geomagnetic field model from measurements of the magnetic field collected by low Earth orbit satellites, which in addition to the internal field provides estimates of the typical current system in the polar ionosphere and successfully accounts for previously unmodelled ionospheric signals in field model residuals. To resolve the ambiguity between the internal and ionospheric fields when using satellite data alone, we impose regularization. We find that the time derivative of the estimated internal field is less contaminated by the polar currents, which is mostly visible in the zonal and near-zonal terms at high spherical harmonic degrees. Distinctive patches of strong secular variation at the core-mantle boundary, which have important implications for core dynamics, persist. Relaxing the temporal regularization reveals annual oscillations, which could indicate remaining ionospheric field or related induced signals in the internal field model. Using principal component analysis we find that the annual oscillations mostly affect the zonal low-degree spherical harmonics of the internal field.