Articles | Volume 21, issue 4
Research article
23 Feb 2021
Research article |  | 23 Feb 2021

Assimilating aerosol optical properties related to size and absorption from POLDER/PARASOL with an ensemble data assimilation system

Athanasios Tsikerdekis, Nick A. J. Schutgens, and Otto P. Hasekamp

A data assimilation system for aerosol, based on an ensemble Kalman filter, has been developed for the ECHAM – Hamburg Aerosol Model (ECHAM-HAM) global aerosol model and applied to POLarization and Directionality of the Earth's Reflectances (POLDER)-derived observations of optical properties. The advantages of this assimilation system is that the ECHAM-HAM aerosol modal scheme carries both aerosol particle numbers and mass which are both used in the data assimilation system as state vectors, while POLDER retrievals in addition to aerosol optical depth (AOD) and the Ångström exponent (AE) also provide information related to aerosol absorption like aerosol absorption optical depth (AAOD) and single scattering albedo (SSA). The developed scheme can simultaneously assimilate combinations of multiple variables (e.g., AOD, AE, SSA) to optimally estimate mass mixing ratio and number mixing ratio of different aerosol species. We investigate the added value of assimilating AE, AAOD and SSA, in addition to the commonly used AOD, by conducting multiple experiments where different combinations of retrieved properties are assimilated. Results are evaluated with (independent) POLDER, Moderate Resolution Imaging Spectroradiometer (MODIS) Dark Target, MODIS Deep Blue and Aerosol Robotic Network (AERONET) observations. The experiment where POLDER AOD, AE and SSA are assimilated shows systematic improvement in mean error, mean absolute error and correlation for AOD, AE, AAOD and SSA compared to the experiment where only AOD is assimilated. The same experiment reduces the global ME against AERONET from 0.072 to 0.001 for AOD, from 0.273 to 0.009 for AE and from 0.012 to 0.002 for AAOD. Additionally, sensitivity experiments reveal the benefits of assimilating AE over AOD at a second wavelength or SSA over AAOD, possibly due to a simpler observation covariance matrix in the present data assimilation framework. We conclude that the currently available AE and SSA do positively impact data assimilation.

1 Introduction

Atmospheric aerosol is a key factor that modifies the effects and intensity of climate change, due to its participation in numerous atmospheric processes that may alter the radiative budget of the planet (Boucher et al., 2013). The diverse aerosol size and chemical composition that affect their transport and removal mechanisms, the complex atmospheric aging processes, the in-cloud condensation growth and the limited as well as indirect information regarding their emission flux and sources make their simulation a challenging task (Huneeus et al., 2011; Kinne et al., 2006; Schutgens and Stier, 2014; Textor et al., 2006) and the aerosol direct, semi-direct and indirect radiative effects very hard to estimate (Carslaw et al., 2013; Fan et al., 2016; Hasekamp et al., 2019b; Koch and Del Genio, 2010; Myhre et al., 2013; Nabat et al., 2014; Tsikerdekis et al., 2017, 2019; Yumimoto and Takemura, 2011).

Data assimilation systems have been employed in the past in order to either adjust the aerosol mixing ratio (Benedetti et al., 2009; Dai et al., 2014, 2019; Schutgens et al., 2010a, b; Di Tomaso et al., 2017; Yumimoto et al., 2007, 2016) or estimate new aerosol emission fluxes (Chen et al., 2018, 2019; Escribano et al., 2017; Huneeus et al., 2012; Pope et al., 2016; Schutgens et al., 2012; Sekiyama et al., 2010; Xu et al., 2013).

Retrieved aerosol optical depth (AOD) from the Moderate Resolution Imaging Spectroradiometer Dark Target algorithm (MODIS-DT) has been extensively assimilated (Benedetti et al., 2009; Dai et al., 2014; Escribano et al., 2017; Huneeus et al., 2012; Schutgens et al., 2012; Di Tomaso et al., 2017; Xu et al., 2013; Yumimoto and Takemura, 2011) or used for validation as independent observations (Dai et al., 2019; Schutgens et al., 2010a, b) in past studies. Studies focused on dust AOD or dust source regions assimilated AOD retrievals from the MODIS Deep Blue algorithm (MODIS-DB) (Escribano et al., 2017; Di Tomaso et al., 2017), while other studies assimilated AOD from the ground-based network of stations Aerosol Robotic Network (AERONET) (Schutgens et al., 2012, 2010a, b) or AOD from the Himawari-8 (Yumimoto et al., 2016, 2018). Retrieved AOD from the POLarization and Directionality of the Earth's Reflectances (POLDER) Generalized Retrieval of Atmosphere and Surface Properties (GRASP) algorithm (Dubovik et al., 2011) was also assimilated by Chen et al. (2019).

Almost all aerosol assimilation systems assimilate AOD. AOD is a quantity that describes the aerosol extinction (scattering plus absorption) in the total column of the atmosphere and is, if all microphysical properties stay the same, related to the amount of aerosols. AOD is affected by the size and absorption of aerosol particles, but assimilating just AOD in one wavelength does not disentangle fine from coarse particles and absorbing from non-absorbing particles. The assimilation of other satellite retrieval products, along with AOD, like the Ångström exponent (AE) and single scattering albedo (SSA), which are more closely linked to size and absorption of aerosol particles, may have a positive impact on data assimilation.

The importance of assimilating total and fine-mode AOD separately was highlighted by Generoso et al. (2007). MODIS total and fine-mode fraction AODs were assimilated by Dubovik et al. (2008) and Huneeus et al. (2012), while recently total and fine-mode fraction AODs from MODIS were assimilated simultaneously for dust-only simulations by Escribano et al. (2017). Assimilating species-specific observations, like dust AOD from the LIdar climatology of Vertical Aerosol Structure for space-based lidar simulation studies (LIVAS) (Amiridis et al., 2013), may also address dust-related biases in the Monitoring Atmospheric Composition and Climate (MACC) reanalysis (Georgoulias et al., 2018). Benefits of aerosol size correction were demonstrated also by simultaneously assimilating AODs in two wavelengths or preferably AOD and AE (Schutgens et al., 2010a, b). In addition, even for remote sensing measurements with relatively high uncertainty on the light-absorbing properties, particle-size-related information and the particle-absorbing properties proved to be highly beneficial for a better representation of aerosol composition (Chen et al., 2019). However, while the AOD is retrieved for at least one wavelength from all satellites, the other observational parameters can be retrieved only by a limited number of remote sensing instruments (e.g., POLDER, Dubovik et al., 2011, Hasekamp et al., 2011; the Ozone Monitoring Instrument (OMI), Torres et al., 2007).

Remote sensing is the best way of obtaining aerosol observations over large regions of the Earth. Satellite instruments do not directly measure aerosol-related information but rather properties of light such as intensity, color and polarization state (Benedetti et al., 2018). Although the direct assimilation of clear-sky radiance has been attempted in the past (Weaver et al., 2007), a typical aerosol data assimilation system assimilates aerosol optical properties, which are obtained using retrieval algorithms that use the satellite clear-sky measurements as input. Satellite sensors and retrieval algorithms introduce uncertainties in these estimates due to satellite radiometric calibrations, aerosol properties assumptions, cloud contamination and surface albedo diverse characteristics (Li et al., 2009). In a data assimilation system, the uncertainty of observations should be defined and given as input.

Global aerosol simulations have shown that the aerosol atmospheric load and microphysical properties between different climate models, or even within the same model with altered parameterizations, are quite diverse (Huneeus et al., 2011; Miller et al., 2006). An ensemble-based data assimilation system uses an ensemble of perturbed simulations to define the uncertainty in states of aerosol in the atmosphere (Schutgens et al., 2010a). This ensemble of simulations can be then adjusted based on aerosol retrievals to derive a new and better estimate of aerosol state, which is represented by the ensemble mean.

A widely used data assimilation method that combines an ensemble of perturbed simulations (a priori or background) and provides a new and better estimate (a posteriori or analysis) is the local ensemble transform Kalman filter (LETKF) (Hunt et al., 2007; Miyoshi and Yamane, 2007). In order to get a new and better estimate of aerosol state, which in our case is defined as the atmospheric aerosol mass mixing ratio and number mixing ratio, LETKF requires two main ingredients: (1) an estimate of the background state and the associated uncertainty; (2) observations and the associated uncertainty. The observations must be related to the state vector (e.g., aerosol mixing ratio); the state vector is converted to simulated observations by a model (in our case, ECHAM-HAM), while the background estimate and its uncertainty may be represented through an ensemble of simulations.

In the present study, we use LETKF to assimilate aerosol optical properties (AOD550, AOD865, AE550−865, AAOD550, SSA550; subscript denotes wavelength in nm) from multi-angle photopolarimetric POLDER measurements retrieved by the algorithm developed at the Netherlands Institute for Space Research (SRON). In Sect. 2, we present POLDER retrievals and the corresponding uncertainty estimate, the ECHAM-HAM global climate–aerosol model, as well as other observational data used as independent observations. Section 3 presents the data assimilation system, the method to produce the perturbed ensemble necessary for LETKF and the experimental setup. Finally, in Sect. 4, the results are partitioned into four distinct segments. The first highlights the importance of combining aerosol optical properties (e.g., AOD, AE and SSA) to acquire a more robust representation of the atmospheric aerosol state, the second evaluates the core experiments with independent observations, the third discusses the preference of some observations over others for the assimilation, and lastly a number of sensitivity experiments are presented regarding some of the parameters in LETKF.

2 Data

2.1 Observational data


AERONET is a global network of ground-based stations of Sun–sky radiometers (Holben et al., 1998) that provides high-quality direct-Sun AOD estimates at various wavelengths from 340 to 1020 nm (Holben et al., 2001) and aerosol inversion SSA estimates (Dubovik and King, 2000). Due to its design, the instrument can provide useful measurements under cloud-free conditions during the day. The AOD uncertainty is estimated to be < 0.02 (Dubovik et al., 2000; Eck et al., 1999) and SSA < 0.03 for AOD440> 0.4 and solar zenith angle greater than 50 (Dubovik et al., 2002; Holben et al., 2006). Thus, AERONET is commonly used as the “ground truth” for the validation of aerosol optical properties, in particular AOD, in both satellite and model studies.

For the evaluation of the assimilated experiments, the V3 L2.0 datasets of the direct-Sun and aerosol inversion datasets are used. Recently, AOD monthly differences between AERONET V3 and V2 datasets were less than 0.002 ± 0.002, highlighting the stability of the network independent of the version (Giles et al., 2019). It is noted that for aerosol absorption properties, AERONET L2.0 aerosol inversion data only include cases with AOD440> 0.4; thus, the evaluation of ECHAM-HAM absorbing properties is restricted to high-AOD cases only. To define the POLDER uncertainty, V3 L1.5 inversion datasets were used because they provide more data points at the cost of accuracy. AERONET L1.5 contains cases of low AOD550 where the retrieval accuracy for aerosol absorption optical depth (AAOD) and SSA is low.

2.1.2 POLDER

The POLDER-3 instrument on the Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) microsatellite, which has been active between 2004 and 2013, has the unique capability of measuring light intensity and polarization properties at multiple viewing angles (up to 16) and multiple wavelengths (0.44 to 1.02 µm). The multi-wavelength and multi-viewing-angle photopolarimetric measurements make better use of the information content of scattered solar radiation in comparison to single-viewing measurements (Hasekamp and Landgraf, 2007; Mishchenko and Travis, 1997); hence, POLDER is an ideal tool for obtaining accurate aerosol microphysical and optical properties. Its native horizontal resolution is 6 km× 6 km but in this study “medium-resolution” data have been used that correspond to 18 km× 18 km. The retrieval algorithm developed at SRON fits a radiative transfer model (Hasekamp, 2005; Schepers et al., 2014) to the multi-angle photopolarimetric measurements of POLDER to derive aerosol optical properties corresponding to a bimodal aerosol size distribution. The retrieved properties for two modes for fine and coarse particles are the effective radius and effective variance, the column number concentration and the real and imaginary parts of the refractive index for each mode, and for the coarse mode additionally the fraction of spherical particles is retrieved (Hasekamp et al., 2011; Lacagnina et al., 2015; Wu et al., 2015). Using the above-mentioned aerosol parameters, for the two modes, AOD, AE, AAOD and SSA can be calculated. The multi-angle multi-wavelength photopolarimetric measurements of POLDER also have the ability to differentiate scattering of cloud droplets from aerosol particles, making it possible to exclude cloud contaminated pixels (Stap et al., 2015). Recently, the algorithm was extended to an arbitrary number of modes (Fu and Hasekamp, 2018), but the present paper uses the bimodal product. In the present study, aggregated (1× 1) POLDER data are used in the assimilation. Global aerosol retrievals from POLDER-3 by the SRON algorithm are available for the year 2006.

The aerosol optical properties of POLDER retrievals demonstrate good agreement with either ground-based (AERONET) or satellite (OMI) retrievals for the year 2006 (Hasekamp et al., 2011; Lacagnina et al., 2015, 2017; Stap et al., 2015). A global evaluation against the AERONET inversion dataset for POLDER ocean and land pixels showed similar results, with absolute differences for three AOD wavelengths and SSA of ± 0.05 (Lacagnina et al., 2015, 2017). Performance over land and ocean pixels is similar for SSA, contrary to AOD, where retrievals over ocean pixels were better (Lacagnina et al., 2017). POLDER AOD agreement with AERONET descends considerably for values below 0.07 over land, which results in some spatial bias patterns, where in low-AOD regions like North America (relative) errors were large and in high-AOD regions like the Sahara errors were small. The largest relative SSA discrepancies are found in Europe and North America (Lacagnina et al., 2017). The recent multi-model version of the algorithm (10 modes instead of two modes) achieved higher accuracy for AOD and similar performance for SSA when compared to AERONET for retrievals over land (Fu and Hasekamp, 2018).

Figure 1The uncertainty and the relative uncertainty of POLDER for AOD550 (a, b) AOD865 (c, d) and AAOD550 (e, f) and the uncertainty for AE550−865 (g) and SSA550 (h) averaged over the period 20 July to 28 August 2006. The global mean of each variable is denoted in the bottom-right corner for each case.

In the present study, the observational uncertainties of POLDER are assessed by evaluating the retrieval against the dataset of AERONET. This approach provides a parameterization of observational uncertainty based on the real errors of POLDER retrievals. Obviously, since AERONET is a spatially sparse ground-based network of stations, some generalization had to be made as far as the performance of POLDER in remote areas. Furthermore, the AERONET retrieval errors, which typically are smaller than the errors of any remote sensing retrieval, at least for AOD, were not taken into account. A detailed description of the POLDER uncertainty estimation is presented in Appendix A.

Global maps of the resulting absolute and relative uncertainties of POLDER averaged over 40 d (the period of our assimilation experiment) for the five variables used in the present study are illustrated in Fig. 1. Also, the AERONET stations used are selected for POLDER retrievals over land, where the errors are higher in comparison to POLDER retrievals over ocean (Lacagnina et al., 2017). For the above-mentioned reasons, in cases of very low AOD550, for example, in the first two AOD550 bins of Fig. A1 and predominantly over ocean pixels, the POLDER uncertainty estimates are probably too conservative. The POLDER uncertainty estimation is based on spatiotemporal POLDER retrievals of a 18 km× 18 km grid with AERONET. The estimated uncertainty was afterwards applied to a coarser 1× 1 grid.

The global mean absolute uncertainties for AOD in the two wavelengths (550 and 865 nm) are quite similar (0.08, 0.06), although AOD865 is lower since the absolute values of AOD865 are lower too. The AE550−865 global mean absolute uncertainty is 0.50, for SSA550 it is 0.085, and for AAOD550 it is 0.012. It can be clearly seen that the error on AE550−865 and SSA550 strongly depends on AOD. For example, the error on SSA is  0.03 in regions with high aerosol loading due to biomass burning, dust outbreaks or industrial pollution but can be  0.10 over the remote ocean.


For a broader spatial coverage, a comparison with the independent observations MODIS Collection 5 Dark Target (MODIS-DT) and MODIS Collection 6 Deep Blue (MODIS-DB) retrievals of AOD550 (Sayer et al., 2014) was conducted. This allows assessment of the performance of the assimilated experiments over ocean and other remote regions, away from AERONET sites. In the case of MODIS-DT, a distinctive version designed specifically for assimilation purposes was used, which was corrected using 4 years of collocated data with AERONET as a basis (Hyer et al., 2011; Shi et al., 2011; Zhang and Reid, 2006). It is noted that both MODIS-DT and MODIS-DB come with their own uncertainties; thus, the comparison with the assimilated experiments cannot be considered a validation within the strict definition of the term. The uncertainties of these products are discussed in Di Tomaso et al. (2017).

2.2 Model simulations

The atmospheric global coupled climate–aerosol modeling system ECHAM-HAM is used as the forward model to generate an ensemble of short-term forecasts. ECHAM is the general circulation part of the modeling system and it simulates the meteorological conditions of the atmosphere in a Gaussian grid, while HAM is the aerosol module that utilizes the meteorological (e.g., wind, turbulence, convection, precipitation) and surface (e.g., surface roughness, bare and vegetated surface fractions) variables of ECHAM to solve the physical and chemical aerosol particle processes.

2.2.1 The ECHAM6-HAM2 aerosol climate model

This study uses the sixth generation of the general circulation model (ECHAM6.3), which was developed at the Max Planck Institute for Meteorology (MPI-M) in Hamburg, Germany (Stevens et al., 2013). The adiabatic processes in the model are based on a spectral-transform dynamical core that simulates some essential meteorological parameters (temperature, surface pressure, vorticity and divergence) while a collection of physical schemes parameterizes the diabatic processes like convection, diffusion, turbulence and gravity waves (Schultz et al., 2018).

We employ ECHAM-HAM with a grid resolution of T63L47 (1.875× 1.875, with 47 mostly tropospheric levels based on a hybrid sigma–pressure coordinate). The radiative transfer calculations in the model are performed by the Rapid Radiative Transfer Model for Global modeling (RRTM-G; Iacono et al., 2008). ECHAM can optionally force the essential meteorological parameters of its adiabatic dynamical core to approach a prescribed field by applying a relaxation technique with time-varying weights (Schultz et al., 2018). Typically, the prescribed field consists of a reanalysis database, like ERA-Interim. It is noted that the physics of the model are not directly influenced by the external nudging data; therefore, ECHAM is still the main driver of the dynamics that are just “nudged” towards a prescribed trajectory that describes the 3-D temperature, surface pressure, vorticity and divergence (Rast et al., 2015). Nudging timescales are 24 h for temperature and surface pressure, 48 h for divergence and 6 h for vorticity.

HAM simulates the physical and chemical processes of aerosol in the atmosphere (Stier et al., 2005; Zhang et al., 2012). The most recent version (HAM2.3) includes new emission schemes for aerosol and aerosol precursors and modified aerosol–cloud interactions that are summarized in Tegen et al. (2019). The M7 aerosol model used in HAM2.3 considers five groups of aerosol species: desert dust (DU), sea salt (SS), organic carbon (OC), black carbon (BC) and sulfates (SO4) (Vignati et al., 2004). Nitrate aerosol particles (NO3), which may be produced by gas-phase nitrate (HNO3) and ammonia (NH3) reactions, are not implemented in HAM2.3. Aerosols are partitioned into seven unimodal lognormal particle size distributions, called modes, separated into two hygroscopic classes (hydrophobic and hydrophilic). Six of these modes consist of an internal mix of various aerosol types. For the nucleation-mode radius, r< 0.005 µm; for the Aitken mode, it is 0.005 µm<r< 0.05 µm; for the accumulation mode, it is 0.05 µm<r< 0.5 µm; and for the coarse mode, it is r> 0.5 µm (Vignati et al., 2004). The cloud and aerosol optical properties are computed using Mie theory for each band of the RRTM-G and organized in lookup tables (Tegen et al., 2019). Absorption and scattering of aerosol particles in the ECHAM-HAM are calculated using the prognostic concentrations of aerosol tracers (Schultz et al., 2018).

All aerosol species are emitted, transported and deposited. Depending on their physical and chemical properties, they can take part in a number of other processes like aerosol–radiation interactions (scattering and absorption), as well as other aerosol microphysical processes (e.g., nucleation, coagulation, aerosol water uptake and cloud activation). The purely natural emitted aerosol types (DU, SS) are introduced to the atmosphere by utilizing the simulated information of ECHAM, mainly the wind and some surface and ocean characteristics. Aerosols that can be emitted or formed chemically by both natural and anthropogenic sources (OC, BC, SO4) are introduced using predefined emission inventories (Zhang et al., 2012).

Sea salt emissions are parameterized using the wind velocity at 10 m as the dominant driver for aerosol particle production, while sea surface temperature (SST) influences mostly the emitted particles' size (Long et al., 2011; Sofiev et al., 2011). In cases where the SST is low, sea salt emissions are lower and the emitted particles are smaller compared to when SST is higher (Sofiev et al., 2011). Sea salt particles are emitted only in the soluble accumulation and coarse modes. Natural emissions of dimethyl sulfide (DMS) over the ocean, which is an aerosol precursor, are calculated online based on the 10 m wind velocity (Nightingale et al., 2000) and the prescribed concentration of DMS on the surface of the ocean (Lana et al., 2011).

Dust emissions are based on the dust source scheme developed by Tegen et al. (2002). Improvements were made in terms of the surface aerodynamic roughness length, soil moisture and soil properties specifically over East Asia by Cheng et al. (2008). Also there were updates regarding the representation of Saharan dust sources using infrared dust index from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) instrument aboard the Meteosat Second Generation satellite by Heinold et al. (2016). Wind velocity at 10 m is the main driver of dust aerosol particle production, while soil properties are taken into account. Saltation processes are simulated following Marticorena and Bergametti (1995). The surface roughness length is fixed globally to a value of 0.001, and the minimum friction velocity threshold for dust mobilization is set to 21 cm s−1. The threshold friction velocity depends on the soil size distribution, vegetation cover and soil moisture (Cheng et al., 2008). The preferential dust emission sources include arid or low vegetated areas and are predefined in accordance with Tegen et al. (2002). Dust particles are initially emitted in the insoluble accumulation and coarse modes, but aging processes like the condensation of sulfuric acid (H2SO4) onto insoluble dust particles lead to soluble dust in the accumulation and coarse modes (Zhang et al., 2012).

The emissions for the remaining aerosol types and aerosol precursors are defined using emission inventories organized into 14 sectors, with each sector corresponding to the emission flux of one or more aerosol types or aerosol precursors (Schultz et al., 2018; Tegen et al., 2019). The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) dataset has been used for the anthropogenic aerosol and aerosol precursor emissions, which consists of monthly mean estimates at a horizontal resolution 0.5× 0.5 (Lamarque et al., 2010). The first version of Global Fire Assimilation System (GFAS) was employed for the representation of biomass burning emissions coming from grass and forest fires. GFAS is a gridded daily product with 0.5× 0.5 horizontal resolution based on the fire radiative power measurements of the MODIS instrument (Kaiser et al., 2012). The daily temporal resolution of GFAS and its accurate spatial representation of fires are the critical characteristics that make it ideal for the daily assimilation cycle applied in this study. Using standard GFAS emissions, several studies reported underestimated AOD; thus, a fire emission factor equal to 3.4 has been proposed (Kaiser et al., 2012; Tegen et al., 2019; Veira et al., 2015). This factor only corrects for the AOD bias and not the AAOD bias, and therefore the GFAS emission rescaling factor of 3.4 has not been adopted in the present study.

Most species are emitted at the lowest level of the model which represents the surface of the Earth. Aerosols related to energy production and ships are emitted directly to the second lowest level of the model. In HAM2.3, 75 % of biomass burning emissions are equally distributed in the planetary boundary layer (PBL), while 17 % and 8 % are emitted in the first and second levels above the PBL, respectively. Val Martin et al. (2010) reported that most biomass burning (BB) emissions occur within the PBL. However, it is noted that the wildfire emission height has a limited effect on global AOD distribution when compared to emission fluxes and dry or wet deposition processes (Veira et al., 2015). The partitioning of aerosol emission to the M7 modes is described in detail in Schutgens and Stier (2014). It is noted that all the experiments of this study use rescaled emissions based on the analysis of Sect. 3.4.

3 Methods

3.1 Local ensemble transform Kalman filter

In this study, we use the LETKF. The Kalman equation (Rodgers, 2000) involves the forecast, also called a priori or background state (xb) of the system, the analysis, also called a posteriori or assimilated state (xa) of the system, the observational data (y), the observational operator (H) and the two error covariance matrices that describe the uncertainties and correlations in the background estimates (P) and the observation estimates (R):

(1) x a = x b + G ( y - H x b ) .

G=PaHTR-1 is called the Kalman gain. The subscripts a and b indicate the analysis and background states, respectively, and T denotes the transpose operator. Equation (1) states that the analysis state (xa) is computed based on the background state (xb) plus the product of the Kalman gain and the difference between the observations and the simulated observations of the background state (y-Hxb) called innovation. Kalman gain is a matrix of weights, which corresponds to every ensemble member and adjusts the impact of innovation on the new analysis state. If Kalman gain is equal to zero, this means that either the model covariance is zero, the observational covariance error matrices is infinite, or there is no dependence of the measurements on the state vector elements, which implies that observations hold no useful information, and xa is set equal to xb. The analysis covariance error (Pα) can be calculated from the background state of the ensemble:

(2) P a = P b I + P b H T R - 1 H ,

where I is the identity matrix. Equations (1) and (2) can be also expressed as the minimization of the cost function:

(3) Ψ ( x a ) = ( x b - x a ) T P - 1 ( x b - x a ) + ( y - H x a ) T R - 1 ( y - H x a ) ,

where the two components on the right side of the equation describe the difference between background and analysis states and the difference between observations and analysis state. In a nutshell, the minimization of cost function (Eq. 3) defines a new and better estimate of the state vector based on observational data, a background estimate and taking into account the errors in both the model and observations.

The comprehensive mathematical formulation of LETKF can be found at Hunt et al. (2007). LETKF was implemented by Schutgens et al. (2010a) for an aerosol application based on previous work by Miyoshi and Yamane (2007). The code was modified to operate with ECHAM-HAM and POLDER in this paper.

Table 1State vector of the assimilation system composed by the mass mixing ratio of all the modes and species and number mixing ratio (PN) for all modes in ECHAM-HAM. Aerosols in nucleation mode are not considered in the assimilation since they contribute a very small fraction in the optical depth of aerosols.

Download Print Version | Download XLSX

The state vector of the system (x) is in our case composed of simulated ECHAM-HAM aerosol mass mixing ratios for every mode and species and number mixing ratio for every mode (23 in total; Table 1).

The observation vector y consists of aerosol optical properties retrieved by POLDER. The observation operator H transforms the aerosol mixing ratio into simulated aerosol optical properties and uses the optical property routines in ECHAM-HAM. The perturbed ensemble of ECHAM-HAM (Sect. 3.2) embodies the model error covariance matrix P, while the calculated POLDER errors (Sect. 2.1.2) are utilized in the observation error covariance matrix R. R assumes uncorrelated between the assimilated variables; hence, the off-diagonal elements of the matrix are set to zero. In reality, correlations do exist between variables, and thus they can affect the assimilation results. The off-diagonal elements of the R matrix can be estimated by constructing a data-derived R matrix (Liu et al., 2019), but it is out of the scope of the present study.

3.2 Model uncertainties and ensemble perturbation

When creating the ensemble to represent the model error covariance matrix, we consider two sources of uncertainty: namely in the aerosol emissions and in the wind speed/direction. Global climate models or climate transport models estimates of aerosol concentration and aerosol optical properties in the atmosphere are diverse and uncertain. The most prominent cause of these uncertainties mainly originates from the emission of natural aerosols; depending on the type, global estimates may differ by a factor of 4 to 16 (Grythe et al., 2014; Huneeus et al., 2011; Lewis and Schwartz, 2004; Miller et al., 2006; Pan et al., 2019; Textor et al., 2006). A recent study highlighted that the total global emissions of OC and BC differ by a factor of 4 based on six biomass burning emission inventories (Pan et al., 2019). Furthermore, a multi-model study from AeroCom phase I (Huneeus et al., 2011) and single-model studies using different DU emission schemes (Miller et al., 2006) indicate that DU global emission fluxes may differ by up to a factor of 6. SS global emission fluxes show the highest uncertainties among all other aerosol species mostly due to the differences in the simulated particles size (Textor et al., 2006) and the differences in the sea-spray-aerosol function of the models (Grythe et al., 2014). A wide range of SS fluxes has been reported depending on the sea-spray-aerosol function used (Grythe et al., 2014), but the most well-accepted range for SS emission flux is between 1.2 and 20 Pg yr−1 (Lewis and Schwartz, 2004), which implies a factor of approximately 16 SS emission flux difference between the highest and the lowest estimates.

Contrary to natural aerosol emissions, the anthropogenic aerosol emissions (OC, BC, SO4) and their precursors (SO2) are better constrained. For eastern China, the diversity (highest to lowest emission inventory) is lower than 1.35 for OC, BC and SO2 (Chang et al., 2015), although for other regions like South America diversity was estimated up to 3 (Granier et al., 2011). It is noted that most anthropogenic emission inventories are used on a monthly basis; hence, the day-to-day variability is not accounted in model simulations. It stands to reason that the uncertainty in daily anthropogenic emissions should be higher than the aforementioned values derived for monthly emissions. Other studies have noted that deposition parameterization causes large uncertainties in the direct and indirect aerosol radiative effects (Lee et al., 2016; Regayre et al., 2018) and hence may be considered in our future studies.

In addition to the emission uncertainties, meteorological factors may affect the emission, transport, deposition and hence the atmospheric lifetime and radiative effect of aerosol particles. In our simulations, the surface pressure vorticity and divergence of ECHAM are nudged towards the ERA-Interim reanalysis. Although ERA-Interim is a reanalysis product that provides a better estimate of the meteorological conditions, it contains uncertainties that may affect the life cycle of aerosol particles. Several studies using wind measurements from ground stations, radiosondes, buoys, cruises ships or even satellite estimations revealed ERA-Interim errors of up to ± 3 m s−1 (Bao and Zhang, 2013; Bromwich et al., 2016; Brunke et al., 2011; Campos and Guedes Soares, 2017; Stopa and Cheung, 2014).

Considering these sources of aerosol uncertainty, originating from aerosol emissions and wind speed/direction, our ECHAM-HAM forecast ensemble was assembled by multiplying the standard aerosol emissions by spatially correlated perturbations to obtain the emission for each member and by nudging each member to a slightly different version of ERA-Interim reanalysis in terms of wind speed.

A typical approach for the emission perturbation of the ensemble in past studies is to use global emission perturbations unique for each member and constant through time (Dai et al., 2014; Schutgens et al., 2010a). Another approach is to use spatiotemporal independent emission perturbations where random numbers are assigned for each grid cell, and often these numbers change in time. The latter case proved to produce very low spread in the ensemble, making the assimilation of observations impractical (Dai et al., 2014; Schutgens et al., 2010a). An intermediate approach is adopted in the present study by using spatially correlated emission perturbations, where the changes from grid to grid are not abrupt but smooth, perturbing emissions with positive or negative numbers over large areas. This technique has been successfully used in the past to derive soil erodibility factors under observation system simulation experiments (OSSEs) using an ensemble adjustment Kalman filter (Khade et al., 2013).

Spatially correlated perturbation fields are generated by first creating an ensemble (32 members in our standard setup) global grid at the resolution of the model (1.875× 1.875) filled with random values, sampled from a Gaussian distribution. This spatial field of perturbations does not contain any spatial correlation between the neighboring grid cells (the spatial field for one member is shown in Fig. S1a in the Supplement). In the next step, the global grid is smoothed by averaging, for each grid box, the surrounding grid boxes that are within a two-grid-box distance (Fig. S1b). The last step is repeated five times in total (Fig. S1b–f). Next, an exponential function is applied for each grid cell separately, making the values positive and the distribution positively skewed (close to lognormal) (Fig. S1g). The final step standardizes the numbers (V) to a mean (MVNEW) equal to the rescaled factors in Fig. S5 in the Supplement (for more information, see Sect. 3.4) and a standard deviation (SD) (SVNEW), equal to 0.65 for each grid (Fig. S1h), in the form of

(4) V NEW = SV NEW V - MV SV + MV NEW .

The variogram model for the spatial field of a member for each of the steps is shown in Fig. S2. According to it, in the final step (Fig. S2h), the variogram model flattens at the distance of 30, which indicates that grid cells are spatially correlated up to a distance of 30, whereas locations further than that are not.

The aforementioned methodology produces 32 different spatially correlated maps, one for each ensemble member. Each aerosol species has its own unique randomly generated set of 32 spatially correlated maps. The aerosol emission fluxes of the model are multiplied (perturbed) by these maps while the model runs. The SD of these spatially correlated perturbations (hence, the uncertainty of the emissions) for each species is equal to 0.65. This simplified approach assumes that the natural and anthropogenic emissions will have the same level of uncertainty.

Similar steps were followed for the creation of the zonal and meridional components of the wind spatially correlated perturbations. Contrary to the emission fluxes where numbers should be strictly positive, the wind vector sign indicates direction and the values can be negative; thus, the exponential function was not applied. The values were standardized with a mean and a SD of 0 and 0.8, respectively, using Eq. (4). The final wind perturbations are a selection of 32 members, different for the two wind components. These perturbations are added to the wind fields of the ERA-Interim reanalysis dataset, creating different perturbed branches that each member of the ensemble is nudged to. That approach can account partially for the uncertainties on the wind. A variety of other wind perturbation methods were tested, like altering the nudging relaxation time or restarting the model from different initial conditions, but were not adopted in the final experiments because either the produced ensemble spread was too small or it shrank close to zero after some days of simulation.

Figure 2An example of the spatially correlated perturbations maps generated for an ensemble member regarding (a) dust emissions and (b) the U component of the wind.

The distribution of perturbations for the dust emission and U component of the wind is shown in Fig. S3. An example of the emission and wind spatially correlated perturbation maps generated for a single member is shown in Fig. 2. In this example, the DU emissions of this member are higher than the default emission parameterization over the Arabian Peninsula and the central-eastern Sahara, while the western part of the Sahara will have lower DU emissions (Fig. 2a). Similarly, the U component (positive eastward, negative westward) of the nudging data for wind over the desert will be 1 m s−1 higher, which close to the surface may increase the emission of dust particles but higher up can also reduce the westward outflow of dust towards the North Atlantic (Fig. 2b). Other members may have a different combination of these parameters in the same area, which may cover other possible scenarios that may match reality.

Figure 3The uncertainty and the relative uncertainty of ECHAM-HAM for AOD550 (a, b) AOD865 (c, d) and AAOD550 (e, f) and the uncertainty for AE550−865 (g) and SSA550 (h) averaged over the period 20 July to 28 August 2006. The uncertainty is estimated using the SD of an ensemble (32 members) where each member used different spatial correlated perturbations on aerosol emissions and wind. The relative uncertainty is defined as the ratio of SD to mean. The global mean of each variable is denoted in the bottom-right corner for each case.

Figure 4On the left side, the grid-cell-scale collocation of the assimilated grid and the surrounding observations is illustrated. The highlighted blue area denotes the assimilated grid cell, the red circle around it represents the local patch size (Lx), and black dots represent the observations. It is noted that in LETKF the Lx distance is represented in grid cell units, but for illustrative purposes here it is represented as a circle. On the right side, the daily assimilation cycle is depicted for 3 d.


The derived model uncertainty, shown in Fig. 3 as the SD of this perturbed ensemble, can be conceptually compared with the POLDER uncertainty (Fig. 1). In all cases, the global mean of model uncertainty is lower than the POLDER uncertainty. This does not mean that the simulated aerosol observations of the model are more accurate in comparison to POLDER. The global mean uncertainty of the model is lower since low aerosol regions (majority of points globally) are far away from emission sources. Hence, the emission perturbation has a very limited effect and the members within the ensemble are quite similar to each other in the low aerosol regions. Moreover, the global mean of the model includes polar regions where the absolute uncertainty is very low. In the remote regions, POLDER microphysical retrievals are uncertain in an absolute sense and POLDER AOD retrievals are uncertain in a relative sense. On the other hand, the exact opposite is observed close to emission sources. For example, the SSA550 model uncertainty over the southern part of Africa can be up to 0.1 (Fig. 3h), while for the same region the POLDER uncertainty is only  0.04 (Fig. 1h). A notable difference between the model and POLDER uncertainties can be also observed for AOD550 in any outflow region over the ocean (e.g., North Atlantic Ocean, South Atlantic Ocean, South Indian Ocean) (Figs. 1a and 3a). Furthermore, it is noted also that other uncertainty factors related to aerosol physical and chemical processes (deposition, aging, vertical convective transport) are not taken into account in our model uncertainty estimation.

3.3 Data assimilation system

Our system consists of three phases, namely the spinup, the perturbation and the data assimilation phase. The first phase includes a single simulation that runs from January 2006 to the end of May 2006 and is a spinup simulation for the ECHAM-HAM meteorological and aerosol state. The second phase consists of an ensemble of runs from June to 28 August 2006, where each member's aerosol emission (DU, SS, OC, BC and SO4) and wind (U zonal and V meridional components) are perturbed using random spatially correlated maps (Sect. 3.2), which serves as the representation of model error in the assimilation phase. The third stage is the assimilation of observations, solved in a daily cycle from 20 July to 28 August 2006 (40 d).

The daily cycle of data assimilation involves daily forecasts (Dayt 00:00 UTC to Dayt+1 00:00 UTC) of all perturbed ensemble members. Upon completion of these simulations, the LETKF code is called, which performs a spatial collocation of the simulated (ECHAM-HAM) and the retrieved (POLDER) observations for four temporal time steps (00:00, 06:00, 12:00, 18:00 UTC). Subsequently, LETKF computes a new analysis state vector (ECHAM-HAM aerosol mixing ratio) on Dayt+1 at 00:00 UTC, which serves as the initial conditions for the next day's forecast. The process is repeated until the end of the data assimilation experiment.

At the grid cell scale, LETKF inflates the observation errors depending on their distance from the assimilated grid. This method, known in literature as “observation localization”, aims to reduce error covariance between distant points which is caused by sampling errors due to the limited ensemble size (Miyoshi and Yamane, 2007). More specifically, the local patch size (Lx) defines the distance between the analyzed grid cell and the observations that will be taken into account for assimilation. The observational error (E) of the observations that are within Lx is adjusted (EA) according to their distance (D) in grid cell units from the assimilated grid using a horizontal correlation length (Ly):

(5) E A = E exp D / L y 2 .

Distant observations get higher errors and thus have a smaller contribution in the changes of the assimilated grid. In our experiments, Lx and Ly are set to 4 and 2, respectively, in grid cell units. Consequently, observations that are up to four grids (7.5) away from the assimilated grid may affect the assimilated grid cell, although these more distant observations are accounted with a 2.7 times greater observational error than normal. An illustration of the daily assimilation cycle and the grid-cell-scale collocation during assimilation is depicted in Fig. 4.

3.4 Rescaling the aerosol emissions

Aerosol models often may under- or overestimate AOD550 or AAOD550 close to the sources. Most of these biases close to sources may be attributed to inaccurate emissions of one or two aerosol species. By rescaling these emissions by region and species, the model-simulated observations will be closer to the real observations. This rescaling of emissions can benefit the data assimilation since the Kalman filter assumes that the model is unbiased. Therefore, a simulation for 2006 was conducted with ECHAM-HAM in order to identify yearly AOD550 and AAOD550 biases by evaluating it against POLDER and MODIS-DT retrievals.

Initially, the mean yearly bias of the 2006 simulation against MODIS-DT AOD550, POLDER AOD550 and POLDER AAOD550, as well as the yearly emission fluxes for all species were plotted (Fig. S4). The results indicate that a positive bias of the model against MODIS-DT AOD550 is most probably driven by an overestimation of SS emission fluxes, since these two spatial patterns match (Fig. S4a and e). For the same reason, the over- or underestimation of the model when compared to POLDER AAOD550 over wildfire or anthropogenic polluted regions may be mainly attributed to BC emission fluxes (Fig. S4c and f). Furthermore, the biases of POLDER AOD550 over the desert are, for the most part, associated with DU emission fluxes (Fig. S4b and d). Emission rescaling factors (RFs) for SO2, SO4, SS and OC were based on biases against MODIS AOD550; for DU, they were based on biases against POLDER AOD550; and for BC, they were based on biases against POLDER AAOD550 (Fig. S5). The emission RFs were calculated as the ratio of observations (OBS) to model (MOD) AOD550 or AAOD550 for each region (r):

(6) RF r = OBS r MOD r .

Afterwards, the global grid was smoothed to avoid spatially steep changes of aerosol emission fluxes by averaging each grid cell with the surrounding values at a distance of two grid cells. This approach does not account for observation errors, it does not consider the interannual and intra-annual biases, and it is noted that the aerosol emission adjustment is based on AOD550 and AAOD550 biases, which may not be directly related to aerosol emissions.

An obvious limitation of this method occurs over the Southern Hemisphere major fire sources, where AOD550 is underestimated and AAOD550 is overestimated by ECHAM-HAM. Following the above-mentioned methodology, the resulting emission factors for fire source regions of the Southern Hemisphere would result in a reduction of BC emission fluxes by 10 %, due to the underestimation of AAOD550, and an increase of OC emission fluxes by more than 50 %, due to the underestimation of AOD550. The outcome would provide an improvement in terms of AOD550 but not for AAOD550, since part of the AAOD550 (10 %–20 %) emerges from other species, like OC. Thus, in order to get a simultaneous improvement in both AOD550 and AAOD550, OC emission fluxes were not adjusted in wildfire regions in the Southern Hemisphere major fire regions (Fig. S5c). The emission perturbation of the ensemble that is used in the core experiments is up to 3 times greater than the rescaling factors applied in this stage; thus, the high and low values of the original unadjusted yearly simulation are still represented in the perturbed ensemble.

3.5 Experimental setup

The experiments are focused on the summer of 2006, and assimilation is performed for the period of 20 July to 28 August 2006. The year was selected according to the availability of POLDER SRON retrievals, while the summer season was chosen due to the high peak of forest fires in the tropical band and the pronounced dust plume over the Atlantic, where AAOD and AE in the model can benefit from the assimilation process. The model simulations were bilinearly interpolated when necessary to a 1× 1 global grid for a direct comparison with the gridded satellite retrievals.

The core experiments presented in Sect. 4.1 assess the potential added value of assimilating aerosol information related to the size and absorption. In the CONTROL experiment, there was no assimilation of any kind of observations. In the MASS experiment, only AOD550 was assimilated. In the SIZE1 and SIZE2 experiments, either AOD550 and AOD865 or AOD550 and AE550−865 were assimilated, respectively. In the ABSORB1 and ABSORB2 experiments, either AOD550 and AAOD550 or AOD550 and SSA550 were assimilated, respectively. Finally, the TOTAL experiment assimilates AOD550 and AE550−865 and SSA550.

Table 2The name along with the assimilated aerosol optical properties and the physical meaning of each experiment. All the assimilated parameters are retrievals of POLDER. All of the experiments used the rescaled emission factors of Sect. 3.4.

Download Print Version | Download XLSX

Sensitivity experiments described in Sect. 4.3 intend to explore LETKF's sensitivity to three main parameters: ensemble size (nens), local patch size (Lx) along with horizontal correlation length (Ly) and inflation factor (ρ) following an analogous analysis conducted by Schutgens et al. (2010b). The nens is essentially the number of members used in ensemble and is connected to the accuracy and diversity of the model-predicted covariant error. The Lx represents the distance around a model grid that defines whether an observation would be considered in the assimilation process, while ρ relates to a technique that multiplies the error covariance matrix of the ensemble to increase the ensemble spread and ensure that the assimilation of new observations is possible in the next assimilation step (because otherwise the ensemble spread decreases during the assimilation process). Table 2 presents a summary of the core and the sensitivity experiments.

4 Results

4.1 Comparison to POLDER observations

We first evaluate the impact of data assimilation by evaluating the daily forecast, starting from the latest analysis, with POLDER data not yet assimilated. The benefit of this sort of evaluation is that biases in POLDER retrievals are effectively removed (because the observations used for either assimilation or evaluation come from the same POLDER dataset), and one can study the merits of data assimilation without the added issue of observational biases. The drawback is that such an evaluation cannot determine if the assimilation of POLDER retrievals actually yields improved aerosol simulation. The evaluation of the system with completely independent observations is presented in Sect. 4.2. Each experiment consists of an ensemble of simulations. The ensemble mean of each experiment is the best estimate for the state vector and the simulated observations. Thus, only the ensemble mean of each experiment is presented in the results. Throughout the results, we use the forecast run of the most recent analysis.

Figure 5Aerosol Optical Depth at 550 nm for (a) POLDER, the (b) CONTROL, (c) MASS, (d) SIZE2, (e) ABSORB2, (f) TOTAL experiments and their differences (model – POLDER; g–k). All fields are spatiotemporally collocated to the available measurements of POLDER for the period 20 July to 28 August 2006. Grey-filled grid cells indicate the absence of any valid POLDER measurements for the study period.

Figure 6Aerosol optical depth at 550 nm of POLDER against the core experiments, regarding all the available points of POLDER for the period 20 July to 28 August 2006. In each subplot, the bold line indicates the perfect model (y=x), while the two thinner lines confine the 200 % and 200 % bias boundaries, respectively. The shade depicts the density of points. N represents the total number of points.


4.1.1 Aerosol optical depth

Figures 5 and 6 show the AOD550 fields for the different data assimilation experiments of Table 2 and their agreement with the POLDER AOD550. Looking at the CONTROL experiment, we see that the model underestimates AOD550 over land in most regions with high AOD550, such as the Sahara (dust), tropical Africa (biomass burning) and Asia (industrial). On the other hand, over ocean the model has the tendency to overestimate AOD550. Overall, the global mean AOD550 of the model (0.146) is substantially lower than that of POLDER (0.228).

When assimilating AOD550 (MASS experiment), the global bias virtually disappears, which indicates the data assimilation system is capable of using the AOD550 information provided by POLDER measurements. When more properties than just AOD550 are assimilated (SIZE2, ABSORB2, TOTAL), the bias in AOD550 gets a bit larger than when only assimilating AOD550 (0.013 to 0.020) but is still much smaller than for the CONTROL experiment. The reason that the agreement with POLDER AOD550 gets slightly worse when assimilating other properties in addition to AOD550 (AE and/or SSA) is that the system needs to find the best compromise in fitting all properties simultaneously. So, in some situations the only way to get a better fit to AE is to degrade the fit to AOD550 (because of different assumptions used in the model).

Overall, local biases are decreasing after assimilation; however, over certain areas, biases can be increased, for example, over the South Atlantic Ocean in the MASS experiment (Fig. 5h). The assimilated AOD550 over Africa in the tropics increases the aerosol mixing ratio over land, which is then transported westward over the South Atlantic. The assimilation of AOD550 over the South Atlantic should compensate some of that effect by decreasing the aerosol mixing ratio but evidently not sufficiently. The assimilation of other aerosol optical properties like AE550−865 and SSA550 reduces the South Atlantic positive AOD550 bias, especially in the case of the TOTAL experiment (Fig. 5k), indicating that the simultaneous assimilation of multiple variables can improve the simulated AOD550 spatial representation in some places.

The AOD550 scatterplots for all core experiments are depicted Fig. 6. The averaged global mean error (ME) is reduced from 0.071 in the CONTROL to values that range between 0.002 to 0.020 in the assimilated experiments. Similarly, the mean absolute error (MAE) is reduced from 0.109 in the CONTROL to 0.078 in the TOTAL experiment. Pearson's correlation (R) increases from 0.668 to approximately 0.8 for all assimilated experiments. The consistent improvement of AOD550 in the assimilation experiments demonstrates that the combination of assimilated observations does not negatively affect AOD550.

Figure 7Similar to Fig. 5 but for the Ångström exponent at 550–865 nm.

Figure 8Similar to Fig. 6 but for the Ångström exponent at 550–865 nm. The two thinner lines confine the 30 % and 30 % bias boundaries, respectively.


4.1.2 Ångström exponent

Figures 7 and 8 depict AE550−865 for POLDER and the data assimilation experiments. The CONTROL experiment overestimates the AE550−865 in most cases (Pacific and Indian oceans, Australia, Siberia) and underestimates it in the western Sahara. Globally, the mean AE550−865 of CONTROL (1.25) is higher than that of POLDER (0.95), which indicates that the model overestimates the ratio of fine- to coarse-mode particles.

The MASS experiment has lower global AE550−865 (1.12), which matches slightly better but still remains higher than POLDER. In the MASS experiment, only AOD550 is assimilated; thus, any information regarding size or chemical composition will be related to a combination of transport and previous cycles that assimilated AOD550 close to sources. More specifically, particle size information may indirectly be introduced into the assimilation system by adjusting the AOD of two neighboring regions with different dominant particle size distributions, like northern Africa and Europe. For example, in Europe, the spatial averaged bias of AE550−865 in the CONTROL experiment was 0.20, while in MASS it is 0.01 (Fig. S6).

The other experiments (SIZE2, ABSORB2 and TOTAL) reduce the global mean of AE550−865 even more (1.06, 1.16 and 0.98, respectively). The combined information of more observations related to mass, size and absorption reduces the local biases of AE550−865 as well as the global ME and MAE while increasing R (Fig. 8g). Also, in the TOTAL experiment, the AE550−865 is improved substantially even in areas where the POLDER uncertainty of AE550−865 is quite high. For example, the POLDER uncertainty in Australia is approximately 0.5 (Fig. 1g), but CONTROL overestimates AE550−865 by 0.8 over land (Fig. 7g).

Figure 9Similar to Fig. 5 but for absorption optical depth at 550 nm.

Figure 10Similar to Fig. 5 but for single scattering albedo at 550 nm.

Figure 11Similar to Fig. 6 but for aerosol absorption optical depth at 550 nm.


4.1.3 Aerosol absorption

Figures 9 and 10 show the agreement in AAOD and SSA, respectively, between the different experiments and the POLDER observations. The CONTROL experiment shows lower global mean aerosol absorption than the POLDER observations; i.e., the AAOD is lower and the SSA is higher. Regionally, the differences are more complex. For example, over the southern part of Africa, with a lot of biomass burning, the AAOD is low compared to POLDER, and the SSA is as well. This means that the difference in AAOD is mostly caused by too-low total AOD but that the aerosols are more absorbing “per particle” in the CONTROL experiment than in POLDER. A similar pattern is observed over northern Eurasia. On the other hand, over the Middle East and over North America, a low bias in AAOD is observed together with a high bias is SSA, while for South America a high bias in AAOD is observed together with a low bias in SSA. So, in these regions, the differences between POLDER and the CONTROL experiment are caused by differences in absorption properties of the aerosols.

Figure 12Similar to Fig. 6 but for single scattering albedo at 550 nm. The two thinner lines show the 10 % and 10 % bias boundaries, respectively.


It is interesting to note that in the MASS experiment the AOD550 is improved (Fig. 5h) but SSA550 and AAOD550 not so much, especially in regions like South America, Africa and the Atlantic Ocean (Figs. 9h and 10h). The reason behind that is easiest to explain over South America, where AOD550 is underestimated and AAOD (SSA) is overestimated (underestimated) in CONTROL. The assimilation of AOD550 (MASS) will increase the aerosol mixing ratio of all aerosols based on their extinction, but it will not be based on their absorption. Thus, AAOD will be increased along with AOD since more aerosols will be in the atmosphere. Specifically, in the Amazon basin, SSA550 of the MASS experiment decreases by 0.032 in comparison to CONTROL, since the BC column burden becomes 4 times higher (Fig. S7b), while the difference of SSA550 between POLDER and the model (spatiotemporal collocated points only) increases from 0.084 to 0.117 (Fig. S7c).

Table 3Global column burden (Tg) of all aerosol species regarding the CONTROL experiment and the induced percentage changes due to assimilation for the MASS, SIZE2, ABSORB1 and TOTAL experiments.

Download Print Version | Download XLSX

Scatterplots of all experiments for AAOD550 and SSA550 are depicted in Figs. 11 and 12, respectively. The global SSA550 ME gets worse from 0.013 to 0.025 for the MASS experiments, the MAE increases also from 0.049 to 0.058, while R decreases from 0.243 to 0.162 (Fig. 12a and b). These results reveal the limitations of an aerosol data assimilation system that assimilates AOD in only one wavelength. Assimilating size information does not really improve the agreement in AAOD or SSA, but as expected, when assimilating SSA (ABSORB2 experiment), the agreement between model and POLDER data significantly improves for both AAOD and SSA. This improvement is maintained for the TOTAL experiment, assimilating AOD, AE and SSA together. The SSA550 of ABSORB2 experiment is still slightly higher over the Amazon basin and lower in the Middle East (Fig. 10j), but overall the simulated absorbing properties are significantly better in comparison to the CONTROL experiment (Fig. 10g) and much better in comparison to the MASS experiment (Fig. 10h).

4.1.4 Aerosol column burden changes

The aerosol column burden changes for each experiment are depicted in Table 3. All experiments reveal that ECHAM-HAM underestimates aerosol column burden globally; thus, the changes for all species are positive. In the MASS experiment, BC column burden is almost 5 times greater than the CONTROL experiment. BC AOD550 contribution to AOD550 is less than 10 % in most regions over the globe; thus, the assimilation of only AOD550 should not affect the BC column burden significantly. On the other hand, OC AOD550 contribution to AOD550 is between 50 % and 90 % in the tropical fire and outflow areas, and thus the assimilation of AOD550 is expected to significantly affect the OC column burden. Although BC and OC emissions are perturbed differently, correlations in these two species will still persist, since both BC and OC are emitted from the same location (but not with the same magnitude) and are following similar transport paths in each member. Thus, we conclude that the large BC column burden increase is related also to the correlations between OC and BC AOD550. The ABSORB1 experiment constrains the BC increase to +27 % in comparison to CONTROL. The TOTAL experiment increases of the column burden range between +20 % and +95 % for all the species. To understand better how these changes are made, we explore further regionally by isolating the effect of assimilating AE550−865 and AAOD550 for Australia and South America, respectively.

Figure 13An overview of the changes due to the addition of POLDER AE550−865 in the assimilation for Australia. Panel (a) depicts the mean values of five aerosol optical properties, while the height of bars indicates the regions relative amount of each property in comparison to the global mean. Panel (b) illustrates the column burden of five aerosol species as simulated in the MASS experiment. The percentage in the column burdens of each species specifies the relative contribution of each species in the total column burden. Panel (c) illustrates the absolute changes caused in aerosol optical properties, due to the assimilation of AE550−865 (SIZE2 – MASS), while the height of bars indicates if the change is just positive or negative. Panel (d) shows the relative changes caused in the mixing ratio for each species due to the assimilation of AE550−865 (SIZE2 – MASS). Panel (e) displays the aerosol optical properties' bias in comparison to POLDER for the MASS and SIZE2 experiments.

Figure 13 depicts the aerosol optical properties and aerosol column burden changes over Australia caused by the assimilation of AE550−865 by subtracting the value for the MASS experiment from SIZE2 (SIZE2 – MASS). Australia is a mixed aerosol area with fairly low aerosol content and many uncertain emission sources within the continent, while satellite retrievals are quite diverse (Schutgens et al., 2020) due to the complex surface albedo of the continent. MASS clearly overestimates the amount of fine particles over the Australian continent compared to POLDER (Fig. 7h); hence, the addition of AE550−865 in the assimilation increases the column burden on rather coarser aerosol groups (DU and SS) and decreases the column burden the modes corresponding to fine particles, where OC, BC and SO4 are dominant. (Fig. 13c). Consequently, the AE550−865 bias against POLDER is reduced from 0.282 in the MASS experiments to 0.166 in the SIZE2 experiment. The bias is decreased also in some other parameters, like AOD550, AAOD550 and SSA550. It is noted that the aerosol mixing ratio changes on the assimilation system are conducted for every aerosol tracer separately (Table 1), but since some aerosol groups contain higher column burden in coarser aerosol modes (DU and SS) and others in finer aerosol modes (OC, BC, SO4), aerosol adjustments due to assimilation are regularly consistent (either positive or negative) for an aerosol species as a whole, like in the case of Australia.

Figure 14Similar to Fig. 13 but for South America. The depicted experiments are MASS and ABSORB1; hence, the changes in panels (c) and (d) are based on the addition of AAOD550 to the assimilation.

Similarly, Fig. 14 depicts the aerosol optical properties and aerosol column burden changes over South America caused by the assimilation of AAOD550 by subtracting the MASS experiment from ABSORB1 (ABSORB1 – MASS). South America, and specifically the Amazon basin, is a major active burning area of the globe in July and August with significant emissions of absorbing aerosol particles. When assimilating AOD550 and AAOD550 (ABSORB1), the absorption optical properties are improved (Fig. 14c), while the AOD550, AOD865 and AE550−865 performance slightly deteriorates. The changes in the simulated absorbing properties are mainly driven by a 81 % reduction of the BC mixing ratio, which reduces the AAOD550 by 0.024 and increases the SSA550 by 0.078 in comparison to the MASS experiment (Fig. 14b). DU and OC changes may in the ABSORB1 experiment affect a small fraction of AAOD550 changes too, but predominantly the other species adjust to match the other assimilated parameter (AOD550) as well as possible.

4.2 Comparison with the independent observations

In this subsection, we evaluate the aerosol fields from the different data assimilation experiments using independent observations. AERONET is the most important data source for this purpose given its high accuracy, especially for AOD from its direct-Sun product. However, spatial coverage by AERONET sites is sparse and entirely absent over the ocean (except for a few islands and coastal stations). MODIS-DT and -DB on the other hand provide close-to-global coverage and have been extensively evaluated. Here, we present the comparisons of our POLDER assimilation experiments with these independent datasets.

Figure 15Comparison between the AERONET (direct-Sun V3 L2) and the CONTROL, TOTAL and POLDER experiments for AOD550 (a–c) and AE550−865 (d–f).


4.2.1 Evaluation with AERONET

In Fig. 15, the background simulation (CONTROL) and the total aerosol assimilated experiment (TOTAL) are evaluated against AERONET direct-Sun V3 L2 dataset for AOD and AE. The background simulation (CONTROL) shows a clear negative bias with AERONET in AOD of 0.072. Assimilating POLDER AOD, AE and SSA simultaneously (TOTAL) removes this bias almost entirely (remaining bias of 0.001). The reduction in MAE is much more moderate (from 0.127 to 0.11). In terms of AE550−865, the ME is reduced from 0.273 to 0.084, and the MAE is reduced from 0.419 to 0.353 (Fig. 15d and e). The spatiotemporal collocated points between POLDER and AERONET indicate that POLDER has a very small ME in AOD550 and AE550−865 (Fig. 15c and f); thus, the assimilation of these variables from POLDER will converge the ensemble mean to a very low bias in comparison to AERONET (Fig. 15b and e).

Figure 16Comparison between the AERONET (aerosol inversion V3 L2) and the CONTROL, MASS, ABSORB2 and TOTAL experiments for AAOD550 (a–d) and SSA550 (e–h).


In Fig. 16, the CONTROL, MASS, ABSORB2 and TOTAL experiments are evaluated using the AERONET aerosol inversion V3 L2 dataset for AAOD550 and SSA550. Both properties deteriorate in the MASS experiment and improve in the ABSORB2 experiment (Fig. 16b and e), while in TOTAL experiment only AAOD550 improves (Fig. 16c and f). It is noted that the AERONET aerosol inversion dataset that provides AAOD550 and SSA550 contains far fewer stations and hence less collocated points with the model (N= 772), in comparison to the AERONET direct-Sun dataset (N= 11 832). The spatiotemporal collocated points of POLDER retrievals with the AERONET inversion V3 L2 dataset were very few (N= 31) for the study period and thus are not shown.

Figure 17The AOD550 of MODIS-DT over land against the CONTROL, TOTAL and POLDER experiments (a–c). Similar setup for MODIS-DB AOD550 (d–f).


Figure 18The AOD550 of MODIS-DT over ocean against the CONTROL, TOTAL and POLDER experiments (a–c).


Figure 19ME and MAE for all the spatiotemporal collocated points between POLDER (assimilated observations) and the experiments, averaged over the whole study period.


Figure 20ME and MAE for all the spatiotemporal collocated points between AERONET (independent observations) and the experiments. AOD550, AOD865 and AE550–865 are calculated using the AERONET direct-Sun V3 L2 dataset, while AAOD550 and SSA550 are calculated using the AERONET aerosol inversion V3 L2 dataset.


4.2.2 Comparison with MODIS-DT and MODIS-DB

The AOD550 of the CONTROL and the TOTAL experiments is compared with MODIS-DT and MODIS-DB over land in Fig. 17. In both cases, the negative ME is reduced in the assimilated experiment from 0.066 to 0.002 (MODIS-DT over land) and 0.103 to 0.029 (MODIS-DB); the MAE is reduced as well and the correlation increases. The assimilation of POLDER observations brings the analysis closer to MODIS-DT and MODIS-DB over land, although the three datasets use different retrieval algorithms. In contrast, the comparison with MODIS-DT over ocean reveals that the assimilation slightly increases the ME and MAE (Fig. 18a and b). The reason for this is that the assimilation of POLDER measurements increases the AOD550 over land, which may have an effect on the AOD550 error over ocean in outflow regions (e.g., South Atlantic; Fig. 5k). These results indicate that more observations are needed over the source area of the outflow region (Africa) or more observations over the outflow region (South Atlantic) in order to constrain the ocean AOD550. Furthermore, the over-ocean AOD550 overestimation of the assimilated experiment is more prominent against MODIS-DT, since POLDER has a higher global mean AOD (by 0.032) than MODIS-DT (Fig. 18c).

4.3 Sensitivity experiments

In practice, many assimilation systems assume uncorrelated observational uncertainties by setting the off-diagonal elements of the observation error covariance matrix R to zero. The SIZE2 experiment results in a somewhat better agreement in AE with POLDER than SIZE1 (Fig. 8c and d), despite the fact that the observations contain the same information. This is most likely an artifact of existing correlations in the observation uncertainties. The nature of error correlations for these POLDER variables is illustrated in Fig. S8. Pearson's correlation for the first group of variables (AOD550, AOD865) is 0.92, while for the second group of variables (AOD550, AE550−865) it is 0.22. When assimilating AOD at two different wavelengths, it becomes important to specify the off-diagonal elements in the R matrix. Ignoring this prevents the LETKF from optimally using the information in the observations.

Similar conclusions can be drawn by the ABSORB1 (assimilation of AOD550 and AAOD550) and ABSORB2 (assimilation of AOD550 and SSA550) experiments. Both experiments are improving the simulated absorbing optical properties of the model (Figs. 11e and f and 12e and f). From the results, it seems that AAOD550 is more efficient in reducing the difference to POLDER in high-aerosol-content situations, like in the Amazon basin, where SSA550 has a larger effect over remote areas with lower aerosol content, like the Pacific Ocean (Fig. S9). Thus, depending on the area of interest, especially in studies that may use a similar assimilation system with a regional climate model, different combinations of variables may be necessary to adequately adjust the simulated absorbing aerosol properties of the model. In our case, SSA550 seems to have more impact at a global scale in comparison to AAOD550. Globally, ABSORB2 has better ME and MAE than ABSORB1. The most likely explanation lies within the correlations in the observational uncertainty of the assimilated variables. In Fig. S10, the estimated POLDER errors show that the correlation between errors in AOD550 and AAOD550 (0.684) is higher than the correlation between errors in AOD550 and SSA550 (0.192), indicating that the latter combination of variables can provide better results in the data assimilation framework of this study.

4.4 LETKF sensitivity experiments

In this subsection, we investigate the sensitivity of the data assimilation system to the number of ensemble members, the localization scale and the inflation parameter. More ensemble members (higher nens) provide a more accurate and detailed description of the model uncertainty at the cost of computational resources required. Previous studies with successful assimilation experiments have used nens values that ranged between 12 and 80 members (Dai et al., 2019; Lin et al., 2008; Rubin et al., 2016; Schutgens et al., 2010b; Sekiyama et al., 2010; Di Tomaso et al., 2017). In all cases, doubling nens did not significantly improve the assimilation results (Schutgens et al., 2010b; Di Tomaso et al., 2017). Rubin et al. (2016) showed that by quadrupling the nens (20 to 80) the root mean square error (RMSE) was improved for most of the AERONET sites and especially for sites affected by spatially large plumes of aerosol events. Specifically, for the Sede Boker AERONET site located in southern Israel, the bias and RMSE were reduced by 50 % and 35 %, respectively. All of the core experiments in this study consist of nens= 32. Two additional sensitivity experiments where conducted with nens= 16 (SMALL) and nens= 64 (LARGE).

Figure 19 indicates that the 64-member-ensemble size experiment (LARGE) managed to get a bit closer to the assimilated observations in all variables. The system managed to use the new ensemble correlations to simultaneously match the three assimilated observations a bit better. In contrast, the 16-member-ensemble size experiment (SMALL) difference between the assimilated observations was higher in comparison to both the 32- and the 64-member-ensemble size experiments for the aforementioned reason. On the other hand, the comparison to AERONET in Fig. 20 does not reveal consistent improvements for any of the variables when increasing the ensemble size. Obviously, the downside of the 64-member-ensemble size experiment is that it had to use double the resources and take more than twice the time to complete in comparison to the 32-member-ensemble size experiment, while the improvements were limited and apparent only when compared to the assimilated observations.

In the same plots, the spatiotemporal collocation between the assimilated grid cell and the observations that affect it in LETKF is tested with the Lx and Ly factors. Theoretically, a larger ensemble size may be able to benefit from information coming from more distant observations (Miyoshi and Yamane, 2007; Schutgens et al., 2010b). Two additional experiments have been conducted, using similar ensemble sizes (32 and 64) to those in the previous sensitivity experiments but with higher values (Lx= 6 and Ly= 3) than the default (Lx= 4 and Ly= 2). By comparing the TOTAL and LOCAL1 experiments as well as LARGE and LOCAL2 (same ensemble size, different localization factors), we can assess the effect of Lx and Ly factors. In both cases, the evaluation against AERONET shows that the TOTAL and LARGE experiments (core experiments) are superior in terms of global ME and MAE, but the evaluation against POLDER shows the opposite. By comparing LOCAL1 and LOCAL2, we can conclude that higher localization factors can benefit from higher ensemble size based on the comparison with POLDER, but yet again the evaluation against AERONET has contradicting results depending on the variable.

The final experiments (INFLATE1 and INFLATE2) test the inflation parameter (ρ) of the LETKF, which is multiplied with the background covariance matrix and prevents the ensemble spread of becoming too small. In each assimilation cycle, the ensemble spread of the analysis decreases, since all the ensemble members are converging to the same assimilated observations. This can create a background uncertainty that may be unrepresentative of the real background uncertainty, which will lead to an assimilation failure (Schutgens et al., 2010b). Here, we test if our spatially correlated perturbation methodology, which was used to describe the model uncertainty, can keep the ensemble spread big enough for the assimilation of POLDER observations. INFLATE1 experiment uses ρ= 1, which basically disables the inflation feature, INFLATE2 uses ρ= 1.5, while the rest of the experiments are using ρ= 1.1. For a direct comparison of the inflation impact, INFLATE1 and INFLATE2 should be compared with TOTAL. Both against POLDER and AERONET, TOTAL is in most cases a bit better than INFLATE1, thus concluding that inflation (ρ= 1.1) is a positive feature for the current assimilation framework. Additionally, when comparing INFLATE2 with TOTAL, the INFLATE2 mean error for POLDER AOD550, AOD865 and ANG550–865 is slightly smaller than that in TOTAL, but in all other cases (variables, statistics, observations) TOTAL is better.

5 Discussion and conclusions

We have presented the development of the first assimilation system for the ECHAM-HAM global aerosol climate model and demonstrated successful assimilation of multiple aerosol optical property retrievals from the POLDER SRON algorithm using an ensemble Kalman filter. The assimilation system uses an ensemble of perturbed simulations to define model uncertainty. The ensemble is created by perturbing the emission fluxes of all aerosol species and wind field of each ensemble member with spatially correlated perturbations. The uncertainty of POLDER observations was defined by evaluating the satellite retrievals with AERONET. The forecast output based on the most recent analysis of all the experiments is compared with POLDER (not yet assimilated), AERONET and MODIS observations.

The experiment in which only POLDER AOD550 is assimilated demonstrates considerable improvement in AOD550 against all observations (POLDER, MODIS-DT over land, MODIS-DB, AERONET) except MODIS-DT over ocean. Furthermore, the AOD550 correction also improves the simulated size representation through the reduction of AE550−865 global ME by 0.13, against POLDER observations. In contrast, it is noted that in the same experiment, both the AAOD550 and SSA550 deteriorate in terms of global ME and MAE. In certain regions, AAOD550 was also dramatically overestimated over Africa, South America and the Atlantic Ocean. These results reveal that AOD-only assimilation may lead to large discrepancies of the simulated aerosol absorbing optical properties.

Several other experiments that assimilated a combination of AOD550 with AOD865, AE550–865, AAOD550 and SSA550 showed consistent improvement in the assimilated variables in comparison to the no-assimilation experiment. The experiment where AOD550, AE550–865 and SSA550 were assimilated simultaneously was the most promising. The difference between model fields and assimilated observations decreased for virtually all aerosol optical properties in comparison to the experiment where only AOD550 was assimilated. The evaluation against AERONET showed that for all variables (except SSA550), both the global ME and MAE were improved in comparison to the CONTROL experiment, demonstrating that our data assimilation system can successfully constrain the simulated aerosol burden, size and absorption properties simultaneously. Our results suggest that it is very important to consider including AE and SSA, or other properties related to aerosol size and absorption, in future operational assimilation applications and especially in reanalysis simulations. Otherwise, aerosol size and absorption may be misrepresented.

Sensitivity experiments on the type and combination of the assimilated observations have been conducted. Assimilating AOD and AE instead of two AODs in different wavelengths reduces AE bias more, while assimilating AOD and SSA instead of AOD and AAOD decreases the bias of SSA and AAOD more. These results are most likely related to the correlations of the assimilated variables that the data assimilation system does not account for. Furthermore, the LETKF sensitivity experiments indicated that there was only a limited effect when varying the ensemble size, localization and inflation. Keeping the values of these parameters within the reported range of past literature provides similar results, and it is safe to say that they do not significantly affect the assimilation performance.

This work concludes that it is crucial to assimilate AE and SSA along with AOD in order to accurately correct for the burden, size and absorption of aerosol particles. The assimilation of other observational retrievals like the effective particle radii (size), column number, the fraction of spheres and the refractive index could improve the aerosol representation in the model further. Accurate measurements of these properties are expected from the SPEXone instrument on the NASA PACE mission (Hasekamp et al., 2019a; Werdell et al., 2019). In addition, the assimilation of AE at different wavelengths could also be an interesting experiment since AE at high wavelengths (e.g., 865 nm) is sensitive to the fine-mode fraction and not the effective radius of aerosol particles, and vice versa for AE at low wavelengths (e.g., 440 nm) (Schuster et al., 2006).

Appendix A

POLDER retrievals with horizontal resolution 18 km× 18 km were spatiotemporally collocated with the AERONET V3 L1.5 aerosol inversion dataset within an hour for the period 2006 to 2009. The AERONET L1.5 data were used in order to acquire more collocated points between POLDER and AERONET. Undoubtedly, this choice may cause an overestimation of POLDER uncertainties for AAOD550 and SSA550, since AERONET L1.5 data include retrievals of aerosol absorbing optical properties in cases of AOD440< 0.4, which are less accurate (Lacagnina et al., 2015). The AOD and AAOD of AERONET have been converted to POLDER wavelengths (550 and 865 nm) using AE:

(A1) AE λ 1 - λ 2 = log ( AOD λ 1 / AOD λ 2 ) log ( λ 1 / λ 2 ) ,

while SSA was calculated by combining AOD and AAOD:


where λ1 and λ2 are the wavelengths in nm. Afterwards, the relative errors (POLDER – AERONET / POLDER) for AOD550, AOD865 and AAOD550 were plotted against POLDER AOD550, AOD865 and AAOD550, respectively (Fig. A1). In a similar fashion, the errors of AE550–865 and SSA550 (POLDER – AERONET) where plotted against POLDER AOD550. In each case, the XX' axis was partitioned in six bins and the SD of the errors was calculated for each bin. Lastly, the SD of the relative differences for each bin and the variables AOD550, AOD865 and AAOD550 was multiplied with the POLDER 1× 1 gridded dataset to represent the POLDER uncertainty. In the case of AE550−865 and SSA550, the SD of the absolute errors for each bin was set as POLDER uncertainty. This fairly simple representation of POLDER uncertainties may be imprecise for some specific areas (e.g., the high-albedo arid areas) or some remote areas (e.g., over ocean) but carries the advantage of being based on high-quality independent observations.

Figure A1POLDER AOD550, AOD865 and AAOD550 uncertainty estimation using the relative errors of POLDER – AERONET and POLDER AE550−865 and SSA550 uncertainty estimation using the differences of POLDER – AERONET. Each variable is partitioned in six bins along the XX' axis. The grey shade indicates the SD and the purple shade illustrates the distribution for each bin.


Appendix B

where sim and obs are the simulated and observed values, while N is the population.

Code and data availability

The data assimilation experiments and the POLDER observations are available from the authors upon request. The AERONET, MODIS-DT and MODIS-DB data are freely available.


The supplement related to this article is available online at:

Author contributions

AT design the experiments, with the help of NAJS and OPH, and carried them out. AT performed the analysis and prepared the manuscript with contributions from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.


We thank the principal investigators, co-investigators and their staff for establishing and maintaining the AERONET sites used in this investigation. The authors would like to thank Paul Snijder for providing the spatiotemporal collocated data between POLDER (18 km× 18 km) and AERONET. This work was carried out on the Dutch national e-infrastructure with the support of SURF Cooperative.

Financial support

This research has been supported by the Dutch Research Council (NWO) and Netherlands Space Office (NSO) (grant no. 2017.008). Athanasios Tsikerdekis is funded by a NWO/NSO project “AEROSOURCE: Estimation of Aerosol Emissions from Polarization Data” (grant no. ALWGO.2017.008).

Review statement

This paper was edited by Chul Han Song and reviewed by two anonymous referees.


Amiridis, V., Wandinger, U., Marinou, E., Giannakaki, E., Tsekeri, A., Basart, S., Kazadzis, S., Gkikas, A., Taylor, M., Baldasano, J., and Ansmann, A.: Optimizing CALIPSO Saharan dust retrievals, Atmos. Chem. Phys., 13, 12089–12106,, 2013. 

Bao, X. and Zhang, F.: Evaluation of NCEP–CFSR, NCEP–NCAR, ERA-Interim, and ERA-40 Reanalysis Datasets against Independent Sounding Observations over the Tibetan Plateau, J. Climate, 26, 206–214,, 2013. 

Benedetti, A., Morcrette, J.-J., Boucher, O., Dethof, A., Engelen, R. J., Fisher, M., Flentje, H., Huneeus, N., Jones, L., Kaiser, J. W., Kinne, S., Mangold, A., Razinger, M., Simmons, A. J., and Suttie, M.: Aerosol analysis and forecast in the European Centre for Medium-Range Weather Forecasts Integrated Forecast System: 2. Data assimilation, J. Geophys. Res., 114, D13205,, 2009. 

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. 

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S. K., Sherwood, S., Stevens, B., and Zhang, X. Y.: chap. Clouds and Aerosols, in: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V. and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2013. 

Bromwich, D. H., Wilson, A. B., Bai, L.-S., Moore, G. W. K., and Bauer, P.: A comparison of the regional Arctic System Reanalysis and the global ERA-Interim Reanalysis for the Arctic, Q. J. Roy. Meteor. Soc., 142, 644–658,, 2016. 

Brunke, M. A., Wang, Z., Zeng, X., Bosilovich, M., and Shie, C. L.: An assessment of the uncertainties in ocean surface turbulent fluxes in 11 reanalysis, satellite-derived, and combined global datasets, J. Climate, 24, 5469–5493,, 2011. 

Campos, R. M. and Guedes Soares, C.: Assessment of three wind reanalyses in the North Atlantic Ocean, J. Oper. Oceanogr., 10, 30–44,, 2017. 

Carslaw, K. S., Lee, L. A., Reddington, C. L., Pringle, K. J., Rap, A., Forster, P. M., Mann, G. W., Spracklen, D. V., Woodhouse, M. T., Regayre, L. A., and Pierce, J. R.: Large contribution of natural aerosols to uncertainty in indirect forcing, Nature, 503, 67–71,, 2013. 

Chang, W., Liao, H., Xin, J., Li, Z., Li, D., and Zhang, X.: Uncertainties in anthropogenic aerosol concentrations and direct radiative forcing induced by emission inventories in eastern China, Atmos. Res., 166, 129–140,, 2015. 

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. 

Chen, C., Dubovik, O., Henze, D. K., Chin, M., Lapyonok, T., Schuster, G. L., Ducos, F., Fuertes, D., Litvinov, P., Li, L., Lopatin, A., Hu, Q., and Torres, B.: Constraining global aerosol emissions using POLDER/PARASOL satellite remote sensing observations, Atmos. Chem. Phys., 19, 14585–14606,, 2019. 

Cheng, T., Peng, Y., Feichter, J., and Tegen, I.: An improvement on the dust emission scheme in the global aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 8, 1105–1117,, 2008. 

Dai, T., Schutgens, N. A. J., Goto, D., Shi, G., and Nakajima, T.: Improvement of aerosol optical properties modeling over Eastern Asia with MODIS AOD assimilation in a global non-hydrostatic icosahedral aerosol transport model, Environ. Pollut., 195, 319–329,, 2014. 

Dai, T., Cheng, Y., Suzuki, K., Goto, D., Kikuchi, M., Schutgens, N. A. J., Yoshida, M., Zhang, P., Husi, L., Shi, G., and Nakajima, T.: Hourly Aerosol Assimilation of Himawari-8 AOT Using the Four-Dimensional Local Ensemble Transform Kalman Filter, J. Adv. Model. Earth Sy., 11, 680–711,, 2019. 

Di Tomaso, E., Schutgens, N. A. J., Jorba, O., and Pérez García-Pando, C.: Assimilation of MODIS Dark Target and Deep Blue observations in the dust aerosol component of NMMB-MONARCH version 1.0, Geosci. Model Dev., 10, 1107–1129,, 2017. 

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., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D., Tanré, D., and Slutsker, I.: Variability of Absorption and Optical Properties of Key Aerosol Types Observed in Worldwide Locations, J. Atmos. Sci., 59, 590–608,<0590:VOAAOP>2.0.CO;2, 2002. 

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. 

Eck, T. F., Holben, B. N., Reid, J. S., Dubovik, O., Smirnov, A., O'Neill, N. T., Slutsker, I., and Kinne, S.: Wavelength dependence of the optical depth of biomass burning, urban, and desert dust aerosols, J. Geophys. Res.-Atmos., 104, 31333–31349,, 1999. 

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. 

Fan, J., Wang, Y., Rosenfeld, D., and Liu, X.: Review of Aerosol–Cloud Interactions: Mechanisms, Significance, and Challenges, J. Atmos. Sci., 73, 4221–4252,, 2016. 

Fu, G. and Hasekamp, O.: Retrieval of aerosol microphysical and optical properties over land using a multimode approach, Atmos. Meas. Tech., 11, 6627–6650,, 2018. 

Generoso, S., Bréon, F. M., Chevallier, F., Balkanski, Y., Schulz, M., and Bey, I.: Assimilation of POLDER aerosol optical thickness into the LMDz-INCA model: Implications for the Arctic aerosol burden, J. Geophys. Res.-Atmos., 112, 1–15,, 2007. 

Georgoulias, A. K., Tsikerdekis, A., Amiridis, V., Marinou, E., Benedetti, A., Zanis, P., Alexandri, G., Mona, L., Kourtidis, K. A., and Lelieveld, J.: A 3-D evaluation of the MACC reanalysis dust product over Europe, northern Africa and Middle East using CALIOP/CALIPSO dust satellite observations, Atmos. Chem. Phys., 18, 8601–8620,, 2018. 

Giles, D. M., Sinyuk, A., Sorokin, M. G., Schafer, J. S., Smirnov, A., Slutsker, I., Eck, T. F., Holben, B. N., Lewis, J. R., Campbell, J. R., Welton, E. J., Korkin, S. V., and Lyapustin, A. I.: Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements, Atmos. Meas. Tech., 12, 169–209,, 2019. 

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. 

Grythe, H., Ström, J., Krejci, R., Quinn, P., and Stohl, A.: A review of sea-spray aerosol source functions using a large global set of sea salt aerosol concentration measurements, Atmos. Chem. Phys., 14, 1277–1297,, 2014. 

Hasekamp, O. P.: Linearization of vector radiative transfer with respect to aerosol properties and its use in satellite remote sensing, J. Geophys. Res., 110, D04203,, 2005. 

Hasekamp, O. P. and Landgraf, J.: Linearization of vector radiative transfer by means of the forward-adjoint perturbation theory and its use in atmospheric remote sensing, in: Light Scattering Reviews 2, edited by: Kokhanovsky, A. A., Springer Berlin Heidelberg, Berlin, Heidelberg, 159–204, 2007. 

Hasekamp, O. P., Litvinov, P., and Butz, A.: Aerosol properties over the ocean from PARASOL multiangle photopolarimetric measurements, J. Geophys. Res., 116, D14204,, 2011. 

Hasekamp, O. P., Fu, G., Rusli, S. P., Wu, L., Di Noia, A., Brugh, J. aan de, Landgraf, J., Martijn Smit, J., Rietjens, J., and van Amerongen, A.: Aerosol measurements by SPEXone on the NASA PACE mission: expected retrieval capabilities, J. Quant. Spectrosc. Ra., 227, 170–184,, 2019a. 

Hasekamp, O. P., Gryspeerdt, E., and Quaas, J.: Analysis of polarimetric satellite measurements suggests stronger cooling due to aerosol-cloud interactions, Nat. Commun., 10, 5405,, 2019b. 

Heinold, B., Tegen, I., Schepanski, K., and Banks, J. R.: New developments in the representation of Saharan dust sources in the aerosol–climate model ECHAM6-HAM2, Geosci. Model Dev., 9, 765–777,, 2016. 

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. 

Holben, B. N., Tanré, D., Smirnov, A., Eck, T. F., Slutsker, I., Abuhassan, N., Newcomb, W. W., Schafer, J. S., Chatenet, B., Lavenu, F., Kaufman, Y. J., Vande Castle, J., Setzer, A., Markham, B., Clark, D., Frouin, R., Halthore, R., Karneli, A., O'Neill, N. T., Pietras, C., Pinker, R. T., Voss, K., and Zibordi, G.: An emerging ground-based aerosol climatology: Aerosol optical depth from AERONET, J. Geophys. Res.-Atmos., 106, 12067–12097,, 2001. 

Holben, B. N., Eck, T. F., Slutsker, I., Smirnov, A., Sinyuk, A., Schafer, J., Giles, D., and Dubovik, O.: Aeronet's Version 2.0 quality assurance criteria, edited by: Tsay, S.-C., Nakajima, T., Singh, R. P., and Sridharan, R., 64080Q,, 2006. 

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. 

Hunt, B. R., Kostelich, E. J., and Szunyogh, I.: Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter, Physica D, 230, 112–126,, 2007. 

Hyer, E. J., Reid, J. S., and Zhang, J.: An over-land aerosol optical depth data set for data assimilation by filtering, correction, and aggregation of MODIS Collection 5 optical depth retrievals, Atmos. Meas. Tech., 4, 379–408,, 2011. 

Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res., 113, D13103,, 2008. 

Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones, L., Morcrette, J.-J., Razinger, M., Schultz, M. G., Suttie, M., and van der Werf, G. R.: Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosciences, 9, 527–554,, 2012. 

Khade, V. M., Hansen, J. A., Reid, J. S., and Westphal, D. L.: Ensemble filter based estimation of spatially distributed parameters in a mesoscale dust model: experiments with simulated and real data, Atmos. Chem. Phys., 13, 3481–3500,, 2013. 

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. 

Koch, D. and Del Genio, A. D.: Black carbon semi-direct effects on cloud cover: review and synthesis, Atmos. Chem. Phys., 10, 7685–7696,, 2010. 

Lacagnina, C., Hasekamp, O. P., Bian, H., Curci, G., Myhre, G., Noije, T., Schulz, M., Skeie, R. B., Takemura, T., and Zhang, K.: Aerosol single-scattering albedo over the global oceans: Comparing PARASOL retrievals with AERONET, OMI, and AeroCom models estimates, J. Geophys. Res.-Atmos., 120, 9814–9836,, 2015. 

Lacagnina, C., Hasekamp, O. P., and Torres, O.: Direct radiative effect of aerosols based on PARASOL and OMI satellite observations, J. Geophys. Res.-Atmos., 122, 2366–2388,, 2017. 

Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039,, 2010. 

Lana, A., Bell, T. G., Simó, R., Vallina, S. M., Ballabrera-Poy, J., Kettle, A. J., Dachs, J., Bopp, L., Saltzman, E. S., Stefels, J., Johnson, J. E., and Liss, P. S.: An updated climatology of surface dimethlysulfide concentrations and emission fluxes in the global ocean, Global Biogeochem. Cy., 25, GB1004,, 2011. 

Lee, L. A., Reddington, C. L., and Carslaw, K. S.: On the relationship between aerosol model uncertainty and radiative forcing uncertainty, P. Natl. Acad. Sci. USA, 113, 5820–5827,, 2016. 

Lewis, R. and Schwartz, E.: Sea Salt Aerosol Production: Mechanisms, Methods, Measurements and Models – A Critical Review, American Geophysical Union, Washington, D.C., 299–344, 2004. 

Li, Z., Zhao, X., Kahn, R., Mishchenko, M., Remer, L., Lee, K.-H., Wang, M., Laszlo, I., Nakajima, T., and Maring, H.: Uncertainties in satellite remote sensing of aerosols and impact on monitoring its long-term trend: a review and perspective, Ann. Geophys., 27, 2755–2770,, 2009. 

Lin, C., Wang, Z., and Zhu, J.: An Ensemble Kalman Filter for severe dust storm data assimilation over China, Atmos. Chem. Phys., 8, 2975–2983,, 2008. 

Liu, Y. A., Li, Z., and Huang, M.: Towards a data-derived observation error covariance matrix for satellite measurements, Remote Sens.-Basel, 11, 1770,, 2019. 

Long, M. S., Keene, W. C., Kieber, D. J., Erickson, D. J., and Maring, H.: A sea-state based source function for size- and composition-resolved marine aerosol production, Atmos. Chem. Phys., 11, 1203–1216,, 2011. 

Marticorena, B. and Bergametti, G.: Modeling the atmospheric dust cycle: 1. Design of a soil-derived dust emission scheme, J. Geophys. Res., 100, 16415,, 1995. 

Miller, R. L., Cakmur, R. V., Perlwitz, J., Geogdzhayev, I. V., Ginoux, P., Koch, D., Kohfeld, K. E., Prigent, C., Ruedy, R., Schmidt, G. A., and Tegen, I.: Mineral dust aerosols in the NASA Goddard Institute for Space Sciences ModelE atmospheric general circulation model, J. Geophys. Res., 111, D06208,, 2006. 

Mishchenko, M. I. and Travis, L. D.: Satellite retrieval of aerosol properties over the ocean using measurements of reflected sunlight: Effect of instrumental errors and aerosol absorption, J. Geophys. Res.-Atmos., 102, 13543–13553,, 1997. 

Miyoshi, T. and Yamane, S.: Local Ensemble Transform Kalman Filtering with an AGCM at a T159/L48 Resolution, Mon. Weather Rev., 135, 3841–3861,, 2007. 

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. 

Nabat, P., Somot, S., Mallet, M., Sevault, F., Chiacchio, M., and Wild, M.: Direct and semi-direct aerosol radiative effect on the Mediterranean climate variability using a coupled regional climate system model, Clim. Dynam., 44, 1127–1155,, 2014. 

Nightingale, P. D., Malin, G., Law, C. S., Watson, A. J., Liss, P. S., Liddicoat, M. I., Boutin, J., and Upstill-Goddard, R. C.: In situ evaluation of air-sea gas exchange parameterizations using novel conservative and volatile tracers, Global Biogeochem. Cy., 14, 373–387,, 2000. 

Pan, X., Ichoku, C., Chin, M., Bian, H., Darmenov, A., Colarco, P., Ellison, L., Kucsera, T., da Silva, A., Wang, J., Oda, T., and Cui, G.: Six global biomass burning emission datasets: intercomparison and application in one global aerosol model, Atmos. Chem. Phys., 20, 969–994,, 2020. 

Pope, R. J., Marsham, J. H., Knippertz, P., Brooks, M. E., and Roberts, A. J.: Identifying errors in dust models from data assimilation, Geophys. Res. Lett., 43, 9270–9279,, 2016. 

Rast, S., Brokopf, R., Esch, M., Gayler, V., Kirchner, I., and Kornblueh, L.: User manual for ECHAM6, PhD, Max Planck Institute of Meteorology, Hamburg, 2015. 

Regayre, L. A., Johnson, J. S., Yoshioka, M., Pringle, K. J., Sexton, D. M. H., Booth, B. B. B., Lee, L. A., Bellouin, N., and Carslaw, K. S.: Aerosol and physical atmosphere model parameters are both important sources of uncertainty in aerosol ERF, Atmos. Chem. Phys., 18, 9975–10006,, 2018. 

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding, World Scientific Publishing Co. Pte. Ltd., Singapore, USA, London,, 2000. 

Rubin, J. I., Reid, J. S., Hansen, J. A., Anderson, J. L., Collins, N., Hoar, T. J., Hogan, T., Lynch, P., McLay, J., Reynolds, C. A., Sessions, W. R., Westphal, D. L., and Zhang, J.: Development of the Ensemble Navy Aerosol Analysis Prediction System (ENAAPS) and its application of the Data Assimilation Research Testbed (DART) in support of aerosol forecasting, Atmos. Chem. Phys., 16, 3927–3951,, 2016. 

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 se ts, and usage recommendations, J. Geophys. Res.-Atmos., 119, 13965–13989,, 2014. 

Schepers, D., aan de Brugh, J. M. J., Hahne, P., Butz, A., Hasekamp, O. P., and Landgraf, J.: LINTRAN v2.0: A linearised vector radiative transfer model for efficient simulation of satellite-born nadir-viewing reflection measurements of cloudy atmospheres, J. Quant. Spectrosc. Ra., 149, 347–359,, 2014. 

Schultz, M. G., Stadtler, S., Schröder, S., Taraborrelli, D., Franco, B., Krefting, J., Henrot, A., Ferrachat, S., Lohmann, U., Neubauer, D., Siegenthaler-Le Drian, C., Wahl, S., Kokkola, H., Kühn, T., Rast, S., Schmidt, H., Stier, P., Kinnison, D., Tyndall, G. S., Orlando, J. J., and Wespes, C.: The chemistry–climate model ECHAM6.3-HAM2.3-MOZ1.0, Geosci. Model Dev., 11, 1695–1723,, 2018. 

Schuster, G. L., Dubovik, O., and Holben, B. N.: Angstrom exponent and bimodal aerosol size distributions, J. Geophys. Res., 111, D07207,, 2006. 

Schutgens, N., Nakata, M., and Nakajima, T.: Estimating Aerosol Emissions by Assimilating Remote Sensing Observations into a Global Transport Model, Remote Sens.-Basel, 4, 3528–3543,, 2012. 

Schutgens, N., Sayer, A. M., Heckel, A., Hsu, C., Jethva, H., de Leeuw, G., Leonard, P. J. T., Levy, R. C., Lipponen, A., Lyapustin, A., North, P., Popp, T., Poulsen, C., Sawyer, V., Sogacheva, L., Thomas, G., Torres, O., Wang, Y., Kinne, S., Schulz, M., and Stier, P.: An AeroCom–AeroSat study: intercomparison of satellite AOD datasets for aerosol model evaluation, Atmos. Chem. Phys., 20, 12431–12457,, 2020. 

Schutgens, N. A. J. and Stier, P.: A pathway analysis of global aerosol processes, Atmos. Chem. Phys., 14, 11657–11686,, 2014. 

Schutgens, N. A. J., Miyoshi, T., Takemura, T., and Nakajima, T.: Applying an ensemble Kalman filter to the assimilation of AERONET observations in a global aerosol transport model, Atmos. Chem. Phys., 10, 2561–2576,, 2010a. 

Schutgens, N. A. J., Miyoshi, T., Takemura, T., and Nakajima, T.: Sensitivity tests for an ensemble Kalman filter for aerosol assimilation, Atmos. Chem. Phys., 10, 6583–6600,, 2010b. 

Sekiyama, T. T., Tanaka, T. Y., Shimizu, A., and Miyoshi, T.: Data assimilation of CALIPSO aerosol observations, Atmos. Chem. Phys., 10, 39–49,, 2010. 

Shi, Y., Zhang, J., Reid, J. S., Holben, B., Hyer, E. J., and Curtis, C.: An analysis of the collection 5 MODIS over-ocean aerosol optical depth product for its implication in aerosol assimilation, Atmos. Chem. Phys., 11, 557–565,, 2011. 

Sofiev, M., Soares, J., Prank, M., de Leeuw, G., and Kukkonen, J.: A regional-to-global model of emission and transport of sea salt particles in the atmosphere, J. Geophys. Res.-Atmos., 116, D21302,, 2011. 

Stap, F. A., Hasekamp, O. P., and Röckmann, T.: Sensitivity of PARASOL multi-angle photopolarimetric aerosol retrievals to cloud contamination, Atmos. Meas. Tech., 8, 1287–1301,, 2015. 

Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric component of the MPI-M earth system model: ECHAM6, J. Adv. Model. Earth Sy., 5, 146–172,, 2013. 

Stier, P., Feichter, J., Kinne, S., Kloster, S., Vignati, E., Wilson, J., Ganzeveld, L., Tegen, I., Werner, M., Balkanski, Y., Schulz, M., Boucher, O., Minikin, A., and Petzold, A.: The aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 5, 1125–1156,, 2005. 

Stopa, J. E. and Cheung, K. F.: Intercomparison of wind and wave data from the ECMWF Reanalysis Interim and the NCEP Climate Forecast System Reanalysis, Ocean Model., 75, 65–83,, 2014. 

Tegen, I., Harrison, S. P., Kohfeld, K., Prentice, I. C., Coe, M., and Heimann, M.: Impact of vegetation and preferential source areas on global dust aerosol: Results from a model study, J. Geophys. Res.-Atmos., 107, AAC 14-1–AAC 14-27,, 2002. 

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. 

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. 

Tsikerdekis, A., Zanis, P., Steiner, A. L., Solmon, F., Amiridis, V., Marinou, E., Katragkou, E., Karacostas, T., and Foret, G.: Impact of dust size parameterizations on aerosol burden and radiative forcing in RegCM4, Atmos. Chem. Phys., 17, 769–791,, 2017. 

Tsikerdekis, A., Zanis, P., Georgoulias, A. K., Alexandri, G., Katragkou, E., Karacostas, T., and Solmon, F.: Direct and semi-direct radiative effect of North African dust in present and future regional climate simulations, Clim. Dynam., 53, 4311–4336,, 2019. 

Val Martin, M., Logan, J. A., Kahn, R. A., Leung, F.-Y., Nelson, D. L., and Diner, D. J.: Smoke injection heights from fires in North America: analysis of 5 years of satellite observations, Atmos. Chem. Phys., 10, 1491–1510,, 2010. 

Veira, A., Kloster, S., Schutgens, N. A. J., and Kaiser, J. W.: Fire emission heights in the climate system – Part 2: Impact on transport, black carbon concentrations and radiation, Atmos. Chem. Phys., 15, 7173–7193,, 2015. 

Vignati, E., Wilson, J., and Stier, P.: M7: An efficient size-resolved aerosol microphysics module for large-scale aerosol transport models, J. Geophys. Res.-Atmos., 109, 1–17,, 2004. 

Weaver, C., da Silva, A., Chin, M., Ginoux, P., Dubovik, O., Flittner, D., Zia, A., Remer, L., Holben, B., and Gregg, W.: Direct Insertion of MODIS Radiances in a Global Aerosol Transport Model, J. Atmos. Sci., 64, 808–827,, 2007. 

Werdell, P. J., Behrenfeld, M. J., Bontempi, P. S., Boss, E., Cairns, B., Davis, G. T., Franz, B. A., Gliese, U. B., Gorman, E. T., Hasekamp, O., Knobelspiesse, K. D., Mannino, A., Martins, J. V., McClain, C. R., Meister, G., and Remer, L. A.: The Plankton, Aerosol, Cloud, Ocean Ecosystem Mission: Status, Science, Advances, B. Am. Meteorol. Soc., 100, 1775–1794,, 2019. 

Wu, L., Hasekamp, O., van Diedenhoven, B., and Cairns, B.: Aerosol retrieval from multiangle, multispectral photopolarimetric measurements: importance of spectral range and angular resolution, Atmos. Meas. Tech., 8, 2625–2638,, 2015. 

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. 

Yumimoto, K. and Takemura, T.: Direct radiative effect of aerosols estimated using ensemble-based data assimilation in a global aerosol climate model, Geophys. Res. Lett., 38, L21802,, 2011. 

Yumimoto, K., Uno, I., Sugimoto, N., Shimizu, A., and Satake, S.: Adjoint inverse modeling of dust emission and transport over East Asia, Geophys. Res. Lett., 34, 1–6,, 2007. 

Yumimoto, K., Nagao, T. M., Kikuchi, M., Sekiyama T. T., Murakami, H., Tanaka, T. Y., Ogi, A., Irie, H., Khatri, P., Okumura, H., Arai, K., Morino, I., Uchino, O., and Maki, T.: Aerosol data assimilation using data from Himawari-8, a next-generation geostationary meteorological satellite, Geophys. Res. Lett., 43, 5886–5894,, 2016. 

Yumimoto, K., Tanaka, T. Y., Yoshida, M., Kikuchi, M., Nagao, T. M., Murakami, H., and Maki, T.: Assimilation and Forecasting Experiment for Heavy Siberian Wildfire Smoke in May 2016 with Himawari-8 Aerosol Optical Thickness, J. Meteorol. Soc. Jpn. Ser. II, 96, 133–149,, 2018.  

Zhang, J. and Reid, J. S.: MODIS aerosol product analysis for data assimilation: Assessment of over-ocean level 2 aerosol optical thickness retrievals, J. Geophys. Res., 111, D22207,, 2006. 

Zhang, K., O'Donnell, D., Kazil, J., Stier, P., Kinne, S., Lohmann, U., Ferrachat, S., Croft, B., Quaas, J., Wan, H., Rast, S., and Feichter, J.: The global aerosol-climate model ECHAM-HAM, version 2: sensitivity to improvements in process representations, Atmos. Chem. Phys., 12, 8911–8949,, 2012. 

Short summary
Accurate representation of aerosols in the atmosphere is hard to achieve due to their complex microphysical and optical properties and uncertain emissions. In our work, we employ a data assimilation method which integrates model simulations with satellite observation related to the amount, size and the light absorption of aerosol. The use of these observations in an experiment improves aerosol representation and it is recommended for utilization in future data assimilation practices.
Final-revised paper