Showing 1–50 of 50 results for author: Walsh, K J

BYORP and Dissipation in Binary Asteroids: Lessons from DART

Authors: Matija Ćuk, Harrison Agrusa, Rachel H. Cueva, Fabio Ferrari, Masatoshi Hirabayashi, Seth A. Jacobson, Jay McMahon, Patrick Michel, Paul S�nchez, Daniel J. Scheeres, Stephen Schwartz, Kevin J. Walsh, Yun Zhang

Abstract: The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary.… ▽ More The Near-Earth binary asteroid Didymos was the target of a planetary defense demonstration mission DART in September 2022. The smaller binary component, Dimorphos, was impacted by the spacecraft in order to measure momentum transfer in kinetic impacts into rubble piles. DART and associated Earth-based observation campaigns have provided a wealth of scientific data on the Didymos-Dimorphos binary. DART revealed a largely oblate and ellipsoidal shape of Dimorphos before the impact, while the post-impact observations suggest that Dimorphos now has a prolate shape. Here we add those data points to the known properties of small binary asteroids and propose new paradigms of the radiative binary YORP (BYORP) effect as well as tidal dissipation in small binaries. We find that relatively spheroidal bodies like Dimorphos made of small debris may experience a weaker and more size-dependent BYORP effect than previously thought. This could explain the observed values of period drift in several well-characterized binaries. We also propose that energy dissipation in small binaries is dominated by relatively brief episodes of large-scale movement of (likely surface) materials, rather than long-term steady-state tidal dissipation. We propose that one such episode was triggered on Dimorphos by the DART impact. Depending on the longevity of this high-dissipation regime, it is possible that Dimorphos will be more dynamically relaxed in time for the Hera mission than it was in the weeks following the impact. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: Accepted for PSJ

Shaking and Tumbling: Short- and Long-Timescale Mechanisms for Resurfacing of Near-Earth Asteroid Surfaces from Planetary Tides and Predictions for the 2029 Earth Encounter by (99942) Apophis

Authors: R. -L. Ballouz, H. Agrusa, O. S. Barnouin, K. J. Walsh, Y. Zhang, R. P. Binzel, V. J. Bray, D. N. DellaGiustina, E. R. Jawin, J. V. DeMartini, A. Marusiak, P. Michel, N. Murdoch, D. C. Richardson, E. G. Rivera-Valent�n, A. S. Rivkin, Y. Tang

Abstract: Spectral characterization of near-Earth asteroids (NEAs) has revealed a continuum of space-weathered states for the surfaces of S-complex NEAs, with Q-class NEAs, an S-complex subclass, most closely matching the un-weathered surfaces of ordinary chondrite meteorites. Dynamical calculations of the orbital evolution of S-complex NEAs revealed that Q-class NEAs tend to have close encounters with terr… ▽ More Spectral characterization of near-Earth asteroids (NEAs) has revealed a continuum of space-weathered states for the surfaces of S-complex NEAs, with Q-class NEAs, an S-complex subclass, most closely matching the un-weathered surfaces of ordinary chondrite meteorites. Dynamical calculations of the orbital evolution of S-complex NEAs revealed that Q-class NEAs tend to have close encounters with terrestrial planets, suggesting that planetary tides may play a role in refreshing NEA surfaces. However, the exact physical mechanism(s) that drive resurfacing through tidal encounters and the encounter distance at which these mechanisms are effective, has remained unclear. Through the lens of the upcoming (99942) Apophis encounter with Earth in 2029, we investigate the potential for surface mobilization through tidally-driven seismic shaking over short-timescales during encounter and subsequent surface slope evolution over longer-timescales driven by tumbling. We perform multi-scale numerical modeling and find that the 2029 encounter will induce short-term tidally-driven discrete seismic events that lead to high-frequency (>0.1 Hz) surface accelerations that reach magnitudes similar to Apophis' gravity, and that may be detectable by modern seismometers. It is still unclear if the shaking we model translates to widespread particle mobilization and/or lofting. We also find there will be a significant change in Apophis' tumbling spin state that could lead to longer-term surface refreshing in response to tumbling-induced surface slope changes. We propose that through these mechanisms, space-weathered S-class asteroid surfaces may become refreshed through the exposure of unweathered underlying material. These results will be tested by the future exploration of Apophis by NASA's OSIRIS-APEX. △ Less

Submitted 4 October, 2024; v1 submitted 7 June, 2024; originally announced June 2024.

Comments: 39 pages, 11 figures, Accepted by The Planetary Science Journal

EL meteorites do date the giant planet instability

Authors: Chrysa Avdellidou, Marco Delbo, David Nesvorny, Kevin J. Walsh, Alessandro Morbidelli

Abstract: In our recent work, we combined dynamical simulations, meteoritic data and thermal models as well as asteroid observations to argue that the current parent body of the EL meteorites was implanted into the asteroid belt not earlier than 60 Myr after the beginning of the Solar System and that the most likely capture mechanism was the giant planet orbital instability. In the study "The link between A… ▽ More In our recent work, we combined dynamical simulations, meteoritic data and thermal models as well as asteroid observations to argue that the current parent body of the EL meteorites was implanted into the asteroid belt not earlier than 60 Myr after the beginning of the Solar System and that the most likely capture mechanism was the giant planet orbital instability. In the study "The link between Athor and EL meteorites does not constrain the timing of the giant planet instability" that appeared in arXiv, Izidoro and collaborators argue that the implantation of Athor into the asteroid belt does not necessarily require that the giant planet orbital instability occurred at the implantation time. Here we provide further arguments that, in the end, the giant planet instability is still the most likely dynamical process to implant asteroid Athor into the asteroid main belt between 60 and 100 Myr after the beginning of the Solar System. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 4 pages, 1 figure; comment on arXiv:2404.10828

Accretion and Uneven Depletion of the Main Asteroid Belt

Authors: Rogerio Deienno, David Nesvorny, Matthew S. Clement, William F. Bottke, Andre Izidoro, Kevin J. Walsh

Abstract: The main asteroid belt (MAB) is known to be primarily composed of objects from two distinct taxonomic classes, generically defined here as S- and C-complex. The former probably originated from the inner solar system (interior to Jupiter's orbit), while the latter probably from the outer solar system. Following this definition, (4) Vesta, a V-type residing in the inner MAB (a < 2.5 au), is the sole… ▽ More The main asteroid belt (MAB) is known to be primarily composed of objects from two distinct taxonomic classes, generically defined here as S- and C-complex. The former probably originated from the inner solar system (interior to Jupiter's orbit), while the latter probably from the outer solar system. Following this definition, (4) Vesta, a V-type residing in the inner MAB (a < 2.5 au), is the sole D > 500 km object akin to S-complex that potentially formed in-situ. This provides a useful constraint on the number of D > 500 km bodies that could have formed, or grown, within the primordial MAB. In this work we numerically simulate the accretion of objects in the MAB region during the time when gas in the protoplanetary disk still existed, while assuming different MAB primordial masses. We then accounted for the depletion of that population happening after gas disk dispersal. In our analysis, we subdivided the MAB into five sub-regions and showed that the depletion factor varies throughout the MAB. This results in uneven radial- and size-dependent depletion of the MAB. We show that the MAB primordial mass has to be $\lesssim$ 2.14$\times$10$^{-3}$ Earth masses. Larger primordial masses would lead to the accretion of tens-to-thousands of S-complex objects with D > 500 km in the MAB. Such large objects would survive depletion even in the outer sub-regions (a > 2.5 au), thus being inconsistent with observations. Our results also indicate that S-complex objects with D > 200-300 km, including (4) Vesta, are likely to be terrestrial planetesimals implanted into the MAB rather than formed in-situ. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: Accepted for publication in PSJ, 23 pages, 7 figures, 4 tables

Tidal disruption of near-Earth asteroids during close encounters with terrestrial planets

Authors: Mikael Granvik, Kevin J. Walsh

Abstract: Numerical modeling has long suggested that gravitationally-bound (or so-called rubble-pile) near-Earth asteroids (NEAs) can be destroyed by tidal forces during close and slow encounters with terrestrial planets. However, tidal disruptions of NEAs have never been directly observed nor have they been directly attributed to any families of NEAs. Here we show population-level evidence for the tidal di… ▽ More Numerical modeling has long suggested that gravitationally-bound (or so-called rubble-pile) near-Earth asteroids (NEAs) can be destroyed by tidal forces during close and slow encounters with terrestrial planets. However, tidal disruptions of NEAs have never been directly observed nor have they been directly attributed to any families of NEAs. Here we show population-level evidence for the tidal disruption of NEAs during close encounters with the Earth and Venus. Debiased model distributions of NEA orbits and absolute magnitudes based on observations by the Catalina Sky Survey during 2005--2012 underpredict the number of NEAs with perihelion distances coinciding with the semimajor axes of Venus and the Earth. A detailed analysis of the orbital distributions of the excess NEAs shows that their characteristics agree with the prediction for tidal disruptions, and they cannot be explained by observational selection effects or orbital dynamics. Accounting for tidal disruptions in evolutionary models of the NEA population partly bridges the gap between the predicted rate of impacts by asteroids with diameters of tens of meters and observed statistics of fireballs in the same size range. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 8 pages, 5 figures, accepted for publication in ApJL

Journal ref: ApJL 960, 2, L9, 2024

An efficient numerical approach to modeling the effects of particle shape on rubble-pile dynamics

Authors: Julian C. Marohnic, Joseph V. DeMartini, Derek C. Richardson, Yun Zhang, Kevin J. Walsh

Abstract: We present an approach for the inclusion of non-spherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model non-spherical components. Though the method may be applied more generally, we detail our implementation in the existing N-body code pkdgrav. It has long been acknowledged that non-spherical grains confer ad… ▽ More We present an approach for the inclusion of non-spherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model non-spherical components. Though the method may be applied more generally, we detail our implementation in the existing N-body code pkdgrav. It has long been acknowledged that non-spherical grains confer additional shear strength and resistance to flow when compared with spheres. As a result, we expect that rubble-pile asteroids will also exhibit these properties and may behave differently than comparable rubble piles composed of idealized spheres. Since spherical particles avoid some significant technical challenges, most DEM gravity codes have used only spherical particles, or have been confined to relatively low resolutions. We also discuss the work that has gone into improving performance with non-spherical grains, building on pkdgrav's existing leading-edge computational efficiency among DEM gravity codes. This allows for the addition of non-spherical shapes while maintaining the efficiencies afforded by pkdgrav's tree implementation and parallelization. As a test, we simulated the gravitational collapse of 25,000 non-spherical bodies in parallel. In this case, the efficiency improvements allowed for an increase in speed by nearly a factor of three when compared with the naive implementation. Without these enhancements, large runs with non-spherical components would remain prohibitively expensive. Finally, we present the results of several small-scale tests: spinup due to the YORP effect, tidal encounters, and the Brazil-nut Effect. In all cases, we find that the inclusion of non-spherical constituents has a measurable impact on simulation outcomes. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 35 pages, 8 figures

Gaia view of primitive inner-belt asteroid families: Searching for the origins of asteroids Bennu and Ryugu

Authors: Marco Delbo, Chrysa Avdellidou, Kevin J. Walsh

Abstract: Near-Earth asteroids Ryugu and Bennu, were visited, characterised, and sampled by the Hayabusa2 and OSIRIS-REx missions: remote sensing data and sample return analysis showed that both asteroids have primitive, hydrated and organic-rich compositions. The dark families of the inner main belt (IMB) that belong to the spectroscopic C-complex have been claimed to be the sources of both Ryugu and Bennu… ▽ More Near-Earth asteroids Ryugu and Bennu, were visited, characterised, and sampled by the Hayabusa2 and OSIRIS-REx missions: remote sensing data and sample return analysis showed that both asteroids have primitive, hydrated and organic-rich compositions. The dark families of the inner main belt (IMB) that belong to the spectroscopic C-complex have been claimed to be the sources of both Ryugu and Bennu. Hence, there has been large effort to characterise them. Here we used the Gaia Data Release 3 (DR3) asteroid reflectance spectra to investigate the 11 known IMB C-complex families (Chaldaea, Chimaera, Clarissa, Erigone, Eulalia, Klio, Polana, Primordial, Sulamitis, Svea, Tamara). For each family, we extracted the family members that have known geometric visible albedo values and Gaia DR3 data and we created an average reflectance spectrum per family between 370 and 950 nm. The average DR3 reflectance spectra of each family were compared with the previous literature data and to Bennu's and Ryugu's spectra. We found that DR3 reflectance spectra of the IMB C-complex families are in general consistent with previous findings with the only exception of the Svea family. We also showed that the Polana and the Eulalia families can be distinguished in the wavelength region 370 - 500 nm. Among all the IMB C-complex families, we determined that the average reflectance spectra of the Eulalia and Polana families are the most similar to those of Bennu and Ryugu, respectively. In particular, Eulalia family's average spectrum is a good match to Bennu's in the wavelength range 450 - 800 nm, while beyond 800 nm the spectrum of Bennu is bluer than that of Eulalia. Moreover, the spectrum of the Polana family has the smallest discrepancy against the spectrum of Ryugu, although this match is formally unsatisfactory (reduced chi^2 ~ 1.9). △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy and Astrophysics

MSC Class: 85.05

Journal ref: A&A 680, A10 (2023)

SPH simulations of high-speed collisions between asteroids and comets

Authors: J. Rozehnal, M. Brož, D. Nesvorný, K. J. Walsh, D. D. Durda, D. C. Richardson, E. Asphaug

Abstract: We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter… ▽ More We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations. △ Less

Submitted 28 May, 2022; originally announced May 2022.

Journal ref: Icarus 383, 115064, 2022

Exomoons as sources of white dwarf pollution

Authors: Isabella L Trierweiler, Alexandra E Doyle, Carl Melis, Kevin J Walsh, Edward D Young

Abstract: Polluted white dwarfs offer a unique way to study the bulk compositions of exoplanetary material, but it is not always clear if this material originates from comets, asteroids, moons, or planets. We combine N-body simulations with an analytical model to assess the prevalence of extrasolar moons as white dwarf (WD) polluters. Using a sample of observed polluted white dwarfs we find that the extrapo… ▽ More Polluted white dwarfs offer a unique way to study the bulk compositions of exoplanetary material, but it is not always clear if this material originates from comets, asteroids, moons, or planets. We combine N-body simulations with an analytical model to assess the prevalence of extrasolar moons as white dwarf (WD) polluters. Using a sample of observed polluted white dwarfs we find that the extrapolated parent body masses of the polluters are often more consistent with those of many solar system moons, rather than solar-like asteroids. We provide a framework for estimating the fraction of white dwarfs currently undergoing observable moon accretion based on results from simulated white dwarf planetary and moon systems. Focusing on a three-planet white dwarf system of Super-Earth to Neptune-mass bodies, we find that we could expect about one percent of such systems to be currently undergoing moon accretions as opposed to asteroid accretion. △ Less

Submitted 18 May, 2022; v1 submitted 16 May, 2022; originally announced May 2022.

Comments: Submitted to AAS journals, comments welcome

Chemical Habitability: Supply and Retention of Life's Essential Elements During Planet Formation

Authors: Sebastiaan Krijt, Mihkel Kama, Melissa McClure, Johanna Teske, Edwin A. Bergin, Oliver Shorttle, Kevin J. Walsh, Sean N. Raymond

Abstract: Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimatel… ▽ More Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 39 pages, 7 figures; under review for publication as a chapter in Protostars and Planets VII (Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura)

Barrel Instability in Binary Asteroids

Authors: Matija Ćuk, Seth A. Jacobson, Kevin J. Walsh

Abstract: Most close-in planetary satellites are in synchronous rotation, which is usually the stable end-point of tidal despinning. Saturn's moon Hyperion is a notable exception by having a chaotic rotation. Hyperion's dynamical state is a consequence of its high eccentricity and its highly prolate shape (Wisdom et al. 1984). As many binary asteroids also have elongated secondaries, chaotic rotation is exp… ▽ More Most close-in planetary satellites are in synchronous rotation, which is usually the stable end-point of tidal despinning. Saturn's moon Hyperion is a notable exception by having a chaotic rotation. Hyperion's dynamical state is a consequence of its high eccentricity and its highly prolate shape (Wisdom et al. 1984). As many binary asteroids also have elongated secondaries, chaotic rotation is expected for moons in eccentric binaries (Ćuk and Nesvorný, 2010), and a minority of asteroidal secondaries may be in that state (Pravec et al. 2016). The question of secondary rotation is also important for the action of the BYORP effect, which can quickly evolve orbits of synchronous (but not non-synchronous) secondaries (Ćuk and Burns, 2005). Here we report results of a large set of short numerical simulations which indicate that, apart from synchronous and classic chaotic rotation, close-in irregularly-shaped asteroidal secondaries can occupy an additional, intermediate rotational state. In this ``barrel instability'' the secondary slowly rolls along its long axis, while the longest axis is staying largely aligned with the primary-secondary line. This behavior may be more difficult to detect through lightcurves than a fully chaotic rotation, but would likewise shut down BYORP. We show that the binary's eccentricity, separation measured in secondary's radii and the secondary's shape are all important for determining whether the system settles in synchronous rotation, chaotic tumbling, or barrel instability. We compare our results for synthetic asteroids with known binary pairs to determine which of these behaviors may be present in the Near-Earth Asteroid binary population. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: Accepted for PSJ

Internal rubble properties of asteroid (101955) Bennu

Authors: P. Tricarico, D. J. Scheeres, A. S. French, J. W. McMahon, D. N. Brack, J. M. Leonard, P. Antreasian, S. R. Chesley, D. Farnocchia, Y. Takahashi, E. M. Mazarico, D. Rowlands, D. Highsmith, K. Getzandanner, M. Moreau, C. L. Johnson, L. Philpott, E. B. Bierhaus, K. J. Walsh, O. S. Barnouin, E. E. Palmer, J. R. Weirich, R. W. Gaskell, M. G. Daly, J. A. Seabrook , et al. (2 additional authors not shown)

Abstract: Exploration of asteroid (101955) Bennu by the OSIRIS-REx mission has provided an in-depth look at this rubble-pile near-Earth asteroid. In particular, the measured gravity field and the detailed shape model of Bennu indicate significant heterogeneities in its interior structure, compatible with a lower density at its center. Here we combine gravity inversion methods with a statistical rubble-pile… ▽ More Exploration of asteroid (101955) Bennu by the OSIRIS-REx mission has provided an in-depth look at this rubble-pile near-Earth asteroid. In particular, the measured gravity field and the detailed shape model of Bennu indicate significant heterogeneities in its interior structure, compatible with a lower density at its center. Here we combine gravity inversion methods with a statistical rubble-pile model to determine the density and size-frequency distribution (SFD) index of the rubble that constitutes Bennu. The best-fitting models indicate that the SFD of the interior is consistent with that observed on the surface, with a cumulative SFD index of approximately $-2.9$. The rubble bulk density is approximately $1.35$ g/cm$^3$, corresponding to a $12$% macro-porosity. We find the largest rubble particle to be approximately $145$ m, whereas the largest void is approximately $10$ m. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 14 pages, 5 figures, accepted to Icarus

Efficiency characterization of the V-shape asteroid family detection method

Authors: Rogerio Deienno, Kevin J. Walsh, Marco Delbo

Abstract: Following the break up of a parent body, the Yarkovsky effect causes asteroid family members to spread in orbital semimajor axis with a rate often inversely proportional to their diameter. This size dependent semimajor axis drift causes family members to form structures in the semimajor axis vs inverse diameter plane that have the shape of the letter V. The V-shape method consists in finding the b… ▽ More Following the break up of a parent body, the Yarkovsky effect causes asteroid family members to spread in orbital semimajor axis with a rate often inversely proportional to their diameter. This size dependent semimajor axis drift causes family members to form structures in the semimajor axis vs inverse diameter plane that have the shape of the letter V. The V-shape method consists in finding the borders of such V-shapes of unknown center and opening. Although successfully employed to find some very old families in the inner main asteroid belt, the V-shape searching method is very sensitive to many parameters. In this work, we first created and evolved a synthetic asteroid family over billions of years. Then, by adding uncertainties to semimajor axis and diameter of the evolved synthetic family components, we randomly generated additional 99 similar, but not perfectly V-shaped, family clones. We chose a fairly low initial velocity dispersion of 20 m/s for our family. Thus, we can more easily relate the spreading in semimajor axis with the family's age (slope of the evolving V). A synthetic background with an initially randomly distributed components was also created and evolved for 100 Myr. Thus, by setting different levels of ratio of the synthetic family and background asteroids, we derived a detection efficiency map for the V-shape method and determined how sensitive the results can be based on signal-to-noise levels. We also determined optimal parameter values for the method's efficiency. We found that, families older than ~3 Gyr are likely undetectable, with a method efficiency of 50% or less, whereas younger families (0.5--2.5 Gyr) are more easily detected by the method, with an efficiency of >80%. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 21 pages, 16 figures, Accepted for publication in Icarus

Spin-driven evolution of asteroids' top-shapes at fast and slow spins seen from (101955) Bennu and (162173) Ryugu

Authors: Masatoshi Hirabayashi, Ryota Nakano, Eri Tatsumi, Kevin J. Walsh, Olivier S. Barnouin, Patrick Michel, Christine M. Hartzell, Daniel T. Britt, Seiji Sugita, Sei-ichiro Watanabe, William F. Bottke, Daniel J. Scheeres, Ronald-Louis Ballouz, Yuichiro Cho, Tomokatsu Morota, Ellen S. Howell, Dante S. Lauretta

Abstract: Proximity observations by OSIRIS-REx and Hayabusa2 provided clues on the shape evolution processes of the target asteroids, (101955) Bennu and (162173) Ryugu. Their oblate shapes with equatorial ridges, or the so-called top shapes, may have evolved due to their rotational conditions at present and in the past. Different shape evolution scenarios were previously proposed; Bennu's top shape may have… ▽ More Proximity observations by OSIRIS-REx and Hayabusa2 provided clues on the shape evolution processes of the target asteroids, (101955) Bennu and (162173) Ryugu. Their oblate shapes with equatorial ridges, or the so-called top shapes, may have evolved due to their rotational conditions at present and in the past. Different shape evolution scenarios were previously proposed; Bennu's top shape may have been driven by surface processing, while Ryugu's may have been developed due to large deformation. These two scenarios seem to be inconsistent. Here, we revisit the structural analyses in earlier works and fill a gap to connect these explanations. We also apply a semi-analytical technique for computing the cohesive strength distribution in a uniformly rotating triaxial ellipsoid to characterize the global failure of top-shaped bodies. Assuming that the structure is uniform, our semi-analytical approach describes the spatial variations in failed regions at different spin periods; surface regions are the most sensitive at longer spin periods, while interiors fail structurally at shorter spin periods. This finding suggests that the shape evolution of a top shape may vary due to rotation and internal structure, which can explain the different evolution scenarios of Bennu's and Ryugu's top shapes. We interpret our results as the indications of top shapes' various evolution processes. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 27 pages, 1 table, 7 figures, In Press in Icarus (2020)

Journal ref: Icarus, Volume 352, 113946 (2020)

Charon: A brief history of tides

Authors: Alyssa Rose Rhoden, Helle L. Skjetne, Wade G. Henning, Terry A. Hurford, Kevin J. Walsh, S. A. Stern, C. B. Olkin, J. R. Spencer, H. A. Weaver, L. A. Young, K. Ennico, the New Horizons Team

Abstract: In 2015, the New Horizons spacecraft flew past Pluto and its moon Charon, providing the first clear look at the surface of Charon. New Horizons images revealed an ancient surface, a large, intricate canyon system, and many fractures, among other geologic features. Here, we assess whether tidal stresses played a significant role in the formation of tensile fractures on Charon. Although presently in… ▽ More In 2015, the New Horizons spacecraft flew past Pluto and its moon Charon, providing the first clear look at the surface of Charon. New Horizons images revealed an ancient surface, a large, intricate canyon system, and many fractures, among other geologic features. Here, we assess whether tidal stresses played a significant role in the formation of tensile fractures on Charon. Although presently in a circular orbit, most scenarios for the orbital evolution of Charon include an eccentric orbit for some period of time and possibly an internal ocean. Past work has shown that these conditions could have generated stresses comparable in magnitude to other tidally fractured moons, such as Europa and Enceladus. However, we find no correlation between observed fracture orientations and those predicted to form due to eccentricity-driven tidal stress. It thus seems more likely that the orbit of Charon circularized before its ocean froze, and that either tidal stresses alone were insufficient to fracture the surface or subsequent resurfacing remove these ancient fractures. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 12 pages, 5 figures, 2 tables

Potential Themis Family Asteroid Contribution to the Jupiter-Family Comet Population

Authors: Henry H. Hsieh, Bojan Novakovic, Kevin J. Walsh, Norbert Schorghofer

Abstract: Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family me… ▽ More Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 au and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 12 pages, 6 figures; accepted for publication in AJ

Planetesimals to Terrestrial Planets: collisional evolution amidst a dissipating gas disk

Authors: Kevin J. Walsh, Harold F. Levison

Abstract: We present numerical simulations of terrestrial planet formation that examine the growth continuously from planetesimals to planets in the inner Solar System. Previous studies show that the growth will be inside-out, but it is still common practice to assume that the entire inner disk will eventually reach a bi-modal distribution of embryos and planetesimals. For the combinations of disk mass, ini… ▽ More We present numerical simulations of terrestrial planet formation that examine the growth continuously from planetesimals to planets in the inner Solar System. Previous studies show that the growth will be inside-out, but it is still common practice to assume that the entire inner disk will eventually reach a bi-modal distribution of embryos and planetesimals. For the combinations of disk mass, initial planetesimal radius and gas disk lifetime explored in this work the entire disk never reaches a simple bi-modal mass distribution. We find that the inside-out growth is amplified by the combined effects of collisional evolution of solid bodies and interactions with a dissipating gas disk. This leads to oligarchic growth never being achieved in different places of the disk at the same time, where in some cases the disk can simultaneoulsy support chaotic growth and giant impacts inside 1 au and runaway growth beyond 2 au. The planetesimal population is efficiently depleted in the inner disk where embryo growth primarily advances in the presence of a significant gas disk. Further out in the disk growth is slower relative to the gas disk dissipation, resulting in more excited planetesimals at the same stage of growth and less efficient accretion. This same effect drives mass loss due to collisional grinding strongly altering the surface density of the accreted planets relative to the initial mass distribution. This effect decreases the Mars-to-Earth mass ratios compared to previous works with no collisional grinding. Similar to some previous findings utilizing vastly different growth scenarios these simulations produce a first generation of planetary embryos that are stable for 10-20 Myr, or 5-10 e-folding times of the gas dissipation timescale, before having an instability and entering the chaotic growth stage. △ Less

Submitted 2 August, 2019; originally announced August 2019.

Comments: 17 figures

Journal ref: Icarus 329,1 (2019)

The early instability scenario: terrestrial planet formation during the giant planet instability, and the effect of collisional fragmentation

Authors: Matthew S. Clement, Nathan A. Kaib, Sean N. Raymond, John E. Chambers, Kevin J. Walsh

Abstract: The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations… ▽ More The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the "Early Instability" scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets' orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system's terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system. △ Less

Submitted 18 December, 2018; originally announced December 2018.

Comments: 16 pages, 7 figures, 4 tables, accepted for publication in Icarus

Overcoming the Challenges Associated with Image-based Mapping of Small Bodies in Preparation for the OSIRIS-REx Mission to (101955) Bennu

Authors: D. N. DellaGiustina, C. A. Bennett, K. Becker, D. R Golish, L. Le Corre, D. A. Cook, K. L. Edmundson, M. Chojnacki, S. S. Sutton, M. P. Milazzo, B. Carcich, M. C. Nolan, N. Habib, K. N. Burke, T. Becker, P. H. Smith, K. J. Walsh, K. Getzandanner, D. R. Wibben, J. M. Leonard, M. M. Westermann, A. T. Polit, J. N. Kidd Jr., C. W. Hergenrother, W. V. Boynton , et al. (16 additional authors not shown)

Abstract: The OSIRIS-REx Asteroid Sample Return Mission is the third mission in NASA's New Frontiers Program and is the first U.S. mission to return samples from an asteroid to Earth. The most important decision ahead of the OSIRIS-REx team is the selection of a prime sample-site on the surface of asteroid (101955) Bennu. Mission success hinges on identifying a site that is safe and has regolith that can re… ▽ More The OSIRIS-REx Asteroid Sample Return Mission is the third mission in NASA's New Frontiers Program and is the first U.S. mission to return samples from an asteroid to Earth. The most important decision ahead of the OSIRIS-REx team is the selection of a prime sample-site on the surface of asteroid (101955) Bennu. Mission success hinges on identifying a site that is safe and has regolith that can readily be ingested by the spacecraft's sampling mechanism. To inform this mission-critical decision, the surface of Bennu is mapped using the OSIRIS-REx Camera Suite and the images are used to develop several foundational data products. Acquiring the necessary inputs to these data products requires observational strategies that are defined specifically to overcome the challenges associated with mapping a small irregular body. We present these strategies in the context of assessing candidate sample-sites at Bennu according to a framework of decisions regarding the relative safety, sampleability, and scientific value across the asteroid's surface. To create data products that aid these assessments, we describe the best practices developed by the OSIRIS-REx team for image-based mapping of irregular small bodies. We emphasize the importance of using 3D shape models and the ability to work in body-fixed rectangular coordinates when dealing with planetary surfaces that cannot be uniquely addressed by body-fixed latitude and longitude. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 31 pages, 10 figures, 2 tables

Rubble Pile Asteroids

Authors: Kevin J. Walsh

Abstract: The moniker rubble pile is typically applied to all solar system bodies with Diameter between 200m and 10km - where in this size range there is an abundance of evidence that nearly every object is bound primarily by self-gravity with significant void space or bulk porosity between irregularly shaped constituent particles. The understanding of this population is derived from wide-ranging population… ▽ More The moniker rubble pile is typically applied to all solar system bodies with Diameter between 200m and 10km - where in this size range there is an abundance of evidence that nearly every object is bound primarily by self-gravity with significant void space or bulk porosity between irregularly shaped constituent particles. The understanding of this population is derived from wide-ranging population studies of derived shape and spin, decades of observational studies in numerous wavelengths, evidence left behind from impacts on planets and moons and the in situ study of a few objects via spacecraft flyby or rendezvous. The internal structure, however, which is responsible for the name rubble pile, is never directly observed, but belies a violent history. Many or most of the asteroids on near-Earth orbits, and the ones most accessible for rendezvous and in situ study, are likely byproducts of the continued collisional evolution of the Main Asteroid Belt. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 40 pages, 14 figures; Annual Review of Astronomy and Astrophysics, 2018, Vol. 56:593-624

Journal ref: Annual Review of Astronomy and Astrophysics Vol. 56:593-624 (Volume publication date September 2018) Walsh, K. J., 2018, Annual Review of Astronomy and Astrophysics, Vol. 56:593-624

Secular Resonance Between Iapetus and the Giant Planets

Authors: Matija Ćuk, Luke Dones, David Nesvorný, Kevin J. Walsh

Abstract: Using numerical integrations, we find that the orbital eccentricity of Saturn's moon Iapetus undergoes prominent multi-Myr oscillations. We identify the responsible resonant argument to be $\varpi-\varpi_{g5}+Ω-Ω_{eq}$, with the terms being the longitudes of pericenter of Iapetus and planetary secular mode $g_5$, Iapetus's longitude of the node and Saturn's equinox. We find that this argument curr… ▽ More Using numerical integrations, we find that the orbital eccentricity of Saturn's moon Iapetus undergoes prominent multi-Myr oscillations. We identify the responsible resonant argument to be $\varpi-\varpi_{g5}+Ω-Ω_{eq}$, with the terms being the longitudes of pericenter of Iapetus and planetary secular mode $g_5$, Iapetus's longitude of the node and Saturn's equinox. We find that this argument currently (on a $10^7$~yr timescale) appears to librate with a very large amplitude. On longer timescales, the behavior of this resonant angle is strongly dependent on the resonant interaction between Saturn's spin axis and the planetary mode $f_8$, with long-term secular resonance being possible if Saturn's equinox is librating relative to the node of the $f_8$ eigenmode. We present analytical estimates of the dependence of the resonant argument on the orbital elements of Iapetus. We find that this Iapetus-$g_5$ secular resonance could have been established only after the passage of Iapetus through the 5:1 mean-motion resonance with Titan, possibly in the last Gyr. Using numerical simulations, we show that the capture into the secular resonance appears to be a low-probability event. While the Iapetus-$g_5$ secular resonance can potentially help us put new constraints on the past dynamics of the Saturnian system, uncertainties in both the spin axis dynamics of Saturn and the tidal evolution rate of Titan make it impossible to make any firm conclusions about the resonance's longevity and origin. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: Accepted for MNRAS

Mars' Growth Stunted by an Early Giant Planet Instability

Authors: Matthew S. Clement, Nathan A. Kaib, Sean N. Raymond, Kevin J. Walsh

Abstract: Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial plan… ▽ More Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial planet formation. While current simulations of terrestrial accretion can reproduce many observed qualities of the solar system, replicating the small mass of Mars requires modification to standard planet formation models. Here we use 800 dynamical simulations to show that an early instability in the outer solar system strongly influences terrestrial planet formation and regularly yields properly sized Mars analogs. Our most successful outcomes occur when the terrestrial planets evolve an additional 1-10 million years (Myr) following the dispersal of the gas disk, before the onset of the giant planet instability. In these simulations, accretion has begun in the Mars region before the instability, but the dynamical perturbation induced by the giant planets' scattering removes large embryos from Mars' vicinity. Large embryos are either ejected or scattered inward toward Earth and Venus (in some cases to deliver water), and Mars is left behind as a stranded embryo. An early giant planet instability can thus replicate both the inner and outer solar system in a single model. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 18 pages, 9 figures, accepted for publication in Icarus

Asteroid Ryugu Before the Hayabusa2 Encounter

Authors: Koji Wada, Matthias Grott, Patrick Michel, Kevin J. Walsh, Antonella M. Barucci, Jens Biele, J�rgen Blum, Carolyn M. Ernst, Jan T. Grundmann, Bastian Gundlach, Axel Hagermann, Maximilian Hamm, Martin Jutzi, Myung-Jin Kim, Ekkehard K�hrt, Lucille Le Corre, Guy Libourel, Roy Lichtenheldt, Alessandro Maturilli, Scott R. Messenger, Tatsuhiro Michikami, Hideaki Miyamoto, Stefano Mottola, Akiko M. Nakamura, Thomas M�ller , et al. (9 additional authors not shown)

Abstract: Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical… ▽ More Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: 62 pages, 10 figures, 8 tables. Submitted to Space Science Reviews, on 6 April 2018

Size-dependent modification of asteroid family Yarkovsky V-shapes

Authors: Bryce T. Bolin, Alessandro Morbidelli, Kevin J. Walsh

Abstract: The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $Γ$. Following the proposed relationship between $Γ$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size depe… ▽ More The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $Γ$. Following the proposed relationship between $Γ$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size dependence than previous thought, leading to a curved family V-shape boundary in semi-major axis, a, vs. 1/$D$ space. This implies that asteroids are drifting faster at larger sizes than previously considered decreasing on average the known ages of asteroid families. The V-Shape curvature is determined for $>$25 families located throughout the Main Belt to quantify the Yarkovsky size-dependent drift rate. We find that there is no correlation between family age and V-shape curvature. In addition, the V-shape curvature decreases for asteroid families with larger heliocentric distances suggesting that the relationship between $Γ$ and $D$ is weaker in the outer MB possibly due to homogenous surface roughness among family members. △ Less

Submitted 4 March, 2018; v1 submitted 11 October, 2017; originally announced October 2017.

Comments: 68 Pages, 69 Figures, 4 Tables, accepted for publication in Astronomy & Astrophysics. arXiv admin note: text overlap with arXiv:1708.03662

Journal ref: A&A 611, A82 (2018)

Initial velocity V-shapes of young asteroid families

Authors: Bryce T. Bolin, Kevin J. Walsh, Alessandro Morbidelli, Marco Delbo

Abstract: Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape… ▽ More Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape identification technique to separate family shapes that are due to the initial ejection velocity field and those that are due to the Yarkovsky effect. $<$20 Myr-old asteroid families provide an opportunity to study the velocity field of family fragments before they become too dispersed. Only the Karin family's initial velocity field has been determined and scales inversely with diameter, $D^{-1}$. We have applied the V-shape identification technique to constrain young families' initial ejection velocity fields by measuring the curvature of their fragments' V-shape correlation in semi-major axis, $a$, vs. $D^{-1}$ space. Curvature from a straight line implies a deviation from a scaling of $D^{-1}$. We measure the V-shape curvature of 11 young asteroid families including the \FYnospace, Aeolia, Brangane, Brasilia, Clarissa, Iannini, Karin, Konig, Koronis(2), Theobalda and Veritas asteroid families. We find that the majority of asteroid families have initial ejection velocity fields consistent with $\sim D^{-1}$ supporting laboratory impact experiments and computer simulations of disrupting asteroid parent bodies. △ Less

Submitted 26 September, 2017; v1 submitted 11 August, 2017; originally announced August 2017.

Comments: 68 Pages, 25 Figures, 2 Tables, accepted for publication in MNRAS

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

Authors: D. S. Lauretta, S. S. Balram-Knutson, E. Beshore, W. V. Boynton, C. Drouet dAubigny, D. N. DellaGiustina, H. L. Enos, D. R. Gholish, C. W. Hergenrother, E. S. Howell, C. A. Johnson, E. T. Morton, M. C. Nolan, B. Rizk, H. L. Roper, A. E. Bartels, B. J. Bos, J. P. Dworkin, D. E. Highsmith, D. A. Lorenz, L. F. Lim, R. Mink, M. C. Moreau, J. A. Nuth, D. C. Reuter , et al. (23 additional authors not shown)

Abstract: 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 Jan. 1, 2019… ▽ More 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 Jan. 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 August 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 Bennus 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. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Comments: 89 pages, 39 figures, submitted to Space Science Reviews - OSIRIS-REx special issue

Is the Grand Tack model compatible with the orbital distribution of main belt asteroids?

Authors: Rogerio Deienno, Rodney S. Gomes, Kevin J. Walsh, Alessandro Morbidelli, David Nesvorny

Abstract: The Asteroid Belt is characterized by the radial mixing of bodies with different physical properties, a very low mass compared to Minimum Mass Solar Nebula expectations and has an excited orbital distribution. Models of the evolution of the Asteroid Belt show that the origin of its structure is strongly linked to the process of terrestrial planet formation. The Grand Tack model presents a possible… ▽ More The Asteroid Belt is characterized by the radial mixing of bodies with different physical properties, a very low mass compared to Minimum Mass Solar Nebula expectations and has an excited orbital distribution. Models of the evolution of the Asteroid Belt show that the origin of its structure is strongly linked to the process of terrestrial planet formation. The Grand Tack model presents a possible solution to the conundrum of reconciling the small mass of Mars with the properties of the Asteroid Belt, including the mass depletion, radial mixing and orbital excitation. However, while the inclination distribution produced in the Grand Tack model is in good agreement with the one observed, the eccentricity distribution is skewed towards values larger than those found today. Here, we evaluate the evolution of the orbital properties of the Asteroid Belt from the end of the Grand Tack model (at the end of the gas nebula phase when planets emerge from the dispersing gas disk), throughout the subsequent evolution of the Solar System including an instability of the Giant Planets approximately 400 My later. Before the instability, the terrestrial planets were modeled on dynamically cold orbits with Jupiter and Saturn locked in a 3:2 mean motion resonance. The model continues for an additional 4.1 Gy after the giant planet instability. Our results show that the eccentricity distribution obtained in the Grand Tack model evolves towards one very similar to that currently observed, and the semimajor axis distribution does the same. The inclination distribution remains nearly unchanged with a slight preference for depletion at low inclination; this leads to the conclusion that the inclination distribution at the end of the Grand Tack is a bit over-excited. △ Less

Submitted 10 January, 2017; originally announced January 2017.

Comments: published on Icarus

Journal ref: Icarus, 272, 114 (2016)

Terrestrial Planet Formation from an Annulus

Authors: Kevin J. Walsh, Hal F. Levison

Abstract: It has been shown that some aspects of the terrestrial planets can be explained, particularly the Earth/Mars mass ratio, when they form from a truncated disk with an outer edge near 1.0 au (Hansen 2009). This has been previously modeled starting from an intermediate stage of growth utilizing pre-formed planetary embryos. We present simulations that were designed to test this idea by following the… ▽ More It has been shown that some aspects of the terrestrial planets can be explained, particularly the Earth/Mars mass ratio, when they form from a truncated disk with an outer edge near 1.0 au (Hansen 2009). This has been previously modeled starting from an intermediate stage of growth utilizing pre-formed planetary embryos. We present simulations that were designed to test this idea by following the growth process from km-sized objects located between 0.7 to 1.0 au up to terrestrial planets. The simulations explore initial conditions where the solids in the disk are planetesimals with radii initially between 3 and 300 km, alternately including effects from a dissipating gaseous solar nebula and collisional fragmentation. We use a new Lagrangian code known as LIPAD (Levison et al. 2012), which is a particle-based code that models the fragmentation, accretion and dynamical evolution of a large number of planetesimals, and can model the entire growth process from km-sizes up to planets. A suite of large (Mars mass) planetary embryos is complete in only 1 Myr, containing most of the system mass. A quiescent period then persists for 10-20 Myr characterized by slow diffusion of the orbits and continued accretion of the remaining planetesimals. This is interupted by an instability that leads to embryos crossing orbits and embyro-embryo impacts that eventually produce the final set of planets. While this evolution is different than that found in other works exploring an annulus, the final planetary systems are similar, with roughly the correct number of planets and good Mars-analogs. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: 23 pages, 15 figs

Journal ref: Astronomical Journal, 152, 3 (2016)

Yarkovsky V-shape identification of asteroid families

Authors: Bryce T. Bolin, Marco Delbo, Alessandro Morbidelli, Kevin J. Walsh

Abstract: There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity,… ▽ More There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 74 pages, 17 figures, accepted for publication in Icarus

Orbital Simulations on Deflecting Near-Earth Objects by Directed Energy

Authors: Qicheng Zhang, Kevin J. Walsh, Carl Melis, Gary B. Hughes, Philip M. Lubin

Abstract: Laser ablation of a Near-Earth Object (NEO) on a collision course with Earth produces a cloud of ejecta which exerts a thrust on the NEO, deflecting it from its original trajectory. Ablation may be performed from afar by illuminating an Earth-targeting asteroid or comet with a stand-off "DE- STAR" system consisting of a large phased-array laser in Earth orbit. Alternatively, a much smaller stand-o… ▽ More Laser ablation of a Near-Earth Object (NEO) on a collision course with Earth produces a cloud of ejecta which exerts a thrust on the NEO, deflecting it from its original trajectory. Ablation may be performed from afar by illuminating an Earth-targeting asteroid or comet with a stand-off "DE- STAR" system consisting of a large phased-array laser in Earth orbit. Alternatively, a much smaller stand-on "DE-STARLITE" system may travel alongside the target, slowly deflecting it from nearby over a long period. This paper presents orbital simulations comparing the effectiveness of both systems across a range of laser and NEO parameters. Simulated parameters include magnitude, duration and, for the stand-on system, direction of the thrust, as well as the type, size and orbital characteristics of the target NEO. These simulations indicate that deflection distance is approximately proportional to the magnitude of thrust and to the square of the duration of ablation, and is inversely proportional to the mass. Furthermore, deflection distance shows strong dependence on thrust direction with the optimal direction of thrust varying with the duration of laser activity. As one example, consider a typical 325 m asteroid: beginning 15 yr in advance, just 2 N of thrust from a ~20 kW stand-on DE-STARLITE system is sufficient to deflect the asteroid by 2 R_e. Numerous scenarios are discussed as is a practical implementation of such a system consistent with current launch vehicle capabilities. △ Less

Submitted 14 January, 2016; originally announced January 2016.

Comments: 10 pages, 11 figures. Accepted by PASP

Journal ref: Publ.Astron.Soc.Pac. 128 (2016) 045001

Formation and Evolution of Binary Asteroids

Authors: Kevin J. Walsh, Seth A. Jacobson

Abstract: 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… ▽ More 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 system's 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. △ Less

Submitted 22 June, 2015; originally announced June 2015.

Comments: 42 pages, 4 figures, contribution to the Asteroids 4 book

The Compositional Structure of the Asteroid Belt

Authors: Francesca E. DeMeo, Conel M. O'D. Alexander, Kevin J. Walsh, Clark R. Chapman, Richard P. Binzel

Abstract: The past decade has brought major improvements in large-scale asteroid discovery and characterization with over half a million known asteroids and over 100,000 with some measurement of physical characterization. This explosion of data has allowed us to create a new global picture of the Main Asteroid Belt. Put in context with meteorite measurements and dynamical models, a new and more complete pic… ▽ More The past decade has brought major improvements in large-scale asteroid discovery and characterization with over half a million known asteroids and over 100,000 with some measurement of physical characterization. This explosion of data has allowed us to create a new global picture of the Main Asteroid Belt. Put in context with meteorite measurements and dynamical models, a new and more complete picture of Solar System evolution has emerged. The question has changed from "What was the original compositional gradient of the Asteroid Belt?" to "What was the original compositional gradient of small bodies across the entire Solar System?" No longer is the leading theory that two belts of planetesimals are primordial, but instead those belts were formed and sculpted through evolutionary processes after Solar System formation. This article reviews the advancements on the fronts of asteroid compositional characterization, meteorite measurements, and dynamical theories in the context of the heliocentric distribution of asteroid compositions seen in the Main Belt today. This chapter also reviews the major outstanding questions relating to asteroid compositions and distributions and summarizes the progress and current state of understanding of these questions to form the big picture of the formation and evolution of asteroids in the Main Belt. Finally, we briefly review the relevance of asteroids and their compositions in their greater context within our Solar System and beyond. △ Less

Submitted 15 June, 2015; originally announced June 2015.

Comments: Accepted chapter in Asteroids IV in the Space Science Series to be published Fall 2015

Formation and evolution of Pluto's small satellites

Authors: Kevin J. Walsh, Harold F. Levison

Abstract: Pluto's system of 5 known satellites are in a puzzling orbital configuration. Each of the four small satellites are on low-eccentricity and low-inclination orbits situated near a mean motion resonance with the largest satellite Charon. The Pluto-Charon binary likely formed as a result of a giant impact and so the simplest explanation for the small satellites is that they accreted from debris of th… ▽ More Pluto's system of 5 known satellites are in a puzzling orbital configuration. Each of the four small satellites are on low-eccentricity and low-inclination orbits situated near a mean motion resonance with the largest satellite Charon. The Pluto-Charon binary likely formed as a result of a giant impact and so the simplest explanation for the small satellites is that they accreted from debris of that collision. The Pluto-Charon binary has evolved outward since its formation due to tidal forces, which drove them into their current doubly synchronous state. Meanwhile, leftover debris from the formation of Charon was not initially distant enough from Pluto-Charon to explain the orbits of the current small satellites. The outstanding problems of the system are the movement of debris outward and the small satellites location near mean motion resonances with Charon. This work explores the dynamical behavior of collisionally interacting debris orbiting the Pluto-Charon system. While this work specifically tests initial disk and ring configurations designed to mimic the aftermath of the disruption of satellites by heliocentric impactors, we generally find that collisional interactions can help move material outwards and keep otherwise unstable material dynamically bound to the Pluto-Charon system. These processes can produce rings of debris whose orbits evolve rapidly due to collisional processes, with increasing pericenters and decreasing semimajor axes. While these rings and disks of debris eventually build satellites significantly further out than the initial locations of a disrupted satellite, they do not show a strong preference for building satellites in or near mean motion resonances with Charon under a wide array of tested conditions. △ Less

Submitted 5 May, 2015; originally announced May 2015.

Comments: 18 pages, 8 figures. Accepted in AJ

Highly Siderophile Elements in the Earth's Mantle as a Clock for the Moon-forming Impact

Authors: Seth A. Jacobson, Alessandro Morbidelli, Sean N. Raymond, David P. O'Brien, Kevin J. Walsh, David C. Rubie

Abstract: According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the 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 Sola… ▽ More According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the 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 My) 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 it 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 since 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 My after condensation. △ Less

Submitted 6 April, 2015; originally announced April 2015.

Comments: Includes Supplemental Material

Journal ref: Nature, Volume 508, Issue 7494, pp. 84-87 (2014)

Earth and Terrestrial Planet Formation

Authors: Seth A. Jacobson, Kevin J. Walsh

Abstract: The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for th… ▽ More The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of planetary embryos. Many of the popular models for the later stages of planet formation are presented. The classical models with the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets. In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The late accreted mass samples the disk nearly evenly. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: Accepted for publication in Early Earth an AGU Monograph edited by James Badro and Michael J. Walter

Journal ref: In: The early Earth: accretion and differentiation (Eds. J. Badro & M. J. Walter) AGU/John Wiley & Sons, Inc., Geophysical Monograph Series 212 (2015) 49-70

The Dynamical Evolution of the Asteroid Belt

Authors: Alessandro Morbidelli, Kevin J. Walsh, David P. O'Brien, David A. Minton, William F. Bottke

Abstract: The asteroid belt is the leftover of the original planetesimal population in the inner solar system. However, currently the asteroids have orbits with all possible values of eccentricities and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal orbits should have been quasi-circular and almost co-planar. The total mass in the asteroid population is a small… ▽ More The asteroid belt is the leftover of the original planetesimal population in the inner solar system. However, currently the asteroids have orbits with all possible values of eccentricities and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal orbits should have been quasi-circular and almost co-planar. The total mass in the asteroid population is a small fraction of that existing primordially. Also, asteroids with different chemical/mineralogical properties are not ranked in an orderly manner with mean heliocentric distance as one could expect from the existence of a radial gradient of the temperature in the proto-planetary disk, but they are partially mixed. These properties show that the asteroid belt has been severely sculpted by one or a series of processes during its lifetime. This paper reviews the processes that have been proposed so far, discussing the properties that they explain and the problems that they are confronted with. Emphasis is paid to the interplay between the dynamical and the collisional evolution of the asteroid population, which allows the use of the size distribution to constrain the dynamical models. We divide the asteroid belt evolution into three phases. The first phase started during the lifetime of the gaseous proto-planetary disk, when the giant planets formed and presumably experienced large-scale migrations, and continued after the removal of the gas, during the build-up of the terrestrial planets. The second phase occurred after the removal of the gaseous proto-planetary disk and it became particularly lively for the asteroid belt when the giant planets suddenly changed their orbits, as a result of a mutual dynamical instability and the interaction with the trans-Neptunian planetesimal disk. The third phase covers the aftermath of the giant planet instability, until today. △ Less

Submitted 25 January, 2015; originally announced January 2015.

Comments: Accepted for publication in Asteroid IV

Water Delivery and Giant Impacts in the 'Grand Tack' Scenario

Authors: David P. O'Brien, Kevin J. Walsh, Alessandro Morbidelli, Sean N. Raymond, Avi M. Mandell

Abstract: A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanis… ▽ More A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk--a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies--the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as `primitive' bodies). We find here that the planets will accrete on order 1-2% of their total mass from primitive planetesimals scattered onto planet-crossing orbits during the formation of the planets. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. While the radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system, we find that the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event. △ Less

Submitted 11 July, 2014; originally announced July 2014.

Journal ref: Published Open Access in Icarus, Volume 239, pp. 74-84 (2014)

Excitation of the Orbital Inclination of Iapetus during Planetary Encounters

Authors: David Nesvorny, David Vokrouhlicky, Rogerio Deienno, Kevin J. Walsh

Abstract: Saturn's moon Iapetus has an orbit in a transition region where the Laplace surface is bending from the equator to the orbital plane of Saturn. The orbital inclination of Iapetus to the local Laplace plane is ~8 deg, which is unexpected, because the inclination should be ~0 if Iapetus formed from a circumplanetary disk on the Laplace surface. It thus appears that some process has pumped up Iapetus… ▽ More Saturn's moon Iapetus has an orbit in a transition region where the Laplace surface is bending from the equator to the orbital plane of Saturn. The orbital inclination of Iapetus to the local Laplace plane is ~8 deg, which is unexpected, because the inclination should be ~0 if Iapetus formed from a circumplanetary disk on the Laplace surface. It thus appears that some process has pumped up Iapetus's inclination while leaving its eccentricity near zero (e=0.03 at present). Here we examined the possibility that Iapetus's inclination was excited during the early solar system instability when encounters between Saturn and ice giants occurred. We found that the dynamical effects of planetary encounters on Iapetus's orbit sensitively depend on the distance of the few closest encounters. In four out of ten instability cases studied here, the orbital perturbations were too large to be plausible. In one case, Iapetus's orbit was practically unneffected. In the remaining five cases, the perturbations of Iapetus's inclination were adequate to explain its present value. In three of these cases, however, Iapetus's eccentricity was excited to 0.1-0.25, and it is not clear whether it could have been damped to its present value (0.03) by some subsequent process (e.g., tides and dynamical friction from captured irregular satellites do not seem to be strong enough). Our results therefore imply that only 2 out of 10 instability cases (~20%) can excite Iapetus's inclination to its present value (~30% of trials lead to >5 deg) while leaving its orbital eccentricity low. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: Astronomical Journal

Properties and evolution of NEO families created by tidal disruption at Earth

Authors: Eva Schunová, Robert Jedicke, Kevin J. Walsh, Mikael Granvik, Richard J. Wainscoat, Nader Haghighipour

Abstract: We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for… ▽ More We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\barτ_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\barτ_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population. △ Less

Submitted 16 May, 2014; originally announced May 2014.

Comments: 60 pages, 16 figures

MSC Class: 85

Terrestrial Planet Formation at Home and Abroad

Authors: Sean N. Raymond, Eiichiro Kokubo, Alessandro Morbidelli, Ryuji Morishima, Kevin J. Walsh

Abstract: We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the… ▽ More We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations. △ Less

Submitted 5 September, 2014; v1 submitted 5 December, 2013; originally announced December 2013.

Comments: Reduced to 24 pages, 9 figures, 2 tables. Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning. (Updated references)

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

Authors: Matthew M. Knight, Kevin J. Walsh

Abstract: 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 comet's response to such a close passage by the Sun is important for planning observational campaigns and for inferring ISON's physical properties. We present new numerical simulations and interpret them in context with the historical track reco… ▽ More 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 comet's response to such a close passage by the Sun is important for planning observational campaigns and for inferring ISON's 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 ISON's 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 ISON's 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. △ Less

Submitted 9 September, 2013; originally announced September 2013.

Comments: Accepted by ApJL; 12 pages of text (pre-print style), 1 figure, 1 table

Numerical simulations of granular dynamics. I. Hard-sphere discrete element method and tests

Authors: Derek C. Richardson, Kevin J. Walsh, Naomi Murdoch, Patrick Michel

Abstract: We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres… ▽ More We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate. △ Less

Submitted 11 June, 2013; originally announced June 2013.

Comments: 54 manuscript pages, 8 figures including 4 in colour (online version only)

Journal ref: Icarus 212 (2011) 427-437

Introducing the Eulalia and new Polana asteroid families: re-assessing primitive asteroid families in the inner Main Belt

Authors: Kevin J. Walsh, Marco Delb�, William F. Bottke, David Vokrouhlick�, Dante S. Lauretta

Abstract: The so-called Nysa-Polana complex of asteroids is a diverse and widespread group. It appears to be two overlapping families of different asteroid taxonomies: (44) Nysa is an E-type asteroid with the lowest number in the midst of a predominantly S-type cluster and (142) Polana is a B-type asteroid near the low-albedo B- and C-type cluster. Using the data from the Wide-field Infrared Survey Explor… ▽ More The so-called Nysa-Polana complex of asteroids is a diverse and widespread group. It appears to be two overlapping families of different asteroid taxonomies: (44) Nysa is an E-type asteroid with the lowest number in the midst of a predominantly S-type cluster and (142) Polana is a B-type asteroid near the low-albedo B- and C-type cluster. Using the data from the Wide-field Infrared Survey Explorer (WISE) mission we have re-analyzed the region around the Nysa-Polana complex in the inner Main Belt, focusing on the low-albedo population. (142) Polana does not appear to be a member of the family of low-albedo asteroids in the Nysa-Polana complex. Rather, the largest is asteroid (495) Eulalia. This asteroid has never before been linked to this complex for an important dynamical reason: it currently has a proper eccentricity slightly below the range of most of the family members. However, its orbit is very close to the 3:1 mean motion resonance with Jupiter and is in a weak secular resonance. We show that its osculating eccentricity varies widely on short timescales and the averaged value diffuses over long timescales. The diffusive orbit, low-albedo, taxonomic similarity and semimajor axis strongly suggests that despite its current proper eccentricity, (495) Eulalia could have recently been at an orbit very central to the family. Hierarchical Clustering Method tests confirm that at an eccentricity of e=0.15, (495) Eulalia could be the parent of the family. The ``Eulalia family'' was formed between 900--1500 Myr ago, and likely resulted from the breakup of a 100--160 km parent body. There is also compelling evidence for an older and more widespread primitive family in the same region of the asteroid belt parented by asteroid (142) Polana. △ Less

Submitted 13 May, 2013; originally announced May 2013.

Comments: 49 pages, 22 figures, Accepted to Icarus

Building Terrestrial Planets

Authors: Alessandro Morbidelli, Jonathan I. Lunine, David P. O`brien, Sean N. Raymond, Kevin J. Walsh

Abstract: This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given… ▽ More This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order 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 planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth. △ Less

Submitted 23 August, 2012; originally announced August 2012.

Journal ref: Annual Reviews of Earth and Planetary Sciences, vol. 40, issue 1, pp. 251-275 (2012)

Physical characterization and origin of binary near-Earth asteroid (175706) 1996 FG3

Authors: Kevin J. Walsh, Marco Delbo, Michael Mueller, Richard P. Binzel, Francesca E. DeMeo

Abstract: The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-DeltaV 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 ob… ▽ More The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-DeltaV 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 (IR) wavelengths and with near-IR reflection spectroscopy. We derive an area-equivalent system diameter of 1.90 \pm 0.28 km (corresponding to approximate component diameters of 1.83 km and 0.51 km, respectively) and a geometric albedo of 0.039 \pm 0.012. 1996 FG3 was previously classified as a C-type asteroid, though the combined 0.4--2.5 micron spectrum with thermal correction indicates classification as B-type; both are consistent with the low measured albedo. Dynamical studies show that 1996 FG3 has 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-DeltaV B-type NEA, (101955) 1999 RQ36. A similar origin for 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and Mars 5:9 mean motion resonances rather than the nu-6 resonance, and we find this to be a low probability, but possible, origin. △ Less

Submitted 21 March, 2012; originally announced March 2012.

Comments: Published in ApJ

Journal ref: Walsh et al. 2011, ApJ, 748, 104

A low mass for Mars from Jupiter's early gas-driven migration

Authors: Kevin J. Walsh, Alessando Morbidelli, Sean N. Raymond, David P. O'Brien, Avi M. Mandell

Abstract: Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens &… ▽ More Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens & Nelson 2008). The terrestrial planets finished accreting much later (Klein et al. 2009), and their characteristics, including Mars' small mass, are best reproduced by starting from a planetesimal disk with an outer edge at about one astronomical unit from the Sun (Wetherill 1978, Hansen 2009) (1 AU is the Earth-Sun distance). Here we report simulations of the early Solar System that show how the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, lead to a planetesimal disk truncated at 1 AU; the terrestrial planets then form from this disk over the next 30-50 million years, with an Earth/Mars mass ratio consistent with observations. Scattering by Jupiter initially empties but then repopulates the asteroid belt, with inner-belt bodies originating between 1 and 3 AU and outer-belt bodies originating between and beyond the giant planets. This explains the significant compositional differences across the asteroid belt. The key aspect missing from previous models of terrestrial planet formation is the substantial radial migration of the giant planets, which suggests that their behaviour is more similar to that inferred for extrasolar planets than previously thought. △ Less

Submitted 24 January, 2012; originally announced January 2012.

Comments: 12 pages, 4 figures + Supplementary Material 46 pages, 10 figures

Journal ref: Nature 475, 206-209 (14 July 2011)

Ridge Formation and De-Spinning of Iapetus via an Impact-Generated Satellite

Authors: H. F. Levison, K. J. Walsh, A. C. Barr, L. Dones

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 ca… ▽ More 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. △ Less

Submitted 9 May, 2011; originally announced May 2011.

Comments: 19 pages, 3 figures; Icarus, in press

Journal ref: Icarus v. 214, p. 773-778, 2011

Stability analysis of the Martian obliquity during the Noachian era

Authors: Ramon Brasser, Kevin J. Walsh

Abstract: We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities… ▽ More We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities below 30$^\circ$ and above 60$^\circ$; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets. △ Less

Submitted 4 February, 2011; originally announced February 2011.

The Effect of an Early Planetesimal-Driven Migration of the Giant Planets on Terrestrial Planet Formation

Authors: Kevin J. Walsh, Alessandro Morbidelli

Abstract: 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… ▽ More 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. △ Less

Submitted 19 January, 2011; originally announced January 2011.

Comments: Accepted for publication in Astronomy and Astrophysics

Collisional and Rotational Disruption of Asteroids

Authors: Kevin J. Walsh, Patrick Michel, Derek C. Richardson

Abstract: Asteroids are leftover pieces from the era of planet formation that help us understand conditions in the early Solar System. Unlike larger planetary bodies that were subject to global thermal modification during and subsequent to their formation, these small bodies have kept at least some unmodified primordial material from the solar nebula. However, the structural properties of asteroids have b… ▽ More Asteroids are leftover pieces from the era of planet formation that help us understand conditions in the early Solar System. Unlike larger planetary bodies that were subject to global thermal modification during and subsequent to their formation, these small bodies have kept at least some unmodified primordial material from the solar nebula. However, the structural properties of asteroids have been modified considerably since their formation. Thus, we can find among them a great variety of physical configurations and dynamical histories. In fact, with only a few possible exceptions, all asteroids have been modified or completely disrupted many times during the age of the Solar System. This picture is supported by data from space mission encounters with asteroids that show much diversity of shape, bulk density, surface morphology, and other features. Moreover, the gravitational attraction of these bodies is so small that some physical processes occur in a manner far removed from our common experience on Earth. Thus, each visit to a small body has generated as many questions as it has answered. In this review we discuss the current state of research into asteroid disruption processes, focusing on collisional and rotational mechanisms. We find that recent advances in modeling catastrophic disruption by collisions have provided important insights into asteroid internal structures and a deeper understanding of asteroid families. Rotational disruption, by tidal encounters or thermal effects, is responsible for altering many smaller asteroids, and is at the origin of many binary asteroids and oddly shaped bodies. △ Less

Submitted 23 June, 2009; originally announced June 2009.

Comments: Accepted for publication to Advanced Science Letters, Special Issue on Computational Astrophysics, edited by Lucio Mayer