Skip to main content Accessibility help
×
Hostname: page-component-586b7cd67f-rdxmf Total loading time: 0 Render date: 2024-11-28T06:06:28.496Z Has data issue: false hasContentIssue false

2 - Origin and Early Differentiation of Carbon and Associated Life-Essential Volatile Elements on Earth

Published online by Cambridge University Press:  03 October 2019

Beth N. Orcutt
Affiliation:
Bigelow Laboratory for Ocean Sciences, Maine
Isabelle Daniel
Affiliation:
Université Claude-Bernard Lyon I
Rajdeep Dasgupta
Affiliation:
Rice University, Houston

Summary

This chapter reviews what is known about the fate of carbon during early differentiation of inner solar system planets. It reviews the nature of carbon fractionation in a magma ocean as compared to the core, mantle, and atmosphere, and how this may have varied between planetary bodies in the solar system. It discusses whether magma ocean processes could have established the present-day budget of carbon in Earth’s bulk silicate, and also reviews possibilities for the early temporal evolution of the mantle carbon budget through core formation, later veneer addition, and magma ocean crystallization processes.

Type
Chapter
Information
Deep Carbon
Past to Present
, pp. 4 - 39
Publisher: Cambridge University Press
Print publication year: 2019
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - SA
This content is Open Access and distributed under the terms of the Creative Commons Attribution licence CC-BY-NC-SA 4.0 https://creativecommons.org/cclicenses/

2.1 Introduction

Earth’s unique status as the only life-harboring planet in the inner Solar System is often linked to its orbit being in the habitable zone that resulted in the stabilization of liquid water on its surface. However, in addition to surface liquid water, the presence of carbon and other life-essential volatile elements (LEVEs) such as nitrogen and sulfur in the surface environment is also critical for the habitability of rocky planets. Earth’s long-term equable climate and chemically habitable surface environment are therefore results of well-tuned fluxes of carbon and other LEVEs, involving the deep and the surficial Earth. Given that Earth developed into a volcano-tectonically active planet with both outgassing and ingassing mechanisms, the surface inventory of LEVEs over million- to billion-year timescales is maintained by interactions of the ocean–atmosphere system with the silicate fraction of the planet. Yet it remains uncertain how and when the initial inventory of LEVEs for the surface reservoir plus the silicate fraction (bulk silicate Earth (BSE) altogether), which set the initial boundary conditions for subsequent planet-scale volatile cycles, got established. The answers to these questions lie within Earth’s formative years – in the building blocks, in the process of Earth’s accretion, and in the early differentiation that made the major reservoirs that constitute the core, the mantle, and the crust–atmosphere.

Earth is a differentiated planet with a central metallic core and an outer shell of dominantly silicate rocks divided into the mantle and crust, which is overlain by the fluid envelope of the ocean–atmosphere system. Direct constraints on carbon and the other life-essential volatile abundance of Earth’s core are lacking (Chapter 3, this volume), although significant concentrations of carbon, sulfur, hydrogen, and nitrogen can be in Earth’s core given that the outer core is 10% lighter than pure Fe–Ni alloy liquid.Reference Birch1Reference Poirier3 More is known about the carbon content of the present-day silicate fraction of the planet and its fluid envelope because Earth’s mantle reservoir is sampled to some extent via mantle-derived melts (e.g. Refs. Reference Javoy, Pineau and Allègre4Reference Dasgupta and Hirschmann9 and Chapter 9, this volume). Yet the estimates of bulk mantle carbon abundance from concentrations in basalts are uncertain partly because most basalts are partially degassed and basalt generation only directly samples typically the top 100–200 km of the mantle. Nonetheless, largely based on the carbon content of basalts (e.g. C/3He, CO2/Nb, CO2/Ba, and C/40Ar ratios), the BSE abundance of carbon is estimated to be ~100–530 ppm.Reference Hirschmann and Dasgupta7, Reference Dasgupta and Hirschmann9Reference Marty13 Among this range of estimated C budget for the BSE, most estimatesReference Hirschmann and Dasgupta7, Reference Dasgupta and Hirschmann9, Reference Bergin, Blake, Ciesla, Hirschmann and Li10Reference Halliday12 have converged to a concentration of 90–130 ppm, and this will be used in this chapter. It can be surmised that this BSE budget of C is an end product of Earth’s accretion and early differentiation,Reference Dasgupta8, Reference Dasgupta and Walker14, Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 but the uncertainty in the C budget of BSE makes testing such hypotheses quantitatively less rigorous. Further insight could be gained, however, from C to other LEVE ratios (i.e. C/X, where X = N, S, or H; e.g. Refs. Reference Hirschmann and Dasgupta7, Reference Bergin, Blake, Ciesla, Hirschmann and Li10, Reference Hirschmann11, Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18). Specifically, these volatile ratios for the terrestrial reservoirs in general and BSE in particular could be compared with the same ratios in various planetary building blocks and evaluated against the known processes of planet formation and early differentiation.

The key questions focus on what planetary building blocks were responsible for delivering C and other associated LEVEs to Earth and when such delivery processes took place. If the C abundance and the C/X ratios of the BSE as well as the isotopic compositions of all LEVEs are similar to any of the undifferentiated meteorites, it may be concluded that such building blocks were responsible for bringing these volatiles to Earth and no postdelivery process altered their ratios (Figure 2.1). On the other hand, if the C abundance and C/X ratios of the BSE are distinct from these known building blocks, one needs to investigate whether any postdelivery processes could have altered these quantities such that the BSE LEVE geochemistry could be established. For example, if core formation processes fractionate C and/or affect the C/X ratios, then could equilibrium core formation scenarios be responsible for establishing carbon, oxygen, hydrogen, nitrogen, and sulfur in BSE (Figure 2.1)? The timing of volatile delivery is also important. If the LEVEs were delivered during the main stage of accretion, then early differentiation processes like core–mantle segregation and magma ocean (MO) degassing would have played an important role in setting up the relative budgets of the LEVEs in the primary reservoirs (Figure 2.1).Reference Dasgupta8, Reference Hirschmann11 The delivery mechanism and the timing of delivery may be intimately linked; however, less interrogation has taken place thus far regarding the delivery mechanism. For example, there is a growing recognition that following an initial period of accretion of undifferentiated planetesimals, differentiated planetary embryos with masses ranging from ~0.01 to 0.1ME (where ME is the mass of present-day Earth) were abundant in the inner Solar System;Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19 therefore, accretion of these differentiated bodies significantly contributed to the growth of larger planets like Earth.Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19, Reference O’Brien, Walsh, Morbidelli, Raymond and Mandell20 Did the impacts or merger of differentiated bodies deliver the LEVEs to the BSE and establish the C abundance and the C/X ratiosReference Li, Dasgupta, Tsuno, Monteleone and Shimizu16, Reference Tsuno, Grewal and Dasgupta17 or did post-core formation addition of undifferentiated or differentiated bodiesReference Marchi, Canup and Walker21 establish the BSE LEVE budget? And what were the relative contributions of the impacts of large differentiated and smaller undifferentiated bodies to bringing in LEVEs versus causing a depletion of LEVEs via atmospheric mass loss?Reference Schlichting and Mukhopadhyay22Reference Ahrens24

Figure 2.1 Plausible initial distributions of C, H, N, and S between Earth’s major reservoirs at the end of accretion and core formation. (a) LEVEs are available during accretion and all LEVEs are initially sequestered in the core and the MO. (b) LEVE delivery during accretion is temporally separate from the main episodes of core formation and thus the core is effectively devoid of LEVEs. (c) Available LEVEs during accretion are effectively sequestered in the core or lost to space (not shown) and the BSE acquires LEVEs via later additions. (d) Core–mantle differentiation takes place in the absence of LEVEs and LEVEs are added via late additions. These distributions eventually evolved into the setup seen in the modern world (e), where these volatile elements are thought to exchange between Earth’s mantle and the crust–ocean–atmosphere system (ingassing and outgassing) regulates the long-term habitability of the planet.

This chapter will discuss the state of the art on what is known about the origin of C relative to the other LEVEs in the BSE, highlighting open questions and gaps in knowledge. We will discuss constraints based on geochemical comparison of the potential terrestrial building blocks with Earth reservoirs and interrogate each key process of planet formation and early differentiation, such as core formation, MO crystallization, MO degassing, and atmospheric loss in shaping the absolute and relative abundances of C and the other LEVEs. The current understanding of the fractionation of C and other LEVEs during the core formation process will also be used to put constraints on the C budget of the core.

2.2 Constraints on the Compositions of Terrestrial Building Blocks

2.2.1 Constraints from Isotopes of Refractory Elements

Geochemical models based on elemental abundances and isotopic constraints predict the accretion of bulk Earth primarily from chondritic meteorites.Reference Ringwood25, Reference Anders and Ebihara26 Because elemental and isotopic fractionation trends for bulk silicate samples and chondritic meteorites have well-defined relationships, it is widely assumed that the building blocks of Earth belong to a class of primitive objects that are sampled by chondritic meteorites.Reference Allègre, Manhès and Lewin27 Close similarity of the abundances of refractory lithophile elements between the BSE and CI led to an early development of a CI model for Earth.Reference Ringwood25 Allègre et al.Reference Allègre, Manhès and Lewin27 further consolidated this model by arguing that volatile depletion trends in the BSE are matched most closely by CI chondrites. However, a distinct disparity of oxygen (second most abundant element in the BSE by mass) isotopes (i.e. 16O, 17O, and 18O) between CI chondrites and the BSE, which could not be explained by any fractionation or differentiation models,Reference Clayton, Onuma and Mayeda28, Reference Clayton29 was one of the first observed isotopic variations that challenged the CI model for Earth. Over the past decades, isotopic dissimilarity between CI chondrites and the BSE has been extended from Δ17O (deviation from a reference mass-dependent fractionation line in a triple isotope of oxygen plotReference Young30) to several radiogenic as well as non-radiogenic isotopic anomalies for a wide suite of refractory elements with contrasting geochemical behaviors: lithophile elements (ε48Ca (Ref. Reference Dauphas31), ε50Ti, μ142Nd), moderately siderophile elements (ε54Cr, ε64Ni, ε92Mo), and highly siderophile elements (HSEs; ε100Ru) (Refs. Reference Dauphas32, Reference Fischer-Gödde and Kleine33 and references therein). Interestingly, the BSE isotopic compositions for all of these elements are almost indistinguishable from enstatite chondrites (E-chondrites). This has led to the development of an E-chondrite model for Earth.Reference Dauphas32, Reference Javoy34, Reference Javoy35 Tracking the isotopic evolution of the growing BSE by taking into account the alloy–silicate partitioning behavior of several lithophile, moderately siderophile, and HSEs during core–mantle equilibration, DauphasReference Dauphas32 proposed that half of the first 60% of accreted mass was composed of E-type chondrites, while the remaining 40% was all E-type chondritic materials. However, there are a few outstanding issues with an E-chondritic model for Earth, namely: (1) the Ca/Mg ratio of the BSE (0.118Reference McDonough and Sun36) is distinctly higher than that of E-chondrites (0.080Reference Wasson and Kallemeyn37), suggesting that the BSE has an enriched refractory elemental abundance relative to E-chondrites;Reference Fitoussi, Bourdon and Wang38 (2) E-chondrites have an Mg/Si ratio that is lower than that of the BSE sampled by upper-mantle rocks,Reference Palme, O’Neill and Carlson39 which can only be explained by either incorporation of Si into the core or via an unsampled lower mantle with a low Mg/Si ratio;Reference Murakami, Ohishi, Hirao and Hirose40 and (3) δ30SiBSE is 0.3–0.4‰ higher than δ30Si of E-chondrites and aubrites.Reference Fitoussi and Bourdon41 Experimentally determined Si isotopic fractionation parametersReference Shahar42, Reference Shahar43 rule out the possibility of explaining the Si isotopic composition as well as the Mg/Si ratio of the BSE via incorporation of Si in the core if Earth was primarily accreted from E-chondrites.Reference Fitoussi and Bourdon41 Although there are studies that argue for a lower mantle or at least a portion of the lower mantle being richer in silica compared to Earth’s upper mantle,Reference Murakami, Ohishi, Hirao and Hirose40, Reference Ballmer, Houser, Hernlund, Wentzcovitch and Hirose44 such findings are debated. To circumvent these issues, Dauphas et al.Reference Dauphas, Poitrasson, Burkhardt, Kobayashi and Kurosawa45 proposed that the bulk Earth and the E-chondrites sample an isotopically similar reservoir with their subsequent chemical evolution offset by ensuing nebular fractionation and planetary differentiation processes. In summary, a large number of geochemical systematics point to the bulk Earth being made from E-chondrite-type materials, although questions on the major element and volatile element abundance mismatch remain.

2.2.2 Constraints from Isotopes of Highly Volatile Elements

In contrast to constraints from refractory elements, constraints from highly volatile elements (i.e. C, N and H (water)) paint a different picture on the potential sources of LEVEs. Although ordinary chondrites (OCs), which sample S-type asteroids in the inner belt (2.1–2.8 AU), have similar amounts of bulk water as the BSE, the D/H ratio of the BSE is not similar to that of OCs,Reference McNaughton, Borthwick, Fallick and Pillinger46 but rather is strikingly similar to that of carbonaceous chondrites, especially CI chondrites,Reference Robert, Merlivat and Javoy47 which are representative of C-type asteroids (beyond 2.8 AU) (Figure 2.2a). Similarity of 15N/14N between the BSE and CI chondrites also points toward a carbonaceous chondritic origin of volatiles on Earth (Figure 2.3b).Reference Marty13, Reference Alexander48, Reference Pearson, Sephton, Franchi, Gibson and Gilmour49 13C/12C of the BSE cannot distinguish between a chondritic or cometary origin of C (Figure 2.2b).Reference Marty, Alexander and Raymond50, Reference Sarafian, Nielsen, Marschall, McCubbin and Monteleone51 However, a clear distinction of the 15N/14N ratio of the BSE from the corresponding cometary values would necessitate C delivery by CI chondrites as well, assuming all the major volatiles (i.e. C, N, and water) were sourced from a similar parent body.

Figure 2.2 Comparison of isotopic compositions of (a) hydrogen and nitrogen and (b) carbon and nitrogen for several Solar System objects and reservoirs. The solar reservoir is depleted while the cometary reservoir is enriched in the heavier isotopes of nitrogen and hydrogen relative to all classes of meteorites as well as the BSE. Although δ13C alone cannot distinguish between the meteoritic and cometary sources, the isotopic compositions of all major LEVEs suggest that the BSE had a similar parent reservoir as carbonaceous chondrites.

Although thus far direct matching of the isotopes of C, H, and N of the BSE with undifferentiated meteorites has been the approachReference Marty and Yokochi59 to determining the source of the LEVEs on Earth, there is a growing realization that there may be fractionation of stable isotopes during early terrestrial differentiation, such as devolatilization, MO degassing, and core–mantle fractionation. Indeed, such fractionation processes may be able to reconcile some of the observed differences between the LEVE isotopic compositions of the BSE and the chondrites.Reference Wood, Li and Shahar60 For example, carbon isotopic compositions of many CI chondritic materials are distinctly lighter (δ13C, approximately –15 to –7‰Reference Marty, Alexander and Raymond50) than the average carbon isotope composition of Earth’s mantle (δ13C, approximately –5‰Reference Kerridge61, Reference Deines62). Similarly, the sulfur isotope composition of Earth’s mantle has also been argued to be non-chondritic.Reference Labidi, Cartigny and Moreira63 Graphite–Fe–carbide melt carbon isotope fractionation experimentsReference Satish-Kumar, So, Yoshino, Kato and Hiroi64 showed that 12C preferentially incorporates into the metallic phase, leaving a 13C-enriched signature in the graphite or diamond. If similar fractionation behavior applies between silicate MO and core-forming alloy liquid, then core formation could explain how Earth’s mantle evolved to have a heavier carbon isotope value compared to the chondritic building blocks.

2.2.3 Constraints from Theoretical Modeling

Because isotopic compositions of refractory elements suggest bulk Earth to be composed of materials similar to E-chondrites while volatile elements favor their delivery via CI chondritic material, the timing of admixing of C-rich material into the accreting zone of Earth becomes key. Such admixing could have taken place either during the primary stage of Earth’s growth or after the main phase of accretion. Similarly, delivery of C-bearing materials may have occurred via accretion of undifferentiated bodies similar to chondrites or via amalgamation of differentiated bodies with C abundances in relevant reservoirs of those bodies set by core–mantle differentiation and atmospheric losses.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 Theoretical models have been used to simulate the early evolution of the Solar System by accounting for dynamics of planetary accretion along with geochemical, cosmochemical, and chronological relationships between accreting and resulting bodies.Reference Walsh, Morbidelli, Raymond, O’Brien and Mandell65, Reference Raymond, O’Brien, Morbidelli and Kaib66 Constrains from Hf/W chronometry predict that Mars-sized planetary embryos were formed from mostly E-type chondritic material sourced from the inner part of the Solar System in the initial ~5 Ma.Reference Dauphas and Pourmand67 Recent advancements like pebble accretion models predict that gas giants like Jupiter also grew synchronously to their present-day mass via trapping of proto-solar nebular gas within the lifetime of the protoplanetary disk.Reference Lambrechts and Johansen68, Reference Johansen, Low, Lacerda and Bizzarro69 However, the growth of gas giants is thought to rapidly deplete the asteroid belt, which would mean minimal interaction of C-rich carbonaceous chondritic material that condensed in the outer parts of the Solar System with the growth zone of terrestrial planet accretion within 2 AU.Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19 With the consideration that the presence of giant planets in the system affects the mixing and delivery of LEVEs from the outer regions to the terrestrial planet-forming region in the disk, the migration of giant planets has been argued to be necessary to perturb the orbits of cometary and asteroid-like bodies, bringing a fraction of these to the inner regions, where they can collide with terrestrial planets, delivering additional volatile elements to the rocky planets.Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19 A specific planetary migration model, the Grand Tack scenario,Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19, Reference O’Brien, Walsh, Morbidelli, Raymond and Mandell20, Reference Walsh, Morbidelli, Raymond, O’Brien and Mandell65 postulated that the inward and outward movement of Jupiter likely repopulated the asteroid belt in the inner parts of the Solar System with C-type chondritic material. However, these models cannot constrain the net mass of volatile-rich material being delivered, as they allow for a total influx of volatile-rich material ≥0.5–2.0% of present-day Earth’s mass, assuming volatile-rich CI chondrites have water in the range of 5–10 wt.%.Reference Morbidelli, Lunine, O’Brien, Raymond and Walsh19, Reference O’Brien, Walsh, Morbidelli, Raymond and Mandell20 Delivery of a C budget that is larger than the present day in the BSE would require offsetting via later differentiation processes, while C delivery being restricted to the present-day budget of the BSE would have to be unaffected by later differentiation processes. A model such as Grand Tack does not directly constrain the relative timing of acquiring the C- and water-rich chondritic materials in the terrestrial planet-forming regions with respect to specific giant impact events. However, with the Moon formation being a specific event, similarity or a lack thereof between C and other LEVE budgets and volatile isotope signatures of Earth and Moon may shed light on when LEVE-rich materials were brought in with respect to the Moon-forming event. Little is known about the C budget of the Moon or the lunar mantle;Reference Wetzel, Hauri, Saal and Rutherford70, Reference Chi, Dasgupta, Duncan and Shimizu71 however, if water and carbon were co-delivered, the appreciable water content of lunar glasses and melt inclusionsReference Saal72, Reference Hauri, Weinreich, Saal, Rutherford and Van Orman73 and the E-chondritic nature of late veneerReference Dauphas32 can only be reconciled if volatiles were delivered to the Earth–Moon system beforeReference Saal, Hauri, Van Orman and Rutherford74, Reference Greenwood75 or during the Moon-forming impact.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 However, if C and other LEVEs were delivered during the main stage of Earth’s accretion, then differentiation processes like core–mantle separation, as well as the delivery mechanism of volatiles – either by smaller, undifferentiated planetesimals composed of primarily CI-chondritic material or via relatively large, differentiated planetary embryos heterogeneously composed of material sourced from both inner and outer parts of the Solar System – can provide additional information on the origin of LEVEs in the BSE.

2.3 C and Other Volatiles: Abundances, Ratios, and Forms in Various Classes of Meteorites and Comparison with the BSE

The present-day BSE is estimated to be depleted in C and other highly volatile elements relative to all classes of undifferentiated meteorites by at least an order of magnitude and by as much as two orders of magnitude compared to the more primitive carbonaceous chondrites (Table 2.1).Reference Marty13 For example, carbonaceous chondrites are C rich,Reference Alexander48, Reference Alexander, McKeegan and Altwegg54, Reference Palme, Lodders, Jones, Holland and Turekian76 while enstatite and OCs are relatively C poor.Reference Alexander, McKeegan and Altwegg54, Reference Kerridge61 Among various carbonaceous chondrites, CI chondrites are the most C rich,Reference Alexander, McKeegan and Altwegg54, Reference Palme, Lodders, Jones, Holland and Turekian76Reference Palme, Jones and Davis79 and they are also richer in other LEVEs such as N (0.19–0.32 wt.%Reference Alexander, McKeegan and Altwegg54, Reference Palme, Lodders, Jones, Holland and Turekian76, Reference Anders and Grevesse77) and H (1.55–2.02 wt.%Reference Alexander, McKeegan and Altwegg54, Reference Anders and Grevesse77) compared to other types of carbonaceous chondrites such as CV, CO, CM, and CR (Table 2.1). Severe depletion of C in the BSE relative to CI chondrites can be explained either by accretion of extremely small quantities of CI chondrites as late additions or via accretion of larger quantities of CI chondrites, but subject to subsequent loss during accretion or differentiation. Such loss of C from carbonaceous building blocks may be expected given that primary carbon is present as soluble as well as insoluble organic molecules in carbonaceous chondrites, which are expected to be unstable at the high temperatures occurring during inner-Solar System processes. Similarly, the carbonates that are products of secondary alterations in carbonaceous chondrites are also expected to volatilize during delivery and accretion in the inner-Solar systems bodies. Even though E-chondrites may have C abundances that are not much greater than that in the BSE, such concentrations are more likely to survive the delivery to the inner-Solar System bodies given their presence in the form of more refractory graphite, diamonds, and carbides.Reference Grady and Wright80 It remains unclear, however, in what form and abundance C was in E-chondrite parent bodies before thermal metamorphism. It cannot be ruled out that unprocessed E-chondritic protoliths were more C rich and may have even contained organics. Despite the survival potential of refractory C-bearing phases in E-chondrites in the inner-Solar System processes, it is generally thought that the E-chondrites had much less water than the present-day BSE,Reference Hutson and Ruzicka81 although the H budget of the BSE has significant uncertainty owing to poor constraints on the water budget of the mantle.Reference Hirschmann82

Table 2.1 C content and C/S, C/N, and C/H weight ratios of major terrestrial reservoirs and chondritic building blocks

ReservoirsC (wt.%)C/SC/NC/H
Terrestrial reservoirs
BSE0.011 ± 0.0020.49 ± 0.1440.00 ± 8.001.13 ± 0.20
Mantle0.008 ± 0.0020.36 ± 0.1072.73 ± 37.722.00 ± 0.70
Crust–ocean–atmosphere0.002 ± 0.0008.89 ± 2.2212.5 ± 3.220.51 ± 0.04
Core0.50 (?)0.32 (?)85.00 (?)8.40 (?)
Carbonaceous chondrites
CO0.63 ± 0.240.32 ± 0.1221.74 ± 19.1610.40 ± 0.00
CV0.92 ± 0.430.42 ± 0.2023.96 ± 24.1610.20 ± 0.00
CM1.89 ± 0.480.58 ± 0.1521.01 ± 8.353.20 ± 0.00
CI4.24 ± 0.770.72 ± 0.1319.66 ± 11.694.80 ± 0.00
Enstatite chondrites
EH0.36 ± 0.130.06 ± 0.0213.73 ± 7.2611.00 ± 0.00
EL0.48 ± 0.160.15 ± 0.0524.43 ± 9.4011.00 ± 0.00
Ordinary chondrites
H0.110.06>104.002.50 ± 0.70
L0.090.0446.10 ± 1.003.00 ± 1.00
LL0.090.0551.22 ± 1.202.60 ± 0.50

The mantle, the crust-ocean-atmosphere, and the BSE averages and 1σ standard deviations are from Hirschmann and DasguptaReference Hirschmann and Dasgupta7 for C/H, Bergin et al.Reference Bergin, Blake, Ciesla, Hirschmann and Li10 for C/N, and HirschmannReference Hirschmann11 for C/S. The carbonaceous chondrite data for C, H, and N are from Figure 2.3, E-chondrite C/N data are from Grady et al.Reference Grady, Wright, Carr and Pillinger52 and Alexander et al.Reference Alexander, Swan and Prombo53 OC C/N data are from Bergin et al.Reference Bergin, Blake, Ciesla, Hirschmann and Li10 and references therein. The C/S data for all chondrites are from Wasson and Kallemeyn,Reference Wasson and Kallemeyn37 and the C/H data for enstatite and OCs are from HirschmannReference Hirschmann11 and references therein. The core compositional estimates are using bulk Earth concentrations from McDonough,Reference McDonough and Carlson83 and Dxalloy/silicate values are from Figure 2.5 and assuming complete core–mantle equilibration. The data for C contents are from Bergin et al.Reference Bergin, Blake, Ciesla, Hirschmann and Li10 for terrestrial reservoirs, from Alexander et al.Reference Alexander48 for carbonaceous chondrites, from Grady et al.Reference Grady, Wright, Carr and Pillinger52 for E-chondrites, and from Wasson and KallemeynReference Wasson and Kallemeyn37 for OCs.

In contrast to absolute abundances, the relative abundance of C with respect to other volatile elements in the BSE (i.e. C/H, C/N and C/S) is a powerful tracer to constrain the C abundance of Earth and provides additional constraints on the processes that must have fractionated the geochemical reservoirs relative to the cosmochemical reservoirs (Figure 2.3).

Figure 2.3 Comparison of the abundances of C, N, S, and H in the BSE normalized to their respective abundances in CI chondrites. In the BSE, C is enriched relative to N, has a similar abundance to S, and is depleted relative to H in comparison to its abundance in CI chondrites. Data from Refs. Reference Bergin, Blake, Ciesla, Hirschmann and Li10, Reference Hirschmann11, and Reference Marty13. Also noted in the figure is the main geochemical information used to estimate the BSE concentrations of the respective elements in the original studies. MORB = mid-ocean ridge basalt.

Due to the comparable 50% condensation temperatures of C and N as well as their abundances being positively correlated in various classes of chondrites, comparison of the C/N ratio of the BSE with that of the building blocks can help track the accretion and differentiation history of C in the early Earth. Recently, Bergin et al.Reference Bergin, Blake, Ciesla, Hirschmann and Li10 estimated the BSE C/N ratio to be 40 ± 8 by calculating the mantle C content using CO2/Nb and CO2/Ba ratios of possibly undegassed ocean island basalts (Figure 2.3).Reference Rosenthal, Hauri and Hirschmann84 This C/N ratio is similar to the estimate of Halliday,Reference Halliday12 but lower than the estimate of Marty.Reference Marty13 We will adopt the C/N ratio from the latest studies.Reference Hirschmann11, Reference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Dalou, Hirschmann, von der Handt, Mosenfelder and Armstrong85 The C/N ratios of volatile-rich CI and CM as well as slightly volatile-poor CO and CV chondrites are lower than the C/N ratio of the BSE by approximately a factor 2, while E-chondrites also have low C/N ratios (5–25; Figure 2.4a). Although OCs show large scatter in C/N ratios due to uncertainties in determining extremely low C and N contents, their average C/N ratios are higher than that of the BSE (Figure 2.4a). Because the average C/N ratio of all classes of chondrites lies in the range of 5–30, the superchondritic C/N ratio of the BSE has been postulated to be a legacy of either preferential incorporation of N into the core or loss of an early N-rich atmosphere.Reference Marty13 However, ureilites, a class of achondrites and assumed to be the mantle restite of an asteroidal body, are extremely C rich (average 3 wt.% C) with intergranular veins of graphite and traces of diamonds.Reference Nabiei86, Reference Kana, Masayuki and Jun-ichi87 It is interesting to note that the C/N ratio of ureilites is much higher than that of the BSE;Reference Hirschmann11, Reference Barrat, Sansjofre, Yamaguchi, Greenwood and Gillet88 therefore, if the C/N ratios of this achondrite are representative of an ureilite parent body and such a parent body significantly contributed to Earth’s accretion,Reference Schiller, Bizzarro and Fernandes89 then preferential loss of N as a necessary condition for bulk silicate reservoirs in Earth is precluded.

Figure 2.4 Comparison of bulk C/N (a), C/H (b), and C/S (c) weight ratios between different terrestrial and chondritic reservoirs. The data are from Table 2.1. CI = Ivuna-type Carbonaceous chondrite; CM = Mighei-type Carbonaceous chondrite; CV = Vigarano-type Carbonaceous chondrite; CO = Ornans-type Carbonaceous chondrite; LL = low total iron, low metal ordinary chondrite; H = highest total iron ordinary chondrite; L = low total iron ordinary chondrite; EL = (low enstatite) enstatite chondrite; EH = (high enstatite) enstatite chondrite.

Assuming a similar C inventory in the mantle as discussed earlier, the C/H ratio of the BSE has been estimated to be 1.13 ± 0.20.Reference Hirschmann and Dasgupta7, Reference Hirschmann11 The C/H ratio of BSE is lower than all classes of chondritic meteorites (Figure 2.4). OCs have a C/H ratio that is two to five times greater than that of the BSE, while E-chondrites, which are extremely water poor, have a 20‑times higher C/H ratio. Carbonaceous chondrites, which have the closest match to the terrestrial D/H ratio, also have a higher C/H ratio; CI and CM chondrites, which are C and water rich, have a C/H ratio that is four to six times higher; while CO and CV chondrites, which are relatively C and water poor, have a C/H ratio that is 12–15 times higher. Therefore, in contrast to a superchondritic C/N ratio, a subchondritic C/H ratio can be explained if C was preferentially incorporated into the core or lost to space relative to H during terrestrial differentiation.

Sulfur has a much higher solar nebula condensation temperature than C, N, and H; therefore, it is not a highly volatile element. Sulfur abundance in the BSE, especially in the mantle, is also better knownReference Ding and Dasgupta90Reference Rehkämper92 relative to that of other highly volatile elements of interest (Figure 2.3). Therefore, the C/S ratio can also act as an important tracer to track the evolution of C in the early Earth. The C/S ratio of the BSE is higher than that of C-depleted E-chondrites, but it is similar to or lower than that of C-rich carbonaceous chondrites (Figure 2.4c). Therefore, a near-chondritic C/S ratio in the BSE stipulates that either both C and S showed similar loss/gain behavior during early differentiation processes or they were primarily delivered in post-differentiation processes via undifferentiated, chondritic materials.

Because the C/H ratio of the BSE is subchondritic, the C/S ratio is near-chondritic, while the C/N ratio is superchondritic, H, N, and S have widely different geochemical behaviors relative to C during core formation or MO degassing and/or they track different stages of terrestrial accretion. The following sections will explore whether there are known processes that can explain the difference between the LEVE geochemistry of the BSE and those of carbonaceous chondrites.

2.4 Establishing LEVE Budgets of the BSE After Core Formation?

2.4.1 The Role of Late Accretion

The theory of the addition of undifferentiated meteorites after the completion of metal–silicate differentiation chiefly stems from the approximately chondritic relative abundance of HSE (e.g. Re, Os, Ir, Ru, Pt, Rh, Pd, and Au) concentrations in Earth’s mantle and the BSE, plus the fact that the BSE is not as depleted of HSEs as it would have been if the entire HSE inventory of the bulk Earth was available before equilibrium core formation took place.Reference Walker93, Reference Mann, Frost, Rubie, Becker and Audétat94 Assuming that core formation originally left the silicate Earth virtually HSE free and later additions of chondritic materials elevated the HSE abundances to the present level, the chondritic mass added is estimated to be 0.5% of Earth’s mass (0.005 ME). WalkerReference Walker93 noted that this estimate may be ~1.5 times higher or a factor of four lowerReference Morgan, Walker, Brandon and Horan95 depending on the exact concentrations of HSEs in the chondritic materials being added. Willbold et al.Reference Willbold, Elliott and Moorbath96 also calculated that in order to offset the pre-late accretion 182W-excess composition of the BSE, 0.08 ME chondritic mass needs to be added. Given that many carbonaceous chondrites are rich in LEVEs, many studies have surmised that late accretion of these materials brought all of the LEVEs to the BSE.Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Albarede97Reference Dasgupta, Chi, Shimizu, Buono and Walker99 If this was the mechanism, then the core formation would have no influence on the relative abundance of various LEVEs. Indeed, the addition of 0.0030–0.0075 ME chondritic materials with 3.48–3.65 wt.% C to the BSE would set the BSE C budget to ~100–250 ppm C (or ~360–900 ppm CO2 in the mantleReference Dasgupta8). Although this C budget may be sufficient to match even the enriched basalt source regions, as briefly discussed by Dasgupta,Reference Dasgupta8 there are several challenges in invoking CI-type carbonaceous chondrite as the chief or sole source of the BSE C and other LEVEs. Section 2.2 already discussed which elemental and isotopic compositions of the BSE are not satisfied if the BSE is chiefly made up of carbonaceous chondrites. Section 2.2 also outlined how the elemental abundance of C and the C/S ratio of the BSE can be satisfied if CI-type carbonaceous chondrite accreted to Earth after core formation was complete, but the C/N and C/H ratios and possibly δ13C and δ34S of the BSE are difficult to explain. Here, we specifically discuss what challenges exist in invoking other known meteorites as the late-accreting materials to explain the entire inventory of carbon and other LEVEs in the BSE.

Given E-chondrites (e.g. (low enstatite) enstatite chondrite (EL) and (high enstatite) enstatite chondrite (EH)) are depleted in water and also poor in C and N,Reference Grady, Wright, Carr and Pillinger52 bringing LEVEs in general and water in particular to Earth exclusively by E-chondrite is not thought to be a realistic mechanism. EH and EL chondrites also have distinctly lower C/N and C/S ratios compared to the BSE (Figure 2.4a and c). Hence, unprocessed E-chondrite, known from the studied samples, could not have brought LEVEs to Earth after core formation in the right proportions. This is at odds with recent propositions that, based on Ru isotopic composition, argued for the late veneers having E-chondritic character.Reference Dauphas32, Reference Fischer-Gödde and Kleine33 One issue in evaluating E-chondrites as plausible terrestrial building blocks based on volatile element geochemistry is that all known E-chondrites studied thus far are heavily metamorphosed; hence, before thermal metamorphism, these chondrites might have been more LEVE rich and C could have been originally organic.Reference Alexander, Fogel, Yabuta and Cody100

OCs are depleted in volatiles with respect to carbonaceous chondrites chiefly because the volatile-rich matrix volume percentage is less than that of carbonaceous chondrites.Reference Alexander, McKeegan and Altwegg54 In fact, bulk C content of most primitive OCs can be perfectly explained by its matrix material being entirely CI-type material.Reference Alexander48 OC matrices are drier, however, with bulk H of OCs falling below the concentration of what would be expected if the entire OC matrix was made up of CI-type materials.Reference Alexander, McKeegan and Altwegg54 Therefore, the C/H ratio of OCs is higher than that of the BSE (Figure 2.4b). Similarly, the C/S ratios of OCs are lower and the C/N ratios of OCs are distinctly higher compared to the estimates for the BSE (Figure 2.4a and c). Hence after core formation, OC delivery cannot be considered as the chief process of origin of Earth’s LEVEs.

Among all stony meteorites, ureilites are the only group of primitive achondrites that are rich in C and depleted in N, leading to higher C/N ratios.Reference Vdovykin101Reference Downes103 Carbon in ureilites is also present in relatively non-labile phases (i.e. as graphites, nanodiamonds, and lonsdaleite) and hence is likely to survive high-temperature accretional processes. Primarily motivated by their high C/N ratios, recent studiesReference Hirschmann11, Reference Barrat, Sansjofre, Yamaguchi, Greenwood and Gillet88 have suggested that C-rich achondrites such as ureilites may be responsible for establishing the abundances of LEVEs in the BSE through post-core formation addition. While the C/N ratio of the BSE could be heavily shaped by the addition of a ureilite-like achondrite, the C/S ratio of ureilite is significantly higher than that estimated for the BSE. Hence, any attempt to match the C/S ratios of the BSE through addition of a ureilite-like achondrite-rich late veneer requires a very unique consideration of the composition and extent of core–mantle equilibration of proto-Earth.Reference Hirschmann11 No model of ureilite-like late-veneer addition can explain the C/N and C/S ratios simultaneously. For example, scenarios of C-rich late-silicate addition in which the superchondritic C/N ratio can be established also result in a superchondritic C/S ratio.Reference Hirschmann11

Overall, late accretion alone does not seem to account for the C budget and C to other LEVE ratios of the present-day BSE (Figure 2.5). Hence, other early Earth processes need to be considered.

Figure 2.5 Application of alloy–silicate partition coefficients and solubilities of LEVEs in the silicate melts to examine the effect of core formation, with varying degrees of alloy–silicate equilibration, with or without loss of an early atmosphere formed via MO degassing on the remnant abundances of LEVEs in the bulk silicate reservoir. (a) LEVEs, when delivered as 0.015 ME late-accreting materials (i.e. 0% alloy–silicate equilibration), cause the volatile abundance to be higher than the present-day BSE. Core formation with increasing degrees of alloy–silicate equilibration increasingly depletes the remnant MO in all LEVEs, with C being much more depleted than other LEVEs, leading to subchondritic C/N, C/H, and C/S ratios. (b) Combining early atmospheric loss with core formation cannot offset C loss to the core due to the lower solubility of C relative to the other LEVEs in the silicate MOs. Bulk Earth volatile abundance data are from McDonough,Reference McDonough and Carlson83 while the alloy–silicate partition coefficients in a deep MO (P = 50 GPa, T = 3500 K; e.g. Siebert et al.Reference Siebert, Badro, Antonangeli and Ryerson104) for C, N, S, and H are from the parametrized relationships of Chi et al.,Reference Chi, Dasgupta, Duncan and Shimizu71 Grewal et al.,Reference Grewal, Dasgupta, Holmes, Costin and Li105 Boujibar et al.,Reference Boujibar106 and Clesi et al.,Reference Clesi107 respectively. Solubility constant data for C, N, S, and H in the silicate melt are from Armstrong et al.,Reference Armstrong, Hirschmann, Stanley, Falksen and Jacobsen108 Roskosz et al.,Reference Roskosz, Bouhifd, Jephcoat, Marty and Mysen109 O’Neill and Mavrogenes,Reference O’Neill and Mavrogenes110 and Hirschmann et al.Reference Hirschmann, Withers, Ardia and Foley111

2.4.2 The Role of Post-core Formation Sulfide Segregation

The segregation of sulfide melt has been shown to be a necessary process for explaining the HSE geochemistry of the BSE.Reference Rubie112 This would also deplete the BSE in terms of its S inventory, but owing to the low solubility of C in sulfide-rich melts (Figure 2.6a)Reference Tsuno, Grewal and Dasgupta17, Reference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Corgne, Wood and Fei113Reference Zhang, Hastings, Von der Handt and Hirschmann116 caused by strong non-ideality of mixing between C- and S-rich components in Fe-rich alloy melts, such sulfide segregation should also elevate the C/S ratio of the BSE. Hence, if late accretion establishes a chondritic C/S ratio, sulfide segregation should alter such a ratio. N also is not compatible in S-rich Fe-alloy melts,Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 hence the S/N ratio of the BSE would also drop after sulfide melt segregation. Similarly, it is also unclear what the effects of late-stage sulfide melt segregation might be on the H budget of the BSE. If H is stored in nominally anhydrous silicate minerals in solid, post-core formation BSE, partitioning of H between segregating sulfide melts and solid silicate minerals could influence the H inventory and hence the C/H ratio of the BSE. If H partitions preferentially into sulfide liquid over relevant nominally anhydrous silicate minerals, then the C/H ratio may also increase above the chondritic value.

Figure 2.6 Experimental data showing the effects of S and Si contents in the Fe-rich alloy melt on C solubility. (a) Carbon solubility decreases monotonically with increases in sulfur content in the alloy melt. Si- and N-bearing alloys have lower C solubility for a given S content in the alloy melt. (b) Carbon solubility decreases linearly with increases in Si content in the alloy or decreases in oxygen fugacity (fO2) of silicate–alloy systems. The presence of S in the alloy does not have a major effect on carbon solubility in silicon-bearing alloys.

2.4.3 The Role of MO–Atmosphere Interactions and Atmospheric Loss

Because C and other LEVEs are atmophile elements, they were likely heavily concentrated in the proto-atmospheres at various stages of MO evolution and terrestrial accretion. Hence, it is worth considering under what circumstances silicate magma–atmosphere interactions and possible atmospheric loss through impact-driven blow-offReference Schlichting, Sari and Yalinewich23 could influence the inventory of LEVEs in the BSE. The timing of LEVE delivery is of critical importance in determining the role of silicate–atmosphere interactions (and atmospheric blow-off). For example, if a large fraction of core formation takes place before the LEVEs were available to the growing Earth,Reference O’Brien, Walsh, Morbidelli, Raymond and Mandell20, Reference Raymond, Quinn and Lunine118 then such core formation would have had negligible influence on the BSE LEVE budget, with MO–atmosphere processes being more important. On the other hand, if LEVEs were delivered throughout the terrestrial accretion history,Reference Marty13, Reference Sarafian, Nielsen, Marschall, McCubbin and Monteleone51, Reference Hallis119 then both atmosphere–MO interactions and metal–silicate equilibration could have played roles in establishing the LEVE budget of the BSE.

If the BSE C–N–S–H budgets were established chiefly by loss of proto-atmosphere overlying MO, then N needed to be lost preferentially to C and C needed to be lost preferentially to H. In other words, the proto-atmosphere overlying MO has to be N rich relative to C, C rich relative to H, and with similar abundances of C and S. Indeed, the N-depleted nature of Earth has been explained previously by loss of an N-rich atmosphere.Reference Tucker and Mukhopadhyay120 But evaluation of whether an atmosphere with a concentration of LEVEs in the order N > C~S > H can be generated overlying a silicate magma reservoir depends on the mixed-volatile solubility in shallow MO as a function of oxygen fugacity (fO2; effective partial pressure of oxygen). HirschmannReference Hirschmann11 considered a simple model evaluating the fractionation of C, N, H, and S between silicate MO, overlying atmosphere, and equilibrating core-forming liquid alloy with three distinct fO2 conditions at the MO surface (i.e. IW – 3.5 (reduced), IW – 2 (intermediate fO2), and IW + 1 (oxidized), where IW refers to the log fO2 set by equilibrium of metallic iron and iron oxide (FeO)). Among these, the most oxidizing condition is not very realistic during or soon after core formation because, for a well-mixed MO with a near-constant Fe3+/FeT ratio, the fO2 gradient is such that the shallow MO is always more reduced.Reference Kress and Carmichael121Reference Zhang, Hirschmann, Cottrell and Withers123 Thus, if a well-mixed MO is close to equilibrium with an Fe-rich core at its base, its surface should be more reduced and therefore, for MO–atmosphere interactions, the relevant fO2 would be less than IW. The only way this may be different is if at higher pressures the fO2 gradient reverses (reduction with increasing depth) and iron disproportionation and subsequent Fe–metal segregation leads to oxidation of metal-free MO.Reference Schaefer and Elkins-Tanton124 Solubility data for most LEVEs are lacking or sparse for peridotitic or ultramafic silicate melts due to difficulty in obtaining a glass and reliable chemical analyses for such compositions.Reference Duncan, Dasgupta and Tsuno125 Based on the data available for silicic to mafic silicate liquids, for log fO2 < IW – 0.5 and with decreasing fO2, N becomes much more soluble (reaching the weight percentage level chiefly as N–H speciesReference Dalou, Hirschmann, von der Handt, Mosenfelder and Armstrong85, Reference Libourel, Marty and Humbert126, Reference Kadik127) than C (tens of ppm by weight mostly as CO32– and to some extent as C–H and C–N speciesReference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Dasgupta, Chi, Shimizu, Buono and Walker99, Reference Armstrong, Hirschmann, Stanley, Falksen and Jacobsen108, Reference Li, Dasgupta and Tsuno117, Reference Duncan, Dasgupta and Tsuno125, Reference Holloway, Pan and Gudmundsson128, Reference Li, Dasgupta and Tsuno129), with the latter mostly decreasing with decreasing fO2 and then leveling off with stabilization of C–H species. Similar to N, H also remains quite soluble in silicate melt in the form of OH and H2O even in fairly reducing conditions such as log fO2 of IW – 1 to IW – 2, with molecular hydrogen, H2, contributing to some degree.Reference Hirschmann, Withers, Ardia and Foley111 Therefore, at log fO2 < IW – 1.5, atmospheric blow-off could lower the C/H ratio from a chondritic value, but this would result in an MO C/N ratio that is subchondritic. Similarly, S is much more soluble in mafic–ultramafic silicate melts (several thousands of ppmReference O’Neill and Mavrogenes110, Reference Ding, Hough and Dasgupta130, Reference Baker and Moretti131) compared to C and hence loss of an atmosphere would result in a residual silicate MO with a C/S ratio distinctly lower than that of chondritic value. Therefore, no condition exists where a proto-atmosphere overlying an MO can attain all of the compositional characteristics necessary to leave an MO with C/N, C/S, and C/H ratios of the modern-day BSE (Figure 2.5). Hence, late accretion of chondritic materials to an alloy-free MO followed by an atmospheric loss cannot be the chief origin of BSE LEVEs.

An alternative suggestionReference Tucker and Mukhopadhyay120 is for atmospheric loss to help define the BSE C/N and C/H ratios, but at conditions where C would be retained in the form of carbonate or bicarbonate ions either as dissolved aqueous species or in the crust, mediated by the presence of a liquid water ocean. According to this model, a superchondritic C/N ratio and a subchondritic C/H ratio could have been attained via loss of an N-rich atmosphere overlying an ocean. This mode of atmospheric loss is supported by the dynamical model of a giant impact that shows that such an impact would preferentially remove the atmosphere relative to the ocean.Reference Genda and Abe132 However, it is unclear what the water chemistry would be in these oceans, which, if they existed, is expected to be highly anoxic and thus may not sequester significant dissolved carbonates.

2.5 Establishing the Volatile Budget of the BSE through Equilibrium Accretion and MO Differentiation

The discussion on the origin of LEVEs in the BSE thus far considered building blocks and processes after the core formation was mostly complete. This approach would be appropriate if all of the LEVEs were delivered through late accretions. The discovery of a solar component (e.g. neon) in deep plume-derived magma,Reference Honda, McDougall, Patterson, Doulgeris and Clague133, Reference Sarda, Staudacher and Allègre134 however, suggests that Earth might have trapped nebular gas early within the first few millions of years of the formation of the Solar System. Indeed, the possibility of acquiring volatile elements such as primordial noble gases directly from nebular gas has been proposed in a number of studies.Reference Mizuno, Nakazawa and Hayashi135, Reference Yokochi and Marty136 Although nebular gas can persist for ∼10 Myr, its median lifetime (i.e. the time at which half of all systems lost their nebular gas disks) is ∼2.5 Myr.Reference Lynne137 Given the mean time of terrestrial accretion (i.e. the time to have grown to 50% of the final mass) is ∼11 MyrReference Yin138 while the mean accretion timescale of Mars is 1.90.8+1.7 Myr,Reference Dauphas and Pourmand67, Reference Tang and Dauphas139 it is likely that although the entire proto-Earth mass may not have equilibrated with the solar nebula, at least a few Mars-sized embryos that contributed to the initial proto-Earth mass likely inherited nebular volatiles. These constraints suggest that early accreting material may have been LEVE poor but likely not LEVE free. However, the LEVE compositions do not suggest direct Solar contribution to be important (Figure 2.2).

The 182Hf/182W ratio of iron meteorites suggests metal–silicate separation in protoplanetary bodies occurred as early as ~1 Myr after the formation of the Solar System. Accordingly, core formation would influence the initial distribution of LEVEs among the planetary reservoirs.Reference Dasgupta8, Reference Hirschmann11, Reference Dasgupta and Walker14, Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16, Reference Tsuno, Grewal and Dasgupta17, Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Dalou, Hirschmann, von der Handt, Mosenfelder and Armstrong85, Reference Dasgupta, Chi, Shimizu, Buono and Walker99, Reference Li, Dasgupta and Tsuno117 Over the past decade, a significant body of research has tried to understand the fate of C and other LEVEs (S, N, and H) during core formation processes and associated alloy–silicate equilibration.Reference Dasgupta and Walker14Reference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Dalou, Hirschmann, von der Handt, Mosenfelder and Armstrong85, Reference Dasgupta, Chi, Shimizu, Buono and Walker99, Reference Boujibar106, Reference Roskosz, Bouhifd, Jephcoat, Marty and Mysen109, Reference Li, Dasgupta and Tsuno117, Reference Rose-Weston, Brenan, Fei, Secco and Frost140Reference Suer, Siebert, Remusat, Menguy and Fiquet142 These investigations chiefly focused on determining the partition coefficients of LEVEs between Fe-rich alloy melts and coexisting silicate melts (i.e. Dxalloy/silicate (= concentration of an element in the alloy melt divided by concentration of the same element in the equilibrium silicate melt), where x = C, N, S, etc.). If Dxalloy/silicate is >1 and if a given LEVE, x, is available during core formation, then it should be preferentially sequestered in the metallic core. On the other hand, if Dxalloy/silicate<1, the LEVE, x, present during core–mantle fractionation should preferentially concentrate in the silicate Earth. Dxalloy/silicate values provide the framework to understand how a given LEVE is expected to be distributed between the metallic and silicate portions of a differentiated planet if alloy and silicate melts equilibrated before separation.

The most robust constraints available on the alloy–silicate partitioning of LEVEs are obtained through high pressure–temperature experiments. Determinations of DCalloy/silicate at high pressure–temperature exist chiefly for graphite/diamond-saturated conditions (i.e. experiments conducted in graphite capsules). Concentrations of C in Fe-rich alloy melts at graphite/diamond saturation and in the absence of any other light nonmetals show small dependence on pressure and temperature (not shown)Reference Dasgupta and Walker14, Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 and are ~5–7 wt.%. However, C solubility in Fe-rich alloy melts decreases strongly with increasing S and Si content, and at S of content 22–36 wt.% or Si content >10–12 wt.%, C solubility is <1 wt.% (Figure 2.6a and b). Carbon solubility in Fe-rich alloy melts also diminishes with increasing Ni content in the S-free or S-poor alloy.Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Tsuno and Dasgupta115, Reference Zhang, Hastings, Von der Handt and Hirschmann116, Reference Tsymbulov and Tsmekhman143 Unlike in alloy melts, C solubility in equilibrium silicate melts is tens to hundreds of ppm and mostly decreases with decreasing fO2 (Figure 2.7). The MOs of Earth and other inner-Solar System bodies are expected to be ultramafic to mafic; that is, relatively poor in silica and rich in MgO. Because of the difficulty of rapidly cooling such depolymerized silicate liquids from well-constrained high pressures–temperatures to glasses at ambient conditions, which is necessary for reliable quantitative analyses of carbon and hydrogen, C solubility data in such ultramafic–mafic liquids at core-forming reduced conditions are sparse.Reference Duncan, Dasgupta and Tsuno125 To address this, experimental studies have been conducted with a silicate melt of varying polymerization so that a meaningful extrapolation of C solubility and DCalloy/silicate to desired silicate melt composition can be made.

Figure 2.7 Experimental data showing the effects of fO2 and water content in the silicate melt on C solubility. Carbon solubility decreases with decreasing log fO2 until ~ΔIW of –2 to –1 followed by an increase below ~IW – 2 depending upon bulk water content in the silicate melt. This change is likely caused by the change in contributions of C as CO32−CO32 at higher fO2 and C as C–H species such as methyl groups at lower fO2. The pressure–temperature–composition (expressed as NBO/T) condition of the experimental glasses are given in the legend. The shaded region represents predicted solubility of C as CO32−CO32 at graphite saturation based on the model of Eguchi and DasguptaReference Eguchi and Dasgupta146 for the entire pressure, temperature, and XH2O space of the experimental data, while the NBO/T parameter is extrapolated for peridotite-like silicate melts (NBO/T = 2.7).

2.5.1 Carbon Speciation in MO

DCalloy/silicate and C solubility in silicate melts are affected by the incorporation mechanism of C in the silicate melt structure. Irrespective of the saturating phase such as CO2-bearing vapor or graphite/diamond, carbon can dissolve in silicate melt as molecular CO2 (CO2mol.CO2mol.Reference Stolper, Fine, Johnson and Newman144Reference Duncan and Dasgupta148), as carbonate anions (CO32−CO32Reference Duncan, Dasgupta and Tsuno125, Reference Stolper and Holloway149, Reference Pan, Holloway and Hervig150) bonded to cation modifiers, possibly as molecular CO or carbonyl,Reference Yoshioka, McCammon, Shcheka and Keppler151, Reference Wetzel, Rutherford, Jacobsen, Hauri and Saal152 and as possible C–H, C–N complexes.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Li, Dasgupta and Tsuno117 For mafic–ultramafic melts, the chief species of interest is CO32−CO32, which complexes most strongly with Ca2+ cations and also with Na+, Mg2+, and maybe K+.Reference Duncan, Dasgupta and Tsuno125, Reference Eguchi and Dasgupta146, Reference Dasgupta, Hirschmann and Smith153Reference Ghiorso and Gualda156 Most experiments on carbon dissolution at reducing conditions employed mafic compositions of variable degrees of polymerization (i.e. tholeiitic to alkali basalts or even more silicic melts). The most common approach to extrapolating C solubility and DCalloy/silicate to systems containing peridotitic compositions is to use a compositional parameter NBO/T, where “NBO” is the total concentration of non-bridging oxygen in the silicate melt and “T” is the total number of tetrahedral cations. Although this method is thought to give reasonable predictions of ultramafic melt C solubility as CO32−CO32, Duncan et al.Reference Duncan, Dasgupta and Tsuno125 showed that the NBO/T approach may overpredict the extent of CO32−CO32 dissolution in silicate melts if an increase of NBO/T is correlated with an increase in CaO rather than MgO in melts. This is because Ca2+ cations complex with CO32−CO32 much more than Mg2+,Reference Eguchi and Dasgupta146 although peridotitic melts are still predicted to have higher CO32−CO32 dissolution capacities than basalts. Figure 2.7 shows CO32−CO32 dissolution as a function of fO2 and shallow MO conditions. Also shown in Figure 2.7 are other possible species such as C–H that may become important at low fO2.

2.5.2 DCalloy/silicate and Its Impact on Carbon Distribution between BSE versus Core in Various Scenarios of Equilibrium Core Formation

Figure 2.8 provides the available data of DCalloy/silicate, showing their dependence on alloy composition and fO2 at shallow MO conditions. Shallow MO DCalloy/silicate at log fO2 of IW – 2 to IW – 0.5 are in the order of ~103, with the values increasing with increasing pressure and decreasing fO2 and decreasing with increasing temperature and melt depolymerization.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16, Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Li, Dasgupta and Tsuno117 In more reducing conditions where the Si content becomes significant in the equilibrium Fe-rich alloy melts and where decreasing C dissolution as CO32−CO32 is offset to some extent by the increasing fraction of C–H complexes (Figure 2.7), DCalloy/silicate may be somewhat lower and may not monotonically increase with decreasing fO2 (Figure 2.8b).Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16, Reference Li, Dasgupta and Tsuno117 Application of empirical parameterizations of DCalloy/silicate as a function of temperature, pressure, fO2, silicate melt composition, and alloy composition (e.g. Ni content) to conditions that many studies considered relevant for the deep terrestrial MO (e.g. pressure = 45–65 GPa, temperature = 3500–4000 K, log fO2 = IW – 2) yields values in the order of 103–104. The existing constraints on DCalloy/silicate in graphite-saturated conditions therefore suggest that carbon is highly siderophile. If core formation was an equilibrium process (i.e. if the entire terrestrial core mass equilibrated with a global silicate MO), post-core formation silicate Earth would be severely depleted in carbon (≪10 ppm), with the extent of depletion being dependent upon the available bulk carbon that participates in such core–mantle fractionation.Reference Dasgupta8, Reference Tsuno, Grewal and Dasgupta17, Reference Chi, Dasgupta, Duncan and Shimizu71, Reference Dasgupta, Chi, Shimizu, Buono and Walker99 This analysis does not imply that the core formation process was a single-stage alloy–silicate equilibration event, but rather what the C distribution between Earth’s metallic and nonmetallic portions would be for an average condition of core–mantle separation.

Figure 2.8 Experimental data showing the effects of S and Si in Fe-rich alloy melts on DCalloy/silicate. (a) DCalloy/silicate decreases with increases in S content in the alloy melt. (b) DCalloy/silicate decreases with increases in Si content in the alloy or indirect decreases in fO2 of the alloy–silicate systems. The presence of S in Si-bearing alloy melts generally decreases DCalloy/silicate in comparison to S-free systems.

Figure 2.9 shows more realistic scenarios of alloy–silicate equilibration and subsequent core segregation in terrestrial MO where the depth, temperature, and fO2 evolve with the extent of accretion.Reference Wade and Wood158 This style of core–mantle fractionation is described as multistage, equilibrium core formation. Competing arguments exist in the literature as to whether, with continuing accretion, proto-Earth’s MO evolved from a relatively reduced to an oxidized conditionReference Wade and Wood158, Reference Rubie159 or a relatively oxidized to a reducing condition.Reference Badro, Brodholt, Piet, Siebert and Ryerson160 Figure 2.9a therefore shows various plausible fO2 paths of MO evolution of a growing Earth and computes the change of DCalloy/silicate along those paths, following Li et al.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 It must be noted that the current data sets of DCalloy/silicate do not show any effects of pressure–temperature in relatively reducing conditions (e.g. log fO2 < IW – 1.5). In such reducing conditions, the controlling variables appear to be fO2, the Si content of the metal, and the water content of the silicate melt. At ~IW – 2 < log fO2 < IW, however, the key variables controlling DCalloy/silicate are pressure, temperature, fO2, and silicate melt composition. Figure 2.9b, however, shows that irrespective of MO condition, DCalloy/silicate remains ≫1, and in fact deeper MO, toward the later stage of accretion, leads to even higher DCalloy/silicate predictions than those experimentally measured at shallow MO conditions. Therefore, the multistage core formation model also leads to C being stripped off highly efficiently from the silicate MO, leaving behind the BSE that is <1–10 ppm for a bulk Earth with 1000 ppm C (Figure 2.9c).

Figure 2.9 The effects of different accretion scenarios on the core–mantle partitioning behavior of carbon. (a) Four possible evolution paths for the FeO content of the mantle with its growth. The “very reduced” and “reduced” paths simulate the accretion of progressively oxidized material with time, and the “oxidized” path simulates accretion of progressively reduced material with time. All of the models converge at log fO2 of IW – 1.8, or present-day FeO content (8 wt.%) of the primitive upper mantle. (b) DCalloy/silicateat each step of accretion (i.e. accretion of 1% ME), calculated using parametrized relationships from Li et al.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 For the “very reduced” and “reduced” paths, where methyl (or alkyl) groups would be the predominant dissolved carbon species in the silicate melt, eq. (4) is chosen from Li et al.,Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 whereas for the “intermediate” and “oxidized” paths, where carbonate groups would be the predominant dissolved carbon species in the silicate melt, eq. (3) is chosen from Li et al.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16 (c) Post-core formation MO resulting from a bulk system with 1000 ppm C and assuming 100% alloy–silicate equilibration at each accretional step. All of the models lead to MO with distinctly less C that those estimated for the BSE. (d) Carbon content in Earth’s core converges at ~3300 ppm at the end of all accretional models.

2.5.3 Comparison of DCalloy/silicate with Alloy–Silicate Melt Partitioning of Other LEVEs

The preceding section argues that equilibrium core formation either via homogeneous or heterogeneous accretion results in a nonmetallic Earth that is too depleted in C compared to the present-day BSE. Explaining the origin of LEVEs in the BSE through equilibrium core formation is exacerbated when DCalloy/silicate is compared with the alloy–silicate partition coefficients of other LEVEs such as S, N, and H.

The partition coefficient of sulfur between alloy and silicate melts, DSalloy/silicate, has been experimentally constrained in a number of studies,Reference Boujibar106, Reference Rose-Weston, Brenan, Fei, Secco and Frost140, Reference Suer, Siebert, Remusat, Menguy and Fiquet142 but only a few recent studies measured DSalloy/silicate and DCalloy/silicate simultaneously.Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 These studies show that over the entire range of shallow MO conditions explored experimentally thus far, DCalloy/silicateDSalloy/silicate, with DCalloy/silicate/DSalloy/silicate varying between ~1000 and 20 as log fO2 varies from IW – 4.5 to IW – 0.5. Therefore, under any conditions of alloy–silicate equilibration, C is more siderophilic than S and equilibrium core formation leaves a subchondritic C/S ratio for the silicate MO.

Experiments also exist that constrained DNalloy/silicate.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Roskosz, Bouhifd, Jephcoat, Marty and Mysen109, Reference Kadik127 Among these, the recent studies that constrain both DCalloy/silicate and DNalloy/silicate from the same experimentsReference Grewal, Dasgupta, Sun, Tsuno and Costin18, Reference Dalou, Hirschmann, von der Handt, Mosenfelder and Armstrong85 show that nitrogen is only mildly siderophile in relatively oxidizing conditions and possibly even lithophile in relatively reducing conditions (e.g. log fO2 < IW – 3). Therefore, DCalloy/silicate/DNalloy/silicate is ~6–120 at 1–7 GPa, 1400–1800°C, and fO2 of IW – 3.5 to IW – 0.5, and hence any models of equilibrium core formation would also produce a subchondritic C/N ratio for the post-core formation BSE.

Alloy–silicate partitioning studies for H are limited, and two available studiesReference Clesi107, Reference Okuchi161 provide contrasting data. OkuchiReference Okuchi161 estimates the DHalloy/silicate value at a fixed pressure of 7.5 GPa and as a function of temperature, but does not provide any pressure effect. Based on the experiments of Okuchi,Reference Okuchi161 H is mildly siderophile. On the other hand, Clesi et al.Reference Clesi107 conducted experiments at 5, 10, and 20 GPa with C-rich alloy melts and measured DHalloy/silicate of 0.04–0.77. Therefore, it remains unclear whether H is retained mostly in the core of differentiated planets or in the silicate fraction. Despite the discrepancy in the available experimental data on DHalloy/silicate, the BSE could attain a subchondritic C/H ratio via an equilibrium core formation process where much more C would be sequestered to the metallic core than H (Figure 2.4). However, if the DHalloy/silicate data of Clesi et al.Reference Clesi107 are appropriate, then equilibrium core formation would lead to a BSE C/H ratio that is much lower than what is estimated for the present-day BSE.

2.6 C and Other LEVE Budgets of the BSE: A Memory of Multistage Accretion and Core Formation Process with Partial Equilibrium?

The discussion presented thus far suggests that no single process of volatile delivery such as via late accretion of any known class of meteorites or comets, through loss of an atmosphere overlying a MO,Reference Schlichting, Sari and Yalinewich23 being subjected to equilibrium core formation processes, or sulfide segregation to the coreReference Rubie112, Reference O’Neill162 can simultaneously explain carbon’s absolute as well as relative abundance compared to other LEVEs. This does not mean that these processes did not impact the volatile inventory of the BSE, and it is plausible that all of these processes acted in tandem. However, the nature of the fractionation of C with respect to H, N, and S for each of these processes is such that even if they acted together, not all of the attributes of the BSE volatile budget can be satisfied. This brings us to consider other styles of terrestrial accretion that may have been essential in establishing the LEVE budget of modern Earth.

2.6.1 Disequilibrium Core Formation

Models of protoplanetary disk dynamics, dust and planetesimal accretion, and planetary migration show that although the first 50% of proto-Earth mass growth was continuous, rapid, and completed within the first few million years of Solar System formation, the protracted growth of Earth spanning ~100 MyrReference Kleine163 involved few abrupt increases in the accreted mass via several giant impacts of Mars-sized or perhaps even larger bodies.Reference Rubie164 The final such impact is thought to be the Moon-forming giant impact.Reference Canup and Asphaug165, Reference Ćuk and Stewart166 A key aspect of accretion via these impacts is that proto-Earth collided with differentiated planetary embryos. Merger of proto-Earth with large differentiated bodies is thought to have taken place via near-disequilibrium merger of the core of the impactor with proto-Earth’s core owing to limited emulsification of the impactor’s core and the lack of complete mixing of the core with the entrained silicate melt.Reference Deguen, Olson and Cardin167, Reference Deguen, Landeau and Olson168 Therefore, core formation involving giant impacts may leave a silicate fraction with a very different LEVE budget than that predicted by equilibrium core formation models. One critical question is whether any differentiated, large impactor could bring LEVE attributes that would help proto-Earth to evolve to the present-day BSE. A full-parameter space involving not only the LEVEs but also other nonvolatile, highly siderophile, and moderately siderophile elements, major elements, and stable isotope anomalies, as well as various accretion and core formation model assumptions, is not yet explored. However, recent studiesReference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 showed how having an S-rich or Si-rich Fe-alloy in the cores of planetary embryos allow attaining high C/N and C/S ratios in the mantles of such embryos. This would occur as the cores of these embryos reach the C solubility limit and expel graphite or diamond to the overlying silicate fraction (Figure 2.10).Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 Therefore, a merger of the silicate reservoirs of such embryos with the proto-silicate Earth reservoir could give rise to the BSE LEVE inventory. Using inverse Monte Carlo simulations, Grewal et al.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 showed that if the merger of the impactor is in perfect disequilibrium, depending on the core size of the impactor and bulk S content, the impactor could be similar in size as Mars (~0.1 ME; Figure 2.10). It is also showed that the bulk C content of such an impactor would not need to be enriched in carbon and could have a bulk C concentration (~1000–5000 ppm C)Reference Tsuno, Grewal and Dasgupta17, Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 similar to those of E-chondrites. These model results are intriguing because these potentially circumvent the long-standing challenge of reconciling the volatile element isotope constraintsReference Sarafian, Nielsen, Marschall, McCubbin and Monteleone51, Reference Marty and Yokochi59 and constraints from the isotopic anomalies of lithophile and siderophile elementsReference Dauphas31, Reference Dauphas32 in the major building blocks of Earth. Some of these latest studies on the origins of terrestrial volatiles through late-stage accretion of a differentiated planetary embryo are also conceptually similar to models that call for the addition of a highly reduced rocky body to the growing Earth to potentially explain the superchondritic Sm/Nd and 142Nd/144Nd anomalies observed in Earth’s present mantle.Reference Wohlers and Wood169, Reference Wohlers and Wood170 The concept of bringing carbon and other LEVEs via accretion of the mantle of a differentiated embryo with little post-impact equilibration is also similar to the suggestion that Earth’s Nb/Ta may be explained by accretion of a differentiated asteroid with a reduced S-rich core that efficiently sequestered Nb.171 Grewal et al.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 also posited that Earth’s volatile-delivering impactor may have been the Moon-forming giant impactor, given the geochemical similarity, including volatiles, between the BSE and the Moon. This is again broadly consistent with Wade and Wood,Reference Wade and Wood172 who demonstrated that near-disequilibrium accretion of a reduced planetary embryo of 0.1–0.2 ME could explain the isotopic constraints, FeO contents, and Nb/Ta and Hf/W ratios of the silicate Earth and Moon.

Figure 2.10 A schematic depiction of the late-stage accretion for a large rocky planet such as Earth where a large fraction of mass increase takes place via giant impacts of small planets and planetary embryos (0.1–0.2 ME).Reference Li, Dasgupta, Tsuno, Monteleone and Shimizu16Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 This mode of accretion creates a unique possibility of establishing the LEVE inventory in the degassable portion of the rocky planets.

2.7 Carbon as a Light Element in the Core

One by-product of our discussion thus far is constraint on the concentration of C in the core. One of the main uncertainties in constraining the C budget of Earth’s core is the bulk C content in the growing Earth through time in the first ~100–120 Myr of the Solar System. It is not realistic to conceive that Earth’s building blocks and therefore the bulk C content did not change through the accretion history, but for simplicity here we assume a fixed C budget for the entire history of core formation or MO equilibration. If core formation was achieved by a single-stage equilibrium event at a deep MO, although not a physically realistic process, C content in the core as shown by DasguptaReference Dasgupta8 would be ~0.2 wt.% for a bulk Earth C content of ~730–1000 ppm. For a multistage core formation, with evolving depth–temperature–fO2 conditions of alloy–silicate equilibration, the core would finally contain ~0.34 wt.% C (Figure 2.9d). Finally, in the model where the C–N–S budget of the BSE is attained via late-stage merger of a large differentiated planetary embryo, the impactor’s S-rich core contributes ~100–700 ppm C to Earth’s core.Reference Tsuno, Grewal and Dasgupta17, Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 Therefore, depending on the volatile budget of the growing Earth leading to the giant impact, from a completely C-free proto-Earth to a proto-Earth with 1000 ppm bulk C, the final core C content would be between 700 ppm and 0.4 wt.%. The upper bound derived here is similar to the previous constraintReference Dasgupta and Walker14 placed using the mantle C abundance. Wood et al.Reference Wood, Li and Shahar60 also attempted to place constraint on the C content of Earth’s core based on the combined effects of C and S in the alloy on alloy–silicate partitioning of W and Mo, with the chief effect being W and Mo becoming more siderophile with increasing C in the alloy melt. By coupling S and C delivery in the final 13% of accretion, Wood et al.Reference Wood, Li and Shahar60 showed that in order to obtain the known Mo and W abundances of the mantle, the bulk C content of the core needs to be <1 wt.%. Badro et al.Reference Badro, Brodholt, Piet, Siebert and Ryerson160 also reached a similar conclusion that, in order to simultaneously satisfy the geochemical and geophysical constraints, Earth’s core composition cannot contain C as one of the major light elements. Therefore, all of the constraints thus far point to Earth’s core being the largest C reservoir, yet C being a minor light element in the core.

2.8 Conclusion

The origins of Earth’s LEVEs, including C, remain debated. Such debate stems in part from the uncertainty regarding the whole-Earth budget of C and other volatile elements and in part from the fate of these elements in putative planetary building blocks throughout their journey in the protoplanetary disk and planet formation process. While the abundance of each LEVE in the BSE is poorly constrained, the ratios of various LEVEs (along with their isotopic compositions) are more reliable geochemical indices to isolate key processes. Some of the key processes to consider are the accretion of undifferentiated chondritic meteorites, compound-specific condensation of the LEVEs, fractionation of the LEVEs during alloy–silicate differentiation, and MO–proto-atmosphere partitioning aided by redox-controlled solubility. Probing each of these processes in the light of the evolution of the LEVE ratios suggests that no single process or combination of processes can result in the currently estimated abundance of LEVEs in the BSE. However, a new paradigm is emerging where the late-stage growth of Earth via the impact and accretion of differentiated planetary embryos (0.1–0.2 ME) may not only explain the abundances of C–N–S, but also satisfy the broad geochemical similarity of Earth and the Moon. Despite the similarity in the arguments put forward by many recent studies that call for the accretion of one or more large planetary embryos with little post-impact equilibration of the impactor’s core with the terrestrial mantle/MO, differences remain in their details. For example, Grewal et al.Reference Grewal, Dasgupta, Sun, Tsuno and Costin18 proposed accretion of a relatively oxidized planet with an S-rich core, whereas Wade and WoodReference Wade and Wood172 called for accretion of a reduced planet with S-rich core. Future studies, therefore, will have to investigate how to reconcile these differences while fulfilling all of the geochemical constraints, both for non-volatile lithophile and siderophile elements and LEVEs. Nonetheless, it appears inevitable that a large planet like Earth’s LEVE budgets were shaped by its protracted growth history with episodes of planetesimals, planetary embryos, and planet accretion, whereas other smaller planets that grew quickly lacked some of these possible supply and modification mechanisms of key ingredients. Whatever the accretion mechanism may be, as long as mostly volatile-depleted materials accreted on Earth throughout its core formation history, C as well as N and H are going to be minor elements in the bulk core. As the search continues to find other rocky, habitable worlds beyond our Solar System, the origins of C–H–N–S on Earth provide a guiding framework.

2.9 Limits of Knowledge and Unknowns

The limits of knowledge on the origin of Earth’s C and other LEVEs derive from current unknowns about many fundamental questions regarding the formation and evolution of Earth and other rocky planets in the Solar System, as well as about the compositions and processes in Earth’s deep interior. The condensation temperature of C that dictates the formation of solid C-bearing compounds from gases is thought to be extremely low,Reference Krasnokutski173 leading to the idea that most C brought to the inner Solar System at the early stage of planet formation is largely lost from the building blocks. However, it remains unknown whether some of the refractory C-bearing phases such as graphite and carbide could have contributed to the C budget of Earth during its early growth stage.

Experimental data on alloy–silicate elemental partitioning for mixed-volatile systems, especially in conditions approaching the deep MO of Earth, are still lacking and are needed in order to refine some of the current ideas based on extrapolated values from relatively low pressure–temperature experiments. Similar experimental data are needed for volatile isotopic systems as well, which are even more restricted,Reference Li, Marty, Shcheka, Zimmermann and Keppler174, Reference Labidi175 and for many systems, such as C,Reference Satish-Kumar, So, Yoshino, Kato and Hiroi64 H, and S,Reference Labidi175 they are either nonexistent or do not extend to the pressure–temperature and compositional space directly relevant for MOs.

Acknowledgments

The authors thank many former and current members of Rice University’s experimental petrology team, including Han Chi, Yuan Li, Megan Duncan, Alexandra Holmes, Kyusei Tsuno, and James Eguchi, whose research results influenced the views expressed in this chapter. The authors thank Jonathan Tucker and coeditor Isabelle Daniel for useful reviews. This work was made possible by funding from NASA grants 80NSSC18K0828 and 80NSSC18K1314, the Deep Carbon Observatory, and a Packard Fellowship for Science and Engineering to RD.

Questions for the Classroom

1 How well do the chondritic meteorites (undifferentiated materials; i.e. mixtures or metals, sulfides, and silicates), derived from the main asteroid belt, represent the building blocks of Earth and other rocky planets and planetary embryos in the inner Solar System? In other words, should the plausible compositions of building blocks including their LEVE budgets remain restricted to the known classes of meteorites?

2 How could we improve our knowledge on the concentration of C, N, S, and H in the BSE? What other approaches can be used to determine the mantle inventory of LEVEs?

3 How can we improve our understanding of the relative contributions of giant impacts versus impacts of planetesimals on the acquisition of LEVEs for large, rocky planets?

4 What information from our knowledge on the origins of LEVEs on Earth can be applied to build a framework for exoplanetary chemical habitability?

References

Birch, F., Elasticity and constitution of the Earth’s interior. J. Geophys. Res. 57, 227286 (1952).CrossRefGoogle Scholar
Birch, F., Density and composition of mantle and core. J. Geophys. Res. 69, 43774388 (1964).Google Scholar
Poirier, J.-P., Light elements in the Earth’s outer core: a critical review. Phys. Earth Planet. Int. 85, 319337 (1994).CrossRefGoogle Scholar
Javoy, M., Pineau, F., & Allègre, C.J., Carbon geodynamic cycle. Nature 300, 171173 (1982).Google Scholar
Sleep, N.H. & Zahnle, K., Carbon dioxide cycling and implications for climate on ancient earth. J. Geophys. Res. 106, 13731399 (2001).Google Scholar
Hayes, J.F. & Waldbauer, J.R., The carbon cycle and associated redox processes through time. Phil. Trans. Royal Soc. London B361, 931950 (2006).Google Scholar
Hirschmann, M.M. & Dasgupta, R., The H/C ratios of Earth’s near-surface and deep reservoirs, and consequences for deep Earth volatile cycles. Chem. Geol. 262, 416 (2009).Google Scholar
Dasgupta, R., Ingassing, storage, and outgassing of terrestrial carbon through geologic time. Rev. Min. Geochem. 75, 183229 (2013).CrossRefGoogle Scholar
Dasgupta, R. & Hirschmann, M.M., The deep carbon cycle and melting in Earth’s interior. Earth Planet. Sci. Lett. 298, 113 (2010).Google Scholar
Bergin, E.A., Blake, G.A., Ciesla, F., Hirschmann, M.M., & Li, J., Tracing the ingredients for a habitable earth from interstellar space through planet formation. Proc. Natl Acad. Sci. 112, 89658970 (2015).CrossRefGoogle ScholarPubMed
Hirschmann, M.M., Constraints on the early delivery and fractionation of Earth’s major volatiles from C/H, C/N, and C/S ratios. Am. Mineral. 101, 540553 (2016).Google Scholar
Halliday, A.N., The origins of volatiles in the terrestrial planets. Geochim. Cosmochim. Acta 105, 146171 (2013).CrossRefGoogle Scholar
Marty, B., The origins and concentrations of water, carbon, nitrogen and noble gases on Earth. Earth Planet. Sci. Lett. 313 –314, 5666 (2012).Google Scholar
Dasgupta, R. & Walker, D., Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle. Geochim. Cosmochim. Acta 72, 46274641 (2008).CrossRefGoogle Scholar
Kuramoto, K. & Matsui, T., Partitioning of H and C between the mantle and core during the core formation in the Earth: its implications for the atmospheric evolution and redox state of early mantle. J. Geophys. Res. 101, 1490914932 (1996).Google Scholar
Li, Y., Dasgupta, R., Tsuno, K., Monteleone, B., & Shimizu, N., Carbon and sulfur budget of the silicate Earth explained by accretion of differentiated planetary embryos. Nat. Geosci. 9, 781785 (2016).Google Scholar
Tsuno, K., Grewal, D.S., & Dasgupta, R., Core–mantle fractionation of carbon in Earth and Mars: the effects of sulfur. Geochim. Cosmochim. Acta 238, 477495 (2018).CrossRefGoogle Scholar
Grewal, D.S., Dasgupta, R., Sun, C., Tsuno, K., & Costin, G., Delivery of carbon, nitrogen and sulfur to the silicate Earth by a giant impact. Science Adv. 5, eaau3669 (2019).Google Scholar
Morbidelli, A., Lunine, J.I., O’Brien, D.P., Raymond, S.N., & Walsh, K.J., Building terrestrial planets. Ann. Rev. Earth Planet. Sci. 40, 251275 (2012).CrossRefGoogle Scholar
O’Brien, D.P., Walsh, K.J., Morbidelli, A., Raymond, S.N., & Mandell, A.M., Water delivery and giant impacts in the “Grand Tack” scenario. Icarus 239, 7484 (2014).Google Scholar
Marchi, S., Canup, R.M., & Walker, R.J., Heterogeneous delivery of silicate and metal to the Earth by large planetesimals. Nat. Geosci. 11, 7781 (2018).Google Scholar
Schlichting, H.E. & Mukhopadhyay, S., Atmosphere impact losses. Space Sci. Rev. 214, 34 (2018).CrossRefGoogle Scholar
Schlichting, H.E., Sari, R.E., & Yalinewich, A., Atmospheric mass loss during planet formation: the importance of planetesimal impacts. Icarus 247, 8194 (2015).Google Scholar
Ahrens, T.J., Impact erosion of terrestrial planetary atmospheres. Ann. Rev. Earth Planet. Sci. 21, 525555 (1993).CrossRefGoogle Scholar
Ringwood, A.E., Chemical evolution of the terrestrial planets. Geochim. Cosmochim. Acta 30, 41104 (1966).CrossRefGoogle Scholar
Anders, E. & Ebihara, M., Solar-system abundances of the elements. Geochim. Cosmochim. Acta 46, 23632380 (1982).Google Scholar
Allègre, C., Manhès, G., & Lewin, É., Chemical composition of the Earth and the volatility control on planetary genetics. Earth Planet. Sci. Lett. 185, 4969 (2001).CrossRefGoogle Scholar
Clayton, R.N., Onuma, N., & Mayeda, T.K., A classification of meteorites based on oxygen isotopes. Earth Planet. Sci. Lett. 30, 1018 (1976).Google Scholar
Clayton, R.N., Oxygen isotopes in meteorites. Ann. Rev. Earth Planet. Sci. 21, 115149 (1993).CrossRefGoogle Scholar
Young, E.D. et al., Oxygen isotopic evidence for vigorous mixing during the Moon-forming giant impact. Science 351, 493496 (2016).Google Scholar
Dauphas, N. et al., Calcium-48 isotopic anomalies in bulk chondrites and achondrites: evidence for a uniform isotopic reservoir in the inner protoplanetary disk. Earth Planet. Sci. Lett. 407, 96108 (2014).CrossRefGoogle Scholar
Dauphas, N., The isotopic nature of the Earth’s accreting material through time. Nature 541, 521 (2017).CrossRefGoogle ScholarPubMed
Fischer-Gödde, M. & Kleine, T., Ruthenium isotopic evidence for an inner Solar System origin of the late veneer. Nature 541, 525 (2017).CrossRefGoogle ScholarPubMed
Javoy, M. et al., The chemical composition of the Earth: enstatite chondrite models. Earth Planet. Sci. Lett. 293, 259268 (2010).CrossRefGoogle Scholar
Javoy, M., The integral enstatite chondrite model of the Earth. Geophys. Res. Lett. 22, 22192222 (1995).Google Scholar
McDonough, W.F. & Sun, S.-S., The composition of the Earth. Chem. Geol. 120, 223253 (1995).Google Scholar
Wasson, J.T. & Kallemeyn, G.W., Compositions of chondrites. Phil. Trans. Royal Soc. A. 325, 535544 (1988).Google Scholar
Fitoussi, C., Bourdon, B., & Wang, X., The building blocks of Earth and Mars: a close genetic link. Earth Planet. Sci. Lett. 434, 151160 (2016).CrossRefGoogle Scholar
Palme, H. & O’Neill, H.S.C., Cosmochemical estimates of mantle composition, in The Mantle and Core, Vol. 2, ed. Carlson, R.W. (Oxford: Pergamon Press, 2007), pp. 138.Google Scholar
Murakami, M., Ohishi, Y., Hirao, N., & Hirose, K., A perovskitic lower mantle inferred from high-pressure, high-temperature sound velocity data. Nature 485, 90 (2012).Google Scholar
Fitoussi, C. & Bourdon, B., Silicon isotope evidence against an enstatite chondrite Earth. Science 335, 14771480 (2012).Google Scholar
Shahar, A. et al., Experimentally determined Si isotope fractionation between silicate and Fe metal and implications for Earth’s core formation. Earth Planet. Sci. Lett. 288, 228234 (2009).Google Scholar
Shahar, A. et al., High-temperature Si isotope fractionation between iron metal and silicate. Geochim. Cosmochim. Acta 75, 76887697 (2011).CrossRefGoogle Scholar
Ballmer, M.D., Houser, C., Hernlund, J.W., Wentzcovitch, R.M., & Hirose, K., Persistence of strong silica-enriched domains in the Earth’s lower mantle. Nat. Geosci. 10, 236 (2017).Google Scholar
Dauphas, N., Poitrasson, F., Burkhardt, C., Kobayashi, H., & Kurosawa, K., Planetary and meteoritic Mg/Si and δ30Si variations inherited from solar nebula chemistry. Earth Planet. Sci. Lett. 427, 236248 (2015).CrossRefGoogle Scholar
McNaughton, N.J., Borthwick, J., Fallick, A.E., & Pillinger, C.T., Deuterium/hydrogen ratios in unequilibrated ordinary chondrites. Nature 294, 639 (1981).Google Scholar
Robert, F., Merlivat, L., & Javoy, M., Deuterium concentration in the early Solar System: hydrogen and oxygen isotope study. Nature 282, 785 (1979).CrossRefGoogle Scholar
Alexander, C.M.O.D. et al., The provenances of asteroids, and their contributions to the volatile inventories of the terrestrial planets. Science 337, 721723 (2012).Google Scholar
Pearson, V.K., Sephton, M.A., Franchi, I.A., Gibson, J.M., & Gilmour, I., Carbon and nitrogen in carbonaceous chondrites: elemental abundances and stable isotopic compositions. Meteor. Planet. Sci. 41, 18991918 (2006).CrossRefGoogle Scholar
Marty, B., Alexander, C.M.O.D., & Raymond, S., Primordial origins of Earth’s carbon. Rev. Mineral. Geochem. 75, 149181 (2013).CrossRefGoogle Scholar
Sarafian, A.R., Nielsen, S.G., Marschall, H.R., McCubbin, F.M., & Monteleone, B.D., Early accretion of water in the inner solar system from a carbonaceous chondrite-like source. Science 346, 623626 (2014).Google Scholar
Grady, M.M., Wright, I.P., Carr, L.P., & Pillinger, C.T., Compositional differences in enstatite chondrites based on carbon and nitrogen stable isotope measurements. Geochim. Cosmochim. Acta 50, 27992813 (1986).CrossRefGoogle Scholar
Alexander, C.M.O.D., Swan, P., & Prombo, C.A., Occurrence and implications of silicon nitride in enstatite chondrites. Meteor. Planet. Sci. 29, 7985 (1994).Google Scholar
Alexander, C.M.O.D., McKeegan, K.D., & Altwegg, K., Water reservoirs in small planetary bodies: meteorites, asteroids, and comets. Space Sci. Rev. 214, 36 (2018).Google Scholar
Sarafian, A.R. et al., Early accretion of water and volatile elements to the inner Solar System: evidence from angrites. Phil. Trans. Royal Soc. A. 375, 20160209 (2017).Google Scholar
Geiss, J. & Gloeckler, G., Isotopic composition of H, HE and NE in the protosolar cloud. Space Sci. Rev. 106, 318 (2003).Google Scholar
Ko, H., Marc, C., Bernard, M., & Kentaro, T., Protosolar carbon isotopic composition: implications for the origin of meteoritic organics. Astrophys. J. 600, 480 (2004).Google Scholar
Marty, B., Chaussidon, M., Wiens, R.C., Jurewicz, A.J.G., & Burnett, D.S., A 15N-poor isotopic composition for the Solar System as shown by genesis solar wind samples. Science 332, 15331536 (2011).Google Scholar
Marty, B. & Yokochi, R., Water in the early Earth. Rev. Mineral. Geochem. 62, 421450 (2006).Google Scholar
Wood, B.J., Li, J., & Shahar, A., Carbon in the core: its influence on the properties of core and mantle. Rev. Mineral. Geochem. 75, 231250 (2013).CrossRefGoogle Scholar
Kerridge, J.F., Carbon, hydrogen and nitrogen in carbonaceous chondrites: abundances and isotopic compositions in bulk samples. Geochim. Cosmochim. Acta 49, 17071714 (1985).CrossRefGoogle ScholarPubMed
Deines, P., The carbon isotope geochemistry of mantle xenoliths. Earth Sci. Rev. 58, 247278 (2002).CrossRefGoogle Scholar
Labidi, J., Cartigny, P., & Moreira, M., Non-chondritic sulphur isotope composition of the terrestrial mantle. Nature 501, 208 (2013).Google Scholar
Satish-Kumar, M., So, H., Yoshino, T., Kato, M., & Hiroi, Y., Experimental determination of carbon isotope fractionation between iron carbide melt and carbon: 12C-enriched carbon in the Earth’s core? Earth Planet. Sci. Lett. 310, 340348 (2011).Google Scholar
Walsh, K.J., Morbidelli, A., Raymond, S.N., O’Brien, D.P., & Mandell, A.M., A low mass for Mars from Jupiter’s early gas-driven migration. Nature 475, 206209 (2011).Google Scholar
Raymond, S.N., O’Brien, D.P., Morbidelli, A., & Kaib, N.A., Building the terrestrial planets: constrained accretion in the inner Solar System. Icarus 203, 644662 (2009).CrossRefGoogle Scholar
Dauphas, N. & Pourmand, A., Hf–W–Th evidence for rapid growth of Mars and its status as a planetary embryo. Nature 473, 489492 (2011).Google Scholar
Lambrechts, M. & Johansen, A., Rapid growth of gas-giant cores by pebble accretion. Astronom. Astrophys. 544, A32 (2012).Google Scholar
Johansen, A., Low, M.-M.M., Lacerda, P., & Bizzarro, M., Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion. Sci. Adv. 1, e1500109 (2015).CrossRefGoogle ScholarPubMed
Wetzel, D.T., Hauri, E.H., Saal, A.E., & Rutherford, M.J., Carbon content and degassing history of the lunar volcanic glasses. Nat. Geosci. 8, 755758 (2015).Google Scholar
Chi, H., Dasgupta, R., Duncan, M., & Shimizu, N., Partitioning of carbon between Fe-rich alloy melt and silicate melt in a magma ocean – implications for the abundance and origin of volatiles in Earth, Mars, and the Moon. Geochim. Cosmochim. Acta 139, 447471 (2014).Google Scholar
Saal, A.E. et al., Volatile content of lunar volcanic glasses and the presence of water in the Moon’s interior. Nature 454, 192195 (2008).CrossRefGoogle ScholarPubMed
Hauri, E.H., Weinreich, T., Saal, A.E., Rutherford, M.C., & Van Orman, J.A., High pre-eruptive water contents preserved in lunar melt inclusions. Science 333, 213215 (2011).Google Scholar
Saal, A.E., Hauri, E.H., Van Orman, J.A., & Rutherford, M.J., Hydrogen isotopes in lunar volcanic glasses and melt inclusions reveal a carbonaceous chondrite heritage. Science 340, 13171320 (2013).Google Scholar
Greenwood, R.C. et al., Oxygen isotopic evidence for accretion of Earth’s water before a high-energy Moon-forming giant impact. Sci. Adv. 4, eaao5928 (2018).Google Scholar
Palme, H., Lodders, K., & Jones, A., Solar system abundances of the elements, in Treatise on Geochemistry (2nd edn.), eds. Holland, H.D. & Turekian, K.K. (Oxford: Elsevier, 2014), pp. 1536.Google Scholar
Anders, E. & Grevesse, N., Abundances of the elements: meteoritic and solar. Geochim. Cosmochim. Acta 53, 197214 (1989).CrossRefGoogle Scholar
Lodders, K., Solar system abundances of the elements, in Principles and Perspectives in Cosmochemistry: Lecture Notes of the Kodai School on “Synthesis of Elements in Stars”, eds. Goswami, A. & Reddy, B. Eswar (Berlin: Springer-Verlag, 2010), pp. 379417.CrossRefGoogle Scholar
Palme, H. & Jones, A., Solar system abundances of the elements, in Meteorites, Comets, and Planets, ed. Davis, A.M. (Amsterdam: Elsevier Ltd., 2003), pp. 4161.Google Scholar
Grady, M.M. & Wright, I.P., Elemental and isotopic abundances of carbon and nitrogen in meteorites. Space Sci. Rev. 106, 231248 (2003).Google Scholar
Hutson, M. & Ruzicka, A., A multi-step model for the origin of E3 (enstatite) chondrites. Meteor. Planet. Sci. 35, 601608 (2000).Google Scholar
Hirschmann, M.M., Water, melting, and the deep Earth H2O cycle. Ann. Rev. Earth Planet. Sci. 34, 629653 (2006).Google Scholar
McDonough, W.F., Compositional model for the Earth’s core, in The Mantle and Core, Vol. 2, ed. Carlson, R.W. (Oxford: Elsevier-Pergamon, 2003), pp. 547568.Google Scholar
Rosenthal, A., Hauri, E.H., & Hirschmann, M.M., Experimental determination of C, F, and H partitioning between mantle minerals and carbonated basalt, CO2/Ba and CO2/Nb systematics of partial melting, and the CO2 contents of basaltic source regions. Earth Planet. Sci. Lett. 412, 7787 (2015).Google Scholar
Dalou, C., Hirschmann, M.M., von der Handt, A., Mosenfelder, J., & Armstrong, L.S., Nitrogen and carbon fractionation during core–mantle differentiation at shallow depth. Earth Planet. Sci. Lett. 458, 141151 (2017).CrossRefGoogle Scholar
Nabiei, F. et al., A large planetary body inferred from diamond inclusions in a ureilite meteorite. Nat. Commun. 9, 1327 (2018).Google Scholar
Kana, N., Masayuki, N., & Jun-ichi, M., Raman spectroscopic study of diamond and graphite in ureilites and the origin of diamonds. Meteor. Planet. Sci. 47, 17281737 (2012).Google Scholar
Barrat, J.-A., Sansjofre, P., Yamaguchi, A., Greenwood, R.C., & Gillet, P., Carbon isotopic variation in ureilites: evidence for an early, volatile-rich inner Solar System. Earth Planet. Sci. Lett. 478, 143149 (2017).Google Scholar
Schiller, M., Bizzarro, M., & Fernandes, V.A., Isotopic evolution of the protoplanetary disk and the building blocks of Earth and the Moon. Nature 555, 507 (2018).Google Scholar
Ding, S. & Dasgupta, R., The fate of sulfide during decompression melting of peridotite – implications for sulfur inventory of the MORB-source depleted upper mantle. Earth Planet. Sci. Lett. 459, 183195 (2017).Google Scholar
Saal, A.E., Hauri, E., Langmuir, C.H., & Perfit, M.R., Vapour undersaturation in primitive mid-ocean-ridge basalt and the volatile content of Earth’s upper mantle. Nature 419, 451455 (2002).Google Scholar
Rehkämper, M. et al., Ir, Ru, Pt, and Pd in basalts and komatiites: new constraints for the geochemical behavior of the platinum-group elements in the mantle. Geochim. Cosmochim. Acta 63, 39153934 (1999).Google Scholar
Walker, R.J., Highly siderophile elements in the Earth, Moon and Mars: update and implications for planetary accretion and differentiation. Geochemisty 69, 101125 (2009).Google Scholar
Mann, U., Frost, D.J., Rubie, D.C., Becker, H., & Audétat, A., Partitioning of Ru, Rh, Pd, Re, Ir and Pt between liquid metal and silicate at high pressures and high temperatures – implications for the origin of highly siderophile element concentrations in the Earth’s mantle. Geochim. Cosmochim. Acta 84, 593613 (2012).Google Scholar
Morgan, J.W., Walker, R.J., Brandon, A.D., & Horan, M.F., Siderophile elements in Earth’s upper mantle and lunar breccias: data synthesis suggests manifestations of the same late influx. Meteor. Planet. Sci. 36, 12571275 (2001).Google Scholar
Willbold, M., Elliott, T., & Moorbath, S., The tungsten isotopic composition of the Earth’s mantle before the terminal bombardment. Nature 477, 195 (2011).Google Scholar
Albarede, F., Volatile accretion history of the terrestrial planets and dynamic implications. Nature 461, 12271233 (2009).Google Scholar
Albarede, F. et al., Asteroidal impacts and the origin of terrestrial and lunar volatiles. Icarus 222, 4452 (2013).Google Scholar
Dasgupta, R., Chi, H., Shimizu, N., Buono, A., & Walker, D., Carbon solution and partitioning between metallic and silicate melts in a shallow magma ocean: implications for the origin and distribution of terrestrial carbon. Geochim. Cosmochim. Acta 102, 191212 (2013).CrossRefGoogle Scholar
Alexander, C.M.O.D., Fogel, M., Yabuta, H., & Cody, G.D., The origin and evolution of chondrites recorded in the elemental and isotopic compositions of their macromolecular organic matter. Geochimi. Cosmochim. Acta 71, 43804403 (2007).Google Scholar
Vdovykin, G.P., Ureilites. Space Sci. Rev. 10, 483510 (1970).Google Scholar
Grady, M.M., Wright, I.P., Swart, P.K., & Pillinger, C.T., The carbon and nitrogen isotopic composition of ureilites: implications for their genesis. Geochim. Cosmochim. Acta 49, 903915 (1985).Google Scholar
Downes, H. et al., Isotopic composition of carbon and nitrogen in ureilitic fragments of the Almahata Sitta meteorite. Meteor. Planet. Sci. 50, 255272 (2015).Google Scholar
Siebert, J., Badro, J., Antonangeli, D., & Ryerson, F.J., Metal–silicate partitioning of Ni and Co in a deep magma ocean. Earth Planet. Sci. Lett. 321 –322, 189197 (2012).Google Scholar
Grewal, D.S., Dasgupta, R., Holmes, A.K., Costin, G., & Li, Y., The fate of nitrogen during core-mantle separation on Earth. Geochim. Cosmochim. Acta 251, 87115 (2019) .Google Scholar
Boujibar, A. et al., Metal–silicate partitioning of sulphur, new experimental and thermodynamic constraints on planetary accretion. Earth Planet. Sci. Lett. 391, 4254 (2014).Google Scholar
Clesi, V. et al., Low hydrogen contents in the cores of terrestrial planets. Science Adv. 4, e1701876 (2018).Google Scholar
Armstrong, L.S., Hirschmann, M.M., Stanley, B.D., Falksen, E.G., & Jacobsen, S.D., Speciation and solubility of reduced C–O–H–N volatiles in mafic melt: implications for volcanism, atmospheric evolution, and deep volatile cycles in the terrestrial planets. Geochim. Cosmochim. Acta 171, 283302 (2015).Google Scholar
Roskosz, M., Bouhifd, M.A., Jephcoat, A.P., Marty, B., & Mysen, B.O., Nitrogen solubility in molten metal and silicate at high pressure and temperature. Geochim. Cosmochim. Acta 121, 1528 (2013).Google Scholar
O’Neill, H.S.C. & Mavrogenes, J.A., The sulfide capacity and the sulfur content at sulfide saturation of silicate melts at 1400°C and 1 bar. J. Petrol. 43, 10491087 (2002).Google Scholar
Hirschmann, M.M., Withers, A.C., Ardia, P., & Foley, N.T., Solubility of molecular hydrogen in silicate melts and consequences for volatile evolution of terrestrial planets. Earth Planet. Sci. Lett. 345, 3848 (2012).Google Scholar
Rubie, D.C. et al., Highly siderophile elements were stripped from Earth’s mantle by iron sulfide segregation. Science 353, 11411144 (2016).Google Scholar
Corgne, A., Wood, B.J., & Fei, Y., C- and S-rich molten alloy immiscibility and core formation of planetesimals. Geochim. Cosmochim. Acta 72, 24092416 (2008).Google Scholar
Dasgupta, R., Buono, A., Whelan, G., & Walker, D., High-pressure melting relations in Fe–C–S systems: implications for formation, evolution, and structure of metallic cores in planetary bodies. Geochim. Cosmochim. Acta 73, 66786691 (2009).CrossRefGoogle Scholar
Tsuno, K. & Dasgupta, R., Fe–Ni–Cu–C–S phase relations at high pressures and temperatures – the role of sulfur in carbon storage and diamond stability at mid- to deep-upper mantle. Earth Planet. Sci. Lett. 412, 132142 (2015).Google Scholar
Zhang, Z., Hastings, P., Von der Handt, A., & Hirschmann, M.M., Experimental determination of carbon solubility in Fe–Ni–S melts. Geochim. Cosmochim. Acta 225, 6679 (2018).CrossRefGoogle Scholar
Li, Y., Dasgupta, R., & Tsuno, K., The effects of sulfur, silicon, water, and oxygen fugacity on carbon solubility and partitioning in Fe-rich alloy and silicate melt systems at 3 GPa and 1600 °C: implications for core–mantle differentiation and degassing of magma oceans and reduced planetary mantles. Earth Planet. Sci. Lett. 415, 5466 (2015).Google Scholar
Raymond, S.N., Quinn, T., & Lunine, J.I., High-resolution simulations of the final assembly of Earth-like planets. 2. Water delivery and planetary habitability. Astrobiology 7, 6684 (2007).Google Scholar
Hallis, L.J. et al., Evidence for primordial water in Earth’s deep mantle. Science 350, 795797 (2015).Google Scholar
Tucker, J.M. & Mukhopadhyay, S., Evidence for multiple magma ocean outgassing and atmospheric loss episodes from mantle noble gases. Earth Planet. Sci. Lett. 393, 254265 (2014).Google Scholar
Kress, V.C. & Carmichael, I.S.E., The compressibility of silicate liquids containing Fe2O3 and the effect of composition, temperature, oxygen fugacity and pressure on their redox states. Contrib. Mineral. Petrol. 108, 8292 (1991).Google Scholar
O’Neill, H.S.C. et al., An experimental determination of the effect of pressure on the Fe3+/∑Fe ratio of an anhydrous silicate melt to 3.0 GPa. Am. Mineral. 91, 404412 (2006).Google Scholar
Zhang, H.L., Hirschmann, M.M., Cottrell, E., & Withers, A.C., Effect of pressure on Fe3+/ΣFe ratio in a mafic magma and consequences for magma ocean redox gradients. Geochim. Cosmochim. Acta 204, 83103 (2017).Google Scholar
Schaefer, L. & Elkins-Tanton, L.T., Magma oceans as a critical stage in the tectonic development of rocky planets. Phil. Trans. Royal Soc. A. 376, 20180109 (2018).Google Scholar
Duncan, M.S., Dasgupta, R., & Tsuno, K., Experimental determination of CO2 content at graphite saturation along a natural basalt–peridotite melt join: implications for the fate of carbon in terrestrial magma oceans. Earth Planet. Sci. Lett. 466, 115128 (2017).Google Scholar
Libourel, G., Marty, B., & Humbert, F., Nitrogen solubility in basaltic melt. Part I. Effect of oxygen fugacity. Geochim. Cosmochim. Acta 67, 41234135 (2003).Google Scholar
Kadik, A.A. et al., Influence of oxygen fugacity on the solubility of nitrogen, carbon, and hydrogen in FeO–Na2O–SiO2–Al2O3 melts in equilibrium with metallic iron at 1.5 GPa and 1400°C. Geochem. Int. 49, 429438 (2011).Google Scholar
Holloway, J.R., Pan, V., & Gudmundsson, G., High-pressure fluid-absent melting experiments in the presence of graphite; oxygen fugacity, ferric/ferrous ratio and dissolved CO2. Eur. J. Mineral. 4, 105114 (1992).Google Scholar
Li, Y., Dasgupta, R., & Tsuno, K., Carbon contents in reduced basalts at graphite saturation: implications for the degassing of Mars, Mercury, and the Moon. J. Geophys. Res. Planets 122, doi:10.1002/2017JE005289 (2017).Google Scholar
Ding, S., Hough, T., & Dasgupta, R., New high pressure experiments on sulfide saturation of high-FeO* basalts with variable TiO2 contents – implications for the sulfur inventory of the lunar interior. Geochim. Cosmochim. Acta 222, 319339 (2018).CrossRefGoogle Scholar
Baker, D.R. & Moretti, R., Modeling the solubility of sulfur in magmas: a 50-year old geochemical challenge. Rev. Mineral. Geochem. 73, 167213 (2011).Google Scholar
Genda, H. & Abe, Y., Enhanced atmospheric loss on protoplanets at the giant impact phase in the presence of oceans. Nature 433, 842844 (2005).Google Scholar
Honda, M., McDougall, I., Patterson, D.B., Doulgeris, A., & Clague, D.A., Possible solar noble-gas component in Hawaiian basalts. Nature 349, 149 (1991).Google Scholar
Sarda, P., Staudacher, T., & Allègre, C.J., Neon isotopes in submarine basalts. Earth Planet. Sci. Lett. 91, 7388 (1988).Google Scholar
Mizuno, H., Nakazawa, K., & Hayashi, C., Dissolution of the primordial rare gases into the molten Earth’s material. Earth Planet. Sci. Lett. 50, 202210 (1980).Google Scholar
Yokochi, R. & Marty, B., A determination of the neon isotopic composition of the deep mantle. Earth Planet. Sci. Lett. 225, 7788 (2004).Google Scholar
Lynne, A.H., Disk-dispersal and planet-formation timescales. Phys. Scipta 2008, 014024 (2008).Google Scholar
Yin, Q. et al., A short timescale for terrestrial planet formation from Hf–W chronometry of meteorites. Nature 418, 949 (2002).Google Scholar
Tang, H. & Dauphas, N., 60Fe–60Ni chronology of core formation in Mars. Earth Planet. Sci. Lett. 390, 264274 (2014).Google Scholar
Rose-Weston, L., Brenan, J.M., Fei, Y., Secco, R.A., & Frost, D.J., Effect of pressure, temperature, and oxygen fugacity on the metal-silicate partitioning of Te, Se, and S: implications for earth differentiation. Geochim. Cosmochim. Acta 73, 45984615 (2009).Google Scholar
Zhang, Y. & Yin, Q.-Z., Carbon and other light element contents in the Earth’s core based on first-principles molecular dynamics. Proc. Natl Acad. Sci. 109, 1957919583 (2012).Google Scholar
Suer, T.-A., Siebert, J., Remusat, L., Menguy, N., & Fiquet, G., A sulfur-poor terrestrial core inferred from metal–silicate partitioning experiments. Earth Planet. Sci. Lett. 469, 8497 (2017).Google Scholar
Tsymbulov, L.B. & Tsmekhman, L.S., Solubility of carbon in sulfide melts of the system Fe–Ni-S. Russ. J. Appl. Chem. 74, 925929 (2001).Google Scholar
Stolper, E., Fine, G., Johnson, T., & Newman, S., Solubility of carbon dioxide in albitic melt. Am. Mineral. 72, 10711085 (1987).Google Scholar
Duncan, M.S. & Dasgupta, R., Pressure and temperature dependence of CO2 solubility in hydrous rhyolitic melt – implications for carbon transfer to mantle source of volcanic arcs via partial melt of subducting crustal lithologies. Contrib. Mineral. Petrol. 169, 119 (2015).Google Scholar
Eguchi, J. & Dasgupta, R., A CO2 solubility model for silicate melts from fluid saturation to graphite or diamond saturation. Chem. Geol. 487, 2338 (2018).Google Scholar
Behrens, H., Ohlhorst, S., Holtz, F., & Champenois, M., CO2 solubility in dacitic melts equilibrated with H2O–CO2 fluids: implications for modeling the solubility of CO2 in silicic melts. Geochim. Cosmochim. Acta 68, 46874703 (2004).Google Scholar
Duncan, M.S. & Dasgupta, R., Rise of Earth’s atmospheric oxygen controlled by efficient subduction of organic carbon. Nat. Geosci. 10, 387392 (2017).Google Scholar
Stolper, E. & Holloway, J.R., Experimental determination of the solubility of carbon dioxide in molten basalt at low pressure. Earth Planet. Sci. Lett. 87, 397408 (1988).Google Scholar
Pan, V., Holloway, J.R., & Hervig, R.L., The pressure and temperature dependence of carbon dioxide solubility in tholeiitic basalt melts. Geochim. Cosmochim. Acta 55, 15871595 (1991).CrossRefGoogle Scholar
Yoshioka, T., McCammon, C.A., Shcheka, S., & Keppler, H., The speciation of carbon monoxide in silicate melts and glasses. Am. Mineral. 100, 16411644 (2015).Google Scholar
Wetzel, D.T., Rutherford, M.J., Jacobsen, S.D., Hauri, E.H., & Saal, A.E., Degassing of reduced carbon from planetary basalts. Proc. Natl Acad. Sci. 110, 80108013 (2013).Google Scholar
Dasgupta, R., Hirschmann, M.M., & Smith, N.D., Partial melting experiments of peridotite + CO2 at 3 GPa and genesis of alkalic ocean island basalts. J. Petrol. 48, 20932124 (2007).Google Scholar
Mallik, A. & Dasgupta, R., Effect of variable CO2 on eclogite-derived andesite-lherzolite reaction at 3 GPa – implications for mantle source characteristics of alkalic ocean island basalts. Geochem. Geophys. Geosyst. 15, 15331557 (2014).Google Scholar
Vetere, F., Holtz, F., Behrens, H., Botcharnikov, R.E., & Fanara, S., The effect of alkalis and polymerization on the solubility of H2O and CO2 in alkali-rich silicate melts. Contrib. Mineral. Petrol. 167, 1014 (2014).Google Scholar
Ghiorso, M.S. & Gualda, G.A.R., An H2O–CO2 mixed fluid saturation model compatible with rhyolite-MELTS. Contrib. Mineral. Petrol. 169, 53 (2015).Google Scholar
Stanley, B.D., Hirschmann, M.M., & Withers, A.C., Solubility of C–O–H volatiles in graphite-saturated martian basalts. Geochim. Cosmochim. Acta 129, 5476 (2014).Google Scholar
Wade, J. & Wood, B.J., Core formation and the oxidation state of the Earth. Earth Planet. Sci. Lett. 236, 7895 (2005).Google Scholar
Rubie, D.C. et al., Heterogeneous accretion, composition and core–mantle differentiation of the Earth. Earth Planet. Sci. Lett. 301, 3142 (2011).Google Scholar
Badro, J., Brodholt, J.P., Piet, H., Siebert, J., & Ryerson, F.J., Core formation and core composition from coupled geochemical and geophysical constraints. Proc. Natl Acad. Sci. 112, 1231012314 (2015).Google Scholar
Okuchi, T., Hydrogen partitioning into molten iron at high pressure: implications for Earth’s core. Science 278, 17811784 (1997).Google Scholar
O’Neill, H.S.C., The origin of the moon and the early history of the earth – a chemical model. Part 2: the earth. Geochim. Cosmochim. Acta 55, 11591172 (1991).Google Scholar
Kleine, T. et al., Hf–W chronology of the accretion and early evolution of asteroids and terrestrial planets. Geochim. Cosmochim. Acta 73, 51505188 (2009).Google Scholar
Rubie, D.C. et al., Accretion and differentiation of the terrestrial planets with implications for the compositions of early-formed Solar System bodies and accretion of water. Icarus 248, 89108 (2015).Google Scholar
Canup, R.M. & Asphaug, E., Origin of the Moon in a giant impact near the end of the Earth’s formation. Nature 412, 708 (2001).Google Scholar
Ćuk, M. & Stewart, S.T., Making the Moon from a fast-spinning Earth: a giant impact followed by resonant despinning. Science 338, 10471052 (2012).Google Scholar
Deguen, R., Olson, P., & Cardin, P., Experiments on turbulent metal–silicate mixing in a magma ocean. Earth Planet. Sci. Lett. 310, 303313 (2011).Google Scholar
Deguen, R., Landeau, M., & Olson, P., Turbulent metal–silicate mixing, fragmentation, and equilibration in magma oceans. Earth Planet. Sci. Lett. 391, 274287 (2014).Google Scholar
Wohlers, A. & Wood, B.J., A Mercury-like component of early Earth yields uranium in the core and high mantle 142Nd. Nature 520, 337340 (2015).Google Scholar
Wohlers, A. & Wood, B.J., Uranium, thorium and REE partitioning into sulfide liquids: implications for reduced S-rich bodies. Geochim. Cosmochim. Acta 205, 226244 (2017).Google Scholar
Münker, C., Fonseca, R.O.C., & Schulz, T., Silicate Earth’s missing niobium may have been sequestered into asteroidal cores. Nat. Geosci. 10, 822 (2017).Google Scholar
Wade, J. & Wood, B.J., The oxidation state and mass of the Moon-forming impactor. Earth Planet. Sci. Lett. 442, 186193 (2016).Google Scholar
Krasnokutski, S.A. et al., Low-temperature condensation of carbon. Astrophys. J. 847, 89 (2017).Google Scholar
Li, Y., Marty, B., Shcheka, S., Zimmermann, L., & Keppler, H., Nitrogen isotope fractionation during terrestrial core–mantle separation. Geochem. Perspect. Lett. 2, 138147 (2016).Google Scholar
Labidi, J. et al., Experimentally determined sulfur isotope fractionation between metal and silicate and implications for planetary differentiation. Geochim. Cosmochim. Acta 175, 181194 (2016).Google Scholar
Figure 0

Figure 2.1 Plausible initial distributions of C, H, N, and S between Earth’s major reservoirs at the end of accretion and core formation. (a) LEVEs are available during accretion and all LEVEs are initially sequestered in the core and the MO. (b) LEVE delivery during accretion is temporally separate from the main episodes of core formation and thus the core is effectively devoid of LEVEs. (c) Available LEVEs during accretion are effectively sequestered in the core or lost to space (not shown) and the BSE acquires LEVEs via later additions. (d) Core–mantle differentiation takes place in the absence of LEVEs and LEVEs are added via late additions. These distributions eventually evolved into the setup seen in the modern world (e), where these volatile elements are thought to exchange between Earth’s mantle and the crust–ocean–atmosphere system (ingassing and outgassing) regulates the long-term habitability of the planet.

Figure 1

Figure 2.2 Comparison of isotopic compositions of (a) hydrogen and nitrogen and (b) carbon and nitrogen for several Solar System objects and reservoirs. The solar reservoir is depleted while the cometary reservoir is enriched in the heavier isotopes of nitrogen and hydrogen relative to all classes of meteorites as well as the BSE. Although δ13C alone cannot distinguish between the meteoritic and cometary sources, the isotopic compositions of all major LEVEs suggest that the BSE had a similar parent reservoir as carbonaceous chondrites.

Data sources: Carbonaceous chondrites, Refs. 48, 49; E-chondrites, Refs. 52, 53; comets, Ref. 54; angrites and Vesta, Refs. 51, 55; and solar, Refs. 56–58.
Figure 2

Table 2.1 C content and C/S, C/N, and C/H weight ratios of major terrestrial reservoirs and chondritic building blocks

Figure 3

Figure 2.3 Comparison of the abundances of C, N, S, and H in the BSE normalized to their respective abundances in CI chondrites. In the BSE, C is enriched relative to N, has a similar abundance to S, and is depleted relative to H in comparison to its abundance in CI chondrites. Data from Refs. 10, 11, and 13. Also noted in the figure is the main geochemical information used to estimate the BSE concentrations of the respective elements in the original studies. MORB = mid-ocean ridge basalt.

Figure 4

Figure 2.4 Comparison of bulk C/N (a), C/H (b), and C/S (c) weight ratios between different terrestrial and chondritic reservoirs. The data are from Table 2.1. CI = Ivuna-type Carbonaceous chondrite; CM = Mighei-type Carbonaceous chondrite; CV = Vigarano-type Carbonaceous chondrite; CO = Ornans-type Carbonaceous chondrite; LL = low total iron, low metal ordinary chondrite; H = highest total iron ordinary chondrite; L = low total iron ordinary chondrite; EL = (low enstatite) enstatite chondrite; EH = (high enstatite) enstatite chondrite.

Figure 5

Figure 2.5 Application of alloy–silicate partition coefficients and solubilities of LEVEs in the silicate melts to examine the effect of core formation, with varying degrees of alloy–silicate equilibration, with or without loss of an early atmosphere formed via MO degassing on the remnant abundances of LEVEs in the bulk silicate reservoir. (a) LEVEs, when delivered as 0.015 ME late-accreting materials (i.e. 0% alloy–silicate equilibration), cause the volatile abundance to be higher than the present-day BSE. Core formation with increasing degrees of alloy–silicate equilibration increasingly depletes the remnant MO in all LEVEs, with C being much more depleted than other LEVEs, leading to subchondritic C/N, C/H, and C/S ratios. (b) Combining early atmospheric loss with core formation cannot offset C loss to the core due to the lower solubility of C relative to the other LEVEs in the silicate MOs. Bulk Earth volatile abundance data are from McDonough,83 while the alloy–silicate partition coefficients in a deep MO (P = 50 GPa, T = 3500 K; e.g. Siebert et al.104) for C, N, S, and H are from the parametrized relationships of Chi et al.,71 Grewal et al.,105 Boujibar et al.,106 and Clesi et al.,107 respectively. Solubility constant data for C, N, S, and H in the silicate melt are from Armstrong et al.,108 Roskosz et al.,109 O’Neill and Mavrogenes,110 and Hirschmann et al.111

Figure 6

Figure 2.6 Experimental data showing the effects of S and Si contents in the Fe-rich alloy melt on C solubility. (a) Carbon solubility decreases monotonically with increases in sulfur content in the alloy melt. Si- and N-bearing alloys have lower C solubility for a given S content in the alloy melt. (b) Carbon solubility decreases linearly with increases in Si content in the alloy or decreases in oxygen fugacity (fO2) of silicate–alloy systems. The presence of S in the alloy does not have a major effect on carbon solubility in silicon-bearing alloys.

Data for N-free alloys are from Dasgupta et al.,99,114 Chi et al.,71 Armstrong et al.,108 and Li et al.,16,117 while data for N-bearing alloys are from Dalou et al.85 and Grewal et al.18
Figure 7

Figure 2.7 Experimental data showing the effects of fO2 and water content in the silicate melt on C solubility. Carbon solubility decreases with decreasing log fO2 until ~ΔIW of –2 to –1 followed by an increase below ~IW – 2 depending upon bulk water content in the silicate melt. This change is likely caused by the change in contributions of C as CO32−CO32− at higher fO2 and C as C–H species such as methyl groups at lower fO2. The pressure–temperature–composition (expressed as NBO/T) condition of the experimental glasses are given in the legend. The shaded region represents predicted solubility of C as CO32−CO32−at graphite saturation based on the model of Eguchi and Dasgupta146 for the entire pressure, temperature, and XH2O space of the experimental data, while the NBO/T parameter is extrapolated for peridotite-like silicate melts (NBO/T = 2.7).

Data for N-free silicate melts are from Refs. 16, 17, 71, 99, 108, 117, 125, 129, 152, and 157.
Figure 8

Figure 2.8 Experimental data showing the effects of S and Si in Fe-rich alloy melts on DCalloy/silicate. (a) DCalloy/silicate decreases with increases in S content in the alloy melt. (b) DCalloy/silicate decreases with increases in Si content in the alloy or indirect decreases in fO2 of the alloy–silicate systems. The presence of S in Si-bearing alloy melts generally decreases DCalloy/silicate in comparison to S-free systems.

Data for N-free alloys are from Dasgupta et al.,99,114 Chi et al.,71 Armstrong et al.,108 Li et al.,16,117 and Tsuno et al.,17 while data for N-bearing alloys are from Dalou et al.85 and Grewal et al.18
Figure 9

Figure 2.9 The effects of different accretion scenarios on the core–mantle partitioning behavior of carbon. (a) Four possible evolution paths for the FeO content of the mantle with its growth. The “very reduced” and “reduced” paths simulate the accretion of progressively oxidized material with time, and the “oxidized” path simulates accretion of progressively reduced material with time. All of the models converge at log fO2 of IW – 1.8, or present-day FeO content (8 wt.%) of the primitive upper mantle. (b) DCalloy/silicateat each step of accretion (i.e. accretion of 1% ME), calculated using parametrized relationships from Li et al.16 For the “very reduced” and “reduced” paths, where methyl (or alkyl) groups would be the predominant dissolved carbon species in the silicate melt, eq. (4) is chosen from Li et al.,16 whereas for the “intermediate” and “oxidized” paths, where carbonate groups would be the predominant dissolved carbon species in the silicate melt, eq. (3) is chosen from Li et al.16 (c) Post-core formation MO resulting from a bulk system with 1000 ppm C and assuming 100% alloy–silicate equilibration at each accretional step. All of the models lead to MO with distinctly less C that those estimated for the BSE. (d) Carbon content in Earth’s core converges at ~3300 ppm at the end of all accretional models.

Figure 10

Figure 2.10 A schematic depiction of the late-stage accretion for a large rocky planet such as Earth where a large fraction of mass increase takes place via giant impacts of small planets and planetary embryos (0.1–0.2 ME).1618 This mode of accretion creates a unique possibility of establishing the LEVE inventory in the degassable portion of the rocky planets.

Save book to Kindle

To save this book to your Kindle, first ensure [email protected] is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about saving to your Kindle.

Note you can select to save to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

Available formats
×

Save book to Dropbox

To save content items to your account, please confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account. Find out more about saving content to Dropbox.

Available formats
×

Save book to Google Drive

To save content items to your account, please confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account. Find out more about saving content to Google Drive.

Available formats
×