Dynamical regimes of giant planet polar vortices

We present a numerical model that reveals a mechanism governing the polar atmospheric dynamics of Jupiter, Saturn, Uranus and Neptune. Exploration of the polar regions of the gas giants has produced surprisingly diverse results, with Cassini finding a single, intense, compact polar cyclone precisely centered on each pole of Saturn, and Voyager data and ground-based observations suggesting Uranus and Neptune have dominant, single polar cyclones as well. The Juno spacecraft at Jupiter finds several tightly packed cyclones surrounding a central cyclone offset from the poles. These discoveries raise questions about the mechanism that differentiates these polar atmospheric dynamics regimes. To help determine what physical mechanisms control these differences, we use the Explicit Planetary Isentropic Coordinate (EPIC) model to carry out forced-turbulence shallow-water simulations in a gamma-plane configuration, i.e. a Cartesian grid with a pole placed at the center. The model is forced by small-scale stochastic mass pulses that parametrically represent cumulus storms. The effects of three parameters, the planetary Burger number, Bu =(L-d/a)(2) (L-d is the Rossby deformation radius, a is the planetary radius), input storm strength, s, and proportion of cyclonic and anticyclonic storms injected into the domain, alpha, are systematically investigated. Bu emerges to be the most important, able to distinguish between four distinct dynamical regimes, matching those of the giant planets, which from large to small Bu, are: i) a large cyclonic polar vortex (i.e., Uranus/Neptune-like), ii) a compact intense cyclonic polar vortex (Saturn-like), iii) two large vortices or one vortex offset from the pole (transitional), and iv) meandering jets with no centrally dominant vortex, or with multiple circumpolar cyclones (Jupiter-like). The boundaries of these regimes are found to be only slightly modulated by the values of s and alpha. By applying this correlation with respect to Bu in reverse, an observation of a particular polar regime could in principle be used to constrain L-d.