We present a new model of Enceladus' internal structure based on the recent shape model by Tajeddine et al. (2017) and the gravity model by Iess et al. (2014). Assuming that the rocky core is homogeneous and in hydrostatic equilibrium, we derive a set of structural models that accurately reproduce the main characteristics of Enceladus' gravity field.
In order to restrict the range of acceptable models we analyze the degree of compensation (ratio of the bottom to the surface load) at spherical harmonic degrees well constrained by the gravity data. In agreement with previous studies, we find that Enceladus' ice shell is close to equilibrium, with the degree of compensation approaching one for models with a hydrostatic core having a radius between 190 km and 195 km.
By computing the flow of ice driven by variations in hydrostatic pressure on the ice water/interface, we demonstrate that the ice shell is in steady state, as suggested by the gravity and shape data, only if the viscosity of ice at the melting temperature is equal to or higher than 3 x 10(14) Pa s, corresponding to diffusion creep with a grain size of 1 mm or larger. The topographic anomalies are maintained by phase changes at the ice/water interface with a melting/freezing rate of a few mm/yr or smaller.
This process is controlled by the heat flux at the top of the ocean, characterized by a strong degree-2 zonal component with amplitudes exceeding 60 rnW/m(2) in both south and north polar regions.