Kevin J. Walsh

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.

THE LINK BETWEEN PLANETARY SYSTEMS, DUSTY WHITE DWARFS, AND METAL-POLLUTED WHITE DWARFS

It has long been suspected that metal-polluted white dwarfs (types DAZ, DBZ, and DZ) and white dwarfs with dusty disks possess planetary systems, but a specific physical mechanism by which planetesimals are perturbed close to a white dwarf has not yet been fully posited. In this paper, we demonstrate that mass loss from a central star during post-main-sequence evolution can sweep planetesimals into interior mean motion resonances with a single giant planet. These planetesimals are slowly removed through chaotic excursions of eccentricity that in time create radial orbits capable of tidally disrupting the planetesimal. Numerical N-body simulations of the solar system show that a sufficient number of planetesimals are perturbed to explain white dwarfs with both dust and metal pollution, provided other white dwarfs have more massive relic asteroid belts. Our scenario requires only one Jupiter-sized planet and a sufficient number of asteroids near its 2:1 interior mean motion resonance. Finally, we show that once a planetesimal is perturbed into a tidal crossing orbit, it will become disrupted after the first pass of the white dwarf, where a highly eccentric stream of debris forms the main reservoir for dust-producing collisions. These simulations, in concert with observations of white dwarfs, place interesting limits on the frequency of planetary systems around main-sequence stars, the frequency of planetesimal belts, and the probability that dust may obscure future terrestrial planet finding missions.

Highly siderophile elements in Earth’s mantle as a clock for the Moon-forming impact

According to the generally accepted scenario, the last giant impact on Earth formed the Moon and initiated the final phase of core formation by melting Earth’s mantle. A key goal of geochemistry is to date this event, but different ages have been proposed. Some argue for an early Moon-forming event, approximately 30 million years (Myr) after the condensation of the first solids in the Solar System, whereas others claim a date later than 50 Myr (and possibly as late as around 100 Myr) after condensation. Here we show that a Moon-forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We use a large number of N-body simulations to demonstrate a relationship between the time of the last giant impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This relationship exists within both the classical scenario and the Grand Tack scenario of terrestrial planet formation, and holds across a wide range of disk conditions. The concentration of highly siderophile elements (HSEs) in Earth’s mantle constrains the mass of chondritic material added to Earth during Late Accretion. Using HSE abundance measurements, we determine a Moon-formation age of 95 ± 32 Myr after condensation. The possibility exists that some late projectiles were differentiated and left an incomplete HSE record in Earth’s mantle. Even in this case, various isotopic constraints strongly suggest that the late-accreted mass did not exceed 1 per cent of Earth’s mass, and so the HSE clock still robustly limits the timing of the Moon-forming event to significantly later than 40 Myr after condensation.

Growing the terrestrial planets from the gradual accumulation of submeter-sized objects

Significance The fact that Mars is so much smaller than both Earth and Venus has been a long-standing puzzle of terrestrial planet formation. Here we show that a new mode of planet formation known as “Viscous Stirred Pebble Accretion,” which has recently been shown to produce the giant planets, also naturally explains the small size of Mars and the low mass of the asteroid belt. Thus there is a unified model that can be used to explain all of the basic properties of our solar system. Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 astronomical units (AU) should roughly be the same mass as Earth. Recently, a new model called Viscously Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants. This model envisions that the cores of the giant planets formed from 100- to 1,000-km bodies that directly accreted a population of pebbles—submeter-sized objects that slowly grew in the protoplanetary disk. Here we apply this model to the terrestrial planet region and find that it can reproduce the basic structure of the inner solar system, including a small Mars and a low-mass asteroid belt. Our models show that for an initial population of planetesimals with sizes similar to those of the main belt asteroids, VSPA becomes inefficient beyond ∼ 1.5 AU. As a result, Mars’s growth is stunted, and nothing large in the asteroid belt can accumulate.

The effect of an early planetesimal-driven migration of the giant planets on terrestrial planet formation

The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is assumed to have occurred immediately after the dissipation of the nebular disk – i.e. “early” with respect to the timing of the late heavy bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer between the planets and scattered planetesimals. Our model forms the terrestrial planets from a disk of material which stretchs from 0.3–4.0 AU, evenly split in mass between planetesimals and planetary embryos. Jupiter and Saturn are initially at 5.4 and 8.7 AU respectively, on orbits with eccentricities comparable to the current ones, and migrate to 5.2 and 9.4 AU with an e-folding time of 5 Myr. Unfortunately, the terrestrial planets formed in the simulations are not good analogs for the current solar system, with Mars typically being much too massive. Moreover, the final distribution of the planetesimals remaining in the asteroid belt is inconsistent with the observed distribution of asteroids. This argues that, even if giant planet migration had occurred early, the real evolution of the giant planets would have to have been of the “jumping-Jupiter” type, i.e. the increase in orbital separation between Jupiter and Saturn had to be dominated by encounters between Jupiter and a third, Neptune-mass planet. This result was already demonstrated for late migrations occuring at the LHB time by previous work, and this paper shows those conclusions hold for early migration as well.

OSIRIS-REx: Sample Return from Asteroid (101955) Bennu

In May of 2011, NASA selected the Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer (OSIRIS-REx) asteroid sample return mission as the third mission in the New Frontiers program. The other two New Frontiers missions are New Horizons, which explored Pluto during a flyby in July 2015 and is on its way for a flyby of Kuiper Belt object 2014 MU69 on January 1, 2019, and Juno, an orbiting mission that is studying the origin, evolution, and internal structure of Jupiter. The spacecraft departed for near-Earth asteroid (101955) Bennu aboard an United Launch Alliance Atlas V 411 evolved expendable launch vehicle at 7:05 p.m. EDT on September 8, 2016, on a seven-year journey to return samples from Bennu. The spacecraft is on an outbound-cruise trajectory that will result in a rendezvous with Bennu in November 2018. The science instruments on the spacecraft will survey Bennu to measure its physical, geological, and chemical properties, and the team will use these data to select a site on the surface to collect at least 60 g of asteroid regolith. The team will also analyze the remote-sensing data to perform a detailed study of the sample site for context, assess Bennu’s resource potential, refine estimates of its impact probability with Earth, and provide ground-truth data for the extensive astronomical data set collected on this asteroid. The spacecraft will leave Bennu in 2021 and return the sample to the Utah Test and Training Range (UTTR) on September 24, 2023.

Will Comet ISON (C/2012 S1) Survive Perihelion?

On 2013 November 28 Comet ISON (C/2012 S1) will pass by the Sun with a perihelion distance of 2.7 solar radii. Understanding the possible outcomes for the comets response to such a close passage by the Sun is important for planning observational campaigns and for inferring ISONs physical properties. We present new numerical simulations and interpret them in context with the historical track record of comet disruptions and of sungrazing comet behavior. Historical data suggest that sizes below ~200?m are susceptible to destruction by sublimation driven mass loss, while we find that for ISONs perihelion distance, densities lower than 0.1?g?cm?3 are required to tidally disrupt a retrograde or non-spinning body. Such low densities are substantially below the range of the best-determined comet nucleus densities, though dynamically new comets such as ISON have few measurements of physical properties. Disruption may occur for prograde rotation at densities up to 0.7?g?cm?3, with the chances of disruption increasing for lower density, faster prograde rotation, and increasing elongation of the nucleus. Given current constraints on ISONs nucleus properties and the typically determined values for these properties among all comets, we find tidal disruption to be unlikely unless other factors (e.g., spin-up via torquing) affect ISON substantially. Whether or not disruption occurs, the largest remnant must be big enough to survive subsequent mass loss due to sublimation in order for ISON to remain a viable comet well after perihelion.

Ridge formation and de-spinning of Iapetus via an impact-generated satellite

Abstract We present a scenario for building the equatorial ridge and de-spinning Iapetus through an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material down to the surface of Iapetus, and then itself tidally evolves outward, thereby helping to de-spin Iapetus. This scenario can de-spin Iapetus an order of magnitude faster than when tides due to Saturn act alone, almost independently of its interior geophysical evolution. Eventually, the satellite is stripped from its orbit by Saturn. The range of satellite and impactor masses required is compatible with the estimated impact history of Iapetus.

Formation and evolution of binary asteroids

Satellites of asteroids have been discovered in nearly every known small body population, and a remarkable aspect of the known satellites is the diversity of their properties. They tell a story of vast differences in formation and evolution mechanisms that act as a function of size, distance from the Sun, and the properties of their nebular environment at the beginning of Solar System history and their dynamical environment over the next 4.5 Gyr. The mere existence of these systems provides a laboratory to study numerous types of physical processes acting on asteroids and their dynamics provide a valuable probe of their physical properties otherwise possible only with spacecraft. Advances in understanding the formation and evolution of binary systems have been assisted by: 1) the growing catalog of known systems, increasing from 33 to nearly 250 between the Merline et al. (2002) Asteroids III chapter and now, 2) the detailed study and long-term monitoring of individual systems such as 1999 KW4 and 1996 FG3, 3) the discovery of new binary system morphologies and triple systems, 4) and the discovery of unbound systems that appear to be end-states of binary dynamical evolutionary paths. Specifically for small bodies (diameter smaller than 10 km), these observations and discoveries have motivated theoretical work finding that thermal forces can efficiently drive the rotational disruption of small asteroids. Long-term monitoring has allowed studies to constrain the systems dynamical evolution by the combination of tides, thermal forces and rigid body physics. The outliers and split pairs have pushed the theoretical work to explore a wide range of evolutionary end-states.

Physical Characterization and Origin of Binary Near-Earth Asteroid (175706) 1996 FG3

The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-?v heliocentric orbit. The orbit of its satellite has provided valuable information about its mass density while its albedo and colors suggest it is primitive or part of the C-complex taxonomic grouping. We extend the physical characterization of this object with new observations of its emission at mid-infrared wavelengths and with near-infrared reflection spectroscopy. We derive an area-equivalent system diameter of 1.90 ? 0.28 km (corresponding to approximate component diameters of 1.83 km and 0.51 km, respectively) and a geometric albedo of 0.039 ? 0.012. (175706) 1996 FG3 was previously classified as a C-type asteroid, though the combined 0.4-2.5??m spectrum with thermal correction indicates classification as B-type; both are consistent with the low measured albedo. Dynamical studies show that (175706) 1996 FG3 most probably originated in the inner main asteroid belt. Recent work has suggested the inner Main Belt (142) Polana family as the possible origin of another low-?v B-type NEA, (101955) 1999 RQ36. A similar origin for (175706) 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and Mars 5:9 mean motion resonances rather than the ?6, and we find this to be a low probability, but possible, origin.