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
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.
more » « less- Award ID(s):
- 2054605
- NSF-PAR ID:
- 10499911
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 237
- Issue:
- 3
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 1251-1274
- Size(s):
- p. 1251-1274
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core‐mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co‐constrain core‐mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports
Q CMB∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long‐lived basal mantle molten layer through much of Earth's history. -
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
-
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
-
more » « less
SUMMARY Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (Sdiff) are potentially strongly influenced by lowermost mantle anisotropy. Sdiff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. Sdiff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SHdiff–SVdiff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of Sdiff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which Sdiff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on Sdiff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SHdiff waves. In particular, we investigate how SHdiff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to Sdiff splitting that can be clearly distinguished from other waveform effects. These results enable us to lay out a strategy for the analysis of Sdiff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of Sdiff and the influence of source-side anisotropy. We demonstrate our Sdiff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained Sdiff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.
