A Full-Coverage Daily Average PM2.5 Retrieval Method with Two-Stage IVW Fused MODIS C6 AOD and Two-Stage GAM Model

Current PM2.5 retrieval maps have many missing values, which seriously hinders their performance in real applications. This paper presents a framework to map full-coverage daily average PM2.5 concentrations from MODIS C6 aerosol optical depth (AOD) products and fill missing pixels in both the AOD and PM2.5 maps. First, a two-stage inversed variance weights (IVW) algorithm was adopted to fuse the MODIS C6 Terra and Aqua AOD products, which fills missing data in MODIS standard AOD data and obtains a high coverage daily average. After that, using the fused MODIS daily average AOD and ground-level PM2.5 in all grid cells, a two-stage generalized additive model (GAM) was implemented to obtain the full-coverage PM2.5 concentrations. Experiments on the Yangtze River Delta (YRD) in 2013−2016 were carefully designed to validate the performance of our proposed framework. The results show that the two-stage IVW could not only improve the spatial coverage of MODIS AOD against the original standard product by 230%, but could also keep its data accuracy. When compared with the ground-level measurements, the two-stage GAM can obtain accurate PM2.5 concentration estimates (R2 = 0.78, RMSE = 19.177 μg/m3, and RPE = 28.9%). Moreover, our method performs better than the inverse distance weighted method and kriging methods in mapping full-coverage daily PM2.5 concentrations. Therefore, the proposed framework provides a good methodology for retrieving full-coverage daily average PM2.5 concentrations from MODIS standard AOD products.