Fluorescence lidar observations of wildfire smoke inside cirrus: a contribution to smoke–cirrus interaction research
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.
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 fluorescence channel into a multiwavelength aerosol lidar (Reichardt 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.
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 Licel 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 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 β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 β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 βF if the relative sensitivity of the channels is calibrated and the nitrogen Raman scattering differential cross-section σR is known. In our research we used cm2 sr−1 at 355 nm from Venable et al. (2011). All β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 β532 is characterized by the fluorescence capacity .
For most of atmospheric particles βF is proportional to the volume of dry matter, while dependence of β532 on particle size is more complicated. As a result, GF depends not only on aerosol type but also on particle size and the relative humidity RH. Uncertainty of βF calculation depends on the chosen value of σ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 βF and GF 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 GF, which we measured during observations in cloudy conditions in the lower troposphere, was about 2 × 10−8. Thus, at least in the measurements with GF above this value, the contribution of optics fluorescence can be ignored.
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 (βF> 1.2 × 10−4 Mm−1 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 β532 and fluorescence βF backscattering coefficients for the period from 02:00–05:30 UTC are shown in Fig. 2a. The fluorescence capacity GF in the center of smoke layer (not shown) was about 4.5 × 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 ± 4 % at 355 nm (δ355) to 4 ± 1 % at 1064 (δ1064). The tropopause height Htr 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 (βF < 0.3 × 10−4 Mm−1 sr−1) and ascended up to 14 500 m, which is above the tropopause. The corresponding vertical profiles of β532 and βF are shown in Fig. 2b. The fluorescence capacity in the center of the layer is about 4.5 × 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 βF = 0.015 × 10−4 Mm−1 sr−1 at 5000 m (near the cloud base) to a minimum value of 0.01 × 10−4 Mm−1 sr−1 at 7000 m inside the cirrus layer. Above the tropopause the fluorescence backscattering increased strongly and reached the maximum (about 0.3 × 10−4 Mm−1 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 × 10−4 to 4.5 × 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, GF was normally close to 4.5 × 10−4. This relatively low range of GF variations allows the estimation of the backscattering coefficient attributed to the smoke particles from fluorescence measurements as follows:
Figure 3d shows the smoke backscattering coefficient , calculated from βF for GF = 4.5 × 10−4, together with β532. The dynamical range of β532 variations is high. To make smoke backscattering visible above Htr, β532 is plotted in expanded scale in Fig. 3d. The values, though strongly oscillating above the tropopause, match the β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β+2α 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 CN, CS and CV, introduced as
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β+2α 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β+2α) are shown in Fig. 5a, while the particle depolarization ratios δ355, δ532 and δ1064 and the lidar ratios at 355 and 532 nm (LR355, LR532) are presented in Fig. 5b. The depolarization ratio decreases with wavelength from 9 ± 1.5 % at 355 nm to 1.5 ± 0.3 % at 1064 nm, and the lidar ratio at 532 nm significantly exceeds the corresponding value at 355 nm (80 ± 12 and 50 ± 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 reff, number, surface area and volume concentrations for seven height bins inside the smoke layer. The effective radius reff in Fig. 5c increases through the layer from 0.15–0.2 µm simultaneously with the increase in the fluorescence capacity GF from 2.8 × 10−4–3.6 × 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 CN = 88 × 104 cm−3 (Mm−1 sr−1)−1, CS = 35 × 104 µm2 cm−3 (Mm−1 sr−1)−1 and CV = 2.4 × 104 µm3 cm−3 (Mm−1 sr−1)−1 (thus, when βF is given in Mm−1 sr−1, the calculated values of N, S and V are given in cm−3, µm2cm−3 and µm3cm−3 respectively). Fluorescence backscattering is proportional to the particle volume concentration, so CV is not sensitive to the effective radius variation. The conversion factors CN and CS, on the contrary, depend on the particle size. Figure 5d shows the profiles of N, S and V obtained by the inversion of 3β+2α observations (symbols), together with corresponding values (NS, SS, VS) obtained from βF, using the mean conversion factors for the seven height bins considered. The volume concentrations V and VS agree well for all seven height bins. For the surface area concentrations the agreement is still good, but for N and NS 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 of the volume concentration V in Fig. 5d to the extinction coefficient α at 532 nm in Fig. 5a, as well as the ratio , 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).
The conversion factors depend on the smoke composition. To estimate the variation range of CN, CS and CV, several smoke episodes were analyzed, and corresponding results are presented in Table 1. The table provides the fluorescence capacity GF and the conversion factors at the heights where 3β+2α data could be calculated. Mean values of 〈CN〉, 〈CS〉 and 〈CV〉 derived for these episodes and corresponding standard deviations are
Table 1 shows also the volume and surface area concentrations of the smoke particles obtained from the inversion of 3β+2α observations (V, S) and calculated from βF (VS, SS) using the conversion factors in Eq. (3). Standard deviations of VS and SS from corresponding values of V3β+2α and S3β+2α 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 β532 and fluorescence βF backscattering coefficients. Fluorescence backscattering shows a maximum at 3.5 km, but it is detected even inside the cloud. The smoke backscattering coefficient , computed from βF for GF = 3.6 × 10−4, agrees well with β532 inside the 2.0–10.0 km height range (Fig. 6b). The height profile of the surface area concentration of the smoke particle SS, calculated from βF using the respective conversion factor in Eq. (3), is shown in Fig. 6c. In the smoke layer, SS is up to 60 µm2 cm−3, while in the center of the cloud at 12–13 km height the average value of SS is 1.6 ± 0.4 µm2 cm−3. Corresponding values of number and volume concentrations in the cloud center are 3.5 ± 1.8 cm−3 and 0.13 ± 0.013 µm3 cm−3.
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 µm2 cm−3 as indicated in Fig. 6c. The obtained INP concentrations of 1–10 L−1 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−1 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−1, 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.
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 β355, β532 and β1064, as well the particle depolarization ratios δ355, δ532 and δ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 δ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 GF = 4.5 × 10−4 at 5.5 km, and it decreases to 3.5 × 10−4 by 8 km. The profile of shown in Fig. 8c is calculated assuming GF = 4.0 × 10−4, and it matches well the profile of β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 β532, while (or the respective fluorescence backscatter coefficient β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 βF, is about 10 µm2 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 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 calculated for GF = 3.5 × 10−4 matches well with β532 below 6.9 km and above 8.0 km (Fig. 10c), but inside the 7.0–8.0 km height range , meaning that GF was decreased. The depolarization ratio in the 7.0–8.0 km height range shows some enhancement (Fig. 10b): in particular, δ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 δ1064, δ532 and δ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 µm2 cm−3 (see Fig. 10d, thin blue line), smoke INP concentrations of 1–10 L−1 are 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).
The operation of a fluorescence channel in the LILAS lidar during strong smoke events in the summer and autumn seasons 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 GF of smoke particles during this period varied within a relatively small range: 2.5–4.5 × 10−4; thus the use of the mean value of GF 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 (CN, CS, CV) 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 〈CN〉, 〈CS〉 and 〈CV〉 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.
Lidar measurements are available upon request (firstname.lastname@example.org).
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.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
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.
This paper was edited by Eduardo Landulfo and reviewed by Alexandros D. Papayannis and one anonymous referee.
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 concentrationin 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., Stachlewska, 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 MOSAiC: 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 severalwavelengths 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., Pettersson, 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., Papayannis, 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 Eastern 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.