physics.ao-phMar 30, 2026

Downscaling weather forecasts from Low- to High-Resolution with Diffusion Models

Authors: Joffrey Dumont Le Brazidec, Simon Lang, Martin Leutbecher, Baudouin Raoult, Gert Mertes, Florian Pinault, Aristofanis Tsiringakis, Pedro Maciel, 

Organizations: European Centre for Medium-Range Weather Forecasts (ECMWF)

Read PDF
Ask AI
Cite

Copy citation

DOINot available for this paper

Abstract

We introduce a probabilistic diffusion-based method for global atmospheric downscaling implemented within the Anemoi framework. The approach transforms low-resolution ensemble forecasts into high-resolution ensembles by learning the conditional distribution of finer-scale residuals, defined as the difference between the high-resolution fields and the interpolated low-resolution inputs. The system is trained on reforecast pairs from ECMWF IFS, using coarse fields at 100 km to reconstruct fine-scale variability at 30 km resolution. The bulk of the training focuses on recovering small-scale structures, while fine-tuning in high-noise regimes enables the generation of extremes. Evaluation against the medium-range IFS ensemble target shows that the model increases probabilistic skill (FCRPS) for surface variables, reproduces target power spectra at small scales, captures physically consistent multivariate relationships such as wind-pressure coupling, and generates extreme values consistent with those of the target ensemble in tropical cyclones.

Figures & tables

Figure 1: Hurricanes Idalia and Franklin as seen in the subseasonal ensemble forecast (eefo, O96 resolution) and its AI-based downscaling to O320 resolution. The background field in each panel shows the low-resolution O96 input, while the inset (black frame) shows the corresponding downscaled O320 prediction from the diffusion model. Top: mean sea level pressure (hPa). Bottom: 10-metre wind speed (m/s). Initialisation date: 29 August 2023, +24 h lead time (valid 30 August 2023), ensemble member 6. The downscaled fields resolve the tight inner-core structure of both tropical cyclones — including the pressure minima and eyewall wind maxima — that is absent from the coarse input.
Field / Weather variableLevelInput/TargetNormalisation
Geopotential, horizontal and vertical wind components, temperature, specific humidityPressure levels: 50, 100, 200, 300, 400, 500, 700, 850, 925, 1000 hPaInput: low-res Target: high-resStandardised, except geopotential (max-scaled)
10m wind (u, v), 2m temperature, 2m dewpoint temperature, surface pressure, skin temperature, mean sea level pressure, total column waterSurfaceInput: low-res Target: high-resStandardised
Land–sea mask, orography, cosine/sine of local time of day, day of year, latitude, longitude, insolation,SurfaceInput: high-resNot normalised
Table 1: Meteorological fields, levels, and normalisation used for training Anemoi-D 2 .
Figure 2: Schematic of the diffusion-based downscaling model. At each denoising step, the model takes as input on the O320 grid: (i) the current noisy residual estimate, (ii) static and temporal forcing fields, and (iii) the coarse O96 atmospheric fields interpolated to O320 ( D96→320D_{96\to 320} ). These features are processed by an encoder–processor–decoder graph neural network: the encoder projects information from O320 to the hidden grid, the processor applies 16 successive graph-transformer layers on the hidden graph, and the decoder projects the result back to O320. In the lower panel, the pink boxes represent these three consecutive components. The model output is a prediction of the residual on O320, which is added to the interpolated coarse field to reconstruct the final high-resolution field.
Figure 3: Receptive field of the GNN processor for a randomly selected node (red). The figure shows how many nodes receive its information after 2, 10, and 16 message-passing layers—that is, across the corresponding 2-, 10-, and 16-hop neighborhoods. The graph structure is fixed; only the range of information propagation increases.
StageStepsNoise setupObjective
110610^{6}EDM lognormal noise distribution ( Karras et al., 2022 ) with mean −1.2-1.2 and standard deviation 1.21.2Base training
210510^{5}Higher-noise fine-tuning with a Karras-style schedule, using ρ=7\rho=7 , σmax⁡=10000\sigma_{\max}=10000 , and Smax⁡=10000S_{\max}=10000Improve generation from strongly degraded states and better capture extremes
Table 2: Summary of the two-stage training procedure.
Figure 4: Fair Continuous Ranked Probability Score (FCRPS) for 2-metre temperature (left) and 10-metre wind speed (right) over the Northern Hemisphere extratropics, from forecast day 1 to 10 (20240201–20240225), evaluated against SYNOP station observations. Comparison between the downscaled ensemble forecast, the IFS medium-range ensemble forecast, and the subseasonal ensemble forecast (input).
Figure 5: Subseasonal ensemble forecast first member (first column), medium-range ensemble forecast first member (second column), and downscaling with the diffusion model of two subseasonal members (third and fourth columns, with the third corresponding to downscaling of the first column). Results are shown for the Amazon rainforest at lead time 24h. Rows display 10u, 10v, 2t, and T850. The diffusion model generates stochastic and physically consistent small-scale features.
Figure 6: Subseasonal ensemble forecast first member (first column), medium-range ensemble forecast first member (second column), and downscaling with the diffusion model of two subseasonal members (third and fourth columns, with the third corresponding to downscaling of the first column). Results are shown for the Himalayas at lead time 24h. Rows display 10u, 10v, 2t, and T850. The diffusion model generates physically consistent small-scale features.
Figure 7: Comparison of global amplitude spectra of weather variables from the subseasonal forecasts, the medium-range forecasts, and the downscaled subseasonal forecasts. The spectra are evaluated over the full globe at forecast day 6 and are based on 31 forecast initialisation dates and 10 ensemble members per date, corresponding to a sample of 310 forecasts.
Figure 8: Normalised pdfs of mean sea level pressure (980–1015 hPa) and 10-meter wind speed (0–30 m/s) for Hurricane Idalia across ensemble members, initialized on 26-30 August 2023, with a 1 to 5-days lead time. The pdfs represent the distributions of gridpoint values within the cyclone region, normalised by the corresponding distribution of the operational analysis. They are shown for the target IFS medium-range ensemble forecasts, the input IFS subseasonal ensemble forecasts, and the downscaled reconstructions (before and after finetuning). The accompanying map illustrates Idalia’s track from 27 to 31 August 2023, marking stages from tropical depression to major hurricane. Idalia cyclone track is adapted from NOAA/NHC https://www.nhc.noaa.gov/data/tcr/index.php?season=2023&basin=atl .
Figure 9: Mean sea level pressure (985–1015 hPa) and wind speed (0–30 m/s) for tropical cyclone Idalia for three ensemble members at 24-hour lead time, initialized on 28 August 2023. The first column shows the target IFS medium-range ensemble forecast MSLP, the second column displays the input IFS subseasonal ensemble forecast MSLP, and the third and fourth columns present the MSLP and wind fields, respectively, reconstructed through downscaling of the IFS subseasonal ensemble forecast.
Figure 10: Normalised pdfs of mean sea level pressure (960–1015 hPa) and 10-meter wind speed (0–30 m/s) for Hurricane Hilary across ensemble members, initialized on 16-20 August 2023, with a 1 to 5-days lead time. The pdfs represent the distributions of gridpoint values within the cyclone region (approximately the dashed orange line), normalised by the corresponding distribution of the operational analysis. They are shown for the target IFS medium-range ensemble forecasts, the input IFS subseasonal ensemble forecasts, and the downscaled reconstructions (before and after finetuning). The map illustrates the cyclone’s track from 16 to 20 August 2023, marking stages from tropical depression to major hurricane. Hilary cyclone track is adapted from NOAA/NHC https://www.nhc.noaa.gov/data/tcr/index.php?season=2023&basin=epac .
Figure 11: Local fields of mean sea level pressure (960–1015 hPa) and wind speed (0–30 m/s) for Tropical Cyclone Hilary across ensemble members at 24-hour lead time, initialized on 17 August 2023. The first column shows the target IFS medium-range ensemble forecast MSLP, the second column displays the input IFS subseasonal ensemble forecast MSLP, and the third and fourth columns present the MSLP and wind fields, respectively, reconstructed through downscaling of the IFS subseasonal ensemble forecast.
More in this topic

Weather Forecasting

Weather forecasting predicts future atmospheric conditions, from local storms to global patterns, using physical models, observations, and machine learning.

See all 175 papers

Explore similar work

CardsList
Jul 19, 2026

Apeliotes: A Diffusion-Based Modeling Framework for km-scale Multi-Level Atmospheric Fields

High-resolution atmospheric data are required to resolve mesoscale and localized meteorological structures, however such datasets remain limited in many regions of the world. Existing high-resolution weather products are typically produced through dynamical downscaling, which is computationally expensive and difficult to scale across locations, variables, and forecast scenarios. These limitations motivate machine-learning-based downscaling systems that can generate multiple weather variables stochastically while producing new high-resolution fields directly. In this paper we present Apeliotes, a framework for high-resolution weather forecasting. Built on the global re-analysis atmospheric data, a pre-trained global weather foundation model, and a regionally trained generative diffusion model, Apeliotes not only provides accurate kilometer-scale weather variables, but also multi-level atmospheric fields which are not directly available in the existing global atmospheric data. Our comprehensive evaluation demonstrates that Apeliotes achieves highly competitive performance. The model predicts vertical wind profile with less than 3% error between truth and predicted fields, achieving correlations of 0.91 for 10-m wind speed and 0.99 for 2-m temperature, with NRMSE values of 0.42 and 0.17, respectively.

Evangelia Rafaela Frastali, Achyut Paudel, Maryam Golbazi +1

Old Dominion University · Norfolk, VA 23529, USA

Diffusion ModelsWeather Forecasting
AskPDF
Sep 15, 2026

Efficient Analog-Initialized Latent Transport for Spatially Coherent Probabilistic Downscaling from Global to Kilometer Scales

Probabilistic downscaling must represent kilometer-scale structure left unresolved by a deterministic regional prediction. We introduce an analog-initialized latent transport method that uses historical regional residuals as a meteorologically informed empirical source. A frozen graph neural network predicts the deterministic state from the Integrated Forecasting System (IFS). For each input, residuals from 20 training dates with similar IFS conditions initialize the ensemble. A conditional flow-matching model transports these fields in a carefully designed 65-factor, 100,000-node latent space, and a grid decoder reconstructs 81 atmospheric channels including precipitation at 2.5km spacing. On a season- and cycle-balanced 2025 test, ensemble CRPS improves over deterministic mean absolute error for all 81 variables. Matched analogs improve point and spatial errors over Gaussian initialization and reduce spectral error compared with covariance matching. Comparisons with released CorrDiff, a local CorrDiff-mini model and full grid EDM provide context on the 2023 validation sample. Total training of the downscaling model takes approximately 12h on one NVIDIA H200, and inference takes about 12s per 81-field member. The frozen networks also generate fields over France when supplied with local static data and a high-resolution IFS-derived analog bank. The method combines empirical initialization with whole-domain latent transport for computationally practical regional ensembles.

Ophélia Miralles

Norwegian Meteorological Institute, Oslo, Norway

Probabilistic ForecastingWeather ForecastingPrecipitation DownscalingFlow Matching
AskPDF
Mar 23, 2026

IPSL-AID: Generative Diffusion Models for Climate Downscaling from Global to Regional Scales

Effective adaptation and mitigation strategies for climate change require high-resolution projections to inform strategic decision-making. Conventional global climate models, which typically operate at resolutions of 150 to 200 kilometers, lack the capacity to represent essential regional processes. IPSL-AID is a global to regional downscaling tool based on a denoising diffusion probabilistic model designed to address this limitation. Trained on ERA5 reanalysis data, it generates 0.25 degree resolution fields for temperature, wind, and precipitation using coarse inputs and their spatiotemporal context. It also models probability distributions of fine-scale features to produce plausible scenarios for uncertainty quantification. The model accurately reconstructs statistical distributions, including extreme events, power spectra, and spatial structures. This work highlights the potential of generative diffusion models for efficient climate downscaling with uncertainty

Kishanthan Kingston, Olivier Boucher, Freddy Bouchet +5

Climate Modeling Center, Sorbonne University, CNRS, IPSL, Paris, France · LMD, CNRS, ENS/PSL, Sorbonne University, IPSL, Paris, France · Climate Modeling and Analysis, LLC, Superior, CO, USA +1

Climate DownscalingPrecipitation DownscalingDiffusion Models
AskPDF