Convective processes are crucial in shaping exoplanetary atmospheres but are computationally expensive to simulate directly. A novel technique of simulating moist convection, especially on tidally locked exoplanets such as those orbiting TRAPPIST-1, is to use a 3D general circulation model (GCM) with a global stretched mesh. This allows us to locally refine the model resolution to a km-scale and resolve deep convection without relying on parameterization. We explore the impact of explicit vs parameterized convection on the climate of TRAPPIST-1e, a confirmed rocky exoplanet in the habitable zone and a primary candidate for atmospheric characterization. We show allowing for explicit convection in a stretched-mesh simulation results primarily in changes in cloud distribution and precipitation on a planetary scale. Nevertheless, the overall climate state is close to that produced with parameterized convection and a non-stretched mesh. Additionally, these novel simulations shed more light on the bistability of the atmospheric circulation on TRAPPIST-1e. Our methodology opens an exciting and computationally feasible avenue for improving our understanding of fine-scale 3D mixing in exoplanetary atmospheres.

The moist shallow water equations offer a promising route for advancing understanding of the coupling of physical parametrisations and dynamics in numerical atmospheric models, an issue known as 'physics-dynamics coupling'. Without moist physics, the traditional shallow water equations are a simplified form of the atmospheric equations of motion and so are computationally cheap, but retain many relevant dynamical features of the atmosphere. Introducing physics into the shallow water model in the form of moisture provides a tool to experiment with numerical techniques for physics-dynamics coupling in a simple dynamical model. In this paper, we compare some of the different moist shallow water models by writing them in a general formulation. The general formulation encompasses three existing forms of the moist shallow water equations and also a fourth, previously unexplored formulation. The equations are coupled to a three-state moist physics scheme that interacts with the resolved flow through source terms and produces two-way physics-dynamics feedback. We present a new compatible finite element discretisation of the equations and apply it to the different formulations of the moist shallow water equations in three test cases. The results show that the models capture generation of cloud and rain and physics-dynamics interactions, and demonstrate some differences between moist shallow water formulations and the implications of these different modelling choices.

Convective processes are crucial in shaping exoplanetary atmospheres but are computationally expensive to simulate directly. A novel technique of simulating moist convection on tidally locked exoplanets is to use a global 3D model with a stretched mesh. This allows us to locally refine the model resolution to 4.7 km and resolve fine-scale convective processes without relying on parameterizations. We explore the impact of mesh stretching on the climate of a slowly rotating TRAPPIST-1e-like planet, assuming it is 1:1 tidally locked. In the stretched-mesh simulation with explicit convection, the climate is 5 K colder and 25% drier than that in the simulations with parameterized convection(with both stretched and quasi-uniform meshes). This is due to the increased cloud reflectivity—because of an increase in low-level cloudiness—and exacerbated by the diminished greenhouse effect due to less water vapor. At the same time, our stretched-mesh simulations reproduce the key characteristics of the global climate of tidally locked rocky exoplanets, without any noticeable numerical artifacts. Our methodology opens an exciting and computationally feasible avenue for improving our understanding of 3D mixing in exoplanetary atmospheres. Our study also demonstrates the feasibility of a global stretched-mesh configuration for LFRic-Atmosphere, the next-generation Met Office climate and weather model.

Local conservation of mass and entropy are becoming increasingly desirable properties for modern numerical weather and climate models. This work presents a Flux-Form Semi-Lagrangian (FFSL) transport scheme, called SWIFT, that facilitates this conservation for tracer variables, whilst maintaining other vital properties such as preservation of a constant, monotonicity and positivity. Importantly, these properties all hold for large Courant numbers and multi-dimensional flow, making the scheme appropriate for use within a dynamical core which takes large time steps. The SWIFT scheme presented here can be seen as an evolution of the FFSL methods of Leonard et al and Lin and Rood. Two-dimensional and three-dimensional schemes consist of a splitting into a sequence of one-dimensional calculations. The new SWIFT splitting presented here allows monotonic and positivity properties from the one-dimensional calculations to be inherited by the multi-dimensional scheme. These one-dimensional calculations involve separating the mass flux into terms that correspond to integer and fractional parts of the Courant number. Key to achieving conservation is coupling the transport of tracers to the transport of the fluid density, through re-use of the discrete mass flux that was calculated from the fluid density in the transport of the tracers. This work also describes how these properties can still be attained when the tracer is vertically-staggered from the density in a Charney-Phillips grid.

The reformulation of the Met Office's dynamical core for weather and climate prediction previously described by the authors is extended to spherical domains using a cubed-sphere mesh. This paper updates the semi-implicit mixed finite-element formulation to be suitable for spherical domains. In particular the finite-volume transport scheme is extended to take account of non-uniform, non-orthogonal meshes and uses an advective-then-flux formulation so that increment from the transport scheme is linear in the divergence. The resulting model is then applied to a standard set of dry dynamical core tests and compared to the existing semi-implicit semi-Lagrangian dynamical core currently used in the Met Office's operational model.

We demonstrate that LFRic-Atmosphere, a model built using the Met Office's GungHo dynamical core, is able to reproduce idealised large-scale atmospheric circulation patterns specified by several widely-used benchmark recipes. This is motivated by the rapid rate of exoplanet discovery and the ever-growing need for numerical modelling and characterisation of their atmospheres. Here we present LFRic-Atmosphere's results for the idealised tests imitating circulation regimes commonly used in the exoplanet modelling community. The benchmarks include three analytic forcing cases: the standard Held-Suarez test, the Menou-Rauscher Earth-like test, and the Merlis-Schneider Tidally Locked Earth test. Qualitatively, LFRic-Atmosphere agrees well with other numerical models and shows excellent conservation properties in terms of total mass, angular momentum and kinetic energy. We then use LFRic-Atmosphere with a more realistic representation of physical processes (radiation, subgrid-scale mixing, convection, clouds) by configuring it for the four TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) scenarios. This is the first application of LFRic-Atmosphere to a possible climate of a confirmed terrestrial exoplanet. LFRic-Atmosphere reproduces the THAI scenarios within the spread of the existing models across a range of key climatic variables. Our work shows that LFRic-Atmosphere performs well in the seven benchmark tests for terrestrial atmospheres, justifying its use in future exoplanet climate studies.

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds. The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. An approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

The Charney–Phillips grid, used in many numerical models of the atmosphere, involves vertically staggering the nodes of the density variable with the nodes of the entropy‐type variable. When moisture is included in such a model, it is either co‐located with density so that moisture can be transported conservatively and consistently with dry mass, or with the entropy‐type variable so that the coupling between moisture and temperature can be represented well. Both properties are desirable, yet at first it appears difficult to obtain both simultaneously. Here, we present a framework to resolve this problem, by co‐locating the moisture mixing ratio with potential temperature but formulating its transport as that of a density on a vertically shifted mesh. Within this framework, particular choices of the operators involved provide the desired conservation and consistency properties of the moisture transport. The framework is described in the context of a finite‐element approach. We also present an explicit Runge–Kutta time‐stepping scheme that is appropriate for use within this framework. This approach is then illustrated through numerical tests, which demonstrate that it does indeed have the desired conservation and consistency properties.

A famous feature of the Camassa–Holm equation is its admission of peaked soliton solutions known as peakons. We investigate this equation under the influence of stochastic transport. Noting that peakons are weak solutions of the equation, we present a finite-element discretization for it, which we use to explore the formation of peakons. Our simulations using this discretization reveal that peakons can still form in the presence of stochastic perturbations. Peakons can emerge both through wave breaking, as the slope turns vertical, and without wave breaking as the inflection points of the velocity profile rise to reach the summit.

A promising development of the last decade in the numerical modelling of geophysical fluids has been the compatible finite‐element framework. Indeed, this will form the basis for the next‐generation dynamical core of the Met Office. For this framework to be useful for numerical weather prediction models, it must be able to handle descriptions of unresolved and diabatic processes. These processes offer a challenging test for any numerical discretisation, and have not yet been described within the compatible finite‐element framework. The main contribution of this article is to extend a discretisation using this new framework to include moist thermodynamics. Our results demonstrate that discretisations within the compatible finite‐element framework can be robust enough also to describe moist atmospheric processes.We describe our discretisation strategy, including treatment of moist processes, and present two configurations of the model using different sets of function spaces with different degrees of finite element. The performance of the model is demonstrated through several test cases. Two of these test cases are new cloudy‐atmosphere variants of existing test cases: inertia–gravity waves in a two‐dimensional vertical slice and a three‐dimensional rising thermal.

We present a new compatible finite element advection scheme for the compressible Euler equations. Unlike the discretisations described in Cotter and Kuzmin (2016) and Shipton et al (2018), the discretisation uses the lowest-order family of compatible finite element spaces, but still retains second-order numerical accuracy. This scheme obtains this second-order accuracy by first `recovering' the function in higher-order spaces, before using the discontinuous Galerkin advection schemes of Cotter and Kuzmin (2016). As well as describing the scheme, we also present its stability properties and a strategy for ensuring boundedness. We then demonstrate its properties through some numerical tests, before presenting its use within a model solving the compressible Euler equations.

A framework of variational principles for stochastic fluid dynamics was presented by Holm (2015), and these stochastic equations were also derived by Cotter et al. (2017). We present a conforming finite element discretisation for the stochastic quasi-geostrophic equation that was derived from this framework. The discretisation preserves the first two moments of potential vorticity, i.e. the mean potential vorticity and the enstrophy. Following the work of Dubinkina and Frank (2007), who investigated the statistical mechanics of discretisations of the deterministic quasi-geostrophic equation, we investigate the statistical mechanics of our discretisation of the stochastic quasi-geostrophic equation. We compare the statistical properties of our discretisation with the Gibbs distribution under assumption of these conserved quantities, finding that there is agreement between the statistics under a wide range of set-ups.

Refractivity of water vapor at THz frequencies was calculated using three different models. These did not agree. A direct measurement was made using time-domain spectroscopy and the correct form of the refractivity model was determined.