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.     

Non-Newtonian topographic relaxation on Europa

Models of topographic support on Europa by lateral shell thickness variations have previously assumed a Newtonian ice viscosity. Here I show that using a more realistic stress-dependent viscosity gives relaxation times which can be significantly different. Topography of wavelength 100 km cannot be supported by lateral shell thickness variations for ∼50 Myr, unless the shell thickness is <10 km or the ice grain size >10 mm. Shorter wavelength topography would require even thinner shells, but may be supported elastically. Global-scale variations in shell thickness, however, can be supported for geological timescales if the shell thickness is O(10 km).     

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.     

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.     

Modeling stresses on satellites due to nonsynchronous rotation and orbital eccentricity using gravitational potential theory

The tidal stress at the surface of a satellite is derived from the gravitational potential of the satellite's parent planet, assuming that the satellite is fully differentiated into a silicate core, a global subsurface ocean, and a decoupled, viscoelastic lithospheric shell. We consider two types of time variability for the tidal force acting on the shell: one caused by the satellite's eccentric orbit within the planet's gravitational field (diurnal tides), and one due to nonsynchronous rotation (NSR) of the shell relative to the satellite's core, which is presumed to be tidally locked. In calculating surface stresses, this method allows the Love numbers h and ?, describing the satellite's tidal response, to be specified independently; it allows the use of frequency-dependent viscoelastic rheologies (e.g. a Maxwell solid); and its mathematical form is amenable to the inclusion of stresses due to individual tides. The lithosphere can respond to NSR forcing either viscously or elastically depending on the value of the parameter , where μ and η are the shear modulus and viscosity of the shell respectively, and ω is the NSR forcing frequency. Δ is proportional to the ratio of the forcing period to the viscous relaxation time. When Δ?1 the response is nearly fluid; when Δ?1 it is nearly elastic. In the elastic case, tensile stresses due to NSR on Europa can be as large as ∼3.3 MPa, which dominate the ∼50 kPa stresses predicted to result from Europa's diurnal tides. The faster the viscous relaxation the smaller the NSR stresses, such that diurnal stresses dominate when Δ?100. Given the uncertainty in current estimates of the NSR period and of the viscosity of Europa's ice shell, it is unclear which tide should be dominant. For Europa, tidal stresses are relatively insensitive both to the rheological structure beneath the ice layer and to the thickness of the icy shell. The phase shift between the tidal potential and the resulting stresses increases with Δ. This shift can displace the NSR stresses longitudinally by as much as 45° in the direction opposite of the satellite's rotation.     

Elastic ice shells of synchronous moons: Implications for cracks on Europa and non-synchronous rotation of Titan

A number of synchronous moons are thought to harbor water oceans beneath their outer ice shells. A subsurface ocean frictionally decouples the shell from the interior. This has led to proposals that a weak tidal or atmospheric torque might cause the shell to rotate differentially with respect to the synchronously rotating interior. Applications along these lines have been made to Europa and Titan. However, the shell is coupled to the ocean by an elastic torque. As a result of centrifugal and tidal forces, the ocean would assume an ellipsoidal shape with its long axis aligned toward the parent planet. Any displacement of the shell away from its equilibrium position would induce strains thereby increasing its elastic energy and giving rise to an elastic restoring torque. In the investigation reported on here, the elastic torque is compared with the tidal torque acting on Europa and the atmospheric torque acting on Titan.Regarding Europa, it is shown that the tidal torque is far too weak to produce stresses that could fracture the ice shell, thus refuting an idea that has been widely advocated. Instead, it is suggested that the cracks arise from time-dependent stresses due to non-hydrostatic gravity anomalies from tidally driven, episodic convection in the satellite’s interior.Two years of Cassini RADAR observations of Titan’s surface have been interpreted as implying an angular displacement of ∼0.24° relative to synchronous rotation. Compatibility of the amplitude and phase of the observed non-synchronous rotation with estimates of the atmospheric torque requires that Titan’s shell be decoupled from its interior. We find that the elastic torque balances the seasonal atmospheric torque at an angular displacement ?0.05°, effectively coupling the shell to the interior. Moreover, if Titan’s surface were spinning faster than synchronous, the tidal torque tending to restore synchronous rotation would almost certainly be larger than the atmospheric torque. There must either be a problem with the interpretation of the radar observations, or with our basic understanding of Titan’s atmosphere and/or interior.     

Probing Europa's interior with natural sound sources

Europa's interior structure may be determined by relatively simple and robust seismo-acoustic echo sounding techniques. The strategy is to use ice cracking events or impacts that are hypothesized to occur regularly on Europa's surface as sources of opportunity. A single passive geophone on Europa's surface may then be used to estimate the thickness of its ice shell and the depth of its ocean by measuring the travel time of seismo-acoustic reflections from the corresponding internal strata. Quantitative analysis is presented with full-field seismo-acoustic modeling of the Europan environment. This includes models for Europan ambient noise and conditions on signal-to-noise ratio necessary for the proposed technique to be feasible. The possibility of determining Europa's ice layer thickness by surface wave and modal analysis with a single geophone is also investigated.     

Energy Distributions for Desorption of Sodium and Potassium from Ice: The Na/K Ratio at Europa

The sputtering and decomposition of the surface of Europa by fast ions and electrons lead to the production of an atomosphere containing sodium and potassium atoms. Here time-of-flight energy distributions are measured for Na and K sputtered from a vapor-deposited ice by 200-eV electrons. These data are then used in a Monte Carlo simulation for alkalis in Europa's atmosphere. Na/K ratios versus distance from Europa are calculated and compared to the recent observations in the range 6 to 18 Europan radii from the surface. Normalizing to the observations, the Na/K ratio for the loss rates is ∼27 and the ratio for the average surface source rates is ∼20. These ratios are very different from the Na/K ratio at Io and are larger than the Na/K ratio suggested for Europa's putative subsurface ocean, consistent with fractionation on freezing and upwelling of ocean material.     

Flexure of Europa's lithosphere due to ridge-loading

Using photoclinometry, topographic profiles across europan ridges have been produced. These profiles allow the identification of bulges in the terrain adjacent to the ridges. The bulges are assumed to have been produced by flexure of the elastic lithosphere due to the load of the ridges, which lie along cracks in the crust. The distance from the crack to these “fore-bulges” depends on the thickness of the elastic plate being flexed. Based on a survey of ridges in Galileo images with resolution <300 m/pixel, the thickness of the elastic lithosphere has been determined by this method at a wide variety of sites along the leading and trailing hemispheres of Europa. The average thickness is about 200 m. The elastic lithosphere underneath smooth dilational bands tends to be thicker than plains morphology, an effect that is pronounced at Thynia Linea and Astypalaea Linea. Among the ridges investigated here, more recent loading correlates with a thicker elastic lithosphere, which may either reflect an intrinsically thicker layer, or less viscous relaxation over the shorter time period.     

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 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.     

Strength of the lithosphere of the Galilean satellites

Several approaches have been used to estimate the ice shell thickness on Callisto, Ganymede, and Europa. Here we develop a method for placing a strict lower bound on the thickness of the strong part of the shell (lithosphere) using measurements of topography. The minimal assumptions are that the strength of faults in the brittle lithosphere is controlled by lithostatic pressure according to Byerlee's law and the shell has relatively uniform density and thickness. Under these conditions, the topography of the ice provides a direct measure of the bending moment in the lithosphere. This topographic bending moment must be less than the saturation bending moment of the yield strength envelope derived from Byerlee's law. The model predicts that the topographic amplitude spectrum decreases as the square of the topographic wavelength. This explains why Europa is rugged at shorter wavelengths (∼10 km) but extremely smooth, and perhaps conforming to an equipotential surface, at longer wavelengths (>100 km). Previously compiled data on impact crater depth and diameter [Schenk, P.M., 2002. Nature 417, 419-421] on Europa show good agreement with the spectral decrease predicted by the model and require a lithosphere thicker than 2.5 km. A more realistic model, including a ductile lower lithosphere, requires a thickness greater than 3.5 km. Future measurements of topography in the 10-100 km wavelength band will provide tight constraints on lithospheric strength.     

Numerical modeling of endogenic thermal anomalies on Europa

A variety of recent resurfacing features have been observed on Europa, which may produce thermal anomalies detectable by a future mission. However, the likelihood of such a detection depends on their size and lifetimes. The results of this numerical study suggest that the lifetime of a thermal anomaly associated with the emplacement of 100 m of water onto the surface of Europa is several hundred years, and ∼10 years for 10 m of water. If warm ice is emplaced on the surface instead of liquid water, these lifetimes decrease by up to a factor of two. Exploration of model parameters indicates that a thin insulating surface layer can double thermal anomaly lifetimes, anomalies emplaced at a latitude of 80° can remain detectable nearly a factor of two longer than those at equatorial latitudes, and anomalies on the night side can remain detectable for up to ∼20% longer than those on the day side. High temperatures are very short-lived as the surface ice cools very rapidly to below 200 K due to sublimation cooling. Assuming steady-state resurfacing, the number of detectable thermal anomalies associated with the emplacement of 100 m of water would be on the order of 10 if the typical resurfacing area is 15 km². If recent resurfacing is dominated by chaos regions with typical areas of 100 to 1000 km² and lifetimes of 1000 to 4000 years, the number of detectable thermal anomalies would be on the order of 1 to 10.     

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.     

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.     

Nonsynchronous Rotation Evidence and Fracture History in the Bright Plains Region, Europa

A geologic map for the Bright Plains in the Conamara Chaos region of Europa is presented and is used to unravel a detailed fracture sequence using cross-cutting relationships and fracture mechanics principles. Fracture orientations in the Bright Plains region rotated with time, consistently in a clockwise sense. This conclusion agrees with the observations of other researchers' northern Europan hemisphere investigations and points strongly toward the fracture sequence being controlled by the effect of nonsynchronous rotation, whereby the outer ice crust of Europa rotates slightly faster than the satellite's interior. This is convincing evidence that Europa's crust has been decoupled from the interior, possibly due to the presence of a liquid ocean beneath the crust.Tidal stresses induced in the ice crust by the combined effects of nonsynchronous rotation and diurnal tidal flexing can be calculated using the assumption that the crust behaves elastically over relatively short time scales (i.e., no viscous relaxation of stresses). The fracture orientations in the Bright Plains area were compared to a global scale tidal stress field to determine the longitudes at which each fracture set developed. The fracture sequence points strongly to the Bright Plains region of the crust having rotated at least 720° (and perhaps up to 900°) with respect to the satellite's interior during the visible fracture history. This amount exceeds previously published estimates of nonsynchronous rotation. The orientations of the most recent surface fractures are incompatible with the current state of stress in the Bright Plains region, implying a period of a few thousand years since the most recent fracturing events based on existing nonsynchronous rotation rate estimates.     

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.     

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.