pith. sign in

arxiv: 2507.03538 · v6 · pith:MNTUY3SZnew · submitted 2025-07-04 · 🌌 astro-ph.EP · physics.geo-ph

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

Souvik Naskar , Jonathan E. Mound , Christopher J. Davies , Andrew T. Clarke This is my paper

Pith reviewed 2026-05-19 06:11 UTC · model grok-4.3

classification 🌌 astro-ph.EP physics.geo-ph
keywords outer core convectionthermochemical convectioncore-mantle boundarystable regionsregional inversion lensesheterogeneous heat fluxseismic detectabilitygeomagnetic signatures
0
0 comments X

The pith

Outer core convection models with heterogeneous CMB heat flux produce diverse local and global stable regions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines whether buoyancy forces from heat and light elements at the inner core boundary drive full convection in Earth's outer core or permit stable layers near the core-mantle boundary. Simulations of thermal, chemical, and thermochemical convection at fixed Ekman number explore a range of thermal and compositional forcing strengths, showing that stable regions form locally near the poles or as broader layers depending on the balance of those forces. Heterogeneous heat flux at the boundary generates thermally stratified zones that coexist with convective areas. A sympathetic reader would care because such regions could alter interpretations of seismic data from the core and influence patterns in the geomagnetic field generated by core flow.

Core claim

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. Simulations incorporating CMB heat flux heterogeneities produce locally stable regional inversion lenses rather than a global layer, allowing stable and unstable regions to coexist. Purely chemical simulations accumulate light elements below the CMB, forming locally stable regions near the poles or global layers depending on compositional forcing strength. These chemically stratified regions persist in thermochemical cases, a

What carries the argument

Regional inversion lenses formed by heterogeneous core-mantle boundary heat flux, which create locally stable zones amid thermochemical convection and allow stable and convective regions to coexist.

If this is right

  • Stable regions form in a range of locations, thicknesses, and strengths that depend on the relative thermal and compositional forcing.
  • Chemically stratified stable zones persist even when thermal buoyancy is strongly destabilizing.
  • Heterogeneous heat flux produces thermally stratified local stable regions that coexist with convection.
  • These regions reach thicknesses and strengths that could be detected seismically and might imprint on geomagnetic observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Seismic data could be used to map the specific locations and extents of these stable regions in the real core.
  • The diversity of stable morphologies suggests that single global-layer models may oversimplify core dynamics.
  • Coupling these convection models to evolving mantle heat flux patterns could test whether stable regions migrate over geological time.

Load-bearing premise

The fixed Ekman number, Prandtl numbers, limited Rayleigh number ranges, and imposed heterogeneous heat flux pattern are assumed to capture essential Earth-like core behavior without full mantle coupling or realistic diffusivities.

What would settle it

Seismic observations showing no thick, strong stable layers near the core-mantle boundary or geomagnetic records lacking predicted signatures from regional stable zones would rule out the simulated stable regions.

Figures

Figures reproduced from arXiv: 2507.03538 by Andrew T. Clarke, Christopher J. Davies, Jonathan E. Mound, Souvik Naskar.

Figure 1
Figure 1. Figure 1: Regimes of force balance in homogeneous models. The symbol shapes indicate purely thermal [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Scalar fields and their gradients in homogeneous simulations. The equatorial, meridional and [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Scalar fields and their gradients in heterogeneous simulations. The interpretations of the contour [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Variation of N2/4Ω2 with depth below CMB near the top of the core for homogeneous (a,b) and heterogeneous (c,d) models, averaged around various locations (see Figure A.8 in SI) for models with compositionally dominated convection (a,c) at RafT = 90 and Rafξ = 30000, and thermally dominated convection (b,d) at RafT = 1200 and Rafξ = 300. 17 [PITH_FULL_IMAGE:figures/full_fig_p019_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Regimes of thermochemical stability in the heterogeneous models. The symbols are colored by [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Scaling of thickness (a,b) and strength (c,d) of the stable regions averaged around Africa (a,c) [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Scalar fields and their gradients in a heterogeneous model at [PITH_FULL_IMAGE:figures/full_fig_p027_7.png] view at source ↗
read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports numerical simulations of thermal, chemical, and thermochemical convection in a spherical shell model of Earth's outer core at fixed Ekman number E=10^{-5} and Prandtl numbers Pr_T=1, Pr_ξ=10. With thermal and chemical flux Rayleigh numbers scanned over modest ranges and with imposed heterogeneous CMB heat flux, the authors find that chemically stratified regions form and persist even under destabilizing thermal buoyancy, while heterogeneous thermal forcing produces regional inversion lenses (RILs). The work concludes that these stable regions exhibit diverse locations, thicknesses, and morphologies that could be seismically detectable and carry geomagnetic signatures.

Significance. If the reported stable-region properties prove robust, the results would advance understanding of possible coexisting stable and convective regions in the outer core by explicitly separating thermal and compositional buoyancy sources under heterogeneous boundary conditions. The systematic exploration of purely chemical, thermal, and combined cases, together with the identification of RILs rather than global layers, provides a useful parameter-space map. Credit is due for the clear demonstration that chemically stable regions survive in thermochemical runs and for the discussion of observable signatures.

major comments (2)
  1. [Abstract and Simulation Setup] Abstract and parameter choices: the central claim that stable regions 'can have a seismically detectable thickness and strength' and 'might also have a signature in geomagnetic observations' rests on simulations performed exclusively at E=10^{-5} (Earth value ~10^{-15}), Pr_T=1, Pr_ξ=10, and limited Ra_T, Ra_ξ ranges. No sensitivity tests or scaling arguments are presented to show that RIL thickness, location, or persistence survive at lower Ekman number or altered diffusivity ratios, where the relative importance of thermal versus compositional buoyancy and the penetration of downwellings can change substantially.
  2. [Methods / Numerical Methods] Numerical validation: the manuscript reports outcomes from simulations at the stated parameter values but supplies no information on grid resolution, time-step convergence, or benchmark comparisons against known solutions for rotating convection or double-diffusive cases. Without these checks, it is impossible to confirm that the reported layer thicknesses and morphologies are free of numerical artifacts, which directly affects the reliability of the seismically detectable strength claim.
minor comments (2)
  1. [Introduction / Governing Equations] The notation for the flux Rayleigh numbers (widetilde{Ra}_T and widetilde{Ra}_ξ) is introduced in the abstract but should be defined explicitly with reference to the governing equations in the main text.
  2. [Results / Figures] Figures illustrating stable-region morphology would benefit from quantitative insets or tables reporting measured layer thickness and buoyancy frequency for each Ra combination.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the significance of our results on coexisting stable and convective regions in thermochemical core convection and for the constructive comments. We address each major comment below and will revise the manuscript to incorporate the requested clarifications and discussion.

read point-by-point responses

  1. Referee: Abstract and parameter choices: the central claim that stable regions 'can have a seismically detectable thickness and strength' and 'might also have a signature in geomagnetic observations' rests on simulations performed exclusively at E=10^{-5} (Earth value ~10^{-15}), Pr_T=1, Pr_ξ=10, and limited Ra_T, Ra_ξ ranges. No sensitivity tests or scaling arguments are presented to show that RIL thickness, location, or persistence survive at lower Ekman number or altered diffusivity ratios, where the relative importance of thermal versus compositional buoyancy and the penetration of downwellings can change substantially.

    Authors: We agree that E=10^{-5} is larger than the estimated Earth value and that the manuscript would benefit from explicit discussion of this limitation. This Ekman number was selected to enable a broad scan of thermal and chemical Rayleigh numbers within available computational resources while still capturing rotational effects. The persistence of chemically stable regions under destabilizing thermal buoyancy is driven primarily by the diffusivity contrast (Pr_ξ=10 versus Pr_T=1), which favors compositional stratification; this mechanism is expected to remain qualitatively robust at lower E, although quantitative layer thicknesses may decrease. In the revised manuscript we will add a dedicated subsection in the Discussion that references existing scaling relations for stable-layer thickness with Ekman number and Prandtl ratio, and that estimates how the reported RIL properties might change at more extreme parameters. This addition will qualify the detectability claims without requiring new simulations. revision: yes

  2. Referee: Numerical validation: the manuscript reports outcomes from simulations at the stated parameter values but supplies no information on grid resolution, time-step convergence, or benchmark comparisons against known solutions for rotating convection or double-diffusive cases. Without these checks, it is impossible to confirm that the reported layer thicknesses and morphologies are free of numerical artifacts, which directly affects the reliability of the seismically detectable strength claim.

    Authors: We acknowledge that the original manuscript omitted explicit numerical validation details. All simulations were performed at resolutions sufficient to resolve the Ekman boundary layers and the radial structure of the stable regions (typically 128–256 radial points and spherical harmonic degrees up to 128–256, with time steps satisfying a Courant number <0.5). Internal resolution-doubling tests on representative cases confirmed that layer thicknesses and morphologies converged to within a few percent. In the revised manuscript we will insert a new paragraph in the Methods section that reports the grid parameters, time-stepping criteria, and convergence tests. We will also add a short appendix comparing our purely thermal convection diagnostics (e.g., Nusselt numbers and flow morphology) against published benchmark solutions for rotating spherical-shell convection at E=10^{-5}. These additions will directly address concerns about numerical artifacts affecting the seismically detectable strength claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity: results are direct numerical outputs

full rationale

The paper reports outcomes from direct numerical integration of the thermochemical convection equations at fixed Ekman number E=10^{-5}, Prandtl numbers Pr_T=1 and Pr_xi=10, and scanned ranges of flux Rayleigh numbers. Stable-region locations, thicknesses, and morphologies are computed results under the chosen boundary conditions and parameter values; they are not obtained by fitting a parameter to a subset of the same data and then relabeling it a prediction, nor by any self-referential definition that equates an output to an input by construction. No load-bearing analytical step reduces the reported diversity or detectability claims to quantities defined inside the paper itself. The work is therefore self-contained as a computational survey.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard fluid-dynamics approximations and a set of dimensionless parameters chosen to reach computationally accessible regimes that are hoped to be representative of the core.

free parameters (3)
  • Ekman number E = 10^{-5}
    Fixed at 10^{-5} to enable simulation at low viscosity; actual core value is orders of magnitude smaller.
  • Prandtl numbers Pr_T and Pr_xi = 1 and 10
    Set to 1 and 10 respectively to represent differing thermal and chemical diffusivities.
  • Thermal and chemical flux Rayleigh numbers = 30-4000 and 30-100000
    Ranges 30-4000 and 30-100000 chosen to span stable-to-unstable regimes.
axioms (2)
  • standard math Boussinesq approximation for buoyancy-driven flow
    Density variations are treated only in the buoyancy term; invoked throughout the model setup.
  • domain assumption Imposed heterogeneous CMB heat flux boundary condition
    Heat flux pattern is prescribed rather than emerging from coupled mantle convection.

pith-pipeline@v0.9.0 · 5817 in / 1546 out tokens · 104116 ms · 2026-05-19T06:11:07.346739+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.