Thermochemical models of outer core convection with heterogeneous core-mantle boundary heat flux

TL;DR

This study uses advanced simulations to explore how thermal and chemical heterogeneities at Earth's core boundary influence convection, stable layer formation, and potential observable signatures.

Contribution

It introduces comprehensive thermochemical convection models at realistic parameters, revealing diverse stable region behaviors influenced by thermal and compositional buoyancy.

Findings

Chemically stratified regions form below the CMB depending on $ ilde{Ra}_\xi$.

Thermally stratified regional inversion lenses (RILs) persist even with destabilising compositional buoyancy.

Stable regions vary in location, morphology, and seismic detectability based on key parameters.

Abstract

Convection in Earth's outer core is driven by the release of heat and light elements at the inner core boundary. A key question is whether these buoyancy sources drive convection throughout the core, or whether a stable layer exists just below the core-mantle boundary (CMB). Recent simulations incorporating CMB heat flux heterogeneities propose locally stable ``regional inversion lenses'' (RILs) rather than a global layer, allowing stable and unstable regions to coexist. However, these simulations combine thermal and compositional anomalies, ignoring differences in diffusivities and boundary conditions. Here we simulate thermal, chemical, and thermochemical convection at Ekman number $E=10^{-5}$, with thermal and chemical flux Rayleigh numbers $\widetilde{Ra}_T=30-4000$ and $\widetilde{Ra}_\xi=30-100000$, and Prandtl numbers $Pr_T=1$ and $Pr_\xi=10$. Purely chemical simulations accumulate light elements below the CMB, forming locally stable regions near the poles or global layers, depending on $\widetilde{Ra}_\xi$. These chemically stratified regions persist in thermochemical simulations even when thermal forcing is destabilising. Introducing heterogeneous CMB heat flux produces thermally stratified RILs even with strongly destabilising compositional buoyancy. Our simulations reveal a diverse range of locations, properties, and morphologies of stable regions depending on $\widetilde{Ra}_T$ and $\widetilde{Ra}_\xi$, they can have a seismically detectable thickness and strength and might also have a signature in geomagnetic observations.