We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the ...large-scale atmospheric flow in hot Jupiter atmospheres. We coupled a flexible chemical kinetics scheme to the Met Office Unified Model, which enables the study of the interaction of chemistry, radiative transfer, and fluid dynamics. We used a newly-released “reduced” chemical network, comprising 30 chemical species, that was specifically developed for its application in 3D atmosphere models. We simulated the atmospheres of the well-studied hot Jupiters HD 209458b and HD 189733b which both have dayside–nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH
4
), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. Three-dimensional mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D – and even 2D – chemical mixing. We find that CH
4
, carbon dioxide (CO
2
), and ammonia (NH
3
) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, which are detectable with near-future instruments.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat ...overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as global circulation models (GCMs) it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P–T profiles in our models.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
ABSTRACT
We present four new secondary eclipse observations for the ultrahot Jupiter WASP-121b acquired using the Hubble Space Telescope Wide Field Camera 3. The eclipse depth is measured to a median ...precision of 60 ppm across 28 spectroscopic channels spanning the 1.12–$1.64\, \mu {\rm m}$ wavelength range. This is a considerable improvement to the 90 ppm precision we achieved previously for a single eclipse observation using the same observing set-up. Combining these data with those reported at other wavelengths, a blackbody spectrum for WASP-121b is ruled out at >6σ confidence and we confirm the interpretation of previous retrieval analyses that found the data are best explained by a dayside thermal inversion. The updated spectrum clearly resolves the water emission band at 1.3–$1.6\, \mu {\rm m}$, with higher signal-to-noise than before. It also fails to reproduce a bump in the spectrum at $1.25\, \mu {\rm m}$ derived from the first eclipse observation, which had tentatively been attributed to VO emission. We conclude that the latter was either a statistical fluctuation or a systematic artefact specific to the first eclipse data set.
To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office Global Circulation Model (GCM), the Unified Model (UM), to these ...exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable agreement between observations and both our days side emission and hot spot offset, however, our night side emissions is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009, ApJ, 699, 564) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the “vertical coherency” found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
In this paper we present 3D atmospheric simulations of the hot Jupiter HD 189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved ...wind. Our model consistently couples the treatment of dynamics, radiative transfer, and chemistry, completing the feedback cycle between these three important processes. The effect of wind-driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD 209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD 189733b, with changes in both the temperature and wind velocities reaching ∼10%. We also present the contribution function, diagnosed from our simulations, and show that wind-driven chemistry has a significant impact on its 3D structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind-driven chemistry on the thermal emission, particularly within the 3.6 m Spitzer/IRAC channel.
The treatment of radiation transport in global circulation models (GCMs) is crucial for correctly describing 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 the 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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Identification of habitable planets beyond our solar system is a key goal of current and future space missions. Yet habitability depends not only on the stellar irradiance, but equally on constituent ...parts of the planetary atmosphere. Here we show, for the first time, that radiatively active mineral dust will have a significant impact on the habitability of Earth-like exoplanets. On tidally-locked planets, dust cools the day-side and warms the night-side, significantly widening the habitable zone. Independent of orbital configuration, we suggest that airborne dust can postpone planetary water loss at the inner edge of the habitable zone, through a feedback involving decreasing ocean coverage and increased dust loading. The inclusion of dust significantly obscures key biomarker gases (e.g. ozone, methane) in simulated transmission spectra, implying an important influence on the interpretation of observations. We demonstrate that future observational and theoretical studies of terrestrial exoplanets must consider the effect of dust.
We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an “Earth-like” ...atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the planetary conditions using a range of eccentricities guided by the observational constraints. Overall, our results are in agreement with previous studies in suggesting Proxima Centauri B may well have surface temperatures conducive to the presence of liquid water. Moreover, we have expanded the parameter regime over which the planet may support liquid water to higher values of eccentricity (≳0.1) and lower incident fluxes (881.7 W m-2) than previous work. This increased parameter space arises because of the low sensitivity of the planet to changes in stellar flux, a consequence of the stellar spectrum and orbital configuration. However, we also find interesting differences from previous simulations, such as cooler mean surface temperatures for the tidally-locked case. Finally, we have produced high-resolution planetary emission and reflectance spectra, and highlight signatures of gases vital to the evolution of complex life on Earth (oxygen, ozone and carbon dioxide).
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK