Articles | Volume 22, issue 15
Research article
01 Aug 2022
Research article |  | 01 Aug 2022

Retrieving ice-nucleating particle concentration and ice multiplication factors using active remote sensing validated by in situ observations

Jörg Wieder, Nikola Ihn, Claudia Mignani, Moritz Haarig, Johannes Bühl, Patric Seifert, Ronny Engelmann, Fabiola Ramelli, Zamin A. Kanji, Ulrike Lohmann, and Jan Henneberger

Understanding the evolution of the ice phase within mixed-phase clouds (MPCs) is necessary to reduce uncertainties related to the cloud radiative feedback in climate projections and precipitation initiation. Both primary ice formation via ice-nucleating particles (INPs) and secondary ice production (SIP) within MPCs are unconstrained, not least because of the lack of atmospheric observations. In the past decades, advanced remote sensing methods have emerged which provide high-resolution data of aerosol and cloud properties and could be key in understanding microphysical processes on a global scale. In this study, we retrieved INP concentrations and ice multiplication factors (IMFs) in wintertime orographic clouds using active remote sensing and in situ observations obtained during the RACLETS campaign in the Swiss Alps. INP concentrations in air masses dominated by Saharan dust and continental aerosol were retrieved from a polarization Raman lidar and validated with aerosol and INP in situ observations on a mountaintop. A calibration factor of 0.0204 for the global INP parameterization by DeMott et al. (2010) is derived by comparing in situ aerosol and INP measurements, improving the INP concentration retrieval for continental aerosols. Based on combined lidar and radar measurements, the ice crystal number concentration and ice water content were retrieved and validated with balloon-borne in situ observations, which agreed with the balloon-borne in situ observations within an order of magnitude. For seven cloud cases the ice multiplication factors (IMFs), defined as the quotient of the ice crystal number concentration to the INP concentration, were calculated. The median IMF was around 80, and SIP was active (defined as IMFs > 1) nearly 85 % of the time. SIP was found to be active at all observed temperatures (−30 to −5C), with the highest IMFs between −20 and −5C. The introduced methodology could be extended to larger datasets to better understand the impact of SIP not only over the Alps but also at other locations and for other cloud types.

1 Introduction and background

The increase in the Earth's global mean temperature in recent years is unequivocal, yet the extent of a cloud cooling effect remains most uncertain (IPCC2021). The radiative feedback of a cloud is a strong function of the hydrometeor phase (Sun and Shine1994). Mixed-phase clouds (MPCs) consisting of both water phases (liquid and solid) contribute strongly to uncertainties in the radiative feedback. They are thermodynamically unstable (see e.g., Korolev et al.2017) because of the lower vapor pressure with respect to ice than with respect to liquid water. This causes the ice crystals to grow at the expense of the evaporation of cloud droplets, which is referred to as the Wegener–Bergeron–Findeisen process (Wegener1911; Bergeron1935; Findeisen1938). The complexity of phase partitioning adds to difficulties simulating MPCs in models (e.g., McCoy et al.2016; Matus and L’Ecuyer2017) where experimental observations could reduce uncertainties (see e.g., Baumgardner et al.2011; Mahrt et al.2019). In the evolution of an MPC, the ice phase plays an important role as it controls precipitation initiation and consequently cloud lifetime (see e.g., Field and Heymsfield2015; Heymsfield et al.2020). Thus, the ice phase determines not only how a cloud impacts the radiative budget but also for how long. The first ice crystals in a cloud either formed within the cloud by homogeneous or heterogeneous ice nucleation or can be externally introduced by, for example, sedimentation from a higher cloud (seeder–feeder process, e.g., Proske et al.2021; Ramelli et al.2021a, and references therein) or levitated from the ground (blowing snow, especially in mountainous terrain, e.g., Beck et al.2018; Walter et al.2020, and references therein). Approximately below −38C, supercooled droplets can freeze homogeneously or heterogeneously. Above that temperature, heterogeneous ice nucleation is favored on sparsely abundant aerosols called ice-nucleating particles (INPs; see e.g.,Wegener1911; Vali1971; Pruppacher and Klett2010; Murray et al.2012). Until today, a variety of aerosol particles acting as INPs in the atmosphere have been identified such as desert and soil dust, organics from biomass burning, marine or terrestrial biogenic particles, atmospherically aged soot, bacteria, and others (Kanji et al.2017; Huang et al.2021). However, due to the geospatial variability of INP sources, atmospheric INP concentrations feature a high spatiotemporal variability, complicating their quantitative assessment (see e.g., DeMott et al.2010; Kanji et al.2017; Murray et al.2021, and references therein). Understanding heterogeneous ice formation and subsequent ice crystal growth is key to understanding the link between aerosols and precipitation formation and is therefore an important step towards constraining weather and climate predictions (Ansmann et al.2019a; Bühl et al.2019). After first ice crystals are found in a cloud, secondary ice production (SIP) can enhance the ice crystal number concentration (ICNC) by various processes, e.g., fragmentation during ice–ice collisions (e.g., Vardiman1978; Takahashi et al.1995) or splintering during riming (Hallett–Mossop process, active between −8 and −3C, Hallett and Mossop1974), given favorable environmental conditions (Korolev and Leisner2020). The Hallett–Mossop process requires the presence of supercooled liquid droplets. Thus, the Hallett–Mossop process can play an important role, especially over orographic terrain, given that the orographic-forcing-induced updrafts could sustain the availability of liquid water droplets (Lohmann et al.2016). Various field studies have observed ICNC exceeding the ambient INP concentration by several orders of magnitude (see e.g., Koenig1963; Auer et al.1969; Hobbs and Rangno1985, 1990, 1998; Gayet et al.2009; Crosier et al.2011; Stith et al.2011; Crawford et al.2012; Heymsfield and Willis2014; Lawson et al.2015; Lasher-Trapp et al.2016; Ladino et al.2017; Mignani et al.2019; Lauber et al.2021; Pasquier et al.2022). However, the underlying SIP processes are weakly constrained (Korolev and Leisner2020). Therefore, the prevalence of SIP in the atmosphere as well as the environmental conditions for SIP to be active remains uncertain.

Uncertainties in predicting atmospheric INP concentrations and ice multiplication are (partly) related to spatiotemporally limited field observations, because of the large effort needed to obtain field data and their point-like characteristic. Remote sensing techniques can be a suitable solution to overcome the issue by providing continuous data in time and at least one spatial dimension. Starting in the 1930s, remote sensing techniques have evolved to very advanced and reliable instruments today (Wandinger2005). Using lidar (light detection and ranging) and radar (radio detection and ranging) instruments, a broad variety of aerosol and cloud properties can be retrieved (e.g., Ansmann et al.2019a). Lidar systems are key instruments for the investigation of aerosol optical properties and their vertical layering in the atmosphere. The retrieved backscatter and extinction coefficients allow for the determination of physical properties of aerosol particles in the atmosphere, such as size and particle number concentration (see e.g., Müller et al.1999; Veselovskii et al.2005; Mamouri and Ansmann2016). Advanced lidar systems also deliver other products such as extinction-to-backscatter ratio (lidar ratio), linear depolarization ratio (LDR), or Ångström exponent from multiwavelength observations of the backscatter and extinction coefficients that allow the determination of more detailed aerosol information, e.g., aerosol type (see e.g., Burton et al.2012; Groß et al.2013; Baars et al.2017). The accuracy of lidar retrievals has been validated in many studies, with a focus on different properties comparing lidar measurements with in situ observations using aircraft or unmanned aerial vehicles (UAVs) (Ferrare et al.1998; Wandinger et al.2002; Sakai et al.2003; Müller et al.2012; Sawamura et al.2017; Schrod et al.2017; Marinou et al.2019; Haarig et al.2019; Düsing et al.2021). Furthermore, the LDR was utilized to determine the contribution of dust to the observed aerosol (Shimizu et al.2004; Tesche et al.2009; Mamouri and Ansmann2016; Haarig et al.2017). Precise knowledge about the aerosol type is essential to use the lidar retrieval to estimate INP concentrations, owed to the fact that the INP concentration is deduced only from a physical property such as particle number concentration and surface area concentration. However, the parameterizations are often proposed for specific aerosol types such as dust (e.g., Niemand et al.2012; DeMott et al.2015; Ullrich et al.2017; Harrison et al.2019), marine aerosol (e.g., McCluskey et al.2018), soot (e.g., Ullrich et al.2017), and global aerosol (e.g., DeMott et al.2010). The feasibility of predicting INP concentration from lidar profiles has been shown by Mamouri and Ansmann (2016). In lidar INP studies, typically three aerosol type categories are used: mineral dust (Ansmann et al.2003; Mamouri and Ansmann2016; Schrod et al.2017; Haarig et al.2019; Marinou et al.2019), marine aerosol (Mamouri and Ansmann2016; Haarig et al.2019), and continental aerosol (Mamouri and Ansmann2016; Schrod et al.2017; Düsing et al.2018; Marinou et al.2019; Düsing et al.2021). If necessary, these categories can be extended by wildfire smoke (Ansmann et al.2021; Engelmann et al.2021) or volcanic ash (Ansmann et al.2011). Whereas the aerosol sources within the first two categories can be comparably narrowed down by source, continental aerosol sources are much more diverse, featuring (among others) biological, organic, or lifted soil particles and also anthropogenic emitted particles from biomass burning or combustion for example (Kanji et al.2017; Huang et al.2021, and references therein), complicating the retrieval of continental INP concentrations. Thus, finding the most ideal INP parameterization for a given aerosol type is key for minimizing uncertainties in predicted INP concentration. Schrod et al. (2017) proposed calibration factors to optimize the INP parameterization by DeMott et al. (2015) for the retrieval of INP concentration in condensation and deposition mode from dust-dominated air masses. Marinou et al. (2019) showed that separating the dust-carrying air masses into a dust and a continental component and applying the parameterizations of DeMott et al. (2015, 2010), respectively, yields the best agreement between lidar retrieval and in situ observations (basing on Schrod et al.2017). However, Marinou et al. (2019) stated further that “additional measurements are required in order to define the optimum INP parameterizations for non-dust atmospheric conditions (e.g., continental, marine, smoke)”, suggesting that an optimal INP parameterization for the retrieval of INP concentration from continental air masses is yet to be proposed.

Remote sensing instrumentation has recently increasingly been used for cloud observations and is becoming a promising method for retrieving ice-crystal-related parameters (see e.g., Seifert et al.2010; Bühl et al.2013; Zhang et al.2014; Ansmann et al.2019a). It allows for a continuous monitoring of the vertical cloud structure and therefore permits the observation of spatiotemporal cloud evolution (Seifert et al.2010). This is a crucial addition to in situ measurements, which only allow for measurements in a constrained time and height frame (Seifert et al.2010). The ice water content (IWC) can be retrieved solely based on radar measurements (Hogan et al.2006), and ICNC can be estimated by linking the IWC to an assumed particle shape and size distribution (see e.g., Bühl et al.2019). Complementing radar retrievals with size information from a collocated lidar improves the accuracy of cloud property retrieval further (Delanoë et al.2013). In previous studies, remote sensing was used to study (primary) ice nucleation in the atmosphere (see e.g., Sakai et al.2003; Ansmann et al.2008; Eidhammer et al.2010; Seifert et al.2010, 2011; Ansmann et al.2019a; Engelmann et al.2021) and to estimate SIP (see e.g., Auer et al.1969; Luke et al.2021; Sotiropoulou et al.2021). The identification of a specific SIP process from atmospheric in situ and remote sensing observations is a challenge. Lauber et al. (2021) recently showed that recirculation of melted ice crystals can enhance the concentration of small ice crystals in the proximity above the melting layer. Whereas deepening the understanding of individual processes seems best achievable in controlled conditions during laboratory studies (Korolev and Leisner2020), the contribution and magnitude of SIP still needs to be assessed in the (complex interacting) atmosphere. From field observations, the ice multiplication factor (IMF) defined as the ICNC divided by the INP concentration can be utilized to quantify SIP. Whereas the IMF allows quantification of the excess of ice crystals, it cannot be used to directly identify a specific SIP process due to the potential (temporal and spatial) difference between the origin of the first ice crystals and SIP to occur.

In this study we combined a suite of in situ and remote sensing instruments to understand cloud formation and evolution in orographic terrain. In February and March 2019, we performed an intensive field campaign in the Swiss Alps in the region of Davos. In a high valley, a combined lidar–radar system was employed along with ground-based in situ aerosol (including INP) observation (Wieder et al.2022b) and balloon-borne in situ cloud observations (Ramelli et al.2021a, b). A second in situ aerosol site was located on a nearby mountaintop (height difference 1.1 kmMignani et al.2021; Wieder et al.2022b). The near collocation of the lidar beam and the mountaintop site (horizontal displacement 3.65km) is ideally suited for aerosol and thus INP closure as also previously shown by Bedoya-Velásquez et al. (2018), who studied the hygroscopicity of aerosol particles. Based on aerosol data collected over 8 weeks and seven observed cloud events, we address the following. First, we validate the lidar retrieval of aerosol number concentration and surface area concentration, which is the basis for the INP concentration retrieval (Sect. 3.1). Second, based on INP concentrations measured during a Saharan dust event, we assess the accuracy of different dust parameterizations and evaluate their performance (Sect. 3.2.1). Third, we evaluate two INP parameterizations (DeMott et al.2010; Ullrich et al.2017) used in the lidar community for the retrieval of INP concentration from continental air masses (Sect. 3.2.2). Based on our in situ measurements, we propose a calibration factor to optimize one INP parameterization and validate the tuning with the lidar observations (Sect. 3.2.3). Fourth, we validate radar-retrieved IWC and ICNC with balloon-borne in situ observations (Sect. 3.3). Lastly, using the tuned parameterization we present a methodology to obtain estimates of IMFs in MPCs (Sect. 3.4).

Figure 1Overview of the measurement locations: (a) Davos-Wolfgang (black star) in the east of Switzerland (map source: Federal Office of Topography). (b) Measurement locations with the local topography around. The OCEANET container and the lower aerosol monitoring site were located in a high valley at Wolfgangpass (blue dot, WOP, 1631 ma.s.l.). A second aerosol measurement site was located on the mountaintop Weissfluhjoch (orange dot, WFJ, 2693 ma.s.l.), which is located approximately 1.1 km higher than WOP and horizontally displaced by 3.65 km. Lidar product retrieval heights for the comparison with in situ measurements at WFJ and WOP are indicated in orange and blue, respectively, right of the lidar beam (green). The topography was extracted from the digital height model DHM200 from the Federal Office of Topography swisstopo. (c) View from WFJ in direction of WOP. The approximate location of the lidar beam is indicated in green. (d) Placement of the aerosol trailer right of the OCEANET container's lidar and radar at WOP. The tethered balloon used for in situ cloud observations is seen in the front right.

2 Measurement setup and methodology

The RACLETS (Role of Aerosols and CLouds Enhanced by Topography on Snow) campaign took place in the region of Davos, Switzerland (Fig. 1a), in February and March 2019, where an extensive set of aerosol, cloud, precipitation, and snow measurements were conducted (Envidat2019; Walter et al.2020; Mignani et al.2021; Ramelli et al.2021a, b; Lauber et al.2021; Georgakaki et al.2021; Wieder et al.2022b). The two main measurement sites were located on a saddle at the entrance of a high valley (Wolfgangpass, 1631 ma.s.l., hereafter referred to as WOP) and the other on a mountaintop (Weissfluhjoch, 2693 ma.s.l., hereafter referred to as WFJ) as shown in Fig. 1b. At both locations, two similarly equipped aerosol measurement sites were set up. In addition, the remote sensing instrumentation was installed at WOP.

2.1 In situ aerosol and INP measurements

The aerosol measurement sites at WFJ and WOP were previously presented (Wieder et al.2022b; Georgakaki et al.2021; Mignani et al.2021). At both sites, ambient air was sampled through a 46 C heated inlet. The heating was a preventive measure to avoid icing of the outside inlet parts, to evaporate activated cloud droplets, and to sublimate ice crystals. The evaporation of volatile compounds of the aerosol cannot fully be excluded. However, the effect is expected to be minor given the high flow rate through the inlet (300 L min−1), such that the temperature of sampled air was likely below 46 C. Furthermore, the degradation of relevant INPs (mostly biological) should only occur at temperatures above 46 C (Kanji et al.2017; Huang et al.2021) and is hence regarded as unlikely (see Wieder et al.2022b, for further details). Downstream, an aerodynamic particle sizer spectrometer (APS; model 3321, TSI Inc., US) and a scanning mobility particle sizer spectrometer (SMPS; model 3938, TSI Inc., US) recorded aerosol size distributions between approximately 10 nm (electric mobility) and 20 µm (aerodynamic diameter). In this study, electric mobility (SMPS) and aerodynamic diameter (APS) were converted to physical diameter assuming a shape factor χ=1.2 and assuming a particle density ρ=2 g cm−3 (Thomas and Charvet2017). After conversion, the observed size range of the SMPS and the APS covered particles with physical diameters between 10 and 400 nm and between 400 and 15 µm, respectively. The total surface area concentration (hereafter referred to as s) of all aerosols was calculated from the entire size distribution utilizing both SMPS and APS data. The number concentrations of particles with radii 250 nm (hereafter referred to as n250) were obtained by integrating the APS distributions for physical diameters 500 nm. The error in the obtained concentrations is given by ±10 %. As described above, we cannot fully exclude the loss of volatile compounds which could shift the size distributions to smaller sizes, thus underestimating in situ predicted INP concentrations when using size-dependent aerosol number concentrations as an input parameter (Sect. 3.2.3). However, as we discussed above the effect is likely negligible.

Ambient aerosol was collected over a time span of 20 min into pure water (15 mL, W4502-1L, Sigma-Aldrich, US) for offline INP analysis using high-flow-rate impingers (Coriolis® μ, Bertin Instruments, France, 300 L min−1) attached to the end of the inlets (Wieder et al.2022b; Mignani et al.2021). The (immersion mode) INP analysis was done on site utilizing the drop-freezing apparatuses LINDA (Stopelli et al.2014) at WFJ and DRINCZ (David et al.2019) at WOP. In short, the liquid sample is split into small aliquots of the same volume (Va) and cooled down in a cryostat. During cooling a camera mounted above the cryostat takes pictures of the droplets in which frozen droplets appear more opaque than unfrozen ones. After all droplets froze, the fraction of frozen droplets as a function of the cryostat temperature (FF(T)) is retrieved in post-processing from the pictures. Consequently, the INP concentration is calculated according to Vali (1971) as

(1) n INP ( T ) = - ln [ 1 - FF ( T ) ] V a × C ,

with C as the conversion factor from INP concentration found in the sample liquid (in mL−1) to the concentration in ambient air (in StdL−1). The error in the obtained INP concentration is temperature dependent and varies per sample, but extends on average to ±50 %. Further details about the aerosol setups as well as the INP data processing can be found in Wieder et al. (2022b).

2.2 Lidar and radar measurements

During RACLETS, a 35 GHz (Ka-band) cloud radar and the multi-wavelength polarization Raman lidar PollyXT-OCEANET (Engelmann et al.2016, hereafter referred to as PollyXT) of the Leibniz Institute for Tropospheric Research (TROPOS) were deployed at the WOP site (Fig. 1d) as part of the mobile observation platform OCEANET-Atmosphere (Griesche et al.2020). The instruments provided a continuous overview on the temporal evolution of the vertical distribution of aerosol particles and clouds from 8 February to 16 March 2019 and are used in the following for remote-sensing-based retrievals of aerosol and cloud microphysical properties.

PollyXT measured vertical profiles of the particle backscatter coefficient (at 355, 532, and 1064 nm), the particle extinction coefficient with the Raman method (nighttime, 355 and 532 nm), and the particle linear depolarization ratio (LDR, 355 and 532 nm). Based on the PollyXT observations, INP number concentration was calculated following the POLIPHON method as described in Mamouri and Ansmann (2016). Here, only a short description of POLIPHON will be provided. In a first step, the dust and non-dust contributions to the backscatter coefficient are separated using the particle LDR at 532 nm. The backscatter contributions are transferred to extinction contributions using an extinction-to-backscatter ratio (lidar ratio) of 55 sr for dust and 50 sr for continental particles (rural background, pollution). These extinction coefficients σi at height z are converted for each aerosol type i (dust or continental aerosol) to n250,i and si following the relation described in Mamouri and Ansmann (2016).


The additional factor C(z) defined as p0/p(z)×T(z)/T0 converts the retrieved concentrations to standard conditions (Stdcm−3 for n250 and µm2 Stdcm−3 for s) at standard pressure p0 and standard temperature T0 according to the ambient pressure p(z) and temperature T(z) at height z. The meteorological data were taken from the COSMO-1 reanalysis data (see Sect. 2.4). The needed conversion factors (c250,i,cs,i) are aerosol type dependent and are obtained from long-term sun-photometer (Aerosol Robotic Network, AERONET,  Holben et al.1998) datasets from the Sahara (dust, Ansmann et al.2019b) and the Alpine site of Davos (continental aerosol, calculated for the present study) following the method described in Mamouri and Ansmann (2016). The applied extinction-to-number-concentration (particle radii 250 nm) conversion factors are c250,dust=0.19Mm cm−3 for Saharan dust and c250,cont.=0.0828Mm cm−3 for continental aerosol (Alps), and the extinction-to-surface-area conversion factors are cs,dust=2.4×10-12Mm m2 cm−3 for Saharan dust and cs,cont.=2.48×10-12Mm m2 cm−3 for continental aerosol (Alps). The uncertainties are 30 % for n250 and 30 %–50 % for s. Further assessments of the uncertainties are provided in Mamouri and Ansmann (2016) and Haarig et al. (2019). With this method, the vertical profiles of the basic input parameters in various INP parametrizations are derived (see Sect. 3.2). The uncertainty in the obtained INP concentrations is about a factor of 3 (Haarig et al.2019).

The 35 GHz cloud radar was of the type Mira-35 (Görsdorf et al.2015). During RACLETS, the radar was operated in vertical-stare mode. Pulses with a length of 208ns were emitted at a repetition frequency of 6000 Hz, resulting in a vertical resolution of 31.17 m and a maximum unambiguous velocity range of 25.6 m s−1, which spans from −12.8 to 12.8 m s−1. The return signals of the emitted linearly polarized pulses were detected separately in the co- and cross-polarized planes. For both channels, Doppler spectra are derived from Fourier transformations of the return signals from a series of 512 consecutive pulses, corresponding to a Doppler-velocity resolution of 0.05 m s−1. The final temporal resolution of the acquired cloud radar dataset of 10 s is obtained from incoherent averaging of 100 consecutive Doppler spectra.

From the cloud radar's reflectivity the IWC was derived according to Hogan et al. (2006) as

(4) IWC = 10 0.000242 Z T + 0.0699 Z - 0.0186 T - 1.63

with the radar reflectivity (Z) and the ambient temperature (T; see Sect. 2.4). For this approach, the expected error in the retrieval is temperature dependent up to a factor of 2. Furthermore, ICNCs were retrieved from the cloud radar observations with the method described in Bühl et al. (2019). The procedure to derive ICNC from the RACLETS cloud radar observations was described by Ramelli et al. (2021a) and is only explained briefly in the following. Observations of LDR, radar reflectivity, Doppler velocity, and Doppler spectral width are compared with a lookup table of forward modeled cloud radar spectra holding microphysical and observable quantities. In this way, ICNC and the distribution of maximum particle diameter are estimated. The uncertainty of the retrieved values is estimated by varying the measurement values within their measurement errors and thereby retrieving many retrieval results at once. The most common value in the resulting distribution of retrieved values is considered the result and the width of the distribution represents its uncertainty. Within this study, the uncertainty in ICNC is a factor of around 4 (see Ramelli et al.2021a).

2.3 In situ cloud measurements

During the campaign, cloud properties were also measured in situ using the tethered balloon system HoloBalloon (Ramelli et al.2021a, b). On HoloBalloon, a HOLographic Imager for Microscopic Objects (HOLIMO) is installed which images an ensemble of cloud particles in the size range from 6 µm and 2 mm in a well-defined sample volume to obtain phase-resolved cloud properties (Ramelli et al.2020). The captured particles are classified into cloud droplets, ice crystals, and artifacts using a convolutional neural network (Touloupas et al.2020) for particles larger than 25 µm and a decision tree for particles smaller than 25 µm. The differentiation between cloud droplets and ice crystals is only done for particles larger than 25 µm based on the particle shape (circular versus non-circular), whereby all particles smaller than 25 µm are classified as cloud droplets. All ice crystals predicted by the neural network are manually confirmed or reclassified to ensure a high classification accuracy of ice crystals. Following this approach, the ICNC can be estimated with an uncertainty of ±5 % for ice crystals larger than 100 µm and ±15 % for ice crystals smaller than 100 µm (Beck2017). The IWC (in kg m−3) is calculated using the mass–diameter relationship given in Cotton et al. (2013):

(5) IWC = i π 6 ρ ice D i , max 3 × V cloud - 1 for  D max 70 µ m i 0.026 D i , max 2 × V cloud - 1 for  D max > 70 µ m ,

where Dmax is the maximum dimension of the ice crystals (in m), ρice is the effective ice density of small ice (700 kg m−3), and Vcloud is the cloud sampling volume (in m3). Due to the underlying assumptions involved in the calculation of the IWC (e.g., mass–diameter relationship, uncertainties in the effective ice crystal density), the uncertainty in the IWC can be up to a factor of 2 (Beck2017). In a comparison to radar observations (Sect. 3.3), we use IWCs and ICNCs measured with HoloBalloon on 8 March 2019, for which further details on the measurements and the synoptic situation are provided in Ramelli et al. (2021a).

2.4 Cloud temperature data

Temperatures at cloud top and within the clouds were retrieved from the COSMO-1 analysis data at corresponding heights above WOP. The COSMO-1 data were provided by MeteoSchweiz for the processing of the remote sensing data.

2.5 Investigated INP parameterizations

DeMott et al. (2015)Ullrich et al. (2017)Harrison et al. (2019)DeMott et al. (2010)Ullrich et al. (2017)

Table 1Overview of INP parameterizations for immersion freezing mode considered in this study along with the aerosol type they are applicable to, the aerosol property they are based on (n250: aerosol number concentration of particles with radii 250 nm; s: aerosol surface area concentration), the applicable temperature range, and the explicit form equation number in this paper. This overview was adapted from the more extensive summary given by Marinou et al. (2019, Table 1).

 Note that D10 was not primarily developed for predicting continental INP concentrations and includes samples from different field observations in North America, the Amazon, and the Pacific also featuring dust-carrying air masses. It was shown in recent publications (Mamouri and Ansmann2016; Marinou et al.2019) that it is suitable to deduce continental INP concentration from lidar observations.

Download Print Version | Download XLSX

INP concentrations can be obtained from remote sensing data by applying a suitable INP parameterization to retrieved aerosol properties. By now, a wide range of INP parameterizations specific to a certain aerosol type (e.g., dust, continental, marine, soot) and a freezing mode (e.g., immersion or deposition) exist. Here, we compare lidar-estimated INP concentrations with in situ observations in immersion mode at temperatures -20C during times of a Saharan dust event and otherwise continental aerosol over the region of Davos. In Table 1, we present the investigated existing parameterizations for immersion freezing which are applicable to air masses during our observations. The overview is based on the summary of Marinou et al. (2019), which gives detailed information about the parameterizations (Marinou et al.2019, Sect. 2). In the presence of Saharan-dust-carrying air masses, we investigate the predictions of INP concentration (given in StdL−1) from the parameterizations of DeMott et al. (2015) (D15) defined as

(6) n INP , D 15 ( T , n 250 ) = cf × ( n 250 ) 1.25 × e 0.46 × ( 273.16 - T ) - 11.6 ,

with cf being a calibration factor (unity by default); Ullrich et al. (2017) (dust parameterization, immersion freezing mode, U17d) defined as

(7) n INP , U 17 d ( T , s ) = 10 - 9 × s × e - 0.517 × T + 150.577 ;

and Harrison et al. (2019) (H19) defined as

(8) n INP , H 19 ( T , s ) = f K × 10 - 5 × s × 10 i = 0 5 c i × ( T - 273.15 ) i

where fK is the fraction of K-feldspar on the particle surface, c0=-3.25, c1=-0.793, c2=-6.91×10-2, c3=-4.17×10-3, c4=-1.05×10-4, and c5=-9.08×10-7. The latter parameterization is a more recent dust parameterization scaling with the relative contribution of K-feldspar, which was found to correlate most strongly with the INP concentration in their study. Note that temperature T is given in kelvin, and both n250 (given in Stdcm−3) and s (given in µm2 Stdcm−3) refer to dry particle dimensions. In the other cases, which are considered air masses carrying continental aerosol, we investigate the INP concentration prediction based on DeMott et al. (2010) (global parameterization, D10) defined as

(9) n INP , D 10 ( T , n 250 ) = 5.94 × 10 - 5 × ( 273.16 - T ) 3.33 × ( n 250 ) 0.0264 × ( 273.16 - T ) + 0.0033

and Ullrich et al. (2017) (soot parameterization, immersion freezing mode, U17s) defined as

(10) n INP , U 17 s ( T , s ) = 7.463 × 10 - 9 × s × e 0.7667 - 0.8525 × ( T - 273.15 ) - 0.0101 × ( T - 273.15 ) 2 .

The latter was developed on soot aerosol, thus making it suitable to specifically capture anthropogenic contributions to continental aerosol. Notably, the application ranges of all parameterizations presented in Table 1 are applicable mainly at temperatures -15C. Therefore, we extrapolate the parameterizations to −5C in accordance with Marinou et al. (2019).

3 Results and discussion

Atmospheric INP concentrations and IMF are estimated from the lidar and radar measurements. Three main uncertainties are involved in the estimation: (i) uncertainties linked to the measurement of the aerosol extinction coefficient and its conversion to number or surface area concentration, (ii) uncertainties linked to the INP parameterization itself, and (iii) using a parameterization not suitable for the dominant aerosol constituent. In the following, we first validate the aerosol properties as input and consequently determine the most suitable INP parameterization for dust and continental aerosol during RACLETS.

3.1 Aerosol concentration comparison

Figure 2Comparison of aerosol properties observed in situ and with lidar. (a) Time series of aerosol number concentration with radii ≥250nm (n250) measured in situ (black) and with lidar (orange) at WFJ. (b) Same as (a) but measured at WOP (lidar data in blue). (c) Time series of surface area concentration (s) measured in situ (black) and with lidar (orange) at WFJ. (d) Same as (c) but measured at WOP (lidar data in blue). Right of each time series panel, a scatter plot compares the in situ to the lidar observations of the data presented on the left. Uncertainties are indicated by the shading for the time series and at representative data points of the scatter plots in black. Note the lidar retrieval height at WOP was 2052 ma.s.l. (see Fig. 1).


In the following, we investigate the accuracy of the lidar retrieval of aerosol properties with in situ observations. For the comparison of lidar observations to in situ observations at WFJ (mountaintop site) and WOP (high valley site), lidar retrievals were taken from the closest height bins at 2695 and 2052 ma.s.l., respectively (Fig. 1b). In the case of WOP, the lowest bin with complete overlap of the lidar (see e.g., Wandinger and Ansmann2002) was taken for the comparison, which is still around 400 m above the in situ site. The incomplete overlap of the emitted and received beam is a general issue in the comparison of lidar measurements to ground-based in situ observations. The comparison to WFJ, more than 1000 m above the lidar site where complete overlap is surely reached, is a great advantage of the present study.

Figure 2a, b show n250 from in situ observations and lidar retrievals for WFJ and WOP, respectively. Note that the accuracy will be discussed only qualitatively. A quantitative assessment and improvement of the lidar aerosol retrieval is beyond the scope of this publication. For both locations in situ and lidar observations agree qualitatively well, in particular for higher concentrations (dusty conditions). At relatively low aerosol concentrations and in the presence of continental aerosol, a plateauing of the lidar-retrieved aerosol concentrations was observed. The clear atmosphere over the Alps with very low values of the extinction coefficient (<10Mm−1) could be responsible for deviations from the assumed linear relationship of extinction to n250 and s, respectively, (Mamouri and Ansmann2016). The larger diurnal variability between the in situ observation and lidar retrieval at WOP (high valley site) compared to WFJ (mountaintop site) can be explained by the diurnal changes of aerosol concentration near the ground (in situ observations) not affecting the air masses on the lidar retrieval height (height difference approx. 400 m; see Fig. 1b). This difference in height and therefore air mass commonly limits a quantitative conclusion between ground-based in situ observations and remote sensing instruments as the well-mixed boundary layer could at times not extend up to the lowest retrieval height. For the retrieval of s (Fig. 2c, d) the aforementioned observations hold equally true, which is not surprising as the surface area relates to the square of particle radius. However, comparing the retrieval accuracy at WFJ (Fig. 2a, c), a stronger bias of the lidar-retrieved surface area concentrations is apparent, which has also been previously reported by Haarig et al. (2019), based on similar observations carried out at Barbados.

Despite the fact that the comparison in this section was purely qualitative, two important insights are gained: (i) due to the better agreement of in situ observations and lidar retrieval for n250 compared to s, INP parameterizations based on n250 should be preferred. (ii) Comparing ground-based in situ observations (collocated with the lidar instrument) does not allow for a comparison even to the lowest lidar retrieval height. Consequently, we restrict the comparison to in situ INP concentrations in the following to observations made at WFJ (mountaintop site), where the lidar retrieval overlaps the in situ observations.

Figure 3Overview of the Saharan dust event over the Davos region starting on 19 February 2019 and the continental aerosol conditions before. (a) Lidar attenuated backscatter coefficient (β1064nmatt.) measured above WOP (high valley site). Note that the gap in data on 20 February is due to technical maintenance of the lidar system. Retrieval height of WFJ (mountaintop site) is indicated in dashed white. (b) Lidar attenuated backscatter coefficient (β1064nmatt.) retrieved from the lidar height bin closest to WFJ (2695 ma.s.l., Fig. 1b). (c) In situ aerosol number concentrations for particles with radii 250 nm (n250, left axis) at WFJ (dashed orange) and WOP (dashed blue) and INP concentrations at −13C (nINP,-13C, right axis) measured in situ at WFJ. INP samples during continental background are indicated with crosses, and samples during Saharan dust are indicated with hexagons. Uncertainties in INP concentration are indicated for representative samples in black.


3.2 INP concentration comparison

Throughout the campaign, the region around Davos was mostly exposed to continental aerosol (e.g., 14–18 February 2019 as seen in Fig. 3). The situation changed 19–22 February 2019, when the synoptic wind situation promoted transport of Saharan dust from North Africa towards Davos. The arrival of the dust plume is clearly visible in the attenuated backscatter coefficient (β1064nmatt.) of the lidar (Fig. 3a, b). Additionally, an increase by more than an order of magnitude in n250 occurred at both sites on 19 February (in situ observations, Fig. 3c). Before 19 February, the higher variability of n250 at WOP (high valley site) compared to WFJ (mountaintop site) can be explained by local sources and accumulation of aerosol beneath a nighttime inversion along with valley and mountain breezes (Wieder et al.2022b). Over the course of 19 February, n250 at WOP followed the same trend and magnitude as at WFJ and continued doing so in the following days. This drastic change in trend underlines that the region was affected and dominated by the long-range-transported Saharan dust. INP concentrations at −13C measured in situ at WFJ also showed an increase of approximately 1 order of magnitude already on 18 February (a day before the strong Saharan dust signal in Fig. 3a, b). INP concentrations could have already been higher on 17 February; however, no in situ measurements are available on that day. It is conceivable that air masses carrying only a smaller fraction of Saharan dust had already influenced the region of Davos. Analysis of the air mass properties and history prompted us to distinguish continental aerosol and Saharan dust for 19 February 11:30 UTC. The FLEXPART dispersion model (Pisso et al.2019) indicated that air masses before that time originated from eastern Europe and Italy (not shown). Between 11:00 and 12:00 UTC (time corresponding to measurements at Davos) the flow over southern Italy was affected by intrusion of air masses from North Africa. To identify when exactly air masses from North Africa carrying Saharan dust entered the flow, we utilized the relative humidity measured at WFJ (Fig. A1a in the Appendix). A local minimum in relative humidity was identified on 19 February at 11:30 UTC. Therefore, we classify samples before that time as still continental. Nonetheless, the Jungfraujoch Research Station (approximately 140 km WSW of Weissfluhjoch) released a Saharan dust event warning only at 13:37 UTC. Brunner et al. (2021) reported that the determination of a Saharan dust event varies based on the underlying parameter used as proxy. The three samples with higher INP concentrations (orange crosses on 18 and 19 February in Fig. 3c) coincided with local maxima in relative humidity and attenuated backscatter (indicated by 1 and 2 in Fig. A1a, b, respectively). Therefore, the higher INP concentrations observed at WFJ at the end of the continental aerosol period could be a result of a stronger exchange with air masses at lower height resulting in uptake of biogenic and soil particles from fertile lands over southern Italy and eastern Europe (Conen et al.2015).

Table 2Overview of the number of samples (N) available for INP closure in dust-dominated and continental air mass cases. Additionally, the number of days over which the samples are taken is given.

 Note that due to technical maintenance of the lidar instrument only 9 of the 11 total samples during the Saharan dust event shown in Fig. 3 could be used.

Download Print Version | Download XLSX

In the following we will investigate the performance of INP parameterizations for dust-carrying and continental aerosol air masses. For the dust comparison (Sect. 3.2.1) we use 9 of the 11 in situ samples indicated by hexagons in Fig. 3c. For the two samples in the morning of 20 February, no lidar data are available for a comparison. For the comparison with continental aerosol (Sect. 3.2.2) we use the 14 samples presented with crosses in Fig. 3c and additional samples that were collected individually during clear-sky conditions before and after cloud events between 23 February and 12 March 2019, resulting in a total of 23 samples collected on 14 d (see Table 2). The grouping into Saharan dust and continental samples is consistent with Mignani et al. (2021) (see Fig. S1 in the Supplement of Mignani et al., 2021).

In the in situ and lidar observations, four sources of uncertainty are present.

  • The first source is the measurement uncertainty of the in situ observation (drop freezing technique), which depends on various parameters such as freezing temperature and number of droplets. Vali (1971) states the uncertainty to be between a factor of 2 and 4 (based on an assay of 150 droplets).

  • The second source is the retrieval uncertainty of the lidar that are passed on to the INP concentration through the applied parameterization. The retrieval uncertainty of the INP concentration from the lidar measurements is given by a factor of 3–10 (Haarig et al.2019; Mamouri and Ansmann2016).

  • The third source is air mass differences due to the horizontal distance between in situ measurements at WFJ (mountaintop site) and the lidar beam, which was 3.65 km (Fig. 1b). However, being located at mountaintop height, it is conceivable to assume no drastic difference between the two locations. Thus, the uncertainties caused by the spatial distance are negligible compared to the other sources of uncertainties.

  • Ultimately, the natural variability of INP is typically referred to as 1 order of magnitude (Kanji et al.2017).

The latter being the largest source of uncertainty, we consequently view data points of in situ measurements and lidar-based observations within 1 order of magnitude as acceptable.

Figure 4Comparison of in situ INP concentrations (nINP,obs.) measured at WFJ (mountaintop site, Fig. 1) to lidar-retrieved INP concentrations (nINP,lidar) during the Saharan dust event (19–21 February 2019) for three dust parameterizations: (a) D15 (DeMott et al.2015), (b) D15 with a calibration factor (cf=0.086) applied as proposed by Schrod et al. (2017), (c) H19 (Harrison et al.2019), and (d) U17d (Ullrich et al.2017). Notably, we compare INP concentrations at temperatures between −20C (dark orange) and −5C (light orange). The 1:1 line is shown in dashed gray. The dashed red lines show a linear regression fit through the logarithmically transformed data points, including the 95 % confidence interval in dotted cyan. Shown in every panel are the mean absolute error (MAE), the fraction of data points falling within a factor of 2 and 5 with respect to the 1:1 line (f2 and f5), and the slope of the fit (including the standard error). The gray shaded area indicates a deviation of 1 order of magnitude from the 1:1 line, which is the upper limit uncertainty of the lidar-based INP retrieval (Mamouri and Ansmann2016). For D15 (a), the data fit is also presented, applying the calibration factor from (b) for comparison (black dashed). For H19 (c), the data fit is also presented, assuming a K-feldspar fraction of 10 % (black dashed) and 1 % (black dotted), respectively. Uncertainties in INP concentration are indicated for representative samples in black.


3.2.1 Retrieval during times of Saharan dust

In situ INP observations are compared against lidar-retrieved INP concentrations using different INP parameterizations to evaluate their performance during Saharan dust presence. In Fig. 4, we present comparisons to the parameterizations of D15, H19, and U17d (Table 1) along with performance measures such as the mean absolute error (MAE), the fraction of data points falling within a factor of 2 and 5 with respect to the 1:1 line (f2 and f5), and the slope of a linear regression performed on the logarithmically transformed data. Among all investigated parameterizations, D15 (Fig. 4a) results in a slope of the linear regression fit that is closest to the desired unity. However, D15 overestimates the INP concentration, which could partially be attributed to the slight overestimation in n250 retrieval (Fig. 2a). Schrod et al. (2017) proposed a calibration factor of cf=0.086 based on a lidar comparison with Saharan dust samples collected with UAVs over Cyprus. Applying this calibration factor to the data pushes the data beyond the 1:1 line, reducing the MAE and increasing f2 and f5 (Fig. 4b). Previously, Mignani et al. (2021) compared INP observations at −15C of the same samples to predictions of D15 using in situ aerosol data and applying the calibration factor of Schrod et al. (2017). In this study, we extend the analysis to the whole temperature spectrum (−20 to −5C) and relied on aerosol retrievals from the lidar. For concentrations between 10−1 and 101StdL−1, which is the concentration range Schrod et al. (2017) used to determine the calibration factor, the fit cuts the 1:1 line (black dashed extension of the red dashed data fit). Note that Schrod et al. (2017) also proposed a slope correction; however, we did not correct the slope because all data points already fall within the uncertainty range (gray shaded area). H19 (Fig. 4c) results in a higher f2 than D15 (using cf=0.086) but with a more shallow slope leading to a lower f5. For the retrieval we assumed a K-feldspar fraction of 20 % in accordance with Kandler et al. (2011), who investigated the relative K-feldspar fractions on Cabo Verde west of central Africa. H19 scales log-linearly with the percentage of K-feldspar. During measurements at Barbados Kandler et al. (2018) found that Saharan dust featured only 1 % K-feldspar. However, using 1 % K-feldspar results in an underestimation of INP concentration (Fig. 4c). A fraction of approximately 10 % K-feldspar would match our observations best, yet the slope of the relation would be too small (0.51). U17d (Fig. 4d) overestimates the INP concentration the most strongly, leading to the highest MAE (1.93) among the investigated parameterizations.

Figure 5Same as for Fig. 4 but during continental background (see Fig. 3) for the two parameterizations: (a) D10 (DeMott et al.2010), and (b) U17s (Ullrich et al.2017).


Our comparison suggests D15 in combination with the calibration factor (cf=0.086) proposed by Schrod et al. (2017) as the preferable choice for INP concentration retrieval from dust air masses (slope closest to unity, lowest MAE). Moreover, D15 should be preferred due to its dependence on n250 instead of s, which is associated with lower lidar retrieval uncertainty (see Sect. 3.1).

3.2.2 Retrieval during times of background aerosol

Aside from a Saharan dust event, the measurement region was not affected by a long-range-transported dominating aerosol species. Thus, the samples collected outside the Saharan dust episode represent a general mix of continental aerosols. In Fig. 5, we compare the in situ observed INP concentration of these samples with lidar retrievals using the parameterizations (a) D10 for global aerosol and (b) U17s for soot, which are the best parameterizations available to mimic continental aerosol (Table 1). D10 results in a large MAE (1.84) and with no data points agreeing within a factor of 2 or 5 (f2=f5=0). Yet, the slope is close to unity (1.05). INP concentrations based on U17s result in smaller MAE (0.66) and higher f2 (0.27) and f5 (0.53), but the slope is more shallow (0.47). Purely based on the error assessment of the parameterizations in their original form, U17s would be the parameterization of choice for continental INP concentration lidar retrieval. The slope of the D10 comparison being close to unity indicates that D10 captures the composition of a continental aerosol mix. Together with the large MAE, the comparison suggests that the expected fraction of INPs among the aerosol mix is overestimated. In the next section, we propose a tuning of D10 to improve the retrieval of continental INP concentration from lidar measurements.

3.2.3 Tuning of D10 based on in situ INP and aerosol measurements

Figure 6Same as for Fig. 5a but comparing in situ INP concentrations (nINP,obs.) measured at WFJ (mountaintop site, Fig. 1) to predicted INP concentrations (nINP,APS) based on D10 now applied to in situ aerosol concentrations (from APS) during continental background (see Fig. 3). (a) D10 applied to in situ aerosol concentrations (from APS) in its standard form. (b) D10 shifted towards the 1:1 line, maximizing the fraction of data points within 1 order of magnitude around the 1:1 line (see Eq. 11).


The slope in the comparison between in situ observations and lidar-retrieved INP concentrations suggests that D10 captures the feature of a continental INP mix (Fig. 5a). The large overprediction by nearly 2 orders of magnitude on average (MAE = 1.84) indicates an overestimation in INP active aerosol fraction. D10 was developed based on nine datasets from different locations (e.g., Alaska, Pacific, Amazon Basin) and thus contains various aerosol types (DeMott et al.2010). Although a large number of non-dust samples contributed to the parameterization, dust samples, e.g., from the Pacific Dust Experiment (PACDEX, Stith et al.2009), were included in the D10 parameterization. Mineral dust is more ice active compared to (non-biogenic) continental aerosol. Therefore, it is not surprising that D10 overestimates the INP concentration when used for non-dust continental aerosol. Using the collocated in situ observations at WFJ (mountaintop site) of n250 (from APS) and INP concentration (Sect. 2.1), we could derive a tuned D10 parameterization optimized for continental aerosol. The optimization is based purely on the in situ observations in order to be independent of the lidar data (which are the test data) as we cannot fully rule out a difference in air mass. Analogously to the lidar retrieval, INP concentration is predicted by n250 (now taken from the APS) and D10 applied to it. The comparison of the obtained INP concentrations (Fig. 6a) exhibits very similar performance measures as in the lidar comparison (Fig. 5a), underlining that the difference in air mass was minuscule. Furthermore, as discussed in Sect. 2.1, the use of a heated inlet could potentially lead to the loss of volatile compounds. This could subsequently decrease the obtained n250, causing an underestimation of in situ predicted INP concentrations. We note that this effect does not apply to lidar-retrieved aerosol number concentrations and therefore lidar-retrieved INP concentrations. Therefore, the similarity of the performance measures also indicates the validity of the assumption that evaporation of volatile compounds is inconsequential. We optimize D10 with a multiplicative calibration factor (cf) as a free parameter (similarly to DeMott et al.2015; Schrod et al.2017) for maximizing the number of data points within an order of magnitude around the 1:1 line (gray shaded area in Fig. 6b). Based on the performance parameters it is apparent how the prediction has improved while the slope remained unchanged. To obtain the tuned INP concentrations (nINP,D10,tuned) from INP concentrations using D10 (nINP,D10, as defined in Eq. 9), we calculated

(11) n INP , D 10 , tuned = cf × n INP , D 10 ,

with the calibration factor cf=0.0204 and the INP concentrations given in StdL−1. The calibration factor indicates that under clean continental aerosol conditions just 2 % of the INPs predicted with D10 are present. The improvement in prediction that was found using the in situ aerosol data is also found for the lidar retrieval when applying the calibration factor (Fig. 7). The now tuned D10 parameterization outperforms U17s (compare Figs. 5b and 7b). Thus, we propose the tuned D10 parameterization (Eq. 11) to retrieve atmospheric INP concentrations from wintertime continental air masses.

Figure 7Same as for Fig. 4 but comparing in situ INP concentrations (nINP,obs.) measured at WFJ (mountaintop site, Fig. 1) to lidar-retrieved INP concentrations (nINP,lidar) using (a) the standard D10 parameterization and (b) the tuned D10 parameterization (see Eq. 11) for continental aerosol. Note that (a) is the same as Fig. 5a.


Figure 8Comparison of the remotely retrieved and in situ measured (a) ice water content (IWC) and (b) ice crystal number concentration (ICNC, nice as variable) on 8 March 2019. The IWC data were retrieved from the lowest size bin (1783 ma.s.l.). Due to measurement gaps, ICNC between 1900 and 2500 ma.s.l. was averaged (median). Left: time series of available remotely retrieved data (blue) and in situ observed data (black). Right: scatter plot of data during times the balloon was at an elevation higher than 1750 ma.s.l. (gray shading in the time series plots, maximal elevation: 1910 ma.s.l.). The gray shading in the scatter plots indicates an error of 1 order of magnitude. The values f10 and f3 indicate the fraction of data points within a factor of 10 and 3, respectively. Uncertainties are indicated by the shading for the time series and at representative data points of the scatter plots in black.


3.3 Validation of remotely retrieved cloud properties

In this section, remotely retrieved ice water contents (IWCs) and ice crystal number concentrations (ICNCs) are compared with balloon-borne in situ observations to validate the accuracy of the retrievals. On 8 March 2019, the balloon performed vertical profiles between 1650 and 1900 ma.s.l. next to the remote sensing instruments (see overview plot (Fig. 6b) in Ramelli et al.2021a). Since the lowest data bin of the radar and lidar was at 1783 ma.s.l., we considered in situ data only for times when the balloon was at an elevation of 1750 ma.s.l. and higher (gray shading in time series plots of Fig. 8, maximum elevation 1910 ma.s.l.). IWC is deduced solely from radar observations and could continuously be taken from the lowest height bin (1783 ma.s.l.). The retrieval of ICNCs needs parallel measurements of the radar and lidar. As the lidar featured data gaps, it was necessary to average (median) the ICNC data between 1900 and 2500 ma.s.l. The uncertainty of remotely retrieved cloud properties is typically given by a factor of 3 (Bühl et al.2019). Within a factor of 3, 68 % and 46 % of the data points agreed for IWC (Fig. 8a) and ICNC (Fig. 8b), respectively. For a factor of 10, the agreement was 92 % and 88 % for IWC (Fig. 8a) and ICNC (Fig. 8b), respectively. Given the inhomogeneous distribution of hydrometeors within an MPC (Korolev et al.2017) and that the balloon could have been horizontally displaced up to 200 m from the remote sensing beams, we assess the comparability to be satisfactory. In the following analysis we will consider remotely retrieved ICNC associated with an uncertainty of 1 order of magnitude.

3.4 IMF estimation in orographic MPCs

Here we propose a method to assess secondary ice production (SIP) on a single cloud basis by calculating ice multiplication factor (IMF) histograms using combined lidar and radar retrievals. First, the procedure is explained for one observed cloud case (Sect. 3.4.1). Second, the effect of sublimation below cloud base on ICNC and the temporal evolution of ICNC profiles is presented. Observed IMF histograms with height and temperature are presented and discussed with regard to the ICNC input data (Sect. 3.4.2). Third, the limitations and caveats of the method are discussed (Sect. 3.4.3). Finally, the results of this analysis are compared to observations of previous field measurements (Sect. 3.4.4).

Figure 9Exemplary illustration of data processing workflow to obtain in-cloud ice multiplication factors (IMFs) for the cloud observed on 12 March 2019. (a) Reflectivity of the cloud radar is given in shades of gray. Ice crystal number concentration (nice) is overlaid in color. The limit of used data for the IMF calculation with height is indicated in red (height span hd, time span td). Additionally, the ambient temperature at the heights of the upper and lower bounds of selected data is given in red (see Tdt and Tdb in Table D1). The INP concentration (nINP,ct) was taken from the last lidar retrieval (solid green shading) at cloud top (dashed blue line). (b) Attenuated backscatter coefficient at 1064 nm (β1064nmatt.) from the lidar in color. Cloud base estimates determined from the attenuated backscatter coefficient at 532 nm (β532nmatt.) are indicated by red dots. Regions where the Cloudnet algorithm detected liquid water are indicated in bright green. The height of cloud top is indicated by the dashed blue line. The region within the yellow rectangle was classified as below cloud base and the region within the cyan rectangle as above cloud base. Overviews of all individual cloud events can be found in Figs. B1 and C1.


3.4.1 Procedure description

Locally in a cloud, the magnitude of ice crystal enhancement can be estimated with the IMF defined as ICNC divided by INP concentration. Note that the effect of blowing snow is neglected, and sedimenting ice crystals from a cloud layer aloft are ruled out based on manual inspection of cloud radar observations. For this analysis, only ICNCs in ice saturated regions were considered, as otherwise ICNCs can decrease due to sublimation. The INP concentration is indicative of the number of primary ice crystals formed in the cloud that could subsequently evoke SIP. Taking a too low INP concentration results in an overestimation of the IMF. The (cumulative) INP concentration increases with decreasing temperature, resulting in the highest INP concentrations at cloud top (assuming that the entire cloud formed by the same ascending air mass). ICNCs at cloud top are generally created due to primary ice formation (see e.g., Crosier et al.2014). At sufficient size, the primary ice crystals will no longer be levitated by the updrafts and start to sediment to lower heights. Thus, using INP concentrations at cloud top is an upper estimate for the primary ice crystals in the cloud and avoids an overestimation of the IMF. In the next section, we determine IMFs for seven individual cloud events based on the following procedure.

  1. Based on the aforementioned considerations, we retrieve the INP concentration at cloud top height from the lidar retrieval closest to the cloud event (Fig. 9a), using the tuned D10 parameterization (Eq. 11), which is the best estimate using ground-based lidar. ICNCs were taken over a height hd (Fig. 9a), excluding heights near the radar echo boundaries to avoid retrieval errors and to avoid zones of cloud top entrainment. Additionally, data were trimmed to a time td (Fig. 9a) to account for echo-free regions at the beginning and end of the clouds. Echo-free regions in the selected data frame (e.g., at a height of approximately 3.6 kma.s.l. at around 13:40 UTC in Fig. 9a) are thought to be negligible for the consequent statistical analysis.

  2. Ice crystals sedimenting out of the cloud will start to sublimate upon entering the unsaturated environment below cloud base. Thus, in a second step, lidar variables are used to determine the cloud base height to discard data from regions where sublimation could occur. In our study, the classification was done manually by inspecting (i) the attenuated backscatter coefficient at 1064 nm (β1064nmatt.), (ii) a cloud base estimate derived from the attenuated backscatter coefficient at 532 nm (β532nmatt.), and (iii) the presence of liquid water flag by the Cloudnet algorithm (Fig. 9b). In most cases criteria (ii) and (iii) were used to determine cloud base because these estimates are deemed to be more reliable than criterion (i). However, in cases where these estimates provided frequent jumps in cloud base, for example, due to snowfall blurring the backscattered signal (see e.g., 15:00–15:30 UTC in Fig. 9b), we used (i) as a proxy for the cloud base. As the cloud base height in general changes with cloud evolution, the data were divided into time slices of approximately 1 h for which individually averaged cloud base heights were determined. The individual sizes of the time slices were set to match distinctive regime changes (e.g., the lowering of the cloud base in Fig. 9b). Consequently, each time slice was vertically divided and classified into below cloud base and above cloud base. Strong precipitation can hamper the lidar signal close to the instrument, disallowing cloud base determination above. Time slices during strong precipitation or without lidar information are therefore classified as unknown.

  3. The IMF at a given height h and time t is consequently calculated as

    (12) IMF ( h , t ) = n ice ( h , t ) n INP , ct ,

    with nice(h,t) being the ICNC at height h and time t and nINP,ct being the INP concentration at cloud top (as illustrated in Fig. 9a).

Figure 10Median ice crystal number concentrations (ICNCs) as a function of altitude for each cloud and for time slices of approximately 1 h. (a) Temporal evolution of the ICNC (nice) profiles: the faintest colors correspond to the first slice of data, and the strongest colors correspond to the last slice of data. ((b) Cloud regime classification for each time slice: ICNCs (nice) taken above cloud base (solid), below cloud base (dashed), and for an unknown cloud base (faint color). For cloud 3, the cloud base could not be determined.


3.4.2 IMF observations of seven orographic MPCs

IMFs were calculated for seven cloud events according to the procedure described above. An overview of the ICNC data and the structure of each cloud is presented in Appendix B. The time slice division and classification are provided in Appendix C. Our methodology uses a fixed INP concentration for all IMFs (Eq. 12). Thus, changes in the IMFs are driven by the evolution of the ICNCs. In a generalized case of cloud evolution, one would expect the ICNCs to increase with time due to SIP as long as conditions for SIP prevail. At cloud top primary ice crystals will nucleate, start growing, and eventually sediment. Due to the sedimentation, ICNC increases downward. If the environmental conditions permit, SIP processes start to enhance the ICNC. With the presence of more ice crystals and warmer temperatures, aggregation becomes more likely, reducing the ICNC towards cloud base and at higher temperatures. With ice crystals sedimenting below cloud base, not only aggregation but also sublimation of entire ice crystals could reduce the ICNC. In this simple picture, the recirculation of ice crystals back into the cloud due to updrafts is not considered. In reality it is not possible a priori to know where in the cloud ICNCs increase or decrease as this strongly depends on the environmental conditions and hence varies from cloud to cloud and throughout a cloud's lifetime.

Figure 10 gives an overview of the median ICNC profiles for each of the seven observed clouds in this study as time slices of approximately 1 h. A general trend of increasing ICNC with time manifests itself across different clouds (e.g., clouds 2, 3, 5, 6, and 7). For clouds 5 and 6 ICNCs increase from cloud top downwards before they start to decrease towards the ground. In contrast, the ICNC profiles of the other clouds (clouds 1, 2, 3, 4, and 7) generally decrease with decreasing height. Note that the top of the profile is below cloud top (see discussion in Sect. 3.4.1). Above cloud base, aggregation can reduce the ICNC. Below cloud base, sublimation can reduce the ICNC in addition to aggregation (Fig. 10b). Despite no cloud base determination, ICNC profiles of clouds 5 and 6 during strong precipitation increase with time near the surface (especially for cloud 6), suggesting that the relative humidity was close to 100 % with respect to ice everywhere below cloud. Note that for cloud 5 in the last hour, there could be ice crystals originating from homogeneous freezing sedimenting from higher up (see Fig. B1d). Cloud base for cloud 3 could not be determined due to a liquid layer in a lower cloud fully attenuating the lidar signal (see Fig. A1b). It is worth mentioning that the ICNCs do not differ strongly between precipitating and non-precipitating clouds.

Figure 11Probability density functions (PDFs) of obtained ice multiplication factors (IMFs, Eq. 12) for ice crystal number concentration data above cloud base (solid black), below cloud base (dashed black), and of unknown cloud base (solid gray). The solid red line indicates unity. The median IMF and the percentage of IMF larger unity for each PDF are given in the legend.


Following the definition of the IMF (Eq. 12), those below unity are linked to regions in the cloud where the INP concentration at cloud top exceeds the local ICNC. Potential regions of IMF close to and lower than unity are (i) in the beginning stage of cloud formation and (ii) near cloud top, where primary ice crystals have just formed and SIP processes are not active yet. Furthermore, ICNC and thus IMF could be reduced (iii) towards cloud base due to aggregation and (iv) below cloud base due to sublimation of ice crystals. Note that sublimation could also act as a SIP process as ice crystals could fragment during sublimation (Korolev et al.2020). The obtained IMF distributions for the three categories (above cloud base, below cloud base, and unknown cloud base) for all clouds combined are presented in Fig. 11. Above cloud base, IMFs >1 (and thus active SIP) were observed nearly 84 % of the time, with a median IMF of around 80 and an interquartile range of 2 orders of magnitude (see Table D1). As expected, IMFs below cloud base are generally found at lower values along with a wider spread than above cloud base. In contrast, the distribution of IMF where no cloud base could be determined is shifted to even higher values than the IMF distribution above cloud base (median IMF higher by approximately a factor of 3). Additionally, the shape resembles the above-cloud-base distribution. These observations suggest that ICNCs in the unknown category were potentially not strongly affected by sublimation and could be classified as above cloud base. Table D1 summarizes the IMF calculations per cloud as well as individual key numbers of the clouds (e.g., concentration of INP and aerosol at cloud top and temperatures at different heights). The individual IMF distributions per cloud and all calculated IMFs combined are presented in Fig. D1.

Figure 12Frequency histogram of all observed ice multiplication factors (IMFs) above cloud base combined with altitude. Frequencies are normalized to the number of observations (right axis) per height bin. Gray background indicates where no ice crystal number concentration data were observed. The dashed white line indicates an IMF of unity. The median IMF of all data is indicated in solid white.


In Fig. 11, it is conspicuous that the IMF distribution (above cloud base) features a heavy tail below unity signifying higher INP concentration at cloud top than local ICNCs. These IMFs can be explained when looking at the IMFs as a function of height (Fig. 12). IMFs below unity are found in two height ranges (one between 3.5 and 4 kma.s.l. and one between 4.5 and 6 kma.s.l.). The two regions belong to the early stages of clouds 5 and 6 (see Fig. 10a) where the cloud was still developing. This suggests that regions of IMFs smaller than unity could be a proxy for regions where primary ice generation dominates. After the first hour, median ICNCs of clouds 5 and 6 increased by over 2 orders of magnitude and 1 order of magnitude, respectively. Accordingly, the IMFs also increased at the corresponding heights above unity, potentially resembling a change from a primary ice-dominated regime towards the onset of SIP.

Figure 13As for Fig. 12, but with the altitude converted to temperature using data from the COSMO-1 analysis.


Despite all uncertainties related to SIP, the effectiveness of the underlying processes varies with temperature (e.g., the Hallett–Mossop process is thought to be active only between −8 and −3C). To elucidate the obtained IMF with regard to ambient temperatures, retrieval height was converted to ambient temperature using the COSMO-1 analysis. The temperatures at the top and base of the data frame were averaged over the time of available ICNC data (as listed in Table D1 and indicated in Fig. B1) and interpolated linearly with height. The resulting histogram of IMFs from all cloud events is presented in Fig. 13. Individual histograms per cloud can be found in Fig. E1. SIP was active in all observed clouds (IMFs>1) and at all observed temperatures, highlighting that SIP processes are not only restricted to warm subfreezing temperatures, where the Hallett–Mossop process prevails. High frequencies of IMF above the median IMF were observed between −10 and −5C. Hanna et al. (2008) found that for cloud tops at temperatures between −20 and −10C frequently precipitation is initiated. Strong SIP leading to high ICNC could promote precipitation initiation given a favorable environment. However, no substantial difference in IMF distributions between precipitating and non-precipitating clouds was found (highest IMFs were even found for a non-precipitating cloud, cloud 7 in Fig. E1). Such a link of the observed IMFs to the ambient environmental conditions is not made since our measurements are restricted to only seven cloud events, with the number of observations biased towards temperatures warmer than −15C.

3.4.3 Uncertainties and caveats of the proposed methodology

As stated in the introduction, using the IMF (as defined in Eq. 12) itself does not allow the inference of the underlying SIP process due to the potential (temporal and spatial) displacement of first ice crystal formation and occurrence of SIP. For the same reason, it is also not possible to quantitatively describe the uncertainty of the calculated IMFs. Three main factors add uncertainty to this method:

  • (i) uncertainties from the ICNC retrieval, which were shown to be within an order of magnitude (Fig. 8),

  • (ii) uncertainties from the INP concentration retrieval at cloud top, which are also within an order of magnitude (Fig. 7b), and

  • (iii) uncertainties from assuming a constant INP concentration with time throughout the cloud.

Based on the retrieved INP concentrations for the seven clouds (Table D1) ranging from approximately 10−3 to 10−1StdL−1, we estimate this uncertainty to also be of 1 order of magnitude. The combined influence of the uncertainties on the calculated IMF is not possible to determine due to the missing spatial (and thus logical) link between cloud top INP concentration and in-cloud ICNC following the considerations above. Our method has additional caveats.

  • All ICNC data above cloud base are assumed to not be affected by sublimation. However, it is not necessarily true that above the detected cloud base no sublimation occurred as dry air layers could have been entrained.

  • In most SIP processes, small ice splinters are created that could then start growing in favorable cloud conditions. This small ice would create a small but sharp peak at small velocities in the Doppler velocity spectrum (see e.g., Fig. 1c in Li et al.2021). In the ICNC retrieval, the entire Doppler velocity spectrum is fit to a unimodal distribution for an assumed ice crystal shape (Bühl et al.2019). Consequently, the secondary ice would only considerably affect the retrieved ICNC when the ice splinters have grown to a size when their distribution approaches and ultimately merges with the background ice distribution.

  • Excluding ICNC found in sublimation conditions could exclude contributions of the SIP process of sublimation fragmentation (Korolev et al.2020). However, the produced secondary ice crystals will not fully sublimate if they re-enter the cloud again where they would contribute to the ICNC above cloud base.

  • In our approach, a constant cloud top height and thus a constant INP concentration at cloud top were assumed for an entire cloud event. We tested that the effect of individual cloud heights did not substantially affect the obtained IMFs.

  • Calculating the ICNC from the remote sensing observations requires a sufficiently strong detectable lidar echo, limiting our method to optically thinner clouds (where the signal is not attenuated), which could induce a bias. Thus, all quantitative findings above need to be treated with caution. The statement that can be conveyed with higher certainty is the observation of SIP being active at all temperatures.

The data of each cloud were selected manually (see red data frames in Fig. B1). A robustness test of the data selection is provided in Appendix F.

3.4.4 Comparison to previous field observations

Figure 14Comparison of obtained ice multiplication factors (IMFs) in this study (gray shading indicating the 10th and 90th percentiles of observations with the median indicated by vertical lines per temperature bin) to previous field observations. The polygons envelop the observed IMFs per study. Circles indicate observations of one IMF at one temperature. Note that for the study of Crawford et al. (2012), no lower limit in IMF could be determined. Note that the retrieved IMFs from Ladino et al. (2017) tend towards 109 for temperatures towards 0 C. The dashed black line and the solid black line indicate an IMF of unity and the observed median IMF in this study, respectively.


Atmospheric ICNC and INP concentrations have been studied for more than 50 years, most often utilizing aircraft as a measurement platform. In the beginning of the millennium it was brought to the cloud physics community's attention that previous measurements could have recorded artificially enhanced ice crystal concentration due to shattering of ice crystals on the tips of the cloud probes (see e.g., Field et al.2006; Korolev et al.2011). As a response, correction algorithms were developed (e.g., Field et al.2006; Lawson2011), and the probe design was adapted (e.g., Korolev et al.2011, 2013). Thus, we compare our findings only to field observations of selected studies from the past decade in which the effect of shattering was mitigated. Details about the individual studies and the procedure of obtaining IMFs from the ice crystal and INP data are described in Appendix G.

Figure 14 comprises the range of observed IMFs in this study along with derived IMFs from previous studies which observed maritime convective systems (Stith et al.2011), continental clouds (Crawford et al.2012; Taylor et al.2016), tropical clouds (Lasher-Trapp et al.2016; Ladino et al.2017), orographic MPCs (Mignani et al.2019; Lauber et al.2021), and Arctic MPCs (Pasquier et al.2022). The IMFs derived from the previous studies not only exhibit great variability among each other but also span over a considerable value range in each individual study. The ranges of IMF from previous studies generally coincide with the observations from our study, strengthening the applicability and validity of our method. Most studies observed clouds at temperatures warmer than −15C, and some acknowledge the Hallett–Mossop process as the most likely mechanism producing secondary ice concentrations around −5C (Crawford et al.2012; Taylor et al.2016; Ladino et al.2017). Independent of temperature we observed IMFs between 101 and 103 (see Fig. 11). At temperatures warmer than −15C, previous studies consistently confirm this observation and suggest that even IMFs of 104 are frequently observed in the atmosphere. The observations of Stith et al. (2011) and Pasquier et al. (2022) suggest SIP also occurring at temperatures much lower than −15C, supporting our findings. Especially at these colder temperatures, the contribution of SIP processes remains uncertain due to the lack of available measurements. At these low temperatures, which are often encountered higher up in the troposphere, our method is a promising tool to investigate IMFs for many clouds and associate ranges of IMFs with regimes of environmental conditions.

4 Conclusions and outlook

In this study we retrieved atmospheric ice-nucleating particle (INP) concentrations in dust-dominated and continental air masses and ice multiplication factors (IMFs) in wintertime orographic mixed-phase clouds (MPCs) using active remote sensing and in situ observations obtained during the RACLETS field campaign in the Swiss Alps in February and March 2019. Using in situ aerosol and INP observations on a mountaintop in proximity to a lidar located in a nearby high valley, we investigated the remote sensing retrieval performance of aerosol properties and different INP parameterizations for dust and continental-aerosol-loaded air masses. In addition, we validated the retrievals of ice crystal number concentration (ICNC) and ice water content (IWC) derived from the combination of radar and lidar data with in situ cloud observations. We proposed a methodology to estimate IMFs from collocated lidar and radar data using a tuned INP parameterization for the continental aerosol. We conclude with the following key findings.

  1. The aerosol retrieval of the lidar could be validated. We found that the retrieval of particle number concentration n250 (number concentration of particles with physical radii 250 nm) is less prone to bias than the retrieval of surface area concentration.

  2. The dust parameterization by DeMott et al. (2015) in combination with the calibration factor derived by Schrod et al. (2017) was found to predict the INP concentrations from Saharan-dust-carrying air masses best (within a factor of 10). Although the calibration factor was derived for condensation and deposition mode INP concentrations at temperatures -20C and below, we could extend its applicability for immersion mode INP concentrations at temperatures -20C.

  3. We found that the INP parameterization by DeMott et al. (2010) (D10) in its original form overpredicts continental INP concentration in the Swiss Alps. We propose a multiplicative calibration factor (cf=0.0204) to significantly improve its performance, making it applicable to wintertime continental aerosol.

  4. Retrievals of ICNC and IWC from collocated lidar and radar measurements were found to agree with in situ cloud observations within an order of magnitude.

  5. We presented a methodology to estimate IMF from collocated lidar and radar measurements. Secondary ice production (SIP) was active nearly 84 % of the time and at all observed temperatures roughly between −30 and −5C. Overall, the median IMF found in wintertime MPCs over the Swiss Alps was around 80, with an interquartile range of 2 orders of magnitude. High frequencies of IMF above 80 were observed between −10 and −5C based on the restricted number of observed clouds.

  6. The IMFs observed with the proposed method showed substantial overlap to IMFs retrieved from previous field observations, especially at temperatures warmer than −15C, where IMFs between 101 and 104 were frequently observed.

Our study profited from an ideal setting for the closure of continental INP concentration retrieval from the lidar observations combined with in situ observations. The proposed calibration factor should be validated during future field campaigns or on large datasets from past field campaigns with suitable setup. Future development of INP parameterizations intended for lidar applications should be based on aerosol number concentration (n250) due to the reduced bias in the lidar retrieval. Different assumptions are included in the proposed methodology to determine IMF from collocated lidar and radar observations. However, the applicability of the methodology is supported by the consistency to the results from previous field observations. The tuning of the D10 parameterization was essential to reduce errors induced by the INP concentration retrieval. As pointed out by Field et al. (2017), reliable quantification of primary ice formation is crucial to constrain SIP. Without the tuning, the INP concentration would have been overestimated by nearly 2 orders of magnitude, and IMFs would have consequently been underestimated by 2 orders of magnitude. The retrieval of ICNC could be further improved by varying the assumed particle shapes. Ultimately, the efficiency and occurrence of proposed processes that cause SIP vary based on the environmental conditions (Field et al.2017; Korolev and Leisner2020). It is thus important to link the observed IMFs back to the environmental conditions the cloud formed in. The application of our methodology to large datasets of ground-based or space-borne aerosol and cloud observations could help to assess IMFs in different cloud types and constrain the environmental conditions for different SIP regimes. Machine learning tools such as random forests or gradient-boosted decision trees could be valuable tools in constraining environmental regimes impacting SIP. This could not only accelerate our understanding of ice multiplication in the atmosphere but also be of immense help to motivate and plan laboratory experiments that investigate individual ice multiplication processes.

Appendix A: Categorization of continental background and Saharan dust samples

Shortly before the Saharan dust event on 19 February 2019, the Davos region was influenced by air masses already featuring high INP concentrations and a higher aerosol content suggested by an increased signal in the attenuated backscatter (see events 1 and 2 in Fig. A1b). The fact that these air masses do not necessarily belong to the following plumes of Saharan dust is further supported by the peaks observed in relative humidity before the phase of Saharan dust (Fig. A1a).

Figure A1(a) Relative humidity observed during the comparison period (14–22 February 2019) on 19 February 2019 at WFJ (mountaintop site). The two peaks in relative humidity during the transition time are indicated with 1 and 2. (b) Observed lidar attenuated backscatter coefficient β1064nmatt. at the height of WFJ (2687 ma.s.l.). The two characteristic time periods for continental background aerosol and Saharan dust are indicated. Two perturbations in β1064nmatt. during the transition time are indicated with 1 and 2, respectively.


Appendix B: Overview of individual clouds

The exemplary procedure for obtaining IMFs was described for cloud 7 observed on 12 March 2019 in Fig. 9. Figure B1 shows the overview for all observed clouds. For cloud 4 (Fig. B1c) the last lidar retrieval was more than 1 d before the ICNC data were available. After the ICNC retrieval time the next lidar profile was also more than a day later. However, the profile used (as indicated in Fig. B1c) and the next lidar profile only differed by a factor of 2; thus a strong deviation during the available ICNC data times seems unlikely. For cloud 5 the last lidar profile prior to the event was only available 2 d before. Thus, the closest profile in the morning of the next day was used. Small parts of cloud 5 reach temperatures below −38C where potentially homogeneous freezing could occur and additional ice crystals could sediment into higher temperature levels.

Figure B1As in Fig. 9a for each panel except for the parameter labeling of hd, td, and nINP,ct. The green shading indicates the period during which the INP concentration at cloud top (dashed blue line) was retrieved. Calculated and retrieved parameters for each cloud are presented in Table D1. The blue dotted line in (d) indicates the height above which homogeneous freezing is possible (T<-38C).


Appendix C: Cloud regime classification

The exemplary procedure for obtaining IMFs was described for cloud 7 observed on 12 March 2019 in Fig. 9. Figure C1 shows the classification for all observed clouds. Cloud 3 (Fig. C1b) was entirely classified as unknown due to the liquid water layer at around 2500 m.

Figure C1Time slicing and classification of the ice crystal number concentration data indicated in Fig. B1 (red rectangle encompasses rectangles presented in this plot). White background represents regions where the lidar either was fully attenuated or underwent automatic maintenance. The attenuated backscatter coefficient at 1064 nm (β1064nmatt.) is presented in the background. A (liquid) cloud base estimate derived from the attenuated backscatter coefficient at 532 nm (β532nmatt.) is indicated by red circles. The automated liquid water classification by the Cloudnet algorithm is presented by lime tiles. Empty rectangles in cyan, yellow, and magenta indicate regions classified as above cloud base, below cloud base, and unknown cloud base, respectively. Cloud top derived from radar observations (Fig. B1) is indicated by the dashed blue line. The dotted light blue line in (d) indicates the height above which homogeneous freezing is possible (T<-38C).


Appendix D: IMF summary and additional cloud information

Table D1 provides a summary of obtained IMF for all ICNC data (red rectangles in Fig. B1) and using only ICNC data above cloud base (cyan rectangles in Fig. C1) to demonstrate the effect of the selection. Generally, the selection does not change the median IMF (cloud 1) or increases it (clouds 2, 4, and 7). In the cases of clouds 5 and 6 the median IMF decreased due to the selection. This is explained by the exclusion of high ICNCs that were present during strong precipitation when the lidar was attenuated close to the ground and no cloud base could be determined. Figure D1 shows IMF histograms per cloud and for all clouds combined using all ICNC data.

Table D1Overview of the secondary ice production (SIP) analysis per investigated cloud and environmental information for each cloud. For each cloud the total number of ice crystal number concentration (ICNC) observations (Nobs.), and the median ICNC (ñice), the times when SIP was active (SIP active, i.e., percentage of IMFs>1), and the median IMF (IMF̃) and its interquartile range (IQRIMF, in orders of magnitude) are presented for all ICNC data available and ICNC data above a detectable cloud base. Further, the concentration of particles with radii 250 nm and INPs at cloud top before the event (n250,ct, nINP,ct), the cloud height (hct), the height (hd) and time (td) over which ICNC was obtained, and the temperatures at cloud top (Tct) and top (Tdt) and bottom (Tdb) of available ICNC data are given. In the last column, the total number of observations (Nobs.) and averages (weighted by number of observations) of the variables are given combined for all clouds.

Download Print Version | Download XLSX

Figure D1Probability density functions (PDFs) of obtained ice multiplication factors (IMFs, Eq. 12) for each of the seven observed clouds and all IMFs combined (thick solid red) using all ice crystal number concentration data (not only above cloud base). The dashed black line indicates unity. The median IMF and the percentage of IMF larger than unity for the combined histogram are given in the legend.


Appendix E: IMF temperature histogram per cloud

IMF histograms with temperature for each cloud for ICNC data above cloud base are presented in Fig. E1. Note that for cloud 3 no cloud base could be determined (see Fig. C1b).

Figure E1Frequency histogram of observed ice multiplication factors (IMFs) above cloud base with temperature per cloud. Frequencies are normalized to the number of observations (right axis) per temperature bin. Radar height was converted to temperature using data from the COSMO-1 analysis. Gray background indicates where no data were observed. The dashed white line indicates unity. The median IMF of all data is indicated in solid white.


Appendix F: Sensitivity of obtained IMF based on data selection

To assess the robustness of our methodology we performed a sensitivity study by changing the upper and lower heights of the selected ICNC by 200 m and the starting and ending times by 20 % of the entire time span towards the center (see hd and td in Fig. 9a). Note that this test was done on the entire ICNC data available, i.e., the data frames indicated in red in Fig. B1. As seen in Table F1 the obtained SIP parameters did not drastically change upon changes in the data selection. The fraction of time when SIP was active did not change by more than five percentage points (highest upon changing the starting time by 20 %, ts+20 %), and the IMF-related parameters did not change by more than a factor of 1.5 (highest upon changing the starting time by 20 %, ts+20 %).

Table F1Benchmarks of sensitivity studies of obtained secondary ice production (SIP) and ice multiplication factor (IMF) based on data selection for both estimates. The data selection limits in height and time were changed by ±200m and by ±20 % of the time span, respectively.

Download Print Version | Download XLSX

Appendix G: Extraction of IMFs from previous publications

Airborne ice crystal measurements could have overestimated the ICNC due to shattering of ice crystals on the cloud probes (see e.g., Field et al.2006; Korolev et al.2011). The retrieval of ICNC from radar observations has also become more reliable with the further technical development of the instruments. Hence, we considered publications from the last decade to compare our obtained IMF to them (Table G1). However, most studies provide information of ambient temperature, INP concentration, and ICNC separately. Additionally, INP concentrations are typically given at the same location or temperature as the observed ICNC – in contrast to our approach. In the following we briefly describe how IMFs were retrieved from each study. From Stith et al. (2011) Fig. 8, maximum and minimum IMF were obtained for the two legs 05:22 and 05:42 UTC. The ambient temperature ranges during the legs were taken from Stith et al. (2011) Fig. 5. Crawford et al. (2012) estimate an upper limit of 0.01 L−1 in INP concentration based on source considerations. From Crawford et al. (2012) Fig. 3, maximal ICNC with temperature were taken. A minimal ICNC could not be distinguished from zero such that IMFs smaller than 1 are possible, but not further determinable. Lasher-Trapp et al. (2016) Fig. 5 presents the excess of ICNC over INP concentration combined as a function of temperature. We considered ICNC data that were not prone to possible side contamination as the authors write. Taylor et al. (2016) Fig. 13 provides ICNC with estimates on INP concentration with temperature. Note that we consider Run 11.1 as a point measurement. Ladino et al. (2017) Fig. 4 encompasses an extensive overview of their ICNC measurements and predicted INP concentration with temperature. IMFs were calculated using the envelope of all ICNC observations and the INP concentration at corresponding altitude. Mignani et al. (2019) provide IMF between 4 and 25 (Fig. 3) at temperatures between −17 and −12C (Fig. 2). Lauber et al. (2021) Fig. 6 provides a median ICNC concentration close to the melting layer. INP concentrations below the detection at this temperature are estimated to be smaller than the lower detection limit of the drop freezing instrument. For measurements close to the melting layer, a temperature of −0.5C was used for plotting. IMFs from Pasquier et al. (2022) were provided by personal communication. Li et al. (2021) do not provide correlated data of ICNC, INP concentration, and temperature. They state that “ice number concentrations tend to be 1–3 orders of magnitude higher than expected INP concentrations for clouds at temperatures of  C or warmer”.

Stith et al. (2011)Crawford et al. (2012)Lasher-Trapp et al. (2016)Taylor et al. (2016)Ladino et al. (2017)Mignani et al. (2019)Lauber et al. (2021)Li et al. (2021)Pasquier et al. (2022)

Table G1Overview of the studies considered in the ice multiplication factor comparison and the corresponding cloud type sampled in the study.

Download Print Version | Download XLSX

Code and data availability

Evaluation scripts used in this study are available upon request. An overview of all performed measurements during RACLETS can be found on the campaign's website The in situ aerosol data can be accessed via (last access: 6 March 2022, and Rösch2020) and (last access: 6 March 2022, et al.2020). The primary lidar output (e.g., backscatter) can be accessed via (last access: 4 June 2022). The Cloudnet products can be accessed via (last access: 4 June 2022). The in situ cloud microphysical properties, lidar-retrieved aerosol properties, lidar cloud base estimates, and remotely retrieved cloud microphysical properties can be accessed via (last access: 4 June 2022, Wieder et al.2022a).

Author contributions

JW and JH conceived the study. NI performed the initial analysis. JW performed the final analysis and prepared the figures for the manuscript. JW and CM performed the aerosol measurements during the campaign. MH provided the processed lidar data. MH performed manual optimization of the data quality to lower the detection limit of the lidar in the clean environment. JB and PS operated the radar and provided the processed radar data. RE performed the lidar measurements. FR provided the in situ cloud data. JW, NI, CM, MH, ZAK, UL, and JH interpreted the data. JW wrote the manuscript with contributions from all co-authors. All authors reviewed the manuscript. JH supervised the project.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors express their gratitude to the whole RACLETS campaign team for their support and many fruitful discussions. In particular, we thank Michael Lehning (WSL/SLF, EPFL) and his team for their support in realizing the RACLETS campaign. A big thanks to Mario Schär and Lucie Roth for their help with the aerosol measurements. We thank Paul Fopp for providing his land at Wolfgangpass and Martin Genter for logistical support at Weissfluhjoch. Our deepest appreciation to Michael Rösch and Marco Vecellio for technical support. We express our deepest gratitude to Nora Els (University of Innsbruck) for providing us with a second impinger. We thank MeteoSchweiz for the meteorological observations and the COSMO-1 analysis data. We thank Annett Skupin and Hannes Griesche for setting up the lidar at Davos and performing the lidar measurements. In addition, special thanks are due to Albert Ansmann for helping in reprocessing the lidar data. The authors thank Franz Conen, Robert David, Julie T. Pasquier, Maxim Samarin, and Colin Tully for discussions and suggestions improving the manuscript. Annika Lauber is thanked for providing the photograph of the setup at Wolfgangpass. We acknowledge the use of the scientific color maps provided by Crameri et al. (2020). Furthermore, the authors want to thank the two anonymous reviewers, for their comments and suggestions, which improved the paper.

Financial support

This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant nos. 200021_169620 and 200021_175824).

Review statement

This paper was edited by Anne Perring and reviewed by two anonymous referees.


Ansmann, A., Bösenberg, J., Chiakovsky, A., Comerón, A., Eckhardt, S., Eixmann, R., Freudenthaler, V., Ginoux, P., Komguem, L., Linné, H., Ángel Lopez Marquez, M., Matthias, V., Mattis, I., Mitev, V., Müller, D., Music, S., Nickovic, S., Pelon, J., Sauvage, L., Sobolewsky, P., Srivastava, M. K., Stohl, A., Torres, O., Vaughan, G., Wandinger, U., and Wiegner, M.: Long-range transport of Saharan dust to northern Europe: The 11-16 October 2001 outbreak observed with EARLINET, J. Geophys. Res.-Atmos., 108, 4783,, 2003. a

Ansmann, A., Tesche, M., Althausen, D., Müller, D., Seifert, P., Freudenthaler, V., Heese, B., Wiegner, M., Pisani, G., Knippertz, P., and Dubovik, O.: Influence of Saharan dust on cloud glaciation in southern Morocco during the Saharan Mineral Dust Experiment, J. Geophys. Res.-Atmos., 113, D04210,, 2008. a

Ansmann, A., Tesche, M., Seifert, P., Groß, S., Freudenthaler, V., Apituley, A., Wilson, K. M., Serikov, I., Linné, H., Heinold, B., Hiebsch, A., Schnell, F., Schmidt, J., Mattis, I., Wandinger, U., and Wiegner, M.: Ash and fine-mode particle mass profiles from EARLINET-AERONET observations over central Europe after the eruptions of the Eyjafjallajökull volcano in 2010, J. Geophys. Res.-Atmos., 116, D00U02,, 2011. a

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,, 2019a. a, b, c, d

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,, 2019b. a

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,, 2021. a

Auer, A. H., Veal, D. L., and Marwitz, J. D.: Observations of Ice Crystal and Ice Nuclei Concentrations in Stable Cap Clouds, J. Atmos. Sci., 26, 1342–1343,<1342:OOICAI>2.0.CO;2, 1969. a, b

Baars, H., Seifert, P., Engelmann, R., and Wandinger, U.: Target categorization of aerosol and clouds by continuous multiwavelength-polarization lidar measurements, Atmos. Meas. Tech., 10, 3175–3201,, 2017. a

Baumgardner, D., Brenguier, J., Bucholtz, A., Coe, H., DeMott, P., Garrett, T., Gayet, J., Hermann, M., Heymsfield, A., Korolev, A., Krämer, M., Petzold, A., Strapp, W., Pilewskie, P., Taylor, J., Twohy, C., Wendisch, M., Bachalo, W., and Chuang, P.: Airborne instruments to measure atmospheric aerosol particles, clouds and radiation: A cook's tour of mature and emerging technology, Atmos. Res., 102, 10–29,, 2011. a

Beck, A.: Observing the Microstructure of Orographic Clouds with HoloGondel, PhD thesis, ETH Zurich,, 2017. a, b

Beck, A., Henneberger, J., Fugal, J. P., David, R. O., Lacher, L., and Lohmann, U.: Impact of surface and near-surface processes on ice crystal concentrations measured at mountain-top research stations, Atmos. Chem. Phys., 18, 8909–8927,, 2018. a

Bedoya-Velásquez, A. E., Navas-Guzmán, F., Granados-Muñoz, M. J., Titos, G., Román, R., Casquero-Vera, J. A., Ortiz-Amezcua, P., Benavent-Oltra, J. A., de Arruda Moreira, G., Montilla-Rosero, E., Hoyos, C. D., Artiñano, B., Coz, E., Olmo-Reyes, F. J., Alados-Arboledas, L., and Guerrero-Rascado, J. L.: Hygroscopic growth study in the framework of EARLINET during the SLOPE I campaign: synergy of remote sensing and in situ instrumentation, Atmos. Chem. Phys., 18, 7001–7017,, 2018. a

Bergeron, T.: On the physics of clouds and precipitation, Proc. 5th Assembly UGGI, Lisbon, Portugal, 1935, 156–180, 1935. a

Brunner, C., Brem, B. T., Collaud Coen, M., Conen, F., Hervo, M., Henne, S., Steinbacher, M., Gysel-Beer, M., and Kanji, Z. A.: The contribution of Saharan dust to the ice-nucleating particle concentrations at the High Altitude Station Jungfraujoch (3580 m a.s.l.), Switzerland, Atmos. Chem. Phys., 21, 18029–18053,, 2021. a

Burton, S. P., Ferrare, R. A., Hostetler, C. A., Hair, J. W., Rogers, R. R., Obland, M. D., Butler, C. F., Cook, A. L., Harper, D. B., and Froyd, K. D.: Aerosol classification using airborne High Spectral Resolution Lidar measurements – methodology and examples, Atmos. Meas. Tech., 5, 73–98,, 2012. a

Bühl, J., Ansmann, A., Seifert, P., Baars, H., and Engelmann, R.: Toward a quantitative characterization of heterogeneous ice formation with lidar/radar: Comparison of CALIPSO/CloudSat with ground-based observations, Geophys. Res. Lett., 40, 4404–4408,, 2013. a

Bühl, J., Seifert, P., Radenz, M., Baars, H., and Ansmann, A.: Ice crystal number concentration from lidar, cloud radar and radar wind profiler measurements, Atmos. Meas. Tech., 12, 6601–6617,, 2019. a, b, c, d, e

Conen, F., Rodríguez, S., Hülin, C., Henne, S., Herrmann, E., Bukowiecki, N., and Alewell, C.: Atmospheric ice nuclei at the high-altitude observatory Jungfraujoch, Switzerland, Tellus B, 67, 25014,, 2015. a

Cotton, R. J., Field, P. R., Ulanowski, Z., Kaye, P. H., Hirst, E., Greenaway, R. S., Crawford, I., Crosier, J., and Dorsey, J.: The effective density of small ice particles obtained from in situ aircraft observations of mid-latitude cirrus, Q. J. Roy. Meteor. Soc., 139, 1923–1934,, 2013. a

Crameri, F., Shephard, G. E., and Heron, P. J.: The misuse of colour in science communication, Nat. Commun., 11, 5444,, 2020. a

Crawford, I., Bower, K. N., Choularton, T. W., Dearden, C., Crosier, J., Westbrook, C., Capes, G., Coe, H., Connolly, P. J., Dorsey, J. R., Gallagher, M. W., Williams, P., Trembath, J., Cui, Z., and Blyth, A.: Ice formation and development in aged, wintertime cumulus over the UK: observations and modelling, Atmos. Chem. Phys., 12, 4963–4985,, 2012. a, b, c, d, e, f, g

Crosier, J., Bower, K. N., Choularton, T. W., Westbrook, C. D., Connolly, P. J., Cui, Z. Q., Crawford, I. P., Capes, G. L., Coe, H., Dorsey, J. R., Williams, P. I., Illingworth, A. J., Gallagher, M. W., and Blyth, A. M.: Observations of ice multiplication in a weakly convective cell embedded in supercooled mid-level stratus, Atmos. Chem. Phys., 11, 257–273,, 2011. a

Crosier, J., Choularton, T. W., Westbrook, C. D., Blyth, A. M., Bower, K. N., Connolly, P. J., Dearden, C., Gallagher, M. W., Cui, Z., and Nicol, J. C.: Microphysical properties of cold frontal rainbands, Q. J. Roy. Meteor. Soc., 140, 1257–1268,, 2014. a

David, R. O., Cascajo-Castresana, M., Brennan, K. P., Rösch, M., Els, N., Werz, J., Weichlinger, V., Boynton, L. S., Bogler, S., Borduas-Dedekind, N., Marcolli, C., and Kanji, Z. A.: Development of the DRoplet Ice Nuclei Counter Zurich (DRINCZ): validation and application to field-collected snow samples, Atmos. Meas. Tech., 12, 6865–6888,, 2019. a

Delanoë, J., Protat, A., Jourdan, O., Pelon, J., Papazzoni, M., Dupuy, R., Gayet, J. F., and Jouan, C.: Comparison of airborne in situ, airborne radar-lidar, and spaceborne radar-lidar retrievals of polar ice cloud properties sampled during the POLARCAT campaign, J. Atmos. Ocean. Tech., 30, 57–73,, 2013. a

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,, 2010. a, b, c, d, e, f, g, h, i, j

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,, 2015. a, b, c, d, e, f, g, h

Düsing, S., Wehner, B., Seifert, P., Ansmann, A., Baars, H., Ditas, F., Henning, S., Ma, N., Poulain, L., Siebert, H., Wiedensohler, A., and Macke, A.: Helicopter-borne observations of the continental background aerosol in combination with remote sensing and ground-based measurements, Atmos. Chem. Phys., 18, 1263–1290,, 2018. a

Düsing, S., Ansmann, A., Baars, H., Corbin, J. C., Denjean, C., Gysel-Beer, M., Müller, T., Poulain, L., Siebert, H., Spindler, G., Tuch, T., Wehner, B., and Wiedensohler, A.: Measurement report: Comparison of airborne, in situ measured, lidar-based, and modeled aerosol optical properties in the central European background – identifying sources of deviations, Atmos. Chem. Phys., 21, 16745–16773,, 2021. a, b

Eidhammer, T., DeMott, P. J., Prenni, A. J., Petters, M. D., Twohy, C. H., Rogers, D. C., Stith, J., Heymsfield, A., Wang, Z., Pratt, K. A., Prather, K. A., Murphy, S. M., Seinfeld, J. H., Subramanian, R., and Kreidenweis, S. M.: Ice initiation by aerosol particles: Measured and predicted ice nuclei concentrations versus measured ice crystal concentrations in an orographic wave cloud, J. Atmos. Sci., 67, 2417–2436,, 2010. a

Engelmann, R., Kanitz, T., Baars, H., Heese, B., Althausen, D., Skupin, A., Wandinger, U., Komppula, M., Stachlewska, I. S., Amiridis, V., Marinou, E., Mattis, I., Linné, H., and Ansmann, A.: The automated multiwavelength Raman polarization and water-vapor lidar PollyXT: the neXT generation, Atmos. Meas. Tech., 9, 1767–1784,, 2016. a

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,, 2021. a, b

Envidat: RACLETS field campaign, (last access: 12 June 2022), 2019. a

Ferrare, R. A., Melfi, S. H., Whiteman, D. N., Evans, K. D., Poellot, M., and Kaufman, Y. J.: Raman lidar measurements of aerosol extinction and backscattering 2. Derivation of aerosol real refractive index, single-scattering albedo, and humidification factor using Raman lidar and aircraft size distribution measurements, J. Geophys. Res.-Atmos., 103, 19673–19689,, 1998. a

Field, P. R. and Heymsfield, A. J.: Importance of snow to global precipitation, Geophys. Res. Lett., 42, 9512–9520,, 2015. a

Field, P. R., Heymsfield, A. J., and Bansemer, A.: Shattering and Particle Interarrival Times Measured by Optical Array Probes in Ice Clouds, J. Atmos. Ocean. Tech., 23, 1357–1371,, 2006. a, b, c

Field, P. R., Lawson, R. P., Brown, P. R. A., Lloyd, G., Westbrook, C., Moisseev, D., Miltenberger, A., Nenes, A., Blyth, A., Choularton, T., Connolly, P., Buehl, J., Crosier, J., Cui, Z., Dearden, C., DeMott, P., Flossmann, A., Heymsfield, A., Huang, Y., Kalesse, H., Kanji, Z. A., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Chapter 7. Secondary Ice Production - current state of the science and recommendations for the future, Meteorol. Monogr., 58, 1–20,, 2017. a, b

Findeisen, W.: Kolloid-meteorologische Vorgänge bei Neiderschlags-bildung, Meteor. Z, 55, 121–133, 1938. a

Gayet, J. F., Treffeisen, R., Helbig, A., Bareiss, J., Matsuki, A., Herber, A., and Schwarzenboeck, A.: On the onset of the ice phase in boundary layer Arctic clouds, J. Geophys. Res.-Atmos., 114, D19201,, 2009. a

Georgakaki, P., Bougiatioti, A., Wieder, J., Mignani, C., Ramelli, F., Kanji, Z. A., Henneberger, J., Hervo, M., Berne, A., Lohmann, U., and Nenes, A.: On the drivers of droplet variability in alpine mixed-phase clouds, Atmos. Chem. Phys., 21, 10993–11012,, 2021. a, b

Griesche, H. J., Seifert, P., Ansmann, A., Baars, H., Barrientos Velasco, C., Bühl, J., Engelmann, R., Radenz, M., Zhenping, Y., and Macke, A.: Application of the shipborne remote sensing supersite OCEANET for profiling of Arctic aerosols and clouds during Polarstern cruise PS106, Atmos. Meas. Tech., 13, 5335–5358,, 2020. a

Groß, S., Esselborn, M., Weinzierl, B., Wirth, M., Fix, A., and Petzold, A.: Aerosol classification by airborne high spectral resolution lidar observations, Atmos. Chem. Phys., 13, 2487–2505,, 2013. a

Görsdorf, U., Lehmann, V., Bauer-Pfundstein, M., Peters, G., Vavriv, D., Vinogradov, V., and Volkov, V.: A 35-GHz Polarimetric Doppler Radar for Long-Term Observations of Cloud Parameters-Description of System and Data Processing, J. Atmos. Ocean. Tech., 32, 675–690,, 2015. a

Haarig, M., Ansmann, A., Althausen, D., Klepel, A., Groß, S., Freudenthaler, V., Toledano, C., Mamouri, R.-E., Farrell, D. A., Prescod, D. A., Marinou, E., Burton, S. P., Gasteiger, J., Engelmann, R., and Baars, H.: Triple-wavelength depolarization-ratio profiling of Saharan dust over Barbados during SALTRACE in 2013 and 2014, Atmos. Chem. Phys., 17, 10767–10794,, 2017. a

Haarig, M., Walser, A., Ansmann, A., Dollner, M., Althausen, D., Sauer, D., Farrell, D., and Weinzierl, B.: Profiles of cloud condensation nuclei, dust mass concentration, and ice-nucleating-particle-relevant aerosol properties in the Saharan Air Layer over Barbados from polarization lidar and airborne in situ measurements, Atmos. Chem. Phys., 19, 13773–13788,, 2019. a, b, c, d, e, f, g

Hallett, J. and Mossop, S. C.: Production of secondary ice particles during the riming process, Nature, 249, 26–28,, 1974. a

Hanna, J. W., Schultz, D. M., and Irving, A. R.: Cloud-top temperatures for precipitating winter clouds, J. Appl. Meteorol. Clim., 47, 351–359,, 2008. a

Harrison, A. D., Lever, K., Sanchez-Marroquin, A., Holden, M. A., Whale, T. F., Tarn, M. D., McQuaid, J. B., and Murray, B. J.: The ice-nucleating ability of quartz immersed in water and its atmospheric importance compared to K-feldspar, Atmos. Chem. Phys., 19, 11343–11361,, 2019. a, b, c, d

Heymsfield, A. and Willis, P.: Cloud conditions favoring secondary ice particle production in tropical maritime convection, J. Atmos. Sci., 71, 4500–4526,, 2014. a

Heymsfield, A. J., Schmitt, C., Chen, C.-C.-J., Bansemer, A., Gettelman, A., Field, P. R., and Liu, C.: Contributions of the Liquid and Ice Phases to Global Surface Precipitation: Observations and Global Climate Modeling, J. Atmos. Sci., 77, 2629–2648,, 2020. a

Hobbs, P. V. and Rangno, A. L.: Ice Particle Concentrations in Clouds, J. Atmos. Sci., 42, 2523–2549,<2523:IPCIC>2.0.CO;2, 1985. a

Hobbs, P. V. and Rangno, A. L.: Rapid Development of High Ice Particle Concentrations in Small Polar Maritime Cumuliform Clouds, J. Atmos. Sci., 47, 2710–2722,<2710:RDOHIP>2.0.CO;2, 1990. a

Hobbs, P. V. and Rangno, A. L.: Microstructures of low and middle-level clouds over the Beaufort Sea, Q. J. Roy. Meteor. Soc., 124, 2035–2071,, 1998. a

Hogan, R. J., Mittermaier, M. P., and Illingworth, A. J.: The Retrieval of Ice Water Content from Radar Reflectivity Factor and Temperature and Its Use in Evaluating a Mesoscale Model, J. Appl. Meteorol. Clim., 45, 301–317,, 2006. a, b

Holben, B., Eck, T., Slutsker, I., Tanré, D., Buis, J., Setzer, A., Vermote, E., Reagan, J., Kaufman, Y., Nakajima, T., Lavenu, F., Jankowiak, I., and Smirnov, A.: AERONET – Federated Instrument Network and Data Archive for Aerosol Characterization, Remote Sens. Environ., 66, 1–16,, 1998. a

Huang, S., Hu, W., Chen, J., Wu, Z., Zhang, D., and Fu, P.: Overview of biological ice nucleating particles in the atmosphere, Environ. Int., 146, 106197,, 2021. a, b, c

IPCC, 2021: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, in press.,, 2021. a

Kandler, K., Schütz, L., Jäckel, S., Lieke, K., Emmel, C., Müller-Ebert, D., Ebert, M., Scheuvens, D., Schladitz, A., Šegvić, B., Wiedensohler, A., and Weinbruch, S.: Ground-based off-line aerosol measurements at Praia, Cape Verde, during the Saharan Mineral Dust Experiment: Microphysical properties and mineralogy, Tellus B, 63, 459–474,, 2011. a

Kandler, K., Schneiders, K., Ebert, M., Hartmann, M., Weinbruch, S., Prass, M., and Pöhlker, C.: Composition and mixing state of atmospheric aerosols determined by electron microscopy: method development and application to aged Saharan dust deposition in the Caribbean boundary layer, Atmos. Chem. Phys., 18, 13429–13455,, 2018. a

Kanji, Z. A., Ladino, L. A., Wex, H., Boose, Y., Burkert-Kohn, M., Cziczo, D. J., and Krämer, M.: Overview of Ice Nucleating Particles, Meteorol. Monogr., 58, 11–33,, 2017. a, b, c, d, e

Koenig, L. R.: The Glaciating Behavior of Small Cumulonimbus Clouds, J. Atmos. Sci., 20, 29–47,<0029:TGBOSC>2.0.CO;2, 1963. a

Korolev, A. and Leisner, T.: Review of experimental studies of secondary ice production, Atmos. Chem. Phys., 20, 11767–11797,, 2020. a, b, c, d

Korolev, A., McFarquhar, G., Field, P. R., Franklin, C., Lawson, P., Wang, Z., Williams, E., Abel, S. J., Axisa, D., Borrmann, S., Crosier, J., Fugal, J., Krämer, M., Lohmann, U., Schlenczek, O., Schnaiter, M., and Wendisch, M.: Mixed-Phase Clouds: Progress and Challenges, Meteorol. Monogr., 58, 1–50,, 2017. a, b

Korolev, A., Heckman, I., Wolde, M., Ackerman, A. S., Fridlind, A. M., Ladino, L. A., Lawson, R. P., Milbrandt, J., and Williams, E.: A new look at the environmental conditions favorable to secondary ice production, Atmos. Chem. Phys., 20, 1391–1429,, 2020. a, b

Korolev, A. V., Emery, E. F., Strapp, J. W., Cober, S. G., Isaac, G. A., Wasey, M., and Marcotte, D.: Small Ice Particles in Tropospheric Clouds: Fact or Artifact?, Airborne Icing Instrumentation Evaluation Experiment, B. Am. Meteorol. Soc., 92, 967–973,, 2011. a, b, c

Korolev, A. V., Emery, E. F., Strapp, J. W., Cober, S. G., and Isaac, G. A.: Quantification of the effects of shattering on airborne ice particle measurements, J. Atmos. Ocean. Tech., 30, 2527–2553,, 2013. a

Ladino, L. A., Korolev, A., Heckman, I., Wolde, M., Fridlind, A. M., and Ackerman, A. S.: On the role of ice-nucleating aerosol in the formation of ice particles in tropical mesoscale convective systems, Geophys. Res. Lett., 44, 1574–1582,, 2017. a, b, c, d, e, f

Lasher-Trapp, S., Leon, D. C., DeMott, P. J., Villanueva-Birriel, C. M., Johnson, A. V., Moser, D. H., Tully, C. S., and Wu, W.: A multisensor investigation of rime splintering in tropical maritime cumuli, J. Atmos. Sci., 73, 2547–2564,, 2016. a, b, c, d

Lauber, A., Henneberger, J., Mignani, C., Ramelli, F., Pasquier, J. T., Wieder, J., Hervo, M., and Lohmann, U.: Continuous secondary-ice production initiated by updrafts through the melting layer in mountainous regions, Atmos. Chem. Phys., 21, 3855–3870,, 2021. a, b, c, d, e, f

Lawson, R. P.: Effects of ice particles shattering on the 2D-S probe, Atmos. Meas. Tech., 4, 1361–1381,, 2011. a

Lawson, R. P., Woods, S., and Morrison, H.: The microphysics of ice and precipitation development in tropical cumulus clouds, J. Atmos. Sci., 72, 2429–2445,, 2015. a

Li, H., Möhler, O., Petäjä, T., and Moisseev, D.: Two-year statistics of columnar-ice production in stratiform clouds over Hyytiälä, Finland: environmental conditions and the relevance to secondary ice production, Atmos. Chem. Phys., 21, 14671–14686,, 2021. a, b, c

Lohmann, U., Henneberger, J., Henneberg, O., Fugal, J. P., Bühl, J., and Kanji, Z. A.: Persistence of orographic mixed-phase clouds, Geophys. Res. Lett., 43, 10512–10519,, 2016. a

Luke, E. P., Yang, F., Kollias, P., Vogelmann, A. M., and Maahn, M.: New insights into ice multiplication using remote-sensing observations of slightly supercooled mixed-phase clouds in the Arctic, P. Natl. Acad. Sci. USA, 118, e2021387118,, 2021. a

Mahrt, F., Wieder, J., Dietlicher, R., Smith, H. R., Stopford, C., and Kanji, Z. A.: A high-speed particle phase discriminator (PPD-HS) for the classification of airborne particles, as tested in a continuous flow diffusion chamber, Atmos. Meas. Tech., 12, 3183–3208,, 2019. a

Mamouri, R.-E. and Ansmann, A.: Potential of polarization lidar to provide profiles of CCN- and INP-relevant aerosol parameters, Atmos. Chem. Phys., 16, 5905–5931,, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m, n

Marinou, E., Tesche, M., Nenes, A., Ansmann, A., Schrod, J., Mamali, D., Tsekeri, A., Pikridas, M., Baars, H., Engelmann, R., Voudouri, K.-A., Solomos, S., Sciare, J., Groß, S., Ewald, F., and Amiridis, V.: Retrieval of ice-nucleating particle concentrations from lidar observations and comparison with UAV in situ measurements, Atmos. Chem. Phys., 19, 11315–11342,, 2019. a, b, c, d, e, f, g, h, i, j

Matus, A. V. and L’Ecuyer, T. S.: The role of cloud phase in Earth’s radiation budget, J. Geophys. Res., 122, 2559–2578,, 2017. a

McCluskey, C. S., Ovadnevaite, J., Rinaldi, M., Atkinson, J., Belosi, F., Ceburnis, D., Marullo, S., Hill, T. C., Lohmann, U., Kanji, Z. A., O'Dowd, C., Kreidenweis, S. M., and DeMott, P. J.: Marine and Terrestrial Organic Ice-Nucleating Particles in Pristine Marine to Continentally Influenced Northeast Atlantic Air Masses, J. Geophys. Res.-Atmos, 123, 6196–6212,, 2018. a

McCoy, D. T., Tan, I., Hartmann, D. L., Zelinka, M. D., and Storelvmo, T.: On the relationships among cloud cover, mixed-phase partitioning, and planetary albedo in GCMs, J. Adv. Model. Earth Sys., 8, 650–668,, 2016. a

Mignani, C., Creamean, J. M., Zimmermann, L., Alewell, C., and Conen, F.: New type of evidence for secondary ice formation at around −15C in mixed-phase clouds, Atmos. Chem. Phys., 19, 877–886,, 2019. a, b, c, d

Mignani, C., Wieder, J., Sprenger, M. A., Kanji, Z. A., Henneberger, J., Alewell, C., and Conen, F.: Towards parameterising atmospheric concentrations of ice-nucleating particles active at moderate supercooling, Atmos. Chem. Phys., 21, 657–664,, 2021. a, b, c, d, e, f

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. a

Müller, D., Lee, K.-H., Gasteiger, J., Tesche, M., Weinzierl, B., Kandler, K., Müller, T., Toledano, C., Otto, S., Althausen, D., and Ansmann, A.: Comparison of optical and microphysical properties of pure Saharan mineral dust observed with AERONET Sun photometer, Raman lidar, and in situ instruments during SAMUM 2006, J. Geophys. Res.-Atmos., 117, D07211,, 2012. a

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,, 2012. a

Murray, B. J., Carslaw, K. S., and Field, P. R.: Opinion: Cloud-phase climate feedback and the importance of ice-nucleating particles, Atmos. Chem. Phys., 21, 665–679,, 2021. a

Niemand, M., Möhler, O., Vogel, B., Vogel, H., Hoose, C., Connolly, P., Klein, H., Bingemer, H., Demott, P., Skrotzki, J., and Leisner, T.: A particle-surface-area-based parameterization of immersion freezing on desert dust particles, J. Atmos. Sci., 69, 3077–3092,, 2012. a

Pasquier, J. T., Henneberger, J., Ramelli, F., Lauber, A., David, R. O., Wieder, J., Carlsen, T., Gierens, R., Maturilli, M., and Lohmann, U.: Conditions favorable for secondary ice production in Arctic mixed-phase clouds, Atmos. Chem. Phys. Discuss. [preprint],, in review, 2022. a, b, c, d, e

Pisso, I., Sollum, E., Grythe, H., Kristiansen, N. I., Cassiani, M., Eckhardt, S., Arnold, D., Morton, D., Thompson, R. L., Groot Zwaaftink, C. D., Evangeliou, N., Sodemann, H., Haimberger, L., Henne, S., Brunner, D., Burkhart, J. F., Fouilloux, A., Brioude, J., Philipp, A., Seibert, P., and Stohl, A.: The Lagrangian particle dispersion model FLEXPART version 10.4, Geosci. Model Dev., 12, 4955–4997,, 2019. a

Proske, U., Bessenbacher, V., Dedekind, Z., Lohmann, U., and Neubauer, D.: How frequent is natural cloud seeding from ice cloud layers (< −35C) over Switzerland?, Atmos. Chem. Phys., 21, 5195–5216,, 2021. a

Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, 2nd edn., Springer, Dordrecht, 954 pp.,, 2010. a

Ramelli, F., Beck, A., Henneberger, J., and Lohmann, U.: Using a holographic imager on a tethered balloon system for microphysical observations of boundary layer clouds, Atmos. Meas. Tech., 13, 925–939,, 2020. a

Ramelli, F., Henneberger, J., David, R. O., Bühl, J., Radenz, M., Seifert, P., Wieder, J., Lauber, A., Pasquier, J. T., Engelmann, R., Mignani, C., Hervo, M., and Lohmann, U.: Microphysical investigation of the seeder and feeder region of an Alpine mixed-phase cloud, Atmos. Chem. Phys., 21, 6681–6706,, 2021a. a, b, c, d, e, f, g, h

Ramelli, F., Henneberger, J., David, R. O., Lauber, A., Pasquier, J. T., Wieder, J., Bühl, J., Seifert, P., Engelmann, R., Hervo, M., and Lohmann, U.: Influence of low-level blocking and turbulence on the microphysics of a mixed-phase cloud in an inner-Alpine valley, Atmos. Chem. Phys., 21, 5151–5172,, 2021b. a, b, c

Sakai, T., Nagai, T., Nakazato, M., Mano, Y., and Matsumura, T.: Ice clouds and Asian dust studied with lidar measurements of particle extinction-to-backscatter ratio, particle depolarization, and water-vapor mixing ratio over Tsukuba, Appl. Opt., 42, 7103–7116,, 2003. a, b

Sawamura, P., Moore, R. H., Burton, S. P., Chemyakin, E., Müller, D., Kolgotin, A., Ferrare, R. A., Hostetler, C. A., Ziemba, L. D., Beyersdorf, A. J., and Anderson, B. E.: HSRL-2 aerosol optical measurements and microphysical retrievals vs. airborne in situ measurements during DISCOVER-AQ 2013: an intercomparison study, Atmos. Chem. Phys., 17, 7229–7243,, 2017. a

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,, 2017. a, b, c, d, e, f, g, h, i, j, k, l, m

Seifert, P., Ansmann, A., Mattis, I., Wandinger, U., Tesche, M., Engelmann, R., Müller, D., Pérez, C., and Haustein, K.: Saharan dust and heterogeneous ice formation: Eleven years of cloud observations at a central European EARLINET site, J. Geophys. Res.-Atmos., 115, D20201,, 2010. a, b, c, d

Seifert, P., Ansmann, A., Groß, S., Freudenthaler, V., Heinold, B., Hiebsch, A., Mattis, I., Schmidt, J., Schnell, F., Tesche, M., Wandinger, U., and Wiegner, M.: Ice formation in ash-influenced clouds after the eruption of the Eyjafjallajkull volcano in April 2010, J. Geophys. Res.-Atmos., 116, D00U04,, 2011. a

Shimizu, A., Sugimoto, N., Matsui, I., Arao, K., Uno, I., Murayama, T., Kagawa, N., Aoki, K., Uchiyama, A., and Yamazaki, A.: Continuous observations of Asian dust and other aerosols by polarization lidars in China and Japan during ACE-Asia, J. Geophys. Res.-Atmos, 109, D19S17,, d19S17, 2004. a

Sotiropoulou, G., Ickes, L., Nenes, A., and Ekman, A. M. L.: Ice multiplication from ice–ice collisions in the high Arctic: sensitivity to ice habit, rimed fraction, ice type and uncertainties in the numerical description of the process, Atmos. Chem. Phys., 21, 9741–9760,, 2021. a

Stith, J. L., Ramanathan, V., Cooper, W. A., Roberts, G. C., DeMott, P. J., Carmichael, G., Hatch, C. D., Adhikary, B., Twohy, C. H., Rogers, D. C., Baumgardner, D., Prenni, A. J., Campos, T., Gao, R., Anderson, J., and Feng, Y.: An overview of aircraft observations from the Pacific Dust Experiment campaign, J. Geophys. Res.-Atmos., 114, 6229–6243,, 2009. a

Stith, J. L., Twohy, C. H., DeMott, P. J., Baumgardner, D., Campos, T., Gao, R., and Anderson, J.: Observations of ice nuclei and heterogeneous freezing in a Western Pacific extratropical storm, Atmos. Chem. Phys., 11, 6229–6243,, 2011. a, b, c, d, e, f

Stopelli, E., Conen, F., Zimmermann, L., Alewell, C., and Morris, C. E.: Freezing nucleation apparatus puts new slant on study of biological ice nucleators in precipitation, Atmos. Meas. Tech., 7, 129–134,, 2014. a

Sun, B. Z. and Shine, K. P.: Studies of the radiative properties of ice and mixed-phase clouds, Q. J. Roy. Meteor. Soc., 120, 111–137,, 1994. a

Takahashi, T., Nagao, Y., and Kushiyama, Y.: Possible High Ice Particle Production during Graupel-Graupel Collisions, J. Atmos. Sci., 52, 4523–4527,<4523:PHIPPD>2.0.CO;2, 1995. a

Taylor, J. W., Choularton, T. W., Blyth, A. M., Liu, Z., Bower, K. N., Crosier, J., Gallagher, M. W., Williams, P. I., Dorsey, J. R., Flynn, M. J., Bennett, L. J., Huang, Y., French, J., Korolev, A., and Brown, P. R. A.: Observations of cloud microphysics and ice formation during COPE, Atmos. Chem. Phys., 16, 799–826,, 2016. a, b, c, d

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

Thomas, D. and Charvet, A.: An Introduction to Aerosols, Aerosol Fi, Elsevier,, 2017. a

Touloupas, G., Lauber, A., Henneberger, J., Beck, A., and Lucchi, A.: A convolutional neural network for classifying cloud particles recorded by imaging probes, Atmos. Meas. Tech., 13, 2219–2239,, 2020. a

Ullrich, R., Hoose, C., Möhler, O., Niemand, M., Wagner, R., Höhler, K., Hiranuma, N., Saathoff, H., and Leisner, T.: A new ice nucleation active site parameterization for desert dust and soot, J. Atmos. Sci., 74, 699–717,, 2017. a, b, c, d, e, f, g, h, i

Vali, G.: Quantitative Evaluation of Experimental Results an the Heterogeneous Freezing Nucleation of Supercooled Liquids, J. Atmos. Sci., 28, 402–409,<0402:QEOERA>2.0.CO;2, 1971. a, b, c

Vardiman, L.: The Generation of Secondary Ice Particles in Clouds by Crystal–Crystal Collision, J. Atmos. Sci., 35, 2168–2180,<2168:TGOSIP>2.0.CO;2, 1978. a

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

Walter, B., Huwald, H., Gehring, J., Bühler, Y., and Lehning, M.: Radar measurements of blowing snow off a mountain ridge, The Cryosphere, 14, 1779–1794,, 2020. a, b

Wandinger, U.: Introduction to Lidar, in: Lidar. Springer Series in Optical Sciences, vol 102, edited by: Weitkamp, C., Springer, New York, NY, 1–18,, 2005. a

Wandinger, U. and Ansmann, A.: Experimental determination of the lidar overlap profile with Raman lidar, Appl. Opt., 41, 511–514,, 2002. a

Wandinger, U., Müller, D., Böckmann, C., Althausen, D., Matthias, V., Bösenberg, J., Weiß, V., Fiebig, M., Wendisch, M., Stohl, A., and Ansmann, A.: Optical and microphysical characterization of biomass- burning and industrial-pollution aerosols from- multiwavelength lidar and aircraft measurements, J. Geophys. Res.-Atmos., 107, 7–20,, 2002.  a

Wegener, A.: Thermodynamik der Atmosphäre, J. A. Barth Verlag, Leipzig, Germany, 331 pp., 1911. a, b

Wieder, J. and Rösch, C.: Aerosol Data Davos Wolfgang, EnviDat [data set],, 2020. a

Wieder, J., Mignani, C., and Rösch, C.: Aerosol Data Weissfluhjoch, EnviDat [data set],, 2020. a

Wieder, J., Ihn, N., Mignani, C., Haarig, M., Bühl, J., Seifert, P., Engelmann, R., Ramelli, F., Kanji, Z. A., Lohmann, U., and Henneberger, J.: Data for the publication “Retrieving ice-nucleating particle concentration and ice multiplication factors using active remote sensing validated by in situ observations”,, 2022a. a

Wieder, J., Mignani, C., Schär, M., Roth, L., Sprenger, M., Henneberger, J., Lohmann, U., Brunner, C., and Kanji, Z. A.: Unveiling atmospheric transport and mixing mechanisms of ice-nucleating particles over the Alps, Atmos. Chem. Phys., 22, 3111–3130,, 2022b. a, b, c, d, e, f, g, h

Zhang, D., Wang, Z., Heymsfield, A., Fan, J., and Luo, T.: Ice concentration retrieval in stratiform mixed-phase clouds using cloud radar reflectivity measurements and 1D ice growth model simulations, J. Atmos. Sci., 71, 3613–3635,, 2014. a

Short summary
Ice formation and its evolution in mixed-phase clouds are still uncertain. We evaluate the lidar retrieval of ice-nucleating particle concentration in dust-dominated and continental air masses over the Swiss Alps with in situ observations. A calibration factor to improve the retrieval from continental air masses is proposed. Ice multiplication factors are obtained with a new method utilizing remote sensing. Our results indicate that secondary ice production occurs at temperatures down to −30 °C.
Final-revised paper