Geophysical and Astrophysical Fluid Dynamics

Abstract:

This archive contains the Python code used to analyse and plot the data in Hammond & Lewis 2021, "The rotational and divergent components of atmospheric circulation on tidally locked planets", as well as the data from the "terrestrial" simulation of the atmosphere of a rocky planet using the general circulation model ExoFMS. It contains three files: 1) HL21_plotter.ipynb This is a Jupyter notebook containing Python code. It reads the data from the ExoFMS simulation and finds its rotational and divergent parts. It then plots the figures used in Hammond & Lewis 2021. 2) data/rotdiv-terr-control-1000-2000_atmos_average_interp.nc The "terrestrial" simulation output, interpolated to uniform pressure levels. This is used to plot quantities such as velocity at a constant pressure. 3) data/rotdiv-terr-control-1000-2000_atmos_average.nc The "terrestrial" simulation output, on the raw model sigma-pressure levels. This is used to calculate the dry static energy budget. The paper also uses a "Hot Jupiter" simulation from the THOR GCM. This is from "THOR 2.0: Major Improvements to the Open-Source General Circulation Model" (Deitrick et al. 2020). The data is available on request to Russell Deitrick (russell.deitrick@csh.unibe.ch). The same analysis can be made using HL21_plotter.ipynb, with small modifications due to the different grid in THOR.

Abstract:

Images of the giant planets Jupiter and Saturn show highly turbulent storms and swirling clouds that reflect the intensity of turbulence in their atmospheres. Quantifying planetary turbulence is inaccessible to conventional tools, however, since they require large quantities of spatially and temporally resolved data. Here we show, using experiments, observations, and simulations, that potential vorticity (PV) is a straightforward and universal diagnostic that can be used to estimate turbulent energy transfer in a stably stratified atmosphere. We use the conservation of PV to define a length scale, LM, representing a typical distance over which PV is mixed by planetary turbulence. LM increases as the turbulent intensity increases and can be estimated from any latitudinal PV profile. Using this principle, we estimate LM within Jupiter's and Saturn's tropospheres, showing for the first time that turbulent energy transfer in Saturn's atmosphere is four times less intense than Jupiter's.

Abstract:

The wintertime Southern Hemisphere extratropical circulation exhibits considerable zonal asymmetries. We investigate the roles of various surface boundary conditions in shaping the mean state using a semi-realistic, atmosphere-only climate model. We find, in agreement with previous literature, that tropical sea surface temperature (SST) patterns are an important contributor to the mean state, while midlatitude SSTs and sea ice extent play a smaller role. Our main finding is that Antarctic orography has a first-order effect on the structure of the midlatitude circulation. In the absence of Antarctic orography, equatorward eddy momentum fluxes associated with the orography are removed and hence convergence of eddy momentum in midlatitudes is reduced. This weakens the Indian Ocean jet, making Rossby wave propagation downstream to the South Pacific less favorable. Consequently, the flow stagnates over the mid- to high-latitude South Pacific and the characteristic split jet pattern is destroyed. Removing Antarctic orography also results in a substantial warming over East Antarctica partly because transient eddies are able to penetrate farther poleward, enhancing poleward heat transport. However, experiments in which a high-latitude cooling is applied indicate that these temperature changes are not the primary driver of circulation changes in the midlatitudes. Instead, we invoke a simple barotropic mechanism in which the orographic slope creates an effective potential vorticity gradient that alters the eddy momentum flux.

Abstract:

Continuous structural parameterization (CSP) is a proposed method for approximating different numerical model parameterizations of the same process as functions of the same grid‐scale variables. This allows systematic comparison of parameterizations with each other and observations or resolved simulations of the same process. Using the example of two convection schemes running in the Met Office Unified Model (UM), we show that a CSP is able to capture concisely the broad behavior of the two schemes, and differences between the parameterizations and resolved convection simulated by a high resolution simulation. When the original convection schemes are replaced with their CSP emulators within the UM, basic features of the original model climate and some features of climate change are reproduced, demonstrating that CSP can capture much of the important behavior of the schemes. Our results open the possibility that future work will estimate uncertainty in model projections of climate change from estimates of uncertainty in simulation of the relevant physical processes.