Generative Atmospheric Super-Resolution from Heterogeneous In Situ Observations through Composable Interfaces
Organizations: School of Mechanical Engineering, Purdue University, West Lafayette, Indiana, USA · College of Information Sciences and Technology, The Pennsylvania State University, University Park, Pennsylvania, USA
Abstract
Atmospheric observations are sparse, heterogeneous, and unevenly distributed, whereas many generative atmospheric models learn distributions over regularly gridded multivariate states. Once pretrained, diffusion models can supply atmospheric priors that can be combined with observation-derived likelihood factors in a Bayesian formulation. However, these observation sources differ substantially in geometry and sampling density, complicating the consistent use of their observations within a common inference framework. Here, we formulate this reconstruction problem as generative atmospheric super-resolution and introduce composable observation interfaces for conditioning a single pretrained 13-variable atmospheric diffusion model. The interfaces convert sparse radiosonde (R), clustered aircraft (A), and dense irregular surface-station (S) observations into source-specific likelihood factors that specify where observations constrain the gridded state, how residuals are counted under uneven sampling, and how strongly each source guides posterior sampling. We developed the aircraft and surface observation interfaces using 2019 observations and evaluated the selected interfaces throughout 2020 without further tuning. Compared with reconstructions conditioned only on radiosonde observations, the composed R+A+S interface reduces RMSE evaluated against ERA5 by across all 13 state variables over the CONUS domain. The aircraft and surface factors provide complementary improvements in upper-air and surface variables. The R+A+S combination also lowers the Continuous Ranked Probability Score (CRPS), while evaluations at held-out aircraft and surface-station observations show reduced prediction errors. Together, these results demonstrate a modular route for conditioning a pretrained atmospheric generative prior on heterogeneous in situ observations without retraining the underlying model.
Figures & tables
| Symbol | Meaning |
|---|---|
| , | Observation-source and state-variable index sets. |
| , | Clean standardized state and its noise-level- counterpart. |
| , | Current noisy sampler state and denoised clean-state estimate at reverse step . |
| , | Generic diffusion-noise level and its value at reverse step . |
| , , | Final likelihood targets, observation operator, and residual-counting weights for source . |
| , | Likelihood weight and noise-dependent variance term for source . |
| Source | Product and sampling geometry | Variables | Spatial domain |
|---|---|---|---|
| Radiosonde (R) | IGRA sparse vertical profiles | , , , , , , , , , , , , | Global |
| Aircraft (A) | MADIS ABO reports clustered along flight routes | , , , , , | CONUS domain |
| Surface station (S) | MADIS METAR surface stations | , , | CONUS domain |
| Source | Retained observations | Constrained variables | State mapping and residual counting | Likelihood parameters |
|---|---|---|---|---|
| Radiosonde (R) | IGRA profiles | , , , , , , , , , , , , | bilinear interpolation; equal target and channel weights | |
| Aircraft (A) | MADIS ABO reports from ACARS direct, MDCRS/ARINC, and Canadian AMDAR; -hPa windows | , , , , , | cell means; equal cell and channel weights | |
| Surface station (S) | MADIS METAR reports | , , | cell means; equal cell and channel weights |
| Conditioning | All 13 | Surface-targeted | Upper-air temperature and wind |
|---|---|---|---|
| CONUS domain | |||
| R+A | |||
| R+S | |||
| R+A+S | |||
| Global | |||
| R+A | |||
| Region | Variable group | (%) | 95% interval (%) |
|---|---|---|---|
| CONUS domain | All 13 variables | ||
| Surface-targeted | |||
| Upper-air temperature and wind | |||
| Global | All 13 variables | ||
| Surface-targeted | |||
| Upper-air temperature and wind |
| Evaluated source | Conditioning comparison | Mean change | 95% interval | Lower-RMSE cases |
|---|---|---|---|---|
| Surface station (S) | Retained 80% of S added to R+A | |||
| Aircraft (A) | Retained 80% of A added to R+S |
| Variable group | RMSE change | CRPS change [95% interval] |
|---|---|---|
| All 13 variables | ||
| Surface-targeted variables | ||
| Upper-air temperature and wind variables |
Appendix figures & tables12 assets
Supplementary material from the paper’s appendix.
Appendix
| Development stage | Candidates compared | Selection criterion | Selected result |
| A interface design | 24 combinations of reporting systems, pressure-matching windows, and residual representations | Mean RMSE changes across all 13 variables and across the upper-air temperature and wind variables | ACARS direct, MDCRS/ARINC, and Canadian AMDAR; -hPa windows; equal-cell mean residuals |
| S interface design | 3 residual representations | Mean RMSE changes across all 13 variables and across the surface-targeted variables | Equal-cell mean residuals |
| A likelihood parameters | Staged parameter comparisons using R+A | RMSE across the upper-air temperature and wind variables and across all 13 variables | |
| S likelihood parameters | Staged parameter comparisons using R+S | Surface-targeted and all-variable RMSE | |
| R+A+S selection | 9 combinations of candidate A and S likelihood-parameter settings | Mean RMSE change across all 13 variables as the primary score, with the surface-targeted group and the upper-air temperature and wind group as secondary checks | The selected A and S settings above are combined, producing an mean RMSE reduction across all 13 variables in 2019 |
| Variable | Unit | R-only RMSE | R+A+S RMSE |
|---|---|---|---|
| 1.725 | 1.380 | ||
| 1.605 | 1.414 | ||
| 1.697 | 1.507 | ||
| 112.4 | 101.8 | ||
| 84.5 | 76.1 | ||
| 2.912 | 2.618 |
| Variable | (%) | 95% interval (%) | Variable | (%) | 95% interval (%) |
|---|---|---|---|---|---|
| Variable | R-only | R+A+S | Variable | R-only | R+A+S |
|---|---|---|---|---|---|
| 0.721 | 0.793 | 0.752 | 0.768 | ||
| 0.667 | 0.664 | 0.701 | 0.709 | ||
| 0.642 | 0.631 | 0.704 | 0.724 | ||
| 0.577 | 0.639 | 0.677 | 0.678 | ||
| 0.613 | 0.680 | 0.856 | 0.858 | ||
| 0.679 | 0.682 | 0.734 | 0.777 |