Aerosol Optical Depth Retrievals Research Articles

Intercomparison of aerosol optical depth retrievals from GAW-PFR and SKYNET sun photometer networks and the effect of calibration

Abstract. In this study, we assess the homogeneity of aerosol optical depth (AOD) between two sun photometer networks, the Global Atmosphere Watch-Precision Filter Radiometer (GAW-PFR) and the European SKYNET radiometers network (ESR), at the common wavelengths of their main instruments (500 and 870 nm). The main focus of this work is to evaluate the effect of the improved Langley plot calibration method (ILP) used by SKYNET and to investigate the factors affecting its performance. We used data from three intercomparison campaigns that took place during 2017–2021. Each campaign was organized at two locations (mountainous rural – Davos, Switzerland; urban – Rome, Italy). Our analysis shows that differences in AOD due to post-processing and instrument differences are minor. The main factor leading to AOD differences is the calibration method. We found a systematic underestimation of AOD in ESR compared to in GAW-PFR due to underestimation of the calibration constant calculated with the ILP method compared to the calibration transfers using the PFR as a reference. The calibration and AOD differences are smaller in Davos, where the traceability criteria are satisfied at 870 nm and where the median differences are below 0.01 at 500 nm. In Rome, the AOD median differences at 500 nm were in the 0.015–0.034 range. We conducted a sensitivity study, which shows that part of the difference can potentially be explained by errors in the assumed surface albedo and instrument solid-view angle provided as inputs to the ILP code (based on Skyrad pack 4.2). Our findings suggest that the ILP method is mainly sensitive to the measured sky radiance. The underestimation in calibration is probably caused by an error in the retrieved scattering AOD (sc-AOD) through the sky radiance inversion. Using an alternative retrieval method (Skyrad MRI pack version 2) to derive sc-AOD and to recalibrate the instruments with the ILP method, we found no significant differences between the retrieved sc-AOD and no systematic increase in the ILP-derived calibration constant when using the MRI pack for sc-AOD inversion instead of the Skyrad 4.2. The potential error may be a result of the model assumptions used for the sky radiance simulations. In conclusion, the on-site calibration of sun photometers has several advantages, including the fact that instrument shipments and data gaps can be avoided. However, it has also the disadvantages of a larger uncertainty and significant systematic differences compared to the traditional Langley calibration performed under low- and constant-AOD conditions at high-altitude sites. The larger uncertainty of the ILP method can be attributed to the required modelling and input parameters.

Aerosol optical depth retrieval using scaled digital number (DN) values of multi-spectral satellite and a generating adversarial model based on deep learning application

ABSTRACT The current quantitative retrieval of Aerosol Optical Depth (AOD) typically uses Top-of-atmosphere (TOA) reflectance data obtained by the scale and offset terms of radiometric calibration. Errors can be introduced during the conversion of Digital Number (DN) values to TOA reflectance, restraining the accurate retrieval of AOD. Particularly when surface reflectance is high, conversion errors can significantly distort the retrieval of Aerosol Optical Depth (AOD), given the minimal impact of aerosols on the radiation detected by satellite sensors. To counteract this, we introduce an innovative retrieval algorithm that harnesses a Generative Adversarial Network (GANN), a deep learning (DL) model, to accurately derive AOD from the scaled digital number (DN) values of multi-spectral satellite imagery. The DL technology enables the use of scaled DN values for AOD retrieval, bypassing the conversion errors associated with TOA reflectance since it relies on a sample learning approach instead of radiative transfer calculations. To ensure the number and representativeness of training samples, a total of 37,103 samples from different underlying surfaces from 2014 to 2018 were obtained using a space-time matching strategy. Our K-cross validation demonstrates that employing DN values for AOD retrieval promises to yield higher accuracy compared to using TOA reflectance. Independent 2019 data were used to estimate GANN, MCD19A2 and MOD04_3K aerosol products. Our model showed superior performance compared to the other products in terms of a higher correlation (R = 0.8184) with AERONET AOD and obtaining more reliable retrievals (7601). The sensitivity experiment suggests that GANN exhibits limitations in areas characterized by high aerosol loads and heterogeneity. Our study can give a novel insight into AOD retrieval based on multi-spectral satellite and DL models.

Post-process correction improves the accuracy of satellite PM2.5 retrievals

Abstract. Estimates of PM2.5 levels are crucial for monitoring air quality and studying the epidemiological impact of air quality on the population. Currently, the most precise measurements of PM2.5 are obtained from ground stations, resulting in limited spatial coverage. In this study, we consider satellite-based PM2.5 retrieval, which involves conversion of high-resolution satellite retrieval of aerosol optical depth (AOD) into high-resolution PM2.5 retrieval. To improve the accuracy of the AOD-to-PM2.5 conversion, we employ the machine-learning-based post-process correction to correct the AOD-to-PM conversion ratio derived from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) reanalysis model data. The post-process-correction approach utilizes a fusion and downscaling of satellite observation and retrieval data, MERRA-2 reanalysis data, various high-resolution geographical indicators, meteorological data, and ground station observations for learning a predictor for the approximation error in the AOD-to-PM2.5 conversion ratio. The corrected conversion ratio is then applied to estimate PM2.5 levels given the high-resolution satellite AOD retrieval data derived from Sentinel-3 observations. The region of study is central Europe during the year 2019. Our model produces PM2.5 estimates with a spatial resolution of 100 m at satellite overpass times with R2 = 0.55 and RMSE = 6.2 µg m−3. The corresponding metrics for monthly averages are R2 = 0.72 and RMSE = 3.7 µg m−3. Additionally, we have incorporated an ensemble of neural networks to provide error envelopes for machine-learning-related uncertainty in the PM2.5 estimates. The proposed approach can produce accurate high-resolution PM2.5 data that can be very useful for air quality monitoring, emission regulation, and epidemiological studies.

Assessment of the impact of NO2 contribution on aerosol-optical-depth measurements at several sites worldwide

Abstract. This work aims at investigating the effect of NO2 absorption on aerosol-optical-depth (AOD) measurements and Ångström exponent (AE) retrievals of sun photometers by the synergistic use of accurate NO2 characterization for optical-depth estimation from co-located ground-based measurements. The analysis was performed for ∼ 7 years (2017–2023) at several sites worldwide for the AOD measurements and AE retrievals by Aerosol Robotic Network (AERONET) sun photometers which use OMI (Ozone Monitoring Instrument) climatology for NO2 representation. The differences in AOD and AE retrievals by NO2 absorption are accounted for using high-frequency columnar NO2 measurements by a co-located Pandora spectroradiometer belonging to the Pandonia Global Network (PGN). NO2 absorption affects the AOD measurements in UV-Vis (visible) range, and we found that the AOD bias is the most affected at 380 nm by NO2 differences, followed by 440, 340, and 500 nm, respectively. AERONET AOD was found to be overestimated in half of the cases, while also underestimated in other cases as an impact of the NO2 difference from “real” (PGN NO2) values. Overestimations or underestimations are relatively low. About one-third of these stations showed a mean difference in NO2 and AOD (at 380 and 440 nm) above 0.5 × 10−4 mol m−2 and 0.002, respectively, which can be considered a systematic contribution to the uncertainties in the AOD measurements that are reported to be of the order of 0.01. However, under extreme NO2 loading scenarios (i.e. 10 % highest differences) at highly urbanized/industrialized locations, even higher AOD differences were observed that were at the limit of or higher than the reported 0.01 uncertainty in the AOD measurement. PGN NO2-based sensitivity analysis of AOD difference suggested that for PGN NO2 varying between 2 × 10−4 and 8 × 10−4 mol m−2, the median AOD differences were found to rise above 0.01 (even above 0.02) with the increase in NO2 threshold (i.e. the lower limit from 2 × 10−4 to 8 × 10−4 mol m−2). The AOD-derivative product, AE, was also affected by the NO2 correction (discrepancies between the AERONET OMI climatological representation of NO2 values and the real PGN NO2 measurements) on the spectral AOD. Normalized frequency distribution of AE (at 440–870 and 340–440 nm wavelength pair) was found to be narrower for a broader AOD distribution for some stations, and vice versa for other stations, and a higher relative error at the shorter wavelength (among the wavelength pairs used for AE estimation) led to a shift in the peak of the AE difference distribution towards a higher positive value, while a higher relative error at a lower wavelength shifted the AE difference distribution to a negative value for the AOD overestimation case, and vice versa for the AOD underestimation case. For rural locations, the mean NO2 differences were found to be mostly below 0.50 × 10−4 mol m−2, with the corresponding AOD differences being below 0.002, and in extreme NO2 loading scenarios, it went above this value and reached above 1.00 × 10−4 mol m−2 for some stations, leading to higher AOD differences but below 0.005. Finally, AOD and AE trends were calculated based on the original AERONET AOD (based on AERONET OMI climatological NO2), and its comparison with the mean differences in the AERONET and PGN NO2-corrected AOD was indicative of how NO2 correction could potentially affect realistic AOD trends.

Impact of Model Spatial Resolution on Global Geophysical Satellite-Derived Fine Particulate Matter.

Global geophysical satellite-derived ambient fine particulate matter (PM2.5) inference relies upon a geophysical relationship (η) from a chemical transport model to relate satellite retrievals of aerosol optical depth (AOD) to surface PM2.5. The resolution dependence of simulated η warrants further investigation. In this study, we calculate geophysical PM2.5 with simulated η from the GEOS-Chem model in its high-performance configuration (GCHP) at cubed-sphere resolutions of C360 (∼25 km) and C48 (∼200 km) and satellite AOD at 0.01° (∼1 km). Annual geophysical PM2.5 concentrations inferred from satellite AOD and GCHP simulations at ∼25 km and ∼200 km resolutions exhibit remarkable similarity (R 2 = 0.96, slope = 1.03). This similarity in part reflects opposite resolution responses across components with population-weighted normalized mean difference (PW-NMD) increasing by 5% to 11% for primary species while decreasing by -30% to -5% for secondary species at finer resolution. Despite global similarity, our results also identify larger resolution sensitivities of η over isolated pollution sources and mountainous regions, where spatial contrast of aerosol concentration and composition is better represented at fine resolution. Our results highlight the resolution dependence of representing near-surface concentrations and the vertical distribution of chemically different species with implications for inferring ground-level PM2.5 from columnar AOD.

Estimating the Impact of a 2017 Smoke Plume on Surface Climate Over Northern Canada With a Climate Model, Satellite Retrievals, and Weather Forecasts

AbstractIn August 2017, a smoke plume from wildfires in British Columbia and the Northwest Territories recirculated and persisted over northern Canada for over two weeks. We compared a full‐factorial set of NASA Goddard Institute for Space Studies ModelE simulations of the plume to satellite retrievals of aerosol optical depth and carbon monoxide, finding that ModelE performance was dependent on the model configuration, and more so on the choice of injection height approach, aerosol scheme and biomass burning emissions estimates than to the choice of horizontal winds for nudging. In particular, ModelE simulations with free‐tropospheric smoke injection, a mass‐based aerosol scheme and comparatively high fire NOx emissions led to unrealistically high aerosol optical depth. Using paired simulations with and without fire emissions, we estimated that for 16 days over an 850,000 km2 region, the smoke decreased planetary boundary layer heights by between 253 and 547 m, decreased downward shortwave radiation by between 52 and 172 Wm−2, and decreased surface temperature by between 1.5°C and 4.9°C, the latter spanning an independent estimate from operational weather forecasts of a 3.7°C cooling. The strongest surface climate effects were for ModelE configurations with more detailed aerosol microphysics that led to a stronger first indirect effect.

Comparisons of aerosol types and optical characters over Shouxian Area China observed from ground- and space-based systems

This study evaluates the performance of moderate-resolution Imaging spectroradiometer (MODIS) in aerosol optical depth(AOD) and Ångström exponent(AE) retrievals under high aerosol loading conditions across various aerosol types, utilizing ground-based and space-borne aerosol measurements in Shouxian, China. The intercomparison reveals cloud-aerosol LiDAR with orthogonal polarization's (CALIOP) efficacy in detecting significant aerosol layers and the refinement of sunphotometer-based aerosol type classification through CALIPSO, achieving approximately 80% accuracy. Analysis of 2016-2017 data indicates substantial aerosol presence in Shouxian, with monthly mean AODs ranging from 0.35 to 0.72 at 550 nm, significantly above the global average. The predominant aerosol types were mixed-type (54.8%), desert dust (21.2%), urban/industrial(15.5%), biomass-burning aerosol (6.4%), and continental aerosol (12.1%), with frequent observations of elevated long-range transported aerosol layers. MODIS AOD retrievals generally align with sunphotometer measurements but exhibit higher biases, especially with increasing AOD magnitudes. However, there is a notable difference between MODIS and sunphotometer aerosol AE measurements, with MODIS accurately assessing BBA but showing varied performance across other aerosol types. The combination of AOD and AE of the DD aerosol type is the most accurate. Further analysis showed that MODIS AOD biases and AE biases are negatively correlated, these negative bias correlations show strong aerosol type sensitivities. Monthly analysis of MODIS and sunphotometer comparisons highlights varying performance, particularly during normalized difference vegetation index (NDVI) transitions, suggesting that local vegetation cycles and associated surface spectral reflectance changes significantly impact MODIS aerosol retrieval accuracy under high aerosol loading conditions.

An improved meteorological variables-based aerosol optical depth estimation method by combining a physical mechanism model with a two-stage model

A two-stage model integrating a spatiotemporal linear mixed effect (STLME) and a geographic weight regression (GWR) model is proposed to improve the meteorological variables-based aerosol optical depth (AOD) retrieval method (Elterman retrieval model-ERM). The proposed model is referred to as the STG-ERM model. The STG-ERM model is applied over the Beijing-Tianjin-Hebei (BTH) region in China for the years 2019 and 2020. The results show that data coverage increased by 39.0% in 2019 and 40.5% in 2020. Cross-validation of the retrieval results versus multi-angle implementation of atmospheric correction (MAIAC) AOD shows the substantial improvement of the STG-ERM model over earlier meteorological models for AOD estimation, with a determination coefficient (R2) of daily AOD of 0.86, root mean squared prediction error (RMSE) and the relative prediction error (RPE) of 0.10 and 36.14% in 2019 and R2 of 0.86, RMSE of 0.12 and RPE of 37.86% in 2020. The fused annual mean AOD indicates strong spatial variation with high value in south plain and low value in northwestern mountainous areas of the BTH region. The overall spatial seasonal mean AOD ranges from 0.441 to 0.586, demonstrating strongly seasonal variation. The coverage of STG-ERM retrieved AOD, as determined in this exercise by leaving out part of the meteorological data, affects the accuracy of fused AOD. The coverage of the meteorological data has smaller impact on the fused AOD in the districts with low annual mean AOD of less than 0.35 than that in the districts with high annual mean AOD of greater than 0.6. If available, continuous daily meteorological data with high spatiotemporal resolution can improve the model performance and the accuracy of fused AOD. The STG-ERM model may serve as a valuable approach to provide data to fill gaps in satellite-retrieved AOD products.

Towards long-term, high-accuracy, and continuous satellite total and fine-mode aerosol records: Enhanced Land General Aerosol (e-LaGA) retrieval algorithm for VIIRS

Long-term, accurate, stable, and continuous aerosol records from space are a major requirement for climate and atmospheric environment research. Due to the limited period span of the single satellite mission, solving the problem requires the combination of multiple satellite missions. In this study, a novel aerosol retrieval algorithm, named enhanced Land General Aerosol Retrieval (e-LaGA), for Visible/Infrared Imager Radiometer Suite (VIIRS) was developed to continue Moderate-resolution Imaging Spectro-radiometer (MODIS) aerosol retrieval. e-LaGA utilizes a Surface Reflectance (SR) relationship model map to more accurately estimate the SR parameter, optimizes the priori aerosol-type map, and improves the multi-band retrieval strategy using the residual-interpolation approach. Additionally, e-LaGA expands the retrieval ability of the traditional Dark Target (DT) algorithm over bright surfaces, and can more accurately retrieve Aerosol Optical Depth (AOD), Fine-mode AOD (AODF), and Fraction (FMF). The validation using global 533 AERONET sites shows that the volume of matchups of AOD (AODF) is approximately 80000, the correlation coefficient (R) is 0.902 (0.879) and the fraction of meeting the expected error envelope (±0.05 ± 0.15τ) is 0.806 (0.804). The inter-comparison indicates that the accuracy of e-LaGA AOD retrievals is comparable to seven commonly used AOD products. The validation performance and spatial distribution pattern of e-LaGA AODF and FMF are in good agreement with the POLDER (Polarization and Directionality of the Earth's Reflectances) products. e-LaGA is also applied to MODIS sensors. The VIIRS e-LaGA AOD, AODF, and FMF retrievals have high consistency with MODIS e-LaGA. Their average bias of AOD is 0.006, which is smaller than 0.024 for the Deep Blue and 0.011 for the DT. These preliminary results demonstrate the robustness of the e-LaGA algorithm and its potential for establishing a long-term climate record of total and fine-mode aerosol by combining multiple satellite missions, which is expected to reduce the uncertainty of climate change research.

DustSCAN: A Five Year (2018-2022) Hourly Dataset of Dust Plumes From SEVIRI

Airborne mineral dust significantly impacts air quality, human health, and the global climate. Due to sparse ground sensors, particularly in source regions, dust monitoring relies mainly on remote sensing through Aerosol Optical Depth (AOD) retrievals from polar-orbiting satellite optical instruments. These are valuable but lack the temporal resolution for precise plume tracking and source characterization. We introduce DustSCAN, a five-year, hourly dust plume dataset derived from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) images on geostationary-orbit Meteosat satellites. Using multi-channel infrared images, we detect atmospheric dust and track hourly dust-affected pixels. These are clustered into discrete plumes using the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm. DustSCAN includes 9950 discrete plumes over 2018-2022 across the Sahara, the Arabian Desert, and Western and Central Asia. It complements existing resources and provides a framework for detailed analysis of dust sources, trajectories, and impacts. Its distinctive event-based and spatio-temporal detail offers an advancement in unraveling the complexities of dust storm dynamics.

A Bayesian Framework to Quantify Uncertainty in Aerosol Optical Model Selection Applied to TROPOMI Measurements

This article presents a method within a Bayesian framework for quantifying uncertainty in satellite aerosol remote sensing when retrieving aerosol optical depth (AOD). By using a Bayesian model averaging technique, we take into account uncertainty in aerosol optical model selection and also obtain a shared inference about AOD based on the best-fitting optical models. In particular, uncertainty caused by forward-model approximations has been taken into account in the AOD retrieval process to obtain a more realistic uncertainty estimate. We evaluated a model discrepancy, i.e., forward-model uncertainty, empirically by exploiting the residuals of model fits and using a Gaussian process to characterise the discrepancy. We illustrate the method with examples using observations from the TROPOspheric Monitoring Instrument (TROPOMI) on the Sentinel-5 Precursor satellite. We evaluated the results against ground-based remote sensing aerosol data from the Aerosol Robotic Network (AERONET).

Evaluation and comparison of MODIS and MISR aerosol products with ground-based monitoring stations in the Amazon Basin

This study evaluates the long-term accuracy of six daily aerosol optical depth (AOD) datasets derived from the Moderate Resolution Imaging Spectroradiometers (MODIS) Dark Target (DT) at 3 and 10 km resolution, Deep Blue (DB), combined DT and DB datasets (DTB) and MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithms and the Multi-angle Imaging SpectroRadiometer (MISR) Version 23, by comparing them with ground-based AOD observations in 18 Aerosol Robotic Network (AERONET) sites over the Amazon basin between 2000 and 2022. The AOD retrieval accuracy was also assessed under varying seasons, aerosol loading, particle size, elevation and land cover type. Overall results showed that Terra and Aqua-MODIS DB and MISR algorithms perform better, with about 88% of AOD retrievals within the expected error range. However, MISR presented a systematic high negative bias. MODIS DT 3 km overestimated AOD values and indicated unsatisfactory retrieval accuracy in almost all scenarios. MODIS and MISR retrievals showed improved performance under low aerosol loading, whereas Aqua-MODIS MAIAC under intermediate conditions (AOD>0.2). MISR indicated a large dependence on a coarse mode-dominated environment with considerable negative bias. Terra-MODIS MAIAC retrievals exhibited a smaller dependence for aerosol particle size distribution, while Terra-MODIS DB performed better for fine mode and Aqua-MODIS DB for coarse dominance. MAIAC also performed well under polluted scenarios in the dry season in contrast to the majority of satellite-based AOD retrievals, which showed good agreement under clean background conditions in the wet season. Almost all satellite-based retrievals showed higher uncertainty under extremely high altitudes (>3000 m) followed by non-satisfactory correlation and significant biases. Additionally, almost all satellite-based AOD retrievals indicated poor performance as vegetation coverage increases. On the contrary, MISR showed higher accuracy under forest type. In general, MODIS and MISR retrieval accuracy demonstrated distinct levels of dependencies that require further improvement in the Amazon Basin.

Retrievals of aerosol optical depth over the western North Atlantic Ocean during ACTIVATE

Abstract. Aerosol optical depth was retrieved from two airborne remote sensing instruments, the Research Scanning Polarimeter (RSP) and Second Generation High Spectral Resolution Lidar (HSRL-2), during the National Aeronautics and Space Administration (NASA) Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). The field campaign offers a unique opportunity to evaluate an extensive 3-year dataset under a wide range of meteorological conditions from two instruments on the same platform. However, a long-standing issue in atmospheric field studies is that there is a lack of reference datasets for properly validating field measurements and estimating their uncertainties. Here we address this issue by using the triple collocation method, in which a third collocated satellite dataset from the Moderate Resolution Imaging Spectroradiometer (MODIS) is introduced for comparison. HSRL-2 is found to provide a more accurate retrieval than RSP over the study region. The error standard deviation of HSRL-2 with respect to the ground truth is 0.027. Moreover, this approach enables us to develop a simple, yet efficient, quality control criterion for RSP data. The physical reasons for the differences in two retrievals are determined to be cloud contamination, aerosols near the surface, multiple aerosol layers, absorbing aerosols, non-spherical aerosols, and simplified retrieval assumptions. These results demonstrate the pathway for optimal aerosol retrievals by combining information from both lidars and polarimeters for future airborne and satellite missions.

Retrieval of high-resolution aerosol optical depth (AOD) using Landsat 8 imageries over different LULC classes over a city along Indo-Gangetic Plain, India.

Aerosol optical depth (AOD) serves as a crucial indicator for assessing regional air quality. To address regional and urban pollution issues, there is a requirement for high-resolution AOD products, as the existing data is of very coarse resolution. To address this issue, we retrieved high-resolution AOD over Kanpur (26.4499°N, 80.3319°E), located in the Indo-Gangetic Plain (IGP) region using Landsat 8 imageries and implemented the algorithm SEMARA, which combines SARA (Simplified Aerosol Retrieval Algorithm) and SREM (Simplified and Robust Surface Reflectance Estimation). Our approach leveraged the green band of the Landsat 8, resulting in an impressive spatial resolution of 30m of AOD and rigorously validated with available AERONET observations. The retrieved AOD is in good agreement with high correlation coefficients (r) of 0.997, a low root mean squared error of 0.035, and root mean bias of - 4.91%. We evaluated the retrieved AOD with downscaled MODIS (MCD19A2) AOD products across various land classes for cropped and harvested period of agriculture cycle over the study region. It is noticed that over the built-up region of Kanpur, the SEMARA algorithm exhibits a stronger correlation with the MODIS AOD product compared to vegetation, barren areas and water bodies. The SEMARA approach proved to be more effective for AOD retrieval over the barren and built-up land categories for harvested period compared with the cropping period. This study offers a first comparative examination of SEMARA-retrieved high-resolution AOD and MODIS AOD product over a station of IGP.

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.

Local and Regional Diurnal Variability of Aerosol Properties Retrieved by DSCOVR/EPIC UV Algorithm

AbstractThe hour‐to‐hour variability of 388 nm aerosol optical depth (AOD) and single scattering albedo (SSA) derived from near UV observations by the Earth Polychromatic Imaging Camera (EPIC) on the Deep Space Climate Observatory has been evaluated at multiple locations around the world. AOD retrievals by the EPIC near UV algorithm (EPICAERUV) have been compared to ground based AOD measurements at 16 Aerosol Robotic Network (AERONET) stations representative of the most commonly observed aerosol types over geographic regions in three continents. Obtained results show that, in general, the EPICAERUV algorithm reproduces closely the hour‐to‐hour AOD variability reported by AERONET ground‐truth observations. Although most sites in the analysis show high correlation between the AOD hourly measurements by the ground‐based and space‐borne measuring techniques. Best algorithm performance is observed in the presence of carbonaceous and desert dust aerosols. The diurnal cycle of the retrieved SSA product was also analyzed. Although, a direct comparison of hourly EPICAERUV retrievals to equivalent ground‐based observations was not possible, the satellite result shows that diurnal SSA variability as large as 0.05 can be observed mostly associated with carbonaceous aerosols. EPICAERUV observed diurnal cycle of retrieved AOD on a regional basis was examined for the unusually active seasons of aerosol production of Saharan desert dust aerosols in 2020, and during the 2023 Canadian wildfires. Results presented in this study confirm the EPIC near UV aerosol product is well suited for observing diurnal variability of aerosols and, therefore, it is an important resource for climate and air quality studies.

Parameterizing spectral surface reflectance relationships for the Dark Target aerosol algorithm applied to a geostationary imager

Abstract. Originally developed for the moderate resolution imaging spectroradiometer (MODIS) in polar, sun-synchronous low earth orbit (LEO), the Dark Target (DT) aerosol retrieval algorithm relies on the assumption of a surface reflectance parameterization (SRP) over land surfaces. Specifically for vegetated and dark-soiled surfaces, values of surface reflectance in blue and red visible-wavelength bands are assumed to be nearly linearly related to each other and to the value in a shortwave infrared (SWIR) wavelength band. This SRP also includes dependencies on scattering angle and a normalized difference vegetation index computed from two SWIR bands (NDVISWIR). As the DT retrieval algorithm is being ported to new sensors to continue and expand the aerosol data record, we assess whether the MODIS-assumed SRP can be used for these sensors. Here, we specifically assess SRP for the Advanced Baseline Imager (ABI) aboard the Geostationary Operational Environmental Satellite (GOES)-16/East (ABIE). First, we find that using MODIS-based SRP leads to higher biases and artificial diurnal signatures in aerosol optical depth (AOD) retrievals from ABIE. The primary reason appears to be that the geostationary orbit (GEO) encounters an entirely different set of observation geometry than does LEO, primarily with regard to solar angles coupled with fixed-view angles. Therefore, we have developed a new SRP for GEO that draws the angular shape of the surface bidirectional reflectance. We also introduce modifications to the parameterization of both red–SWIR and blue–red spectral relationships to include additional information. The revised red–SWIR SRP includes the solar zenith angle, NDVISWIR, and land-type percentage from an ancillary database. The blue–red SRP adds dependencies on the scattering angle and NDVISWIR. The new SRPs improve the AOD retrieval of ABIE in terms of overall less bias and mitigation of the overestimation around local noon. The average bias of the DT AOD compared to the Aerosol Robotic Network (AERONET) AOD shows a reduction from 0.08 to 0.03, while the bias of local solar noon decreases from 0.12 to 0.03. The agreement between the DT and AERONET AOD is established through a regression slope of 1.06 and a y intercept of 0.01 with a correlation coefficient of 0.74. By using the new SRP, the percentage of data falling within the expected error range (±0.05 % + 15 %) is notably increased from 54 % to 78 %.

Time series retrieval of Multi-wavelength Aerosol optical depth by adapting Transformer (TMAT) using Himawari-8 AHI data

Accurate aerosol optical depth (AOD) with high temporal resolution is required for dynamic monitoring changes of aerosol properties. The Himawari-8 Advanced Himawari imager (AHI) data with 10-min temporal resolution has been widely used for AOD retrieval but usually involves radiative transfer models and assumptions on surface reflection and aerosol properties that are hard to satisfy. This paper introduces a data-driven algorithm without any priori assumptions on surface or aerosol properties called Time series retrieval of Multi-wavelength AOD by adapting Transformer (TMAT). TMAT directly retrieves AOD at eight-wavelengths every 10-min throughout the day, using AHI top of atmosphere (TOA) reflectances in 6 wavebands and ancillary data. TMAT explores diurnal temporal change patterns of AOD and surface reflection using time series data. TMAT was trained and evaluated using samples derived by matching a year (2019) of AHI full-disk data and AERONET/SONET observations with 12,768 diurnal time series, each of which includes up to 81 10-min 6-band AHI TOA observations and 81 (time) × 8 (wavelength) AOD values per day. There are up to 273,413 (<12,768 × 81) TOA observations in 6 wavebands and up to 1,420,716 (<12,768 × 81 × 8) AOD values in all the time series samples, with filled values (−9999) to represent no data due to cloud contamination and incomplete wavelength configurations of AOD stations. The masking mechanism in Transformer is leveraged to exclude the contribution of the AHI filled values in time series AOD retrieval and exclude the contribution of AOD filled values in training. TMAT was trained with a physical constraint that the retrieved AOD decreases toward longer wavelengths. TMAT AOD was evaluated by separating training and validation samples using site-specific leave-one-station-out validation. TMAT AOD shows extremely high accuracy at all wavelengths, with RMSE = 0.101, R = 0.949 and 81.69% of retrievals within ±0.05+15%∗AODsites for the AOD at 550 nm. We also validated the accuracy of the NASA multiangle implementation of atmospheric correction (MAIAC) AHI AOD over these sites with RMSE = 0.186, R = 0.832 and 52.10% of retrievals within ±0.05+15%∗AODsites. The Ångström exponent (AE) derived from the multi-wavelength TMAT AOD also shows good agreement with AERONET/SONET AE, with a correlation coefficient of 0.620 and an RMSE of 0.291. The TMAT effectiveness to capture diurnal temporal change patterns of surface reflection and AOD was evidenced by the significant reduction of AOD accuracy when using a non-time series deep learning method. The effectiveness of multi-wavelength joint retrieval was evidenced by the significant reduction of AE accuracy when using single wavelength retrieval models. Evidence shows the training and evaluation split methods have a large impact on evaluation accuracy and this study further used ground network independent training and evaluation split to confirm the leave-one-station-out validation accuracy. Despite this, future study will focus on more spatially independent validation and more comprehensive comparison between the MAIAC and TMAT results.

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.

Improving Dust Aerosol Optical Depth (DAOD) Retrieval from the GEOKOMPSAT-2A (GK-2A) Satellite for Daytime and Nighttime Monitoring.

The Advanced Meteorological Image (AMI) onboard GEOKOMPSAT 2A (GK-2A) enables the retrieval of dust aerosol optical depth (DAOD) from geostationary satellites using infrared (IR) channels. IR observations allow the retrieval of DAOD and the dust layer altitude (24 h) over surface properties, particularly over deserts. In this study, dust events in northeast Asia from 2020 to 2021 were investigated using five GK-2A thermal IR bands (8.7, 10.5, 11.4, 12.3, and 13.3 μm). For the dust cloud, the brightness temperature differences (BTDs) of 10.5 and 12.3 μm were consistently negative, while the BTD of 8.7 and 10.5 μm varied based on the dust intensity. This study exploited these optical properties to develop a physical approach for DAOD lookup tables (LUTs) using IR channels to retrieve the DAOD. To this end, the characteristics of thermal radiation transfer were simulated using the forward model; dust aerosols were explained by BTD (10.5, 12.3 μm)-an intrinsic characteristic of dust aerosol. The DAOD and dust properties were gained from a brightness temperature (BT) of 10.5 μm and BTD of 10.5, 12.3 μm. Additionally, the cumulative distribution function (CDF) was employed to strengthen the continuity of 24-h DAOD. The CDF was applied to the algorithm by calculating the conversion value coefficient for the DAOD error correction of the IR, with daytime visible aerosol optical depth as the true value. The results show that the DAOD product can be successfully applied during the daytime and nighttime to continuously monitor the flow of yellow dust from the GK-2A satellite in northeast Asia. In particular, the validation results for IR DAOD were similar to the active satellite product (CALIPSO/CALIOP) results, which exhibited a tendency similar to that for IR DAOD at night.