Geophysical and Astrophysical Fluid Dynamics

Abstract:

We investigate the dynamics of Jupiter's upper troposphere and lower stratosphere using a General Circulation Model that includes two-stream radiation and optional heating from below. Based on the MITgcm dynamical core, this is a new generation of the Oxford Jupiter model [Zuchowski, L.C. et al., 2009. Plan. Space Sci., 57, 1525--1537, doi:10.1016/j.pss.2009.05.008]. We simulate Jupiter's atmosphere at up to 0.7 degree horizontal resolution with 33 vertical levels down to a pressure of 18 bar, in configurations with and without a 5.7 W/m2 interior heat flux. Simulations ran for 130000-150000 days to allow the deep atmosphere to come into radiative equilibrium. Baroclinic instability generates alternating, eddy-driven, midlatitude jets in both cases. With interior heating the zonal jets migrate towards the equator and become barotropically unstable. This generates Rossby waves that radiate away from the equator, depositing westerly momentum there via eddy angular momentum flux convergence and spinning up a super-rotating 20 m/s equatorial jet throughout the troposphere. There are 30-35 zonal jets with latitudinal separation comparable with the real planet, and there is strong eddy activity throughout. Without interior heating the jets do not migrate and a divergent eddy angular momentum flux at the equator spins up a broad, 50 m/s sub-rotating equatorial jet with weak eddy activity at low latitudes.

Abstract:

We examine the ammonia and water cycles in Jupiter's upper troposphere and lower stratosphere during spin-up of a multiple zonal jet circulation using the Oxford Jupiter GCM. Jupiter's atmosphere is simulated at 512 x 256 horizontal resolution with 33 vertical levels between 0.01 and 18 bar, putting the lowest level well below the expected water cloud base. Simulations with and without a 5.7 W/m2 interior heat source were run for 130000-150000d to allow the deep atmosphere to come into radiative-convective-dynamical equilibrium, with variants on the interior heating case including varying the initial tracer distribution, particle condensate diameter, and cloud process timescales. The cloud scheme includes simple representations of the ammonia and water cycles. Ammonia vapour changes phase to ice, and reacts with hydrogen sulphide to produce ammonium hydrosulphide. Water changes phases between vapour, liquid, and ice depending on local environmental conditions, and all condensates sediment at their respective Stokes velocities. With interior heating, clouds of ammonia ice, ammonium hydrosulphide ice, and water ice form with cloud bases around 0.4 bar, 1.5 bar, and 3 bar respectively. Without interior heating the ammonia cloud base forms in the same way, but the ammonium hydrosulphide and water clouds sediment to the bottom of the domain. The liquid water cloud is either absent or extremely sparse. Zonal structures form that correlate regions of strong latitudinal shear with regions of constant condensate concentration, implying that jets act as barriers to the mixing. Regions with locally high and low cloud concentrations also correlated with regions of upwelling and downwelling, respectively. Shortly after initialisation, the ammonia vapour distribution up to the cloud base resembles the enhanced concentration seen in Juno observations, due to strong meridional mean circulation at the equator. The resemblance decays rapidly over time, but suggests that at least some of the relevant physics is captured by the model. The comparison should improve with additional microphysics and better representation of the deep ammonia reservoir.

Abstract:

The energetics of possible global atmospheric circulation patterns in an Earth-like atmosphere are explored using a simplified GCM based on the University of Hamburg’s Portable University Model for the Atmosphere (designated here as PUMA-S), forced by linear relaxation towards a prescribed temperature field and subject to Rayleigh surface drag and hyperdiffusive dissipation. Results from a series of simulations, obtained by varying planetary rotation rate Ω with an imposed equator-to-pole temperature difference, were analysed to determine the structure and magnitude of the heat transport and other contributions to the energy budget for the time-averaged, equilibrated flow. These show clear trends with rotation rate, with the most intense Lorenz energy cycle for an Earth-sized planet occurring with a rotation rate around half that of the present day Earth (i.e. Ω* = Ω/ΩE = 1/2, where ΩE is the rotation rate of the Earth). KE and APE spectra, EK(n) and EA(n) (where n is total spherical wavenumber), also show clear trends with rotation rate, with n^-3 enstrophy-dominated spectra around Ω* = 1 and steeper (~ n^-5) slopes in the zonal mean flow with little evidence for the n^-5/3 spectrum anticipated for an inverse KE cascade. Instead, both KE and APE spectra become almost flat at scales larger than the internal Rossby radius, Ld, and exhibit near-equipartition at high wavenumbers. At Ω* << 1, the spectrum becomes dominated by KE with EK(n) ~ (2 - 3)EA(n) at most wavenumbers and a slope that tends towards n^-5/3 across most of the spectrum. Spectral flux calculations show that enstrophy and APE are almost always cascaded downscale, regardless of rotation rate. KE cascades are more complicated, however, with downscale transfers across almost all wavenumbers, dominated by horizontally divergent modes, for Ω* ≲ 1/4. At higher rotation rates, transfers of KE become increasingly dominated by rotational (horizontally non-divergent) components with strong upscale transfers (dominated by eddy-zonal flow interactions) for scales larger than Ld and weaker downscale transfers for scales smaller than Ld.