Earth’s magnetic field is generated by turbulent motion in its fluid outer core. Although the bulk of the outer core is vigorously convecting and well mixed, some seismic, geomagnetic and geodynamic evidence suggests that a global stably stratified layer exists at the top of Earth’s core. Such a layer would strongly influence thermal, chemical and momentum exchange across the core–mantle boundary and thus have important implications for the dynamics and evolution of the core. Here we argue that the relevant scenario is not global stratification, but rather regional stratification arising solely from the lateral variations in heat flux at the core–mantle boundary. Using our extensive suite of numerical simulations of the dynamics of the fluid core with het- erogeneous core–mantle boundary heat flux, we predict that thermal regional inversion layers extend hundreds of kilometres into the core under anomalously hot regions of the lowermost mantle. Although the majority of the outermost core remains actively convecting, sufficiently large and strong regional inversion layers produce a one-dimensional temperature profile that mimics a globally stratified layer below the core–mantle boundary—an apparent thermal stratification despite the average heat flux across the core–mantle boundary being strongly superadiabatic.
more »
« less
Changes in core-mantle boundary heat flux patterns throughout the supercontinent cycle
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 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 core-mantle boundary heat flux on Earth and its changes over time using compressible global 3-D mantle convection models in the geodynamic modeling 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 core-mantle boundary 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 core-mantle boundary 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% 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.
more »
« less
- Award ID(s):
- 2054605
- PAR ID:
- 10498733
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- ISSN:
- 0956-540X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Large igneous provinces (LIPs) have been linked to both surface and deep mantle processes. During the formation, tenure, and breakup of the supercontinent Pangea, there is an increase in emplacement events for both continental and oceanic LIPs. There is currently no clear consensus on the origin of LIPs, but a hypothesis relates their formation to crustal emplacement of hot plume material originating in the deep mantle. The interaction of subducted slabs with the lowermost mantle thermal boundary and subsequent return-flow is a key control on such plume generation. This mechanism has been explored for LIPs below the interior of a supercontinent (i.e., continental LIPs). However, a number of LIPs formed exterior to Pangea (e.g., Ontong Java Plateau), with no consensus on their formation mechanism. Here, we consider the dynamics of supercontinent processes as predicted by numerical models of mantle convection, and analyse whether circum-supercontinent subduction could generate both interior (continental) and exterior (oceanic) deep-mantle plumes. Our numerical models show that subduction related to the supercontinent cycle can reproduce the location and timing of the Ontong Java Plateau, Caribbean LIP, and potentially the Shatsky Rise, by linking the origin of these LIPs to the return-flow that generated deep mantle exterior plumes.more » « less
-
Abstract Determining the fate of subducted oceanic crust is critical for understanding material cycling through Earth’s deep interior and sources of mantle heterogeneity. A key control on the distribution of subducted slabs over long timescales is the bridgmanite to post-perovskite phase transition in the lowermost mantle, thought to cause rheological weakening. Using high-resolution computational models, we show that the ubiquitous presence of weak post-perovskite at the core-mantle boundary can facilitate or prevent the accumulation of basaltic oceanic crust, depending on the amount of weakening and the crustal thickness. Moderately weak post-perovskite ( ~ 10–100× weaker) facilitates segregation of crust from subducted slabs, increasing basalt accumulation in dense piles. Conversely, very weak post-perovskite (more than 100× weaker) promotes vigorous plumes that entrain more crustal material, decreasing basalt accumulation. Our results reconcile the contradicting conclusions of previous studies and provide insights into the accumulation of subducted crust in the lowermost mantle throughout Earth’s history.more » « less
-
Abstract The high-pressure melting curve of FeO controls key aspects of Earth’s deep interior and the evolution of rocky planets more broadly. However, existing melting studies on wüstite were conducted across a limited pressure range and exhibit substantial disagreement. Here we use an in-situ dual-technique approach that combines a suite of >1000 x-ray diffraction and synchrotron Mössbauer measurements to report the melting curve for Fe1-xO wüstite to pressures of Earth’s lowermost mantle. We further observe features in the data suggesting an order-disorder transition in the iron defect structure several hundred kelvin below melting. This solid-solid transition, suggested by decades of ambient pressure research, is detected across the full pressure range of the study (30 to 140 GPa). At 136 GPa, our results constrain a relatively high melting temperature of 4140 ± 110 K, which falls above recent temperature estimates for Earth’s present-day core-mantle boundary and supports the viability of solid FeO-rich structures at the roots of mantle plumes. The coincidence of the defect order-disorder transition with pressure-temperature conditions of Earth’s mantle base raises broad questions about its possible influence on key physical properties of the region, including rheology and conductivity.more » « less
-
There has been longstanding controversy about whether the influence of lateral variations in core-mantle boundary heat flow can be detected in paleomagnetic records of geomagnetic field behavior. Their signature is commonly sought in globally distributed records of virtual geomagnetic pole (VGP) paths that have been claimed to exhibit specific longitudinal preferences during polarity transitions and excursions. These preferences have often been linked to thermal effects from large low seismic velocity areas (LLVPs) in the lowermost mantle, but the results have been contested because of potential sensitivity to sparse temporal and spatial sampling. Recently developed time varying global paleofield models spanning various time intervals in 1–100 ka, three of which include excursions, allow us to complement assessments of spatial distributions of transitional VGP paths with distributions of minimum field intensity. Robustness of the results is evaluated using similar products from four distinct numerical dynamo simulations with and without variable thermal boundary conditions and including stable geomagnetic polarity, excursions and reversals. We determine that VGP distributions are less useful than minimum field intensity in linking the influences of thermal CMB structure to geographical variations in actual paleofield observables, because VGP correlations depend strongly on good spatial sampling of a sufficient number of relatively rare events. These results provide a basis for evaluating comparable observations from four paleofield models. The distribution of VGP locations provide unreliable results given the restricted time span and available data locations. Rough correlations of global distributions of minimum intensity with areas outside the LLVPs give some indications of mantle control during excursions, although the results for the eastern hemisphere are complex, perhaps highlighting uncertainties about the hemispheric balance between thermal and compositional variations in the lowermost mantle. However, access to other geomagnetic properties (such as intensity and radial field at the CMB) provides a strong argument for using extended and improved global paleofield models to resolve the question of mantle influence on the geodynamo from the observational side.more » « less
An official website of the United States government

