# Characterization of smoke and dust episode over West Africa: comparison of MERRA-2 modeling with multiwavelength Mie-Raman lidar observations

Igor Veselovskii, Philippe Goloub, Thierry Podvin, Didier Tanre, Arlindo da Silva, Peter Colarco, Patricia Castellanos, Mikhail Korenskiy, Qiaoyun Hu, David Whiteman, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Igor Veselovskii, Philippe Goloub, Thierry Podvin, Didier Tanre, Arlindo da Silva, et al. Characterization of smoke and dust episode over West Africa: comparison of MERRA-2 modeling with multiwavelength Mie-Raman lidar observations. Atmospheric Measurement Techniques, 2018, 11 (2), pp.949-969. 10.5194/amt-11-949-2018. hal-04290021

HAL Id: hal-04290021

https://ulco.hal.science/hal-04290021v1

Submitted on 20 Nov 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Characterization of smoke and dust episode over West Africa: comparison of MERRA-2 modeling with multiwavelength Mie–Raman lidar observations

Igor Veselovskii<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Philippe Goloub<fnref n="4" />, Thierry Podvin<fnref n="4" />, Didier Tanre<fnref n="4" />, Arlindo da Silva<fnref n="3" />, Peter Colarco<fnref n="3" />, Patricia Castellanos<fnref n="3" />,<fnref n="5" />, Mikhail Korenskiy<fnref n="1" />, Qiaoyun Hu<fnref n="3" />, David N. Whiteman<fnref n="4" />, Daniel Pérez-Ramírez<fnref n="6" />, Patrick Augustin<fnref n="7" />, Marc Fourmentin<fnref n="7" />, and Alexei Kolgotin<fnref n="1" />

**Correspondence:** Igor Veselovskii (iveselov@hotmail.com)

Received: 21 September 2017 – Discussion started: 5 October 2017

Revised: 13 January 2018 – Accepted: 16 January 2018 – Published: 16 February 2018

**Abstract.** Observations of multiwavelength Mie–Raman lidar taken during the SHADOW field campaign are used to analyze a smoke–dust episode over West Africa on 24–27 December 2015. For the case considered, the dust layer extended from the ground up to approximately 2000 m while the elevated smoke layer occurred in the 2500–4000 m range. The profiles of lidar measured backscattering, extinction coefficients, and depolarization ratios are compared with the vertical distribution of aerosol parameters provided by the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The MERRA-2 model simulated the correct location of the near-surface dust and elevated smoke layers. The values of modeled and observed aerosol extinction coefficients at both 355 and 532 nm are also rather close. In particular, for the episode reported, the mean value of difference between the measured and modeled extinction coefficients at 355 nm is 0.01 km⁻<fnref n="1" /> with SD of 0.042 km⁻<fnref n="1" />. The model predicts significant concentration of dust particles inside the elevated smoke layer, which is supported by an increased depolarization ratio of 15 % observed in the center of this layer. The modeled at 355 nm the lidar ratio of 65 sr in the near-surface dust layer is close to the observed value (70 ± 10) sr. At 532 nm, however, the simulated lidar ratio (about 40 sr) is lower than measurements

(55 ± 8 sr). The results presented demonstrate that the lidar and model data are complementary and the synergy of observations and models is a key to improve the aerosols characterization.

## 1 Introduction

Atmospheric aerosols are an important factor influencing the Earth's radiative budget, though its impact is still highly uncertain due largely to the complicated mechanisms of aerosol–cloud interaction. Aerosol particles serve as cloud condensation nuclei and ice-nucleating particles, creating a strong impact on cloud and precipitation formation. However, different aerosol types differ significantly in their ability to initiate drop and ice crystal nucleation. There is thus a clear need for a better knowledge on vertically resolved optical, physical, and chemical aerosol properties. Lidar is a recognized instrument for vertical profiling of aerosol properties, and the possibility to invert lidar observations at several wavelengths to aerosol microphysical properties has been extensively studied both theoretically and experimentally over the past two decades (e.g., Müller et al., 1999, 2016; Veselovskii et al., 2002; Böckman et al., 2005). These

<footnote n="1">¹Physics Instrumentation Center of GPI, Troitsk, Moscow, Russia</footnote>

<footnote n="2">²Joint Center for Earth Systems Technology, UMBC, Baltimore, USA</footnote>

<footnote n="3">³NASA Goddard Space Flight Center, Greenbelt, USA</footnote>

<footnote n="4">⁴Laboratoire d'Optique Atmosphérie, Université de Lille-CNRS, Villeneuve d'Ascq, France</footnote>

<footnote n="5">⁵Universities Space Research Association, Columbia, Maryland, USA</footnote>

<footnote n="6">⁶Applied Physics Department, University of Granada, Spain</footnote>

<footnote n="7">⁷Laboratoire de Physico-Chimie de l'Atmosphère, Université du Littoral Côte d'Opale, France</footnote>
studies revealed the importance of using Raman or HSRL (high spectral resolution lidar) systems, which allow independent measurements of aerosol extinction and backscattering coefficients to be made. At present, the most practical configuration of Raman (HSRL) lidar is based on a triple wavelength Nd:YAG laser. Such a lidar provides the so-called $3\beta + 2\alpha$ set of observations, including three backscattering (355, 532, 1064 nm) and two extinction (355, 532 nm) coefficients.

However, the problem of inversion of $3\beta + 2\alpha$ observations is underdetermined (Chemyakin et al., 2016; Alexandrov and Mishchenko, 2017). As a result, instead of a unique solution, a family of solutions should be considered, leading to an increase in retrieval uncertainties. Still the estimation of volume density (V) and effective radius ($r_{eff}$) with uncertainty below 30% is possible, especially when the fine mode in the particle size distribution (PSD) is predominant (e.g., Veselovskii et al., 2004; Müller et al., 2005, 2016; Pérez-Ramírez et al., 2013). The refractive index (RI) can be also estimated from the measurements, although the uncertainty of such estimation is significant: for the real part ($m_R$) of RI the uncertainty is normally about $\pm$ 0.05 and for the imaginary part ($m_I$) it is about 50% when $m_I > 0.01$ (Veselovskii et al., 2004; Müller et al., 2016). Proposed improvements of inversion schemes were considered in recent publications (Chemyakin et al., 2014; Kolgotin et al., 2016), but these improvements are not able to resolve the fundamental issue: the information content of $3\beta + 2\alpha$ observations is insufficient to support exact solution of the problem and additional information should be used in retrievals to improve the accuracy of the retrieved products (Veselovskii et al., 2005; Burton et al., 2016; Alexandrov and Mishchenko, 2017; Kahnert and Andersson, 2017).

We should recall also that in the inversion schemes considered, the RI is normally assumed to be spectrally and size independent, which is generally not the case in the atmosphere. The irregularity of the particles shape can be also a significant error source. Moreover, the volume density and effective radius obtained from $3\beta + 2\alpha$ observations are attributed to the whole size distribution, which is of limited practical use because of the importance of characterizing the particle properties separately for the fine and coarse modes. Considering these issues makes the inverse problem even more underdetermined, emphasizing the need for additional input information.

One opportunity to get this additional information is by combining the lidar observations with aerosol transport models (Kahnert and Andersson, 2017). Models provide the vertical distribution of mass mixing ratios of chemical aerosol components, which can be used as “initial guess” in the inversion scheme. Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), offers a unique opportunity to provide such an “initial guesses” of the vertical structure of aerosol chemical composition. MERRA-2 is produced with NASA’s global Earth system

model, GEOS-5 (Goddard Earth Observing System version 5) (Gelaro et al., 2017) and includes an online coupling with the Goddard Chemistry, Aerosol, Radiation and Transport model (GOCART), which allows for assimilation of aerosol optical depth (AOD) from spaceborne and surface instruments such as MODIS, AVHRR, MISR, and AERONET (Randles et al., 2017). The fundamental data that MERRA-2 provides are vertical profiles of the mass mixing ratios of five aerosol components: dust, sea salt (SS), black and organic carbon (BC and OC), and sulfate (SU) aerosols. The main optical parameters related to lidar measurements, such as aerosol extinction and backscattering coefficients can be calculated basing on these data. The principal question arising, however, is how well the reanalysis reproduces independent observations and thus can provide a realistic initial guess for a lidar inversion scheme. Buchard et al. (2017) and Randles et al. (2017) extensively validated MERRA-2 with independent surface and aircraft observations of particulate matter (PM<sub>2.5</sub>) and AOD, as well as space-based observations of absorption aerosol optical depth and aerosol index. The extinction profiles derived from airborne HSRL measurements were also compared with modeling, finding generally good agreement between the observations and MERRA-2.

For global validation of the aerosol vertical distribution, the modeled profiles of attenuated backscatter were compared to spaceborne Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) observations (Winker et al., 2009), and a good consistency between simulations and observations was reported (Nowotnick et al., 2015; Buchard et al., 2017). Additional opportunities for model validation are provided by ground-based multiwavelength Raman or HSRL systems. Such lidars by their nature have limited spatial coverage but are well suited for the characterization of the vertical distribution of particle properties at a chosen location.

In our paper, we consider Raman lidar observations taken during a smoke–dust episode over West Africa in December 2015 during the SHADOW (SaHAran Dust Over West Africa) campaign (Veselovskii et al., 2016) and compare the vertical profiles of particle parameters with MERRA-2. The simultaneous presence of dust and smoke layers in the atmosphere provides an opportunity to test the ability of the model to reproduce the vertical structure of aerosol properties over the observation site.

## 2 Measurement setup and data analysis

### 2.1 Observation site

The observation site is located at the Institute for Research and Development Center (IRDC), Mbour, Senegal (14° N, 17° W). Information about the SHADOW campaign and instruments at the IRDC site can be found in the recent publication by Veselovskii et al. (2016). During the SHADOW campaign data from three lidar instruments were available:
organic carbon is based on the 100% brown carbon case from Hammer et al. (2016) and it is implemented as described in Colarco et al. (2017).

The sources of aerosols in the model include wind-speed-
based emissions of dust and sea salt, fossil fuel combustion,
biomass burning, biofuel consumption, biogenic particulate
organic matter, and oxidation of dimethyl sulfide and SO₂,
which includes volcanic sources. Aerosol sinks include con-
vective scavenging, dry deposition, and wet removal, where
aerosol hygroscopic growth is considered in the calculation
of particle fall velocity and deposition velocity. The model
resolution is 0.5° × 0.625° latitude by longitude with 72
hybrid-eta layers from the surface to 0.01 hPa. Additional de-
tails of the simulation can be found in Randles et al. (2017)
and Buchard et al. (2017).

In MERRA-2, aerosol and meteorological observations
are jointly assimilated within GEOS-5. Aerosols are as-
similated by means of analysis splitting and the local dis-
placement ensemble methodology (Buchard et al., 2015,
2016). The system assimilated MODIS, AVHRR, MISR, and
AERONET 550 nm AOD. AERONET measurements are in-
terpolated to 550 nm using the Ångström relationship and the
closest available channels, generally 500 and 675 nm. The
assimilation determines an AOD increment, which corrects
the model AOD in a way that minimizes the differences be-
tween the model and observations. The AOD increment both
corrects for misplaced aerosol plumes and scales the aerosol
mass mixing ratio to match the observations. The 2-D AOD
increment does not contain enough information to correct ei-
ther the vertical distribution of aerosols or the aerosol com-
position. Thus, the model determines the aerosol speciation,
optical properties, and vertical distributions, while the AOD
increments modulate the aerosol mass. Thus, the assimilated
aerosol distributions and physical and optical properties arise
from the forecast model assumptions and the formulation of
the aerosol data assimilation algorithm.

3 Experimental results

The smoke layers from forest fires near the Equator were regularly observed over the instrumentation site during the wintertime measurement sessions made in December 2015–January 2016. In our study we will focus on a strong smoke episode that occurred on 24–27 December 2015. Air mass back trajectories over Mbour on 25 December 2015 at 04:00 UTC are shown in Fig. 1 together with map of fires on 20 December 2015 (https://worldview.earthdata.nasa.gov).

The air masses below 1000 m (red line) are transported
over the desert and are strongly loaded by dust, while air
masses at 3000 m (green line) arrive from the south and pass
over the regions of forest fires and thus can transport smoke
particles. The Cimel MPL operated continuously through the
period of 24–27 December and thus monitored the arrival and
evolution of the smoke layer, as shown in Fig. 2. An elevated

smoke layer appears on 24 December around 00:00 UTC. The aerosol layer becomes thicker during the day but remains confined to the height interval of 2.5–4.0 km and stays well separated from the dust layer, which extends from the ground to approximately 2.0 km. This structure of the layers is preserved throughout 25 December, as well. Cirrus clouds appear at 08:00 UTC on 24 December at a height of 8 km, soon after the smoke layer arrival (Fig. 2a) and persist throughout the smoke episode. After 12:00 UTC on 24 December the clouds start descending and by 07:00 UTC on 25 December the cloud base is below 6 km (Fig. 2b). On 26 December strong precipitation of ice particles occurs (Fig. 2c) and, finally, on 27 December the cloud is located at the top of the smoke–dust layer (Fig. 2d).

Multiwavelength Raman lidar observations are available for the 23–25 December period only. The height–temporal evolution of the particle backscattering coefficient β<sub>532</sub>, depolarization ratio δ<sub>532</sub>, and water vapor mixing ratio w measured by Raman lidar on the nights 23–24 and 24–25 December 2015 are shown in Fig. 3. Due to the geometrical overlap factor the extinction data can be processed starting from approximately 750 m, and thus plots of all parameters start at this height. The depolarization ratios of pure dust observed during SHADOW are in the 30–35 % range (Veselovskii et al., 2016), while the depolarization ratio of smoke at 532 nm normally is below 10 % (e.g., Tesche et al., 2011; Burton et al., 2015). Hence depolarization measurements provide a convenient way to separate the aerosols into dust and smoke components.

On the night of 23–24 December the dust layer extends up to 2500 m, but with a high depolarization ratio (> 30 %), which is usually associated with pure dust, is observed only below 1000 m, meaning that in 1000–2500 m range the dust is probably mixed with smoke. The optical depth of the elevated smoke layer is rather small on 23–24 December (0.1 at 05:00 UTC), but on 24–25 December it increases up to 0.25, making possible the calculation of extinction coefficients from the Raman lidar signals. For analyzing the vertical distribution of smoke and dust particle parameters, we focus on the nighttime measurements of 24–25 December 2015.

Figure 4 shows the horizontal wind direction and speed measured by the wind lidar on 24–25 December. The range-corrected signal of the wind lidar can be evaluated starting from 100 m height, and the corresponding height–temporal image is shown in Fig. 5. The wind speed was measured in the dust layer (< 1500 m) for the whole period, but inside the smoke layer the backscatter signal is lower, so the measurements were possible only in the period of 16:00–22:00 UTC on 24 December. During 24–25 December 2015, the wind in the low troposphere (< 1500 m) is mainly dominated by the easterly Harmattan continental trades. Deceleration and acceleration of the lower part of the Harmattan (< 1000 m) are observed, respectively, in the beginning of the afternoon and during the night. The vertical profile of the wind speed
NOAA HYSPLIT MODEL
Backward trajectories ending at 04:00 UTC 25 Dec 15
GDAS meteorological data

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Five-day backward trajectories for the air mass in Mbour at altitudes 750, 1500, and 3500 m, on 25 December 2015 at 04:00 UTC, together with the map of forest fires on 20 December 2015.</figcaption></figure>

demonstrates the presence of a low-level jet (LLJ), where the maximum wind speed (jet speed) is located at a height of 350 m (LLJ height) at 01:00 UTC. LLJs are known to contribute to regional horizontal aerosol transport and to increase vertical mixing. Indeed, the LLJ occurrence at 01:00 UTC increases the aerosol loading by transporting desert dust. The corresponding increase of backscattering due to the LLJ at 01:00 UTC on 25 December can also be seen in Fig. 3.

The vertical profiles of temperature *T*, potential temperature *Θ*, wind direction and speed, together with relative humidity (RH) and WVMR from radiosonde launched from Dakar at 00:00 UTC on 25 December 2015, are shown in Fig. 6. The profile of wind speed and wind direction obtained from the sonde confirms that the LLJ observed with lidar at Mbour is not a local phenomenon, because it is also observed at Dakar. The vertical profile of the potential temperature suggests that the nocturnal boundary layer top corresponds to the LLJ height. Above 3000 m, the lidar and sonde depict southerly winds which transport the smoke plume. The water vapor mixing ratio increases above 2500 m; as a result the RH in the smoke layer reaches 75 % while in the dust layer RH is below 30 %.

To quantify the vertical distribution of particle parameters, Fig. 7 shows the profiles of backscattering ($\beta_{355}$, $\beta_{532}$, $\beta_{1064}$), extinction ($\alpha_{355}$, $\alpha_{532}$) coefficients and the particle depolarization ratio ($\delta_{532}$) derived from Raman lidar measurements for three temporal intervals on the night of 24–25 December: 19:00–23:00, 01:00–04:00, and 04:00–07:00 UTC. For the profiles presented, the uncertainty of both $\beta$ and $\alpha$ computations is estimated to be below 10 % for the Raman technique and to be below 20 % for $\beta_{1064}$ computation by the Klett method. The relative uncertainty of depolarization measurements is below 15 %. The extinction and backscattering Ångström exponents $A_{355/532}^{\alpha}$, $A_{355/532}^{\beta}$, and $A_{532/1064}^{\beta}$ are given by Fig. 8.

For the first temporal interval (Fig. 7a) dust and smoke layers are well separated. Extinction coefficients $\alpha_{355}$ and $\alpha_{532}$ differ in the smoke layer ($\alpha_{355} > \alpha_{532}$), but inside the near-surface dust layer (below 1750 m) the extinction values are nearly the same. The depolarization ratio is $\delta_{532} = 35 \pm 5 \%$ at 750 m and it gradually decreases with height to $27 \pm 4 \%$ at 1750 m. Above that height $\delta_{532}$ decreases quickly, indicating an increase in the contribution of smoke particles. For the second and third temporal intervals the dust and smoke lay-
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Range-corrected lidar signal (in arbitrary units) of Cimel MPL for 24–27 December 2015.</figcaption></figure>

ers appear to mix, leading to layering in the backscattering coefficient in the 1000–2000 m range. The EAE in this range is increased up to 0.5 (Fig. 8c), indicating that these layers may contain significant amounts of smoke.

The EAE of pure dust observed during SHADOW is slightly negative $A_{355/532}^{\alpha} \approx -0.1$ (Veselovskii et al., 2016). In Fig. 8a the EAE below 1500 m is about $0.2 \pm 0.2$, so the dust likely contains some amount of smoke. Values of EAE close to zero are observed in Fig. 8b and c below 1000, where the depolarization ratio increases up to $35 \pm 5\%$. Inside the dust layer $\beta_{355} < \beta_{532}$, so the corresponding backscattering Ångström exponent is negative. The negative values of $A_{355/532}^{\beta}$ have been already reported by Veselovskii et al. (2016), where negative BAE was attributed to an in-

crease of the imaginary part of the complex RI at 355 nm compared to 532 nm. In the center of the elevated layer at 3100 m $\delta_{532} = 14 \pm 3\%$, while at the top of this layer $\delta_{532}$ decreases to $6 \pm 1.5\%$ (Fig. 8a), indicating a possible presence of dust particles in the center of the elevated layer. The loading of elevated layer with dust particles is supported also by the profiles of $A_{355/532}^{\beta}$: for all three temporal intervals $A_{355/532}^{\beta}$ demonstrates the dip in the center of the elevated layer, while $A_{355/532}^{\alpha}$ and $A_{532/1064}^{\beta}$ do not decrease in the 2500–4000 m range. As mentioned, for pure dust $A_{355/532}^{\beta}$ is negative, so presence of dust in the center of smoke layer should decrease the backscattering Ångström exponent. The presence of dust in the smoke layer is not surprising, because
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Height-temporal distributions of the backscattering coefficient and particle depolarization ratio at 532 nm together with the water vapor mixing ratio derived from the Raman lidar measurements on the nights of 23–24 (a) and 24–25 December 2015 (b).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Time-height section of horizontal wind direction (arrows) and wind speed (color map) deduced from Doppler lidar during 24–25 December 2015. Leftward and downward arrows represent, respectively, easterly wind and northerly wind.</figcaption></figure>

<caption>**Figure 5.** Time-height section of the logarithmic range-corrected lidar signal (in arbitrary units) deduced from the Doppler lidar measurements during the night of 24–25 December 2015.</caption>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Vertical profiles of (a) temperature *T*, potential temperature Θ, (b) wind direction and speed, and (c) relative humidity (RH) and water vapor mixing ratio (WVMR) measured by the radiosonde in Dakar at 00:00 UTC on 25 December 2015. Solid line in plot (a) shows the aerosol backscattering coefficient at 532 nm in arbitrary units measured by the Raman lidar at 21:00 UTC on 24 December.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Vertical profiles of the aerosol backscattering (β<sub>355</sub>, β<sub>532</sub>, β<sub>1064</sub>) and extinction (α<sub>355</sub>, α<sub>532</sub>) coefficients together with the particle depolarization ratio (δ<sub>532</sub>) for three temporal intervals: 19:00–23:00, 01:00–04:00, and 04:00–07:00 UTC on 24–25 December 2015. The values of δ<sub>532</sub> are multiplied by factor 0.02.</figcaption></figure>

upwelling airflows in forest fires region can lift a signifi-
cant amount of dust together with biomass burning products
(Nisantzi et al., 2014). We should also mention that the spec-
tral dependence of the imaginary part (and thus A<sup>β</sup><sub>355/532</sub>) de-
pends on the dust origin. In particular, no negative values
of A<sup>β</sup><sub>355/532</sub> of dust were reported during the SAMUM cam-
paign, so the lidar ratios at 532 nm and 355 nm were close
(Tesche et al., 2011).

Lidar ratio profiles at 355 and 532 nm, for the same tempo-
ral intervals as in Fig. 7, are shown in Fig. 9. The lidar ratios
in the dust layer at 532 and 355 nm for 19:00–23:00 UTC

period are LR<sub>532</sub> = 55 ± 8 sr and LR<sub>355</sub> = 70 ± 10 sr, respec-
tively. At the top of the elevated layer, where the smoke par-
ticles are predominant, the lidar ratios for the same period
are higher: LR<sub>532</sub> = 65 ± 10 sr and LR<sub>355</sub> = 75 ± 11 sr. Due
to the presence of dust in the center of the elevated layer,
the height dependence of lidar ratios shows a decrease, with
a minimum at approximately 3000 m for all three temporal
intervals. The decrease is more pronounced at 532 nm be-
cause the difference between smoke and dust lidar ratios is
larger at this wavelength. The lidar ratios below 2000 m at
01:00–04:00 and 04:00–07:00 UTC become strongly oscil-
<figure><img src="image_11.png" /><figcaption>Figure 8. Extinction ($A_{355/532}^{\alpha}$) and backscattering ($A_{355/532}^{\beta}$, $A_{532/1064}^{\beta}$) Ångström exponents together with backscattering coefficient $\beta_{532}$ for the same three temporal intervals as in Fig. 7.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>Figure 9. Lidar ratios at 355 nm (solid lines) and 532 nm (dash lines) for three temporal intervals from Fig. 7.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Figure 10. Particle depolarization ratio as a function of the extinction Ångström exponent derived from data shown in Figs. 7 and 8.</figcaption></figure>

lating because of high gradients of backscattering and extinction coefficients at low altitudes and are not shown due to high uncertainties.

Figure 10 shows the dependence of the particle depolarization ratio $\delta_{532}$ on the extinction Ångström exponent derived from data in Figs. 7 and 8. The depolarization ratio monotonically decreases while EAE rises from 0 to 0.9. Thus observed high values of the depolarization ratio are attributed to big dust particles with EAE close to zero, while small smoke particles are characterized by low depolarization (below 10%). If depolarization ratios of smoke $\delta^s$ and dust $\delta^d$ are known, the contributions of smoke and dust particles to the total backscattering can be separated $\beta = \beta^s + \beta^d$ (Sugi-

moto and Lee, 2006; Tesche et al., 2009; Miffler et al., 2012; David et al., 2013; Burton et al., 2014). Assuming that the depolarization ratios of dust and smoke particles do not change with height the contributions $\beta^d$ and $\beta^s$ can be calculated as suggested by Tesche et al. (2009):

$$ \beta^d = \beta \frac{(\delta - \delta^s)(1 + \delta^d)}{(\delta^d - \delta^s)(1 + \delta)} \quad \text{and} \quad \beta^s = \beta - \beta^d. \qquad (6) $$

In our computations we used values $\delta^d = 35\%$ and $\delta^s = 7\%$.

The results of the decomposition of $\beta_{532}$ for $\beta_{532}^d$ and $\beta_{532}^s$ components for the same three temporal intervals as in Fig. 7 are shown in Fig. 11. This figure presents the total backscattering coefficient $\beta_{532}$ together with the particle
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Contributions of dust and smoke to the total backscattering coefficient $\beta_{532}$ together with particle depolarization ratio $\delta_{532}$ for three temporal intervals on 24–25 December 2015. Magenta and grey regions correspond to dust and smoke contribution to total scattering $\beta_{532} = \beta_{532}^d + \beta_{532}^s$.</figcaption></figure>

depolarization ratio $\delta_{532}$. The dust contribution to backscattering is marked with magenta, while the residual backscattering $\beta_{532} - \beta_{532}^d$ is attributed to the smoke and is marked with grey. For the height regions with low backscattering the uncertainty of $\beta_{532}$ is high, so the decomposition for these regions is not shown. The dust is predominant below 1700 m for 19:00–23:00 UTC period, but even the elevated layer contains a significant amount of dust: at 3100 m $\beta_{532}^d \approx 0.3\beta_{532}$. After 01:00 UTC the smoke layers descend (Fig. 3) and their contribution to backscattering becomes significant down to 1000 m height.

## 4 Comparison of lidar measurements with MERRA-2

MERRA-2 provides the vertical distribution of mass mixing ratios of five aerosol components, so for each of these components the extinction, backscattering coefficients, and depolarization ratios can be calculated. The vertical profiles of extinction coefficient of dust, BC, OC, SS, and SU, together with total extinction $\alpha_{532}$, are shown in Fig. 12 for 03:00 UTC and 21:00 UTC on 24 December 2015. At 03:00 UTC the aerosol is localized below 3000 m. Dust extinction is predominant, but contribution of OC to the total extinction coefficient rises with height, reaching maxima at 2250 m. The presence of a significant amount of OC agrees with the low values of depolarization ratio above 1500 m for this temporal interval in Fig. 3.

At 21:00 UTC an elevated layer with a maximum of extinction at 3150 m is observed (Fig. 12b). In this layer OC and dust provide similar contributions to extinction (about 40% at 3150 m height). From the results shown in Fig. 11a

we can estimate the contribution of dust to $\alpha_{532}$ in the center of the elevated layer as 30% (by assuming the dust lidar ratio $LR_{532} = 55$ sr), so the measured and simulated dust contributions are in good agreement. Below 1750 m the dust is the main contributor to the extinction coefficient providing 88% of $\alpha_{532}$ at 1000 m (Fig. 12b). The observed dust contribution to $\alpha_{532}$ at the same height is about 90% (Fig. 11a), which again shows good agreement between the model and measurements. Total contribution of BC and SU to extinction is below 20% in the elevated smoke layer, and in the near-surface dust layer their contribution is negligible. The extinction coefficients can be recalculated to the backscattering using model lidar ratios of the aerosol components. Figure 12c shows the profiles of backscattering coefficients at 532 nm computed for the same temporal interval as in Fig. 12b. The simulation of the backscattering coefficient is more challenging than that of extinction, because backscatter depends more strongly on the particle morphology and RI. A detailed comparison of measured and modeled profiles of backscattering coefficients will be performed later in this section.

As mentioned, the comparison of model and observed values is more straightforward for extinction coefficients. Figure 13 shows the time series of extinction profiles at 355 and 532 nm modeled for the night of 24–25 December 2015 at 18:00, 21:00, 00:00, 03:00, and 06:00 UTC. The profiles are shifted relative to each other by $0.2 \text{ km}^{-1}$. For comparison, the same figure presents the profiles of extinction coefficients derived from Raman lidar measurements. The model reproduces well the location of the elevated smoke layer as well as the top of the near-surface dust layer. However, the model
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Vertical profiles of extinction coefficients at (a) 03:00 UTC, (b) 21:00 UTC and (c) backscattering coefficients at 21:00 UTC on 24 December 2015 from MERRA-2 model at 532 nm. Profiles are given for five aerosol components: dust, black carbon (BC), organic carbon (OC), sea salt (SS), and sulfates (SU), together with total extinction $\alpha_{532}$ and backscattering $\beta_{532}$.</figcaption></figure>

does not resolve the oscillations of extinction profile below 2000 m at 03:00 and 06:00 UTC on 25 January.

To quantify the difference between the measured ($\alpha^{\text{meas}}$) and modeled ($\alpha^{\perp}$) extinction coefficients the difference $\Delta\alpha = \alpha^{\text{meas}} - \alpha^{\perp}$ was calculated. The statistical analysis of the frequency distribution of $\Delta\alpha$ for all five profiles in Fig. 13 shows that at 355 nm the mean value of $\Delta\alpha$ is -0.01 km$^{-1}$ and SD of 0.042 km$^{-1}$. With typical values of extinction coefficient in elevated smoke layer and near-surface dust layer being on the order of 0.2 km$^{-1}$, the relative difference of modeled and measured extinction is estimated to be below 25% for the time period considered. The results of statistical analysis for $\alpha_{532}$ nm are similar.

To analyze how well the model reproduces the temporal variations of aerosol optical depth, Fig. 14 presents AODs at 355 nm on 23–24 December 2015 for two height intervals: 750–2000 and 2500–4500 m. The first interval corresponds to the near-surface dust layer, while the second interval corresponds to the elevated smoke layer. The AOD is calculated from the Raman backscatter channel, and in the daytime measurements could be processed only in the dust layer due to enhanced background noise. Thus daytime measurements in the elevated smoke layer are not plotted. The time of the appearance of the smoke layer is well represented in the model results (about 00:00 UTC on 24 December), but the lidar-derived AOD of this layer increases rapidly from the first appearance of the layer, while in the model the rapid increase in AOD growth starts approximately 5 h later. The model predicts that the maximum value of AOD in the smoke layer (0.27) is reached at 20:00–24:00 UTC interval, which reasonably agrees with observations: mean value of measured AOD for this interval is 0.23 ± 0.02. After midnight the

modeled AOD of the smoke layer decreases quickly, while lidar measured AOD stays about 0.25. The measured AOD of the near-surface layer agrees with the model. The observed AOD exceeds the model values in the beginning (at 00:00 UTC on 24 December measured and modeled AODs are 0.24 and 0.175, respectively), but after 10:00 UTC the values are in better agreement. Thus, we can conclude that the model reproduces the temporal variability of AOD in the dust and smoke layers.

The agreement between modeled and observed extinction profiles provides an opportunity to test how well the backscattering coefficients can be modeled. Simulation of backscattering coefficients is especially challenging for dust for several reasons. First of all, we are not confident in the accuracy of the presumed scattering phase function in the backward direction. Second, the backscattering coefficient strongly depends on the particle RI, in particular on the imaginary part, which may vary over a wide range depending on dust origin. The in situ ground measurements in West Africa, performed during the SAMUM field campaign, demonstrate that the mean value of $m_1$ for dust episodes is about 0.003 at 532 nm and 0.02 at 355 nm. However, deviation from these mean values for every individual measurement can be significant (Müller et al., 2009; Kandler et al., 2011; Ansmann et al., 2011). The imaginary part of RI of dust in the model is assumed to be 0.007 at 355 nm, following previous OMI data analysis (Torres et al., 2007), and 0.0025 at 532 nm.

Figure 15 shows measured and modeled backscattering coefficients at 355 and 532 nm for the same five temporal intervals as in Fig. 13. At 355 nm the modeled and measured values agree for both the smoke and dust layers. However, at 532 nm the aerosol backscattering coefficients agree
<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Comparison of extinction profiles at (a) 355 nm and (b) 532 nm derived from Raman lidar measurements (line) and modeled by MERRA-2 (line + symbols) on the night of 24–25 December 2015. Model profiles are provided at 18:00, 21:00, 00:00, 03:00, and 06:00 UTC. The lidar measurements are given for temporal intervals centered at 19:00, 21:00, 00:00, 03:00, and 06:00 UTC. For each profile, 2 h of measurements are averaged. The profiles are shifted relatively to each other by 0.2 km⁻¹.</figcaption></figure>

only inside elevated layer, while below 1750 m the modeled $\beta_{532}$ significantly exceeds the measured values. As mentioned, the modeled lidar ratio $LR_{532}$ for the mixture is close to 40 sr at 1000 m, while the measured lidar ratio in the near-surface dust layer is 55 ± 8 sr. The reason for this disagreement could be that the assumed imaginary part of the RI for dust (0.0025 at 532 nm) is too low. Recall, however, that we cannot determine the imaginary part of the RI for dust by simply adjusting the modeled lidar ratio to the measured one, because the lidar ratio depends on several factors besides $m_1$, such as the PSD and the aspect ratio of the ellipsoids used in the model. It is possible that

<figure><img src="image_18.png" /><figcaption>**Figure 14.** Aerosol optical depth at 355 nm on 23–24 December 2015 obtained from MERRA-2 (line + symbols) and from the Raman lidar measurements (solid lines). The results are given for two height intervals: 750–2000 m (red) and 2500–4500 m (black). Zero of timescale corresponds to 00:00 UTC on 24 December.</figcaption></figure>

the PSD in the model is weighted too much toward fine-mode dust. The modeled and measured particle intensive parameters, such as extinction $A_{355-532}^{\alpha}$ and backscattering $A_{355-532}^{\beta}$ Ångström exponents and the particle depolarization ratio $\delta_{532}$, are shown in Fig. 16. The measurements are averaged over 19:00–23:00 UTC interval while modeled values are given for 21:00 UTC. The model reproduces well the observed vertical distribution of $A_{355-532}^{\alpha}$ in both the dust and the elevated layer. As follows from Fig. 7a, inside the dust layer $\beta_{355} < \beta_{532}$, so the corresponding $A_{355-532}^{\beta}$ is negative with a minimum value of about −0.4. The model predicts values of $A_{355-532}^{\beta}$ as low as −1.4. The modeled BAE is sensitive to the choice of the imaginary part of RI at 355 and 532 nm and, as mentioned, the chosen $m_1(532) = 0.0025$ may be too low for this episode. In the elevated layer the modeled $A_{355-532}^{\beta}$ is close to the observed one. The modeled BAE has no minimum in the center of the elevated layer, because the modeled ratio of dust and OC aerosol concentrations shows only a small variation throughout the elevated layer.

The model reproduces reasonably well the depolarization in the elevated layer, but inside the dust layer the modeled $\delta_{532}$ is significantly lower than what is observed (22 % compared to 35 %). This problem is well known: the spheroidal model underestimates the depolarization ratio when typical dust PSD and complex RI are used (Veselovskii et al., 2010; Wiegner et al., 2009; Müller et al., 2013; Nowottnick et al., 2015).

One of the MERRA-2 data products is WVMR, which helps to identify atmospheric parcels, is critically important for determining atmospheric stability, and serves as the source of water for aerosol hygroscopic growth. Figure 17 shows five model profiles of WVMR together with the results of Raman lidar measurements for the same temporal
<figure><img src="image_19.png" /></figure>

**Figure 15.** Backscattering coefficients at (a) 355 nm and (b) 532 nm measured by Raman lidar (solid line) and modeled by MERRA-2 (line + symbols) on the night of 24–24 December 2015. Profiles are shifted relatively to each other by 0.0025 km⁻¹ sr⁻¹. The temporal intervals are the same as in Fig. 13.

intervals as in Fig. 13. The model reproduces rather well the WVMR profile inside the elevated layer (2500–4500 m) on 24 December, though on 25 December the modeled values in this range are lower than the observations. In the near-surface dust layer, the deviation of modeled values from the measurements is larger. Statistical analysis of the deviation of modeled values from lidar measurements for all five profiles shows that mean difference is 0.04 g kg⁻¹ with SD of differences of 1.6 g kg⁻¹. Thus in the elevated layer, where WVMR is approximately 8 g kg⁻¹, the agreement is quite good, but in the dust layer, which is characterized by low water vapor content (below 4 g kg⁻¹), the difference may be up to 40 %.

<figure><img src="image_20.png" /></figure>

**Figure 16.** Extinction ($A_{355/532}^{\alpha}$) and backscattering ($A_{355/532}^{\beta}$) Ångström exponents together with the particle depolarization ratio $\delta_{532}$ obtained from lidar measurements (line) and from MERRA-2 modeling (line + symbols). Lidar data are averaged over the 19:00–23:00 UTC period while model data are given for 21:00 UTC.

<figure><img src="image_21.png" /></figure>

**Figure 17.** Water vapor mixing ratio derived from Raman lidar measurements (solid line) and obtained from the model (line + symbols) on the night of 24–25 December 2015. Temporal intervals are the same as in Fig. 13. The profiles are shifted relatively each other by 5 g kg⁻¹.

## 5 Inversion of lidar measurements to particle microphysical properties

In the previous section, as validation of the model output we compared the modeled aerosol optical parameters, such as extinction, backscattering coefficients, and depolarization ratio with the values derived from lidar measurements in a straightforward way. The comparison of particle micro-
physical properties such as volume, effective radius, and complex RI, however, is not straightforward, since it needs inversion of the measurements and requires additional assumptions. In the case of dust particles the inversion becomes especially challenging for the following reasons:

- The size distribution of dust contains a strong coarse mode with particle radii extending up to ~15 µm, and the estimation of properties for such big particles is difficult since measurements are only performed in the wavelength range 0.355–1.064 µm.

- The inversions have to consider the RI as spectrally independent. In fact, the imaginary part of the dust RI is spectrally dependent with a strong enhancement at 355 nm compared to 532 nm.

- The dust particles are not spherical and so the application of Mie formulas for the forward modeling results in errors in computing the scattering phase function.

Regarding the shape issue, one of the ways to mimic the scattering properties of dust particles is to use the model of randomly oriented spheroids (Mishchenko et al., 1997; Dubovik et al., 2006). The implementation of this model for inversion of dust lidar measurements is described in Veselovskii et al. (2010, 2016) and Müller et al. (2013). This algorithm was used also for inversion of our $3\beta + 2\alpha$ observations. The range of particle radius in the inversion has been set to a minimum and maximum of 0.075 and 15 µm, respectively. The real part of RI was allowed to vary in the range 1.35–1.65, while the imaginary part varied in the range 0–0.02. The RI was assumed to be spectrally independent. The effects of a possible spectral dependence of the imaginary part of RI were considered in Veselovskii et al. (2016).

Profiles of the effective radius, volume density, and real part of the RI retrieved from optical measurements in Fig. 7a are shown in Fig. 18. The inversion was performed for two cases, with the assumption of all spherical particles or all spheroids. A realistic solution (for the mixture of spherical and non-spherical particles) should be closer to spheroids in the dust layer, while in the elevated layer (where depolarization ratio is below 15 %) it should be closer to the results obtained with spheres. The model results provided by MERRA-2 are shown on the same plot. The effective radius and volume density obtained in assumption of spherical particles are always higher than the values obtained with spheroids. The modeled effective radius at 1100 m height is 1.1 µm, which is close to $r_{eff} = 0.95 \pm 0.3$ µm obtained from lidar measurements using the spheroids model. Comparing the lidar retrievals with model in the dust layer, we should keep in mind that inside 1500–2000 m height range the dust particles are mixed with biomass burning products; thus the use of only spheroids in retrieval underestimates the effective radius and volume. Moreover, accounting for the spectral dependence of the imaginary part of the dust may additionally increase the

retrieved values of V and $r_{eff}$ by factor 1.2–1.3 (Veselovskii et al., 2016).

Lidar-derived effective radius in the elevated layer at 3000 m is approximately 0.4 and 0.5 µm when spheroids and spheres are used, respectively, while the modeled value is 0.3 µm. The reason for the lower value of modeled effective radius is the contribution of black carbon, which is characterized by small size and relatively low hygroscopic growth. Recall that in the inversion of lidar measurements, the smallest radius considered is 0.075 µm. Modeled values of the volume density agree well with lidar retrievals in both dust and elevated layers.

The estimation of the real part of RI from lidar measurements is sensitive to the type of kernel functions chosen for retrieval. In the regularization algorithm the treatment of dust particles as spheres strongly underestimates $m_R$ (Veselovskii et al., 2010), so results obtained with spheres in the dust layer are not shown in Fig. 18c. At 1000 m the $m_R$ retrieved with spheroids is $1.52 \pm 0.05$, which agrees well with the modeled value. Inside the elevated smoke layer, where fine-mode particles predominate, the application of spheroids overestimates $m_R$. The lidar-derived real part of RI at 3000 m is $1.43 \pm 0.05$ for spheres and $1.51 \pm 0.05$ for spheroids, so we expect that the true value would lie within this. The simulated value of $m_R = 1.50$ in the elevated layer is quite high, which is again the result of BC contribution.

The single scattering albedo (SSA) is one of the key parameters to be retrieved and conclusions about the potential of the multiwavelength lidar method strongly rely on its ability to profile SSA. Figure 19 shows SSA at 355, 532, and 1064 nm. As mentioned, the spectral dependence of $m_I$ was not accounted for and the algorithm retrieves an average value of the imaginary part over the interval of 355–1064 nm. In particular, for dust and OC the imaginary part is underestimated at 355 nm and overestimated at 532 and 1064 nm. As a result, in the dust layer the retrieved SSA exceeds the model values at 355 nm, while at 532 nm and 1064 the situation is opposite. Still at a height of 1000 m, the difference between modeled and lidar-derived SSAs is below 0.04 for all wavelengths. In the elevated layer, where the spectral dependence of $m_I$ is less pronounced, the simulated and retrieved SSAs agree well with a corresponding difference of less than 0.02.

## 6 Summary and conclusion

The synergy of lidar observations with the aerosol transport model has a great potential to improve the characterization of aerosols properties, and as a first step in such synergy one has to demonstrate how well observations and models agree and describe the same aerosol scenario. For that we have considered a smoke–dust episode over West Africa to compare the vertical profiles of particle parameters modeled by MERRA-2 and retrieved from Raman lidar measurements. In the case selected, the simultaneous presence of the dust and smoke
<figure><img src="image_22.png" /><figcaption>**Figure 18.** Profiles of (a) effective radius, (b) particle volume, and (c) real part of the refractive index on 24 December 2015 retrieved from $3\beta+2\alpha$ lidar measurements shown in Fig. 7a (solid symbols) and provided by MERRA-2 for 21:00 UTC (open symbols). Inversion of lidar measurements was performed in assumption of spherical particles (s) and using the model of spheroids (ns).</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Figure 19.** The single scattering albedo at 355 nm (blue), 532 nm (green), and 1064 nm (red) on 24 December 2015 retrieved from $3\beta + 2\alpha$ lidar measurements shown in Fig. 7a (solid symbols) and provided by the MERRA-2 model for 21:00 UTC (line + open symbols). For inversion of lidar data the spheroids (ns) were used below 2000 m and spheres (s) above 2000 m.</figcaption></figure>

layers resulted in significant height variation of particle parameters, providing a good opportunity to test the models' capability to reproduce complicated vertical structure. Modeled and observed vertical profiles of $\alpha_{355}$ and $\alpha_{532}$ show good similarity: MERRA-2 provides the correct location of both the near-surface and elevated layers.

The modeling of the dust lidar ratio is challenging due to irregularity of the particles shape and due to the spectral dependence of the imaginary part of the RI. The $m_1$ can change significantly for dust of different origin and this variability may be accounted for in future model developments. The modeled at 355 nm the lidar ratio of 65 sr in the near-surface dust layer is close to the observed value ($70 \pm 10$ sr). At 532 nm, however, the simulated dust lidar ratio (about 40 sr) is lower than measurements ($55 \pm 7$ sr). This discrepancy may be an indication that $m_1$ of dust during the episode considered is higher than the value assumed in the model. Another possible explanation is that the model PSD is too much weighted toward fine-mode dust. The measured lidar ratios at the top of the elevated layer, where smoke particles are predominant, are $LR_{355} = 75 \pm 11$ sr and $LR_{532} = 70 \pm 10$ sr, which is close to the corresponding model values for organic carbon of 71 and 66 sr, respectively.

MERRA-2 predicts the existence of a significant amount of dust in the elevated smoke layer, and the high values of observed depolarization ratio agree with this prediction. The existence of minima of $A_{355/532}^\beta$ in the center of the elevated layer, characterized by the highest $\delta_{532}$, also supports this finding. Moreover, the lidar ratios at both 355 and 532 nm also have a minima in the center of the layer because the lidar ratio of dust is lower than that of smoke. The contributions of dust and smoke particles to the aerosol backscattering and extinction coefficient at 532 nm evaluated from particle depolarization ratio agree with the values provided by the model. Of course an analysis of only one episode is not sufficient for broad conclusions regarding how well the model reproduces the vertical distribution of particle properties. More measurements at different locations are needed. However, the results presented here demonstrate that observations and the
MERRA-2 model contribute in a complementary way, allow-
ing the separation of the contributions of different chemical
component of the aerosol mixture.

The motivation for this work is to show that the aerosol
transport model has sufficient skill to serve as an additional
constraint in inversion of 3β + 2α lidar observations and de-
velopment of such constrained inversion is in progress. As-
similation of lidar measured parameters in the model is the

subject of our future efforts.

*Data availability.* MPL observations for SHADOW campaign are available at http://loa-ptfi.univ-lille1.fr/lidar/calendars/cal_2015_mbour.php (LOA, 2018).
Appendix A: Optical properties of aerosol components
in MERRA-2 model

<figure><img src="image_24.png" /></figure>

**Figure A1.** Dependence of the growth factor of organic carbon, black carbon, sulfates, and sea salt on relative humidity (RH) used in MERRA-2. For the sea salt the results are given for five size bins from Table A1. The growth factor increases with increase of bin number.

Table A1 summarizes the main characteristics of five
aerosol components: dust, sea salt, black carbon, organic car-
bon, and sulfates used in MERRA-2 model. For dust and sea
salt five size bins are considered. All values are given for
the relative humidity of 0. Thus OC, BC, and SU with the
effective radii of 0.09, 0.04, and 0.157 µm, respectively, are
presented by the fine fraction only, while dust and sea salt
contribute to both fine and coarse fractions.

The dust particles are assumed to be hydrophobic, but
other aerosol components may present significant hygro-
scopic growth. To account for the effect of RH, the growth
factor *g*, which is the ratio of particle radius at current RH
to the dry particle radius, is introduced. Figure A1 shows de-
pendence of the growth factor of different aerosol compo-
nents on RH. For sea salt the results are given for five size
bins from Table 1. Each bin has a different growth factor: *g*
increases with increase of particle radius. Relative humidity
modifies also the particle complex refractive index. Depen-
dence of the real and the imaginary part of particle compo-
nents on RH is shown in Fig. A2. For dry sea salt particles RI
is supposed to be the same for all size bins. However, in the
process of hygroscopic growth the RI of different bins be-
haves differently: both *m*<sub>R</sub> and *m*<sub>I</sub> decrease with bin number
 radius) increasing.
<table><caption>**Table A1.** Parameters of the aerosol components, such as minimal radius ($r_{min}$), maximal radius ($r_{max}$), effective radius ($r_{eff}$), and real ($m_R$) and imaginary ($m_I$) part of the refractive index at 355, 532, and 1064 nm used in MERRA-2 model. For dust and sea salt, five size bins are considered. All values are given for RH = 0.</caption><thead><tr><th>Component</th><th></th><th>r<sub>min</sub><br>μm</th><th>r<sub>max</sub><br>μm</th><th>r<sub>eff</sub><br>μm</th><th>m<sub>R355</sub></th><th>m<sub>R532</sub></th><th>m<sub>R1064</sub></th><th>m<sub>I355</sub></th><th>m<sub>I532</sub></th><th>m<sub>I1064</sub></th></tr></thead><tbody><tr><td rowspan="5">Dust</td><td>Bin 1</td><td>0.1</td><td>1.0</td><td>0.64</td><td>1.53</td><td>1.53</td><td>1.53</td><td>0.007</td><td>0.0026</td><td>0.0022</td></tr><tr><td>Bin 2</td><td>1</td><td>1.5</td><td>1.32</td><td>1.53</td><td>1.53</td><td>1.53</td><td>0.007</td><td>0.0026</td><td>0.0022</td></tr><tr><td>Bin 3</td><td>1.5</td><td>3.0</td><td>2.30</td><td>1.53</td><td>1.53</td><td>1.53</td><td>0.007</td><td>0.0026</td><td>0.0022</td></tr><tr><td>Bin 4</td><td>3.0</td><td>7.0</td><td>4.17</td><td>1.53</td><td>1.53</td><td>1.53</td><td>0.007</td><td>0.0026</td><td>0.0022</td></tr><tr><td>Bin 5</td><td>7.0</td><td>10.0</td><td>7.67</td><td>1.53</td><td>1.53</td><td>1.53</td><td>0.007</td><td>0.0026</td><td>0.0022</td></tr><tr><td rowspan="5">Sea salt</td><td>Bin 1</td><td>0.03</td><td>0.1</td><td>0.08</td><td>1.51</td><td>1.50</td><td>1.47</td><td>2.9E-7</td><td>1.2E-8</td><td>1.97E-4</td></tr><tr><td>Bin 2</td><td>0.1</td><td>0.5</td><td>0.27</td><td>1.51</td><td>1.50</td><td>1.47</td><td>2.9E-7</td><td>1.2E-8</td><td>1.97E-4</td></tr><tr><td>Bin 3</td><td>0.5</td><td>1.5</td><td>1.07</td><td>1.51</td><td>1.50</td><td>1.47</td><td>2.9E-7</td><td>1.2E-8</td><td>1.97E-4</td></tr><tr><td>Bin 4</td><td>1.5</td><td>5</td><td>2.55</td><td>1.51</td><td>1.50</td><td>1.47</td><td>2.9E-7</td><td>1.2E-8</td><td>1.97E-4</td></tr><tr><td>Bin 5</td><td>5</td><td>10</td><td>7.3</td><td>1.51</td><td>1.50</td><td>1.47</td><td>2.9E-7</td><td>1.2E-8</td><td>1.97E-4</td></tr><tr><td>OC</td><td></td><td>0.01</td><td>0.29</td><td>0.09</td><td>1.53</td><td>1.53</td><td>1.52</td><td>0.048</td><td>0.009</td><td>0.016</td></tr><tr><td>BC</td><td></td><td>0.01</td><td>0.29</td><td>0.04</td><td>1.75</td><td>1.75</td><td>1.75</td><td>0.46</td><td>0.44</td><td>0.44</td></tr><tr><td>SU</td><td></td><td>0.01</td><td>0.29</td><td>0.157</td><td>1.45</td><td>1.43</td><td>1.42</td><td>1E-8</td><td>1E-8</td><td>2.9E-6</td></tr></tbody></table>

<figure><img src="image_25.png" /><figcaption>**Figure A2.** Dependence of the real and imaginary part of the refractive index of organic carbon (OC), black carbon (BC), sulfates (SU), and sea salt (SS) on relative humidity (RH) used in the MERRA-2 model. For the sea salt the results are given for five size bins from Table 1. Both $m_R$ and $m_I$ decrease with bin number increasing.</figcaption></figure>
*Competing interests.* The authors declare that they have no conflict of interest.

*Acknowledgements.* The authors are very grateful to IRD-Dakar (Institut de Recherche pour le Développement) for their welcome and efficient support and also thank Labex CaPPA for supporting this campaign. The CaPPA project (Chemical and Physical Properties of the Atmosphere) is funded by the French National Research Agency (ANR) through the PIA (Programme d'Investissement d'Avenir) under contract "ANR-11-LABX-0005-01" and by the Regional Council "Nord-Pas de Calais" and the "European Funds for Regional Economic Development" (FEDER). Development of lidar retrieval algorithms was supported by Russian Science Foundation (project 16-17-10241).

Edited by: Omar Torres

Reviewed by: three anonymous referees

## References

Alexandrov, M. and Mishchenko, M.: Information content of bistatic lidar observations of aerosols from space, Opt. Express, 25, 134–150, 2017.

Ansmann, A., Wandinger, U., Riebesell, M., Weitkamp, C., and Michaelis, W.: Independent measurement of extinction and backscatter profiles in cirrus clouds by using a combined Raman elastic-backscatter lidar, Appl. Optics, 31, 7113–7131, 1992.

Ansmann, A., Petzold, A., Kandler, K., Tegen, I., Wendisch, M., Müller, D., Weinzierl, B., Müller, T., and Heintzenberg, J.: Saharan mineral dust experiments SAMUM-1 and SAMUM-2: what have we learned?, Tellus B, 63, 403–429, 2011.

Böckmann, C., Miranova, I., Müller, D., Scheidenbach, L., and Nessler, R.: Microphysical aerosol parameters from multiwavelength lidar, J. Opt. Soc. Am. A, 22, 518–528, 2005.

Buchard, V., da Silva, A. M., Colarco, P. R., Darmenov, A., Randles, C. A., Govindaraju, R., Torres, O., Campbell, J., and Spurr, R.: Using the OMI aerosol index and absorption aerosol optical depth to evaluate the NASA MERRA Aerosol Reanalysis, Atmos. Chem. Phys., 15, 5743–5760, https://doi.org/10.5194/acp-15-5743-2015, 2015.

Buchard, V., da Silva, A. M., Randles, C. A., Colarco, P., Ferrare, R., Hair, J., Hostetler, C., Tackett, J., and Winker, D.: Evaluation of the surface PM2.5 in Version 1 of the NASA MERRA Aerosol Reanalysis over the United States, Atmos. Environ., 125, 100–111, 2016.

Buchard, V., Randles, C. A., da Silva, A. M., Darmenov, A., Colarco, P. R., Govindaraju, R., Ferrare, R., Hair, J., Beyersdorf, A. J., Ziemba, L. D., and Yu, H.: The MERRA-2 aerosol reanalysis, 1980 onward. Part II: Evaluation and case studies, J. Climate, 30, 6851–6872, 2017.

Burton, S. P., Vaughan, M. A., Ferrare, R. A., and Hostetler, C. A.: Separating mixtures of aerosol types in airborne High Spectral Resolution Lidar data, Atmos. Meas. Tech., 7, 419–436, https://doi.org/10.5194/amt-7-419-2014, 2014.

Burton, S. P., Hair, J. W., Kahnert, M., Ferrare, R. A., Hostetler, C. A., Cook, A. L., Harper, D. B., Berkoff, T. A., Seaman, S. T., Collins, J. E., Fenn, M. A., and Rogers, R. R.: Observations of the spectral dependence of linear particle depolarization ratio of aerosols using NASA Langley airborne High Spectral Resolution Lidar, Atmos. Chem. Phys., 15, 13453–13473, https://doi.org/10.5194/acp-15-13453-2015, 2015.

Burton, S. P., Chemyakin, E., Liu, X., Knobelspiesse, K., Stamnes, S., Sawamura, P., Moore, R. H., Hostetler, C. A., and Ferrare, R. A.: Information content and sensitivity of the $3\beta + 2\alpha$ lidar measurement system for aerosol microphysical retrievals, Atmos. Meas. Tech., 9, 5555–5574, https://doi.org/10.5194/amt-9-5555-2016, 2016.

Chemyakin, E., Müller, D., Burton, S., Kolgotin, A., Hostetler, C., and Ferrare, R.: Arrange and average algorithm for the retrieval of aerosols parameters from multiwavelength HSRL/Raman lidar data, Appl. Optics, 53, 7252–7266, 2014.

Chemyakin, E., Burton, S., Kolgotin, A., Müller, D., Hostetler, C., and Ferrare, R.: Retrieval of aerosol parameters from multiwavelength lidar: investigation of the underlying inverse mathematical problem, Appl. Optics, 55, 2188–2202, 2016.

Chin, M., Ginoux, P., Kinne, S., Torres, O., Holben, B., Duncan, B., Martin, R. V., Logan, J. A., Higurashi, A., and Nakajima, T.: Tropospheric aerosol optical thickness from the GOCART model and comparisons with satellite and Sun photometer measurements, J. Atmos. Sci., 59, 461–483, 2002.

Colarco, P., da Silva, A., Chin, M., and Diehl, T.: Online simulations of global aerosol distributions in the NASA GEOS-4 model and comparisons to satellite and ground-based aerosol optical depth, J. Geophys. Res., 115, D14207, https://doi.org/10.1029/2009JD012820, 2010.

Colarco, P. R., Nowotnick, E. P., Randles, C. A., Yi, B., Yang, P., Kim, K.-M., Smith, J., and Bardeen, C. G.: Impact of radiatively interactive dust aerosols in the NASA GEOS-5 climate model: sensitivity to dust particle shape and refractive index, J. Geophys. Res.-Atmos., 119, 753–786, https://doi.org/10.1002/2013JD020046, 2014.

Colarco, P. R., Gassó, S., Ahn, C., Buchard, V., da Silva, A. M., and Torres, O.: Simulation of the Ozone Monitoring Instrument aerosol index using the NASA Goddard Earth Observing System aerosol reanalysis products, Atmos. Meas. Tech., 10, 4121–4134, https://doi.org/10.5194/amt-10-4121-2017, 2017.

David, G., Thomas, B., Nousiainen, T., Miffre, A., and Rairoux, P.: Retrieving simulated volcanic, desert dust and sea-salt particle properties from two/three-component particle mixtures using UV-VIS polarization lidar and T matrix, Atmos. Chem. Phys., 13, 6757–6776, https://doi.org/10.5194/acp-13-6757-2013, 2013.

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

Freudenthaler, V., Esselborn, M., Wiegner, M., Heese, B., Tesche, M., Ansmann, A., Müller, D., Althausen, D., Wirth, M., Fix, A., Ehret, G., Knippertz, P., Toledano, C., Gasteiger, J., Garhammer, M., and Seefeldner, M.: Depolarization ratio profiling at several wavelengths in pure Saharan dust during SAMUM 2006, Tellus B, 61, 165–179, 2009.
Gelaro, R., McCarty, W., Suarez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., Da Silva, A. M., Gu, W., Kim, G. K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454 2017.

Gerber, H. E.: Relative-humidity parameterization of the Navy Aerosol Model (NAM), Tech. Rep. NTIS ADA1632090, Naval Research Laboratory, Washington, DC, 1985.

Hammer, M. S., Martin, R. V., van Donkelaar, A., Buchard, V., Torres, O., Ridley, D. A., and Spurr, R. J. D.: Interpreting the ultraviolet aerosol index observed with the OMI satellite instrument to understand absorption by organic aerosols: implications for atmospheric oxidation and direct radiative effects, Atmos. Chem. Phys., 16, 2507–2523, https://doi.org/10.5194/acp-16-2507-2016, 2016.

Hess, M., Koepke, P., and Schult, I.: Optical properties of aerosols and clouds: the software package OPAC, B. Am. Meteorol. Soc., 79, 831–844, 1998.

Kahnert, M. and Andersson, E.: How much information do extinction and backscattering measurements contain about the chemical composition of atmospheric aerosol?, Atmos. Chem. Phys., 17, 3423–3444, https://doi.org/10.5194/acp-17-3423-2017, 2017.

Kandler, K., Lieke, K., Benker, N., Emmel, C., Küpper, M., Müller-Ebert, D., Ebert, M., Scheuvens, D., Schladitz, A., Schütz, L., and Weinbruch, S.: Electron microscopy of particles collected at Praia, Cape Verde, during the Saharan Mineral Dust Experiment: particle chemistry, shape, mixing state and complex refractive index, Tellus B, 63, 475–496, 2011.

Klett, J. D.: Stable analytical inversion solution for processing lidar returns, Appl. Optics, 20, 211–220, 1981.

Kolgotin, A., Müller, D., Chemyakin, E., and Romanov, A.: Improved identification of the solution space of aerosol microphysical properties derived from the inversion of profiles of lidar optical data, part 1: theory, Appl. Optics, 55, 9839–9849, 2016.

LOA: MPL observations for SHADOW campaign, available at: http://loa-ptfi.univ-lille1.fr/lidar/calendars/cal_2015_mbour.php, last access: 15 February 2018.

Mifflre, A., David, G., Thomas, B., Rairoux, P., M. Fjaeraa, A., Kristiansen, N. I., and Stohl, A.: Volcanic aerosol optical properties and phase partitioning behavior after long-range advection characterized by UV-Lidar measurements, Atmos. Environ., 48, 76–84, 2012.

Mishchenko, M. I., Travis, L. D., Kahn, R. A., and West, R. A.: Modeling phase functions for dustlike tropospheric aerosols using a mixture of randomly oriented polydisperse spheroids, J. Geophys. Res., 102, 16831–16847, 1997.

Müller, D., Wandinger, U., and Ansmann, A.: Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: theory, Appl. Optics, 38, 2346–2357, 1999.

Müller, D., Mattis, I., Wandinger, U., Ansmann, A., Althausen, D., and Stohl, A.: Raman lidar observations of aged Siberian and Canadian forest fire smoke in the free troposphere over Germany in 2003: microphysical

particle characterization, J. Geophys. Res., 110, D17201, https://doi.org/10.1029/2004JD005756, 2005.

Müller, D., Veselovskii, I., Kolgotin, A., Tesche, M., Ansmann, A., and Dubovik, O.: Vertical profiles of pure dust (SAMUM-1) and mixed smoke-dust plumes (SAMUM-2) inferred from inversion of multiwavelength Raman/polarization lidar data and comparison to AERONET retrievals and in-situ observations, Appl. Optics, 52, 3178–3202, 2013.

Müller, D., Böckmann, C., Kolgotin, A., Schneidenbach, L., Chemyakin, E., Rosemann, J., Znak, P., and Romanov, A.: Microphysical particle properties derived from inversion algorithms developed in the framework of EARLINET, Atmos. Meas. Tech., 9, 5007–5035, https://doi.org/10.5194/amt-9-5007-2016, 2016.

Müller, T., Schladitz, A., Massling, A., Kaaden, N., Kandler, K., and Wiedensohler, A.: Spectral absorption coefficients and imaginary parts of refractive indices of Saharan dust during SAMUM-1, Tellus B, 61, 79–95, 2009.

Nisantzi, A., Mamouri, R. E., Ansmann, A., and Hadjimitsis, D.: Injection of mineral dust into the free troposphere during fire events observed with polarization lidar at Limassol, Cyprus, Atmos. Chem. Phys., 14, 12155–12165, https://doi.org/10.5194/acp-14-12155-2014, 2014.

Nowotnick, E. P., Colarco, P. R., Welton, E. J., and da Silva, A.: Use of the CALIOP vertical feature mask for evaluating global aerosol models, Atmos. Meas. Tech., 8, 3647–3669, https://doi.org/10.5194/amt-8-3647-2015, 2015.

Pérez-Ramírez, D., Whiteman, D. N., Veselovskii, I., Kolgotin, A., Korenskiy, M., and Alados-Arboledas, L.: Effects of systematic and random errors on the retrieval of particle microphysical properties from multiwavelength lidar measurements using inversion with regularization, Atmos. Meas. Tech., 6, 3039–3054, https://doi.org/10.5194/amt-6-3039-2013, 2013.

Randles, C. A., Da Silva, A. M., Buchard, V., Colarco, P. R., Darmenov, A., Govindaraju, R., Smirnov, A., Holben, B., Ferrare, R., Hair, J., Shinozuka, Y., and Flynn, J.: The MERRA-2 aerosol reanalysis, 1980 onward. Part I: System description and data assimilation evaluation, J. Climate, 30, 6823–6850, 2017.

Sugimoto, N. and Lee, C. H.: Characteristics of dust aerosols inferred from lidar depolarization measurements at two wavelength, Appl. Optics, 45, 7468–7474, 2006.

Tesche, M., Ansmann, A., Müller, D., Althausen, D., Engelmann, R., Freudenthaler, V., and Groß, S.: Vertically resolved separation of dust and smoke over Cape Verde using multiwavelength Raman and polarization lidars during Saharan Mineral Dust Experiment 2008, J. Geophys. Res., 114, D13202, https://doi.org/10.1029/2009JD011862, 2009.

Tesche, M., Groß, S., Ansmann, A., Müller, D., Althausen, D., Freudenthaler, V., and Esselborn, M.: Profiling of Saharan dust and biomass-burning smoke with multiwavelength polarization Raman lidar at Cape Verde, Tellus B, 63, 649–676, https://doi.org/10.1111/j.1600-0889.2011.00548.x, 2011.

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, https://doi.org/10.1029/2007JD008809, 2007.

Veselovskii, I., Kolgotin, A., Griaznov, V., Müller, D., Wandinger, U., and Whiteman, D.: Inversion with regularization for the retrieval of tropospheric aerosol parameters from
multi-wavelength lidar sounding, Appl. Optics, 41, 3685–3699, 2002.

Veselovskii, I., Kolgotin, A., Griaznov, V., Müller, D., Franke, K., and Whiteman, D. N.: Inversion of multiwavelength Raman lidar data for retrieval of bimodal aerosol size distribution, Appl. Optics, 43, 1180–1195, 2004.

Veselovskii, I., Kolgotin, A., Müller, D., and Whiteman, D. N.: Information content of multiwavelength lidar data with respect to microphysical particle properties derived from eigenvalue analysis, Appl. Optics, 44, 5292–5303, 2005.

Veselovskii, I., O. Dubovik, A. Kolgotin, T. Lapyonok, P. Di Girolamo, D. Summa, D. N. Whiteman, M. Mishchenko, and Tanré, D.: Application of randomly oriented spheroids for retrieval of dust particle parameters from multiwavelength lidar measurements, J. Geophys. Res., 115, D21203, https://doi.org/10.1029/2010JD014139, 2010.

Veselovskii, I., Goloub, P., Podvin, T., Bovchaliuk, V., Derimian, Y., Augustin, P., Fourmentin, M., Tanre, D., Korenskiy, M., Whiteman, D. N., Diallo, A., Ndiaye, T., Kolgotin, A., and Dubovik, O.: Retrieval of optical and physical properties of African dust from multiwavelength Raman lidar measurements during the SHADOW campaign in Senegal, Atmos. Chem. Phys., 16, 7013–7028, https://doi.org/10.5194/acp-16-7013-2016, 2016.

Whiteman, D., Melfi, S., and Ferrare, R.: Raman lidar system for measurement of water vapor and aerosols in the Earth’s atmosphere, Appl. Optics, 31, 3068–3082, 1992.

Wiegner, M., Gasteiger, J., Kandler, K., Weinzierl, B., Rasp, K., Esselborn, M., Freudenthaler, V., Heese, B., Toledano, C., Tesche, M., and Althausen, D.: Numerical simulations of optical properties of Saharan dust aerosols with emphasis on lidar applications, Tellus B, 61, 180–194, 2009.

Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y., and Powell, K. A.: Overview of the CALIPSO mission and CALIOP data processing algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323, 2009.