Abstract:Impact-induced erosion of the Earth's early crust during accretion of terrestrial bodies can significantly modify the primordial chemical composition of the Bulk Silicate Earth (BSE, that is, the composition of the crust added to the present-day mantle). In particular, it can be particularly efficient in altering the abundances of elements having a strong affinity for silicate melts (i.e. incompatible elements) as the early differentiated crust was preferentially enriched in those. Here, we further develop an erosion model (EROD) to quantify the effects of collisional erosion on the final composition of the BSE. Results are compared to the present-day BSE composition models and constraints on Earth's accretion processes are provided. The evolution of the BSE chemical composition resulting from crustal stripping is computed for entire accretion histories of about 50 Earth analogs in the context of the Grand Tack model. The chosen chemical elements span a wide range of incompatibility degrees. We find that a maximum loss of 40wt% can be expected for the most incompatible lithophile elements such as Rb, Th or U in the BSE when the crust is formed from low partial melting rates. Accordingly, depending on both the exact nature of the crust-forming processes during accretion and the accretion history itself, Refractory Lithophile Elements (RLE) may not be in chondritic relative proportions in the BSE. In that case, current BSE estimates may need to be corrected as a function of the geochemical incompatibility of these elements. Alternatively, if RLE are indeed in chondritic relative proportions in the BSE, accretion scenarios that are efficient in affecting the BSE chemical composition should be questioned.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the impact of collisional erosion on the composition of the Bulk Silicate Earth (BSE) during the formation of Earth - like planets, especially in the context of the "Grand Tack" model. Specifically, the study focuses on how impact - induced erosion of the early crust during the accretion of planetary bodies can significantly change the original chemical composition of the BSE, especially the abundance changes of incompatible elements (such as Rb, Th or U) that have a strong affinity for silicate melts. By developing an erosion model (EROD), the impact of collisional erosion on the final composition of the BSE is quantified and compared with models of the modern BSE composition to provide constraints on the Earth's accretion process.
### Research Background
The chemical compositions of the Earth and other inner - solar - system rocky bodies are inherited from nebular and post - nebular processes. Traditionally, the compositional gradients of solar - system bodies are mainly explained by the condensation sequence of elements under nebular conditions in the early stages of solar - system formation. However, post - nebular processes (such as evaporation, differentiation, collisional erosion) may lead to further chemical fractionation among planetary bodies. These processes are considered to occur in the final stages of planet formation, including the so - called giant - planet - collision stage, which begins 2 - 3 million years after the formation of CAIs and lasts until 200 million years. This stage begins before the gas in the disk disappears, if the gas disappears approximately 3 million years later.
### Research Objectives
This paper aims to quantitatively simulate the impact of collisional erosion on the final composition of the BSE. The study selects a series of chemical elements with different geochemical properties, including different degrees of incompatibility, to explore the potential fractionation of these elements in the BSE. The study specifically focuses on the following points:
- Tracking the impact - induced fractionation of different RLEs pairs (such as Lu/Hf, Sm/Nd, etc.) in Earth - like planets within about 100 - 200 million years.
- Estimating the loss of individual RLEs relative to the original material.
- In - depth study of the popular "Grand Tack" scenario, especially the impact of the date of the Moon - formation event.
- Considering the impact of small planetary bodies.
- Calculating the coupling between all embryos to consider their own fractionation and their further impact in impacts.
- Exploring the impact of the partial melting rate in the crust - formation history.
### Methods
Use N - body numerical simulations to track the collision history of each embryo, and combine with analytical models to calculate the mass transfer in each impact. In addition, perform geochemical mass balance to recalculate the composition of the embryo after each impact. By integrating the entire collision history, the final fractionation of various elements in the crust and mantle of each embryo can be quantified.
### Results
The research results show that for a selected series of elements, the degree of fractionation caused by collisional erosion varies. For example, the fractionation of Rb/Sr is the most obvious, with an average fractionation value of \(\epsilon=- 0.16\pm0.06\); while the fractionation of Th/U is the smallest, with an average fractionation value of \(\epsilon = - 0.008\pm5\times10^{-4}\). These results are consistent with the ratio of liquid - solid partition coefficients of these element pairs. A negative value indicates the preferential loss of the numerator element, and a positive value indicates the preferential loss of the denominator element.
### Conclusions
This study provides a new perspective for understanding the role of collisional erosion in the formation process of Earth - like planets, especially the possible impact on the proportion of RLEs in the BSE. The research results are of great significance for evaluating the rationality of current BSE model assumptions and exploring the Earth - formation scenarios in which these assumptions remain unchanged.