This study presents a novel methodology for the remote monitoring of aerosol
components over large spatial and temporal domains. The concept is realized
within the GRASP (Generalized Retrieval of ...Aerosol and Surface Properties)
algorithm to directly infer aerosol components from the measured radiances.
The observed aerosols are assumed to be mixtures of hydrated soluble
particles embedded with black carbon, brown carbon, iron oxide, and other
(non-absorbing) insoluble inclusions. The complex refractive indices of the
dry components are fixed a priori (although the refractive index of the
soluble host is allowed to vary with hydration), and the complex refractive
indices of the mixture are computed using mixing rules. The volume fractions
of these components are derived along with the size distribution and the
fraction of spherical particles, as well as the spectral surface reflectance in
cases when the satellite data are inverted. The retrieval is implemented as a
statistically optimized fit in a continuous space of solutions. This
contrasts with most conventional approaches in which the type of aerosol is
either associated with a pre-assumed aerosol model that is included in a set
of look-up tables, or determined from the analysis of the retrieved aerosol
optical parameters (e.g., single scattering albedo, refractive index, among others, provided by the AERONET retrieval algorithm); here, we retrieve the aerosol
components explicitly. The approach also bridges directly to the quantities
used in global chemical transport models. We first tested the approach with
synthetic data to estimate the uncertainty, and then applied it to real
ground-based AERONET and spaceborne POLDER/PARASOL observations; thus, the
study presents a first attempt to derive aerosol components from satellite observations specifically tied to global chemical transport model quantities. Our results
indicate aerosol optical characteristics that are highly consistent with
standard products (e.g., R of ∼0.9 for aerosol optical
thickness) and demonstrate an ability to separate intrinsic optical
properties of fine- and coarse-sized aerosols. We applied our method to
POLDER/PARASOL radiances on the global scale and obtained spatial and
temporal patterns of the aerosol components that agree well with existing
knowledge on aerosol sources and transport features. Finally, we discuss
limitations and perspectives of this new technique.
Quantitative estimations of atmospheric aerosol absorption are rather uncertain due to the lack of reliable information about the global distribution. Because the information about aerosol properties ...is commonly provided by single-viewing photometric satellite sensors that are not sensitive to aerosol absorption. Consequently, the uncertainty in aerosol radiative forcing remains one of the largest in the Assessment Reports of the Intergovernmental Panel on Climate Change (IPCC AR5 and AR6). Here, we use multi-angular polarimeters (MAP) to provide constraints on emission of absorbing aerosol species and estimate global aerosol absorption optical depth (AAOD) and its climate effect. Our estimate of modern-era mid-visible AAOD is 0.0070 that is higher than IPCC by a factor of 1.3-1.8. The black carbon instantaneous direct radiative forcing (BC DRF) is +0.33 W/m
+0.17, +0.54. The MAP constraint narrows the 95% confidence interval of BC DRF by a factor of 2 and boosts confidence in its spatial distribution.
We invert global black carbon (BC), organic carbon (OC) and desert dust (DD) aerosol emissions from POLDER/PARASOL spectral aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) ...using the GEOS-Chem inverse modeling framework. Our inverse modeling framework uses standard a priori emissions to provide a posteriori emissions that are constrained by POLDER/PARASOL AODs and AAODs. The following global emission values were retrieved for the three aerosol components: 18.4 Tg/yr for BC, 109.9 Tg/yr for OC and 731.6 Tg/yr for DD for the year 2010. These values show a difference of +166.7 %, +184.0 % and −42.4 %, respectively, with respect to the a priori values of emission inventories used in “standard” GEOS-Chem runs. The model simulations using a posteriori emissions (i.e., retrieved emissions) provide values of 0.119 for global mean AOD and 0.0071 for AAOD at 550 nm, which are +13.3 % and +82.1 %, respectively, higher than the AOD and AAOD obtained using the a priori values of emissions. Additionally, the a posteriori model simulation of AOD, AAOD, single scattering albedo, Ångström exponent and absorption Ångström exponent show better agreement with independent AERONET, MODIS and OMI measurements than the a priori simulation. Thus, this study suggests that using satellite-constrained global aerosol emissions in aerosol transport models can improve the accuracy of simulated global aerosol properties.
Proven by multiple theoretical and practical studies, multi-angular spectral polarimetry is ideal for comprehensive retrieval of properties of aerosols. Furthermore, a large number of advanced space ...polarimeters have been launched recently or planned to be deployed in the coming few years (Dubovik et al., 2019). Nevertheless, at present, practical utilization of aerosol products from polarimetry is rather limited, due to the relatively small number of polarimetric compared to photometric observations, as well as challenges in making full use of the extensive information content available in these complex observations. Indeed, while in recent years several new algorithms have been developed to provide enhanced aerosol retrievals from satellite polarimetry, the practical value of available aerosol products from polarimeters yet remains to be proven. In this regard, this paper presents the analysis of aerosol products obtained by the Generalized Retrieval of Atmosphere and Surface Properties (GRASP) algorithm from POLDER/PARASOL observations. After about a decade of development, GRASP has been adapted for operational processing of polarimetric satellite observations and several aerosol products from POLDER/PARASOL observations have been released. These updated PARASOL/GRASP products are publicly available (e.g., http://www.icare.univ-lille.fr, last access: 16 October 2018, http://www.grasp-open.com/products/, last access: 28 March 2020); the dataset used in the current study is registered under https://doi.org/10.5281/zenodo.3887265 (Chen et al., 2020).
Biomass burning is the main global source of fine primary carbonaceous aerosols in the form of organic carbon (OC) and black carbon (BC). We present an approach to estimate biomass burning aerosol ...emissions based on the measurement of radiative energy released during combustion. We make use of both Aqua and Terra MODIS observations to estimate the fire radiative energy using a simple model to parameterize the fire diurnal cycle based on the long‐term ratio between Terra and Aqua MODIS FRP. The parameterization is developed using cases of frequent (up to 12 times daily) MODIS observations, geostationary data from SEVIRI, and precessing observations from TRMM VIRS. FRE‐based emission coefficients for the organic and black carbon (OCBC) component of fine mode aerosols are computed from multiple regions encompassing grassland/savanna, tropical forest, and extratropical forest biomes using OCBC emission estimates derived from the MODIS fine mode aerosol product and an inverse aerosol transport model. The values of emission coefficients for OCBC retrieved were 2.7 ± 0.3 g/MJ for grassland/savanna, 8.6 ± 0.8 g/MJ for tropical forest, and 14.4 ± 0.8 g/MJ for extratropical forest. The FRE monthly data are then used to estimate OCBC emissions from biomass burning on a global basis. For 2001 to 2007, our annual estimates are comparable to previously published values. According to our estimate, the OCBC emissions are the largest for 2003 (18.8 Tg), roughly 20% above average and primarily driven by wildland fires in the Lake Baikal region (Russia).
The possibility of using shape mixtures of randomly oriented spheroids for modeling desert dust aerosol light scattering is discussed. For reducing calculation time, look‐up tables were simulated for ...quadrature coefficients employed in the numerical integration of spheroid optical properties over size and shape. The calculations were done for 25 bins of the spheroid axis ratio ranging from ∼0.3 (flattened spheroids) to ∼3.0 (elongated spheroids) and for 41 narrow size bins covering the size parameter range from ∼0.012 to ∼625. The look‐up tables were arranged into a software package, which allows fast, accurate, and flexible modeling of scattering by randomly oriented spheroids with different size and shape distributions. In order to evaluate spheroid model and explore the possibility of aerosol shape identification, the software tool has been integrated into inversion algorithms for retrieving detailed aerosol properties from laboratory or remote sensing polarimetric measurements of light scattering. The application of this retrieval technique to laboratory measurements by Volten et al. (2001) has shown that spheroids can closely reproduce mineral dust light scattering matrices. The spheroid model was utilized for retrievals of aerosol properties from atmospheric radiation measured by AERONET ground‐based Sun/sky‐radiometers. It is shown that mixtures of spheroids allow rather accurate fitting of measured spectral and angular dependencies of observed intensity and polarization. Moreover, it is shown that for aerosol mixtures with a significant fraction of coarse‐mode particles (radii ≥ ∼1 μm), the nonsphericity of aerosol particles can be detected as part of AERONET retrievals. The retrieval results indicate that nonspherical particles with aspect ratios ∼1.5 and higher dominate in desert dust plumes, while in the case of background maritime aerosol spherical particles are dominant. Finally, the potential of using AERONET derived spheroid mixtures for modeling the effects of aerosol particle nonsphericity in other remote sensing techniques is discussed. For example, the variability of lidar measurements (extinction to backscattering ratio and signal depolarization ratio) is illustrated and analyzed. Also, some potentially important differences in the sensitivity of angular light scattering to parameters of nonspherical versus spherical aerosols are revealed and discussed.
We have developed a method for evaluating the fidelity of the Aerosol Robotic Network (AERONET) retrieval algorithms by mimicking atmospheric extinction and radiance measurements in a laboratory ...experiment. This enables radiometric retrievals that use the same sampling volumes, relative humidities, and particle size ranges as observed by other in situ instrumentation in the experiment. We use three Cavity Attenuated Phase Shift (CAPS) monitors for extinction and University of Maryland Baltimore County’s (UMBC) three-wavelength Polarized Imaging Nephelometer (PI-Neph) for angular scattering measurements. We subsample the PI-Neph radiance measurements to angles that correspond to AERONET almucantar scans, with simulated solar zenith angles ranging from 50 ∘ to 77 ∘ . These measurements are then used as input to the Generalized Retrieval of Aerosol and Surface Properties (GRASP) algorithm, which retrieves size distributions, complex refractive indices, single-scatter albedos, and bistatic LiDAR ratios for the in situ samples. We obtained retrievals with residuals less than 8% for about 90 samples. Samples were alternately dried or humidified, and size distributions were limited to diameters of less than 1.0 or 2.5 μ m by using a cyclone. The single-scatter albedo at 532 nm for these samples ranged from 0.59 to 1.00 when computed with CAPS extinction and Particle Soot Absorption Photometer (PSAP) absorption measurements. The GRASP retrieval provided single-scatter albedos that are highly correlated with the in situ single-scatter albedos, and the correlation coefficients ranged from 0.916 to 0.976, depending upon the simulated solar zenith angle. The GRASP single-scatter albedos exhibited an average absolute bias of +0.023–0.026 with respect to the extinction and absorption measurements for the entire dataset. We also compared the GRASP size distributions to aerodynamic particle size measurements, using densities and aerodynamic shape factors that produce extinctions consistent with our CAPS measurements. The GRASP effective radii are highly correlated (R = 0.80) and biased under the corrected aerodynamic effective radii by 1.3% (for a simulated solar zenith angle of θ ∘ = 50 ∘ ); the effective variance indicated a correlation of R = 0.51 and a relative bias of 280%. Finally, our apparatus was not capable of measuring backscatter LiDAR ratios, so we measured bistatic LiDAR ratios at a scattering angle of 173 degrees. The GRASP bistatic LiDAR ratios had correlations of 0.71 to 0.86 (depending upon simulated θ ∘ ) with respect to in situ measurements, positive relative biases of 2–10%, and average absolute biases of 1.8–7.9 sr.
We present a method for simultaneously retrieving aerosol and surface parameters from ground-based and satellite observations collocated in space and time. We show that a combination of down and ...up-looking observations provides sufficient measurement constraints for characterizing both aerosol and surface properties with minimal assumptions. In order to employ this concept in AERONET processing, the standard inverse algorithm Dubovik, O. & King, M. D. (2000), A flexible inversion algorithm for retrieval of aerosol optical properties from sun and sky radiance measurements.
Journal of Geophysical Research, 105, 20673–20696 has been modified to retrieve surface reflectance in addition to aerosol parameters when co-incident satellite measurements are available.
The method was applied to observations of smoke and desert dust over the Mongu (Zambia) and Solar Village (Saudi Arabia) AERONET sites respectively. The AERONET data were complemented by available observations from the MISR, MODIS, and POLDER-2 satellite sensors. The retrieved bidirectional reflectance factor (BRF) and surface albedo comparison shows good agreement between results obtained using observations from different satellites.
The robustness of the method is tested by analyzing surface albedo time series retrieved during periods of high aerosol optical depth variability and low seasonal changes in surface reflectance. The analysis shows that the performance of retrieval algorithm is stable under different aerosol loadings. It is shown that much of the observed surface albedo temporal variability could be attributed mostly to the combined uncertainty in satellite radiometric calibration and aerosol vertical distribution for Mongu and to differences in satellite angular sampling on different days for Solar Village.
The sensitivity of surface retrievals to assumptions on aerosol vertical distribution and aerosol particle shape are analyzed. It is found that the maximum error in retrieved surface albedo at 0.44 μm is 0.035 for aerosol optical depth 0.85 at 0.44 μm. For aerosol optical depths lower than ∼
0.7 the error in retrieved surface albedo is less than 0.02. Analysis of particle shape assumptions on surface retrievals showed that aerosol particle non-sphericity significantly affects the angular shape of BRF, but not the surface albedo.
Finally, the sensitivity of AERONET aerosol retrievals to uncertainty in assumed surface reflectance is analyzed by comparing aerosol retrievals obtained with different surface assumptions. It is found that the uncertainty in surface reflectance model employed in the version 1 AERONET operational algorithm is larger than was previously assumed in Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., & Slutsker, I. (2000), Accuracy assessment of aerosol optical properties retrieved from AERONET sun and sky radiance measurements.
Journal of Geophysical Research, 105, 9791–9806 and may have more significant effect on the retrieved aerosol properties than was documented in that work. In particular, larger errors were encountered for the real part of the refractive index (∼
0.05–0.07 increase) and maximum of the particle size distribution (∼
20% decrease) retrievals for the Mongu case, when the aerosol optical depth was relatively small (∼
0.4 at 0.44 μm). The retrieved single scattering albedo uncertainties were within the error bars (0.03) estimated in Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., & Slutsker, I. (2000), Accuracy assessment of aerosol optical properties retrieved from AERONET sun and sky radiance measurements.
Journal of Geophysical Research, 105, 9791–9806, with the exception of the 0.44 μm retrievals for the desert dust case when they increased by ∼
0.09 and 0.07 for low and high aerosol loadings respectively.
The Ocean and Land Color Instrument (OLCI) onboard the Copernicus Sentinel-3A satellite is a medium-resolution and multi-spectral push-broom imager acquiring radiance in 21 spectral bands covering ...from the visible to the far near-infrared. These measurements are primary dedicated to land & ocean color applications, but actually include also reliable information for atmospheric aerosol and surface brightness characterization. In the framework of the EUMETSAT funded study to support the Copernicus Program, we describe the retrieval of aerosol and surface properties from OLCI single-viewing multi-spectral Top-Of-Atmosphere (TOA) radiances based on the Generalized Retrieval of Atmosphere and Surface Properties (GRASP) algorithm. The high potential of the OLCI/GRASP configuration stems from the attempt to retrieve both aerosol load and surface reflectance simultaneously using a globally consistent high-level approach. For example, both over land and ocean surfaces OLCI/GRASP uses 9 spectral channels (albeit with different weights), strictly the same prescribed aerosol models and globally the same a priori constraints (though with some differences for observations over land and ocean). Due to the lack of angular multi-viewing information, the directional properties of underlying surface are strongly constrained in the retrieval: over ocean the Fresnel reflection together with foam/whitecap albedo are exclusively computed using a priori wind speed; over land, the Bidirectional Reflectance Distribution Function (BRDF) is slightly adjusted from a priori values of climatological Ross-Li volumetric and geometric terms. Meanwhile, the isotropic reflectance is retrieved globally under mild spectral smoothness constraints. It should be noticed that OLCI/GRASP configuration employs innovative multi-pixel concept (Dubovik et al., 2011) that enhance retrieval by simultaneously inverting large group of pixels. The concept allows for benefiting from knowledge about natural variability of the retrieved parameters.
The obtained OLCI/GRASP products were validated with the Aerosol Robotic Network (AERONET) and Maritime Aerosol Network (MAN) and intercompared with the Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol and surface products. The overall performance is quite comparable to the community-referenced MODIS. Over ocean the OLCI/GRASP results are encouraging with 67% of the AOD (550 nm) satisfying the Global Climate Observing System (GCOS) requirement using AERONET coastal sites and 74% using MAN deep ocean measurements, and an AOD (550 nm) bias 0.01 with AERONET and nearly zero bias with MAN. Over land, 48% of OLCI/GRASP AOD (550 nm) satisfy the GCOS requirement and a bias within ±0.01 for total and AOD < 0.2. Key challenges are identified and discussed: adequate screening of cloud contaminations, retrieval of aerosol over bright surfaces and in the regions containing complex mixtures of aerosol.
•We retrieve aerosol and surface from OLCI/Sentinel-3A based on GRASP algorithm.•The directional BRDF are constrained using wind speed and climatology.•OLCI/GRASP products are validated with AERONET, MAN and intercompared with MODIS.
Cosmic dust is a polydisperse mixture of irregular, often aggregated, particles. Previous attempts have tried to simulate polarimetric properties of this dust using aggregate dust models, but it has ...not been possible to consider particle sizes larger than a couple of microns due to limitations of computer memory and processing power. Attempts have also been made to replace aggregates by polydisperse regular particles (spheres, spheroids, cylinders), but those models could not consistently reproduce the observed photopolarimetric characteristics. In this study, we introduce to the astronomical community the software package developed by Dubovik et al. (2006) for modeling light scattering by a polydisperse mixture of randomly oriented smooth and rough spheroids of a variety of aspect ratios. The roughness of spheroids is defined by a normal distribution of the surface slopes, and its degree depends on the standard deviation of the distribution (which is zero for smooth surface and greater than zero for rough surface). The pre-calculated kernels in the software package allow for fast, accurate, and flexible modeling of different size and shape distributions. We present our results of a systematic investigation of polarization obtained with the rough and smooth spheroid models; we study differences in their phase angle dependence and how those differences change with the particle size distribution. We found that the difference between smooth and rough particles increases with increasing effective size parameter and affects mainly the value and position of the maximum polarization. Negative polarization was found to be typical only for silicate-like refractive indexes and only when the particles have size parameters within 2.5–25. As an example of an application of the rough spheroid model, we made computations for rough spheroids that have a size distribution and composition typical for cometary dust. We found that a mixture of porous rough spheroids made of absorbing material compositionally similar to comet Halley’s dust and solid silicate spheroids, dominated by particles of size parameter 5<x<20, can reproduce angular and spectral characteristics of the brightness and polarization observed for cometary dust.
•A package to compute light scattering by polydisperse rough spheroids is described.•The package is used to study polarimetric differences between rough and smooth spheroids.•We model polarization for materials typical for cosmic dust, including porous particles.•A model, consistent with photopolarimetric observations of comet dust, is presented.