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.
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.
Abstract
We present results from one-dimensional atmospheric simulations investigating the effect of varying the carbon-to-oxygen (C/O) ratio on the thermal structure, chemical composition, and ...transmission and emission spectra, for irradiated hydrogen-dominated atmospheres. We find that each of these properties of the atmosphere is strongly dependent on the individual abundances (relative to hydrogen) of carbon and oxygen. We confirm previous findings that different chemical equilibrium compositions result from sets of element abundances with the same C/O ratio but with different individual abundances of carbon and oxygen. We investigate the effect of this difference in composition on the thermal structure and simulated spectra. We also simulate observations using the PandExo tool and show that these differences are observationally significant with current (i.e. Hubble Space Telescope) and future (i.e. James Webb Space Telescope) instruments. We conclude that it is important to consider the full set of individual element abundances, with respect to hydrogen, rather than the ratios of only two elements, such as the C/O ratio, particularly when comparing model predictions with observed transmission and emission spectra.
ABSTRACT We present a primary transit observation for the ultra-hot (Teq ∼ 2400 K) gas giant expolanet WASP-121b, made using the Hubble Space Telescope Wide Field Camera 3 in spectroscopic mode ...across the 1.12-1.64 m wavelength range. The 1.4 m water absorption band is detected at high confidence (5.4 ) in the planetary atmosphere. We also reanalyze ground-based photometric light curves taken in the B, r′, and z′ filters. Significantly deeper transits are measured in these optical bandpasses relative to the near-infrared wavelengths. We conclude that scattering by high-altitude haze alone is unlikely to account for this difference and instead interpret it as evidence for titanium oxide and vanadium oxide absorption. Enhanced opacity is also inferred across the 1.12-1.3 m wavelength range, possibly due to iron hydride absorption. If confirmed, WASP-121b will be the first exoplanet with titanium oxide, vanadium oxide, and iron hydride detected in transmission. The latter are important species in M/L dwarfs and their presence is likely to have a significant effect on the overall physics and chemistry of the atmosphere, including the production of a strong thermal inversion.
Aims. To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics ...and thermal profile. Methods. We coupled the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self-consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption, and crucially, scattering. We used fluxes directly obtained from the UM to produce synthetic spectral energy distributions and phase curves. Results. Our simulations show extensive cloud formation in both HD 209458 b and HD 189733 b. However, cooler temperatures in the latter result in higher cloud particle number densities. Large particles, reaching 1 μm in diameter, can form due to high particle growth velocities, and sub-μm particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, results in enhanced cloud coverage for latitudes above 30° and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet, for HD 209458 b and HD 189733 b, as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, the simulated optical phase curve for HD 209458 b shows no “offset”, in contrast to observations. Efficient scattering of stellar irradiation by cloud particles results in a local maximum cooling of up to 250 K in the upper atmosphere, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements, which have been employed to interpret observations and determine thermal profiles for these planets. This suggests that readers should be cautious of interpretations neglecting such cloud feedback and scattering, and that the subject merits further study.
ABSTRACT
The effect of enhanced UV irradiation associated with stellar flares on the atmospheric composition and temperature of gas giant exoplanets was investigated. This was done using a 1D ...radiative-convective-chemical model with self-consistent feedback between the temperature and the non-equilibrium chemistry. It was found that flare-driven changes to chemical composition and temperature give rise to prolonged trends in evolution across a broad range of pressure levels and species. Allowing feedback between chemistry and temperature plays an important role in establishing the quiescent structure of these atmospheres, and determines their evolution due to flares. It was found that cooler planets are more susceptible to flares than warmer ones, seeing larger changes in composition and temperature, and that temperature–chemistry feedback modifies their evolution. Long-term exposure to flares changes the transmission spectra of gas giant atmospheres; these changes differed when the temperature structure was allowed to evolve self-consistently with the chemistry. Changes in spectral features due to the effects of flares on these atmospheres can be associated with changes in composition. The effects of flares on the atmospheres of sufficiently cool planets will impact observations made with JWST. It is necessary to use self-consistent models of temperature and chemistry in order to accurately capture the effects of flares on features in the transmission spectra of cooler gas giants, but this depends heavily on the radiation environment of the planet.
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).
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.
ABSTRACT
Transport-induced quenching, i.e. the homogenization of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical ...modelling of this process exists, the three-dimensional nature of transport-induced quenching is underexplored. Here, we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b, HD 189733b, HD 209458b, and WASP-17b including coupled hydrodynamics, radiative transfer, and chemistry. Our simulations were performed with two chemical schemes: a chemical kinetics scheme, which is capable of capturing transport-induced quenching, and a simpler, more widely used chemical equilibrium scheme. We find that transport-induced quenching is predicted to occur in atmospheres of all planets in our sample; however, the extent to which it affects their synthetic spectra and phase curves varies from planet to planet. This implies that there is a ‘sweet spot’ for the observability of signatures of transport-induced quenching, which is controlled by the interplay between the dynamics and chemistry.