First retrieval of daily 160 m aerosol optical depth over urban areas using Gaofen-1/6 synergistic observations: Algorithm development and validation

The satellite-based aerosol optical depth (AOD), which can provide continuous spatial observations of aerosol loadings, is widely adopted to estimate atmospheric environmental quality and evaluate its risk for human health. However, current satellite-retrieved AOD products characterized by a comparatively coarse spatial resolution (≥1 km) can hardly analyze the structure of atmospheric pollution or its correlation with urban landscapes over populated urban areas. Existed studies have tried to address this deficiency by retrieving high-resolution AOD using Landsat images, whose long revisit period (16 days nominally), however, largely limits its applications related to urban air pollution research. To achieve both high spatial resolution and expected revisit period from satellite observation, in this study, a comprehensive AOD retrieval framework is developed for Wide-Field-of-View (WFV) satellite sensors to yield a daily AOD dataset over urban areas with 160 m spatial resolution, based on the Gaofen-1 and Gaofen-6 synergistic observations. To address the crucial challenge that the high spatial resolution and complex urban landscape both contribute a dramatic variation of land surface bidirectional reflectance, a Simulated-Annealing-coupled Semiempirical Modified Rahman-Pinty-Verstraete (SAS-MRPV) scheme is proposed to model and estimate the land surface reflectance (LSR) in the AOD retrieval framework, where the SAS-MRPV scheme is implemented using simulated annealing iteration initialized by the bidirectional reflectance distribution function (BRDF) products from the Moderate Resolution Imaging Spectrometer (MODIS) as a priori knowledge. The validation results demonstrate that the Gaofen WFV AOD retrievals exhibit good agreement with the ground-based AErosol RObotic NETwork (AERONET) and SONET (Sun-sky radiometer Observation NETwork) AOD, demonstrating the correlation coefficient and the expected error (EE)±(0.05 + 0.15AODAERONET/SONET) ratio respectively reaching up to 0.97 and over 80 %, and only slight bias fluctuations are found under different land cover types, aerosol loading, and seasonal conditions. Moreover, the validation results of Gaofen WFV AOD retrievals, with LSR estimated based on the SAS-MRPV scheme, the MODIS BRDF Products, and the Minimum Reflectance Technique (MRT) method, demonstrate better accuracy and completeness of AOD retrievals using the SAS-MPRV scheme over both natural and impervious land surfaces, indicating the stability and reliability of proposed SAS-MRPV LSR determination scheme in WFV AOD retrieval. In addition to that, the error analyses and quality control results demonstrate that the SAS-MRPV scheme is effective and necessary to guarantee the reliability and accuracy of land surface bidirectional reflectance estimation by identifying and excluding the land surface pixels unsuitable for LSR modeling, which mitigates the uncertainty and yield better accuracy in the final WFV AOD products. Furthermore, the inter-comparison between WFV-derived AOD against the operational MODIS Deep Blue (10 km), Dark Target (3 km), and Multi-Angle Implementation of Atmospheric Correction (MAIAC, 1 km) AOD products demonstrates that the 160 m WFV AOD retrievals, on the basis of obtaining similar aerosol overall descriptions as MODIS AOD retrievals, possess the capability to characterize the spatial distribution of atmospheric pollution at a finer scale with smoother variation under both clean and polluted conditions. With outstanding accuracy and reliable performance, the developed comprehensive high-resolution AOD retrieval framework exhibits substantial potential for operational AOD retrieval of Gaofen WFV satellite observations, which could be directly applied to other urban areas supporting further studies related to air pollution emission management and health risk assessment at extra-fine spatial scale.