Coupled convection and tidal dissipation in Europa’s ice shell using non-Newtonian grain-size-sensitive (GSS) creep rheology

We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10⁻⁵), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10⁻⁵ on Europa.     

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.     

Effects of low-viscous layers and a non-zero obliquity on surface stresses induced by diurnal tides and non-synchronous rotation: The case of Europa

In this study we present a semi-analytical Maxwell-viscoelastic model of the variable tidal stress field acting on Europa’s surface. In our analysis, we take into account surface stresses induced by the small eccentricity of Europa’s orbit, the non-zero obliquity of Europa’s spin axis - both acting on a diurnal 3.55-days timescale - and the reorientation of the ice shell as a result of non-synchronous rotation (NSR). We assume that Europa’s putative ocean is covered by an ice shell, which we subdivide in a low-viscous and warm lower ice layer (asthenosphere, viscosity 10¹²-10¹⁷ Pa s), and a high-viscous and cold upper ice layer (lithosphere, viscosity 10²¹ Pa s).Viscoelastic relaxation influences surface stresses in two ways: (1) through viscoelastic relaxation in the lithosphere and (2) through the viscoelastic tidal response of Europa’s interior. The amount of relaxation in the lithosphere is proportional to the ratio between the period of the forcing mechanism and the Maxwell time of the high-viscous lithosphere. As a result, this effect is only relevant to surface stresses caused by the slow NSR mechanism. On the other hand, the importance of the viscoelastic response on surface stresses is proportional to the ratio between the relaxation time (τ_j) of a given viscoelastic mode j and the period of the forcing function. On a diurnal timescale the fast relaxation of transient modes related to the low viscosity of the asthenosphere can alter the magnitude and phase shift of the diurnal stress field at Europa’s surface. The effects are largest, up to 20% in magnitude and 7° in phase for ice rigidities lower than 3.487 GPa, when the relaxation time of the aforementioned transient modes approaches the inverse of the average angular rate of Europa’s orbit. On timescales relevant for NSR (>10⁴ years) the magnitude and phase shift of NSR surface stresses can be affected by viscoelastic relaxation of the ocean-ice boundary. This effect, however, becomes only important when the behavior of the lithosphere w.r.t. NSR approaches the fluid limit, i.e. for strong relaxation in the lithosphere. The combination of NSR and diurnal stresses for different amounts of viscoelastic relaxation of NSR stresses in the lithosphere leads to a large variety of global stress fields that can explain the formation of the large diversity of lineament morphologies observed on Europa’s surface. Variation of the amount of relaxation in the lithosphere is likely due to changes in the spin rate of Europa and/or the rheological properties of the surface.In addition, we show that a small obliquity(<1°) can have a considerable effect on Europa’s diurnal stress field. A non-zero obliquity breaks the symmetric distribution of stress patterns with respect to the equator, thereby affecting the magnitude and orientation of the principal stresses at the surface. As expected, increasing the value of Europa’s obliquity leads to larger diurnal stresses at the surface, especially when Europa is located 90° away from the nodes formed by the intersection of its orbital and equatorial planes.     

Heat flow, lenticulae spacing, and possibility of convection in the ice shell of europa

Two opposing models to explain the geological features observed on Europa’s surface have been proposed. The thin-shell model states that the ice shell is only a few kilometers thick, transfers heat by conduction only, and can become locally thinner until it exposes an underlying ocean on the satellite’s surface. According to the thick-shell model, the ice shell may be several tens of kilometers thick and have a lower convective layer, above which there is a cold stagnant lid that dissipates heat by conduction. Whichever the case, from magnetic data there is strong support for the presence of a layer of salty liquid water under the ice. The present study was performed to examine whether the possibility of convection is theoretically consistent with surface heat flows of ∼100-200 mW m⁻², deduced from a thin brittle lithosphere, and with the typical spacing of 15-23 km proposed for the features usually known as lenticulae. It was obtained that under Europa’s ice shell conditions convection could occur and also account for high heat flows due to tidal heating of the convective (nearly isothermal) interior, but only if the dominant water ice rheology is superplastic flow (with activation energy of 49 kJ mol⁻¹; this is the rheology thought dominant in the warm interior of the ice shell). In this case the ice shell would be ∼15-50 km thick. Furthermore, in this scenario explaining the origin of the lenticulae related to convective processes requires ice grain size close to 1 mm and ice thickness around 15-20 km.     

Thermal-orbital evolution of Io and Europa

Coupled thermal-orbital evolution models of Europa and Io are presented. It is assumed that Io, Europa, and Ganymede evolve in the Laplace resonance and that tidal dissipation of orbital energy is an internal heat source for both Io and Europa. While dissipation in Io occurs in the mantle as in the mantle dissipation model of Segatz et al. (1988, Icarus 75, 187), two models for Europa are considered. In the first model dissipation occurs in the silicate mantle while in the second model dissipation occurs in the ice shell. In the latter model, ice shell melting and variations of the shell thickness above an ocean are explicitly included. The rheology of both the ice and the rock is cast in terms of a viscoelastic Maxwell rheology with viscosity and shear modulus depending on the average temperature of the dissipating layer. Heat transfer by convection is calculated using a parameterization for strongly temperature-dependent viscosity convection. Both models are consistent with the present orbital elements of Io, Europa, and Ganymede. It is shown that there may be phases of quasi-steady evolution with large or small dissipation rates (in comparison with radiogenic heating), phases with runaway heating or cooling and oscillatory phases during which the eccentricity and the tidal heating rate will oscillate. Europa's ice thickness varies between roughly 3 and 70 km (dissipation in the silicate layer) or 10 and 60 km (dissipation in the ice layer), suggesting that Europa's ocean existed for geological timescales. The variation in ice thickness, including both convective and purely conductive phases, may be reflected in the formation of different geological surface features on Europa. Both models suggest that at present Europa's ice thickness is several tens of km thick and is increasing, while the eccentricity decreases, implying that the satellites evolve out of resonance. Including lithospheric growth in the models makes it impossible to match the high heat flux constraint for Io. Other heat transfer processes than conduction through the lithosphere must be important for the present Io.     

A Melt-through Model for Chaos Formation on Europa

The character of chaotic terrain on Europa is consistent with its formation by the melting of a thin conducting ice shell from below. Tidal dissipation can provide adequate energy for such a process. For example, only a few percent of Europa's predicted tidal heat, spread over a region 200 km in diameter, can lead to large melt regions within a few tens of thousands of years. Stronger, more localized concentrations result in melt-through in significantly shorter times (i.e., a few hundred years). The time scale for melt-through is shorter than the time scale for the solid-state viscous inflow of ice by several orders of magnitude. In general, modest concentrations of tidal heat can melt ice away faster than viscous inflow, leading to melt-through. A mechanism to transmit these heat concentrations through the ocean is required for this model. Such heat transport could be the result of convective plumes in the ocean driven by seafloor volcanism or by the destabilization of a stratified ocean.     

Thermal Equilibrium States of Europa's Ice Shell: Implications for Internal Ocean Thickness and Surface Heat Flow

Ice-shell thickness and ocean depth are calculated for steady state models of tidal dissipation in Europa's ice shell using the present-day values of the orbital elements. The tidal dissipation rate is obtained using a viscoelastic Maxwell rheology for the ice, the viscosity of which has been varied over a wide range, and is found to strongly increase if an (inviscid) internal ocean is present. To determine steady state values, the tidal dissipation rate is equated to the heat-transfer rate through the ice shell calculated from a parameterized model of convective heat transfer or from a thermal conduction model, if the ice layer is found to be stable against convection. Although high dissipation rates and heat fluxes of up to 300 mWm⁻² are, in principle, possible for Europa, these values are unrealistic because the states for which they are obtained are thermodynamically unstable. Equilibrium models have surface heat flows around 20 mWm⁻² and ice-layer thicknesses around 30 km, which is significantly less than the total thickness of the H₂O-layer. These results support models of Europa with ice shells a few tens of kilometers thick and around 100-km-thick subsurface oceans.     

Constraints on Europa’s rotational dynamics from modeling of tidally-driven fractures

Cycloids, arcuate features observed on Europa’s surface, have been interpreted as tensile cracks that form in response to diurnal tidal stress caused by Europa’s orbital eccentricity. Stress from non-synchronous rotation may also contribute to tidal stress, and its influence on cycloid shapes has been investigated as well. Obliquity, fast precession, and physical libration would contribute to tidal stress but have often been neglected because they were expected to be negligibly small. However, more sophisticated analyses that include the influence of Jupiter’s other large satellites and the state of Europa’s interior indicate that perhaps these rotational parameters are large enough to alter the tidal stress field and the formation of tidally-driven fractures. We test tidal models that include obliquity, fast precession, stress due to non-synchronous rotation, and physical libration by comparing how well each model reproduces observed cycloids. To do this, we have designed and implemented an automated parameter-searching algorithm that relies on a quantitative measure of fit quality, which we use to identify the best fits to observed cycloids. We then apply statistical techniques to determine the tidal model best supported by the data. By incorporating obliquity, fits to observed southern hemisphere cycloids improve, and we can reproduce equatorial and equator-crossing cycloids. Furthermore, we find that obliquity plus physical libration is the tidal model best supported by the data. With this model, the obliquities range from 0.32° to 1.35°. The libration amplitudes are 0.72–2.44°, and the libration phases are −6.04° to 17.72° with one outlier at 84.5°. The variability in obliquity is expected if Europa’s ice shell is mechanically decoupled from the interior, and the libration amplitudes are plausible in the presence of a subsurface ocean. Indeed, the presence of a decoupling ocean may result in feedbacks that cause all of these rotational parameters to become time-variable.     

Seas under ice: Stability of liquid-water oceans within icy worlds

The present-day existence of internal oceans under the outer ice shell of several icy satellites of the Solar System has been recently proposed. The presence of antifreeze substances decreasing ice’s melting point (and tidal heating in Europa’s case) has been generally believed to allow the stability of such oceans; limited cooling of the water (ice plus liquid) layer, due to stability against convection or to stagnant lid convection in the icy shell, have been also considered. Here we propose that even pure liquid-water oceans could survive today within several icy worlds, and we consider some factors affecting thermal modeling in these bodies. So, the existence of such oceans would be a natural consequence of the physical properties of water ice, independently from the addition of antifreeze substances or any other special conditions. The inclusion of these substances would contribute to expand the conditions for water to stay liquid and to increase ocean’s volume.     

Equilibrium Convection on a Tidally Heated and Stressed Icy Shell of Europa for a Composite Water Ice Rheology

Water ice I rheology is a key factor for understanding the thermal and mechanical state of the outer shell of the icy satellites. Ice flow involves several deformation mechanisms (both Newtonian and non-Newtonian), which contribute to different extents depending on the temperature, grain size, and applied stress. In this work I analyze tidally heated and stressed equilibrium convection in the ice shell of Europa by considering a composite viscosity law which includes diffusion creep, basal slip, grain boundary sliding and dislocation creep, and. The calculations take into account the effect of tidal stresses on ice flow and use grain sizes between 0.1 and 100 mm. An Arrhenius-type relation (useful for parameterized convective models) is found then by fitting the calculated viscosity between 170 and 273 K to an exponential regression, which can be expressed in terms of pre-exponential constant and effective activation energy. I obtain convective heat flows between ~40 and ~60 mW m^?2, values lower than those usually deduced (~100 mW m^?2) from geological indicators of lithospheric thermal state, probably indicating heterogeneous tidal heating. On the other hand, for grain sizes larger than ~0.3 mm the thicknesses of the ice shell and convective sublayer are ~20–30 km and ~5–20 km respectively, values in good agreement with the available information for Europa. So, some fundamental geophysical characteristics of the ice shell of Europa could be arising from the properties of the composite water ice rheology.     

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.     

A whole-moon thermal history model of Europa: Impact of hydrothermal circulation and salt transport

A whole-moon numerical model of Europa is developed to simulate its thermal history. The thermal evolution covers three phases: (i) an initial, roughly 0.5 Gyr-long period of radiogenic heating and differentiation, (ii) a long period from 0.5 Gyr to 4 Gyr with continuing radiogenic heating but no tidal dissipative heating (TDH), and (iii) a final period covering the last 0.5 Gyr until the present, during which TDH is active. Hydrothermal plumes develop after the initial period of heating and differentiation and transport heat and salt from Europa’s silicate mantle to its ice shell. We find that, even without TDH, vigorous hydrothermal convection in the rocky mantle can sustain flow in an ocean layer throughout Europa’s history. When TDH becomes active, the ice shell melts quickly to a thickness of about 20 km, leaving an ocean 80 km or more deep. Parameterized convection in the ice shell is non-uniform spatially, changes over time, and is tied to the deeper ocean–mantle dynamics. We also find that the dynamics are affected by salt concentrations. An initially non-uniform salt distribution retards plume penetration, but is homogenized over time by turbulent diffusion and time-dependent flow driven by initial thermal gradients. After homogenization, the uniformly distributed salt concentrations are no longer a major factor in controlling plume transport. Salt transport leads to the formation of a heterogeneous brine layer and salt inclusions at the bottom of the ice shell; the presence of salt in the ice shell could strongly influence convection in that layer.     

Possibility of Convection for Diffusion (Newtonian) Viscosity in the Ice Shell of Europa?

Recent investigations into convection in Europa’s ice shell have been based on non-Newtonian (stress-dependent) or Newtonian (stress-independent) viscosity for water ice. However, despite the wide use of Newtonian convection, experimentally observed water ice flow is non-Newtonian, and analysis of stability against convection of the ice shell using updated flow laws has been only performed for non-Newtonian rheologies. Here we use the flow law proposed for diffusion creep to analyze the possibility of the onset of convection for Newtonian viscosity in relation to the thermal state of Europa. Our findings indicate that for diffusion creep convection might have started, but that significantly lower heat flows (and equivalently higher shell thicknesses) and/or grain sizes are required than for superplastic flow, which is the most probable flow mechanism if the ice shell is convective.     

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 orbital-thermal evolution and global expansion of Ganymede

The tectonically and cryovolcanically resurfaced terrains of Ganymede attest to the satellite's turbulent geologic history. Yet, the ultimate cause of its geologic violence remains unknown. One plausible scenario suggests that the Galilean satellites passed through one or more Laplace-like resonances before evolving into the current Laplace resonance. Passage through such a resonance can excite Ganymede's eccentricity, leading to tidal dissipation within the ice shell. To evaluate the effects of resonance passage on Ganymede's thermal history we model the coupled orbital-thermal evolution of Ganymede both with and without passage through a Laplace-like resonance. In the absence of tidal dissipation, radiogenic heating alone is capable of creating large internal oceans within Ganymede if the ice grain size is 1 mm or greater. For larger grain sizes, oceans will exist into the present epoch. The inclusion of tidal dissipation significantly alters Ganymede's thermal history, and for some parameters (e.g. ice grain size, tidal Q of Jupiter) a thin ice shell (5 to 20 km) can be maintained throughout the period of resonance passage. The pulse of tidal heating that accompanies Laplace-like resonance capture can cause up to 2.5% volumetric expansion of the satellite and contemporaneous formation of near surface partial melt. The presence of a thin ice shell and high satellite orbital eccentricity would generate moderate diurnal tidal stresses in Ganymede's ice shell. Larger stresses result if the ice shell rotates non-synchronously. The combined effects of satellite expansion, its associated tensile stress, rapid formation of near surface partial melt, and tidal stress due to an eccentric orbit may be responsible for creating Ganymede's unique surface features.     

Shell thickness variations and the long-wavelength topography of Titan

The long-wavelength topography of Titan has an amplitude larger than that expected from tidal and rotational distortions at its current distance from Saturn. This topography is associated with small gravity anomalies, indicating a high degree of compensation. Both observations can be explained if Titan has a floating, isostatically-compensated ice shell with a spatially-varying thickness. The spatial variations arise because of laterally-variable tidal heating within the ice shell. Models incorporating shell thickness variations result in an improved fit to the observations and a degree-two tidal Love number h_2t consistent with expectations, without requiring Titan to have moved away from Saturn. Our preferred models have a mean shell thickness of ≈100 km in agreement with the observed gravity anomalies, and a heat flux appropriate to a chondritic Titan. Shell thickness variations are eliminated by convection; we therefore conclude that Titan’s ice shell is not convecting at the present day.     

Europa’s ridged plains and smooth low albedo plains: Distinctive compositions and compositional gradients at the leading side-trailing side boundary

This investigation uses linear mixture modeling employing cryogenic laboratory reference spectra to estimate surface compositions and water ice grain sizes of Europa’s ridged plains and smooth low albedo plains. Near-infrared spectra for 23 exposures of ridged plains materials are analyzed along with 11 spectra representing low albedo plains. Modeling indicates that these geologic units differ both in the relative abundance of non-ice hydrated species and in the abundance and grain sizes of water ice. The background ridged plains in our study area appear to consist predominantly of water ice (∼46%) with approximately equal amounts (on average) of hydrated sulfuric acid (∼27%) and hydrated salts (∼27%). The solutions for the smooth low albedo plains are dominated by hydrated salts (∼62%), with a relatively low mean abundance of water ice (∼10%), and an abundance of hydrated sulfuric acid similar to that found in ridged plains (∼27%). The model yields larger water ice grain sizes (100 μm versus 50-75 μm) in the ridged plains. The 1.5-μm water ice absorption band minimum is found at shorter wavelengths in the low albedo plains deposits than in the ridged plains (1.498 ± .003 μm versus 1.504 ± .001 μm). The 2.0-μm band minimum in the low albedo plains exhibits a somewhat larger blueshift (1.964 ± .006 μm versus 1.983 ± .006 μm for the ridged plains).The study area spans longitudes from 168° to 185°W, which includes Europa’s leading side-trailing side boundary. A well-defined spatial gradient of sulfuric acid hydrate abundance is found for both geologic units, with concentrations increasing in the direction of the trailing side apex. We associate this distribution with the exogenic effects of magnetospheric charged particle bombardment and associated chemical processing of surface materials (the radiolytic sulfur cycle). However, one family of low albedo plains exposures exhibits sulfuric acid hydrate abundances up to 33% lower than found for adjacent exposures, suggesting that these materials have undergone less processing, thus implying that these deposits may have been emplaced more recently.Modeling identifies high abundances (to 30%) of magnesium sulfate brines in the low albedo plains exposures. Our investigation marks the first spectroscopic identification of MgSO₄ brine on Europa. We also find significantly higher abundances of sodium-bearing species (bloedite and mirabilite) in the low albedo plains. The results illuminate the role of radiolytic processes in modifying the surface composition of Europa, and may provide new constraints for models of the composition of Europa’s putative subsurface ocean.     

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.     

Convective-conductive transitions and sensitivity of a convecting ice shell to perturbations in heat flux and tidal-heating rate: Implications for Europa

We investigate the response of conductive and convective ice shells on Europa to variations of heat flux and interior tidal-heating rate. We present numerical simulations of convection in Europa's ice shell with Newtonian, temperature-dependent viscosity and tidal heating. Modest variations in the heat flux supplied to the base of a convective ice shell, ΔF, can cause large variations of the ice-shell thickness Δδ. In contrast, for a conductive ice shell, large ΔF involves relatively small Δδ. We demonstrate that, for a fluid with temperature-dependent viscosity, the heat flux undergoes a finite-amplitude jump at the critical Rayleigh number Racr. This jump implies that, for a range of heat fluxes relevant to Europa, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. We show that, as a result, modest variations in heat flux near the critical Rayleigh number can force the ice shell to switch between the thin, conductive and thick, convective configurations over a ∼10⁷-year interval, with thickness changes of up to ∼10-30 km. Depending on the orbital and thermal history, such switches might occur repeatedly. However, existing evolution models based on parameterized-convection schemes have to date not allowed these transitions to occur. Rapid thickening of the ice shell would cause radial expansion of Europa, which could produce extensional tectonic features such as fractures or bands. Furthermore, based on interpretations for how features such as chaos and ridges are formed, several authors have suggested that Europa's ice shell has recently undergone changes in thickness. Our model provides a mechanism for such changes to occur.