Atmospheric processes

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:

Marine plastic debris that is discarded into the ocean and eventually beaches is an acute problem for small island nations such as the Seychelles. To address this problem and enable anticipatory action, the sources of marine plastic debris need to be identified. Recent observations suggest that not all of the debris that arrives at the Seychelles is from terrestrial input. However, there is currently a lack of quantitative attribution of maritime debris sources.Therefore, we investigate the potential shipping and fishing vessel plastic debris sources in the southwestern Indian Ocean that beach at the Seychelles. We use a 2D Lagrangian particle-tracking model, based on OceanParcels, with particles that are advected by currents from a 1/50° (~2 km) regional ocean model. We combine this with satellite-tracked fishing and shipping data to inform particle starting locations and weightings. This model resolution allows us to resolve island to sub-island accumulation patterns.Spatially, model results suggest that debris beaching at the Seychelles from shipping vessels is concentrated along a limited number of high-activity shipping routes. The port destinations of these routes are consistent with the origin of plastic bottles inferred from labels in a previous study. Model results also show that 50-66% of fishing debris that beaches at the Seychelles is discarded within its own exclusive economic zone, depending on the island group.Temporally, the season during which debris is discarded strongly impacts the likelihood of beaching. This seasonal profile varies in amplitude and phase between islands due to wind-driven surface current changes. Additionally, we find debris beaching patterns can vary substantially between islands and on a sub-island scale. This highlights the importance of higher-resolution models for investigating plastic accumulation at kilometre-scale islands. We also assess the model against observations including bottle drifters and seasonal accumulation data at the Aldabra Atoll, the latter of which is consistent with sub-island scale accumulation seasonality from the model.Given that most marine debris originates from a few major shipping routes and from within the Seychelles exclusive economic zone, we suggest targeted enforcement of MARPOL Annex V could tackle the source of the issue.

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.