Surface reflectance (SR) estimation is the most critical preprocessing step for deriving geophysical parameters in multi-sensor remote sensing. Most state-of-the-art SR estimation methods, such as the vector version of the Second Simulation of the Satellite Signal in the Solar Spectrum (6SV) radiative transfer (RT) model, depend on accurate information on aerosol and atmospheric gases. In this study, a Simplified and Robust Surface Reflectance Estimation Method (SREM) based on the equations from 6SV RT model, without integrating information of aerosol particles and atmospheric gasses, is proposed and tested using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper plus (ETM+), and Landsat 8 Operational Land Imager (OLI) data from 2000 to 2018. For evaluation purposes, (i) the SREM SR retrievals are validated against in situ SR measurements collected by Analytical Spectral Devices (ASD) from the South Dakota State University (SDSU) site, USA; (ii) cross-comparison between the SREM and Landsat spectral SR products, i.e., Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Landsat 8 Surface Reflectance Code (LaSRC), are conducted over 11 urban (2013–2018), 13 vegetated (2013–2018), and 11 desert/arid (2000 to 2018) sites located over different climatic zones at a global scale; (iii) the performance of the SREM spectral SR retrievals for low to high aerosol loadings is evaluated; (iv) spatio-temporal cross-comparison is conducted for six Landsat paths/rows located in Asia, Africa, Europe, and the United States of America from 2013 to 2018 to consider a large variety of land surfaces and atmospheric conditions; (v) cross-comparison is also performed for the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), and the Soil Adjusted Vegetation Index (SAVI) calculated from both the SREM and Landsat SR data; (vi) the SREM is also applied to the Sentinel-2A and Moderate Resolution Imaging Spectrometer (MODIS) data to explore its applicability; and (vii) errors in the SR retrievals are reported using the mean bias error (MBE), root mean squared deviation (RMSD), and mean systematic error (MSE). Results depict significant and strong positive Pearson’s correlation (r), small MBE, RMSD, and MSE for each spectral band against in situ ASD data and Landsat (LEDAPS and LaSRC) SR products. Consistency in SREM performance against Sentinel-2A (r = 0.994, MBE = −0.009, and RMSD = 0.014) and MODIS (r = 0.925, MBE = 0.007, and RMSD = 0.014) data suggests that SREM can be applied to other multispectral satellites data. Overall, the findings demonstrate the potential and promise of SREM for use over diverse surfaces and under varying atmospheric conditions using multi-sensor data on a global scale.

In this study, Terra-MODIS (Moderate Resolution Imaging Spectroradiometer) Collections 6 and 6.1 (C6 & C6.1) aerosol optical depth (AOD) retrievals with the recommended high-quality flag (QF = 3) were retrieved from Dark-Target (DT), Deep-Blue (DB) and merged DT and DB (DTB) level–2 AOD products for verification against Aerosol Robotic Network (AERONET) Version 3 Level 2.0 AOD data obtained from 2004–2014 for three sites located in the Beijing-Tianjin-Hebei (BTH) region. These are: Beijing, located over mixed bright urban surfaces, XiangHe located over suburban surfaces, and Xinglong located over hilly and vegetated surfaces. The AOD retrievals were also validated over different land-cover types defined by static monthly NDVI (Normalized Difference Vegetation Index) values obtained from the Terra-MODIS level-3 product (MOD13A3). These include non-vegetated surfaces (NVS, NDVI < 0.2), partially vegetated surfaces (PVS, 0.2 ≤ NDVI ≤ 0.3), moderately vegetated surfaces (MVS, 0.3 < NDVI < 0.5) and densely vegetated surfaces (DVS, NDVI ≥ 0.5). Results show that the DT, DB, and DTB-collocated retrievals achieve a high correlation coefficient of ~ 0.90–0.97, 0.89–0.95, and 0.86–0.95, respectively, with AERONET AOD. The DT C6 and C6.1 collocated retrievals were comparable at XiangHe and Xinglong, whereas at Beijing, the percentage of collocated retrievals within the expected error (↔EE) increased from 21.4% to 35.5%, the root mean square error (RMSE) decreased from 0.37 to 0.24, and the relative percent mean error (RPME) decreased from 49% to 27%. These results suggest significant relative improvement in the DT C6.1 product. The percentage of DB-collocated AOD retrievals ↔EE was greater than 70% at Beijing and Xinglong, whereas less than 66% was observed at XiangHe. Similar to DT AOD, DTB AOD retrievals performed well at XiangHe and Xinglong compared with Beijing. Regionally, DB C6 and C6.1-collocated retrievals performed better than DT and DTB in terms of good quality retrievals and relatively small errors. For diverse vegetated surfaces, DT-collocated retrievals reported small errors and good quality retrievals only for NVS and DVS, whereas larger errors were reported for PVS. MVS. DB contains good quality AOD retrievals over PVS, MVS, and DVS compared with NVS. DTB C6.1 collocated retrievals were better than C6 over NVS, PVS, and DVS. C6.1 is substantially improved overall, compared with C6 at local and regional scales, and over diverse vegetated surfaces.

In urban areas, fine particulate matter (PM_2.5) associated with local vehicle emissions can cause respiratory and cardiorespiratory disease and increased mortality rates, but less so in rural areas. However, Hong Kong may be a special case, since the whole territory often suffers from regional haze from nearby mainland China, as well as local sources. Therefore, to understand which areas of Hong Kong may be affected by damaging levels of fine particulates, PM_2.5 data were obtained from March 2005 to February 2009 for urban, suburban, and rural air quality monitoring stations; namely Central (city area, commercial area, and urban populated area), Tsuen Wan (city area, commercial area, urban populated, and residential area), Tung Chung (suburban and residential area), Yuen Long (urban and residential area), and Tap Mun (remote rural area). To evaluate the relative contributions of regional and local pollution sources, the study aimed to test the influence of weather conditions on PM_2.5 concentrations. Thus, meteorological parameters including temperature, relative humidity, wind speed, and wind directions were obtained from the Hong Kong Observatory. The results showed that Hong Kong’s air quality is mainly affected by regional aerosol emissions, either transported from the land or ocean, as similar patterns of variations in PM_2.5 concentrations were observed over urban, suburban, and rural areas of Hong Kong. Only slightly higher PM_2.5 concentrations were observed over urban sites, such as Central, compared to suburban and rural sites, which could be attributed to local automobile emissions. Results showed that meteorological parameters have the potential to explain 80% of the variability in daily mean PM_2.5 concentrations—at Yuen Long, 77% at Tung Chung, 72% at Central, 71% at Tsuen Wan, and 67% at Tap Mun, during the spring to summer part of the year. The results provide not only a better understanding of the impact of regional long-distance transport of air pollutants on Hong Kong’s air quality but also a reference for future regional-scale collaboration on air quality management.

The SEMARA approach, an integration of the Simplified and Robust Surface Reflectance Estimation (SREM) and Simplified Aerosol Retrieval Algorithm (SARA) methods, was used to retrieve aerosol optical depth (AOD) at 550 nm from a Landsat 8 Operational Land Imager (OLI) at 30 m spatial resolution, a Terra-Moderate Resolution Imaging Spectroradiometer (MODIS) at 500 m resolution, and a Visible Infrared Imaging Radiometer Suite (VIIRS) at 750 m resolution over bright urban surfaces in Beijing. The SEMARA approach coupled (1) the SREM method that is used to estimate the surface reflectance, which does not require information about water vapor, ozone, and aerosol, and (2) the SARA algorithm, which uses the surface reflectance estimated by SREM and AOD measurements obtained from the Aerosol Robotic NETwork (AERONET) site (or other high-quality AOD) as the input to estimate AOD without prior information on the aerosol optical and microphysical properties usually obtained from a look-up table constructed from long-term AERONET data. In the present study, AOD measurements were obtained from the Beijing AERONET site. The SEMARA AOD retrievals were validated against AOD measurements obtained from two other AERONET sites located at urban locations in Beijing, i.e., Beijing_RADI and Beijing_CAMS, over bright surfaces. The accuracy and uncertainties/errors in the AOD retrievals were assessed using Pearson’s correlation coefficient (r), root mean squared error (RMSE), relative mean bias (RMB), and expected error (EE = ± 0.05 ± 20%). EE is the envelope encompassing both absolute and relative errors and contains 68% (±1σ) of the good quality retrievals based on global validation. Here, the EE of the MODIS Dark Target algorithm at 3 km resolution is used to report the good quality SEMARA AOD retrievals. The validation results show that AOD from SEMARA correlates well with AERONET AOD measurements with high correlation coefficients (r) of 0.988, 0.980, and 0.981; small RMSE of 0.08, 0.09, and 0.08; and small RMB of 4.33%, 1.28%, and −0.54%. High percentages of retrievals, i.e., 85.71%, 91.53%, and 90.16%, were within the EE for Landsat 8 OLI, MODIS, and VIIRS, respectively. The results suggest that the SEMARA approach is capable of retrieving AOD over urban areas with high accuracy and small errors using high to medium spatial resolution satellite remote sensing data. This approach can be used for aerosol monitoring over bright urban surfaces such as in Beijing, which is frequently affected by severe dust storms and haze pollution, to evaluate their effects on public health.

Numerous studies (hereafter GA: general approach studies) have been made to classify aerosols into desert dust (DD), biomass-burning (BB), clean continental (CC), and clean maritime (CM) types using only aerosol optical depth (AOD) and Ångström exponent (AE). However, AOD represents the amount of aerosol suspended in the atmospheric column while the AE is a qualitative indicator of the size distribution of the aerosol estimated using AOD measurements at different wavelengths. Therefore, these two parameters do not provide sufficient information to unambiguously classify aerosols into these four types. Evaluation of the performance of GA classification applied to AErosol Robotic NETwork (AERONET) data, at sites for situations with known aerosol types, provides many examples where the GA method does not provide correct results. For example, a thin layer of haze was classified as BB and DD outside the crop burning and dusty seasons respectively, a thick layer of haze was classified as BB, and aerosols from known crop residue burning events were classified as DD, CC, and CM by the GA method. The results also show that the classification varies with the season, for example, the same range of AOD and AE were observed during a dust event in the spring (20^th March 2012) and a smog event in the autumn (2nd November 2017). The results suggest that only AOD and AE cannot precisely classify the exact nature (i.e., DD, BB, CC, and CM) of aerosol types without incorporating more optical and physical properties. An alternative approach, AEROsol generic classification using a novel Satellite remote sensing Approach (AEROSA), is proposed to provide aerosol amount and size information using AOD and AE, respectively, from the Terra-MODIS (MODerate resolution Imaging Spectroradiometer) Collection 6.1 Level 2 combined Dark Target and Deep Blue (DTB) product and AERONET Version 3 Level 2.0 data. Although AEROSA is also based on AOD and AE, it does not claim the nature of aerosol types, instead providing information on aerosol amount and size. The purpose is to introduce AEROSA for those researchers who are interested in the generic classification of aerosols based on AOD and AE, without claiming the exact aerosol types such as DD, BB, CC, and CM. AEROSA not only provides 9 generic aerosol classes for all observations but can also accommodate variations in location and season, which GA aerosol types do not.

This letter reports uncertainties in the Aqua-Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2 dark target (DT), deep blue (DB), and multiangle implementation of atmospheric correction (MAIAC) aerosol optical depth (AOD) during the COVID-19 lockdown period (February–May 2020) compared to the pre-COVID-19 period (February–May 2019). Validation of AOD retrievals was conducted against AErosol RObotic NETwork (AERONET) Version 3 Level 1.5 AOD data obtained from three sites located in urban (Beijing_CAMS and Beijing_RADI) and suburban (XiangHe) areas of China. The results show the poor performance of the DT and DB algorithms compared to the MAIAC algorithm, which performed better during the lockdown period. Overall, all MODIS algorithms overestimated the AOD and showed higher positive bias under high aerosol loading conditions during lockdown than during prelockdown. This is mainly attributed to the overestimation of the aerosol single-scattering albedo (SSA), which was found higher during lockdown than during the same period in 2019.