Geoscientific Model Development (Oct 2023)

A Regional multi-Air Pollutant Assimilation System (RAPAS v1.0) for emission estimates: system development and application

  • S. Feng,
  • F. Jiang,
  • F. Jiang,
  • Z. Wu,
  • H. Wang,
  • H. Wang,
  • W. He,
  • Y. Shen,
  • L. Zhang,
  • Y. Zheng,
  • C. Lou,
  • Z. Jiang,
  • W. Ju,
  • W. Ju

DOI
https://doi.org/10.5194/gmd-16-5949-2023
Journal volume & issue
Vol. 16
pp. 5949 – 5977

Abstract

Read online

Top-down atmospheric inversion infers surface–atmosphere fluxes from spatially distributed observations of atmospheric composition in order to quantify anthropogenic and natural emissions. In this study, we developed a Regional multi-Air Pollutant Assimilation System (RAPAS v1.0) based on the Weather Research and Forecasting–Community Multiscale Air Quality (WRF–CMAQ) modeling system model, the three-dimensional variational (3D-Var) algorithm, and the ensemble square root filter (EnSRF) algorithm. This system can simultaneously assimilate hourly in situ CO, SO2, NO2, PM2.5, and PM10 observations to infer gridded emissions of CO, SO2, NOx, primary PM2.5 (PPM2.5), and coarse PM10 (PMC) on a regional scale. In each data assimilation window, we use a “two-step” scheme, in which the emissions are inferred first and then input into the CMAQ model to simulate initial conditions (ICs) of the next window. The posterior emissions are then transferred to the next window as prior emissions, and the original emission inventory is only used in the first window. Additionally, a “super-observation” approach is implemented to decrease the computational costs, observation error correlations, and influence of representative errors. Using this system, we estimated the emissions of CO, SO2, NOx, PPM2.5, and PMC in December and July 2016 over China using nationwide surface observations. The results show that compared to the prior emissions (2016 Multi-resolution Emission Inventory for China – MEIC 2016)), the posterior emissions of CO, SO2, NOx, PPM2.5, and PMC in December 2016 increased by 129 %, 20 %, 5 %, 95 %, and 1045 %, respectively, and the emission uncertainties decreased by 44 %, 45 %, 34 %, 52 %, and 56 %, respectively. With the inverted emissions, the RMSE of simulated concentrations decreased by 40 %–56 %. Sensitivity tests were conducted with different prior emissions, prior uncertainties, and observation errors. The results showed that the two-step scheme employed in RAPAS is robust in estimating emissions using nationwide surface observations over China. This study offers a useful tool for accurately quantifying multi-species anthropogenic emissions at large scales and in near-real time.