Articles | Volume 19, issue 23
Atmos. Chem. Phys., 19, 14585–14606, 2019
Atmos. Chem. Phys., 19, 14585–14606, 2019

Research article 03 Dec 2019

Research article | 03 Dec 2019

Constraining global aerosol emissions using POLDER/PARASOL satellite remote sensing observations

Constraining global aerosol emissions using POLDER/PARASOL satellite remote sensing observations
Cheng Chen1,5, Oleg Dubovik1, Daven K. Henze2, Mian Chin3, Tatyana Lapyonok1, Gregory L. Schuster4, Fabrice Ducos1, David Fuertes5, Pavel Litvinov5, Lei Li1,6, Anton Lopatin5, Qiaoyun Hu1, and Benjamin Torres1 Cheng Chen et al.
  • 1Laboratoire d'Optique Atmosphérique (LOA), UMR8518 CNRS, Université de Lille, Villeneuve-d'ascq, 59655, France
  • 2Department of Mechanical Engineering, University of Colorado, Boulder, Colorado 80309, USA
  • 3NASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA
  • 4NASA Langley Research Center, Hampton, Virginia 23681, USA
  • 5GRASP-SAS, Remote Sensing Developments, Université de Lille, Villeneuve-d'ascq, 59655, France
  • 6State Key Laboratory of Severe Weather (LASW) and Institute of Atmospheric Composition, Chinese Academy of Meteorological Sciences, CMA, Beijing, 100081, China

Correspondence: Cheng Chen ( and Oleg Dubovik (


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−1 for BC, 109.9 Tg yr−1 for OC and 731.6 Tg yr−1 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.

1 Introduction

Atmospheric aerosol emission inventories are often used to drive chemical transport model (CTM) simulations of aerosol distributions on regional and global scales (Boucher, 2015; Brasseur and Jacob, 2017; Granier et al., 2011). Satellite-retrieved columnar aerosol optical depth (AOD) is directly related to light extinction due to the presence of aerosols; hence, satellite-retrieved columnar AOD is widely used to evaluate the spatial and temporal variability of aerosols simulated from CTMs (e.g., Chin et al., 2002; Ginoux et al., 2006; Kinne et al., 2006, 2003; Liu et al., 2012; Ocko and Ginoux, 2017; Pozzer et al., 2015; Schulz et al., 2006; Tegen et al., 2019). A general agreement has been shown for columnar AOD between model simulations and satellite observations in the Aerosol Comparisons between Observations and Models (AeroCom) “Experiment A” multimodel assessments (Kinne et al., 2006). However, this study also revealed a large model diversity of species-specific AOD and aerosol absorption optical depth (AAOD), which encourages research to harmonize and improve the emissions of individual aerosol species and aerosol precursors and the representation of aerosol absorption and other elements (Kinne et al., 2006; Samset et al., 2018). Accurate knowledge of spatial and temporal distribution of species-specific aerosol emissions is also useful for numerical weather prediction (NWP) (Benedetti et al., 2018; Xian et al., 2019).

There have been several efforts to improve aerosol emission inventories by using satellite observations and inverse modeling (Dubovik et al., 2008; Huneeus et al., 2012, 2013; Wang et al., 2012; Xu et al., 2013; Zhang et al., 2015, 2005; Escribano et al., 2016, 2017). However, most of them are regional studies that focus on a single aerosol or gas species, or they use predefined regions to reduce the size of state vector (Huneeus et al., 2012; Zhang et al., 2015, 2005). There are only a few studies that use detailed satellite information to simultaneously retrieve emissions of multiple aerosol components at the native spatial resolution of a forward CTM model. For example, we have previously developed an inverse modeling framework for retrieving aerosol emissions of black carbon (BC), organic carbon (OC) and desert dust (DD) components within the GEOS-Chem model (Chen et al., 2018). Specifically, the emissions of BC, OC and DD were simultaneously derived from satellite retrievals of spectral AOD and AAOD that were provided by the PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) products. Chen et al. (2018) successfully used this method at the regional scale over all of Africa and the Arabian Peninsula.

This work expands the inversion algorithm of Chen et al. (2018) to the global scale. We also refine an assumption that defines the observation error covariance matrix for the recently released PARASOL Level 3 AOD and AAOD aerosol products generated by the GRASP (Generalized Retrieval of Atmosphere and Surface Properties) algorithm. The method is applied to derive global BC, OC and DD aerosol emissions for the year 2010, using the updated Hemispheric Transport of Air Pollution (HTAP) Phase 2 emission data (HTAP, 2010; Janssens-Maenhout et al., 2015) for the initial estimate of anthropogenic emissions. Then the satellite-derived global a posteriori aerosol emissions are intensively evaluated in comparison to the a priori emission inventories and other independent measurements.

Figure 1General concept of satellite remote sensing of global aerosol emissions.

2 Methods and data

Figure 1 demonstrates the general concept of satellite remote sensing of the global distribution and strength of aerosol emissions. A priori emissions are used with the GEOS-Chem model to create a simulated aerosol loading, which are then compared to observed spectral AOD and AAOD from PARASOL. The inverse modeling based on adjoint GEOS-Chem iteratively optimizes a priori emissions to minimize the differences between observed and modeled AOD and AAOD. A posteriori emissions are equivalent to retrieved or optimized emissions, which are then input to the GEOS-Chem model simulation, and the results are verified with independent AERONET, MODIS and OMI aerosol products.

2.1 POLDER/PARASOL aerosol dataset generated by the GRASP algorithm

The GRASP algorithm implements statistically optimized fitting of diverse observations using the multi-term LSM (least square method) (Dubovik, 2004). The basic concept of this approach was introduced and implemented in the AERONET algorithm developed for aerosol characterization from ground-based radiometric observations (Dubovik and King, 2000; Dubovik et al., 2000; King and Dubovik, 2013). Dubovik et al. (2011, 2014) have adapted and extended this concept in GRASP algorithm, which is designed to retrieve aerosol and surface properties from satellite and other observations. As a new inversion development, a multipixel retrieval concept was implemented in the GRASP (Dubovik et al., 2011). Using this concept, the satellite retrieval is implemented as a statistically optimal simultaneous fitting of observations over a large number of pixels. This approach allows for improved accuracy of the retrieval by applying known a priori constraints on temporal or/and spatial variability of the derived parameters. In addition, GRASP is a highly versatile algorithm that has been applied for a large variety of different types of satellite, ground-based and airborne remote-sensing measurements by photometers, lidars, satellite sensors, nephelometers, sky cameras, etc. (Benavent-Oltra et al., 2017; Dubovik et al., 2019; Espinosa et al., 2017; Hu et al., 2019; Li et al., 2019; Lopatin et al., 2013; Román et al., 2017, 2018; Torres et al., 2017; Tsekeri et al., 2017).

The spaceborne multidirectional, multispectral polarized POLDER-3 (Polarization and Directionality of the Earth's Reflectances) imager on board PARASOL can measure the global angular distribution of intensity and polarization of solar radiation reflected to space by the earth–atmosphere system (Deschamps et al., 1994; Deuzé et al., 1999, 2001; Tanré et al., 2011). Throughout this article, we use “PARASOL” to denote POLDER/PARASOL observations. The GRASP algorithm inverts PARASOL comprehensive measurements to derive aerosol properties (e.g., extinction, absorption, size and composition) and surface BRDF (bidirectional reflectance distribution function) and BPDF (bidirectional polarization distribution function) properties. The development of the GRASP algorithm is described in Dubovik et al. (2011, 2014), and a description of some PARASOL/GRASP aerosol products can be found in Chen et al. (2018), Kokhanovsky et al. (2015), Popp et al. (2016) and Sayer et al. (2018). The accuracy of the GRASP algorithm configured for AERONET measurements has been evaluated with a laboratory experiment by Schuster et al. (2019).

We use PARASOL/GRASP Level 3 AOD and AAOD at six wavelengths (443, 490, 565, 670, 865 and 1020 nm) in this study. The Level 3 PARASOL/GRASP aerosol products are rescaled in 1×1 spatial resolution, and have been archived at the AERIS/ICARE Data and Services Center (, last access: 15 March 2019). Information about PARASOL/GRASP Level 3 products is also available from the GRASP-OPEN website (, last access: 15 March 2019). In order to colocate with model 2 (latitude) ×2.5 (longitude) spatial resolution the PARASOL/GRASP aerosol data are aggregated into 2×2.5 model grid boxes. Any grid box with less than two PARASOL/GRASP retrievals is omitted. The global distribution of PARASOL/GRASP 2×2.5 spectral AOD and AAOD at 443, 490, 565, 670, 865 and 1020 nm used in the inversion is presented in Fig. S1 in the Supplement.

2.2 GEOS-Chem inverse modeling framework

The GEOS-Chem chemical transport model simulates the spatial and temporal mass distribution of each aerosol species by modeling transport processes (e.g., advection, convection, diffusion, deposition) and source injection (Bey et al., 2001; Brasseur and Jacob, 2017; Jacob, 1999). Once the global distribution of aerosol mass is known, it is generally converted to distributions of aerosol extinction (AOD) and absorption (AAOD) by modeling the aerosol microphysical and optical properties (Martin et al., 2003). For this study, we simulate five major aerosol components, including sulfate (SU), BC, OC, DD (seven bins, with effective radii of re=0.14, 0.24, 0.45, 0.80, 1.40, 2.40 and 4.50 µm), and sea salt (SS; re=0.80µm for the accumulation mode and re=5.73µm for the coarse mode). The microphysics of each aerosol species used in the simulation are given in Chen et al. (2018). Previously, we adopted two prevalent assumptions of BC refractive indices (Case 1: m=1.75-0.45i; Case 2: m=1.95-0.79i). As suggested by the sensitivity test by Chen et al. (2018), using these two different assumptions leads to an additional factor of ∼2.0 differences for total retrieved BC emissions because of the differences in the mass absorption efficiency (βa; Case 1: βa=4.5 m2 g−1; Case 2: βa=6.3 m2 g−1 at 565 nm). However, the recommended value for uncoated BC at 550 nm is 7.5±1.2 m2 g−1 (Bond and Bergstrom, 2006). Hence, in this study, we retrieve BC emissions using βa from the Case 2 assumption. The details of microphysical properties of BC, OC, DD, SU and SS used in the inversion are described in Chen et al. (2018). The assumption of external mixing of spherical particles is adopted in our inversion, as it is commonly done in most CTMs. It should be noted, however, that the particle morphologies and mixing states could have strong affects on scattering and absorption properties, thus affecting mass to optical conversion (Liu and Mishchenko, 2018). For example, the “lensing effect” of less absorbing components coated on BC could amplify total aerosol absorption (Lesins et al., 2002). The absorption enhancement due to coating is estimated to be ∼1.5 (Bond and Bergstrom, 2006). A recent study by Curci et al. (2019) implemented partial internal mixing for regional simulation and found it could improve simulation of total absorption while the spectral dependence cannot be well reproduced. Therefore, in this approach, as well as generally in CTMs, there is some intrinsic ambiguity in assumptions influencing efficiency of scattering and absorption of aerosol particles. This ambiguity is certainly among the major factors affecting accuracy of derived emissions in the current approach.

The inverse modeling framework used here was originally presented by Chen et al. (2018), wherein BC, OC and DD aerosol emissions were estimated simultaneously using spectral AOD and AAOD observations. The framework uses an adjoint of the GEOS-Chem model that was developed by Henze et al. (2007, 2009) and Wang et al. (2012). Here, SU and SS emissions are kept fixed (similar to Chen et al., 2018); in future studies, we plan to retrieve SU and SS emissions together with other emissions by using additional spatial and temporal (smoothness) constraints (Dubovik et al., 2008).

Our inverse modeling method iteratively seeks adjustments of aerosol emissions that can minimize the cost function J(S)

(1) J S = 1 2 H S - f obs T C obs - 1 H S - f obs + 1 2 γ r S - S a T C a - 1 S - S a .

In Eq. (1), the forward model H is a CTM (e.g., GEOS-Chem). fobs is the vector of observed parameters used for inversion, Cobs is the error covariance matrix of fobs. The vector S describes the 4-D distribution of emissions, and Sa is the a priori estimates of emissions. Ca is the error covariance matrix of Sa. γr is a regularization parameter.

The inversion is initialized using a priori model emissions plus a spatially uniform value of 10−4 for DD, 10−6 for BC and 5×10-6 for OC over land grid boxes where Sa=0; this allows for the detection of new sources and performs satisfactorily even at locations where the a priori knowledge of aerosol emission is poor (Chen et al., 2018). In this study, the a priori DD emission is based upon the mineral dust entrainment and deposition (DEAD) scheme (Zender et al., 2003) and the GOCART dust source function (Ginoux et al., 2001) that was implemented in GEOS-Chem model by Fairlie et al. (2007). We adopted the improved fine mode dust contribution scheme by Zhang et al. (2013) for both inversion and simulation. We used anthropogenic emissions from the Hemispheric Transport of Air Pollution (HTAP) Phase 2, with biomass-burning emissions from version 4s of the Global Fire Emissions Database (GFED; Randerson et al., 2012; van der Werf et al., 2017). The GEOS-Chem sulfate module was developed by Park et al. (2004) and carbonaceous aerosol simulations by Park et al. (2003), both of which are based upon the GOCART model scheme (Chin et al., 2002). The sea salt simulation is based upon Jaeglé et al. (2011).

The inversion system derives daily total BC, OC and DD aerosol emissions for each grid box. The daily ratio between biomass burning and anthropogenic contribution for BC and OC and the proportion of DD bins for each grid box are kept as the a priori GEOS-Chem assumption. Distinguishing the anthropogenic contribution from total emission is crucial for climate effects evaluation. Here, we propose a simple method to estimate daily anthropogenic BC (EBC_INV_AN) and OC (EOC_INV_AN) emissions from our retrieved total emission (EBC_INV) by using daily proportion of anthropogenic emissions over each grid box from a priori emission database:

(2) E BC _ INV _ AN ( x , y , t ) = E BC _ AN ( x , y , t ) E BC _ AN ( x , y , t ) + E BC _ BB ( x , y , t ) E BC _ INV ( x , y , t ) , E OC _ INV _ AN ( x , y , t ) = E OC _ AN ( x , y , t ) E OC _ AN ( x , y , t ) + E OC _ BB ( x , y , t ) E OC _ INV ( x , y , t ) ,

where EBC_INV(x,y,t) and EOC_INV(x,y,t) represent retrieved total BC and OC emission, respectively. EBC_INV_AN(x,y,t) and EOC_INV_AN(x,y,t) are derived anthropogenic BC and OC emissions from retrieved total BC and OC emission database. EBC_AN(x,y,t) and EOC_AN(x,y,t) represent anthropogenic BC and OC emissions from HTAP v2 database, and EBC_BB(x,y,t) and EOC_BB(x,y,t) are BC and OC emitted from biomass burning adapted from GFED v4s database. x,y and t indicate indexes of longitude, latitude and time, respectively.

2.3 Error covariance matrix

The error covariance matrix Cobs of observations fobs needs to be prescribed in our inversion system. In our previous study, we tested and adopted a spectral weighting schemeto fit the absolute value of AOD and AAOD at different wavelengths (Chen et al., 2018). Specifically, PARASOL/GRASP AOD and AAOD were rigorously evaluated against the AERONET dataset. The current understanding of the accuracy of PARASOL/GRASP products is ∼0.05 for absolute AOD. The AAOD accuracy is unknown but expected to be positively correlated with the AOD value; therefore, the AAOD accuracy is assumed to be 0.05∕AOD. In this work, we assume the covariance error matrices (Cobs and Ca) are diagonal, and we use these accuracy estimates for AOD and AAOD as the diagonal terms of the observation covariance matrix Cobs.

The error covariance matrix Ca for the a priori emission dataset Sa is not accurately known. As a result, this matrix is often defined using rather simple strategies that are based upon uncertainty estimates found in the literature (e.g., Escribano et al., 2016, 2017; Huneeus et al., 2012). In our previous study, we adopted a very small regularization parameter γr (e.g., 1.0×10-4) to force the inversion to rely upon the observations (Chen et al., 2018). Additionally, our previous simulations using the a posteriori emissions have shown good agreement with independent measurements, so we decided to use the same strategy for our global inversions in this study.

3 Results

3.1 Emissions

We applied our method to retrieve BC, OC and DD daily emissions on a global basis using PARASOL spectral AODs and AAODs for the year 2010. The retrieved emissions dataset is publicly available at the LOA website (, last access: 3 April 2019). Table 1 summarizes the annual a priori and a posteriori emissions of five aerosol components (DD, BC, OC, SU and SS). The a posteriori DD, BC and OC emissions are retrieved from PARASOL spectral AOD and AAOD, while the a posteriori SS emission and the SU sources (from SO2 oxidation) remain the same as the a priori simulation. We break the 1-year global inversion into 12 months, and the maximum of iterations is set to 40 for each individual retrieval run.

The retrieved aerosol emissions are 18.4 Tg yr−1 for BC, 109.9 Tg yr−1 for OC and 731.6 Tg yr−1 for DD in the year 2010. These new emissions indicate an increase of 166.7 % for BC and 184.0 % for OC, while a decrease of 42.4 % for DD with respect to the a priori GEOS-Chem emission database is seen. Table 2 compares the retrieved annual anthropogenic and biomass-burning BC and OC emissions with a priori emission database in GEOS-Chem. The method used to separate anthropogenic from total emission is described in Sect. 2.2. The retrieved anthropogenic emissions are 14.8 Tg yr−1 for BC and 85.6 Tg yr−1 for OC, representing an increase of 217.3 % for BC and 357.8 % for OC. Meanwhile, the retrieved biomass-burning emissions of BC and OC are 3.6 and 24.3 Tg yr−1, corresponding to increases of 56.5 % and 21.5 %, respectively, with respect to the a priori emission database. The comparison of the spatial distribution of anthropogenic and biomass-burning emissions of BC and OC is presented in Figs. S5 and S6 in the Supplement.

Table 1Total source strengths (unit: Tg yr−1) of five major aerosol components in the year 2010.

a Estimated based upon a priori GEOS-Chem model emission database. b Retrieved emission database that is used for a posteriori GEOS-Chem model simulation. c The SO2 and SS emissions are kept the same as a priori GEOS-Chem emission database in this study. d Particle radius ranging from 0.1 to 6.0 µm. e Unit: Tg S yr−1

Download Print Version | Download XLSX

Table 2Anthropogenic (AN) and biomass-burning (BB) emissions (unit: Tg yr−1) of BC and OC in the year 2010.

Download Print Version | Download XLSX

Spatial distributions of the a priori and a posteriori annual emissions of BC, OC and DD, and their differences (a posteriori minus a priori emissions), are shown in Fig. 2a, b and c for the year 2010. The a posteriori BC emissions are generally greater than the a priori BC emissions throughout the globe; the a posteriori increases are particularly significant over certain regions, such as Southeast Asia and central and northwest China. However, there are also notable decreases of BC a posteriori emissions (with respect to a priori emissions) observed in several grid boxes over South America and north China. The largest increase in OC emissions is in southern Africa, where biomass burning is the predominant source. Consistent with BC, there is also a slight decrease in OC emissions over the high-emission grid boxes in South America. In contrast, the a posteriori DD emissions are reduced throughout the global desert regions. As a reference, the seasonal cycle of the a posteriori BC, OC and DD emissions are presented in the Supplement illustrations in Figs. S2–S4.

Notably, the a posteriori BC and OC emission distributions are more homogenous than the a priori emission inventories. This phenomenon is probably because the emissions are reported for 2×2.5 grid boxes, and this is too coarse to characterize cities with high anthropogenic activities. The spread of a posteriori BC and OC emission distributions leads to more grid boxes being influenced by absorbing aerosols over India and China, which can be supported by evaluation with AERONET retrieval of AAOD (see Sect. 3.3.1). Briefly, there is a set of AERONET data of moderate-to-high aerosol absorption (AAOD >0.07 at 550 nm) where the a priori simulation was very close to zero, and the a posteriori simulation adjusted them (see Sect. 3.3.1).

Figure 2Global distribution of emissions for 2010 for BC (left column), OC (middle column) and DD (right column) based on (a) a priori, (b) a posteriori emission datasets, and the differences between a posteriori and a priori emissions (c).

Figure 3Comparison of monthly emissions for a priori and a posteriori datasets: (a) BC emissions, (b) OC emissions and (c) DD emissions.


Comparisons of monthly global total BC, OC and DD emissions between a priori and a posteriori emission inventories for 2010 are shown in Fig. 3. Both a priori and a posteriori BC and OC emission inventories show a maximum in August and September, while the second peak in March observed in the a priori database shift to April and May in the a posteriori database. However, a posteriori BC and OC emissions are higher than the a priori emissions throughout the year. The a posteriori to a priori ratio for monthly BC emissions is up to 4.4 in April and down to 1.1 in July. Meanwhile, the a posteriori to a priori ratio for monthly OC emissions is up to 4.6 in October and down to 1.4 in July. In contrast, the a posteriori DD emissions capture seasonal variations that are similar to the a priori DD emissions. The a posteriori DD emissions are reduced consistently throughout the year. The posteriori and priori monthly DD emissions ratio slightly varies between 0.51 (December) and 0.78 (June).

3.2 Simulation of AOD and AAOD

In a numerical modeling experiment for the year 2010, we used the a priori and a posteriori emissions inventories as inputs to the GEOS-Chem model. Then we compared the simulated AODs and AAODs that were generated from the a priori emissions inventory to the AODs and AAODs that were likewise generated from the a posteriori emissions inventory. The results are shown in Figs. 4 and 5 for AOD and AAOD at 550 nm.

Both simulations with the a priori (Fig. 4a) and a posteriori (Fig. 4b) emissions show that high values of AOD appear over the Sahara in north Africa and the Taklimakan and Gobi deserts in Asia, which are strong dust source regions. Figure 4 also indicates high values of AOD in East Asia, where anthropogenic aerosols (e.g., BC, OC and SU) are the predominant components.

One of the major differences (Fig. 4c) between a priori and a posteriori AOD is that the a posteriori AOD over desert regions is reduced, especially over the Sahara; meanwhile, the a posteriori AOD increases over industrial and biomass-burning regions. The strongest increase of a posteriori AOD occurs in southern Africa, where it is associated with biomass-burning emissions.

Figure 5 shows the comparison of global distribution of a priori (Fig. 5a) and a posteriori (Fig. 5b) AAOD in 2010. Figure 5c clearly reveals that the a posteriori AAOD is higher than the a priori AAOD throughout the globe, except over the Sahara. The differences are high over southern Africa, India, Southeast Asia and central China, where they are associated with biomass-burning and anthropogenic emissions.

The statistics of global annual mean AOD and AAOD at 550 nm for the five major aerosol components using a priori and a posteriori simulations are shown in Table 2. The a priori GEOS-Chem global mean AOD at 550 nm is 0.105, while the a posteriori simulation showed a slightly increased value of 0.119. The dust AOD decreases from 0.031 to 0.019 and its relative contribution to total AOD decreases from 29.9 % to 16.1 %, owing to the reduction of global DD emission from 1269.4 Tg yr−1 (a priori) to 731.6 Tg yr−1 (a posteriori). On the other hand, the a posteriori simulation indicates significant increase in carbonaceous AOD from a priori 0.014 (BC: 0.003; OC: 0.011) to a posteriori 0.040 (BC: 0.008; OC: 0.032). The resulting OC emissions increase from 38.7 to 109.9 Tg yr−1, and OC becomes the largest contributor to the total AOD. With the a posteriori simulation, OC accounts for 26.8 % of the total AOD, whereas the a priori simulation indicated DD as the largest contributor to total AOD (at 29.9 %).

In general, the global mean AAOD at 550 nm significantly increases from 0.0039 (a priori) to 0.0071 (a posteriori), i.e., by a factor of 1.8. In particular, DD and BC are the two major components to the total aerosol absorption that collectively account for 90.1 % (a priori) and 88.2 % (a posteriori) of total AAOD. The a posteriori DD AAOD decreases by 42.8 % relative to the a priori simulation, from 0.0014 to 0.0008. The BC AAOD increases from 0.0020 (a priori) to 0.0054 (a posteriori), a factor of ∼2.7, and the relative contribution to total AAOD increases from 52.7 % to 76.9 %.

In comparison to the AeroCom Phase II assessments of global AOD and AAOD (Myhre et al., 2013), our GEOS-Chem a posteriori global mean AOD is ∼8 % smaller. However, our a posteriori GEOS-Chem global mean AAOD (0.0071) is ∼69 % larger than the AeroCom multimodel mean of 0.0042±0.0019 (1 SD).

Figure 4Comparison of GEOS-Chem simulation of global aerosol optical depth in the year 2010 at 550 nm based upon a priori and a posteriori emission datasets: (a) a priori AOD, (b) a posteriori AOD and (c) a posteriori minus a priori AOD.

Figure 5Same as Fig. 4 but for AAOD.

3.3 Evaluation with independent measurements

To achieve a more robust evaluation, the a priori and a posteriori aerosol properties simulated with the GEOS-Chem model are evaluated with other independent measurements that are not used in our emission inversion. The definition of the statistics used in the comparison, including correlation coefficient (R), root mean square error (RMSE), mean absolute error (MAE), normalized mean bias (NMB), and mean bias (MB), can be found below:


where M represents the model results; O represents the observations; M and O are the mean values for the simulations and observations, respectively; and N is the number of data points.

3.3.1 Comparison with AERONET measurements

The Aerosol Robotic Network (AERONET; Holben et al., 1998) has provided comprehensive and accurate aerosol data from a worldwide ground-based sun photometer network for more than 2 decades. The data products include measurements of multiple-wavelength AODs and Ångström exponents (AExp). Additionally, the AERONET products also include retrievals of AAOD, single scattering albedo (SSA), absorption Ångström exponent (AAExp), size distribution, and complex refractive index (Dubovik and King, 2000, Dubovik et al., 2000, 2002, 2006). The AERONET aerosol dataset has been commonly used for satellite product validation and model evaluation for a wide range of aerosol research topics. In this section, the AERONET version 2 daily level 2.0 (Smirnov et al., 2000) AOD, AAOD, SSA at 550 nm, AExp (440–870 nm), and AAExp (440–870 nm) products are used to evaluate the a priori and a posteriori GEOS-Chem model aerosol simulation. We convert spectral AOD and AAOD to SSA, AExp and AAExp using the following equations:


We evaluate the GEOS-Chem model-simulated daily aerosol properties with the AERONET dataset. Globally, 282 AERONET sites that have available data for the year 2010 are all taken into account. We should note that AERONET aerosol datasets are averaged of all available data during the daytime (lunar photometer aerosol products will probably be available in the future; e.g., Barreto et al., 2016, and Berkoff et al., 2011) and that the GEOS-Chem daily aerosol properties are averaged based upon daytime and nighttime simulations. Additionally, the use of AERONET data comparing against corresponding 2×2.5 grid box model simulation could bring some sampling issues that cannot be ignored (Schutgens et al., 2016).

Table 3Global mean AOD and AAOD at 550 nm of five major aerosol components. Numbers in parenthesis are relative contributions of different components to AOD and AAOD at 550 nm. Statistics are based upon a priori and a posteriori GEOS-Chem simulations for the year 2010.

Download Print Version | Download XLSX

Figure 6Comparison of a priori and a posteriori GEOS-Chem simulated AOD, AAOD, SSA at 550 nm, AExp (440–870 nm), and AAExp (440–870 nm) with AERONET Level 2 daily aerosol products. The data for density plots are all aggregated into 100 bins for both x and y directions spanning from minimum to maximum values on the axes. The correlation coefficient (R), mean absolute error (MAE), normalized mean bias (NMB), and mean bias (MB) are also provided in the panels.


Figure 6 shows the comparison of a priori (red) and a posteriori (blue) GEOS-Chem model simulated AOD, AAOD, SSA, AExp and AAExp with the AERONET dataset. The a posteriori GEOS-Chem simulation has a much better agreement with AERONET data for AAOD. The correlation coefficient (R) of a priori GEOS-Chem AOD with AERONET equals 0.59, mean absolute error (MAE) is 0.13, and normalized mean bias (NMB) is 8.15 %. The a posteriori simulation improves R to 0.66 and MAE to 0.12. The NMB from a posteriori simulation (21.49 %) is higher than a priori simulation, which indicates the a posteriori GEOS-Chem AOD is generally ∼20 % higher than the AERONET dataset. This is possibly due to the use of AERONET point measurements to evaluate model simulation at the 2×2.5 grid box. The statistics for simulated AAOD show significant improvement from the a priori (R=0.05, MAE =0.12, NMB =-78.45 %, MB =-0.012) to the a posteriori (R=0.49, MAE =0.02, NMB =0.71 %, MB =0.000) simulations. The a priori GEOS-Chem simulated AAOD is lower than AERONET data (NMB =-78.45 %, MB =-0.012), while the a posteriori AAOD is much closer to AERONET (NMB =-0.71 %, MB =0.000).

The inversion framework derives aerosol emissions by fitting spectral AOD and AAOD from PARASOL. The SSA, AExp and AAExp are the derived products from spectral AOD and AAOD (Eqs. 9–11). The GEOS-Chem model SSA correlations with AERONET improve slightly from R=0.34 for the a priori to R=0.47 for the a posteriori simulation. The a priori and a posteriori simulations of AExp show similar performance when evaluating with AERONET data. However, the a posteriori R=0.70 for AExp is slightly worse than the a priori (R=0.74), which is likely related to the ∼20 % higher bias of the a posteriori AOD. Consistent with AAOD, the GEOS-Chem AAExp shows improvements from the a priori simulation (R=0.01, MAE =0.71, NMB =6.84 %, MB =0.15) to the a posteriori simulation (R=0.62, MAE =0.35, NMB =-7.35 %, MB =-0.11), which indicates a better representation of the spectral dependence of aerosol absorption (Russell et al., 2010; Schuster et al., 2016).

Figure 7The differences between a posteriori and a priori GEOS-Chem simulated AOD, AAOD, SSA, AExp and AAExp correlation coefficients (R) with AERONET daily aerosol products over all sites with colocated data in 2010.

To better understand the regional performance of a priori and a posteriori GEOS-Chem simulations, we conducted the comparison of daily aerosol products of AOD, AAOD, SSA, AExp and AAExp over all AERONET sites with colocated data in 2010. Figure 7 shows the differences in correlation coefficients between the a posteriori and a priori simulations. The red circles indicate sites where the a posteriori simulation has higher correlation with AERONET than the a priori simulation. Alternatively, the blue circles indicate sites where the a priori simulation shows better correlation than the a posteriori simulations. There are 202 out of a total of 282 sites that show improved correlation for AOD using the a posteriori emission data. Additionally, 176 out of 272 sites show improvement for AExp, 97 of 162 sites show improvement for SSA, 89 of 167 sites show improvement for AAOD, and 84 of 167 sites have a better correlation for AAExp. The a posteriori simulation loses correlation with AERONET for AAOD in central Eurasia and the western United States, which needs further investigation in future studies.

Intensive wild fire events over central Russia during the Summer of 2010 have been reported in several studies (e.g., Chubarova et al., 2012; Gorchakova and Mokhov, 2012; Huijnen et al., 2012; Péré et al., 2014; R'Honi et al., 2013). An increase in daily AOD is observed at the AERONET site Moscow_MSU_MO (55.707 N, 37.522 E) from late July to early August, with a maximum on 7 August 2010. The ground-based observation of this short-term fire event is suitable to evaluate the retrieved emissions used for daily simulation results. The comparison of the time series of the daily AOD and AAOD at 550 nm has been conducted for site Moscow_MSU_MO, and the results are shown in Fig. 8. The geolocations of Moscow_MSU_MO site are apparent in Fig. 9 as a red star. The a priori GEOS-Chem simulation underestimates aerosol extinction and absorption during the fire event, which may be associated with underestimation of biomass-burning emissions. The a posteriori GEOS-Chem simulation has much better agreement with observations in terms of temporal variation and intensity using retrieved daily BC and OC emissions.

Figure 8Time serial plot of AOD and AAOD at 550 nm from AERONET (blue star), PARASOL/GRASP (red circles), a priori (black line) and a posteriori (green line) GEOS-Chem simulations at the Moscow_MSU_MO (55.707 N, 37.522 E) site.


3.3.2 Comparison with MODIS aerosol products

In this section, we evaluate the agreement of a priori and a posteriori GEOS-Chem simulated AOD at 550 nm with the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 (C6) Level 3 merged aerosol products (Sayer et al., 2014). Dark Target (DT) and Deep Blue (DB) are two well known retrieval algorithms developed for processing MODIS atmospheric aerosol products (Hsu et al., 2004, 2013; Kaufman et al., 1997; Levy et al., 2013; Remer et al., 2005; Sayer et al., 2013; Tanré et al., 1997). The MODIS/Aqua merged aerosol products combing the DB with DT land/ocean data provide more gap-filled retrievals, which are suitable for model evaluation. In order to colocate with model data, the 1 level 3 products are then aggregated into model 2×2.5 grid boxes, and any grid box with less than 2 MODIS AOD is omitted.

Figure 9Comparison of a priori and a posteriori GEOS-Chem simulated seasonal AOD at 550 nm with MODIS C6 Dark Target and Deep Blue merged products.

Figure 10Differences between a priori and a posteriori GEOS-Chem seasonal AOD with MODIS C6 merged products.

In order to assess the aerosol seasonal cycle and some peaks of aerosol loading, we focus on AOD seasonal pattern from MODIS, a priori and a posteriori GEOS-Chem simulations (Fig. 9), and their differences (Fig. 10). The a priori GEOS-Chem simulation strongly overestimates aerosol loading over the Sahara throughout the year; the aerosol loading over Sahara is reduced in the a posteriori GEOS-Chem simulation, while the simulated AOD is still slightly higher than the MODIS AOD. Over India, the a priori simulation underestimates aerosol loading during December–January–February (DJF), June–July–August (JJA) and September–October–November (SON), and the a posteriori simulation improves the consistency between model and observation especially in DJF and SON. The a posteriori model simulated AOD in JJA is still lower than MODIS aerosol products over India. High aerosol loading occurred over eastern China throughout the year, which can be inferred both from observations and simulations. While the a posteriori simulated AOD over eastern China is lower than MODIS data in JJA, over biomass-burning regions (e.g., southern Africa and South America) the model simulation shows a consistent seasonal variability with MODIS data and the peaks are in JJA and SON. One of the major discrepancies between the a priori model simulation and MODIS over biomass-burning regions is that the a priori model simulated AOD is lower than observations. The biomass-burning peak over South America is observed in SON by MODIS; however, the a priori model simulated AOD in JJA is higher than that in SON. The a posteriori simulation using retrieved emissions reduces this bias. The 1–2 month delayed biomass-burning peak inferred from observations has also been reported over Africa in a recent study by Zheng et al. (2018). Overall, the a posteriori model simulated AOD shows a better agreement with independent MODIS observations over southern Africa and South America, where the aerosol is associated with biomass-burning emissions. Over central Europe, MODIS observed a high aerosol loading event in JJA (wild fires events over central Russia in summer 2010; see Fig. 8), which is not well reproduced by the a priori GEOS-Chem simulation. The retrieved emission data help to improve the a posteriori simulation reporting high aerosol loading there; however, the a posteriori AOD is still somewhat lower than MODIS.

The statistics of a priori and a posteriori GEOS-Chem simulated AOD evaluated with MODIS AOD are shown in Table 3. The evaluation was conducted at 550 nm for daily products in the year 2010. The GEOS-Chem simulation shows a better agreement with independent MODIS AOD from a priori (R=0.49, RMSE =0.19, MAE =0.054, NMB =-18.6 %) to a posteriori (R=0.59, RMSE =0.15, MAE =0.045, NMB =-13.2 %) simulations.

Table 4Statistics for evaluation of a priori and a posteriori daily GEOS-Chem simulation with MODIS AOD and OMI AAOD in the year 2010.

Download Print Version | Download XLSX

3.3.3 Comparison with OMI

In this section, we discuss the seasonal variability of AAOD from the Ozone Monitoring Instrument (OMI) near-UV algorithm (OMAERUV) and the a priori and a posteriori GEOS-Chem simulation, as well as their differences. The OMAERUV algorithm uses two UV wavelengths to derive columnar AAOD (Torres et al., 2007). The climatology of Atmospheric Infrared Sounder (AIRS) carbon monoxide (CO) observations and aerosol layer height information from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) are adopted in the latest OMAERUV algorithm (Torres et al., 2013), and the assessment of the OMI/OMAERUV aerosol products are described in Ahn et al. (2014) and Jethva et al. (2014). The OMI/OMAERUV Level 3 aerosol products with 1 spatial resolution are used for model evaluation. To aggregate to model resolution, any grid box with less than 2 OMI/OMAERUV retrievals is omitted.

Figure 11Comparison of a priori and a posteriori GEOS-Chem simulated seasonal AAOD at 500 nm with OMI/OMAERUV products.

Figure 12Differences between a priori and a posteriori GEOS-Chem seasonal AAOD with OMI/OMAERUV products.

Figure 11 shows seasonal variability of AAOD at 500 nm for OMI/OMAERUV (left column), a priori GEOS-Chem (middle column), and a posteriori GEOS-Chem (right column). The differences between a priori and a posteriori GEOS-Chem AAOD with OMI are presented in Fig. 12. Here, the GEOS-Chem 500 nm AAOD is interpolated using AAOD at two wavelengths based on absorbing Ångström exponent, αa443-865=lnτa443τa865/ln(865443). Consistent with the simulations of AOD, the a priori GEOS-Chem simulated AAOD over the Sahara is higher than OMI AAOD products throughout the year. However, the a priori AAOD is lower than OMI data in the westward grid boxes over the Atlantic Ocean, which suggests that the removal processes in the model may be too rapid (Ridley et al., 2012, 2016). The a posteriori GEOS-Chem AAOD over the Sahara source region decreases (with respect to a priori AAOD) and shows a good match with OMI AAOD; nevertheless, the discrepancy for transported dust over ocean grid boxes still exists.

The a posteriori GEOS-Chem simulation indicates moderate aerosol absorption (AAOD500nm=0.05) over central and northern China in SON, whereas the a priori model simulation and OMI data show less aerosol absorption (AAOD500nm=0.02). The a posteriori GEOS-Chem simulation shows good agreement with OMI over South America in JJA, but the a priori AAOD is higher than both a posteriori AAOD and OMI. Seasonal biomass-burning events can cause high AAOD in the OMI observations as well as the a priori and a posteriori GEOS-Chem simulations over southern Africa and South America; however, the a posteriori AAOD is higher than the a priori and OMI AAODs. Over India, OMI-observed aerosol absorption shows a peak season during March–April–May (MAM), and the a posteriori simulation successfully reproduces this seasonal peak. However, the a posteriori AAOD is slightly higher than OMI in DJF and SON over India. This is partially due to the fact that the inversion framework was fixed to SU as a priori emission database, which could result in the propagation of the uncertainty from a priori SU into a posteriori BC and OC emissions over regions with predominant SU source.

The statistics of a priori and a posteriori GEOS-Chem daily AAOD evaluated with OMI AAOD are shown in Table 3. The a priori AAOD shows a similar correlation coefficient (R=0.39) to the a posteriori AAOD (R=0.38) and similar mean absolute error (a priori: MAE =0.015; a posteriori: MAE =0.012). However, the a posteriori AAOD presents a better root mean square error (a priori: RMSE =0.026; a posteriori: RMSE =0.014). Meanwhile, the normalized mean bias is significant for both the a priori (NMB =-17.8 %) and a posteriori (NMB =19.2 %) simulations, which indicates that both GEOS-Chem AAOD simulations differ from OMI by ∼20 %.

3.3.4 Comparison with IMPROVE surface measurements of BC concentration

This section describes an evaluation that we conducted for simulated surface concentrations of BC mass. The Interagency Monitoring of Protected Visual Environments (IMPROVE) is a network of in situ aerosol measurement sites located in US national parks (Malm et al., 1994, 2003). We obtained the surface concentration measurements of BC at IMPROVE sites and evaluated the a priori and a posteriori GEOS-Chem simulation of surface BC at these sites for the year 2010. Results for the annual mean surface concentration of BC are shown in Fig. 13. Here, we see that the a priori GEOS-Chem surface BC concentration is lower than the IMPROVE data, with NMB =-26.20 %, MB =-0.051, which is consistent with previous work (Chin et al., 2014; Wang et al., 2014). The a posteriori GEOS-Chem simulation has higher values and therefore shows better agreement with the high surface BC concentrations over the eastern US than the a priori simulation. However, the a posteriori GEOS-Chem simulation overestimates BC concentrations over most of the central and western US sites where surface BC concentrations are low; this leads to an ensemble NMB =+18.29 % and MB =+0.035. In terms of correlation coefficient, the a posteriori simulation slightly improves it from 0.43 to 0.47.

Figure 13Evaluation a priori (a) and a posteriori (b) GEOS-Chem simulated annual mean surface BC concentration in the Unite States with observations (circles) from the IMPROVE network for the year 2010. The correlation coefficient (R) and normalized mean bias (NMB) are provided in the bottom right corner.

4 Discussion and conclusions

In this study, we have used PARASOL spectral AOD and AAOD generated by the GRASP algorithm to retrieve global BC, OC and DD aerosol emissions based upon the development of the GEOS-Chem inverse modeling framework. Specifically, PARASOL/GRASP AOD and AAOD at six wavelengths (443, 490, 565, 670, 865 and 1020 nm) were used to correct the aerosol emission fields using the inverse modeling framework developed by Chen et al. (2018). This resulted in improved global daily aerosol emissions at the spatial resolution of transport models.

The retrieved global annual DD emission of 731.6 Tg yr−1 is 42.4 % less than the a priori GEOS-Chem DD model value of 1269.4 Tg yr−1. The retrieved DD annual emission is also near the low end of the estimates provided by the AeroCom project (700–4000 Tg yr−1; Huneeus et al., 2011; Kinne et al., 2006; Textor et al., 2006). This is partially due to the exclusion of super coarse mode dust particles in the retrieval, as only particle radii ranging from 0.1 to 6.0 µm are taken into account. An overestimation of simulated global DD emission was also reported by Huneeus et al. (2012), who assimilated MODIS AOD into a global aerosol model. However, the overestimation of dust lifetime, poor characterization of soil data used for dust mobilization calculation and the overestimation of small dust particles (Kok et al., 2017) are also non-negligible.

We used the GFED v4s and HTAP v2 emission inventories for the a priori GEOS-Chem simulation, where the global BC emission is 6.9 Tg yr−1. The retrieved global BC emission is 18.4 Tg yr−1, which is 166.7 % higher than the a priori model inventories. Our estimation of global BC emission is close to the other top-down studies by Huneeus et al. (2012) and Cohen and Wang (2014). The former study estimates global BC emission of 15 Tg yr−1, whereas the latter study gives the emission range from 14.6 to 22.2 Tg yr−1. In comparison, the best estimate from bottom-up inventory methods is 7.5 Tg yr−1, with an uncertainty range of 2.0 to 29.0 Tg yr−1 (Bond et al., 2013). Correspondingly, the retrieved global annual OC emission flux is 109.9 Tg yr−1, which is 184.0 % higher than a priori emission inventories (38.7 Tg yr−1). Huneeus et al. (2012) estimated global organic matter (OM) emission as 119 Tg yr−1, which is equivalent to ∼85 Tg yr−1 OC based upon a conversion factor of 1.4 for OM  OC.

We introduced a method to separate anthropogenic and biomass-burning BC and OC from retrieved total BC and OC emissions by using a priori daily proportion. The retrieved anthropogenic BC emissions are 14.8 and 4.6 Tg yr−1, which is 217.3 % higher than a priori inventory adopted from HTAP v2. The retrieved biomass-burning BC emission is 4.6 Tg yr−1, showing an increase of 56.5 % with respect to a priori GFED v4s emission inventory in the GEOS-Chem model. The retrieved total OC emission is split into anthropogenic OC of 85.6 Tg yr−1 and biomass-burning OC of 24.3 Tg yr−1, which are 357.8 % and 21.5 % higher than the a priori emission inventories, respectively.

The resulting GEOS-Chem a posteriori annual mean AOD and AAOD using the retrieved emission data are 0.119 and 0.0071, respectively, at the 550 nm wavelength. These calculations indicate a decrease of 8 % for AOD and an increase of 69 % for AAOD with respect to the AeroCom Phase II multimodel assessment.

The fidelity of the results is confirmed by evaluating the a posteriori simulations of aerosol properties with independent measurements. Namely, in order to validate the retrieved emissions, the a posteriori model simulations of AOD, AAOD, SSA, AExp and AAExp were compared to independent measurements from AERONET, MODIS and OMI. We also note that the AERONET dataset is temporally more frequent than the PARASOL observations that we used to obtain the a posteriori emissions. The a posteriori GEOS-Chem daily AODs and AAODs show a better agreement (higher correlation coefficients and lower biases) with AERONET values than the a priori simulation. In addition, the a posteriori SSA, AExp and AAExp also show good agreement with AERONET data; this indicates that PARASOL provided sufficient constraints for fitting spectral AOD and AAOD, and that the retrieved emission dataset can provide reliable model simulations. Besides, a posteriori GEOS-Chem AOD and AAOD exhibit a similar seasonal pattern with the MODIS AOD and OMI AAOD, respectively, during all seasons, which indicates that the retrieved emissions are capable of capturing the major events (e.g., dust hot spots, biomass burning, anthropogenic activities). However, the a posteriori simulation overestimates AOD and AAOD over that Sahara dust source region, while underestimating AOD and AAOD over grid boxes located downwind over the Atlantic Ocean. This is probably caused by an overestimation of the retrieved DD emission over the Sahara combined with a GEOS-Chem removal process that may be too rapid (Ridley et al., 2012, 2016). During biomass-burning seasons (e.g., JJA and SON over South America and southern Africa), the a posteriori AOD shows good agreement with MODIS; meanwhile, the a posteriori AAOD is slightly higher than the OMI AAOD. This could be caused by light-absorbing OC that is not included in the simulation or the inversion. Light-absorbing OC is also known as brown carbon (BrC), and it is characterized by absorption that decreases from UV to mid-visible wavelengths (Feng et al., 2013; Lack et al., 2012). The lack of BrC in our framework could cause the retrieval to generate more BC in order to capture the observed aerosol absorption (since BC has strong absorption throughout the visible and near-infrared wavelengths; Sato et al., 2003). As a consequence, the a posteriori AAOD may be overestimated from mid-visible to near-infrared wavelengths when BrC is not included.

Our evaluation of BC surface concentrations with the IMPROVE network over the United States indicates the possibility to improve the aerosol mass simulation based on inversion of satellite-derived columnar spectral aerosol extinction and absorption. However, the a posteriori simulation shows overestimation of surface BC concentration over sites with low levels of BC, which is probably due to an overestimation of a posteriori BC emissions over low-loading regions or a modeled BC lifetime that is too long (Lund et al., 2018; Wang et al., 2014).

The result of this study is the satellite-based aerosol emission database that is commonly used as GEOS-Chem default input and was adjusted to tune the global POLDER/PARASOL observations of spectral AOD and AAOD. The analysis in the paper, as well as previous work by Chen et al. (2018), shows that the aerosol distribution modeled with the satellite-based aerosol emission improves the agreement of modeling results with the independent AERONET, MODIS, OMI and IMPROVE data. These validation results support the validity of the identified corrections of the emissions. Therefore, if this database is used for initialization of not only the GEOS-Chem model but also any other aerosol transport or GCM (global climate model), then the effects of these suggested significant corrections to the amount of mass of DD, BC and OC, and their spatiotemporal effects on the climate and environment, can be studied.

To recapitulate, we derived global BC, OC and DD aerosol emission fields in a GEOS-Chem modeling framework that was constrained with PARASOL/GRASP spectral AODs and AAODs. Our study shows that this method can be useful for improving global aerosol simulations with CTMs. In the future, we plan to use the entire PARASOL dataset to generate a satellite-based aerosol emission database; this is expected to improve multiyear aerosol simulations of AOD, AAOD, SSA, AExp and AAExp in CTMs. In addition, the efforts to better understand the aerosol life cycle at the process level (Textor et al., 2007) are essential to inversion of aerosol emission, aerosol prediction and aerosol climate effect evaluation.

Data availability

The AERONET version 2 data are available at (Holben et al., 1998). The MODIS C6 merged AOD and OMI/OMAERUV AAOD products are publicly available at NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC,, Levy et al., 2013; Sayer et al., 2014). The IMPROVE surface BC concentration data are available at (Malm et al., 1994, 2003). The retrieved daily DD, BC and OC aerosol emissions for the year 2010 are publicly available at (Chen et al., 2018, 2019).


The supplement related to this article is available online at:

Author contributions

CC, OD, DKH and TL contributed to the inversion algorithm development. DF, PL and AL prepared the aerosol dataset from POLDER/PARASOL generated by the GRASP algorithm. CC, TL and FD carried out the data processing. CC analyzed the results with contribution from OD, MC, DKH and GLS; LL, QH, PL, AL and BT participated in scientific discussions. CC and OD wrote the paper with inputs from all authors.

Competing interests

The authors declare that they have no conflict of interest.


This work is supported by the Laboratory of Excellence CaPPA – Chemical and Physical Properties of the Atmosphere – project, which is funded by the French National Research Agency (ANR) under contract ANR-ll-LABX-0005-01. We would like to thank the AERONET, ACTRIS project infrastructure, MODIS and IMPROVE teams for sharing the data used in this study.

Financial support

This research has been supported by the French National Research Agency (grant no. ANR-ll-LABX-0005-01).

Review statement

This paper was edited by Johannes Quaas and reviewed by two anonymous referees.


Ahn, C., Torres, O., and Jethva, H.: Assessment of OMI near-UV aerosol optical depth over land, J. Geophys. Res.-Atmos., 119, 2457–2473,, 2014. 

Barreto, Á., Cuevas, E., Granados-Muñoz, M.-J., Alados-Arboledas, L., Romero, P. M., Gröbner, J., Kouremeti, N., Almansa, A. F., Stone, T., Toledano, C., Román, R., Sorokin, M., Holben, B., Canini, M., and Yela, M.: The new sun-sky-lunar Cimel CE318-T multiband photometer – a comprehensive performance evaluation, Atmos. Meas. Tech., 9, 631–654,, 2016. 

Benavent-Oltra, J. A., Román, R., Granados-Muñoz, M. J., Pérez-Ramírez, D., Ortiz-Amezcua, P., Denjean, C., Lopatin, A., Lyamani, H., Torres, B., Guerrero-Rascado, J. L., Fuertes, D., Dubovik, O., Chaikovsky, A., Olmo, F. J., Mallet, M., and Alados-Arboledas, L.: Comparative assessment of GRASP algorithm for a dust event over Granada (Spain) during ChArMEx-ADRIMED 2013 campaign, Atmos. Meas. Tech., 10, 4439–4457,, 2017. 

Benedetti, A., Reid, J. S., Knippertz, P., Marsham, J. H., Di Giuseppe, F., Rémy, S., Basart, S., Boucher, O., Brooks, I. M., Menut, L., Mona, L., Laj, P., Pappalardo, G., Wiedensohler, A., Baklanov, A., Brooks, M., Colarco, P. R., Cuevas, E., da Silva, A., Escribano, J., Flemming, J., Huneeus, N., Jorba, O., Kazadzis, S., Kinne, S., Popp, T., Quinn, P. K., Sekiyama, T. T., Tanaka, T., and Terradellas, E.: Status and future of numerical atmospheric aerosol prediction with a focus on data requirements, Atmos. Chem. Phys., 18, 10615–10643,, 2018. 

Berkoff, T. A., Sorokin, M., Stone, T., Eck, T. F., Hoff, R., Welton, E., Holben, B., Berkoff, T. A., Sorokin, M., Stone, T., Eck, T. F., Hoff, R., Welton, E., and Holben, B.: Nocturnal Aerosol Optical Depth Measurements with a Small-Aperture Automated Photometer Using the Moon as a Light Source, J. Atmos. Ocean. Tech., 28, 1297–1306,, 2011. 

Bey, I., Jacob, D. J., Yantosca, R. M., Logan, J. A., Field, B. D., Fiore, A. M., Li, Q., Liu, H. Y., Mickley, L. J., and Schultz, M. G.: Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res.-Atmos., 106, 23073–23095,, 2001. 

Bond, T. C. and Bergstrom, R. W.: Light Absorption by Carbonaceous Particles: An Investigative Review, Aerosol Sci. Tech., 40, 27–67,, 2006. 

Bond, T. C., Doherty, S. J., Fahey, D. W., Forster, P. M., Berntsen, T., DeAngelo, B. J., Flanner, M. G., Ghan, S., Kärcher, B., Koch, D., Kinne, S., Kondo, Y., Quinn, P. K., Sarofim, M. C., Schultz, M. G., Schulz, M., Venkataraman, C., Zhang, H., Zhang, S., Bellouin, N., Guttikunda, S. K., Hopke, P. K., Jacobson, M. Z., Kaiser, J. W., Klimont, Z., Lohmann, U., Schwarz, J. P., Shindell, D., Storelvmo, T., Warren, S. G., and Zender, C. S.: Bounding the role of black carbon in the climate system: A scientific assessment, J. Geophys. Res.-Atmos., 118, 5380–5552,, 2013. 

Boucher, O.: Atmospheric aerosols?: properties and climate impacts, 1st Edn., Springer Netherlands, 2015. 

Brasseur, G. and Jacob, D. J.: Modeling of atmospheric chemistry, Cambridge University Press, 2017. 

Chen, C., Dubovik, O., Henze, D. K., Lapyonak, T., Chin, M., Ducos, F., Litvinov, P., Huang, X., and Li, L.: Retrieval of desert dust and carbonaceous aerosol emissions over Africa from POLDER/PARASOL products generated by the GRASP algorithm, Atmos. Chem. Phys., 18, 12551–12580,, 2018. 

Chin, M., Ginoux, P., Kinne, S., Torres, O., Holben, B. N., Duncan, B. N., Martin, R. V., Logan, J. A., Higurashi, A., and Nakajima, T.: Tropospheric Aerosol Optical Thickness from the GOCART Model and Comparisons with Satellite and Sun Photometer Measurements, J. Atmos. Sci., 59, 461–483,<0461:TAOTFT>2.0.CO;2, 2002. 

Chin, M., Diehl, T., Tan, Q., Prospero, J. M., Kahn, R. A., Remer, L. A., Yu, H., Sayer, A. M., Bian, H., Geogdzhayev, I. V., Holben, B. N., Howell, S. G., Huebert, B. J., Hsu, N. C., Kim, D., Kucsera, T. L., Levy, R. C., Mishchenko, M. I., Pan, X., Quinn, P. K., Schuster, G. L., Streets, D. G., Strode, S. A., Torres, O., and Zhao, X.-P.: Multi-decadal aerosol variations from 1980 to 2009: a perspective from observations and a global model, Atmos. Chem. Phys., 14, 3657–3690,, 2014. 

Chubarova, N., Nezval', Ye., Sviridenkov, I., Smirnov, A., and Slutsker, I.: Smoke aerosol and its radiative effects during extreme fire event over Central Russia in summer 2010, Atmos. Meas. Tech., 5, 557–568,, 2012. 

Cohen, J. B. and Wang, C.: Estimating global black carbon emissions using a top-down Kalman Filter approach, J. Geophys. Res.-Atmos., 119, 307–323,, 2014. 

Curci, G., Alyuz, U., Barò, R., Bianconi, R., Bieser, J., Christensen, J. H., Colette, A., Farrow, A., Francis, X., Jiménez-Guerrero, P., Im, U., Liu, P., Manders, A., Palacios-Peña, L., Prank, M., Pozzoli, L., Sokhi, R., Solazzo, E., Tuccella, P., Unal, A., Vivanco, M. G., Hogrefe, C., and Galmarini, S.: Modelling black carbon absorption of solar radiation: combining external and internal mixing assumptions, Atmos. Chem. Phys., 19, 181–204,, 2019. 

Deschamps, P.-Y., Breon, F.-M., Leroy, M., Podaire, A., Bricaud, A., Buriez, J.-C., and Seze, G.: The POLDER mission: instrument characteristics and scientific objectives, IEEE T. Geosci. Remote, 32, 598–615,, 1994. 

Deuzé, J. L., Herman, M., Goloub, P., Tanré, D., and Marchand, A.: Characterization of aerosols over ocean from POLDER/ADEOS-1, Geophys. Res. Lett., 26, 1421–1424,, 1999. 

Deuzé, J. L., Bréon, F. M., Devaux, C., Goloub, P., Herman, M., Lafrance, B., Maignan, F., Marchand, A., Nadal, F., Perry, G., and Tanré, D.: Remote sensing of aerosols over land surfaces from POLDER-ADEOS-1 polarized measurements, J. Geophys. Res.-Atmos., 106, 4913–4926,, 2001. 

Dubovik, O.: Optimization of Numerical Inversion in Photopolarimetric Remote Sensing, in: Photopolarimetry in Remote Sensing, Kluwer Academic Publishers, Dordrecht, 65–106, 2004. 

Dubovik, O. and King, M. D.: A flexible inversion algorithm for retrieval of aerosol optical properties from Sun and sky radiance measurements, J. Geophys. Res.-Atmos., 105, 20673–20696,, 2000. 

Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., and Slutsker, I.: Accuracy assessments of aerosol optical properties retrieved from Aerosol Robotic Network (AERONET) Sun and sky radiance measurements, J. Geophys. Res.-Atmos., 105, 9791–9806,, 2000. 

Dubovik, O., Holben, B. N., Lapyonok, T., Sinyuk, A., Mishchenko, M. I., Yang, P., and Slutsker, I.: Non-spherical aerosol retrieval method employing light scattering by spheroids, Geophys. Res. Lett., 29, 54-1–54–4,, 2002. 

Dubovik, O., Sinyuk, A., Lapyonok, T., Holben, B. N., Mishchenko, M., Yang, P., Eck, T. F., Volten, H., Muñoz, O., Veihelmann, B., van der Zande, W. J., Leon, J.-F., Sorokin, M., and Slutsker, I.: Application of spheroid models to account for aerosol particle nonsphericity in remote sensing of desert dust, J. Geophys. Res., 111, D11208,, 2006. 

Dubovik, O., Lapyonok, T., Kaufman, Y. J., Chin, M., Ginoux, P., Kahn, R. A., and Sinyuk, A.: Retrieving global aerosol sources from satellites using inverse modeling, Atmos. Chem. Phys., 8, 209–250,, 2008. 

Dubovik, O., Herman, M., Holdak, A., Lapyonok, T., Tanré, D., Deuzé, J. L., Ducos, F., Sinyuk, A., and Lopatin, A.: Statistically optimized inversion algorithm for enhanced retrieval of aerosol properties from spectral multi-angle polarimetric satellite observations, Atmos. Meas. Tech., 4, 975–1018,, 2011. 

Dubovik, O., Lapyonok, T., Litvinov, P., Herman, M., Fuertes, D., Ducos, F., Torres, B., Derimian, Y., Huang, X., Lopatin, A., Chaikovsky, A., Aspetsberger, M., and Federspiel, C.: GRASP: a versatile algorithm for characterizing the atmosphere, SPIE Newsroom,, 2014. 

Dubovik, O., Li, Z., Mishchenko, M. I., Tanré, D., Karol, Y., Bojkov, B., Cairns, B., Diner, D. J., Espinosa, W. R., Goloub, P., Gu, X., Hasekamp, O., Hong, J., Hou, W., Knobelspiesse, K. D., Landgraf, J., Li, L., Litvinov, P., Liu, Y., Lopatin, A., Marbach, T., Maring, H., Martins, V., Meijer, Y., Milinevsky, G., Mukai, S., Parol, F., Qiao, Y., Remer, L., Rietjens, J., Sano, I., Stammes, P., Stamnes, S., Sun, X., Tabary, P., Travis, L. D., Waquet, F., Xu, F., Yan, C., and Yin, D.: Polarimetric remote sensing of atmospheric aerosols: Instruments, methodologies, results, and perspectives, J. Quant. Spectrosc. Ra.., 224, 474–511,, 2019. 

Escribano, J., Boucher, O., Chevallier, F., and Huneeus, N.: Subregional inversion of North African dust sources, J. Geophys. Res.-Atmos., 121, 8549–8566,, 2016. 

Escribano, J., Boucher, O., Chevallier, F., and Huneeus, N.: Impact of the choice of the satellite aerosol optical depth product in a sub-regional dust emission inversion, Atmos. Chem. Phys., 17, 7111–7126,, 2017. 

Espinosa, W. R., Remer, L. A., Dubovik, O., Ziemba, L., Beyersdorf, A., Orozco, D., Schuster, G., Lapyonok, T., Fuertes, D., and Martins, J. V.: Retrievals of aerosol optical and microphysical properties from Imaging Polar Nephelometer scattering measurements, Atmos. Meas. Tech., 10, 811–824,, 2017. 

Fairlie, T. D., Jacob, D. J., and Park, R. J.: The impact of transpacific transport of mineral dust in the United States, Atmos. Environ., 41, 1251–1266,, 2007. 

Feng, Y., Ramanathan, V., and Kotamarthi, V. R.: Brown carbon: a significant atmospheric absorber of solar radiation?, Atmos. Chem. Phys., 13, 8607–8621,, 2013. 

Ginoux, P., Chin, M., Tegen, I., Prospero, J. M., Holben, B., Dubovik, O., and Lin, S.-J.: Sources and distributions of dust aerosols simulated with the GOCART model, J. Geophys. Res.-Atmos., 106, 20255–20273,, 2001. 

Ginoux, P., Horowitz, L. W., Ramaswamy, V., Geogdzhayev, I. V., Holben, B. N., Stenchikov, G., and Tie, X.: Evaluation of aerosol distribution and optical depth in the Geophysical Fluid Dynamics Laboratory coupled model CM2.1 for present climate, J. Geophys. Res., 111, D22210,, 2006. 

Gorchakova, I. A. and Mokhov, I. I.: The radiative and thermal effects of smoke aerosol over the region of Moscow during the summer fires of 2010, Izv. Atmos. Ocean. Phys., 48, 496–503,, 2012. 

Granier, C., Bessagnet, B., Bond, T., D'Angiola, A., Denier van der Gon, H., Frost, G. J., Heil, A., Kaiser, J. W., Kinne, S., Klimont, Z., Kloster, S., Lamarque, J.-F., Liousse, C., Masui, T., Meleux, F., Mieville, A., Ohara, T., Raut, J.-C., Riahi, K., Schultz, M. G., Smith, S. J., Thompson, A., van Aardenne, J., van der Werf, G. R., and van Vuuren, D. P.: Evolution of anthropogenic and biomass burning emissions of air pollutants at global and regional scales during the 1980–2010 period, Climatic Change, 109, 163–190,, 2011. 

Henze, D. K., Hakami, A., and Seinfeld, J. H.: Development of the adjoint of GEOS-Chem, Atmos. Chem. Phys., 7, 2413–2433,, 2007. 

Henze, D. K., Seinfeld, J. H., and Shindell, D. T.: Inverse modeling and mapping US air quality influences of inorganic PM2.5 precursor emissions using the adjoint of GEOS-Chem, Atmos. Chem. Phys., 9, 5877–5903,, 2009. 

Holben, B. N., Eck, T. F., Slutsker, I., Tanré, D., Buis, J. P., Setzer, A., Vermote, E., Reagan, J. A., Kaufman, Y. J., Nakajima, T., Lavenu, F., Jankowiak, I., and Smirnov, A.: AERONET – A Federated Instrument Network and Data Archive for Aerosol Characterization, Remote Sens. Environ., 66, 1–16,, 1998. 

Hsu, N. C., Tsay, S.-C., King, M. D., and Herman, J. R.: Aerosol Properties Over Bright-Reflecting Source Regions, IEEE T. Geosci. Remote, 42, 557–569,, 2004. 

Hsu, N. C., Jeong, M.-J., Bettenhausen, C., Sayer, A. M., Hansell, R., Seftor, C. S., Huang, J., and Tsay, S.-C.: Enhanced Deep Blue aerosol retrieval algorithm: The second generation, J. Geophys. Res.-Atmos., 118, 9296–9315,, 2013. 

HTAP: UNECE: Hemispheric Transport of Air Pollution 2010: Part A: Ozone and Particulate Matter, Air Pollution Studies No. 17, edited by: Dentener, F., Keating, T., and Akimoto, H., UNECE Information Service, Geneva, Switzerland, 2010. 

Hu, Q., Goloub, P., Veselovskii, I., Bravo-Aranda, J.-A., Popovici, I. E., Podvin, T., Haeffelin, M., Lopatin, A., Dubovik, O., Pietras, C., Huang, X., Torres, B., and Chen, C.: Long-range-transported Canadian smoke plumes in the lower stratosphere over northern France, Atmos. Chem. Phys., 19, 1173–1193,, 2019. 

Huijnen, V., Flemming, J., Kaiser, J. W., Inness, A., Leitão, J., Heil, A., Eskes, H. J., Schultz, M. G., Benedetti, A., Hadji-Lazaro, J., Dufour, G., and Eremenko, M.: Hindcast experiments of tropospheric composition during the summer 2010 fires over western Russia, Atmos. Chem. Phys., 12, 4341–4364,, 2012. 

Huneeus, N., Schulz, M., Balkanski, Y., Griesfeller, J., Prospero, J., Kinne, S., Bauer, S., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Fillmore, D., Ghan, S., Ginoux, P., Grini, A., Horowitz, L., Koch, D., Krol, M. C., Landing, W., Liu, X., Mahowald, N., Miller, R., Morcrette, J.-J., Myhre, G., Penner, J., Perlwitz, J., Stier, P., Takemura, T., and Zender, C. S.: Global dust model intercomparison in AeroCom phase I, Atmos. Chem. Phys., 11, 7781–7816,, 2011. 

Huneeus, N., Chevallier, F., and Boucher, O.: Estimating aerosol emissions by assimilating observed aerosol optical depth in a global aerosol model, Atmos. Chem. Phys., 12, 4585–4606,, 2012. 

Huneeus, N., Boucher, O., and Chevallier, F.: Atmospheric inversion of SO2 and primary aerosol emissions for the year 2010, Atmos. Chem. Phys., 13, 6555–6573,, 2013. 

Jacob, D. J.: Introduction to Atmospheric Chemistry, Princeton University Press, 1999. 

Jaeglé, L., Quinn, P. K., Bates, T. S., Alexander, B., and Lin, J.-T.: Global distribution of sea salt aerosols: new constraints from in situ and remote sensing observations, Atmos. Chem. Phys., 11, 3137–3157,, 2011. 

Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Dentener, F., Muntean, M., Pouliot, G., Keating, T., Zhang, Q., Kurokawa, J., Wankmüller, R., Denier van der Gon, H., Kuenen, J. J. P., Klimont, Z., Frost, G., Darras, S., Koffi, B., and Li, M.: HTAP_v2.2: a mosaic of regional and global emission grid maps for 2008 and 2010 to study hemispheric transport of air pollution, Atmos. Chem. Phys., 15, 11411–11432,, 2015. 

Jethva, H., Torres, O., and Ahn, C.: Global assessment of OMI aerosol single-scattering albedo using ground-based AERONET inversion, J. Geophys. Res.-Atmos., 119, 9020–9040,, 2014. 

Kaufman, Y. J., Tanré, D., Remer, L. A., Vermote, E. F., Chu, A., and Holben, B. N.: Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer, J. Geophys. Res., 102, 17051–17067,, 1997. 

King, M. D. and Dubovik, O.: Determination of aerosol optical properties from inverse methods, in: Aerosol Remote Sensing, Springer Berlin Heidelberg, Berlin, Heidelberg, 101–136, 2013. 

Kinne, S., Lohmann, U., Feichter, J., Schulz, M., Timmreck, C., Ghan, S., Easter, R., Chin, M., Ginoux, P., Takemura, T., Tegen, I., Koch, D., Herzog, M., Penner, J., Pitari, G., Holben, B., Eck, T., Smirnov, A., Dubovik, O., Slutsker, I., Tanré, D., Torres, O., Mishchenko, M., Geogdzhayev, I., Chu, D. A., and Kaufman, Y.: Monthly averages of aerosol properties: A global comparison among models, satellite data, and AERONET ground data, J. Geophys. Res., 108, 4634,, 2003. 

Kinne, S., Schulz, M., Textor, C., Guibert, S., Balkanski, Y., Bauer, S. E., Berntsen, T., Berglen, T. F., Boucher, O., Chin, M., Collins, W., Dentener, F., Diehl, T., Easter, R., Feichter, J., Fillmore, D., Ghan, S., Ginoux, P., Gong, S., Grini, A., Hendricks, J., Herzog, M., Horowitz, L., Isaksen, I., Iversen, T., Kirkevåg, A., Kloster, S., Koch, D., Kristjansson, J. E., Krol, M., Lauer, A., Lamarque, J. F., Lesins, G., Liu, X., Lohmann, U., Montanaro, V., Myhre, G., Penner, J., Pitari, G., Reddy, S., Seland, O., Stier, P., Takemura, T., and Tie, X.: An AeroCom initial assessment – optical properties in aerosol component modules of global models, Atmos. Chem. Phys., 6, 1815–1834,, 2006. 

Kok, J. F., Ridley, D. A., Zhou, Q., Miller, R. L., Zhao, C., Heald, C. L., Ward, D. S., Albani, S., and Haustein, K.: Smaller desert dust cooling effect estimated from analysis of dust size and abundance, Nat. Geosci., 10, 274–278,, 2017. 

Kokhanovsky, A. A., Davis, A. B., Cairns, B., Dubovik, O., Hasekamp, O. P., Sano, I., Mukai, S., Rozanov, V. V., Litvinov, P., Lapyonok, T., Kolomiets, I. S., Oberemok, Y. A., Savenkov, S., Martin, W., Wasilewski, A., Di Noia, A., Stap, F. A., Rietjens, J., Xu, F., Natraj, V., Duan, M., Cheng, T., and Munro, R.: Space-based remote sensing of atmospheric aerosols: The multi-angle spectro-polarimetric frontier, Earth-Sci. Rev., 145, 85–116,, 2015. 

Lack, D. A., Langridge, J. M., Bahreini, R., Cappa, C. D., Middlebrook, A. M., and Schwarz, J. P.: Brown carbon and internal mixing in biomass burning particles, P. Natl. Acad. Sci. USA, 109, 14802–14807,, 2012. 

Lesins, G., Chylek, P., and Lohmann, U.: A study of internal and external mixing scenarios and its effect on aerosol optical properties and direct radiative forcing, J. Geophys. Res.-Atmos., 107, AAC 5-1–AAC 5-12,, 2002. 

Levy, R. C., Mattoo, S., Munchak, L. A., Remer, L. A., Sayer, A. M., Patadia, F., and Hsu, N. C.: The Collection 6 MODIS aerosol products over land and ocean, Atmos. Meas. Tech., 6, 2989–3034,, 2013. 

Li, L., Dubovik, O., Derimian, Y., Schuster, G. L., Lapyonok, T., Litvinov, P., Ducos, F., Fuertes, D., Chen, C., Li, Z., Lopatin, A., Torres, B., and Che, H.: Retrieval of aerosol components directly from satellite and ground-based measurements, Atmos. Chem. Phys., 19, 13409–13443,, 2019. 

Liu, L. and Mishchenko, M.: Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study, Remote Sens., 10, 1634,, 2018. 

Liu, X., Easter, R. C., Ghan, S. J., Zaveri, R., Rasch, P., Shi, X., Lamarque, J.-F., Gettelman, A., Morrison, H., Vitt, F., Conley, A., Park, S., Neale, R., Hannay, C., Ekman, A. M. L., Hess, P., Mahowald, N., Collins, W., Iacono, M. J., Bretherton, C. S., Flanner, M. G., and Mitchell, D.: Toward a minimal representation of aerosols in climate models: description and evaluation in the Community Atmosphere Model CAM5, Geosci. Model Dev., 5, 709–739,, 2012. 

Lopatin, A., Dubovik, O., Chaikovsky, A., Goloub, P., Lapyonok, T., Tanré, D., and Litvinov, P.: Enhancement of aerosol characterization using synergy of lidar and sun-photometer coincident observations: the GARRLiC algorithm, Atmos. Meas. Tech., 6, 2065–2088,, 2013. 

Lund, M. T., Samset, B. H., Skeie, R. B., Watson-Parris, D., Katich, J. M., Schwarz, J. P., and Weinzierl, B.: Short Black Carbon lifetime inferred from a global set of aircraft observations, npj Clim. Atmos. Sci., 1, 31,, 2018. 

Malm, W. C., Sisler, J. F., Huffman, D., Eldred, R. A., and Cahill, T. A.: Spatial and seasonal trends in particle concentration and optical extinction in the United States, J. Geophys. Res., 99, 1347,, 1994. 

Malm, W. C., Day, D. E., Kreidenweis, S. M., Collett, J. L., and Lee, T.: Humidity-dependent optical properties of fine particles during the Big Bend Regional Aerosol and Visibility Observational Study, J. Geophys. Res.-Atmos., 108, 4279,, 2003. 

Martin, R. V., Jacob, D. J., Yantosca, R. M., Chin, M., and Ginoux, P.: Global and regional decreases in tropospheric oxidants from photochemical effects of aerosols, J. Geophys. Res.-Atmos., 108, 4097,, 2003. 

Myhre, G., Samset, B. H., Schulz, M., Balkanski, Y., Bauer, S., Berntsen, T. K., Bian, H., Bellouin, N., Chin, M., Diehl, T., Easter, R. C., Feichter, J., Ghan, S. J., Hauglustaine, D., Iversen, T., Kinne, S., Kirkevåg, A., Lamarque, J.-F., Lin, G., Liu, X., Lund, M. T., Luo, G., Ma, X., van Noije, T., Penner, J. E., Rasch, P. J., Ruiz, A., Seland, Ø., Skeie, R. B., Stier, P., Takemura, T., Tsigaridis, K., Wang, P., Wang, Z., Xu, L., Yu, H., Yu, F., Yoon, J.-H., Zhang, K., Zhang, H., and Zhou, C.: Radiative forcing of the direct aerosol effect from AeroCom Phase II simulations, Atmos. Chem. Phys., 13, 1853–1877,, 2013. 

Ocko, I. B. and Ginoux, P. A.: Comparing multiple model-derived aerosol optical properties to spatially collocated ground-based and satellite measurements, Atmos. Chem. Phys., 17, 4451–4475,, 2017. 

Park, R. J., Jacob, D. J., Chin, M., and Martin, R. V.: Sources of carbonaceous aerosols over the United States and implications for natural visibility, J. Geophys. Res., 108, 4355,, 2003. 

Park, R. J., Jacob, D. J., Field, B. D., Yantosca, R. M., and Chin, M.: Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: Implications for policy, J. Geophys. Res., 109, D15204,, 2004. 

Péré, J. C., Bessagnet, B., Mallet, M., Waquet, F., Chiapello, I., Minvielle, F., Pont, V., and Menut, L.: Direct radiative effect of the Russian wildfires and its impact on air temperature and atmospheric dynamics during August 2010, Atmos. Chem. Phys., 14, 1999–2013,, 2014. 

Popp, T., de Leeuw, G., Bingen, C., Brühl, C., Capelle, V., Chedin, A., Clarisse, L., Dubovik, O., Grainger, R., Griesfeller, J., Heckel, A., Kinne, S., Klüser, L., Kosmale, M., Kolmonen, P., Lelli, L., Litvinov, P., Mei, L., North, P., Pinnock, S., Povey, A., Robert, C., Schulz, M., Sogacheva, L., Stebel, K., Stein Zweers, D., Thomas, G., Tilstra, L., Vandenbussche, S., Veefkind, P., Vountas, M., and Xue, Y.: Development, Production and Evaluation of Aerosol Climate Data Records from European Satellite Observations (Aerosol_cci), Remote Sens., 8, 421,, 2016. 

Pozzer, A., de Meij, A., Yoon, J., Tost, H., Georgoulias, A. K., and Astitha, M.: AOD trends during 2001–2010 from observations and model simulations, Atmos. Chem. Phys., 15, 5521–5535,, 2015. 

Randerson, J. T., Chen, Y., van der Werf, G. R., Rogers, B. M., and Morton, D. C.: Global burned area and biomass burning emissions from small fires, J. Geophys. Res.-Biogeo., 117, G04012,, 2012. 

Remer, L. A., Kaufman, Y. J., Tanré, D., Mattoo, S., Chu, D. A., Martins, J. V., Li, R.-R., Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., Holben, B. N., Remer, L. A., Kaufman, Y. J., Tanré, D., Mattoo, S., Chu, D. A., Martins, J. V., Li, R.-R., Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., and Holben, B. N.: The MODIS Aerosol Algorithm, Products, and Validation, J. Atmos. Sci., 62, 947–973,, 2005. 

R'Honi, Y., Clarisse, L., Clerbaux, C., Hurtmans, D., Duflot, V., Turquety, S., Ngadi, Y., and Coheur, P.-F.: Exceptional emissions of NH3 and HCOOH in the 2010 Russian wildfires, Atmos. Chem. Phys., 13, 4171–4181,, 2013. 

Ridley, D. A., Heald, C. L., and Ford, B.: North African dust export and deposition: A satellite and model perspective, J. Geophys. Res.-Atmos., 117, D02202,, 2012. 

Ridley, D. A., Heald, C. L., Kok, J. F., and Zhao, C.: An observationally constrained estimate of global dust aerosol optical depth, Atmos. Chem. Phys., 16, 15097–15117,, 2016. 

Román, R., Torres, B., Fuertes, D., Cachorro, V. E., Dubovik, O., Toledano, C., Cazorla, A., Barreto, A., Bosch, J. L., Lapyonok, T., González, R., Goloub, P., Perrone, M. R., Olmo, F. J., de Frutos, A., and Alados-Arboledas, L.: Remote sensing of lunar aureole with a sky camera: Adding information in the nocturnal retrieval of aerosol properties with GRASP code, Remote Sens. Environ., 196, 238–252,, 2017. 

Román, R., Benavent-Oltra, J. A., Casquero-Vera, J. A., Lopatin, A., Cazorla, A., Lyamani, H., Denjean, C., Fuertes, D., Pérez-Ramírez, D., Torres, B., Toledano, C., Dubovik, O., Cachorro, V. E., de Frutos, A. M., Olmo, F. J., and Alados-Arboledas, L.: Retrieval of aerosol profiles combining sunphotometer and ceilometer measurements in GRASP code, Atmos. Res., 204, 161–177,, 2018. 

Russell, P. B., Bergstrom, R. W., Shinozuka, Y., Clarke, A. D., DeCarlo, P. F., Jimenez, J. L., Livingston, J. M., Redemann, J., Dubovik, O., and Strawa, A.: Absorption Angstrom Exponent in AERONET and related data as an indicator of aerosol composition, Atmos. Chem. Phys., 10, 1155–1169,, 2010. 

Samset, B. H., Stjern, C. W., Andrews, E., Kahn, R. A., Myhre, G., Schulz, M., and Schuster, G. L.: Aerosol Absorption: Progress Towards Global and Regional Constraints, Curr. Clim. Chang. Reports, 4, 65–83,, 2018. 

Sato, M., Hansen, J., Koch, D., Lacis, A., Ruedy, R., Dubovik, O., Holben, B., Chin, M., and Novakov, T.: Global atmospheric black carbon inferred from AERONET, P. Natl. Acad. Sci. USA, 100, 6319–6324,, 2003. 

Sayer, A. M., Hsu, N. C., Bettenhausen, C., and Jeong, M.-J.: Validation and uncertainty estimates for MODIS Collection 6 “Deep Blue” aerosol data, J. Geophys. Res.-Atmos., 118, 7864–7872,, 2013. 

Sayer, A. M., Munchak, L. A., Hsu, N. C., Levy, R. C., Bettenhausen, C., and Jeong, M.-J.: MODIS Collection 6 aerosol products: Comparison between Aqua's e-Deep Blue, Dark Target, and “merged” data sets, and usage recommendations, J. Geophys. Res.-Atmos., 119, 13965–13989,, 2014. 

Sayer, A. M., Hsu, N. C., Lee, J., Kim, W. V., Dubovik, O., Dutcher, S. T., Huang, D., Litvinov, P., Lyapustin, A., Tackett, J. L., and Winker, D. M.: Validation of SOAR VIIRS Over-Water Aerosol Retrievals and Context Within the Global Satellite Aerosol Data Record, J. Geophys. Res.-Atmos., 123, 2018JD029465,, 2018. 

Schulz, M., Textor, C., Kinne, S., Balkanski, Y., Bauer, S., Berntsen, T., Berglen, T., Boucher, O., Dentener, F., Guibert, S., Isaksen, I. S. A., Iversen, T., Koch, D., Kirkevåg, A., Liu, X., Montanaro, V., Myhre, G., Penner, J. E., Pitari, G., Reddy, S., Seland, Ø., Stier, P., and Takemura, T.: Radiative forcing by aerosols as derived from the AeroCom present-day and pre-industrial simulations, Atmos. Chem. Phys., 6, 5225–5246,, 2006. 

Schuster, G. L., Dubovik, O., Arola, A., Eck, T. F., and Holben, B. N.: Remote sensing of soot carbon – Part 2: Understanding the absorption Ångström exponent, Atmos. Chem. Phys., 16, 1587–1602,, 2016. 

Schuster, G. L., Espinosa, W. R., Ziemba, L. D., Beyersdorf, A. J., Rocha-Lima, A., Anderson, B. E., Martins, J. V., Dubovik, O., Ducos, F., Fuertes, D., Lapyonok, T., Shook, M., Derimian, Y., and Moore, R. H.: A Laboratory Experiment for the Statistical Evaluation of Aerosol Retrieval (STEAR) Algorithms, Remote Sens., 11, 498,, 2019. 

Schutgens, N. A. J., Gryspeerdt, E., Weigum, N., Tsyro, S., Goto, D., Schulz, M., and Stier, P.: Will a perfect model agree with perfect observations? The impact of spatial sampling, Atmos. Chem. Phys., 16, 6335–6353,, 2016. 

Smirnov, A., Holben, B. N., Eck, T. F., Dubovik, O., and Slutsker, I.: Cloud-Screening and Quality Control Algorithms for the AERONET Database, Remote Sens. Environ., 73, 337–349,, 2000. 

Tanré, D., Kaufman, Y. J., Herman, M., and Mattoo, S.: Remote sensing of aerosol properties over oceans using the MODIS/EOS spectral radiances, J. Geophys. Res.-Atmos., 102, 16971–16988,, 1997. 

Tanré, D., Bréon, F. M., Deuzé, J. L., Dubovik, O., Ducos, F., François, P., Goloub, P., Herman, M., Lifermann, A., and Waquet, F.: Remote sensing of aerosols by using polarized, directional and spectral measurements within the A-Train: the PARASOL mission, Atmos. Meas. Tech., 4, 1383–1395,, 2011. 

Tegen, I., Neubauer, D., Ferrachat, S., Siegenthaler-Le Drian, C., Bey, I., Schutgens, N., Stier, P., Watson-Parris, D., Stanelle, T., Schmidt, H., Rast, S., Kokkola, H., Schultz, M., Schroeder, S., Daskalakis, N., Barthel, S., Heinold, B., and Lohmann, U.: The global aerosol–climate model ECHAM6.3–HAM2.3 – Part 1: Aerosol evaluation, Geosci. Model Dev., 12, 1643–1677,, 2019. 

Textor, C., Schulz, M., Guibert, S., Kinne, S., Balkanski, Y., Bauer, S., Berntsen, T., Berglen, T., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Feichter, H., Fillmore, D., Ghan, S., Ginoux, P., Gong, S., Grini, A., Hendricks, J., Horowitz, L., Huang, P., Isaksen, I., Iversen, I., Kloster, S., Koch, D., Kirkevåg, A., Kristjansson, J. E., Krol, M., Lauer, A., Lamarque, J. F., Liu, X., Montanaro, V., Myhre, G., Penner, J., Pitari, G., Reddy, S., Seland, Ø., Stier, P., Takemura, T., and Tie, X.: Analysis and quantification of the diversities of aerosol life cycles within AeroCom, Atmos. Chem. Phys., 6, 1777–1813,, 2006. 

Textor, C., Schulz, M., Guibert, S., Kinne, S., Balkanski, Y., Bauer, S., Berntsen, T., Berglen, T., Boucher, O., Chin, M., Dentener, F., Diehl, T., Feichter, J., Fillmore, D., Ginoux, P., Gong, S., Grini, A., Hendricks, J., Horowitz, L., Huang, P., Isaksen, I. S. A., Iversen, T., Kloster, S., Koch, D., Kirkevåg, A., Kristjansson, J. E., Krol, M., Lauer, A., Lamarque, J. F., Liu, X., Montanaro, V., Myhre, G., Penner, J. E., Pitari, G., Reddy, M. S., Seland, Ø., Stier, P., Takemura, T., and Tie, X.: The effect of harmonized emissions on aerosol properties in global models – an AeroCom experiment, Atmos. Chem. Phys., 7, 4489–4501,, 2007. 

Torres, B., Dubovik, O., Fuertes, D., Schuster, G., Cachorro, V. E., Lapyonok, T., Goloub, P., Blarel, L., Barreto, A., Mallet, M., Toledano, C., and Tanré, D.: Advanced characterisation of aerosol size properties from measurements of spectral optical depth using the GRASP algorithm, Atmos. Meas. Tech., 10, 3743–3781,, 2017. 

Torres, O., Tanskanen, A., Veihelmann, B., Ahn, C., Braak, R., Bhartia, P. K., Veefkind, P., and Levelt, P.: Aerosols and surface UV products from Ozone Monitoring Instrument observations: An overview, J. Geophys. Res., 112, D24S47,, 2007. 

Torres, O., Ahn, C., and Chen, Z.: Improvements to the OMI near-UV aerosol algorithm using A-train CALIOP and AIRS observations, Atmos. Meas. Tech., 6, 3257–3270,, 2013. 

Tsekeri, A., Lopatin, A., Amiridis, V., Marinou, E., Igloffstein, J., Siomos, N., Solomos, S., Kokkalis, P., Engelmann, R., Baars, H., Gratsea, M., Raptis, P. I., Binietoglou, I., Mihalopoulos, N., Kalivitis, N., Kouvarakis, G., Bartsotas, N., Kallos, G., Basart, S., Schuettemeyer, D., Wandinger, U., Ansmann, A., Chaikovsky, A. P., and Dubovik, O.: GARRLiC and LIRIC: strengths and limitations for the characterization of dust and marine particles along with their mixtures, Atmos. Meas. Tech., 10, 4995–5016,, 2017. 

Wang, J., Xu, X., Henze, D. K., Zeng, J., Ji, Q., Tsay, S.-C., and Huang, J.: Top-down estimate of dust emissions through integration of MODIS and MISR aerosol retrievals with the GEOS-Chem adjoint model, Geophys. Res. Lett., 39, L08802,, 2012. 

Wang, Q., Jacob, D. J., Spackman, J. R., Perring, A. E., Schwarz, J. P., Moteki, N., Marais, E. A., Ge, C., Wang, J., and Barrett, S. R. H.: Global budget and radiative forcing of black carbon aerosol: Constraints from pole-to-pole (HIPPO) observations across the Pacific, J. Geophys. Res.-Atmos., 119, 195–206,, 2014. 

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720,, 2017. 

Xian, P., Reid, J. S., Hyer, E. J., Sampson, C. R., Rubin, J. I., Ades, M., Asencio, N., Basart, S., Benedetti, A., Bhattacharjee, P. S., Brooks, M. E., Colarco, P. R., da Silva, A. M., Eck, T. F., Guth, J., Jorba, O., Kouznetsov, R., Kipling, Z., Sofiev, M., Perez Garcia-Pando, C., Pradhan, Y., Tanaka, T., Wang, J., Westphal, D. L., Yumimoto, K. and Zhang, J.: Current state of the global operational aerosol multi-model ensemble: An update from the International Cooperative for Aerosol Prediction (ICAP), Q. J. R. Meteor. Soc., 145, 176–209,, 2019. 

Xu, X., Wang, J., Henze, D. K., Qu, W., and Kopacz, M.: Constraints on aerosol sources using GEOS-Chem adjoint and MODIS radiances, and evaluation with multisensor (OMI, MISR) data, J. Geophys. Res.-Atmos., 118, 6396–6413,, 2013. 

Zender, C. S., Bian, H., and Newman, D.: Mineral Dust Entrainment and Deposition (DEAD) model: Description and 1990s dust climatology, J. Geophys. Res., 108, 4416,, 2003. 

Zhang, L., Kok, J. F., Henze, D. K., Li, Q., and Zhao, C.: Improving simulations of fine dust surface concentrations over the western United States by optimizing the particle size distribution, Geophys. Res. Lett., 40, 3270–3275,, 2013.  

Zhang, L., Henze, D. K., Grell, G. A., Carmichael, G. R., Bousserez, N., Zhang, Q., Torres, O., Ahn, C., Lu, Z., Cao, J., and Mao, Y.: Constraining black carbon aerosol over Asia using OMI aerosol absorption optical depth and the adjoint of GEOS-Chem, Atmos. Chem. Phys., 15, 10281–10308,, 2015. 

Zhang, S., Penner, J. E., and Torres, O.: Inverse modeling of biomass burning emissions using Total Ozone Mapping Spectrometer aerosol index for 1997, J. Geophys. Res., 110, D21306,, 2005. 

Zheng, B., Chevallier, F., Ciais, P., Yin, Y., and Wang, Y.: On the Role of the Flaming to Smoldering Transition in the Seasonal Cycle of African Fire Emissions, Geophys. Res. Lett., 45, 11998–12007,, 2018. 

Short summary
Global BC, OC and DD aerosol emissions are inverted from POLDER/PARASOL observations for the year 2010 based on the GEOS-Chem inverse modeling framework. The retrieved emissions are 18.4 Tg yr−1 BC, 109.9 Tg yr−1 OC and 731.6 Tg yr−1 DD, which indicate an increase of 166.7 % for BC and 184.0 % for OC, while a decrease of 42.4 % for DD with respect to GEOS-Chem a priori emission inventories is seen. Global annul mean AOD and AAOD resulting from retrieved emissions are 0.119 and 0.0071 at 550 nm.
Final-revised paper