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.
- NSF-PAR ID:
- 10089007
- Date Published:
- Journal Name:
- Geophysical research letters
- Volume:
- 46
- ISSN:
- 0094-8276
- Page Range / eLocation ID:
- 89-98
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Density of Fe‐Ni‐C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon
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 (
v p ) of Fe‐Ni‐C liquids are higher than that of pure Fe‐Ni alloy, promoting carbon as a possible candidate to explain the elevatedv p in the Earth's outer core. However, the values and slopes of both density andv p of 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 highv p of Fe‐Ni‐C liquids does not match the presumptivev p of the lunar outer core well, indicating that carbon is less likely to be its dominant light element. -
Abstract Hexagonal close‐packed (hcp) structured Fe‐Ni alloy is believed to be the dominant phase in the Earth's inner core. This phase is expected to contain 4%–5% light elements, such as Si and H. While the effects of individual light element candidates on the equation of state (EoS) of the hcp Fe metal have been studied, their combined effects remain largely unexplored. In this study, we report the equations of state for two hcp‐structured Fe‐Si‐H alloys, namely Fe0.83Si0.17H0.07and Fe0.83Si0.17H0.46, using synchrotron X‐ray diffraction measurements up to 125 GPa at 300 K. These alloys were synthesized by cold compression of Fe‐9wt%Si in either pure H2or Ar‐H2mixture medium in diamond‐anvil cells. The volume increase caused by a H atom in hcp Fe‐Si‐H alloys is approximately eight times greater than that by a Si atom. We used the improved data set to develop a composition‐dependent EoS that covers a wide range of compositions. Our calculated density and bulk sound velocity of hcp Fe‐Si‐H alloys suggest a large trade‐off between Si and H contents in fitting the seismic properties of the inner core. Combining our new EoS with geophysical and geochemical constraints, we propose 1.6–3 wt% Si and 0.15–0.6 wt% H in the Earth's inner core.
-
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
-
Abstract Knowledge of the sound velocity of core materials is essential to explain the observed anomalously low shear wave velocity (
V S) and high Poisson’s ratio (σ ) in the solid inner core. To date, neitherV Snorσ of Fe and Fe-Si alloy have been measured under core conditions. Here, we presentV Sandσ 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 lowV Sand highσ of the inner core. A possible solution would add carbon (C) into the solid inner core that could further decreaseV Sand 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.