The zonal winds on the surfaces of giant planets vary with latitude. Jupiter and Saturn, for example, have several bands of alternating eastward (prograde) and westward (retrograde) jets relative to the angular velocity of their global magnetic fields. These surface wind profiles are likely manifestations of the variations in depth and latitude of angular velocity deep within the liquid interiors of these planets. Two decades ago it was proposed that this differential rotation could be maintained by vortex stretching of convective fluid columns that span the interiors of these planets from the northern hemisphere surface to the southern hemisphere surface. This now classic mechanism explains the differential rotation seen in laboratory experiments and in computer simulations of, at best, weakly turbulent convection in rotating constant-density fluid spheres. However, these experiments and simulations are poor approximations for the density-stratified strongly-turbulent interiors of giant planets. The long thin global convective columns predicted by the classic geostrophic theory for these planets would likely not develop. Here we propose a much more robust mechanism for maintaining differential rotation in radius based on the local generation of vorticity as rising plumes expand and sinking plumes contract. Our high-resolution two-dimensional computer simulations demonstrate how this mechanism could maintain either prograde or retrograde surface winds in the equatorial region of a giant planet depending on how the density scale height varies with depth.

Local generation of vorticity occurs in rotating density-stratified fluids as fluid parcels move radially, expanding or contracting with respect to the background density stratification. Thermal convection in rotating 2D equatorial simulations demonstrates this mechanism. The convergence of the vorticity into zonal flow structures as a function of radius depends on the shape of the density profile, with the prograde jet forming in the region of the disk where the greatest number of density scale heights occurs. The number of stable jets that form in the fluid increases with decreasing Ekman number and decreases with increasing thermal driving. This local form of vorticity generation via the density stratification is likely to be of great importance in bodies that are quickly rotating, highly turbulent, and have large density changes, such as Jovian planets. However, it is likely to be of lesser importance in the interiors of planets such as the Earth, which have smaller density stratifications and are less turbulent.

We study how the pattern of thermal convection and differential rotation in the interior of a giant gaseous planet is affected by the presence of a small solid core as a function of the planetary rotation rate. We show, using 2D anelastic, hydrodynamic simulations, that the presence of a small solid core results in significantly different flow structure relative to that of a fully convective interior only if there is little or no planetary rotation.