Theory of the rotation of Janus and Epimetheus

The saturnian coorbital satellites Janus and Epimetheus present a unique dynamical configuration in the Solar System, because of high-amplitude horseshoe orbits, due to a mass ratio of order unity. As a consequence, they swap their orbits every 4 years, while their orbital periods is about 0.695 days. Recently, Tiscareno et al. (Tiscareno, M.S., Thomas, P.C., Burns, J.A. [2009]. Icarus 204, 254-261) got observational informations on the shapes and the rotational states of these satellites. In particular, they detected an offset in the expected equilibrium position of Janus, and a large libration of Epimetheus.We here propose to give a three-dimensional theory of the rotation of these satellites in using these observed data, and to compare it to the observed rotations. We consider the two satellites as triaxial rigid bodies, and we perform numerical integrations of the system in assuming the free librations as damped.The periods of the three free librations we get, associated with the three dimensions, are respectively 1.267, 2.179 and 2.098 days for Janus, and 0.747, 1.804 and 5.542 days for Epimetheus. The proximity of 0.747 days to the orbital period causes a high sensitivity of the librations of Epimetheus to the moments of inertia. Our theory explains the amplitude of the librations of Janus and the error bars of the librations of Epimetheus, but not an observed offset in the orientation of Janus.     

Tectonic patterns on reoriented and despun planetary bodies

Several processes may produce global tectonic patterns on the surface of a planetary body. The stresses associated with distortions of biaxial figures due to despinning or reorientation were first calculated by Vening Meinesz [Vening Meinesz, F.A., 1947. Trans. Am. Geophys. Union 28 (1), 1-23]. We adopt a mathematically equivalent, but physically more meaningful treatment for distortions associated with rotation. The new approach allows us to find analytic solutions for the general case of stresses associated with distortions of biaxial or triaxial planetary figures. Distortions of biaxial figures may be driven by variations in rotation rate, rotation axis orientation, or the combination of both. Distortions of triaxial figures may be driven by the same mechanisms and/or variations in tidal axis orientation for tidally deformed satellites. While the magnitude of the resulting stresses depends on the adopted elastic and physical parameters, the expected tectonic pattern is independent of these parameters for these mechanisms. Reorientation of the rotation/tidal axis alone is expected to produce normal/thrust faulting provinces enclosing the initial rotation/tidal poles, and thrust/normal faulting provinces enclosing the final rotation/tidal poles. Reorientation of both the rotation and tidal axis results in a wide variety of tectonic patterns for different reorientation geometries. On Europa, the tidal axis reorientation which generally accompanies rotation axis reorientations may provide an alternative explanation for tectonic features that have been interpreted as evidence for non-synchronous rotation. The observed tectonic pattern on Enceladus is more easily explained by a large reorientation (∼90°) of the rotation axis, than by rotation rate variations.     

Triaxiality of satellites and small bodies in the Solar System

The triaxial figures are very common shape of most of planetary satellites as well as of smaell bodies as asteroids. There are 21 satellites in the Solar System triaxial figures of which were detected in situ evidently (Davies et al., 1995). However, the total number of triaxiaxial satellites in the Solar System should be in fact larger. In this paper the general theory of triaxiality due to tidal forces is discussed in regard to the very recent numerical data. Since they orbit synchronously, as a rule: their orbital periods are equal to the rotational periods, the tidal forces may be responsible for their triaxial figures. On the other hand the origin of triaxiality of asteroids due to another process and the of their figures cannot be axplained by the tidal effects.     

The equilibrium of rubble-pile satellites: The Darwin and Roche ellipsoids for gravitationally held granular aggregates

Many new small moons of the giant planets have been discovered recently. In parallel, satellites of several asteroids, e.g., Ida, have been found. Strikingly, a majority of these new-found planetary moons are estimated to have very low densities, which, along with their hypothesized accretionary origins, suggests a rubble internal structure. This, coupled to the fact that many asteroids are also thought to be particle aggregates held together principally by self-gravity, motivates the present investigation into the possible ellipsoidal shapes that a rubble-pile satellite may achieve as it orbits an aspherical primary. Conversely, knowledge of the shape will constrain the granular aggregate's orbit—the closer it gets to a primary, both primary's tidal effect and the satellite's spin are greater. We will assume that the primary body is sufficiently massive so as not to be influenced by the satellite. However, we will incorporate the primary's possible ellipsoidal shape, e.g., flattening at its poles in the case of a planet, and the proloidal shape of asteroids. In this, the present investigation is an extension of the first classical Darwin problem to granular aggregates. General equations defining an ellipsoidal rubble pile's equilibrium about an ellipsoidal primary are developed. They are then utilized to scrutinize the possible granular nature of small inner moons of the giant planets. It is found that most satellites satisfy constraints necessary to exist as equilibrated granular aggregates. Objects like Naiad, Metis and Adrastea appear to violate these limits, but in doing so, provide clues to their internal density and/or structure. We also recover the Roche limit for a granular satellite of a spherical primary, and employ it to study the martian satellites, Phobos and Deimos, as well as to make contact with earlier work of Davidsson [Davidsson, B., 2001. Icarus 149, 375–383]. The satellite's interior will be modeled as a rigid-plastic, cohesion-less material with a Drucker–Prager yield criterion. This rheology is a reasonable first model for rubble piles. We will employ an approximate volume-averaging procedure that is based on the classical method of moments, and is an extension of the virial method [Chandrasekhar, S., 1969. Ellipsoidal Figures of Equilibrium. Yale Univ. Press, New Haven] to granular solid bodies.     

Analytical description of physical librations of saturnian coorbital satellites Janus and Epimetheus

Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every 4 years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (Tiscareno, M.S., Thomas, P.C., Burns, J.A. [2009]. Icarus 204, 254-261) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to Cassini ISS images. By a quasi-periodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long- and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results obtained for the short-period librations are studied with the help of an adiabatic-like approach.     

Orbital evolution of small binary asteroids

About 15% of both near-Earth and main-belt asteroids with diameters below 10 km are now known to be binary. These small asteroid binaries are relatively uniform and typically contain a fast-spinning, flattened primary and a synchronously rotating, elongated secondary that is 20-40% as large (in diameter) as the primary. The principal formation mechanism for these binaries is now thought to be YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect induced spin-up of the primary followed by mass loss and accretion of the secondary from the released material. It has previously been suggested (?uk, M. [2007]. Astrophys. J. 659, L57-L60) that the present population of small binary asteroids is in a steady state between production through YORP and destruction through binary YORP (BYORP), which should increase or decrease secondary’s orbit, depending on the satellite’s shape. However, BYORP-driven evolution has not been directly modeled until now. Here we construct a simple numerical model of the binary’s orbital as well the secondary’s rotational dynamics which includes BYORP and selected terms representing main solar perturbations. We find that many secondaries should be vulnerable to chaotic rotation even for relatively low-eccentricity mutual orbits. We also find that the precession of the mutual orbit for typical small binary asteroids might be dominated by the perturbations from the prolate and librating secondary, rather than the oblate primary. When we evolve the mutual orbit by BYORP we find that the indirect effects on the binary’s eccentricity (through the coupling between the orbit and the secondary’s spin) dominate over direct ones caused by the BYORP acceleration. In particular, outward evolution causes eccentricity to increase and eventually triggers chaotic rotation of the secondary. We conclude that the most likely outcome will be reestablishing of the synchronous lock with a “flipped” secondary which would then evolve back in. For inward evolution we find an initial decrease of eccentricity and secondary’s librations, to be followed by later increase. We think that it is likely that various forms of dissipation we did not model may damp the secondary’s librations close to the primary, allowing for further inward evolution and a possible merger. We conclude that a merger or a tidal disruption of the secondary are the most likely outcomes of the BYORP evolution. Dissociation into heliocentric pairs by BYORP alone should be very difficult, and satellite loss might be restricted to the minority of systems containing more than one satellite at the time.     

Impact cratering and break up of the small bodies of the Solar System

We consider the largest impact craters observed on small satellites and asteroids and the impact disruption of such bodies. Observational data are considered from 21 impact-like structures on 13 satellites and 8 asteroids (target body radii in the range 0.7-265 km). If the radius of the target body is R and the diameter of the largest crater observed on this body D, the ratio D/R is then the main observational parameter of interest. This is found on the observed bodies and compared to data obtained in the laboratory. Taking the largest observed value for D/R as a proxy for the ratio Dc/R (where Dc is the diameter of the largest crater that can be formed on a body without shattering it) it was found that for the observed icy satellites Dc,icy≈1.2R and for the asteroids and the rocky satellites Dc,rocky≈1.6R. In laboratory experiments with ice targets at impactor speeds of 1 to 3 km s⁻¹ we obtained Dc,icy≈1.64R.     

A steady-state model of NEA binaries formed by tidal disruption of gravitational aggregates

We present results of a simulation of a steady-state binary near-Earth asteroid (NEA) population. This study combines previous work on tidal disruption of gravitational aggregates [Walsh, K.J., Richardson, D.C., 2006. Icarus 180, 201-216] with a Monte Carlo simulation of NEA planetary encounters. Evolutionary effects include tidal evolution and binary disruption from close planetary encounters. The results show that with the best known progenitor (small Main Belt asteroids) shape and spin distributions, and current estimates of NEA lifetime and encounter probabilities, that tidal disruption should account for approximately 1-2% of NEAs being binaries. Given the best observed estimate of a ∼15% binary NEA fraction, we conclude that there are other formation mechanisms that contribute significantly to this population. We also present the expected distribution of binary orbital and physical properties for the steady-state binary NEAs formed by tidal disruption. We discuss the effects on binary fraction and properties due to changes in the least constrained parameters, and other possible effects on our model that could account for differences between the presented results and the observed binary population. Finally, we model possible effects of a significant population of binaries migrating to the near-Earth population from the Main Belt.     

Tidal disruptions: A continuum theory for solid bodies

Although the theory of Roche 1847 for the tidal disruption limits of orbiting satellites assumes a fluid body, a length to diameter of exactly 2.07:1, and a particular body orientation, the theory is commonly applied to the satellites of the Solar System and to small asteroids and comets passing nearby a planet. Clearly these bodies are neither fluid nor generally are that elongated, so a more appropriate theory is needed. Here we present exact analytical results for the distortion and disruption limits of solid spinning ellipsoidal bodies subjected to tidal forces, using the Drucker-Prager strength model with zero cohesion. It is the appropriate model for dry granular materials such as sands and rocks, for rubble-pile asteroids and comets, and for larger satellites, asteroids and comets where the cohesion can be ignored. This study uses the same approach as the studies of spin limits for solid ellipsoidal bodies given in [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303]. It is a static theory that predicts conditions for breakup and predicts the nature of the deformations at the limit state, but does not track the dynamics of the body as it comes apart. The strength is characterized by a single material parameter associated with an angle of friction, which can range from zero to 90°. The case with zero friction angle has no shear strength whatsoever, so it is then the model of a fluid or gas. The case of 90° represents a material that cannot fail in shear, but still has zero tensile strength. Typical dry soils have angles of friction of 30°-40°. Since the static fluid case is included in the theory as a special case, the classical results of Roche [Roche, E.A., 1847. Acad. Sci. Lett. Montpelier. Mem. Section Sci. 1, 243-262] and Jeans [Jeans, J.H., 1917. Mem. R. Astron. Soc. London 62, 1-48] are included and re-derived in their entirety; but the general solid case has much more variety and applicability. We consider both the spin-locked case, appropriate for most satellites of the Solar System; and the zero spin case, a possible case for a passing stray body. Detailed plots of many special cases are presented, in terms of shape, orientation and mass densities. A very typical result gives a closest approach d=1.5(ρ/ρP)^1/3R in terms of the planet radius R, and the satellite and planet mass densities ρ and ρP. We also use the theory to distinguish between conditions allowing global shape changes leading to new equilibrium states, or those leading to complete disruption. We apply the theory to the potentially hazardous Asteroid 99942 Apophis due to pass very near the Earth in 2029, and conclude it is extremely unlikely to experience any tidal readjustments during its passage. The states of many of the satellites of the Solar System are compared to the theory, and we find that all are well within their tidal disruption limits for expected values of the internal friction.     

Internal stresses in phobos and other triaxial bodies

The unusual dynamical behavior of Phobos, its strange appearance, and its mysterious network of grooves all make it an intriguing object. Geophysical studies, though, have been hampered by the lack of suitable theories applicable to nonspherical bodies. In this paper the Martian satellites are modeled as homogeneous, elastic triaxial ellipsoids subject to tidal, rotational, and self-gravitational stresses. A novel semianalytical treatment then gives the stress and strain fields throughout their interiors. Yield phenomena and their possible surface expressions are also investigated. The results indicate that Phobos and Deimos have always been stable with respect to tidal fracture or disruption, but that Phobos will probably break up before colliding with Mars. Applications of the new formulation to other nonspherical bodies in the solar system are also discussed.     

A survey for satellites of Venus

We present a systematic survey for satellites of Venus using the Baade-Magellan 6.5 m telescope and IMACS wide-field CCD imager at Las Campanas observatory in Chile. In the outer portions of the Hill sphere the search was sensitive to a limiting red magnitude of about 20.4, which corresponds to satellites with radii of a few hundred meters when assuming an albedo of 0.1. In the very inner portions of the Hill sphere scattered light from Venus limited the detection to satellites of about a kilometer or larger. Although several main belt asteroids were found, no satellites (moons) of Venus were detected.     

Tidal evolution of Dysnomia, satellite of the dwarf planet Eris

The past tidal evolution of the satellite Dysnomia of the dwarf planet Eris can be inferred from the current physical and orbital properties of the system. Preliminary considerations, which assumed a circular orbit for the satellite, suggested that the satellite formed close to the planet, perhaps as a result of a giant impact, and that it is thus unlikely that smaller satellites lie further out. However, if the satellite's orbit is eccentric, even if the eccentricity is very small, a qualitatively different past tidal evolution may be indicated. Early in the Solar System's history, the satellite may have been on a highly eccentric orbit much farther from the planet than it is now, suggestive of a capture origin. Additional satellites farther out cannot be ruled out.     

Plasma, plumes and rings: Saturn system dynamics as recorded in global color patterns on its midsize icy satellites

New global maps of the five inner midsize icy saturnian satellites, Mimas, Enceladus, Tethys, Dione, and Rhea, have been constructed in three colors (UV, Green and near-IR) at resolutions of 1 km/pixel. The maps reveal prominent global patterns common to several of these satellites but also three major color features unique to specific satellites or satellite subgroups. The most common features among the group are first-order global asymmetries in color properties. This pattern, expressed on Tethys, Dione and Rhea, takes the form of a ∼1.4-1.8 times enhancement in redness (expressed as IR/UV ratio) of the surface at the center of the trailing hemisphere of motion, and a similar though significantly weaker IR/UV enhancement at the center of the leading hemisphere. The peak in redness on the trailing hemisphere also corresponds to a known decrease in albedo. These double hemispheric asymmetries are attributable to plasma and E-ring grain bombardment on the trailing and leading hemispheres, respectively, for the outer three satellites Tethys, Dione and Rhea, whereas as E-ring bombardment may be focused on the trailing hemisphere of Mimas due to its orbital location interior to Enceladus. The maps also reveal three major deviations from these basic global patterns. We observe the previously known dark bluish leading hemisphere equatorial band on Tethys but have also discovered a similar band on Mimas. Similar in shape, both features match the surface patterns expected for irradiation of the surface by incident MeV electrons that drift in a direction opposite to the plasma flow. The global asymmetry on Enceladus is offset ∼40° to the west compared to the other satellites. We do not consider Enceladus in detail here, but the global distribution of bluish material can be shown to match the deposition pattern predicted for plume fallback onto the surface (Kempf, S., Beckmann, U., Schmidt, S. [2010]. Icarus 206, 446-457. doi:10.1016/j.icarus.2009.09.016). E-ring deposition on Enceladus thus appears to mask or prevent the formation of the lenses and hemispheric asymmetries we see on the other satellites. Finally, we observe a chain of discrete bluish splotches along the equator of Rhea. Unlike the equatorial bands of Tethys and Mimas, these splotches form a very narrow great circle ?10-km wide (north-to-south) and appear to be related to surface disruption, exposing fresh, bluish ice on older crater rims. This feature is unique to Rhea and may have formed by impact onto its surface of orbiting material.     

Rotational fission of contact binary asteroids

The energetics and dynamics of contact binary asteroids as they approach and pass the rotational fission limit is studied. We presume that the asteroids are subject to an external torque, such as from the YORP effect, that increases their angular momentum. Furthermore, we assume the asteroids can be described by a fairly minimal model comprised of a sphere and ellipsoid resting on each other. The minimum energy configurations for contact binary asteroids at different levels of angular momentum are computed and discussed. We find distinct transitions between different configurations as the angular momentum of the system is increased. These indicate that rapidly rotating contact binary asteroids may seek out clearly different relative configurations than slowly rotating systems. We find a single end state of the systems prior to rotational fission, and distinct dynamical outcomes as a function of mass distribution and shape when the rotational fission limit is exceeded. Our theoretical results agree qualitatively with observed properties of near-Earth asteroids, and can be used to help explain the spin-rate barrier, contact binaries, and the observed morphology of most NEO binaries.     

Modeling global impact effects on middle-sized icy bodies: applications to Saturn's moons

Disrupted terrains that form as a consequence of giant impacts may help constrain the internal structures of planets, asteroids, comets and satellites. As shock waves and powerful seismic stress waves propagate through a body, they interact with the internal structure in ways that may leave a characteristic impression upon the surface. Variations in peak surface velocity and tensile stress, related to landform degradation and surface rupture, may be controlled by variations in core size, shape and density. Caloris Basin on Mercury and Imbrium Basin on the Moon have disturbed terrain at their antipodes, where focusing is most intense for an approximately symmetric spheroid. Although, the icy saturnian satellites Tethys, Mimas, and Rhea possess giant impact structures, it is not clear whether these structures have correlated disrupted terrains, antipodal or elsewhere. In anticipation of high-resolution imagery from Cassini, we investigate antipodal focusing during giant impacts using a 3D SPH impact model. We first investigate giant impacts into a fiducial 1000 km diameter icy satellite with a variety of core radii and compositions. We find that antipodal disruption depends more on core radius than on core density, suggesting that core geometry may express a surface signature in global impacts on partially differentiated targets. We model Tethys, Mimas, and Rhea according to their image-derived shapes (triaxial for Tethys and Mimas and spherical for Rhea), varying core radii and densities to give the proper bulk densities. Tethys shows greater antipodal values of peak surface velocity and peak surface tensile stress, indicating more surface damage, than either Mimas or Rhea. Results for antipodal and global fragmentation and terrain rupture are inconclusive, with the hydrocode itself producing global disruption at the limits of model resolution but with peak fracture stresses never exceeding the strength of laboratory ice.     

The rotation of Janus and Epimetheus

Epimetheus, a small moon of Saturn, has a rotational libration (an oscillation about synchronous rotation) of 5.9°±1.2°, placing Epimetheus in the company of Earth’s Moon and Mars’ Phobos as the only natural satellites for which forced rotational libration has been detected. The forced libration is caused by the satellite’s slightly eccentric orbit and non-spherical shape.Detection of a moon’s forced libration allows us to probe its interior by comparing the measured amplitude to that predicted by a shape model assuming constant density. A discrepancy between the two would indicate internal density asymmetries. For Epimetheus, the uncertainties in the shape model are large enough to account for the measured libration amplitude. For Janus, on the other hand, although we cannot rule out synchronous rotation, a permanent offset of several degrees between Janus’ minimum moment of inertia (long axis) and the equilibrium sub-Saturn point may indicate that Janus does have modest internal density asymmetries.The rotation states of Janus and Epimetheus experience a perturbation every 4 years, as the two moons “swap” orbits. The sudden change in the orbital periods produces a free libration about synchronous rotation that is subsequently damped by internal friction. We calculate that this free libration is small in amplitude (<0.1°) and decays quickly (a few weeks, at most), and is thus below the current limits for detection using Cassini images.     

Aqueous alteration affecting the irregular outer planets satellites: Evidence from spectral reflectance

The surface reflectance properties of the irregular outer planets satellites are probed for evidence for the presence of aqueous alteration products on their surfaces using the strong correlation between the 3.0-μm water of hydration absorption feature and the 0.7-μm Fe²⁺ → Fe³⁺ oxidized iron feature seen in low-albedo asteroid reflectances, in an effort to expand our understanding of the composition of the precursor bodies from which the dynamical satellite clusters are derived. Equations converting Johnson V and Kron-Cousins RI photometry to Eight Color Asteroid Survey v (0.550 μm), w (0.701 μm), and x (0.853 μm) photometry are derived from relationships defined by Howell (1995, Ph.D. thesis), and coupled with an algorithm previously defined to detect the presence of the 0.7-μm absorption feature in ECAS asteroid photometry [Vilas, F., 1994. Icarus 111, 456-467]. Broadband VRI photometry of Ch-class Asteroid 19 Fortuna acquired during 2004 confirms the efficacy of this method of identifying the presence of the 0.7-μm feature. Photometric observations of many recently discovered irregular outer jovian, saturnian, uranian, and neptunian satellites, coupled with limited asteroid spectroscopy, were examined for the presence of aqueous alteration. The dynamical clusters of outer irregular jovian satellites are mixed between objects that do and do not show this absorption feature. Multiple observations of some objects test both positively and negatively, similar to the surface variegation that has been observed among many C-class asteroids in the main asteroid belt. Evidence for aqueous alteration on these jovian satellites augers for an origin in or near the same location as the asteroids now occupying the aqueous alteration zone (2.6-3.5 AU), at heliocentric distances internal to Jupiter's orbit. Among the saturnian irregular satellites, only S IX Phoebe shows limited evidence of aqueous alteration from ground-based observations. The other satellites show no sign of this feature, and have general reflectance properties very similar to the D-class asteroids, supporting an origin for their precursor bodies in the outer Solar System, perhaps the Centaur region. Only two uranian satellites were tested: U XVII Caliban tests positively for the feature. The differences in surface reflectance properties support the idea that Caliban and U XVI Sycorax derive from separate parent bodies. One observation of neptunian satellite N II Nereid shows no sign of this absorption feature.     

Three-body capture of irregular satellites: Application to Jupiter

We investigate a new theory of the origin of the irregular satellites of the giant planets: capture of one member of a ∼100-km binary asteroid after tidal disruption. The energy loss from disruption is sufficient for capture, but it cannot deliver the bodies directly to the observed orbits of the irregular satellites. Instead, the long-lived capture orbits subsequently evolve inward due to interactions with a tenuous circumplanetary gas disk.We focus on the capture by Jupiter, which, due to its large mass, provides a stringent test of our model. We investigate the possible fates of disrupted bodies, the differences between prograde and retrograde captures, and the effects of Callisto on captured objects. We make an impulse approximation and discuss how it allows us to generalize capture results from equal-mass binaries to binaries with arbitrary mass ratios.We find that at Jupiter, binaries offer an increase of a factor of ∼10 in the capture rate of 100-km objects as compared to single bodies, for objects separated by tens of radii that approach the planet on relatively low-energy trajectories. These bodies are at risk of collision with Callisto, but may be preserved by gas drag if their pericenters are raised quickly enough. We conclude that our mechanism is as capable of producing large irregular satellites as previous suggestions, and it avoids several problems faced by alternative models.     

Constraints from deuterium on the formation of icy bodies in the Jovian system and beyond

We consider the role of deuterium as a potential marker of location and ambient conditions during the formation of small bodies in our Solar system. We concentrate in particular on the formation of the regular icy satellites of Jupiter and the other giant planets, but include a discussion of the implications for the Trojan asteroids and the irregular satellites. We examine in detail the formation of regular planetary satellites within the paradigm of a circum-Jovian subnebula. Particular attention is paid to the two extreme potential subnebulae—“hot” and “cold”. In particular, we show that, for the case of the “hot” subnebula model, the D:H ratio in water ice measured from the regular satellites would be expected to be near-Solar. In contrast, satellites which formed in a “cold” subnebula would be expected to display a D:H ratio that is distinctly over-Solar. We then compare the results obtained with the enrichment regimes which could be expected for other families of icy small bodies in the outer Solar system—the Trojan asteroids and the irregular satellites. In doing so, we demonstrate how measurements by Laplace, the James Webb Space Telescope, HERSCHEL and ALMA will play an important role in determining the true formation locations and mechanisms of these objects.     

Orbital resonances in the inner neptunian system: II. Resonant history of Proteus, Larissa, Galatea, and Despina

We investigate the orbital history of the small neptunian satellites discovered by Voyager 2. Over the age of the Solar System, tidal forces have caused the satellites to migrate radially, bringing them through mean-motion resonances with one another. In this paper, we extend our study of the largest satellites Proteus and Larissa [Zhang, K., Hamilton, D.P., 2007. Icarus 188, 386-399] by adding in mid-sized Galatea and Despina. We test the hypothesis that these moons all formed with zero inclinations, and that orbital resonances excited their tilts during tidal migration. We find that the current orbital inclinations of Proteus, Galatea, and Despina are consistent with resonant excitation if they have a common density . Larissa's inclination, however, is too large to have been caused by resonant kicks between these four satellites; we suggest that a prior resonant capture event involving either Naiad or Thalassa is responsible. Our solution requires at least three past resonances with Proteus, which helps constrain the tidal migration timescale and thus Neptune's tidal quality factor: 9000<QN<36,000. We also improve our determination of Qs for Proteus and Larissa, finding 36<QP<700 and 18<QL<200. Finally, we derive a more general resonant capture condition, and work out a resonant overlap criterion relevant to satellite orbital evolution around an oblate primary.