Tidal dynamical considerations constrain the state of an ocean on Enceladus

In previous work, solutions to the non-dissipative Laplace Tidal Equations (LTE) were used to provide bounds on the heat generated by the response of a subsurface ocean on Enceladus to an obliquity component of tidal forces. Here we improve these bounds using solutions from the LTE with a generic dissipation term explicitly added. We find solutions for a wider range of ocean tidal responses that include both unstratified (barotropic) and stratified (baroclinic) flow responses to obliquity as well as eccentricity components of the tidal forces. We consolidate the results in three ocean tidal scenarios on Enceladus that can explain the high heat fluxes (∼7 mW/m2 globally averaged) inferred from measurements by the Cassini spacecraft: (1) a deep (1–50 km) barotropic ocean responding to obliquity tidal forces, where obliquity is at least 0.1°; (2) a shallow (∼360 m) barotropic ocean responding to eccentricity tidal forces; (3) a stratified (baroclinic) ocean responding to eccentricity tidal forces where the density-weighted “equivalent depth” (typically much smaller than the ocean’s physical depth) is near 360 m. The ocean is assumed to be global, but extensions for a semi-global case are also described. A more general result which is independent of the specific scenarios proposed is that an ocean attempting to freeze (with an associated decrease in its liquid depth, which affects the ocean’s dynamical response to the tidal forcing) must first pass through resonant configurations with a greatly increased generation of ocean tidal heat (exceeding 1 W/m2 to 1 kW/m2) that would act to halt further freezing and stagnate the ocean state in this configuration so long as there is still orbital energy to provide the tidal forces. With an additional assumption that the ocean has evolved from a more energetic state where the depth of the liquid ocean was greater, we obtain the three scenarios proposed.