skip to main content

Title: Multidisciplinary Constraints on the Thermal‐Chemical Boundary Between Earth's Core and Mantle

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 supportsQCMB∼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.

more » « less
Award ID(s):
1554554 1664332
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Earth's lowermost mantle is characterized by two large low shear velocity provinces (LLSVPs). The regions outside the LLSVPs have been suggested to be strongly influenced by subducted slabs and, therefore, much colder than the LLSVPs. However, localized low‐velocity seismic anomalies have been detected in the subduction‐influenced regions, whose origin remains unclear. Here, three‐dimensional geodynamic calculations are performed, and they show that linear, ridge‐like hot thermal anomalies, or thermal ridges, form in the relatively cold, downwelling regions of the lowermost mantle. Like the formation of Richter rolls due to sublithosphere small‐scale convection (SSC), the thermal ridges form as a result of SSC from the basal thermal boundary layer and they extend in directions parallel to the surrounding mantle flow. The formation of thermal ridges in subduction regions of the lowermost mantle is very sensitive to the thermal structures of the subducted materials, and thermal heterogeneities brought to the bottom of the mantle by subducting slabs greatly promote the formation of thermal ridges. The formation of thermal ridges is also facilitated by the increase of core‐mantle boundary heat flux and vigor of lowermost mantle convection. The thermal ridges may explain the low‐velocity seismic anomalies outside of the LLSVPs in the lowermost mantle. The results suggest that the relatively cold, subduction‐influenced regions of the Earth's lowermost mantle may contain localized hot anomalies.

    more » « less
  2. Abstract

    An internally generated magnetic field once existed on the Moon. This field reached high intensities (∼10–100μT, perhaps intermittently) from ∼4.3 to 3.6 Gyr ago and then weakened to ≲5μT before dissipating by ∼1.9–0.8 Gyr ago. While the Moon’s metallic core could have generated a magnetic field via a dynamo powered by vigorous convection, models of a core dynamo often fail to explain the observed characteristics of the lunar magnetic field. In particular, the core alone may not contain sufficient thermal, chemical, or radiogenic energy to sustain the high-intensity fields for >100 Myr. A recent study by Scheinberg et al. suggested that a dynamo hosted in electrically conductive, molten silicates in a basal magma ocean (BMO) may have produced a strong early field. However, that study did not fully explore the BMO’s coupled evolution with the core. Here we show that a coupled BMO–core dynamo driven primarily by inner core growth can explain the timing and staged decline of the lunar magnetic field. We compute the thermochemical evolution of the lunar core with a 1D parameterized model tied to extant simulations of mantle evolution and BMO solidification. Our models are most sensitive to four parameters: the abundances of sulfur and potassium in the core, the core’s thermal conductivity, and the present-day heat flow across the core–mantle boundary. Our models best match the Moon’s magnetic history if the bulk core contains ∼6.5–8.5 wt% sulfur, in agreement with seismic structure models.

    more » « less
  3. Abstract

    Temporal trends in the paleomagnetic dipole moment exhibit the property of positive skewness. On average, positive trends are larger and occur less frequently than negative trends over timescales of several tens of kyr. We explore the origin of this property using numerical geodynamo models. A suite of models reveals that skewness arises for a restricted set of boundary conditions. Models driven by heat flow at the top and bottom boundaries exhibit very little skewness, whereas models driven solely by heat flow on the lower boundary produce significant positive skewness. Further increases in skewness occur in the presence of thermal stratification at the top of the core. The level of skewness in the geodynamo models is correlated with estimates of upwelling near the core‐mantle boundary. Sustained upwelling is expected to increase magnetic‐flux expulsion, contributing to higher levels of skewness. Similar behavior is recovered from stochastic models in which the dipole is generated by a random series of cyclonic convection events. Skewness in the stochastic models is quantitatively similar to estimates from the geodynamo models when the average recurrence time of the convection events is 100 years. Extending the stochastic models to the paleomagnetic field implies a longer recurrence time of 1,000 years or more. We interpret this recurrence time in terms of the timing of flux‐expulsion events rather than individual convective events. Abrupt increases in the dipole moment from flux expulsion can produce skewed trends on timescales of tens of kyr.

    more » « less
  4. 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.

    more » « less
  5. Abstract

    Debate continues on the amount and distribution of radioactive heat producing elements (i.e., U, Th, and K) in the Earth, with estimates for mantle heat production varying by an order of magnitude. Constraints on the bulk‐silicate Earth's (BSE) radiogenic power also places constraints on overall BSE composition. Geoneutrino detection is a direct measure of the Earth's decay rate of Th and U. The geoneutrino signal has contributions from the local (40%) and global (35%) continental lithosphere and the underlying inaccessible mantle (25%). Geophysical models are combined with geochemical data sets to predict the geoneutrino signal at current and future geoneutrino detectors. We propagated uncertainties, both chemical and physical, through Monte Carlo methods. Estimated total signal uncertainties are on the order of20%, proportionally with geophysical and geochemical inputs contributing30% and70%, respectively. We find that estimated signals, calculated using CRUST2.0, CRUST1.0, and LITHO1.0, are within physical uncertainty of each other, suggesting that the choice of underlying geophysical model will not change results significantly, but will shift the central value by up to15%. Similarly, we see no significant difference between calculated layer abundances and bulk crustal heat production when using these geophysical models. The bulk crustal heat production is calculated as 7  2 TW, which includes an increase of 1 TW in uncertainty relative to previous studies. Combination of our predicted lithospheric signal with measured signals yield an estimated BSE heat production of 21.5  10.4 TW. Future improvements, including uncertainty attribution and near‐field modeling, are discussed.

    more » « less