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

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.     

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.     

Measurements of the elastic thickness under ancient lunar terrain

The part of a planet's gravity that is coherent with its topography provides information about the deflection of its crust after loading, and hence the resistance of its lithosphere to bending at the time the load was emplaced. We used observed line of sight accelerations from Lunar Prospector, together with the accelerations we would have expected if anomalies in the gravity field were only caused by topography, to estimate the admittance and coherence between topography and gravity over several regions of the lunar nearside. We then compared our estimates to what we would expect if the lithosphere behaved as a floating elastic plate or shell, assuming a linear relationship between topography and gravity. We found in the region surrounding the southern highland crater Clavius that the data can be modeled using a thin plate with Te=12±5 km and uncorrelated loads at the surface and base of the upper crust. A spherical shell model with surface loading is less satisfactory: to fit the admittance adequately requires topography with wavelengths over 400 km to be formed when Te≈ 1 km and the remainder when Te≈7.5 km. By contrast, the apparent lack of compensation around the youngest giant impact basins requires a plate with Te>80 km or a shell with Te>25 km. Our results indicate the thickness of the lunar lithosphere increased from ≈12 km in the pre-Nectarian to >25 km in the Nectarian.     

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.     

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.     

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.     

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.     

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.     

The great thickness debate: Ice shell thickness models for Europa and comparisons with estimates based on flexure at ridges

Estimates of the thickness of the ice shell of Europa range from <1 to >30 km. The higher values are generally assumed to be estimates of the entire ice shell thickness, which may include a lower ductile layer of ice, whereas many of the smaller thickness estimates are based on analyses that only consider that portion of the ice layer that behaves elastically at a particular strain rate. One example of the latter is flexure analysis, in which the elastic ice layer is modeled as a plate or sphere that is flexed under the weight of a surface load. We present calculations based on flexure analysis in which we model the elastic ice layer as flexing under a line-load caused by ridges. We use precisely located, parallel flanking cracks as indicators of the location of greatest tensile stress induced by flexure. Our elastic thickness results are spatially variable: ∼500-2200 m (two sites) and ∼200-1000 m (one site). Thorough analysis of Europan flexure studies performed by various researchers shows that the type of model selected causes the greatest variability in the thickness results, followed by the choice of Young's modulus, which is poorly constrained for the Europan ice shell. Comparing our results to those of previously published flexure analyses for Europa, we infer spatial variability in the elastic ice thickness (at the time of load emplacement), with smooth bands having the thinnest elastic ice thickness of all areas studied. Because analysis of flexure-induced fracturing can only reveal the elastic thickness at the time of load emplacement, calculated thickness variability between features having different ages may also reflect a temporal variability in the thickness of Europa's ice shell.     

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.     

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.     

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.     

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.     

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.     

Deciphering plume–lithosphere interactions beneath Europe from topographic signatures

When a mantle plume rises and impinges on the base of the lithosphere, it expectably produces variations in surface topography. Taking into consideration a realistic mantle rheology, plume ascent rates can reach tens to hundreds of metres per year, whereupon the impingement of the plume head at the base of the lithosphere can be considered as an “impact". Recent numerical experiments based on tectonically realistic formulations for the lithosphere and a representative mantle rheology, have shown that plume-induced undulations exhibit temporal successions of uplift and subsidence at various wavelengths. From spectral (Fourier) analyses of the undulations would appear that two groups of wavelengths (200–400 km and 60–100 km) predominate. Interestingly, a spectral analysis of Europe's topography also reveals two dominant groups. In the present study, we have used a spectral analysis with a wavelet formulation in order to discriminate between tectonically-induced undulations (uni-directional deformation) and plume-induced undulations (omni-directional deformation). The European lithosphere is well-suited for this approach since it has been suggested that two mantle plumes (the Massif Central and the Eifel area) underlie Western Europe, where Alpine compression has folded the lithosphere over several hundreds of kilometres. The wavelet analysis of Europe's surface topography confirms that the energy distribution of the topographic undulations outside the two main volcanic provinces is homogeneous, thus contrasting with the coexistence of both large-scale and medium-scale high-energy features that are obtained for the Massif Central and Eifel areas. Similar signatures are also found beneath the northern Sudetes area. The wavelet approach dedicated to the analysis of plume-related topographic signatures needs, however, detailed theoretical grounds and would probably benefit from two-dimensional analyses.     

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.     

Models of the Internal Structure of Callisto

Models of the internal structure of Callisto were constructed and the extent of its differentiation was determined based on geophysical information from the Galileo spacecraft (the mass, the radius, the mean density, and the moment of inertia), geochemical data (the chemical composition of meteorites), and the equations of state of water, ices, and meteoritic material. The thickness and the phase state of the water-ice shell were defined as well as the ice concentrations in the rock-ice mantle and the bulk concentration of H₂O. The constraints on the density distribution in the mantle and the size of the rock-iron core were derived. We considered models of the internal structure of Callisto in which the presence of a continuous ice shell was assumed (models without ocean) and models with an internal ocean. We demonstrated that it is possible to apply three-layer models with an icy shell up to 320 km in thickness and a rock-iron core in different combinations with a rock-ice mantle. These models do not reject a two-layer structure of Callisto (an ice lithosphere plus a rock-ice mantle or a rock-ice mantle plus a rock-iron core) and a one-layer model of the satellite composed only of a rock-ice mantle with an ice concentration that is variable in depth. Taking into account the chemically bound water, the bulk content of H₂O in the satellite is found to be 49–55 wt %. For the model with an internal ocean, the geophysically allowed thickness of the water-ice shell of Callisto was estimated to be 270–315 km with thicknesses of the icy crust and the underlying water layer of 135–150 and 120–180 km, respectively. The results of reconstruction of the composition and structure of the regular satellites of Jupiter allow us to conclude that they were possibly formed from material whose composition was close to ordinary L/LL chondrites at relatively low temperatures, lower than the temperature of evaporation of iron and Fe-Mg silicates.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 321–341.Original Russian Text Copyright © 2005 by Kuskov, Kronrod.     

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.     

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.