91̽»¨

Abstract:

Machine learning hold the promise of unlocking more accurate and realistic parameterizations of atmospheric processes, but brings its own set of challenges and drawbacks. Among top issues are generalization, stability and interpretability. Here we present a parameter-efficient neural network parameterization which aims to address these issues by incorporating physical knowledge to a high degree. By predicting fluxes and microphysical process rates instead of total tendencies, the conservation of water can be hardcoded, which is shown to improve online performance. Furthermore, a physically motivated architecture based on vertically recurrent neural networks enables high computational efficiency and a low number of parameters. The models are trained and evaluated using a superparameterization setup with real orography. The impact of incorporating stochasticity is also discussed.

Abstract:

Representing and quantifying uncertainty in physical parameterisations is a central challenge in weather and climate modelling, and approaches are often developed separately for different timescales. Here, we consider the separation of uncertainty by source using machine learning frameworks for subgrid-scale parameterisations. In this context, aleatoric uncertainty arises from internal variability in the training data, and epistemic uncertainty, arises from poorly constrained parameters during training. Using the Lorenz 1996 system as a testbed for simplified chaotic dynamics, we deal with uncertainties through a unified framework using Bayesian Neural Networks, to explore how the different sources of uncertainty evolve over different prediction timescales.

Abstract:

In the Model Uncertainty-Model Intercomparison Project (MUMIP) we compare parameterisation packages from different modelling centres using their single-column modelling (SCM) frameworks. We will showcase the dataset from an Indian Ocean experiment at a 0.2 degrees grid covering one month, with about 10 million simulations of each model. These parametrised models are compared against a convection-permitting benchmark from DYAMOND under common dynamical constraints. We will show differences and similarities in precipitation patterns and physics tendencies among four models and show how these differences can be generalised. Following earlier works, we find that at coarse grids that do not resolve convection, parameterisation packages tend to produce overconfident tendencies compared to the convection-permitting benchmark. Furthermore, we test several hypotheses on the MUMIP dataset to explain the differences. We use the data to explore the foundations of stochastic physical parametrisations. Would stochastic physics effectively overcome the overconfidence for good reasons? May the stochastic perturbations actually have a physically meaningful quantitative interpretation? Can stochastic physics be used to partially overcome truncation and grid spacing limitations?

Abstract:

The Atlantic Multidecadal Variability (AMV) and the North Atlantic Oscillation (NAO) are the dominant modes of oceanic and atmospheric variability in the North Atlantic, respectively, and are key sources of predictability from seasonal to decadal timescales. However, the physical processes and feedback mechanisms linking the AMV and NAO, and the role of diabatic processes in these feedbacks, remain debated. We present a data-driven dynamical modelling framework which captures coupled decadal variability in AMV, NAO, and North Atlantic precipitation. Applying equation discovery methods to observational data, we identify deterministic low-order dynamical models consisting of three coupled ordinary differential equations. These models reproduce observed North Atlantic decadal variability and show robust out-of-sample predictive skill on multi-annual to decadal lead times. The resulting model dynamics include a distinct quasi-periodic 20-year oscillation consistent with a damped oceanic mode of variability. Notably, precipitation-related terms feature prominently in the low-order models, suggesting an important role for latent heat release and freshwater fluxes in mediating ocean–atmosphere interactions. We propose new feedback mechanisms between North Atlantic sea surface temperature and the NAO, with precipitation acting as a dynamical bridge. By linearising the low-order models and computing finite-time Lyapunov exponents, we find that North Atlantic precipitation is more predictable in a positive AMV phase. We then analyse several decadal prediction ensemble experiments based on initialised hindcasts and find comparable state-dependent predictability of precipitation. Overall, these results illustrate how data-driven equation discovery can provide mechanistic hypotheses and new insight beyond conventional analyses of observations and climate model simulations.

Abstract:

Uncertainty in projections of future regional climate change remains large, driven by structural differences among Earth System Models and the influence of internal climate variability. Existing uncertainty-reduction approaches, including emergent constraints and Bayesian variants, primarily focus on forced climate responses derived from simple aggregate metrics, thereby requiring strong assumptions and exploiting only low-dimensional climate information. Here we propose a data-driven deep-learning framework that directly forecasts spatially and monthly resolved decadal mean climatologies of surface temperature anomalies from the 2030s to the 2090s, using only recent monthly trajectories spanning 1980-2025. The training ensemble contains 265 historical+SSP2-4.5 simulations, distributed across 40 ESMs from 25 different families (i.e., modelling centers) over which the cross validation is performed. The architecture couples pluri-annual to multi-decadal temporal convolutions with a spatial U-Net encoder-decoder and is evaluated on CMIP6 simulations using a leave-one-model-family-out cross-validation (LOMFO-CV) design to ensure generalisation across separately developed ESMs. Predictive uncertainty is quantified via LOMFO-CV errors, yielding conservative and reliable ranges that incorporate irreducible internal variability and systematic model shifts.To further evaluate the predictive capacity beyond the CMIP6 distribution, we evaluated the network on historical+SSP2-4.5 simulations from a recent HadGEM3-GC5 model hierarchy developed within the European Eddy-Rich ESMs (EERIE) project, the European contribution to HighResMIP2 for CMIP7. In particular, the eddy-rich GC5-HH configuration explicitly simulates mesoscale ocean dynamics that are absent in CMIP6-type models, providing a rigorous test of generalisation to richer and more realistic physical representations. Despite these substantial differences, the network successfully reproduces warming trajectories and future climate patterns for all three model configurations (GC5-LL, GC5-MM, GC5-HH), with forecast errors largely contained within empirically calibrated uncertainty bounds from the LOMFO-CV, both globally and locally. These results, notably for GC5-HH and its more realistic physics, strengthens confidence in the applicability of the framework to real-world data.When applied to observations, the extracted end-of-century global-mean surface temperature and its uncertainty range are consistent with prior estimates from Bayesian frameworks. At local scales, the network reduces uncertainty by 40% (2030s) to 30% (2090s) on average, and by up to 75% in some regions for all future decades. Importantly, these uncertainty estimates account not only for uncertainty in the forced response (as emergent constraint methods do), but also for errors associated with predicting different realisations of internal variability, providing a physically meaningful reduction of local and global climate uncertainty.