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.

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 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.

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.