Climate Simulations of Hot Jupiters: Developing and Applying an Accurate Radiation Scheme

To date more than 1500 exoplanets have been discovered. A large number of these are hot Jupiters, Jupiter-sized planets orbiting < 0.1 au from their parent stars, due to limitations in observational techniques making them easier to detect than smaller planets in wider orbits. This is also, for the same reasons, the class of exoplanets with the most observational constraints. Due to the very large interaction between these planets and their parent stars they are believed to be tidally locked, causing a large temperature contrast between the permanently hot day side and colder night side. There are still many open questions about these planets. Many are observed to have inflated radii, i.e. the observed radius is larger for a given mass than evolutionary models predict. A mechanism that can transport some of the stellar heating into the interior of the planet may be able to explain this. The presence of hazes or clouds has been inferred on some planets, but their composition and distribution remain unknown. According to chemical equilibrium models TiO and VO should be present on the day side of the hottest of these planets, but these molecules have not yet been detected. Cold traps, where these molecules condense out on the night side, have been suggested to explain this. The efficiency of the heat redistribution from the day side to the night side has been found to vary significantly between different planets; the mechanism behind this is still unknown. To begin to answer many of these questions we need models capturing the three-dimensional nature of the atmospheres of these planets. General circulation models (GCMs) do this by solving the equations of fluid dynamics for the atmosphere coupled to a radiative transfer scheme. GCMs have previously been applied to several exoplanets, but many solve simplified fluid equations (shallow water or primitive equations) or highly parametrised radiation schemes (temperature-forcing, gray or band-averaged opacities). We here present an adaptation of the Met Office Unified Model (UM), a GCM used for weather predictions and climate studies for the Earth, to hot Jupiters. The UM solves the full 3D Euler equations for the fluid, and the radiation scheme uses the two-stream approximation and correlated-k method, which are state of the art for both Earth and exoplanet GCMs. This makes it ideally suited for the study of hot Jupiters. An important part of this work is devoted to the adaptation of the radiation scheme of the UM to hot Jupiters. This includes calculation of opacities for the main absorbers in these atmospheres from state-of-the-art high temperature line lists, the calculation of k-coefficients from these opacities, and making sure all aspects of the scheme perform satisfactorily at high temperatures and pressures. We have tested approximations made in previous works such as the two-stream approximation, use of band-averaged opacities and different treatments of gaseous overlap. Uncertainties in current models, such as the lack of high temperature line broadening parameters for these atmospheres, are discussed. We couple the adapted radiation scheme to the UM dynamical core, which has been tested independently. Our first application is devoted to one of the most well-observed hot Jupiters, HD 209458b. Differences between previous modelling works and our model are discussed, and we compare results from the full coupled model with results obtained using a temperature-forcing scheme. We have also developed a tool to calculate synthetic phase curves, and emission and transmission spectra from the output of our 3D model. This enables us to directly compare our model results to observations and test the effect of various parameters and model choices on observable quantities

A uniform analysis of HD209458b Spitzer/IRAC lightcurves with Gaussian process models

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.Comment: 18 pages, 5 figures, 6 tables. Accepted by MNRA

Results from a set of three-dimensional numerical experiments of a hot Jupiter atmosphere

We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.Comment: 26 pages, 22 Figures. Accepted for publication in Astronomy and Astrophysic

Habitable Climate Scenarios for Proxima Centauri b With a Dynamic Ocean

The nearby exoplanet Proxima Centauri b will be a prime future target for characterization, despite questions about its retention of water. Climate models with static oceans suggest that an Earth-like Proxima b could harbor a small dayside region of surface liquid water at fairly warm temperatures despite its weak instellation. We present the first 3-dimensional climate simulations of Proxima b with a dynamic ocean. We find that an ocean-covered Proxima b could have a much broader area of surface liquid water but at much colder temperatures than previously suggested, due to ocean heat transport and depression of the freezing point by salinity. Elevated greenhouse gas concentrations do not necessarily produce more open ocean area because of possible dynamic regime transitions. For an evolutionary path leading to a highly saline present ocean, Proxima b could conceivably be an inhabited, mostly open ocean planet dominated by halophilic life. For an ocean planet in 3:2 spin-orbit resonance, a permanent tropical waterbelt exists for moderate eccentricity. Simulations of Proxima Centauri b may also be a model for the habitability of planets receiving similar instellation from slightly cooler or warmer stars, e.g., in the TRAPPIST-1, LHS 1140, GJ 273, and GJ 3293 systems.Comment: Submitted to Astrobiology; 38 pages, 12 figures, 5 table

A Hybrid Line List for CH4 and Hot Methane Continuum

Molecular line lists (a catalogue of transition frequencies and line strengths) are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. Methods. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. Results. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34 to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power R = 1,000,000 to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections

Accuracy tests of radiation schemes used in hot Jupiter global circulation models

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.European Community’s Seventh Framework Programme FP7/2007-2013Science & Technology Facilities Council (STFC)Royal Societ

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a 3-Dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of Solar System and exoplanetary terrestrial planets. Its parent model, known as ModelE2 (Schmidt et al. 2014), is used to simulate modern and 21st Century Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (slowly rotating to more rapidly rotating than modern Earth, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the Solar System such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn's moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents we can then expand its capabilities further to those exoplanetary rocky worlds that have been discovered in the past and those to be discovered in the future. We discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.Comment: Revisions since previous draft. Now submitted to Astrophysical Journal Supplement Serie