Interstellar comets offer direct samples of volatiles from distant protoplanetary disks. 2I/Borisov is the first notably active interstellar comet discovered in our solar system[1]. Comets are condensed samples of the gas, ice, and dust that were in a star's protoplanetary disk during the formation of its planets and inform our understanding on how chemical compositions and abundances vary with distance from the central star. Their orbital migration moves volatiles[2], organic material, and prebiotic chemicals in their host system[3]. In our solar system, hundreds of comets have been observed remotely, and a few have been studied up close by space missions[4]. However, knowledge of extrasolar comets has been limited to what could be gleaned from distant, unresolved observations of cometary regions around other stars, with only one detection of carbon monoxide[5]. Here we report that the coma of 2I/Borisov contains significantly more CO than H2O gas, with abundances of at least 173%, more than three times higher than previously measured for any comet in the inner (<2.5 au) solar system[4]. Our ultraviolet observations of 2I/Borisov provide the first glimpse into the ice content and chemical composition of the protoplanetary disk of another star that is substantially different from our own.

The discovery of the second interstellar object 2I/Borisov on 2019 August 30 raises the question of whether it was ejected recently from a nearby stellar system. Here we compute the asymptotic incoming trajectory of 2I/Borisov, based on both recent and pre-discovery data extending back to December 2018, using a range of force models that account for cometary outgassing. From Gaia DR2 astrometry and radial velocities, we trace back in time the Galactic orbits of 7.4 million stars to look for close encounters with 2I/Borisov. The closest encounter we find took place 910kyr ago with the M0V star Ross 573, at a separation of 0.068pc (90% confidence interval of 0.053-0.09pc) with a relative velocity of 23km/s. This encounter is nine times closer than the closest past encounter identified for the first interstellar object 1I/'Oumuamua. Ejection of 2I/Borisov via a three-body encounter in a binary or planetary system is possible, although such a large ejection velocity is unlikely to be obtained and Ross 573 shows no signs of binarity. We also identify and discuss some other recent close encounters, recognizing that if 2I/Borisov is more than about 10Myr old, our search would be unlikely to find its parent system.

Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two.

Near-Earth Objects (NEOs) that orbit the Sun on or within Earth's orbit are tricky to detect for Earth-based observers due to their proximity to the Sun in the sky. These small bodies hold clues to the dynamical history of the inner solar system as well as the physical evolution of planetesimals in extreme environments. Populations in this region include the Atira and Vatira asteroids, as well as Venus and Earth co-orbital asteroids. Here we present a twilight search for these small bodies, conducted using the 1.2-m Oschin Schmidt and the Zwicky Transient Facility (ZTF) camera at Palomar Observatory. The ZTF twilight survey operates at solar elongations down to $35^\circ$ with limiting magnitude of $r=19.5$. During a total of 40 evening sessions and 62 morning sessions conducted between 2018 November 15 and 2019 June 23, we detected 6 Atiras, including 2 new discoveries 2019 AQ$_3$ and 2019 LF$_6$, but no Vatiras or Earth/Venus co-orbital asteroids. NEO population models show that these new discoveries are likely only the tip of the iceberg, with the bulk of the population yet to be found. The population models also suggest that we have only detected 5--$7\%$ of the $H<20$ Atira population over the 7-month survey. Co-orbital asteroids are smaller in diameters and require deeper surveys. A systematic and efficient survey of the near-Sun region will require deeper searches and/or facilities that can operate at small solar elongations.

The Pioneer 10 and 11 and Voyager 1 and 2 spacecraft, launched in the 1970s, are heading out of the solar system. Using the astrometric and radial velocity data from the second Gaia data release, we integrate the trajectories of 7.4 million stars, and the spacecraft, through a Galactic potential in order to identify those stars the spacecraft will pass closest to. The closest encounters for all spacecraft take place at separations between 0.2 and 0.5 pc within the next million years. The closest encounter will be by Pioneer 10 with the K8 dwarf HIP 117795, at 0.23 pc in 90 kyr at a high relative velocity of 291 km/s.

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.

Near-Earth asteroid (410777) 2009 FD is a potentially hazardous asteroid with potential impacts on Earth at the end of the 22nd century. The astrometry collected during the 2019 apparition provides information on the trajectory of (410777) by constraining the Yarkovsky effect, which is the main source of uncertainty for future predictions, and informing the impact hazard assessment. We included the Yarkovsky effect in the force model and estimated its magnitude from the fit to the (410777) optical and radar astrometric data. We performed the (410777) hazard assessment over 200 years by using two independent approaches: the NEODyS group adopted a generalisation of the Line Of Variations method in a 7-dimensional space, whereas the JPL team resorted to the Multilayer Clustered Sampling technique. We obtain a 4-sigma detection of the Yarkovsky effect acting on (410777), which corresponds to a semimajor axis drift of (3.6 +/- 0.9) x 10^{-3} au/Myr. As for the hazard results of both teams, the main impact possibility in 2185 is ruled out and the only remaining one is in 2190, but with a probability smaller than 10^{-8}.

Though pairs of dynamically associated asteroids in the Main Belt have been identified and studied for over a decade, very few pair systems have been identified in the near-Earth asteroid population. We present data and analysis that supports the existence of two genetically related pairs in near-Earth space. The members of the individual systems, 2015 EE7 -- 2015 FP124 and 2017 SN16 -- 2018 RY7, are found to be of the same spectral taxonomic class, and both pairs are interpreted to have volatile-poor compositions. In conjunction with dynamical arguments, this suggests that these two systems formed via YORP spin-up and/or dissociation of a binary precursor. Backwards orbital integrations suggest a separation age of <10 kyr for the pair 2017 SN16 -- 2018 RY7, making these objects amongst the youngest multiple asteroid systems known to date. A unique separation age was not realized for 2015 EE7 -- 2015 FP124 due to large uncertainties associated with these objects' orbits. Determining the ages of such young pairs is of great value for testing models of space weathering and asteroid spin-state evolution. As the NEO catalog continues to grow with current and future discovery surveys, it is expected that more NEO pairs will be found, thus providing an ideal laboratory for studying time dependent evolutionary processes that are relevant to asteroids throughout the Solar System.

We present a photometric and dynamical study of comet C/2010 U3 (Boattini), which was seen active in prediscovery data as early as 2005 November at a new inbound record heliocentric distance $r_{\rm H} = 25.8$ au. Two outburst events around 2009 and 2017 were observed. The coma and tail of the comet consist of dust grains of $\sim$10 $\mu$m in radius, ejected protractedly at speeds $\lesssim$50 m s$^{-1}$ near the subsolar point, and are subjected to the Lorentz force, solar gravitation and radiation pressure force altogether. The prolonged activity indicates that sublimation of supervolatiles (e.g., CO, CO$_2$) is at play, causing a net mass-loss rate $\gtrsim$1 kg s$^{-1}$. To sustain the mass loss, the nucleus radius has to be $\gtrsim$0.1 km. The color of the cometary dust, similar to other long-period comets, is redder than the solar colors, but we also observed potential color variations when the comet was at $10 < r_{\rm H} < 15$ au, concurrent with the onset of crystallisation of amorphous water ice, if at all. Using publicly available and our refined astrometric measurements, we estimated the precise trajectory of the comet and propagated it backward to its previous perihelion. We found that the comet has visited the planetary region $1.96 \pm 0.04$ Myr ago, with barycentric perihelion distance $q = 8.364 \pm 0.004$ au. Thus, C/2010 U3 (Boattini) is almost certainly a dynamically old comet from the Oort cloud, and the observed activity cannot be caused by retained heat from the previous apparition.

The recent close approach of the NEA (3200) Phaethon offered a rare opportunity to obtain high-quality observational data. We used the newly obtained optical light curves to improve the spin and shape model of Phaethon and to determine its surface physical properties derived by thermophysical modeling. We also used the available astrometric observations of Phaethon, including those obtained by the Arecibo radar and the Gaia spacecraft, to constrain the secular drift of the orbital semimajor axis. This constraint allowed us to estimate the bulk density by assuming that the drift is dominated by the Yarkovsky effect. We used the convex inversion model to derive the 3D shape model of Phaethon, and a detailed numerical approach for an accurate analysis of the Yarkovsky effect. We obtained a unique solution for Phaethon's pole orientation at $(318,-47)^{\circ}$ ecliptic longitude and latitude (uncertainty of $5^{\circ}$), and confirm the previously reported thermophysical properties ($D=5.1\pm0.2$ km, $\Gamma=600\pm200$ SI). Phaethon has a top-like shape with possible north-south asymmetry. The characteristic size of the regolith grains is 1-2 cm. The orbit analysis reveals a secular drift of the semimajor axis of $-(6.9\pm1.9)\times 10^{-4}$ au Myr$^{-1}$. With the derived volume-equivalent size of 5.1~km, the bulk density $\rho$ is $1.67\pm0.47$ g cm$^{-3}$. If the size is slightly larger $\sim5.7$ km, as suggested by radar data, $\rho$ would decrease to $1.48\pm0.42$ g cm$^{-3}$. We further investigated the suggestion that Phaethon may be in a cluster with asteroids (155140) 2005 UD and (225416) 1999 YC that was formed by rotational fission of a critically spinning parent body. Phaethon's $\rho$ is consistent with typical values for large ($>100$ km) C-complex asteroids and supports its association with asteroid (2) Pallas. These findings render a cometary origin unlikely for Phaethon.

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.

The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity.

We studied the Near-Earth Asteroid (NEA) candidates posted on the Minor Planet Center's Near-Earth Object Confirmation Page (NEOCP) between years 2013 and 2016. Out of more than 17,000 NEA candidates, while the majority became either new discoveries or were associated with previously known objects, about 11% were unable to be followed-up or confirmed. We further demonstrate that of the unconfirmed candidates, 926+/-50 are likely to be NEAs, representing 18% of discovered NEAs in that period. Only 11% (~93) of the unconfirmed NEA candidates were large (having absolute magnitude H<22). To identify the reasons why these NEAs were not recovered, we analyzed those from the most prolific asteroid surveys: Pan-STARRS, the Catalina Sky Survey, the Dark Energy Survey, and the Space Surveillance Telescope. We examined the influence of plane-of-sky positions and rates of motion, brightnesses, submission delays, and computed absolute magnitudes, as well as correlations with the phase of the moon and seasonal effects. We find that delayed submission of newly discovered NEA candidate to the NEOCP drove a large fraction of the unconfirmed NEA candidates. A high rate of motion was another significant contributing factor. We suggest that prompt submission of suspected NEA discoveries and rapid response to fast moving targets and targets with fast growing ephemeris uncertainty would allow better coordination among dedicated follow-up observers, decrease the number of unconfirmed NEA candidates, and increase the discovery rate of NEAs.

We present an updated set of near-Earth asteroids with a Yarkovsky-related semimajor axis drift detected from the orbital fit to the astrometry. We find 87 reliable detections after filtering for the signal-to-noise ratio of the Yarkovsky drift estimate and making sure the estimate is compatible with the physical properties of the analyzed object. Furthermore, we find a list of 24 marginally significant detections, for which future astrometry could result in a Yarkovsky detection. A further outcome of the filtering procedure is a list of detections that we consider spurious because unrealistic or not explicable with the Yarkovsky effect. Among the smallest asteroids of our sample, we determined four detections of solar radiation pressure, in addition to the Yarkovsky effect. As the data volume increases in the near future, our goal is to develop methods to generate very long lists of asteroids with reliably detected Yarkovsky effect, with limited amounts of case by case specific adjustments. Furthermore, we discuss the improvements this work could bring to impact monitoring. In particular, we exhibit two asteroids for which the adoption of a non-gravitational model is needed to make reliable impact predictions.

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

We imaged Comet 252P/2000 G1 (LINEAR) (hereafter 252P) with the Hubble Space Telescope and both 252P and P/2016 BA$_{14}$ (PanSTARRS) (hereafter BA$_{14}$) with the Discovery Channel Telescope in March and April 2016, surrounding its close encounter to Earth. The r'-band $Af\rho$ of 252P in a 0.2"-radius aperture were $16.8\pm0.3$ and $57\pm1$ cm on March 14 and April 4, respectively, and its gas production rates were: $Q$(OH) = $(5.8\pm0.1)\times10^{27}$ s$^{-1}$, and $Q$(CN) = $(1.25\pm0.01)\times10^{25}$ s$^{-1}$ on April 17. The r'-band upper limit $Af\rho$ of BA1$_{14}$ was $0.19\pm0.01$ cm in a 19.2"-radius aperture, and $Q$(CN) = $(1.4\pm0.1)10^{22}$ s$^{-1}$ on April 17, 2017. 252P shows a bright and narrow jet of a few hundred kilometers long in the sunward direction, changing its projected position angle in the sky with a periodicity consistent with 7.24 hours. However, its photometric lightcurve is consistent with a periodicity of 5.41 hours. We suggest that the nucleus of 252P is likely in a non-principal axis rotation. The nucleus radius of 252P is estimated to be about $0.3\pm0.03$ km, indicating an active fraction of 40% to >100% in its 2016 apparition. Evidence implies a possible cloud of slow-moving grains surrounding the nucleus. The activity level of 252P in the 2016 apparition increased by two orders of magnitude from its previous apparitions, making this apparition unusual. On the other hand, the activity level of BA14 appears to be at least three orders of magnitude lower than that of 252P, despite its ten times or larger surface area.

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.

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

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.

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.

Although the nucleus of comet C/2013 A1 (Siding Spring) will safely pass Mars in October 2014, the dust in the coma and tail will more closely approach the planet. Using a dynamical model of comet dust, we estimate the impact fluence. Based on our nominal model no impacts are expected at Mars. Relaxing our nominal model's parameters, the fluence is no greater than ~10^-7 grains/m^2 for grain radii larger than 10 {\mu}m. Mars orbiting spacecraft are unlikely to be impacted by large dust grains, but Mars may receive as many as ~10^7 grains, or ~100 kg of total dust. We also estimate the flux of impacting gas molecules commonly observed in comet comae.

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 have a close encounter with Mars on October 19, 2014. We model the dynamical evolution of dust grains from the time of their ejection from the comet nucleus to the Mars close encounter, and determine the flux at Mars. Constraints on the ejection velocity from Hubble Space Telescope observations indicate that the bulk of the grains will likely miss Mars, although it is possible that a few-percent of grains with higher velocities will reach Mars, peaking approximately 90--100 minutes after the close approach of the nucleus, and consisting mostly of millimeter-radius grains ejected from the comet nucleus at a heliocentric distance of approximately 9~AU or larger. At higher velocities, younger grains from sub-millimeter to several millimeter can reach Mars too, although an even smaller fraction of grains is expected have these velocities, with negligible effect on the peak timing. Using NEOWISE observations of the comet, we can estimate that the maximum fluence will be of the order of $10^{-7}$ grains/m$^2$. We include a detailed analysis of how the expected fluence depends on the grain density, ejection velocity, and size-frequency distribution, to account for current model uncertainties and in preparation of possible refined model values in the near future.

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.

We study the time evolution of the impact probability for synthetic but realistic impacting and close approaching asteroids detected in a simulated all-sky survey. We use the impact probability to calculate the impact warning time as the time interval between when an object reaches a Palermo Scale value of -2 and when it impacts Earth. A simple argument shows that warning time ~ D^x with the exponent in the range [1.0,1.5] and our derived value was $x=1.3\pm0.1$. The low-precision astrometry from the single simulated all-sky survey could require many days or weeks to establish an imminent impact for asteroids >100m diameter that are discovered far from Earth. Most close approaching asteroids are quickly identified as not being impactors but a size-dependent percentage, even for those >50m diameter, have a persistent impact probability of >10^{-6} on the day of closest approach. Thus, a single all-sky survey can be of tremendous value in identifying Earth impacting and close approaching asteroids in advance of their closest approach but it can not solve the problem on its own: high-precision astrometry from other optical or radar systems is necessary to rapidly establish an object as an impactor or close approacher. We show that the parallax afforded by surveying the sky from two sites is only of benefit for a small fraction of the smallest objects detected within a couple days before impact: probably not enough to justify the increased operating costs of a 2-site survey. Finally, the survey cadence within a fixed time span is relatively unimportant to the impact probability calculation. We tested three different reasonable cadences and found that one provided ~10x higher (better) value for the impact probability on the discovery night for the smallest (10m diameter) objects but the consequences on the overall impact probability calculation is negligible.

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.

The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1,000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.

We consider a network of telescopes capable of scanning all the observable sky each night and targeting Near-Earth objects (NEOs) in the size range of the Tunguska-like asteroids, from 160 m down to 10 m. We measure the performance of this telescope network in terms of the time needed to discover at least 50% of the impactors in the considered population with a warning time large enough to undertake proper mitigation actions. The warning times are described by a trimodal distribution and the telescope network has a 50% probability of discovering an impactor of the Tunguska class with at least one week of advance already in the first 10 yr of operations of the survey. These results suggest that the studied survey would be a significant addition to the current NEO discovery efforts.

We present the results of a large scale simulation, reproducing the behavior of a data center for the build-up and maintenance of a complete catalog of space debris in the upper part of the low Earth orbits region (LEO). The purpose is to determine the performances of a network of advanced optical sensors, through the use of the newest orbit determination algorithms developed by the Department of Mathematics of Pisa (DM). Such a network has been proposed to ESA in the Space Situational Awareness (SSA) framework by Carlo Gavazzi Space SpA (CGS), Istituto Nazionale di Astrofisica (INAF), DM, and Istituto di Scienza e Tecnologie dell'Informazione (ISTI-CNR). The conclusion is that it is possible to use a network of optical sensors to build up a catalog containing more than 98% of the objects with perigee height between 1100 and 2000 km, which would be observable by a reference radar system selected as comparison. It is also possible to maintain such a catalog within the accuracy requirements motivated by collision avoidance, and to detect catastrophic fragmentation events. However, such results depend upon specific assumptions on the sensor and on the software technologies.

Given a set of astrometric observations of the same object, the problem of orbit determination is to compute the orbit and to assess its uncertainty and reliability. For the next generation surveys, with much larger number density of observed objects, new algorithms or substantial revisions of the classical ones are needed. The problem has three main steps, preliminary orbit, least squares orbit, and quality control. The classical theory of preliminary orbits was incomplete: the consequences of the topocentric correction had not been fully studied. We show that it is possible to account for this correction, possibly with an increase in the number of preliminary solutions, without impairing the overall orbit determination performance. We have developed modified least squares orbit determination algorithms that can be used to improve the reliability of the procedure. We have tested the complete procedure on two simulations with number densities comparable to that expected from the next generation surveys such as Pan-STARRS and LSST. To control the problem of false identifications we have introduced a quality control on the fit residuals based on an array of metrics and a procedure to remove duplications and contradictions in the output. The results confirm that large sets of discoveries can be obtained with good quality orbits and very high success rate losing only 0.6 to 1.3% of objects and a false identification rate in the range 0.02 to 0.06%.