The deep ocean (~100 km) of Europa, Jupiter's moon, is covered by a thick (tens of km) icy shell, and is one of the most probable places in the solar system to find extraterrestrial life. Yet, its ocean dynamics and its interaction with the ice cover have so far received little attention. Previous studies suggested that Europa's ocean is turbulent, yet neglected to take into account the effects of ocean salinity and appropriate boundary conditions for the ocean's temperature. Here, the ocean dynamics of Europa is studied using global ocean models that include non-hydrostatic effects, a full Coriolis force, consistent top and bottom heating boundary conditions, and including the effects of melting and freezing of ice on salinity. The density is found to be dominated by salinity effects and the ocean is very weakly stratified. The ocean exhibits strong transient vertical convection, eddies, low latitude zonal jets and Taylor columns parallel to Europa's axis of rotation. In the equatorial region, the Taylor columns do not intersect the ocean bottom and propagate equatorward, while off the equator, the Taylor columns are static. The meridional oceanic heat transport is intense enough to result in a nearly uniform ice thickness, that is expected to be observable in future missions.

Planets with high obliquity receive more radiation in the polar regions than at low latitudes, and thus, assuming an ocean-covered surface with sufficiently high heat capacity, their meridional temperature gradient was shown to be reversed for the entire year. The objective of this work is to investigate the drastically different general circulation of such planets, with an emphasis on the tropical Hadley circulation and the mid-latitude baroclinic eddy structure. We use a 3D dry dynamic core model, accompanied by an eddy-free configuration and a generalized 2D Eady model. When the meridional temperature gradient $T_y$ is reversed, the Hadley cell remains in the same direction, because the surface wind pattern and hence the associated meridional Ekman transport are not changed, as required by the baroclinic eddy momentum transport. The Hadley cell under reversed $T_y$ also becomes much shallower and weaker, even when the magnitude of the gradient is the same as in the normal case. The shallowness is due to the bottom-heavy structure of the baroclinic eddies in the reverse case, and the weakness is due to the weak wave activity. We propose a new mechanism to explain the mid-latitude eddy structure for both cases, and verify it using the generalized Eady model. With seasonal variations included, the annual mean circulation resembles that under perpetual annual mean setup. Approaching the solstices, a strong cross-equator Hadley cell forms in both cases, and about 2/3 of the Hadley circulation is driven by eddies, as shown by eddy-free simulations and using a decomposition of the Hadley cell.

The nature of the flow below the cloud level on Jupiter and Saturn is still unknown. Relating the flow on these planets to perturbations in their density field is key to the analysis of the gravity measurements expected from both the Juno (Jupiter) and Cassini (Saturn) spacecrafts during 2016-17. Both missions will provide latitude-dependent gravity fields, which in principle could be inverted to calculate the vertical structure of the observed cloud-level zonal flow on these planets. Theories to date connecting the gravity field and the flow structure have been limited to potential theories under a barotropic assumption, or estimates based on thermal wind balance that allow analyzing baroclinic wind structures, but have made simplifying assumptions. Those include the effects of the deviations from spherical symmetry, the centrifugal force due to density perturbations, and self-gravitational effects of the density perturbations. Recent studies attempted to include some effects but not in a self-consistent manner. The present study introduces such a self-consistent perturbation approach to the thermal wind balance that incorporates all physical effects, and applies it to several example wind structures, both barotropic and baroclinic. The contribution of each term is analyzed, and the results are compared in the barotropic limit to those of potential theory. It is found that the dominant balance involves the original simplified thermal wind approach. This balance produces a good order-of-magnitude estimate of the gravitational moments, and is able, therefore, to address the order one question of how deep the flows are given measurements of gravitational moments. The additional terms are significantly smaller and none of these terms is dominant, so any approximation attempting to improve over the simplified thermal wind approach needs to include all other terms.