- Award ID(s):
- 1853196
- NSF-PAR ID:
- 10137699
- Date Published:
- Journal Name:
- Nature geoscience
- ISSN:
- 1752-0894
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Light elements in Earth’s core play a key role in driving convection and influencing geodynamics, both of which are crucial to the geodynamo. However, the thermal transport properties of iron alloys at high-pressure and -temperature conditions remain uncertain. Here we investigate the transport properties of solid hexagonal close-packed and liquid Fe-Si alloys with 4.3 and 9.0 wt % Si at high pressure and temperature using laser-heated diamond anvil cell experiments and first-principles molecular dynamics and dynamical mean field theory calculations. In contrast to the case of Fe, Si impurity scattering gradually dominates the total scattering in Fe-Si alloys with increasing Si concentration, leading to temperature independence of the resistivity and less electron–electron contribution to the conductivity in Fe-9Si. Our results show a thermal conductivity of ∼100 to 110 W⋅m −1 ⋅K −1 for liquid Fe-9Si near the topmost outer core. If Earth’s core consists of a large amount of silicon (e.g., > 4.3 wt %) with such a high thermal conductivity, a subadiabatic heat flow across the core–mantle boundary is likely, leaving a 400- to 500-km-deep thermally stratified layer below the core–mantle boundary, and challenges proposed thermal convection in Fe-Si liquid outer core.more » « less
-
SUMMARY The Earth’s magnetic field is generated by a dynamo in the outer core and is crucial for shielding our planet from harmful radiation. Despite the established importance of the core–mantle boundary (CMB) heat flux as driver for the dynamo, open questions remain about how heat flux heterogeneities affect the magnetic field. Here, we explore the distribution of the CMB heat flux on Earth and its changes over time using compressible global 3-D mantle convection models in the geodynamic modelling software ASPECT. We discuss the use of the consistent boundary flux method as a tool to more accurately compute boundary heat fluxes in finite element simulations and the workflow to provide the computed heat flux patterns as boundary conditions in geodynamo simulations. Our models use a plate reconstruction throughout the last 1 billion years—encompassing the complete supercontinent cycle—to determine the location and sinking speed of subducted plates. The results show how mantle upwellings and downwellings create localized heat flux anomalies at the CMB that can vary drastically over Earth’s history and depend on the properties and evolution of the lowermost mantle as well as the surface subduction zone configuration. The distribution of hot and cold structures at the CMB changes throughout the supercontinent cycle in terms of location, shape and number, indicating that these structures fluctuate and might have looked very differently in Earth’s past. We estimate the resulting amplitude of spatial heat flux variations, expressed by the ratio of peak-to-peak amplitude to average heat flux, q*, to be at least 2. However, depending on the material properties and the adiabatic heat flux out of the core, q* can easily reach values >30. For a given set of material properties, q* generally varies by 30–50 per cent over time. Our results have implications for understanding the Earth’s thermal evolution and the stability of its magnetic field over geological timescales. They provide insights into the potential effects of the mantle on the magnetic field and pave the way for further exploring questions about the nucleation of the inner core and the past state of the lowermost mantle.
-
We numerically and theoretically investigate the Boussinesq Eady model, where a rapidly rotating density-stratified layer of fluid is subject to a meridional temperature gradient in thermal wind balance with a uniform vertically sheared zonal flow. Through a suite of numerical simulations, we show that the transport properties of the resulting turbulent flow are governed by quasigeostrophic (QG) dynamics in the rapidly rotating strongly stratified regime. The ‘vortex gas’ scaling predictions put forward in the context of the two-layer QG model carry over to this fully three-dimensional system: the functional dependence of the meridional flux on the control parameters is the same, the two adjustable parameters entering the theory taking slightly different values. In line with the QG prediction, the meridional heat flux is depth-independent. The vertical heat flux is such that turbulence transports buoyancy along isopycnals, except in narrow layers near the top and bottom boundaries, the thickness of which decreases as the diffusivities go to zero. The emergent (re)stratification is set by a simple balance between the vertical heat flux and diffusion along the vertical direction. Overall, this study demonstrates how the vortex-gas scaling theory can be adapted to quantitatively predict the magnitude and vertical structure of the meridional and vertical heat fluxes, and of the emergent stratification, without additional fitting parameters.more » « less
-
Abstract Convection in planetary mantles is in the so‐called mixed heating mode; it is driven by heating from below, due to a hotter core, as well as heating from within, due to radiogenic heating and secular cooling. Thus, in order to model the thermal evolution of terrestrial planets, we require the parameterization of heat flux for mixed heated convection in particular. However, deriving such a parameterization from basic principles is an elusive task. While scaling laws for purely internal heating and purely basal heating have been successfully determined using the idea that thermal boundary layers are marginally stable, recent theoretical analyses have questioned the applicability of this idea to convection in the mixed heating mode. Here, we present a scaling approach that is rooted in the physics of convection, including the boundary layer stability criterion. We show that, as long as interactions between thermal boundary layers are properly accounted for, this criterion succeeds in describing relationships between thermal boundary layer (TBL) properties for mixed heated convection. The surface heat flux of a convecting fluid is locally determined by the properties of the upper TBL, as opposed to globally determined. Our foundational scaling approach can be readily extended to nearly any complexity of convection within planetary mantles.
-
Abstract Understanding Jupiter's present‐day interior structure and dynamics is key to constraining planetary accretion models. In particular, the extent of stable stratification (i.e., non‐convective regions) in the planet strongly influences long‐term cooling processes, and may record primordial heavy element gradients from early in a planet's formation. Because the Galileo entry probe measured a subsolar helium abundance, Jupiter interior models often invoke an outer stably stratified region due to helium rain. Additionally, Juno gravity data suggest a deeper, potentially stratified dilute core extending halfway through the planet. However, fits to Jupiter's gravitational data are non‐unique, and outstanding uncertainty over the equations of state for hydrogen and helium remain. Here, we use high‐resolution numerical magnetohydrodynamic simulations of Jupiter's magnetic field to place constraints on the extent of stable stratification within the planet. We find that compared to traditional interior models, an upper stably stratified layer between 0.9 and 0.95 Jupiter radii (
R J ) helps to explain both Jupiter's dipolar magnetic field and zonal winds. In contrast, an extended dilute core that is entirely stably stratified (no convective layers) yields significantly worse fits to both. However, our models with extended deep stratification still generate dipolar magnetic fields if an upper stratified region is also present. Overall, we find that a planet with a dilute core i.e., strongly stably stratified is increasingly challenging to reconcile with Jupiter's magnetic field and winds. Thus if a dilute core is present, alternative modalities such as a fully convective dilute core, a complex multilayered interior structure, or double diffusive convection may be required.