Abstract
We investigated the effect of variations in ice shell thickness and of the tiger stripe fractures crossing Enceladus' south polar terrain on the moon's tidal deformation by performing finite element calculations in three-dimensional geometry. The combination of thinning in the polar region and the presence of faults has a synergistic effect that leads to an increase of both the displacement and stress in the south polar terrain by an order of magnitude compared to that of the traditional model with a uniform shell thickness and without faults. Assuming a simplified conductive heat transfer and neglecting the heat sources below the ice shell, we computed the global heat budget of the ice shell. For the inelastic properties of the shell described by a Maxwell viscoelastic model, we show that unrealistically low average viscosity of the order of 10(13) Pa s is necessary for preserving the volume of the ocean, suggesting the important role of the heat sources in the deep interior. Similarly, low viscosity is required to predict the observed delay of the plume activity, which hints at other delaying mechanisms than just the viscoelasticity of the ice shell. The presence of faults results in large spatial and temporal heterogeneity of geysering activity compared to the traditional models without faults. Our model contributes to understanding the physical mechanisms that control the fault activity, and it provides potentially useful information for future missions that will sample the plume for evidence of life. Key Words: Enceladus-Tidal deformation-Faults-Variable ice shell thickness-Tidal heating-Plume activity and timing. Astrobiology 17, 941-954.