Comparison and validation of band residual difference algorithm and principal component analysis algorithm for retrievals of atmospheric SO2 columns from satellite observations

Abstract

Remote sensing technology provides an unprecedented tool for the continuous and real-time monitoring of atmospheric SO2 from volcanic eruption and anthropogenic emission. The Global Ozone Monitoring Experiment (GOME), SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY), and Ozone Monitoring Instrument (OMI) have high SO2 monitoring capability. The OMI, which was launched on the EOS/Aura platform in July 2004, has the same hyperspectral measurements as the GOME and SCIAMACHY, but offers the improved spatial resolution at nadir (1324 km2) and daily global coverage for short-lifetime SO2. For OMI operational SO2 planetary boundary layer (PBL) retrieval, the previous band residual difference (BRD) algorithm has been replaced by principal component analysis (PCA) algorithm, which effectively reduces the systematic biases in SO2 column retrievals. However, there are few studies on the evaluations and validations of PCA SO2 retrievals over China, and the long-term comparisons with BRD SO2 retrievals also need to be conducted. In this study, the accuracies of PCA and BRD SO2 retrievals are validated by using ground-based multi axis differential optical absorption spectroscopy (MAX-DOAS) located in Beijing, and regional atmospheric modeling system, community multi-scale air quality (RAMS-CMAQ) modeling system model which can simulate the vertical distribution of atmospheric SO2. Moreover, BRD and PCA SO2 retrievals from oceanic area, eastern China and Reunion volcanic eruption are compared to find the long-term trend and spatiotemporal differences between SO2 columns. Finally, the uncertainty of SO2 retrieval, caused by measurement errors, band selection and input parameter errors in radiative transfer model, are analysed to understand the limitations of BRD and PCA algorithms. Results show that both PCA and BRD SO2 retrievals over Beijing are lower than ground-based MAX-DOAS measurements of SO2. PCA and BRD SO2 retrievals over eastern China are lower than the simulated SO2 columns from RAMS-CMAQ in winter 2008, but in July and August BRD SO2 columns are higher than RAMS-CMAQ simulations. The values of SO2 columns from BRD over China are more consistent with those from ground-based MAX-DOAS and RAMS-CMAQ model than from PCA. Although PCA algorithm effectively reduces the noise in SO2 column retrieval, SO2 columns from PCA over China are lower than those from BRD. For oceanic area where SO2 amount is nearly zero, the standard deviation of PCA results is lower than that of BRD, but the absolute value of averaged PCA SO2 column is larger than that of BRD. In the case of Reunion volcanic eruption with SO2 columns larger than 25 DU, the BRD SO2 columns are lower than PCA retrievals. Meanwhile, with the increase of SO2 column, the difference between BRD and PCA SO2 retrievals increases. Detailed uncertainty analysis shows the influences of measurement errors, band selection and inputs of radiative transfer model on the retrieval results. This study is important for developing the retrieval algorithm, and can also improve the application of OMI SO2 products.