Europa’s disk-resolved ultraviolet spectra: Relationships with plasma flux and surface terrains

The full set of high-resolution observations from the Galileo Ultraviolet Spectrometer (UVS) is analyzed to look for spectral trends across the surface of Europa. We provide the first disk-resolved map of the 280 nm SO₂ absorption feature and investigate its relationship with sulfur and electron flux distributions as well as with surface features and relative surface ages. Our results have implications for exogenic and endogenic sources. The large-scale pattern in SO₂ absorption band depth is again shown to be similar to the pattern of sulfur ion implantation, but with strong variations in band depth based on terrain. In particular, the young chaos units show stronger SO₂ absorption bands than expected from the average pattern of sulfur ion flux, suggesting a local source of SO₂ in those regions, or diapiric heating that leads to a sulfur-rich lag deposit.While the SO₂ absorption feature is confined to the trailing hemisphere, the near UV albedo (300-310 nm) has a global pattern with a minimum at the center of the trailing hemisphere and a maximum at the center of the leading hemisphere. The global nature of the albedo pattern is suggestive of an exogenic source, and several possibilities are discussed. Like the SO₂ absorption, the near UV albedo also has local variations that depend on terrain type and age.     

Far ultraviolet spectral properties of Saturn’s rings from Cassini UVIS

Spectra taken by the Cassini Ultraviolet Imaging Spectrograph (UVIS) of Saturn’s C ring, B ring, Cassini Division, and A ring have been analyzed in order to characterize ring particle surface properties and water ice abundance in the rings. UVIS spectra sense the outer few microns of the ring particles. Spectra of the normalized reflectance (I/F) in all four regions show a characteristic water ice absorption feature near 165 nm. Our analysis shows that the fractional abundance of surface water ice is largest in the outer B ring and decreases by over a factor of 2 across the inner C ring. We calculate the mean path length of UV photons through icy ring particle regolith and the scattering asymmetry parameter using a Hapke reflectance model and a Shkuratov reflectance model to match the location of the water ice absorption edge in the data. Both models give similar retrieved values of the photon mean length, however the retrieved asymmetry (g) values are different. The photon mean path lengths are nearly uniform across the B and A rings. Shortward of 165 nm the rings exhibit a slope that turns up towards shorter wavelengths, while the UV slope of 180/150 nm (reflectance outside the water absorption ratioed to that inside the absorption band) tracks I/F with maxima in the outer B ring and in the central A ring. Retrieved values of the scattering asymmetry parameter show the regolith grains to be highly backscattering in the FUV spectral regime.     

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.     

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.     

Monte Carlo simulations of the water vapor plumes on Enceladus

Monte Carlo simulations are used to model the July 14, 2005 UVIS stellar occultation observations of the water vapor plumes on Enceladus. These simulations indicate that the observations can be best fit if the water molecules ejected along the Tiger Stripes in the South Polar region of Enceladus have a vertical surface velocity of 300-500 m/s at the surface. The high surface velocity suggests that the plumes on Enceladus originate from some depth beneath the surface. The total escape rate of water molecules is 4-6×¹⁰²⁷ s⁻¹, or 120-180 kg/s, consistent with previous works, and more than 100 times the estimated mass escape rate for ice particles. The average deposition rate in the South Polar region is on the order of 10¹¹ cm⁻² s⁻¹, yielding a resurfacing rate as high as 3×¹⁰⁻⁴ cm/yr. The globally averaged deposition rate of water molecules is about one order of magnitude lower.     

Enceladus: The likely dominant nitrogen source in Saturn's magnetosphere

The spatial distribution of N⁺ in Saturn's magnetosphere obtained from Cassini Plasma Spectrometer (CAPS) data can be used to determine the spatial distribution and relative importance of the nitrogen sources for Saturn's magnetosphere. We first summarize CAPS data from 15 orbits showing the spatial and energy distribution of the nitrogen component of the plasma. This analysis re-enforces our earlier discovery [Smith, H.T., Shappirio, M., Sittler, E.C., Reisenfeld, D., Johnson, R.E., Baragiola, R.A., Crary, F.J., McComas, D.J., Young, D.T., 2005. Geophys. Res. Lett. 32 (14). L14S03] that Enceladus is likely the dominant nitrogen source for Saturn's inner magnetosphere. We also find a sharp enhancement in the nitrogen ion to water ion ratio near the orbit of Enceladus which, we show, is consistent with the presence of a narrow Enceladus torus as described in [Johnson, R.E., Liu, M., Sittler Jr., E.C., 2005. Geophys. Res. Lett. 32. L24201]. The CAPS data and the model described below indicate that N⁺ ions are a significant fraction of the plasma in this narrow torus. We then simulated the combined Enceladus and Titan nitrogen sources using the CAPS data as a constraint. This simulation is an extension of the model we employed earlier to describe the neutral tori produced by the loss of nitrogen from Titan [Smith, H.T., Johnson, R.E., Shematovich, V.I., 2004. Geophys. Res. Lett. 31 (16). L16804]. We show that Enceladus is the principal nitrogen source in the inner magnetosphere but Titan might account for a fraction of the observed nitrogen ions at the largest distances discussed. We also show that the CAPS data is consistent with Enceladus being a molecular nitrogen source with a nitrogen to water ratio roughly consistent with INMS [Waite, J.H., and 13 colleagues, 2006. Science 311 (5766), 1419-1422], but out-gassing of other nitrogen-containing species, such as ammonia, cannot be ruled out.     

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.     

Near-infrared spectra of the leading and trailing hemispheres of Enceladus

We present individual spectra 0.8-2.5 μm of the leading and trailing hemispheres of Enceladus obtained with the CorMASS spectrograph on the 1.8 m Vatican Advanced Technology Telescope (VATT) at the Mount Graham International Observatory. While the absorption bands of water ice dominate the spectrum of both hemispheres, most of these bands are stronger on the leading hemisphere than the trailing hemisphere. In addition, longward of 1 μm, the continuum slope is greater on the leading hemisphere than the trailing hemisphere. These differences could be produced by the presence of particles on the trailing side that are smaller and/or microstructurally more complex than those on the leading side, consistent with the preferential erosion or structural degradation of regolith particle grains on the trailing side by magnetospheric sweeping. We also explore compositional differences between the two hemispheres by applying Hapke spectrophotometric mixture models to the spectra whose components include water ice and ammonia hydrate (1% NH₃⋅H₂O). We find that spectral models which include as much as 25% by weight ammonia hydrate intimately mixed with water ice and covering 80% of the illuminated area of the satellite fit the observed spectrum of both the leading and trailing hemispheres. Areal (checkerboard) mixing models of ammonia hydrate and water ice fit the leading hemisphere with 15% of the surface comprised of ammonia hydrate and the trailing hemisphere with 10% ammonia hydrate. Therefore, while these spectral data do not contain an unambiguous detection of ammonia hydrate on Enceladus, our spectral models do not preclude the presence of a modest amount of 1% NH₃⋅H₂O on both hemispheres. We examine spectral differences and similarities between both hemispheres and the tenuous E ring within which Enceladus orbits. The spectral resolution (R=λ/Δλ) of these CorMASS data (R∼300) is comparable to but nevertheless higher than that of the Visual-Infrared Mapping Spectrometer (VIMS) (R=225) onboard the Cassini spacecraft.     

The oxidation state of hydrothermal systems on early Enceladus

The discovery of CO₂, CH₄, and N₂ in a plume at Enceladus provides useful clues about the chemistry and evolution of this moon of Saturn. Here, we use chemical equilibrium and kinetic calculations to estimate the oxidation state of hydrothermal systems on early Enceladus, with the assumption that the plume's composition was inherited from early hydrothermal fluids. Chemical equilibrium calculations are performed using the CO₂/CH₄ ratio in the plume, and kinetic calculations are conducted using equations from fluid dynamics and chemical kinetics. Our results suggest that chemical equilibrium between CO₂ and CH₄ would have been reachable at temperatures above ∼200 °C in hydrothermal systems. The oxidation state of the hydrothermal systems would have been close to the pyrrhotite-pyrite-magnetite (PPM) or fayalite-magnetite-quartz (FMQ) redox buffer (i.e., terrestrial-like) if the plume's CO₂ and CH₄ equilibrated in hydrothermal systems long ago. As for minerals, we suggest that iron metal would have been oxidized to magnetite by the escape of H₂ from the early satellite. Our calculations also indicate that, assuming CO₂ and CH₄ reached chemical equilibrium, magnetite would not have been oxidized to hematite in hydrothermal systems, perhaps due to insufficient H₂ escape. It is shown that, if Enceladus accreted as much NH₃ as comets contain, the presence of N₂ and deficiency of NH₃ in the plume can be understood in the context of chemical equilibrium in the C-N-O-H system. We conclude by proposing an evolutionary hypothesis in which the fairly oxidized nature of the plume can be explained by a brief episode of oxidation caused by short-lived radioactivity. These suggestions can be rigorously tested by acquiring gravity and isotopic data in the future.     

Ultraviolet observations of Phoebe from the Cassini UVIS

Observations of Saturn's distant moon Phoebe were made at far-ultraviolet (FUV) (1100-1900 Å) and extreme-ultraviolet (EUV) (600-1100 Å) wavelengths by the Cassini Ultraviolet Imaging Spectrograph (UVIS) during the Cassini spacecraft flyby on June 11, 2004. These are the first UV spectra of Phoebe and the first detection of water ice on a Solar System surface using FUV wavelengths. The characteristics of water ice in the FUV are presented, and Hapke models are used to interpret the spectra in terms of composition and grain size; the use of both areal and intimate mixing models is explored. Non-ice species used in these models include carbon, ice tholin, Triton tholin, poly-HCN and kerogen. Satisfactory disk-integrated fits are obtained for intimate mixtures of ∼10% H₂O plus a non-ice species. Spatially resolved regions of higher (∼20%) and lower (∼5%) H₂O ice concentrations are also detected. Phoebe does not display any evidence of volatile activity. Upper limits on atomic oxygen and carbon are 5×¹⁰¹¹ and 2×¹⁰¹² atoms/cm², respectively, for solar photon scattering. The UVIS detection of water ice on Phoebe, and the ice amounts detected, are consistent with IR measurements and contribute to the evidence for a Phoebe origin in the outer Solar System rather than in the main asteroid belt.     

A model for the temperature-dependence of tidal dissipation in convective plumes on icy satellites: Implications for Europa and Enceladus

To explain the formation of surface features on Europa, Enceladus, and other satellites, many authors have postulated the spatial localization of tidal heating within convective plumes. However, the concept that enhanced tidal heating can occur within a convective plume has not been rigorously tested. Most models of this phenomenon adopt a tidal heating with a temperature-dependence derived for an incompressible, homogeneous (zero-dimensional) Maxwell material, but it is unclear whether this formulation is relevant to the heterogeneous situation of a warm plume surrounded by cold ice. To determine whether concentrated dissipation can occur in convective plumes, we develop a two-dimensional model to compute the volumetric dissipation rate for an idealized, vertically oriented, isolated convective plume obeying a Maxwellian viscoelastic compressible rheology. We apply the model to the Europa and Enceladus ice shells, and we investigate the consequences for partial melting and resurfacing processes on these bodies. We find that the tidal heating is strongly temperature dependent in a convective ice plume and could produce elevated temperatures and local partial melting in the ice shells of Europa and Enceladus. Our calculation provides the first quantitative verification of the hypothesis by Sotin et al. [Sotin, C., Head, J.W., Tobie, G., 2002. Geophys. Res. Lett. 29. 74-1] and others that the tidal dissipation rate is a strong function of temperature inside a convective plume. On Europa, such localized heating could help allow the formation of domes and chaos terrains by convection. On Enceladus, localized tidal heating in a thermal plume could explain the concentrated activity at the south pole and its associated heat transport of 2-7 GW.     

Arecibo radar observations of Rhea, Dione, Tethys, and Enceladus

We have measured the bulk radar reflectance properties of the mid-size saturnian satellites Rhea, Dione, Tethys, and Enceladus with the Arecibo Observatory's 13 cm wavelength radar system during the 2004 through 2007 oppositions of the Saturn system. Comparing to the better studied icy Galilean satellites, we find that the total reflectivities of Rhea and Tethys are most similar to Ganymede while Dione is most similar to Callisto. Enceladus' reflectivity falls between those of Ganymede and Europa. The mean circular polarization ratios of the saturnian satellites range from ∼0.8 to 1.2, and are on average lower than those of the icy Galilean satellites at this wavelength although still larger than expected for single reflections off the surface. The ratio for the trailing hemisphere of Enceladus may be the exception with a value ?0.56. The 13 cm wavelength radar albedos and polarization ratios may be systematically lower than similar results from the Cassini orbiter's RADAR instrument at 2.2 cm wavelength [Ostro, S.J., and 19 colleagues, 2006. Icarus 183, 479-490]. Overall, these reflectivities and polarization properties, together with the shapes of the echo spectra, suggest subsurface multiple scattering to be the dominant reflection mechanism although operating less efficiently than on the large icy moons of Jupiter. All these saturnian moons and icy jovian moons are atmosphere-less, low temperature water ice surfaces, and any differences in radar properties may be indicative of differences in composition or the effects of various processes that modify the regolith structure. The degree of variation in radar properties with wavelength on each satellite may constrain the thickness and efficiency of the scattering layer.     

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.     

Regolith depth growth on an icy body orbiting Saturn and evolution of bidirectional reflectance due to surface composition changes

Using a Markov chain model, we consider the regolith growth on a small body in orbit around Saturn, subject to meteoritic bombardment, and assuming all impact ejecta are re-collected. We calculate the growth of regolith and the fractional pollution, assuming an initial pure ice body and amorphous carbon as a pollutant. We extend the meteorite flux of Cuzzi and Estrada (Cuzzi, J., Estrada, P. [1998]. Icarus 132, 1-35) to larger sizes to consider the effect of disruption of the moonlet on other moonlets in the ensemble. This is a relatively small effect, completely negligible for moonlets of 1 m radius. For the given impact model, fractional pollution reaches 22% for 1 m bodies, but only 3% for 10 m bodies, 1.7% for 20 m bodies, and 1% for 30 m bodies after 4 byr. By considering an ensemble of moonlets, which have identical cross-sections for releasing and capturing ejecta, this analysis can be extended to a model of particles in Saturn’s rings, where the calculated spectra can be compared to observed ring spectra. The measured spectral reflectance of Saturn’s rings from Cassini observations therefore constrains the size and age of the ring particles. The comparison between 1 m, 10 m, 20 m, and 30 m particles confirms that for larger ring mass, the current rings would be less polluted; for the largest particles, we expect negligible changes in the UV spectrum after 4 byr of meteoritic bombardment. We consider two end members for mixing of the meteoritic material: areal and intimate. Given the uncertainties in the actual mixing of the meteoritic infall and in its composition (as a worst case, we assume the meteoritic material is 100% amorphous carbon, intimately mixed) initially pure ice 30 m ring particles would darken after 4 byr of exposure by 15%.     

Photometric and spectral analysis of the distribution of crystalline and amorphous ices on Enceladus as seen by Cassini

Photometric and spectral analysis of data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) has yielded significant results regarding the properties and composition of the surface of Saturn's satellite Enceladus. We have obtained spectral cubes of this satellite, containing both spatial and spectral information, with a wavelength distribution in the infrared far more extensive than from any previous observations and at much higher spatial resolution. Using a composite mosaic of the satellite, we map the distribution of crystalline and amorphous ices on the surface of Enceladus according to a “crystallinity factor” and also the depth of the temperature- and structure-dependent 1.65 micron water-ice band. These maps show the surface of Enceladus to be mostly crystalline, with a higher degree of crystallinity at the “tiger-stripe” cracks and a larger amorphous signature between these stripes. These results suggest recent geological activity at the “tiger stripe” cracks and an intriguing atmospheric environment over the south pole where amorphous ice is produced either through intense radiative bombardment, flash-freezing of cryovolcanic liquid, or rapid condensation of water vapor particles on icy microspherules or on the surface of Enceladus.     

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

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

Tidal evolution of Mimas, Enceladus, and Dione

The tidal evolution through several resonances involving Mimas, Enceladus, and/or Dione is studied numerically with an averaged resonance model. We find that, in the Enceladus-Dione 2:1 e-Enceladus type resonance, Enceladus evolves chaotically in the future for some values of k₂/Q. Past evolution of the system is marked by temporary capture into the Enceladus-Dione 4:2 ee^′-mixed resonance. We find that the free libration of the Enceladus-Dione 2:1 e-Enceladus resonance angle of 1.5° can be explained by a recent passage of the system through a secondary resonance. In simulations with passage through the secondary resonance, the system enters the current Enceladus-Dione resonance close to tidal equilibrium and thus the equilibrium value of tidal heating of 1.1(18,000/QS) GW applies. We find that the current anomalously large eccentricity of Mimas can be explained by passage through several past resonances. In all cases, escape from the resonance occurs by unstable growth of the libration angle, sometimes with the help of a secondary resonance. Explanation of the current eccentricity of Mimas by evolution through these resonances implies that the Q of Saturn is below 100,000. Though the eccentricity of Enceladus can be excited to moderate values by capture in the Mimas-Enceladus 3:2 e-Enceladus resonance, the libration amplitude damps and the system does not escape. Thus past occupancy of this resonance and consequent tidal heating of Enceladus is excluded. The construction of a coherent history places constraints on the allowed values of k₂/Q for the satellites.     

The opposition surge of Enceladus: HST observations 338-1022 nm

Hubble Space Telescope (HST) Wide-Field Planetary Camera (WFPC2) observations at phase angles in the range α=0.26°-6.4° obtained at every opposition and near quadrature between October 1996 and December 2002 reveal the opposition effect of Enceladus. We present a photometric analysis of nearly 200 images obtained through the five broadband UVBRI filters (F336W, F439W, F555W, F675W, and F814W) and the F785LP and F1042M filters from which we generate mutually consistent solar and rotational phase curves. Our solar phase curves reveal a dramatic, sharp increase in the albedo (from 0.11 mag in the F675W filter to 0.17 mag in the F785LP filter) as phase angles decrease from 2° to 0.26°. A slight opposition effect is evident in data from the F1042M filter (λ_eff=1022 nm); however, the smallest phase angle currently available for observations from this filter is α=0.58°. With the addition of data from the F255W filter we demonstrate the wavelength dependence of the albedo of the trailing hemisphere from 275 to 1022 nm. Our rotation curves show that the trailing hemisphere is ∼0.06 mag brighter than the leading when observed at wavelengths between 338 and 868 nm and 0.11 mag brighter than the leading at 1022 nm. We have supplemented the phase curve from the F439W filter (λ_eff=434 nm) with Voyager clear filter (λ_eff=480 nm) observations made at larger phase angles (α=13°-43°) to produce a phase curve with the most extensive phase angle coverage possible to date. This newly expanded range of phase angles enhances the ability of the Hapke photometric model (Hapke B., 2002, Icarus 157, 523-534) to relate physical characteristics of the surface of Enceladus to the manner in which incident light is reflected from it. We present Hapke 2002 model fits to solar phase curves from each UVBRI filter as well as from the F785LP and F1042M filters. Geometric albedos derived from these model fits range from p=0.92±0.01 at 1022 nm to p=1.41±0.03 at 549 nm, necessitating an increase of about 20% from previously derived values. Our Hapke fits demonstrate that the opposition surge of Enceladus is best described by a model which combines both moderate shadow-hiding and narrow coherent backscattering components.     

Tidal heating in Enceladus

The heating in Enceladus in an equilibrium resonant configuration with other saturnian satellites can be estimated independently of the physical properties of Enceladus. We find that equilibrium tidal heating cannot account for the heat that is observed to be coming from Enceladus. Equilibrium heating in possible past resonances likewise cannot explain prior resurfacing events.     

Episodic volcanism on Enceladus: Application of the Ojakangas-Stevenson model

The main equations in the paper “Episodic volcanism of tidally heated satellites with application to Io” by Ojakangas and Stevenson [Icarus 66, 341-358] are presented; numerical integration of these equations confirms the results of Ojakangas and Stevenson [Icarus 66, 341-358] for Io. Application to Enceladus is considered. It is shown that Enceladus does not oscillate about the tidal equilibrium in this model by both new nonlinear stability analysis and numerical integration of the model equations.