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.

The majority of potentially habitable exoplanets detected orbit stars cooler than the Sun, and therefore are irradiated by a stellar spectrum peaking at longer wavelengths than that incident on Earth. Here, we present results from a set of simulations of tidally-locked terrestrial planets orbiting three different host stars to isolate the effect of the stellar spectra on the simulated climate. Specifically, we perform simulations based on TRAPPIST-1e, adopting an Earth-like atmosphere and using the UK Met Office Unified Model in an idealised 'aqua-planet' configuration. Whilst holding the planetary parameters constant, including the total stellar flux (900 W/m$^2$) and orbital period (6.10 Earth days), we compare results between simulations where the stellar spectrum is that of a quiescent TRAPPIST-1, Proxima Centauri and the Sun. The simulations with cooler host stars had an increased proportion of incident stellar radiation absorbed directly by the troposphere compared to the surface. This, in turn, led to an increase in the stability against convection, a reduction in overall cloud coverage on the dayside (reducing scattering), leading to warmer surface temperatures. The increased direct heating of the troposphere also led to more efficient heat transport from the dayside to the nightside and, therefore, a reduced day-night temperature contrast. We inferred that planets with an Earth-like atmosphere orbiting cooler stars had lower dayside cloud coverage, potentially allowing habitable conditions at increased orbital radii, compared to similar planets orbiting hotter stars for a given planetary rotation rate.

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 couple 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 use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate 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. 3D 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, detectable with near-future instruments.

Using a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidally-locked Earth-like exoplanets, Trappist-1e and Proxima Centauri b, to the choice of a convection parameterization. Compared to a mass-flux convection parameterization, a simplified convection adjustment parameterization leads to a $>$60% decrease of the cloud albedo, increasing the mean day-side temperature by $\approx$10 K. The representation of convection also affects the atmospheric conditions of the night side, via a change in planetary-scale wave patterns. As a result, using the convection adjustment scheme makes the night-side cold traps warmer by 17-36 K for the planets in our simulations. The day-night thermal contrast is sensitive to the representation of convection in 3D GCM simulations, so caution should be taken when interpreting emission phase curves. The choice of convection treatment, however, does not alter the simulated climate enough to result in a departure from habitable conditions, at least for the atmospheric composition and planetary parameters used in our study. The near-surface conditions both in the Trappist-1e and Proxima b cases remain temperate, allowing for an active water cycle. We further advance our analysis using high-resolution model experiments, in which atmospheric convection is simulated explicitly. Our results suggest that in a hypothetical global convection-permitting simulation the surface temperature contrast would be higher than in the coarse-resolution simulations with parameterized convection. In other words, models with parameterized convection may overestimate the inter-hemispheric heat redistribution efficiency.

We use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman mechanism and the hydrogen oxide (HO$_x$, describing the sum of OH and HO$_2$) catalytic cycle. We find that Proxima Centauri radiates with sufficient UV energy to initialize the Chapman mechanism. The result is a thin but stable ozone layer that peaks at 0.75 parts per million at 25 km. The quasi-stationary distribution of atmospheric ozone is determined by photolysis driven by incoming stellar radiation and by atmospheric transport. Ozone mole fractions are smallest in the lowest 15 km of the atmosphere at the sub-stellar point and largest in the nightside gyres. Above 15 km the ozone distribution is dominated by an equatorial jet stream that circumnavigates the planet. The nightside ozone distribution is dominated by two cyclonic Rossby gyres that result in localized ozone hotspots. On the dayside the atmospheric lifetime is determined by the HO$_x$ catalytic cycle and deposition to the surface, with nightside lifetimes due to chemistry much longer than timescales associated with atmospheric transport. Surface UV values peak at the substellar point with values of 0.01 W/m$^2$, shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds.

We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism can not be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed $f_{\textrm{sed}}$ in EddySed), but also the temperature-pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parameterised cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments, is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission and 4.5 $\mu$m Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.

We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called `shallow-fluid' and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the day and night sides. The implications for the atmospheric chemistry and, therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets.

In this paper we present three-dimensional 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 $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional 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 \textmu m Spitzer/IRAC channel.