The conventional least-squares asteroid mass determination algorithm allows us to solve for the mass of a large subject asteroid that is perturbing the trajectory of a smaller test asteroid. However, ...this algorithm is necessarily a first approximation, ignoring the possibility that the subject asteroid may itself be perturbed by the test asteroid, or that the encounter's precise geometry may be entangled with encounters involving other asteroids. After reviewing the conventional algorithm, we use it to calculate the masses of 30 main-belt asteroids. Compared to our previous results, we find new mass estimates for eight asteroids (11 Parthenope, 27 Euterpe, 51 Neimausa, 76 Freia, 121 Hermione, 324 Bamberga, 476 Hedwig, and 532 Herculina) and significantly more precise estimates for six others (2 Pallas, 3 Juno, 4 Vesta, 9 Metis, 16 Psyche, and 88 Thisbe). However, we also find that the conventional algorithm yields questionable results in several gravitationally coupled cases. To address such cases, we describe a new algorithm that allows the epoch state vectors of the subject asteroids to be included as solve-for parameters, allowing for the simultaneous solution of the masses and epoch state vectors of multiple subject and test asteroids. We then apply this algorithm to the same 30 main-belt asteroids and conclude that mass determinations resulting from current and future high-precision astrometric sources (such as Gaia) should conduct a thorough search for possible gravitational couplings and account for their effects.
Abstract We characterize asteroid (16) Psyche using high-precision astrometry, including the recent Gaia Focused Product Release. The gravitational perturbations of Psyche on other asteroids can be ...observable in the case of mutual encounters. Using a least squares approach, we estimate the mass of Psyche by fitting astrometric data of asteroids that come within 0.05 au of Psyche. Combining the resulting individual estimates, we find GM = 1.601 ± 0.017 km 3 s −2 . This result is robust against variations in the orbit determination setup and specific data set used. The volume and equivalent radius of Psyche are currently constrained by occultations and radar and optical imaging to (5.75 ± 0.19) × 10 6 km 3 and 111 − 0.5 + 2 km, respectively. Given the volume of Psyche, our mass estimate corresponds to a bulk density of 4172 ± 145 kg m −3 , which is compatible with an M-type taxonomic classification. Finally, the phase-dependent photocenter offset is visible in the residuals of Gaia astrometric observations of Psyche. This effect is consistent with the size of Psyche.
The solid, central part of a comet--its nucleus--is subject to destructive processes, which cause nuclei to split at a rate of about 0.01 per year per comet. These destructive events are due to a ...range of possible thermophysical effects; however, the geophysical expressions of these effects are unknown. Separately, over two-thirds of comet nuclei that have been imaged at high resolution show bilobate shapes, including the nucleus of comet 67P/Churyumov-Gerasimenko (67P), visited by the Rosetta spacecraft. Analysis of the Rosetta observations suggests that 67P's components were brought together at low speed after their separate formation. Here, we study the structure and dynamics of 67P's nucleus. We find that sublimation torques have caused the nucleus to spin up in the past to form the large cracks observed on its neck. However, the chaotic evolution of its spin state has so far forestalled its splitting, although it should eventually reach a rapid enough spin rate to do so. Once this occurs, the separated components will be unable to escape each other; they will orbit each other for a time, ultimately undergoing a low-speed merger that will result in a new bilobate configuration. The components of four other imaged bilobate nuclei have volume ratios that are consistent with a similar reconfiguration cycle, pointing to such cycles as a fundamental process in the evolution of short-period comet nuclei. It has been shown that comets were not strong contributors to the so-called late heavy bombardment about 4 billion years ago. The reconfiguration process suggested here would preferentially decimate comet nuclei during migration to the inner solar system, perhaps explaining this lack of a substantial cometary flux.
Abstract
Near-Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and ...comets, which starts with their removal from the respective source regions—the main belt and trans-Neptunian scattered disk—and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the solar system. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main-belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the
ν
6
and 3:1 resonances producing ≃30% of NEOs with absolute magnitudes
H
= 15 and ≃80% of NEOs with
H
= 25. Hence, the large and small NEOs have different orbital distributions. The inferred flux of
H
< 18 bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate (≃2 × 10
−4
au Myr
−1
for diameter
D
= 1 km and semimajor axis
a
= 2.5 au). This implies obliquities
θ
≃ 0° for
a
< 2.5 au and
θ
≃ 180° for
a
> 2.5 au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model.
As an application of our recent observational error model, we present the astrometric masses of 26 main-belt asteroids. We also present an integrated ephemeris of 300 large asteroids, which was used ...in the mass determination algorithm to model significant perturbations from the rest of the main belt. After combining our mass estimates with those of other authors, we study the bulk porosities of over 50 main-belt asteroids and observe that asteroids as large as 300 km in diameter may be loose aggregates. This finding may place specific constraints on models of main-belt collisional evolution. Additionally, we observe that C-group asteroids tend to have significantly higher macroporosity than S-group asteroids.
Abstract
The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multiband wide-field synoptic survey will transform our view of the ...solar system, with the discovery and monitoring of over five million small bodies. The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor planets and passing interstellar objects. Creating an inventory of the solar system is one of the four main LSST science drivers. The LSST observing cadence is a complex optimization problem that must balance the priorities and needs of all the key LSST science areas. To design the best LSST survey strategy, a series of operation simulations using the Rubin Observatory scheduler have been generated to explore the various options for tuning observing parameters and prioritizations. We explore the impact of the various simulated LSST observing strategies on studying the solar system’s small body reservoirs. We examine what are the best observing scenarios and review what are the important considerations for maximizing LSST solar system science. In general, most of the LSST cadence simulations produce ±5% or less variations in our chosen key metrics, but a subset of the simulations significantly hinder science returns with much larger losses in the discovery and light-curve metrics.
In this paper, we discuss the detection of systematic biases in star positions of the USNO A1.0, A2.0, and B1.0 catalogs, as deduced from the residuals of numbered asteroid observations. We present a ...technique for the removal of these biases, and validate this technique by illustrating the resulting improvements in numbered asteroid residuals, and by establishing that debiased orbits predict omitted observations more accurately than do orbits derived from non-debiased observations. We also illustrate the benefits of debiasing to high-precision astrometric applications such as asteroid mass determination and collision analysis, including a refined prediction of the impact probability of 99942 Apophis. Specifically, we find the IP of Apophis to be lowered by nearly an order of magnitude to 4.5
×
10
−6 for the 2036 close approach.
•Asteroid Bennu has a well-determined orbit due primarily to 12years of radar ranging.•The Yarkovsky effect on Bennu causes a semimajor axis drift of 284±1.5m/year.•We estimate Bennu’s bulk density ...at 1260±70kg/m3 and macroporosity 40±10%.•Bennu has a 1 in 2700 chance of an Earth impact late in the 22nd century.
The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ36), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999 to 2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate da/dt=(-19.0±0.1)×10-4au/Myr or 284±1.5m/year. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with 3σ significance in da/dt.
The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of ρ=1260±70kg/m3, which indicates a macroporosity in the range 40±10% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is (7.8±0.9)×1010kg and GM=5.2±0.6m3/s2.
Bennu’s Earth close approaches are deterministic over the interval 1654–2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and numerous potential impacts in subsequent years, from 2175 to 2196. The highest individual impact probability is 9.5×10-5 in 2196, and the cumulative impact probability is 3.7×10-4, leading to a cumulative Palermo Scale of −1.70.
Asteroid (410777) 2009 FD could hit 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 seven-dimensional space that includes orbital elements and the parameter characterizing the Yarkovsky effect. The highest impact probability is 2.7 × 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 × 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.
This paper analyzes the current population of known near-Earth asteroids (NEAs) and presents statistics on the recoverability of NEAs with both targeted observation campaigns and all-sky surveys. For ...an asteroid to be observable at a future apparition, given the right geometry, the plane-of-sky uncertainty must be small enough to be covered by a telescope's field of view and the asteroid must be brighter than the detector's limiting magnitude. Since recoverability is a telescope-dependent property, we select two representative instruments that span a wide range of capability and availability: the 1.0 m I52 telescope of the Catalina Sky Survey and the Hyper Suprime-Cam of the 8.2 m Subaru telescope. Based on this choice, we classify asteroids as recoverable, potentially recoverable, and not recoverable depending on whether they could be detected with an I52-class telescope, only with a Subaru-class telescope, or with neither, respectively. Using these definitions, we find that the majority (90%) of NEAs with H < 22 and most (93%) potentially hazardous asteroids are recoverable or potentially recoverable in the next 50 yr. When considering fainter asteroids down to H ≤ 28, about two-thirds of the NEA population and half of the low minimum-orbit intersection distance (MOID) asteroids (MOID ≤ 0.05 au) are either recoverable or potentially recoverable. As of 2019 October 13, the Sentry risk list includes 193 objects with an impact probability greater than 10−6 that are not recoverable. The fraction of NEAs and low-MOID NEAs that are not recoverable can be reduced by up to 47% and 43%, respectively, when incorporating statistical estimates of serendipitous recoveries by all-sky surveys.