Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Solid tidal friction above a liquid water reservoir as the origin of the south pole hotspot on Enceladus

Earth, Jupiter's moon Io and Saturn's tiny moon Enceladus are the only solid objects in the Solar System to be sufficiently geologically active for their internal heat to be detected by remote sensing. Interestingly, the endogenic activity on Enceladus is only located on a specific region at the south pole, from which jets of water vapor and ice particles have been observed [Spencer, J.R., and 9 colleagues, 2006. Science 311, 1401-1405; Porco, C.C., and 24 colleagues, 2006. Science 311, 1393-1401]. The current polar location of the thermal anomaly can possibly be explained by diapir-induced reorientation of the satellite [Nimmo, F., Pappalardo, R.T., 2006. Nature 441, 614-616], but the thermal anomaly triggering and the heat power required to sustain it over geological timescales remain problematic. Using a three-dimensional viscoelastic numerical model simulating the response of Enceladus to tidal forcing, we explore the effect of a low-viscosity anomaly in the ice shell, localized to the south polar region, on the tidal dissipation patterns. We demonstrate that only interior models with a liquid water layer at depth can explain the observed magnitude of dissipation rate and its particular location at the south pole. Moreover, we show that tidally-induced heat must be generated over a relatively broad region in the southern hemisphere, and it is then transferred toward the south pole where it is episodically released during relatively short resurfacing events. As large tidal dissipation and internal melting cannot be induced in the south polar region in the absence of a pre-existing liquid decoupling layer, we propose that liquid water must have been present in the interior for a very long period of time, and possibly since the satellite formation. Owing to the orbital equilibrium requirement [Meyer, J., Wisdom, J., 2007. Icarus 188, 535-539], sustaining some liquid water at depth is impossible if heat is continuously emitted at a rate of 4-8 GW at the south pole. Based on that requirement, we propose that the current thermal emission is not in equilibrium with the heat production, and that the thermal emission rate is abnormally high at present time. Alternatively, continuous dissipation at a rate of 1-2 GW in the ice shell at the south pole should be sufficient to induce internal melting and it could sustain a layer of liquid water at depth over geologic timescales.     

A model for the temperature-dependence of tidal dissipation in convective plumes on icy satellites: Implications for Europa and Enceladus

To explain the formation of surface features on Europa, Enceladus, and other satellites, many authors have postulated the spatial localization of tidal heating within convective plumes. However, the concept that enhanced tidal heating can occur within a convective plume has not been rigorously tested. Most models of this phenomenon adopt a tidal heating with a temperature-dependence derived for an incompressible, homogeneous (zero-dimensional) Maxwell material, but it is unclear whether this formulation is relevant to the heterogeneous situation of a warm plume surrounded by cold ice. To determine whether concentrated dissipation can occur in convective plumes, we develop a two-dimensional model to compute the volumetric dissipation rate for an idealized, vertically oriented, isolated convective plume obeying a Maxwellian viscoelastic compressible rheology. We apply the model to the Europa and Enceladus ice shells, and we investigate the consequences for partial melting and resurfacing processes on these bodies. We find that the tidal heating is strongly temperature dependent in a convective ice plume and could produce elevated temperatures and local partial melting in the ice shells of Europa and Enceladus. Our calculation provides the first quantitative verification of the hypothesis by Sotin et al. [Sotin, C., Head, J.W., Tobie, G., 2002. Geophys. Res. Lett. 29. 74-1] and others that the tidal dissipation rate is a strong function of temperature inside a convective plume. On Europa, such localized heating could help allow the formation of domes and chaos terrains by convection. On Enceladus, localized tidal heating in a thermal plume could explain the concentrated activity at the south pole and its associated heat transport of 2-7 GW.     

Viscosity of the lithosphere of Enceladus

High resolution Voyager II images of Enceladus reveal that some regions on its surface are highly cratered; the most heavily cratered surfaces probably date back to a period of heavy bombardment. The forms of many of the craters on Enceladus are similar to those of fresh lunar craters, but many of the craters are much shallower in depth, and the floors of some craters are bowed up. The flattering of craters and bowing up of the floors are indicative of viscous relaxation of the topography. Analysis of the forms of the flattened craters suggests that the viscosity at the top of the lithosphere, in the most heavily cratered regions, is between 10²⁴ and 10²⁵ P. The exact time scale for the collapse of the craters is not known, but probably was between 100 my and 4 gy. The flattened craters are located in distinct zones that are adjacent to zones, of similar age, where craters have not flattened. The zones where flattened craters occur possibly are regions in which the heat flow was (or is) higher than that in the adjacent terrains. Because the temperature at the top of the lithosphere of Enceladus would be less than or equal to that of Ganymede and Callisto, if it is covered by a thick regolith, and because the required viscosity on Enceladus is one to two orders of magnitude less than that for Ganymede and Callisto, it can be concluded that the lithospheric material on Enceladus is different from that of Ganymede and Callisto. Enceladus probably has a mixture of ammonia ice and water ice in the lithosphere, whereas the lithospheres of Ganymede and Callisto are composed primarily of water ice.     

Enceladus: Present internal structure and differentiation by early and long-term radiogenic heating

Pre-Cassini images of Saturn's small icy moon Enceladus provided the first indication that this satellite has undergone extensive resurfacing and tectonism. Data returned by the Cassini spacecraft have proven Enceladus to be one of the most geologically dynamic bodies in the Solar System. Given that the diameter of Enceladus is only about 500 km, this is a surprising discovery and has made Enceladus an object of much interest. Determining Enceladus' interior structure is key to understanding its current activity. Here we use the mean density of Enceladus (as determined by the Cassini mission to Saturn), Cassini observations of endogenic activity on Enceladus, and numerical simulations of Enceladus' thermal evolution to infer that this satellite is most likely a differentiated body with a large rock-metal core of radius about 150 to 170 km surrounded by a liquid water-ice shell. With a silicate mass fraction of 50% or more, long-term radiogenic heating alone might melt most of the ice in a homogeneous Enceladus after about 500 Myr assuming an initial accretion temperature of about 200 K, no subsolidus convection of the ice, and either a surface temperature higher than at present or a porous, insulating surface. Short-lived radioactivity, e.g., the decay of ²⁶Al, would melt all of the ice and differentiate Enceladus within a few million years of accretion assuming formation of Enceladus at a propitious time prior to the decay of ²⁶Al. Long-lived radioactivity facilitates tidal heating as a source of energy for differentiation by warming the ice in Enceladus so that tidal deformation can become effective. This could explain the difference between Enceladus and Mimas. Mimas, with only a small rock fraction, has experienced relatively little long-term radiogenic heating; it has remained cold and stiff and less susceptible to tidal heating despite its proximity to Saturn and larger eccentricity than Enceladus. It is shown that the shape of Enceladus is not that of a body in hydrostatic equilibrium at its present orbital location and rotation rate. The present shape could be an equilibrium shape corresponding to a time when Enceladus was closer to Saturn and spinning more rapidly, or more likely, to a time when Enceladus was spinning more rapidly at its present orbital location. A liquid water layer on Enceladus is a possible source for the plume in the south polar region assuming the survivability of such a layer to the present. These results could place Enceladus in a category similar to the large satellites of Jupiter, with the core having a rock-metal composition similar to Io, and with a deep overlying ice shell similar to Europa and Ganymede. Indeed, the moment of inertia factor of a differentiated Enceladus, C/MR², could be as small as that of Ganymede, about 0.31.     

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.     

Tidal heating and the long-term stability of a subsurface ocean on Enceladus

Tidal dissipation has been suggested as the heat source for the south polar thermal anomaly on Enceladus. We find that under present-day conditions and assuming Maxwellian behavior, tidal dissipation is negligible in the silicate core. Dissipation may be significant in the ice shell if the shell is decoupled from the silicate core by a subsurface ocean. We have run a series of self-consistent convection and conduction models in 2D axisymmetric and 3D spherical geometry in which we include the spatially-variable tidal heat production. We find that in all cases, the shell removes more heat from the interior than can be produced in the core by radioactive decay, resulting in cooling of the interior and the freezing of any ocean. Under likely conditions, a 40-km thick ocean made of pure water would freeze solid on a ∼30 Ma timescale. An ocean containing other chemical components will have a lower freezing point, but even a water-ammonia eutectic composition will only prolong the freezing, not prevent it. If the eccentricity of Enceladus were higher (e?0.015) in the past, the increased dissipation in the ice shell may have been sufficient to maintain a liquid layer. We cannot therefore rule out the presence of a transient ocean, as a relic of an earlier era of greater heating. If the eccentricity is periodically pumped up, the ocean may have thickened and thinned on a similar timescale as the orbital evolution, provided the ocean never froze completely. We conclude that the current heat flux of Enceladus and any possible subsurface ocean is not in steady-state, and is the remnant of an epoch of higher eccentricity and tidal dissipation.     

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.     

Tidal synchronization of close-in satellites and exoplanets. III. Tidal dissipation revisited and application to Enceladus

This paper deals with a new formulation of the creep tide theory (Ferraz-Mello in Celest Mech Dyn Astron 116:109, 2013—Paper I) and with the tidal dissipation predicted by the theory in the case of stiff bodies whose rotation is not synchronous but is oscillating around the synchronous state with a period equal to the orbital period. We show that the tidally forced libration influences the amount of energy dissipated in the body and the average perturbation of the orbital elements. This influence depends on the libration amplitude and is generally neglected in the study of planetary satellites. However, they may be responsible for a 27% increase in the dissipation of Enceladus. The relaxation factor necessary to explain the observed dissipation of Enceladus (\(\gamma =1.2{-}3.8\times 10^{-7}\ \mathrm{s}^{-1}\)) has the expected order of magnitude for planetary satellites and corresponds to the viscosity \(0.6{-}1.9 \times 10^{14}\) Pa s, which is in reasonable agreement with the value recently estimated by Efroimsky (Icarus 300:223, 2018) (\(0.24 \times 10^{14}\) Pa s) and with the value adopted by Roberts and Nimmo (Icarus 194:675, 2008) for the viscosity of the ice shell (\(10^{13}{-}10^{14}\) Pa s). For comparison purposes, the results are extended also to the case of Mimas and are consistent with the negligible dissipation and the absence of observed tectonic activity. The corrections of some mistakes and typos of paper II (Ferraz-Mello in Celest Mech Dyn Astron 122:359, 2015) are included at the end of the paper.     

Distribution of icy particles across Enceladus' surface as derived from Cassini-VIMS measurements

The surface of Enceladus consists almost completely of water ice. As the band depths of water ice absorptions are sensitive to the size of particles, absorptions can be used to map variations of icy particles across the surface. The Visual and Infrared Mapping Spectrometer (VIMS) observed Enceladus with a high spatial resolution during three Cassini flybys in 2005 (orbits EN 003, EN 004 and EN 011). Based on these data we measured the band depths of water ice absorptions at 1.04, 1.25, 1.5, and 2 μm. These band depths were compared to water ice models that represent theoretically calculated reflectance spectra for a range of particle diameters between 2 μm and 1 mm. The agreement between the experimental (VIMS) and model values supports the assumption that pure water ice characterizes the surface of Enceladus and therefore that variations in band depth correspond to variations in water ice particle diameters. Our measurements show that the particle diameter of water ice increases toward younger tectonically altered surface units with the largest particles exposed in relatively “fresh” surface material. The smallest particles were generally found in old densely cratered terrains. The largest particles (∼0.2 mm) are concentrated in the so called “tiger stripes” at the south pole. In general, the particle diameters are strongly correlated with geologic features and surface ages, indicating a stratigraphic evolution of the surface that is caused by cryovolcanic resurfacing and impact gardening.     

Enceladus' south polar sea

Recent observations of the south pole of Saturn's moon Enceladus by the Cassini spacecraft have revealed an active world, powered by internal heat. In this paper, we propose that localized subsurface melting on Enceladus has produced an internal south polar sea. Evidence for this localized sea comes from the shape of Enceladus, which does not match a differentiated body at its current orbital position. We show that melting induced by the observed heat flow at the south pole produces a large enough pit to match the shape of Enceladus with a differentiated rock and ice interior. Numerical modeling of melting and ice flow shows that the sea produced beneath the south pole is stable against inflow of ductile ice from its surroundings for the duration of the heating. The shape modification due to melting also produces a negative degree-two gravity anomaly, which can reorient the spin axis of Enceladus in order to place the sea at the pole.     

Tidally-induced melting events as the origin of south-pole activity on Enceladus

The intense activity at the south pole of Enceladus hints at an internal water reservoir. However, there is no direct evidence of liquid water at present and its long-term stability in the interior remains problematic. By modeling heat production and transfer in the ice shell in a spherical geometry, we show that tidal heating naturally leads to a concentration of convective hot upwellings in the south polar region, favoring the preservation of liquid water at depth. We show that large volumes of water are produced within the ice shell at the south pole during periods of elevated orbital eccentricity (3–5 times the present-day value). Strong lateral variations in the melt production and crystallization rates result in stress concentration in the south polar region, thus providing an explanation for the tectonic activity observed today. We predict that an internal ocean may be sustained over the long term as the consequence of repeated periods with elevated orbital eccentricity, leading to episodic melting and resurfacing events.     

Tidally Induced Volcanism

The dissipation of tidal energy causes the ongoing silicate volcanism on Jupiter's satellite, Io, and cryovolcanism almost certainly has resurfaced parts of Saturn's satellite, Enceladus, at various epochs distributed over the latter's history. The maintenance of tidal dissipation in Io and the occurrence of the same on Enceladus depends crucially on the maintenance of the respective orbital eccentricities by the existence of mean motion resonances with nearby satellites. A formation of the resonances among the Galilean satellites by differential expansion of the satellite orbits from tides raised on Jupiter by the satellites means the onset of the volcanism on Io could be relatively recent. If, on the other hand, the resonances formed by differential migration from resonant interactions of the satellites with the disk of gas and particles from which they formed, Io would have been at least intermittently volcanically active throughout its history. Either means of assembling the Galilean satellite resonances lead to the same constraint on the dissipation function of Jupiter Q _J 10⁶, where the currently high heat flux from Io seems to favor episodic heating as Io's eccentricity periodically increases and decreases. Either of the two models might account for sufficient tidal dissipation in the icy satellite Enceladus to cause at least occasional cryovolcanism over much of its history. However, both models are assumption-dependent and not secure, so uncertainty remains on how tidal dissipation resurfaced Enceladus.     

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.     

Release of N2, CH4, CO2, and H2O from surface ices on Enceladus

We vapor deposit at 20 K a mixture of gases with the specific Enceladus plume composition measured in situ by the Cassini INMS [Waite, J.H., Combi, M.R., Ip, W.H., Cravens, T.E., McNutt, R.L., Kasprzak, W., Yelle, R., Luhmann, J., Niemann, H., Gell, D., Magee, B., Fletcher, G., Lunine, J., Tseng, W.L., 2006. Science 311, 1419-1422] to form a mixed molecular ice. As the sample is slowly warmed, we monitor the escaping gas quantity and composition with a mass spectrometer. Pioneering studies [Schmitt, B., Klinger, J., 1987. Different trapping mechanisms of gases by water ice and their relevance for comet nuclei. In: Rolfe, E.J., Battrick, B. (Eds.), Diversity and Similarity of Comets. SP-278. ESA, Noordwijk, The Netherlands, pp. 613-619; Bar-Nun, A., Kleinfeld, I., Kochavi, E., 1988. Phys. Rev. B 38, 7749-7754; Bar-Nun, A., Kleinfeld, I., 1989. Icarus 80, 243-253] have shown that significant quantities of volatile gases can be trapped in a water ice matrix well above the temperature at which the pure volatile ice would sublime. For our Enceladus ice mixture, a composition of escaping gases similar to that detected by Cassini in the Enceladus plume can be generated by the sublimation of the H₂O:CO₂:CH₄:N₂ mixture at temperatures between 135 and 155 K, comparable to the high temperatures inferred from the CIRS measurements [Spencer, J.R., Pearl, J.C., Segura, M., Flasar, F.M., Mamoutkine, A., Romani, P., Buratti, B.J., Hendrix, A.R., Spilker, L.J., Lopes, R.M.C., 2006. Science 311, 1401-1405] of the Enceladus “tiger stripes.” This suggests that the gas escape phenomena that we measure in our experiments are an important process contributing to the gases emitted from Enceladus. A similar experiment for ice deposited at 70 K shows that both the processes of volatile trapping and release are temperature dependent over the temperature range relevant to Enceladus.     

Saturn's icy satellites: Thermal and structural models

Thermal histories of the small icy Saturnian satellites Mimas, Tethys, Dione, Rhea, and Iapetus are constructed by assuming that they formed as homogeneous ice-silicate mixtures. The models include effects of radiogenic and accretional heating, conductive and subsolidus convective heat transfer, and lithospheric growth. Accretional heating is unlikely to have melted the water ice in the interiors of these bodies and solid state creep of the predominately ice material precludes melting by radiogenic heating. Mimas is so small that its thermal evolution is essentially purely conductive; at present it is a cold, nearly isothermal body. Any subsolidus convection or thermal activity in Mimas would have been confined to a brief period in its early history and would have been due to a warm formation. The four largest satellites are big enough and contain sufficient heat-producing silicates that solid state convection beneath a rigid lithosphere is inevitable independent of initial conditions. Dione and Rhea have convective interiors for most of their thermal histories, while Tethys and Iapetus have mainly conductive thermal histories with early periods of convective 0activity. The thermal histories of the five satellites for the last 4 by are independent of initial conditions; at present they have cold, conductive interiors. The model thermal histories are qualitatively consistent with the appearances of these satellites: Mimas has an ancient heavily cratered surface, Tethys and probably Iapetus have both heavily cratered and more lightly cratered areas, and Dione and Rhea have extensively modified surfaces. Because of their similar sizes and densities, Mimas and Enceladus are expected to have similar surfaces and thermal histories, but instead Enceladus has the most modified surface of all the small icy Saturnian satellites. Our results suggest a heat source for Enceladus, in addition to radiogenic and accretional heating; tidal dissipation is a possibility. Because the water ice in these bodies does not melt, resurfacing must be accomplished by the melting of a low-melting-temperature minor component such as ammonia hydrate.     

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.     

MAIC-2, a latitudinal model for the Martian surface temperature, atmospheric water transport and surface glaciation

The Mars Atmosphere-Ice Coupler MAIC-2 is a simple, latitudinal model, which consists of a set of parameterisations for the surface temperature, the atmospheric water transport and the surface mass balance (condensation minus evaporation) of water ice. It is driven directly by the orbital parameters obliquity, eccentricity and solar longitude (L_s) of perihelion. Surface temperature is described by the Local Insolation Temperature (LIT) scheme, which uses a daily and latitude-dependent radiation balance. The evaporation rate of water is calculated by an expression for free convection, driven by density differences between water vapor and ambient air, the condensation rate follows from the assumption that any water vapour which exceeds the local saturation pressure condenses instantly, and atmospheric transport of water vapour is approximated by instantaneous mixing. Glacial flow of ice deposits is neglected. Simulations with constant orbital parameters show that low obliquities favour deposition of ice in high latitudes and vice versa. A transient scenario driven by a computed history of orbital parameters over the last 10 million years produces essentially monotonically growing polar ice deposits during the most recent 4 million years, and a very good agreement with the observed present-day polar layered deposits. The thick polar deposits sometimes continue in thin ice deposits which extend far into the mid latitudes, which confirms the idea of “ice ages” at high obliquity.     

Tidal dissipation,orbital evolution,and the nature of Saturn's inner satellites

Estimates of tidal damping times of the orbital eccentricities of Saturn's inner satellites place constraints on some satellite rigidities and dissipation functions Q. These constraints favor rock-like rather than ice-like properties for Mimas and probably Dione. Photometric and other observational data are consistent with relatively higher densities for these two satellites, but require lower densities for Tethys, Enceladus, and Rhea. This leads to a nonmonotonic density distribution for Saturn's inner satellites, apparently determined by different mass fractions of rocky materials. In spite of the consequences of tidal dissipation for the orbital eccentricity decay and implications for satellite compositions, tidal heating is not an important contributor to the thermal history of any Saturnian satellite.     

Parallel Proton Heating in Solar Wind Using Generalized (r, q) Distribution Function

Kinetic theory is used to calculate the power dissipated by obliquely propagating Alfvén waves to heat the solar wind protons, using the Generalized (r, q) distribution function. The evolution of power dissipation of protons with increasing heliocentric distance is subsequently determined. Comparison between theoretical and observational results with data shows good agreement, especially for the slow solar wind streams. Previous results where a Maxwellian distribution function was used to calculate the power dissipated did not match well with observations.