AirMPA: A Meteorology-to-Pollution Adapter for Global Air Quality Forecasting
Organizations: Fudan University · Shanghai Academy of AI for Science · University of Shanghai for Science and Technology
Abstract
Accurate real-time forecasting of atmospheric pollutants is essential for reducing exposure risks and improving public health. Traditional numerical forecasting systems suffer from incomplete parameterizations, uncertain inputs, and high computational cost. Existing deep learning approaches improve efficiency, but often entangle meteorological and pollutant variables within unified forecasting frameworks, limiting flexibility and pollutant-specific representation learning. To address these challenges, we propose AirMPA, a decoupled autoregressive framework for meteorology-to-pollution forecasting. AirMPA predicts 13 atmospheric pollutants at 0.4° spatial resolution and 12-h temporal resolution using meteorological fields, emission priors, and historical pollutant concentrations. The model is built on an improved Swin-Transformer architecture with zonal release attention and dual-branch fusion, enabling more effective modeling of global transport continuity and heterogeneous pollutant dynamics. During training, AirMPA uses the fifth-generation atmospheric reanalysis from the European Centre for Medium-Range Weather Forecasts (ERA5) as meteorological forcing, while during inference it directly ingests forecast fields from external meteorological models, enabling stable autoregressive prediction up to 5 days ahead. Experiments show that AirMPA consistently outperforms Aurora and surpasses the Copernicus Atmosphere Monitoring Service (CAMS) in approximately 90% of pollutant-lead-time evaluation cases over the 5-day forecast horizon, while outperforming CAMS on all variables beyond 48 h. These results demonstrate the potential of decoupled meteorology-to-pollution forecasting for global air quality prediction.
Figures & tables
| Metric | 1d | 3d | 5d | Overall |
| Better cases | 11/13 | 11/13 | 10/13 | 32/39 |
| Mean RMSE |
| Input | 1d | 3d | 5d | Overall |
| ERA5 | ||||
| FuXi | ||||
| HRES |
| Model | PM 2.5 | O 3 | TC NO 2 | Average RMSE | ||||
| 1d | 5d | 1d | 5d | 1d | 5d | 1d | 5d | |
| Swin Transformer (Baseline) | 4.31 | 6.08 | 2.595 | 4.435 | 0.651 | 0.870 | 17.68 | 32.20 |
| w/ Emission Input | 4.29 | 6.03 | 2.590 | 4.630 | 0.648 | 0.872 | 16.84 | 30.08 |
| w/ Cube Embedding ∗ | 4.25 | 6.13 | 2.558 | 4.425 | 0.619 | 0.845 | 16.76 | 30.00 |
| w/ Zonal Release Attention ∗ | 4.28 | 6.03 | 2.543 | 4.420 | 0.627 | 0.864 | 16.80 | 30.05 |
| w/ Dual-branch Fusion ∗ | 4.12 | 5.79 | 2.540 | 4.361 | 0.622 | 0.826 | 16.28 | 28.95 |
Appendix figures & tables14 assets
Supplementary material from the paper’s appendix.
Appendix
| Data Type | Source | Variables | Spatial Res. | Temporal Res. | Period |
| Pollutants | CAMS-RA | pm1, pm2p5, pm10, tcco, tcno2, tc_no, gtco3, tcso2, co_1000, no_1000, no2_1000, so2_1000, go3_1000 | 12 h | 2010.01–2021.12 | |
| CAMS-FC | Same as above | 12 h | 2017.10–2022.05 | ||
| Meteorology | ERA5 | U10m, v10m, d2m, t2m, msl, tcc, tp; z, u, v, t (at 13 levels: 50, 100, 150, 200, 250, 300, 400, 500, 600, 700, 850, 925, 1000 hPa) | 12 h | 2010.01–2023.12 | |
| FuXi | U10m, v10m, d2m, t2m, msl, tcc, tp; z, u, v, t (at 13 levels) | 12 h | 2022.06–2022.11 | ||
| Emissions | CAMS emission | ammonia, black carbon, carbon monoxide, nitrogen oxides, NMVOCs, sulfur dioxide | 12 h | 2010.01–2023.12 |
| Stage / Split | Dataset | Purpose | Period |
| Stage 1 | CAMS-RA | Global pre-training | 2010.01–2021.12 |
| Stage 2 | CAMS-FC | Operational fine-tuning | 2017.10–2022.05 |
| Test | CAMS-FC | Final evaluation | 2022.06–2022.11 |
| Stage | Objective | Learning rate | Loss settings |
| Stage 1 | Global pre-training on CAMS reanalysis | Hybrid loss: , , | |
| Stage 2 | Operational fine-tuning on CAMS forecast data | Hybrid loss: , , | |
| Stage 3 | Extreme-event optimization | continued from Stage 2 | Weighted RMSE loss: overall RMSE , extreme-event RMSE |
| Variable | PM 1 | PM 2.5 | PM 10 | TC CO | TC NO | TC NO 2 | TC SO 2 | TC O 3 |
| Training extreme threshold | 8.0 | 25.0 | 50.0 | 900 | 1.5 | 3.5 | 3.0 | 7500 |
| Threshold (90%) | 5.2323 | 20.8496 | 42.1983 | 811.4433 | 1.4372 | 3.2538 | 1.7426 | 7227.6772 |
| Variable (at 1000 hpa) | CO | NO | NO 2 | SO 2 | O 3 | |||
| Training extreme threshold | 135 | 0.05 | 0.8 | 0.5 | 45 | |||
| Threshold (90%) | 110.1245 | 0.0265 | 0.4042 | 0.2312 | 38.8720 | |||
| Model | Zonal | Cube | Dual-Branch | Emission |
| Release Attention | Embedding | Fusion | Input | |
| Swin-Transformer (backbone) | ||||
| w/ Emission | ✓ | |||
| w/ Zonal Release Attention | ✓ | ✓ | ||
| w/ Cube Embedding | ✓ | ✓ | ||
| w/ Dual-Branch Fusion | ✓ | ✓ |
| Variable | Day 1 | Day 3 | Day 5 |
| PM 1 | 19.36 | 13.88 | 11.79 |
| PM 2.5 | 21.89 | 15.17 | 8.10 |
| PM 10 | 33.46 | 30.39 | 24.90 |
| TC CO | 15.38 | 13.58 | 18.14 |
| TC NO 2 | 20.18 | 13.84 | 8.10 |
| TC NO | 29.33 | 24.76 | 18.13 |
| Day 1 | Day 3 | Day 5 | ||||
| Variable | ERA5 | HRES | ERA5 | HRES | ERA5 | HRES |
| PM 1 | 1.53 | -0.08 | 3.79 | 0.17 | 6.94 | 0.38 |
| PM 2.5 | 2.46 | 0.13 | 8.29 | 0.48 | 16.12 | 0.82 |
| PM 10 | 2.99 | 0.32 | 10.44 | 0.48 | 19.66 | 1.42 |
| TC CO | 2.15 | -0.25 | 9.22 | -0.53 | 17.06 | -1.57 |
| TC NO 2 | 3.40 | 0.16 | 11.52 | 1.90 | 19.92 | 3.56 |
Air Quality Prediction
Air quality prediction estimates future pollutant levels and maps local exposure to support health warnings, planning, and pollution control.
See all 42 papers