skip to main content


Title: Earth’s core could be the largest terrestrial carbon reservoir
Abstract

Evaluating carbon’s candidacy as a light element in the Earth’s core is critical to constrain the budget and planet-scale distribution of this life-essential element. Here we use first principles molecular dynamics simulations to estimate the density and compressional wave velocity of liquid iron-carbon alloys with ~4-9 wt.% carbon at 0-360 gigapascals and 4000-7000 kelvin. We find that for an iron-carbon binary system, ~1-4 wt.% carbon can explain seismological compressional wave velocities. However, this is incompatible with the ~5-7 wt.% carbon that we find is required to explain the core’s density deficit. When we consider a ternary system including iron, carbon and another light element combined with additional constraints from iron meteorites and the density discontinuity at the inner-core boundary, we find that a carbon content of the outer core of 0.3-2.0 wt.%, is able to satisfy both properties. This could make the outer core the largest reservoir of terrestrial carbon.

 
more » « less
Award ID(s):
1753125
NSF-PAR ID:
10383861
Author(s) / Creator(s):
; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Communications Earth & Environment
Volume:
2
Issue:
1
ISSN:
2662-4435
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon‐rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. Here, we report the density and elasticity of Fe90Ni10‐3 wt.% C and Fe90Ni10‐5 wt.% C liquid alloys using synchrotron‐based X‐ray absorption experiments and first‐principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe90Ni10liquid by ∼2.9–3.1% at 2 GPa, and ∼3.4–3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe‐Ni‐C liquids are similar to or slightly higher than those of Fe‐Ni liquids. Thus, the calculated compressional velocities (vp) of Fe‐Ni‐C liquids are higher than that of pure Fe‐Ni alloy, promoting carbon as a possible candidate to explain the elevatedvpin the Earth's outer core. However, the values and slopes of both density andvpof the studied two Fe‐Ni‐C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The highvpof Fe‐Ni‐C liquids does not match the presumptivevpof the lunar outer core well, indicating that carbon is less likely to be its dominant light element.

     
    more » « less
  2. Using first principles molecular dynamic simulations, we explore the effects of nitrogen (N) on the density and sound velocity of liquid iron and evaluate its potential as a light element in the Earth’s outer core. Our results suggest that Fe-N melt cannot simultaneously explain the density and seismic velocity of the Earth’s outer core. Although ~2.0 wt.% N can explain the bulk sound velocity of the outer core, such N content only lowers the density of liquid Fe by ~3%. Matching both the velocity and density by the other light elements limits the N in the core to ≪2.0 wt.%. Our finding suggests that nitrogen is a minor to trace element in the Earth’s core and is consistent with the geochemical mass balance with terrestrial abundance of N and alloy-silicate partitioning data, which suggest that there cannot be significant N in the core. 
    more » « less
  3. Abstract

    Knowledge of the sound velocity of core materials is essential to explain the observed anomalously low shear wave velocity (VS) and high Poisson’s ratio (σ) in the solid inner core. To date, neitherVSnorσof Fe and Fe-Si alloy have been measured under core conditions. Here, we presentVSandσderived from direct measurements of the compressional wave velocity, bulk sound velocity, and density of Fe and Fe-8.6 wt%Si up to ~230 GPa and ~5400 K. The new data show that neither the effect of temperature nor incorporation of Si would be sufficient to explain the observed lowVSand highσof the inner core. A possible solution would add carbon (C) into the solid inner core that could further decreaseVSand increaseσ. However, the physical property-based Fe-Si-C core models seemingly conflict with the partitioning behavior of Si and C between liquid and solid Fe.

     
    more » « less
  4. null (Ed.)
    Understanding the effect of carbon on the density of hcp (hexagonal-close-packed) Fe-C alloys is essential for modeling the carbon content in the Earth’s inner core. Previous studies have focused on the equations of state of iron carbides that may not be applicable to the solid inner core that may incorporate carbon as dissolved carbon in metallic iron. Carbon substitution in hcp-Fe and its effect on the density have never been experimentally studied. We investigated the compression behavior of Fe-C alloys with 0.31 and 1.37 wt % carbon, along with pure iron as a reference, by in-situ X-ray diffraction measurements up to 135 GPa for pure Fe, and 87 GPa for Fe-0.31C and 109 GPa for Fe-1.37C. The results show that the incorporation of carbon in hcp-Fe leads to the expansion of the lattice, contrary to the known effect in body-centered cubic (bcc)-Fe, suggesting a change in the substitution mechanism or local environment. The data on axial compressibility suggest that increasing carbon content could enhance seismic anisotropy in the Earth’s inner core. The new thermoelastic parameters allow us to develop a thermoelastic model to estimate the carbon content in the inner core when carbon is incorporated as dissolved carbon hcp-Fe. The required carbon contents to explain the density deficit of Earth’s inner core are 1.30 and 0.43 wt % at inner core boundary temperatures of 5000 K and 7000 K, respectively. 
    more » « less
  5. 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