During its initial orbital phase in early 2019, the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) asteroid sample return mission detected small particles apparently emanating from the surface of the near-Earth asteroid (101955) Bennu in optical navigation images. Identification and characterization of the physical and dynamical properties of these objects became a mission priority in terms of both spacecraft safety and scientific investigation. Traditional techniques for particle identification and tracking typically rely on manual inspection and are often time-consuming. The large number of particles associated with the Bennu events and the mission criticality rendered manual inspection techniques infeasible for long-term operational support. In this work, we present techniques for autonomously detecting potential particles in monocular images and providing initial correspondences between observations in sequential images, as implemented for the OSIRIS-REx mission.

Astrometric positions of moving objects in the Solar System have been measured using a variety of star catalogs in the past. Previous work has shown that systematic errors in star catalogs can affect the accuracy of astrometric observations. That, in turn, can influence the resulting orbit fits for minor planets. In order to quantify these systematic errors, we compare the positions and proper motion of stellar sources in the most utilized star catalogs to the second release of the Gaia star catalog. The accuracy of Gaia astrometry allows us to unambiguously identify local biases and derive a scheme that can be used to correct past astrometric observations of solar system objects. Here we provide a substantially improved debiasing scheme for 26 astrometric catalogs that were extensively used in minor planet astrometry. Revised corrections near the galactic center eliminate artifacts that could be traced back to reference catalogs used in previous debiasing schemes. Median differences in stellar positions between catalogs now tend to be on the order of several tens of milliarcseconds (mas) but can be as large as 175 mas. Median stellar proper motion corrections scatter around 0.3 mas/yr and range from 1 to 4 mas/yr for star catalogs with and without proper motion, respectively. The tables in this work are meant to be applied to existing optical observations. They are not intended to correct new astrometric measurments as those should make use of the Gaia astrometric catalog. Since previous debiasing schemes already reduced systematics in past observations to a large extent, corrections beyond the current work may not be needed in the foreseeable future.

We observed the near-Earth asteroid (101955) Bennu from the ground in 1999 and 2005, and with the Hubble Space Telescope in 2012, to constrain its rotation rate. The data reveal an acceleration of $2.64 \pm 1.05 \times 10^{-6} \mathrm{deg\ day}^{-2}$, which could be due to a change in the moment of inertia of Bennu or to spin up from the YORP effect or other source of angular momentum. The best solution is within 1 sigma of the period determined by Nolan et al. (2013). The OSIRIS-REx mission will determine the rotation state independently in 2019. Those measurements should show whether the change in rotation rate is a steady increase (due, for example, to the YORP effect) or some other phenomenon. The precise shape and surface properties measured by the OSIRIS-REx science team will allow for a better understanding of variations in rotation rate of small asteroids.

Making an inventory of the Solar System is one of the four fundamental science requirements for the Large Synoptic Survey Telescope (LSST). The current baseline footprint for LSST's main Wide-Fast-Deep (WFD) Survey observes the sky below 0$^\circ$ declination, which includes only half of the ecliptic plane. Critically, key Solar System populations are asymmetrically distributed on the sky: they will be entirely missed, or only partially mapped, if only the WFD occurs. We propose a Northern Ecliptic Spur (NES) mini survey, observing the northern sky up to +10$^\circ$ ecliptic latitude, to maximize Solar System science with LSST. The mini survey comprises a total area of $\sim$5800 deg$^2$/604 fields, with 255 observations/field over the decade, split between g,r, and z bands. Our proposed survey will 1) obtain a census of main-belt comets; 2) probe Neptune's past migration history, by exploring the resonant structure of the Kuiper belt and the Neptune Trojan population; 3) explore the origin of Inner Oort cloud objects and place significant constraints on the existence of a hypothesized planet beyond Neptune; and 4) enable precise predictions of KBO stellar occultations. These high-ranked science goals of the Solar System Science Collaboration are only achievable with this proposed northern survey.

1I/`Oumuamua is the first confirmed interstellar body in our Solar System. Here we report on observations of `Oumuamua made with the Spitzer Space Telescope on 2017 November 21--22 (UT). We integrated for 30.2~hours at 4.5 micron (IRAC channel 2). We did not detect the object and place an upper limit on the flux of 0.3 uJy (3sigma). This implies an effective spherical diameter less than [98, 140, 440] meters and albedo greater than [0.2, 0.1, 0.01] under the assumption of low, middle, or high thermal beaming parameter eta, respectively. With an aspect ratio for `Oumuamua of 6:1, these results correspond to dimensions of [240:40, 341:57, 1080:180] meters, respectively. We place upper limits on the amount of dust, CO, and CO2 coming from this object that are lower than previous results; we are unable to constrain the production of other gas species. Both our size and outgassing limits are important because `Oumuamua's trajectory shows non-gravitational accelerations that are sensitive to size and mass and presumably caused by gas emission. We suggest that `Oumuamua may have experienced low-level post-perihelion volatile emission that produced a fresh, bright, icy mantle. This model is consistent with the expected eta value and implied high albedo value for this solution, but, given our strict limits on CO and CO2, requires another gas species --- probably H2O --- to explain the observed non-gravitational acceleration. Our results extend the mystery of `Oumuamua's origin and evolution.

We present lightcurves and derive periods and amplitudes for a subset of 38 near earth objects (NEOs) observed at 4.5 microns with the IRAC camera on the the Spitzer Space Telescope, many of them having no previously reported rotation periods. This subset was chosen from about 1800 IRAC NEO observations as having obvious periodicity and significant amplitude. For objects where the period observed did not sample the full rotational period, we derived lower limits to these parameters based on sinusoidal fits. Lightcurve durations ranged from 42 to 544 minutes, with derived periods from 16 to 400 minutes. We discuss the effects of lightcurve variations on the thermal modeling used to derive diameters and albedos from Spitzer photometry. We find that both diameters and albedos derived from the lightcurve maxima and minima agree with our previously published results, even for extreme objects, showing the conservative nature of the thermal model uncertainties. We also evaluate the NEO rotation rates, sizes, and their cohesive strengths.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

We perform high fidelity simulations of a wide-field telescopic survey searching for Near-Earth Objects (NEO) larger than 140\,m, focusing on the observational and detection model, detection efficiency and accuracy. As a test survey we select the Large Synoptic Survey Telescope. We use its proposed pointings for a 10-year mission and model the detection of near-Earth objects in the fields. We discuss individual model parameters for magnitude losses, vignetting, fading, asteroid rotation and colors, fill factor, limiting magnitude, rate of motion, field shape and rotation and survey patterns and we assess results in terms of the cumulative completeness of the detected population as a function of size and time. Additionally, we examine the sources of modeling uncertainty and derive the overall NEO population completeness for the baseline LSST survey to be $55\pm5$% for NEOs with absolute magnitude brighter than 22. Including already discovered objects and ongoing surveys, the population completeness at the end of the LSST baseline survey should reach $\sim 77\%$.

We have conducted a detailed simulation of LSST's ability to link near-Earth and main belt asteroid detections into orbits. The key elements of the study were a high-fidelity detection model and the presence of false detections in the form of both statistical noise and difference image artifacts. We employed the Moving Object Processing System (MOPS) to generate tracklets, tracks and orbits with a realistic detection density for one month of the LSST survey. The main goals of the study were to understand whether a) the linking of Near-Earth Objects (NEOs) into orbits can succeed in a realistic survey, b) the number of false tracks and orbits will be manageable, and c) the accuracy of linked orbits would be sufficient for automated processing of discoveries and attributions. We found that the overall density of asteroids was more than 5000 per LSST field near opposition on the ecliptic, plus up to 3000 false detections per field in good seeing. We achieved 93.6% NEO linking efficiency for H<22 on tracks composed of tracklets from at least three distinct nights within a 12-day interval. The derived NEO catalog was comprised of 96% correct linkages. Less than 0.1% of orbits included false detections, and the remainder of false linkages stemmed from main belt confusion, which was an artifact of the short time span of the simulation. The MOPS linking efficiency can be improved by refined attribution of detections to known objects and by improved tuning of the internal kd-tree linking algorithms.

We performed a statistical analysis of the astrometric errors for the major asteroid surveys. We analyzed the astrometric residuals as a function of observation epoch, observed brightness and rate of motion, finding that astrometric errors are larger for faint observations and some stations improved their astrometric quality over time. Based on this statistical analysis we develop a new weighting scheme to be used when performing asteroid orbit determination. The proposed weights result in ephemeris predictions that can be conservative by a factor as large as 1.5. However, the new scheme is more robust with respect to outliers and better handles faint detections.

This report describes the methodology and results of an assessment study of the performance of the Large Synoptic Survey Telescope (LSST) in its planned efforts to detect and catalog near-Earth objects (NEOs). The baseline LSST survey approach is designed to make only two visits to a given field in a given night, leading to two possible NEO detections per night. These nightly pairs must be linked across nights to derive orbits of moving objects. However, the presence of false detections in the data stream leads to the possibility of high rates of false tracklets, and the ensuing risk that the resulting orbit catalog may be contaminated by false orbits. To understand the challenges stemming from a two-visit-per-night cadence, we conducted high-fidelity linkage tests on a full-density simulated LSST detection stream. We also sought to quantify the overall performance of LSST as an NEO discovery system under the hypothesis that detections from the baseline cadence can be successfully linked. Our simulations revealed that in 10 years LSST would catalog ~60% of NEOs with absolute magnitude H<22, which is a proxy for 140m and larger objects. This results neglects linking losses and the contribution of any other NEO surveys. Including our worst-case linking efficiency we reach a overall performance assessment of 55% completeness of NEOs with H<22. We estimate that survey mis-modeling could account for systematic errors of up to 5%. We find that restricting the evaluation metric to so-called Potentially Hazardous Asteroids (PHAs) increases the completeness by 3-4%, and that including the benefits of past and expected future NEO survey activity increases completeness at the end of the baseline LSST survey by 15-20%. Assembling these results leads to a projection that by the end of the baseline LSST survey the NEO catalog will be 80+/-5% complete for PHAs with H < 22.

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.

Near Earth Objects (NEOs) are small Solar System bodies whose orbits bring them close to the Earth's orbit. We are carrying out a Warm Spitzer Cycle 11 Exploration Science program entitled NEOSurvey --- a fast and efficient flux-limited survey of 597 known NEOs in which we derive diameter and albedo for each target. The vast majority of our targets are too faint to be observed by NEOWISE, though a small sample has been or will be observed by both observatories, which allows for a cross-check of our mutual results. Our primary goal is to create a large and uniform catalog of NEO properties. We present here the first results from this new program: fluxes and derived diameters and albedos for 80 NEOs, together with a description of the overall program and approach, including several updates to our thermal model. The largest source of error in our diameter and albedo solutions, which derive from our single band thermal emission measurements, is uncertainty in eta, the beaming parameter used in our thermal modeling; for albedos, improvements in Solar System absolute magnitudes would also help significantly. All data and derived diameters and albedos from this entire program are being posted on a publicly accessible webpage at nearearthobjects.nau.edu .

Near-Earth asteroid 2014 AA entered the Earth's atmosphere on 2014 January 2, only 21 hours after being discovered by the Catalina Sky Survey. In this paper we compute the trajectory of 2014 AA by combining the available optical astrometry, seven ground-based observations over 69 minutes, and the International Monitoring system detection of the atmospheric impact infrasonic airwaves in a least-squares orbit estimation filter. The combination of these two sources of observations results in a tremendous improvement in the orbit uncertainties. The impact time is 3:05 UT with a 1-sigma uncertainty of 6 min, while the impact location corresponds to a west longitude of 44.7 deg and a latitude of 13.1 deg with a 1-sigma uncertainty of 140 km. The minimum impact energy estimated from the infrasound data and the impact velocity result in an estimated minimum mass of 22.6 t. By propagating the trajectory of 2014 AA backwards we find that the only window for finding precovery observations is for the three days before its discovery.

The Mars flyby of C/2013 A1 (Siding Spring) represented a unique opportunity for imaging a long-period comet and resolving its nucleus and rotation period. Because of the small encounter distance and the high relative velocity, the goal of successfully observing C/2013 A1 from the Mars orbiting spacecrafts posed strict accuracy requirements on the comet's ephemerides. These requirements were hard to meet, as comets are known for being highly unpredictable: astrometric observations can be significantly biased and nongravitational perturbations affect comet trajectories. Therefore, even prior to the encounter, we remeasured a couple of hundred astrometric images obtained with ground-based and Earth-orbiting telescopes. We also observed the comet with the Mars Reconnaissance Orbiter's High Resolution Imaging Science Experiment (HiRISE) camera on 2014 October 7. In particular, these HiRISE observations were decisive in securing the trajectory and revealed that out-of-plane nongravitational perturbations were larger than previously assumed. Though the resulting ephemeris predictions for the Mars encounter allowed observations of the comet from the Mars orbiting spacecrafts, post-encounter observations show a discrepancy with the pre-encounter trajectory. We reconcile this discrepancy by employing the Rotating Jet Model, which is a higher fidelity model for nongravitational perturbations and provides an estimate of C/2013 A1's spin pole.

We describe systematic ranging, an orbit determination technique especially suitable to assess the near-term Earth impact hazard posed by newly discovered asteroids. For these late warning cases, the time interval covered by the observations is generally short, perhaps a few hours or even less, which leads to severe degeneracies in the orbit estimation process. The systematic ranging approach gets around these degeneracies by performing a raster scan in the poorly-constrained space of topocentric range and range rate, while the plane of sky position and motion are directly tied to the recorded observations. This scan allows us to identify regions corresponding to collision solutions, as well as potential impact times and locations. From the probability distribution of the observation errors, we obtain a probability distribution in the orbital space and then estimate the probability of an Earth impact. We show how this technique is effective for a number of examples, including 2008 TC3 and 2014 AA, the only two asteroids to date discovered prior to impact.

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.

Asteroid 2009 FD could impact Earth between 2185 and 2196. The long term propagation to the possible impacts and the intervening planetary encounters make 2009 FD one of the most challenging asteroids in terms of hazard assessment. To compute accurate impact probabilities we model the Yarkovsky effect by using the available physical characterization of 2009 FD and general properties of the Near Earth Asteroid population. We perform the hazard assessment with two independent methods: the first method is a generalization of the standard impact monitoring algorithms in use by NEODyS and Sentry, while the second one is based on a Monte Carlo approach. Both methods generate orbital samples in a 7 dimensional space that includes orbital elements and the parameter characterizing the Yarkovsky effect. The highest impact probability is $2.7 \times 10^{-3}$ for an impact during the 2185 Earth encounter. Impacts after 2185 corresponding to resonant returns are possible, the most relevant being in 2190 with a probability of $3 \times 10^{-4}$. Both numerical methods can be used in the future to handle similar cases. The structure of resonant returns and the list of the possible keyholes on the Target Plane of the scattering encounter in 2185 can be predicted by an analytic theory.

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.

We provide a scheme to correct asteroid astrometric observations for star catalog systematic errors due to inaccurate star positions and proper motions. As reference we select the most accurate stars in the PPMXL catalog, i.e., those based on 2MASS astrometry. We compute position and proper motion corrections for 19 of the most used star catalogs. The use of these corrections provides better ephemeris predictions and improves the error statistics of astrometric observations, e.g., by removing most of the regional systematic errors previously seen in Pan-STARRS PS1 asteroid astrometry. The correction table is publicly available at ftp://ssd.jpl.nasa.gov/pub/ssd/debias/debias_2014.tgz and can be freely used in orbit determination algorithms to obtain more reliable asteroid trajectories.

We report on observations of near-Earth asteroid 2011 MD with the Spitzer Space Telescope. We have spent 19.9 h of observing time with channel 2 (4.5 {\mu}m) of the Infrared Array Camera and detected the target within the 2{\sigma} positional uncertainty ellipse. Using an asteroid thermophysical model and a model of nongravitational forces acting upon the object we constrain the physical properties of 2011 MD, based on the measured flux density and available astrometry data. We estimate 2011 MD to be 6 (+4/-2) m in diameter with a geometric albedo of 0.3 (+0.4/-0.2) (uncertainties are 1{\sigma}). We find the asteroid's most probable bulk density to be 1.1 (+0.7/-0.5) g cm^{-3}, which implies a total mass of (50-350) t and a macroporosity of >=65%, assuming a material bulk density typical of non-primitive meteorite materials. A high degree of macroporosity suggests 2011 MD to be a rubble-pile asteroid, the rotation of which is more likely to be retrograde than prograde.

Comet C/2013 A1 (siding Spring) will experience a high velocity encounter with Mars on October 19, 2014 at a distance of 135,000 km +- 5000 km from the planet center. We present a comprehensive analysis of the trajectory of both the comet nucleus and the dust tail. The nucleus of C/2013 A1 cannot impact on Mars even in the case of unexpectedly large nongravitational perturbations. Furthermore, we compute the required ejection velocities for the dust grains of the tail to reach Mars as a function of particle radius and density and heliocentric distance of the ejection. A comparison between our results and the most current modeling of the ejection velocities suggests that impacts are possible only for millimeter to centimeter size particles released more than 13 au from the Sun. However, this level of cometary activity that far from the Sun is considered extremely unlikely. The arrival time of these particles spans a 20-minute time interval centered at October 19, 2014 at 20:09 TDB, i.e., around the time that Mars crosses the orbital plane of C/2013 A1. Ejection velocities larger than currently estimated by a factor >2 would allow impacts for smaller particles ejected as close as 3 au from the Sun. These particles would reach Mars from 43 to 130 min after the nominal close approach epoch of the purely gravitational trajectory of the nucleus.

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m.

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.

In this paper we perform an assessment of the 2880 Earth impact risk for asteroid (29075) 1950 DA. To obtain reliable predictions we analyze the contribution of the observational dataset and the astrometric treatment, the numerical error in the long-term integration, and the different accelerations acting on the asteroid. The main source of uncertainty is the Yarkovsky effect, which we statistically model starting from 1950 DA's available physical characterization, astrometry, and dynamical properties. Before the release of 2012 radar data, this modeling suggests that 1950 DA has 99% likelihood of being a retrograde rotator. By using a 7-dimensional Monte Carlo sampling we map 1950 DA's uncertainty region to the 2880 close approach b-plane and find a 5 x 10^-4 impact probability. With the recently released 2012 radar observations, the direct rotation is definitely ruled out and the impact probability decreases to 2.5 x 10^-4.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

We assess the risk of an Earth impact for asteroid (99942) Apophis by means of a statistical analysis accounting for the uncertainty of both the orbital solution and the Yarkovsky effect. We select those observations with either rigorous uncertainty information provided by the observer or a high established accuracy. For the Yarkovsky effect we perform a Monte Carlo simulation that fully accounts for the uncertainty in the physical characterization, especially for the unknown spin orientation. By mapping the uncertainty information onto the 2029 b-plane and identifying the keyholes corresponding to subsequent impacts we assess the impact risk for future encounters. In particular, we find an impact probability greater than 10^-6 for an impact in 2068. We analyze the stability of the impact probability with respect to the assumptions on Apophis' physical characterization and consider the possible effect of the early 2013 radar apparition.

We seek evidence of the Yarkovsky effect among Near Earth Asteroids (NEAs) by measuring the Yarkovsky-related orbital drift from the orbital fit. To prevent the occurrence of unreliable detections we employ a high precision dynamical model, including the Newtonian attraction of 16 massive asteroids and the planetary relativistic terms, and a suitable astrometric data treatment. We find 21 NEAs whose orbital fits show a measurable orbital drift with a signal to noise ratio (SNR) greater than 3. The best determination is for asteroid (101955) 1999 RQ36, resulting in the recovery of one radar apparition and an orbit improvement by two orders of magnitude. In addition, we find 16 cases with a lower SNR that, despite being less reliable, are good candidates for becoming stronger detections in the future. In some cases it is possible to constrain physical quantities otherwise unknown by means of the detected orbital drift. Furthermore, the distribution of the detected orbital drifts shows an excess of retrograde rotators that can be connected to the delivery mechanism from the most important NEA feeding resonances and allows us to infer the distribution for NEAs obliquity. We discuss the implications of the Yarkovsky effect for impact predictions. In particular, for asteroid (29075) 1950 DA our results favor a retrograde rotation that would rule out an impact in 2880.

We have identified and quantified semi-major axis drifts in Near-Earth Asteroids (NEAs) by performing orbital fits to optical and radar astrometry of all numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the most reliable and strongest drift rates. Our selection criteria include a Yarkovsky sensitivity metric that quantifies the detectability of semi-major axis drift in any given data set, a signal-to-noise metric, and orbital coverage requirements. In 42 cases, the observed drifts (~10^-3 AU/Myr) agree well with numerical estimates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is the dominant non-gravitational process affecting these orbits, and allows us to derive constraints on asteroid physical properties. In 12 cases, the drifts exceed nominal Yarkovsky predictions, which could be due to inaccuracies in our knowledge of physical properties, faulty astrometry, or modeling errors. If these high rates cannot be ruled out by further observations or improvements in modeling, they would be indicative of the presence of an additional non-gravitational force, such as that resulting from a loss of mass of order a kilogram per second. We define the Yarkovsky efficiency f_Y as the ratio of the change in orbital energy to incident solar radiation energy, and we find that typical Yarkovsky efficiencies are ~10^-5.

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.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We have performed a simulation of a next generation sky survey's (Pan-STARRS 1) efficiency for detecting Earth-impacting asteroids. The steady-state sky-plane distribution of the impactors long before impact is concentrated towards small solar elongations (Chesley and Spahr, 2004) but we find that there is interesting and potentially exploitable behavior in the sky-plane distribution in the months leading up to impact. The next generation surveys will find most of the dangerous impactors (>140m diameter) during their decade-long survey missions though there is the potential to miss difficult objects with long synodic periods appearing in the direction of the Sun, as well as objects with long orbital periods that spend much of their time far from the Sun and Earth. A space-based platform that can observe close to the Sun may be needed to identify many of the potential impactors that spend much of their time interior to the Earth's orbit. The next generation surveys have a good chance of imaging a bolide like 2008TC3 before it enters the atmosphere but the difficulty will lie in obtaining enough images in advance of impact to allow an accurate pre-impact orbit to be computed.

The potentially hazardous asteroid (101955) 1999 RQ36 has the possibility of collision with the Earth in the latter half of the 22nd century, well beyond the traditional 100-year time horizon for routine impact monitoring. The probabilities accumulate to a total impact probability of approximately 10E-3, with a pair of closely related routes to impact in 2182 comprising more than half of the total. The analysis of impact possibilities so far in the future is strongly dependent on the action of the Yarkovsky effect, which raises new challenges in the careful assessment of longer term impact hazards. Even for asteroids with very precisely determined orbits, a future close approach to Earth can scatter the possible trajectories to the point that the problem becomes like that of a newly discovered asteroid with a weakly determined orbit. If the scattering takes place late enough so that the target plane uncertainty is dominated by Yarkovsky accelerations then the thermal properties of the asteroid,which are typically unknown, play a major role in the impact assessment. In contrast, if the strong planetary interaction takes place sooner, while the Yarkovsky dispersion is still relatively small compared to that derived from the measurements, then precise modeling of the nongravitational acceleration may be unnecessary.

Near Earth Objects (NEOs) are fragments of remnant primitive bodies that date from the era of Solar System formation. At present, the physical properties and origins of NEOs are poorly understood. We have measured thermal emission from three NEOs -- (6037) 1988 EG, 1993 GD, and 2005 GL -- with Spitzer's IRAC instrument at 3.6, 4.5, 5.8, and 8.0 microns (the last object was detected only at 5.8 and 8.0 microns). The diameters of these three objects are 400 m, 180 m, and 160 m, respectively, with uncertainties of around 20% (including both observational and systematic errors). For all three the geometric albedos are around 0.30, in agreement with previous results that most NEOs are S-class asteroids. For the two objects detected at 3.6 and 4.5 microns, diameters and albedos based only on those data agree with the values based on modeling the data in all four bands. This agreement, and the high sensitivity of IRAC, show the promise of the Spitzer Warm Mission for determining the physical parameters for a large number of NEOs.