Distribution of icy particles across Enceladus' surface as derived from Cassini-VIMS measurements

The surface of Enceladus consists almost completely of water ice. As the band depths of water ice absorptions are sensitive to the size of particles, absorptions can be used to map variations of icy particles across the surface. The Visual and Infrared Mapping Spectrometer (VIMS) observed Enceladus with a high spatial resolution during three Cassini flybys in 2005 (orbits EN 003, EN 004 and EN 011). Based on these data we measured the band depths of water ice absorptions at 1.04, 1.25, 1.5, and 2 μm. These band depths were compared to water ice models that represent theoretically calculated reflectance spectra for a range of particle diameters between 2 μm and 1 mm. The agreement between the experimental (VIMS) and model values supports the assumption that pure water ice characterizes the surface of Enceladus and therefore that variations in band depth correspond to variations in water ice particle diameters. Our measurements show that the particle diameter of water ice increases toward younger tectonically altered surface units with the largest particles exposed in relatively “fresh” surface material. The smallest particles were generally found in old densely cratered terrains. The largest particles (∼0.2 mm) are concentrated in the so called “tiger stripes” at the south pole. In general, the particle diameters are strongly correlated with geologic features and surface ages, indicating a stratigraphic evolution of the surface that is caused by cryovolcanic resurfacing and impact gardening.     

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.     

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.     

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.     

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.     

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.     

The three-dimensional structure of Saturn’s E ring

Saturn’s diffuse E ring consists of many tiny (micron and sub-micron) grains of water ice distributed between the orbits of Mimas and Titan. Various gravitational and non-gravitational forces perturb these particles’ orbits, causing the ring’s local particle density to vary noticeably with distance from the planet, height above the ring-plane, hour angle and time. Using remote-sensing data obtained by the Cassini spacecraft in 2005 and 2006, we investigate the E-ring’s three-dimensional structure during a time when the Sun illuminated the rings from the south at high elevation angles (>15°). These observations show that the ring’s vertical thickness grows with distance from Enceladus’ orbit and its peak brightness density shifts from south to north of Saturn’s equator plane with increasing distance from the planet. These data also reveal a localized depletion in particle density near Saturn’s equatorial plane around Enceladus’ semi-major axis. Finally, variations are detected in the radial brightness profile and the vertical thickness of the ring as a function of longitude relative to the Sun. Possible physical mechanisms and processes that may be responsible for some of these structures include solar radiation pressure, variations in the ambient plasma, and electromagnetic perturbations associated with Saturn’s shadow.     

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.     

Thermal inertia and bolometric Bond albedo values for Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus as derived from Cassini/CIRS measurements

Spectra taken by Cassini’s Composite Infrared Spectrometer (CIRS) between 10 and 600 cm⁻¹ (17-1000 μm) of surface thermal emission of Mimas, Enceladus, Tethys, Dione, Rhea and Iapetus have been used to derive the thermal inertia and bolometric Bond albedo values. Only an upper limit for the bolometric Bond albedo of Iapetus’ dark leading side could be determined due to the insensitivity of the thermal model to albedo when albedos are very low. The thermal inertia in this region however is better constrained. The CIRS coverage of Enceladus is extensive enough that the latitudinal variation in these values from 60°S to 70°N has been determined in 10° wide bins. The bolometric Bond albedos determined here are consistent with literature values which show the surface of the saturnian icy moons to be covered in ice contaminated to varying degrees. The thermal inertia of the moons is shown to be in the range 9-, approximately 2-6 times lower than that of the Galilean satellites, implying a less well consolidated and more porous surface. The thermal inertias of Iapetus and Phoebe are somewhat higher, suggesting that the very low thermal inertias of satellites from Rhea inwards may be related to their probable coating of E-ring material. Latitudinal variations on the surface of Enceladus show that the bolometric Bond albedo and thermal inertia increase towards the active plume source at the south pole.     

Geysers of Enceladus: Quantitative analysis of qualitative models

Aspects of two qualitative models of Enceladus’ dust plume—the so-called “Cold Faithful” [Porco, C.C., et al., 2006. Cassini observes the active south pole of Enceladus. Science 311, 1393-1401; Ingersoll, A.P., et al., 2006. Models of the Enceladus plumes. In: Bulletin of the American Astronomical Society, vol. 38, p. 508] and “Frigid Faithful” [Kieffer, S.W., et al., 2006. A clathrate reservoir hypothesis for Enceladus’ south polar plume. Science 314, 1764; Gioia, G., et al., 2007. Unified model of tectonics and heat transport in a Frigid Enceladus. Proc. Natl. Acad. Sci. 104, 13578-13591] models—are analyzed quantitatively. The former model assumes an explosive boiling of subsurface liquid water, when pressure exerted by the ice crust is suddenly released due to an opening crack. In the latter model the existence of a deep shell of clathrates below Enceladus’ south pole is conjectured; clathrates can decompose explosively when exposed to vacuum through a fracture in the outer icy shell. For the Cold Faithful model we estimate the maximal velocity of ice grains, originating from water splashing in explosive boiling. We find that for water near the triple point this velocity is far too small to explain the observed plume properties. For the Frigid Faithful model we consider the problem of momentum transfer from gas to ice particles. It arises since any change in the direction of the gas flow in the cracks of the shell requires re-acceleration of the entrained grains. While this effect may explain the observed speed difference of gas and grains if the gas evaporates from triple point temperature (273.15 K) [Schmidt, J., et al., 2008. Formation of Enceladus dust plume. Nature 451, 685], the low temperatures of the Frigid Faithful model imply a too dilute vapor to support the observed high particle fluxes in Enceladus’ plume.     

Tidally-induced melting events as the origin of south-pole activity on Enceladus

The intense activity at the south pole of Enceladus hints at an internal water reservoir. However, there is no direct evidence of liquid water at present and its long-term stability in the interior remains problematic. By modeling heat production and transfer in the ice shell in a spherical geometry, we show that tidal heating naturally leads to a concentration of convective hot upwellings in the south polar region, favoring the preservation of liquid water at depth. We show that large volumes of water are produced within the ice shell at the south pole during periods of elevated orbital eccentricity (3–5 times the present-day value). Strong lateral variations in the melt production and crystallization rates result in stress concentration in the south polar region, thus providing an explanation for the tectonic activity observed today. We predict that an internal ocean may be sustained over the long term as the consequence of repeated periods with elevated orbital eccentricity, leading to episodic melting and resurfacing events.     

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.     

Photochemistry of methane-water ices

We report a study on the broadband ultraviolet photolysis of methane-water ice mixtures, at low methane concentrations and temperatures relevant to the icy satellites of the outer Solar System. The photochemistry of these mixtures is dominated by the action of hydroxyl radicals on methane and the resulting products. This implies that, given sufficient exposure time, the methane will eventually be completely oxidized to carbon dioxide. The presence of methane inhibits the formation of hydrogen peroxide by serving as a trap for hydroxyl radicals. The distribution of photochemical products is broadly similar to that previously conducted using ion and electron sources, with some differences possibly attributable to the difference in radiation source. The results are applicable to a variety of icy bodies in the Solar System. On Enceladus, where methane mixed with water is measured in the plumes, methane in the surface ices is subject to oxidation and will eventually be converted to CO₂. The CH stretch feature detected in the VIMS spectra of the Enceladus surface ice suggests that methane is currently being supplied to the surface ice, likely from re-condensation of the plume gas.     

A redetermination of the ice/vapor ratio of Enceladus’ plumes: Implications for sublimation and the lack of a liquid water reservoir

The discovery of plumes of H₂O vapor and ice particles erupting from the south pole of Enceladus, the tiny frigid satellite of Saturn, sparked controversy over whether these plumes are produced by boiling, or by sublimation with subsequent recondensation of the sublimated vapor [Porco, C.C., Helfenstein, P., Thomas, P.C., Ingersoll, A.P., Wisdom, J., West, R., Neukum, G., Denk, T., Wagner, R., Roatsch, T., Kieffer, S., Turtle, E., McEwen, A., Johnson, T.V., Rathbun, J., Veverka, J., Wilson, D., Perry, J., Spitale, J., Brahic, A., Burns, J.A., DelGenio, A.D., Dones, L., Murray, C.D., Squyres, S., 2006. Science 311, 1393-1401]. Porco et al.’s analysis that the masses of ice (I) and vapor (V) in the plume were comparable was taken to argue against the occurrence of sublimation and recondensation, leading to the hypothesis that the reservoir was boiling water, possibly as close as 7 m to the surface. Thus, it has been advocated that Enceladus should be a target for astrobiology exploration. Here we show, with recalculations using the original data and methodologies, as well as with new sensitivity studies, that the mass of ice in the column is significantly less than the mass of water vapor, and that by considering three additional effects, I/V is likely to be <0.2-0.1. This means that the plume is dominated by vapor that the thermodynamics permits to be easily produced by sublimation with recondensation. The low I/V ratio provides no compelling criterion for consideration of a liquid water reservoir. The uncertainties on the I/V ratio have not previously been discussed in the literature. Although the I/V ratio is sensitive to particle sizes and size distributions, the masses of ice (I) and vapor (V) are not comparable in any scenario constrained by available observations. We thus discuss the implications of sublimation from a thermodynamic point of view in a context that has not been presented previously. Constraints on I/V ratio from future spacecraft measurements of the plume, in conjunction with consideration of the total plume composition and multicomponent analysis, can help constrain source conditions for the plume.     

The E-ring in the vicinity of Enceladus: II. Probing the moon's interior—The composition of E-ring particles

The population of Saturn's outermost tenuous E-ring is dominated by tiny water ice particles. Active volcanism on the moon Enceladus, embedded in the E-ring, has since late 2005 been known to be a major source of particles replenishing the ring. Therefore particles in the vicinity of Enceladus may provide crucial information about the dynamical and chemical processes occurring below the moon's icy surface. Here we present a statistical evaluation of more than 2000 impact ionisation mass spectra of Saturn's E-ring particles, with sizes predominantly below 1 μm, detected by the Cosmic Dust Analyser onboard the Cassini spacecraft. We focus on the identification of non-water features in spectra otherwise dominated by water ice signatures. Here we specify the categorisation of two different spectrum types, which probably represent two particle populations. Type I spectra imply pure water ice particles, whereas in Type II spectra organic compounds and/or silicate minerals are identified as impurities within the icy particles. This finding supports the hypothesis of a dynamic interaction of Enceladus' rocky core with liquid water.     

Radar imagery of Mercury’s putative polar ice: 1999-2005 Arecibo results

We present an updated survey of Mercury’s putative polar ice deposits, based on high-resolution (1.5-km) imaging with the upgraded Arecibo S-band radar during 1999-2005. The north pole has now been imaged over a full range of longitude aspects, making it possible to distinguish ice-free areas from radar-shadowed areas and thus better map the distribution of radar-bright ice. The new imagery of the south pole, though derived from only a single pair of dates in 2005, improves on the pre-upgrade Arecibo imagery and reveals many additional ice features. Some medium-size craters located within 3° of the north pole show near-complete ice coverage over their floors, central peaks, and southern interior rim walls and little or no ice on their northern rim walls, while one large (90 km) crater at 85°N shows a sharp ice-cutoff line running across its central floor. All of this is consistent with the estimated polar extent of permanent shading from direct sunlight. Some craters show ice in regions that, though permanently shaded, should be too warm to maintain unprotected surface ice owing to indirect heating by reflected and reradiated sunlight. However, the ice distribution in these craters is in good agreement with models invoking insulation by a thin dust mantle. Comparisons with Goldstone X-band radar imagery indicate a wavelength dependence that could be consistent with such a dust mantle. More than a dozen small ice features have been found at latitudes between 67° and 75°. All of this low-latitude ice is probably sheltered in or under steep pole-facing crater rim walls, although, since most is located in the Mariner-unimaged hemisphere, confirmation must await imaging by the MESSENGER orbiter. These low-latitude features are concentrated toward the “cold longitudes,” possibly indicating a thermal segregation effect governed by indirect heating. The radar imagery places the corrected locations of the north and south poles at 7°W, 88.35°N and 90°W, 88.7°S, respectively, on the original Mariner-based maps.     

On the origin of south polar folds on Enceladus

Recent high-resolution Cassini images of the south polar terrain of Enceladus reveal regions of short-wavelength deformation, inferred to be compressional folds between the Baghdad and Damascus tiger stripes (Spencer, J.R., Barr, A.C., Esposito, L.W., Helfenstein, P., Ingersoll, A.P., Jaumann, R., McKay, C.P., Nimmo, F., Waite, J.H. [2009a]. Enceladus: An active cryovolcanic satellite. In: Saturn after Cassini-Huygens. Springer, New York, pp. 683-722). Here, we use Fourier analysis of the bright/dark variations to show that the folds have a dominant wavelength of 1.1 ± 0.4 km. We use the simple model of lava flow folding from Fink (Fink, J. [1980]. Geology 8, 250-254) to show that the folds could form in an ice shell with an upper high-viscosity boundary layer of thickness <400 m, with a driving stress of 40-80 kPa, and strain rate between 10⁻¹⁴ s⁻¹ and 10⁻¹² s⁻¹. Such deformation rates imply resurfacing of the SPT in 0.05-5 Myr, consistent with its estimated surface age. Measurements of fold topography and more sophisticated numerical modeling can narrow down the conditions of fold formation and provide valuable constraints on the thermal structure of the ice shell on 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.     

Ionization chemistry in H<Subscript>2</Subscript>O-dominated atmospheres of icy moons

The processes of the formation and dynamics of tenuous gaseous envelopes of icy moons in giant-planet systems are considered. Tenuous exospheres with relatively dense surface layers are likely to form around more massive icy satellites, such as, for example, the Galilean satellites Europa and Ganymede in the Jovian system. Escaping exospheres are formed in the case of low-mass icy moons, as happens for the icy satellite Enceladus in the Saturnian system. The main parent component of such gaseous envelopes is water vapor, which enters into the atmosphere as a result of thermal degassing processes, nonthermal radiolysis, and other active processes and phenomena on the icy surface of a satellite. A numerical kinetic model has been developed to study on a molecular level the processes of the formation, chemical evolution, and dynamics of tenuous gaseous envelopes dominated mainly by H₂O. The ionization processes in such tenuous gaseous envelopes are caused by solar ultraviolet (UV) radiation and solar-wind and/or magnetospheric plasma. The primary processes when ultraviolet solar photons and plasma electrons affect the tenuous gas of the H₂O-dominated atmosphere are responsible for the chemical diversity of the gaseous envelopes of icy moons. Ionization chemistry, including ion-molecular reactions, dissociative recombination of molecular ions, and the reactions of the charge exchange with magnetospheric ions, is important for the formation of chemical diversity in gaseous envelopes of icy satellites. The model considered in the study was used to numerically simulate the formation and development of chemical diversity in the tenuous gaseous envelope of Enceladus. The numerical results were compared to the direct Cassini measurements during its close flyby near Enceladus.