Chemical Convection and Stratification in the Earth's Outer Core

Convection in the Earth's outer core is driven by buoyancy sources of both thermal and compositional origin. The thermal and compositional molecular diffusivities differ by several orders of magnitude, which can affect the dynamics in various ways. So far, the large majority of numerical simulations have been performed within the codensity framework that consists in combining temperature and composition, assuming artificially enhanced diffusivities for both variables. In this study, we use a particle-in-cell method implemented in a 3D dynamo code to conduct a first qualitative exploration of pure compositional convection in a rotating spherical shell. We focus on the end-member case of infinite Schmidt number by totally neglecting the compositional diffusivity. We show that compositional convection has a very rich physics that deserves several more focused and quantitative studies. We also report, for the first time in numerical simulations, the self-consistent formation of a chemically stratified layer at the top of the shell caused by the accumulation of chemical plumes and blobs emitted at the bottom boundary. When applied to likely numbers for the Earth's core, some (possibly simplistic) physical considerations suggest that a stratified layer formed in such a scenario would be probably weakly stratified and may be compatible with magnetic observations.