Tungsten isotope evolution during Earth’s formation and new constraints on the viability of accretion simulations

1D.C. Rubie, 2K.I. Dale, 3G. Nathan, 4M. Nakajima, 5E.S. Jennings, 1G.J. Golabek, 3S.A. Jacobson, 2,6A. Morbidelli
Earth and Planetary Science Letters 651, 119139 Open Access Link to Article [https://doi.org/10.1016/j.epsl.2024.119139]
1Bayerisches Geoinstitut, University of Bayreuth, Universitätsstrasse 30, 95440 Bayreuth, Germany
2Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, Nice, France
3Department of Earth & Environmental Sciences, Michigan State University, 288 Farm Lane, East Lansing, MI 48823, USA
4Department of Earth and Environmental Sciences, University of Rochester, 227 Hutchison Hall, Rochester, NY 14627, USA
5School of Natural Sciences, Birkbeck, University of London, Malet Street, London WC1E 7HX, UK
6Collège de France, CNRS, PSL Univ., Sorbonne Univ., Paris, 75014, France
Copyright Elsevier

The Hf-W isotopic system is the reference chronometer for determining the chronology of Earth’s accretion and differentiation. However, its results depend strongly on uncertain parameters, including the extent of metal-silicate equilibration and the siderophility of tungsten. Here we show that a multistage core-formation model based on N-body accretion simulations, element mass balance and metal-silicate partitioning, largely eliminates these uncertainties. We modified the original model of Rubie et al. (2015) by including (1) smoothed particle hydrodynamics estimates of the depth of melting caused by giant impacts and (2) the isotopic evolution of 182W. We applied two metal-silicate fractionation mechanisms: one when the metal delivered by the cores of large impactors equilibrates with only a small fraction of the impact-induced magma pond and the other when metal delivered by small impactors emulsifies in global magma oceans before undergoing progressive segregation. The latter is crucial for fitting the W abundance and 182W anomaly of Earth’s mantle. In addition, we show, for the first time, that the duration of magma ocean solidification has a major effect on Earth’s tungsten isotope anomaly. We re-evaluate the six Grand Tack N-body simulations of Rubie et al. (2015). Only one reproduces ε182W=1.9 ± 0.1 of Earth’s mantle, otherwise accretion is either too fast or too slow. Depending on the characteristics of the giant impacts, results predict that the Moon formed either 143–183 Myr or 53–62 Myr after the start of the solar system. Thus, independent evaluations of the Moon’s age provide an additional constraint on the validity of accretion simulations.

Discuss