# Fluorescence lidar observations of wildfire smoke inside cirrus: a contribution to smoke-cirrus interaction research

Igor Veselovskii, Qiaoyun Hu, Albert Ansmann, Philippe Goloub, Thierry Podvin, Mikhail Korenskiy

► To cite this version:

Igor Veselovskii, Qiaoyun Hu, Albert Ansmann, Philippe Goloub, Thierry Podvin, et al.. Fluorescence lidar observations of wildfire smoke inside cirrus: a contribution to smoke-cirrus interaction research. Atmospheric Chemistry and Physics, 2022, 22, pp.5209-5221. 10.5194/acp-22-5209-2022. insu-03686318

HAL Id: insu-03686318

https://insu.hal.science/insu-03686318

Submitted on 3 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

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

# Fluorescence lidar observations of wildfire smoke inside cirrus: a contribution to smoke–cirrus interaction research

Igor Veselovskii¹, Qiaoyun Hu², Albert Ansmann³, Philippe Goloub², Thierry Podvin², and
Mikhail Korenskiy¹

¹Prokhorov General Physics Institute of the Russian Academy of Sciences, Moscow, Russia

²Univ. Lille, CNRS, UMR 8518 – LOA – Laboratoire d'Optique Atmosphérique, Lille, 59000, France

³Leibniz Institute for Tropospheric Research, Leipzig, Germany

**Correspondence:** Qiaoyun Hu (qiaoyun.hu@univ-lille.fr)

Received: 6 December 2021 – Discussion started: 8 December 2021

Revised: 1 February 2022 – Accepted: 23 February 2022 – Published: 21 April 2022

**Abstract.** A remote sensing method, based on fluorescence lidar measurements, that allows us to detect and to quantify the smoke content in the upper troposphere and lower stratosphere (UTLS) is presented. The unique point of this approach is that smoke and cirrus properties are observed in the same air volume simultaneously. In this article, we provide results of fluorescence and multiwavelength Mie–Raman lidar measurements performed at ATOLL (ATmospheric Observation at liLLe) observatory from Laboratoire d'Optique Atmosphérique, University of Lille, during strong smoke episodes in the summer and autumn seasons of 2020. The aerosol fluorescence was induced by 355 nm laser radiation, and the fluorescence backscattering was measured in a single spectral channel, centered at 466 nm and having 44 nm width. To estimate smoke particle properties, such as number, surface area and volume concentration, the conversion factors, which link the fluorescence backscattering and the smoke microphysical properties, are derived from the synergy of multiwavelength Mie–Raman and fluorescence lidar observations. Based on two case studies, we demonstrate that the fluorescence lidar technique provides the possibility to estimate the smoke surface area concentration within freshly formed cirrus layers. This value was used in the smoke ice nucleating particle (INP) parameterization scheme to predict ice crystal number concentrations in cirrus generation cells.

## 1 Introduction

Aerosol particles in the height regime of the upper-troposphere and lower-stratosphere (UTLS) play an important role in processes of heterogeneous ice formation; however, our current understanding of these processes is still insufficient for a trustworthy implementation in numerical weather and climate prediction models. The ability of aerosol particles to act as ice nucleating particles (INPs) depends on meteorological factors such as temperature and ice supersaturation (as a function of vertical velocity), as well as on the aerosol type in the layer in which cirrus developed (Kanji et al., 2017). Heterogeneous ice nucleation initiated by insoluble inorganic materials such as mineral dust has been studied

for a long time (e.g., DeMott et al., 2010, 2015; Hoose and Möhler, 2012; Murray et al., 2012; Boose et al., 2016; Schrod et al., 2017; Ansmann et al., 2019b), while the potential of omnipresent organic particles, especially of frequently occurring, aged and long-range-transported wildfire smoke particles, to act as INPs is less well explored and thus not well understood (Knopf et al., 2018). Wildfire smoke can reach the lower stratosphere via pyrocumulonimbus (pyroCb) convection (Fromm et al., 2010; Peterson et al., 2018, 2021; Ansmann et al., 2018; Hu et al., 2019; Khaykin et al., 2020) or via self-lifting processes (Boers et al., 2010; Ohneiser et al., 2021). It is widely assumed that the ability of smoke particles to serve as INPs mainly depends on the organic material (OM) in the shell of the coated smoke particles (Knopf
et al., 2018) but may also depend on mineral components in the smoke particles (Jahl et al., 2021). The ice nucleation efficiency may increase with increasing duration of the long-range transport as Jahl et al. (2021) suggested. Disregarding the progress made in this atmospheric research field during the last years, the link between ice nucleation efficiency and the smoke particle chemical and morphological properties is still largely unresolved (Umo et al., 2015; Grawe et al., 2016; China et al., 2017; Knopf et al., 2018).

To contribute to the field of smoke–cirrus interaction research, we present a laser remote sensing method that allows us simultaneously to detect and quantify the smoke particle amount inside of cirrus layers, together with cirrus properties, and to provide INP estimates in regions close to the cloud top where ice formation usually begins. The unique point of our approach is that, for the first time, smoke and cirrus properties are observed in the same air volume simultaneously. Recently, a first attempt (closure study) was performed to investigate the smoke impact on High Arctic cirrus formation (Engelmann et al., 2021). However, the aerosol measurements had to be performed outside the cloud layers, and then an assumption was needed in which the estimated aerosol (and estimated INP) concentration levels also hold inside the cirrus layers. Now, we propose a method to directly determine INP-relevant smoke parameters inside the cirrus layer during ice nucleation events. This also offers the opportunity to illuminate whether an INP reservoir can be depleted in cirrus evolution processes or not. Furthermore, this new lidar detection method permits a clear discrimination between, for example, smoke and mineral dust INPs.

Multiwavelength Mie–Raman lidars or high-spectral-resolution lidars (HSRLs) are favorable instruments to provide the vertical profiles of the physical properties of tropospheric aerosol particles. In particular, the inversion of the so-called 3β + 2α lidar observations, based on the measurement of height profiles of three aerosol backscatter coefficients at 355, 532 and 1064 nm and two extinction coefficients at 355 and 532 nm, allows us to estimate smoke microphysical properties (Müller et al., 1999, 2005; Veselovskii et al., 2002, 2015). However, the aerosol content in the UTLS height range can be low so that particle extinction coefficients cannot be determined with sufficient accuracy and are thus not available in the lidar inversion data analysis. To resolve this issue Ansmann et al. (2019a, 2021) used the synergy of polarization lidar measurements and Aerosol Robotic Network (AERONET) sunphotometer observations (Holben et al., 1998) to derive conversion factors (to convert backscatter coefficients into microphysical particle properties) and to estimate INP concentrations for dust and smoke aerosols with the retrieved aerosol surface area concentration as aerosol input.

Dust particles are very efficient ice nuclei in contrast to wildfire smoke particles. In this context, the following question arises: how can we unambiguously discriminate smoke from dust particles? This is realized by integrating a fluo-

rescence channel into a multiwavelength aerosol lidar (Richard et al., 2017; Richardson et al., 2019; Veselovskii et al., 2020, 2021). The fluorescence capacity of smoke (ratio of fluorescence backscattering to the overall aerosol backscattering) significantly exceeds corresponding values for other types of aerosol, such as dust or anthropogenic particles (Veselovskii et al., 2020, 2021), and thus allows us to discriminate smoke from other aerosol types. The fluorescence technique provides therefore the unique opportunity to monitor ice formation in well-identified wildfire smoke layers and thus to create a good basis for long-term investigations of smoke–cirrus interaction.

In this article, we present results of fluorescence and multiwavelength Mie–Raman lidar measurements performed at ATOLL (ATmospheric Observation at LiLLe) of the Laboratoire d'Optique Atmosphérique, University of Lille, during strong smoke episodes in the summer and autumn seasons of 2020. The results demonstrate that the fluorescence lidar is capable of monitoring the smoke in the UTLS height range and inside the cirrus clouds formed at or below the tropopause. We start with a brief description of the experimental setup in Sect. 2. In the first part of the result section (Sect. 3.1 and 3.2), it is explained how smoke optical properties can be quantified by using fluorescence backscattering information and how we can estimate smoke microphysical properties (volume, surface area and number concentration) from measured fluorescence backscatter coefficients. In this approach, multiwavelength Mie–Raman aerosol lidar observations are used in addition. The retrieved values of the smoke particle surface area concentration are then the aerosol input in the smoke INP estimation. A case study is discussed in Sect. 3.2. Two case studies are then presented in Sect. 3.3 to demonstrate the capability of a fluorescence lidar to monitor ice formation in extended smoke layers and to provide detailed information on aerosol microphysical properties and smoke-relate INP concentration levels.

## 2 Experimental setup

The multiwavelength Mie–Raman lidar LILAS (Lille Lidar AtmosphereS) is based on a tripled Nd:YAG laser with a 20 Hz repetition rate and pulse energy of 70 mJ at 355 nm. Backscattered light is collected by a 40 cm aperture Newtonian telescope, and the lidar signals are digitized with Lical transient recorders of 7.5 m range resolution, allowing simultaneous detection in the analog and photon counting mode. The system is designed for simultaneous detection of elastic and Raman backscattering, allowing the so called 3β + 2α + 3δ data configuration, including three particle backscattering coefficients (β355, β532, β1064), two extinction coefficients (α355, α532) and three particle depolarization ratios (δ355, δ532, δ1064). The particle depolarization ratio, determined as a ratio of cross- and co-polarized components of the particle backscattering coefficient, was calculated and
calibrated in the same way as described in Freudenthaler et al. (2009). The aerosol extinction and backscattering coefficients at 355 and 532 nm were calculated from Mie–Raman observations (Ansmann et al., 1992), while $\beta_{1064}$ was derived by the Klett method (Fernald, 1984; Klett, 1985). Additional information about atmospheric parameters was available from radiosonde measurements performed at Herstmonceux (UK) and Beauvechain (Belgium) stations, located 160 and 80 km away from the observation site respectively.

This lidar system is also capable of performing aerosol fluorescence measurements. A part of the fluorescence spectrum is selected by a wideband interference filter of 44 nm width centered at 466 nm (Veselovskii et al., 2020, 2021). The strong sunlight background at daytime restricts the fluorescence observations to nighttime hours. To characterize the fluorescence properties of aerosol, the fluorescence backscattering coefficient $\beta_F$ is calculated from the ratio of fluorescence and nitrogen Raman backscatters, as described in Veselovskii et al. (2020). This approach allows us to evaluate the absolute values of $\beta_F$ if the relative sensitivity of the channels is calibrated and the nitrogen Raman scattering differential cross-section $\sigma_R$ is known. In our research we used $\sigma_R = 2.744 \times 10^{-30} \text{ cm}^2 \text{ sr}^{-1}$ at 355 nm from Venable et al. (2011). All $\beta_F$ profiles presented in this work were smoothed with the Savitzky–Golay method using a second-order polynomial with 21 points in the window. The efficiency of fluorescence backscattering with respect to elastic backscattering $\beta_{532}$ is characterized by the fluorescence capacity $G_F = \frac{\beta_F}{\beta_{532}}$.

For most of atmospheric particles $\beta_F$ is proportional to the volume of dry matter, while dependence of $\beta_{532}$ on particle size is more complicated. As a result, $G_F$ depends not only on aerosol type but also on particle size and the relative humidity RH. Uncertainty of $\beta_F$ calculation depends on the chosen value of $\sigma_R$ and on the relative transmission of optical elements in fluorescence and nitrogen channels. These system parameters do not change with time. The relative sensitivity of photomultiplier tubes, however, may change. Regular calibration of the channels' relative sensitivity (Veselovskii et al., 2020) demonstrates that corresponding uncertainty can be up to 10%. At high altitudes the statistical uncertainty becomes predominant. We recall also that only a part of the fluorescence spectra was selected by the interference filter in the receiver, so provided values of $\beta_F$ and $G_F$ are specific to the filter used. Analyzing the fluorescence measurements we should keep in mind that the sensitivity of this technique can be limited by the fluorescence of optics in the lidar receiver. The minimal value of $G_F$, which we measured during observations in cloudy conditions in the lower troposphere, was about $2 \times 10^{-8}$. Thus, at least in the measurements with $G_F$ above this value, the contribution of optics fluorescence can be ignored.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Range-corrected lidar signal at 1064 nm, volume depolarization ratio at 1064 nm and fluorescence backscattering coefficient (in $10^{-4} \text{ Mm}^{-1} \text{ sr}^{-1}$) on 4–5 November 2020.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Vertical profiles of aerosol backscattering coefficient $\beta_{532}$ and fluorescence backscattering coefficient $\beta_F$ on 5 November 2020 for the periods (a) 02:00–05:30 UTC and (b) 18:00–22:15 UTC. Open symbols show the temperature profile measured by the radiosonde launched at Herstmonceux (UK).</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Smoke fluorescence in the presence of clouds on 24–25 November 2020. (a, b) Spatiotemporal variations in the range-corrected lidar signal and volume depolarization at 1064 nm. (c) Vertical profiles of the aerosol $\beta_{532}$ and fluorescence $\beta_F$ backscattering coefficients. (d) Aerosol backscattering coefficient $\beta_{532}$, together with smoke backscattering coefficient $\beta_{532}^s$, computed from $\beta_F$ for $G_F = 4.5 \times 10^{-4}$. Open symbols show temperature profile measured by the radiosonde at Herstmonceux.</figcaption></figure>

# 3 Results of the measurements

## 3.1 Observation of smoke particles in UTLS

Smoke particles produced by intensive fires and transported across the Atlantic are regularly observed in the UTLS height range over Europe (Müller et al., 2005; Hu et al., 2019; Baars et al., 2019, 2021). One such event, observed over Lille in the night of 4–5 November 2020, is shown in Fig. 1. The figure provides height-time displays of the range-corrected lidar signal and the volume depolarization ratio at 1064 nm, together with the fluorescence backscattering coefficient. A narrow smoke layer occurred in the upper troposphere in the period from 23:00–06:00 UTC. The smoke was detected at heights above 12 km after midnight. The particles caused a low volume depolarization ratio ($< 5 \%$) at 1064 nm and strong fluorescence backscattering ($\beta_F > 1.2 \times 10^{-4} \text{ Mm}^{-1}\text{sr}^{-1}$). The backward trajectory analysis indicated that the aerosol layer was transported over the Atlantic and contained products of North American wild fires.

Vertical profiles of aerosol $\beta_{532}$ and fluorescence $\beta_F$ backscattering coefficients for the period from 02:00–05:30 UTC are shown in Fig. 2a. The fluorescence capacity $G_F$ in the center of smoke layer (not shown) was about $4.5 \times 10^{-4}$. The depolarization ratio of aged smoke in the UTLS height range usually shows a strong spectral dependence (Haarig et al., 2018; Hu et al., 2019). For the case presented in Fig. 2a the particle depolarization ratio in the center of the smoke layer decreased from $16 \pm 4 \%$ at 355 nm ($\delta_{355}$) to $4 \pm 1 \%$ at 1064 ($\delta_{1064}$). The tropopause height $H_{tr}$ was at about 13 000 m; thus the main part of the smoke layer was below the tropopause. By the end of day the smoke layer became weaker ($\beta_F < 0.3 \times 10^{-4} \text{ Mm}^{-1}\text{sr}^{-1}$) and ascended up

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Range-corrected lidar signal at 1064 nm on 23–24 June 2020, revealing a smoke layer between 4500 and 5200 m height.</figcaption></figure>

to 14 500 m, which is above the tropopause. The corresponding vertical profiles of $\beta_{532}$ and $\beta_F$ are shown in Fig. 2b. The fluorescence capacity in the center of the layer is about $4.5 \times 10^{-4}$, which is close to the value observed during the 02:00–05:30 UTC period.

An important advantage of the fluorescence lidar technique is the ability to monitor smoke particles inside cirrus clouds. The results of smoke observations in the presence of ice clouds are shown in Fig. 3. Cirrus clouds occurred during the whole night in the height range from 6.0–10.0 km. To quantify the fluorescence backscattering inside the cloud (which was rather weak in this case), the lidar signals were averaged over the full 18:00–06:00 UTC time interval in Fig. 3a. The fluorescence backscatter coefficient shown in Fig. 3c decreased from $\beta_F = 0.015 \times 10^{-4} \text{ Mm}^{-1}\text{sr}^{-1}$ at 5000 m (near the cloud base) to a minimum value of $0.01 \times 10^{-4} \text{ Mm}^{-1}\text{sr}^{-1}$ at 7000 m inside the cirrus layer. Above the tropopause the fluorescence backscattering increased strongly and reached the maximum (about $0.3 \times 10^{-4} \text{ Mm}^{-1}\text{sr}^{-1}$) at 11 000–13 000 m height.
The analysis of fluorescence measurements performed during strong smoke episodes in the summer and autumn of 2020, when smoke layers from North American fires frequently reached Europe, demonstrates that the fluorescence capacity varied within the range of $2.5 \times 10^{-4}$ to $4.5 \times 10^{-4}$. The variations are a function of smoke composition, relative humidity and particle size. However, in the upper troposphere where relative humidity is low, $G_F$ was normally close to $4.5 \times 10^{-4}$. This relatively low range of $G_F$ variations allows the estimation of the backscattering coefficient attributed to the smoke particles from fluorescence measurements as follows:

$$ \beta_{532}^s = \frac{\beta_F}{G_F}. \qquad (1) $$

Figure 3d shows the smoke backscattering coefficient $\beta_{532}^s$, calculated from $\beta_F$ for $G_F = 4.5 \times 10^{-4}$, together with $\beta_{532}$. The dynamical range of $\beta_{532}$ variations is high. To make smoke backscattering visible above $H_{lr}$, $\beta_{532}$ is plotted in expanded scale in Fig. 3d. The $\beta_{532}^s$ values, though strongly oscillating above the tropopause, match the $\beta_{532}$, indicating that the smoke contribution to backscattering was predominant.

## 3.2 Estimation of smoke particle content based on fluorescence measurements

The possibility to detect fluorescence backscattering inside the cirrus clouds reveals also the opportunity for a quantitative characterization of the smoke content. This can be realized by a synergistic use of fluorescence and multiwavelength Mie-Raman lidar observations. The flow chart, summarizing the main steps of this procedure, is presented in Appendix A. For the smoke layers with sufficient optical depth, the number $N$, surface area $S$ and volume $V$ concentrations can be evaluated by inverting the $3\beta + 2\alpha$ observations consisting of three backscatter coefficients (355, 532, 1064 nm) and two extinction coefficients (355, 532 nm) (Müller et al., 1999; Veselovskii et al., 2002; Pérez-Ramírez et al., 2013). The conversion factors $C_N$, $C_S$ and $C_V$, introduced as

$$ C_N = \frac{N}{\beta_F}, \quad C_S = \frac{S}{\beta_F}, \quad C_V = \frac{V}{\beta_F}, \qquad (2) $$

allow the estimation of smoke particle concentration inside the clouds from fluorescence backscattering, assuming that smoke contribution to the fluorescence is predominant. Moreover, it allows the estimation of the particle concentration in weak smoke layers in UTLS, where $3\beta + 2\alpha$ observations are normally not available.

On 23–24 June 2020, a strong smoke layer was observed at 4500–5500 m height during the whole night (Fig. 4). The vertical profiles of the aerosol backscattering and extinction coefficients ($3\beta + 2\alpha$) are shown in Fig. 5a, while the particle depolarization ratios $\delta_{355}$, $\delta_{532}$ and $\delta_{1064}$ and the lidar ratios at 355 and 532 nm (LR$_{355}$, LR$_{532}$) are presented in Fig. 5b.

<figure><img src="image_8.png" /><figcaption>Figure 5. Smoke layer on 23–24 June 2020. (a) Vertical profiles of backscattering ($\beta_{355}$, $\beta_{532}$, $\beta_{1064}$) and extinction ($\alpha_{355}$, $\alpha_{532}$) coefficients. (b) Particle depolarization ratios ($\delta_{355}$, $\delta_{532}$, $\delta_{1064}$) and lidar ratios (LR$_{355}$, LR$_{532}$). (c) Fluorescence backscattering ($\beta_F$), fluorescence capacity ($G_F$) and the particle effective radius ($r_{eff}$). (d) Number ($N$, $N^S$), surface area ($S$, $S^S$) and volume ($V$, $V^S$) concentrations obtained by inversion of $3\beta + 2\alpha$ observations (symbols) and calculated from the fluorescence backscattering (lines) by using the mean conversion factors defined in Eq. (2).</figcaption></figure>

The depolarization ratio decreases with wavelength from $9 \pm 1.5\%$ at 355 nm to $1.5 \pm 0.3\%$ at 1064 nm, and the lidar ratio at 532 nm significantly exceeds the corresponding value at 355 nm ($80 \pm 12$ and $50 \pm 7.5$ sr respectively), which is typical for aged smoke (Müller et al., 2005). The multiwavelength observations were inverted to determine the particle effective radius $r_{eff}$, number, surface area and volume concentrations for seven height bins inside the smoke layer. The effective radius $r_{eff}$ in Fig. 5c increases through the layer from 0.15–0.2 µm simultaneously with the
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Observation of smoke fluorescence on 12–13 September 2020 at 21:00–03:00 UTC. (a) Vertical profiles of the aerosol backscattering coefficient $\beta_{532}$ and fluorescence backscattering coefficient $\beta_F$. (b) Aerosol backscattering coefficient $\beta_{532}$, together with smoke backscattering coefficient $\beta_{532}^s$, computed from $\beta_F$ for $G_F = 3.6 \times 10^{-4}$. (c) Surface area concentration of the smoke particles calculated from $\beta_F$ using the respective conversion factor from Eq. (3). Open symbols show the temperature profile measured by the radiosonde at Herstmonceux.</figcaption></figure>

increase in the fluorescence capacity $G_F$ from $2.8 \times 10^{-4}$–$3.6 \times 10^{-4}$. Retrieved values of N, S and V were used for the calculation of the conversion factors (Eq. 2) for each height bin. In the center of the smoke layer (at 4.9 km) the factors are $C_N = 88 \times 10^4 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}$, $C_S = 35 \times 10^4 \text{ µm}^2 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}$ and $C_V = 2.4 \times 10^4 \text{ µm}^3 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}$ (thus, when $\beta_F$ is given in $\text{Mm}^{-1} \text{ sr}^{-1}$, the calculated values of N, S and V are given in $\text{cm}^{-3}$, $\text{µm}^2 \text{ cm}^{-3}$ and $\text{µm}^3 \text{ cm}^{-3}$ respectively). Fluorescence backscattering is proportional to the particle volume concentration, so $C_V$ is not sensitive to the effective radius variation. The conversion factors $C_N$ and $C_S$, on the contrary, depend on the particle size. Figure 5d shows the profiles of N, S and V obtained by the inversion of $3\beta + 2\alpha$ observations (symbols), together with corresponding values ($N^S$, $S^S$, $V^S$) obtained from $\beta_F$, using the mean conversion factors for the seven height bins considered. The volume concentrations V and $V^S$ agree well for all seven height bins. For the surface area concentrations the agreement is still good, but for N and $N^S$ the difference is up to 30%. We need to emphasize that the conversion factors presented are specific for our lidar system (for the interference filter installed in the fluorescence channel). It is worthwhile to mention that the ratio $V/\alpha_{532}$ of the volume concentration V in Fig. 5d to the extinction coefficient $\alpha$ at 532 nm in Fig. 5a, as well as the ratio $S/\alpha_{532}$, is very close to respective extinction-to-volume and extinction-to-surface-area-concentration conversion factors presented for aged wildfire smoke by Ansmann et al. (2021).

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Formation of ice particles at heights above 10 km inside a smoke layer on 11–12 September 2020. Spatiotemporal variations in range-corrected lidar signal at 1064 nm, volume depolarization ratio at 1064 nm and fluorescence backscattering coefficient (in $10^{-4} \text{ Mm}^{-1} \text{ sr}^{-1}$).</figcaption></figure>
<table><caption>**Table 1.** Conversion factors $C_N$, $C_S$ and $C_V$ and fluorescence capacity $G_F$ at height $H$ for five smoke episodes. Volume and surface area concentrations of smoke particles, obtained by the inversion of $\beta + 2\alpha$ lidar observations ($V, S$), are given, together with values calculated from fluorescence measurements ($V^S, S^S$) and using conversion factors (Eq. 3).</caption><thead><tr><th rowspan="2">Date<br>dd/mm/yy</th><th rowspan="2">H<br>km</th><th rowspan="2">C<sub>N</sub><br>10<sup>4</sup> cm<sup>-3</sup><br>(Mm<sup>-1</sup> sr<sup>-1</sup>)<sup>-1</sup></th><th rowspan="2">C<sub>S</sub><br>10<sup>4</sup> µm<sup>2</sup> cm<sup>-3</sup><br>(Mm<sup>-1</sup> sr<sup>-1</sup>)<sup>-1</sup></th><th rowspan="2">C<sub>V</sub><br>10<sup>4</sup> µm<sup>3</sup> cm<sup>-3</sup><br>(Mm<sup>-1</sup> sr<sup>-1</sup>)<sup>-1</sup></th><th rowspan="2">G<sub>F</sub><br>10<sup>-4</sup></th><th colspan="2">V (µm<sup>3</sup> cm<sup>-3</sup>)</th><th colspan="2">S (µm<sup>2</sup> cm<sup>-3</sup>)</th></tr><tr><th>V</th><th>V<sup>S</sup></th><th>S</th><th>S<sup>S</sup></th></tr></thead><tbody><tr><td>23/06/20</td><td>4.9</td><td>88</td><td>35</td><td>2.4</td><td>3.5</td><td>21 ± 4</td><td>19 ± 4</td><td>306 ± 75</td><td>237 ± 60</td></tr><tr><td>11/09/20</td><td>7.5</td><td>75</td><td>28</td><td>2.0</td><td>3.9</td><td>7.6 ± 1.6</td><td>8.7 ± 1.6</td><td>111 ± 25</td><td>111 ± 25</td></tr><tr><td>14/09/20</td><td>6.0</td><td>90</td><td>34</td><td>2.3</td><td>3.7</td><td>6.4 ± 1.3</td><td>6.1 ± 1.3</td><td>94 ± 25</td><td>78 ± 20</td></tr><tr><td>17/09/20</td><td>6.8</td><td>21</td><td>21</td><td>2.3</td><td>2.9</td><td>8.0 ± 1.6</td><td>7.8 ± 1.6</td><td>73 ± 18</td><td>100 ± 25</td></tr><tr><td>20/09/20</td><td>4.9</td><td>33</td><td>22</td><td>2.0</td><td>4.3</td><td>2.7 ± 0.5</td><td>2.9 ± 0.6</td><td>31 ± 8</td><td>37 ± 9</td></tr></tbody></table>

The conversion factors depend on the smoke composition. To estimate the variation range of $C_N$, $C_S$ and $C_V$, several smoke episodes were analyzed, and corresponding results are presented in Table 1. The table provides the fluorescence capacity $G_F$ and the conversion factors at the heights where $3\beta + 2\alpha$ data could be calculated. Mean values of $\langle C_N \rangle$, $\langle C_S \rangle$ and $\langle C_V \rangle$ derived for these episodes and corresponding standard deviations are

$$ \begin{aligned} \langle C_N \rangle &= (61 \pm 32) \times 10^4 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}, \\ \langle C_S \rangle &= (28 \pm 6.4) \times 10^4 \text{ µm}^2 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}, \\ \langle C_V \rangle &= (2.2 \pm 0.2) \times 10^4 \text{ µm}^3 \text{ cm}^{-3} (\text{Mm}^{-1} \text{ sr}^{-1})^{-1}. \end{aligned} \quad (3) $$

Table 1 shows also the volume and surface area concentrations of the smoke particles obtained from the inversion of $3\beta + 2\alpha$ observations ($V, S$) and calculated from $\beta_F$ ($V^S, S^S$) using the conversion factors in Eq. (3). Standard deviations of $V^S$ and $S^S$ from corresponding values of $V_{3\beta+2\alpha}$ and $S_{3\beta+2\alpha}$ are 10% and 25% respectively.

The mean conversion factors in Eq. (3) are now used to estimate the smoke microphysical properties inside the cloud, assuming in addition that the predominant contribution to the fluorescence is provided by the smoke. Table 2 summarizes the number, surface area and volume concentrations of smoke particles inside the ice clouds, estimated from fluorescence measurements for four episodes considered in this paper. On 12–13 September 2020, the smoke layer with high fluorescence and low depolarization ratio at 1064 nm (below 4%) was observed during the whole night inside the 2.0–5.0 km height range. The cirrus cloud occurred above 11 000 m also during the whole night. Figure 6a presents vertical profiles of the aerosol $\beta_{532}$ and fluorescence $\beta_F$ backscattering coefficients. Fluorescence backscattering shows a maximum at 3.5 km, but it is detected even inside the cloud. The smoke backscattering coefficient $\beta_{532}^S$, computed from $\beta_F$ for $G_F = 3.6 \times 10^{-4}$, agrees well with $\beta_{532}$ inside the 2.0–10.0 km height range (Fig. 6b). The height profile of the surface area concentration of the smoke particle $S^S$, calculated from $\beta_F$ using the respective conversion factor in Eq. (3), is shown in Fig. 6c. In the smoke layer, $S^S$ is up to

60 µm² cm⁻³, while in the center of the cloud at 12–13 km height the average value of $S^S$ is 1.6 ± 0.4 µm² cm⁻³. Corresponding values of number and volume concentrations in the cloud center are 3.5 ± 1.8 cm⁻³ and 0.13 ± 0.013 µm³ cm⁻³.

The temperature in the cloud ranged from about –50 to almost –70 °C and was –68 °C at cirrus top in Fig. 6b where ice nucleation usually starts. We applied the immersion freezing INP parameterization of Knopf and Alpert (2013) for leonardite (a standard humic acid surrogate material) and assume that this humic compound represents the amorphous organic coating of smoke particles. The INP parameterization for smoke particles is summarized for lidar applications in Ansmann et al. (2021). The selected parameterization allows the estimation of the INP concentration as a function of ambient air temperature (freezing temperature), ice supersaturation, particle surface area and time period for which a certain level of ice supersaturation is given. We simply assume a constant ice supersaturation of around 1.45 during a time period of 600 s (upwind phase of a typical gravity wave in the upper troposphere). The temperature at cirrus top height is set to –68 °C and the aerosol surface area concentration to 2.0 µm² cm⁻³ as indicated in Fig. 6c. The obtained INP concentrations of 1–10 L⁻¹ for these meteorological and aerosol environmental conditions can be regarded as the predicted number concentration of ice crystals nucleated in the cirrus top region. Ice crystal number concentrations of 1–10 L⁻¹ are typical values in cirrus layers when heterogeneous ice nucleation dominates (typical values of INP concentrations and supersaturation are discussed, for example, in Sullivan et al., 2016; Ansmann et al., 2019b, 2021; Engelmann et al., 2021). It should be mentioned that the required very high ice supersaturation levels of close to 1.5 (ice supersaturation of 1.1–1.2 is sufficient in the case of mineral dust particles) are still lower than the threshold supersaturation level of > 1.5 at which homogeneous freezing starts to dominate. At low updraft velocities around 10–25 cm s⁻¹, as usually given in gravity waves in the upper troposphere (Barahona et al., 2017), heterogeneous ice nucleation very likely dominates the ice production when cirrus evolves in detected aerosol layers.
<table><caption>**Table 2.** Number $N^S$, surface area $S^S$ and volume $V^S$ concentrations of smoke particles inside the ice cloud at height $H$ estimated from fluorescence measurements by applying the conversion factors in Eq. (3) for four measurement sessions.</caption><thead><tr><th>Date<br>dd/mm/yy</th><th>Time<br>UTC</th><th>H<br>km</th><th>10<sup>-4</sup> Mm<sup>-1</sup> sr<sup>-1</sup></th><th>β<sub>F</sub><br>cm<sup>-3</sup></th><th>N<sup>S</sup><br>μm<sup>2</sup> cm<sup>-3</sup></th><th>S<sup>S</sup><br>μm<sup>3</sup> cm<sup>-3</sup></th><th>V<sup>S</sup><br>μm<sup>3</sup> cm<sup>-3</sup></th></tr></thead><tbody><tr><td>12/09/20</td><td>01:20–01:45</td><td>10.5</td><td>0.32</td><td>20 ± 10</td><td>9 ± 2.3</td><td>0.7 ± 0.15</td><td></td></tr><tr><td>12–13/09/20</td><td>21:00–03:00</td><td>12.5</td><td>0.06</td><td>3.5 ± 1.8</td><td>1.6 ± 0.4</td><td>0.13 ± 0.03</td><td></td></tr><tr><td>17/09/20</td><td>22:45–23:45</td><td>8.5</td><td>6.5</td><td>400 ± 200</td><td>180 ± 45</td><td>14 ± 3</td><td></td></tr><tr><td>24–25/11/20</td><td>18:00–06:00</td><td>8.0</td><td>0.013</td><td>0.8 ± 0.4</td><td>0.36 ± 0.09</td><td>0.03 ± 0.006</td><td></td></tr></tbody></table>

## 3.3 Ice formation inside the smoke layers

During September 2020 we observed several episodes with ice cloud formation inside of smoke layers. One such episode occurred on 11–12 September 2020 and is shown in Fig. 7. The height-time display of the fluorescence backscattering coefficient reveals the smoke layer in the 5.0–10.0 km height range. Inside this layer, we can observe a short time interval of 15 min with a strongly increased depolarization ratio around 10.5 km height (red spots), indicating ice cloud formation. Figure 8 shows vertical profiles of the aerosol backscattering coefficients $\beta_{355}$, $\beta_{532}$ and $\beta_{1064}$, as well the particle depolarization ratios $\delta_{355}$, $\delta_{532}$ and $\delta_{1064}$, for two temporal intervals. The first interval (23:00–00:30 UTC) is prior to ice cloud formation, and the second one (01:20–01:45 UTC) covers the ice occurrence period. The depolarization ratios at all three wavelengths were < 5% below 6 km height. Above that height $\delta_{355}$ significantly increased reaching the value of 10% at 7 km (Fig. 8b), which is indicative of a change in the particle shape (from spherical to irregular shape). The fluorescence capacity also changed with height, being about $G_F = 4.5 \times 10^{-4}$ at 5.5 km, and it decreases to $3.5 \times 10^{-4}$ by 8 km. The profile of $\beta_{532}^S$ shown in Fig. 8c is calculated assuming $G_F = 4.0 \times 10^{-4}$, and it matches well the profile of $\beta_{532}$ for the whole height range. The aerosol layer at 10.5 km is thus a pure smoke layer. Ice formation at 10.5 km (Fig. 8d–f) leads to a significant increase in $\beta_{532}$, while $\beta_{532}^S$ (or the respective fluorescence backscatter coefficient $\beta_F$) remains low and at the same level as observed below the cirrus layer, i.e., below 10 km height. The depolarization ratios at all three wavelengths increases to typical cirrus values around 40%. The temperature at 10.5 km is about –50 °C, and the surface area concentration of the smoke particles inside the cloud, estimated from $\beta_F$, is about $10 \mu m^2 cm^{-3}$ (see Fig. 8f, thin blue line). For these temperature and aerosol conditions, smoke INP concentrations of $1–10 L^{-1}$ are yielded for ice supersaturation values even below 1.4 (1.38–1.4) and updraft duration of 600 s. When comparing Fig. 8c and f at cirrus level it seems to be that ice nucleation on the smoke particles widely depleted the smoke INP reservoir.

Another case of ice formation in the smoke layer was observed on 17–18 September 2020. Strong smoke layers occurred in the 5.0–9.0 km height range as shown in Fig. 9. During the period from 22:30–00:00 UTC, the depolarization

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Formation of ice particles at 10–11 km height inside a smoke layer on 11–12 September 2020. Vertical profiles of (a, d) the aerosol backscattering coefficients $\beta_{355}$, $\beta_{532}$ and $\beta_{1064}$; (b, e) the particle depolarization ratios $\delta_{355}$, $\delta_{532}$ and $\delta_{1064}$; and (c, f) $\beta_{532}$, together with backscattering coefficient of smoke $\beta_{532}^S$, calculated from fluorescence backscattering $\beta_F$ assuming $G_F = 4.0 \times 10^{-4}$. Panel (f) shows also the smoke surface area concentration $S^S$ of the smoke particles calculated from $\beta_F$ by applying the respective conversion factor in Eq. (3). Results are given for the time intervals 23:00–00:30 and 01:20–01:45 UTC, which are prior to and during ice cloud formation at 10.5 km height. The temperature profile measured by the radiosonde at Herstmonceux is shown with open symbols in panel (c).</figcaption></figure>
<figure><img src="image_12.png" /></figure>

**Figure 9.** Formation of ice particles at heights above 8 km inside the smoke layer on 17–18 September 2020. Spatiotemporal variations in range-corrected lidar signal at 1064 nm, volume depolarization ratio at 1064 nm and fluorescence backscattering coefficient (in $10^{-4}$ Mm$^{-1}$ sr$^{-1}$).

tion increased at 8.5 km height, indicating ice formation. Vertical profiles of the particle parameters prior to and during ice formation are shown in Fig. 10. The $\beta_{532}^s$ calculated for $G_F = 3.5 \times 10^{-4}$ matches well with $\beta_{532}$ below 6.9 km and above 8.0 km (Fig. 10c), but inside the 7.0–8.0 km height range $\beta_{532} > \beta_{532}^s$, meaning that $G_F$ was decreased. The depolarization ratio in the 7.0–8.0 km height range shows some enhancement (Fig. 10b): in particular, $\delta_{532}$ increased from 10%–12%. Cloud formation at 8.5 km (Fig. 10d) led to a significantly smaller increase in the depolarization ratio compared to the case on 11–12 September. Prior to the cloud formation the values of $\delta_{1064}$, $\delta_{532}$ and $\delta_{355}$ at 8.5 km were 3%, 10% and 13% respectively (Fig. 10b), and in the cloud corresponding depolarization ratios increase up to 9%, 15% and 20%. The reason is probably that the signal averaging period from 22:45–23:45 UTC includes a cloud-free section. Three gravity waves obviously crossed the lidar field site and triggered ice nucleation just before 23:00 UTC, 15–30 min after 23:00 UTC and around midnight (00:00 UTC). The temperature at cloud top at about 8.5–8.6 km height was close to –35 °C. For this high temperature and the high particle surface area concentration of $200 \ \mu m^2$ cm$^{-3}$ (see Fig. 10d, thin blue line), smoke INP concentrations of $1–10$ L$^{-1}$ are

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

**Figure 10.** Formation of ice particles at 8.5–8.6 km height inside a smoke layer on 17 September 2020. Vertical profiles of (a) the aerosol backscattering coefficients $\beta_{355}$, $\beta_{532}$ and $\beta_{1064}$; (b) the particle depolarization ratios $\delta_{355}$, $\delta_{532}$ and $\delta_{1064}$; and (c, d) $\beta_{532}$, together with backscattering coefficient of smoke $\beta_{532}^s$, calculated from fluorescence backscattering $\beta_F$ assuming $G_F = 3.5 \times 10^{-4}$. Results are given for the time intervals (a–c) 21:30–22:30 UTC and (d) 22:45–23:45 UTC, which are prior to and during ice formation at 8.5 km height. Panel (d) shows also the surface area concentration of the smoke particle $S^s$ calculated from $\beta_F$ by applying the respective conversion factor from Eq. (3). The temperature profile measured by the radiosonde at Herstmonceux is shown with open symbols in panel (c).

yielded for a relatively low ice supersaturation of 1.30–1.33 and an updraft period of 600 s. Again, a depletion of the INP reservoir is visible after the formation of the cirrus layer (see Fig. 6c and d around and above 8.5 km height).

## 4 Conclusion

The operation of a fluorescence channel in the LILAS lidar during strong smoke events in the summer and autumn sea-
sons of 2020 has demonstrated the ability of the fluorescence lidar technique to discriminate ice from smoke particles in atmospheric layers in the UTLS height range in great detail. The fluorescence capacity $G_F$ of smoke particles during this period varied within a relatively small range: 2.5–4.5 × 10⁻⁴; thus the use of the mean value of $G_F$ allows us to estimate the contribution of smoke to the total particle backscattering coefficient. The fluorescence lidar technique makes it possible to estimate smoke parameters, such as number, surface area and volume concentration, in the UTLS height range in a quantitative way by applying conversion factors ($C_N$, $C_S$, $C_V$) which link the fluorescence backscattering and the smoke microphysical properties. These factors, derived from the synergy of multiwavelength Mie–Raman and fluorescence lidar observations, show some variation from episode to episode; however, the use of mean values ($C_N$), ($C_S$) and ($C_V$) allows the estimation of smoke properties in the UTLS height regime with reasonable accuracy. Based on two case studies, we demonstrated that the fluorescence lidar technique provides the unique possibility to characterize the smoke particles and their amount inside cirrus cloud layers. The smoke input parameter (surface area concentration) in smoke INP parameterization schemes that are used to predict ice crystal number concentrations in cirrus generation cells can now be estimated within freshly formed cirrus layers.

The smoke parameters such as fluorescence capacity and conversion factors were derived from observations of aged wildfire smoke transported over the Atlantic in 2020. However, smoke composition depends on many factors, such as burning materials type, flame temperature and environmental conditions; thus the smoke fluorescence properties may also vary. Hence, it is important to perform the measurements for different locations and seasons. The fluorescence backscattering in the UTLS height range is quite weak, so to perform measurements with higher temporal resolution, more powerful lidar systems are needed. A dedicated high-power lidar, LIFE (laser induced fluorescence explorer), will be designed and operated at ATOLL in the frame of OBS4CLIM/ACTRIS-France.

## Appendix A: Estimation of smoke parameters from Mie–Raman and fluorescence lidar measurements

<figure><img src="image_14.png" /><figcaption>**Figure A1.** Flow chart showing the main steps of the procedure of smoke parameter estimation from multiwavelength Mie–Raman and fluorescence lidar measurements. Procedure includes the following steps. (i) For a strong smoke layer the $3\beta + 2\alpha$ data set, derived from multiwavelength Mie–Raman lidar observations, is inverted to the particle number $N$, surface $S$ and volume $V$ density. (ii) Conversion factors $C_N$, $C_S$ and $C_V$ are calculated from Eq. (2) by using the fluorescence backscattering coefficient $\beta_F$. (iii) Different smoke events are analyzed to get mean values of conversion factors $\langle C_N \rangle$, $\langle C_S \rangle$ and $\langle C_V \rangle$. These mean values are used to estimate smoke concentration in weak layers in UTLS and inside cirrus clouds in regular observations. The mean value of smoke fluorescence capacity $\langle G_F \rangle$ allows the estimation of smoke contribution $\beta_{532}^s$ to the total backscattering coefficient $\beta_{532}$.</figcaption></figure>

**Data availability.** Lidar measurements are available upon request (philippe.goloub@univ-lille.fr).

**Author contributions.** IV processed the data and wrote the paper. QH and TP performed the measurements. AA analyzed results of fluorescence measurements and ice formation in smoke layers. PG supervised the project and helped with paper preparation. MK developed software for data processing.

**Competing interests.** The contact author has declared that neither they nor their co-authors have any competing interests.
**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** We acknowledge funding from the CaPPA project funded by the ANR through the PIA under contract ANR-11-LABX-0005-01, the "Hauts de France" Regional Council (project CLIMIBIO) and the European Regional Development Fund (FEDER). The development of the algorithm for analysis of fluorescence observations was supported by Russian Science Foundation (project 21-17-00114). The ESA/QA4EO program is greatly acknowledged for support of observation activities at LOA.

**Financial support.** This research has been supported by the Agence Nationale de la Recherche (grant no. ANR-11-LABX-0005-01), the European Regional Development Fund (FEDER), the "Hauts de France" Regional Council (project CLIMIBIO), and the Russian Science Foundation (project 21-17-00114).

Publisher's note: the article processing charges for this publication were not paid by a Russian or Belarusian institution.

**Review statement.** This paper was edited by Eduardo Landulfo and reviewed by Alexandros D. Papayannis and one anonymous referee.

## References

Ansmann, A., Riebesell, M., Wandinger, U., Weitkamp, C., Voss, E., Lahmann, W., and Michaelis, W.: Combined Raman elastic-backscatter lidar for vertical profiling of moisture, aerosols extinction, backscatter, and lidar ratio, Appl. Phys. B, 55, 18–28, 1992.

Ansmann, A., Baars, H., Chudnovsky, A., Mattis, I., Veselovskii, I., Haarig, M., Seifert, P., Engelmann, R., and Wandinger, U.: Extreme levels of Canadian wildfire smoke in the stratosphere over central Europe on 21–22 August 2017, Atmos. Chem. Phys., 18, 11831–11845, https://doi.org/10.5194/acp-18-11831-2018, 2018.

Ansmann, A., Mamouri, R.-E., Hofer, J., Baars, H., Althausen, D., and Abdullaev, S. F.: Dust mass, cloud condensation nuclei, and ice-nucleating particle profiling with polarization lidar: updated POLIPHON conversion factors from global AERONET analysis, Atmos. Meas. Tech., 12, 4849–4865, https://doi.org/10.5194/amt-12-4849-2019, 2019a.

Ansmann, A., Mamouri, R.-E., Bühl, J., Seifert, P., Engelmann, R., Hofer, J., Nisantzi, A., Atkinson, J. D., Kanji, Z. A., Sierau, B., Vrekoussis, M., and Sciare, J.: Ice-nucleating particle versus ice crystal number concentration in altocumulus and cirrus layers embedded in Saharan dust: a closure study, Atmos. Chem. Phys., 19, 15087–15115, https://doi.org/10.5194/acp-19-15087-2019, 2019b.

Ansmann, A., Ohneiser, K., Mamouri, R.-E., Knopf, D. A., Veselovskii, I., Baars, H., Engelmann, R., Foth, A., Jimenez, C., Seifert, P., and Barja, B.: Tropospheric and stratospheric

wildfire smoke profiling with lidar: mass, surface area, CCN, and INP retrieval, Atmos. Chem. Phys., 21, 9779–9807, https://doi.org/10.5194/acp-21-9779-2021, 2021.

Baars, H., Ansmann, A., Ohneiser, K., Haarig, M., Engelmann, R., Althausen, D., Hanssen, I., Gausa, M., Pietruczuk, A., Szkop, A., Stachlewска, I. S., Wang, D., Reichardt, J., Skupin, A., Mattis, I., Trickl, T., Vogelmann, H., Navas-Guzmán, F., Haefele, A., Acheson, K., Ruth, A. A., Tatarov, B., Müller, D., Hu, Q., Podvin, T., Goloub, P., Veselovskii, I., Pietras, C., Haeffelin, M., Fréville, P., Sicard, M., Comerón, A., Fernández García, A. J., Molero Menéndez, F., Córdoba-Jabonero, C., Guerrero-Rascado, J. L., Alados-Arboledas, L., Bortoli, D., Costa, M. J., Dionisi, D., Liberti, G. L., Wang, X., Sannino, A., Papagiannopoulos, N., Boselli, A., Mona, L., D'Amico, G., Romano, S., Perrone, M. R., Belegante, L., Nicolae, D., Grigorov, I., Gialitaki, A., Amiridis, V., Soupiona, O., Papayannis, A., Mamouri, R.-E., Nisantzi, A., Heese, B., Hofer, J., Schechner, Y. Y., Wandinger, U., and Pappalardo, G.: The unprecedented 2017–2018 stratospheric smoke event: decay phase and aerosol properties observed with the EARLINET, Atmos. Chem. Phys., 19, 15183–15198, https://doi.org/10.5194/acp-19-15183-2019, 2019.

Baars, H., Radenz, M., Floutsi, A. A., Engelmann, R., Althausen, D., Heese, B., Ansmann, A., Flament, T., Dabas, A., Trapon, D., Reitebuch, O., Bley, S., and Wandinger, U.: Californian wildfire smoke over Europe: A first example of the aerosol observing capabilities of Aeolus compared to ground-based lidar, Geophys. Res. Lett., 48, e2020GL092194. https://doi.org/10.1029/2020GL092194, 2021.

Barahona, D., Molod, A., and Kalesse, H.: Direct estimation of the global distribution of vertical velocity within cirrus clouds, Sci. Rep., 7, 206840, https://doi.org/10.1038/s41598-017-07038-6, 2017.

Boers, R., de Laat, A. T., Zweers, D. C. S., and Dirksen, R. J.: Lifting Potential of Solar-Heated Aerosol Layers, Geophys. Res. Lett., 37, L24802, https://doi.org/10.1029/2010GL045171, 2010.

Boose, Y., Sierau, B., García, M. I., Rodríguez, S., Alastuey, A., Linke, C., Schnaiter, M., Kupiszewski, P., Kanji, Z. A., and Lohmann, U.: Ice nucleating particles in the Saharan Air Layer, Atmos. Chem. Phys., 16, 9067–9087, https://doi.org/10.5194/acp-16-9067-2016, 2016.

China, S., Alpert, P. A., Zhang, B., Schum, S., Dzepina, K., Wright, K., Owen, R. C., Fialho, P., Mazzoleni, L. R., Mazzoleni, C., and Knopf, D. A.: Ice cloud formation potential by free tropospheric particles from long-range transport over the Northern Atlantic Ocean, J. Geophys. Res.-Atmos., 122, 3065–3079, https://doi.org/10.1002/2016JD025817, 2017.

DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D., Twohy, C. H., Richardson, M. S., Eidhammer, T., and Rogers, D. C.: Predicting global atmospheric ice nuclei distributions and their impacts on climate, P. Natl. Acad. Sci. USA, 107, 11217–11222, https://doi.org/10.1073/pnas.0910818107, 2010.

DeMott, P. J., Prenni, A. J., McMeeking, G. R., Sullivan, R. C., Petters, M. D., Tobo, Y., Niemand, M., Möhler, O., Snider, J. R., Wang, Z., and Kreidenweis, S. M.: Integrating laboratory and field data to quantify the immersion freezing ice nucleation activity of mineral dust particles, Atmos. Chem. Phys., 15, 393–409, https://doi.org/10.5194/acp-15-393-2015, 2015.

Engelmann, R., Ansmann, A., Ohneiser, K., Griesche, H., Radenz, M., Hofer, J., Althausen, D., Dahlke, S., Maturilli, M.,
Veselovskii, I., Jimenez, C., Wiesen, R., Baars, H., Bühl, J., Gebauer, H., Haarig, M., Seifert, P., Wandinger, U., and Macke, A.: Wildfire smoke, Arctic haze, and aerosol effects on mixed-phase and cirrus clouds over the North Pole region during MO-SAiC: an introduction, Atmos. Chem. Phys., 21, 13397–13423, https://doi.org/10.5194/acp-21-13397-2021, 2021.

Fernald, F. G.: Analysis of atmospheric lidar observations: some comments, Appl. Opt., 23, 652–653, doi.org/10.1364/AO.23.000652, 1984.

Freudenthaler, V., Esselborn, M., Wiegner, M., Heese, B., Tesche, M., and co-authors: Depolarization ratio profiling at several wavelengths in pure Saharan dust during SAMUM 2006, Tellus, 61B, 165–179, 2009.

Fromm, M., Lindsey, D. T., Servranckx, R., Yue, G., Trickl, T., Sica, R., Doucet, P., and Godin-Beekmann, S. E.: The untold story of pyrocumulonimbus, B. Am. Meteorol. Soc., 91, 1193–1209, https://doi.org/10.1175/2010bams3004.1, 2010.

Grawe, S., Augustin-Bauditz, S., Hartmann, S., Hellner, L., Petersson, J. B. C., Prager, A., Stratmann, F., and Wex, H.: The immersion freezing behavior of ash particles from wood and brown coal burning, Atmos. Chem. Phys., 16, 13911–13928, https://doi.org/10.5194/acp-16-13911-2016, 2016.

Haarig, M., Ansmann, A., Baars, H., Jimenez, C., Veselovskii, I., Engelmann, R., and Althausen, D.: Depolarization and lidar ratios at 355, 532, and 1064 nm and microphysical properties of aged tropospheric and stratospheric Canadian wildfire smoke, Atmos. Chem. Phys., 18, 11847–11861, https://doi.org/10.5194/acp-18-11847-2018, 2018.

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.

Hoose, C. and Möhler, O.: Heterogeneous ice nucleation on atmospheric aerosols: a review of results from laboratory experiments, Atmos. Chem. Phys., 12, 9817–9854, https://doi.org/10.5194/acp-12-9817-2012, 2012.

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

Jahl, L. G., Brubaker, T. A., Polen, M. J., Jahn, L. G., Cain, K. P., Bowers, B. B., Fahy, W. D., Graves, S., and Sullivan, R. C.: Atmospheric aging enhances the ice nucleation ability of biomass-burning aerosol, Sci. Adv., 7, eabd3440, https://doi.org/10.1126/sciadv.abd3440, 2021.

Kanji, Z. A., Ladino, L. A., Wex, H., Boose, Y., Burkert-Kohn, M., Cziczo, D. J., and Krämer, M.: Chapter 1: Overview of ice nucleating particles, Meteor. Monogr. Am. Meteorol. Soc., 58, 1.1–1.33, https://doi.org/10.1175/amsmonographs-d-16-0006.1, 2017.

Khaykin, S., Legras, B., Bucci, S., Sellitto, P., Isaksen, L., Tencé, F., Bekki, S., Bourassa, A., Rieger, L., Tawada, D., Jumelet, J., and Godin-Beekmann, S.: The 2019/20 Australian wildfires generated a persistent smoke-charged vortex rising up to 35 km altitude, Commun. Earth Environ., 1, 22, https://doi.org/10.1038/s43247-020-00022-5, 2020.

Klett, J. D.: Lidar inversion with variable backscatter/extinction ratios, Appl. Opt., 24, 1638–1643, 1985.

Knopf, D. A. and Alpert, P. A.: A water activity based model of heterogeneous ice nucleation kinetics for freezing of water and aqueous solution droplets, Farad. Discuss., 165, 513–534, https://doi.org/10.1039/c3fd00035d, 2013.

Knopf, D. A., Alpert, P. A., and Wang, B.: The role of organic aerosol in atmospheric ice nucleation: a review, ACS Earth and Space Chemistry, 2, 168–202, 2018.

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

Müller, D., Mattis, I., Wandinger, U., Ansmann, A., Althausen, A., 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.

Murray, B. J., O'Sullivan, D., Atkinson, J. D., and Webb, M. E.: Ice nucleation by particles immersed in supercooled cloud droplets, Chem. Soc. Rev., 41, 6519–6554, https://doi.org/10.1039/c2cs35200a, 2012.

Ohneiser, K., Ansmann, A., Chudnovsky, A., Engelmann, R., Ritter, C., Veselovskii, I., Baars, H., Gebauer, H., Griesche, H., Radenz, M., Hofer, J., Althausen, D., Dahlke, S., and Maturilli, M.: The unexpected smoke layer in the High Arctic winter stratosphere during MOSAiC 2019–2020, Atmos. Chem. Phys., 21, 15783–15808, https://doi.org/10.5194/acp-21-15783-2021, 2021.

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.

Peterson, D. A., Campbell, J. R., Hyer, E. J., Fromm, M. D., Kablick, G. P., Cossuth, J. H., and DeLand, M. T.: Wildfire-driven thunderstorms cause a volcano-like stratospheric injection of smoke, npj Clim. Atmos. Sci., 1, 30, https://doi.org/10.1038/s41612-018-0039-3, 2018.

Peterson, D. A., Fromm, M. D., McRae, R. H. D., Campbell, J. R., Hyer, E. J., Taha, G., Camacho, C. P., Kablick III, G. P., Schmidt, C. C., and DeLand, M. T.: Australia's Black Summer pyrocumulonimbus super outbreak reveals potential for increasingly extreme stratospheric smoke events, npj Clim. Atmos. Sci., 4, 38, https://doi.org/10.1038/s41612-021-00192-9, 2021.

Reichardt, J., Leinweber, R., and Schwebe, A.: Fluorescing aerosols and clouds: investigations of co-existence, Proceedings of the 28th ILRC, Bucharest, Romania, 25–30 June 2017, https://doi.org/10.1051/epjconf/201817605010, 2017.

Richardson, S. C., Mytilinaios, M., Foskinis, R., Kyrou, C., Payannis, A., Pyrri, I., Giannoutsou, E., and Adamakis, I. D. S.: Bioaerosol detection over Athens, Greece using the laser induced fluorescence technique, Sci. Total Environ., 696, 133906, https://doi.org/10.1016/j.scitotenv.2019.133906, 2019.

Schrod, J., Weber, D., Drücke, J., Keleshis, C., Pikridas, M., Ebert, M., Cvetković, B., Nickovic, S., Marinou, E., Baars, H., Ansmann, A., Vrekoussis, M., Mihalopoulos, N., Sciare, J., Curtius, J., and Bingemer, H. G.: Ice nucleating particles over the East-
ern Mediterranean measured by unmanned aircraft systems, Atmos. Chem. Phys., 17, 4817–4835, https://doi.org/10.5194/acp-17-4817-2017, 2017.

Sullivan, S. C., Morales Betancourt, R., Barahona, D., and Nenes, A.: Understanding cirrus ice crystal number variability for different heterogeneous ice nucleation spectra, Atmos. Chem. Phys., 16, 2611–2629, https://doi.org/10.5194/acp-16-2611-2016, 2016.

Umo, N. S., Murray, B. J., Baeza-Romero, M. T., Jones, J. M., Lea-Langton, A. R., Malkin, T. L., O'Sullivan, D., Neve, L., Plane, J. M. C., and Williams, A.: Ice nucleation by combustion ash particles at conditions relevant to mixed-phase clouds, Atmos. Chem. Phys., 15, 5195–5210, https://doi.org/10.5194/acp-15-5195-2015, 2015.

Venable, D. D., Whiteman, D. N., Calhoun, M. N., Dirisu, A. O., Connell, R. M., and Landulfo, E.: Lamp mapping technique for independent determination of the water vapor mixing ratio calibration factor for a Raman lidar system, Appl. Opt., 50, 4622–4632, 2011.

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. Opt., 41, 3685–3699, 2002.

Veselovskii, I., Whiteman, D. N., Korenskiy, M., Suvorina, A., Kolgotin, A., Lyapustin, A., Wang, Y., Chin, M., Bian, H., Kucsera, T. L., Pérez-Ramírez, D., and Holben, B.: Characterization of forest fire smoke event near Washington, DC in summer 2013 with multi-wavelength lidar, Atmos. Chem. Phys., 15, 1647–1660, https://doi.org/10.5194/acp-15-1647-2015, 2015.

Veselovskii, I., Hu, Q., Goloub, P., Podvin, T., Korenskiy, M., Pujol, O., Dubovik, O., and Lopatin, A.: Combined use of Mie-Raman and fluorescence lidar observations for improving aerosol characterization: feasibility experiment, Atmos. Meas. Tech., 13, 6691–6701, https://doi.org/10.5194/amt-13-6691-2020, 2020.

Veselovskii, I., Hu, Q., Goloub, P., Podvin, T., Choël, M., Visez, N., and Korenskiy, M.: Mie-Raman-fluorescence lidar observations of aerosols during pollen season in the north of France, Atmos. Meas. Tech., 14, 4773–4786, https://doi.org/10.5194/amt-14-4773-2021, 2021.