Metasomatic clathrate xenoliths as a possible source for the south polar plumes of Enceladus

The composition and abundance of volatile gases observed in the jets emanating from fissures near the south pole of Saturn's moon Enceladus are strongly indicative of outgassing from clathrate hydrates which formed as a result of hydrothermal activity rather than nebula condensation. I suggest that fluids must be able to permeate the ice shell, extensively metasomatising the mantle by emplacement of clathrates along fractures and grain boundaries, which subsequently are entrained in rising cryomagmas as xenoliths. These are carried upwards to the point where they dissociate, releasing their gas load into the magma and promoting the vigorous ice fountaining observed—a direct analogue of terrestrial basaltic fire fountains caused by volatile exsolution. This clathrate xenolith model can explain the measured volatile abundances, eruption velocities, the ice to vapour ratio in the plumes, and the vent temperatures.     

Tidal heating and the long-term stability of a subsurface ocean on Enceladus

Tidal dissipation has been suggested as the heat source for the south polar thermal anomaly on Enceladus. We find that under present-day conditions and assuming Maxwellian behavior, tidal dissipation is negligible in the silicate core. Dissipation may be significant in the ice shell if the shell is decoupled from the silicate core by a subsurface ocean. We have run a series of self-consistent convection and conduction models in 2D axisymmetric and 3D spherical geometry in which we include the spatially-variable tidal heat production. We find that in all cases, the shell removes more heat from the interior than can be produced in the core by radioactive decay, resulting in cooling of the interior and the freezing of any ocean. Under likely conditions, a 40-km thick ocean made of pure water would freeze solid on a ∼30 Ma timescale. An ocean containing other chemical components will have a lower freezing point, but even a water-ammonia eutectic composition will only prolong the freezing, not prevent it. If the eccentricity of Enceladus were higher (e?0.015) in the past, the increased dissipation in the ice shell may have been sufficient to maintain a liquid layer. We cannot therefore rule out the presence of a transient ocean, as a relic of an earlier era of greater heating. If the eccentricity is periodically pumped up, the ocean may have thickened and thinned on a similar timescale as the orbital evolution, provided the ocean never froze completely. We conclude that the current heat flux of Enceladus and any possible subsurface ocean is not in steady-state, and is the remnant of an epoch of higher eccentricity and tidal dissipation.     

Saturn in hot water: Viscous evolution of the Enceladus torus

The detection of outgassing water vapor from Enceladus is one of the great breakthroughs of the Cassini mission. The fate of this water once ionized has been widely studied; here we investigate the effects of purely neutral-neutral interactions within the Enceladus torus. We find that, thanks in part to the polar nature of the water molecule, a cold (∼180 K) neutral torus would undergo rapid viscous heating and spread to the extent of the observed hydroxyl cloud, before plasma effects become important. We investigate the physics behind the spreading of the torus, paying particular attention to the competition between heating and rotational line cooling. A steady-state torus model is constructed, and it is demonstrated that the torus will be observable in the millimeter band with the upcoming Herschel satellite. The relative strength of rotational lines could be used to distinguish between physical models for the neutral cloud.     

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.     

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.     

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.     

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.     

Negative ions in the Enceladus plume

During Cassini’s Enceladus encounter on 12th March 2008, the Cassini Electron Spectrometer, part of the CAPS instrument, detected fluxes of negative ions in the plumes from Enceladus. It is thought that these ions include negatively charged water group cluster ions associated with the plume and forming part of the ‘plume ionosphere’. In this paper we present our observations, argue that these are negative ions, and present preliminary mass identifications. We also suggest mechanisms for production and loss of the ions as constrained by the observations. Due to their short lifetime, we suggest that the ions are produced in or near the water vapour plume, or from the extended source of ice grains in the plume. We suggest that Enceladus now joins the Earth, Comet Halley and Titan as locations in the Solar System where negative ions have been directly observed although the ions observed in each case have distinctly different characteristics.     

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.     

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

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.     

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.     

Enceladus: Present internal structure and differentiation by early and long-term radiogenic heating

Pre-Cassini images of Saturn's small icy moon Enceladus provided the first indication that this satellite has undergone extensive resurfacing and tectonism. Data returned by the Cassini spacecraft have proven Enceladus to be one of the most geologically dynamic bodies in the Solar System. Given that the diameter of Enceladus is only about 500 km, this is a surprising discovery and has made Enceladus an object of much interest. Determining Enceladus' interior structure is key to understanding its current activity. Here we use the mean density of Enceladus (as determined by the Cassini mission to Saturn), Cassini observations of endogenic activity on Enceladus, and numerical simulations of Enceladus' thermal evolution to infer that this satellite is most likely a differentiated body with a large rock-metal core of radius about 150 to 170 km surrounded by a liquid water-ice shell. With a silicate mass fraction of 50% or more, long-term radiogenic heating alone might melt most of the ice in a homogeneous Enceladus after about 500 Myr assuming an initial accretion temperature of about 200 K, no subsolidus convection of the ice, and either a surface temperature higher than at present or a porous, insulating surface. Short-lived radioactivity, e.g., the decay of ²⁶Al, would melt all of the ice and differentiate Enceladus within a few million years of accretion assuming formation of Enceladus at a propitious time prior to the decay of ²⁶Al. Long-lived radioactivity facilitates tidal heating as a source of energy for differentiation by warming the ice in Enceladus so that tidal deformation can become effective. This could explain the difference between Enceladus and Mimas. Mimas, with only a small rock fraction, has experienced relatively little long-term radiogenic heating; it has remained cold and stiff and less susceptible to tidal heating despite its proximity to Saturn and larger eccentricity than Enceladus. It is shown that the shape of Enceladus is not that of a body in hydrostatic equilibrium at its present orbital location and rotation rate. The present shape could be an equilibrium shape corresponding to a time when Enceladus was closer to Saturn and spinning more rapidly, or more likely, to a time when Enceladus was spinning more rapidly at its present orbital location. A liquid water layer on Enceladus is a possible source for the plume in the south polar region assuming the survivability of such a layer to the present. These results could place Enceladus in a category similar to the large satellites of Jupiter, with the core having a rock-metal composition similar to Io, and with a deep overlying ice shell similar to Europa and Ganymede. Indeed, the moment of inertia factor of a differentiated Enceladus, C/MR², could be as small as that of Ganymede, about 0.31.     

Solid tidal friction above a liquid water reservoir as the origin of the south pole hotspot on Enceladus

Earth, Jupiter's moon Io and Saturn's tiny moon Enceladus are the only solid objects in the Solar System to be sufficiently geologically active for their internal heat to be detected by remote sensing. Interestingly, the endogenic activity on Enceladus is only located on a specific region at the south pole, from which jets of water vapor and ice particles have been observed [Spencer, J.R., and 9 colleagues, 2006. Science 311, 1401-1405; Porco, C.C., and 24 colleagues, 2006. Science 311, 1393-1401]. The current polar location of the thermal anomaly can possibly be explained by diapir-induced reorientation of the satellite [Nimmo, F., Pappalardo, R.T., 2006. Nature 441, 614-616], but the thermal anomaly triggering and the heat power required to sustain it over geological timescales remain problematic. Using a three-dimensional viscoelastic numerical model simulating the response of Enceladus to tidal forcing, we explore the effect of a low-viscosity anomaly in the ice shell, localized to the south polar region, on the tidal dissipation patterns. We demonstrate that only interior models with a liquid water layer at depth can explain the observed magnitude of dissipation rate and its particular location at the south pole. Moreover, we show that tidally-induced heat must be generated over a relatively broad region in the southern hemisphere, and it is then transferred toward the south pole where it is episodically released during relatively short resurfacing events. As large tidal dissipation and internal melting cannot be induced in the south polar region in the absence of a pre-existing liquid decoupling layer, we propose that liquid water must have been present in the interior for a very long period of time, and possibly since the satellite formation. Owing to the orbital equilibrium requirement [Meyer, J., Wisdom, J., 2007. Icarus 188, 535-539], sustaining some liquid water at depth is impossible if heat is continuously emitted at a rate of 4-8 GW at the south pole. Based on that requirement, we propose that the current thermal emission is not in equilibrium with the heat production, and that the thermal emission rate is abnormally high at present time. Alternatively, continuous dissipation at a rate of 1-2 GW in the ice shell at the south pole should be sufficient to induce internal melting and it could sustain a layer of liquid water at depth over geologic timescales.     

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.     

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.     

Enceladus' south polar sea

Recent observations of the south pole of Saturn's moon Enceladus by the Cassini spacecraft have revealed an active world, powered by internal heat. In this paper, we propose that localized subsurface melting on Enceladus has produced an internal south polar sea. Evidence for this localized sea comes from the shape of Enceladus, which does not match a differentiated body at its current orbital position. We show that melting induced by the observed heat flow at the south pole produces a large enough pit to match the shape of Enceladus with a differentiated rock and ice interior. Numerical modeling of melting and ice flow shows that the sea produced beneath the south pole is stable against inflow of ductile ice from its surroundings for the duration of the heating. The shape modification due to melting also produces a negative degree-two gravity anomaly, which can reorient the spin axis of Enceladus in order to place the sea at the pole.     

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.     

Thermal convection in ice-I shells of Titan and Enceladus

Cassini-Huygens observations have shown that Titan and Enceladus are geologically active icy satellites. Mitri and Showman [Mitri, G., Showman, A.P., 2005. Icarus 177, 447-460] and McKinnon [McKinnon, W.B., 2006. Icarus 183, 435-450] investigated the dynamics of an ice shell overlying a pure liquid-water ocean and showed that transitions from a conductive state to a convective state have major implications for the surface tectonics. We extend this analysis to the case of ice shells overlying ammonia-water oceans. We explore the thermal state of Titan and Enceladus ice-I shells, and also we investigate the consequences of the ice-I shell conductive-convective switch for the geology. We show that thermal convection can occur, under a range of conditions, in the ice-I shells of Titan and Enceladus. Because the Rayleigh number Ra scales with δ³/ηb, where δ is the thickness of the ice shell and ηb is the viscosity at the base of the ice-I shell, and because ammonia in the liquid layer (if any) strongly depresses the melting temperature of the water ice, Ra equals its critical value for two ice-I shell thicknesses: for relatively thin ice shell with warm, low-viscosity base (Onset I) and for thick ice shell with cold, high-viscosity base (Onset II). At Onset I, for a range of heat fluxes, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. Switches between these states can cause large, rapid changes in the ice-shell thickness. For Enceladus, we demonstrate that an Onset I transition can produce tectonic stress of ∼500 bars and fractures of several tens of km depth. At Onset II, in contrast, we demonstrate that zero equilibrium states exist for a range of heat fluxes. For a mean heat flux within this range, the satellite experiences oscillations in surface heat flux and satellite volume with periods of ∼50-800 Myr even when the interior heat production is constant or monotonically declining in time; these oscillations in the thermal state of the ice-I shell would cause repeated episodes of extensional and compressional tectonism.