Viscoelastic relaxation of Enceladus's ice shell

We propose a dynamic model relating the non-hydrostatic long-wavelength topography on Enceladus to loading at the base of the ice shell. In dynamic equilibrium, the surface topography induced by a bottom load is very close to that inferred for Airy isostasy. However, due to the small size of the moon the relaxation to equilibrium is significantly longer than the maximum Maxwell relaxation time. During the relaxation, the upper load/bottom load ratio can be either smaller than 1 or larger than 1, depending on the initial state and the loading history.