Showing 1–50 of 60 results for author: Farnocchia, D

The Dynamical Origins of the Dark Comets and a Proposed Evolutionary Track

Authors: Aster G. Taylor, Jordan K. Steckloff, Darryl Z. Seligman, Davide Farnocchia, Luke Dones, David Vokrouhlicky, David Nesvorny, Marco Micheli

Abstract: So-called 'dark comets' are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubbl… ▽ More So-called 'dark comets' are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the $ν_6$ resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 17 pages, 7 figures, 12 supplementary figures & pages. Accepted for publication in Icarus

NEOMOD 3: The Debiased Size Distribution of Near Earth Objects

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Carson Fuls, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Marco Delbo

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

Submitted 29 April, 2024; originally announced April 2024.

Comments: Icarus, in press

NEOMOD 2: An Updated Model of Near-Earth Objects from a Decade of Catalina Sky Survey Observations

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Eric Christensen, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Mikael Granvik

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

Submitted 14 December, 2023; originally announced December 2023.

Comments: Icarus, in press. arXiv admin note: text overlap with arXiv:2306.09521

Seasonally Varying Outgassing as an Explanation for Dark Comet Accelerations

Authors: Aster G. Taylor, Davide Farnocchia, David Vokrouhlicky, Darryl Z. Seligman, Jordan K. Steckloff, Marco Micheli

Abstract: Significant nonradial, nongravitational accelerations with magnitudes incompatible with radiation-driven effects have been reported in seven small, photometrically inactive near-Earth objects. Two of these objects exhibit large transverse accelerations (i.e., within the orbital plane but orthogonal to the radial direction), and six exhibit significant out-of-plane accelerations. Here, we find that… ▽ More Significant nonradial, nongravitational accelerations with magnitudes incompatible with radiation-driven effects have been reported in seven small, photometrically inactive near-Earth objects. Two of these objects exhibit large transverse accelerations (i.e., within the orbital plane but orthogonal to the radial direction), and six exhibit significant out-of-plane accelerations. Here, we find that anisotropic outgassing resulting from differential heating on a nucleus with nonzero spin-pole obliquity, averaged over an eccentric orbit, can explain these accelerations for most of the objects. This balanced outgassing model depends on three parameters -- the spin pole orientation (R.A. and Dec.) and an acceleration magnitude. For these "dark comets" (excepting 2003 RM), we obtain parameter values that reproduce the observed nongravitational accelerations. We derive formulae for the component accelerations under certain assumptions for the acceleration scaling over heliocentric distance. Although we lack estimates of these objects' spin axes to confirm our values, this mechanism is nevertheless a plausible explanation for the observed accelerations, and produces accurate perturbations to the heliocentric motions of most of these objects. This model may also be applied to active objects outside of the dark comets group. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 17 pages, 4 figures. Published in Icarus

Detectability of the Yarkovsky Effect in the Main Belt

Authors: Denise Hung, David J. Tholen, Davide Farnocchia, Federica Spoto

Abstract: We attempt to a detect signal of Yarkovsky-related acceleration in the orbits of 134 main belt asteroids (MBAs) we observed with the University of Hawai'i 88 inch telescope, supplemented with observations publicly available from the Minor Planet Center and Gaia Data Release 3. We estimated expected Yarkovsky acceleration values based on parameters derived through thermophysical modeling, but we we… ▽ More We attempt to a detect signal of Yarkovsky-related acceleration in the orbits of 134 main belt asteroids (MBAs) we observed with the University of Hawai'i 88 inch telescope, supplemented with observations publicly available from the Minor Planet Center and Gaia Data Release 3. We estimated expected Yarkovsky acceleration values based on parameters derived through thermophysical modeling, but we were not able to find any reliable detections of Yarkovsky in our sample. Through tests with synthetic observations however, we estimated the minimum observational arc length needed to detect the Yarkovsky effect for all of our sample MBAs, which in nearly every case exceeded the current arc length of the existing observations. We find that the Yarkovsky effect could be detectable within a couple of decades of discovery for a 100 m MBA assuming 0.1" astrometric accuracy, which is at the size range detectable by the upcoming Vera Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 2 November, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 28 pages, 6 figures, 6 tables, accepted by PSJ

Constraints on fifth forces and ultralight dark matter from OSIRIS-REx target asteroid Bennu

Authors: Yu-Dai Tsai, Davide Farnocchia, Marco Micheli, Sunny Vagnozzi, Luca Visinelli

Abstract: Using the OSIRIS-REx mission and ground-based tracking data for the asteroid Bennu, we derive new constraints on fifth forces and ultralight dark matter. The bounds we obtain are strongest for mediator masses $m \sim 10^{-18} - 10^{-17}\,{\rm eV}$, where we currently achieve the tightest bounds. Our limits can be translated to a wide class of models leading to Yukawa-type fifth forces, and we demo… ▽ More Using the OSIRIS-REx mission and ground-based tracking data for the asteroid Bennu, we derive new constraints on fifth forces and ultralight dark matter. The bounds we obtain are strongest for mediator masses $m \sim 10^{-18} - 10^{-17}\,{\rm eV}$, where we currently achieve the tightest bounds. Our limits can be translated to a wide class of models leading to Yukawa-type fifth forces, and we demonstrate how they apply to $U(1)_B$ dark photons and baryon-coupled scalars. Our results demonstrate the potential of asteroid tracking in probing well-motivated extensions of the Standard Model and ultralight dark matter satisfying the fuzzy dark matter constraints. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 1 figure, 7 pages plus references

Report number: UCI-HEP-TR-2023-04; FERMILAB-PUB-23-538-T-V; CA21106; CA21136

Interstellar Object Uncertainty Evolution and Effect on Fast Flyby Delivery and Required Delta-V

Authors: Declan Mages, Davide Farnocchia, Benjamin Donitz

Abstract: Interstellar objects (ISOs) are small bodies that can travel through our solar system from other star systems. When present in our solar system, they represent an opportunity to study the properties and origins of these objects, as well as the potential for cross-pollination of material between star systems. With current propulsion technology, rendezvous with these objects is likely infeasible, an… ▽ More Interstellar objects (ISOs) are small bodies that can travel through our solar system from other star systems. When present in our solar system, they represent an opportunity to study the properties and origins of these objects, as well as the potential for cross-pollination of material between star systems. With current propulsion technology, rendezvous with these objects is likely infeasible, and thus the maximum science return results from a rapid response flyby and impactor. However, while trajectories to ISOs may be feasible, their potentially high ephemeris uncertainties and high-speed hyperbolic orbits present significant challenges to navigation. In this paper we assess these challenges by modeling the uncertainties of reachable synthetic ISOs as a function of time, as derived by measurements from ground observatories and an approaching spacecraft. From these uncertainties we derive the final delivery accuracy of fast flyby spacecraft to the ISO and required statistical delta-v for navigation. We find that these two challenges can lead to hundreds of meters-per-second or even kilometers-per-second of required statistical delta-v for navigation, reduce delivery accuracy to hundreds of kilometers, and make autonomous navigation a requirement. △ Less

Submitted 21 July, 2023; originally announced July 2023.

The hazardous km-sized NEOs of the next thousands of years

Authors: Oscar Fuentes-Muñoz, Daniel J. Scheeres, Davide Farnocchia, Ryan S. Park

Abstract: The catalog of km-sized near-Earth objects (NEOs) is nearly complete. Typical impact monitoring analyses search for possible impacts over the next 100 years and none of the km-sized objects represent an impact threat over that time interval. Assessing the impact risk over longer time scales is a challenge since orbital uncertainties grow. To overcome this limitation we analyze the evolution of the… ▽ More The catalog of km-sized near-Earth objects (NEOs) is nearly complete. Typical impact monitoring analyses search for possible impacts over the next 100 years and none of the km-sized objects represent an impact threat over that time interval. Assessing the impact risk over longer time scales is a challenge since orbital uncertainties grow. To overcome this limitation we analyze the evolution of the Minimum Orbit Intersection Distance (MOID), which bounds the closest possible encounters between the asteroid and the Earth. The evolution of the MOID highlights NEOs that are in the vicinity of the Earth for longer periods of time, and we propose a method to estimate the probability of a deep Earth encounter during these periods. This metric is used to rank the km-sized catalog in terms of their long-term impact hazard to identify targets of potential interest for additional observation and exploration. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 20 pages, 17 figures, Accepted for Publication in The Astronomical Journal

ASSIST: An Ephemeris-Quality Test Particle Integrator

Authors: Matthew J. Holman, Arya Akmal, Davide Farnocchia, Hanno Rein, Matthew J. Payne, Robert Weryk, Daniel Tamayo, David M. Hernandez

Abstract: We introduce ASSIST, a software package for ephemeris-quality integrations of test particles. ASSIST is an extension of the REBOUND framework and makes use of its IAS15 integrator to integrate test particle trajectories in the field of the Sun, Moon, planets, and 16 massive asteroids, with the positions of the masses coming from the JPL DE441 ephemeris and its associated asteroid perturber file. T… ▽ More We introduce ASSIST, a software package for ephemeris-quality integrations of test particles. ASSIST is an extension of the REBOUND framework and makes use of its IAS15 integrator to integrate test particle trajectories in the field of the Sun, Moon, planets, and 16 massive asteroids, with the positions of the masses coming from the JPL DE441 ephemeris and its associated asteroid perturber file. The package incorporates the most significant gravitational harmonics and general relativistic corrections. ASSIST also accounts for position- and velocity-dependent non-gravitational effects. The first order variational equations are included for all terms to support orbit fitting and covariance mapping. This new framework is meant to provide an open-source package written in a modern language to enable high-precision orbital analysis and science by the small body community. ASSIST is open source, freely distributed under the GNU General Public license, version 3. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to PSJ (27 January 2023), Revised (28 March 2023)

(523599) 2003 RM: The asteroid that wanted to be a comet

Authors: Davide Farnocchia, Darryl Z. Seligman, Mikael Granvik, Olivier Hainaut, Karen J. Meech, Marco Micheli, Robert Weryk, Steven R. Chesley, Eric J. Christensen, Detlef Koschny, Jan T. Kleyna, Daniela Lazzaro, Michael Mommert, Richard Wainscoat

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

Submitted 17 January, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 14 pages, 9 figures, 2 tables, accepted by the Planetary Science Journal

Dark Comets? Unexpectedly Large Nongravitational Accelerations on a Sample of Small Asteroids

Authors: Darryl Z. Seligman, Davide Farnocchia, Marco Micheli, David Vokrouhlický, Aster G. Taylor, Steven R. Chesley, Jennifer B. Bergner, Peter Vereš, Olivier R. Hainaut, Karen J. Meech, Maxime Devogele, Petr Pravec, Rob Matson, Sam Deen, David J. Tholen, Robert Weryk, Edgard G. Rivera-Valentín, Benjamin N. L. Sharkey

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

Submitted 17 January, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 21 pages, 9 figures, 3 tables, accepted for publication in the Planetary Science Journal

Anticipating the DART impact: Orbit estimation of Dimorphos using a simplified model

Authors: Shantanu P. Naidu, Steven R. Chesley, Davide Farnocchia, Nick Moskovitz, Petr Pravec, Petr Scheirich, Cristina Thomas, Andrew S. Rivkin

Abstract: 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 cau… ▽ More 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σ$) 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σ$). The formal $3σ$ 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. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: Accepted to PSJ

New Constraints on Dark Matter and Cosmic Neutrino Profiles through Gravity

Authors: Yu-Dai Tsai, Joshua Eby, Jason Arakawa, Davide Farnocchia, Marianna S. Safronova

Abstract: We derive purely gravitational constraints on dark matter and cosmic neutrino profiles in the solar system using asteroid (101955) Bennu. We focus on Bennu because of its extensive tracking data and high-fidelity trajectory modeling resulting from the OSIRIS-REx mission. We find that the local density of dark matter is bound by… ▽ More We derive purely gravitational constraints on dark matter and cosmic neutrino profiles in the solar system using asteroid (101955) Bennu. We focus on Bennu because of its extensive tracking data and high-fidelity trajectory modeling resulting from the OSIRIS-REx mission. We find that the local density of dark matter is bound by $ρ_{\rm DM}\lesssim 3.3\times 10^{-15}\;\rm kg/m^3 \simeq 6\times10^6\,\barρ_{\rm DM}$, in the vicinity of $\sim 1.1$ au (where $\barρ_{\rm DM}\simeq 0.3\;\rm GeV/cm^3$). We show that high-precision tracking data of solar system objects can constrain cosmic neutrino overdensities relative to the Standard Model prediction $\bar{n}_ν$, at the level of $η\equiv n_ν/\bar{n}_ν\lesssim 1.7 \times 10^{11}(0.1 \;{\rm eV}/m_ν)$ (Saturn), comparable to the existing bounds from KATRIN and other previous laboratory experiments (with $m_ν$ the neutrino mass). These local bounds have interesting implications for existing and future direct-detection experiments. Our constraints apply to all dark matter candidates but are particularly meaningful for scenarios including solar halos, stellar basins, and axion miniclusters, which predict or allow overdensities in the solar system. Furthermore, introducing a DM-SM long-range fifth force with a strength $\tildeα_D$ times stronger than gravity, Bennu can set a constraint on $ρ_{\rm DM}\lesssim \barρ_{\rm DM}\left(6 \times 10^6/\tildeα_D\right)$. These constraints can be improved in the future as the accuracy of tracking data improves, observational arcs increase, and more missions visit asteroids. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 2 figures, 2 tables, 9 pages plus references

Report number: UCI-HEP-TR-2022-11

Journal ref: JCAP 02 (2024) 029

Modelling the water and carbon dioxide production rates of Comet 67P/Churyumov-Gerasimenko

Authors: Björn J. R. Davidsson, Nalin H. Samarasinha, Davide Farnocchia, Pedro J. Gutiérrez

Abstract: The European Space Agency Rosetta/Philae mission to Comet 67P/Churyumov-Gerasimenko in 2014-2016 is the most complete and diverse investigation of a comet carried out thus far. Yet, many physical and chemical properties of the comet remain uncertain or unknown, and cometary activity is still not a well-understood phenomenon. We here attempt to place constraints on the nucleus abundances and sublim… ▽ More The European Space Agency Rosetta/Philae mission to Comet 67P/Churyumov-Gerasimenko in 2014-2016 is the most complete and diverse investigation of a comet carried out thus far. Yet, many physical and chemical properties of the comet remain uncertain or unknown, and cometary activity is still not a well-understood phenomenon. We here attempt to place constraints on the nucleus abundances and sublimation front depths of H2O and CO2 ice, and to reconstruct how the nucleus evolved throughout the perihelion passage. We employ the thermophysical modelling code 'Numerical Icy Minor Body evolUtion Simulator', or NIMBUS, to search for conditions under which the observed H2O and CO2 production rates are simultaneously reproduced before and after perihelion. We find that the refractories to water-ice mass ratio of relatively pristine nucleus material is mu~1, that airfall material has mu~2, and that the molar abundance of CO2 relative H2O near 30 per cent. The dust mantle thickness is typically < 2 cm. The average CO2 sublimation front depths near aphelion were ~3.8 m and ~1.9 m on the northern and southern hemispheres, respectively, but varied substantially with time. We propose that airfall material is subjected to substantial fragmentation and pulverisation due to thermal fatigue during the aphelion passage. Sub-surface compaction of material due to CO2 activity near perihelion seems to have reduced the diffusivity in a measurable way. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 22 pages, 19 figures. This is a pre--copyedited, author--produced PDF of an article accepted for publication in MNRAS following peer review

Journal ref: Mon. Not. R. Astron. Soc. 509, 3065-3085 (2022)

Internal rubble properties of asteroid (101955) Bennu

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

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

Submitted 23 August, 2021; originally announced August 2021.

Comments: 14 pages, 5 figures, accepted to Icarus

A novel approach to asteroid impact monitoring and hazard assessment

Authors: Javier Roa, Davide Farnocchia, Steven R. Chesley

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

Submitted 6 August, 2021; originally announced August 2021.

Comments: 19 pages, 13 figures. To be published in the Astronomical Journal

The impact and recovery of asteroid 2018 LA

Authors: Peter Jenniskens, Mohutsiwa Gabadirwe, Qing-Zhu Yin, Alexander Proyer, Oliver Moses, Tomas Kohout, Fulvio Franchi, Roger L. Gibson, Richard Kowalski, Eric J. Christensen, Alex R. Gibbs, Aren Heinze, Larry Denneau, Davide Farnocchia, Paul W. Chodas, William Gray, Marco Micheli, Nick Moskovitz, Christopher A. Onken, Christian Wolf, Hadrien A. R. Devillepoix, Quanzhi Ye, Darrel K. Robertson, Peter Brown, Esko Lyytinen , et al. (41 additional authors not shown)

Abstract: The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to… ▽ More The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. 23 meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as a HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of about 156 cm diameter with high bulk density about 2.85 g/cm3, a relatively low albedo pV about 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of about 0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth-impacting orbit via the nu_6 resonance. The impact that ejected 2018 LA in an orbit towards Earth occurred 22.8 +/- 3.8 Ma ago. Zircons record a concordant U-Pb age of 4563 +/- 11 Ma and a consistent 207Pb/206Pb age of 4563 +/- 6 Ma. A much younger Pb-Pb phosphate resetting age of 4234 +/- 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: Meteoritics & Planetary Science (2021)

(6478) Gault: Physical characterization of an active main-belt asteroid

Authors: Maxime Devogèle, Marin Ferrais, Emmanuel Jehin, Nicholas Moskovitz, Brian A. Skiff, Stephen E. Levine, Annika Gustafsson, Davide Farnocchia, Marco Micheli, Colin Snodgrass, Galin Borisov, Jean Manfroid, Youssef Moulane, Zouhair Benkhaldoun, Artem Burdanov, Francisco J. Pozuelos, Michael Gillon, Julien de Wit, Simon F. Green, Philippe Bendjoya, Jean-Pierre Rivet, Luy Abe, David Vernet, Colin Orion Chandler, Chadwick A. Trujillo

Abstract: In December 2018, the main-belt asteroid (6478)~Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observation… ▽ More In December 2018, the main-belt asteroid (6478)~Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observations over two apparitions to monitor its activity, rotation period, composition, and possible non-gravitational orbital evolution. We find that Gault has a rotation period of $P = 2.4929 \pm 0.0003$ hours with a lightcurve amplitude of $0.06$ magnitude. This short rotation period close to the spin barrier limit is consistent with Gault having a density no smaller than $ρ= 1.85$~g/cm$^3$ and its activity being triggered by the YORP spin-up mechanism. Analysis of the Gault phase curve over phase angles ranging from $0.4^{\circ}$ to $23.6^{\circ}$ provides an absolute magnitude of $H = 14.81 \pm 0.04$, $G1=0.25 \pm 0.07$, and $G2= 0.38 \pm 0.04$. Model fits to the phase curve find the surface regolith grain size constrained between 100-500 $\rmμ$m. Using relations between the phase curve and albedo we determine that the geometrical albedo of Gault is $p_{\rm v} = 0.26 \pm 0.05$ corresponding to an equivalent diameter of $D = 2.8^{+0.4}_{-0.2}$ km. Our spectroscopic observations are all consistent with an ordinary chondrite-like composition (S, or Q-type in the Bus-DeMeo taxonomic classification). A search through archival photographic plate surveys found previously unidentified detections of Gault dating back to 1957 and 1958. Only the latter had been digitized, which we measured to nearly double the observation arc of Gault. Finally, we did not find any signal of activity during the 2020 apparition or non-gravitational effects on its orbit. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 15 pages, 11 figures; Accepted in MNRAS

Hayabusa2 Extended Mission: New Voyage to Rendezvous with a Small Asteroid Rotating with a Short Period

Authors: M. Hirabayashi, Y. Mimasu, N. Sakatani, S. Watanabe, Y. Tsuda, T. Saiki, S. Kikuchi, T. Kouyama, M. Yoshikawa, S. Tanaka, S. Nakazawa, Y. Takei, F. Terui, H. Takeuchi, A. Fujii, T. Iwata, K. Tsumura, S. Matsuura, Y. Shimaki, S. Urakawa, Y. Ishibashi, S. Hasegawa, M. Ishiguro, D. Kuroda, S. Okumura , et al. (19 additional authors not shown)

Abstract: Hayabusa2 is the Japanese Asteroid Return Mission and targeted the carbonaceous asteroid Ryugu, conducted by the Japan Aerospace Exploration Agency (JAXA). The goal of this mission was to conduct proximity operations including remote sensing observations, material sampling, and a Small Carry-On Impact experiment, as well as sample analyses. As of September 2020, the spacecraft is on the way back t… ▽ More Hayabusa2 is the Japanese Asteroid Return Mission and targeted the carbonaceous asteroid Ryugu, conducted by the Japan Aerospace Exploration Agency (JAXA). The goal of this mission was to conduct proximity operations including remote sensing observations, material sampling, and a Small Carry-On Impact experiment, as well as sample analyses. As of September 2020, the spacecraft is on the way back to Earth with samples from Ryugu with no critical issues after the successful departure in November 2019. Here, we propose an extended mission in which the spacecraft will rendezvous with a small asteroid with ~30 m - ~40 m in diameter that is rotating at a spin period of ~10 min after an additional ~10-year cruise phase. We introduce that two scenarios are suitable for the extended mission. In the first scenario, the spacecraft will perform swing-by maneuvers at Venus once and Earth twice to arrive at asteroid 2001 AV43. In the second scenario, it will perform swing-by maneuvers at Earth twice to reach asteroid 1998 KY26. In both scenarios, the mission will continue until the early 2030s. JAXA recently released the decision that the spacecraft will rendezvous with 1998 KY26. This paper focuses on our scientific assessments of the two scenarios but leaves the decision process to go to 1998 KY26 for future reports. Rendezvous operations will be planned to detail the physical properties and surrounding environments of the target, one of the smallest elements of small planetary bodies. By achieving the planned operations, the mission will provide critical hints on the violent histories of collisions and accumulations of small bodies in the solar system. Furthermore, the established scientific knowledge and techniques will advance key technologies for planetary defense. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Comments: 51 pages, 9 figures, and 10 tables. The article is in press in Advances in Space Research

FUV Observations of the Inner Coma of 46P/Wirtanen

Authors: John W. Noonan, Walter M. Harris, Steven Bromley, Davide Farnocchia, Jian-Yang Li, Kathleen E. Mandt, Joel Wm. Parker, Kumar Venkataramani, Dennis Bodewits

Abstract: Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical componen… ▽ More Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical components for proper interpretation and analysis of comets. The close approach of comet 46P/Wirtanen in 2018-2019 provided a unique opportunity to study the inner most parts of a cometary coma with the Hubble Space Telescope Cosmic Origins Spectrograph, rarely accessible with remote observations, at length scales (100's of km) and wavelengths (900-1430 Angstroms) previously probed only by the European Space Agency's Rosetta spacecraft. Our observations show a complex picture for the inner coma; atomic production rates for H and O that show water is the dominant source of both, an abundance of atomic sulfur that is difficult to explain with the lifetimes of common sulfur parent molecules, and a density distribution that is poorly fit with both Haser and vectorial models. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 16 pages, 9 figures. Accepted in A'Hearn Special Issue of Planetary Science Journal

Impact Probability Under Aleatory And Epistemic Uncertainties

Authors: C. Tardioli, D. Farnocchia, M. Vasile, S. R. Chesley

Abstract: We present an approach to estimate an upper bound for the impact probability of a potentially hazardous asteroid when part of the force model depends on unknown parameters whose statistical distribution needs to be assumed. As case study we consider Apophis' risk assessment for the 2036 and 2068 keyholes based on information available as of 2013. Within the framework of epistemic uncertainties, un… ▽ More We present an approach to estimate an upper bound for the impact probability of a potentially hazardous asteroid when part of the force model depends on unknown parameters whose statistical distribution needs to be assumed. As case study we consider Apophis' risk assessment for the 2036 and 2068 keyholes based on information available as of 2013. Within the framework of epistemic uncertainties, under the independence and non-correlation assumption, we assign parametric families of distributions to the physical properties of Apophis that define the Yarkovsky perturbation and in turn the future orbital evolution of the asteroid. We find ${\rm IP}\leq 5\times 10^{-5}$ for the 2036 keyhole and ${\rm IP}\leq 1.6\times 10^{-5}$ for the 2068 keyhole. These upper bounds are largely conservative choices due to the rather wide range of statistical distributions that we explored. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Journal ref: Celestial Mechanics and Dynamical Astronomy (2020)

Establishing Earth's Minimoon Population through Characterization of Asteroid 2020 CD$_3$

Authors: Grigori Fedorets, Marco Micheli, Robert Jedicke, Shantanu P. Naidu, Davide Farnocchia, Mikael Granvik, Nicholas Moskovitz, Megan E. Schwamb, Robert Weryk, Kacper Wierzchoś, Eric Christensen, Theodore Pruyne, William F. Bottke, Quanzhi Ye, Richard Wainscoat, Maxime Devogèle, Laura E. Buchanan, Anlaug Amanda Djupvik, Daniel M. Faes, Dora Föhring, Joel Roediger, Tom Seccull, Adam B. Smith

Abstract: We report on our detailed characterization of Earth's second known temporary natural satellite, or minimoon, asteroid 2020 CD3. An artificial origin can be ruled out based on its area-to-mass ratio and broadband photometry, which suggest that it is a silicate asteroid belonging to the S or V complex in asteroid taxonomy. The discovery of 2020 CD3 allows for the first time a comparison between know… ▽ More We report on our detailed characterization of Earth's second known temporary natural satellite, or minimoon, asteroid 2020 CD3. An artificial origin can be ruled out based on its area-to-mass ratio and broadband photometry, which suggest that it is a silicate asteroid belonging to the S or V complex in asteroid taxonomy. The discovery of 2020 CD3 allows for the first time a comparison between known minimoons and theoretical models of their expected physical and dynamical properties. The estimated diameter of 1.2+0.4-0.2 m and geocentric capture approximately a decade after the first known minimoon, 2006 RH120, are in agreement with theoretical predictions. The capture duration of 2020 CD3 of at least 2.7 yr is unexpectedly long compared to the simulation average, but it is in agreement with simulated minimoons that have close lunar encounters, providing additional support for the orbital models. 2020 CD3's atypical rotation period, significantly longer than theoretical predictions, suggests that our understanding of meter-scale asteroids needs revision. More discoveries and a detailed characterization of the population can be expected with the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 22 pages, 6 Figures, 5 Tables, to appear in the Astronomical Journal

The carbon monoxide-rich interstellar comet 2I/Borisov

Authors: D. Bodewits, J. W. Noonan, P. D. Feldman, M. T. Bannister, D. Farnocchia, W. M. Harris, J. -Y. Li, K. E. Mandt, J. Wm. Parker, Z. Xing

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

Submitted 19 April, 2020; originally announced April 2020.

Journal ref: Nature Astronomy, 2020

A search for the origin of the interstellar comet 2I/Borisov

Authors: Coryn A. L. Bailer-Jones, Davide Farnocchia, Quanzhi Ye, Karen J. Meech, Marco Micheli

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

Submitted 7 January, 2020; v1 submitted 21 December, 2019; originally announced December 2019.

Comments: 7 pages; accepted to A&A

Journal ref: A&A 634, A14 (2020)

A Twilight Search for Atiras, Vatiras and Co-orbital Asteroids: Preliminary Results

Authors: Quanzhi Ye, Frank J. Masci, Wing-Huen Ip, Thomas A. Prince, George Helou, Davide Farnocchia, Eric C. Bellm, Richard Dekany, Matthew J. Graham, Shrinivas R. Kulkarni, Thomas Kupfer, Ashish Mahabal, Chow-Choong Ngeow, Daniel J. Reiley, Maayane T. Soumagnac

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

Submitted 12 December, 2019; originally announced December 2019.

Comments: AJ accepted

Future stellar flybys of the Voyager and Pioneer spacecraft

Authors: Coryn A. L. Bailer-Jones, Davide Farnocchia

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

Submitted 7 December, 2019; originally announced December 2019.

Comments: Extended version of our article published on 3 April 2019 in RNAAS

Journal ref: Research Notes of the American Astronomical Society, 3, 59 (2019)

Pre-discovery Activity of New Interstellar Comet 2I/Borisov Beyond 5 AU

Authors: Quanzhi Ye, Michael S. P. Kelley, Bryce T. Bolin, Dennis Bodewits, Davide Farnocchia, Frank J. Masci, Karen J. Meech, Marco Micheli, Robert Weryk, Eric C. Bellm, Eric Christensen, Richard Dekany, Alexandre Delacroix, Matthew J. Graham, Shrinivas R. Kulkarni, Russ R. Laher, Ben Rusholme, Roger M. Smith

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

Submitted 23 December, 2019; v1 submitted 13 November, 2019; originally announced November 2019.

Comments: AJ in press

Star catalog position and proper motion corrections in asteroid astrometry II: The Gaia era

Authors: Siegfried Eggl, Davide Farnocchia, Alan B. Chamberlin, Steven R. Chesley

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

Submitted 10 September, 2019; originally announced September 2019.

Comments: submitted to Icarus

Yarkovsky effect detection and updated impact hazard assessment for Near-Earth Asteroid (410777) 2009 FD

Authors: Alessio Del Vigna, Javier Roa, Davide Farnocchia, Marco Micheli, Dave Tholen, Francesca Guerra, Federica Spoto, Giovanni Battista Valsecchi

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

Submitted 10 July, 2019; v1 submitted 12 June, 2019; originally announced June 2019.

Comments: 8 pages, 2 figures

Journal ref: A&A 627, L11 (2019)

A common origin for dynamically associated near-Earth asteroid pairs

Authors: Nicholas Moskovitz, Petr Fatka, Davide Farnocchia, Maxime Devogele, David Polishook, Cristina A. Thomas, Michael Mommert, Louis D. Avner, Richard P. Binzel, Brian Burt, Eric Christensen, Francesca DeMeo, Mary Hinkle, Joseph L. Hora, Mitchell Magnusson, Robert Matson, Michael Person, Brian Skiff, Audrey Thirouin, David Trilling, Lawrence H. Wasserman, Mark Willman

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

Submitted 28 May, 2019; originally announced May 2019.

Comments: 8 figures, 2 tables; Accepted to Icarus

C/2010 U3 (Boattini): A Bizarre Comet Active at Record Heliocentric Distance

Authors: Man-To Hui, Davide Farnocchia, Marco Micheli

Abstract: 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 $μ$m in radius, ejected protractedly at speeds $\lesssim$50 m s… ▽ More 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 $μ$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. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Comments: Accepted by AJ

(3200) Phaethon: Bulk density from Yarkovsky drift detection

Authors: Josef Hanus, David Vokrouhlicky, Marco Delbo, Davide Farnocchia, David Polishook, Petr Pravec, Kamil Hornoch, Hana Kucakova, Peter Kusnirak, Robert Stephens, Brian Warner

Abstract: 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… ▽ More 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, $Γ=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 $ρ$ is $1.67\pm0.47$ g cm$^{-3}$. If the size is slightly larger $\sim5.7$ km, as suggested by radar data, $ρ$ 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 $ρ$ 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. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 620, L8 (2018)

Spitzer Observations of Interstellar Object 1I/`Oumuamua

Authors: DE Trilling, M Mommert, JL Hora, D Farnocchia, P Chodas, J Giorgini, HA Smith, S Carey, CM Lisse, M Werner, A McNeill, SR Chesley, JP Emery, G Fazio, YR Fernandez, A Harris, M Marengo, M Mueller, A Roegge, N Smith, HA Weaver, K Meech, M Micheli

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

Submitted 19 November, 2018; originally announced November 2018.

Journal ref: Published in the Astronomical Journal, 156, 261 (2018)

Plausible home stars of the interstellar object 'Oumuamua found in Gaia DR2

Authors: C. A. L. Bailer-Jones, D. Farnocchia, K. J. Meech, R. Brasser, M. Micheli, S. Chakrabarti, M. W. Buie, O. R. Hainaut

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

Submitted 24 September, 2018; originally announced September 2018.

Comments: Accepted to The Astronomical Journal

Detecting the Yarkovsky effect among near-Earth asteroids from astrometric data

Authors: Alessio Del Vigna, Laura Faggioli, Andrea Milani, Federica Spoto, Davide Farnocchia, Benoit Carry

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

Submitted 15 May, 2018; originally announced May 2018.

Comments: 36 pages, 1 figure

Journal ref: A&A 617, A61 (2018)

Unconfirmed Near-Earth Objects

Authors: Peter Vereš, Matthew J. Payne, Matthew J. Holman, Davide Farnocchia, Gareth V. Williams, Sonia Keys, Ian Boardman

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

Submitted 30 May, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 35 pages, 19 figures, 8 tables

Journal ref: 2018, Astronomical Journal, Vol 156, Number 1

YORP and Yarkovsky effects in asteroids (1685) Toro, (2100) Ra-Shalom, (3103) Eger, and (161989) Cacus

Authors: J. Durech, D. Vokrouhlicky, P. Pravec, J. Hanus, D. Farnocchia, Yu. N. Krugly, V. R. Ayvazian, P. Fatka, V. G. Chiorny, N. Gaftonyuk, A. Galad, R. Groom, K. Hornoch, R. Y. Inasaridze, H. Kucakova, P. Kusnirak, M. Lehky, O. I. Kvaratskhelia, G. Masi, I. E. Molotov, J. Oey, J. T. Pollock, V. G. Shevchenko, J. Vrastil, B. D. Warner

Abstract: 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… ▽ More 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} ω/ \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}ω/ \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σ$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}ω/ \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}ω/ \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3σ$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Journal ref: A&A 609, A86 (2018)

The Unusual Apparition of Comet 252P/2000 G1 (LINEAR) and Comparison with Comet P/2016 BA14 (PanSTARRS)

Authors: Jian-Yang Li, Michael S. P. Kelley, Nalin H. Samarasinha, Davide Farnocchia, Max J. Mutchler, Yanqiong Ren, Xiaoping Lu, David J. Tholen, Tim Lister, Marco Micheli

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

Submitted 17 August, 2017; originally announced August 2017.

Comments: 31 pages, 15 figures, 4 tables, accepted by AJ

Statistical Analysis of Astrometric Errors for the Most Productive Asteroid Surveys

Authors: Peter Vereš, Davide Farnocchia, Steven R. Chesley, Alan B. Chamberlin

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

Submitted 1 June, 2017; v1 submitted 9 March, 2017; originally announced March 2017.

Comments: 25 pages, 11 figures, 5 tables

The splitting of double-component active asteroid P/2016 J1 (PANSTARRS)

Authors: Fernando Moreno, Francisco Pozuelos, Bojan Novakovic, Javier Licandro, Antonio Cabrera-Lavers, Bryce Bolin, Robert Jedicke, Brett Gladman, Michele Bannister, Stephen Gwyn, Peter Veres, Kenneth Chambers, Serge Chastel, Larry Denneau, Heather Flewelling, Mark Huber, Eva Schunova-Lilly, Eugene Magnier, Richard Wainscoat, Christopher Waters, Robert Weryk, Davide Farnocchia, Marco Micheli

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

Submitted 13 February, 2017; originally announced February 2017.

Comments: Accepted by ApJ Letters, Feb. 17, 2017

The atmospheric impact trajectory of asteroid 2014 AA

Authors: D. Farnocchia, S. R. Chesley, P. G. Brown, P. W. Chodas

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

Submitted 13 January, 2016; originally announced January 2016.

High precision comet trajectory estimates: the Mars flyby of C/2013 A1 (Siding Spring)

Authors: D. Farnocchia, S. R. Chesley, M. Micheli, A. Delamere, R. S. Heyd, D. J. Tholen, J. D. Giorgini, W. M. Owen, L. K. Tamppari

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

Submitted 7 July, 2015; originally announced July 2015.

Systematic ranging and late warning asteroid impacts

Authors: D. Farnocchia, S. R. Chesley, M. Micheli

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

Submitted 31 March, 2015; originally announced April 2015.

A Study of Dust and Gas at Mars from Comet C/2013 A1 (Siding Spring)

Authors: Michael S. P. Kelley, Tony L. Farnham, Dennis Bodewits, Pasquale Tricarico, Davide Farnocchia

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

Submitted 13 August, 2014; v1 submitted 12 August, 2014; originally announced August 2014.

Comments: Accepted for publication in Astrophysical Journal Letters. 7 pages, 3 figures, 2 tables

Non-gravitational Perturbations and Virtual Impactors: the case of asteroid 2009 FD

Authors: F. Spoto, A. Milani, D. Farnocchia, S. R. Chesley, M. Micheli, G. B. Valsecchi, D. Perna, O. Hainaut

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

Submitted 22 October, 2014; v1 submitted 4 August, 2014; originally announced August 2014.

Star catalog position and proper motion corrections in asteroid astrometry

Authors: D. Farnocchia, S. R. Chesley, A. B. Chamberlin, D. J. Tholen

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

Submitted 31 July, 2014; originally announced July 2014.

Physical Properties of Near-Earth Asteroid 2011 MD

Authors: M. Mommert, D. Farnocchia, J. L. Hora, S. R. Chesley, D. E. Trilling, P. W. Chodas, M. Mueller, A. W. Harris, H. A. Smith, G. G. Fazio

Abstract: 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 μm) of the Infrared Array Camera and detected the target within the 2σ 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 M… ▽ More 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 μm) of the Infrared Array Camera and detected the target within the 2σ 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σ). 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. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: 20 pages, 4 figures

Journal ref: M. Mommert et al. 2014, ApJL, 789, L22

Delivery of dust grains from comet C/2013 A1 (Siding Spring) to Mars

Authors: Pasquale Tricarico, Nalin H. Samarasinha, Mark V. Sykes, Jian-Yang Li, Tony L. Farnham, Michael S. P. Kelley, Davide Farnocchia, Rachel Stevenson, James M. Bauer, Robert E. Lock

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

Submitted 28 April, 2014; originally announced April 2014.

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

Trajectory analysis for the nucleus and dust of comet C/2013~A1 (Siding Spring)

Authors: Davide Farnocchia, Steven R. Chesley, Paul W. Chodas, Pasquale Tricarico, Michael S. P. Kelley, Tony L. Farnham

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

Submitted 17 April, 2014; originally announced April 2014.

Constraining the Physical Properties of Near-Earth Object 2009 BD

Authors: M. Mommert, J. L. Hora, D. Farnocchia, S. R. Chesley, D. Vokrouhlický, D. E. Trilling, M. Mueller, A. W. Harris, H. A. Smith, G. G. Fazio

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

Submitted 30 March, 2014; originally announced March 2014.

Comments: 28 pages, 8 figures, accepted for publication in ApJ