Research ArticleSPACE SCIENCES

Origin of uranium isotope variations in early solar nebula condensates

See allHide authors and affiliations

Science Advances  04 Mar 2016:
Vol. 2, no. 3, e1501400
DOI: 10.1126/sciadv.1501400


High-temperature condensates found in meteorites display uranium isotopic variations (235U/238U), which complicate dating the solar system’s formation and whose origin remains mysterious. It is possible that these variations are due to the decay of the short-lived radionuclide 247Cm (t1/2 = 15.6 My) into 235U, but they could also be due to uranium kinetic isotopic fractionation during condensation. We report uranium isotope measurements of meteoritic refractory inclusions that reveal excesses of 235U reaching ~+6% relative to average solar system composition, which can only be due to the decay of 247Cm. This allows us to constrain the 247Cm/235U ratio at solar system formation to (7.0 ± 1.6) × 10−5. This value provides new clues on the universality of the nucleosynthetic r-process of rapid neutron capture.

  • nucleosynthesis
  • meteorites
  • uranium isotopes
  • curium-247
  • r-process
  • fine-grained CAI
  • group II REE pattern


All elements beyond the iron peak (above ~70 atomic mass units) are the products of three main processes of stellar nucleosynthesis: the s- (slow neutron capture), r- (rapid neutron capture), and p-process (proton process) (1, 2). Unlike the s- and p-process, which are relatively well understood [neutron capture in asymptotic giant branch (AGB) stars for the s-process and photodisintegration of seed nuclei in supernovae for the p-process] (3, 4), little is known regarding the astrophysical conditions under which r-process nuclides are produced (58). The r-process label may comprise some diversity as it was suggested that up to three r-processes were responsible for producing light r-nuclides (A < 140), heavy r-nuclides (A > 140), and actinides (for example, U and Th). The existence of an “actinide” production site is motivated by two main observations: (i) the ages of some old stars inferred from their Th/Eu ratios are negative, meaning that they must have formed with a Th/Eu ratio that was higher than that relevant to the solar system (SS) (9), and (ii) the abundance of the short-lived radionuclide (SLR) 244Pu [t1/2 = 79.3 million years (My)] in meteorites is low compared to another nominal r-process radionuclide, 182Hf (7). However, 244Pu has a long half-life and its stellar yield is uncertain (10), which makes it insensitive to the history of nucleosynthesis before SS formation and whether or not multiple r-process sites contributed to the synthesis of the actinides. In contrast, 247Cm, which decays into 235U, has a much shorter half-life (t1/2 = 15.6 My) and would be very well suited to address this question. Unfortunately, its abundance in the early solar system (ESS) has been the subject of some debate that boils down to knowing whether the uranium isotope variations (that is, 235U/238U ratio) measured in early-formed nebular condensates [calcium- and aluminum-rich inclusions (CAIs)] are due to the decay of 247Cm or isotopic fractionation during condensation. The two cannot be easily distinguished because U has only two long-lived isotopes and both mechanisms would produce similar correlations between light U isotope enrichments and U concentrations.

Recently, excesses in 235U of up to 3.5‰ were documented in four fine-grained CAIs (11). To demonstrate that those variations are due to 247Cm decay, one must show that the δ235U values correlate with the Cm/U parent-to-daughter ratios. Because Cm has no long-lived isotopes, another element must be used as a proxy for Cm in isochron diagrams. On the basis of their near-identical valence states, ionic radii, and volatilities, the light rare earth elements (REEs) are thought to behave similarly to Pu and Cm during nebular processes (1214). This conclusion is supported by the observed coherent behavior of Pu, Nd, and Sm in pyroxene/melt and phosphate/melt partitioning experiments (15) and during magmatic differentiation in achondrites (12, 14). Although Th has sometimes been used as a proxy for Cm (11, 16), coherent behavior for the Th-Cm and Th-Pu pairs is neither expected (13) nor observed (12), and the light REEs (for example, Nd or Sm) are therefore taken as more reliable proxies of Cm behavior during condensation.

The 235U excesses observed by Brennecka et al. (11) correlate with Nd/U ratios, which the authors interpreted as evidence of live 247Cm in the ESS at a level of (247Cm/235U)ESS = (1.1 to 2.4) × 10−4. However, such 235U enrichments could also reflect mass-dependent fractionation during condensation of solid CAIs from nebular gas (1719). Indeed, the kinetic theory of gases predicts that the lighter isotope (235U) should condense faster than the heavier isotope (238U), resulting in fractionations that can reach Embedded Image for condensation of atomic U in a low-pressure gas. In both cases of 247Cm decay and fractionation during condensation, one would expect to find correlations between U isotopic variations and the degree of U depletion. It is therefore presently undecided which mechanism is responsible for the U isotope variations documented in CAIs.

Definitive evidence of live 247Cm therefore awaits the discovery of meteoritic material that is highly depleted in uranium and displays 235U excess outside of the ±6‰ window allowed by fractionation at condensation. Because the abundance of 247Cm in the ESS is expected to be low (10), such large excesses of 235U from 247Cm decay will only be resolvable in phases with 144Nd/238U atomic ratios exceeding ~2100 (atomic ratios will be used hereafter). Therefore, numerous slabs of the Allende meteorite were examined, and 15 CAIs (12 fine-grained and 3 coarse-grained) were extracted and digested with acids. Many of these CAIs revealed group II REE patterns and large U depletions indicative of incomplete condensation of refractory lithophile elements (see the Supplementary Materials). A method was developed to measure accurately the U isotopic compositions of samples with low U contents (0.1 to 0.01 ng; see the Supplementary Materials).


As in previous studies, most samples display low 144Nd/238U ratios (that is, <900), and their δ235U values are within 6‰ of the bulk SS value at δ235U = +0.31‰ (relative to CRM-112a) (19, 20). These samples display a trend between δ235U and 144Nd/238U similar to that described previously on similar samples (11) and which had been taken as evidence that 247Cm was alive in the ESS. The samples display significant scatter around the best-fit line that cannot be entirely explained by analytical uncertainty [the mean square weighted deviation (MSWD) of the regression for samples with 144Nd/238U ratios lower than 900 is 46]. Even coarse-grained CAIs that are nondepleted in U (low 144Nd/238U ratios) have δ235U that vary between −0.16 and +1.35‰, well outside of analytical uncertainty.

One sample (named Curious Marie) has an extremely high 144Nd/238U ratio (~22,640) and a 235U excess of +58.9 ± 1.9‰ (equivalent to an absolute 238U/235U ratio of 130.17 ± 0.27; Fig. 1 and Table 1). For comparison, the highest 144Nd/238U ratio and 235U excess measured in CAIs prior to this work were 794 and +3.43‰ (238U/235U = 137.37), respectively (11). Considerable effort was expended to confirm this result (see the Supplementary Materials). The measurement was triplicated using different sample purification schemes and various measurement setups. The tests yielded δ235U values of +52.79 ± 14.91‰ after one purification step with 235U measured on Faraday, +59.12 ± 2.80‰ after two purifications with 235U on an electron multiplier, and +58.97 ± 2.72‰ after three purifications with 235U on an electron multiplier. Extensive testing was also done to ensure that no interferences or matrix effects were affecting the measurement by combining the matrix cut from Curious Marie with CRM-112a, purifying the U by column chemistry, and finding that the measured U isotopic composition is correct (that is, identical to pure CRM-112a) after purification.

Fig. 1 δ235U plotted as a function of the 144Nd/238U atomic ratio in meteoritic samples.

Open circles, previous studies (11, 1719, 33, 34); blue circles, Allende CAIs from this work; light-blue square, bulk Allende from this work]. The +59‰ δ235U value observed in the Curious Marie CAI is well outside the range of variations expected from fractionation during condensation (gray rectangle) and is thus interpreted as definitive evidence for live 247Cm in the ESS. The scatter in the data (for example, at very low Nd/U ratios) suggests that stable isotopic fractionation during evaporation/condensation also influenced the U isotopic composition of CAIs. The slope of the two-point isochron between Curious Marie and the rest of the samples translates into a 247Cm/235U of (5.6 ± 0.3) × 10−5 at the time of Nd/U fractionation in Curious Marie. Accounting for a possible delay between this fractionation event (possibly related to the extensive alteration of this CAI) and the formation of the SS of 5 ± 5 My, the inferred 247Cm/235U at SS formation is (7.0 ± 1.6) × 10−5 (red line).

Table 1 Type, REE pattern, mass, U content, 144Nd/238U atomic ratio, and U isotopic composition of the samples analyzed in this study.

n is the number of replicate analyses for each sample (starting from digested sample). Nd/U ratios and isotopic compositions are corrected for blank contributions (see table S2). δ235U = [(235U/238U)sample/(235U/238U)CRM-112a − 1] × 103. Error bars are 95% confidence intervals. Nd/U ratio was calculated using the U concentration from the double-spike technique, and Nd concentration from the standard addition technique (see the Supplementary Materials).

View this table:


Evidence for 247Cm

The finding of such a large excess 235U in a normal CAI [in opposition to fractionated and unknown nuclear (FUN) CAIs, which is a group of refractory inclusions that display FUN effects; see the Supplementary Materials] can only be explained by decay of 247Cm into 235U. Below, we examine other processes that could potentially lead to U isotope variations in ESS materials and show that they all suffer from serious shortcomings:

(i) Isotopic fractionation during secondary processes (for example, aqueous alteration and/or redox processes), either on the meteorite parent body or on Earth, can impart large isotopic fractionation to light elements (such as Li). However, the degree of fractionation generally decreases as the mass of the element increases, and for U on Earth, the variations are limited to ~1.5 to 2‰ (21). Although this process could lead to some scatter in the data around the isochron, it cannot explain a +59‰ anomaly.

(ii) Large isotopic variations of nucleosynthetic origin have been documented for refractory elements in CAIs and are usually readily identified as departures from mass-dependent fractionation. For uranium, which has only two stable isotopes, it is impossible to distinguish between 247Cm decay and nucleosynthetic anomalies. Nevertheless, considering that anomalies on the order of a few permil in heavy elements (for example, Ba, Sm, and Nd) are only found in FUN CAIs (normal CAIs have more subdued anomalies of a few tenths of permil) and given that Curious Marie (δ238U ~ 59‰) is not a FUN CAI (see the Supplementary Materials), the large 235U excess found here cannot be reasonably ascribed to the presence of nucleosynthetic anomalies.

(iii) Finally, U isotopic fractionation could be the result of fractionation during evaporation/condensation processes. During condensation, the light isotope of U condenses faster, leading to large 235U depletion in the condensing gas and the instantaneous solid, but 235U excesses limited to +6‰ (see the Supplementary Materials). Similarly, during evaporation, the highest 235U excess predicted by the kinetic theory of gases is limited to ~+6‰. Consequently, evaporation/condensation processes cannot explain the +59‰ 235U excess observed in Curious Marie, which leaves 247Cm decay as the only possible explanation.

Closure time in Curious Marie and 247Cm/235U ratio in the ESS

The large excess of 235U found in the Curious Marie CAI is definitive evidence that 247Cm was alive in the ESS. An initial (247Cm/235U) ratio of ~5.6 × 10−5 at the time of closure can be calculated using the slope of the isochron in Fig. 1. Even though the samples with low 144Nd/238U ratios define a trend with δ235U, the slope of the isochron in Fig. 1 is mainly leveraged by Curious Marie (144Nd/238U ratio ~ 22,640). The initial 247Cm/235U ratio of (5.6 ± 0.3) × 10−5 thus corresponds to the time when this CAI acquired its high 144Nd/238U ratio. Terrestrial alteration is ruled out because the Allende meteorite is an observed fall that did not experience much terrestrial weathering. The extreme uranium depletion in the Curious Marie CAI is thus most likely due to solar nebula condensation and/or nebular/parent body alteration. All fine-grained CAIs in this study display a typical group II REE pattern, thought to represent a snapshot in time and space of the condensation sequence (22, 23). The most refractory REEs (heavy REEs except Tm and Yb) are depleted in these CAIs because they were sequestered in ultrarefractory dust such as perovskite or hibonite, whereas the more volatile REEs (Eu and Yb) and uranium are also depleted because they stayed in the gas phase when those CAIs formed. In most samples, U and Yb present similar levels of depletion relative to solar composition and the abundance of other refractory lithophile elements (Fig. 2), indicating that those two elements have similar behaviors during evaporation/condensation processes under solar nebula conditions. However, in Curious Marie, the U/Nd ratio is 1000 times lower than solar composition, whereas the Yb/Nd ratio is only depleted by a factor of 50. If U and Yb have similar behaviors during condensation, one would expect the U/Nd ratio to be 50 times lower than solar, not 1000 times as is observed. Compared to other fine-grained CAIs analyzed, Curious Marie is peculiar because it is extremely altered, which is manifested by the extensive replacement of high-T phases by low-T alteration products such as nepheline and sodalite [its Na2O is 15 weight percent (wt %) when other CAIs are all lower than 6.7 wt %; table S1]. Such alteration may have mobilized U, producing a 20-fold U depletion on top of the 50-fold depletion associated with condensation.

Fig. 2 Correlation between U and Yb abundances relative to solar composition and the abundance of Nd (a refractory lithophile element).

Circles denote data from different studies (3537) and this study, and triangles denote data from Brennecka et al. (11). This correlation indicates that U and Yb have similar behaviors during evaporation/condensation processes under solar nebula conditions. The Curious Marie CAI plots off the 1:1 line, with a CI-normalized U/Nd ratio 20 times lower than that of Yb/Nd. This CAI is extremely altered (see main text and table S1), and such alteration may have mobilized U, producing a 20-fold U depletion on top of the 50-fold depletion associated with condensation.

Some of these alteration products could have formed (i) in the nebula, (ii) on an earlier generation of water-rich asteroid, or (iii) during aqueous alteration on Allende itself [for example, Ross et al. (24) and Russell and MacPherson (25) and references therein]. Regardless of the location, dating of aqueous alteration products on meteorite parent bodies with extinct radionuclides 36Cl (t1/2 = 0.301 My), 26Al (t1/2 = 0.717 My), 53Mn (t1/2 = 3.74 My), or 129I (t1/2 = 15.7 My) suggests that it took place no later than 10 My after SS formation. This is a conservative upper limit because secondary alteration phases in fine-grained inclusions may have formed very early in the nebula. The half-life of 247Cm is 15.6 My, meaning that a time span of ~5 ± 5 My between SS formation and CAI alteration (and possible U depletion) translates into a correction for the initial 247Cm/235U ratio of 25 ± 25% (the uncertainty on this factor takes into account the uncertainty on the closure age). Our present best estimate of the initial SS 247Cm/235U ratio is thus (7.0 ± 1.6) × 10−5, which is equivalent to (247Cm/238U)ESS = (2.2 ± 0.5) × 10−5 and (247Cm/232Th)ESS = (9.7 ± 2.2) × 10−6. This value is in agreement with the 247Cm/235U ratio of (1.1 to 2.4) × 10−4 obtained by Brennecka et al. (11) based on CAI measurements and an upper limit of ~4 × 10−3 inferred from earlier meteoritic measurements (26). It is also in line with the lower estimate derived from modeling of galactic chemical evolution (GCE) (10), which predicts an initial ratio of (5.0 ± 2.5) × 10−5.

Fractionation of U isotopes during evaporation/condensation

Some of the samples with low 144Nd/238U ratios show significant scatter around the isochron (see inset of Fig. 1), outside of analytical uncertainties. In particular, samples with 144Nd/238U ratio as low as ~30 span a range of δ235U values of 3.5‰, a feature that led previous workers (17, 18) to question the conclusion of Brennecka et al. (11) that 247Cm was responsible for U isotope variations. This scatter is most likely due to isotopic fractionation during condensation, suggesting that better isochronous behavior could be obtained if such fractionation could be corrected.

Implications for the r-process

In addition to the dating implications that large 235U/238U variations have on the Pb-Pb ages of the CAIs (11, 21), the existence of 247Cm in the ESS has implications for the nucleosynthesis of r-process elements. The simplest GCE model that successfully reproduces the metallicity distribution of G-dwarfs requires infall of gas onto the galactic disk (10, 27, 28). In such a model, the abundance ratio of an SLR normalized to a stable nuclide isEmbedded Image(1)where N is the abundance, P is the production ratio, τ is the mean life of the SLR (τ = 1/λ), T is the age of the galaxy at the time when the ratio NSLR/NStable is to be calculated, and k is a constant that distinguishes closed-box (k = 0) versus infall models (typically, k = 2 ± 1) (8, 2729). If T is taken as the presolar age of the galaxy, T* (= TGTSS = 8.7 Gy), then the ratio NSLR/NStable in the interstellar medium (ISM) at the time of isolation of the protosolar molecular cloud from fresh nucleosynthetic input can be calculated. This value can be compared to the NSLR/NStable ratio in the ESS as obtained from meteoritic measurements. The difference between the two values is often interpreted as a “free-decay interval” (Δ) between the last nucleosynthetic event that produced the SLR and the formation of the SSEmbedded Image(2)At present, the meteoritic abundances of only three short-lived r-nuclides (129I, 182Hf, and 244Pu) have been estimated, yielding three different Δ values (100 ± 7 My, ~35 My, and 158 ± 85 My, respectively) (7). The limit of our understanding of the r-process is illustrated by the fact that there have been as many r-processes proposed as short-lived r-nuclides investigated.

Using the value of (247Cm/232Th)ESS = (9.7 ± 2.2) × 10−6 obtained in this study and the open nonlinear GCE model of Dauphas et al. (28) with k = 1.7, we obtained a free-decay interval of Δ = 98 ± 14 My. This value is in agreement with the Δ value of ~100 My derived from 129I and 244Pu but is much longer than the value of ~35 My obtained from 107Pd and 182Hf (Fig. 3). If 107Pd and 182Hf are indeed pure r-process isotopes, then a multiplicity of r-process environments is needed to explain the inconsistent Δ values (5, 7). This is evident when looking at Fig. 3: 107Pd and 129I have similar (NSLR/NStable)/(PSLR/PStable) ratios but different half-lives, making it impossible for both isotopes to be produced in a single event/process, no matter what model is used for the evolution of the SLR abundances in the giant molecular cloud parental to the SS [for example, free-decay interval versus three-phase mixing ISM (8, 30)]. However, recent nucleosynthetic models have reconsidered the origin of 107Pd and 182Hf and find a significant s-process contribution (70 to 80%) for both isotopes (6, 31). In such a framework, the initial abundance of all r-process SLR in the ESS can be explained by a single r-process making environment, which last injected material into the protosolar molecular cloud ~100 My before SS formation.

Fig. 3 Meteoritic abundance ratios of extinct radionuclides to stable nuclides produced by the same process (for example, 129I/127Ir), normalized to stellar production ratios versus mean lives (τ = t1/2 /ln2).

The superscript “r” refers to the r-process component of the cosmic abundance [obtained after subtracting the s-process contribution from solar abundances (3)]. When the normalizing isotope is not stable (for example, X/238U or X/232Th), the R/P ratio [that is, (NSLR/NStable)/(PSLR/PStable)] is corrected for the decay of the long-lived isotope by multiplication by the N/P ratio of the normalizing isotope in the ESS [0.71 for 238U and 0.89 for 232Th, values from Nittler and Dauphas (10)]. Triangles denote p-process isotopes, diamonds denote r-process nuclides with possible large s-process contributions, and circles denote r-process isotopes. The data are compared to model steady-state abundances in the ISM, using the model of Dauphas et al. (28) with k = 1.7 and a presolar age of the galaxy of 8.7 Gy. Dotted curves show model abundances assuming free-decay intervals of 10, 35, and 100 My, respectively. The abundances of all r-process nuclides can be explained by a single r-process environment that was active throughout the history of the galaxy and from which the solar system parent material was isolated about 100 My before SS formation, provided that 182Hf and 107Pd in meteorites originate from the s-process (6, 31). See the Supplementary Materials for source data.

Although some low metallicity halo stars point to the existence of at least three distinct r-processes, these stars formed from a gas that had been enriched in r-process nuclides by synthesis in very low metallicity stars formed early in the history of the galaxy. Meteorite evidence suggests that such r-process multiplicity may only be relevant to exotic conditions that prevailed in the earliest generation of stars in the history of the galaxy and that a single r-process may still be relevant to long-term models of the chemical evolution of the galaxy.


Twelve fine-grained and three coarse-grained CAIs were selected for this study (fig. S1). One of the coarse-grained CAIs (TS32) was obtained in powder form directly and was described in a previous publication (32). All other CAIs were identified in Allende slabs and extracted with clean stainless steel dental tools before digestion by acids. A small chip of each CAI (all but TS32) was extracted using clean stainless steel dental tools under a stereoscopic zoom microscope, and mounted in epoxy for characterization.

All samples were mapped using a JEOL JSM-5800LV scanning electron microscope (SEM). Images of a selected field of view of each CAI, secondary electron, backscattered electron, and false-color RGB (Mg/Ca/Al) are shown in figs. S2 to S15. The REE patterns of the samples were determined using a Quadrupole LA-ICPMS (laser ablation inductively coupled plasma mass spectrometry) at the Field Museum. Diagnostic group II REE patterns were identified in all 12 fine-grained CAIs, indicating that the Nd/U ratio (a proxy for the Cm/U ratio) of these samples was high and that these samples were therefore well suited to search for 235U excesses coming from 247Cm decay.

After digestion but before column chemistry, a small aliquot of each sample (~2 to 3%) was used to determine U and REE concentrations. The samples were spiked with IRMM-3636 U double spike (233U-236U) and processed through U/TEVA resin (Eichrom) column chromatography, using Optima-grade and sub-boiling double-distilled acids.

Isotopic analyses were performed on a ThermoFinnigan Neptune MC-ICPMS at the Origins Lab (University of Chicago), equipped with an OnToolBooster 150 Jet Pump (Pfeiffer) and using Jet sample cones and X skimmer cones. An Aridus II desolvating nebulizer was used for sample introduction, and the measurements were done using a static cup configuration [see (21)]. When a ~+55‰ 235U excess was discovered in the Curious Marie CAI, a new measurement setup was developed that allowed characterization of the U isotopic composition of 0.1 to 0.01 ng of U with a precision of ±2 to 3‰ (see the Supplementary Materials). Extensive testing was done to ensure that no systematic bias was affecting the data, the details of which can be found in the Supplementary Materials.


Supplementary material for this article is available at

Materials and Methods

Curious Marie: not a FUN CAI.

The 247Cm-235U chronometer and the initial abundance of 247Cm in the ESS.

GCE model

Fig. S1. Photos of typical fine-grained and coarse-grained CAIs.

Figs. S2 to S15. Secondary electron, backscattered electron, and false-color RGB maps of all samples.

Fig. S16. REE and U-Th abundance patterns of all 12 fine-grained CAIs analyzed in this study.

Fig. S17. Results of the standard addition measurements conducted on the Curious Marie CAI.

Fig. S18. U blank from new U/TEVA resin as a function of the volume of 0.05 M HCl passed through the column.

Fig. S19. Results of precision tests of U isotopic measurements done using various instrumental setups.

Fig. S20. Comparison of the δ235U determined on the CAIs using 80% (x axis) and 20% (y axis) of the sample.

Fig. S21. Flowchart of the tests conducted on the Curious Marie CAI.

Fig. S22. Evolution of the U isotopic composition of the gas, instantaneous solid, and cumulative solid, as a function of the fraction of U condensed.

Table S1. Results of SEM analysis on small chips of CAIs mounted in epoxy (wt %).

Table S2. Summary of U isotopic compositions and concentrations of CAIs and geostandards.

Table S3. Specifics of U isotopic measurements on MC-ICPMS for low U samples.

Table S4. Compilation of chemistry blanks and effect on U “stable” isotope ratio measurements.

Table S5. Summary of the Ti data obtained on geostandards and the Curious Marie CAI.

Table S6. Production ratios of selected SLRs produced by the s-, r-, and p-process and present in the ESS, normalized to a stable isotope produced in the same or similar nucleosynthetic process.

Table S7. Compilation of Cm/U isochron data and free-decay interval data from experimental and theoretical studies.

Table S8. Selected extinct radionuclides produced by the s-, r-, and p-process.

References (3878)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank P. Heck, J. Holstein, and the Robert A. Pritzker Center for Meteoritics and Polar Studies at the Field Museum for providing some of the CAI samples; I. Steele for help in operating the SEM at the University of Chicago; and L. Dussubieux for help in operating the LA-ICPMS at the Field Museum. We are grateful to A. M. Davis and B. S. Meyer for discussions. Funding: This work was supported by grants from NASA (Laboratory Analysis of Returned Samples, NNX14AK09G; Cosmochemistry, OJ-30381-0036A and NNX15AJ25G) and NSF (Petrology and Geochemistry, EAR144495; Cooperative Studies of the Earth’s Deep Interior, EAR150259) to N.D. and a NASA grant (Cosmochemistry, NNX13AE73G) to L.G. Author contributions: N.D. and F.L.H.T. designed the research; F.L.H.T. and L.G. selected the samples; F.L.H.T. performed the research; F.L.H.T., N.D., and L.G. wrote the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data are available from the authors upon request. This is Origins Lab contribution number 91.

Stay Connected to Science Advances

Navigate This Article