Search

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., Comment: 26 Pages, 8 Figures, 1 Table

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., Comment: 34 pages, 9(12) figures; Submitted to Geoscientific Model Development; Comments are welcome (see Discussion tab on the journal's website: https://egusphere.copernicus.org/preprints/2023/egusphere-2023-647)

We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow., Comment: 26 pages, 22 Figures. Accepted for publication in Astronomy and Astrophysics

Geophysical models of the atmosphere and ocean invariably involve parameterizations. These represent two distinct areas: Subgrid processes that the model cannot resolve, and diabatic sources in the equations, due to radiation for example. Hence, coupling between these physics parameterizations and the resolved fluid dynamics and also between the dynamics of the air and water, is necessary. In this paper weather and climate models are used to illustrate the problems. Nevertheless the same applies to other geophysical models. This coupling is an important aspect of geophysical models. However, often model development is strictly segregated into either physics or dynamics. As a consequence, this area has many unanswered questions. Recent developments in the design of dynamical cores, extended process physics and predicted future changes of the computational infrastructure are increasing complexity. This paper reviews the state-of-the-art of the physics-dynamics coupling in geophysical models, surveys the analysis techniques, and illustrates open questions in this field. This paper focuses on two objectives: To illustrate the phenomenology of the coupling problem with references to examples in the literature and to show how the problem can be analysed. Proposals are made on how to advance the understanding and upcoming challenges with emerging modeling strategies. This paper is of interest to model developers who aim to improve the models and have to make choices on and test new implementations, to users who have to understand choices presented to them and finally users of outputs, who have to distinguish physical features from numerical problems in the model data.

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., Comment: 24 pages, 11 figures. Accepted for publication in Astronomy and Astrophysics

We demonstrate that both the current (New Dynamics), and next generation (ENDGame) dynamical cores of the UK Met Office global circulation model, the UM, reproduce consistently, the long-term, large-scale flows found in several published idealised tests. The cases presented are the Held-Suarez test, a simplified model of Earth (including a stratosphere), and a hypothetical tidally locked Earth. Furthermore, we show that using simplifications to the dynamical equations, which are expected to be justified for the physical domains and flow regimes we have studied, and which are supported by the ENDGame dynamical core, also produces matching long-term, large-scale flows. Finally, we present evidence for differences in the detail of the planetary flows and circulations resulting from improvements in the ENDGame formulation over New Dynamics., Comment: 34 Pages, 23 Figures. Accepted for publication in Geoscientific Model Development (pre-proof version)