How different are deterministic physics suites when coupled to fixed model dynamics and why?

Journal of the European Meteorological Society Elsevier 5 (2026) 100041

Authors:

Edward Groot, Hannah Christensen, Xia Sun, Kathryn Newman, Wahiba Lfarh, Romain Roehrig, Lisa Bengtsson, Julia Simonson

Abstract:

It is often difficult to attribute uncertainty and errors in atmospheric models to designated model components. This is because sub-grid parameterised processes interact strongly with the large-scale transport represented by the explicit model dynamics. We carry out experiments with prescribed large-scale dynamics and different sub-grid physics suites. This dataset has been constructed for the Model Uncertainty Model Intercomparison Project (MUMIP), in which each suite forecasts sub-grid tendencies at a 22km grid. The common dynamics is derived from a convection-permitting benchmark: an ICON DYAMOND experiment (2.5km grid). We compare four different physics suites for atmospheric models in an Indian Ocean experiment. We analyse their joint PDFs of precipitation and associated physics tendencies for a full month, where precipitation is used to diagnose uncertainty of convective activity. We find that all physics suites produce very similar precipitation amounts, with very high correlations between models, i.e.,  > 0.95 at the native grid. However, the convection-permitting benchmark is more dissimilar from each of the physics suites, with correlations of  ≈ 0.80. Similarly, we show that the vertically averaged physics tendencies in the free-troposphere are highly similar between the four physics suites, yet different if reconstructed for the benchmark. The water vapour sink is very closely linked with precipitation in the four physics suites. This suggests that the coarse-grid models are overconfident. We hypothese is that variation in unresolved convective structures can lead to variation in the dynamics, following a given amount of latent heating at fine grids, but not in our physics suites. The difference appears to be caused by the explicit interactions between gravity waves, occurring at fine grids only Groot et al. (2024) We assess whether their non-linear feedback from convective precipitation systems explains our joint PDFs of precipitation. The slightly exponential curve supports the interaction mechanism. These findings are further evidence for a non-linear feedback between convective organisation/aggregation and dynamics. This feedback has been studied earlier in a real-case study with ICON by looking from the fixed-physics rather than the fixed-dynamics perspective. Our current results may indicate that sub-grid physics with stochastic physics perturbations emulate convective organisation effects.

The impact of stochastic sea ice perturbations on seasonal forecasts

Weather and Climate Dynamics 7:3 (2026) 1593-1618

Authors:

K Strommen, M Mayer, A Storto, J Spaeth, S Tietsche

Abstract:

Sea ice ensemble forecasts can be highly underdispersive, meaning that the ensemble spread is notably lower than the average forecast error. One common strategy to address underdispersion is to add stochastic perturbations to the forecasts. We detail the implementation of a stochastically perturbed parameterisation (SPP) scheme for SI3, the sea ice component used by the Integrated Forecast System (IFS), the forecast model used and developed by the European Centre for Medium-Range Weather Forecasts (ECMWF). We then evaluate its impact on seasonal forecasts of Northern Hemisphere summer and winter. The inclusion of SPP is found to enhance ensemble spread for sea ice concentration (SIC) and sea ice thickness (SIT) forecasts by around 10 % relative to a forecast with no SPP, which results in a better calibrated probabilistic forecast. Some small but robust changes to the mean state are also found, including a general decrease in the mean SIC and a redistribution of the winter ice from the central Arctic to the ice edge. These changes reduce or increase the mean bias depending on the region. Changes to the mean and spread of the sea ice result in changes to the mean and spread of air temperature up to at least 850 hPa, altering the mean air temperature biases of the model. An apparent consequence of this is a significant increase in the anomaly correlation coefficients of 500 hPa geopotential height (Z500) over the Euro-Atlantic domain in winter, which partially projects onto the North Atlantic Oscillation. We conclude that sea ice stochastic perturbations can be a valuable contribution to increased reliability of seasonal forecasts of the sea ice itself and can impact seasonal forecasts of the atmosphere at high and mid latitudes.

How Do AI Climate Models Respond to Warming Across Climate Zones?

(2026)

Authors:

Charlotte C Merchant, Milan Klöwer, Bradley Stanley-Clamp, Maren Höver, Simon LL Michel, Edward Groot, Hannah M Christensen

New insights into decadal climate variability in the North Atlantic revealed by data-driven dynamical models

Earth Syst. Dynam., 17, 1061–1079, 2026

Authors:

Andrew J. Nicoll, Hannah M. Christensen, Chris Huntingford, and Doug Smith

Abstract:

The Atlantic Multidecadal Variability (AMV) and the North Atlantic Oscillation (NAO) are the dom-
inant 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 low-order
models consisting of three coupled ordinary differential equations. These models reproduce observed 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 pro-
pose new feedback mechanisms between North Atlantic sea surface temperature and the NAO, with precipitation
acting as a dynamical bridge. Overall, these results illustrate how equation discovery can provide mechanistic
hypotheses and new insight beyond conventional analyses of observations and climate model simulations.

New insights into decadal climate variability in the North Atlantic revealed by data-driven dynamical models

Earth System Dynamics Copernicus Publications 17:4 (2026) 1061-1079

Authors:

Andrew J Nicoll, Hannah M Christensen, Chris Huntingford, Doug Smith

Abstract:

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 low-order models consisting of three coupled ordinary differential equations. These models reproduce observed 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. Overall, these results illustrate how equation discovery can provide mechanistic hypotheses and new insight beyond conventional analyses of observations and climate model simulations.