A contraction model for the flattening and equatorial ridge of Iapetus

Others have explained the excess flattening of Iapetus by a model in which the moon formed at a high spin rate, achieved isostatic equilibrium by very rapid interior heating caused by short-lived radioactive isotopes (SLRI), and subsequently cooled, locking in the excess flattening with respect to an equilibrium shape at its present spin rate. Here we propose an alternate model that does not require an unusually high initial spin rate or the SLRI. The initial formation of Iapetus results in a slightly oblate spheroid with porosity >10%. Radioactive heating by long-lived isotopes warms the interior to about 200 K, at which point it becomes ductile and the interior compacts by 10%, while the 120 km-thick exterior shell remains strong. The shell must deform to match the reduced volume of the ductile interior, and we propose that this deformation occurs along the equator, perhaps focused by a thinner equatorial shell. The final shape of the collapsed sphere matches the observed shape of Iapetus today, described as an oblate ellipse, except along the equator where strain concentration forms a broad ridge. To maintain this non-equilibrium shape, the thickness of the shell must exceed 120 km. Testing the equatorial focusing hypothesis will require a model that includes non-linear processes to account for the finite yield strength of the thick lithosphere. Nevertheless, we show that the stress in the lithosphere generated by the contraction of the interior is about 3 times greater than the stress needed to deform the lithosphere, so some type of lithospheric deformation is expected.     

Impact basin relaxation at Iapetus

We investigate impact basin relaxation on Iapetus by combining a 3D thermal evolution model (Robuchon, G., Choblet, G., Tobie, G., Cadek, O., Sotin, C., Grasset, O. [2010]. Icarus 207, 959-971) with a spherical axisymmetric viscoelastic relaxation code (Zhong, S., Paulson, A., Wahr, J. [2003]. Geophys. J. Int. 155, 679-695). Due to the progressive cooling of Iapetus, younger basins relax less than older basins. For an ice reference viscosity of 10¹⁴ Pa s, an 800 km diameter basin relaxes by 30% if it formed in the first 50 Myr but by 10% if it formed at 1.2 Gyr. Bigger basins relax more rapidly than smaller ones, because the inferred thickness of the ice shell exceeds the diameter of all but the largest basins considered. Stereo topography shows that all basins 600 km in diameter or smaller are relaxed by 25% or less. Our model can match the relaxation of all the basins considered, within error, by assuming a single basin formation age (4.36 Ga for our nominal viscosity). This result is consistent with crater counts, which show no detectable age variation between the basins examined.     

Revisiting the thermal inertia of Iapetus: Clues to the thickness of the dark material

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m⁻² K⁻¹ s^−1/2 while the light material varies between 15 and 25 J m⁻² K⁻¹ s^−1/2. Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.     

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.     

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.     

Deciphering the origin of the regular satellites of gaseous giants - Iapetus: The Rosetta ice-moon

Ever since their discovery the regular satellites of Jupiter and Saturn have held out the promise of providing an independent set of observations with which to test theories of planet formation. Yet elucidating their origins has proven elusive. Here we show that Iapetus can serve to discriminate between satellite formation models. Its accretion history can be understood in terms of a two-component gaseous subnebula, with a relatively dense inner region, and an extended tail out to the location of the irregular satellites, as in the SEMM model of Mosqueira and Estrada (2003a,b) (Mosqueira, I., Estrada, P.R. [2003a]. Icarus 163, 198-231; Mosqueira, I., Estrada, P.R. [2003b]. Icarus 163, 232-255). Following giant planet formation, planetesimals in the feeding zone of Jupiter and Saturn become dynamically excited, and undergo a collisional cascade. Ablation and capture of planetesimal fragments crossing the gaseous circumplanetary disks delivers enough collisional rubble to account for the mass budgets of the regular satellites of Jupiter and Saturn. This process can result in rock/ice fractionation as long as the make up of the population of disk crossers is non-homogeneous, thus offering a natural explanation for the marked compositional differences between outer solar nebula objects and those that accreted in the subnebulae of the giant planets. For a given size, icy objects are easier to capture and to ablate, likely resulting in an overall enrichment of ice in the subnebula. Furthermore, capture and ablation of rocky fragments become inefficient far from the planet for two reasons: the gas surface density of the subnebula is taken to drop outside the centrifugal radius, and the velocity of interlopers decreases with distance from the planet. Thus, rocky objects crossing the outer disks of Jupiter and Saturn never reach a temperature high enough to ablate either due to melting or vaporization, and capture is also greatly diminished there. In contrast, icy objects crossing the outer disks of each planet ablate due to the melting and vaporization of water-ice. Consequently, our model leads to an enhancement of the ice content of Iapetus, and to a lesser degree those of Titan, Callisto and Ganymede, and accounts for the (non-stochastic) compositions of these large, low-porosity outer regular satellites of Jupiter and Saturn. For this to work, the primordial population of planetesimals in the Jupiter-Saturn region must be partially differentiated, so that the ensuing collisional cascade produces an icy population of ?1 m size fragments to be ablated during subnebula crossing. We argue this is likely because the first generation of solar nebula ∼10 km planetesimals in the Jupiter-Saturn region incorporated significant quantities of ²⁶Al. This is the first study successfully to provide a direct connection between nebula planetesimals and subnebulae mixtures with quantifiable and observable consequences for the bulk properties of the regular satellites of Jupiter and Saturn, and the only explanation presently available for Iapetus’ low density and ice-rich composition.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.     

Hydrocarbons on Saturn's satellites Iapetus and Phoebe

Material of low geometric albedo (pV?0.1) is found on many objects in the outer Solar System, but its distribution in the saturnian satellite system is of special interest because of its juxtaposition with high-albedo ice. In the absence of clear, diagnostic spectral features, the composition of this low-albedo (or “dark”) material is generally inferred to be carbon-rich, but the form(s) of the carbon is unknown. Near-infrared spectra of the low-albedo hemisphere of Saturn's satellite Iapetus were obtained with the Visible-Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft at the fly-by of that satellite of 31 December 2004, yielding a maximum spatial resolution on the satellite's surface of ∼65 km. The spectral region 3-3.6 μm reveals a broad absorption band, centered at 3.29 μm, and concentrated in a region comprising about 15% of the low-albedo surface area. This is identified as the CH stretching mode vibration in polycyclic aromatic hydrocarbon (PAH) molecules. Two weaker bands attributed to CH₂ stretching modes in aliphatic hydrocarbons are found in association with the aromatic band. The bands most likely arise from aromatic and aliphatic units in complex macromolecular carbonaceous material with a kerogen- or coal-like structure, similar to that in carbonaceous meteorites. VIMS spectra of Phoebe, encountered by Cassini on 11 June 2004, also show the aromatic hydrocarbon band, although somewhat weaker than on Iapetus. The origin of the PAH molecular material on these two satellites is unknown, but PAHs are found in carbonaceous meteorites, cometary dust particles, circumstellar dust, and interstellar dust.     

Sublimation’s impact on temporal change of albedo dichotomy on Iapetus

Iapetus, one of the saturnian moons, has an extreme albedo contrast between the leading and trailing hemispheres. The origin of this albedo dichotomy has led to several hypotheses, however it remains controversial. To clarify the origin of the dichotomy, the key approach is to investigate the detailed distribution of the dark material. Recent studies of impact craters and surface temperature from Cassini spacecraft data implied that sublimation of H₂O ice can occur on Iapetus’ surface. This ice sublimation can change the albedo distribution on the moon with time.In this study, we evaluate the effect of ice sublimation and simulate the temporal change of surface albedo. We assume the dark material and the bright ice on the surface to be uniformly mixed with a certain volume fraction, and the initial albedo distribution to incorporate the dark material deposits on the surface. That is, the albedo at the apex is lowest and concentrically increases in a sinusoidal pattern. This situation simulates that dark materials existed around the Iapetus’ orbit billions of years ago, and the synchronously rotating Iapetus swept the material and then deposited it on its surface. The evolution of the surface albedo during 4.0 Gyr is simulated by estimating the surface temperature from the insolation energy on Iapetus including the effect of Saturn’s eccentricity and Iapetus’ obliquity precession, and evaluating the sublimation rate of H₂O ice from the Iapetus’ surface.As a result, we found that the distribution of the surface albedo changed dramatically after 4.0 Gyr of evolution. The sublimation has three important effects on the resultant surface albedo. First, the albedo in the leading hemisphere has significantly decreased to approach the minimum value. Second, the albedo distribution has been elongated along the equator. Third, the edge of the low albedo region has become clear. Considering the effect of ice sublimation, the current albedo distribution can be reconstructed from the sinusoidal albedo distribution, suggesting the apex-antapex cratering asymmetry as a candidate for the origin of the albedo dichotomy. From the model analysis, we obtained an important aspect that the depth of the turn-over layer where the darkening process proceeded for 4 Gyr should be an order of 10 cm, which is consistent with evaluation from the Cassini radar observations.     

Tectonics on Iapetus: Despinning, respinning, or something completely different?

Saturn’s moon Iapetus is unique in that it has apparently despun while retaining a substantial equatorial bulge. Stresses arising from such a non-hydrostatic shape should in principle cause surface deformation (tectonics). As part of a search for such a tectonic signature, lineaments (linear surface features) on Iapetus were mapped on both its bright and dark hemispheres. Lineament orientations were then compared to model stress patterns predicted for spin-down from a rotation period of 16.5 h (or less) to its present synchronous period, and for a range of lithospheric thicknesses. Many lineaments are straight segments of crater rimwalls, which may be faults or joints reactivated during complex crater collapse. Most striking are several large troughs on the bright, trailing hemisphere. These troughs appear to be extensional and are distinctive on that hemisphere, because the interior floors and walls of the troughs contain dark material. Globally, no specific evidence of strike slip or thrust offsets are seen, but this could be due to the age and degraded nature of any such features. We find that observed lineament orientations do not correlate with predicted patterns due to despinning on either hemisphere (the equatorial ridge was specifically excluded from this analysis, and is considered separately). Modest evidence for preferred orientations ±40° from north could be construed as consistent with respinning, which is not necessarily far-fetched. Assuming the rigidity of unfractured ice, predicted maximum lithospheric differential stresses from despinning range from ∼1 MPa to ∼160 MPa for the elastic spheroid and thin lithosphere limits, respectively (although it is only for thicker elastic lithospheres that we expect a nonhydrostatic state to be maintained over geologic time against lithospheric failure). The tectonic signature of despinning may have been obscured over time because the surface of Iapetus is very ancient, Iapetus’ thick lithosphere may have inhibited the full tectonic expression of despinning, or both. Several prominent lineaments strike E–W, and are thus parallel to the equatorial ridge (though not physically close to it), but a tectonic or volcanic origin for the ridge is highly problematic.     

Solar phase curves and phase integrals for the leading and trailing hemispheres of Iapetus from the Cassini Visual Infrared Mapping Spectrometer

We performed photometry of Cassini Visual Infrared Mapping Spectrometer observations of Iapetus to produce the first phase integrals calculated directly from solar phase curves of Iapetus for the leading hemisphere and to estimate the phase integrals for the trailing hemisphere. We also explored the phase integral dependence on wavelength and geometric albedo. The extreme dichotomy of the brightness of the leading and trailing sides of Iapetus is reflected in their phase integrals. Our phase integrals, which are lower than the results of Morrison et al. (Morrison, D., Jones, T.J., Cruikshank, D.P., Murphy, R.E. [1975]. Icarus 24, 157-171) and Squyres et al. (Squyres, S.W., Buratti, B.J., Veverka, J., Sagan, C. [1984]. Icarus 59, 426-435), have profound implications on the energy balance and volatile transport on this icy satellite.     

Grooves on small saturnian satellites and other objects: Characteristics and significance

High-resolution images from the Cassini Imaging Science Subsystem (ISS) show parallel sets of grooves on Epimetheus and Pandora. Grooves have previously been observed on other satellites and asteroids, including Phobos, Gaspra, Ida, Eros, and minor occurrences on Phoebe. Sets of parallel grooves are so far observed only on satellites known or likely to be subject to significant tidal stresses, such as forced librations. Grooves on asteroids and on satellites not subject to significant forced librations occur in more globally disorganized patterns that may reflect impacts, varying internal structures, or even thermal stresses. The patterns and individual morphologies of grooves on the tidally-affected satellites suggest fracturing in weak materials due to tidal stresses and forced librations.     

Iapetus surface variability revealed from statistical clustering of a VIMS mosaic: The distribution of CO2

We present a detailed study of an Iapetus mosaic of VIMS data with high spatial resolution (0.5 × 0.5° or ∼6.4 km/pixel). The spectra were taken in August 2007 and provide the highest VIMS spatial resolution data for this object during Cassini’s primary mission. We analyze this set of data using a statistical clustering approach to reduce the analysis of a large number of data (∼10⁴ spectra from 0.35 to 5.10 μm) to the study of seven representative groups accounting for 99.6% of the surface covered by the original sample. We analyze the spectral absorption bands in the spectra of the different clusters indicative of different composition over the observed surface. We find coherence between the distribution of the clusters and the geographical features on the surface. We give special attention to the study of the water ice and CO₂ bands. We find that CO₂ is widespread over the entire surface being studied, including the bright and dark areas on Iapetus’ surface, and is probably trapped at the molecular level with other materials. The strength of the CO₂ band in the areas where both, H₂O- and carbon-bearing materials exist, gives support to the hypothesis that this volatile is formed on the surface of Iapetus as a product of irradiation of these two components. Finally, we also compare the Iapetus CO₂ with that on other satellites confirming, that there are evident differences on the center, depth and width of the band on Iapetus and Phoebe, where CO₂ has been suggested to be endogenous.     

The albedo dichotomy of Iapetus measured at UV wavelengths

The dramatic hemispheric dichotomy in albedo displayed by Saturn's moon Iapetus has intrigued astronomers for centuries. Here we report on far-ultraviolet observations of Iapetus' bright and dark terrains from Cassini. We compare the reflectance spectra of Iapetus's dark terrain, Hyperion and Phoebe and find that both Phoebe and Hyperion are richer in water ice than Iapetus' dark terrain. Spectra of the lowest latitudes of the dark terrain display the diagnostic water ice absorption feature; water ice amounts increase within the dark material away from the apex (at 90° W longitude, the center of the dark leading hemisphere), consistent with thermal segregation of water ice. The water ice in the darkest, warmest low latitude regions is not expected to be stable and may be a sign of ongoing or recent emplacement of the dark material from an exogenic source.     

Spectrophotometry of the small satellites of Saturn and their relationship to Iapetus, Phoebe, and Hyperion

We have obtained broadband spectrophotometric observations of four of the recently discovered small satellites of Saturn (Gladman et al., 2001, Nature 412, 163-166). The new data enable an understanding of the provenance, composition, and interrelationships among these satellites and the other satellites of Saturn, particularly Iapetus, Phoebe, and Hyperion. Temporal coverage of one satellite (S21 Tarvos) was sufficient to determine a partial rotational lightcurve. Our major findings include: (1) the satellites are red and similar in color, comparable to D-type asteroids, some KBOs, Iapetus, and Hyperion; (2) none of the satellites, including those from the “Phoebe Group” has any spectrophotometric relationship to Phoebe; and (3) S21 Tarvos exhibits a rotational lightcurve, although the data are not well-constrained and more observations are required to fit a definitive period. Dust created by meteoritic impacts and ejected from these satellites and additional undiscovered ones may be the source of the exogenous material deposited on the low-albedo side of Iapetus. Recent work which states that the small irregular satellites of Saturn have impacted Phoebe at least 6-7 times in the age of the Solar System (Nesvorny et al., 2003, Astron. J. 126, 398-429), suggests that such collisions may have propelled additional material from both Phoebe and the small irregular satellites toward Iapetus. The accretion of material from outer retrograde satellites may be a process that also occurs on Callisto and the uranian satellites.     

The stability and transport of carbon dioxide on Iapetus

Carbon dioxide has been detected associated with Iapetus' dark material by the Cassini spacecraft. This CO₂ may be primordial and/or resulting from ongoing production by photolysis of water-ice in the presence of carbonaceous material [Allamandola, L.J., Sandford, S.A., Valero, G.J., 1988. Icarus 76, 225-252]. Although any primordial CO₂ would likely be complexed with the dark material and thus stable against thermal transport to Iapetus' poles [Buratti, B.J., and 28 colleagues, 2005. Astrophys. J. 622, L149-L152], active production of CO₂ would result in some fraction of the CO₂ being mobile enough to allow the accumulation of CO₂ at Iapetus' poles. We develop a computer model to simulate ballistic transport of CO₂ ice on Iapetus, accounting for Iapetus' gravitational binding energy and polar cold traps. We find that the residence time of CO₂ ice outside the polar regions is very short; a sheet of CO₂ ice near the equator of Iapetus decreases in thickness at a rate of 50 mm year⁻¹. The sublimated CO₂ will ballistically move around Iapetus until it reaches the polar cold traps where it can be sequestered for up to 15 years. If the total surface inventory of CO₂ exceeds 3×¹⁰⁷ kg, the polar ice cap will be permanent. While CO₂ is moving around the surface, a small percentage will eventually reach escape velocity and be lost from the system. As such, a seasonal polar cap is lost at rate of 12% every solar orbit as the CO₂ moves between the two polar cold traps.     

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.     

Production and detection of carbon dioxide on Iapetus

Cassini VIMS detected carbon dioxide on the surface of Iapetus during its insertion orbit. We evaluated the CO₂ distribution on Iapetus and determined that it is concentrated almost exclusively on Iapetus’ dark material. VIMS spectra show a 4.27-μm feature with an absorption depth of 24%, which, if it were in the form of free ice, requires a layer 31 nm thick. Extrapolating for all dark material on Iapetus, the total observable CO₂ would be 2.3 × 10⁸ kg.Previous studies note that free CO₂ is unstable at 10 AU over geologic timescales. Carbon dioxide could, however, be stable if trapped or complexed, such as in inclusions or clathrates. While complexed CO₂ has a lower thermal volatility, loss due to photodissociation by UV radiation and gravitational escape would occur at a rate of 2.6 × 10⁷ kg year⁻¹. Thus, Iapetus’ entire inventory of surface CO₂ could be lost within a few decades.The high loss/destruction rate of CO₂ requires an active source. We conducted experiments that generated CO₂ by UV radiation of simulated icy regolith under Iapetus-like conditions. The simulated regolith was created by flash-freezing degassed water, crushing it into sub-millimeter sized particles, and then mixing it with isotopically labeled amorphous carbon (¹³C) dust. These samples were placed in a vacuum chamber and cooled to temperatures between 50 K and 160 K. The samples were irradiated with UV light, and the products were measured using a mass spectrometer, from which we measured ¹³CO₂ production at a rate of 2.0 × 10¹² mol s⁻¹. Extrapolating to Iapetus and adjusting for the solar UV intensity and Iapetus’ surface area, we calculated that CO₂ production for the entire surface would be 1.1 × 10⁷ kg year⁻¹, which is only a factor of two less than the loss rate. As such, UV photochemical generation of CO₂ is a plausible source of the detected CO₂.     

Dynamics of rotationally fissioned asteroids: Source of observed small asteroid systems

We present a model of near-Earth asteroid (NEA) rotational fission and ensuing dynamics that describes the creation of synchronous binaries and all other observed NEA systems including: doubly synchronous binaries, high-e binaries, ternary systems, and contact binaries. Our model only presupposes the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect, “rubble pile” asteroid geophysics, and gravitational interactions. The YORP effect torques a “rubble pile” asteroid until the asteroid reaches its fission spin limit and the components enter orbit about each other (Scheeres, D.J. [2007]. Icarus 189, 370-385). Non-spherical gravitational potentials couple the spin states to the orbit state and chaotically drive the system towards the observed asteroid classes along two evolutionary tracks primarily distinguished by mass ratio. Related to this is a new binary process termed secondary fission - the secondary asteroid of the binary system is rotationally accelerated via gravitational torques until it fissions, thus creating a chaotic ternary system. The initially chaotic binary can be stabilized to create a synchronous binary by components of the fissioned secondary asteroid impacting the primary asteroid, solar gravitational perturbations, and mutual body tides. These results emphasize the importance of the initial component size distribution and configuration within the parent asteroid. NEAs may go through multiple binary cycles and many YORP-induced rotational fissions during their approximately 10 Myr lifetime in the inner Solar System. Rotational fission and the ensuing dynamics are responsible for all NEA systems including the most commonly observed synchronous binaries.     

Exchange of ejecta between Telesto and Calypso: Tadpoles, horseshoes, and passing orbits

We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn’s major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent’s escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction of launch.The slowest ejecta follow suborbital trajectories and re-impact their source moon in less than one day. Slightly faster debris barely escape their parent’s Hill sphere and are confined to tadpole orbits, librating about Tethys’ triangular Lagrange points L₄ (leading, near Telesto) or L₅ (trailing, near Calypso) with nearly the same orbital semi-major axis as Tethys, Telesto, and Calypso. These ejecta too eventually re-impact their source moon, but with a median lifetime of a few dozen years. Those which re-impact within the first 10 years or so have lifetimes near integer multiples of 348.6 days (half the tadpole period).Still faster debris with azimuthal velocity components ?10 m/s enter horseshoe orbits which enclose both L₄ and L₅ as well as L₃, but which avoid Tethys and its Hill sphere. These ejecta impact either Telesto or Calypso at comparable rates, with median lifetimes of several thousand years. However, they cannot reach Tethys itself; only the fastest ejecta, with azimuthal velocities ?40 m/s, achieve “passing orbits” which are able to encounter Tethys. Tethys accretes most of these ejecta within several years, but some 1% of them are scattered either inward to hit Enceladus or outward to strike Dione, over timescales on the order of a few hundred years.