Implications from Galileo Observations on the Interior Structure and Chemistry of the Galilean Satellites

Data from the recent gravity measurements by the Galileo mission are used to construct wide ranges of interior structure and composition models for the Galilean satellites of Jupiter. These models show that mantle densities of Io and Europa are consistent with an olivine-dominated mineralogy with the ratios of Mg to Fe components depending on mantle temperature for Io and on ice shell thickness for Europa. The mantle density and composition depend relatively little on core composition. The size of the core is largely determined by the core's composition with core radius increasing with the concentration of a light component such as sulfur. For Io, the range of possible core sizes is between 38 and 53% of the satellite's radius. For Europa, there is also a substantial effect of the thickness of the ice layer which is varied between 120 and 170 km on the core size. Core sizes are between 10 and 45% of Europa's radius. The core size of Ganymede ranges between one-quarter and one-third of the surface radius depending on its sulfur content and the thickness of the ice shell. A subset of the Ganymede models is consistent with an olivine-dominated mantle mineralogy. The thickness of the silicate mantle above the core varies between 900 and 1100 km. The outermost ice shell is about 900 km in thickness and is further subdivided by pressure-induced phase transitions into ice I, ice III, ice V, and ice VI layers. Callisto should be differentiated, albeit incompletely. It is proposed that this satellite was never molten at a large scale but differentiated through the convective gradual unmixing of the ice and the metal/rock component. Bulk iron-to-silicon ratios Fe/Si calculated for the inner pair of satellites, Io and Europa, are less than the CI carbonaceous chondrite value of 1.7±0.1, whereas ratios for the outer pair, Ganymede and Callisto, cover a broad range above the chondritic value. Although the ratios are uncertain, in particular for Ganymede and Callisto, the values are sufficiently distinct to suggest a difference in composition between these two pairs of satellites. This may indicate a difference in iron-silicon fractionation during the formation of both classes of satellites in the protojovian nebula.     

Librations of the Galilean satellites: The influence of global internal liquid layers

The four Galilean satellites are thought to harbor one or even two global internal liquid layers beneath their surface layer. The iron core of Io and Ganymede is most likely (partially) liquid and also the core of Europa may be liquid. Furthermore, there are strong indications for the existence of a subsurface ocean in Europa, Ganymede, and Callisto. Here, we investigate whether libration observations can be used to prove the existence of these liquid layers and to constrain the thickness of the overlying solid layers. For Io, the presence of a small liquid core increases the libration of the mantle by a few percent with respect to an entirely solid Io and mantle libration observations could be used to determine the mantle thickness with a precision of several tens of kilometers given that the libration amplitude can be measured with a precision of 1 m. For Europa, Ganymede, and Callisto, the presence of a water ocean close to the surface increases by at least an order of magnitude the ice shell libration amplitude with respect to an entirely solid satellite. The shell libration depends essentially on the shell thickness and to a minor extent on the density difference between the ocean and the ice shell. The possible presence of a liquid core inside Europa and Ganymede has no noticeable influence on their shell libration. For a precision of several meters on the libration measurements, in agreement with the expected accuracy with the NASA/ESA EJSM orbiter mission to Europa and Ganymede, an error on the shell thickness of a few tens kilometers is expected. Therefore, libration measurements can be used to detect liquid layers such as Io’s core or water subsurface oceans in Europa, Ganymede, and Callisto and to constrain the thickness of the overlying solid surface layers.     

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.     

Secondary and sesquinary craters on Europa

We address impact cratering on Io and Europa, with the emphasis on the origin of small craters on Europa as secondary to the primary impacts of comets on Io, Europa, and Ganymede. In passing we also address the origin of secondary craters generated by Zunil, a well-studied impact crater on Mars that is a plausible analog to impact craters on Io. At nominal impact rates, and taking volcanic resurfacing into account, we find that there should be 1.3 impact craters on Io, equally likely to be of any diameter between 100 m and 20 km. The corresponding model age of Europa's surface is between 60 and 100 Ma. This range of ages does not include a factor three uncertainty stemming from the uncertain sizes and numbers of comets. The mass of basaltic impact ejecta from Io to reach Europa is found to meet or exceed the micrometeoroid flux as a source of rock-forming elements to Europa's ice crust. To describe impact ejecta in more detail we adapt models for impact-generated spalls and Grady-Kipp fragments originally developed by Melosh. Our model successfully reproduces the observed size-number distributions of small craters on both Mars and Europa. However, the model predicts that a significant fraction of the 200-500 m diameter craters on Europa are not traditional secondary craters but are instead sesquinary craters caused by impact ejecta from Io that had gone into orbit about Jupiter. This prediction is not supported by observation, which implies that high speed spalls usually break up into smaller fragments that make smaller sesquinary craters. Iogenic basalts are also interesting because they provide stratigraphic horizons on Europa that in principle could be used to track historic motions of the ice, and they provide materials suitable to radiometric dating of Europa's surface.     

The production of Ganymede's magnetic field

One of the great discoveries of NASA's Galileo mission was the presence of an intrinsically produced magnetic field at Ganymede. Generation of the relatively strong (750 nT) field likely requires dynamo action in Ganymede's metallic core, but how such a dynamo has been maintained into the present epoch remains uncertain. Using a one-dimensional, three layer thermal model of Ganymede, we find that magnetic field generation can only occur if the sulfur mass fraction in Ganymede's core is very low (?3%) or very high (?21%), and the silicate mantle can cool rapidly (i.e. it has a viscosity like wet olivine). However, these requirements are not necessarily compatible with cosmochemical and physical models of the satellite. We therefore investigate an alternative scenario for producing Ganymede's magnetic field in which passage through an eccentricity pumping Laplace-like resonance in Ganymede's past enables present day dynamo action in the metallic core. If sufficient tidal dissipation occurs in Ganymede's silicate mantle during resonance passage, silicate temperatures can undergo a runaway which prevents the core from cooling until the resonance passage ends. The rapid silicate and core cooling that follows resonance escape triggers dynamo action via thermal and/or compositional convection. To test the feasibility of this mechanism we couple our thermal model with an orbital evolution model to examine the effects of resonance passage on Ganymede's silicate mantle and metallic core. We find that, contrary to expectations, there are no physically plausible scenarios in which tidal heating in the silicates is sufficient to cause the thermal runaway necessary to prevent core cooling. These findings are robust to variations in the silicate rheology, tidal dissipation factor of Jupiter (QJ), structure of the ice shell, and the inclusion of partial melting in the silicate mantle. Resonance passage therefore appears unlikely to explain Ganymede's magnetic field and we must appeal to the special conditions described above to explain the presence of the field.     

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.     

Thermal equilibrium in Europa's ice shell

Models of tidal-convective equilibrium for Europa's ice shell are computed using a laboratory-derived composite flow law for ice. Volume diffusion creep is found to dominate the flow law at equilibrium, and thus the thickness of the shell is strongly dependent on the poorly known grain size of the ice. This grain size is, however, constrained to be less than a few millimeters if equilibrium is achieved at the current eccentricity. Europa's ice shell cannot be thinner than about 16 km in equilibrium at present, since tidal dissipation cannot generate enough heat in such a thin shell to balance the heat transport. No conductive equilibria are found; this is likely due to the fact that most of a conductive shell must be cold if temperature gradients are to be large enough to carry the heat. A minimum eccentricity of about 0.0025 (about 1/4 the present value) below which there are no equilibria is also found.     

Empirical constraints on the salinity of the europan ocean and implications for a thin ice shell

Induced electrical currents within Europa inferred from Galileo spacecraft magnetometer instrument data have been interpreted as due to a salty europan ocean. Published compositional models for Europa's ocean, based on aqueous leaching of carbonaceous chondrites, range over five orders of magnitude in predicted magnesium sulfate concentrations. We combine the Galileo spacecraft magnetometer-derived oceanic conductivities and radio Doppler data-derived interior models with laboratory conductivity vs concentration data for both magnesium sulfate solutions and terrestrial seawater to determine empirically the range of salt concentrations permitted for Europa's ocean. Solutions for both a three-layer spherical model, and a five-layer half-space model, that satisfy current preferred best fits to magnetometer data imply high, near-saturation salt concentrations and require a europan ice shell of less than 15 km thick, with a best fit at 4 km ice thickness. Adding a conductive core and mantle has a negligible effect on the amplitude when ocean conductivities are greater than a few Siemens per meter. Similarly, we find that including a realistic ionosphere has a negligible effect. We examine the implications of these results for the subsurface habitability of Europa.     

Thermal and topographic tests of Europa chaos formation models from Galileo E15 observations

Stereo topography of an area near Tyre impact crater, Europa, reveals chaos regions characterised by marginal cliffs and domical topography, rising to 100-200 m above the background plains. The regions contain blocks which have both rotated and tilted. We tested two models of chaos formation: a hybrid diapir model, in which chaos topography is caused by thermal or compositional buoyancy, and block motion occurs due to the presence of near-surface (1-3 km) melt; and a melt-through model, in which chaos regions are caused by melting and refreezing of the ice shell. None of the hybrid diapir models tested generate any melt within 1-3 km of the surface, owing to the low surface temperature. A model of ocean refreezing following melt-through gives effective elastic thicknesses and ice shell thicknesses of 0.1-0.3 and 0.5-2 km, respectively. However, for such low shell thicknesses the refreezing model requires implausibly large lateral density contrasts (50-100 kg m⁻³) to explain the elevation of the centres of the chaos regions. Although a global equilibrium ice shell thickness of ≈2 km is possible if Europa's mantle resembles that of Io, it is unclear whether local melt-through events are energetically possible. Thus, neither of the models tested here gives a completely satisfactory explanation for the formation of chaos regions. We suggest that surface extrusion of warm ice may be an important component of chaos terrain formation, and demonstrate that such extrusion is possible for likely ice parameters.     

Internal structure of Europa and Callisto

Models of the internal structure of completely differentiated Europa and partially differentiated Callisto have been constructed on the basis of Galileo gravity measurements, geochemical constraints on composition of ordinary and carbonaceous chondrites, and thermodynamic data on the equations of state of water, high-pressure ices, and meteoritic material. We assume thermal and mechanical equilibrium for the interiors of the satellites. A geophysically and geochemically permissible thickness of Europa's outer water-ice shell lies between 105 and 160 km (6.2-9.2% of total mass). Our results show that the bulk composition of the rock-iron core of Europa may be described by material approaching the L/LL-type chondrites in composition, but cannot be correlated either with the material of CI chondrites or H chondrites. For Europa's L/LL-chondritic models, core radii are estimated to be 470-640 km (5.3-12.5% of total mass). The allowed thickness of Europa's H₂O layer ranges from 115±10 km for a differentiated L/LL-type chondritic mantle with a crust to 135±10 km for an undifferentiated mantle. We show that Callisto must only be partially differentiated into an outer ice-I layer, a water ocean, a rock-ice mantle, and a rock-iron core (mixture of anhydrous silicates and/or hydrous silicates + FeFeS alloy). We accept that the composition of the rock-iron material of Callisto is similar to the bulk composition of L/LL-type chondritic material containing up to 10-15% of iron and iron sulfide. Assuming conductive heat transfer through the ice-I crust [Ruiz, 2001. The stability against freezing of an internal liquid-water ocean on Gallisto. Nature, 412, 409-411], heat flows were estimated and the possibility of the existence of a water ocean in Callisto was evaluated. The liquid phase is stable (not freezing) beneath the ice crust, if the heat flow is between 3.3 and 3.7 mW m⁻², which corresponds to the heat flow from radiogenic sources. The thickness of the ice-I crust is 135-150 km, and that of the underlying water layer, 120-180 km. The results of modeling support the hypothesis that Callisto may have an internal liquid-water ocean. The allowed total (maximum) thickness of the outer water-ice shell is up to 270-315 km. Rock-iron core radii, depending on the presence or absence of hydrous silicates, do not exceed 500-700 km, the thickness of an intermediate ice-rock mantle is not less than 1400 km, and its density is in the range of 1960-2500 kg m⁻³. The surface temperature of Callisto is expected to be 100-112 K. The total amount of H₂O in Callisto is found to be 49-55 wt%. The correspondence between the density and moment of inertia values for bulk ice-free Io, rock-iron core of ice-poor Europa, and rock-iron cores of Ganymede and Callisto shows that their bulk compositions may be, in general, similar and may be described by the composition close to a material of the L/LL-type chondrites with the (Fe_tot/Si) weight ratios ranging from 0.9 to 1.3. Planetesimals composed of these types of ordinary chondrites could be considered as analogues of building material for the rock-iron cores of the Galilean satellites. Similarity of bulk composition of the rock-iron cores of the inner and outer satellites implies the absence of iron-silicon fractionation in the protojovian nebula.     

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.     

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.     

Oceans in the icy Galilean satellites of Jupiter?

Equilibrium models of heat transfer by heat conduction and thermal convection show that the three satellites of Jupiter—Europa, Ganymede, and Callisto—may have internal oceans underneath ice shells tens of kilometers to more than a hundred kilometers thick. A wide range of rheology and heat transfer parameter values and present-day heat production rates have been considered. The rheology was cast in terms of a reference viscosity ν₀ calculated at the melting temperature and the rate of change A of viscosity with inverse homologous temperature. The temperature dependence of the thermal conductivity k of ice I has been taken into account by calculating the average conductivity along the temperature profile. Heating rates are based on a chondritic radiogenic heating rate of 4.5 pW kg⁻¹ but have been varied around this value over a wide range. The phase diagrams of H₂O (ice I) and H₂O + 5 wt% NH₃ ice have been considered. The ice I models are worst-case scenarios for the existence of a subsurface liquid water ocean because ice I has the highest possible melting temperature and the highest thermal conductivity of candidate ices and the assumption of equilibrium ignores the contribution to ice shell heating from deep interior cooling. In the context of ice I models, we find that Europa is the satellite most likely to have a subsurface liquid ocean. Even with radiogenic heating alone the ocean is tens of kilometers thick in the nominal model. If tidal heating is invoked, the ocean will be much thicker and the ice shell will be a few tens of kilometers thick. Ganymede and Callisto have frozen their oceans in the nominal ice I models, but since these models represent the worst-case scenario, it is conceivable that these satellites also have oceans at the present time. The most important factor working against the existence of subsurface oceans is contamination of the outer ice shell by rock. Rock increases the density and the pressure gradient and shifts the triple point of ice I to shallower depths where the temperature is likely to be lower then the triple point temperature. According to present knowledge of ice phase diagrams, ammonia produces one of the largest reductions of the melting temperature. If we assume a bulk concentration of 5 wt% ammonia we find that all the satellites have substantial oceans. For a model of Europa heated only by radiogenic decay, the ice shell will be a few tens of kilometers thinner than in the ice I case. The underlying rock mantle will limit the depth of the ocean to 80-100 km. For Ganymede and Callisto, the ice I shell on top of the H₂O-NH₃ ocean will be around 60- to 80-km thick and the oceans may be 200- to 350-km deep. Previous models have suggested that efficient convection in the ice will freeze any existing ocean. The present conclusions are different mainly because they are based on a parameterization of convective heat transport in fluids with strongly temperature dependent viscosity rather than a parameterization derived from constant-viscosity convection models. The present parameterization introduces a conductive stagnant lid at the expense of the thickness of the convecting sublayer, if the latter exists at all. The stagnant lid causes the temperature in the sublayer to be warmer than in a comparable constant-viscosity convecting layer. We have further modified the parameterization to account for the strong increase in homologous temperature, and therefore decrease in viscosity, with depth along an adiabat. This modification causes even thicker stagnant lids and further elevated temperatures in the well-mixed sublayer. It is the stagnant lid and the comparatively large temperature in the sublayer that frustrates ocean freezing.     

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.     

Heat flow and thickness of a convective ice shell on Europa for grain size-dependent rheologies

Although it is mostly accepted that the lower part of the ice shell of Europa is actively convective, there is still much uncertainty about the flow mechanism dominating the rheology of this convective layer, which largely depends on the grain size of the ice. In this work, we examined thermal equilibrium states in a tidally heated and strained convective shell, for two rheologies sensitive to grain size, grain boundary sliding and diffusion creep. If we take a lower limit of 70 mW m⁻² for the surface heat flow, according to some geological features observed, the ice grain size should be less than 2 or 0.2 mm for grain boundary sliding or diffusion creep respectively. If in addition the thickness of the ice shell is constrained to a few tens of kilometers and it is assumed that the thickness of the convective layer is related to lenticulae spacing, then grain sizes between 0.2 and 2 mm for grain boundary sliding, and between 0.1 and 0.2 mm for diffusion creep are obtained. Also, local convective layer thicknesses deduced from lenticulae spacing are more similar to those here derived for grain boundary sliding. Our results thus favor grain boundary sliding as the dominant rheology for the water ice in Europa's convective layer, since this flow mechanism is able to satisfy the imposed constraints for a wider range of grain sizes.     

Fracture penetration in planetary ice shells

The proposed past eruption of liquid water on Europa and ongoing eruption of water vapor and ice on Enceladus have led to discussion about the feasibility of cracking a planetary ice shell. We use a boundary element method to model crack penetration in an ice shell subjected to tension and hydrostatic compression. We consider the presence of a region at the base of the ice shell in which the far-field extensional stresses vanish due to viscoelastic relaxation, impeding the penetration of fractures towards a subsurface ocean. The maximum extent of fracture penetration can be limited by hydrostatic pressure or by the presence of the unstressed basal layer, depending on its thickness. Our results indicate that Europa's ice shell is likely to be cracked under 1-3 MPa tension only if it is ?2.5 km thick. Enceladus' ice shell may be completely cracked if it is capable of supporting ∼1-3 MPa tension and is less than 25 km thick.     

Internal structures of the Galilean satellites

New models for the interiors of Io, Ganymede, and Callisto are proposed. The model of Io consists of a thin, high-rigidity outer layer separated from a solid interior by a thin, molten or partially molten shell. The modulus of rigidity of the outer layer must be at least 100 times larger than that of the underlying partially molten shell. These layers have thicknesses of order 100 km or less. The near-surface partially molten layer was most likely produced early in Io's history as a consequence of accretional heating; enhanced tidal heating in the outer rigid layer has kept the underlying region partially molten to the present day. The model of Ganymede consists of an ice outer layer, a shell of undifferentiated, primordial ice-silicate mixture, and a rock core. Accretional heating is responsible for melting the ice in the outer layers of Ganymede's initially homogeneous ice-silicate interior. Most of the rock in this outer layer accumulates in a shell on top of Ganymede's early cold and rigid central region; the water in the outer layer quickly refreezes. Heating of the undifferentiated region by the decay of radioactive elements in the silicate fraction would gradually warm it and reduce its viscosity. The rock layer would become gravitationally unstable and sink through the undifferentiated materials to form a rock core. Callisto's heavily cratered surface strongly suggests that relatively little, if any, ice-rock differentiation has occured in its interior.     

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.     

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.     

The heat flow of Europa

The heat flow from Europa has profound implications for ice shell thickness and structure, as well as for the existence of an internal ocean, which is strongly suggested by magnetic data. The brittle-ductile transition depth and the effective elastic thickness of the lithosphere are here used to perform heat flow estimations for Europa. Results give preferred heat flow values (for a typical geological strain rate of ¹⁰⁻¹⁵ s⁻¹) of 70-110 mW m⁻² for a brittle-ductile transition 2 km deep (the usually accepted upper limit for the brittle-ductile transition depth in the ice shell of Europa), 24-35 mW m⁻² for an effective elastic thickness of 2.9 km supporting a plateau near the Cilix impact crater, and >130 mW m⁻² for effective elastic thicknesses of ?0.4 km proposed for the lithosphere loaded by ridges and domes. These values are clearly higher than those produced by radiogenic heating, thus implying an important role for tidal heating. The ?19-25 km thick ice shell proposed from the analysis of size and depth of impact structures suggests a heat flow of ?30-45 mW m⁻² reaching the ice shell base, which in turn would imply an important contribution to the heat flow from tidal heating within the ice shell. Tidally heated convection in the ice shell could be capable to supply ∼100 mW m⁻² for superplastic flow, and, at the Cilix crater region, ∼35-50 mW m⁻² for dislocation creep, which suggests local variations in the dominant flow mechanism for convection. The very high heat flows maybe related to ridges and domes could be originated by preferential heating at special settings.