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.     

Understanding the behavior of Prometheus and Pandora

We revisit the dynamics of Prometheus and Pandora, two small moons flanking Saturn's F ring. Departures of their orbits from freely precessing ellipses result from mutual interactions via their 121:118 mean motion resonance. Motions are chaotic because the resonance is split into four overlapping components. Orbital longitudes were observed to drift away from predictions based on Voyager ephemerides. A sudden jump in mean motions took place close to the time at which the orbits' apses were antialigned in 2000. Numerical integrations reproduce both the longitude drifts and the jumps. The latter have been attributed to the greater strength of interactions near apse antialignment (every 6.2 yr), and it has been assumed that this drift-jump behavior will continue indefinitely. We re-examine the dynamics of the Prometheus-Pandora system by analogy with that of a nearly adiabatic, parametric pendulum. In terms of this analogy, the current value of the action of the satellite system is close to its maximum in the chaotic zone. Consequently, at present, the two separatrix crossings per precessional cycle occur close to apse antialignment. In this state libration only occurs when the potential's amplitude is nearly maximal, and the “jumps” in mean motion arise during the short intervals of libration that separate long stretches of circulation. Because chaotic systems explore the entire region of phase space available to them, we expect that at other times the Prometheus-Pandora system would be found in states of medium or low action. In a low action state it would spend most of the time in libration, and separatrix crossings would occur near apse alignment. We predict that transitions between these different states can happen in as little as a decade. Therefore, it is incorrect to assume that sudden changes in the orbits only happen near apse antialignment.     

Finding the trigger to Iapetus’ odd global albedo pattern: Dynamics of dust from Saturn’s irregular satellites

The leading face of Saturn’s moon Iapetus, Cassini Regio, has an albedo only one tenth that on its trailing side. The origin of this enigmatic dichotomy has been debated for over 40 years, but with new data, a clearer picture is emerging. Motivated by Cassini radar and imaging observations, we investigate Soter’s model of dark exogenous dust striking an originally brighter Iapetus by modeling the dynamics of the dark dust from the ring of the exterior retrograde satellite Phoebe under the relevant perturbations. In particular, we study the particles’ probabilities of striking Iapetus, as well as their expected spatial distribution on the Iapetian surface. We find that, of the long-lived particles (?5 μm), most particle sizes (?10 μm) are virtually certain to strike Iapetus, and their calculated distribution on the surface matches up well with Cassini Regio’s extent in its longitudinal span. The satellite’s polar regions are observed to be bright, presumably because ice is deposited there. Thus, in the latitudinal direction we estimate polar dust deposition rates to help constrain models of thermal migration invoked to explain the bright poles (Spencer, J.R., Denk, T. [2010]. Science 327, 432-435). We also analyze dust originating from other irregular outer moons, determining that a significant fraction of that material will eventually coat Iapetus—perhaps explaining why the spectrum of Iapetus’ dark material differs somewhat from that of Phoebe. Finally we track the dust particles that do not strike Iapetus, and find that most land on Titan, with a smaller fraction hitting Hyperion. As has been previously conjectured, such exogenous dust, coupled with Hyperion’s chaotic rotation, could produce Hyperion’s roughly isotropic, moderate-albedo surface.     

Appearance of Saturn’s F ring azimuthal channels for the anti-alignment configuration between the ring and Prometheus

In this article we explore the aspect of the F ring with respect to the anti-alignment configuration between the ring and Prometheus. We focus our attention on the shape of the F ring’s azimuthal channels which were first reported by Porco et al. (Porco, C.C., Baker, E., Barbara, J., Beurle, K., Brahic, A., Burns, J.A., Charnoz, S., Cooper, N., Dawson, D.D., Del Genio, A.D., Denk, T., Dones, L., Dyudina, U., Evans, M.W., Giese, B., Grazier, K., Helfenstein, P., Ingersoll, A.P., Jacobson, R.A., Johnson, T.V., McEwen, A., Murray, C.D., Neukum, G., Owen, W.M., Perry, J., Roatsch, T., Spitale, J., Squyres, S., Thomas, P., Tiscareno, M., Turtle, E., Vasavada, A.R., Veverka, J., Wagner, R., West, R. [2005] Science, 307, 1226-1236) and numerically explored by Murray et al. (Murray, C.D., Chavez, C., Beurle, K., Cooper, N., Evans, M.W., Burns, J.A., Porco, C.C. [2005] Nature 437, 1326-1329) who found excellent agreement between Cassini’s ISS reprojected images and their numerical model via a direct comparison. We find that for anti-alignment the channels are wider and go deeper inside the ring material. From our numerical model we find a new feature, an island in the middle of the channel. This island is made up of the particles that have been perturbed the most by Prometheus and only appears when this satellite is close to apoapsis. In addition, plots of the anti-alignment configuration for different orbital stages of Prometheus are obtained and discussed here.     

New Cassini RADAR results for Saturn’s icy satellites

Cassini radar tracks on Saturn’s icy satellites through the end of the Prime Mission in 2008 have increased the number of radar albedo estimates from 10 (Ostro et al., 2006) to 73. The measurements sample diverse subradar locations (and for Dione, Rhea, and Iapetus almost always use beamwidths less than half the target angular diameters), thereby constraining the satellites’ global radar albedo distributions. The echoes result predominantly from volume scattering, and their strength is thus strongly sensitive to ice purity and regolith maturity. The combination of the Cassini data set and Arecibo 13-cm observations of Enceladus, Tethys, Dione, Rhea (Black et al., 2007), and Iapetus (Black et al., 2004) discloses an unexpectedly complex pattern of 13-to-2-cm wavelength dependence. The 13-cm albedos are generally smaller than 2-cm albedos and lack the correlation seen between 2-cm and optical geometric albedos. Enceladus and Iapetus are the most interesting cases. We infer from hemispheric albedo variations that the E-ring has a prominent effect on the 13-cm radar “lightcurve”. The uppermost trailing-side regolith is too fresh for meteoroid bombardment to have developed larger-scale heterogeneities that would be necessary to elevate the 13-cm radar albedo, whereas all of Enceladus is clean and mature enough for the 2-cm albedo to be uniformly high. For, Iapetus, the 2-cm albedo is strongly correlated with optical albedo: low for the optically dark, leading-side material and high for the optically bright, trailing-side material. However, Iapetus’ 13-cm albedo values show no significant albedo dichotomy and are several times lower than 2-cm values, being indistinguishable from the weighted mean of 13-cm albedos for main-belt asteroids, 0.15 ± 0.10. The leading side’s optically dark contaminant must be present to depths of at least one to several decimeters, so 2-cm albedos can mimic the optical dichotomy; however, it does not have to extend any deeper than that. The fact that both hemispheres of Iapetus look Asteroid-like at 13 cm means that coherent backscattering itself is not nearly as effective as it is at 2 cm. Since Iapetus’ entire surface is mature regolith, the wavelength dependence must involve composition, not structure. Either the composition is a function of depth everywhere (with electrical loss much greater at depths greater than a decimeter or two), or the intrinsic electrical loss of some pervasive constituent is much higher at 13 cm than at 2 cm. Ammonia is a candidate for such a contaminant. If ammonia’s electrical properties do not depend on frequency, and if ammonia is globally much less abundant within the upper one or two decimeters than at greater depths, then coherent backscattering would effectively be shut down at 13 cm, explaining the Asteroid-like 13-cm albedo.     

A primordial origin for the atmospheric methane of Saturn’s moon Titan

The origin of Titan’s atmospheric methane is a key issue for understanding the origin of the saturnian satellite system. It has been proposed that serpentinization reactions in Titan’s interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan’s planetesimals before its formation. Here, we point out that serpentinization reactions in Titan’s interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan’s water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan’s interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite’s planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan’s interior can be up to ∼1300 times the current mass of atmospheric methane.     

Irregular satellite capture during planetary resonance passage

The passage of Jupiter and Saturn through mutual 1:2 mean-motion resonance has recently been put forward as explanation for their relatively high eccentricities [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] and the origin of Jupiter's Trojans [Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. Additional constraints on this event based on other small-body populations would be highly desirable. Since some outer satellite orbits are known to be strongly affected by the near-resonance of Jupiter and Saturn (“the Great Inequality”; ?uk, M., Burns, J.A., 2004b. Astron. J. 128, 2518-2541), the irregular satellites are natural candidates for such a connection. In order to explore this scenario, we have integrated 9200 test particles around both Jupiter and Saturn while they went through a resonance-crossing event similar to that described by Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]. The test particles were positioned on a grid in semimajor axes and inclinations, while their initial pericenters were put at just 0.01 AU from their parent planets. The goal of the experiment was to find out if short-lived bodies, spiraling into the planet due to gas drag (or alternatively on orbits crossing those of the regular satellites), could have their pericenters raised by the resonant perturbations. We found that about 3% of the particles had their pericenters raised above 0.03 AU (i.e. beyond Iapetus) at Saturn, but the same happened for only 0.1% of the particles at Jupiter. The distribution of surviving particles at Saturn has strong similarities to that of the known irregular satellites. If saturnian irregular satellites had their origin during the 1:2 resonance crossing, they present an excellent probe into the early Solar System's evolution. We also explore the applicability of this mechanism for Uranus, and find that only some of the uranian irregular satellites have orbits consistent with resonant pericenter lifting. In particular, the more distant and eccentric satellites like Sycorax could be stabilized by this process, while closer-in moons with lower eccentricity orbits like Caliban probably did not evolve by this process alone.     

How the Enceladus dust plume feeds Saturn’s E ring

Pre-Cassini models of Saturn’s E ring [Horányi, M., Burns, J., Hamilton, D., 1992. Icarus 97, 248-259; Juhász, A., Horányi, M., 2002. J. Geophys. Res. 107, 1-10] failed to reproduce its peculiar vertical structure inferred from Earth-bound observations [de Pater, I., Martin, S.C., Showalter, M.R., 2004. Icarus 172, 446-454]. After the discovery of an active ice-volcanism of Saturn’s icy moon Enceladus the relevance of the directed injection of particles for the vertical ring structure of the E ring was swiftly recognised [Juhász, A., Horányi, M., Morfill, G.E., 2007. Geophys. Res. Lett. 34, L09104; Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., Economou, T., Schmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. However, simple models for the delivery of particles from the plume to the ring predict a too small vertical ring thickness and overestimate the amount of the injected dust.Here we report on numerical simulations of grains leaving the plume and populating the dust torus of Enceladus. We run a large number of dynamical simulations including gravity and Lorentz force to investigate the earliest phase of the ring particle life span. The evolution of the electrostatic charge carried by the initially uncharged grains is treated selfconsistently. Freshly ejected plume particles are moving in almost circular orbits because the Enceladus orbital speed exceeds the particles’ ejection speeds by far. Only a small fraction of grains that leave the Hill sphere of Enceladus survive the next encounter with the moon. Thus, the flux and size distribution of the surviving grains, replenishing the ring particle reservoir, differs significantly from the flux and size distribution of the particles freshly ejected from the plume. Our numerical simulations reproduce the vertical ring profile measured by the Cassini Cosmic Dust Analyzer (CDA) [Kempf, S., Beckmann, U., Moragas-Klostermeyer, G., Postberg, F., Srama, R., EconoDmou, T., Smchmidt, J., Spahn, F., Grün, E., 2008. Icarus 193, 420-437]. From our simulations we calculate the deposition rates of plume particles hitting Enceladus’ surface. We find that at a distance of 100 m from a jet a 10 m sized ice boulder should be covered by plume particles in 10⁵-10⁶ years.     

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.     

The fate of ejecta from Hyperion

Ejecta from Saturn's moon Hyperion are subject to powerful perturbations from nearby Titan, which control their ultimate fate. We have performed numerical integrations to simulate a simplified system consisting of Saturn (including optical flattening as well as dynamical oblateness), its main ring system (treated as a massless flat annulus), the moons Tethys, Dione, Titan, Hyperion, and Iapetus, and the Sun (treated simply as a massive satellite). At several different points in Hyperion's orbit, 1050 massless particles, more or less evenly distributed over latitude and longitude, were ejected radially outward from 1 km above Hyperion's mean radius at speeds 10% faster than escape speed from Hyperion. Most of these particles were removed within the first few thousand years, but ∼3% of them survived the entire 100,000-year duration of the simulations. Ejecta from Hyperion are much more widely scattered than previously thought, and can cross the orbits of all of Saturn's satellites. About 9% of all the particles escaped from the saturnian system, but Titan accreted ∼78% of the total, while Hyperion reaccreted only ∼5%. This low efficiency of reaccretion may help to account for Hyperion's small size and rugged shape. Only ∼1% of all the particles hit other satellites, and another ∼1% impacted Saturn itself, while ∼3% of them struck its main rings. The high proportion of impacts into Saturn's rings is surprising; these collisions show a broad decline in impact speed with time, suggesting that Hyperion ejecta gradually spread inwards. Additional simulations were used to investigate the dependence of ejecta evolution on launch speed, the mass of Hyperion, and the presence of the Sun. In general, the wide distribution of ejecta from Hyperion suggests that it does contribute to “Population II” craters on the inner satellites of Saturn. Ejecta which escape from a satellite into temporary orbit about its planet, but later reimpact into the same moon or another one produce “poltorary” impacts, intermediate in character between primary and secondary impacts. It may be possible to distinguish poltorary craters from primary and secondary craters on the basis of morphology.     

Total particulate mass in Enceladus plumes and mass of Saturn’s E ring inferred from Cassini ISS images

The eclipse mosaic (PIA08329) of the Saturn system, taken on September 15, 2006 when Cassini was in Saturn’s shadow, contains numerous color images of the Enceladus plume and the E ring at phase angles ranging from 173° to 179°. These forward-scattering observations sample the diffraction peak for particle radii in the 1–5 μm range. The phase angle dependence and total brightness are sensitive indicators of the total mass of solid material in the plume. We fit the data with a variety of particle shapes and size distributions, and find that the median radius of the equivalent-volume sphere is 3.1 μm, with an uncertainty of ±0.5 μm. The total mass of particles in the plume is (1.45 ± 0.5) × 10⁵ kg. We have not considered variations with altitude in the particle size and shape distribution, and we leave that for another paper. We find that the brightness of the E ring varies with position in the orbit, not only because of the viewing geometry, e.g., variations in phase angle, but also because of some unknown intrinsic variability. The total mass of solid material in the E ring is (12 ± 5.5) × 10⁸ kg. For the plume, the production rate of particles – the mass per unit time leaving the vents is 51 ± 18 kg s⁻¹. We estimate that 9% of these particles are escaping from Enceladus, implying lifetimes of ∼8 years for the E ring particles. Based on three comparisons with vapor amounts from ultraviolet spectroscopy, the ice/vapor ratio is in the range 0.35–0.70. This high ratio poses a problem for theories in which particles form by condensation from the gas phase, and could indicate that particles are formed as spray from a liquid reservoir.     

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.     

Chaotic motions of prometheus and pandora

Recent HST images of the saturnian satellites Prometheus and Pandora show that their longitudes deviate from predictions of ephemerides based on Voyager images. Currently Prometheus is lagging and Pandora leading these predictions by somewhat more than 20°. We show that these discrepancies are fully accounted for by gravitational interactions between the two satellites. These peak every 24.8 days at conjunctions and excite chaotic perturbations. The Lyapunov exponent for the Prometheus-Pandora system is of order 0.3 year⁻¹ for satellite masses based on a nominal density of 0.63 g cm⁻³. Interactions are strongest when the orbits come closest together. This happens at intervals of 6.2 years when their apses are antialigned. In this context, we note the sudden changes of opposite signs in the mean motions of Prometheus and Pandora at the end of 2000 occurred around the time their apsidal lines were antialigned.     

Origin of chaos in the Prometheus-Pandora system

We demonstrate that the chaotic orbits of Prometheus and Pandora are due to interactions associated with the 121:118 mean motion resonance. Differential precession splits this resonance into a quartet of components equally spaced in frequency. Libration widths of the individual components exceed the splitting, resulting in resonance overlap which causes the chaos. Mean motions of Prometheus and Pandora wander chaotically in zones of width 1.8 and 3.1 deg yr⁻¹, respectively. A model with 1.5 degrees of freedom captures the essential features of the chaotic dynamics. We use it to show that the Lyapunov exponent of 0.3 yr⁻¹ arises because the critical argument of the dominant member of the resonant quartet makes approximately two separatrix crossings every 6.2 year precessional cycle.     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

Rotational dynamics of planetary satellites: A survey of regular and chaotic behavior

The problem of observability of chaotic regimes in the rotation of planetary satellites is studied. The analysis is based on the inertial and orbital data available for all satellites discovered up to now. The Lyapunov spectra of the spatial chaotic rotation and the full range of variation of the spin rate are computed numerically by integrating the equations of the rotational motion; the initial data are taken inside the main chaotic layer near the separatrices of synchronous resonance in phase space. The model of a triaxial satellite in a fixed elliptic orbit is adopted. A short Lyapunov time along with a large range of variation of the spin rate are used as criteria for observability of the chaotic motion. Independently, analysis of stability of the synchronous state with respect to tilting the axis of rotation provides a test for the physical opportunity for a satellite to rotate chaotically. Finally, a calculation of the times of despinning due to tidal evolution shows whether a satellite's spin could evolve close to the synchronous state. Apart from Hyperion, already known to rotate chaotically, only Prometheus and Pandora, the 16th and 17th satellites of Saturn, pass all these four tests.     

Accretion of Saturn’s mid-sized moons during the viscous spreading of young massive rings: Solving the paradox of silicate-poor rings versus silicate-rich moons

The origin of Saturn’s inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn’s rings is debated. Charnoz et al. [Charnoz, S., Salmon J., Crida A., 2010. Nature 465, 752–754] introduced the idea that the smallest inner moons could form from the spreading of the rings’ edge while Salmon et al. [Salmon, J., Charnoz, S., Crida, A., Brahic, A., 2010. Icarus 209, 771–785] showed that the rings could have been initially massive, and so was the ring’s progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested by Canup [Canup, R., 2010. Nature 468, 943–946]. However, the process driving mid-sized moon accretion from the icy debris disks has not been investigated in details. In particular, Canup’s (2010) model does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea’s mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Q_p < 2000) is required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites’ cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process). These moons may be either mostly icy, or, if they contain a significant amount of rock, already differentiated from the ice without the need for radiogenic heating. The resulting inner mid-sized moons may be significantly younger than the Solar System and a ∼1 Gyr formation delay is possible between Mimas and Rhea. The rings resulting from this process would evolve to a state compatible with current mass estimates of Saturn’s rings, and nearly devoid of silicates, apart from isolated silicate chunks coated with ice, interpreted as today Saturn’s rings’ propellers and ring-moons (like Pan or Daphnis).     

Simulation and analysis of the dynamics of binary near-Earth Asteroid (66391) 1999 KW4

Near-Earth Asteroid (66391) 1999 KW4 was the subject of the recently published first extensive radar imaging, shape and mutual orbit modeling, and physical and dynamical characterization of a binary asteroid. In this paper we present in detail our numerical simulation of KW4 behind that work. Our propagations of the system with some variation in estimated parameters cover the set of KW4's possible current dynamical states consistent with the body models and other information obtained directly from the observations. We also apply our implementation of this simulation capability to address some of the dynamical mechanisms by which KW4 may be moved into the more energetically excited of those possible current states, particularly solar gravity interaction. Through comparison of the results with certain features of the observation data, we conclude that the actual KW4 system is not in the most energetically relaxed configuration but must be moderately excited. The system occupies a generalized Cassini state 2 which is different from that considered in most previously published treatments of Cassini states in that it involves co-precession of the primary's spin frame and the mutual orbit rather than co-precession of a satellite's spin frame and that satellite's orbit about the primary. We present a simple analytical theory describing the system's dynamics, which should be applicable to any other binary systems, of which KW4 is representative, in which a massive, roughly oblate primary is spinning rapidly relative to the rate of its mutual orbit with an on-average synchronous, elongated secondary. We examine separately both the effect of the larger binary component's oblateness, and the effect of the smaller component's roughly triaxial ellipsoid shape. The simple analytical formulae obtained agree with full-detail numerical simulation results, and can be used for remote estimation of binary mass properties from observed system motion.     

Normal modes of synchronous rotation

The dynamics of synchronous rotation and physical librations are revisited in order to establish a conceptually simple and general theoretical framework applicable to a variety of problems. Our motivation comes from disagreements between the results of numerical simulations and those of previous theoretical studies, and also because different theoretical studies disagree on basic features of the dynamics. We approach the problem by decomposing the orientation matrix of the body into perfectly synchronous rotation and deviation from the equilibrium state. The normal modes of the linearized equations are computed in the case of a circular satellite orbit, yielding both the periods and the eigenspaces of three librations. Libration in longitude decouples from the other two, vertical modes. There is a fast vertical mode with a period very close to the average rotational period. It corresponds to tilting the body around a horizontal axis while retaining nearly principal-axis rotation. In the inertial frame, this mode appears as nutation and free precession. The other vertical mode, a slow one, is the free wobble. The effects of the nodal precession of the orbit are investigated from the point of view of Cassini states. We test our theory using numerical simulations of the full equations of the dynamics and discuss the disagreements among our study and previous ones. The numerical simulations also reveal that in the case of eccentric orbits large departures from principal-axis rotation are possible due to a resonance between free precession and wobble. We also revisit the history of the Moon's rotational state and show that it switched from one Cassini state to another when it was at 46.2 Earth radii. This number disagrees with the value 34.2 derived in a previous study.     

Fates of satellite ejecta in the Saturn system

We use conventional numerical integrations to assess the fates of impact ejecta in the Saturn system. For specificity we consider impact ejecta launched from four giant craters on three satellites: Herschel on Mimas, Odysseus and Penelope on Tethys, and Tirawa on Rhea. Speeds, trajectories, and size of the ejecta are consistent with impact on a competent surface (“spalls”) and into unconsolidated regolith. We do not include near-field effects, jetting, or effects peculiar to highly oblique impact. Ejecta are launched at velocities comparable to or exceeding the satellite's escape speed. Most ejecta are swept up by the source moon on time-scales of a few to several decades, and produce craters no larger than 19 km in diameter, with typical craters in the range of a few km. As much as 17% of ejecta reach satellites other than the source moon. Our models generate cratering patterns consistent with a planetocentric origin of most small impact craters on the saturnian icy moons, but the predicted craters tend to be smaller than putative Population II craters. We conclude that ejecta from the known giant craters in the saturnian system do not fully account for Population II craters.