A nucleosynthetic origin for the Earth's anomalous (142)Nd composition

Testing the Ce Limit of Mass Bias Correction Using 145Nd/142Nd as Normalizing Ratio in Radiogenic Neodymium Isotope Analysis by MC-ICP-MS

RATIONALE

Neodymium isotopes are a powerful geochemical tool that has widely been used in terrestrial and extraterrestrial studies. Modern multicollector inductively coupled plasma mass spectrometers (MC-ICP-MS) allow fast, accurate, and precise analysis of the radiogenic Nd isotope ratio 143Nd/144Nd. These analyses comprise relatively high instrumental mass bias that is typically corrected for using the stable 146Nd/144Nd of 0.7219 and an exponential law. The instrument is usually tuned to optimize the operating conditions for isotope analysis, but this tuning is a trade-off primarily between signal intensity, stability, and accuracy. Alternative, more effective approaches for mass bias correction have been proposed that use 145Nd/142Nd as normalizing ratio. However, one drawback of using this ratio is that the efficient removal of Ce from Nd is required to minimize the effect of isobaric interference of 142Ce on 142Nd.

METHODS

Here, we analyzed international Nd and rock reference materials using a Thermo Scientific Neptune Plus MC-ICP-MS to evaluate the sensitivity of 145Nd/142Nd-based mass bias correction to varying Ce/Nd and in comparison with the commonly used 146Nd/144Nd-based correction.

RESULTS

Our results show that the corrected 143Nd/144Nd of Ce-doped JNdi-1 and Ce-containing USGS BCR-2, NOD-A-1, and NOD-P-1 reference materials are insensitive to Ce/Nd of up ~1.

CONCLUSIONS

The correction of instrumental mass bias with 145Nd/142Nd as a normalizing ratio yields, as previously suggested, improved trueness and precision of 143Nd/144Nd data in comparison with 146Nd/144Nd-based corrections, even under high Ce/Nd of to ~1. This allows improved optimization of signal intensity during instrument tuning, which is particularly useful for natural samples with low Nd content.

High-precision Sm isotope analysis by thermal ionisation mass spectrometry for large meteorite samples (>1 g)

This study presents a new procedure for high-precision Sm isotope analysis by thermal ionisation mass spectrometry (TIMS) for geological samples. A four-step chemical separation scheme results in sharp separation of Sm and Nd from the same sample aliquot. The first step utilises anion exchange resin to remove Fe from the sample solution. Two different liquid-liquid extraction resins are then used to isolate rare-earth elements (TRU-Spec) and purify Sm from Nd (DGA). Fractionation occurs on the DGA resin due to the nuclear field shift effect, but this is negligible if yields greater than 70% are achieved. Different analytical setups were tested to ascertain their ionisation efficiencies on TIMS. The effect of activators composed of Pt and Ta was tested on single Re filaments but the conventional double Re filament assembly provided efficient ionisation and more stable ion beams. The determination of nucleosynthetic isotope variations requires high precision for all Sm isotope ratios. We aimed to improve the precision on the scarce 144Sm isotope (3% of all Sm). Static, multistatic and dynamic methods were tested. Isotope ratios were normalised to both 147Sm/152Sm and 152Sm/148Sm for comparison. The dynamic methods failed to provide better precision on ratios involving 144Sm, whereas the multistatic method yielded improved precisions between 13 and 22 ppm (twice the standard deviation, 2 SD) on the 144Sm/152Sm ratio. Synthetic standards have variable Sm isotope compositions, thus requiring systematic and precise characterisation against terrestrial samples. Analyses conducted using this new procedure yielded high-precision values which were consistent with literature data for an array of terrestrial rock standards and the meteorite Allende.

The Ni isotopic composition of Ryugu reveals a common accretion region for carbonaceous chondrites

The isotopic compositions of samples returned from Cb-type asteroid Ryugu and Ivuna-type (CI) chondrites are distinct from other carbonaceous chondrites, which has led to the suggestion that Ryugu/CI chondrites formed in a different region of the accretion disk, possibly around the orbits of Uranus and Neptune. We show that, like for Fe, Ryugu and CI chondrites also have indistinguishable Ni isotope anomalies, which differ from those of other carbonaceous chondrites. We propose that this unique Fe and Ni isotopic composition reflects different accretion efficiencies of small FeNi metal grains among the carbonaceous chondrite parent bodies. The CI chondrites incorporated these grains more efficiently, possibly because they formed at the end of the disk's lifetime, when planetesimal formation was also triggered by photoevaporation of the disk. Isotopic variations among carbonaceous chondrites may thus reflect fractionation of distinct dust components from a common reservoir, implying CI chondrites/Ryugu may have formed in the same region of the accretion disk as other carbonaceous chondrites.

Dating the Solar System's giant planet orbital instability using enstatite meteorites

The giant planets of the Solar System formed on initially compact orbits, which transitioned to the current wider configuration by means of an orbital instability. The timing of that instability is poorly constrained. In this work, we use dynamical simulations to demonstrate that the instability implanted planetesimal fragments from the terrestrial planet region into the asteroid main belt. We use meteorite data to show that the implantation occurred >60 million years (Myr) after the Solar System began to form. Combining this constraint with a previous upper limit derived from Jupiter's trojan asteroids, we conclude that the orbital instability occurred 60 to 100 Myr after the beginning of Solar System formation. The giant impact that formed the Moon occurred within this range, so it might be related to the giant planet instability.

The combined Hf and Nd isotope evolution of the depleted mantle requires Hadean continental formation

The onset and rates of continental growth are first-order indicators of early Earth dynamics, and whether substantial crust existed in the Hadean or much later has long been debated. Here, we present a theoretical analysis of published Hf and Nd isotopic data representing the depleted mantle and demonstrate that continental growth must have started in the early Hadean. Whereas the traditional interpretation of depleted mantle signatures in crustal rocks assumes unrealistic instantaneous mantle mixing, our modeling incorporates the effect of a finite mixing time over which these signatures are recorded in rocks produced through mantle melting. This effect is shown to delay, by as much as 0.65 to 0.75 billion years, the appearance of the earliest depleted mantle signatures in continental crust. Our results suggest that published observations of εHf, ε¹⁴³Nd, and μ¹⁴²Nd require Hadean growth of continental crust, with a minimum of 50% of today's continental volume already existing by the end of Hadean.

Nd isotope variation between the Earth-Moon system and enstatite chondrites

Reconstructing the building blocks that made Earth and the Moon is critical to constrain their formation and compositional evolution to the present. Neodymium (Nd) isotopes identify these building blocks by fingerprinting nucleosynthetic components. In addition, the ¹⁴⁶Sm-¹⁴²Nd and ¹⁴⁷Sm-¹⁴³Nd decay systems, with half-lives of 103 million years and 108 billion years, respectively, track potential differences in their samarium (Sm)/Nd ratios. The difference in Earth's present-day ¹⁴²Nd/¹⁴⁴Nd ratio compared with chondrites^1,2, and in particular enstatite chondrites, is interpreted as nucleosynthetic isotope variation in the protoplanetary disk. This necessitates that chondrite parent bodies have the same Sm/Nd ratio as Earth's precursor materials². Here we show that Earth and the Moon instead had a Sm/Nd ratio approximately 2.4 ± 0.5 per cent higher than the average for chondrites and that the initial ¹⁴²Nd/¹⁴⁴Nd ratio of Earth's precursor materials is more similar to that of enstatite chondrites than previously proposed^1,2. The difference in the Sm/Nd ratio between Earth and chondrites probably reflects the mineralogical distribution owing to mixing processes within the inner protoplanetary disk. This observation simplifies lunar differentiation to a single stage from formation to solidification of a lunar magma ocean³. This also indicates that no Sm/Nd fractionation occurred between the materials that made Earth and the Moon in the Moon-forming giant impact.

Earth's composition was modified by collisional erosion

The samarium-146 (¹⁴⁶Sm)-neodymium-142 (¹⁴²Nd) short-lived decay system (half-life of 103 million years) is a powerful tracer of the early mantle-crust evolution of planetary bodies. However, an increased ¹⁴²Nd/¹⁴⁴Nd in modern terrestrial rocks relative to chondrite meteorites has been proposed to be caused by nucleosynthetic anomalies, obscuring early Earth's differentiation history. We use stepwise dissolution of primitive chondrites to quantify nucleosynthetic contributions on the composition of chondrites. After correction for nucleosynthetic anomalies, Earth and the silicate parts of differentiated planetesimals contain resolved excesses of ¹⁴²Nd relative to chondrites. We conclude that only collisional erosion of primordial crusts can explain such compositions. This process associated with planetary accretion must have produced substantial loss of incompatible elements, including long-term heat-producing elements such as uranium, thorium, and potassium.

Half-life and initial Solar System abundance of 146Sm determined from the oldest andesitic meteorite

¹⁴⁶Sm-¹⁴²Nd radioactive systematics can provide constraints on the timing of early differentiation processes on Earth, Moon, and Mars. The uncertainties related to the initial abundance and half-life of the extinct isotope ¹⁴⁶Sm impede the interpretation of the ¹⁴⁶Sm-¹⁴²Nd systematics of planetary materials. The accurate determinations of Sm, Nd, and Mg isotopic compositions of the oldest “andesitic” achondrite Erg Chech 002 (EC 002) define a crystallization age of 1.8 Myr after the formation of the Solar System and provide the most accurate and reliable initial ratio of ¹⁴⁶Sm/¹⁴⁴Sm for the Solar System at 0.00840 ± 0.00032 using a ¹⁴⁶Sm half-life of 103 Ma, making EC 002 an anchor for ¹⁴⁶Sm-¹⁴²Nd systematics for Earth and planetary materials.

The formation and differentiation of planetary bodies are dated using radioactive decay systems, including the short-lived ¹⁴⁶Sm-¹⁴²Nd (T_½ = 103 or 68 Ma) and long-lived ¹⁴⁷Sm-¹⁴³Nd (T_½ = 106 Ga) radiogenic pairs that provide relative and absolute ages, respectively. However, the initial abundance and half-life of the extinct radioactive isotope ¹⁴⁶Sm are still debated, weakening the interpretation of ¹⁴⁶Sm-¹⁴²Nd systematics obtained for early planetary processes. Here, we apply the short-lived ²⁶Al-²⁶Mg, ¹⁴⁶Sm-¹⁴²Nd, and long-lived ¹⁴⁷Sm-¹⁴³Sm chronometers to the oldest known andesitic meteorite, Erg Chech 002 (EC 002), to constrain the Solar System initial abundance of ¹⁴⁶Sm. The ²⁶Al-²⁶Mg mineral isochron of EC 002 provides a tightly constrained initial δ²⁶Mg* of −0.009 ± 0.005 ‰ and (²⁶Al/²⁷Al)₀ of (8.89 ± 0.09) × 10⁻⁶. This initial abundance of ²⁶Al is the highest measured so far in an achondrite and corresponds to a crystallization age of 1.80 ± 0.01 Myr after Solar System formation. The ¹⁴⁶Sm-¹⁴²Nd mineral isochron returns an initial ¹⁴⁶Sm/¹⁴⁴Sm ratio of 0.00830 ± 0.00032. By combining the Al-Mg crystallization age and initial ¹⁴⁶Sm/¹⁴⁴Sm ratio of EC 002 with values for refractory inclusions, achondrites, and lunar samples, the best-fit half-life for ¹⁴⁶Sm is 102 ± 9 Ma, corresponding to the physically measured value of 103 ± 5 Myr, rather than the latest and lower revised value of 68 ± 7 Ma. Using a half-life of 103 Ma for ¹⁴⁶Sm, the ¹⁴⁶Sm/¹⁴⁴Sm abundance of EC 002 translates into an initial Solar System ¹⁴⁶Sm/¹⁴⁴Sm ratio of 0.00840 ± 0.00032, which represents the most reliable and precise estimate to date and makes EC 002 an ideal anchor for the ¹⁴⁶Sm-¹⁴²Nd clock.

Terrestrial planet formation from lost inner solar system material

Two fundamentally different processes of rocky planet formation exist, but it is unclear which one built the terrestrial planets of the solar system. They formed either by collisions among planetary embryos from the inner solar system or by accreting sunward-drifting millimeter-sized “pebbles” from the outer solar system. We show that the isotopic compositions of Earth and Mars are governed by two-component mixing among inner solar system materials, including material from the innermost disk unsampled by meteorites, whereas the contribution of outer solar system material is limited to a few percent by mass. This refutes a pebble accretion origin of the terrestrial planets but is consistent with collisional growth from inner solar system embryos. The low fraction of outer solar system material in Earth and Mars indicates the presence of a persistent dust-drift barrier in the disk, highlighting the specific pathway of rocky planet formation in the solar system.

Nickel isotopic evidence for late-stage accretion of Mercury-like differentiated planetary embryos

Earth's habitability is closely tied to its late-stage accretion, during which impactors delivered the majority of life-essential volatiles. However, the nature of these final building blocks remains poorly constrained. Nickel (Ni) can be a useful tracer in characterizing this accretion as most Ni in the bulk silicate Earth (BSE) comes from the late-stage impactors. Here, we apply Ni stable isotope analysis to a large number of meteorites and terrestrial rocks, and find that the BSE has a lighter Ni isotopic composition compared to chondrites. Using first-principles calculations based on density functional theory, we show that core-mantle differentiation cannot produce the observed light Ni isotopic composition of the BSE. Rather, the sub-chondritic Ni isotopic signature was established during Earth's late-stage accretion, probably through the Moon-forming giant impact. We propose that a highly reduced sulfide-rich, Mercury-like body, whose mantle is characterized by light Ni isotopic composition, collided with and merged into the proto-Earth during the Moon-forming giant impact, producing the sub-chondritic Ni isotopic signature of the BSE, while delivering sulfur and probably other volatiles to the Earth.

Remnants of early Earth differentiation in the deepest mantle-derived lavas

The noble gas isotope systematics of ocean island basalts suggest the existence of primordial mantle signatures in the deep mantle. Yet, the isotopic compositions of lithophile elements (Sr, Nd, Hf) in these lavas require derivation from a mantle source that is geochemically depleted by melt extraction rather than primitive. Here, this apparent contradiction is resolved by employing a compilation of the Sr, Nd, and Hf isotope composition of kimberlites-volcanic rocks that originate at great depth beneath continents. This compilation includes kimberlites as old as 2.06 billion years and shows that kimberlites do not derive from a primitive mantle source but sample the same geochemically depleted component (where geochemical depletion refers to ancient melt extraction) common to most oceanic island basalts, previously called PREMA (prevalent mantle) or FOZO (focal zone). Extrapolation of the Nd and Hf isotopic compositions of the kimberlite source to the age of Earth formation yields a ¹⁴³Nd/¹⁴⁴Nd-¹⁷⁶Hf/¹⁷⁷Hf composition within error of chondrite meteorites, which include the likely parent bodies of Earth. This supports a hypothesis where the source of kimberlites and ocean island basalts contains a long-lived component that formed by melt extraction from a domain with chondritic ¹⁴³Nd/¹⁴⁴Nd and ¹⁷⁶Hf/¹⁷⁷Hf shortly after Earth accretion. The geographic distribution of kimberlites containing the PREMA component suggests that these remnants of early Earth differentiation are located in large seismically anomalous regions corresponding to thermochemical piles above the core-mantle boundary. PREMA could have been stored in these structures for most of Earth's history, partially shielded from convective homogenization.

Astronomical context of Solar System formation from molybdenum isotopes in meteorite inclusions

Calcium-aluminum-rich inclusions (CAIs) in meteorites are the first solids to have formed in the Solar System, defining the epoch of its birth on an absolute time scale. This provides a link between astronomical observations of star formation and cosmochemical studies of Solar System formation. We show that the distinct molybdenum isotopic compositions of CAIs cover almost the entire compositional range of material that formed in the protoplanetary disk. We propose that CAIs formed while the Sun was in transition from the protostellar to pre-main sequence (T Tauri) phase of star formation, placing Solar System formation within an astronomical context. Our results imply that the bulk of the material that formed the Sun and Solar System accreted within the CAI-forming epoch, which lasted less than 200,000 years.

Potassium isotope anomalies in meteorites inherited from the protosolar molecular cloud

Potassium (K) and other moderately volatile elements are depleted in many solar system bodies relative to CI chondrites, which closely match the composition of the Sun. These depletions and associated isotopic fractionations were initially believed to result from thermal processing in the protoplanetary disk, but so far, no correlation between the K depletion and its isotopic composition has been found. Our new high-precision K isotope data correlate with other neutron-rich nuclides (e.g., 64Ni and 54Cr) and suggest that the observed 41K variations have a nucleosynthetic origin. We propose that K isotope anomalies are inherited from an isotopically heterogeneous protosolar molecular cloud, and were preserved in bulk primitive meteorites. Thus, the heterogeneous distribution of both refractory and moderately volatile elements in chondritic meteorites points to a limited radial mixing in the protoplanetary disk.

Combined mass-dependent and nucleosynthetic isotope variations in refractory inclusions and their mineral separates to determine their original Fe isotope compositions

Calcium-aluminum-rich inclusions (CAIs) are the oldest dated materials that provide crucial information about the isotopic reservoirs present in the early Solar System. For a variety of elements, CAIs have isotope compositions that are uniform yet distinct from later formed solid material. However, despite being the most abundant metal in the Solar System, the isotopic composition of Fe in CAIs is not well constrained. In an attempt to determine the Fe isotopic compositions of CAIs, we combine extensive work from a previously studied CAI sample set with new isotopic work characterizing mass-dependent and mass-independent (nucleosynthetic) signatures in Mg, Ca, and Fe. This investigation includes work on three mineral separates of the Allende CAI Egg 2. For all isotope systems investigated, we find that in general, fine-grained CAIs exhibit light mass-dependent isotopic signatures relative to terrestrial standards, whereas igneous CAIs have heavier isotopic compositions relative to the fine-grained CAIs. Importantly, the mass-dependent Fe isotope signatures of bulk CAIs show a range of both light (fine-grained CAIs) and heavy (igneous CAIs) isotopic signatures relative to bulk chondrites, suggesting that Fe isotope signatures in CAIs largely derive from mass fractionation events such as condensation and evaporation occurring in the nebula. Such signatures show that a significant portion of the secondary alteration experienced by CAIs, particularly prevalent in fine-grained inclusions, occurred in the nebula prior to accretion into their respective parent bodies. Regarding nucleosynthetic Fe isotope signatures, we do not observe any variation outside of analytical uncertainty in bulk CAIs compared to terrestrial standards. In contrast, all three Egg 2 mineral separates display resolved mass-independent excesses in ⁵⁶Fe compared to terrestrial standards. Furthermore, we find that the combined mass-dependent and nucleosynthetic Fe isotopic compositions of the Egg 2 mineral separates are well correlated, likely indicating that Fe indigenous to the CAI is mixed with less anomalous Fe, presumably from the solar nebula. Thus, these reported nucleosynthetic anomalies may point in the direction of the original Fe isotope composition of the CAI-forming region, but they likely only provide a minimum isotopic difference between the original mass-independent Fe isotopic composition of CAIs and that of later formed solids.

Hadean silicate differentiation preserved by anomalous 142Nd/144Nd ratios in the Réunion hotspot source

Active volcanic hotspots can tap into domains in Earth's deep interior that were formed more than two billion years ago. High-precision data on variability in tungsten isotopes have shown that some of these domains resulted from differentiation events that occurred within the first fifty million years of Earth history. However, it has not proved easy to resolve analogous variability in neodymium isotope compositions that would track regions of Earth's interior whose composition was established by events occurring within roughly the first five hundred million years of Earth history. Here we report ¹⁴²Nd/¹⁴⁴Nd ratios for Réunion Island igneous rocks, some of which are resolvably either higher or lower than the ratios in modern upper-mantle domains. We also find that Réunion ¹⁴²Nd/¹⁴⁴Nd ratios correlate with helium-isotope ratios (³He/⁴He), suggesting parallel behaviour of these isotopic systems during very early silicate differentiation, perhaps as early as 4.39 billion years ago. The range of ¹⁴²Nd/¹⁴⁴Nd ratios in Réunion basalts is inconsistent with a single-stage differentiation process, and instead requires mixing of a conjugate melt and residue formed in at least one melting event during the Hadean eon, 4.56 billion to 4 billion years ago. Efficient post-Hadean mixing nearly erased the ancient, anomalous ¹⁴²Nd/¹⁴⁴Nd signatures, and produced the relatively homogeneous ¹⁴³Nd/¹⁴⁴Nd composition that is characteristic of Réunion basalts. Our results show that Réunion magmas tap into a particularly ancient, primitive source compared with other volcanic hotspots, offering insight into the formation and preservation of ancient heterogeneities in Earth's interior.

Radial mixing and Ru-Mo isotope systematics under different accretion scenarios

The Ru-Mo isotopic compositions of inner Solar System bodies may reflect the provenance of accreted material and how it evolved with time, both of which are controlled by the accretion scenario these bodies experienced. Here we use a total of 116 N-body simulations of terrestrial planet accretion, run in the Eccentric Jupiter and Saturn (EJS), Circular Jupiter and Saturn (CJS), and Grand Tack scenarios, to model the Ru-Mo anomalies of Earth, Mars, and Theia analogues. This model starts by applying an initial step function in Ru-Mo isotopic composition, with compositions reflecting those in meteorites, and traces compositional evolution as planets accrete. The mass-weighted provenance of the resulting planets reveals more radial mixing in Grand Tack simulations than in EJS/CJS simulations, and more efficient mixing among late-accreted material than during the main phase of accretion in EJS/CJS simulations. We find that an extensive homogenous inner disk region is required to reproduce Earth's observed Ru-Mo composition. EJS/CJS simulations require a homogeneous reservoir in the inner disk extending to ≥3-4 AU (≥74-98% of initial mass) to reproduce Earth's composition, while Grand Tack simulations require a homogeneous reservoir extending to ≥3-10 AU (≥97-99% of initial mass), and likely to ≥6-10 AU. In the Grand Tack model, Jupiter's initial location (the most likely location for a discontinuity in isotopic composition) is ~3.5 AU; however, this step location has only a 33% likelihood of producing an Earth with the correct Ru-Mo isotopic signature for the most plausible model conditions. Our results give the testable predictions that Mars has zero Ru anomaly and small or zero Mo anomaly, and the Moon has zero Mo anomaly. These predictions are insensitive to wide variations in parameter choices.

Tungsten Isotopes in Planets

The short-lived Hf-W isotope system has a wide range of important applications in cosmochemistry and geochemistry. The siderophile behavior of W, combined with the lithophile nature of Hf, makes the system uniquely useful as a chronometer of planetary accretion and differentiation. Tungsten isotopic data for meteorites show that the parent bodies of some differentiated meteorites accreted within 1 million years after Solar System formation. Melting and differentiation on these bodies took ~1-3 million years and was fueled by decay of 26Al. The timescale for accretion and core formation increases with planetary mass and is ~10 million years for Mars and >34 million years for Earth. The nearly identical 182W compositions for the mantles of the Moon and Earth are difficult to explain in current models for the formation of the Moon. Terrestrial samples with ages spanning ~4 billion years reveal small 182W variations within the silicate Earth, demonstrating that traces of Earth's earliest formative period have been preserved throughout Earth's history.

The isotopic nature of the Earth's accreting material through time

The Earth formed by accretion of Moon- to Mars-size embryos coming from various heliocentric distances. The isotopic nature of these bodies is unknown. However, taking meteorites as a guide, most models assume that the Earth must have formed from a heterogeneous assortment of embryos with distinct isotopic compositions. High-precision measurements, however, show that the Earth, the Moon and enstatite meteorites have almost indistinguishable isotopic compositions. Models have been proposed that reconcile the Earth-Moon similarity with the inferred heterogeneous nature of Earth-forming material, but these models either require specific geometries for the Moon-forming impact or can explain only one aspect of the Earth-Moon similarity (that is, ¹⁷O). Here I show that elements with distinct affinities for metal can be used to decipher the isotopic nature of the Earth's accreting material through time. I find that the mantle signatures of lithophile O, Ca, Ti and Nd, moderately siderophile Cr, Ni and Mo, and highly siderophile Ru record different stages of the Earth's accretion; yet all those elements point to material that was isotopically most similar to enstatite meteorites. This isotopic similarity indicates that the material accreted by the Earth always comprised a large fraction of enstatite-type impactors (about half were E-type in the first 60 per cent of the accretion and all of the impactors were E-type after that). Accordingly, the giant impactor that formed the Moon probably had an isotopic composition similar to that of the Earth, hence relaxing the constraints on models of lunar formation. Enstatite meteorites and the Earth were formed from the same isotopic reservoir but they diverged in their chemical evolution owing to subsequent fractionation by nebular and planetary processes.

Early formation of the Moon 4.51 billion years ago

Establishing the age of the Moon is critical to understanding solar system evolution and the formation of rocky planets, including Earth. However, despite its importance, the age of the Moon has never been accurately determined. We present uranium-lead dating of Apollo 14 zircon fragments that yield highly precise, concordant ages, demonstrating that they are robust against postcrystallization isotopic disturbances. Hafnium isotopic analyses of the same fragments show extremely low initial ¹⁷⁶Hf/¹⁷⁷Hf ratios corrected for cosmic ray exposure that are near the solar system initial value. Our data indicate differentiation of the lunar crust by 4.51 billion years, indicating the formation of the Moon within the first ~60 million years after the birth of the solar system.