Core-surface kinematic control of polarity reversals in advanced geodynamo simulations

Abstract

The geomagnetic field has undergone hundreds of polarity reversals over Earth's history, at a variable pace. In numerical models of Earth's core dynamics, reversals occur with increasing frequency when the convective forcing is increased past a critical level. This transition has previously been related to the influence of inertia in the force balance. Because this force is subdominant in Earth's core, concerns have been raised regarding the geophysical applicability of this paradigm. Reproducing the reversal rate of the past million years also requires forcing conditions that do not guarantee that the rest of the geomagnetic variation spectrum is reproduced. These issues motivate the search for alternative reversal mechanisms. Using a suite of numerical models where buoyancy is provided at the bottom of the core by inner-core freezing, we show that the magnetic dipole amplitude is controlled by the relative strength of subsurface upwellings and horizontal circulation at the core surface. A relative weakening of upwellings brings the system from a stable to a reversing dipole state. This mechanism is purely kinematic because it operates irrespectively of the interior force balance. It is therefore expected to apply at the physical conditions of Earth's core. Subsurface upwellings may be impeded by stable stratification in the outermost core. We show that with weak stratification levels corresponding to a nearly adiabatic core surface heat flow, a single model reproduces the observed geomagnetic variations ranging from decades to millions of years. In contrast with the existing paradigm, reversals caused by this stable top core mechanism become more frequent when the level of stratification increases i.e. when the core heat flow decreases. This suggests that the link between mantle dynamics and magnetic reversal frequency needs to be reexamined.