Tropospheric and stratospheric wildfire smoke profiling with lidar: Mass, surface area, CCN and INP retrieval

We present retrievals of tropospheric and stratospheric height profiles of particle mass, volume, and surface area concentrations in the case of wildfire smoke layers as well as estimates of smoke-related cloud condensation nucleus (CCN) and ice-nucleating particle (INP) concentrations from single-wavelength backscatter lidar measurements at ground and in space. A central role in the data analysis play conversion factors to convert the measured optical into microphysical properties. The set of needed conversion parameters for wildfire smoke are derived from AERONET observations of major smoke events caused 5 by record-breaking wildfires in western Canada in August 2017 and southeastern Australia in January-February 2020. The new smoke analysis scheme is applied to stratospheric CALIPSO observations of fresh smoke plumes over northern Canada in 2017 and New Zealand in January 2020 and to ground-based lidar observation in southern Chile in aged Australian smoke layers in January 2020. These case studies show the potential of spaceborne and ground-based lidars to document large-scale and long-lasting wildfire smoke events in large detail and thus to provide valuable information for climate-, cloud-, and air 10 chemistry modeling efforts performed to investigate the role of wildfire smoke in the atmospheric system.


Introduction
Record-breaking injections of Canadian and Australian wildfire smoke into the upper troposphere and lower stratosphere (UTLS) in 2017 and 2020 caused strong perturbations of stratospheric aerosol conditions in the Northern and Southern Hemisphere. The smoke reached heights up to 23 km (Canadian smoke, 2017) Torres et al., 2020) and more than 30 km (Australian smoke, 2020) Kablick et al., 2020;Khaykin et al., 2020), spread over large parts of the stratosphere, and remained detectable for 6-12 months. Smoke particles influence climate conditions Hirsch and Koren, 2021) by strong absorption of solar radiation and by acting as cloud condensation nuclei (CCN) and ice-nucleating particles (INPs) in cloud evolution processes Knopf et al., 2018). As discussed by Ohneiser et al. (2021), smoke may have even been involved in the complex processes leading to the record-breaking stratospheric ozone-depletion events in the Arctic and Antarctica in 2020 (CAMS, 2021). Recent studies suggest that such major hemispheric perturbations may become more frequent in the fu-A. Ansmann et al.: Lidar-based smoke retrievals ture within a changing global climate with more hot and dry weather conditions (Liu et al., , 2014Kitzberger et al., 2017;Kirchmeier-Young et al., 2019;Dowdy et al., 2019;Jones et al., 2020;Witze, 2020).
Lidars around the world and in space are favorable instruments to monitor and document high-altitude aerosol layers in the troposphere and lower stratosphere over long time periods. This was impressively demonstrated after major volcanic eruptions such as the El Chichón and Mt. Pinatubo events (Jäger, 2005;Trickl et al., 2013;Sakai et al., 2016;Zuev et al., 2019). As main aerosol proxies the measured particle backscatter coefficient and the related column-integrated backscatter are used. These optical quantities allow a precise and detailed study of the decay behavior of stratospheric aerosol perturbations. Furthermore, for volcanic aerosol a conversion technique was introduced to derive climate and air-chemistry-relevant parameters such as particle extinction coefficient and related aerosol optical thickness (AOT), mass, and surface area concentration from the backscatter lidar observations (Jäger and Hofmann, 1991;Jäger et al., 1995;Jäger andDeshler, 2002, 2003). Analogously, such a conversion scheme is needed for the analysis of free-tropospheric and stratospheric wildfire smoke layers but is not available yet. The two major stratospheric smoke events in 2017 and 2020 motivated us to develop a respective smoke-related data analysis concept. The technique covers the retrieval of smoke microphysical properties and the estimation of cloud-relevant aerosol properties such as cloud condensation nuclei (CCN) and ice-nucleating particle (INP) number concentrations. The focus is on backscatter lidar observations at 532 nm, but can easily be extended to 355 and 1064 nm, the other two main laser wavelengths used in atmospheric lidar studies. A preliminary version of the new method was already applied to describe the decay of stratospheric perturbation after the major Canadian smoke injection in the second half of year 2017  and in recent studies of stratospheric smoke observed over the North Pole region with ground-based lidar during the winter half year of 2019-2020 (Ohneiser et al., 2021). The retrieval scheme is easy to handle and applicable to lidar observation from ground and in space and thus can also be used to evaluate measurements acquired by the spaceborne CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) lidar Omar et al., 2009;Kar et al., 2019), CATS (Cloud-Aerosol Transport System aboard the International Space Station, ISS) (Proestakis et al., 2019), and the Aeolus lidar (Reitebuch, 2012;Reitebuch et al., 2020;Baars et al., 2021), which continuously monitor the global aerosol distribution.
For completeness, alternative lidar techniques are available to derive microphysical properties of smoke layers from lidar observations (Müller et al., 1999a(Müller et al., , 2014Veselovskii et al., 2002Veselovskii et al., , 2012. These comprehensive inversion methods were successfully applied to wildfire smoke layers in the troposphere Murayama et al., 2004;Müller et al., 2005;Tesche et al., 2011;Alados-Arboledas et al., 2011;Veselovskii et al., 2015) as well as in the stratosphere (Haarig et al., 2018) and even to a stratospheric volcanic aerosol observation (Mattis et al., 2010). However, this sophisticated approach needs lidar observation at multiple wavelengths of very high quality and is strongly based on directly observed particle extinction coefficient profiles which are not easy to obtain, especially not during the final phase of major stratospheric perturbations. The lidar inversion technique can sporadically provide valuable information about the relationship between the optical and microphysical properties of observed aerosol layers and thus can be used to check the reliability of applied sun-photometer-based conversion factors as shown in Sect. 5.5.
The article is organized as follows. An introduction into the complex chemical, microphysical, morphological, and optical properties of wildfire smoke and the ability of these particles to influence ice formation in clouds is given in Sect. 2. In Sect. 3, we provide an overview of the methodological concept, i.e., the way we derive the microphysical and cloud-relevant smoke properties from height profiles of the particle backscatter coefficient. A central role in the data analysis is played by conversion factors Ansmann, 2016, 2017). The way we determined the smoke conversion factors from Aerosol Robotic Network (AERONET) (Holben et al., 1998) sun photometer observations is described in Sect. 4. Section 5 presents the results of the AERONET correlation analysis and the derived set of conversion parameters for fire smoke as obtained from respective observations with AERONET sun photometers in North America, southern Africa, southern South America, and Antarctica. A summary of the studies and an uncertainty analysis is given in Sect. 6. Case studies of observations of stratospheric Australian smoke with ground-based Raman lidar in Punta Arenas, Chile, in January 2020 and of fresh tropospheric smoke with the CALIPSO lidar over the United States in September 2020 are discussed in Sect. 7. Concluding remarks are given in Sect. 8.

Wildfire smoke characteristics
The development of a smoke-related conversion method is a difficult task because of the complexity of smoke chemical, microphysical, and morphological properties. To facilitate the discussions in the next sections, a good knowledge of smoke characteristics is necessary and provided in this section. The overview is based on the smoke research and discussions presented by Fiebig et al. (2003), Müller et al. (2005Müller et al. ( , 2007a, Dahlkötter et al. (2014), China et al. (2015China et al. ( , 2017, Knopf et al. (2018), and Mishchenko (2018, 2020).

Chemical, physical, and morphological properties
First of all, the types of fires, e.g., flaming versus smoldering combustion, the fuel type (burning material), and the combustion efficiency at given environmental and soil moisture conditions determine the initial chemical composition and size distribution of the smoke particles injected into the atmosphere. Burning of biomass at higher temperatures, during flaming fires, generates smaller particles than smoldering fires (Müller et al., 2005). In forest fires, the flaming stage is usually followed by a longer period of smoldering fires.
Smoke particles from forest fires are largely composed of organic material (organic carbon, OC) and, to a minor degree, of black carbon (BC). The BC mass fraction is typically < 5 % (Dahlkötter et al., 2014;Yu et al., 2019) but may reach values of 10 %-15 % in cases of complex mixtures of anthropogenic haze with domestic, forest, and agricultural fire smoke . Biomass burning aerosol also consists of humic-like substances (HULIS), which represent large macromolecules (Mayol-Bracero et al., 2002;Schmidl et al., 2008a, b;Fors et al., 2010;Graber and Rudich, 2006). The particles and released vapors within biomass burning plumes undergo chemical and physical aging processes during long-range transport. There is strong evidence from lidar observations that smoke particles grow in size during the aging phase (Müller et al., 2007a). Processes that lead to the increase in particle size are hygroscopic growth of the particles, gas-to-particle conversion of inorganic and organic vapors during transport, condensation of large organic molecules from the gas phase in the first few hours of aging, coagulation, and photochemical and cloud-processing mechanisms. The lidar observations are in agreement with modeling studies of Fiebig et al. (2003), who used the theory of particle aging processes described by Reid and Hobbs (1998). Condensational growth dominates the increase in particle size in the first 2 d after emission of a plume. Thereafter coagulation in the increasingly diluted plumes becomes the dominating process. A significant shift of the particle size distribution indicated by an increase in the number median radius from about 0.2 µm shortly after emission to about 0.35 µm after 6 d of travel was found in several cases of Canadian smoke by Müller et al. (2007a). The aging effect has to be considered in the retrieval of smoke conversion factors. We distinguish fresh and aged smoke observations in Sect. 5. Dahlkötter et al. (2014) analyzed aircraft in situ measurements of a smoke layer advected from North America and observed over Germany at 10-12 km height in September 2011 and found, in agreement with many other airborne in situ observations, an almost monomodal size distribution of smoke particles with a pronounced accumulation mode (particles with diameters from roughly 200 to about 1400 to 1800 nm). A distinct coarse mode was absent.
The black-carbon-containing smoke particles showed coating thicknesses of roughly 50-220 nm and shell-to-core diameter ratios of typically 2-3. Dahlkötter et al. (2014) assumed a concentric-spheres core-shell morphology for the strongly-light-absorbing BC core and further assumed purely-light-scattering coating material (i.e., no absorption by the shell) in their analysis of the airborne in situ observations. The authors emphasized that their core-shell model is an idealized scenario because the BC cores of combustion particles are fractal-like or compact aggregates and BC can be mixed with light-scattering material in different ways, including, e.g., surface contact of BC with the light-scattering components, full immersion of BC in the light-scattering component, or immersion of the light-scattering components in the BC aggregate. A process that can produce near-surface BC morphology is coagulation of almost bare BC aggregates with BC-free particles. Condensation of secondary organic or inorganic aerosol components on BC particles can result in particles either with core-shell morphology (concentric or eccentric) or with near-surface BC morphology. All these possible morphology features must be considered in the discussion and estimation of the smoke optical properties and of the potential of smoke particles to serve as INP (Sects. 2.2 and 3.1).
Changes in the morphology (size, shape, and internal structure) of smoke particles and their internal mixing state (e.g., soot particle coating) are ongoing during long-range transport. As China et al. (2015) pointed out, freshly emitted soot particles, i.e., BC particles, are typically hydrophobic, lacy fractal-like aggregates of carbonaceous monomers and become hydrophilic as a result of coating and other aging processes. Lace soot undergoes compaction upon humidification. All these effects lead to an increased ability of smoke particles to serve as CCN with increasing long-range travel time.
Soot compaction (and collapse of the core structures) changes also the scattering and absorption cross sections depending on the refractive index, the monomer diameter, and the structural details. Many publications dealing with the optical properties became available in recent years (China et al., 2015;Mishchenko, 2018, 2020;Kahnert, 2017;Yu et al., 2019;Gialitaki et al., 2020). Liu and Mishchenko (2018) mentioned that their model considers 11 different model morphologies ranging from bare soot to completely embedded soot-sulfate and soot-brown carbon mixtures. In agreement with earlier studies, they found that for the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Thus absorption increases as soot accumulates more coating material during long-range transport. As a general finding of the modeling studies, the absorption enhancement is a complex function of many factors such as the size and shape of the soot aerosols, the mixing state, the location of soot within the host, and the amount and composition of the coating material. All these facts make it necessary to distinguish between fresh smoke (<2.5 d after injection) and aged wildfire smoke (>2.5 d of long-range transport) in our attempt to determine smoke conversion parameters. As already mentioned, smoke particles after long-range transport seem to be favorable CCN because they become increasingly hydrophilic during aging. In contrast to the impact of smoke on cloud droplet formation, the characterization of their influence on ice nucleation is rather difficult. The link between ice nucleation efficiency and particle chemical and morphological properties and the ongoing modifications of the properties during long-range transport is largely unresolved (China et al., 2017). However, it is widely assumed that the ability of smoke particles to serve as INP mainly depends on the organic material (OM) in the shell of the coated smoke particles (Knopf et al., 2018). BC is not considered to be an important contributor to immersion freezing (Möhler et al., 2005;Ullrich et al., 2017;Schill et al., 2020;Kanji et al., 2020), which is assumed to be the preferred heterogeneous ice nucleation mode. Knopf et al. (2018) present a review on the role of organic aerosol (OA) and OM in atmospheric ice nucleation. A unique feature of OA particles is that they can be amorphous and can exist in different phases, including liquid, semisolid, and solid (or glassy) states, in response to changes in temperature (T ) and relative humidity (RH) (Koop et al., 2011;Zobrist et al., 2008;Knopf et al., 2018). At low temperatures, e.g., in the UTLS region, where the atmospheric temperature can be as low as 180 K, it is conceivable to assume that the particles are in a glassy state. Most of the secondary organic aerosol particles are solid above 500 hPa (about 5 km) according to modeling studies and for temperatures < 240 K .
It has been shown that humic and fulvic matter can act as deposition nucleation and immersion freezing INPs (Wang and Knopf, 2011;Rigg et al., 2013;Knopf and Alpert, 2013;Knopf et al., 2018). Furthermore, these macromolecules can undergo amorphous phase transition under typical tropospheric conditions (Wang et al., 2012;Slade et al., 2017) similar to the processes we assume the organic coating of the smoke particles experience.
Aerosol particles serving as INPs usually provide an insoluble, solid surface that can facilitate the freezing of water (Knopf et al., 2018). Deposition ice nucleation is defined as ice formation occurring on the INP surface by water vapor deposition from the supersaturated gas phase. Although, recent studies suggest that deposition ice nucleation can be the result of pore condensation freezing, where homogeneous ice nucleation occurs at lower supersaturation in nanometersized pores (David et al., 2019;Marcolli, 2014). When the supercooled smoke particle takes up water or its shell deliquesces, immersion freezing can proceed, where the INP immersed in an aqueous solution can initiate freezing (Knopf et al., 2018;Berkemeier et al., 2014). Finally, if the smoke particle becomes completely liquid (and no insoluble part within the particle is left), homogeneous freezing will take place at temperatures below 235 K (Koop et al., 2000).
However, in reality, at given air mass lifting conditions, the ice nucleation process can be very complex. The time that solid OM needs for transition to a more liquid state, termed as humidity-induced amorphous deliquescence, can range from several minutes to days at temperatures low enough for ice formation (Mikhailov et al., 2009;Berkemeier et al., 2014;Knopf et al., 2018). Thus the phase change (as function of T and RH) can be longer than typical cloud activation time periods (governed by the updraft velocity), potentially inhibiting full deliquescence and allowing the OA or the organic coating to serve as INP. When amorphous OA or OM are involved in ice nucleation, the condensed-phase diffusion processes within OA particles will most probably govern the ice nucleation pathway (Wang et al., 2012).
The following potential scenarios of atmospheric ice nucleation are uniquely attributable to the presence of amorphous OM. (1) Ice formation in the glassy region may be due to ice nucleation on the solid organic particle, i.e., deposition ice nucleation. (2) During partial deliquescence, a residual solid core is coated by an aqueous shell, and immersion freezing may proceed. (3) At full deliquescence RH, where the particles are completely liquid (and contain no solid soot fragments), homogeneous freezing will occur at temperatures below about 238 K. (4) The presence of a glassy phase in disequilibrium with surrounding water vapor (e.g., cloud activation at fast updrafts as discussed below) may suppress or initiate ice nucleation beyond the homogeneous ice nucleation limit (Berkemeier et al., 2014;Knopf et al., 2018). A slower updraft velocity allows for more time for deliquescence to proceed, potentially resulting in full deliquescence of the OA particle at warmer and drier conditions compared to when a faster updraft is active. Therefore, the same OM can be present in different phase states under the same atmospheric thermodynamic conditions (i.e., T and relative humidity over ice RH i ), resulting in different ice nucleation pathways and corresponding ice nucleation rates. OA particle size or coating thickness can also impact the rate and atmospheric altitude of the organic phase change, as larger particles or thicker coatings require more time to reach full deliquescence (Charnawskas et al., 2017). There are many more peculiarities of amorphous OM that make INP parameterization and prediction efforts very complicated as discussed in detail by Knopf et al. (2018).
Since amorphous smoke OA may take up water and partially deliquesce, resulting in an aqueous solution at possibly subsaturated conditions, we apply the water-activity-based immersion freezing (ABIFM) parameterization Alpert and Knopf, 2016) and homogeneous ice nucleation parameterization by Koop et al. (2000). ABIFM derives the number of INPs per volume of air for a given time period, when T , RH, and particle surface area s are known (see Sect. 3.1.1). A deposition ice nucleation scheme based on classical nucleation theory is outlined in addition (Sect. 3.1.3) to cover the potential pathway of glassy smoke particles to serve as INPs. Again, T , RH, and s are input in the INP estimation.
To demonstrate the prediction or retrieval of smoke INP profiles from lidar observations in Sect. 7, we apply two example OA model systems serving as surrogates of amorphous organic smoke particles. One is based on a macromolecular humic or fulvic acid that undergoes amorphous phase transitions in response to changes in RH and T (Wang et al., 2012) and free-troposphere long-range-transported particles that possess an organic coating acting as INPs (China et al., 2017).

Methodological background: microphysical properties from backscatter coefficients
The goal of the study is to provide a set of conversion parameters that permits the estimation of smoke microphysical properties from particle backscatter coefficients measured at 532 nm. A smoke observation with ground-based lidar at Punta Arenas, in southern Chile, is shown in Fig. 1  . We will use this measurement as a case study in Sect. 7.1 and will apply all conversion procedures to this observation. The methodological background of the conversion of optical into microphysical particle properties is given by Ansmann (2016, 2017). It is out of the scope of this article to present a detailed approach of how an aerosol layer can be unambiguously identified and classified as a smoke layer. In case of single-wavelength backscatter lidars, backward trajectory analysis is the main tool to identify smoke layers and link them to the most probable fire source region. In the case of modern aerosol lidars equipped with polarization-sensitive channels and aerosol and molecular backscatter channels at several wavelengths, favorable conditions are given to identify smoke layers based on the complex set of available information on particle backscatter and extinction coefficients, depolarization ratio, and lidar ratio Müller et al., 2005;Tesche et al., 2011;Burton et al., 2012Burton et al., , 2015Giannakaki et al., 2015;Giannakaki et al., 2016;Prata et al., 2017;Haarig et al., 2018;Hu et al., 2019;Adam et al., 2020;Ohneiser et al., 2020Ohneiser et al., , 2021. However, an unambiguous and accurate quantification of the smoke fraction or contribution to the measured optical backscatter and extinction properties and the separation of smoke and soil dust fractions remains difficult. Soil dust may have been injected together with the smoke by the hot fires.
Regarding the separation of smoke and dust fractions by means of the polarization lidar technique (Tesche et al., 2009(Tesche et al., , 2011Nisantzi et al., 2014), we have to distinguish two branches. As long as the smoke-containing layers occur at low altitudes (in the lower and middle troposphere up to 5-7 km height), we can apply the traditional approach to determine the smoke fraction in dust-smoke mixtures by as-suming a low smoke depolarization ratio of <0.05 and a high mineral dust depolarization ratio of 0.31. In the lower and middle troposphere, aging of the smoke particles is usually fast, including the development of a spherical shape of the aged smoke particles. Furthermore, most of the smoke particles are liquid (at least the shell) at comparably high temperatures and moisture levels. All this leads to a low smoke depolarization ratio at all laser wavelengths from 355 to 1064 nm (Haarig et al., 2018).
However, if the smoke is lifted directly into the upper troposphere and lower stratosphere (UTLS), the smoke properties and aging features may be significantly different. With increasing height, and thus decreasing temperature, water vapor content, and amount of condensable gases, the aging process slows down and the smoke particles become partly glassy. These effects seem to prohibit the development of a perfect spherical shape of the shells. As a consequence, the depolarization ratio can be as high as 0.15-0.2 at 532 nm at greater heights (Burton et al., 2015;Haarig et al., 2018;Hu et al., 2019;Ohneiser et al., 2020). However, we also observed low smoke depolarization ratios in the UTLS region (Ohneiser et al., 2021). Thus, in the case of UTLS smoke observations, the dust-smoke separation technique cannot be used. We have to assume that smoke layers are dominated by smoke (smoke fraction > 0.9) in the UTLS regime, and the soil dust fraction can be neglected at these heights.
To obtain height profiles of smoke in terms of volume concentration v(z), surface area concentration s(z), particle number concentrations n 50 (z), considering all particles with radius > 50 nm, and the large-particle number concentration n 250 (z), considering particles with particle radius > 250 nm, we have the following four basic relationships: (1) s(z) = c s Lβ(z) , (2) n 250 (z) = c 250 Lβ(z) , (3) with the particle backscatter coefficient β(z) at height z and the extinction-to-backscatter or lidar ratio L. The needed conversion factors c v , c s , c 250 , and c 50 and the extinction exponent x for 532 nm are obtained from the analysis of AERONET observations during situations dominated by wildfire smoke. The results of our smoke-related AERONET data analysis are presented in Sect. 5. An important input parameter is the smoke lidar ratio L, required to obtain the smoke extinction coefficient σ = Lβ in the first step of the conversion procedure. As discussed in the review of Adam et al. (2020), the smoke lidar ratio can vary from 25 to 150 sr at 532 nm. However, most studies show that the 532 nm lidar ratio is typically in the range of 70 sr ± 25 sr. For 355 nm, lidar ratios were mostly found around 75 ± 25 sr for fresh smoke and 55 ± 20 sr for aged smoke. Table 1 provides an overview of the large range of smoke lidar ratios. Aged smoke shows a characteristic L ra- tio of L 355 nm /L 532 nm < 1. This feature allows a clear unambiguous identification of smoke layers after long-range transport (Müller et al., 2005;Noh et al., 2009;Nicolae et al., 2013;Ohneiser et al., 2020). The reason for the large spectrum of lidar ratios is the complex smoke properties (size, shape, composition) as discussed in Sect. 2. Extended discussions on smoke lidar ratios can be found in Nicolae et al. (2013), Haarig et al. (2018), and Adam et al. (2020).
We recommend to use a lidar ratio of 55 sr for 355 nm and 70 sr for 532 nm for aged smoke if there is no possibility to obtain actual lidar ratio information from Raman lidar Veselovskii et al., 2015;Haarig et al., 2018;Ohneiser et al., 2020Ohneiser et al., , 2021 or High Spectral Resolution Lidar (HSRL) observations Burton et al., 2015), or in the way Prata et al. (2017) proposed in the case of the CALIPSO lidar to estimate the lidar ratio of smoke layers embedded in clear air. For fresh smoke, an appropriate value for the lidar ratio seems to be 70-80 sr at both wavelengths.
From the obtained values of v, s, and n 50 further relevant parameters can be calculated. The smoke mass concentration m is given by with ρ the density of the smoke particles. Li et al. (2016) investigated different smoke aerosols in the laboratory by burning of different straw types and found densities of 1.1 to 1.4 g cm −3 for the produced smoke particles. For organic particles ρ OM was about 1.05 ± 0.15 g cm −3 , and for ρ EC (elemental carbon) they yielded 1.8 g cm −3 . Chen et al. (2017) reviewed the smoke research in China and concluded that the smoke particle density is 1.0-1.9 g cm −3 . Thus in cases with 2 %-10 % of BC the overall smoke particle density should be in the range of 1.0-1.3 g cm −3 . The particle concentration n 50 is a good aerosol proxy for aerosol particles serving as cloud condensation nuclei (CCN), The CCN concentration is a strong function of updraft speed and thus water supersaturation S w . The number concentration n 50 roughly indicates the CCN concentration for weak updrafts and frequently observed low water supersaturations of S w = 0.2 %. Water supersaturation values may be in the range of 0.4 %-0.7 % in strong updrafts. Then the CCN concentration is a factor of about 2 higher than n 50 . In the case of free-tropospheric and stratospheric smoke, we assume that the relative humidity in the smoke plumes is typically < 60 % so that the derived n 50 values represent the number concentrations for dry aerosol particles, required in the CCN estimation. The estimation of CCN concentration in cases with high relative humidity and corresponding aerosol water-uptake effects is described in Mamouri and Ansmann (2016).
The particle concentration n 250 indicates the reservoir of favorable INPs and is even used as input in dust-INP parameterizations (DeMott et al., 2015). However, in the case of smoke the input parameter in the INP retrieval is the surface area concentration s,  The INP concentration is a function of s, the ice supersaturation S i (which occurs during lifting processes), and temperature T . Details of the complex INP parameterization are given in Sect. 3.1. Finally, information on smoke particle number concentrations (n 50 , n 250 ) and surface area concentration s at stratospheric heights is of interest in studies of heterogeneous formation of polar stratospheric clouds (PSCs). A significant increase in smoke aerosol particle concentration may have a sensitive impact on the evolution of PSCs and their microphysical properties (Voigt et al., 2005;Hoyle et al., 2013;Engel et al., 2013;Zhu et al., 2015).
In order to use the developed smoke retrieval formalism presented here in the case of backscatter lidars operated at single wavelengths of λ = 355 or 1064 nm backscatter lidars, we need to estimate the respective backscatter coefficient at 532 nm in the first step. The 532 nm backscatter profiles within smoke layers may be estimated by using typical smoke color ratios β(532 nm)/β(λ). This aspect is further discussed in Sect. 6.

INP parameterization
As discussed in Sect. 2.2, the estimation of INP concentrations is challenging due to the chemical complexity of the smoke aerosol. The parameterizations introduced in this section cover the OM-related ice nucleation for the temperature range in the upper troposphere (< −40 • C). Only for these low temperatures, organic smoke particles may be able to influence ice nucleation in the atmosphere. In the following, we present procedures to compute INP concentrations for immersion freezing, deposition ice nucleation, and homogeneous freezing.

Immersion freezing
Organic smoke particles that have undergone long-range transport are chemically complex, and INP parameterizations that capture the ice formation rate at upper tropospheric and lower stratospheric conditions (i.e., including subsaturated conditions) are scarce (Knopf et al., 2018). Knopf and Alpert (2013) introduced the water-activity-based immersion freezing model ABIFM, drawn from the water-activity-based homogeneous ice nucleation theory (Koop et al., 2000). Knopf and Alpert (2013) present an ABIFM parameterization for two types of humic compounds based also on experimental data by Rigg et al. (2013) that is valid for saturated and subsaturated atmospheric conditions. For demonstration of our method, we chose to apply the ABIFM for leonardite (a standard humic acid surrogate material) to represent the amorphous organic coating of smoke particles. The ABIFM allows prediction of the ice particle production rate J het,I as a function of ambient air temperature T (freezing temperature), ice supersaturation S i , particle surface area s, and time period t for which a certain level of ice supersaturation S i is given. For demonstration purposes, we simply assume a constant supersaturation period t of 10 min (600 s). Such supersaturation conditions may occur during the upwind phase of a gravity wave.
According to Eqs. (6)-(8) in Alpert and Knopf (2016), we calculate the so-called water activity criterion (Koop et al., 2000) in the first step: The term a w,i in Eq. (8), is the ratio of ice saturation pressure P i to water saturation pressure P w as function of temperature T and can be accurately determined by using Eq. (7) in Koop and Zobrist (2009). When the condensed phase and vapor phase are in equilibrium, the water activity a w is equal to RH w (written as 0.75 if RH w = 75 %) in the air parcel in which ice nucleation takes place (e.g., in a cirrus layer at height z at temperature T ). Relative humidity and temperature values may be available from radiosonde ascents or taken from databases with re-analyzed global atmospheric data. However, the actual RH w and T values during the lifting process (associated with cooling and increase in RH w and decrease in T in the air parcel) remain always unknown and need to be estimated in the studies of a potential smoke impact on cirrus formation. The organic aerosol type leonardite needs a relative humidity over ice RH i of about 130 % or a w = 0.2 at −50 • C to become efficiently activated as INP.
In the next step, the ice crystal nucleation rate coefficient J het,I (in cm −2 s −1 ) is calculated: The particle parameters b and k are determined from laboratory studies for different organic aerosol material. Table 2 contains the parameters for two different natural organic substances (Pahokee peat and leonardite)  which serve as surrogates of the organic coating of the atmospheric smoke particles. Leonardite, an oxidation product of lignite, is a humic-acid-containing soft waxy particle (mineraloid), black or brown in color, and soluble in alkaline solutions. Both substances served as surrogates for humic-like substances (HULIS, Sect. 2.1) in extended immersion freezing laboratory studies Rigg et al., 2013). Organic aerosols containing HULIS are ubiquitous in the atmosphere. We also applied the ABIFM parameterization to aerosol samples representing free-tropospheric aerosol (FTA, China et al., 2017) collected on substrates on the Azores for offline micro-spectroscopic single-particle analysis and ice nucleation experiments. According to backward trajectories, the air masses arriving at the Azores crossed western parts of North America during the main fire season (August-September). FTA showed clear smoke signatures. Note that Eq. (10) delivers strongly fluctuating solutions of J het,I when a w is small, and it delivers robust, less fluctuating J het,I values for a w > 0.1. In the final step, we obtain the number concentration of smoke INP for the immersion freezing mode, with the surface area concentration s of the smoke particles in cm 2 m −3 and the time period t (in seconds) for which constant or almost constant ice supersaturation conditions are given. This can be the time period of a short updraft event (of a few minutes, 120-300 s) or of the lifting period of a gravity wave (>600 s). Long-lasting lifting phases of gravity waves can be up to 20 minutes (1200 s) as our Doppler lidar and radar observations conducted in several field campaigns during the last 10 years indicate.

Homogeneous freezing
Alternatively to smoke particles acting as heterogeneous INPs, we need to consider full deliquescence of smoke particles so that homogeneous freezing comes into play. Following Koop et al. (2000), the ice nucleation rate coefficient for homogeneous freezing is obtained from for 0.26 < a w < 0.34. The INP concentration is then obtained from with the particle volume concentration v in cm 3 m −3 . Homogeneous freezing proceeds at RH i ≈ 150 % at −50 • C (i.e., a w ≈ 0.31), whereas 130 % ( a w = 0.2) is required at −50 • C to activate leonardite-containing particles. Thus at slow ascent conditions heterogeneous ice nucleation on smoke particles may dominate ice formation in cirrus layers.

Deposition nucleation
Wang and Knopf (2011) provide a simplified parameterization of deposition ice nucleation (DIN) based on classical nucleation theory that describes the DIN efficiency of humic and fulvic acid compounds as a function of ambient temperature T and the humidity parameters RH i and S i . An alternative DIN parameterization is provided by, e.g., Hoose et al. (2010). A detailed description of the approach presented here is given in Sect. 3.6 in Wang and Knopf (2011), and thus only a brief introduction is given in the following.
The INP efficiencies are expressed as a function of the contact angle , which describes the relationship of surface free energies among the three involved interfaces including water vapor, ice embryo, and INP. is parameterized as a function of RH i (Eq. 8 in Wang and Knopf, 2011).
The compatibility parameter m = cos( ) (expressing the match between ice embryo and INP) is then used to determine the so-called geometric factor f g (m ) (Eq. 7 in Wang and Knopf, 2011), the free energy of ice embryo formation F g,het (f g , T , S i ) (Eq. 6 in Wang and Knopf, 2011), and finally the ice crystal nucleation rate J het,D (Eq. 5 in Wang and Knopf, 2011) with the Boltzmann constant k B . The final step is then In terms of the contact-angle-based approach, = 180 • represents the case of homogeneous ice nucleation. The smaller , the greater the propensity of the INP to act as deposition nucleation INP.
At the end of this section it remains to be emphasizes that we put together several INP parameterizations in Sect. 3.1 for demonstration purposes. The research on the smoke impact on atmospheric ice formation is ongoing (Knopf et al.,  China et al. (2017) 2018). Presently, uncertainties in the prediction of J het,I and J het,D for organic aerosols are very high (Wang and Knopf, 2011;China et al., 2017). However, the procedures introduced above allow us to estimate INP concentration profiles for organic aerosols and to study the potential impact of wildfire smoke on ice formation in tropospheric mixedphase and ice clouds. In the upcoming years, strong field activities are required, including comparisons of airborne in situ with lidar observations of smoke INP concentrations as successfully performed in the case of Saharan dust (Schrod et al., 2017;Marinou et al., 2019) and so-called cirrus closure experiments as realized in the case of cirrus formation in pronounced Saharan dust layers (Ansmann et al., 2019b) in order to check the applicability of developed smoke INP parameterizations and to quantify the uncertainties in the INP estimates under real-world meteorological, cloud, and aerosol conditions. A first closure study with respect to smoke-cirrus interaction was recently presented by Engelmann et al. (2020).

AERONET sites and data analysis
The AERONET database (AERONET, 2021) contains unique multiyear climatological data sets of spectrally resolved aerosol optical properties and related underlying microphysical properties of aerosol particles (e.g., size distribution, volume, and surface area concentration). These AERONET products are available in the database for purely marine, dust, biomass-burning smoke, and anthropogenic haze conditions as well as for complex mixtures of these basic aerosol types. We used the advantage of the AERONET database already to derive the conversion parameters for marine and Saharan dust conditions Ansmann, 2016, 2017) and extended the dust-related study later on to many desert dust regions around the world (Ansmann et al., 2019a). Now, we apply the methodology to the wildfire aerosol type.

AERONET sites
The smoke conversion parameters c v , c s , c 50 , c 250 , and x, required to solve Eqs.
(1)-(4), were determined from sun photometer observations at nine AERONET stations, distributed over several continents. Figure 2 shows the considered AERONET stations. The observations at these sites cover the full range of smoke scenarios, from fresh to aged  plumes, for different fire types and burning material, and smoke occurrence in the troposphere and stratosphere. Yellowknife (AERONET site: Yellowknife Aurora) and Churchill in Canada were selected because these AERONET sites were located in the outflow region of major smoke plumes which originated from the record-breaking wildfires in British Columbia Torres et al., 2020), Canada, in August 2017. Strong pyrocumulonimbus (pyroCb) towers (Fromm et al., 2010) developed and lifted enormous amounts of wildfire smoke into the upper troposphere and lower stratosphere (UTLS) from 21:00 UTC on 12 August to 00:30 UTC on 13 August 2017 (Peterson et al., 2018). The smoke observation at Yellowknife and Churchill could be thus well assigned to the time after injection and allowed us to study the change in the smoke conversion parameters as a function of time from 12-18 h to about 5 d after injection.
The AERONET stations at Rio Gallegos (CEILAP-RG), Argentina; Punta Arenas (Punta-Arenas-UMAG), Chile, at the southernmost tip of South America; and Marambio in Antarctica were selected because well-aged smoke layers crossed these stations in January and February 2020 A. Ansmann et al.: Lidar-based smoke retrievals . The smoke originated from strong fires in southeastern Australia and traveled the 10 000 km distance within 8-12 d. Strong pyroCb activity lifted the smoke layers up to UTLS heights, and self-lifting processes (Boers et al., 2010) caused further ascent to heights 10-20 km above the tropopause Kablick et al., 2020;Khaykin et al., 2020). The background AOT levels are clearly below 0.05 at 532 nm at these high northern and southern mid-latitudinal stations, far away from industrialized centers, so that the smoke layers could be clearly identified and dominated the sun photometer observations over many days (Yellowknife, Churchill) and weeks (Punta Arenas, Rio Gallegos, Marambio).
In order to consider several centers of biomass burning of global importance we selected six further AERONET stations. Smoke from exceptionally strong forest fires in the western United States and western Canada was observed over Reno (University of Nevada, Reno), Nevada, and Table  Mountain (Table Mountain, CA), California, from the end of August to mid-October 2020 (in close distance to the fire sources) and allowed the determination of conversion parameters for very fresh and mixtures of fresh and aged North American tropospheric smoke layers.
We downloaded long-term observations performed at the AERONET stations Alta Floresta, Brazil (Amazonian forest fires); Mongu, Zambia, in southern Africa; Mukdahan, Thailand; and Singapore in Southeast Asia to consider observations in key fire areas of global importance. The Mongu data sets consists of sun photometer observations at the Mongu site from 1997-2009 and at the Mongu Inn site from 2013-2019. Fairly constant burning conditions are given at Mongu from July to November of each year. The long-term observations in the Amazon region, southern Africa, and Southeast Asia cover smoldering and flaming fires, fresh and aged smoke layers, and agricultural, grassland, savannah, peat, forest, and bush fires. The selection of these AERONET stations in key burning areas was guided by the smoke study of Sayer et al. (2014).
The AERONET smoke studies are supplemented by multiwavelength lidar observations of smoke conversion parameters. These vertically resolved observations were performed at Punta Arenas, Chile ; Manaus, Brazil (Baars et al., 2012); near Washington, DC (Veselovskii et al., 2015); at Cabo Verde; in the outflow regime of central western African smoke (Tesche et al., 2011), at Leipzig and Lindenberg, Germany Haarig et al., 2018); and on the German icebreaker Polarstern drifting through the high Arctic close to the North Pole during the winter half year of 2019-2020 Ohneiser et al., 2021). The lidar results are shown in Sect. 5.5. The retrieval of the microphysical properties was based on backscatter coefficients measured at 355, 532, and 1064 nm and extinction values at 355 and 532 nm (Müller et al., 1999a, b;Veselovskii et al., 2002), except for the smoke observations over Lindenberg in the summer of 1998. Here, particle backscatter coefficients at six wavelength (355, 400, 532, 710, 800, 1064 nm) and extinction coefficients at 355 and 532 nm were available .

AERONET data analysis
We used the version-3 level-2.0 inversion AERONET products (AERONET, 2021) in the case of the long-term observations in the Amazon region, southern Africa, and Southeast Asia and level-1.5 data in the case of the remaining stations. The reason for using level-1.5 data was to significantly increase the number of available observations in our smoke-related studies. Many observations showing high to very high smoke AOTs could not pass the strict criteria of the AERONET data quality checks and were thus removed from the level-2.0 data set. We compared the level-2.0 AERONET products with the corresponding (reduced) level-1.5 products to guarantee that the used level-1.5 data set was of high quality.
In agreement with the AERONET data analysis of Sayer et al. (2014), we used the fine-mode AOTs stored in the AERONET database. Smoke particles form a well-developed accumulation mode (with sizes up to about 1 µm in radius) and the related optical properties are assigned as fine-mode AERONET products (Sayer et al., 2014). However, as will be discussed in Sect. 5.1, a bimodal distribution (accumulation plus coarse mode) was often retrieved from the AERONET sun and sky observations. This was also pointed out by Sayer et al. (2014). The second mode is probably related to soil, road, and desert dust or marine aerosol in the planetary boundary layer. The comparison with respective lidar observations clearly indicates that smoke produces a pronounced accumulation mode only. A coarse mode is absent. Thus, we computed the smoke-related values of s, v, n 50 , and n 250 from the downloaded size distributions by considering the size classes 1-11 only (covering the accumulation mode and thus the radius range up to 0.9-0.95 µm) and correlated these calculated microphysical values with the fine-mode AOT at 532 nm as stored in the AERONET database to finally obtain the conversion parameters. Details of the computation of s, v, n 50 , and n 250 from the AERONET size distributions can be found in Ansmann (2016, 2017).
We begin the discussion of the AERONET results with an overview of the smoke measurements at Yellowknife and Churchill (stratospheric smoke), Reno and Table Mountain (tropospheric smoke), and at Punta Arenas, Rio Gallegos, and Marambio (stratospheric smoke) in Fig. 3. The downloaded AOT data sets (AERONET, 2021) contain values of fine-mode, coarse-mode, and total AOT for 440, 675, 870, and 1020 nm. The AOT τ for 532 nm is obtained from the 440 nm AOT τ 440 and the Ångström exponent a by τ = τ 440 (440/532) a .
The Ångström exponent a is defined as a = ln(τ 440 /τ 675 )/ ln(675/440) with wavelengths λ of 440 and 675 nm. We separately computed 532 nm AOT for fine-mode, coarse-mode, and total aerosol size distributions by using respective fine, coarse, and total aerosol Ångström exponents. In Fig. 3, the total, i.e., fine-mode plus coarse-mode, AOT is shown. In all other figures below, we exclusively used the fine-mode AOT at 532 nm. In cases with a strong smoke occurrence, the fine-mode fraction is usually > 0.9.
The measurements at Yellowknife and Churchill in Fig. 3a were performed 0.5-2.5 d and 2-5 d after injection of smoke into the UTLS height range over British Columbia, Canada, respectively. The injection took place between 21:00 UTC on 12 August 2017 and 00:30 UTC on 13 August 2017 (Peterson et al., 2018). As can be seen, the first smoke plumes arrived over Yellowknife, Canada, already 12-18 h after injection. The 532 nm AOT reached values of almost 2.5. The smoke plumes traveled southeastward and crossed Churchill about 1.5-4 d later. A maximum AOT of 2.7 was measured over Churchill. At clean background conditions the AOT is about 0.025 to 0.05 at these Canadian AERONET stations. To consider all smoke observations over Yellowknife from 13-15 August 2017 (days 225-227) we set the AOT threshold level to 0.45; i.e., we considered cases with total 532 nm AOT of ≥ 0.45, only, in our conversion study.
Rather strong fires occurred in California during the late summer and early autumn of 2020 (Fig. 3b). Mixtures of fresh and aged smoke were observed over Reno and Table Mountain. We increased the 532 nm total AOT threshold level to 0.6 to avoid a significant impact of urban haze on the wildfire smoke observations and derivation of smoke conversion parameters. The haze-related AOT was about 0.1-0.25. The exclusive use of the AERONET fine-mode products further eliminated the potential impact of non-smoke aerosol such as coarse dust and marine particles on the correlation studies. Figure 3c shows the observations of aged Australian wildfire smoke in southern South America and northern Antarctica. The smoke traveled more than 10 000 km within 8-12 d before reaching our combined lidar and AERONET station at Punta Arenas . The diluted smoke caused 532 nm AOTs mostly between 0.05 and 0.3. Maximum values were close to 0.5. At clean background conditions, the AOT is in the range from 0.025-0.035. In our smoke-related AERONET data analysis, we considered all observations with AOT > 0.05 and again carefully checked that all used cases, even those with low AOT, showed clear and dominating smoke signatures (i.e., a pronounced accumulation mode). We selected the low AOT threshold of 0.05 to have sufficient cases in our conversion study for welldefined aged smoke. For each of the shown AOT observation in Fig. 3 we downloaded the required size distributions and computed the respective column-integrated values of s col , v col , n 50,col , and n 250,col (by considering the size classes 1-11).
To obtain the smoke extinction-to-volume conversion factor c v , required to derive volume and mass concentrations with Eqs. (1) and (5), the ratio of the vertically integrated (column) particle volume concentration v col to the fine-mode 532 nm AOT τ was formed for each individual smoke obser- with the layer mean volume concentration v and the layer mean particle extinction coefficient σ . The introduced layer depth D has no impact on the retrieval of the conversion factors and is only introduced to move from column-integrated values and AOT to more lidar-relevant quantities like concentrations and extinction coefficients. In this study, we set D = 1000 m as in the studies before Ansmann, 2016, 2017). For each smoke observation j (from number j = 1 to J ), available in the AERONET database, we computed c v,j and then determined the mean value, which we interpret as a representative smoke conversion factor, In the same way, the conversion factors c 250 , needed to estimate the large-particle number concentration with Eq. (3), and c s , required in the surface area retrieval with Eq.
(2), were computed: It is noteworthy to emphasize again that only the accumulation-mode size range (radius classes 1-11) was considered in the computation of n 250 and s.
In the retrieval of the conversion parameters required to obtain n 50 (Eq. 4), we used a different approach (Mamouri and Ansmann, 2016). Following the procedure suggested by Shinozuka et al. (2015), we applied a log-log regression analysis to the log(n 50,j )-log(σ j ) data field and determined in this way representative values for c 50 and x that fulfill best the relationship, log(n 50 ) = log(c 50 ) + x log(σ ) .

AERONET results
We begin the result section with a discussion of observed smoke size distributions in Sect. 5.1. The continuous growth of smoke particles during the first days after emission is linked to a continuous change in the conversion factors. Therefore, the conversion parameters are significantly different for fresh and aged smoke. In Sect. 5.2, we then present the results of the AERONET-based correlation analysis, starting with the most simple scenarios of well-defined aged smoke observed over the AERONET stations in southern South America and northern Antarctica. Afterwards, we illuminate the link between the microphysical properties v, s, n 50 , and n 250 and the measured light-extinction coefficient σ for mixtures of fresh and aged smoke in North America (Sect. 5.3) and over the subtropical and tropical stations in South America, southern Africa, and Southeast Asia (Sect. 5.4). In addition, in Sect. 5.5, we compare the AERONET findings with lidar observations of smoke conversion factors. The lidarbased approach is an independent method to determine microphysical properties from measured optical effects and thus provides a favorable opportunity to check the relationship between microphysical and optical properties of smoke layers as obtained from the AERONET analysis.

Smoke particle size distributions: from fresh to aged smoke
As emphasized in Sect. 2, the particle size distribution of smoke particles changes with time during the first days after injection into the atmosphere as a result of particle aging processes (chemical processing, particle collisions, and coagulation). The changing size distribution has a strong influence on the microphysical and optical properties as well as the correlation between v, s, n 50 , and n 250 and the smoke extinction coefficient σ . Figure 4 provides insight into the full range of size distributions of atmospheric smoke particles. The smallest particles found at Alta Floresta indicate rather fresh smoke, probably just a few hours after emission. The size distributions for Yellowknife (measured on 13 August 2017, 23:18 UTC) and Churchill were observed about 20 h and 3.5 d after injection of smoke into the UTLS height region, respectively. Aged smoke after long-range transport over more than 1 week was observed at Punta Arenas (8 d after emission) and Lindenberg (10.5 d after emission). It is obvious that the size distribution is shifted towards larger particles with increasing residence time in the atmosphere. All size distributions are normalized so that the integral over each shown size distribution is one. Lidar observation conducted at Leipzig, 180 km to the southwest of Lindenberg (Haarig et al., 2018), and over Punta Arenas  agree well with the respective AERONET size distributions. The lidar observations corroborate that the smoke size distribution is unimodal. Figure 5 shows unimodal as well as bimodal size distribution in cases clearly dominated by smoke. Similar bimodal size distributions were presented in the smoke study of Sayer et al. (2014). The weak coarse mode may result from aerosols in the boundary layer (marine particles, soil, and road dust). The lidar observations do not show this coarse mode.
To consider the changing smoke size distributions shown in Fig. 4 in the smoke data analysis, it would be desirable  to have conversion parameter sets for fresh, weakly aged, and aged smoke particles. However, in all likelihood such an approach would be impractical and/or unreasonably difficult. As will be discussed below in detail, the majority of AERONET smoke observations close to the fire regions indicate that fresh smoke was usually mixed with enhanced levels of background aerosol which, to a large extent, consists of aged smoke. This regional background aerosol obviously builds up over the fire regions during the long-lasting fire seasons. Therefore, we decided to distinguish just between two different measurement scenarios: (a) aged smoke observations (smoke observed after long-range transport over 5 d and more) and (b) measurements of mixtures of fresh and aged smoke (in the near-range to large fire areas). For these two scenarios we developed conversion parameterizations. Figure 6 shows the relationship between (a) the smoke volume concentration v and the smoke-related extinction coefficient σ , (b) particle surface area concentration s and σ , and (c) the particle number concentration of larger smoke particles n 250 and σ for aged Australian smoke. The correlation between the number concentration n 50 and σ is discussed in Sect. 5.6. As a general impression, a clear relationship between v, s, and n 250 and σ is found, at least up to extinction coefficients of 300 Mm −1 (or 0.3 in terms of the fine-mode AOT at 532 nm). The spread in the data reflects variations in the smoke properties (size distribution, refractive index). However, the relatively low scatter in the data is a sign for large similarities in the smoke properties (observed over several weeks). This may be related to the fact that the flamingfire type prevailed, eucalyptus trees were the main burning material, smoke lifting was always linked to strong pyroCb activity and thus similar lifting features, and the size distributions of aged smoke particles after 8-12 d long-range transport are at all very similar.

AERONET results for aged smoke
The mean relationships between v, s, and n 250 and σ are visualized by straight blue lines. The respective mean conversion factors c v , c s , and c 250 are given as numbers in the different panels and also summarized in Table 3. These mean conversion factors were computed from the data in Fig. 6a, b, and c by using the Eqs. (19), (21), and (20), respectively.

AERONET results for mixtures of fresh and aged
North American smoke Figure 7 presents the correlations between the smoke volume concentration v and the smoke extinction coefficient σ (Fig. 7a) and between the smoke surface area concentration s and the smoke extinction coefficient (Fig. 7b) for North American forest fires. The forests in the western United States and Canada mainly consist of pine, fir, aspen, and cedar trees. The flaming-fire type probably prevailed in August 2017 and August-October 2020. The observations in  Table 3.  by aged smoke (larger particles). Nevertheless, a clear relationship between the computed volume and surface area concentrations and the measured smoke extinction coefficient is given. We used the observations at Yellowknife (1-2 d old stratospheric smoke) and Reno (tropospheric smoke, observed a few hours to several days after injection) to compute the conversion parameters and mean relationships visualized by red solid lines in Fig. 7. The mean values of c v and c s , as given in the figures, were calculated with Eqs. (19) and (21). Only the Yellowknife and Reno data in Fig. 7 were considered in this computation. All mean conversion factors are summarized in Table 3.
The Yellowknife data points (fresh smoke) are close the red lines. This may indicate that the respective conversion factors (given as numbers in Fig. 7) describe predominately fresh and weakly aged North American smoke properties. The blue straight lines (for aged Australian smoke) seem to define the lower limit of the range of values in Fig. 7. Many observations taken at Table Mountain, east of Los Angeles (tropospheric smoke), and at Churchill (2-5 d old stratospheric smoke) are close to the blue lines for aged smoke.

AERONET results for mixtures of fresh and aged Amazonian, African, and Southeast Asian smoke
In this section, we switch from short-term observations of record-breaking and major fire episodes to long-term observations (partly over decades) in key burning areas of global importance. We assume that these long-term observations cover the full range of smoke-property-influencing aspects (smoldering and flaming fires, very different fuel types, short-to long-range smoke transport, and related smoke aging effects). Figure 8 presents the correlations between the computed smoke values of the volume concentration v and surface area concentration s and the smoke-related fine-mode extinction coefficient σ at 532 nm for all four selected subtropical and tropical stations. A relatively strong variability is found for the relationship between the surface area concentration and extinction coefficient in Fig. 8b, and even significant differences between the different data sets (Southeast Asian vs. African and Amazonian observations) are visible. In contrast, a quite narrow distribution of all observations is given for the volume-to-extinction relationship in Fig. 8a. The spread in the data is again widely a function of the size distribution and thus of the age of the smoke layers. As in Fig. 7, the upper part of the data fields is strongly influenced by smaller particles and thus fresh smoke, whereas the lower part is controlled by larger particles and thus aged smoke.
The green straight lines show the mean regression lines for the Mongu, Zambia, data set. The computation of the mean conversion factors is performed in the same way as described in the forgoing sections. We included again the mean regression lines for aged Australian smoke (blue lines) and also for comparably fresh North American smoke (red lines) in Fig. 8. It is obvious that the blue lines for aged smoke indicate the lower boundary of the data range in Fig. 8a and b well. On the other hand, the upper boundary of the data field seems to be less well defined. Obviously many of the observed plumes of tropical and subtropical fires, especially over Zambia and the Amazon region, are just a few hours old, and thus the smoke particles were very small. The smoke particles of the Amazon region, southern Africa, and Southeast Asia are frequently considerably smaller than North American smoke particles (represented by the red lines in Fig. 8).
It is noteworthy to mention that Sayer et al. (2014) analyzed the relationship between the column smoke volume concentrations v col and the 550 nm fine-mode AOT τ 550 for a large number of AERONET stations around the world with strong impact of wildfire smoke and found similar mean values for the ratio v col /τ 550 as given for the extinction-tovolume conversion factor c v in our figures and in the summarizing Table 3. The study of Sayer et al. (2014) includes also Russian stations (Moscow, Tomsk, Yakutsk). We may thus conclude that our conversion parameter set well covers main aspects and characteristics of wildfire smoke layers around the world.

AERONET vs. lidar smoke observations
Lidar provides an independent approach to derive microphysical parameters of smoke and thus to determine the link between the retrieved microphysical and measured optical properties of smoke particles (Müller et al., 1999a(Müller et al., , 2014Veselovskii et al., 2002Veselovskii et al., , 2012. This option provides the favorable opportunity to check the quality and robustness of our results obtained by analyzing the AERONET data. One of the main problems of sun photometer observations is that the entire vertical column is observed so that, e.g., boundarylayer aerosols can be a disturbing factor in the study of lofted tropospheric and stratospheric smoke plumes. These problems are absent in the case of profiling techniques such as lidar. In the case of active remote sensing methods, the optical and microphysical properties are exclusively determined for the smoke layers. However, the uncertainties in the lidar retrievals can be large, and thus the obtained data products can scatter over a wide range just as a function of these uncertainties. In Fig. 9, lidar data sets of smoke observations from 53 • S (Punta Arenas) to 86 • N (North Pole range) are considered. Correlation between v and s values and σ for fresh and aged smoke plumes originating from fires in western Canada, eastern Siberia, southeastern Australia, eastern United States, the Amazon Basin, and central western Africa are shown. The AERONET-derived mean relationship between v, s, and n 250 Figure 10. Relationship between smoke extinction coefficient σ (532 nm) and particle number concentration n 50 for the combined Reno and Yellowknife data set (fresh and aged smoke) and the combined South American and Antarctic data set (aged smoke). and σ for aged, fresh, and the long-term African observations as discussed in the foregoing sections are shown again as blue, red, and green lines.
A large scatter in the lidar-based smoke correlation values is visible in Fig. 9 with data points even below the blue lines and above the green lines. This large scatter is partly related to the specific retrieval methodology and data analysis strategy as well as to varying assumptions in the analysis of the different lidar data packages. The most robust results (less sensitive to input errors) are obtained in terms of surface area concentrations when using the inversion algorithm of Müller et al. (1999a, b). This method was applied to the lidar observations at Praia, Cabo Verde; Manaus, Brazil; and Lindenberg, Germany. The other observations taken at Leipzig, Punta Arenas, and the North Pole region were analyzed by applying the analysis scheme of Veselovskii et al. (2002).
In Fig. 9b, it can be seen that most of the smoke layers observed over Praia (smoke from central western Africa) contain aged smoke particles (the data points are close to the blue line), and only a minor part of the observations indicate fresh smoke plumes (these data points are close to the green line). Many smoke layers contained a mixture of fresh and aged smoke. All the lidar data, representing smoke after longrange transport (Lindenberg, Leipzig, North Pole, Punta Arenas), are close to the blue line for aged smoke or even below this line and thus in good agreement with the AERONETbased correlation studies. From the consistency found in the correlations shown, based on AERONET and lidar observations, we can conclude that the AERONET smoke conversion parameters presented here allow trustworthy retrieval of smoke microphysical properties from backscatter lidar observations. Table 3. Smoke conversion parameters required in the conversion of the particle extinction coefficient σ at 532 nm into particle number concentrations n 50 and n 250 , surface area concentration s, and volume concentration v. The mean values and SD for c v , c s , c 250 , c 50 , and x are obtained from the extended AERONET data analysis. Effective radius r eff information is given as well. The conversion factors are derived from the AERONET observations at Yellowknife (Y), Reno (R), Alta Floresta (AF), Punta Arenas (PA), Rio Gallegos (RG), Marambio (Ma), Mongu (Mo), Mukdahan (Mu), and Singapore (S). The conversion parameters for South America (AF), southern Africa (Mo), and Southeast Asia (Mu, S) consider observations with AOT > 0.9 at 532 nm only.

Recommended smoke parameterization
Observations close to fire 0.16 ± 0.02 3.0 ± 0.6 0.18 ± 0.09 100 ± 50 0.75 ± 0.08 regions (fresh + aged smoke) Observations far away from 0.13 ± 0.01 1.75 ± 0.25 0.35 ± 0.08 17 ± 5 0.79 ± 0.08 fire regions (aged smoke) 5.6 AERONET results: n 50 retrieval Figure 10 shows the correlation between the CCN-relevant particle number concentration n 50 and the extinction coefficient σ for two contrasting smoke data sets, i.e., for the observations of aged Australian smoke and, on the other hand, for the observations of fresh smoke (Yellowstone) and mixtures of fresh and aged smoke (Reno). According to the applied regression analysis, fresh smoke plumes contain much more CCN-relevant small particles (roughly a factor of 3 more) than aged plumes. For a given extinction coefficient of σ = 100 Mm −1 , n 50 is 635 cm −3 (for aged Australian smoke over Punta Arenas), 1900 cm −3 (for North American smoke), and 3200 cm −3 (for Mongu, Zambia). The numbers for n 50 and the extinction exponent x (see Eq. 4) in Fig. 10 and Table 3 are obtained by considering the respective data sets shown in the figure or mentioned in the table in the linear regression analysis described in Sect. 4.2.
6 Summary of AERONET-derived conversion parameters and retrieval uncertainties Table 3 provides an overview of the derived mean conversion parameters for the different AERONET observational data sets, discussed in the foregoing section. Since the smoke size distribution widely controls the derived conversion pa-rameters, we added the information on the effective radius, which is the particle-surface-area-weighted mean radius of the smoke accumulation mode and can be regarded as a typical radius of the observed smoke particles. For aged smoke, the effective radius is largest. It is much lower for the mixtures of fresh and aged smoke. We recommended the use of the two conversion parameter sets in the lower part of Table 3 in the analysis of smoke layers observed with backscatter lidars.
In Table 4, the uncertainties in the input parameters and the smoke retrieval products are listed. The uncertainties in the conversion parameters are estimated from the SD values in Table 3. The relative uncertainties in the required smoke lidar ratio L and smoke particle density ρ follow from the discussions in Sect. 3. Three scenarios of lidar backscatter profiling are compared in Table 4. In the case of a Raman lidar or a HSRL, the determination of the particle backscatter coefficient in clearly identified smoke layers is possible with high accuracy (10 % relative uncertainty) as our experience shows Veselovskii et al., 2015;Haarig et al., 2018;Ohneiser et al., 2020Ohneiser et al., , 2021. In addition, the lidar ratio L is measured with a typical relative uncertainty of around 20 %. In the case of a powerful groundbased elastic-backscatter lidar, the smoke lidar ratio must be estimated in the determination of the extinction coefficient. The lidar ratio is even required as input in the basic determination of the backscatter coefficient profiles. The backscatter profile retrieval may be possible with a relative uncertainty of 15 %. In the case of comparably weak backscatter signals measured from space (e.g., with the CALIPSO lidar), we assume an uncertainty of 25 % in Table 4 in the profiling of the backscatter coefficient. Details of the uncertainties in the CALIPSO aerosol backscatter coefficients are given in Young et al. (2013Young et al. ( , 2018. Finally, the relative uncertainties in the smoke microphysical retrieval products are obtained by error propagation applied to Eqs. (1)-(5) in Sect. 3.
As can be seen in Table 4, the retrieval of volume, mass, and surface area concentrations of detected smoke layers is possible with an overall uncertainty of about 25 %-35 % (for both fresh or near-source smoke and for aged smoke after long-range transport) in the case of Raman lidars or HSRLs, when the smoke lidar ratios are measured. The respective uncertainties are 40 %-50 % when smoke profiling is performed with an elastic backscatter lidar so that the lidar ratio L needs to be estimated. The number concentrations n 50 and n 250 can be only roughly estimated with a typical uncertainty of about 50 %-70 %. Again, the retrieval uncertainties are lowest when measurements are performed with a groundbased Raman lidar or a HSRL. The uncertainties are then of the order of 35 %-50 % in the case of aged smoke.
Uncertainties in the estimates of CCN and INP concentrations are not listed in Table 4. Comparisons with airborne in situ observations of CCN profiles suggest that the uncertainty in the lidar-based CCN estimation is around 50 %, and in extreme cases up to a factor of 2 (−50 % to 100 %) (Düsing et al., 2018;Haarig et al., 2019;Genz et al., 2020). In the case of INP estimations, it is too early for an in-depth uncertainty analysis. A considerable number of dedicated field campaigns and further laboratory studies are needed before a trustworthy quantification of uncertainties in the INP estimation is possible (see also the discussion at the end of Sect. 3.1).
At the end of the section it should be mentioned that the developed method (here for 532 nm) can be applied to single-wavelength 355 and 1064 nm backscatter lidar observations as well. We recommend in these cases to estimate the 532 nm backscatter profiles from the measured 355 or 1064 nm backscatter profiles by using properly estimated smoke color ratios β 532 /β 355 and β 532 /β 1064 (the index denotes wavelength in nm). Extended overviews of observed wavelength dependencies of smoke backscatter coefficients can be found in Burton et al. (2012) and Adam et al. (2020). In a follow-on project, we may repeat the procedure presented here for 532 nm for the wavelength of 355 nm to cover spaceborne 355 nm HSRL lidar observations of the European Space Agency. Such an approach was already presented by Ansmann (2016, 2017) in the case of the marine and Saharan dust types.  Fig. 1) in terms of the smoke extinction coefficient σ and particle mass concentration m. Extinction coefficients were obtained by multiplying the respective backscatter coefficients with a lidar ratio of 95 sr. The errors margins (thin dotted) indicate relative uncertainties as given in Table 4 for the Raman lidar option in the case of aged smoke.

Lidar case studies
We applied the new smoke conversion scheme to two contrasting smoke observations. In Fig. 1, an aged stratospheric Australian smoke layer was shown, observed with an advanced multiwavelength Raman lidar (Polly: portable lidar system) (Engelmann et al., 2016) at Punta Arenas, Chile, in January 2020. This case will be further analyzed in Sect. 7.1. As a second contrasting example, we selected a measurement of the spaceborne CALIPSO lidar over North and South Dakota, USA. A comparably fresh tropospheric smoke layer was detected in September 2020. The smoke originated from strong wildfires in the western part of the United States and Canada. This case study is presented in Sect. 7.2.

Aged Australian smoke in the stratosphere observed with ground-based Raman lidar
In the framework of a multiyear measurement campaign, we monitored the stratospheric perturbation caused by the record-breaking Australian bushfires with a polarization Raman lidar Polly over a full year, starting in January 2020 . A measurement example is shown in Fig. 1. The results obtained by applying the conversions scheme in Sect. 3 are presented in Figs. 11-13. In the first step, we calculated the extinction coefficients from the 532 nm backscatter coefficients by using a smoke lidar ratio of L = 95 sr as measured with the Raman lidar Polly . Then we applied the conversion factor c v in Table 3 for aged smoke to obtain the volume concentration v. By assuming a particle density of ρ = 1.15 g cm −3  Retrieval results for 29 January 2020 in terms of surface area s and particle number concentration n 50 (proxy for CCN) with error margins representing the uncertainties as given in Table 4 for the Raman lidar option in the case of aged smoke.
for the smoke particles, we obtain the mass concentration m shown in Fig. 11. Such a high aerosol pollution level of 15 µg m −3 at heights from 20-26 km height has never been observed in the stratosphere before, not even after major volcanic eruptions (Trickl et al., 2013;Sakai et al., 2016). Stratospheric background levels are of the order of 0.01 µg m −3 Taha et al., 2021). As shown in Fig. 1, the particle depolarization ratio was significantly enhanced as a result of fast lifting by pyroCb clouds over the Australian fire regions . The aging process was obviously not fully completed and the particles were probably glassy. This may explain the deviation from the perfect spherical shape of the particles and the enhanced depolarization ratios (Gialitaki et al., 2020). Figure 12 shows the derived surface area concentration s and the particle number concentration n 50 . Information on number concentrations and surface area at stratospheric heights is of interest, e.g., in PSC and ozone research. A record-breaking ozone depletion was observed in the stratosphere over Antarctica starting in September 2020 (CAMS, 2021). PSC particles play a strong role in this context because they permit the activation of chlorine components (on the surfaces of the PSC particles) which subsequently destroy ozone molecules. Even if we assume a strong decay of the stratospheric smoke perturbation by a factor of 10 or 100 in the Southern Hemisphere (at mid-to high latitudes) from January 2020 to September 2020, and thus a reduction in the smoke number concentration from about 500 cm −3 in Fig. 12 (in the height range from 21 to 25.5 km height) to 50 or even 5 cm −3 , such number concentrations are still high and are in the range of particle concentrations typically observed in PSCs (0.01-10 cm −3 ) (Jumelet et al., 2008(Jumelet et al., , 2009). Smoke particles may be able to serve as nuclei in processes of het- Figure 13. Retrieval results for 29 January 2020 in terms of INP concentrations n INP,I and n INP,D and large particle number concentration n 250 (considering particles with radius > 250 nm). See text for more details of the INP computations in the case of immersion freezing (red profiles) and deposition nucleation (olive profiles). We consider leonardite as the organic aerosol substance (see Table 2). The INP concentrations are estimated by assuming an air parcel lifting period of 600 s (period of supersaturation with a w of 0.175 (low INP numbers) and 0.2 (high INP values)) and ice nucleation temperature of −50 • C. erogeneous nucleation of PSC particles Zhu et al., 2018) and thus may influence PSC microphysical properties. On the other hand, smoke surface area concentrations of around 120-130 µm 2 cm −3 in the stratospheric layer in Fig. 12 are extremely high, and even if the smoke concentration was reduced by a factor of 10 to 100 until September 2020, surface area concentrations of around 10 or around 1 µm 2 cm −3 are still high and partly in the same range of typical surface area concentrations in PSC clouds (Jumelet et al., 2008(Jumelet et al., , 2009) so that at least a weak influence on ozone depletion by providing surface areas for chlorine activation cannot be excluded.
The surface area concentration is also an essential aerosol input parameter in the INP parameterization and thus an important quantity in the research field of aerosol-cloud interaction with focus on mixed-phase-cloud and cirrus formation in the troposphere. INP estimates are shown in Fig. 13. We use the aerosol type parameters for leonardite as given in Table 2 in the calculation of immersion freezing INP (n INP,I , Sect. 3.1.1). The calculations start with the computation of the water activation criterion a w (Eq. 8). Ice nucleation is a strong function of the vertical velocity (lifting of moist air parcels), which leads to ice supersaturation and thus determines a w . In the case study here, we assume realistic upper-tropospheric cirrus formation conditions and ignore in this demonstration of INP number estimation that we observed the smoke layer in the dry stratosphere 10-15 km above the local tropopause. We assumed RH w = 79.85 % and 82.35 % and a temperature T of −50 • C. The corresponding RH i values are around 125 % and 130 %. Homogeneous freezing proceeds in significant numbers at about RH i = 150 % at −50 • C. Thus, for slow air lifting, smoke particles potentially acting as INPs have a good chance to sensitively influence cirrus formation. With these input values for RH w and T , we obtain a w = 0.175 and 0.2. The value for the ice melting point a w,i (Eq. 9) is 0.6235 at −50 • C. Afterwards, we calculated the ice nucleation rate J het,I (Eq. 10) and the INP concentration n INP,I (Eq. 11) by assuming a lifting period of 600 s during which ice supersaturation conditions according to a w of 0.175 and 0.2 are given. We also computed deposition nucleation INP solutions (n INP,D , Sect. 3.1.3) by assuming the same T , RH i , RH w , and S i input parameters together with an overall lifting period of 600 s. Figure 13 shows the results of the n INP,I and n INP,D estimations. A strong dependence on relative humidity and ice supersaturation is visible. Obviously a threshold value of ice supersaturation S i has to be reached and exceeded before efficient immersion freezing in the case of leonardite starts. The estimated deposition nucleation INP concentration is much higher than the immersion-freezing INP values for the assumed atmospheric conditions. The obtained high INP numbers are directly correlated to the large amount of smoke particles. The obtained INP number concentrations are not too uncommon. For example, INP number concentrations reached about 10-100 L −1 in a Saharan dust plume (DeMott et al., 2003). Neglecting any radiative heating effects of the smoke layer and microphysical processes such as sedimentation and competition for water vapor, these results clearly indicate that organic smoke particles can impact ice formation processes in the upper troposphere during favorable moisture conditions and gravity wave activity with updraft phases lasting longer than several minutes.
In Fig. 13, also values for n 250 (large particle fraction) are shown. It is usually assumed that particles with diameters > 500 nm can be regarded as the overall reservoir for INPs (DeMott et al., 2015). Number concentrations of 10-100 cm −3 or 10 000-100 000 L −1 indicate that this reservoir of large smoke particles cannot be depleted when n INP is in the range of 0.1 to 100 L −1 , not even during extended cirrus formation processes lasting several hours.
The competitive process to heterogeneous ice nucleation is homogeneous freezing. If ice supersaturation S i reaches sufficient levels, corresponding to a w of 0.29-0.31, n INP,hom (Eq. 13) would be of the order of 600-700 L −1 for v ≈ 10 −8 cm 3 L −1 (mean value of the 20-26 km layer).
As mentioned, the uncertainty in the INP retrieval is large and is widely related to the current status of our knowledge about smoke INP type characteristics. The lidar input parameters s and v could be obtained with low relative errors of 25 %-35 %. The research on the role of wildfire smoke particles in cirrus and PSC formation is one of the key topics in atmospheric research with focus on aerosol-cloud interaction (Knopf et al., 2018).

North American smoke in the troposphere observed with the CALIPSO lidar
Strong fires occurred in the western United States and western Canada during the late summer of 2020. The smoke even reached Europe (Baars et al., 2021). Figure 14 shows an overflight of CALIPSO from North Dakota down to Texas. Two smoke layer complexes were detected between 5 and 10 km height: one over North and South Dakota and another one over Texas. According to the backward trajectory analysis, the plumes over North and South Dakota originated from fires in western Canada and were observed after 1 d of travel from the main fire area to North and South Dakota. The plumes over Texas originated from heavy Californian fires and were observed after 3-5 d of travel time from the Californian smoke sources. Cirrus formed in the neighborhood of the smoke layers. Figure 15 presents the height profiles of smoke extinction coefficient, mass concentration, surface area concentration, and estimated INP concentration for the smoke layers detected over eastern North and South Dakota. We used a lidar ratio of 70 sr to convert the measured smoke backscatter coefficients into extinction values and applied the conversion parameter set for fresh smoke as recommended in Table 3. A potential influence of multiple scattering was ignored. For dense aerosol layers, multiple scattering can introduce substantial unquantified errors into the CALIPSO lidar retrievals of particle backscatter and extinction coefficients Liu et al., 2011). However, multiple scattering effects may be of the order of about 5 %-10 % (underestimation of the true extinction coefficient by 5 %-10 %) in cases Figure 15. CALIPSO smoke observation in the stratosphere over North and South Dakota on 13 September 2020 in terms of (a) particle extinction coefficient σ and mass concentration m and (b) INP concentration estimates n INP,I for T = −40 • C, a w = 0.2, and two different organic substances (leonardite, LEO, and freetropospheric smoke aerosol, FTA; see Table 2). The lidar-derived input parameter is the shown surface area concentration s. The CALIPSO aerosol backscatter coefficients were downloaded and averaged over the range from 45-48.5 • N (CALIPSO, 2020b) and then multiplied with 70 sr to obtain the extinction coefficients. Error margins (thin dotted lines) show the uncertainties in σ , m, and s as given in Table 4 (fifth column).
with smoke layer optical thickness ≤ 0.5 at 532 nm (Prata et al., 2017). The particle depolarization ratio was > 0.1 and thus indicated the presence of nonspherical smoke particles.
According to Table 4 (fifth column), the uncertainties in the lidar products are higher now compared to measurements with ground-based Raman lidar at Punta Arenas. Relative uncertainties of 40 %-45 % in the extinction coefficient are mainly caused by the lidar ratio assumption. The uncertainties in the mass and surface area concentrations are around 50 %.
In the computation of the immersion-freezing INP concentrations in Fig. 15b, we highlight the impact of the selected organic aerosol type (leonardite, LEO, vs. free-tropospheric aerosol, FTA, Table 2). We assumed a water activity criterion of a w = 0.2 or RH i =130 % at T = −40 • C and again set the period during which ice nucleation was possible at these thermodynamic conditions to t = 600 s. As can be seen, the assumed organic substance can have a very sensitive impact on the estimated INP values. The third organic substance in Table 2 (Pahokee peat) leads to similar INP values as obtained for leonardite. In cirrus research, it is thus essential to know the origin of the smoke and a good knowledge of the burning material to be able to properly characterize the aerosol type involved in the cloud formation studies. We presented a new method that permits the retrieval of tropospheric and stratospheric height profiles of smoke particle mass, volume, surface area, and number concentrations as well as first-order estimates of CCN and INP concentrations from single-wavelength backscatter lidar observations. The developed smoke retrieval method is of special importance for spaceborne backscatter lidars such as CALIPSO and the numerous ground-based lidars permitting high-quality observations of height profiles of the particle backscatter coefficient at 532 nm up to stratospheric heights. The method allows us to characterize smoke microphysical and optical properties even at very low smoke pollution levels and thus during the entire decay phase of long-lasting stratospheric perturbations, from thick smoke plumes to aerosol background conditions. Even if advanced multiwavelength Raman or HSRL observations are available for the characterization of pronounced smoke layers so that the lidar inversion procedure can be applied to obtain the microphysical properties, our method based on conversion factors is useful for comparisons to corroborate the quality of the solutions obtained with advanced multiwavelength lidar systems.
The required conversion factors were determined from AERONET observations. In this approach, we distinguished observations of aged smoke and mixtures of fresh and aged smoke. A crucial task is the estimation of smoke INP concentrations because of the complex characteristics of smoke particles. Now, a consistent methodology is available to characterize wildfire smoke plumes in terms of microphysical and cloud-relevant parameters. This will allow us to study smoke-cirrus interaction as well as the potential impact of smoke particles on PSC formation and ozone depletion in large detail. We applied the new smoke analysis scheme to ground-based as well as spaceborne CALIPSO observations to highlight the potential of single-wavelength lidars (on the ground and in space) to significantly contribute to an extended monitoring and microphysical characterization of tropospheric and stratospheric smoke layers and thus to provide valuable information for climate, cloud, and air chemistry modeling efforts.
Author contributions. The paper was written by AA with contributions (data analysis, methodological concepts) from KO, REM, DAK, IV, HB, and AF. The co-authors RE, CJ, PS, and BB were involved in the field observations and took care of all the smoke measurements with the Polly lidars at Punta Arenas and aboard RV Polarstern.
Special issue statement. This article is part of the special issue "EARLINET aerosol profiling: contributions to atmospheric and climate research". It is not associated with a conference.
Financial support. The authors acknowledge support through the European Research Infrastructure for the observation of Aerosol, Clouds and Trace Gases (ACTRIS) under grant agreement no. 654109 and 739530 from the European Union's Horizon 2020 research and innovation program. We thank AERONET-Europe for providing an excellent calibration service. AERONET-Europe is part of the ACTRIS project. R.-E. M. has been financially supported by the SIROCCO project (grant no. EXCELLENCE/1216/0217) co-funded by the Republic of Cyprus and the structural funds of the European Union for Cyprus through the Research and Innovation Foundation. Thanks are also provided to the ERATOSTHENES Centre of Excellence, which was established after receiving funding by the Republic of Cyprus and the EU H2020 Widespread Teaming program with grant agreement no. 857510 (https://excelsior2020. eu/, last access: 20 January 2021). The field observations at Punta Arenas are partly funded by the German Science Foundation (DFG) project PICNICC with project no. 408008112. The development of the lidar inversion algorithm used to analyze Polly data was supported by the Russian Science Foundation (project 16-17-10241). D.K. acknowledges support by the DOE grant DE-SC0021034. The Polarstern Polly data were produced as part of the international Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) with the tag MOSAiC20192020 and project ID AWI_PS122_00.
Review statement. This paper was edited by Lucia Mona and reviewed by two anonymous referees.