  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 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. 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
regularisation 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.
