Seismic observations suggest that the uppermost region of Earth’s liquid outer core is buoyant, with slower velocities than the bulk outer core. One possible mechanism for the formation of a stably stratified layer is immiscibility in molten iron alloy systems, which has yet to be demonstrated at core pressures. We find immisci- bility between liquid Fe-Si and Fe-Si-O persisting to at least 140 GPa through a combination of laser-heated diamond-anvil cell experiments and first-principles molecular dynamics simulations. High-pressure immiscibility in the Fe-Si-O system may explain a stratified layer atop the outer core, complicate differentiation and evolution of the deep Earth, and affect the structure and intensity of Earth’s magnetic field. Our results support silicon and oxy- gen as coexisting light elements in the core and suggest that SiO2 does not crystallize out of molten Fe-Si-O at the core-mantle boundary.
more »
« less
Thermal conductivity of Fe-Si alloys and thermal stratification in Earth’s core
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
- Award ID(s):
- 1901801
- NSF-PAR ID:
- 10315202
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 119
- Issue:
- 1
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases
δ ‐AlOOH (delta phase), MgSiO2(OH)2(phase H), andε ‐FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core‐mantle boundary. We have constrained the equation of state, including the effects of a spin crossover in the Fe3+atoms, of (Al, Fe)‐phase H: Al0.84Fe3+0.07Mg0.02Si0.06OOH, using powder X‐ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)‐phase H andδ ‐(Al, Fe)OOH. The changes in spin state of Fe3+in (Al, Fe)‐phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared withδ ‐(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)‐phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)‐phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid‐lower mantle (∼1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)‐phase H could reproduce some of the seismic signatures of large, low seismic‐velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ ‐AlOOH)‐(MgSiO2(OH)2)‐(ε ‐FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle. -
more » « less
Abstract Silicate liquids are important agents of thermal evolution, yet their thermal conductivity is largely unknown. Here, we determine the thermal conductivity of a silicate liquid by combining the Green‐Kubo method with a machine learning potential of
ab initio quality over the entire pressure regime of the mantle. We find that the thermal conductivity of MgSiO3liquid is 1.1 W m−1 K−1at the 1 bar melting point, and 4.0 W m−1 K−1at core‐mantle boundary conditions. The thermal conductivity increases with compression, while remaining nearly constant on isochoric heating. The pressure dependence arises from the increasing bulk modulus on compression, and the weak temperature dependence arises from the saturation of the phonon mean free path due to structural disorder. The thermal conductivity of silicate liquids is less than that of ambient mantle, a contrast that may be important for understanding melt generation, and heat flux from the core. -
more » « less
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
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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2119001119
