Observational signatures of core-mantle interactions inferred from top-heavy thermochemical geodynamo simulations

The time-averaged geomagnetic field provides crucial insights into deep Earth dynamics and thermal core-mantle interactions. Paleomagnetic observations and numerical dynamo simulations are equivocal regarding the longitudinal structure of the time-averaged field, though the latter have often considered a generic buoyancy source, which may obscure distinct signatures of thermal and chemical buoyancy that arise near the equator and poles, respectively. In this study, we present a new suite of top-heavy geodynamo simulations, varying the relative strengths of thermal and chemical driving and comparing the resultant magnetic signatures to observational field models spanning centuries to tens of thousands of years. None of the spatially-averaged measures of field morphology and variability we tested could robustly distinguish between different levels of chemical driving or the presence of heterogeneous outer boundary heat flux. On the other hand, observational constraints requiring longitudinal variations in time-averaged inclination anomaly are readily matched by simulations with heterogeneous outer boundary thermal forcing, in contrast to those with homogeneous mantle heat flux. Longitudinal field structures are reduced, but not erased, by elevated chemical driving, which also promotes the formation and deepening of polar minima in the radial magnetic field. Our simulations indicate that both the strong heat flux heterogeneity and chemical driving in Earth’s core are likely to result in small but persistent departures from the geocentric axial dipole approximation.