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.     

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.     

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.     

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.     

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.     

Coupled convection and tidal dissipation in Europa’s ice shell

We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 10¹³ Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other 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.     

On convection in ice I shells of outer Solar System bodies, with detailed application to Callisto

It has been argued that the dominant non-Newtonian creep mechanisms of water ice make the ice shell above Callisto's ocean, and by inference all radiogenically heated ice I shells in the outer Solar System, stable against solid-state convective overturn. Conductive heat transport and internal melting (oceans) are therefore predicted to be, or have been, widespread among midsize and larger icy satellites and Kuiper Belt objects. Alternatively, at low stresses (where non-Newtonian viscosities can be arbitrarily large), convective instabilities may arise in the diffusional creep regime for arbitrarily small temperature perturbations. For Callisto, ice viscosities are low enough that convection is expected over most of geologic time above the internal liquid layer for plausible ice grain sizes (?a few mm); the alternative for early Callisto, a conducting shell over a very deep ocean (>450 km), is not compatible with Callisto's present partially differentiated state. Moreover, if convection is occurring today, the stagnant lid would be quite thick (∼100 km) and compatible with the lack of active geology. Nevertheless, Callisto's steady-state heat flow may have fallen below the convective minimum for its ice I shell late in Solar System history. In this case convection ends, the ice shell melts back at its base, and the internal ocean widens considerably. The presence of such an ocean, of order 200 km thick, is compatible with Callisto's moment-of-inertia, but its formation would have caused an ∼0.25% radial expansion. The tectonic effects of such a late, slow expansion are not observed, so convection likely persists in Callisto, possibly subcritically. Ganymede, due to its greater size, rock fraction and full differentiation, has a substantially higher heat flow than Callisto and has not reached this tectonic end state. Titan, if differentiated, and Triton should be more similar to Ganymede in this regard. Pluto, like Callisto, may be near the tipping point for convective shutdown, but uncertainties in its size and rock fraction prevent a more definitive assessment.     

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.     

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.     

Onset of Convection,Heat Flow and Thickness of the Europa‘s ice Shell

The observational evidence given by Galileo spacecraft about Europa supports an icy rigid layer of several kilometers over another ductile layer of ice in convection, which floats over an internal ocean of liquid water. Before the onset of convection, heat is transmitted into the crust by conduction. The heat flow analysis in the potentially convective layer gives values higher than those obtained previously by tidal dissipation models, and suggests that the ice may be limited to a thin layer of ∼4 km total thickness. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

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.     

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.     

Effects of plasticity on convection in an ice shell: Implications for Europa

Europa's surface exhibits numerous pits, uplifts, and disrupted chaos terrains that have been suggested to result from convection in the ice shell. To test this hypothesis, we present numerical simulations of convection in an ice shell including the effects of plasticity, which provides a simple continuum representation for brittle or semibrittle deformation along discrete fractures. Plastic deformation occurs when stresses reach a specified yield stress; at lower stresses, the fluid flow follows a Newtonian, temperature-dependent viscosity. Four distinct modes of behavior can occur. For yield stresses exceeding ∼1 bar, plastic effects are negligible and stagnant-lid convection, with no surface motion and minimal topography, results. At intermediate yield stresses, a stagnant lid forms but deforms plastically, leading to surface velocities up to several millimeters per year. Slightly smaller yield stresses allow episodic, catastrophic overturns of the upper conductive lid, with (transient) stagnant lids forming in between overturn events. The smallest yield stresses allow continual recycling of the upper lid, with simultaneous, gradual ascent of warm ice to the surface and descent of cold, near-surface ice into the interior. The exact yield stresses over which each regime occurs depend on the ice-shell thickness, melting-temperature viscosity, and activation energy for viscous creep. To form hummocky matrix and translate chaos plates by several kilometers, substantial surface strain must accompany chaos formation, and the three plasticity-dominated convection modes described here can provide such deformation. Our simulations suggest that, if yield stresses of ∼0.2-1 bar are relevant to Europa, then convection in Europa's ice shell can produce chaos-like structures at the surface. However, our simulations have difficulty explaining Europa's numerous pits and uplifts. When plasticity forces the upper lid to participate in the convection, dynamic topography of ∼50-100-m amplitude results, but the topographic structures generally have diameters of 30-100 km, an order of magnitude wider than typical pits and uplifts. None of our simulations produced isolated pits or uplifts of any diameter.     

On the origins of band topography, Europa

We use stereo-derived topography of extensional bands on Europa to show that these features can be elevated by 100-150 m with respect to the surroundings, and that the positive topography sometimes extends beyond the band margins. Lateral variations in shell thickness cannot maintain the observed topography for timescales greater than ∼0.1 Myr. Lateral density variations can maintain the observed topography indefinitely; mean density contrasts of 5 and 50 kg m⁻³ are required for shell thicknesses of 20 and 2 km, respectively. Density variations caused by temperature contrasts require either present-day heating or that bands are young features (<1 Myr old). Stratigraphic analyses suggest that these mechanisms are unlikely. The observation that bands form from ridges may be explained by an episode of shear-heating on ridges weakening the ridge area, and leading to strain localization during extension. Fracture porosity is likely to persist over Myr timescales in the top one-third to one-quarter of the conductive part of the ice shell. Lateral variations in this porosity (of order 20%) are the most likely mechanism for producing band topography if the ice shell is thin (≈2 km); porosity variations of 2% or less are required if the shell is thicker (≈20 km). If the ice shell is thick, lateral variations in salt content are a more likely mechanism. Warm ice will tend to lose dense, low-melting temperature phases and be buoyant relative to colder, salt-rich ice. Thus, lateral density variations will arise naturally if bands have been the sites of either localized heating or upwelling of warm ice during extension.     

A model for the interior structure, evolution, and differentiation of Callisto

The recently measured dimensionless moment of inertia (MoI) factor for Callisto of 0.3549±0.0042 (Anderson et al., 2001, Icarus, 153, 157-161) poses a problem: its value cannot be explained by a model in which Callisto is completely differentiated into an ice shell above a rock shell and an iron core such as its neighboring satellite Ganymede nor can it be explained by a model of a homogeneous, undifferentiated ice-rock satellite. We show that Callisto may be incompletely differentiated into an outer ice-rock shell in which the volumetric rock concentration is close to the primordial one at the surface and decreases approximately linearly with depth, an ice mantle mostly depleted of rock, and an about 1800 km rock-ice core in which the rock concentration is close to the close-packing limit. The ice-rock shell thickness depends on uncertain rheology parameters and the heat flow and can be roughly 50 to 150 km thick. We show that if Callisto accreted from a mix of metal bearing rock and ice and if the average size of the rocks was of the order of meters to tens of meters, then Callisto may have experienced a gradual, but still incomplete unmixing of the two components. An ocean in Callisto at a depth of 100-200 km is difficult to obtain if the ice is pure H₂O and if the ice-rock lithosphere is 100 km or more thick; a water ocean is more plausible for ice contaminated by ammonia, methane or salts; or for pure H₂O at a depth of 400-600 km.     

Long-Term Evolution of Objects in the Kuiper Belt Zone—Effects of Insolation and Radiogenic Heating

The Kuiper Belt zone is unique insofar as the major heat sources of objects a few tens of kilometers in size—solar radiation on the one hand and radioactive decay on the other—have comparable power. This leads to unique evolutionary patterns, with heat waves propagating inward from the irradiated surface and outward from the radioactively heated interior. A major radioactive source that is considered in this study is ²⁶Al. The long-term evolution of several models with characteristics typical of Kuiper Belt objects is followed by means of a 1-D numerical code that solves the heat and mass balance equations on a spherically symmetric grid. The free parameters considered are radius (10-500 km), heliocentric distance (30-120 AU), and initial ²⁶Al content (0-5×10⁻⁸ by mass). The initial composition assumed is a porous mixture of ices (H₂O, CO, and CO₂) and dust. Gases released in the interior are allowed to escape to the surface. It is shown that, depending on parameters, the interior may reach quite high temperatures (up to 180 K). The models suggest that Kuiper Belt objects are likely to lose the ices of very volatile species during early evolution; ices of less volatile species are retained in a surface layer, about 1 km thick. The models indicate that the amorphous ice crystallizes in the interior, and hence some objects may also lose part of the volatiles trapped in amorphous ice. Generally, the outer layers are far less affected than the inner part, resulting in a stratified composition and altered porosity distribution. These changes in structure and composition should have significant consequences for the short-period comets, which are believed to be descendants of Kuiper Belt objects.     

The global shape of Europa: Constraints on lateral shell thickness variations

The global shape of Europa is controlled by tidal and rotational potentials and possibly by lateral variations in ice shell thickness. We use limb profiles from four Galileo images to determine the best-fit hydrostatic shape, yielding a mean radius of 1560.8±0.3 km and a radius difference a−c of 3.0±0.9 km, consistent with previous determinations and inferences from gravity observations. Adding long-wavelength topography due to proposed lateral variations in shell thickness results in poorer fits to the limb profiles. We conclude that lateral shell thickness variations and long-wavelength isostatically supported topographic variations do not exceed 7 and 0.7 km, respectively. For the range of rheologies investigated (basal viscosities from 10¹⁴ to ) the maximum permissible (conductive) shell thickness is 35 km. The relative uniformity of Europa's shell thickness is due to either a heat flux from the silicate interior, lateral ice flow at the base of the shell, or convection within the shell.     

The geothermal gradient of Io: Consequences for lithosphere structure and volcanic eruptive activity

We solve numerically the equations describing the transfer of heat through the lithosphere of Io by a mixture of conduction and volcanic advection as proposed by O’Reilly and Davies (O’Reilly, T.C., Davies, G.F. [1981]. Geophys. Res. Lett. 8, 313-316), removing the requirement that average material properties must be used. As expected, the dominance of advective heat transfer by volcanic eruptions means that Io’s geothermal gradient well away from volcanic centres is very small, of order 1 K km⁻¹. This result is independent of any reasonable assumptions about the radiogenic heating rate in the lithosphere. The lithosphere temperature does not increase greatly above the surface temperature until the base of the lithosphere is approached, except in limited areas around shallow magma bodies. As a consequence, solid volatile sulphur compounds mobilized by volcanic processes and re-deposited on the surface of Io commonly remain solid until they reach great depths as they are progressively buried by ongoing activity. For current estimates of the volcanic heat transfer rate, melting of SO₂ does not begin until a depth of ∼20 km and sulphur remains solid to a depth of ∼26 km in a 30 km thick lithosphere. Rising magmas can incorporate fluids from these deep sulphur compound aquifers, and we quantify the major influence that this can have on the bulk density of the magma and hence the resulting possible intrusion and eruption styles.