Planetary Climate 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:

Groundbased and spacecraft telescopic observations, combined with an intensive modeling effort, have greatly enhanced our understanding of hot giant planets and brown dwarfs over the past ten years. Although these objects are all fluid, hydrogen worlds with stratified atmospheres overlying convective interiors, they exhibit an impressive diversity of atmospheric behavior. Hot Jupiters are strongly irradiated, and a wealth of observations constrain the day-night temperature differences, circulation, and cloudiness. The intense stellar irradiation, presumed tidal locking and modest rotation leads to a novel regime of strong day-night radiative forcing. Circulation models predict large day-night temperature differences, global-scale eddies, patchy clouds, and, in most cases, a fast eastward jet at the equator—equatorial superrotation. The warm Jupiters lie farther from their stars and are not generally tidally locked, so they may exhibit a wide range of rotation rates, obliquities, and orbital eccentricities, which, along with the weaker irradiation, leads to circulation patterns and observable signatures predicted to differ substantially from hot Jupiters. Brown dwarfs are typically isolated, rapidly rotating worlds; they radiate enormous energy fluxes into space and convect vigorously in their interiors. Their atmospheres exhibit patchiness in clouds and temperature on regional to global scales—the result of modulation by large-scale atmospheric circulation. Despite the lack of irradiation, such circulations can be driven by interaction of the interior convection with the overlying atmosphere, as well as self-organization of patchiness due to cloud-dynamical-radiative feedbacks. Finally, irradiated brown dwarfs help to bridge the gap between these classes of objects, experiencing intense external irradiation as well as vigorous interior convection. Collectively, these diverse objects span over six orders of magnitude in intrinsic heat flux and incident stellar flux, and two orders of magnitude in rotation rate—thereby placing strong constraints on how the circulation of giant planets (broadly defined) depend on these parameters. A hierarchy of modeling approaches have yielded major new insights into the dynamics governing these phenomena.

Abstract:

M dwarf stars are excellent candidates around which to search for exoplanets, including temperate, Earth-sized planets. To evaluate the photochemistry of the planetary atmosphere, it is essential to characterize the UV spectral energy distribution of the planet's host star. This wavelength regime is important because molecules in the planetary atmosphere such as oxygen and ozone have highly wavelength-dependent absorption cross sections that peak in the UV (900–3200 Å). We seek to provide a broadly applicable method of estimating the UV emission of an M dwarf, without direct UV data, by identifying a relationship between noncontemporaneous optical and UV observations. Our work uses the largest sample of M dwarf star far- and near-UV observations yet assembled. We evaluate three commonly observed optical chromospheric activity indices—Hα equivalent widths and log10 LHα/Lbol, and the Mount Wilson Ca II H&K S and R'HK indices—using optical spectra from the HARPS, UVES, and HIRES archives and new HIRES spectra. Archival and new Hubble Space Telescope COS and STIS spectra are used to measure line fluxes for the brightest chromospheric and transition region emission lines between 1200 and 2800 Å. Our results show a correlation between UV emission-line luminosity normalized to the stellar bolometric luminosity and Ca II R'HK with standard deviations of 0.31–0.61 dex (factors of ~2–4) about the best-fit lines. We also find correlations between normalized UV line luminosity and Hα log10 LHα/Lbol and the S index. These relationships allow one to estimate the average UV emission from M0 to M9 dwarfs when UV data are not available.