Accretion Of Terrestrial Planets Research Articles

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, εRe,Os,orRuSi > εRe,Os,orRuS which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186Os and 187Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

Early bombardment of the moon: Connecting the lunar crater record to the terrestrial planet formation

The lunar crater record features ∼50 basins. The radiometric dating of Apollo samples indicates that the Imbrium basin formed relatively late — from the planet formation perspective — some ≃3.9 Ga. Here we develop a dynamical model for impactors in the inner solar system to provide context for the interpretation of the lunar crater record. The contribution of cometary impactors is found to be insignificant. Asteroids produced most large impacts on the terrestrial worlds in the last ≃3 Gyr. The great majority of early impactors were rocky planetesimals left behind at ∼0.5–1.5 au after the terrestrial planet accretion. The population of terrestrial planetesimals was reduced by disruptive collisions in the first t∼20 Myr after the gas disk dispersal. We estimate that there were ∼4×105 diameter d>10 km bodies when the Moon formed (total planetesimal mass ∼0.015MEarth at t∼50 Myr). The early bombardment of the Moon was intense. To accommodate ∼50 known basins, the lunar basins that formed before ≃4.35–4.41 Ga must have been erased. The late formation of Imbrium occurs with a ∼15–35% probability in our model. About 20 d>10-km bodies were expected to hit the Earth between 2.5 and 3.5 Ga, which is comparable to the number of known spherule beds in the late Archean. We discuss implications of our model for the lunar/Martian crater chronologies, Late Veneer, and noble gases in the Earth atmosphere.

Formation of Lunar Basins from Impacts of Leftover Planetesimals

The Moon holds important clues to the early evolution of the solar system. Some 50 impact basins (crater diameter D > 300 km) have been recognized on the lunar surface, implying that the early impact flux was much higher than it is now. The basin-forming impactors were suspected to be asteroids released from an inner extension of the main belt (1.8–2.0 au). Here we show that most impactors were instead rocky planetesimals left behind at ∼0.5–1.5 au after the terrestrial planet accretion. The number of basins expected from impacts of leftover planetesimals largely exceeds the number of known lunar basins, suggesting that the first ∼200 Myr of impacts are not recorded on the lunar surface. The Imbrium basin formation (age ≃3.92 Gyr; impactor diameter d ≳ 100 km) occurs with a 15%–35% probability in our model. Imbrium must have formed unusually late to have only two smaller basins (Orientale and Schrödinger) forming afterwards. The model predicts ≃20 d > 10 km impacts on the Earth 2.5–3.5 Gyr ago (Ga), which is comparable to the number of known spherule beds in the late Archean.

Early differentiation processes on Mars inferred from silicon isotopes

Accretion of terrestrial planets involved partial or global melting events such as magma oceans or magma ponds. Mars experienced large-scale differentiation very early in its history, as shown by its 146Sm-142Nd and 182Hf-182W record. The broad variations in ε142Nd and ε142W of SNC meteorites highlight the presence of mantle sources that must have remained isolated, at least partly, after the crystallization of a global magma ocean. In this study, we have investigated whether the crystallization of the martian magma ocean could have generated mantle reservoirs characterized by different silicon isotope signatures, as the fractionation of Si isotopes between minerals and melts is known to depend on pressure. Thus, the goal of this study was to investigate whether there were any relationships between magma ocean crystallisation and possible variations in the Si isotope record of SNC meteorites. High resolution silicon isotope measurements were performed on twelve meteorites from the Shergottite, Nakhlite and Chassignite groups using a Neptune Plus MC-ICP-MS in dry plasma mode. The δ30Si values are in good agreement with previous studies but display a narrower range of variations with a mean value at −0.46 ‰ ± 0.07 (2SD). A magma ocean crystallization model shows that the range of δ30Si in SNCs is consistent with that generated by magma ocean crystallisation. In particular, there is a correlation between calculated 147Sm/144Nd for the moderately depleted mantle sources with δ30Si values; this correlation is consistent with the crystallization model if one includes trapped melt in the cumulates. In contrast, enriched shergottites displayed a very homogenous composition in Sm/Nd ratios, despite significant variability in δ30Si. This observation could be related to either fluid-rock interactions or redox effect during magma differentiation. Altogether, silicon isotope compositions of SNC provide new constraints about magma ocean crystallization processes in Mars.

Solar System evolution and terrestrial planet accretion determined by Zr isotopic signatures of meteorites

Nucleosynthetic isotope signatures in meteorites provide key insights into the structure and dynamics of the solar protoplanetary disk and the accretion history of the planets. We present high-precision Zr isotopic data of a comprehensive suite of non-carbonaceous (NC) and carbonaceous (CC) meteorites, and find that various meteorite groups, including enstatite chondrites, exhibit 96Zr enrichments, whereas there is no resolved 91Zr and 92Zr variability. These new Zr isotope data reveal the same fundamental NC-CC dichotomy observed for several other elements, where CC meteorites are more anomalous compared to NC meteorites and are shifted towards the isotopic composition of Ca-Al-rich inclusions (CAIs). For Zr and other elements, the CC composition is reproduced as a mixture of materials with CAI-like and NC-like isotopic compositions in approximately constant proportions, despite these elements exhibiting disparate nucleosynthetic origins or different cosmo- and geochemical behaviors. These constant mixing proportions are inconsistent with an origin of the dichotomy by thermal processing or selective dust-sorting in the disk but indicate mixing of isotopically distinct materials with broadly solar chemical compositions. This corroborates models in which the NC-CC dichotomy reflects time-varied infall from an isotopically heterogeneous molecular cloud. Among NC meteorites, the isotope anomalies in Zr are linearly correlated with those of other elements, which likewise reflects primordial mixing. Lastly, the new Zr isotope data reinforce the notion that Earth incorporated s-process enriched material from the innermost Solar System, which is not represented by known meteorites. By contrast, contributions to Earth and Mars from outer Solar System CC-like materials were limited, indicating that these planets did not form by pebble accretion, which would have led to high CC fractions.

Metal–silicate partitioning of W and Mo and the role of carbon in controlling their abundances in the bulk silicate earth

The liquid metal–liquid silicate partitioning of molybdenum and tungsten during core formation must be well-constrained in order to understand the evolution of Earth and other planetary bodies, in particular because the Hf–W isotopic system is used to date early planetary evolution. The partition coefficients DMo and DW have been suggested to depend on pressure, temperature, silicate and metal compositions, although previous studies have produced varying and inconsistent models. Additionally, the high cationic charges of W and Mo in silicate melts make their partition coefficients particularly sensitive to oxygen fugacity. We combine 48 new high pressure and temperature experimental results with a comprehensive database of previous experiments to re-examine the systematics of Mo and W partitioning, and produce revised partitioning models from the large combined dataset. W partitioning is particularly sensitive to silicate and metallic melt compositions and becomes more siderophile with increasing temperature. We show that W has a 6+ oxidation state in silicate melts over the full experimental fO2 range of ΔIW −1.5 to −3.5. Mo has a 4+ oxidation state, and its partitioning is less sensitive to silicate melt composition but also depends on metallic melt composition. DMo stays approximately constant with increasing depth in Earth. Both W and Mo become more siderophile with increasing C content of the metal: we therefore performed experiments with varying C concentrations and fit epsilon interaction parameters: εCMo = −7.03 ± 0.30 and εCW = −7.38 ± 0.57.W and Mo along with C are incorporated into a combined N-body accretion and core–mantle differentiation model, which already includes the major rock-forming elements as well as S, and moderately and highly siderophile elements. In this model, oxidation and volatility gradients extend through the protoplanetary disk so that Earth accretes heterogeneously. These gradients, as well as the metal–silicate equilibration pressure, are fitted using a least squares optimisation so that the model Earth-like planet reproduces the composition of the bulk silicate Earth (BSE) in terms of 17 simulated element concentrations (Mg, Fe, Si, Ni, Co, Nb, Ta, V, Cr, S, Pt, Pd, Ru, Ir, W, Mo, and C). The effects of the interaction of W and Mo with Si, S, O, and C in metal are included. Using this model with six separate terrestrial planet accretion simulations, we show that W and Mo require the early accreting Earth to be sulfur-depleted and carbon-enriched so that W and Mo are efficiently partitioned into Earth’s core and do not accumulate in the mantle. When this is the case, the produced Earth-like planets possess mantle compositions matching the BSE for all simulated elements. However, there are two distinct groups of estimates of the bulk mantle’s C abundance in the literature: low (∼100 ppm) and high (∼800 ppm), and all six models are consistent with the higher estimated carbon abundance. The low BSE C abundance would be achievable when the effects of the segregation of dispersed metal droplets produced in deep magma oceans by the disproportionation of Fe2+ to Fe3+ plus metallic Fe is included.

Efficient tidal dissipation in Deimos

After its formation Deimos, Mars' second and outer satellite, is expected to have been catastrophically disrupted, possibly multiple times, by planetesimals left over from primary terrestrial planet accretion. The expected last time such a disruption occurred is ca. 4.3 Gyr ago. After this disruption and its subsequent expected coagulation from a thin circum-martian ring, Deimos likely experienced a gradual excitation in its orbital eccentricity from flybys with other leftover planetesimals. The cumulative effect from these encounters indicates that Deimos' eccentricity should have been up to an order of magnitude greater than it is now, with a mean value of 1.5−1.1+2.3×10−3. Based on this eccentricity excitation, Deimos' calculated age, the assumption of synchronous rotation throughout its existence and that the eccentricity is only tidally damped in Deimos and not in Mars, the resulting tidal parameters of Deimos are k2D/QD=3.52−2.89+2.04 × 10−5 (2σ). The large uncertainties in the tidal parameters are caused by the wide eccentricity distribution resulting from encounters with stray planetesimals. Employing a simple tidal model, I calculate that Deimos' shear viscosity is ηD=2.13−0.79+9.77×1012 Pa s, far lower than that of cold rock. Under the assumption that Deimos' interior is akin to a rubble pile its internal rigidity is μrp ~ 120 MPa, while the tidal Love number falls in the range k2D ~ (1 − 7) × 10−5. The low rigidity and viscosity indicate that most of the tidal damping likely occurs in the surface regolith layer. As a consequence, Deimos appears to be near its maximum tidal dissipation efficiency (QD ~ 1). An early episode of chaotic tumbling could have initially lowered Deimos' eccentricity to near zero, but this was likely followed by an excitation of its eccentricity through a past resonance with Phobos ca. 2 Gyr ago. This resonance crossing excited Deimos' eccentricity to 0.002, and subsequent synchronous damping over the next 2.5 Gyr requires k2D/QD ~ 7 × 10−5, which still implies efficient damping in Deimos.

Giant impacts stochastically change the internal pressures of terrestrial planets.

Pressure is a key parameter in the physics and chemistry of planet formation and evolution. Previous studies have erroneously assumed that internal pressures monotonically increase with the mass of a body. Using smoothed particle hydrodynamics and potential field method calculations, we demonstrate that the hot, rapidly rotating bodies produced by giant impacts can have much lower internal pressures than cool, slowly rotating planets of the same mass. Pressures subsequently increase because of thermal and rotational evolution of the body. Using the Moon-forming impact as an example, we show that the internal pressures after the collision could have been less than half that in present-day Earth. The current pressure profile was not established until Earth cooled and the Moon receded, a process that may take up to tens of millions of years. Our work defines a new paradigm for pressure evolution during accretion of terrestrial planets: stochastic changes driven by impacts.

Hydrogen and helium ingassing during terrestrial planet accretion

The amount of atmospheric gas absorbed into a planetary interior during formation depends on the temperature, pressure, composition, and dynamics of its atmosphere, along with the planet mass, its composition, rate of accretion, and the average age of its surface. We develop a boundary layer model for terrestrial planet ingassing during accretion that quantifies the amounts of hydrogen and helium dissolved into a global magma ocean from a hot, gravitationally captured nebular atmosphere in the first few million years of solar system history. Assuming that most of Earth's accretion occurred within the lifetime of its nebular atmosphere, one or more ocean masses of ingassed hydrogen are predicted if the surface pressure is high and the magma surface is young. In addition, Earth would have acquired hundreds of petagrams of helium-3 from the same nebular atmosphere. Our model predicts negligible hydrogen ingassing from a gravitationally captured nebular atmosphere on Mars, but shows that vast amounts of hydrogen ingassing – tens of ocean masses – are possible during the accretion of Super-Earth-sized terrestrial planets.

Mineralogy of volatile‐rich clasts in brecciated meteorites

AbstractMeteoritic breccias are valuable samples as they can contain rare materials from the early solar system as clasts. Volatile‐rich, CI‐ and CM‐like clasts may represent parent body lithologies, which cannot be found as individual meteorites in today's meteorite collections. In order to reveal a better knowledge about the presence and chemical characteristics and variability of volatile (water‐bearing) materials in the early solar system these clasts play an important role. Such materials may have been available as the volatile component during the accretion of terrestrial planets. To understand the distribution of volatile‐rich materials in the solar system, we studied CI‐ and CM‐like clasts in brecciated meteorites including polymict ureilites, HEDs, CR, CB, CH, and ordinary chondrites. CI‐like clasts occur throughout all of the mentioned meteorite groups, whereas the CM‐like clasts have only been identified in HEDs and ordinary chondrites. The abundance of volatile‐rich clasts in general decreases in the order CH > CR > ureilites > HEDs > CB > OC > R. The mineralogy of CI‐like clasts is similar to CI chondrites, but their compositions of phyllosilicates differ. The mineralogy of CM‐like clasts clearly links them to CM chondrites. They must have been delivered to the HED parent body by low‐velocity impacts after differentiation and volcanism, as there is no evidence for high shock and heating processes. Additionally, we propose that CI‐like clasts in the CR, CB, and CH chondrites are a primary component of the appropriate parent bodies (accretionary breccias). Conversely, the CI‐like clasts in polymict ureilites and HEDs represent an infall as (micro)meteorites or as low‐velocity impactors, which happened after the accretion and differentiation of the appropriate parent bodies.

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 homogeneous 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.

Prevalence of chaos in planetary systems formed through embryo accretion

The formation of the solar system’s terrestrial planets has been numerically modeled in various works, and many other studies have been devoted to characterizing our modern planets’ chaotic dynamical state. However, it is still not known whether our planets fragile chaotic state is an expected outcome of terrestrial planet accretion. We use a suite of numerical simulations to present a detailed analysis and characterization of the dynamical chaos in 145 different systems produced via terrestrial planet formation in Kaib and Cowan (2015). These systems were created in the presence of a fully formed Jupiter and Saturn, using a variety of different initial conditions. They are not meant to provide a detailed replication of the actual present solar system, but rather serve as a sample of similar systems for comparison and analysis. We find that dynamical chaos is prevalent in roughly half of the systems we form. We show that this chaos disappears in the majority of such systems when Jupiter is removed, implying that the largest source of chaos is perturbations from Jupiter. Chaos is most prevalent in systems that form 4 or 5 terrestrial planets. Additionally, an eccentric Jupiter and Saturn is shown to enhance the prevalence of chaos in systems. Furthermore, systems in our sample with a center of mass highly concentrated between ∼0.8–1.2 AU generally prove to be less chaotic than systems with more exotic mass distributions. Through the process of evolving systems to the current epoch, we show that late instabilities are quite common in our systems. Of greatest interest, many of the sources of chaos observed in our own solar system (such as the secularly driven chaos between Mercury and Jupiter) are shown to be common outcomes of terrestrial planetary formation. Thus, consistent with previous studies such as Laskar (1996), the solar system’s marginally stable, chaotic state may naturally arise from the process of terrestrial planet formation.

ORBITAL SURVIVAL OF METER-SIZE AND LARGER BODIES DURING GRAVITATIONALLY UNSTABLE PHASES OF PROTOPLANETARY DISK EVOLUTION

A long-standing problem in the collisional accretion of terrestrial planets is the possible loss of m-size bodies through their inward migration onto the protostar as a result of gas drag forces. Such inward migration can be halted, and indeed even reversed, in a protoplanetary disk with local pressure maxima, such as marginally gravitationally unstable (MGU) phases of evolution, e.g., FU Orionis events. Results are presented for a suite of three-dimensional models of MGU disks extending from 1 to 10 AU and containing solid particles with sizes of 1 cm, 10 cm, 1 m, or 10 m, subject to disk gas drag and gravitational forces. These hydrodynamical models show that over disk evolution time scales of years or longer, during which over half the gaseous disk mass is accreted by the protostar, very few 1 and 10 m bodies are lost through inward migration: most bodies survive and orbit stably in the outer disk. A greater fraction of 1 and 10 cm particles are lost to the central protostar during these time periods, as such particles are more closely tied to the disk gas accreting onto the protostar, but even in these cases, a significant fraction survive and undergo transport from the hot inner disk to the cold outer disk, perhaps explaining the presence of small refractory particles in Comet Wild 2. Evidently MGU disk phases offer a means to overcome the m-sized migration barrier to collisional accumulation.

THE ATMOSPHERES OF EARTHLIKE PLANETS AFTER GIANT IMPACT EVENTS

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

No universal minimum-mass extrasolar nebula: evidence againstin situaccretion of systems of hot super-Earths

AbstractIt has been proposed that the observed systems of hot super-Earths formed in situ from high-mass discs. By fitting a disc profile to the entire population of Kepler planet candidates, Chiang & Laughlin constructed a ‘minimum-mass extrasolar nebula’ with surface density profile Σ ∝ r−1.6. Here, we use multiple-planet systems to show that it is inconsistent to assume a universal disc profile. Systems with 3–6 low-mass planets (or planet candidates) produce a diversity of minimum-mass discs with surface density profiles ranging from Σ ∝ r−3.2 to Σ ∝ r0.5 (5th–95th percentile). By simulating the transit detection of populations of synthetic planetary systems designed to match the properties of observed super-Earth systems, we show that a universal disc profile is statistically excluded at high confidence. Rather, the underlying distribution of minimum-mass discs is characterized by a broad range of surface density slopes. Models of gaseous discs can only explain a narrow range of slopes (roughly between r0 and r−1.5). Yet accretion of terrestrial planets in a gas-free environment preserves the initial radial distribution of building blocks. The known systems of hot super-Earths must therefore not represent the structure of their parent gas discs and cannot have predominantly formed in situ. We instead interpret the diversity of disc slopes as the imprint of a process that re-arranged the solids relative to the gas in the inner parts of protoplanetary discs. A plausible mechanism is inward type 1 migration of Mars- to Earth-mass planetary embryos, perhaps followed by a final assembly phase.

Building Terrestrial Planets

This article reviews our current understanding of terrestrial planet formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the review, we combine the results of these theoretical models with geochemical, cosmochemical, and chronological constraints to outline a comprehensive scenario of the early evolution of our solar system. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planet populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and Earth and the origin of water (and other volatiles) on our planet are discussed in depth.

Origin, Internal Structure and Evolution of 4 Vesta

Asteroid 4 Vesta is the only preserved intact example of a large, differentiated protoplanet like those believed to be the building blocks of terrestrial planet accretion. Vesta accreted rapidly from the solar nebula in the inner asteroid belt and likely melted due to heat released due to the decay of 26Al. Analyses of meteorites from the howardite-eucrite-diogenite (HED) suite, which have been both spectroscopically and dynamically linked to Vesta, lead to a model of the asteroid with a basaltic crust that overlies a depleted peridotitic mantle and an iron core. Vesta’s crust may become more mafic with depth and might have been intruded by plutons arising from mantle melting. Constraints on the asteroid’s moments of inertia from the long-wavelength gravity field, pole position and rotation, informed by bulk composition estimates, allow tradeoffs between mantle density and core size; cores of up to half the planetary radius can be consistent with plausible mantle compositions. The asteroid’s present surface is expected to consist of widespread volcanic terrain, modified extensively by impacts that exposed the underlying crust or possibly the mantle. Hemispheric heterogeneity has been observed by poorly resolved imaging of the surface that suggests the possibility of a physiographic dichotomy as occurs on other terrestrial planets. Vesta might have had an early magma ocean but details of the early thermal structure are far from clear owing to model uncertainties and paradoxical observations from the HEDs. Petrological analysis of the eucrites coupled with thermal evolution modeling recognizes two possible mechanisms of silicate-metal differentiation leading to the formation of the basaltic achondrites: equilibrium partial melting or crystallization of residual liquid from the cooling magma ocean. A firmer understanding the plethora of complex physical and chemical processes that contribute to melting and crystallization will ultimately be required to distinguish among these possibilities. The most prominent physiographic feature on Vesta is the massive south polar basin, whose formation likely re-oriented the body axis of the asteroid’s rotation. The large impact represents the likely major mechanism of ejection of fragments that became the HEDs. Observations from the Dawn mission hold the promise of revolutionizing our understanding of 4 Vesta, and by extension, the nature of collisional, melting and differentiation processes in the nascent solar system.

Stable calcium isotopic composition of meteorites and rocky planets

New measurements of mass-dependent calcium isotope effects in meteorites, lunar and terrestrial samples show that Earth, Moon, Mars, and differentiated asteroids (e.g., 4-Vesta and the angrite and aubrite parent bodies) are indistinguishable from primitive ordinary chondritic meteorites at our current analytical resolution (± 0.07‰ SD for the 44Ca/ 40Ca ratio). In contrast, enstatite chondritic meteorites are slightly enriched in heavier calcium isotopes (ca. + 0.5‰) and primitive carbonaceous chondritic meteorites are depleted in heavier calcium isotopes (ca. − 0.5‰). The calcium isotope effects cannot be easily ascribed to evaporation or intraplanetary differentiation processes. The isotopic variations probably survive from the earliest stages of nebular condensation, and indicate that condensation occurred under non-equilibrium (undercooled nebular gas) conditions. Some of this early high-temperature calcium isotope heterogeneity is recorded by refractory inclusions (Niederer and Papanastassiou, 1984) and survived in planetesimals, but virtually none of it survived through terrestrial planet accretion. The new calcium isotope data suggest that ordinary chondrites are representative of the bulk of the refractory materials that formed the terrestrial planets; enstatite and carbonaceous chondrites are not. The enrichment of light calcium isotopes in bulk carbonaceous chondrites implies that their compositions are not fully representative of the solar nebula condensable fraction.

Accretion of the Earth

The origin of the Earth and its Moon has been the focus of an enormous body of research. In this paper I review some of the current models of terrestrial planet accretion, and discuss assumptions common to most works that may require re-examination. Density-wave interactions between growing planets and the gas nebula may help to explain the current near-circular orbits of the Earth and Venus, and may result in large-scale radial migration of proto-planetary embryos. Migration would weaken the link between the present locations of the planets and the original provenance of the material that formed them. Fragmentation can potentially lead to faster accretion and could also damp final planet orbital eccentricities. The Moon-forming impact is believed to be the final major event in the Earth's accretion. Successful simulations of lunar-forming impacts involve a differentiated impactor containing between 0.1 and 0.2 Earth masses, an impact angle near 45 degrees and an impact speed within 10 per cent of the Earth's escape velocity. All successful impacts-with or without pre-impact rotation-imply that the Moon formed primarily from material originating from the impactor rather than from the proto-Earth. This must ultimately be reconciled with compositional similarities between the Earth and the Moon.

Noble gas signature of the Late Heavy Bombardment in the Earth's atmosphere

The Lunar cratering record is consistent with the occurrence of a late heavy bombardment (LHB), which marked the end of terrestrial planet accretion 3.8 billion years ago. However, clear evidence of a LHB on Earth has not yet been identified. Based on a volatile budget of the terres- trial mantle, the atmosphere and hydrosphere we propose that the LHB did indeed occur on Earth and that we are breath- ing its aftermaths. The terrestrial atmosphere and hydro- sphere is enriched in noble gases relative to the abundance of volatiles in the mantle. This enrichment is consistent with the mass delivered to Earth during the LHB, as recently proposed from dynamical modelling (Gomez et al., 2005), if this mate- rial comprised a few Kuiper-belt (cometary) objets (KBOs) mixed in with a population of largely chondritic (i.e. aster- oidal) impactors. The fraction of KBOs necessary to account for the atmospheric composition is, however, much lower (<1%) than the one ( 50%) inferred from modelling.