Geochemical Perspectives Letters
Geochemical
Perspectives Letters
Geochemical
Perspectives
  • Submit here
  • Track your paper
  • For authors
  • e-Alerts
  • Home
  • About
    • About the journal
    • Editorial Board
    • Publication Policy
    • Publication Ethics
  • Submission & Review
    • Copyright & Permissions
    • Information for Authors
    • Information for Reviewers
  • Current issue
  • All issues
  • Submit
Select Page Menu

by admin | Sep 7, 2021 | mainpost, vol18 | 0 comments

C.R.M. Jackson, E. Cottrell, Z. Du, N.R. Bennett, Y. Fei

18

2122

7

March

2021

15

July

2021

7

September

2021

37

42

0

Next article >> << Previous article

High pressure redistribution of nitrogen and sulfur during planetary stratification

C.R.M. Jackson1,

1Department of Earth and Environmental Sciences, Tulane University, 6823 St. Charles Avenue, New Orleans, LA 70118, USA

E. Cottrell2,

2Department of Minerals Sciences, National Museum of Natural History, Smithsonian Institution, 10th and Constitution Ave NW, Washington D.C. 20560, USA

Z. Du3,4,

3State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou, China
4CAS Center for Excellence in Deep Earth Sciences, Guangzhou, 510640, China

N.R. Bennett5,

5Department of Earth Sciences, University of Toronto, 22 Russell St, Toronto, ON M5S 3B1, Canada

Y. Fei6

6Earth and Planets Laboratory, Carnegie Institution for Science, 5241 Broad Branch Rd NW, Washington D.C. 20015, USA

Affiliations | Corresponding Author | Cite as | Funding information

C.R.M. Jackson
Email: cjackson2@tulane.edu

1Department of Earth and Environmental Sciences, Tulane University, 6823 St. Charles Avenue, New Orleans, LA 70118, USA
2Department of Minerals Sciences, National Museum of Natural History, Smithsonian Institution, 10th and Constitution Ave NW, Washington D.C. 20560, USA
3State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou, China
4CAS Center for Excellence in Deep Earth Sciences, Guangzhou, 510640, China
5Department of Earth Sciences, University of Toronto, 22 Russell St, Toronto, ON M5S 3B1, Canada
6Earth and Planets Laboratory, Carnegie Institution for Science, 5241 Broad Branch Rd NW, Washington D.C. 20015, USA

Jackson, C.R.M., Cottrell, E., Du, Z., Bennett, N.R., Fei, Y. (2021) High pressure redistribution of nitrogen and sulfur during planetary stratification. Geochem. Persp. Let. 18, 37–42.

Smithsonian; NSF EAR 1725315; Tulane University, NASA.

Geochemical Perspectives Letters v18 | doi: 10.7185/geochemlet.2122
Received 7 March 2021 | Accepted 15 July 2021 | Published 7 September 2021

Copyright © 2021 The Authors

Published by the European Association of Geochemistry
under Creative Commons License CC BY 4.0

Keywords: nitrogen, sulfur, core formation, accretion

PDF PDF+SI
  • Share this article

  • Article views:
    1,118

    Cumulative count of HTML views and PDF downloads.

  • Download Citation
  • Rights & Permissions


top

Abstract

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information

Nitrogen is essential to life, and yet is also the most depleted element in the Earth relative to gas-rich chondrites. A key expression of Earth’s N depletion is its elevated sulfur-nitrogen (S/N) ratio. Primordial stratification into a core, mantle, and atmosphere is the largest mass transfer process that terrestrial planets experience, but the data required to evaluate how S/N ratios respond to primordial stratification of Earth-sized planets do not exist. We report new metal-silicate partitioning experiments on N up to 26 GPa and 3437 K. Our data indicate that nitrogen becomes more siderophile with increasing pressure and less siderophile with increasing carbon and nickel in the metal phase. We apply our new experiments with literature data for S partitioning to a core formation-primordial atmosphere degassing model. Our model demonstrates that the S/N ratio of the observable Earth can be set during primordial stratification under the same extreme P-T conditions that satisfy refractory siderophile element budgets while also yielding a bulk planet with S contents near that estimated from Earth’s volatility trend.

Figures and Tables

Figure 1 Concentrations of volatile elements (H, C, S, and N) in Earth’s observable reservoirs normalised to CI chondrites. Nitrogen appears depleted in Earth compared to other volatiles. Data are from Hirschmann (2016) and Wasson and Kallemeyn (1988).

Figure 2 Controls on nitrogen partitioning. (a) Oxidising conditions make nitrogen more siderophile. Literature data are recalculated to 1973 K and 1 GPa. (b) Nickel and (c) Sulfur make nitrogen less siderophile. (d) Hotter temperatures make nitrogen less siderophile. Extrapolation of temperature effect at 1 GPa underpredicts values compared to values measured for higher pressure LH-DAC experiments. (e) A 1∶1 comparison of observations versus predictions from Equation 1. The R2 for the literature fit is 0.76. (f) Application of Equation 1 to mantle liquidus geotherm at ΔIW−2. Nitrogen becomes more siderophile with depth.

Figure 3 Models for redistribution of N and S in response to planetary differentiation. (a) An example of a model iteration of S/NOE as function of core formation pressure and atmospheric fO2. Solid lines are model solutions for a given atmospheric fO2 (log unit spacing). Horizontal dotted lines delineate the estimated range of S/NOE for Earth (Hirschmann, 2016). Observable Earth is the sum of the atmosphere and magma ocean. (b) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (c) An example of a model iteration of S/NOE.MO as function of core formation pressure and atmospheric fO2. Horizontal dotted lines delineate the estimated range of S/NOE for Earth. Observable Earth is only the magma ocean. (d) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (e) Predicted bulk planet S as a function of core formation pressure. Dashed lines are 1 sigma that propagate uncertainties from the observable Earth S and parameterised values. The dotted line is the estimate of observable Earth S from the volatility trend (5600 ppm; Dreibus and Palme, 1996).

Table 1 Experimental run conditions and chemical analyses used in parameterisation.

Figure 1 Figure 2 Figure 3 Table 1

View all figures and tables





top

Introduction

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


Nitrogen is the most depleted element in the observable Earth when compared to CI chondrites, materials that most faithfully record the bulk composition of the solar system (observable Earth is the sum of the atmosphere, oceans, crust, and mantle) (e.g., Marty, 2012

Marty, B. (2012) The origins and concentrations of water, carbon, nitrogen and noble gases on earth. Earth and Planetary Science Letters 313–314, 56–66.

; Halliday, 2013

Halliday, A.N. (2013) The origins of volatiles in the terrestrial planets. Geochimica et Cosmochimica Acta 105, 146–171.

). The depletion of N is a fundamental expression of the integrated processes that determine the volatile budget of planets (Fig. 1).


Figure 1 Concentrations of volatile elements (H, C, S, and N) in Earth’s observable reservoirs normalised to CI chondrites. Nitrogen appears depleted in Earth compared to other volatiles. Data are from Hirschmann (2016)

Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.

and Wasson and Kallemeyn (1988)

Wasson, J.T., Kallemeyn, G.W. (1988) Compositions of chondrites. Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences 325, 535–544.

.
Full size image


There are many early acting processes that can deplete nitrogen and other volatile elements, either in a bulk planet or its observable reservoirs. Chondrites have various volatile element patterns that may reflect formation in regions of the Solar System with fractionated volatiles or the operation of incomplete condensation/evaporation reactions (Bergin et al., 2015

Bergin, E.A., Blake, G.A., Ciesla, F., Hirschmann, M.M., Li, J. (2015) Tracing the ingredients for a habitable earth from interstellar space through planet formation. Proceedings of the National Academy of Sciences 112, 8965–8970.

). With accretion and differentiation, terrestrial bodies form atmospheres, and these atmospheres can be lost (Tucker and Mukhopadhyay, 2014

Tucker, J.M., Mukhopadhyay, S. (2014) Evidence for multiple magma ocean outgassing and atmospheric loss episodes from mantle noble gases. Earth and Planetary Science Letters 393, 254–265.

; Schlichting et al., 2015

Schlichting, H.E., Sari, R., Yalinewich, A. (2015) Atmospheric mass loss during planet formation: The importance of planetesimal impacts. Icarus 247, 81–94.

). Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013

Roskosz, M., Bouhifd, M.A., Jephcoat, A.P., Marty, B., Mysen, B.O. (2013) Nitrogen solubility in molten metal and silicate at high pressure and temperature. Geochimica et Cosmochimica Acta 121, 15–28.

; Kadik et al., 2015

Kadik, A., Koltashev, V., Kryukova, E., Plotnichenko, V., Tsekhonya, T., Kononkova, N. (2015) Solubility of nitrogen, carbon, and hydrogen in FeO–Na2O–Al2O3–SiO2 melt and liquid iron alloy: Influence of oxygen fugacity. Geochemistry International 53, 849–868.

; Li et al., 2016

Li, Y., Marty, B., Shcheka, S., Zimmermann, L., Keppler, H. (2016) Nitrogen isotope fractionation during terrestrial core-mantle separation. Geochemical Perspectives Letters 2, 138–147.

; Dalou et al., 2017

Dalou, C., Hirschmann, M.M., von der Handt, A., Mosenfelder, J., Armstrong, L.S. (2017) Nitrogen and carbon fractionation during core–mantle differentiation at shallow depth. Earth and Planetary Science Letters 458, 141–151.

; Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

,b

Grewal, D.S., Dasgupta, R., Sun, C., Tsuno, K., Costin, G. (2019b) Delivery of carbon, nitrogen, and sulfur to the silicate earth by a giant impact. Science Advances 5, eaau3669.

; Speelmanns et al., 2019

Speelmanns, I.M., Schmidt, M.W., Liebske, C. (2019) The almost lithophile character of nitrogen during core formation. Earth and Planetary Science Letters 510, 186–197.

).

The role of pressure and temperature on core formation chemistry become larger for larger planetary bodies, as indicated by the refractory siderophile element concentrations observed for Earth, Mars, and Vesta (Righter and Drake, 1996

Righter, K., Drake, M.J. (1996) Core formation in Earth’s moon, Mars, and Vesta. Icarus 124, 513–529.

). Studies agree that increasing temperature makes nitrogen less siderophile (Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

; Speelmanns et al., 2019

Speelmanns, I.M., Schmidt, M.W., Liebske, C. (2019) The almost lithophile character of nitrogen during core formation. Earth and Planetary Science Letters 510, 186–197.

). Extrapolating the temperature effect on metal-silicate partition coefficients for nitrogen (, atomic) to the average core formation temperature for Earth implies lithophile behaviour (e.g., Speelmanns et al., 2019

Speelmanns, I.M., Schmidt, M.W., Liebske, C. (2019) The almost lithophile character of nitrogen during core formation. Earth and Planetary Science Letters 510, 186–197.

). In the absence of a large effect of pressure, core formation would apparently have little ability to modify the observable budget of nitrogen during the accretion of larger rocky worlds, such as Earth.

Previously published high pressure experiments suggest that pressure favours N incorporation into cores (Roskosz et al., 2013

Roskosz, M., Bouhifd, M.A., Jephcoat, A.P., Marty, B., Mysen, B.O. (2013) Nitrogen solubility in molten metal and silicate at high pressure and temperature. Geochimica et Cosmochimica Acta 121, 15–28.

, Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

); however, the pressure effect cannot be confidently resolved to be different from zero (Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

). Towards this end, we report values from experiments conducted up to 26 GPa and 3437 K, using a laser heated diamond-anvil cell (LH-DAC). We supplement our LH-DAC data, with a systematic series of piston cylinder experiments to enable a robust parameterisation.

top

Methods

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


We conducted lower P-T partitioning experiments (1773–2413 K, 0.95–2.38 GPa) using a piston cylinder (PC) and higher P-T experiments (3046–3437 K, 23.2–25.6 GPa) using a LH-DAC (Tables 1, S-1). Experimental fO2 conditions ranged from ΔIW−6.6 to ΔIW−1.0 (log unit deviations from the iron-wüstite buffer). We quantified the composition of reacted metal-silicate pairs using field emission electron microprobe analysis (Table 1). An example backscatter image is provided in Figure S-1, along with a typical time-temperature path for LH-DAC experiments. We also compare major element partitioning data of LH-DAC experiments to literature data in Figure S-2 as an evaluation of data quality. Table S-2 reports starting materials compositions. Further details on methods are provided in the Supplementary Information.

Table 1 Experimental run conditions and chemical analyses used in parameterisation.
ExperimentNoteSeriesT, KT, ±P, GPaP, ±log(DN), atomiclog(γFe)Ni, atomicS, atomicC, atomicN, atomicΔIW
PC_N_EXP1P, T, S, Ni, N, fO2, IP1973100.950.051.40−0.080.0000.0000.1590.024−1.59
PC_N_EXP2low S metalS1973100.950.051.44−0.060.0000.0140.1330.038−1.62
PC_N_EXP2high S metalS1973100.950.051.170.180.0000.3720.0250.020−1.80
PC_N_EXP3fO21973100.950.050.01−0.010.0000.0000.0530.002−6.56
PC_N_EXP4fO21973100.950.05−0.68−0.040.0000.0000.1090.003−5.28
PC_N_EXP5fO21973100.950.05−0.59−0.060.0000.0000.1360.007−4.65
PC_N_EXP6P1973102.380.051.36−0.050.0000.0000.1190.106−1.61
PC_N_EXP7high S metalS1973102.380.051.43−0.040.0000.1090.0910.086−1.57
PC_N_EXP8fO21973100.950.050.97−0.070.0000.0000.1420.031−2.50
PC_N_EXP9T2373101.050.050.77−0.060.0000.0000.1510.019−1.64
PC_N_EXP10T2403101.050.050.46−0.070.0000.0000.1600.013−1.62
PC_N_EXP11low S metalS2413101.050.051.05−0.060.0000.0120.1430.033−1.69
PC_N_EXP12low S metalS1973101.660.051.37−0.050.0000.0300.1110.074−1.69
PC_N_EXP12high S metalS1973101.660.051.200.110.0000.3010.0260.050−1.78
PC_N_EXP13Ni1973101.000.050.80−0.050.2250.0000.1170.019−1.57
PC_N_EXP14T, N1773100.950.051.56−0.070.0000.0000.1380.039−1.70
PC_N_EXP15Ni1973100.950.050.84−0.070.0120.0000.1430.032−2.48
PC_N_EXP17T, N2173101.050.050.90−0.060.0000.0000.1470.029−1.64
PC_N_EXP18T, N1873101.050.051.47−0.070.0000.0000.1460.033−1.60
PC_N_EXP19IP1973100.950.050.92−0.080.0000.0000.1570.030−1.75
PC_N_EXP22IP1973100.950.050.83−0.060.0000.0000.1370.042−1.71
PC_N_EXP23IP1973100.950.050.82−0.070.0000.0000.1490.036−1.61
PC_N_EXP24Ni1973100.950.051.11−0.070.1040.0000.1480.009−1.85
DAC_N_EXP1 spot 2P31424623.802.351.55−0.060.2560.0000.1600.135−1.18
DAC_N_EXP1 spot 3P343725525.603.901.78−0.100.2150.0000.2190.054−0.96
DAC_N_EXP1 spot 4P30461323.202.101.59−0.060.3030.0000.1540.118−1.15


Series: P, pressure; T, temperature; S, sulfur; Ni, nickel; N-C, nitrogen-carbon; fO2, oxygen fugacity; IP, ionic porosity; Ni, S, C, N, atomic fraction of metal phase.



top

Results

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


We completed independent series of experiments to quantify the P-T-X controls on values (Fig. 2, Table 1). Individual series correlations indicate values depend strongly on fO2 (R2 = 0.99, p value = 0.005, n = 4), temperature (R2 = 0.97, p value = 0.001, n = 6), and pressure (R2 = 0.98, p value < 0.001, n = 9) over other conditions relevant to magma oceans (Fig. 2a,d), whereas , , , and the ionic porosity of silicate liquid have significant (p values < 0.05) but minor effects (Figs. 2b,c, S-4, S-5). Our fO2 and temperature findings accord with previous work (Dalou et al., 2017

Dalou, C., Hirschmann, M.M., von der Handt, A., Mosenfelder, J., Armstrong, L.S. (2017) Nitrogen and carbon fractionation during core–mantle differentiation at shallow depth. Earth and Planetary Science Letters 458, 141–151.

; Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

; Speelmanns et al., 2019

Speelmanns, I.M., Schmidt, M.W., Liebske, C. (2019) The almost lithophile character of nitrogen during core formation. Earth and Planetary Science Letters 510, 186–197.

). Our lower pressure experiments indicate that N-C interactions in Fe alloy make both elements less siderophile and yield an epsilon value () of 7.6 ± 1.4 (R2 = 0.90, p value = 0.004, n = 6; Fig. S-4), qualitatively consistent with the Steelmaking Data Sourcebook (1988)

Steelmaking Data Sourcebook (1988) The Japan Society for the Promotion of Science: The 19th Committee on Steelmaking. Gordon and Breach Science Publishers, New York.

and values collected under C-free conditions (Grewal et al., 2021

Grewal, D.S., Dasgupta, R., Hough, T., Farnell, A. (2021) Rates of protoplanetary accretion and differentiation set nitrogen budget of rocky planets. Nature Geoscience 14, 369–376.

). Our value calculations follow the approach of Ma (2001)

Ma, Z. (2001) Thermodynamic description for concentrated metallic solutions using interaction parameters. Metallurgical and Materials Transactions B 32, 87–103.

.

Our LH-DAC experiments consistently demonstrated siderophile behaviour (Fig. 2d) despite their high temperature, and this requires a large, positive pressure effect. These high P-T experiments contain C and Ni, and were relatively oxidising. Quantifying the pressure effect therefore requires simultaneous consideration of these other parameters, which is enabled by our PC experiments (Equation 1, see below). Our determination of a positive effect of pressure is consistent with previous work, but the larger range of P we investigated results in more precise determination of the P-T coefficients in Equations 1 and 2 (c.f., Roskosz et al., 2013

Roskosz, M., Bouhifd, M.A., Jephcoat, A.P., Marty, B., Mysen, B.O. (2013) Nitrogen solubility in molten metal and silicate at high pressure and temperature. Geochimica et Cosmochimica Acta 121, 15–28.

, Grewal et al., 2019a

Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.

).


Figure 2 Controls on nitrogen partitioning. (a) Oxidising conditions make nitrogen more siderophile. Literature data are recalculated to 1973 K and 1 GPa. (b) Nickel and (c) Sulfur make nitrogen less siderophile. (d) Hotter temperatures make nitrogen less siderophile. Extrapolation of temperature effect at 1 GPa underpredicts values compared to values measured for higher pressure LH-DAC experiments. (e) A 1∶1 comparison of observations versus predictions from Equation 1. The R2 for the literature fit is 0.76. (f) Application of Equation 1 to mantle liquidus geotherm at ΔIW−2. Nitrogen becomes more siderophile with depth.
Full size image


To parameterise our data, we first recalculate experiments to a carbon-free baseline using our newly derived value. We then conduct an equal weight, least squares regression on parameters identified as significant in their individual series (1/T, P/T, , and ΔIW). Our approach yields the following expression (R2 = 0.95, p value < 0.001, n = 22; Fig. 2e):

 Eq. 1



or if the effect of carbon is included in the parameterisation

 Eq. 2



Uncertainties are reported as 1σ. Note that , , and refer to the expanded concentration expression associated with the notation of Ma (2001)

Ma, Z. (2001) Thermodynamic description for concentrated metallic solutions using interaction parameters. Metallurgical and Materials Transactions B 32, 87–103.

, and that positive coefficients indicate a reduction in the value. Application of Equation 1 to a mantle liquidus geotherm at ΔIW−2 indicates a monotonic increase in with depth (Fig. 2f).

We parameterise our data alone because we completed systematic series of experiments to isolate specific effects on partitioning. Inclusion of all published data is accompanied by a large number of free parameters, and prevents resolution of N, S and C effects that are clearly observable in our results (Figs. 2, S-4). Predictions of Equation 1 are compared to literature data in Figures 2e, S-6. Additional details regarding this regression are provided in the Supplementary Information. The covariance matrix for Equation 1 is reported in Table S-3. We also report our fitting of literature data with the effects identified as significant here as Equation S-7.

top

Modelling N Distribution throughout Earth

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


Having established that pressure strongly modulates values, we now calculate the equilibrium distribution of N and S between cores, magma oceans, and atmospheres for a range of plausible stratification conditions. We include S because multiple high pressure studies show agreement and now constrain its metal-silicate partitioning up to extreme P-T conditions. Our parameterisation of values is applied here and is detailed in the Supplementary Information (Fig. S-8).

Mass balance is used to solve for the equilibrium distribution of N and S as required by partitioning and magma solubility constraints (Equation 1; Libourel et al., 2003

Libourel, G., Marty, B., Humbert, F. (2003) Nitrogen solubility in basaltic melt. Part I. Effect of oxygen fugacity. Geochimica et Cosmochimica Acta 67, 4123–4135.

). The mantle is assumed to be completely molten, and the mass fraction of the core is 0.325. Uncertainties are evaluated by varying and according to their covariance matrices and iterating the mass balance model. No S is assumed to be present in the atmosphere given its relatively high solubility in magma. All scenarios assume a bulk planet S/N ratio set by enstatite chondrite (85) and that the magma ocean contains 225 ± 25 ppm S (Hirschmann, 2016

Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.

). We focus on enstatite chondrites because their average S/N ratio is intermediate between gas-rich (29; CI and CM) and ordinary chondrites (404; Wasson and Kallemeyn, 1988

Wasson, J.T., Kallemeyn, G.W. (1988) Compositions of chondrites. Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences 325, 535–544.

). Models that consider S/N ratios for other groups of chondrites are presented in the Supplementary Information (Figs. S-9, S-10).

Our first goal is to quantify how the S/N ratio of observable Earth (S/NOE, 82 ± 29, dotted lines in Fig. 3a,c) varies in response to differentiation conditions. We consider models with core formation fO2 of ΔIW−2 while varying core formation pressure and atmospheric fO2 (Fig. 3). We focus on ΔIW−2 because the FeO content of Earth implies an average fO2 for core formation near this value. Pressures of metal-silicate equilibration range up to 60 GPa, an upper limit in the single stage framework for Earth (e.g., Fischer et al., 2015

Fischer, R.A., Nakajima, Y., Campbell, A.J., Frost, D.J., Harries, D., Langenhorst, F., Miyajima, N., Pollok, K., Rubie, D.C. (2015) High pressure metal–silicate partitioning of Ni, Co, V, Cr, Si, and O. Geochimica et Cosmochimica Acta 167, 177–194.

). Our model takes a single stage approach, and while accretion is a multi-step process, single stage calculations capture average P-T-X conditions that can be readily compared between elements (Siebert et al., 2011

Siebert, J., Corgne, A., Ryerson, F.J. (2011) Systematics of metal–silicate partitioning for many siderophile elements applied to earth’s core formation. Geochimica et Cosmochimica Acta 75, 1451–1489.

). Atmospheric fO2 varies in our model between ΔIW−3 and ΔIW+3. Atmospheric fO2 conditions are a free parameter in light of evidence for the depth dependence of the Fe+3/Fetot ratio of magma in equilibrium with iron (Zhang et al., 2017

Zhang, H., Hirschmann, M., Cottrell, E., Withers, A. (2017) Effect of pressure on Fe3+/Fe ratio in a mafic magma and consequences for magma ocean redox gradients. Geochimica et Cosmochimica Acta 204, 83–103.

; Armstrong et al., 2019

Armstrong, K., Frost, D.J., McCammon, C.A., Rubie, D.C., Ballaran, T.B. (2019) Deep magma ocean formation set the oxidation state of earth’s mantle. Science 365, 903–906.

; Deng etal., 2020

Deng, J., Du, Z., Karki, B., Ghosh D., Lee, K. (2020) A magma ocean origin to divergent redox evolutions of rocky planetary bodies and early atmospheres. Nature Communications 11, 2007.

) and the temperature dependence of Fe+3/Fetot for any fO2 (Sossi et al., 2020

Sossi, P.A., Burnham, A.D., Badro, J., Lanzirotti, A., Newville, M., O’Neill, H.S.C. (2020) Redox state of earth’s magma ocean and its venus-like early atmosphere. Science Advances 6, eabd1387.

).


Figure 3 Models for redistribution of N and S in response to planetary differentiation. (a) An example of a model iteration of S/NOE as function of core formation pressure and atmospheric fO2. Solid lines are model solutions for a given atmospheric fO2 (log unit spacing). Horizontal dotted lines delineate the estimated range of S/NOE for Earth (Hirschmann, 2016

Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.

). Observable Earth is the sum of the atmosphere and magma ocean. (b) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (c) An example of a model iteration of S/NOE.MO as function of core formation pressure and atmospheric fO2. Horizontal dotted lines delineate the estimated range of S/NOE for Earth. Observable Earth is only the magma ocean. (d) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (e) Predicted bulk planet S as a function of core formation pressure. Dashed lines are 1 sigma that propagate uncertainties from the observable Earth S and parameterised values. The dotted line is the estimate of observable Earth S from the volatility trend (5600 ppm; Dreibus and Palme, 1996

Dreibus, G., Palme, H. (1996) Cosmochemical constraints on the sulfur content in the Earth’s core. Geochimica et Cosmochimica Acta 60, 1125–1130.

).
Full size image


Our models demonstrate the strong sensitivity of S/NOE ratios to core formation pressure; at low pressure (<10 GPa), S/NOR ratios are low, but with increasing pressure, S/NOE ratios rise as S becomes less siderophile and N becomes more siderophile. The rise is such that core formation near 30 GPa results in S/NOE ratios that match Earth for more reducing atmospheres. More oxidising atmospheres require higher pressure core formation to satisfy Earth’s S/NOE ratio, as more N remains in the atmosphere, unable to partition into the core (Fig. 3a,b).

It is possible that the atmosphere is partially lost during accretion, leading to preferential depletion of N relative to S. We model this possibility as the end member of the magma ocean being the only contributor of N and S to later-forming observable reservoirs (Fig. 3c,d). In these cases S/NOE,MO (S/N ratio of the magma ocean) also matches Earth near 30 GPa, and model results are relatively insensitive to atmospheric fO2. Cases of intermediate atmospheric loss plot between the two end members for atmospheric contribution considered here (c.f., Fig. 3a–d).

Taken together, our models demonstrate that higher pressure core formation scenarios can satisfy the S/NOE ratio of Earth, assuming a bulk planet S/N ratio similar to enstatite chondrite. Bulk planet S/N ratios closer to gas-rich chondrites (low S/N) require even higher core formation pressures, while bulk planet S/N ratios closer to ordinary chondrites (high S/N) require lower pressure to satisfy the S/NOE ratio constraint (Figs. S-9, S-10). This all serves to highlight the importance of core formation pressure for modulating planetary volatile budgets.

It is significant that higher pressure core formation (>30 GPa) and a core-mantle fO2 of ΔIW−2 can produce S/NOE ratios equal to Earth because these are the same differentiation conditions implied by moderately siderophile element and FeO concentrations in the observable Earth in a single stage framework (e.g., Siebert et al., 2011

Siebert, J., Corgne, A., Ryerson, F.J. (2011) Systematics of metal–silicate partitioning for many siderophile elements applied to earth’s core formation. Geochimica et Cosmochimica Acta 75, 1451–1489.

; Fischer et al., 2015

Fischer, R.A., Nakajima, Y., Campbell, A.J., Frost, D.J., Harries, D., Langenhorst, F., Miyajima, N., Pollok, K., Rubie, D.C. (2015) High pressure metal–silicate partitioning of Ni, Co, V, Cr, Si, and O. Geochimica et Cosmochimica Acta 167, 177–194.

). Importantly, higher pressure core formation also yields a bulk planet with 6650–2050 ppm S (30–60 GPa) (Fig. 3e), and these values compare favourably with estimates of 5600 ppm S for bulk Earth based on the volatility trend (Dreibus and Palme, 1996

Dreibus, G., Palme, H. (1996) Cosmochemical constraints on the sulfur content in the Earth’s core. Geochimica et Cosmochimica Acta 60, 1125–1130.

). The multiple successes of higher pressure models for explaining Earth’s volatile budget are important because they suggest that volatiles are modulated by the same core formation events that modulate refractory elements. This suggestion contrasts with previous hypotheses that smaller, volatile-rich bodies preferentially contributed to Earth’s volatile budget (e.g., Grewal et al., 2019b

Grewal, D.S., Dasgupta, R., Sun, C., Tsuno, K., Costin, G. (2019b) Delivery of carbon, nitrogen, and sulfur to the silicate earth by a giant impact. Science Advances 5, eaau3669.

), decoupling the accretion of life-enabling elements from refractory elements.

It is well established within the Solar System that planetary body size correlates with average core formation pressures (Righter and Drake, 1996

Righter, K., Drake, M.J. (1996) Core formation in Earth’s moon, Mars, and Vesta. Icarus 124, 513–529.

), and our work therefore predicts a direct relationship between planet size and its distribution of volatiles. This link should enable more precise evaluations of exoplanet habitability, worlds for which size remains a central constraint on their geologic evolution.

top

Acknowledgements

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


We acknowledge discussions with Ian Ocampo during a summer internship at NMNH. We additionally acknowledge the support provided by C. Prescher, E. Greenberg, and V. Prakapenka during the experiments conducted at GSECARS and in data reduction. We thank Rob Wardell and Tim Rose for their support of the microscopy work associated with this manuscript. We thank Tim McCoy for discussions regarding the model results.

CRMJ acknowledges support by NSF-EAR grant 1725315, NASA grant 80NSSC21K0377, Smithsonian GVP fellowship, Carnegie Institution for Science (Geophysical Lab) Fellowship, and start up from Tulane University. Experiments were partially supported by NASA grant to YF. ZD thanks the support from Carnegie Fellowship, grants from Chinese Academy of Sciences and State Key Laboratory of Isotope Geochemistry (No. 29Y93301701, 51Y8340107), as well as Strategic Priority Research Program (B) (XDB18030604). Portions of this work were performed at GeoSoilEnviroCARS (The University of Chicago, Sector 13), the Advanced Photon Source, Argonne National Laboratory. GeoSoilEnviroCARS is supported by the National Science Foundation–Earth Sciences (EAR 1128799) and the Department of Energy–GeoSciences (DE-FG02-94ER14466). This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under contract number DE-AC02-06CH11357.

Editor: Cin-Ty Lee

top

References

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information

Armstrong, K., Frost, D.J., McCammon, C.A., Rubie, D.C., Ballaran, T.B. (2019) Deep magma ocean formation set the oxidation state of earth’s mantle. Science 365, 903–906.
Show in context

Atmospheric fO2 conditions are a free parameter in light of evidence for the depth dependence of the Fe+3/Fetot ratio of magma in equilibrium with iron (Zhang et al., 2017; Armstrong et al., 2019; Deng etal., 2020) and the temperature dependence of Fe+3/Fetot for any fO2 (Sossi et al., 2020).
View in article


Bergin, E.A., Blake, G.A., Ciesla, F., Hirschmann, M.M., Li, J. (2015) Tracing the ingredients for a habitable earth from interstellar space through planet formation. Proceedings of the National Academy of Sciences 112, 8965–8970.
Show in context

Chondrites have various volatile element patterns that may reflect formation in regions of the Solar System with fractionated volatiles or the operation of incomplete condensation/evaporation reactions (Bergin et al., 2015).
View in article


Dalou, C., Hirschmann, M.M., von der Handt, A., Mosenfelder, J., Armstrong, L.S. (2017) Nitrogen and carbon fractionation during core–mantle differentiation at shallow depth. Earth and Planetary Science Letters 458, 141–151.
Show in context

Our fO2 and temperature findings accord with previous work (Dalou et al., 2017; Grewal et al., 2019a; Speelmanns et al., 2019).
View in article
Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Deng, J., Du, Z., Karki, B., Ghosh D., Lee, K. (2020) A magma ocean origin to divergent redox evolutions of rocky planetary bodies and early atmospheres. Nature Communications 11, 2007.
Dreibus, G., Palme, H. (1996) Cosmochemical constraints on the sulfur content in the Earth’s core. Geochimica et Cosmochimica Acta 60, 1125–1130.
Show in context

The dotted line is the estimate of observable Earth S from the volatility trend (5600 ppm; Dreibus and Palme, 1996).
View in article
Importantly, higher pressure core formation also yields a bulk planet with 6650–2050 ppm S (30–60 GPa) (Fig. 3e), and these values compare favourably with estimates of 5600 ppm S for bulk Earth based on the volatility trend (Dreibus and Palme, 1996).
View in article


Fischer, R.A., Nakajima, Y., Campbell, A.J., Frost, D.J., Harries, D., Langenhorst, F., Miyajima, N., Pollok, K., Rubie, D.C. (2015) High pressure metal–silicate partitioning of Ni, Co, V, Cr, Si, and O. Geochimica et Cosmochimica Acta 167, 177–194.
Show in context

Pressures of metal-silicate equilibration range up to 60 GPa, an upper limit in the single stage framework for Earth (e.g., Fischer et al., 2015).
View in article
It is significant that higher pressure core formation (>30 GPa) and a core-mantle fO2 of ΔIW−2 can produce S/NOE ratios equal to Earth because these are the same differentiation conditions implied by moderately siderophile element and FeO concentrations in the observable Earth in a single stage framework (e.g., Siebert et al., 2011; Fischer et al., 2015).
View in article


Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., Li, Y., Tsuno, K. (2019a) The fate of nitrogen during core-mantle separation on earth. Geochimica et Cosmochimica Acta 251, 87–115.
Show in context

Studies agree that increasing temperature makes nitrogen less siderophile (Grewal et al., 2019a; Speelmanns et al., 2019).
View in article
Previously published high pressure experiments suggest that pressure favours N incorporation into cores (Roskosz et al., 2013, Grewal et al., 2019a); however, the pressure effect cannot be confidently resolved to be different from zero (Grewal et al., 2019a).
View in article
Previously published high pressure experiments suggest that pressure favours N incorporation into cores (Roskosz et al., 2013, Grewal et al., 2019a); however, the pressure effect cannot be confidently resolved to be different from zero (Grewal et al., 2019a).
View in article
Our fO2 and temperature findings accord with previous work (Dalou et al., 2017; Grewal et al., 2019a; Speelmanns et al., 2019).
View in article
Our determination of a positive effect of pressure is consistent with previous work, but the larger range of P we investigated results in more precise determination of the P-T coefficients in Equations 1 and 2 (c.f., Roskosz et al., 2013, Grewal et al., 2019a).
View in article
Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Grewal, D.S., Dasgupta, R., Sun, C., Tsuno, K., Costin, G. (2019b) Delivery of carbon, nitrogen, and sulfur to the silicate earth by a giant impact. Science Advances 5, eaau3669.
Show in context

This suggestion contrasts with previous hypotheses that smaller, volatile-rich bodies preferentially contributed to Earth’s volatile budget (e.g., Grewal et al., 2019b), decoupling the accretion of life-enabling elements from refractory elements.
View in article
Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Grewal, D.S., Dasgupta, R., Hough, T., Farnell, A. (2021) Rates of protoplanetary accretion and differentiation set nitrogen budget of rocky planets. Nature Geoscience 14, 369–376.
Show in context

Our lower pressure experiments indicate that N-C interactions in Fe alloy make both elements less siderophile and yield an epsilon value () of 7.6 ± 1.4 (R2 = 0.90, p value = 0.004, n = 6; Fig. S-4), qualitatively consistent with the Steelmaking Data Sourcebook (1988) and values collected under C-free conditions (Grewal et al., 2021).
View in article


Halliday, A.N. (2013) The origins of volatiles in the terrestrial planets. Geochimica et Cosmochimica Acta 105, 146–171.
Show in context

Nitrogen is the most depleted element in the observable Earth when compared to CI chondrites, materials that most faithfully record the bulk composition of the solar system (observable Earth is the sum of the atmosphere, oceans, crust, and mantle) (e.g., Marty, 2012; Halliday, 2013).
View in article


Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.
Show in context

All scenarios assume a bulk planet S/N ratio set by enstatite chondrite (85) and that the magma ocean contains 225 ± 25 ppm S (Hirschmann, 2016).
View in article
Horizontal dotted lines delineate the estimated range of S/NOE for Earth (Hirschmann, 2016).
View in article
Data are from Hirschmann (2016) and Wasson and Kallemeyn (1988).
View in article


Kadik, A., Koltashev, V., Kryukova, E., Plotnichenko, V., Tsekhonya, T., Kononkova, N. (2015) Solubility of nitrogen, carbon, and hydrogen in FeO–Na2O–Al2O3–SiO2 melt and liquid iron alloy: Influence of oxygen fugacity. Geochemistry International 53, 849–868.
Show in context

Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Li, Y., Marty, B., Shcheka, S., Zimmermann, L., Keppler, H. (2016) Nitrogen isotope fractionation during terrestrial core-mantle separation. Geochemical Perspectives Letters 2, 138–147.
Show in context

Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Libourel, G., Marty, B., Humbert, F. (2003) Nitrogen solubility in basaltic melt. Part I. Effect of oxygen fugacity. Geochimica et Cosmochimica Acta 67, 4123–4135.
Show in context

Mass balance is used to solve for the equilibrium distribution of N and S as required by partitioning and magma solubility constraints (Eq. 1; Libourel et al., 2003).
View in article


Ma, Z. (2001) Thermodynamic description for concentrated metallic solutions using interaction parameters. Metallurgical and Materials Transactions B 32, 87–103.
Show in context

Our value calculations follow the approach of Ma (2001).
View in article
Note that , , and refer to the expanded concentration expression associated with the notation of Ma (2001), and that positive coefficients indicate a reduction in the value.
View in article


Marty, B. (2012) The origins and concentrations of water, carbon, nitrogen and noble gases on earth. Earth and Planetary Science Letters 313–314, 56–66.
Show in context

Nitrogen is the most depleted element in the observable Earth when compared to CI chondrites, materials that most faithfully record the bulk composition of the solar system (observable Earth is the sum of the atmosphere, oceans, crust, and mantle) (e.g., Marty, 2012; Halliday, 2013).
View in article


Righter, K., Drake, M.J. (1996) Core formation in Earth’s moon, Mars, and Vesta. Icarus 124, 513–529.
Show in context

The role of pressure and temperature on core formation chemistry become larger for larger planetary bodies, as indicated by the refractory siderophile element concentrations observed for Earth, Mars, and Vesta (Righter and Drake, 1996).
View in article
It is well established within the Solar System that planetary body size correlates with average core formation pressures (Righter and Drake, 1996), and our work therefore predicts a direct relationship between planet size and its distribution of volatiles.
View in article


Roskosz, M., Bouhifd, M.A., Jephcoat, A.P., Marty, B., Mysen, B.O. (2013) Nitrogen solubility in molten metal and silicate at high pressure and temperature. Geochimica et Cosmochimica Acta 121, 15–28.
Show in context

Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article
Previously published high pressure experiments suggest that pressure favours N incorporation into cores (Roskosz et al., 2013, Grewal et al., 2019a); however, the pressure effect cannot be confidently resolved to be different from zero (Grewal et al., 2019a).
View in article
Our determination of a positive effect of pressure is consistent with previous work, but the larger range of P we investigated results in more precise determination of the P-T coefficients in Equations 1 and 2 (c.f., Roskosz et al., 2013, Grewal et al., 2019a).
View in article


Schlichting, H.E., Sari, R., Yalinewich, A. (2015) Atmospheric mass loss during planet formation: The importance of planetesimal impacts. Icarus 247, 81–94.
Show in context

With accretion and differentiation, terrestrial bodies form atmospheres, and these atmospheres can be lost (Tucker and Mukhopadhyay, 2014; Schlichting et al., 2015).
View in article


Speelmanns, I.M., Schmidt, M.W., Liebske, C. (2019) The almost lithophile character of nitrogen during core formation. Earth and Planetary Science Letters 510, 186–197.
Show in context

Studies agree that increasing temperature makes nitrogen less siderophile (Grewal et al., 2019a; Speelmanns et al., 2019).
View in article
Extrapolating the temperature effect on metal-silicate partition coefficients for nitrogen (, atomic) to the average core formation temperature for Earth implies lithophile behaviour (e.g., Speelmanns et al., 2019).
View in article
Our fO2 and temperature findings accord with previous work (Dalou et al., 2017; Grewal et al., 2019a; Speelmanns et al., 2019).
View in article
Core formation occurs in parallel to planetary accretion, and because nitrogen can display both siderophile and lithophile behaviour, the effect of core formation on the apparent nitrogen depletion in Earth remains uncertain (Roskosz et al., 2013; Kadik et al., 2015; Li et al., 2016; Dalou et al., 2017; Grewal et al., 2019a,b; Speelmanns et al., 2019).
View in article


Siebert, J., Corgne, A., Ryerson, F.J. (2011) Systematics of metal–silicate partitioning for many siderophile elements applied to earth’s core formation. Geochimica et Cosmochimica Acta 75, 1451–1489.
Show in context

Our model takes a single stage approach, and while accretion is a multi-step process, single stage calculations capture average P-T-X conditions that can be readily compared between elements (Siebert et al., 2011).
View in article
It is significant that higher pressure core formation (>30 GPa) and a core-mantle fO2 of ΔIW−2 can produce S/NOE ratios equal to Earth because these are the same differentiation conditions implied by moderately siderophile element and FeO concentrations in the observable Earth in a single stage framework (e.g., Siebert et al., 2011; Fischer et al., 2015).
View in article


Sossi, P.A., Burnham, A.D., Badro, J., Lanzirotti, A., Newville, M., O’Neill, H.S.C. (2020) Redox state of earth’s magma ocean and its venus-like early atmosphere. Science Advances 6, eabd1387.
Show in context

Atmospheric fO2 conditions are a free parameter in light of evidence for the depth dependence of the Fe+3/Fetot ratio of magma in equilibrium with iron (Zhang et al., 2017; Armstrong et al., 2019; Deng etal., 2020) and the temperature dependence of Fe+3/Fetot for any fO2 (Sossi et al., 2020).
View in article


Steelmaking Data Sourcebook (1988) The Japan Society for the Promotion of Science: The 19th Committee on Steelmaking. Gordon and Breach Science Publishers, New York.
Show in context

Our lower pressure experiments indicate that N-C interactions in Fe alloy make both elements less siderophile and yield an epsilon value () of 7.6 ± 1.4 (R2 = 0.90, p value = 0.004, n = 6; Fig. S-4), qualitatively consistent with the Steelmaking Data Sourcebook (1988) and values collected under C-free conditions (Grewal et al., 2021).
View in article


Tucker, J.M., Mukhopadhyay, S. (2014) Evidence for multiple magma ocean outgassing and atmospheric loss episodes from mantle noble gases. Earth and Planetary Science Letters 393, 254–265.
Show in context

With accretion and differentiation, terrestrial bodies form atmospheres, and these atmospheres can be lost (Tucker and Mukhopadhyay, 2014; Schlichting et al., 2015).
View in article


Wasson, J.T., Kallemeyn, G.W. (1988) Compositions of chondrites. Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences 325, 535–544.
Show in context

Data are from Hirschmann (2016) and Wasson and Kallemeyn (1988).
View in article


Zhang, H., Hirschmann, M., Cottrell, E., Withers, A. (2017) Effect of pressure on Fe3+/Fe ratio in a mafic magma and consequences for magma ocean redox gradients. Geochimica et Cosmochimica Acta 204, 83–103.
Show in context

Atmospheric fO2 conditions are a free parameter in light of evidence for the depth dependence of the Fe+3/Fetot ratio of magma in equilibrium with iron (Zhang et al., 2017; Armstrong et al., 2019; Deng etal., 2020) and the temperature dependence of Fe+3/Fetot for any fO2 (Sossi et al., 2020).
View in article



top

Supplementary Information

Abstract | Introduction | Methods | Results | Modelling N Distribution throughout Earth | Acknowledgements | References | Supplementary Information


The Supplementary Information includes:
  • Experimental Methods
  • Demonstrations of Equilibrium
  • Validation of Predictive Power of Equations 1 and 2
  • Further Details for Nitrogen Distribution Model
  • Model Results for Other Chondrite Groups
  • Tables S-1 to S-3
  • Figures S-1 to S-10
  • Supplementary Information References


Download Tables S-1 to S-3 (Excel).

Download the Supplementary Information (PDF).
top

Figures



Figure 1 Concentrations of volatile elements (H, C, S, and N) in Earth’s observable reservoirs normalised to CI chondrites. Nitrogen appears depleted in Earth compared to other volatiles. Data are from Hirschmann (2016)

Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.

and Wasson and Kallemeyn (1988)

Wasson, J.T., Kallemeyn, G.W. (1988) Compositions of chondrites. Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences 325, 535–544.

.
Back to article


Figure 2 Controls on nitrogen partitioning. (a) Oxidising conditions make nitrogen more siderophile. Literature data are recalculated to 1973 K and 1 GPa. (b) Nickel and (c) Sulfur make nitrogen less siderophile. (d) Hotter temperatures make nitrogen less siderophile. Extrapolation of temperature effect at 1 GPa underpredicts values compared to values measured for higher pressure LH-DAC experiments. (e) A 1∶1 comparison of observations versus predictions from Equation 1. The R2 for the literature fit is 0.76. (f) Application of Equation 1 to mantle liquidus geotherm at ΔIW−2. Nitrogen becomes more siderophile with depth.
Back to article


Figure 3 Models for redistribution of N and S in response to planetary differentiation. (a) An example of a model iteration of S/NOE as function of core formation pressure and atmospheric fO2. Solid lines are model solutions for a given atmospheric fO2 (log unit spacing). Horizontal dotted lines delineate the estimated range of S/NOE for Earth (Hirschmann, 2016

Hirschmann, M.M. (2016) Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. American Mineralogist 101, 540–553.

). Observable Earth is the sum of the atmosphere and magma ocean. (b) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (c) An example of a model iteration of S/NOE.MO as function of core formation pressure and atmospheric fO2. Horizontal dotted lines delineate the estimated range of S/NOE for Earth. Observable Earth is only the magma ocean. (d) Fractional success of the model between iterations for producing a S/NOE ratio that matches Earth. (e) Predicted bulk planet S as a function of core formation pressure. Dashed lines are 1 sigma that propagate uncertainties from the observable Earth S and parameterised values. The dotted line is the estimate of observable Earth S from the volatility trend (5600 ppm; Dreibus and Palme, 1996

Dreibus, G., Palme, H. (1996) Cosmochemical constraints on the sulfur content in the Earth’s core. Geochimica et Cosmochimica Acta 60, 1125–1130.

).
Back to article

  • Contact us
  • |
  • Subscribe
Sign up for our newsletter
  • Connect with us
  • facebook
  • twitter
  • youtube

Geochemical Perspectives Letters is a registered trademark of the European Association of Geochemistry
ISSN 2410-339X (print) | ISSN 2410-3403 (online)
EAG Privacy Policy