Search

This work is dedicated to debias the Near-Earth Objects (NEO) population based on observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS) telescopes. We have applied similar methods used to develop the recently released NEO model generator (NEOMOD), once debiasing the NEO population using data from Catalina Sky Survey (CSS) G96 telescope. ATLAS is composed of four different telescopes. We first analyzed observational data from each of all four telescopes separately and later combined them. Our results highlight main differences between CSS and ATLAS, e.g., sky coverage and survey power at debiasing the NEO population. ATLAS has a much larger sky coverage than CSS, allowing it to find bright NEOs that would be constantly "hiding" from CSS. Consequently, ATLAS is more powerful than CSS at debiasing the NEO population for H $\lesssim$ 19. With its intrinsically greater sensitivity and emphasis on observing near opposition, CSS excels in the debiasing of smaller objects. ATLAS, as an all sky survey designed to find imminent hazardous objects, necessarily spends a significant fraction of time looking at places on the sky where objects do not appear, reducing its power for debiasing the population of small objects. We estimate a NEO population completeness of $\approx$ 88%$^{+3\%}_{-2\%}$ for H $<$ 17.75 and $\approx$ 36%$^{+1\%}_{-1\%}$ for H $<$ 22.25. Those numbers are similar to previous estimates (within error bars for H $<$ 17.75) from CSS, yet, around 3% and 8% smaller at their face values, respectively. We also confirm previous finding that the $\nu_6$ secular resonance is the main source of small and faint NEOs at H = 28, whereas the 3:1 mean motion resonance with Jupiter dominates for larger and brighter NEOs at H = 15., Comment: Accepted for publication in Icarus, 26 pages, 13 figures, 2 tables

Our previous model (NEOMOD2) for the orbital and absolute magnitude distribution of Near Earth Objects (NEOs) was calibrated on the Catalina Sky Survey observations between 2013 and 2022. Here we extend NEOMOD2 to include visible albedo information from the Wide-Field Infrared Survey Explorer. The debiased albedo distribution of NEOs can be approximated by the sum of two Rayleigh distributions with the scale parameters p_V,dark=0.03 and p_V,bright=0.17. We find evidence for smaller NEOs having (on average) higher albedos than larger NEOs; this is likely a consequence of the size-dependent sampling of different main belt sources. These inferences and the absolute magnitude distribution from NEOMOD2 are used to construct the debiased size distribution of NEOs. We estimate 830+/-60 NEOs with diameters D>1 km and 20,000+/-2,000 NEOs with D>140 m. The new model, NEOMOD3, is available via the NEOMOD Simulator -- an easy-to-operate code that can be used to generate user-defined samples (orbits, sizes and albedos) from the model., Comment: Icarus, in press

The NASA Double Asteroid Redirection Test (DART) spacecraft successfully impacted the Didymos-Dimorphos binary asteroid system on 2022 September 26 UTC. We provide an update to its pre-impact mutual orbit and estimate the post-impact physical and orbital parameters, derived using ground-based photometric observations taken from July 2022 to February 2023. We found that the total change of the orbital period was $-33.240 \pm 0.072$ min. (all uncertainties are 3$\sigma$). We obtained the eccentricity of the post-impact orbit to be $0.028 \pm 0.016$ and the apsidal precession rate of $7.3 \pm 2.0$ deg./day from the impact to 2022 December 2. The data taken later in December to February suggest that the eccentricity dropped close to zero or the orbit became chaotic approximately 70 days after the impact. Most of the period change took place immediately after the impact but in a few weeks following the impact it was followed by additional change of $-27^{+19}_{-58}$ seconds or $-19 \pm 18$ seconds (the two values depend on the approach we used to describe the evolution of the orbital period after the impact -- an exponentially decreasing angular acceleration or an assumption of a constant orbital period, which changed abruptly some time after the impact, respectively). We estimate the pre-impact Dimorphos-Didymos size ratio was $0.223 \pm 0.012$ and the post-impact is $0.202 \pm 0.018$, which indicates a marginally significant reduction of Dimorphos' volume by ($9 \pm 9) \%$ as the result of the impact., Comment: 31 pages, 7 figures

Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This resolves previous model inconsistencies related to the population of large NEOs. We estimate there are 936+/-29 NEOs with absolute magnitude H<17.75 (diameter D>1 km for the reference albedo p_V=0.14). The slope of the NEO size distribution for H=25-28 is found to be relatively shallow (cumulative index 2.6) and the number of H<28 NEOs (D>9 m) is determined to be (1.20+/-0.04)x10^7. Small NEOs have a different orbital distribution and higher impact probabilities than large NEOs. We estimate 0.034+/-0.002 impacts of H<28 NEOs on the Earth per year, which is near the low end of the impact flux range inferred from atmospheric bolide observations. Relative to a model where all NEOs are delivered directly from the main belt, the population of small NEOs detected by G96 shows an excess of low-eccentricity orbits with a=1--1.6 au that appears to increase with H. We suggest that the population of very small NEOs is boosted by tidal disruption of large NEOs during close encounters to the terrestrial planets. When the effect of tidal disruption is (approximately) accounted for in the model, we estimate 0.06+/-0.01 impacts of H<28 NEOs on the Earth per year, which is more in line with the bolide data., Comment: Icarus, in press. arXiv admin note: text overlap with arXiv:2306.09521

On 26 September 2022, the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the satellite of binary near-Earth asteroid (65803) Didymos. This demonstrated the efficacy of a kinetic impactor for planetary defense by changing the orbital period of Dimorphos by 33 minutes (Thomas et al. 2023). Measuring the period change relied heavily on a coordinated campaign of lightcurve photometry designed to detect mutual events (occultations and eclipses) as a direct probe of the satellite's orbital period. A total of 28 telescopes contributed 224 individual lightcurves during the impact apparition from July 2022 to February 2023. We focus here on decomposable lightcurves, i.e. those from which mutual events could be extracted. We describe our process of lightcurve decomposition and use that to release the full data set for future analysis. We leverage these data to place constraints on the post-impact evolution of ejecta. The measured depths of mutual events relative to models showed that the ejecta became optically thin within the first ~1 day after impact, and then faded with a decay time of about 25 days. The bulk magnitude of the system showed that ejecta no longer contributed measurable brightness enhancement after about 20 days post-impact. This bulk photometric behavior was not well represented by an HG photometric model. An HG1G2 model did fit the data well across a wide range of phase angles. Lastly, we note the presence of an ejecta tail through at least March 2023. Its persistence implied ongoing escape of ejecta from the system many months after DART impact., Comment: 52 pages, 5 tables, 9 figures, accepted to PSJ

With the successful impact of the NASA DART spacecraft in the Didymos-Dimorphos binary asteroid system, we provide an initial analysis of the post-impact perturbed binary asteroid dynamics. To compare our simulation results with observations, we introduce a set of "observable elements" calculated using only the physical separation of the binary asteroid, rather than traditional Keplerian elements. Using numerical methods that treat the fully spin-orbit-coupled dynamics, we estimate the system's mass and the impact-induced changes in orbital velocity, semimajor axis, and eccentricity. We find that the changes to the mutual orbit depend strongly on the separation distance between Didymos and Dimorphos at the time of impact. If Dimorphos enters a tumbling state after the impact, this may be observable through changes in the system's eccentricity and orbit period. We also find that any DART-induced reshaping of Dimorphos would generally reduce the required change in orbital velocity to achieve the measured post-impact orbit period and will be assessed by the ESA Hera mission in 2027., Comment: Accepted for publication in PSJ

Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $\nu_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 \times 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $\theta \simeq 0^\circ$ for $a<2.5$~au and $\theta \simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model., Comment: AJ

The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of asteroid deflection. Here we report the first determination of the momentum transferred to an asteroid by kinetic impact. Based on the change in the binary orbit period, we find an instantaneous reduction in Dimorphos's along-track orbital velocity component of 2.70 +/- 0.10 mm/s, indicating enhanced momentum transfer due to recoil from ejecta streams produced by the impact. For a Dimorphos bulk density range of 1,500 to 3,300 kg/m$^3$, we find that the expected value of the momentum enhancement factor, $\beta$, ranges between 2.2 and 4.9, depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have equal densities of 2,400 kg/m$^3$, $\beta$= 3.61 +0.19/-0.25 (1 $\sigma$). These $\beta$ values indicate that significantly more momentum was transferred to Dimorphos from the escaping impact ejecta than was incident with DART. Therefore, the DART kinetic impact was highly effective in deflecting the asteroid Dimorphos., Comment: accepted by Nature

The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multi-band wide-field synoptic survey will transform our view of the solar system, with the discovery and monitoring of over 5 million small bodies.The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor planets and passing interstellar objects. Creating an inventory of the solar system is one of the four main LSST science drivers. The LSST observing cadence is a complex optimization problem that must balance the priorities and needs of all the key LSST science areas. To design the best LSST survey strategy, a series of operation simulations using the Rubin Observatory scheduler have been generated to explore the various options for tuning observing parameters and prioritizations. We explore the impact of the various simulated LSST observing strategies on studying the solar system's small body reservoirs. We examine what are the best observing scenarios and review what are the important considerations for maximizing LSST solar system science. In general, most of the LSST cadence simulations produce +/-5% or less variations in our chosen key metrics, but a subset of the simulations significantly hinder science returns with much larger losses in the discovery and light curve metrics., Comment: Accepted to ApJS, 103 pages (including references), 43 figures, 9 Tables. Videos will be available in the online journal formatted and published version of the paper [v2.0 submission corrects the author list metadata from the arxiv initial submission and updates the abstract]

While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority space mission related to asteroid mitigation. NASA's Double Asteroid Redirection Test (DART) mission is the first full-scale test of kinetic impact technology. The mission's target asteroid was Dimorphos, the secondary member of the S-type binary near-Earth asteroid (65803) Didymos. This binary asteroid system was chosen to enable ground-based telescopes to quantify the asteroid deflection caused by DART's impact. While past missions have utilized impactors to investigate the properties of small bodies those earlier missions were not intended to deflect their targets and did not achieve measurable deflections. Here we report the DART spacecraft's autonomous kinetic impact into Dimorphos and reconstruct the impact event, including the timeline leading to impact, the location and nature of the DART impact site, and the size and shape of Dimorphos. The successful impact of the DART spacecraft with Dimorphos and the resulting change in Dimorphos's orbit demonstrates that kinetic impactor technology is a viable technique to potentially defend Earth if necessary., Comment: Accepted by Nature

The Double Asteroid Redirection Test (DART) spacecraft successfully performed the first test of a kinetic impactor for asteroid deflection by impacting Dimorphos, the secondary of near-Earth binary asteroid (65803) Didymos, and changing the orbital period of Dimorphos. A change in orbital period of approximately 7 minutes was expected if the incident momentum from the DART spacecraft was directly transferred to the asteroid target in a perfectly inelastic collision, but studies of the probable impact conditions and asteroid properties indicated that a considerable momentum enhancement ($\beta$) was possible. In the years prior to impact, we used lightcurve observations to accurately determine the pre-impact orbit parameters of Dimorphos with respect to Didymos. Here we report the change in the orbital period of Dimorphos as a result of the DART kinetic impact to be -33.0 +/- 1.0 (3$\sigma$) minutes. Using new Earth-based lightcurve and radar observations, two independent approaches determined identical values for the change in the orbital period. This large orbit period change suggests that ejecta contributed a significant amount of momentum to the asteroid beyond what the DART spacecraft carried., Comment: Accepted by Nature

We report a statistically significant detection of nongravitational acceleration on the sub-kilometer near-Earth asteroid (523599) 2003 RM. Due to its orbit, 2003 RM experiences favorable observing apparitions every 5 years. Thus, since its discovery, 2003 RM has been extensively tracked with ground-based optical facilities in 2003, 2008, 2013, and 2018. We find that the observed plane-of-sky positions cannot be explained with a purely gravity-driven trajectory. Including a transverse nongravitational acceleration allows us to match all observational data, but its magnitude is inconsistent with perturbations typical of asteroids such as the Yarkovsky effect or solar radiation pressure. After ruling out that the orbital deviations are due to a close approach or collision with another asteroid, we hypothesize that this anomalous acceleration is caused by unseen cometary outgassing. A detailed search for evidence of cometary activity with archival and deep observations from Pan-STARRS and the VLT does not reveal any detectable dust production. However, the best-fitting H$_2$O sublimation model allows for brightening due to activity consistent with the scatter of the data. We estimate the production rate required for H$_2$O outgassing to power the acceleration, and find that, assuming a diameter of 300 m, 2003 RM would require Q(H$_2$O)$\sim10^{23}$ molec s$^{-1}$ at perihelion. We investigate the recent dynamical history of 2003 RM and find that the object most likely originated in the mid-to-outer main belt ($\sim86\%$) as opposed to from the Jupiter-family comet region ($\sim11\%$). Further observations, especially in the infrared, could shed light on the nature of this anomalous acceleration., Comment: 14 pages, 9 figures, 2 tables, accepted by the Planetary Science Journal

We report statistically significant detections of non-radial nongravitational accelerations based on astrometric data in the photometrically inactive objects 1998 KY$_{26}$, 2005 VL$_1$, 2016 NJ$_{33}$, 2010 VL$_{65}$, 2016 RH$_{120}$, and 2010 RF$_{12}$. The magnitudes of the nongravitational accelerations are greater than those typically induced by the Yarkovsky effect and there is no radiation-based, non-radial effect that can be so large. Therefore, we hypothesize that the accelerations are driven by outgassing, and calculate implied H$_2$O production rates for each object. We attempt to reconcile outgassing induced acceleration with the lack of visible comae or photometric activity via the absence of surface dust and low levels of gas production. Although these objects are small and some are rapidly rotating, surface cohesive forces are stronger than the rotational forces and rapid rotation alone cannot explain the lack of surface debris. It is possible that surface dust was removed previously, perhaps via outgassing activity that increased the rotation rates to their present day value. We calculate dust production rates of order $\sim10^{-4}$ g s$^{-1}$ in each object assuming that the nuclei are bare, within the upper limits of dust production from a sample stacked image of 1998 KY$_{26}$ of $\dot{M}_{\rm Dust}<0.2$ g s$^{-1}$. This production corresponds to brightness variations of order $\sim0.0025\%$, which are undetectable in extant photometric data. We assess the future observability of each of these targets, and find that the orbit of 1998 KY$_{26}$ -- which is also the target for the extended Hayabusa2 mission -- exhibits favorable viewing geometry before 2025., Comment: 21 pages, 9 figures, 3 tables, accepted for publication in the Planetary Science Journal

We used the times of occultations and eclipses between the components of the 65803 Didymos binary system observed in its lightcurves from 2003-2021 to estimate the orbital parameters of Dimorphos relative to Didymos. We employed a weighted least-squares approach and a modified Keplerian orbit model in order to accommodate the effects from non-gravitational forces such as Binary YORP that could cause a linear change in mean motion over time. We estimate that the period of the mutual orbit at the epoch 2022 September 26.0 TDB, the day of the DART impact, is $11.9214869 \pm 0.000028$~h ($1\sigma$) and that the mean motion of the orbit is changing at a rate of $(5.0\pm 1.0)\times 10^{-18}$~rad s$^{-2}$ $(1\sigma$). The formal $3\sigma$ uncertainty in orbital phase of Dimorphos during the planned Double Asteroid Redirection Test (DART) mission is $5.4^\circ$. Observations from July to September 2022, a few months to days prior to the DART impact, should provide modest improvements to the orbital phase uncertainty and reduce it to about $4.2^\circ$. These results, generated using a relatively simple model, are consistent with those generated using the more sophisticated model of \citet{scheirich22}, which demonstrates the reliability of our method and adds confidence to these mission-critical results., Comment: Accepted to PSJ

NASA's Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos's response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor $\beta$, showing that a particular direction-specific $\beta$ will be directly determined by the DART results, and that a related direction-specific $\beta$ is a figure of merit for a kinetic impact mission. The DART $\beta$ determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos's near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction., Comment: 30 pages, 7 figures. Planetary Science Journal, in press, accepted 2022 September 22

Orbit-determination programs find the orbit solution that best fits a set of observations by minimizing the RMS of the residuals of the fit. For near-Earth asteroids, the uncertainty of the orbit solution may be compatible with trajectories that impact Earth. This paper shows how incorporating the impact condition as an observation in the orbit-determination process results in a robust technique for finding the regions in parameter space leading to impacts. The impact pseudo-observation residuals are the b-plane coordinates at the time of close approach and the uncertainty is set to a fraction of the Earth radius. The extended orbit-determination filter converges naturally to an impacting solution if allowed by the observations. The uncertainty of the resulting orbit provides an excellent geometric representation of the virtual impactor. As a result, the impact probability can be efficiently estimated by exploring this region in parameter space using importance sampling. The proposed technique can systematically handle a large number of estimated parameters, account for nongravitational forces, deal with nonlinearities, and correct for non-Gaussian initial uncertainty distributions. The algorithm has been implemented into a new impact monitoring system at JPL called Sentry-II, which is undergoing extensive testing. The main advantages of Sentry-II over JPL's currently operating impact monitoring system Sentry are that Sentry-II can systematically process orbits perturbed by nongravitational forces and that it is generally more robust when dealing with pathological cases. The runtimes and completeness of both systems are comparable, with the impact probability of Sentry-II for 99% completeness being $3\times10^{-7}$., Comment: 19 pages, 13 figures. To be published in the Astronomical Journal

Vera C. Rubin Observatory will be a key facility for small body science in planetary astronomy over the next decade. It will carry out the Legacy Survey of Space and Time (LSST), observing the sky repeatedly in u, g, r, i, z, and y over the course of ten years using a 6.5 m effective diameter telescope with a 9.6 square degree field of view, reaching approximately r = 24.5 mag (5-{\sigma} depth) per visit. The resulting dataset will provide extraordinary opportunities for both discovery and characterization of large numbers (10--100 times more than currently known) of small solar system bodies, furthering studies of planetary formation and evolution. This white paper summarizes some of the expected science from the ten years of LSST, and emphasizes that the planetary astronomy community should remain invested in the path of Rubin Observatory once the LSST is complete., Comment: White paper submitted to the 2020 Planetary Astronomy Decadal Survey (7 pages, 1 figure)

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., Comment: 23 pages, 10 figures

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., Comment: submitted to Icarus

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., Comment: White Paper submitted in response to the Call for LSST Cadence Optimization White Papers

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., Comment: 16 pages, 4 figures, 3 tables, to appear in the Astrophysical Journal Supplement Series

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., Comment: 7 pages; Feedback welcome

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., Comment: 75 pages. See page 3 for Table of Contents and page 4 for Executive Summary

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., Comment: 8 pages, 7 figures, 2 tables, accepted for publication in MNRAS

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\%$., Comment: 44 pages, 20 figures, 11 tables. arXiv admin note: substantial text overlap with arXiv:1705.06209

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., Comment: 41 pages, 23 figures, 6 tables. arXiv admin note: substantial text overlap with arXiv:1705.06209

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.

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., Comment: 25 pages, 11 figures, 5 tables

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.

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.

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., Comment: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook