Satellite-retrieved Aerosol Optical Depth Research Articles

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.

A deep learning-based imputation method for missing gaps in satellite aerosol products by fusing numerical model data

Satellite-based aerosol optical depth (AOD) products are commonly used in various aerosol-related studies, such as aerosol pollution mapping and aerosol-climate interactions. However, these satellite AOD products often suffer from significant missing gaps due to cloud cover and limitations in the retrieval algorithm. To address this issue, some studies take advantage of real-time seamless simulation of numerical models and successfully fill in these gaps by establishing a regression relationship between satellite AOD and numerical model AOD. However, these previous studies usually use satellite AOD retrievals as the regression target, which limits the accuracy of the imputation results by the original accuracy of satellite AOD retrievals and also consumes a considerable amount of time. To overcome these limitations, this study proposes a spatiotemporal imputation model called Bi-ConvRNN, which combines convolutional neural networks (CNN) and bidirectional recurrent neural networks (Bi-RNN). The model takes both satellite AOD retrievals and numerical model AOD data as input and utilizes the weighted mean squared error (MSE) loss function of multiple AOD datasets, e.g., ground-based data, satellite retrievals, and numerical simulation, as the optimization target to improve the imputation accuracy. The proposed model is evaluated using hourly COMS GOCI AOD products. In the independent test set, the AOD results generated by the Bi-ConvRNN model in the region containing GOCI AOD retrievals can break the accuracy of original GOCI AOD products with the accuracy improved from R2 = 0.70 [RMSE = 0.15] to R2 = 0.84 [RMSE = 0.11], and the filling accuracy, e.g. R2 = 0.79, [RMSE = 0.14], in the region without GOCI AOD retrievals are still better than those of the original GOCI AOD retrievals. Additionally, the Bi-ConvRNN model demonstrates satisfactory filling efficiency, requiring only 0.12 s to fill in the missing gaps of hourly GOCI AOD products per day. These results highlight the efficiency and reliability of the proposed model in filling the gaps in satellite AOD products, and the filled AOD results have great potential for further aerosol-related research.

Estimation of PM2.5 and PM10 Mass Concentrations in Beijing Using Gaofen-1 Data at 100 m Resolution

Due to the advantage of high spatial coverage, using satellite-retrieved aerosol optical depth (AOD) data to estimate PM2.5 and PM10 mass concentrations is a current research priority. Statistical models are the common method of PM estimation currently, which do not require the knowledge of complex chemical and physical interactions. However, the statistical models rely on station data, which results in less accurate PM estimation concentrations in areas where station data are missing. Hence, a new hybrid model, with low dependency on on-site data, was proposed for PM2.5 and PM10 mass concentration estimation. The Gaofen-1 satellite and MODIS data were employed to estimate PM2.5 and PM10 concentrations with 100 m spatial resolution in Beijing, China. Then, the estimated PM2.5/10 mass concentration data in 2020 were employed to conduct a spatio-temporal analysis for the investigation of the particulate matter characteristic in Beijing. The estimation result of PM2.5 was validated by the ground stations with R2 ranging from 0.91 to 0.98 and the root mean square error (RMSE) ranging from 4.51 μg/m3 to 17.04 μg/m3, and that for PM10 was validated by the ground stations with R2 ranging from 0.85 to 0.98 and the RMSE ranging from 6.98 µg/m3 to 29.00 µg/m3. The results showed that the hybrid model has a good performance in PM2.5/10 estimation and can improve the coverage of the results without sacrificing the effectiveness of the model, providing more detailed spatial information for urban-scale studies.

Cloud condensation nuclei concentrations derived from the CAMS reanalysis

Abstract. Determining number concentrations of cloud condensation nuclei (CCN) is one of the first steps in the chain in analysis of cloud droplet formation, the direct microphysical link between aerosols and cloud droplets, and a process key for aerosol–cloud interactions (ACI). However, due to sparse coverage of in situ measurements and difficulties associated with retrievals from satellites, a global exploration of their magnitude, source as well as temporal and spatial distribution cannot be easily obtained. Thus, a better representation of CCN numbers is one of the goals for quantifying ACI processes and achieving uncertainty-reduced estimates of their associated radiative forcing. Here, we introduce a new CCN dataset which is derived based on aerosol mass mixing ratios from the latest Copernicus Atmosphere Monitoring Service reanalysis (CAMSRA) in a diagnostic model that uses CAMSRA aerosol properties and a simplified kappa-Köhler framework suitable for global models. The emitted aerosols in CAMSRA are not only based on input from emission inventories using aerosol observations, they also have a strong tie to satellite-retrieved aerosol optical depth (AOD) as this is assimilated as a constraining factor in the reanalysis. Furthermore, the reanalysis interpolates for cases of poor or missing retrievals and thus allows for a full spatiotemporal quantification of CCN numbers. The derived CCN dataset captures the general trend and spatial and temporal distribution of total CCN number concentrations and CCN from different aerosol species. A brief evaluation with ground-based in situ measurements demonstrates the improvement of the modelled CCN over the sole use of AOD as a proxy for CCN as the overall correlation coefficient improved from 0.37 to 0.71. However, we find the modelled CCN from CAMSRA to be generally high biased and find a particular erroneous overestimation at one heavily polluted site which emphasises the need for further validation. The CCN dataset (https://doi.org/10.26050/WDCC/QUAERERE_CCNCAMS_v1, Block, 2023), which is now freely available to users, features 3-D CCN number concentrations of global coverage for various supersaturations and aerosol species covering the years 2003–2021 with daily frequency. This dataset is one of its kind as it offers lots of opportunities to be used for evaluation in models and in ACI studies.

A Deep-Learning and Transfer-Learning Hybrid Aerosol Retrieval Algorithm for FY4-AGRI: Development and Verification over Asia

The Advanced Geosynchronous Radiation Imager (AGRI) is a mission-critical instrument for the Fengyun series of satellites. AGRI acquires full-disk images every 15 min and views East Asia every 5 min through 14 spectral bands, enabling the detection of highly variable aerosol optical depth (AOD). Quantitative retrieval of AOD has hitherto been challenging, especially over land. In this study, an AOD retrieval algorithm is proposed that combines deep learning and transfer learning. The algorithm uses core concepts from both the Dark Target (DT) and Deep Blue (DB) algorithms to select features for the machine-learning (ML) algorithm, allowing for AOD retrieval at 550 nm over both dark and bright surfaces. The algorithm consists of two steps: ①A baseline deep neural network (DNN) with skip connections is developed using 10 min Advanced Himawari Imager AODs as the target variable, and ②Sunphotometer AODs from 89 ground-based stations are used to fine-tune the DNN parameters. Out-of-station validation shows that the retrieved AOD attains high accuracy, characterized by a coefficient of determination (R2) of 0.70, a mean bias error (MBE) of 0.03, and a percentage of data within the expected error (EE) of 70.7%. A sensitivity study reveals that the top-of-atmosphere reflectance at 650 and 470 nm, as well as the surface reflectance at 650 nm, are the two largest sources of uncertainty impacting the retrieval. In a case study of monitoring an extreme aerosol event, the AGRI AOD is found to be able to capture the detailed temporal evolution of the event. This work demonstrates the superiority of the transfer-learning technique in satellite AOD retrievals and the applicability of the retrieved AGRI AOD in monitoring extreme pollution events.

Addressing Biases in Ambient PM2.5 Exposure and Associated Health Burden Estimates by Filling Satellite AOD Retrieval Gaps over India.

Ambient PM2.5 exposure statistics in countries with limited ground monitors are derived from satellite aerosol optical depth (AOD) products that have spatial gaps. Here, we quantified the biases in PM2.5 exposure and associated health burden in India due to the sampling gaps in AOD retrieved by a Moderate Resolution Imaging Spectroradiometer. We filled the sampling gaps and derived PM2.5 in recent years (2017-2022) over India, which showed fivefold cross-validation R2 of 0.92 and root mean square error (RMSE) of 11.8 μg m-3 on an annual scale against ground-based measurements. If the missing AOD values are not accounted for, the exposure would be overestimated by 19.1%, translating to an overestimation in the mortality burden by 93,986 (95% confidence interval: 78,638-110,597) during these years. With the gap-filled data, we found that the rising ambient PM2.5 trend in India has started showing a sign of stabilization in recent years. However, a reduction in population-weighted exposure balanced out the effect of the increasing population and maintained the mortality burden attributable to ambient PM2.5 for 2022 (991,058:798,220-1,183,896) comparable to the 2017 level (1,014,766:812,186-1,217,346). Therefore, a decline in exposure alone is not sufficient to significantly reduce the health burden attributable to ambient PM2.5 in India.

Diurnal time representation of MODIS, VIIRS, MISR, and AHI over Asia and Oceania

Retrievals of satellite aerosol optical depth (AOD) products are used for studies of air pollution, energy balance, and climate change. However, few studies have evaluated the diurnal time representation (DTR) of satellite AOD products. In this study, we proposed a statistically valid threshold method to investigate the DTR of Moderate Resolution Imaging Spectroradiometer (MODIS), Visible Infrared Imaging Radiometer Suite (VIIRS), Multi-angle Imaging Spectroradiometer (MISR), and Advanced Himawari Imager (AHI) AOD products by comparing their retrievals to the daily AERONET AOD values and explored the ability of AHI and AERONET to capture daily aerosol variability at 54 sites in Asia and Oceania from July 2015 to December 2019. For polar-orbiting satellites, based on collocations with AERONET measurements, MISR-derived AOD showed high consistency with daily average AERONET AOD, with a correlation coefficient (R) of 0.94 and 80.85% matchups falling into the envelopes of expected error (EE, ± (0.05 + 0.15 × AODAERONET)). Compared with single MODIS products, a combination of satellite data from the Terra and Aqua MODIS dark target (DT) algorithm might obtain a dataset with better DTR (R = 0.94 and %Within EE = 72.10). The satellite-retrieved AOD did a good job to represent the daily mean AERONET AOD, because AERONET AOD averaged within a ± 30 min time window of different satellite overpasses closely mirrored the daily mean AERONET AOD. Furthermore, the seasonal average diurnal variations (SADV) of AERONET AOD were changing insignificant, except for the urban coverage area in the Southern Hemisphere and equatorial climate. For geostationary satellite, the DTR of AHI was relatively poorer in comparison to the other sensors (R = 0.84, %Within EE = 58.00), and AHI AOD was seriously underestimated in South Asia and overestimated in Oceania. Besides, the SADVs of AHI AOD presented several patterns summarized as wave type, unimodal type, stable fluctuation type, and other types, which meant that AHI AOD retrievals did not accurately capture the diurnal variability of aerosols. This study provides a reference for the aerosol community to select satellite AOD products with appropriate DTR.

An integrated monitoring system (IMS) for air quality: Observations of a unique ozone-exceedance event in Maryland

Ozone (O3) exceedance events in the mid-Atlantic states of the U.S. are typically observed during the summer months. They frequently result from local diurnal emissions of NO2 and VOCs where photochemical production of O3 is most prevalent. Occasionally, summertime O3-exceedances in the region are driven by long-range transport of pollutants or regional nocturnal dynamics. Discussed herein is the use of multiple observation platforms as an integrated monitoring system (IMS) to characterize the spatiotemporal confluence of long-range transport and nocturnal low-level jet induced O3-exceedances impacting air quality in the mid-Atlantic region. The IMS discussed is comprised of enhanced monitoring of vertically resolved observations of O3, atmospheric backscatter, winds, and state variables (e.g., temperature and relative humidity). This is in combination with surface-O3 observations, satellite-retrieved aerosol optical depth (by the moderate resolution imaging spectroradiometer (MODIS)), and simulations of air mass transport by NOAA's hybrid single-particle Lagrangian trajectory (HYSPLIT) model. Here we discuss an exceedance event that impacted the Mid-Atlantic region on May 20, 2021, where seven air monitoring sites in Maryland were driven by the combination of nocturnal low-level jet downmixing of pollutants and long-range transport of biomass-burning pollutants/O3-precursors. The most notable features of this case are this event's deviation from decadal O3-exceedance seasonality trends, the large spatial context of high surface-O3 impacting traditionally non-exceedance sites (sites that have not exceeded in years), the arrival of a well-mixed polluted air mass into the region, and nocturnal increases in surface-O3 induced by NLLJ.

Data augmentation for bias correction in mapping PM2.5 based on satellite retrievals and ground observations

As most air quality monitoring sites are in urban areas worldwide, machine learning models may produce substantial estimation bias in rural areas when deriving spatiotemporal distributions of air pollutants. The bias stems from the issue of dataset shift, as the density distributions of predictor variables differ greatly between urban and rural areas. We propose a data-augmentation approach based on the multiple imputation by chained equations (MICE-DA) to remedy the dataset shift problem. Compared with the benchmark models, MICE-DA exhibits superior predictive performance in deriving the spatiotemporal distributions of hourly PM2.5 in the megacity (Chengdu) at the foot of the Tibetan Plateau, especially for correcting the estimation bias, with the mean bias decreasing from –3.4 µg/m3 to –1.6 µg/m3. As a complement to the holdout validation, the semi-variance results show that MICE-DA decently preserves the spatial autocorrelation pattern of PM2.5 over the study area. The essence of MICE-DA is strengthening the correlation between PM2.5 and aerosol optical depth (AOD) during the data augmentation. Consequently, the importance of AOD is largely enhanced for predicting PM2.5, and the summed relative importance value of the two satellite-retrieved AOD variables increases from 5.5% to 18.4%. This study resolved the puzzle that AOD exhibited relatively lower importance in local or regional studies. The results of this study can advance the utilization of satellite remote sensing in modeling air quality while drawing more attention to the common dataset shift problem in data-driven environmental research.

The Effect of Abrupt Changes to Sources of PM10 and PM2.5 Concentrations in Three Major Agglomerations in Mexico

In the three major urban agglomerations in Mexico (Mexico City, Monterrey, and Guadalajara), a significant change to anthropogenic sources of air pollution happened in March–May 2020, when policies implemented to stop the spread of the COVID-19 virus in Mexico caused the reduction of some anthropogenic sources of air pollution. We study the effect of these significant changes to air pollution sources using satellite-retrieved aerosol optical depth (AOD) and particulate matter (PM10 and PM2.5) concentrations from ground stations. The Chow test was applied to study trend changes in PM concentrations from 1 January to 30 May 2020. The Mann–Whitney non-parametric test was then used to compare average PM concentrations in April and May pre-lockdown, during lockdown in 2020, and post-lockdown in 2021. The assessment was further performed by evaluating the exceedance of national air quality standard maxima. The trend analysis showed that PM10 concentrations were reduced during lockdown in Mexico City and Monterrey, whereas no change was found for PM10 in Guadalajara and PM2.5 in the three cities. Further analysis showed that in Mexico City and Guadalajara, average PM10 and PM2.5 concentrations decreased by 12% in April and May 2020. However, in Monterrey, average PM10 and PM2.5 concentrations increased by 2.76% and 11.07%, respectively, in April 2021 due to a severe drought that caused dry soils and dust around the city. The results of this research can be used to implement policies for reducing anthropogenic sources to improve the air quality in urban areas.

Meteorological and anthropogenic contributions to changes in the Aerosol Optical Depth (AOD) over China during the last decade

To reduce air pollution over China, a series of measures to decrease concentrations of aerosols and trace gases has been implemented during the last two decades. The effects of these measures on aerosol concentrations, represented here by satellite-retrieved Aerosol Optical Depth (AOD), were evaluated by comparison of observed AOD time series from 2010 to 2021, over five representative areas, with model simulations using emissions fixed to those in 2010, but with actual meteorological data. Thus, model AOD was only influenced by meteorological conditions and the comparison with satellite AOD allows for the separation of meteorological and anthropogenic contributions to the AOD. The model data shows the influence of meteorological factors on the variation of the AOD, between +10% (unfavorable, leading to AOD increase) and −10% (favorable, leading to AOD decrease). Unfavorable meteorological conditions counteract effects of emission control measures. The satellite data show the occurrence of an AOD maximum in 2014 over three areas due to meteorological effects. The decrease of the AOD, varying between 22% and 40%, is mostly due to anthropogenic effects, but meteorological contributions are substantial and vary between the five focus regions. The satellite and model data also show that during 2017–2021 the AOD reduction flattened and AOD fluctuated within 10% of the average value for these years, with a substantial meteorological contribution. This analysis suggests that the Three-year Action Plan for Clean Air was less effective than the Clean Air Action Plan.

Difference between global and regional aerosol model classifications and associated implications for spaceborne aerosol optical depth retrieval

AbstractAerosol model are generally adopted to describe the physical and optical characteristics of aerosols in different regions, and this is important for climate and environmental studies as well as satellite aerosol retrieval. However, current mainstream aerosol retrieval algorithms still lack differentiated and dynamic descriptions of aerosol models, and the aerosol optical depth (AOD) retrieval accuracy is hence largely compromised. Therefore, in this study, we aimed to investigate the difference between global and regional aerosol model classifications and the impact of these differences on satellite AOD retrieval. We developed global and regional (for China) aerosol models through the K-means cluster analysis method, and used the 6SV radiative transfer model based on these global and regional aerosol models to simulate air pollution events. The cluster analysis revealed that both global and regional aerosols can be clustered into five main aerosol types with different optical and physical parameters. The aerosol types include strongly absorbing, moderately absorbing, weakly absorbing, dust and coarse-fine mixed aerosols. These two aerosol models exhibited distinct seasonal variations. Comparing the clustering results obtained with the two models at five Aerosol Robotic Network (AERONET) sites, we found that the proportion of strongly absorbing aerosols in the global results is low. According to the China clustering results, a higher proportion of strongly absorbing aerosols occurred in autumn, winter, and early spring, and dust aerosols attained an increasing proportion in spring. Comparing the aerosol parameters with observational data, the regional aerosol model performed better than the global aerosol model. Simulations of AOD during a pollution episode and in different seasons at the Beijing-CAMS site with the 6SV radiative transfer model further confirmed that the regional aerosol model achieved a better simulation accuracy than the global model. The findings suggest that regional aerosol models can provide a useful reference for satellite remote sensing algorithm, climate change and air pollution assessment, and that the dynamic description of the aerosol model is a key parameter in satellite AOD retrieval algorithms.

Satellite Aerosol Optical Depth Retrieval Based on Fully Connected Neural Network (FCNN) and a Combine Algorithm of Simplified Aerosol Retrieval Algorithm and Simplified and Robust Surface Reflectance Estimation (SREMARA)

Satellite Aerosol Optical Depth Retrieval Based on Fully Connected Neural Network (FCNN) and a Combine Algorithm of Simplified Aerosol Retrieval Algorithm and Simplified and Robust Surface Reflectance Estimation (SREMARA)

Regional PM2.5 Estimation for Southern Ontario through Geographically Weighted Regression

In this study, a geographically weighted regression (GWR) approach was adopted to forecast regional concentration of particulate matter 2.5 (PM2.5) for the southern Ontario based on both in situ meteorological measurement and Satellite retrievals of aerosol optical depth (AOD). The correlation between monitored concentration of PM2.5 and Satellite-retrieved AOD would be quantified. The ground-level PM2.5 for South Ontario area was then predicted using GWR with AOD and meteorological variables considered as inputs. The results indicated that performance of GWR was slightly better than the ordinary least squares (OLS) model, indicating spatial variations between independent and dependent variables. Consequently, the GWR model can help us to predict the PM2.5 concentration in terms of time or region with satellite data, and also help improve satellite data inversion.

Assessment of Satellite AOD during the 2020 Wildfire Season in the Western U.S.

Satellite remote sensing of aerosol optical depth (AOD) is essential for detection, characterization, and forecasting of wildfire smoke. In this work, we evaluate the AOD (550 nm) retrievals during the extreme wildfire events over the western U.S. in September 2020. Three products are analyzed, including the Moderate-resolution Imaging Spectroradiometers (MODIS) Multi-Angle Implementation of Atmospheric Correction (MAIAC) product collections C6.0 and C6.1, and the NOAA-20 Visible Infrared Imaging Radiometer (VIIRS) AOD from the NOAA Enterprise Processing System (EPS) algorithm. Compared with the Aerosol Robotic Network (AERONET) data, all three products show strong linear correlations with MAIAC C6.1 and VIIRS presenting overall low bias (<0.06). The accuracy of MAIAC C6.1 is found to be substantially improved with respect to MAIAC C6.0 that drastically underestimated AOD over thick smoke, which validates the effectiveness of updates made in MAIAC C6.1 in terms of an improved representation of smoke aerosol optical properties. VIIRS AOD exhibits comparable uncertainty with MAIAC C6.1 with a slight tendency of increased positive bias over the AERONET AOD range of 0.5–3.0. Averaging coincident retrievals from MAIAC C6.1 and VIIRS provides a lower root mean square error and higher correlation than for the individual products, motivating the benefit of blending these datasets. MAIAC C6.1 and VIIRS are further compared to provide insights on their retrieval strategy. When gridded at 0.1° resolution, MAIAC C6.1 and VIIRS provide similar monthly AOD distribution patterns and the latter exhibits a slightly higher domain average. On daily scale, over thick plumes near fire sources, MAIAC C6.1 reports more valid retrievals where VIIRS tends to have retrievals designated as low or medium quality, which tends to be due to internal quality checks. Over transported smoke near scattered clouds, VIIRS provides better retrieval coverage than MAIAC C6.1 owing to its higher spatial resolution, pixel-level processing, and less strict cloud masking. These results can be used as a guide for applications of satellite AOD retrievals during wildfire events and provide insights on future improvement of retrieval algorithms under heavy smoke conditions.

Spatiotemporal high-resolution imputation modeling of aerosol optical depth for investigating its full-coverage variation in China from 2003 to 2020

Investigating spatiotemporal variations of atmospheric aerosols is important for climate change and environmental research. Although satellite aerosol optical depth (AOD) retrieved by the MAIAC (Multiangle Implementation of Atmospheric Correct) algorithm provides a unique opportunity to represent global aerosol loading with high spatiotemporal resolution, accurate assessment of long-term aerosol loading countrywide is still challenging due to its non-random missingness. This study aimed to develop an adaptive spatiotemporal high-resolution imputation modeling framework for AOD that incorporates random forest models and multisource data (the simulated AOD, meteorological, and surface condition data) to support full-coverage long- and short-term aerosol studies in China. Aided by the time-stratified approach, the imputation model was constructed for each day, and the MAIAC AOD was used as the target variable. The proposed approach could effectively capture the massive spatiotemporal variability in a large amount of data and deliver full-coverage AODs with high accuracies at a daily timescale (i.e., overall validation R2 against ground-level AOD measurements of 0.77). We then employed the proposed approach to impute the daily MAIAC retrieved AOD towards complete coverage for China for 2003–2020. Due to the complete coverage, the spatial pattern of monthly/seasonal/yearly mean AOD imputations has better representativeness than that of original MAIAC retrievals. Comparison analysis shows that the monthly/seasonal/yearly aerosol loading over most of China tends to be underestimated by temporal aggregates of original satellite-retrieved AODs. Such underestimation is particularly severe in summer and over the North China Plain (the amount of underestimation >0.2). Consequently, our full-coverage AOD imputations can advance scientific research and environmental management by supporting national and local complete pictures of both short-term episodes and long-term trends in atmospheric aerosols.

Spatiotemporal continuous estimates of daily 1 km PM2.5 from 2000 to present under the Tracking Air Pollution in China (TAP) framework

Abstract. High spatial resolution PM2.5 data covering a long time period are urgently needed to support population exposure assessment and refined air quality management. In this study, we provided complete-coverage PM2.5 predictions with a 1 km spatial resolution from 2000 to the present under the Tracking Air Pollution in China (TAP, http://tapdata.org.cn/, last access: 3 October 2022) framework. To support high spatial resolution modeling, we collected PM2.5 measurements from both national and local monitoring stations. To correctly reflect the temporal variations in land cover characteristics that affected the local variations in PM2.5, we constructed continuous annual geoinformation datasets, including the road maps and ensemble gridded population maps, in China from 2000 to 2021. We also examined various model structures and predictor combinations to balance the computational cost and model performance. The final model fused 10 km TAP PM2.5 predictions from our previous work, 1 km satellite aerosol optical depth retrievals, and land use parameters with a random forest model. Our annual model had an out-of-bag R2 ranging between 0.80 and 0.84, and our hindcast model had a by-year cross-validation R2 of 0.76. This open-access, 1 km resolution PM2.5 data product, with complete coverage, successfully revealed the local-scale spatial variations in PM2.5 and could benefit environmental studies and policymaking.

An integrated approach of Belief Rule Base and Convolutional Neural Network to monitor air quality in Shanghai

Accurate monitoring of air quality can reduce its adverse impact on earth. Ground-level sensors can provide fine particulate matter (PM2.5) concentrations and ground images. But, such sensors have limited spatial coverage and require deployment cost. PM2.5 can be estimated from satellite-retrieved Aerosol Optical Depth (AOD) too. However, AOD is subject to uncertainties associated with its retrieval algorithms and constrain the spatial resolution of estimated PM2.5. AOD is not retrievable under cloudy weather as well. In contrast, satellite images provide continuous spatial coverage with no separate deployment cost. Accuracy of monitoring from such satellite images is hindered due to uncertainties of sensor data of relevant enviromental parameters, such as, relative humidity, temperature, wind speed and wind direction. Belief Rule Based Expert System (BRBES) is an efficient algorithm to address these uncertainties. Convolutional Neural Network (CNN) is suitable for image analytics. Hence, we propose a novel model by integrating CNN with BRBES to monitor air quality from satellite images with improved accuracy. We customized CNN and optimized BRBES to increase monitoring accuracy further. An obscure image has been differentiated between polluted air and cloud in our model. Valid environmental data (temperature, wind speed and wind direction) have been adopted to further strengthen the monitoring performance of our proposed model. Three-year observation data (satellite images and environmental parameters) from 2014 to 2016 of Shanghai have been employed to analyze and design our proposed model. The results conclude that the accuracy of our model to monitor PM2.5 of Shanghai is higher than only CNN and other conventional Machine Learning methods. Real-time validation of our model on near real-time satellite images of April-2021 of Shanghai shows average difference between our calculated PM2.5 concentrations and the actual one within ± 5.51.

Joint features random forest (JFRF) model for mapping hourly surface PM2.5 over China

Ambient PM2.5 exerts strong regional pattern for its ability for long-range transport, implying that including the features of surrounding stations may improve the accuracy of machine-learning based model to estimate surface PM2.5 from the satellite-retrieved aerosol optical depth (AOD). However, most of current models either just use single point features, or simply average the observed surface PM2.5 from adjacent stations based on a fixed spatial proportional relationship. The question that how to properly take advantage of the features of surrounding stations for retrieving PM2.5 is still not well addressed. Here we propose an integrated algorithm called joint features random forest (JFRF) model which includes complex feature differences with surrounding stations and the observation of stations to learn the dynamic relations with the PM2.5 of target pixel, rather than the weighted average feature (WAF) only by surface PM2.5 as traditional models (with WAF) used. Results of cross validation suggest better performance of JFRF (R2 = 0.61–0.8; RMSE = 15.97–20.91 μg/m3) than single point feature model (ΔR2 = 0.09–0.3). JFRF also exhibits better performance than traditional models (with WAF) (ΔR2 = 0.05–0.11), particularly in regions with large AOD gradient (accounts for 33% of the total test set), which is of great significance for accurately representing the spatial heterogeneity of PM2.5 (e.g., pollution edging and hot spots areas). And the exclusion of AOD from the features significantly reduced the model performance (ΔR2 = −0.07 ∼ −0.1). Therefore, our study demonstrates the important of the feature differences of surrounding stations and satellite-retrieved AOD in representing the regional pattern of PM2.5 and further helping the machine-learning based model to improve the accuracy in estimating surface PM2.5.

Evaluation and comparison of MODIS and VIIRS aerosol optical depth (AOD) products over regions in the Eastern Mediterranean and the Black Sea

Aerosol loading has direct and indirect impacts on radiative forcing and climate; additionally, exposure to high aerosol concentrations has adverse effects on human health. Several inversion algorithms are used to retrieve aerosol properties, primarily aerosol optical depth (AOD). Satellite aerosol products are available from several moderate spatial resolution sensors and provide information on the spatial and temporal dynamics of aerosol concentration. This study provides an evaluation and inter-comparison of multiple aerosol products from Terra and Aqua Moderate Resolution and Imaging Spectroradiometer (MODIS) and Suomi National Polar-Orbiting Partnership (Suomi NPP) Visible-Infrared Imager Radiometer Suite (VIIRS) between 2014 and 2018 in the Eastern Mediterranean and the Black Sea. This evaluation focuses on four MODIS aerosol products including Dark Target (DT) AOD retrievals at 10-km and 3-km spatial resolutions, Deep Blue (DB) AOD retrievals at10-km spatial resolution, Multi-Angle Implementation of Atmospheric Correction (MAIAC) AOD retrievals at 1-km spatial resolution, and VIIRS DB AOD retrievals at 6-km spatial resolution. Satellite-based AOD retrievals are validated against in-situ AOD measurements from three AERONET ground stations located in urban/land, rural/coastal, and ocean surfaces.The results show that Pearson's correlation coefficient (r) between the satellite AOD retrievals and AERONET AOD measurements range between 0.45 and 0.93 whilst the Root-Mean-Squared Errors (RMSE) varies from 0.047 to 0.159 which highlights the reliability of satellite AOD retrievals' in the region. The best performing products over the urban/land surface are the MODIS MAIAC and VIIRS DB aerosol products with the RMSE of 0.048–0.061 and a higher percentage (81%–89%) of retrievals falling within expected error for land. Over the ocean, the MODIS DT aerosol products have the lowest errors (RMSE 0.047–0.061) and the highest percentage (70%–80%) of retrievals falling within the expected error. The performances of all algorithms over the rural/coastal region are poorer in terms of over- and under-estimations of AOD with a small percentage of retrievals (48%–87%) within expected error. Among the satellite aerosol products, the VIIRS aerosol product performs better over the coastal area in terms of high correlation coefficient (r = 0.88), low error (RMSE = 0.059), and the large percentage of retrievals within expected error (88%). The poor performance found over the more complex coastal regions suggests further improvements are needed in land-water-sediment masks, surface reflectance estimation, and aerosol model assumptions over mixed surface types.The aerosol products show the highest agreement (r ∼ 0.77–0.98) over the ocean surface whilst the coastal site aerosol products verified the greatest dissimilarity in aerosol loading variations (r ∼ 0.50–0.91). Lastly, the good correlation (r > 0.73) between the VIIRS and MODIS AOD products illustrates the potential of VIIRS to provide aerosol data continuity with MODIS.