Atmospheric processes

Abstract:

Mesoscale Convective Systems (MCSs), with length scales of 100 to 1000 km or more, fall into the "grey zone" of global models with grid spacings of 10s of km. Their under-resolved nature leads to model deficiencies in representing MCS latent heating, whose vertical structure critically shapes large-scale circulations. To address this challenge, we use analysis increments—the corrections applied by Data Assimilation (DA) to the model's prior state—from a 10 km Met Office operational forecast model to inform the development of a stochastic parameterization for MCS latent heating. To focus on errors in MCS feedback rather than errors due to a missing MCS, we select analysis increments from 1037 MCS tracks that the model successfully captures at the start of the DA cycle.A Machine Learning–based Gaussian Mixture Model reveals that the vertical structure of temperature analysis increments is probabilistically linked to the atmospheric environment. Bottom-heavy heating increments tend to occur in low Total Column Water Vapor (TCWV) conditions, suggesting that the model underestimates low-level convective heating in relatively dry environments. In contrast, top-heavy heating increments are linked to a moist layer overturning structure—characterized by high TCWV and strong vertical wind shear—indicating model underestimation of upper-level condensate detrainment in such environments. This probabilistic relationship is implemented in the Met Office operational forecast model as part of the MCS: PRIME stochastic scheme, which corrects MCS-related uncertainties during model integration. By enhancing top-heavy heating, the scheme backscatters kinetic energy from the mesoscale to larger scales, improving predictions of Indian seasonal rainfall and the Madden–Julian Oscillation (MJO). Future work will assess its impact on forecast busts and its potential to extend predictability.

Abstract:

A parameterisation suite is the combination of all parameterisation schemes that is used by a numerical model of the atmosphere. These parameterisation (or “physics”) suites are widely seen as the most uncertain components of atmospheric models.  In MUMIP we compare deterministic parameterisation suites from across different modelling centres under common prescribed large-scale dynamics. In the first MUMIP experiment, these dynamical tendencies have been derived by coarse-graining the convection-permitting ICON DYAMOND simulation to 0.2 degree resolution. We use these realistic spatiotemporal dynamical patterns to drive millions  of single column model simulations over the tropical Indian Ocean with prescribed SSTs. We use this data to estimate the uncertainty from their physics across four models, each using their default convection-parametrised physics suites. The models are: IFS, GFS, RAP and ARPEGE. The distributions of precipitation rate, convective available potential energy (CAPE), convective inhibition (CIN) and level of neutral buoyancy are analysed, as well as individual model tendencies and rate of change of CAPE and CIN as a function of lead time and, for instance, the diurnal cycle . We find notable differences across the physics suites and even more strongly between convection-parameterised physics suites and the convection-permitting ICON DYAMOND benchmark. Furthermore, we relate these diagnostics to biases in temperature and specific humidity. We also develop a framework for the detection of statistical relations among diagnostics and/or their change. The framework may for instance be used to quantify the impact of spin-up compared to persistence ("memory") and randomness within a dataset and to identify similarity in the physics across modelling centres. In this contribution some of the early results of the international MUMIP project will be presented and we hope to encourage other researchers to use and/or complement the data of MUMIP. Please refer to https://mumip.web.ox.ac.uk for details of how to get involved.   

Abstract:

Machine learning has the potential to improve the physical realism and/or computational efficiency of parameterizations. A typical approach has been to feed concatenated vertical profiles to a dense neural network. However, feed‐forward networks lack the connections to propagate information sequentially through the vertical column. Here we examine if predictions can be improved by instead traversing the column with recurrent neural networks (RNNs) such as Long Short‐Term Memory (LSTMs). This method encodes physical priors (locality) and uses parameters more efficiently. Firstly, we test RNN‐based radiation emulators in the Integrated Forecasting System. We achieve near‐perfect offline accuracy, and the forecast skill of a suite of global weather simulations using the emulator are for the most part statistically indistinguishable from reference runs. But can radiation emulators provide both high accuracy and a speed‐up? We find optimized, state‐of‐the‐art radiation code on CPU generally faster than RNN‐based emulators on GPU, although the latter can be more energy efficient. To test the method more broadly, and explore recent challenges in parameterization, we also adapt it to data sets from other studies. RNNs outperform reference feed‐forward networks in emulating gravity waves, and when combined with horizontal convolutions, for non‐local unified parameterization. In emulation of moist physics with memory, the RNNs have similar offline accuracy as ResNets, the previous state‐of‐the‐art. However, the RNNs are more efficient, and more stable in autoregressive semi‐prognostic tests. Multi‐step autoregressive training improves performance in these tests and enables a latent representation of convective memory. Recently proposed linearly recurrent models achieve similar performance to LSTMs.

Abstract: