Orbital evolution of the Gefion and Adeona asteroid families: close encounters with massive asteroids and the Yarkovsky effect

Asteroid families are groupings of minor planets identified by clustering in their proper orbital elements; these objects have spectral signatures consistent with an origin in the break-up of a common parent body. From the current values of proper semimajor axes a of family members one might hope to estimate the ejection velocities with which the fragments left the putative break-up event (assuming that the pieces were ejected isotropically). However, the ejection velocities so inferred are consistently higher than N-body and hydro-code simulations, as well as laboratory experiments, suggest. To explain this discrepancy between today’s orbital distribution of asteroid family members and their supposed launch velocities, we study whether asteroid family members might have been ejected from the collision at low speeds and then slowly drifted to their current positions, via one or more dynamical processes. Studies show that the proper a of asteroid family members can be altered by two mechanisms: (i) close encounters with massive asteroids, and (ii) the Yarkovsky non-gravitational effect. Because the Yarkovsky effect for kilometer-sized bodies decreases with asteroid diameter D, it is unlikely to have appreciably moved large asteroids (say those with D > 15 km) over the typical family age (1-2 Gyr).For this reason, we numerically studied the mobility of family members produced by close encounters with main-belt, non-family asteroids that were thought massive enough to significantly change their orbits over long timescales. Our goal was to learn the degree to which perturbations might modify the proper a values of all family members, including those too large to be influenced by the Yarkovsky effect. Our initial simulations demonstrated immediately that very few asteroids were massive enough to significantly alter relative orbits among family members. Thus, to maximize gravitational perturbations in our 500-Myr integrations, we investigated the effect of close encounters on two families, Gefion and Adeona, that have high encounter probabilities with 1 Ceres, by far the largest asteroid in the main belt. Our results show that members of these families spreads in a of less than 5% since their formation. Thus gravitational interactions cannot account for the large inferred escape velocities.The effect of close encounters with massive asteroids is, however, not entirely negligible. For about 10% of the simulated bodies, close encounters increased the “inferred” ejection velocities from sub-100 m/s to values greater than 100 m/s, beyond what hydro-code and N-body simulations suggest are the maximum possible initial ejection velocity for members of Adeona and Gefion with D > 15 km. Thus this mechanism of mobility may be responsible for the unusually high inferred ejection speeds of a few of the largest members of these two families.To understand the orbital evolution of the entire family, including smaller members, we also performed simulations to account for the drift of smaller asteroids caused by the Yarkovsky effect. Our two sets of simulations suggest that the two families we investigated are relatively young compared to larger families like Koronis and Themis, which have estimated ages of about 2 Byr. The Adeona and Gefion families seems to be no more than 600 and 850 Myr old, respectively.     

Prometheus and Pandora: masses and orbital positions during the Cassini tour

Hubble Space Telescope (HST) images of Prometheus and Pandora show longitude discrepancies of about 20° with respect to the Voyager ephemerides, with an abrupt change in mean motion at the end of 2000 (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07). These discrepancies are anti-correlated and arise from chaotic interactions between the two moons, occurring at interval of 6.2 yr, when their apses are anti-aligned (Goldreich and Rappaport, 2003a, Icarus 162, 391-399). This behavior is attributed to the overlap of four 121:118 apse-type mean motion resonances (Goldreich and Rappaport, 2003b, Icarus 166, 320-327). We study the Prometheus-Pandora system using a Radau-type integrator taking into account Saturn's oblateness up to and including terms in J₆, plus the effects of the major satellites. We first confirm the chaotic behavior of Prometheus and Pandora. By fitting the numerical integrations to the HST data (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07), we derive the satellite masses. The resulting GM values (with their standard 3-σ errors) for Prometheus and Pandora are respectively and . Using the nominal shape of the two moons (Thomas, 1989, Icarus 77, 248-274), we derive Prometheus and Pandora's densities, 0.40^+0.03_−0.07 and 0.49^+0.05_−0.09 g cm⁻³, respectively. Our numerical fits also enable us to constrain the time of the latest apse anti-alignment in 2000. Finally, using our fit, we predict the orbital positions of the two satellites during the Cassini tour, and provide a lower limit of the uncertainties due to chaos. These uncertainties amount to about 0.2° in mean longitude at the arrival of the Cassini spacecraft in July 2004, and to about 3° in 2008, at the end of the nominal tour.     

Orbits and masses of Saturn's coorbital and F-ring shepherding satellites

We have obtained numerically integrated orbits for Saturn's coorbital satellites, Janus and Epimetheus, together with Saturn's F-ring shepherding satellites, Prometheus and Pandora. The orbits are fit to astrometric observations acquired with the Hubble Space Telescope and from Earth-based observatories and to imaging data acquired from the Voyager spacecraft. The observations cover the 38 year period from the 1966 Saturn ring plane crossing to the spring of 2004. In the process of determining the orbits we have found masses for all four satellites. The densities derived from the masses for Janus, Epimetheus, Prometheus, and Pandora in units of g cm⁻³ are , , , and , respectively.     

The orbit of the Kappa Cygnid and related meteor streams

A study of the Kappa Cygnid and other minor streams of the August epoch is presented based on a computer search in a sample of 3518 photographic meteoroid orbits. Four different meteoroid streams with radiants in Cygnus, Draco and Lyra, were found. Three of these: the Kappa Cygnids, the Alpha Lyrids and the Zeta Draconids are identified with meteor showers reported by nineteenth century visual observers. The fourth stream, the August Lyrids, consists of six meteors with radiants in Lyra centered on = 277°.6, = 46°.2. No previous visual reports of this stream have been found. It is interesting to note that all four meteoroid streams are coincident in time; their orbits are all of short period and they all have very nearly the same orientation of semi-major axis.     

Effects of thermal radiation on the dynamics of binary NEAs

The Yarkovsky force, produced when thermal radiation is re-emitted asymmetrically, causes significant orbital evolution of asteroids in the 10 m-10 km size range. When acting on a non-spherical body, the momentum carried by this radiation generally produces a torque, called the YORP effect, which may be important in re-orienting asteroidal spins. Here we explore a related effect, the “binary YORP” (BYORP), that can modify the orbit of a synchronously rotating secondary in a binary system. It arises because a locked secondary is effectively an asymmetric appendage of the primary. It should be particularly important for Near-Earth Asteroids (NEAs) owing to their small sizes, proximity to the Sun, and benign collisional environment. To estimate BYORP's strength, we subjected 100 random Gaussian spheroids to the thermal radiation model of Rubincam (2000, Radiative spin-up and spin-down of small asteroids, Icarus, 148, 2-11). For most shapes, a significant torque arose on the secondary's orbit, typically modifying the orbit's size, eccentricity and inclination in less than 10⁵ years, for components of 1 and 0.3 km radii, separated by 2 km, at 1 AU, each of density 1750 kg m⁻³. Together YORP and BYORP are capable of synchronizing secondaries and circularizing orbits, making tidal dissipation unnecessary to explain the evolved state of observed NEA pairs. However, BYORP's rapid timescale poses a problem for the abundance of observed NEA binaries, since their formation rate is thought to be much slower. We consider and reject the following resolutions of this quandary: (i) the approximation using Gaussian spheroids inadequately models YORP; (ii) most secondaries are not synchronous, but inhabit other spin-orbit resonances (very unlikely); (iii) tidal dissipation is much more efficient than previously estimated, and thus capable of stabilizing observed systems; and (iv) moderately close encounters with planets can re-orient secondaries, turning BYORP into a slower, random-walk process. Finally, we speculate that most observed binary NEAs are in a stable state in which the obliquity-changing torques of YORP (acting on the primary) and BYORP cancel out, and that those systems must be close to 55° inclination, where the momentum-changing torques of both YORP and BYORP tend to be very small. Some retrograde systems might develop such that the nodes precess at a Sun-synchronous rate, while some prograde ones might move into the “evection” resonance. All three of these hypotheses can be tested directly by comparison with the i, Ω and ? observed for NEA binaries.     

Asteroid identification at discovery

We present a novel method for the search of linkages among astrometric observations of asteroids, that is, tentative identifications among asteroids observed. Having two different master sets of asteroid observations each containing a number of separate subsets, we define a linkage as a pair of subsets residing in separate master sets that can be tied together with an orbit for given observational errors. To find linkages among a wealth of observations we use an efficient stepwise filtering approach. First, we start with what we call phase-space address comparison. The first step substantially reduces the initially huge amount of pairs by requiring that pairs to be subjected to further analysis have similar geocentric spherical coordinates at common epochs (for example, at three epochs). Second, we search for orbits for each of the selected pairs of subsets. Succeeding in the effort proves that a linkage exists. If there are contradictions among linkages found—for example, a single subset being linked to several mutually exclusive subsets—additional new or archive observations are usually needed to discard erroneous linkages. The new method is built on six-dimensional statistical orbital inversion (Ranging), and is therefore particularly suitable for analyzing objects with the shortest observational arcs, that is, newly discovered asteroids (and comets). Results from extensive and successful tests on simulated survey observations are presented and discussed. Theoretical and empirical scaling results show that the method is applicable to future large-scale surveys that will increase the rate of asteroid discovery by at least two orders of magnitude. The successful linking of faint single-night observation sets obtained with the Very Large Telescope are briefly reviewed.     

Critical Periodic Orbits in the Restricted Three–Body Problem With Oblateness

In a recent paper in this journal series of horizontally critical symmetric periodic orbits of the six basic families of the photogravitational restricted three–body problem were computed (Perdios, 2003). In this paper, such series are determined in the framework of the restricted three–body problem when the more massive primary is an oblate spheroid. The vertical stability of the horizontally critical orbits is also computed.     

Cassini ISS astrometric observations of the inner jovian satellites, Amalthea and Thebe

We present a total of 289 new astrometric observations of the inner jovian satellites, Amalthea and Thebe, obtained using the Cassini ISS narrow angle camera. Observations were made using image sequences from 2000 December 11-12 (inbound) and 2001 January 15-16 (outbound), at phase angles of approximately 2° and 122°, respectively. Target distances were of order 284 R_J, giving a maximum resolution of approximately 100 km/pixel. Centroided line and sample values for 239 observations of Amalthea and 50 of Thebe are provided, together with estimated camera pointing information for each image. Orbit fitting using a uniformly precessing Keplerian ellipse model, taking into account the oblateness of Jupiter up to terms in J₆, gave RMS fit residuals of 0.364 and 0.443 pixel for Amalthea and Thebe, respectively (equivalent to 0.450 and 0.547 arcsec). RMS residuals relative to the JPL JUP230 ephemeris were 0.306 and 0.604 pixel (equivalent to 0.378 and 0.746 arcsec), for Amalthea and Thebe. The fitted orbital parameters confirm the relatively high inclinations of these satellites (0.374°±0.002° and 1.076°±0.003°, respectively), equivalent to maximum vertical displacements above Jupiter's equatorial plane of 1188±6 and 4240±12 km, respectively, consistent with current estimates of the half-thicknesses of the Amalthea and Thebe gossamer rings [Ockert-Bell, M.E., Burns, J.A., Dauber, I.J., Thomas, P.C., Veverka, J., Belton, M.J.S., Klaasen, K.P., 1999. Icarus 138, 188-213].     

Orbit determination with very short arcs: II. Identifications

When the observational data are not enough to compute a meaningful orbit for an asteroid/comet we can represent the data with an attributable, i.e., two angles and their time derivatives. The undetermined variables range and range rate span an admissible region of Solar System orbits, which can be sampled by a set of Virtual Asteroids (VAs) selected by means of an optimal triangulation [Milani, A., Gronchi, G.F., de' Michieli Vitturi, M., Kne?evi?, Z., 2004. Celest. Mech. Dyn. Astron. 90, 59-87]. The attributable 4 coordinates are the result of a fit and they have an uncertainty, represented by a covariance matrix. Two short arcs of observations, represented by two attributables, can be linked by considering for each VA (in the admissible region of the first arc) the covariance matrix for the prediction at the time of the second arc, and by comparing it with the attributable of the second arc with its own covariance. By defining an identification penalty we can select the VAs allowing to fit together both arcs and compute a preliminary orbit. Two attributables may not be enough to compute an orbit with convergent differential corrections. Thus the preliminary orbit is used in a constrained differential correction, providing solutions along the Line Of Variation which can be used as second generation VAs to further predict the observations at the time of a third arc. In general the identification with a third arc will ensure a well determined orbit, to which additional sets of observations can be attributed. To test these algorithms we use a large scale simulation and measure the completeness, the reliability and the efficiency of the overall procedure to build up orbits by accumulating identifications. Under the conditions expected for the next generation asteroid surveys, the methods developed in this and in the preceding papers are efficient enough to be used as primary identification methods, with very good results. One important property is that the completeness in finding the possible identifications is as good for comparatively rare orbits, such as the ones of Near-Earth Objects, as for main belt orbits.     

Evolution of a Keplerian disk of colliding and fragmenting particles: a kinetic model with application to the Edgeworth-Kuiper belt

We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune.