Of profound astrobiological interest, Enceladus appears to have a global saline subsurface ocean, indicating water-rock reaction at present or in the past, an important mechanism in the moon’s potential habitability. Here, we investigate how salinity and the partition of heat production between the silicate core and the ice shell affect ocean dynamics and the associated heat transport—a key factor determining equilibrium ice shell geometry. Assuming steady-state conditions, we show that the meridional overturning circulation of the ocean, driven by heat and salt exchange with the poleward-thinning ice shell, has opposing signs at very low and very high salinities. Regardless of these differing circulations, heat and fresh water converge toward the equator, where the ice is thick, acting to homogenize thickness variations. Among scenarios explored here, the pronounced ice thickness variations observed on Enceladus are most consistent with heating that is predominantly in the ice shell and a salinity of intermediate range.

We numerically explore convection and general circulation of an ocean, encased in a spherical shell of uniform thickness, heated from below by a spatially uniform heat flux, and whose temperature at the upper surface is relaxed to the freezing point of water. The role of salt is not considered. We describe the phenomenology and equilibrium solutions across a broad range of two key non‐dimensional numbers: the natural Rossby number, a measure of the influence of rotation, and the ratio of inner to outer radius of the moon's ocean, a measure of the geometry of the moon's tangent cylinder. Two distinct regimes of circulation are identified, both dominated by Taylor columns aligned with the rotation axis—“plumes” and “rolls” which predominate inside and outside the tangent cylinder, respectively. Inside the tangent cylinder, convective plumes align with Taylor columns that extend from the bottom to the ice shell. The plumes energize geostrophic turbulence which in turn generates a general circulation consisting counter‐rotating zonal jets. Moreover, the plumes are efficient at transferring heat from the bottom to the surface, resulting in loss of heat from the ocean to the polar ice. If the plumes are suppressed rolls outside the tangent cylinder become the dominant mode of heat transfer resulting in equatorial cooling. We conclude that if moons such as Enceladus and Europa were to be predominantly heated from below, they will likely have a “Jovian‐like” circulation: an unstratified, turbulent, geostrophically controlled ocean with strong “Taylor column” behavior and a circulation dominated by counter‐rotating zonal jets.

Of profound astrobiological interest is that not only does Enceladus have a water ocean, but it also appears to be salty, important for its likely habitability. Here, we investigate how salinity affects ocean dynamics and equilibrium ice shell geometry and use knowledge of ice shell geometry and tidal heating rates to help constrain ocean salinity. We show that the vertical overturning circulation of the ocean, driven from above by melting and freezing and the temperature dependence of the freezing point of water on pressure, has opposing signs at very low and very high salinities. In both cases, heat and freshwater converges toward the equator, where the ice is thick, acting to homogenise thickness variations. In order to maintain observed ice thickness variations, ocean heat transport should not overwhelm tidal heating rates within the ice, which are small in equatorial regions. This can only happen when the ocean’s salinity has intermediate values, order 20 psu. In this case polar-sinking driven by meridional temperature variations is largely canceled by equatorial-sinking circulation driven by salinity variations and a consistent ocean circulation, ice shell geometry and tidal heating rate can be achieved.

Of profound astrobiological interest is that not only does Enceladus have a water ocean, but it also appears to be salty, important for its likely habitability. Here, we investigate how salinity affects ocean dynamics and equilibrium ice shell geometry and use knowledge of ice shell geometry and tidal heating rates to help constrain ocean salinity. We show that the vertical overturning circulation of the ocean, driven from above by melting and freezing and the temperature dependence of the freezing point of water on pressure, has opposing signs at very low and very high salinities. In both cases, heat and freshwater converges toward the equator, where the ice is thick, acting to homogenize thickness variations. In order to maintain observed ice thickness variations, ocean heat convergence should not overwhelm heat loss rates through the equatorial ice sheet. This can only happen when the ocean's salinity has intermediate values, order 20 psu. In this case polar-sinking driven by meridional temperature variations is largely canceled by equatorial-sinking circulation driven by salinity variations and a consistent ocean circulation, ice shell geometry and tidal heating rate can be achieved.Since the Cassini and Galileo mission, Enceladus (a satellite of Saturn) and Europa (a satellite of Jupiter) have been revealed to have a high astrobiological potential, satisfying all three 1