Earth and Mars - distinct inner solar system products

Page created by Mildred Salinas
 
CONTINUE READING
Earth and Mars - distinct inner solar system products
Earth and Mars – distinct inner solar system products
                                                                         Takashi Yoshizaki∗1 and William F. McDonough1,2,3
                                               1
                                                   Department of Earth Science, Graduate School of Science, Tohoku University, Sendai,
arXiv:2006.11051v2 [astro-ph.EP] 2 Feb 2021

                                                                                  Miyagi 980-8578, Japan
                                                      2
                                                        Department of Geology, University of Maryland, College Park, MD 20742, USA
                                                   3
                                                     Research Center for Neutrino Science, Tohoku University, Sendai, Miyagi 980-8578,
                                                                                           Japan

                                                                                               February 3, 2021

                                              Abstract
                                              Composition of terrestrial planets records planetary accretion, core-mantle and crust-mantle differ-
                                              entiation, and surface processes. Here we compare the compositional models of Earth and Mars
                                              to reveal their characteristics and formation processes. Earth and Mars are equally enriched in
                                              refractory elements (1.9 × CI), although Earth is more volatile-depleted and less oxidized than
                                              Mars. Their chemical compositions were established by nebular fractionation, with negligible con-
                                              tributions from post-accretionary losses of moderately volatile elements. The degree of planetary
                                              volatile element depletion might correlate with the abundances of chondrules in the accreted ma-
                                              terials, planetary size, and their accretion timescale, which provides insights into composition and
                                              origin of Mercury, Venus, the Moon-forming giant impactor, and the proto-Earth. During its for-
                                              mation before and after the nebular disk’s lifetime, the Earth likely accreted more chondrules and
                                              less matrix-like materials than Mars and chondritic asteroids, establishing its marked volatile de-
                                              pletion. A giant impact of an oxidized, differentiated Mars-like (i.e., composition and mass) body
                                              into a volatile-depleted, reduced proto-Earth produced a Moon-forming debris ring with mostly a
                                              proto-Earth’s mantle composition. Chalcophile and some siderophile elements in the silicate Earth
                                              added by the Mars-like impactor were extracted into the core by a sulfide melt (∼0.5% of the mass
                                              of the Earth’s mantle). In contrast, the composition of Mars indicates its rapid accretion of lesser
                                              amounts of chondrules under nearly uniform oxidizing conditions. Mars’ rapid cooling and early
                                              loss of its dynamo likely led to the absence of plate tectonics and surface water, and the present-
                                              day low surface heat flux. These similarities and differences between the Earth and Mars made the
                                              former habitable and the other inhospitable to uninhabitable.
                                                   ∗
                                                       Corresponding author. E-mail: takashiy@tohoku.ac.jp

                                                                                                         1
Earth and Mars - distinct inner solar system products
Keywords
Earth; Mars; chondrites; solar system; cosmochemistry

1    Introduction
Earth and Mars share many chemical and physical attributes, but are distinct in size and inventory of
volatile elements. Chemical data from surface rocks and meteorites combined with seismological
and geodetic observations (e.g., mass, density, moment of inertia (MOI)) provide a multiply con-
strained compositional model of the terrestrial planets (e.g., Ringwood, 1966; Morgan and Anders,
1980; Wänke, 1981; Longhi et al., 1992). The physicochemical similarities and differences between
Earth and Mars provide useful insights on formation and evolution processes of these bodies, and
potentially provide insights into the present-day properties and origin of other rocky bodies (Venus,
Mercury, and exoplanets).

     By definition, Earth and Mars are located within a habitable zone in which liquid water is avail-
able (Cockell et al., 2016; Ehlmann et al., 2016). However, currently life exists only on Earth, and
it remains unclear if Mars were inhabited or uninhabited in its history. The occurrence of life only
on Earth indicates that there are compositional limits for a habitable planet formation, which can
be revealed by an Earth-Mars comparison. Previous comparative studies of chemical compositions
of the interiors of Earth and Mars described their differences including abundances of volatile or
siderophile elements and redox states (Anders and Owen, 1977; Dreibus and Wänke, 1987; Wänke
and Dreibus, 1994), but the recent advances in compositional modeling of both Earth and Mars will
provide further insights into their characters and origins.

    Isotopic compositions of solar system materials provide strong constraints on the sources of
the terrestrial planets. Mass-independent, nucleosynthetic isotope anomalies in meteorites reveal
a heterogeneous distribution of distinct presolar materials and provide the basis for classifying the
non-carbonaceous (NC) and carbonaceous meteorite (CC) groups (Trinquier et al., 2007, 2009;
Warren, 2011; Kruijer et al., 2017a). The NC meteorites appear to be from the inner solar system,
whereas the CC meteorites are considered to be samples from the outer solar system, including the
Trojan asteroids and possibly beyond Jupiter’s orbit (Walsh et al., 2011; Kruijer et al., 2017a). These
isotopic observations demonstrate a link between Earth and enstatite chondrites, whereas Mars is
thought to be most closely related to ordinary chondrites (e.g., Trinquier et al., 2007, 2009; Javoy
et al., 2010; Warren, 2011; Dauphas, 2017).

    Here we compare the compositional models of Earth (e.g., McDonough and Sun, 1995; Palme
and O’Neill, 2014) and Mars (e.g., Wänke and Dreibus, 1994; Taylor, 2013; Yoshizaki and Mc-
Donough, 2020) to clarify their characteristics and formation processes. We explore the nature of
the building blocks of the terrestrial planets, based on their compositional similarities and differ-
ences with chondritic meteorites. These comparisons provide a basis for insights into the conditions
for habitable planet formation and evolution.

                                                  2
Earth and Mars - distinct inner solar system products
2     Comparisons
2.1   Bulk planet
The geochemical classification of elements consists of four element groups: lithophile (rock-loving),
siderophile (metal-loving), chalcophile (sulfide-loving), and atmophile elements. Lithophile ele-
ments (e.g., Si, Mg, Ca, Al, Ti, Na, K, rare earth elements) are preferentially incorporated into
oxide phases, and thus concentrated in a silicate shell in a differentiated planetary body. There-
fore, their abundances in the mantle can be converted to the bulk composition if we also know a
mass fraction of a metallic core in a planet. In addition, elements are also classified based on their
volatilities in a gas of solar composition at 10 Pa (Lodders, 2003).

    Historically, there have been three families of compositional models for the bulk silicate Earth,
in terms of the proportions of the refractory elements relative to Mg and Si, which, together with
Fe and O, make up ∼90% of rocky planets. These models have low (Javoy, 1995; Javoy et al.,
2010; Warren, 2008; Caro and Bourdon, 2010), medium (Ringwood, 1975; Jagoutz et al., 1979;
Wänke, 1981; Hart and Zindler, 1986; Allègre et al., 1995; McDonough and Sun, 1995; Palme and
O’Neill, 2014), and high (Wasserburg et al., 1963; Turcotte et al., 2001; Turcotte and Schubert, 2014)
proportion of refractory elements (McDonough, 2016 and references therein). Here we estimate a
primitive (i.e., the least melt-depleted) mantle composition based on compositional trends from
basalts and mantle rocks and propose a model with medium refractory element abundance (e.g.,
McDonough and Sun, 1995; Palme and O’Neill, 2014; see also Section A.1). As summarized in
Palme and O’Neill (2014), these geochemical models of Earth are basically similar to each other,
and differences between the models do not affect comparison of Earth and Mars discussed in the
current paper.

    A similar geochemical approach has been applied in constraining compositional models of Mars
(Wänke and Dreibus, 1994; Taylor, 2013; Yoshizaki and McDonough, 2020). Here we adopt the
compositional model by Yoshizaki and McDonough (2020), which is compositionally and miner-
alogically similar to those of Wänke and Dreibus (1994) and Taylor (2013) (Fig. A.1; see also Bertka
and Fei, 1997; Khan et al., 2018; Smrekar et al., 2019 for Martian mantle mineralogy models), but
differ in model development. For example, Wänke and Dreibus (1994) assumed that Mn and more
refractory elements (including Fe, Mg and Si) are in chondritic relative abundances in bulk Mars,
which is denied by Yoshizaki and McDonough (2020). The details of previous Mars models are
summarized by Taylor (2013) and Yoshizaki and McDonough (2020).

    By establishing the planet’s budget of the 36 refractory elements, recognizing that ratios of
refractory elements (e.g., Ca/Al, Th/U) and Fe/Ni are constant (±15% or better) in the chondritic
building blocks (e.g., Wasson and Kallemeyn, 1988; Alexander, 2019a,b), and either knowing the
relative mass of the planet’s core or the mantle’s Mg# (atomic ratio of Mg/(Mg + Fe)), limit the range
of acceptable Mg/Si (i.e., olivine/pyroxene) values for planetary compositional models (Figs. 1
and 2A). The Fe/Mg vs Fe/Si correlation seen in both the NC and CC chondrites is suggestive of

                                                  3
Earth and Mars - distinct inner solar system products
metal-silicate gradients in the solar system (Figs. 1 and 2B); interestingly, Earth and Mars also
follow this correlation.

    The approach used here finds that the bulk Earth contains 1.85 × CI abundances for the refrac-
tory elements (e.g., McDonough and Sun, 1995; Palme and O’Neill, 2014), it has an olivine/pyroxene
proportion equivalent to that of pyrolite (Ringwood, 1966), and is depleted in volatile elements (i.e.,
sub-CI Rb/Sr, K/U, and S; Fig. 3A; Gast, 1960; Wasserburg et al., 1964). Mars also contains ∼ 1.9
× CI abundances for the refractory elements and is depleted in volatile elements, but less so than
Earth (Fig. 3B; Yoshizaki and McDonough, 2020). The net atomic oxygen/(metallic Fe) values of
Earth and Mars (3.7 and 8.7, respectively; Table 1) provide a measure of its average oxidation state.
These compositional models for Earth and Mars provide a time-integrated perspective for materials
available for accretion at 1 and 1.5 AU, respectively. This spatial sampling contrasts with the tem-
poral comparison, as the mean timescales for Mars’ formation (τMars accretion ∼ 2 Myr; Dauphas and

Pourmand, 2011) differs from that for Earth (τEarth
                                                accretion & 30 Myr; Kleine et al., 2009).

2.2     Bulk silicate planet
2.2.1    Geochemistry and mineralogy
Chondritic ratios of refractory lithophile elements in the bulk silicate Earth and Mars (BSE and
BSM, respectively) are further confirmed and constrained by Lu-Hf, Sm-Nd and La-Ce isotopic
systematics of these planetary materials (e.g., Bouvier et al., 2008; Burkhardt et al., 2016; Willig
et al., 2020). The moderately volatile, lithophile elements (lithophile MVE; e.g., alkali metals)
are depleted in Earth and Mars when compared to chondritic meteorites, both planets showing
comparable correlations with their condensation temperatures (Fig. 3). Earth, with K/U of 14000
and Rb/Sr of 0.032 (i.e., MVE/refractory ratios; cf., 68000 and 0.30 in CI chondrites, respectively),
shows such a strong depletion trend in these elements (Fig. 3A). Mars, with a K/U of 20000 and
Rb/Sr of 0.068, also shows this trend (Fig. 3B), albeit less depleted than the Earth’s.

    Siderophile and chalcophile elements are depleted in the BSE and BSM compared to lithophile
elements with similar condensation temperatures, due to their incorporation into the metallic cores
(Fig. 4). The abundances of siderophile elements in the BSE and BSM do not show a sub-parallel
trend with that defined by lithophiles (Fig. 4). The degree of siderophile and chalcophile element
depletion reflects a combination of planetary building block compositions and element fractionation
processes during planetary accretion and differentiation (McDonough, 2014).

    The BSE and BSM both have high concentrations and chondritic relative proportions of highly
siderophile elements (HSE: Re, Os, Ir, Pt, Ru, Rh, Pd and Au; Fig. 4), which is at odds with any
combination of high-pressure and temperature conditions for the partitioning of these elements into
a core forming melt (Brandon et al., 2012; Walker et al., 2015; Day et al., 2016; Tait and Day,
2018). For example, core-mantle equilibration model in Martian interior (Righter et al., 2015)

                                                  4
Earth and Mars - distinct inner solar system products
predicts sub-CI Re/Os and super-CI Ir/Os and Ir/Ru values in silicates, which are inconsistent with
the chondritic HSE pattern in the BSM (Brandon et al., 2012; Tait and Day, 2018). In addition,
Righter et al. (2015) assumes a high S concentration (>10 wt%) in the Martian core, which is
much higher than a value predicted by the Martian volatility trend of lithophile elements (≤7 wt%;
Yoshizaki and McDonough, 2020). These shared features of the HSE abundances in Earth and Mars
are considered as evidence for late accretion of chondritic materials during the final stage (e.g., after
≥98% accretion) of planetary growth in the inner solar system (Kimura et al., 1974; Brandon et al.,
2012; Day et al., 2016). Delivery of these HSE was also accompanied by an addition of volatile
gases and fluids (H, C, N and O; Albarede, 2009; Alexander et al., 2012; Marty, 2012), although the
late-accreted material to Mars might have been volatile-depleted (Wang and Becker, 2017; Righter
et al., 2019).

    On the other hand, the BSE is characterized superchondritic Ru/Ir and possibly Pd/Ir, which
cannot be easily accommodated by the late addition of volatile-rich materials (Becker et al., 2006).
Experimental studies showed that these high ratios in the BSE might reflect P , T , and composition-
dependent changes in partitioning behaviors of Ru and Pd, and proposed no need for the late addition
of these elements (Righter et al., 2008; Wheeler et al., 2011; Laurenz et al., 2016; Righter et al.,
2018). Alternatively, these inconsistencies between the BSE and chondrites show an unrepresenta-
tive sampling of asteroidal materials in our meteorite collections (Walker et al., 2015).

     The average, time-integrated, planet-scale redox condition (i.e., metal-silicate equilibrium) is
recorded in the Mg# of its mantle, with the BSE and BSM having an Mg# of ∼0.89 and 0.75–0.8,
respectively, and its core mass fraction (i.e., Mars ∼20% and Earth ∼ 32%; Section 2.3). These
attributes document Mars being more oxidized than Earth. In addition, Mars lacks depletion in
redox-sensitive, nominally lithophile elements (i.e., V, Cr, and Mn), and its mantle has a lower Hf/W
value (∼3.5; Dauphas and Pourmand, 2011; Yoshizaki and McDonough, 2020), which contrast with
those for Earth (Fig. 4). The more oxidized conditions for the Martian mantle are also indicated by
mineralogy, trace element compositions and valence states of Fe and Eu in Martian meteorites. The
Martian mantle, however, as recorded in Martian meteorites, shows a heterogeneous redox state
(e.g., Herd et al., 2001, 2002; Wadhwa, 2001, 2008; Goodrich et al., 2003; McCanta et al., 2009;
Shearer et al., 2011; Righter et al., 2016; Herd, 2019). Oxygen might be one of the light elements in
the Earth’s core (Ohtani and Ringwood, 1984), but its limited solubility in iron liquids at high P -T
conditions (O’Neill et al., 1998) indicates small effects of the core formation in the BSE’s Mg#, and
supports the more reduced state of the BSE as compared to the BSM.

    Earth and Mars have similar "upper mantle" mineralogies (Fig. 5). Both contain olivine, ortho-
and clino-pyroxenes, and garnet, with the Earth’s mantle being richer in olivine than the Martian
mantle (∼60% and ∼45% in modal proportion, respectively). Differences in composition (Table 1)
and interior thermal gradient (Breuer and Spohn, 2003; Katsura et al., 2010) of Earth and Mars
results in different depths of the olivine-wadsleyite and wadsleyite-ringwoodite phase transitions.
For Earth, these transitions are at 410 km and 520 km depth, respectively, and in Mars, they occur at

                                                   5
1000 km and 1300 km depth, respectively. These phase transitions are sharper in the Earth’s mantle
than the Martian mantle, because of the former’s higher Mg#, greater garnet abundance, and hotter
temperature (Frost, 2003; Filiberto and Dasgupta, 2011; Putirka, 2016). The Martian mantle might
not have a layer of bridgmanite, which is a dominant phase in the Earth’s lower mantle (Fig. 6).

2.2.2   Heat producing elements (HPE)
The rocky planets are powered by both primordial accretion energy (and lesser amounts from core
formation) and radiogenic heat by radioactive decays of heat producing elements (HPE: K, Th,
U). All three of the HPE are incompatible, lithophile elements that have been excluded from their
cores and concentrated into planetary crusts (e.g., Corgne et al., 2007; Blanchard et al., 2017; Wip-
perfurth et al., 2018), which is beneficial for carrying out remote gamma-ray surveys of planetary
surfaces (e.g., Peplowski et al., 2011; Surkov et al., 1987; Prettyman et al., 2015). The Martian
crust is thicker (estimated to be 30–60 km thick with an average value of ∼50 km; Zuber et al.,
2000; McGovern et al., 2002; Wieczorek and Zuber, 2004; Humayun et al., 2013), less enriched
in incompatible elements, and contains ∼50% of the HPE budget of Mars, whereas Earth has a
thinner, more chemically evolved crust, which contains only ∼35% of the planet’s HPE budget (Ta-
ble 1 and Section A.2; Rudnick and Gao, 2014). Mars’ internal heating Rayleigh number of >105
(Section A.3) is consistent with a convecting Martian mantle and a dynamically stabilized Tharsis
bulge (McKenzie et al., 2002; Kiefer, 2003).

    With nearly a factor of 2 greater surface to volume ratio for Mars than Earth, the former cooled
much faster (Filiberto and Dasgupta, 2011; Baratoux et al., 2011, 2013; Filiberto, 2017; Breuer and
Moore, 2015; Putirka, 2016). Mars’ surface heat flux is 19 ± 1 mW/m2 (Parro et al., 2017) (cf.,
Earth’s average is 90 mW/m2 ; Jaupart et al., 2015), equivalent to a global heat flux of 2.75 ± 0.15
TW, of which 2.5 TW comes from radioactive decay (Table 1 and Section A.2). Mars’ planetary
Urey ratio (U r: total radioactive heat production relative to total surface heat loss) is 0.92 ± 0.05,
which is higher than that estimated for Earth (0.43) and its mantle U rMars (planetary U rMars minus
the crustal fraction) is 0.8, more than double of that of the Earth’s mantle (0.33). The high U rMars
demonstrates a minor contribution from secular cooling to the total Martian heat flux. Importantly,
the modeled Martian surface heat flow varies between 14 and 25 mW/m2 (Parro et al., 2017). A
direct measurement the Martian surface heat flux by NASA’s Interior Exploration using Seismic
Investigations, Geodesy and Heat Transport (InSight) mission (Banerdt et al., 2020) will provide
important constraints on its heat production and thermal history.

2.3     Core
The mantles of both Earth and Mars are depleted in siderophile and chalcophile elements due to
core extraction (Fig. 4). Their core compositions are modeled using the following constraints: the
planet’s chondritic Fe/Ni value, temperature, mass, density, MOI, and for Earth, seismic profile (see

                                                  6
McDonough, 2014; Yoshizaki and McDonough, 2020). Also, based on iron meteorite mineralogy
(Scott and Krot, 2014) and cosmic abundance of elements (Lodders, 2020), planetary cores might
contain sulfides, carbides and phosphides, and possibly other phases (e.g., silicides if formed under
highly reduced conditions).

     The Earth’s core (32% by mass) has an outer liquid layer and a solid inner core (∼5 wt%) (Fig. 6).
Estimates of the size of the Martian core ranges from 18–25 wt% and 1500–1900 km, respectively
(Section A.1.3). Data from multiply orbiting satellites provide a precise (±2%) measure of Mars’
Love number solution (k2 ) and document a possible existence of a partially molten core (Yoder et al.,
2003; Genova et al., 2016; Konopliv et al., 2016). This observation indicates significant amounts
of light elements in the core that lowers its solidus temperature.

    Sulfur is a prime candidate for light elements in a planetary core. The volatile element depletion
trends for both Earth and Mars (Fig. 4) constrain the core’s S content of ∼1.8% in Earth and ∼7%
in Mars (Yoshizaki and McDonough, 2020; see also Sections A.1.2 and A.1.3). Since the estimated
S contents of the Earth’s and Mars’ core are not high enough to explain their respected core density
deficits, additional light elements in these cores are needed (e.g., Si, O, H and C). The additional
light elements are likely to be different in these cores because metal-silicate partitioning behavior
of elements depends on conditions of core formation (e.g., timing, P , T , fO2 ) and compositions of
coexisting silicate and metallic phases.

    As described above, a planet’s volatility trend sets expectations for the absolute amount of
siderophile and chalcophile elements in the planet and defines the proportion of these elements
that were partitioned into the core. The degree of depletion below the volatility trend defines an
element’s empirically established metal/silicate partition coefficient (McDonough, 2014; Yoshizaki
and McDonough, 2020). The Martian core is enriched in siderophile and chalcophile MVE (e.g., P,
Ge, S) as compared to the Earth’s, and there appears to be little to no Ga in the cores of terrestrial
planets, which contrasts with that seen in iron meteorite compositions (Figs. 7 and 8). Compared to
the Earth’s core, the Martian core is smaller, it might contain larger amounts of H and O, and it has
a lower Fe/Ni value, with its mantle being enriched in Fe (Table 1). Overall, these factors contribute
to Mars’ lower uncompressed density as compared to the Earth’s (3750 vs 4060 kg/m3 ; Table 1).
There is no unique model for core compositions of Earth and Mars that satisfies the geodetic and
geochemical constraints (McDonough, 2014; Yoshizaki and McDonough, 2020), and further ex-
perimental and theoretical efforts are needed to constrain the light element budget in the planetary
cores.

                                                  7
3      Discussion
3.1     Origin of volatile depletion in terrestrial planets
From Mercury to Mars and beyond to Vesta, the second biggest asteroid, their surface K/Th values
are significantly lower relative to chondrites (Surkov et al., 1987; Taylor et al., 2006; Peplowski et al.,
2011; Prettyman et al., 2015). The origins of this depletion and that of other MVE remain elusive,
with possible explanations including post-nebular volatile loss due to internal or impact-induced
heating (e.g., O’Neill and Palme, 2008; Norris and Wood, 2017) and incomplete condensation of
nebular gas (e.g., Palme and O’Neill, 2014; Siebert et al., 2018). Since planetary K/Th ratios and rel-
ative amounts of more refractory olivine (Mg2 SiO4 ) and less refractory pyroxene (MgSiO3 ) do not
correlate with their heliocentric distances in both the solar and extra-solar systems (e.g., Fig. 2A; van
Boekel et al., 2004; Kessler-Silacci et al., 2006; Bouwman et al., 2008, 2010; Sargent et al., 2008),
these volatile depletions do not solely reflect an outward temperature decrease in the protoplanetary
disk.

    The absence of heavy isotope enrichment in rocks from Earth and Mars for multiple isotope sys-
tems (e.g., K, Zn, Rb, Fe, Cd; Humayun and Clayton, 1995; Nebel et al., 2011; Paniello et al., 2012a;
Sossi et al., 2016, 2018; Pringle and Moynier, 2017; Wombacher et al., 2008) indicates negligible
post-accretionary evaporative loss of the MVE. In contrast, isotopically heavier siderophile or chal-
cophile compositions (e.g., Fe, Ga, Sn) of terrestrial and lunar mantles might reflect (1) core-mantle
differentiation, (2) evaporation of the giant impactor’s core and addition of the metal/sulfur-loving
elements to the terrestrial mantle during the giant impact (see Section 3.3), or (3) planetary surface
processes (e.g., Poitrasson et al., 2004; Kato and Moynier, 2017; Creech and Moynier, 2019). The
53
  Mn- 53Cr and 87Rb- 87Sr isotope systematics of meteorites and terrestrial samples are consistent
with a volatility-dependent, gas-solid fractionation during first few Myr of the solar system (e.g.,
Shukolyukov and Lugmair, 2006; Trinquier et al., 2008; Hans et al., 2013; Moynier et al., 2012),
which might be prior to the dissipation of the nebular gas (∼ t0 + 5 Myr; Wang et al., 2017).

     There is negligible evidence for post-accretionary MVE losses in the chemical compositions of
these planets. The bulk Earth and Mars, together with chondrites, show negligible evidence of evap-
orative losses in ratios of Na, Mn and Ti (Fig. 9 and Section A.1.1; O’Neill and Palme, 2008; Siebert
et al., 2018), which have distinct relative volatilities during condensation and evaporation. Values
of Mn/Na and Na/Ti in the bulk Earth, Mars and chondrites are consistent with an incomplete nebu-
lar condensation, in which earlier condensates are removed from a nebular gas before a completion
of more volatile species (e.g., Larimer, 1967; Larimer and Anders, 1967) (cf. O’Neill and Palme,
2008). The plot of Mn/Na versus Na/Ti indicates that a formation of precursors of Earth and Mars
at higher temperatures than their isotopically linked counterparts (enstatite and ordinary chondrites,
respectively), which is also suggested by the former’s higher Mg/Si (i.e., olivine/pyroxene), Al/Si,
and Al/Mg ratios (Kerridge, 1979; Larimer, 1979; Dauphas et al., 2015; Morbidelli et al., 2020).

      Collectively, the compositions of rocky planets, as compared to their chondritic relatives, likely

                                                    8
reflects volatility-dependent chemical fractionation in the protoplanetary disk, rather than the post-
accretionary losses of MVE. In contrast, small differentiated asteroids (e.g., parent bodies of Eucrite
and Angrite) show clear evidence for preferential loss of MVE (Fig. 9; O’Neill and Palme, 2008)
and heavy isotope enrichment (e.g., Paniello et al., 2012b; Pringle et al., 2017; Tian et al., 2019),
which are the hallmarks of evaporative losses during or after accretion. Volatile elements might
have escaped from these small bodies during and/or after their accretion, due to their weak gravity
field.

3.2   The building blocks of the terrestrial planets
The planetary building blocks appeared to be made up of high-temperature materials, dominantly
chondrules, which are silicate droplets formed by transient heating events within first few Myr of the
solar system evolution (Connelly et al., 2012; Bollard et al., 2017), and are an essential component
of chondrites (e.g., Scott and Krot, 2014; Russell et al., 2018). MVE composition of chondrules
provides a record of incomplete nebular condensation rather than evaporation processes (Humayun
and Clayton, 1995; Alexander et al., 2000, 2008; Galy et al., 2000; Pringle et al., 2017). The re-
fractory element enrichment and MVE depletion in the BSE and BSM are comparable to those ob-
served for chondrules from carbonaceous chondrites (Fig. 10; Hewins and Herzberg, 1996; Mahan
et al., 2018). Significantly, mass-dependent Ca (Huang et al., 2010; Magna et al., 2015; Amsellem
et al., 2017; Simon et al., 2017; Bermingham et al., 2018) and Mg (Bizzarro et al., 2004; Young
and Galy, 2004; Wiechert and Halliday, 2007; Bouvier et al., 2013; Olsen et al., 2016; Hin et al.,
2017) isotopic compositions of chondrules, BSE and BSM also suggest that they inherited fraction-
ated isotopic signatures from similar precursor materials, which have experienced high-temperature
gas-solid fractionation processes.

     Recent theoretical models of planetary growth favor formation of terrestrial planets via accretion
of chondrule-sized pebbles (e.g., Johansen et al., 2015a,b; Levison et al., 2015). These models
predict a rapid accretion of Mars-sized bodies under the presence of a nebular disk. Thus, Mars,
with a τMars
         accretion of ∼2 Myr (Dauphas and Pourmand, 2011; Kruijer et al., 2017b; Bouvier et al.,

2018) (cf. Marchi et al., 2020), formed within the nebular disk lifetime (∼5 Myr; Wang et al., 2017),
                                      accretion ≥30 Myr (e.g., Kleine et al., 2009), which documents
whereas Earth is suggested to have τEarth
its accretion stretched beyond the lifetime of the nebular disk.

    Differences in planetary MVE abundances (Fig. A.3) likely reflect aspects of their accretion
histories and sizes. Chondrites, the least MVE-depleted bodies’ materials, contain fine-grained,
volatile-rich matrix that accounts for the largest fraction of the inventory of volatiles (e.g., Alexan-
der, 2005, 2019a,b; Bland et al., 2005; Zanda et al., 2018). Although a relationship between chon-
drules and matrix remains poorly understood, their coexistence in the protoplanetary disk before
planetesimal accretion is accepted (e.g., Hezel et al., 2018a; Zanda et al., 2018). Some carbonaceous
chondritic asteroids accreted more chondrules and less matrix, resulting in more MVE-depleted
compositions (Figs. 11B and A.4). According to the pebble accretion model (Johansen et al.,

                                                   9
2015a,b; Levison et al., 2015), planetesimal growth starts off by preferentially accreting the smallest
particles, and as the body grows, it prefers to accrete larger and larger size particles. This mechanism
would lead to a growing planetesimal having a larger chondrule/matrix ratio, and becoming more
MVE-depleted as its mass increases. Thus, small chondritic asteroids likely co-accreted chondrules
and matrix, whereas Mars, with its intermediate volatile depletion, size and accretion timescale,
likely accreted a greater fraction of chondrules to matrix (Fig. 11). Finally, the Earth’s prolonged
accretion history was dominated by chondrule accretion, resulting in its significant MVE-depleted
composition. The MVE depletion scales with size of the chondritic parent body and planet (Fig. 11),
implying an accretion driven process from an undepleted nebula (i.e., CI (solar) composition). Ex-
ceptions are small differentiated bodies (e.g., the Moon and Vesta) that experienced post-nebular
volatile loss due to internal or impact-induced heating (Section 3.1). Thus, the chondrule-driven
planetary growth plays a critical role in establishing the planetary MVE-depleted characteristics.

    In contrast, enstatite and ordinary chondrites, whose isotopic composition and redox state are
most related to Earth and Mars, respectively (e.g., Trinquier et al., 2007, 2009; Warren, 2011;
Dauphas, 2017; Fig. 2B), do not show a clear refractory enrichment nor MVE depletion (Figs. 2A,
11 and A.4). In addition, chondrules from these NC chondrites are less depleted in MVE than the
BSE and BSM (Fig. A.5). Thus, the chondrule-driven MVE-depletion scenario discussed above
cannot be applied to the NC chondrites we have in our collections. These observations indicate that
the NC chondrites represent refractory-poor, volatile-rich counterparts of the inner rocky planets.
Morbidelli et al. (2020) showed that the low Mg/Si and Al/Si solids, which are comparable to those
of NC chondrites (Fig. 2A), condense after removal of early-formed, high-temperature condensates
from the system.

    Differences in the timing of planetary accretion might also be important in establishing their rel-
ative abundances of MVE. Chronology of meteorites combined with thermal modeling of asteroids
(Sugiura and Fujiya, 2014; Kruijer et al., 2017a; Zhu et al., 2020) indicates that the undifferentiated
NC planetesimals accreted 1–2 Myr after formation of differentiated NC bodies (i.e., iron meteorite
parent bodies and terrestrial planets). The variation in the absolute ages of chondrules (Connelly
et al., 2012; Bollard et al., 2017), occurrence of relict grains (Jones, 1994; Weisberg et al., 2011;
Tenner et al., 2018) and igneous rims (Krot and Wasson, 1995) in chondrules, and the presence
of compound chondrules (Wasson et al., 1995) indicate that chondrule formation was a repeated
event. Thus, the longer chondrules remained in the accretionary disk, the more opportunities it has
being recycled by later events. Ordinary chondrite chondrules records an admixing of MVE-rich
CC-like materials to MVE-poor chondrule precursors into the NC reservoir (Mahan et al., 2018;
Schiller et al., 2018; Bollard et al., 2019). Additions of MVE-rich materials and repeated chondrule
recycling produce younger chondrules with higher MVE abundance (Mahan et al., 2018). Thus,
terrestrial planets, which are dominated by earlier materials that experienced fewer recycle events
and MVE addition, are more depleted in MVE as compared to the NC chondrites, which accreted
the younger MVE-enriched chondrules.

    In contrast, chondrule formation scenarios predict less frequent chondrule formation/recycle

                                                  10
events in the outer solar system (e.g., Morris et al., 2012; Johnson et al., 2015; Sanders and Scott,
2018; Pilipp et al., 1998). This prediction is consistent with chemical and isotopic signatures of CC
chondrules (e.g., Hewins and Zanda, 2012; Tenner et al., 2018; Mahan et al., 2018). Consequently,
the CC chondrules could have preserved their MVE depletion until the accretion of CC chondritic
asteroids, and thus are chemically comparable to the early NC materials which formed the inner
rocky planets.

    This chondrule-rich accretion model for the terrestrial planets reveals the limited ability to reach
greater levels of enrichment in refractory elements (Figs. 2A and 10B), which is a concern for
planetary compositional models promoting high refractory element abundance (Section 2.1). To
reach higher levels of refractory element enrichment beyond that seen in CV chondrites and their
chondrules, larger amounts of refractory inclusions are needed, which is inconsistent with chondritic
REE ratios of the BSE (Stracke et al., 2012; Dauphas and Pourmand, 2015).

3.3     Formation models of the terrestrial planets
3.3.1    Previous models of the Earth’s accretion and Moon formation
Depletion of moderately to highly siderophile and chalcophile elements in the BSE is consistent
with Earth’s initial accretion from highly reduced, volatile-depleted materials, that was later oxi-
dized by volatile-rich additions (e.g., Wänke, 1981; Wänke et al., 1984; Wänke and Dreibus, 1988;
O’Neill, 1991b; Rubie et al., 2011, 2015). Support for this temporal evolution in volatile accretion
is found in multiple radiogenic isotope systems (e.g., Rb-Sr, U-Pb, Ag-Pd, I-Pu-Xe; Halliday and
Porcelli, 2001; Albarede, 2009; Schönbächler et al., 2010; Mukhopadhyay, 2012; Ballhaus et al.,
2013; Maltese and Mezger, 2020), metal-silicate partitioning behaviors of elements such as W, Mo,
S and C (e.g., Wade et al., 2012; Li et al., 2016; Suer et al., 2017; Tsuno et al., 2018; Ballhaus et al.,
2017), and N-body simulations of planetary formation (e.g., Raymond et al., 2006; Morbidelli et al.,
2012). However, the degree of volatile element depletion in the pre-impact proto-Earth remains un-
constrained.

    The two-component accretion models for the growth of Earth envisage mixing of highly re-
duced, volatile-depleted materials with oxidized, volatile-rich, "CI-chondritic" materials in some
proportion, not well defined, but generally conceived their mass ratios to be in the range of 60:40
to 90:10 (Wänke et al., 1984; Wänke and Dreibus, 1988; O’Neill, 1991b). Some mass estimates
of the late oxidized addition are comparable to a Mars-sized, Moon-forming impactor (e.g., Canup
and Asphaug, 2001; Canup, 2004, 2008). The impactor’s composition has often been assumed to
be CI-chondritic (e.g., O’Neill, 1991a,b; Maltese and Mezger, 2020).

    However, the CI-chondritic impactor model has been multiply challenged. The common MVE
depletion in rocky differentiated bodies in the solar system (e.g., sub-solar K/Th ratios; Section 3.1)
might indicate a similar MVE depletion in the Mars-sized impactor. In addition, the CI-like im-
pactor model requires a significant volatile depletion in the proto-Earth, perhaps at levels seen in

                                                   11
angrites and calcium-aluminum-rich inclusions (CAIs), which show heavy Mg- and Si-isotope en-
richment and significant depletion of moderately volatile elements (e.g., Grossman et al., 2000,
2008; Pringle et al., 2014) due to significant evaporative losses of the major and more volatile ele-
ments by impact-induced or transient nebular heating events (e.g., Stolper and Paque, 1986; Richter
et al., 2002; Pringle et al., 2014; Yoshizaki et al., 2019; Young et al., 2019). However, such isotopic
signatures are not recognized for Earth (Section 3.1). Furthermore, models predicting composi-
tional zoning in the protoplanetary disk have a region of CI-like material accreting at ≥15 AU
(Desch et al., 2018), where it is predicted that the disk mass is too low to form a Mars-sized body.

    Compositional similarities of Earth and Moon in multiple isotope systems set strict constraints
on the nature of both the impactor and proto-Earth, that is, they are derived from a similar NC iso-
topic reservoir that appears to be restricted to inner solar system sources (e.g., Wiechert et al., 2001;
Trinquier et al., 2009; Warren, 2011; Zhang et al., 2012; Greenwood et al., 2018; Dauphas, 2017;
Kruijer et al., 2017a). This observation is consistent with dynamical simulations which predict low
probabilities of a Moon-forming impactor originating from >10 AU (e.g., Jackson et al., 2018).
Thus, these constraints exclude the isotopically distinct CI chondrite as an impactor candidate. In
contrast, recent Mo and O isotopic data from a broad range of NC and CC materials challenge this
exclusion and find support for a CC-like impactor and/or vigorous mixing of proto-Earth and im-
pactor, requesting further constraints to reveal the origin of the Earth-Moon system (Young et al.,
2016; Budde et al., 2019; Cano et al., 2020). Furthermore, the recently-proposed synestia model
(Lock et al., 2018) has the Moon forming in a vapor cloud surrounding Earth. This vapor cloud was
produced by a large impact, resulting in a well-mixed, chemically equilibrated proto-earth and the
vapor cloud.

3.3.2   Mars-like Moon-forming giant impactor model
Here we propose a model for the origin of the Moon. Our model is a modification of the original
Wänke et al. (1984)’s model and later modified by O’Neill (1991a,b). In our version of this model,
we envision the Earth’s formational history in 4 steps and use a differentiated Mars-like composition
(Yoshizaki and McDonough, 2020), instead of CI composition, for the giant impactor. Figs. 12
and 13 show the details of our model: Fig. 12 specifically highlights the differences between our
model and those previously presented. Our model starts with
   1. the accretion of the proto-Earth (reduced and volatile-depleted) accompanied by continuous
      core-mantle differentiation (∼90% of Earth’s mass),

   2. followed by a late-stage Moon-forming giant impact event (30–100 Myr after t0 ; Kleine et al.,
      2009) that adds the final ∼10% mass (oxidized and volatile-enriched) to Earth and forms a
      protolunar accretion disk,

   3. subsequently, the mantle loses a Fe-Ni (±O) sulfide liquid (sulfide matte; O’Neill, 1991b) to
      the core (∼0.5% BSE mass),

                                                   12
4. and finally, the BSE receives the addition of (∼0.5% of the BSE mass) a late accretion com-
      ponent that brings the highly siderophile and chalcophile elements in chondritic proportions
      and highly volatile gases and fluids (see Section A.5 for details of the modeling).

    In this scenario, the bulk proto-Earth already contains ≥80% of the present-day Earth’s budgets
of most of the MVE (e.g., K/Th ∼ 3200; Rb/Sr ∼ 0.026; Table A.2), since a Mars-like impactor
contributes only a limited amount of additional MVE (Fig. 13). The addition of MVE by the im-
pactor leads to sulfur saturation in the magma ocean, generating an Fe-Ni (± O) sulfide liquid (post-
impact sulfide matte; O’Neill, 1991b; Rubie et al., 2016). The sulfide matte precipitates through
a crystallizing mantle into the core due to its high immiscibility, low wetting angle, and high den-
sity (Gaetani and Grove, 1999; Rose and Brenan, 2001). Assuming that the sulfide matte, with a
present-day mantle sulfide-like composition (Lorand and Conquéré, 1983), extracted all S from the
post-impact Earth’s mantle; its mass is estimated to be ∼0.5 wt% of the present-day Earth’s silicate
mantle (Section A.5).

    The siderophile and chalcophile element abundances of the BSE are reproduced when 10–15%
of a planetary mass is added by the impactor and most of the impactor’s core equilibrates with
the proto-Earth’s mantle (Fig. 13 and Section A.5). In this scenario, the bulk proto-Earth contains
≥80% of the present-day Earth’s budget of most of the MVE (e.g., K/Th ∼ 3200; Rb/Sr ∼ 0.026;
Table A.2), with the Mars-like impactor contributing a limited amount of MVE (Fig. 13). We en-
vision the proto-Earth as having a MVE abundance comparable to that of chondrules (Figs. 10
and 11B). Thus, the proto-Earth’s composition is consistent with the chondrule-rich accretion sce-
nario (Section 3.2) and requires no need for a post-accretionary loss of MVE from the proto-Earth
before the Moon-forming event (Section 3.1).

    The mass fraction of the impactor contributing to the lunar composition provides a critical con-
straint on the lunar formation models. The lunar mantle is depleted in nominally lithophile V, Cr
and Mn (Dreibus and Wänke, 1979). Such lunar mantle depletion can be achieved by high-T or
S-rich conditions during lunar core formation, but these conditions seem unlikely (Steenstra et al.,
2016). If a Mars-sized impactor with no V, Cr or Mn-depletion in its mantle contributed >70%
of Moon, as predicted by the canonical giant impact models (Canup and Asphaug, 2001; Canup,
2004, 2008), ≥40% of evaporative losses of Mn and Cr are needed to produce their depletion in
the lunar mantle (Fig. 14). Such significant evaporation of these elements is inconsistent with their
least volatile nature among MVE (Gellissen et al., 2019; Sossi et al., 2019). Thus, the proto-Earth’s
V, Cr and Mn-depleted mantle should be a primary source of the Moon-forming materials, as sup-
ported by their isotopic similarities (e.g., Wiechert et al., 2001; Warren, 2011; Zhang et al., 2012;
Greenwood et al., 2018). A recent geochemical model of Earth’s Moon formation and differentia-
tion prefers a present-day BSE-like MVE composition for the bulk proto-Moon (i.e., lunar source
materials before gas-melt fractionation; Righter, 2019), which is consistent with the proto-Earth’s
MVE abundance and its large contribution to the lunar source materials (Fig. 14 and Table A.2).
The proto-Earth origin of the Moon is also consistent with recent particle hydrodynamic collision

                                                 13
simulations (Hosono et al., 2019) which showed that a terrestrial magma ocean might be a source
of lunar building blocks.

3.3.3    Mars’ accretion and core formation
Our compositional estimate of Mars requires a much simpler formation history as compared to
Earth’s multi-stage formation scenario. The lack of depletion in nominally lithophile elements in
the BSM (V, Cr, and Mn; Fig. 4) is consistent with its accretion under uniformly oxidized and low-T
conditions (e.g., Wood et al., 2009). Mars might have experienced a large impact(s) that produced
its hemispheric crustal dichotomy and possibly Martian moons (e.g., Marinova et al., 2008; Canup
and Salmon, 2018), but estimates of a putative small impactor (∼0.1% of the Mars’ mass; Canup
and Salmon, 2018) would result in negligible changes in its bulk composition. The 0.5–1% mass
addition of chondritic late accretion material provides highly siderophile and chalcophile elements
to the Martian mantle (Tait and Day, 2018; Yoshizaki and McDonough, 2020).

    Righter and Chabot (2011) and Yang et al. (2015) estimated Martian core formation happening
at 14 ± 3 GPa and 2100 ± 200 K based on its previous compositional models. Following the
approaches of Righter and Chabot (2011) and Yang et al. (2015), the Martian siderophile element
distribution of Yoshizaki and McDonough (2020) can be modeled by single stage P -T and redox
conditions (Fig. A.8), whereas there is no unique solution for Martian core formation (Fig. A.9).

     The simple formation history of Mars, combined with its rapid and early accretion (Dauphas and
Pourmand, 2011; Kruijer et al., 2017b; Bouvier et al., 2018; Marchi et al., 2020), is consistent with
its status as a planetary embryo. Given compositional and redox state gradients in the protoplanetary
disk, Mars might record the chemical characteristics of nebular materials in the Mars’ orbit in the
first few years of the solar system, whereas Earth might have incorporated oxidized materials from
greater heliocentric distance in its later accretion stage (Section 3.3.2; Rubie et al., 2015). The Mars’
status as a planetary embryo suggest that a Mars-sized body commonly has a Mars-like composition,
supporting the Mars-like Moon-forming impactor scenario (Section 3.3.2).

3.4     Conditions for a habitable planet formation
The early τMars
            accretion age implies that Mars underwent global-scale silicate melting and rapid core

formation due to heating from short-lived 26Al and 60Fe (Fig. A.12; Dauphas and Pourmand, 2011;
McDonough et al., 2020). The peak radiogenic heating occurs at about 1 to 5 Myr after t0 , well
within the time scale for Mars accretion. With τMars
                                                  accretion = 2 Myr, the radiogenic energy supplied

by Al is comparable to Mars’ gravitational binding energy (∼ 7 × 1029 J and ∼1030 J, respec-
    26

tively). During the first 10 Myr, radiogenic heating ( 26Al, 60Fe and 40K, in order of significance) is
comparable to the planet’s primordial energy and is a major control on its thermal evolution.

                                                   14
Mars has the attributes needed for a rocky planet to be biologically available in its early his-
tory (e.g., Ehlmann et al., 2016), and it has a higher bulk heat production than Earth’s (3.9 vs 3.3
pW/kg; Table 1). Nonetheless, it has rapidly lost much of its primordial energy (i.e., accretion and
core differentiation) due to its larger surface to volume ratio (a factor of 2) and smaller core size
(i.e., reduced bottom heating), and it is in waning stages of limited fuel resources (Parro et al., 2017;
Yoshizaki and McDonough, 2020). Basal heating of the Martian mantle by its core enhances its ther-
mal evolution, while the transfer of hydrogen from the adjacent ringwoodite to the core (Shibazaki
et al., 2009; Yoshizaki and McDonough, 2020) reduces the lifetime of the dynamo (O’Rourke and
Shim, 2019). Collectively, these processes likely contributed to dynamo termination at ∼4 Ga and
loss of the protective magnetosphere (Acuña et al., 1999; Arkani-Hamed, 2004; Lillis et al., 2008).
This magnetosphere shields the planet from atmospheric losses, enhances its surface UV radiation,
and leads to dramatic climate changes (Ehlmann et al., 2016).

    Volatile elements may play a critical role in establishing the amount of light elements in and
solidus of the core. The amount of water and other volatile species in the planet’s interior and
surface may potentially create the appropriate conditions for the initiation of plate tectonics (e.g.,
Albarede, 2009; Ehlmann et al., 2016). Likewise, the heat-producing elements and a reduced core
solidus keep the metallic core convecting and lead to the creation of a magnetic field, which shields
a planet’s surface from cosmic rays. Together, heat-producing and volatile elements regulate a
planet’s cooling history, drive crustal differentiation, and make it habitable (Ehlmann et al., 2016).
The simple formation history of Mars (Section 3.3.3; Rubie et al., 2015) emphasizes the uniqueness
of Earth, the sole habitable planet in our solar system. In turn, the Earth-Mars comparison indicates
that high-temperature nebular chemical processes and timescales of planetary accretion are essential
in making habitable planets.

    Depletion in volatile elements appears to be a common feature of the terrestrial planets (Fig. 3;
Surkov et al., 1987; Peplowski et al., 2011), and may likely be so for rocky exoplanets (Harri-
son et al., 2018, 2021). The relationship between the planetary volatile depletion, size, accretion
timescale and abundance of chondrules (Fig. 11) predict accretion timescales of Venus and Mer-
cury of 30–100 Myr and 2–10 Myr, respectively, based on their Earth- and Mars-like K/Th ratios,
respectively (Surkov et al., 1987; Peplowski et al., 2011). The predicted ages of Mercury and Venus
provide a foundation for future investigation of their thermal history and habitability.

    Venus and Earth are quite similar in their physical and chemical properties (size, bulk K/Th
ratio, Mg# in basalt and possibly MOI; e.g., Surkov et al., 1987; Dumoulin et al., 2017), but their
present-day surface conditions are distinct. Venus does not show a clear evidence for a giant impact,
and its size and K/Th ratio (∼3000; Surkov et al., 1987) are comparable to those of the proto-Earth
(K/Th ∼ 3200; Table A.2). Therefore, the present-day Venus might be comparable to the pre-
impact proto-Earth. Further observational, cosmochemical and theoretical investigations of Venus
may provide useful insights into the pre- and post-formational history of Earth, the only habitable
planet in the today’s solar system.

                                                   15
4    Conclusions
The refractory element enhancement and volatile depletion of Earth and Mars were established
by a nebular chemical fractionation. Post-accretionary losses of moderately volatile elements are
negligible. The degree of volatile element depletion correlates with the abundance of accreted
chondrules, planetary size, and their accretion timescale. The present-day bulk silicate Earth com-
position is consistent with a Moon-forming impactor having a Mars-like size and composition.
Planetary chemistry, which is related to many factors including the building block composition, the
timing, duration and sequence of accretion and its differentiation history, play an essential role in
making a planet habitable.

Acknowledgments
We thank many our colleagues who have listened to various versions of this project and given
feedback, especially Eiji Ohtani, Nick Schmerr, Beda Roskovec, Ondřej Šrámek, Sarah Stewart-
Mukhopadhyay, Kevin Righter, and Henri Samuel. We also thank Attilio Rivoldini for helpful
comments. We greatly appreciate Hugh O’Neill, Bernard Wood, and an anonymous referee for
their constructive reviews, which helped improve the manuscript. We thank the journal editor
Astrid Holzheid for editorial efforts. TY acknowledges supports from the Japanese Society for
the Promotion of Science (JP18J20708) and GP-EES and DIARE research grants. WFM gratefully
acknowledges NSF support (EAR1650365).

Author contributions
TY and WFM proposed and conceived various portions of this study and together calculated the
compositional models of planets. The manuscript was jointly written by TY and WFM and they
read and approved the final manuscript.

Competing interests
The authors declare no competing interests.

Additional information
Correspondence and requests for materials should be addressed to TY.

Data and materials availability
Materials used in this study are available within the paper or supplementary materials.

                                                 16
Table 1: Physical and chemical properties of Earth and Mars. Modeled values are in normal and reference
values in italic fonts.
             Observation                 Unit         Crust           Mantle           Core       Bulk planet         Reference value
              Earth
             Massa,b                       kg      3.12 × 1022      4.00 × 1024    1.94 × 1024   5.972 × 1024    5.97218(60) × 1024
          Mean densitya,b,c              kg/m3        2800             4400           11870          5514              5514(2)
          Moment of inertiaa                –          1%               88%            11%          0.3308           0.330690(9)
                               b,d,e,f
                                          TW           7.3              12.6            0            19.9               46(3)g
 Heat production (K, Th, U)
                                         pW/kg         232               3.1            0             3.3                 –
              Mars
              Massh                        kg      2.56 × 1022      5.01 × 1023    1.17 × 1023   6.419 × 1023      6.417(3) × 1023
           Mean densityi                 kg/m3        3010             3640           6910           3936              3935(1)
          Moment of inertiah                –          7%               89%            4%           0.3638            0.3639(1)
                                j,k,l
                                          TW           1.3               1.3            0             2.5              2.7(2)g
     Heat production (K, Th, U)
                                         pW/kg          49               2.5            0             3.9                 –

                 Mantle + crustd,k                       Cored,k                          Bulk planetd,k,m         Chondriten
 wt%            Earth       Mars           wt%       Earth     Mars            wt%       Earth      Mars        EH     L                CI
 SiO2          44.9          45.5           Fe        85.5         79.5           O       29.7       36.3       30.0       35.8     29.9
 TiO2          0.20          0.17           Ni         5.1          7.4           Fe      32.0       23.7       31.0       22.7     27.7
 Al2 O3        4.44          3.59           O           2           5.2           Mg      15.4       15.3       11.3       15.7     14.3
 MnO            0.14         0.37           S          1.8          6.6           Si      16.1       17.4       17.8       19.5     16.0
  FeO          8.06          14.7           H         0.06          0.9           Ni      1.82        1.4       1.87       1.27     1.58
 MgO            37.8         31.0           Co        0.25         0.33           Ca      1.71       1.69       0.91       1.38     1.36
 CaO            3.54        2.88            P         0.20         0.33           Al      1.59       1.56       0.86       1.29     1.26
 Na2 O         0.36          0.59           Si          4            0            S       0.59        1.2        6.2       2.3       8.0
 K2 O          0.034        0.043           Cr        0.75           0
 P2 O5         0.021         0.17                                              Total      99.0       98.5       100        100      100
  NiO          0.25         0.046          Total      99.7         100.2
 Cr2 O3        0.15          0.88                                              Mg/Si      0.96       0.88       0.63       0.81     0.89
                                                                               Al/Si      0.10       0.09       0.05       0.07     0.08
K (ppm)         280          360                                               Fe/Si       2.0        1.4        1.7        1.2      1.7
Th (ppb)         76          68
U (ppb)         20           18

 Total         100.2        99.9
 a                              b                                 c
   Chambat et al. (2010).         Wipperfurth et al. (2020).        Dziewonski and Anderson (1981).
 d                           e                            f                               g
   McDonough (2014).           Jaupart et al. (2015).       Rudnick and Gao (2014).         Global surface heat loss.
 h                             i                             j                                   k
   Konopliv et al. (2016).       Rivoldini et al. (2011).      Taylor and McLennan (2009).         Yoshizaki and McDonough (2020).
 l                         m                   n
   Parro et al. (2017).      This study.         Alexander (2019a,b) (Volatile-free, normalized to total = 100 wt%).

                                                              17
Fig. 1: (A) Ternary plot and (B) scatter ratio plot of major elements Si, Fe and Mg in bulk Earth (Mc-
Donough, 2014), bulk Mars (Yoshizaki and McDonough, 2020) and chondritic meteorites (Alexan-
der, 2019a,b). The regression line for chondrites (except for sulfide-rich EH) is also shown in the
lower panel. TL–Tagish Lake.

                                                 18
Fig. 2: (A) Weight ratios of Al/Si vs Mg/Si values for planets and chondrites. Earth and Mars
show higher RLE abundances (i.e., Al) and Mg/Si value (i.e., olivine/pyroxene ratio) compared to
enstatite and ordinary chondrites, which are respectively their isotopically identified relatives. (B)
The Urey-Craig diagram (after Urey and Craig, 1953) illustrates relative redox condition for Earth,
Mars, and chondrites. Chondrite classification (non-carbonaceous (NC; red) vs carbonaceous (CC;
blue) groups) follows Warren (2011) and Kruijer et al. (2017a). TL–Tagish Lake. Data sources are
as in Fig. 1.

                                                 19
Fig. 3: Lithophile element abundances in the bulk silicate Earth and Mars normalized to CI abun-
dances and refractory lithophile element Al. The scale of right ordinate shows the CI-normalized
abundances. The values are plotted against 50% condensation temperature (K) of elements at 10 Pa
(open; Wood et al., 2019) and (filled; Lodders, 2003). CI abundance is from Lodders (2020). Other
data sources are as in Fig. 1.

                                               20
Fig. 4: Siderophile (red square) and chalcophile (green triangle) element abundances in the bulk
silicate Earth and Mars plotted against 50% condensation temperature (K) of elements at 10 Pa
(Wood et al., 2019). The values for y axes follow the same convention as Fig. 3. Also shown are
time-integrated planetary volatility trends at 1 AU (Earth) and 1.5 AU (Mars) defined by lithophile
elements (Fig. 3). Data sources are as in Fig. 1.

                                                21
Fig. 5: Mineralogy and physical properties of planetary mantles. (A, C) Phase transitions and (B,
D) P- and S-wave velocities and density from surface to the core-mantle boundary in the Earth’s and
Martian mantle. Red lines show temperature profile in the mantle. Figures for Earth are based on the
mantle compositional model of McDonough (2014) and model geotherm from Katsura et al. (2010).
Figures for Mars are reprinted from T. Yoshizaki and W.F. McDonough (2020) The composition
of Mars, Geochimica et Cosmochimica Acta 273, 137–162, Copyright (2020), with permission
from Elsevier. Abbreviations: Brg–bridgmanite; Ca-pv–Ca-perovskite; Cpx–clinopyroxene; C2/c–
high-pressure clinopyroxene; Fp–ferropericlase; Gt–garnet; Ol–olivine; Opx–orthopyroxene; Pl–
plagioclase; Ring–ringwoodite; Sp–spinel; St–stishovite; Wad–wadsleyite.

                                                22
Fig. 6: Interior structure of Earth and Mars. Abbreviations: Ol–olivine; Px–pyroxene; Gt–garnet;
Wad–wadsleyite; Ring–ringwoodite; Brg–bridgmanite; Cpv–Ca-perovskite; Ppv–post-perovskite;
Fper–ferropericlase. The Mars figure is reprinted from T. Yoshizaki and W.F. McDonough (2020)
The composition of Mars, Geochimica et Cosmochimica Acta 273, 137–162, Copyright (2020),
with permission from Elsevier.

                                              23
Fig. 7: Composition of metallic cores of Earth (McDonough, 2014), Mars (Yoshizaki and Mc-
Donough, 2020) and NC-group iron meteorite parent bodies (Wasson and Richardson, 2001;
Chabot, 2018). (A) Siderophile and chalcophile element and (B) nominally lithophile element abun-
dances. Elemental abundances are normalized to CI chondrite composition and Ni. Elements are
arranged by their 50% nebular condensation temperatures (Wood et al., 2019), while W and Re
are replaced to highlight chondritic highly siderophile element abundance in planetary cores. Type
IIIAB and IVA iron meteorites have fractional crystallization (i.e., magmatic) origin and are isotopi-
cally classified as NC group meteorites (Burkhardt et al., 2011; Kruijer et al., 2017a; Poole et al.,
2017).

                                                 24
Fig. 8: Composition of metallic cores of Earth (McDonough, 2014), Mars (Yoshizaki and Mc-
Donough, 2020) and fractionally crystallized (i.e., magmatic) iron meteorites (Goldstein et al.,
2009). (A) Ge versus Ni. (B) Ir versus Ni. Iron meteorites are classified into non-carbonaceous
(NC: in red) and carbonaceous (CC: in blue) groups based on their Mo and W isotopic compositions
(Burkhardt et al., 2011; Kruijer et al., 2017a; Poole et al., 2017). Silicate-bearing (non-magmatic)
iron groups are not shown. The IIG iron meteorites are also excluded from the plot because isotopic
data to classify them into NC or CC group are not available. CI chondritic ratios are from Lodders
(2020).

                                                25
You can also read