Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere-where temperature increases with altitude-these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5σ confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms.We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H 2 -H 2 , H 2 -He, H 2 O, CO, CO 2 , CH 4 , NH 3 , K, Na, and TiO, VO if they are present in the atmosphere.We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH 3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J − H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.

A correlation between giant-planet mass and atmospheric heavy elemental abundance was first noted in the past century from observations of planets in our own Solar System and has served as a cornerstone of planet-formation theory. Using data from the Hubble and Spitzer Space Telescopes from 0.5 to 5 micrometers, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent HO absorption bands with a maximum base-to-peak amplitude of 525 parts per million in the transmission spectrum. Using the water abundance as a proxy for metallicity, we measured HAT-P-26b's atmospheric heavy element content ([Formula: see text] times solar). This likely indicates that HAT-P-26b's atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals.

We are adapting the global circulation model (GCM) of the UK Met Office, the so-called unified model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large-scale, time-averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non-hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.

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 a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) 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.

We present a primary transit observation for the ultra hot (T eq ∼ 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 lightcurves 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. 8 We quote equilibrium temperature Teq as the blackbody temperature required for planetary thermal emission to balance the stellar irradiation, assuming zero Bond albedo and uniform daynight recirculation. Due to the approximate nature of equilibrium temperature, we round values to the nearest 100 K in this Letter.

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.

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 couple 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 use 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, result in enhanced cloud coverage for latitudes > 30 o 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 both a note of caution of interpretations neglecting such cloud feedback and scattering, and merits further study.