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.     

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.     

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.     

Mass and interior of Enceladus from Cassini data analysis

Gravity results are available from radio Doppler data acquired by the Deep Space Network during the encounter of the Cassini spacecraft with Enceladus in February 2005. We report the mass of Enceladus to be (1.0798±0.0016)×¹⁰²⁰ kg, which implies a density of . For a core made of hydrated silicates with a density of 2500 kg m⁻³ the core radius is ∼190 km and the quadrupole moment C₂₂∼1.4×¹⁰⁻³. If Enceladus is in hydrostatic equilibrium, the larger than previously anticipated density implies that the recently proposed secondary spin-orbit resonance cannot be present. Therefore, the source of endogenic activity of Enceladus remains unexplained.     

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.     

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.     

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.     

Efficiency of neutrino-induced radio measurements to inspect local areas of Enceladus

The Cassini probe observed a young and smooth surface around the south pole of Enceladus, while around the north pole the surface was found to be relatively old and inactive (Porco, C.C. et al. [2006]. Science 311, 1393–1401). This heterogeneous surface implies that the ice thickness of Enceladus is not uniform between the north and south polar regions. Determining the thickness of the icy layer is important to confirm the existence of an internal ocean as well as to reveal the heating mechanism of Enceladus. We show that the measurement of radio waves induced by cosmic neutrinos can be an effective method to constrain the ice thickness of a localized area where conventional gravity or electromagnetic field measurements cannot be used. This method could be used to constrain the thickness of the icy layer on Enceladus even if the ice is a few tens of kilometers thick, measuring over a period of several years, which greatly exceeds the ability of radar sounding, and hence could be used in future orbiter missions.     

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.     

Shapes of the saturnian icy satellites and their significance

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

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.     

Degree-one convection and the origin of Enceladus' dichotomy

Recently, the Cassini spacecraft has detected ongoing geologic activity near the south pole of Saturn's moon Enceladus. In contrast, the satellite's north-polar region is heavily cratered and appears to have been geologically inactive for a long time. We propose that this hemispheric dichotomy is caused by interior dynamics with degree-one convection driving the south-polar activity. We investigate a number of core sizes and internal heating rates for which degree-one convection occurs. The numerical simulations imply that a core radius of less than 100±20 km and an energy input at a rate of 3.0 to 5.5 GW would be required for degree-one convection to prevail. This is within the range of the observed thermal power release near Enceladus' south pole. Provided that Enceladus is not fully differentiated, degree-one convection is found to be a viable mechanism to explain Enceladus' hemispheric dichotomy.     

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.     

On the stability of hypothetical satellites coorbital to Mimas or Enceladus

In the present work, we study the stability of hypothetical satellites that are coorbital with Enceladus and Mimas. We performed numerical simulations of 50 particles around the triangular Lagrangian equilibrium points of Enceladus and Mimas taking into account the perturbation of Mimas, Enceladus, Tethys, Dione, Titan and the oblateness of Saturn. All particles remain on tadpole orbits after 10 000 yr of integration. Since in the past the orbit of Enceladus and Mimas expanded due to the tidal perturbation, we also simulated the system with Enceladus and Mimas at several different values of semimajor axes. The results show that in general the particles remain on tadpole orbits. The exceptions occur when Enceladus is at semimajor axes that correspond to 6:7, 5:6 and 4:5 resonances with Mimas. Therefore, if Enceladus and Mimas had satellites librating around their Lagrangian triangular points in the past, they would have been removed if Enceladus crossed one of these first-order resonances with Mimas.     

Ring torque on janus and the melting of Enceladus

Voyager 2 images show parts of Enceladus' surface to be very smooth, lacking craters down to the resolution limit of 4 km. This absence of craters indicates geologically recent resurfacing, probably due to internal melting. However, calculations of current heating mechanisms, including radioactive decay and tidal heating due to Enceladus' resonance with Dione, yield heating rates too small to cause melting. The orbital mean motion of Janus (1980S1) is slightly less than twice that of Enceladus and, according to theoretical calculations, is currently decreasing as Janus' orbit evolves outward due to resonant torques from Saturn's rings. If Janus were ever locked into a stable 2:1 orbital commensurability with Enceladus, the resulting angular momentum transfer could have sufficiently enhanced the eccentricity of Enceladus' orbit for the ensuing tidal heating to have melted Enceladus' interior. The existence of a Laplace-like three-body resonance including Dione, although unlikely, would have increased heating. If Janus were indeed held in resonance with Enceladus until recently (10⁷–10⁸ years B.P.) when the lock was disrupted by an unspecified event (possibly a catastrophic collision which simultaneously created the coorbital pair, or by the influence of Dione) both the recent internal activity of Enceladus and the proximity of Janus to Saturn's rings may be explained. However, the predicted rapid time scale for ring evolution due to resonant torques from Saturn's inner moons remains a major problem.     

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.     

On the thermal history of Saturn’s satellites Titan and Enceladus

The thermal histories of two geologically active satellites of Saturn—Titan and Enceladus—are discussed. During the Cassini mission, it was found that there are both nitrogen-containing compounds—NH₃ and N₂-and CO₂ and CH₄ in the water plumes of Enceladus; at that, ammonia is the prevailing form. This may testify that during evolution, the material of the satellite was warmed up to T ∼ 500–600 K, when NH₃ (the form of nitrogen capable of being accreted) could only be partly converted into N₂. Contrary to Enceladus, the temperature inside Titan probably reached values higher than 800 K or even higher than 1000 K, since the process of the chemical dissociation of ammonia was completely finished on this satellite and its atmosphere contains only molecular nitrogen. While the internal heating of Titan up to high temperatures can be explained by its large mass, the heating source for Enceladus’ interior is far from evident. Such traditional heating sources as the energy of gravitational differentiation and the radiogenic heating due to shortliving ²⁶Al and ⁶⁰Fe could not be effective. The first one is because of the small size of Enceladus (RE ≈ 250 km), and the inefficiency of the second one is caused by the fact that the satellite was formed not earlier than 8–10 Myr after the formation of calcium and aluminum-enriched inclusions in carbonaceous chondrites (CAI), i.e., after ²⁶Al had completely decayed. In the present paper, we propose other heating mechanisms-the heat of long-living radioactive elements and tidal heat, which could provide the observed chemical composition of the water plumes of Enceladus rather than only the differentiation of its protomatter into the ironstone core and the ice mantle.     

Release of N2, CH4, CO2, and H2O from surface ices on Enceladus

We vapor deposit at 20 K a mixture of gases with the specific Enceladus plume composition measured in situ by the Cassini INMS [Waite, J.H., Combi, M.R., Ip, W.H., Cravens, T.E., McNutt, R.L., Kasprzak, W., Yelle, R., Luhmann, J., Niemann, H., Gell, D., Magee, B., Fletcher, G., Lunine, J., Tseng, W.L., 2006. Science 311, 1419-1422] to form a mixed molecular ice. As the sample is slowly warmed, we monitor the escaping gas quantity and composition with a mass spectrometer. Pioneering studies [Schmitt, B., Klinger, J., 1987. Different trapping mechanisms of gases by water ice and their relevance for comet nuclei. In: Rolfe, E.J., Battrick, B. (Eds.), Diversity and Similarity of Comets. SP-278. ESA, Noordwijk, The Netherlands, pp. 613-619; Bar-Nun, A., Kleinfeld, I., Kochavi, E., 1988. Phys. Rev. B 38, 7749-7754; Bar-Nun, A., Kleinfeld, I., 1989. Icarus 80, 243-253] have shown that significant quantities of volatile gases can be trapped in a water ice matrix well above the temperature at which the pure volatile ice would sublime. For our Enceladus ice mixture, a composition of escaping gases similar to that detected by Cassini in the Enceladus plume can be generated by the sublimation of the H₂O:CO₂:CH₄:N₂ mixture at temperatures between 135 and 155 K, comparable to the high temperatures inferred from the CIRS measurements [Spencer, J.R., Pearl, J.C., Segura, M., Flasar, F.M., Mamoutkine, A., Romani, P., Buratti, B.J., Hendrix, A.R., Spilker, L.J., Lopes, R.M.C., 2006. Science 311, 1401-1405] of the Enceladus “tiger stripes.” This suggests that the gas escape phenomena that we measure in our experiments are an important process contributing to the gases emitted from Enceladus. A similar experiment for ice deposited at 70 K shows that both the processes of volatile trapping and release are temperature dependent over the temperature range relevant to Enceladus.     

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.     

Tidal dissipation in small viscoelastic ice moons: The case of Enceladus

Tidal dissipation is investigated in a viscoelastic homogeneous sphere having the orbital and physical characteristics of the icy inner satellite of Saturn, Enceladus. The dissipated power is calculated for Kelvin-Voigt and Maxwell rheologies, whose dissipation function can be expressed in terms of viscosity. Expressions for the dissipated power as a function of viscosity is calculated in both cases and compared to the expression found for a lossy elastic body. A physical law relating viscosity of water ice to temperature and grain size is introduced and the feedback between dissipated power and temperature is investigated. It is found that tidal dissipation with current orbital eccentricity alone cannot account for the surface activity observed on Enceladus, if it is composed of water ice.