Abstract
At the south pole of Saturn's icy moon Enceladus, eruptions of water vapour and ice emanate from warm tectonic ridges(1-4). Observations in the infrared(5) and visible(6) spectra have shown an orbital modulation of the plume brightness, which suggests that the eruption activity is influenced by tidal forces.
However, the observed activity seems to be delayed by several hours with respect to predictions based on simple tidal models(6,7). Here we simulate the viscoelastic tidal response of Enceladus with a full three-dimensional numerical model(8,9) and show that the delay in eruption activity may be a natural consequence of the viscosity structure in the south-polar region and the size of the putative subsurface ocean.
By systematically comparing simulations of variations in normal stress along faults with plume brightness data, we show that the observed activity is reproduced for two classes of interior models with contrasting thermal histories: a low-viscosity convective region above a polar sea extending about 45 degrees-60 degrees from the south pole at a depth below the surface as small as 30 km, or a convecting ice shell of 60-70 km in thickness above a global ocean. Our analysis further shows that the eruption activity is controlled by the average normal stress applied across the cracks, thus providing a constraint on the eruption mechanism.