The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores (kcore_emb) and the initial embryo masses in N-body simulations, rather than the simulation type, where higher values of kcore_emb and larger initial embryo masses correspond to higher concentrations of Ni, Co, Mo, and W in Earth analog mantles and higher concentrations of Si and O in Earth analog cores. As a result, larger initial embryo masses require smaller values of kcore_emb to match Earth's mantle composition. On the other hand, compositions of Earth analog cores are sensitive to the temperatures of equilibration and fO2 of accreting material. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. We use our compiled simulations to explore the relationship between initial embryo masses and the melting history of Earth analogs, where the complex interplay between the timing between impacts, magma ocean lifetimes, and volatile delivery could affect the compositions of Earth analogs formed from different simulation types. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core (kcore) and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact.
Delivery of significant volumes of howardite-eucrite-diogenite (HED) meteorites from Vesta to Earth is linked to the largest impact events on that protoplanet, especially the giant Rheasilvia basin. Dawn mapping of Rheasilvia reveals a well-preserved impact structure with extensive deformation patterns and moderate superposed cratering. Spiral and radial deformation features across the floor of Rheasilvia are also common to complex craters on Ceres and icy satellites, as well as larger bodies like the Moon, Mars, and Mercury, when floor units are not buried by impact melt and debris. These deformation features predate the solidification of such melt/debris mixtures and are important structural elements in complex craters related to failure in converging floor material during late stages of crater formation. The similarity and commonality of such features and lack of preserved ancient units within Rheasilvia indicate that these patterns formed during the latter stages of the impact event and that the original basin floor and subsequent cratering record are preserved. Updated counts and interpretation of this observed superposed cratering indicate an similar to 0.8-0.9 Ga age of basin formation, supporting a young impact and ejection age. This young age is consistent with recent Ar-Ar ages of feldspar grains in HED (Lindsay et al., 2015), and it may also be consistent with ancient bulk Ar-Ar ages (>similar to 3.2 Ga) for some HEDs because resetting such ages requires prolonged heating from the residual impact heating of rocks left on Vesta during earlier impacts. The complex morphology and recent formation of Rheasilvia basin are also consistent with dynamical estimates for delivery of these meteorites to the Earth-Moon system.
Determining how and when Mars formed has been a long-standing challenge for planetary scientists. The size and orbit of Mars are difficult to reproduce in classical simulations of planetary accretion, and this has inspired models of inner solar system evolution that are tuned to produce Mars-like planets. However, such models are not always coupled to geochemical constraints. Analyses of Martian meteorites using the extinct hafnium-tungsten (Hf-W) radioisotopic system, which is sensitive to the timing of core formation, have indicated that the Martian core formed within a few million years of the start of the solar system itself. This has been interpreted to suggest that, unlike Earth's protracted accretion, Mars grew to its modern size very rapidly. These arguments, however, generally rely on simplified growth histories for Mars. Here, we combine likely accretionary histories from a large number of N-body simulations with calculations of metal-silicate partitioning and Hf-W isotopic evolution during core formation to constrain the range of conditions that could have produced Mars. We find that there is no strong correlation between the final masses or orbits of simulated Martian analogs and their 182W anomalies, and that it is readily possible to produce Mars-like Hf-W isotopic compositions for a variety of accretionary conditions. The Hf-W signature of Mars is very sensitive to the oxygen fugacity (fO2) of accreted material because the metal-silicate partitioning behavior of W is strongly dependent on redox conditions. The average fO2 of Martian building blocks must fall in the range of 1.3-1.6 log units below the iron-wüstite buffer to produce a Martian mantle with the observed Hf/W ratio. Other geochemical properties (such as sulfur content) also influence Martian 182W signatures, but the timing of accretion is a more important control. We find that while Mars must have accreted most of its mass within ~5 million years of solar system formation to reproduce the Hf-W isotopic constraints, it may have continued growing afterwards for over 50 million years. There is a high probability of simultaneously matching the orbit, mass, and Hf-W signature of Mars even in cases of prolonged accretion if giant impactor cores were poorly equilibrated and merged directly with the proto-Martian core.
High-resolution mapping of Ceres, Vesta and the icy satellites of Saturn, Uranus and Pluto reveals a rich variety of well-preserved impact crater morphologies on these low gravity bodies. These objects provide a natural laboratory to study effects of composition on crater formation processes under similar surface gravity conditions (though mean impact velocities vary by several factors). Simple craters occur on all these bodies but subtle differences in morphology on Ceres and Vesta are recognized. Immature complex craters (with large floor mounds but not terraces or conical central peaks) occur on Vesta and while smaller than predicted are consistent with its silicate composition. Asymmetric simple craters (with incomplete scarp development) on all bodies are likely related to differential overburden stresses in the rim, and their occurrence is consistent with lower crustal strength on icy bodies including Ceres. Immature and mature complex craters exhibit increasing degrees of complexity, including spiral floor deformation patterns (related to failure in converging floor material), central peaks, and impact melt. Cerean crater morphologic types and simple-complex transition diameters are smaller than on Vesta but similar to those on icy satellites, indicating a much weaker rheology for Ceres' outer layers under impact conditions. These are consistent with geophysical indications of a low-density water ice and probably clathrate rich outer shell. Fluidized floor deposits (impact melt or melt-solid mixtures) are significant in craters >25 km across on Ceres but absent on Saturn satellites. Central pit craters are common on Ceres (at diameters of similar to 75 to 150 km consistent with gravity scaling from the larger Galilean satellites) but are absent on Saturnian satellites and Charon. The contrasting impact melt and central pit behaviors on Ceres and Saturn's moons is contrary to expectation given the higher impact velocities at Saturn but might be related to lower internal temperatures, or the higher fraction of non-ice material on Ceres. The correlation or scaling of transition diameters to surface gravity is near -0.65 rather than -1, perhaps due to increased porosity on lower gravity bodies. The fundamental similarity of crater morphologies on Ceres and icy satellites, however, indicates that the weaker rheology of water ice results in similar craters even if the non-(ice+clathrate) components are as high as similar to 30 vol%.
Hydrothermal processes in impact environments on water-rich bodies such as Mars and Earth are relevant to the origins of life. Dawn mapping of dwarf planet (1) Ceres has identified similar deposits within Occator crater. Here we show using Dawn high-resolution stereo imaging and topography that Ceres' unique composition has resulted in widespread mantling by solidified water- and salt-rich mud-like impact melts with scattered endogenic pits, troughs, and bright mounds indicative of outgassing of volatiles and periglacial-style activity during solidification. These features are distinct from and less extensive than on Mars, indicating that Occator melts may be less gas-rich or volatiles partially inhibited from reaching the surface. Bright salts at Vinalia Faculae form thin surficial precipitates sourced from hydrothermal brine effusion at many individual sites, coalescing in several larger centers, but their ages are statistically indistinguishable from floor materials, allowing for but not requiring migration of brines from deep crustal source(s).
The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanet systems will allow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine, usingN-body simulations, how the properties of planetary systems are determined during the final stages of assembly, when planets accrete from embryos and planetesimals. While accretion is a chaotic process, trends emerge allowing certain features of an ensemble of planetary systems to provide a memory of the initial distribution of solid mass around a star prior to accretion. We also useepos, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, via either planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes (>1 au) has a profound impact on the observed planetary architectures and thus on the formation and delivery of volatiles to possible habitable worlds.
The planetary building blocks that formed in the terrestrial planet region were likely very dry, yet water is comparatively abundant on Earth. Here we review the various mechanisms proposed for the origin of water on the terrestrial planets. Various in-situ mechanisms have been suggested, which allow for the incorporation of water into the local planetesimals in the terrestrial planet region or into the planets themselves from local sources, although all of those mechanisms have difficulties. Comets have also been proposed as a source, although there may be problems fitting isotopic constraints, and the delivery efficiency is very low, such that it may be difficult to deliver even a single Earth ocean of water this way. The most promising route for water delivery is the accretion of material from beyond the snow line, similar to carbonaceous chondrites, that is scattered into the terrestrial planet region as the planets are growing. Two main scenarios are discussed in detail. First is the classical scenario in which the giant planets begin roughly in their final locations and the disk of planetesimals and embryos in the terrestrial planet region extends all the way into the outer asteroid belt region. Second is the Grand Tack scenario, where early inward and outward migration of the giant planets implants material from beyond the snow line into the asteroid belt and terrestrial planet region, where it can be accreted by the growing planets. Sufficient water is delivered to the terrestrial planets in both scenarios. While the Grand Tack scenario provides a better fit to most constraints, namely the small mass of Mars, planets may form too fast in the nominal case discussed here. This discrepancy may be reduced as a wider range of initial conditions is explored. Finally, we discuss several more recent models that may have important implications for water delivery to the terrestrial planets.
The geologic context of red organic‐rich materials ( ROR ) found across an elongated 200 km region on Ceres is evaluated with spectral information from the multispectral framing camera ( FC ) and the visible and near‐infrared mapping spectrometer ( VIR ) of Dawn. Discrete areas of ROR materials are found to be associated with small fresh craters less than a few hundred meters in diameter. Regions with the highest concentration of discrete ROR areas exhibit a weaker diffuse background of ROR materials. The observed pattern could be consistent with a field of secondary impacts, but no appropriate primary crater has been found. Both endogenic and exogenic sources are being considered for these distinctive organic materials.
Five decades of observations of Ceres suggest that the dwarf planet has a composition similar to carbonaceous meteorites and may have an ice-rich outer shell protected by a silicate layer. NASA's Dawn spacecraft has detected ubiquitous clays, carbonates and other products of aqueous alteration across the surface of Ceres, but surprisingly it has directly observed water ice in only a few areas. Here we use Dawn Framing Camera observations to analyse lobate morphologies on Ceres' surface and we infer the presence of ice in the upper few kilometres of Ceres. We identify three distinct lobate morphologies that we interpret as surface flows: thick tongue-shaped, furrowed flows on steep slopes; thin, spatulate flows on shallow slopes; and cuspate sheeted flows that appear fluidized. The shapes and aspect ratios of these flows are different from those of dry landslides-including those on ice-poor Vesta-but are morphologically similar to ice-rich flows on other bodies, indicating the involvement of ice. Based on the geomorphology and poleward increase in prevalence of these flows, we suggest that the shallow subsurface of Ceres is comprised of mixtures of silicates and ice, and that ice is most abundant near the poles.
Ceres is the largest asteroid belt object, and the Dawn spacecraft observed Ceres since 2015. Dawn observed two morphologically distinct linear features on Ceres's surface: secondary crater chains and pit chains. Pit chains provide unique insights into Ceres's interior evolution. We interpret pit chains called the Samhain Catenae as the surface expression of subsurface fractures. Using the pit chains' spacings, we estimate that the localized thickness of Ceres's fractured, outer layer is approximately ≥58 km, at least ~14 km greater than the global average. We hypothesize that extensional stresses, induced by a region of upwelling material arising from convection/diapirism, formed the Samhain Catenae. We derive characteristics for this upwelling material, which can be used as constraints in future interior modeling studies. For example, its predicted location coincides with Hanami Planum, a high‐elevation region with a negative residual gravity anomaly, which may be surficial evidence for this proposed region of upwelling material.
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.
FOR THEIR ORIGIN. C. M. Pieters, A. Nathues, G. Thangiam, H. Hoffman, C. De Sanctis, E. Ammannito, H. Hiesinger, J. H. Pasckert, D. P. O’Brien, J. C. Castillo-Rogez, O. Ruesch, L. A. McFadden, F. Tosi, F. Zambon, C. A. Raymond, C. T. Russell, Brown Univ., DEEPS, Providence, RI 02912 (Carle_Pieters@brown.edu), Max Planck Institute for Solar System Research, Goettingen, Germany, Istituto di Astrofisica e Planetologia Spaziali, Istituto Nazionale di Astrofisica (INAF), Rome, Italy, University of California, Los Angeles, CA, USA, Westfälische Wilhelms-Universität Münster, Germany, Planetary Science Institute, Tucson, AZ, USA, Jet Propulsion Laboratory, California Inst. of Technology, Pasadena, CA, USA, Goddard Space Flight Center, Greenbelt, MD, USA.
Analysis of Dawn spacecraft Framing Camera image data allows evaluation of the topography and geomorphology of features on the surface of Ceres. The dwarf planet is dominated by numerous craters, but other features are also common. Linear structures include both those associated with impact craters and those that do not appear to have any correlation to an impact event. Abundant lobate flows are identified, and numerous domical features are found at a range of scales. Features suggestive of near-surface ice, cryomagmatism, and cryovolcanism have been identified. Although spectroscopic analysis has currently detected surface water ice at only one location on Ceres, the identification of these potentially ice-related features suggests that there may be at least some ice in localized regions in the crust.