A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit period. They may have been subsequently dissociated by impacts and planetary encounters with at least one binary component retaining its slow rotation. Surviving binaries on Trojan orbits would continue to evolve by tides and spin-changing impacts over 4.5 Gyr. To explain the observed fraction of slow rotators, we find that at least ~15-20% of outer disk bodies with diameters 15<D<50 km would have to form as equal-size binaries with 12<ab/R<30, where ab is the binary semimajor axis and R=D/2. The mechanism proposed here could also explain very slow rotators found in other small body populations.

We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

(modified from published version) Twelve years ago the Catalina Sky Survey discovered Earth's first known natural geocentric object other than the Moon, a few-meter diameter asteroid designated 2006 RH120. Despite significant improvements in ground-based telescope and detector technology in the past decade the asteroid surveys have not discovered another temporarily-captured orbiter (TCO; colloquially known as minimoons). Within a few years the Large Synoptic Survey Telescope (LSST) will either begin to regularly detect TCOs or force a re-analysis of the creation and dynamical evolution of small asteroids in the inner solar system. The first studies of the provenance, properties, and dynamics of Earth's minimoons suggested that there should be a steady state population with about one 1- to 2-meter diameter captured objects at any time. That model was then improved and extended to include the population of temporarily-captured flybys (TCFs), objects that fail to make an entire revolution around Earth while energetically bound to the Earth-Moon system. Several different techniques for discovering TCOs have been considered but their small diameters, proximity, and rapid motion make them challenging targets for existing ground-based optical, meteor, and radar surveys. We expect that if the TCO population is confirmed, and new objects are frequently discovered, they can provide new opportunities for 1) studying the dynamics of the Earth-Moon system, 2) testing models of the production and dynamical evolution of small asteroids from the asteroid belt, 3) rapid and frequent low delta-v missions to multiple minimoons, and 4) evaluating in-situ resource utilization techniques on asteroidal material. Here we review the past decade of minimoon studies in preparation for capitalizing on the scientific and commercial opportunities of TCOs in the first decade of LSST operations.

The orbital distribution of trans-Neptunian objects provides strong evidence for the radial migration of Neptune. The outer planets' orbits are thought to have become unstable during the early stages with Jupiter having scattering encounters with a Neptune-class planet. As a consequence, Jupiter jumped inward by a fraction of an au, as required from inner solar system constraints, and obtained its current orbital eccentricity. The timing of these events is often linked to the lunar Late Heavy Bombardment that ended ~700 Myr after the dispersal of the protosolar nebula (t_0). Here we show instead that planetary migration started shortly after t_0. Such early migration is inferred from the survival of the Patroclus-Menoetius binary Jupiter Trojan. The binary formed at t <~ t_0 within a massive planetesimal disk once located beyond Neptune. The longer the binary stayed in the disk, the greater the likelihood that collisions would strip its components from one another. The simulations of its survival indicate that the disk had to have been dispersed by migrating planets within <100 Myr of t_0. This constraint implies that the planetary migration is unrelated to the formation of the youngest lunar basins.

The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

On April 13, 2029, near-Earth asteroid (NEA) (99942) Apophis will pass at a distance of 6 Earth radii from Earth. This event will provide researchers with a unique opportunity to study the effects of tidal forces experienced by an asteroid during a close encounter with a terrestrial planet. Binzel et al. (2010) predicted that close flybys of terrestrial planets by NEAs would cause resurfacing of their regolith due to seismic shaking. In this work we present the best pre-encounter near-infrared spectra of Apophis obtained so far. These new data were obtained during the 2013 apparition using the NASA Infrared Telescope Facility (IRTF). We found that our spectral data is consistent with previous observations by Binzel et al. (2009) but with a much higher signal-to-noise ratio. Spectral band parameters were extracted from the spectra and were used to determine the composition of the asteroid. Using a naive Bayes classifier, we computed the likelihood of Apophis being an LL chondrite to be >99% based on mol% of Fa vs. Fs. Using the same method, we estimated a probability of 89% for Apophis being an LL chondrite based on ol/(ol+px) and Fs. The results from the dynamical model indicate that the most likely source region for Apophis is the nu6 resonance in the inner main belt. Data presented in this study (especially Band I depth) could serve as a baseline to verify seismic shaking during the 2029 encounter.

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.

The impact cratering record of the Moon and the terrestrial planets provides important clues about the formation and evolution of the Solar System. Especially intriguing is the epoch 3.8-3.9 Gyr ago (Ga), known as the Late Heavy Bombardment (LHB), when the youngest lunar basins such as Imbrium and Orientale formed. The LHB was suggested to originate from a slowly declining impactor flux or from a late dynamical instability. Here we develop a model for the historical flux of large asteroid impacts and discuss how it depends on various parameters, including the time and nature of the planetary migration/instability. We find that the asteroid impact flux dropped by 1 to 2 orders of magnitude during the first 1 Gyr and remained relatively unchanged over the last 3 Gyr. The early impacts were produced by asteroids whose orbits became excited during the planetary migration/instability, and by those originating from the inner extension of the main belt (E-belt; semimajor axis 1.6<a<2.1 au). The profiles obtained for the early and late versions of the planetary instability initially differ, but end up being similar after ~3 Ga. Thus, the time of the instability can only be determined by considering the cratering and other constraints during the first ~1.5 Gyr of the Solar System history. Our absolute calibration of the impact flux indicates that asteroids were probably not responsible for the LHB, independently of whether the instability happened early or late, because the calibrated flux is not large enough to explain Imbrium/Orientale and a significant share of large lunar craters. Comets and leftovers of the terrestrial planet formation provided additional, and probably dominant source of impacts during early epochs.

Small near-Earth asteroids (>20 meters) are interesting because they are progenitors for meteorites in our terrestrial collection. Crucial to our understanding of the effectiveness of our atmosphere in filtering low-strength impactors is the physical characteristics of these small near-Earth asteroids (NEAs). In the past, characterization of small NEAs has been a challenge because of the difficulty in detecting them prior to close Earth flyby. In this study we physically characterized the 2-meter diameter near-Earth asteroid 2015 TC25 using ground-based optical, near-infrared and radar assets during a close flyby of the Earth (distance 69,000 miles) in Oct. 2015. Our observations suggest that its surface composition is similar to aubrites, a rare class of high albedo differentiated meteorites. Aubrites make up only 0.14 % of all know meteorites in our terrestrial meteorite collection. 2015 TC25 is also a very fast rotator with a rotation period of 133 seconds. We compared spectral and dynamical properties of 2015 TC25 and found the best candidate source body in the inner main belt to be the 70-km diameter E-type asteroid (44) Nysa. We attribute difference in spectral slope between the two objects to the lack of regolith on the surface of 2015 TC25. Using the albedo of E-type asteroids (50-60%) we refine the diameter of 2015 TC25 to 2-meters making it one of the smallest NEA ever to be characterized.

Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants. This model envisions that the cores of the giant planets formed from 100 to 1000 km bodies that directly accreted a population of pebbles --- sub-meter sized objects that slowly grew in the protoplanetary disk. Here we apply this model to the terrestrial planet region and find that it can reproduce the basic structure of the inner Solar System, including a small Mars and a low-mass asteroid belt. Our models show that for an initial population of planetesimals with sizes similar to those of the main belt asteroids, VSPA becomes inefficient beyond $\sim\!$1.5 AU. As a result, Mars's growth is stunted and nothing large in the asteroid belt can accumulate.

We perform a search for dormant comets, asteroidal objects of cometary origin, in the near-Earth asteroid (NEA) population based on dynamical and physical considerations. Our study is based on albedos derived within the ExploreNEOs program and is extended by adding data from NEOWISE and the Akari asteroid catalog. We use a statistical approach to identify asteroids on orbits that resemble those of short-period near-Earth comets using the Tisserand parameter with respect to Jupiter, the aphelion distance, and the minimum orbital intersection distance with respect to Jupiter. From the sample of NEAs on comet-like orbits, we select those with a geometric albedo $p_V \leq 0.064$ as dormant comet candidates, and find that only $\sim$50% of NEAs on comet-like orbits also have comet-like albedos. We identify a total of 23 NEAs from our sample that are likely to be dormant short-period near-Earth comets and, based on a de-biasing procedure applied to the cryogenic NEOWISE survey, estimate both magnitude-limited and size-limited fractions of the NEA population that are dormant short-period comets. We find that 0.3-3.3% of the NEA population with $H \leq 21$, and $9^{+2}_{-5}$% of the population with diameters $d \geq 1$ km, are dormant short-period near-Earth comets.

We describe the fall of Annama meteorite occurred in the remote Kola Peninsula (Russia) close to Finnish border on April 19, 2014 (local time). The fireball was instrumentally observed by the Finnish Fireball Network. From these observations the strewnfield was computed and two first meteorites were found only a few hundred meters from the predicted landing site on May 29th and May 30th 2014, so that the meteorite (an H4-5 chondrite) experienced only minimal terrestrial alteration. The accuracy of the observations allowed a precise geocentric radiant to be obtained, and the heliocentric orbit for the progenitor meteoroid to be calculated. Backward integrations of the orbits of selected near-Earth asteroids and the Annama meteoroid showed that they rapidly diverged so that the Annama meteorites are unlikely related to them. The only exception seems to be the recently discovered 2014UR116 that shows a plausible dynamic relationship. Instead, analysis of the heliocentric orbit of the meteoroid suggests that the delivery of Annama onto an Earth-crossing Apollo type orbit occurred via the 4:1 mean motion resonance with Jupiter or the nu6 secular resonance, dynamic mechanisms that are responsible for delivering to Earth most meteorites studied so far.

We explored the statistical and compositional link between Chelyabinsk meteoroid and potentially hazardous asteroid (86039) 1999 NC43 to investigate their proposed relation proposed by Borovi\v{c}ka et al. (2013). Using detailed computation we confirm that the orbit of the Chelyabinsk impactor is anomalously close to 1999 NC43. We find about (1-3) x 10-4 likelihood of that to happen by chance. Taking the standpoint that the Chelyabinsk impactor indeed separated from 1999 NC43 by a cratering or rotational fission event, we run a forward probability calculation, which is an independent statistical test. However, we find this scenario is unlikely at the about (10-3 -10-2) level. We also verified compositional link between Chelyabinska and 1999NC43. Mineralogical analysis of Chelyabinsk (LL chondrite) and (8) Flora (the largest member of the presumed LL chondrite parent family) shows that their olivine and pyroxene chemistries are similar to LL chondrites. Similar analysis of 1999 NC43 shows that its olivine and pyroxene chemistries are more similar to L chondrites than LL chondrites (like Chelyabinsk). We also took photometric observations of 1999 NC43 over 54 nights during two apparitions (2000, 2014). The lightcurve of 1999 NC43 resembles simulated lightcurves of tumblers in Short-Axis Mode with the mean wobbling angle 20-30 deg. While, a mechanism of the non-principal axis rotation excitation is unclear, we can rule out the formation of asteroid in disruption of its parent body as a plausible cause, as it is unlikely that the rotation of an asteroid fragment from catastrophic disruption would be nearly completely halted. Considering all these facts, we find the proposed link between the Chelyabinsk meteoroid and the asteroid 1999 NC43 to be unlikely.

The Yarkovsky effect describes a small but significant force that affects the orbital motion of meteoroids and asteroids smaller than $30-40$ kilometers in diameter. It is caused by sunlight; when these bodies heat up in the Sun, they eventually re-radiate the energy away in the thermal waveband, which in turn creates a tiny thrust. This recoil acceleration is much weaker than solar and planetary gravitational forces, but it can produce measurable orbital changes over decades and substantial orbital effects over millions to billions of years. The same physical phenomenon also creates a thermal torque that, complemented by a torque produced by scattered sunlight, can modify the rotation rates and obliquities of small bodies as well. This rotational variant has been coined the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. During the past decade or so, the Yarkovsky and YORP effects have been used to explore and potentially resolve a number of unsolved mysteries in planetary science dealing with small bodies. Here we review the main results to date, and preview the goals for future work.

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.

The Design Reference Asteroid (DRA) is a compilation of all that is known about the OSIRIS-REx mission target, asteroid (101955) Bennu. It contains our best knowledge of the properties of Bennu based on an extensive observational campaign that began shortly after its discovery, and has been used to inform mission plan development and flight system design. The DRA will also be compared with post-encounter science results to determine the accuracy of our Earth-based characterization efforts. The extensive observations of Bennu in 1999 has made it one of the best-characterized near-Earth asteroids. Many physical parameters are well determined, and span a number of categories: Orbital, Bulk, Rotational, Radar, Photometric, Spectroscopic, Thermal, Surface Analog, and Environment Properties. Some results described in the DRA have been published in peer-reviewed journals while others have been reviewed by OSIRIS-REx Science Team members and/or external reviewers. Some data, such as Surface Analog Properties, are based on our best knowledge of asteroid surfaces, in particular those of asteroids Eros and Itokawa. This public release of the OSIRIS-REx Design Reference Asteroid is a annotated version of the internal OSIRIS-REx document OREX-DOCS-04.00-00002, Rev 9 (accepted by the OSIRIS-REx project on 2014-April-14). The supplemental data products that accompany the official OSIRIS-REx version of the DRA are not included in this release. We are making this document available as a service to future mission planners in the hope that it will inform their efforts.

Vesta has a complex cratering history, with ancient terrains as well as recent large impacts that have led to regional resurfacing. Crater counts can help constrain the relative ages of different units on Vesta's surface, but converting those crater counts to absolute ages requires a chronology function. We present a cratering chronology based on the best current models for the dynamical evolution of asteroid belt, and calibrate it to Vesta using the record of large craters on its surface. While uncertainties remain, our chronology function is broadly consistent with an ancient surface of Vesta as well as other constraints such as the bombardment history of the rest of the inner Solar System and the Ar-Ar age distribution of howardite, eucrite and diogenite (HED) meteorites from Vesta.

We investigated the spectral and compositional properties of Chelyabinsk meteorite to identify its possible parent body in the main asteroid belt. Our analysis shows that the meteorite contains two spectrally distinct but compositionally indistinguishable components of LL5 chondrite and shock blackened/impact melt material. Our X-ray diffraction analysis confirms that the two lithologies of the Chelyabinsk meteorite are extremely similar in modal mineralogy. The meteorite is compositionally similar to LL chondrite and its most probable parent asteroid in the main belt is a member of the Flora family. Intimate mixture of LL5 chondrite and shock blackened/impact melt material from Chelyabinsk provides a spectral match with (8) Flora, the largest asteroid in the Flora family. The Baptistina family and Flora family overlap each other in dynamical space. Mineralogical analysis of (298) Baptistina and 9 small family members shows that their surface compositions are similar to LL chondrites, although their absorption bands are subdued and albedos lower when compared to typical S-type asteroids. A range of intimate mixtures of LL5 chondrite and shock blackened/impact melt material from Chelyabinsk provides spectral matches for all these BAF members. We suggest that the presence of a significant shock/impact melt component in the surface regolith of BAF members could be the cause of lower albedo and subdued absorption bands. The parent asteroid of BAF was either a member of the Flora family or had the same basic composition as the Floras (LL Chondrite). The shock pressures produced during the impact event generated enough impact melt or shock blackening to alter the spectral properties of BAF, but keep the BAF composition largely unchanged.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

Olivine-dominated asteroids are a rare type of objects formed either in nebular processes or through magmatic differentiation. The analysis of meteorite samples suggest that at least 100 parent bodies in the main belt experienced partial or complete melting and differentiation before being disrupted. However, only a few olivine-dominated asteroids, representative of the mantle of disrupted differentiated bodies, are known to exist. Due to the paucity of these objects in the main belt their origin and evolution have been a matter of great debate over the years. In this work we present a detailed mineralogical analysis of twelve olivine-dominated asteroids. Within our sample we distinguish two classes, one that we call monomineralic-olivine asteroids and another referred to as olivine-rich asteroids. For the monomineralic-olivine asteroids the olivine chemistry was found to range from ~ Fo49 to Fo70, consistent with the values measured for brachinites and R chondrites. In the case of the olivine-rich asteroids we determined their olivine and low-Ca pyroxene abundance using a new set of spectral calibrations derived from the analysis of R chondrites spectra. We found that the olivine abundance for these asteroids varies from 0.68 to 0.93, while the fraction of low-Ca pyroxene to total pyroxene ranges from 0.6 to 0.9. A search for dynamical connections between the olivine-dominated asteroids and asteroid families found no genetic link (of the type core-mantel-crust) between these objects.

We propose a new Kepler mission, called NEOKepler, that would survey near Earth's orbit to identify potentially hazardous objects (PHOs). To understand its surveying power, Kepler's large field of view produces an 'etendue' (A*Omega) that is 4.5 times larger than the best survey telescope currently in operation. In this paper, we investigate the feasibility of NEOKepler using a double "fence post" survey pattern that efficiently detects PHOs. In a simulated 12-month survey, we estimate that NEOKepler would detect ~150 new NEOs with absolute magnitudes of less than 21.5, ~50 of which would be new PHOs. This would increase the annual PHO discovery rate by at least 50% and improve upon our goal of discovering 90% of PHOs by the end of 2020. Due to its heliocentric orbit, Kepler would also be sensitive to objects inside Earth's orbit, discovering more objects in its first year than are currently known to exist. Understanding this undersampled sub-population of NEOs will reveal new insights into the actual PHO distribution by further constraining current NEO models. As an alternative science goal, NEOKepler could employ a different observing strategy to discover suitable targets for NASA's Asteroid Redirect Mission.

The NASA Dawn mission has extensively examined the surface of asteroid Vesta, the second most massive body in the main belt. The high quality of the gathered data provides us with an unique opportunity to determine the surface and internal properties of one of the most important and intriguing main belt asteroids (MBAs). In this paper, we focus on the size frequency distributions (SFDs) of sub-kilometer impact craters observed at high spatial resolution on several selected young terrains on Vesta. These small crater populations offer an excellent opportunity to determine the nature of their asteroidal precursors (namely MBAs) at sizes that are not directly observable from ground-based telescopes (i.e., below ~100 m diameter). Moreover, unlike many other MBA surfaces observed by spacecraft thus far, the young terrains examined had crater spatial densities that were far from empirical saturation. Overall, we find that the cumulative power-law index (slope) of small crater SFDs on Vesta is fairly consistent with predictions derived from current collisional and dynamical models down to a projectile size of ~10 m diameter (Bottke et al., 2005a,b). The shape of the impactor SFD for small projectile sizes does not appear to have changed over the last several billions of years, and an argument can be made that the absolute number of small MBAs has remained roughly constant (within a factor of 2) over the same time period. The apparent steady state nature of the main belt population potentially provides us with a set of intriguing constraints that can be used to glean insights into the physical evolution of individual MBAs as well as the main belt as an ensemble.

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.

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.

[Abridged] The Moon is a unique platform from and on which to conduct astrophysical measurements. The Lunar University Network for Astrophysics Research (LUNAR) and the Center for Lunar Origins and Evolution (CLOE) teams within the NASA Lunar Science Institute (NLSI) are illustrating how the Moon can be used as a platform to advance important goals in astrophysics. Of relevance to Astrophysics and aligned with NASA strategic goals, all three of the primary research themes articulated by New Worlds, New Horizons in Astronomy & Astrophysics are being addressed by LUNAR and CLOE, namely Probing Cosmic Dawn, Understanding New Worlds, and Physics of the Universe

We revisit the early evolution of the Moon's bombardment. Our work combines modeling (based on plausible projectile sources and their dynamical decay rates) with constraints from the lunar crater record, radiometric ages of the youngest lunar basins, and the abundance of highly siderophile elements in the lunar crust and mantle. We deduce that the evolution of the impact flux did not decline exponentially over the first billion years of lunar history, but also there was no prominent and "narrow" impact spike some 3.9 Gy ago, unlike that typically envisioned in the lunar cataclysm scenario. Instead, we show the timeline of the lunar bombardment has a sawtooth-like profile, with an uptick in the impact flux near 4.1 Gy ago. The impact flux at the beginning of this weaker cataclysm was 5-10 times higher than the immediately preceding period. The Nectaris basin should have been one of the first basins formed at the sawtooth. We predict the bombardment rate since about 4.1Gy ago declined slowly and adhered relatively close to classic crater chronology models (Neukum and Ivanov (1994)). Overall we expect that the sawtooth event accounted for about 1/4 of the total bombardment suffered by the Moon since its formation. Consequently, considering that about 12-14 basins formed during the sawtooth event, we expect that the net number of basins formed on the Moon was about 45-50. From our expected bombardment timeline, we derived a new and improved lunar chronology suitable for use on Pre-Nectarian surface units. According to this chronology, a significant portion of the oldest lunar cratered terrains has an age of 4.38-4.42 Gyr. Moreover, the largest lunar basin, South Pole Aitken, is older than 4.3Gy, and therefore was not produced during the lunar cataclysm.

NASA's Dawn spacecraft observations of asteroid (4) Vesta reveal a surface with the highest albedo and color variation of any asteroid we have observed so far. Terrains rich in low albedo dark material (DM) have been identified using Dawn Framing Camera (FC) 0.75 {\mu}m filter images in several geologic settings: associated with impact craters (in the ejecta blanket material and/or on the crater walls and rims); as flow-like deposits or rays commonly associated with topographic highs; and as dark spots (likely secondary impacts) nearby impact craters. This DM could be a relic of ancient volcanic activity or exogenic in origin. We report that the majority of the spectra of DM are similar to carbonaceous chondrite meteorites mixed with materials indigenous to Vesta. Using high-resolution seven color images we compared DM color properties (albedo, band depth) with laboratory measurements of possible analog materials. Band depth and albedo of DM are identical to those of carbonaceous chondrite xenolith-rich howardite Mt. Pratt (PRA) 04401. Laboratory mixtures of Murchison CM2 carbonaceous chondrite and basaltic eucrite Millbillillie also show band depth and albedo affinity to DM. Modeling of carbonaceous chondrite abundance in DM (1-6 vol%) is consistent with howardite meteorites. We find no evidence for large-scale volcanism (exposed dikes/pyroclastic falls) as the source of DM. Our modeling efforts using impact crater scaling laws and numerical models of ejecta reaccretion suggest the delivery and emplacement of this DM on Vesta during the formation of the ~400 km Veneneia basin by a low-velocity (<2 km/sec) carbonaceous impactor. This discovery is important because it strengthens the long-held idea that primitive bodies are the source of carbon and probably volatiles in the early Solar System.

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (>10^5 yr at 1 AU) than postulated in the standard collisional models (10^4 yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5)x10^11 km^2 and 4x10^19 g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be 10^4-10^5 kg/s. The input is up to 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate.

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

The size distribution in the cold classical Kuiper belt can be approximated by two idealized power laws: one with steep slope for radii R>R* and one with shallow slope for R<R*, where R*~25-50 km. Previous works suggested that the SFD roll-over at R* can be the result of extensive collisional grinding in the Kuiper belt that led to the catastrophic disruption of most bodies with R<R*. Here we use a new code to test the effect of collisions in the Kuiper belt. We find that the observed roll-over could indeed be explained by collisional grinding provided that the initial mass in large bodies was much larger than the one in the present Kuiper belt, and was dynamically depleted. In addition to the size distribution changes, our code also tracks the effects of collisions on binary systems. We find that it is generally easier to dissolve wide binary systems, such as the ones existing in the cold Kuiper belt today, than to catastrophically disrupt objects with R~R*. Thus, the binary survival sets important limits on the extent of collisional grinding in the Kuiper belt. We find that the extensive collisional grinding required to produce the SFD roll-over at R* would imply a strong gradient of the binary fraction with R and separation, because it is generally easier to dissolve binaries with small components and/or those with wide orbits. The expected binary fraction for R<R* is <0.1. The present observational data do not show such a gradient. Instead, they suggest a large binary fraction of ~0.4 for R=30-40 km. This may indicate that the roll-over was not produced by disruptive collisions, but is instead a fossil remnant of the KBO formation process.

Space missions to NEOs are being planned at all major space agencies, and recently a manned mission to an NEO was announced as a NASA goal. Efforts to find and select suitable targets (plus backup targets) are severely hampered by our lack of knowledge of the physical properties of dynamically favorable NEOs. In particular, current mission scenarios tend to favor primitive low-albedo objects. For the vast majority of NEOs the albedo is unknown. Here we report new constraints on the size and albedo of 65 NEOs with rendezvous deltaV < 7 km/s. Our results are based on thermal-IR flux data obtained in the framework of our ongoing (2009--2011) ExploreNEOs survey (Trilling et al. 2010) using NASA's "Warm Spitzer" space telescope. As of 2010 July 14, we have results for 293 objects in hand (including the 65 low-deltaV NEOs presented here); before the end of 2011 we expect to have measured the size and albedo of ~ 700 NEOs (including probably ~ 160 low-deltaV NEOs). While there are reasons to believe that primitive volatile-rich materials are universally low in albedo, the converse need not be true: the orbital evolution of some dark objects likely has caused them to lose their volatiles by coming too close to the Sun. For all our targets, we give the closest perihelion distance they are likely to have reached (using orbital integrations from Marchi et al. 2009) and corresponding upper limits on the past surface temperature. Low-deltaV objects for which both albedo and thermal history may suggest a primitive composition include (162998) 2001 SK162, (68372) 2001 PM9, and (100085) 1992 UY4.

We have begun the ExploreNEOs project in which we observe some 700 Near Earth Objects (NEOs) at 3.6 and 4.5 microns with the Spitzer Space Telescope in its Warm Spitzer mode. From these measurements and catalog optical photometry we derive albedos and diameters of the observed targets. The overall goal of our ExploreNEOs program is to study the history of near-Earth space by deriving the physical properties of a large number of NEOs. In this paper we describe both the scientific and technical construction of our ExploreNEOs program. We present our observational, photometric, and thermal modeling techniques. We present results from the first 101 targets observed in this program. We find that the distribution of albedos in this first sample is quite broad, probably indicating a wide range of compositions within the NEO population. Many objects smaller than one kilometer have high albedos (>0.35), but few objects larger than one kilometer have high albedos. This result is consistent with the idea that these larger objects are collisionally older, and therefore possess surfaces that are more space weathered and therefore darker, or are not subject to other surface rejuvenating events as frequently as smaller NEOs.

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these `space weathering' processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r*, at which planetary encounters can reset the surface. We found that t_sw ~ 1e6 yr (at 1 AU) and r* ~ 5 R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 1e5 yr would require that r* > 20 R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 1e7 yr. We found that t_sw ~ q^2, where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.

We have searched Microvariability and Oscillations of STars (MOST) satellite photometry obtained in 2004, 2005, and 2007 of the solar-type star HD 209458 for Trojan asteroid swarms dynamically coupled with the system's transiting "hot Jupiter" HD 209458b. Observations of the presence and nature of asteroids around other stars would provide unique constraints on migration models of exoplanetary systems. Our results set an upper limit on the optical depth of Trojans in the HD 209458 system that can be used to guide current and future searches of similar systems by upcoming missions. Using cross-correlation methods with artificial signals implanted in the data, we find that our detection limit corresponds to a relative Trojan transit depth of 1\times10-4, equivalent to ~1 lunar mass of asteroids, assuming power-law Trojan size distributions similar to Jupiter's Trojans in our solar system. We confirm with dynamical interpretations that some asteroids could have migrated inward with the planet to its current orbit at 0.045 AU, and that the Yarkovsky effect is ineffective at eliminating objects of > 1 m in size. However, using numerical models of collisional evolution we find that, due to high relative speeds in this confined Trojan environment, collisions destroy the vast majority of the asteroids in <10 Myr. Our modeling indicates that the best candidates to search for exoTrojan swarms in 1:1 mean resonance orbits with "hot Jupiters" are young systems (ages of about 1 Myr or less). Years of Kepler satellite monitoring of such a system could detect an asteroid swarm with a predicted transit depth of 3\times10-7.

Our goal is to understand primary accretion of the first planetesimals. The primitive meteorite record suggests that sizeable planetesimals formed in the asteroid belt over a period longer than a million years, each composed entirely of an unusual, but homogeneous, mixture of mm-size particles. We sketch a scenario in which primary accretion of 10-100km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted mm-size particles (generically "chondrules"). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100km diameter, which has been identified as a "fossil" property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer solar system planetesimals, which may also have a preferred size in the 100km diameter range. We estimate formation rates of planetesimals and assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For nebula parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the "fossil" sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion (the model naturally helps explain the peculiar size distribution of chondrules within such objects). The scenario also produces 10-100km diameter primary KBOs. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between keplerian and pressure-supported orbital velocities, than "canonical" models of the solar nebula. We discuss several potential explanations for these differences.

The zodiacal cloud is a thick circumsolar disk of small debris particles produced by asteroid collisions and comets. Here, we present a zodiacal cloud model based on the orbital properties and lifetimes of comets and asteroids, and on the dynamical evolution of dust after ejection. The model is quantitatively constrained by IRAS observations of thermal emission, but also qualitatively consistent with other zodiacal cloud observations. We find that 85-95% of the observed mid-infrared emission is produced by particles from the Jupiter-family comets (JFCs) and $<$10% by dust from long period comets. Asteroidal dust is found to be present at $<$10%. We suggest that spontaneous disruptions of JFCs, rather than the usual cometary activity driven by sublimating volatiles, is the main mechanism that librates cometary particles into the zodiacal cloud. Our results imply that JFC particles represent $\sim$85% of the total mass influx at Earth. Since their atmospheric entry speeds are typically low ($\approx$14.5 km s$^{-1}$ mean for D=100-200 $\mu$m with $\approx$12 km s$^{-1}$ being the most common case), many JFC grains should survive frictional heating and land on the Earth's surface. This explains why most micrometeorites collected in antarctic ice have primitive carbonaceous composition. The present mass of the inner zodiacal cloud at $<$5 AU is estimated to be 1-$2\times10^{19}$ g, mainly in D=100-200 $\mu$m particles. The inner zodiacal cloud should have been $>10^4$ times brighter during the Late Heavy Bombardment (LHB) epoch $\approx$3.8 Gyr ago, when the outer planets scattered numerous comets into the inner solar system. The bright debris disks with a large 24-$\mu$m excess observed around mature stars may be an indication of massive cometary populations existing in those systems.

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter D~100km in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ~100km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk ``jumped'' from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100km, probably even 1,000km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size-range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-Neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the `Nice model'. This leads us to predict that the HD of hot population should have the same slope of the HD of the cold population for 6.5 < H < 9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H~9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-Neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. (2004) and Fuentes and Holman (2008). Finally, we show that the constraint on the total mass of the primordial trans-Neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.

The recent discovery of the 2003EL61 collisional family in the Kuiper belt (Brown et al., 2007) is surprising because the formation of such a family is a highly improbable event in today's belt. Assuming Brown et al.'s estimate of the size of progenitors, we find that the probability that a Kuiper belt object was involved in such a collision since primordial times is less than roughly 0.001. In addition, it is not possible for the collision to have occurred in a massive primordial Kuiper belt because the dynamical coherence of the family would not have survived whatever event produced the currently observed orbital excitation. Here we suggest that the family is the result of a collision between two scattered disk objects. We show that the probability that a collision occurred between two such objects with sizes similar to those advocated in Brown et al. (2007) and that the center of mass of the resulting family is on an orbit typical of the Kuiper belt can be as large as 47%. Given the large uncertainties involved in this estimate, this result is consistent with the existence of one such family. If true, this result has implications far beyond the origin of a single collisional family, because it shows that collisions played an important role in shaping the dynamical structure of the small body populations that we see today.