Lunar shape does not record a past eccentric orbit

The Moon has long been known to have an overall shape not consistent with expected past tidal forces. It has recently been suggested (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) that the present lunar moments of inertia indicate a past high-eccentricity orbit and, possibly, a past non-synchronous spin-orbit resonance. Here I show that the match between the lunar shape and the proposed orbital and spin states is much less conclusive than initially proposed. Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) spin and shape evolution scenarios also completely ignore the physics of the capture into such resonances, which require prior permanent deformation, as well as tidal despinning to the relevant resonance. If the early lunar orbit was eccentric, the Moon would have been rotating at an equilibrium non-synchronous rate determined by it eccentricity. This equilibrium supersynchronous rotation would be much too fast to allow a synchronous spin-orbit lock at e = 0.49, while the capture into the 3:2 resonance is possible only for a very constrained lunar eccentricity history and assuming some early permanent lunar tri-axiality. Here I show that large impacts in the early history of the Moon would have frequently disrupted this putative resonant rotation, making the rotation and eccentricity solutions of Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) unstable. I conclude that the present lunar shape cannot be used to support the hypothesis of an early eccentric lunar orbit.     

Librations of the Galilean satellites: The influence of global internal liquid layers

The four Galilean satellites are thought to harbor one or even two global internal liquid layers beneath their surface layer. The iron core of Io and Ganymede is most likely (partially) liquid and also the core of Europa may be liquid. Furthermore, there are strong indications for the existence of a subsurface ocean in Europa, Ganymede, and Callisto. Here, we investigate whether libration observations can be used to prove the existence of these liquid layers and to constrain the thickness of the overlying solid layers. For Io, the presence of a small liquid core increases the libration of the mantle by a few percent with respect to an entirely solid Io and mantle libration observations could be used to determine the mantle thickness with a precision of several tens of kilometers given that the libration amplitude can be measured with a precision of 1 m. For Europa, Ganymede, and Callisto, the presence of a water ocean close to the surface increases by at least an order of magnitude the ice shell libration amplitude with respect to an entirely solid satellite. The shell libration depends essentially on the shell thickness and to a minor extent on the density difference between the ocean and the ice shell. The possible presence of a liquid core inside Europa and Ganymede has no noticeable influence on their shell libration. For a precision of several meters on the libration measurements, in agreement with the expected accuracy with the NASA/ESA EJSM orbiter mission to Europa and Ganymede, an error on the shell thickness of a few tens kilometers is expected. Therefore, libration measurements can be used to detect liquid layers such as Io’s core or water subsurface oceans in Europa, Ganymede, and Callisto and to constrain the thickness of the overlying solid surface layers.     

The fate of primordial lunar Trojans

Multiple large impact basins on the lunar nearside formed in a relatively-short interval around 3.8-3.9 Gyr ago, in what is known as the Lunar Cataclysm (LC; also known as Late Heavy Bombardment). It is widely thought that this impact bombardment has affected the whole Solar System or at least all the inner planets. But with non-lunar evidence for the cataclysm being relatively weak, a geocentric cause of the Lunar Cataclysm cannot yet be completely ruled out [Ryder, G., 1990. Eos 71, 313, 322-323]. In principle, late destabilization of an additional Earth satellite could result in its tidal disruption during a close lunar encounter (cf. [Asphaug, E., Agnor, C.B., Williams, Q., 2006. Nature 439, 155-160]). If the lost satellite had D>500 km, the resulting debris can form multiple impact basins in a relatively short time, possibly explaining the LC. Canup et al. [Canup, R.M., Levison, H.F., Stewart, G.R., 1999. Astron. J. 117, 603-620] have shown that any additional satellites of Earth formed together with (and external to) the Moon would be unable to survive the rapid initial tidally-driven expansion of lunar orbit. Here we explore the fate of objects trapped in the lunar Trojan points, and find that small lunar Trojans can survive the Moon's orbital evolution until they and the Moon reach 38 Earth radii, at which point they are destabilized by a strong solar resonance. However, the dynamics of Trojans containing enough mass to cause the LC (diameters >150 km) is more complex; we find that such objects do not survive the passage through a weaker solar resonance at 27 Earth radii. This distance was very likely reached by the Moon long before the LC, which seems to rule out the disruption of lunar Trojans as a cause of the LC.     

Kinetic modeling of sodium in the lunar exosphere

Our knowledge about the lunar environment is based on a large volume of ground-based, remote, and in situ observations. These observations have been conducted at different times and sampled different pieces of such a complex system as the surface-bound exosphere of the Moon. Numerical modeling is the tool that can link results of these separate observations into a single picture. Being validated against previous measurements, models can be used for predictions and interpretation of future observations results.     

Coupled thermal-orbital evolution of the early Moon

Coupled thermal-orbital histories of early lunar evolution are considered in a simple model. We consider a plagioclase lid, overlying a magma ocean, overlying a solid mantle. Tidal dissipation occurs in the plagioclase lid and heat transport is by conduction and melt migration. We find that large orbital eccentricities can be obtained in this model. We discuss possible consequences of this phase of large eccentricities for the shape of the Moon and geochronology of lunar samples. We find that the orbit can pass through the shape solution of Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M. [2006]. Science 313, 652), but we argue that the shape cannot be maintained against elastic deformation as the orbit continues to evolve.     

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.     

Constraints on the source of lunar cataclysm impactors

Multiple impact basins formed on the Moon about 3.8 Gyr ago in what is known as the lunar cataclysm or Late Heavy Bombardment. Many workers currently interpret the lunar cataclysm as an impact spike primarily caused by main-belt asteroids destabilized by delayed planetary migration. We show that morphologically fresh (class 1) craters on the lunar highlands were mostly formed during the brief tail of the cataclysm, as they have absolute crater number density similar to that of the Orientale basin and ejecta blanket. The connection between class 1 craters and the cataclysm is supported by the similarity of their size-frequency distribution to that of stratigraphically-identified Imbrian craters. Majority of lunar craters younger than the Imbrium basin (including class 1 craters) thus record the size-frequency distribution of the lunar cataclysm impactors. This distribution is much steeper than that of main-belt asteroids. We argue that the projectiles bombarding the Moon at the time of the cataclysm could not have been main-belt asteroids ejected by purely gravitational means.     

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.     

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.     

Iapetus' geophysics: Rotation rate, shape, and equatorial ridge

Iapetus has preserved evidence that constrains the modeling of its geophysical history from the time of its accretion until now. The evidence is (a) its present 79.33-day rotation or spin rate, (b) its shape that corresponds to the equilibrium figure for a hydrostatic body rotating with a period of ∼16 h, and (c) its high, equatorial ridge, which is unique in the Solar System. This paper reports the results of an investigation into the coupling between Iapetus' thermal and orbital evolution for a wide range of conditions including the spatial distributions with time of composition, porosity, short-lived radioactive isotopes (SLRI), and temperature. The thermal model uses conductive heat transfer with temperature-dependent conductivity. Only models with a thick lithosphere and an interior viscosity in the range of about the water ice melting point can explain the observed shape. Short-lived radioactive isotopes provide the heat needed to decrease porosity in Iapetus' early history. This increases thermal conductivity and allows the development of the strong lithosphere that is required to preserve the 16-h rotational shape and the high vertical relief of the topography. Long-lived radioactive isotopes and SLRI raise internal temperatures high enough that significant tidal dissipation can start, and despin Iapetus to synchronous rotation. This occurred several hundred million years after Iapetus formed. The models also constrain the time when Iapetus formed because the successful models are critically dependent upon having just the right amount of heat added by SLRI decay in this early period. The amount of heat available from short-lived radioactivity is not a free parameter but is fixed by the time when Iapetus accreted, by the canonical concentration of ²⁶Al, and, to a lesser extent, by the concentration of ⁶⁰Fe. The needed amount of heat is available only if Iapetus accreted between 2.5 and 5.0 Myr after the formation of the calcium aluminum inclusions as found in meteorites. Models with these features allow us to explain Iapetus' present synchronous rotation, its fossil 16-h shape, and the context within which the equatorial ridge arose.     

Was there ever a Terminal Lunar Cataclysm?: With lunar viscosity arguments

About 30 years ago there was a suggestion by several able scientists at the California Institute of Technology that the Moon had undergone a Terminal Lunar Cataclysm. This meant that most of the early impact cratering had been concentrated strongly at about the time of formation of the Imbrium basin. This solution was discussed in many papers and the idea of a cataclysm gradually faded away. In about 1990 it was again revived by several scientists. The idea of a Terminal Lunar Cataclysm at about the time the Imbrium basin was formed was advanced albeit in a somewhat different manner. The present paper has been written to analyze the various observations and interpretations that have been advanced to permit a cataclysm. It is concluded that the three main proposals, which, if correct, would have permitted a cataclysm to have occurred, are each faulty and not consistent with such a cataclysm. To demonstrate this conclusion it was necessary to determine absolute ages of various lunar features. This meant, in part, determinations of the existence and nature of lunar crustal viscosity consistent with times of formation of six lunar basins. The results of such studies yielded an internally consistent model which requires a long period from the original formation of the Moon at about 4.5 byr to a time slightly earlier than that of the formation of the Imbrium basin at about 3.84 byr. On this model there is no indication of a clustering of events and it is concluded that a Terminal Lunar Cataclysm never occurred.     

Rebuttal to the comment by Malhotra and Strom on “Constraints on the source of lunar cataclysm impactors”

?uk et al. (?uk, M. Gladman, B.J., Stewart, S.T. [2010]. Icarus 207 590-594) concluded that the the lunar cataclysm (late heavy bombardment) was recorded in lunar Imbrian era craters, and that their size distribution is different from that of main belt asteroids (which may have been the dominant pre-Imbrian impactors). This result would likely preclude the asteroid belt as the direct source of lunar cataclysm impactors. Malhotra and Strom (Malhotra, R., Strom, R.G. [2011]. Icarus) maintain that the lunar impactor population in the Imbrian era was the same as in Nectarian and pre-Nectarian periods, and this population had a size distribution identical to that of main belt asteroids. In support of this claim, they present an Imbrian size distribution made from two data sets published by Wilhelms et al. (Wilhelms, D.E., Oberbeck, V.R., Aggarwal, H.R. [1978]. Proc. Lunar Sci. Conf. 9, 3735-3762). However, these two data sets cannot be simply combined as they represent areas of different ages and therefore crater densities. Malhotra and Strom (Malhotra, R., Strom, R.G. [2011]. Icarus) differ with the main conclusion of Wilhelms et al. (Wilhelms, D.E., Oberbeck, V.R., Aggarwal, H.R. [1978]. Proc. Lunar Sci. Conf. 9, 3735-3762) that the Nectarian and Imbrian crater size distributions were different. We conclude that the available data indicate that the lunar Imbrian-era impactors had a different size distribution from the older ones, with the Imbrian impactor distribution being significantly richer in small impactors than that of older lunar impactors or current main-belt asteroids.     

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.     

Polar wander of Mars: Evidence from giant impact basins

We investigate the polar wander of Mars in the last ∼4.2 Ga. We identify two sets of basins from the 20 giant impact basins reported by Frey [Frey, H., 2008. Geophys. Res. Lett. 35, L13203] which trace great circles on Mars, and propose that the great circles were the prevailing equators of Mars at the impact times. Monte Carlo tests are conducted to demonstrate that the two sets of basins are most likely not created by random impacts. Also, fitting 63,771 planes to randomly selected sets of 5, 6, or 7 basins indicated that the identified two sets are unique. We propose three different positions for the rotation pole of Mars, besides the present one. Accordingly, Tharsis bulge was initially formed at ∼50 N and moved toward the equator while rotating counterclockwise due to the influence of the two newly forming volcanic constructs, Alba Patera and Elysium Rise. The formation of the giant impact basins, subsequent mass concentrations (mascons) in Argyre, Isidis, and Utopia basins, and surface masses of volcanic mountains such as Ascraeus, Pavonis, Arsia and Olympus, caused further polar wander which rotated Tharsis bulge clockwise to arrive at its present location. The extensive polar motion of Mars during 4.2-3.9 Ga implies a weak lithosphere on a global scale, deduced from a total of 72,000 polar wander models driven by Tharsis bulge, Alba Patera and Elysium Rise as the major mass perturbations. Different compensation states, 0-100%, are examined for each of the surface loads, and nine different thicknesses are considered for an elastic lithosphere. The lithosphere must have been very weak, with an elastic thickness of less than 5 km, if the polar wander was driven by these mass perturbations.     

Analytical model of the long-period forced longitude librations of Mercury

The shaking of Mercury’s orbit by the planets forces librations in longitude in addition to those at harmonics of the orbital period that have been used to detect Mercury’s molten core. We extend the analytical formulation of Peale et al. (Peale, S.J., Margot, J.L., Yseboodt, M. [2009]. Icarus 199, 1-8) in order to provide a convenient means of determining the amplitudes and phases of the forced librations without resorting to numerical calculations. We derive an explicit relation between the amplitude of each forced libration and the moment of inertia parameter (B-A)/C_m. Far from resonance with the free libration period, the libration amplitudes are directly proportional to (B-A)/C_m. Librations with periods close to the free libration period of ∼12 years may have measurable (∼arcsec) amplitudes. If the free libration period is sufficiently close to Jupiter’s orbital period of 11.86 years, the amplitude of the forced libration at Jupiter’s period could exceed the 35 arcsec amplitude of the 88-day forced libration. We also show that the planetary perturbations of the mean anomaly and the longitude of pericenter of Mercury’s orbit completely determine the libration amplitudes.While these signatures do not affect spin rate at a detectable level (as currently measured by Earth-based radar), they have a much larger impact on rotational phase (affecting imaging, altimetry, and gravity sensors). Therefore, it may be important to consider planetary perturbations when interpreting future spacecraft observations of the librations.     

Central magnetic anomalies of Nectarian-aged lunar impact basins: Probable evidence for an early core dynamo

A re-examination of all available low-altitude LP magnetometer data confirms that magnetic anomalies are present in at least four Nectarian-aged lunar basins: Moscoviense, Mendel-Rydberg, Humboldtianum, and Crisium. In three of the four cases, a single main anomaly is present near the basin center while, in the case of Crisium, anomalies are distributed in a semi-circular arc about the basin center. These distributions, together with a lack of other anomalies near the basins, indicate that the sources of the anomalies are genetically associated with the respective basin-forming events. These central basin anomalies are difficult to attribute to shock remanent magnetization of a shocked central uplift and most probably imply thermoremanent magnetization of impact melt rocks in a steady magnetizing field. Iterative forward modeling of the single strongest and most isolated anomaly, the northern Crisium anomaly, yields a paleomagnetic pole position at 81° ± 19°N, 143° ± 31°E, not far from the present rotational pole. Assuming no significant true polar wander since the Crisium impact, this position is consistent with that expected for a core dynamo magnetizing field. Further iterative forward modeling demonstrates that the remaining Crisium anomalies can be approximately simulated assuming a multiple source model with a single magnetization direction equal to that inferred for the northernmost anomaly. This result is most consistent with a steady, large-scale magnetizing field. The inferred mean magnetization intensity within the strongest basin sources is ∼1 A/m assuming a 1-km thickness for the source layer. Future low-altitude orbital and surface magnetometer measurements will more strongly constrain the depth and/or thicknesses of the sources.     

Tidal dynamical considerations constrain the state of an ocean on Enceladus

In previous work, solutions to the non-dissipative Laplace Tidal Equations (LTE) were used to provide bounds on the heat generated by the response of a subsurface ocean on Enceladus to an obliquity component of tidal forces. Here we improve these bounds using solutions from the LTE with a generic dissipation term explicitly added. We find solutions for a wider range of ocean tidal responses that include both unstratified (barotropic) and stratified (baroclinic) flow responses to obliquity as well as eccentricity components of the tidal forces. We consolidate the results in three ocean tidal scenarios on Enceladus that can explain the high heat fluxes (∼7 mW/m² globally averaged) inferred from measurements by the Cassini spacecraft: (1) a deep (1-50 km) barotropic ocean responding to obliquity tidal forces, where obliquity is at least 0.1°; (2) a shallow (∼360 m) barotropic ocean responding to eccentricity tidal forces; (3) a stratified (baroclinic) ocean responding to eccentricity tidal forces where the density-weighted “equivalent depth” (typically much smaller than the ocean’s physical depth) is near 360 m. The ocean is assumed to be global, but extensions for a semi-global case are also described. A more general result which is independent of the specific scenarios proposed is that an ocean attempting to freeze (with an associated decrease in its liquid depth, which affects the ocean’s dynamical response to the tidal forcing) must first pass through resonant configurations with a greatly increased generation of ocean tidal heat (exceeding 1 W/m² to 1 kW/m²) that would act to halt further freezing and stagnate the ocean state in this configuration so long as there is still orbital energy to provide the tidal forces. With an additional assumption that the ocean has evolved from a more energetic state where the depth of the liquid ocean was greater, we obtain the three scenarios proposed.     

Obliquity tides do not significantly heat Enceladus

Recently, Tyler [Tyler, R.H., 2009. Geophys. Res. Lett. 36, L15205; Tyler, R., 2011. Icarus, 211, 770-779] proposed that the tide due to an obliquity of greater than 0.1° might drive resonant flow in a liquid ocean at Enceladus, and that dissipation of the ocean’s kinetic energy may be an alternate source for the observed global heat flux. While there is currently no measurement of Enceladus’ obliquity, dissipation is expected to drive the spin pole to a Cassini state. Under this assumption, we find that Enceladus should occupy Cassini state 1 and that the obliquity of Enceladus should be less than 0.0015° for values of the degree-2 gravity coefficient C_2,2 between 1.0 × 10⁻³ and 2.5 × 10⁻³. Unless there is a significant free obliquity or the gravity coefficient C_2,2 has been significantly overestimated, it is unlikely that obliquity-driven flow in a subsurface ocean is the source of the extreme heat on Enceladus.     

Shapes of the saturnian icy satellites and their significance

The sizes and shapes of six icy saturnian satellites have been measured from Cassini Imaging Science Subsystem (ISS) data, employing limb coordinates and stereogrammetric control points. Mimas, Enceladus, Tethys, Dione and Rhea are well described by triaxial ellipsoids; Iapetus is best represented by an oblate spheroid. All satellites appear to have approached relaxed, equilibrium shapes at some point in their evolution, but all support at least 300 m of global-wavelength topography. The shape of Enceladus is most consistent with a homogeneous interior. If Enceladus is differentiated, its shape and apparent relaxation require either lateral inhomogeneities in an icy mantle and/or an irregularly shaped core. Iapetus supports a fossil bulge of over 30 km, and provides a benchmark for impact modification of shapes after global relaxation. Satellites such as Mimas that have smoother limbs than Iapetus, and are expected to have higher impact rates, must have relaxed after the shape of Iapetus was frozen.     

Tidal encounters of ellipsoidal granular asteroids with planets

We investigate planetary fly-bys of asteroids using an approximate volume-averaged method that offers a relatively simple, but very flexible, approach to study the rotational dynamics of ellipsoids. The asteroid is considered to be a deformable, prolate ellipsoid, with its interior being modeled as a rigid-granular material. Effects due to the asteroid's rotation, its self-gravity and gravitational interaction with the planet are included. Using a simplified approach allows us to explore in detail the mechanics of asteroid's deformations and disruptions during planetary encounters. We also compare our results with those obtained by Richardson et al. [Richardson, D.C., Bottke Jr., W.F., Love, S.G., 1998. Icarus 134, 47-76] who used a large numerical code. We find that many of the features reported by them can indeed be captured by our rather simple methodology, and we discuss the reasons why some of our results differ from theirs.