Unravelling the role of increased model resolution in simulating surface temperature fields using explainable AI

Artificial Intelligence for the Earth Systems, In press (2026)

Authors:

Simon L. L. Michel, Kristian Strommen, Jürgen Kröger, Malcolm J. Roberts, Jin-Song von Storch, and Hannah M. Christensen

Abstract:

Minimizing climate model biases is essential to reduce uncertainties in future climate projections. Despite recent advances, improvements between generations of Earth System Models (ESMs) remain modest, largely due to continued reliance on subgrid-scale parametrizations. These are required because CMIP6 model resolutions are too coarse to explicitly simulate small-scale processes such as ocean mesoscale eddies and deep atmospheric convection, which strongly influence climate patterns. Recent gains in computational power have enabled higher-resolution ESMs that can resolve some of these processes, thereby reducing the need for parametrization. However, robustly detecting improved modeling with increased resolution remains challenging due to internal variability and model-to-model discrepancies. This study employs a convolutional neural network (CNN) classifier combined with explainable AI (XAI) to assess the role of resolution in simulating winter surface temperature fields in an ensemble of 17 control climate simulations with varying oceanic and atmospheric resolutions. The CNN distinguishes between ESMs of varying resolutions using temperature snapshots, while XAI identifies the regions driving these classifications, offering deeper insight into ESM behavior. The results show that ESMs with similar ocean grid resolutions are more often confused with each other by the CNN than those from different modeling centers, highlighting the central role of ocean resolution, particularly mesoscale eddies, in shaping climate simulations. Although limited to surface air temperature, the approach provides a more nuanced perspective on ESM differences and performance than traditional bias analyses. The framework can be extended to other variables and ESM features, offering a powerful tool for ESM intercomparison and evaluation.

Error in ERA5 2m temperature identified using GraphCast

Quarterly Journal of the Royal Meteorological Society Wiley (2026) e70309

Authors:

Hannah Christensen, Jack Barker, Bobby Antonio, Massimo Bonavita, Mohamed Dahoui, Patricia de Rosnay

Abstract:

Reanalyses such as ERA5 have long been foundational for weather and climate science. They have also found a new use case, as training and verification data for machine-learnt weather prediction (MLWP) models. Here we compare short-lead time (6h) forecasts from the MLWP model GraphCast against ERA5. In doing so, we identify a recurrent, spatially coherent error in 2m Temperature centred on the Ethiopian Highlands, that occurs predominantly at 0600 UTC. We find the same error feature is present in other MLWP models trained on ERA5. We show that these error events arise from errors in ERA5 that are also present in the ECMWF operational analysis. They arise from the 2D optimal interpolation procedure, when surface reports are assimilated that are temporally displaced compared to the background forecast. This produces spuriously warm analysis increments over Ethiopia on approximately 7% of dates at 0600 UTC across the reanalysis record. The spread from the ensemble of data assimilation partially flags these cases but is underdispersive. We assess the impact on a MLWP system trained on ERA5. While the MLWP model can largely ignore these unphysical error events, a small systematic degradation in forecast skill over the region is observed. We discuss implications for using reanalysis as truth in machine learning training and verification, and recommend simple changes to reduce such artefacts in future analyses.

Local temperatures rapidly depart from historical experience beyond 1.5°C

Under review (2026)

Authors:

Simon L. L. Michel, Hannah M. Christensen

Abstract:

Local temperatures increasingly exceed values without precedent in the historical record, but when these excursions become permanent is uncertain due to the underestimated role of internal climate variability. Here we recast the problem in probabilistic terms by asking, for a given global warming level, how likely is it that annual mean local temperature will remain above its 1850-1999 maximum through 2100, given the preceding climate history. We train a neural network on 408 CMIP6 simulations under SSP2-4.5, the emissions trajectory closest to current policy commitments, and obtain calibrated probabilities skilful under out-of-sample validation. Applied to observations, regions where temperatures remain above their historical maximum are confined to just 2.2% of global land area at 1.5°C, but expand rapidly thereafter, reaching nearly half the present-day global population and 35% of aboveground biomass by 2°C, and most of humanity and nearly 80% of aboveground biomass by 3°C. Current policy trajectories will therefore commit a rapidly growing share of humanity and Earth’s ecosystems to twenty-first century local temperatures sustained above their historical maxima.

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