skip to main content

Title: The Lithophile Element Budget of Earth's Core

The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO2, and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

more » « less
Award ID(s):
1947614 1725349
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. Comparing compositional models of the terrestrial planets provides insights into physicochemical processes that produced planet-scale similarities and differences. The widely accepted compositional model for Mars assumes Mn and more refractory elements are in CI chondrite proportions in the planet, including Fe, Mg, and Si, which along with O make up >90% of the mass of Mars. However, recent improvements in our understandings on the composition of the solar photosphere and meteorites challenge the use of CI chondrite as an analog of Mars. Here we present an alternative model composition for Mars that avoids such an assumption and is based on data from Martian meteorites and spacecraft observations. Our modeling method was previously applied to predict the Earth’s composition. The model establishes the absolute abundances of refractory lithophile elements in the bulk silicate Mars (BSM) at 2.26 times higher than that in CI carbonaceous chondrites. Relative to this chondritic composition, Mars has a systematic depletion in moderately volatile lithophile elements as a function of their condensation temperatures. Given this finding, we constrain the abundances of siderophile and chalcophile elements in the bulkMars and its core. The Martian volatility trend is consistent with <7 wt% S in its core, which is significantly lower than that assumed in most core models (i.e., >10 wt% S). Furthermore, the occurrence of ringwoodite at the Martian core-mantle boundary might have contributed to the partitioning of O and H into the Martian core. 
    more » « less
  2. Abstract The Bishop Tuff (BT), erupted from the Long Valley caldera in California, displays two types of geochemical gradients with temperature: one is related to magma mixing, whereas the other is found in the high-SiO2 rhyolite portion of the Bishop Tuff and is characterized by twofold or lower concentration variations in minor and trace elements that are strongly correlated with temperature. It is proposed that the latter zonation, which preceded phenocryst growth, developed as a result of mineral–melt partitioning between interstitial melt and surrounding crystals in a parental mush, from which variable melt fractions were segregated. To test this hypothesis, trends of increasing vs decreasing element concentrations with temperature (as a proxy for melt fraction), obtained from published data on single-clast pumice samples from the high-SiO2 rhyolite portion of the Bishop Tuff, were used to infer their relative degrees of incompatibility vs compatibility between crystals and melt in the parental mush. Relative compatibility values (RCVi) for all elements i, defined as the concentration slope with temperature divided by average concentration, are shown to be linearly correlated with their respective bulk partition coefficients (bulk Di). Mineral–melt partition coefficients from the literature were used to constrain the average stoichiometry of the crystallization/melting reaction in the parental mush: 32 % quartz + 34 % plagioclase + 31 % K-feldspar + 1·60 % biotite + 0·42 % titanomagnetite + 0·34 % ilmenite + 0·093 % allanite + 0·024 % zircon + 0·025 % apatite = 100 % liquid. The proportions of tectosilicates in the reaction (i.e. location of eutectic) are consistent with depths of melt segregation of ~400–550 MPa and an activity of H2O of ~0·4–0·6. Temperatures of <770–780 °C are constrained by allanite in the reaction. Evidence that a fluid phase was present in the parental mush is seen in the decreasing versus increasing H2O and CO2 contents with temperature in the segregated interstitial melt that formed the high-SiO2 rhyolite portion of the Bishop Tuff. The presence of an excess fluid phase, which strongly partitions CO2 relative to the melt, is required to explain the compatible behavior of CO2, whereas the fluid abundance must have been low to explain the incompatible behavior of H2O. Calculated degassing paths for interstitial melts, which segregated from the parental mush and ascended to shallower depths to grow phenocrysts, match published volatile analyses in quartz-hosted melt inclusions and constrain fluid abundances in the mush to be ≤1 wt%. The source of volatiles in the parental mush, irrespective of whether it formed by crystallization or partial melting, must have been primarily from associated basalts, as granitoid crust is too volatile poor. Approximately twice as much basalt as rhyolite is needed to provide the requisite volatiles. The determination of bulk Di for several elements gives the bulk composition of the parental leucogranitic mush and shows that it is distinct from Mesozoic Sierran arc granitoids, as expected. Collectively, the results from this study provide new constraints for models of the complex, multi-stage processes throughout the Plio-Quaternary, involving both mantle-derived basalt and pre-existing crust, that led to the origin of the parental body to the Bishop Tuff. 
    more » « less
  3. Abstract

    Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software “VapoRock” combines the MELTS liquid model with gas-species properties from multiple thermochemistry tables. VapoRock calculates the partial pressures of 34 gaseous species in equilibrium with magmatic liquid in the system Si–Mg–Fe–Al–Ca–Na–K–Ti–Cr–O at desired temperatures and oxygen fugacities (fO2, or partial pressure of O2). Comparison with experiments shows that pressures and melt-oxide activities (which vary over many orders of magnitude) are reproduced to within a factor of ∼3, consistent with measurement uncertainties. We also benchmark the model against a wide selection of igneous rock compositions including bulk silicate Earth, predicting elemental vapor abundances that are comparable to (Na, Ca, and Al) or more realistic than (K, Si, Mg, Fe, and Ti) those of the closed-source MAGMA code (with maximum deviations by factors of 10–300 for K and Si). Vapor abundances depend critically on the activities of liquid components. The MELTS model underpinning VapoRock was calibrated and extensively tested on natural igneous liquids. In contrast, MAGMA’s liquid model assumes ideal mixtures of a limited set of chemically simplified pseudospecies, which only roughly approximates the nonideal compositional interactions typical of many-component natural silicate melts. Finally, we explore how relative abundances of SiO and SiO2provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constrainingfO2in outgassed exoplanetary mantles.

    more » « less
  4. Context . The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R ⊕ and 3.0 M ⊕ . It is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. Aims . To prepare for future studies, we aim to thoroughly characterise the planetary system with new accurate and precise data collected with state-of-the-art photometers from space and spectrometers and interferometers from the ground. Methods . We collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X at the 8.1 m Gemini North telescope and CARMENES at the 3.5 m Calar Alto telescope, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed near-infrared interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes (AstroLAB, LCOGT, OSN, TJO). This extraordinary and rich data set was the input for our comprehensive analysis. Results . From interferometry, we measure a limb-darkened disc angular size of the star Gl 486 at θ LDD = 0.390 ± 0.018 mas. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius R * = 0.339 ± 0.015 R ⊕ . We also measure a stellar rotation period at P rot = 49.9 ± 5.5 days, an upper limit to its XUV (5-920 A) flux informed by new Hubble /STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Moreover, we imposed restrictive constraints on the presence of additional components, either stellar or sub-stellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R p = 1.343 −0.062 +0.063 R ⊕ and M p = 3.00 −0.12 +0.13 M ⊕ , with relative uncertainties of the planet radius and mass of 4.7% and 4.2%, respectively. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope ( Webb ) observations. 
    more » « less
  5. Abstract

    High pressure and temperature experiments were carried out on the oxide mixtures corresponding to the bridgmanite stoichiometry under the hydrous shallow lower mantle conditions (24–25 GPa and 1673–1873 K with 5–10 wt. % of water in the starting material). Oxide mixtures investigated correspond to MgSiO3, (Mg, Fe)SiO3, (Mg, Al, Si)O3, and (Mg, Fe, Al, Si)O3. Melting was observed in all runs. Partitioning of various elements, including Mg, Fe, Si, and H is investigated. Melting under hydrous lower mantle conditions leads to increased (Mg + Fe)O/SiO2in the melt compared to the residual solids. The residual solids often contain a large amount of stishovite, and the melt contains higher (Mg,Fe)O/SiO2ratio than the initial material. (Mg + Fe)O‐rich hydrous melt could explain the low‐velocity anomalies observed in the shallow lower mantle and a large amount of stishovite in the residual solid may be responsible for the scattering of seismic waves in the mid‐lower mantle and may explain the “stishovite paradox. Since stishovite‐rich materials are formed only when silica‐rich source rock (MORB) is melted (not a typical peridotitic rock [bulk silicate Earth]), seismic scattering in the lower mantle provides a clue on the circulation of subducted MORB materials. To estimate hydrogen content, we use a new method of estimating the water content of unquenchable melts, and also propose a new interpretation of the significance of superhydrous phase B inclusions in bridgmanite. The results provide revised values of water partitioning between solid minerals and hydrous melts that are substantially higher than previous estimates.

    more » « less