J.-T. Lin1 , M.-Y. Liu1 , J.-Y. Xin2 , K. F. Boersma3,4 , R. Spurr5 , R. Martin6,7 , and Q. Zhang8
1 Laboratory for Climate and Ocean-Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871, China
2 State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, 100029, PR China
3 Royal Netherlands Meteorological Institute, De Bilt, the Netherlands
4 Meteorology and Air Quality department, Wageningen University, Wageningen, the Netherlands
5 RT Solutions Inc., Cambridge, MA 02138, USA
6 Department of Physics and Atmospheric Science, Dalhousie University, Halifax, Nova Scotia, Canada
7 Atomic and Molecular Physics Division, Harvard – Smithsonian Center for Astrophysics, Cambridge, Massachusetts, USA
8 Ministry of Education Key Laboratory for Earth System Modeling, Center for Earth System Science, Tsinghua University, Beijing, China
Abstract
Satellite retrievals of vertical column densities (VCDs) of tropospheric nitrogen dioxide ( ) normally do not explicitly account for aerosol optical effects and surface reflectance anisotropy that vary with space and time. Here, we conduct an improved retrieval of VCDs over China, called the POMINO algorithm, based on measurements from the Ozone Monitoring Instrument (OMI), and we test the importance of a number of aerosol and surface reflectance treatments in this algorithm. POMINO uses a parallelized LIDORT-driven AMFv6 package to derive tropospheric air mass factors via pixel-specific radiative transfer calculations with no look-up tables, taking slant column densities from DOMINO v2. Prerequisite cloud optical properties are derived from a dedicated cloud retrieval process that is fully consistent with the main retrieval. Aerosol optical properties are taken from GEOS-Chem simulations constrained by MODIS aerosol optical depth (AOD) data. MODIS bidirectional reflectance distribution function (BRDF) data are used for surface reflectance over land. For the present analysis, POMINO level-2 data for 2012 are aggregated into monthly means on a 0.25°long 0.25°lat. grid. POMINO-retrieved annual mean VCDs vary from 15–25 over the polluted North China Plain (NCP) to below over much of western China. Using POMINO to infer Chinese emissions of nitrogen oxides leads to annual anthropogenic emissions of 9.05 , an increase from 2006 (Lin, 2012) by about 19 %. Replacing the MODIS BRDF data with the OMLER v1 monthly climatological albedo data affects VCDs by up to 40% for certain locations and seasons. The effect on constrained NOx emissions is small. Excluding aerosol information from the retrieval process (this is the traditional “implicit” treatment) enhances annual mean VCDs by 15–40% over much of eastern China. Seasonally, VCDs are reduced by 10–20%over parts of the NCP in spring and over northern China in winter, despite the general enhancements in summer and fall. The effect on subsequently constrained annual emissions is between -5 and +30% with large seasonal and spatial dependence. The implicit aerosol treatment also tends to exclude days with high pollution, since aerosols are interpreted as effective clouds and the respective OMI pixels are often excluded by cloud screening; this is a potentially important sampling bias. Therefore an explicit treatment of aerosols is important for space-based retrievals and emission constraints. A comprehensive independent measurement network with sufficient spatial and temporal representativeness is needed to further evaluate the different satellite retrieval approaches.