Microphysical controls on the stratocumulus topped boundary-layer structure during VOCALS-REx

Simulations at a range of resolutions are compared to observations from the South-East Pacific taken during VOCALS-REx. It is found that increased horizontal and vertical resolution make only small improvements to the bulk properties of the simulated cloud and drizzle, but the highest resolution simulation is able to realistically represent mesoscale features in the cloud field. We focus on the highest resolution simulation and demonstrate that a poor representation of the cloud microphysics results in excessive drizzle production. This promotes persistent drizzle induced decoupling of the boundary layer, giving a poor representation of the observed diurnal cycle of stratocumulus. Two simple changes to the microphysics scheme are implemented: a modified autoconversion parametrization and a new representation of the rain drop size distribution. This results in a more realistic simulation of boundary-layer diurnal decoupling, and improvements to the cloud liquid water path and surface drizzle rate

Friction in mid-latitude cyclones: an Ekman-PV mechanism

The mechanism by which the atmospheric boundary layer reduces the intensity of mid-latitude cyclones is investigated. It is demonstrated that two alternative theories, Ekman pumping and the baroclinic potential vorticity (PV) mechanism, in fact act in union to maximize the spin-down. Ekman pumping aids the ventilation of PV from the boundary layer, and shapes the resulting PV anomaly into one of increased static stability. PV inversion techniques are used to demonstrate how this anomaly reduces the coupling between the upper- and lower-levels within the cyclone, reducing the growth rate

The Limits of the Primitive Equations of Dynamics for Warm, Slowly Rotating Small Neptunes and Super Earths (article)

This is the author accepted manuscript. The final version is available from American Astronomical Society / IOP Publishing via the DOI in this record.The dataset associated with this article is located in ORE at: https://doi.org/10.24378/exe.1023We 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.Leverhulme TrustScience and Technology Facilities CouncilEuropean Research Counci

Exploring the climate of Proxima B with the Met Office Unified Model

This is the author accepted manuscript. The final version is available from EDP Sciences via the DOI in this record.The corrigendum to this article is in ORE at: http://hdl.handle.net/10871/34331We 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).I.B., J.M. and P.E. acknowledge the support of a Met Office Academic Partnership secondment. B.D. thanks the University of Exeter for support through a Ph.D. studentship. N.J.M. and J.G.’s contributions were in part funded by a Leverhulme Trust Research Project Grant, and in part by a University of Exeter College of Engineering, Mathematics and Physical Sciences studentship. We acknowledge use of the MONSooN system, a collaborative facility supplied under the Joint Weather and Climate Research Programme, a strategic partnership between the Met Office and the Natural Environment Research Council. This work also used the University of Exeter Supercomputer, a DiRAC Facility jointly funded by STFC, the Large Facilities Capital Fund of BIS and the University of Exeter

Ozone chemistry on tidally locked M dwarf planets

This is the final version. Available from Oxford University Press via the DOI in this recordWe 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 (HOx, describing the sum of OH and HO2) 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 HOx 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/m2 , shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds.Leverhulme TrustScience and Technology Facilities Council (STFC

A modern-day Mars climate in the Met Office Unified Model: dry simulations

We present results from the Met Office Unified Model (UM), a world-leading climate and weather model, adapted to simulate a dry Martian climate. We detail the adaptation of the basic parameterisations and analyse results from two simulations, one with radiatively active mineral dust and one with radiatively inactive dust. These simulations demonstrate how the radiative effects of dust act to accelerate the winds and create a mid-altitude isothermal layer during the dusty season. We validate our model through comparison with an established Mars model, the Laboratoire de Météorologie Dynamique planetary climate model (PCM), finding good agreement in the seasonal wind and temperature profiles but with discrepancies in the predicted dust mass mixing ratio and conditions at the poles. This study validates the use of the UM for a Martian atmosphere, highlights how the adaptation of an Earth general circulation model (GCM) can be beneficial for existing Mars GCMs and provides insight into the next steps in our development of a new Mars climate model.</p

Overcast on Osiris: 3D radiative-hydrodynamical simulations of a cloudy hot Jupiter using the parametrized, phase-equilibrium cloud formation code EDDYSED (article)

This is the final version. Available from OUP via the DOI in this recordThe dataset associated with this article is available in ORE: https://doi.org/10.24378/exe.1483We 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 cannot be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed fsed 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 parametrized 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 μm Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.Leverhulme TrustScience and Technology Facilities Council (STFC

Observable signatures of wind-driven chemistry with a fully consistent three dimensional radiative hydrodynamics model of HD 209458b (article)

This is the final version of the article. Available from American Astronomical Society / IOP Publishing via the DOI in this record.The dataset associated with this article is located in ORE at: http://hdl.handle.net/10871/32579We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD 209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model–observation discrepancy in the 4.5 µm Spitzer/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominent.BD and DKS acknowledge funding from the European Research Council (ERC) under the European Unions Seventh Framework Programme (FP7/2007-2013) / ERC grant agreement no. 336792. NJM is part funded by a Leverhulme Trust Research Project Grant. JM and IAB acknowledge the support of a Met Office Academic Partnership secondment. ALC is funded by an STFC studentship. DSA acknowledges support from the NASA Astrobiology Program through the Nexus for Exoplanet System Science. This work used the DiRAC Complexity system, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility. This equipment is funded by BIS National E-Infrastructure capital grant ST/K000373/1 and STFC DiRAC Operations grant ST/K0003259/1. DiRAC is part of the National E-Infrastructure

Atmospheric convection plays a key role in the climate of tidally-locked terrestrial exoplanets: insights from high-resolution simulations

This is the final version. Available from IOP Piblishing via the DOI in this recordUsing a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidallylocked 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 ≈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.Leverhulme TrustScience and Technology Facilities Council (STFC

Simulating the cloudy atmospheres of HD 209458 b and HD 189733 b with the 3D Met Office Unified Model

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 anddynamics, 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 HD209458 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.PostprintPeer reviewe