skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. Light element alloying in iron is required to explain density deficit and seismic wave velocities in Earth’s core. However, the light element composition of the Earth’s core seems hard to constrain as nearly all light element alloying would reduce the density and sound velocity (elastic moduli). The alloying light elements include oxidizing elements like oxygen and sulfur and reducing elements like hydrogen and carbon, yet their chemical effects in the alloy system are less discussed. Moreover, Fe-X-ray Absorption Near Edge Structure (Fe-XANES) fingerprints have been studied for silicate materials with ferrous and ferric ions, while not many X-ray absorption spectroscopy (XAS) studies have focused on iron alloys, especially at high pressures. To investigate the bonding nature of iron alloys in planetary interiors, we presented X-ray absorption spectroscopy of iron–nitrogen and iron–carbon alloys at high pressures up to 50 GPa. Together with existing literature on iron–carbon, –hydrogen alloys, we analyzed their edge positions and found no significant difference in the degree of oxidation among these alloys. Pressure effects on edge positions were also found negligible. Our theoretical simulation of the valence state of iron, alloyed with S, C, O, N, and P also showed nearly unchanged behavior under pressures up to 300 GPa. This finding indicates that the high pressure bonding of iron alloyed with light elements closely resembles bonding at the ambient conditions. We suggest that the chemical properties of light elements constrain which ones can coexist within iron alloys. 
    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. Abstract Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe‐C‐N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal‐structured Fe7(N0.75C0.25)3up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First‐principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe7(N0.75C0.25)3. Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe‐C‐N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)‐rich planetary cores, thereby influencing the dynamics of such cores. 
    more » « less