Subsequent to the Moon's formation, late accretion to the terrestrial planets modified their silicate crusts and mantles. We combine dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected mass augmentation onto the terrestrial planets from three sources: planetesimals left over from primary accretion, asteroids from the E-belt, and comets arriving from the outer Solar System. We present estimates of the amount of cometary material striking the terrestrial planets in an early episode of planetesimal-driven giant planet migration (Mojzsis et al., 2019). The Moon and Mars suffer proportionally higher cometary accretion than Venus and Earth. The background mass addition from small leftover planetesimals to Earth and Mars is far less than independent estimates based on their mantle abundances of highly-siderophile elements. This supports the theory that both planets were struck by single large bodies that delivered most of their mass augmentation since primary accretion. We present fits to the impact chronologies valid from 4500 Ma to ca. 3700 Ma. The lunar timeline thus obtained is at odds with both the calibrated Neukum (Neukum et al., 2001) and Werner (Werner et al., 2014; Werner, 2019) chronologies. For Mars, the match with its calibrated Werner chronology is no better; by increasing the mass of the E-belt by a factor of four the dynamical lunar and martian chronologies are in line with that of Werner (2019) and match constraints from the current population of Hungaria asteroids. Yet, neither of our dynamical timelines fit well with that of Neukum. The dynamical lunar and martian chronologies are also different from each other. Consequently, the usual extrapolation of such chronologies from one planetary body to the other is technically inappropriate.

The widespread presence of ribonucleic acid (RNA) catalysts and cofactors in Earth's biosphere today suggests that RNA was the first biopolymer to support Darwinian evolution. However, most "path-hypotheses" to generate building blocks for RNA require reduced nitrogen-containing compounds not made in useful amounts in the CO2-N2-H2O atmospheres of the Hadean. We review models for Earth's impact history that invoke a single ~10^23 kg impactor (Moneta) to account for measured amounts of platinum, gold, and other siderophilic ("iron-loving") elements on the Earth and Moon. If it were the last sterilizing impactor, Moneta would have reduced the atmosphere but not its mantle, opening a "window of opportunity" for RNA synthesis, a period when RNA precursors rained from the atmosphere to land holding oxidized minerals that stabilize advanced RNA precursors and RNA. Surprisingly, this combination of physics, geology, and chemistry suggests a time when RNA formation was most probable, ~120 +/- 100 million years after Moneta's impact, or ~4.36 +/- 0.1 billion years ago. Uncertainties in this time are driven by uncertainties in rates of productive atmosphere loss and amounts of sub-aerial land.

The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars.

Immediately after their formation, the terrestrial planets experienced intense impact bombardment by comets, leftover planetesimals from primary accretion, and asteroids. This temporal interval in solar system evolution, termed late accretion, thermally and chemically modified solid planetary surfaces and may have impeded the emergence of life on the Hadean Earth. The sources and tempo of late accretion are, however, vague. Here, we present a timeline that relates variably retentive radiometric ages from asteroidal meteorites, to new dynamical models of late accretion that invokes giant planet migration. Reconciliation of the geochronological data with dynamical models shows that giant planet migration immediately leads to an intense 30 Myr influx of comets to the entire solar system. The absence of whole-sale crustal reset ages after 4450 Ma for the most resilient chronometers from Earth, Moon, Mars, Vesta and various meteorite parent bodies confines the onset of giant planet migration to no later than ca. 4480 Ma. Waning impacts from planetesimals, asteroids (and a minor cometary component) continue to strike the inner planets through a protracted monotonic decline in impactor flux; this is in agreement with predictions from crater chronology. Amended global 3-D thermal analytical bombardment models derived from our new impact mass-production functions show that persistent niches for prebiotic chemistry on the early Hadean Earth could endure late accretion for at least the last 4400 Myr.

Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars.

Radiometric dating indicates that Mars accreted in the first ~4 Myr of solar system formation, which coincides with the formation and possible migration of Jupiter. While nebular gas from the protoplanetary disk was still present, Jupiter may have migrated inwards and tacked at 1.5 AU in a 3:2 resonance with Saturn. This migration excited planetary building blocks in the inner solar system, resulting in extensive mixing and planetesimal removal. Here we evaluate the plausible nature of Mars' building blocks, focusing in particular on how its growth was influenced by the formation and migration of Jupiter. We use a combination of dynamical simulations and an isotopic mixing model that traces the accretion of elements with different affinities for metal. Dynamical simulations show that Jupiter's migration causes the late stages of Earth's and Mars' accretion to be dominated by EC-type (enstatite chondrite) material due to the loss of OC (ordinary chondrite) planetesimals. Our analysis of available isotopic data for SNC meteorites shows that Mars consists of approximately 68%$^{+0}_{-39}$ EC + 32%$^{+35}_{-0}$ OC by mass (2{\sigma}). The large uncertainties indicate that isotopic analyses of martian samples are for the most part too imprecise to definitely test model predictions; in particular it remains uncertain whether or not Mars accreted predominantly EC material in the latter stages of its formation history. Dynamical simulations also provide no definitive constraint on Mars' accretion history due to the great variety of dynamical pathways that the martian embryo exhibits. The present work calls for new measurements of isotopic anomalies in SNC meteorites targeting siderophile elements (most notably Ni, Mo and Ru) to constrain Mars' accretion history and its formation location.

Impact bombardment on Mercury in the solar system late accretion phase (ca. 4.4 to 3.8 Ga) caused considerable mechanical, chemical and thermal reworking of its silicate reservoirs (crust and mantle). Depending on the frequency, size and velocity of such impactors, effects included regional and global scale crustal melting, and thermal perturbations of the mercurian mantle. We use a 3D transient heating model to test the effects of two bombardment scenarios on early (pre Tolstojan) Mercury mantle and crust. Results show that rare impacts by the largest (greater than 100 km diameter) bodies deliver sufficient heat to the shallow mercurian mantle producing high-temperature ultra-magnesian (komatiitic s.s.) melts. Impact heating leading to effusive (flood) volcanism can account for eponymous High-Magnesium Region (HMR) observed during the MErcury Surface, Space Environment, GEochemistry Ranging (MESSENGER) mission. We find that late accretion to Mercury induced volumetrically significant crustal melting (less than 58 vol.percent), mantle heating and melt production, which, combined with extensive resurfacing (less than 100 percent), also explains why its oldest cratering record was effectively erased, consistent with crater-counting statistics.

We conducted zircon (U-Th)/He (ZHe) analysis of lunar impact-melt breccia 14311 with the aim of leveraging radiation damage accumulated in zircon over extended intervals to detect low-temperature or short-lived impact events that have previously eluded traditional isotopic dating techniques. Our ZHe data record a coherent date vs. effective Uranium concentration (eU) trend characterized by >3500 Ma dates from low (<75 ppm) eU zircon grains, and ca. 110 Ma dates for high (>100 ppm) eU grains. A progression between these date populations is apparent for intermediate (75-100 ppm) eU grains. Thermal history modeling constrains permissible temperatures and cooling rates during and following impacts. Modeling shows that the data are most simply explained by impact events at ca. 3950 Ma and ca. 110 Ma, and limits allowable temperatures of heating events between 3950-110 Ma. Modeling of solar cycling thermal effects at the lunar surface precludes this as the explanation for the ca. 110 Ma ZHe dates. We propose a sample history characterized by zircon resetting during the ca. 3950 Ma Imbrium impact event, with subsequent heating during an impact at ca. 110 Ma that ejected the sample to the vicinity of its collection site. Our data show that zircon has the potential to retain 4He over immense timescales (>3950 Myrs), thus providing a valuable new thermochronometer for probing the impact histories of lunar samples, and martian or asteroidal meteorites.

Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is sequestered into iron cores. Here, we analyze the mechanical and chemical fate of the core of such a Late Veneer impactor, and trace how its HSEs are suspended - and thus pollute - the mantle. For the statistically most-likely oblique collision (about 45degree), the impactor's core elongates and thereafter disintegrates into a metallic hail of small particles (about 10 m). Some strike the orbiting Moon as sesquinary impactors, but most re-accrete to Earth as secondaries with further fragmentation. We show that a single oblique impactor provides an adequate amount of HSEs to the primordial terrestrial silicate reservoirs via oxidation of (<m-sized) metal particles with a hydrous, pre-impact, early Hadean Earth.

Once the terrestrial planets had mostly completed their assembly, bombardment continued by planetesimals left-over from accretion. Highly siderophile element (HSE) abundances in Mars' mantle imply its late accretion supplement was 0.8 wt.%; Earth and the Moon obtained an additional 0.7 wt.% and 0.02 wt.%, respectively. The disproportionately high Earth/Moon accretion ratio is explicable by stochastic addition of a few remaining Ceres-sized bodies that preferentially targeted Earth. Here we show that Mars' late accretion budget also requires a colossal impact, a plausible visible remnant of which is the hemispheric dichotomy. The addition of sufficient HSEs to the martian mantle entails an impactor of at least 1200 km in diameter to have struck Mars before ca. 4430 Ma, by which time crust formation was well underway. Thus, the dichotomy could be one of the oldest geophysical features of the martian crust. Ejected debris could be the source material for its satellites.

With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other. However, bulk elemental and isotopic data for martian meteorites demonstrate that key aspects of Mars' composition are markedly different from that of Earth. This suggests that Mars formed outside of the terrestrial feeding zone during primary accretion. It is therefore probable that Mars always remained significantly farther from the Sun than Earth; its growth was stunted early and its mass remained relatively low. Here we identify a potential dynamical pathway that forms Mars in the asteroid belt and keeps it outside of Earth's accretion zone while at the same time accounting for strict age and compositional constraints, as well as mass differences. Our uncommon pathway (approximately 2% probability) is based on the Grand Tack scenario of terrestrial planet formation, in which the radial migration by Jupiter gravitationally sculpts the planetesimal disc at Mars' current location. We conclude that Mars' formation requires a specific dynamical pathway, while this is less valid for Earth and Venus. We further predict that} Mars' volatile budget is most likely different from Earth's and that Venus formed close enough to our planet that it is expected to have a nearly identical composition from common building blocks.

It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a multitude of planetesimals left over from accretion. Here we place limits on the mass and evolution of these planetesimals based on constraints from the highly siderophile element (HSE) budget of the Moon. Outcomes from a combination of N-body and Monte Carlo simulations of planet formation lead us to four key conclusions about the nature of this early epoch. First, matching the terrestrial to lunar HSE ratio requires either that the late veneer on Earth consisted of a single lunar-size impactor striking the Earth before 4.45 Ga, or that it originated from the impact that created the Moon. An added complication is that analysis of lunar samples indicates the Moon does not preserve convincing evidence for a late veneer like Earth. Second, the expected chondritic veneer component on Mars is 0.06 weight percent. Third, the flux of terrestrial impactors must have been low ( <=10^(-6) M_earth/Myr) to avoid wholesale melting of Earth's crust after 4.4~Ga, and to simultaneously match the number of observed lunar basins. This conclusion leads to an Hadean eon which is more clement than assumed previously. Last, after the terrestrial planets had fully formed, the mass in remnant planetesimals was ~10^(-3) M_earth, lower by at least an order of magnitude than most previous models suggest. Our dynamically and geochemically self-consistent scenario requires that future N-body simulations of rocky planet formation either directly incorporate collisional grinding or rely on pebble accretion.

The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here we have investigated how the model fares under a range of initial conditions and migration course-change (`tack') locations. We have run a large number of N-body simulations with a tack location of 1.5 AU and 2 AU and tested initial conditions using equal mass planetary embryos and a semi-analytical approach to oligarchic growth. We make use of a recent model of the protosolar disc that takes account of viscous heating, include the full effect of type 1 migration, and employ a realistic mass-radius relation for the growing terrestrial planets. Results show that the canonical tack location of Jupiter at 1.5 AU is inconsistent with the most massive planet residing at 1 AU at greater than 95% confidence. This favours a tack farther out at 2 AU for the disc model and parameters employed. Of the different initial conditions, we find that the oligarchic case is capable of statistically reproducing the orbital architecture and mass distribution of the terrestrial planets, while the equal mass embryo case is not.