Retrieval of ice-nucleating particle concentrations from lidar observations and comparison with UAV in situ measurements

Aerosols that are efficient ice-nucleating particles (INPs) are crucial for the formation of cloud ice via heterogeneous nucleation in the atmosphere. The distribution of INPs on a large spatial scale and as a function of height determines their impact on clouds and climate. However, in situ measurements of INPs provide sparse coverage over space and time. A promising approach to address this gap is to retrieve INP concentration profiles by combining particle concentration profiles derived by lidar measurements with INP efficiency parameterizations for different freezing mechanisms (immersion freezing, deposition nucleation). Here, we assess the feasibility of this new method for both ground-based and spaceborne lidar measurements, using in situ observations collected with unmanned aerial vehicles (UAVs) and subsequently analyzed with the FRIDGE (FRankfurt Ice nucleation Deposition freezinG Experiment) INP counter from an experimental campaign at Cyprus in April 2016. Analyzing five case studies we calculated the cloud-relevant particle number concentrations using lidar measurements (n250,dry with an uncertainty of 20 % to 40 % and Sdry with an uncertainty of 30 % to 50 %), and we assessed the suitability of the different INP parameterizations with respect to the temperature range and the type of particles considered. Specifically, our analysis suggests that our calculations using the parameterization of Ullrich et al. (2017) (applicable for the temperature range−50 to−33 C) agree within 1 order of magnitude with the in situ observations of nINP; thus, the parameterization of Ullrich et al. (2017) can efficiently address the deposition nucleation pathway in dust-dominated environments. Additionally, our calculations using the combination of the parameterizations of DeMott et al. (2015, 2010) (applicable for the temperature range −35 to −9 C) agree within 2 orders of magnitude with the in situ observations of INP concentrations (nINP) and can thus efficiently address the immersion/condensation pathway of dust and nondust particles. The same conclusion is derived from the compilation of the parameterizations of DeMott et al. (2015) for dust and Ullrich et al. (2017) for soot. Published by Copernicus Publications on behalf of the European Geosciences Union. 11316 E. Marinou et al.: Evaluation of lidar-derived INP concentrations Furthermore, we applied this methodology to estimate the INP concentration profiles before and after a cloud formation, indicating the seeding role of the particles and their subsequent impact on cloud formation and characteristics. More synergistic datasets are expected to become available in the future from EARLINET (European Aerosol Research Lidar Network) and in the frame of the European ACTRISRI (Aerosols, Clouds, and Trace gases Research Infrastructure). Our analysis shows that the developed techniques, when applied on CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) spaceborne lidar observations, are in agreement with the in situ measurements. This study gives us confidence for the production of global 3D products of cloud-relevant particle number concentrations (n250,dry, Sdry and nINP) using the CALIPSO 13-year dataset. This could provide valuable insight into the global heightresolved distribution of INP concentrations related to mineral dust, as well as possibly other aerosol types.

Abstract. Aerosols that are efficient ice-nucleating particles (INPs) are crucial for the formation of cloud ice via heterogeneous nucleation in the atmosphere. The distribution of INPs on a large spatial scale and as a function of height determines their impact on clouds and climate. However, in situ measurements of INPs provide sparse coverage over space and time. A promising approach to address this gap is to retrieve INP concentration profiles by combining particle concentration profiles derived by lidar measurements with INP efficiency parameterizations for different freezing mechanisms (immersion freezing, deposition nucleation). Here, we assess the feasibility of this new method for both ground-based and spaceborne lidar measurements, using in situ observations collected with unmanned aerial vehicles (UAVs) and subsequently analyzed with the FRIDGE (FRankfurt Ice nucleation Deposition freezinG Experiment) INP counter from an experimental campaign at Cyprus in April 2016. Analyzing five case studies we calculated the cloud-relevant particle number concentrations using lidar measurements (n 250,dry with an uncertainty of 20 % to 40 % and S dry with an uncertainty of 30 % to 50 %), and we assessed the suitability of the different INP parameterizations with respect to the temperature range and the type of particles considered. Specifically, our analysis suggests that our calculations using the parameterization of Ullrich et al. (2017) (applicable for the temperature range −50 to −33 • C) agree within 1 order of magnitude with the in situ observations of n INP ; thus, the parameterization of Ullrich et al. (2017) can efficiently address the deposition nucleation pathway in dust-dominated environments. Additionally, our calculations using the combination of the parameterizations of DeMott et al. (2015DeMott et al. ( , 2010 (applicable for the temperature range −35 to −9 • C) agree within 2 orders of magnitude with the in situ observations of INP concentrations (n INP ) and can thus efficiently address the immersion/condensation pathway of dust and nondust particles. The same conclusion is derived from the compilation of the parameterizations of DeMott et al. (2015) for dust and Ullrich et al. (2017) for soot.

Introduction
The interaction of aerosol particles with clouds and the related climatic effects have been in the focus of atmospheric research for several decades. Aerosols can act as cloud condensation nuclei (CCN) in liquid water clouds and as icenucleating particles (INPs) in mixed-phase and ice clouds. Changes in their concentration affect cloud extent, lifetime, particle size and radiative properties (Lohmann and Feichter, 2005;Tao et al., 2012;Altaratz et al., 2014;Rosenfeld et al., 2014). As important these interactions are, they are the source of the highest uncertainty in assessing the anthropogenic climate change (IPCC Fifth Assessment Report, Seinfeld et al., 2016).
All clouds producing ice require, for temperatures above ∼ −35 • C, the presence of INPs. Compared to CCN, INPs are rare (about one particle in a million acts as an INP; Nenes et al., 2014) and become increasingly sparse with increasing temperature (Pruppacher and Klett, 1997;Kanji et al., 2017).
Aerosol species which are identified in the past as potentially important INPs are mineral dust, biological species (pollen, bacteria, fungal spores and plankton), carbonaceous combustion products, soot, volcanic ash and sea spray (Murray et al., 2012;DeMott et al., 2015b). From these aerosol types, mineral dust and soot are efficient INPs at temperatures below −15 to −20 • C (dust) and −40 • C (soot), and they have been studied extensively for their INP properties in field experiments and laboratory studies (Twohy et al., 2009(Twohy et al., , 2017Kamphus et al., 2010;Hoose and Möhler, 2012;Murray et al., 2012;Sullivan et al., 2016;Ullrich et al., 2017). Biological particles are one of the most active INP species; however, their abundance is likely low on a global scale, particularly when compared to other aerosol types such as mineral dust (Morris et al., 2014). It has been suggested that soil and clay particles may act as carriers of biological nanoscale INPs (e.g., proteins), which could potentially contribute to a global/local source of INPs (Schnell and Vali, 1976;O'Sullivan et al., 2014O'Sullivan et al., , 2015O'Sullivan et al., , 2016. Finally, marine aerosols (with possible influence of a biological microlayer close to the surface) are also important INPs in areas where the influence of mineral dust is less pronounced (e.g., Southern Ocean; Wilson et al., 2015;Vergara-Temprado et al., 2017).
There is a variety of pathways for heterogeneous ice nucleation: contact freezing, immersion freezing, condensation freezing and deposition nucleation (Vali et al., 2015). Individual ice nucleation pathways dominate at characteristic temperatures and supersaturation ranges. Observational studies have shown that immersion freezing dominates at temperatures higher than −30 • C, while deposition nucleation dominates below −35 • C Westbrook et al., 2011;de Boer et al., 2011). The factors that regulate the efficiency of heterogeneous ice nucleation are qualitatively understood, but no general theory of heterogeneous ice nucleation exists yet. It has been shown that, in regions not influenced by sea salt aerosol, INP concentrations are strongly correlated with the number of aerosol particles with dry radius greater than 250 nm (n 250,dry ) which form the reservoir of favorable INPs (DeMott et al., 2010. However, we have limited knowledge on how the ice nuclei activity of these particles together with their spatial and vertical distributions depend on cloud nucleation conditions (i.e., temperature (T ) and supersaturation over water (ss w ) and ice (ss i )). Furthermore, field measurements of INP concentrations are very localized in space and time, whilst there are large regions without any data at all (Murray et al., 2012). The lack of data inhibits our quantitative understanding of aerosol-cloud interactions and requires new strategies for obtaining datasets (Seinfeld et al., 2016;Bühl et al., 2016).
Active remote sensing with aerosol lidar and cloud radar provides valuable data for studying aerosol-cloud interaction since it enables observations with high vertical and temporal resolution over long time periods (Ansmann et al., 2005;Illingworth et al., 2007;Seifert et al., 2010;de Boer et al., 2011;Kanitz et al., 2011;Bühl et al., 2016). Lidar measurements can provide profiles of n 250,dry (the number of aerosol particles with dry radius greater than 250 nm) and S dry (the aerosol particle dry surface area concentration) related to mineral dust, continental pollution and marine aerosol, as described in Ansmann (2015, 2016). Their methodology uses lidar-derived optical parameters (i.e., the particle backscatter coefficient, lidar ratio and particle depolarization ratio) to separate the contribution of mineral dust in the lidar profiles ) and subsequently applies sun-photometer-based parameterizations to transform the optical property profiles into profiles of aerosol mass, number and surface area concentration (Ansmann et al., 2012;Ansmann, 2015, 2016). The latter can E. Marinou et al.: Evaluation of lidar-derived INP concentrations 11317 then be used as input to INP parameterizations that have been obtained from laboratory and field measurements (e.g., De-Mott et al., 2010Niemand et al., 2012;Steinke et al., 2015;Ullrich et al., 2017) to derive profiles of INP concentrations (n INP ).
The INP retrieval calculated from the lidar measurements provides a promising insight into atmospheric INP concentrations. To date, there has been no other evaluation of the lidar-derived profiles of n 250,dry , S dry and n INP by means of independent in situ observations apart from one dust case in Schrod et al. (2017). The study presented here compares n 250,dry and n INP as inferred from spaceborne and ground-based lidar observations to findings from airborne in situ measurements using data from the joint experiment INUIT-BACCHUS-ACTRIS (Ice Nuclei Research Unit -Impact of Biogenic versus Anthropogenic emissions on Clouds and Climate: towards a Holistic UnderStanding -Aerosols, Clouds, and Trace gases Research Infrastructure) held in April 2016 in Cyprus (Schrod et al., 2017;Mamali et al., 2018). The paper starts with a review of the different INP parameterizations for mineral dust, soot and continental aerosols in Sect. 2. Section 3 describes the instruments used in this study and the methodology to retrieve INP concentrations from lidar measurements. The results of the intercomparison between the lidar-derived and unmanned aerial vehicle (UAV)-measured n 250,dry and n INP profiles are presented and discussed in Sect. 4 before the paper closes with a summary in Sect. 5.

INP parameterizations
A variety of parameterizations has been proposed to obtain n INP from aerosol concentration measurements. In particular, a global aerosol type-independent n INP parameterization was introduced by DeMott et al.  Ullrich et al. (2017), and soot-specific n INP parameterizations were proposed by Murray et al. (2012) and Ullrich et al. (2017). The aforementioned parameterizations address immersion freezing at or above water saturation and deposition nucleation for ice saturation ratios ranging from unity up to the homogeneous freezing threshold and water saturation. Table 1 provides an overview of the temperature ranges and the freezing mechanisms for which these parameterizations are applicable.
Regarding immersion freezing, the aerosols that are activated to droplets can contribute to ice formation. In turn, the ability of a particle to be activated as a cloud droplet mainly depends on the cloud supersaturation, its diameter, the water adsorption characteristics and the composition of soluble coatings (Levin et al., 2005;Kumar et al., 2011a, b;Garimella et al., 2014;Bègue et al., 2015). Kumar et al. (2011b) showed that all dry-generated dust samples with ra-dius > 50 nm are activated to CCN at a water supersaturation (ss w ) of 0.5 %, while the activation radius increases to > 250 nm when water supersaturation decreases to ss w ≈ 0.1 %. This is the minimum level of ss w required to activate INP for immersion freezing.
For immersion freezing of dust particles, the parameterization of Ullrich et al. (2017) (U17-imm) ( Table 1; Eq. 1) is based on heterogeneous ice nucleation experiments at the cloud chamber AIDA (Aerosol Interaction and Dynamics in the Atmosphere) of the Karlsruhe Institute of Technology. The desert dust ground samples used in this study originated from different desert dust locations around the world (Saharan, Taklamakan Desert, Canary Islands, Israel). The parameterization quantifies the desert dust ice nucleation efficiency as a function of ice-nucleation-active surface-site density n s (T ) and dust dry surface area concentration S d,dry . If the CCN activated fraction is less than 50 %, Eq. (1) for U17imm needs to be scaled to be representative for the CCN activated S dry (Ullrich et al., 2017). In this work, we apply the U17-imm parameterization taking into consideration the total S dry .
Additionally, the parameterization of DeMott et al. (2015) (D15) ( Table 1; Eq. 2) addresses the immersion and condensation freezing activity of natural mineral dust particles based on laboratory studies using the continuous flow diffusion chamber (CFDC) of the Colorado State University (CSU) and field data from atmospheric measurements in Saharan dust layers. D15 quantifies n INP as a function of temperature and the total number concentration of dust particles with dry radii larger than 250 nm (n 250,d,dry ). We note here that the ambient values of measured n INP (p, T ) need to be transferred in standard (std) pressure and temperature conditions (n 250,d,dry (p 0 , T 0 , T )) before the use of (Eq. 2).
For the deposition nucleation of dust particles, the parameterizations of Steinke et al. (2015) and Ullrich et al. (2017) (S15 and U17-dep, respectively) quantify the ice nucleation efficiency as a function of S d,dry and n s (T , S ice ) with S ice the ice saturation ratio. Both were based on AIDA laboratory studies, but they used different dust samples. U17-dep (Table 1; Eq. 3) was based on ground desert dust samples from the Sahara, Taklamakan Desert, Canary Islands and Israel, while S15 (Table 1; Eq. 4) was based on dust samples from Arizona, which were treated (washed, milled, treated with acid) and are much more ice active than natural desert dust particles on average. Although S15 parameterization was based on treated dust samples which usually show an enhanced freezing efficiency, it is used in the NMME-DREAM model (Nonhydrostatic Mesoscale Model on E grid, Janjic et al., 2001;Dust REgional Atmospheric Model, Nickovic et al., 2001;Pérez et al., 2006) for INP concentration estimations (Nickovic et al., 2016). For this reason, it is included in this work.
For the ice activation of soot particles, Ullrich et al. (2017) introduced two parameterizations, one for immersion freezing (Table 1; Eq. 5) and a second one for deposition nucle- Table 1. Overview of INP parameterizations used in this study together with the freezing mode and the temperature range for which they have been developed. The parameterizations of D15 and U17-imm have been extrapolated to the temperature range from −36 to −1 • C. In the equations, n 250,dry is in particles per cubic centimeters (cm −3 ), n INP in particles per liter (L −1 ), T (z) in Kelvin (K) and P in hectopascals (hPa). p 0 and T 0 stand for standard pressure and temperature. The particle concentrations used are derived assuming an extinction coefficient of 50 Mm −1 for each of the different aerosol types (dust, continental, soot). The shaded areas take into account a range of the extinction coefficient from 10 Mm −1 (lower limit) to 200 Mm −1 (upper limit). The error bars mark the error of the respective parameterizations from error propagation using the uncertainties provided in Table 2. Negative error bars that exceed the scale are not shown. In the deposition mode (a), the bold lines correspond to ice supersaturation of 1.15 and the dashed lines to ice supersaturation of 1.05, 1.1, 1.2, 1.3 and 1.4. The black and orange dots indicate the maximum temperatures for which the parameterizations have been developed. In the immersion mode (b), the parameterizations are extrapolated over the immersion freezing temperature range (dashed lines). ation (Table 1; Eq. 6). Both were based on experiments at the AIDA chamber with soot samples generated from four different devices and quantify the soot ice nucleation efficiency as a function of S dry and n s (T ) (for immersion) and n s (T , S ice ) (for deposition).
Finally, the global type-independent n INP parameterization of DeMott et al. (2010) (Table 1; Eq. 7) was based on field data collected during nine field campaigns (in Colorado, eastern Canada, the Amazon, Alaska, and the Pacific basin) and analyzed with the CFDC instrument of the CSU. As the majority of the samples used for D10 were nondesert continental aerosols, this INP parameterization has been considered to be suitable for addressing the immersion and condensation freezing activity of mixtures of anthropogenic haze, biomass burning smoke, biological particles, soil and road dust . From here on these mixtures are addressed as continental aerosols.
The n 250,dry and S dry used in all the aforementioned parameterizations are calculated from the lidar extinction profiles as described in Sect. 3.2 and shown in Figs. A1 and A2 in the Appendix. Figure 1 provides an indication of the relative differences of the observed n INP in nature for immersion (right) and deposition (left) modes and in relation to the different aerosol compositions by showing a summary of the different n INP parameterizations. Specifically, the plot shows the fraction of the ice-activated particles (f i = n INP /n 50,dry ) for desert dust (dark blue, orange, red, light blue), continental (green) and soot (black). The particle concentrations used here are derived assuming an extinction coefficient of 50 Mm −1 for each of the different aerosol types (dust, continental, soot). The shaded areas take into account a range of the extinction coefficient from 10 Mm −1 (lower limit) to 200 Mm −1 (upper limit). The error bars mark the cumulative error in f i that results from the uncertainty in the lidar observations and their conversion to mass concentration as well as from the errors in the respective parameterizations. An overview of the typical values and the uncertainties used for the error estimation in this study is provided in Table 2 Figure 1a shows that, for deposition mode, the dust iceactivated fractions from S15 are several orders of magnitude higher than those of U17-imm (e.g., 4 orders of magnitude at −40 • C and ss i = 1.15 %). Furthermore, the deposition iceactivation fractions of dust and soot (from U17-dep) differ significantly, with soot being more active than dust for T <-38 • C (up to 2 orders of magnitude) and dust being more active than soot for T > −38 • C (up to 4 orders of magnitude). Figure 1a shows that, for immersion mode, the dust iceactivated fractions obtained from D15 are 1 order of magnitude lower than those calculated with U17-imm.

E. Marinou et al.: Evaluation of lidar-derived INP concentrations
Laboratory ice nucleation measurements and corresponding instrument intercomparisons have shown that at a single temperature differences between 2 and 4 orders of magnitude are observed as a result of the natural variability of the INP active fraction (DeMott et al., 2010 or the use of different INP counters (Burkert-Kohn et al., 2017). Hereon, we consider D15 and U17-imm as the lower and upper bounds of the immersed n INP estimations for dust INP populations. Figure 1b illustrates the dust activation increase of up to 6 orders of magnitude within the mixed-phase temperature regime (−15 to −35 • C). For a 5 • C decrease, n d,INP increases by about 1 order of magnitude. Moreover, we see that at T < −18 • C the immersion freezing desert dust ice activation (D15) is higher than the continental one (D10), while this changes at T > −18 • C. On the contrary, soot (U17imm) always has a lower f i than dust (from either D15 or U17-imm). The ice-activated fractions of continental (D10) and soot (U17-imm) aerosols have a relative difference that is always less than 60 % at T < −18 • C. At higher temperatures they diverge with continental f i to exceed the soot one by 1 order of magnitude at T > −11 • C.
Additionally, Fig. 1 provides an indication of the error induced at the lidar-estimated n INP due to errors in the selected values of T and ss i . The right panel shows that, for immersion mode, a 5 • C error in the assumed T can introduce an error of 1 order of magnitude in the dust-related n INP estimations (U17-imm and D15) and 1/2 order of magnitude in the nondust-related estimations of D10. The same error (1/2 order of magnitude) is induced in the U17-imm(soot) (for T < −18 • C). For deposition mode, a 5 • C error in the assumed T can introduce an error of 1/2 order of magnitude in the dust-related n INP estimations (U17-dep(dust) and S15). For the U17-dep(soot) estimates, and at T > −45 • C, the error in the assumed T has a significant impact in the n INP product (e.g., 1 order of magnitude between T = −45 and −40 • C). On the contrary, at T < −45 • C, the error in the assumed T has less impact in the final n INP product (between 100 % and 200 % for 5 • C T error).
Regarding the deposition nucleation, a large variability of the onset saturation ratio is observed in laboratory experiments of different studies, with ss i varying for example at −40 • C between 1 and 1.5 . In Fig. 1a, we see the effect of the ss i on the estimated n INP . In S15, n INP values increase by 1 order of magnitude for a 0.1 increase in the ss i . In U17-dep(dust), a 3-orders-ofmagnitude n INP range is observed at −30 • C for ss i between 1.05 and 1.4. The range is wider at lower temperatures (4 orders at −50 • C). In U17-dep(soot) a 4-orders-of-magnitude n INP range is observed at T < −40 • C for ss i between 1.05 and 1.3. This variability provides an indication of the error induced in the lidar-estimated n INP product due to the error in the selected ss i . In the n INP profiles presented in Fig. 11, ss i = 1.15 is assumed (bold line here).

Instruments and methodology
The INUIT-BACCHUS-ACTRIS campaign in April 2016 was organized within the framework of the projects Ice Nuclei Research Unit (INUIT; https://www.ice-nuclei.de/ the-inuit-project/, last access: 8 August 2019); Impact of Biogenic versus Anthropogenic emissions on Clouds and Climate: towards a Holistic UnderStanding (BACCHUS; http://www.bacchus-env.eu/, last access: 8 August 2019); and Aerosols, Clouds, and Trace gases Research Infrastructure (ACTRIS; https://www.actris.eu/, last access: 8 August 2019) and focused on aerosols, clouds and ice nucleation within dust-laden air over the Eastern Mediterranean. Although dust was the main component observed, other aerosol types were present as well such as soot and continental aerosols.
The atmospheric measurements conducted during the campaign included remote sensing with aerosol lidar and sun photometers as well as in situ particle sampling with two UAVs. The UAV provided observations of the INP abundance in the lower troposphere and they were operated from the airfield of the Cyprus Institute at Orounda (35 • 05 42 N, 33 • 04 53 E; 327 m a.s.l.; about 21 km west of Nicosia) (Schrod et al., 2017). An Aerosol Robotic Network (AERONET, Holben et al., 1998)

Lidar measurements
The EARLINET PollyXT-NOA lidar measurements at 532 nm are used in this study for the derivation of particle optical properties and mass concentration profiles. Quicklooks of all PollyXT measurements can be found on the web page of PollyNet (Raman and polarization lidar network, http://polly.tropos.de, last access: 8 August 2019). Pol-lyXT operates using a Nd:YAG laser that emits light at 355, 532 and 1064 nm. The receiver features 12 channels that enable measurements of elastically (three channels) and Raman scattered light (387 and 607 channels for aerosols, 407 for water vapor) as well the depolarization of the incoming light at 355 and 532 nm. It also performs near-range measurements of two elastic and two Raman channels. More details about the instrument and its measurements are provided in Engelmann et al. (2016) and Baars et al. (2016). In brief, the nighttime backscatter (β p ) and extinction (α p ) coefficient profiles at 532 nm are derived using the Raman method proposed by Ansmann et al. (1992). The volume and particle depolarization ratio profiles are derived using the methodologies described in Freudenthaler et al. (2009) andFreudenthaler (2016). The daytime backscatter and extinction coefficient profiles are derived using the Klett-Fernald method (Klett, 1981;Fernald , 1984), assuming a constant value for the lidar ratio (LR). The daytime Klett profiles in Sect. 4.1 were derived using a lidar ratio of 50 sr on 15 April and of 40 sr on 5, 9, 21 and 22 April as well as a vertical smoothing length using a sliding average of 232.5 m. The integrated extinction coefficient profiles calculated with these LRs agree well with the collocated AERONET aerosol optical depth (AOD) observations. The LR values also are in agreement with the nighttime Raman measurements indicating mixtures of dust and anthropogenic/continental particles at heights between 1 and 3 km. The 2-D backscatter coefficient curtain for Fig. 4 is calculated with the methodology described by Baars et al. (2017).
In this work we also use spaceborne observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite (Winker et al., 2009). During the campaign period CALIPSO passed over Nicosia at a distance of 5 km on 5 and 21 April 2016. Here, we use the CALIPSO L2 version 4 (V4) aerosol profile products of 21 April 2016 and consider only qualityassured retrievals Tackett et al., 2018).

INP retrieval from lidar measurements
We calculated the n INP profiles from the lidar measurements by first separating the lidar backscatter profile into its dust and nondust components using the aerosol-type separation technique introduced by Shimizu et al. (2004) and . For this method we consider a dust particle linear depolarization ratio of δ d = 0.31±0.04 Ansmann et al., 2011a) and a nondust particle linear depolarization ratio of δ nd = 0.05 ± 0.03 (Müller et al., 2007;Groß et al., 2013;Baars et al., 2016;Haarig et al., 2017). The observed particle linear depolarization ratio in between these marginal values is therefore attributed to a mixture of the two aerosol types. The dust extinction coefficient (α d ) is calculated using the mean LR of 45 ± 11 sr for dust transported to Cyprus (Nisantzi et al., 2015). For the nondust component, the extinction coefficient (α c ) is calculated using a LR of 50 ± 25 sr which is representative for nondesert continental mixtures (Mamouri and Ansmann, 2014;Kim et al., 2018). The profiles of n 250,d,dry , S d,dry , n 250,c,dry and S c,dry are calculated from the extinction coefficient profiles using the POLIPHON algorithm (POlarization-LIdar PHOtometer Networking) and AERONET-based parameterizations proposed by Ansmann (2015, 2016). Table 3 provides an overview of the corresponding formulas used for the calculations. Weinzierl et al. (2009) showed that for dust environments the AERONET-derived values of S dry are about 95 % of the total particle surface area concentration (including particles with radius < 50 nm). This assumption has been validated against airborne in situ observations of the particle size distribution during the Saharan Mineral Dust Experiment (SAMUM; Ansmann et al., 2011b) in Morocco. The correlation drops to ∼ 0.85 ± 0.10 for urban environments based on ground-based in situ measurements of particle size distributions at the urban site of Leipzig .
The uncertainty in the products (considering the initial errors provided in Table 2) are as follows: the estimated n 250,d,dry uncertainty is 30 % in well-detected dessert dust layers (δ d = 0.3), 37 % in less pronounced aerosol layers (δ d = 0.2) and exceeds 94 % in aerosol layers with low dust contribution (δ d < 0.1). The uncertainty of the estimated S d,dry values is 38 % in well-detected dessert dust layers, 44 % in less pronounced aerosol layers and exceeds 97 % in aerosol layers with low dust contribution. The overall uncertainties of the combined (dust and continental) n 250,dry and S dry values are between 20 % and 40 % and between 30 % and 50 % respectively. The steps of the procedure for obtaining the profile of n 250,dry and S c,dry , as described here, are illustrated in an example in Fig. 2. In this example, we use the PollyXT measurements at Nicosia between 01:00 and 02:00 UTC on 21 April 2016.
In the final step, the n INP profiles are estimated using the ice nuclei parameterizations presented in Sect. 2 (Eqs. 1-7). For these calculations we are using collocated modeled profiles of the pressure, temperature and humidity fields. Specifically, for the PollyXT-based n INP calculations we use hourly outputs from the Weather Research and Forecasting atmospheric model (WRF; Skamarock et al., 2008) which is operational at the National Observatory of Athens at a mesoscale resolution of 12 km × 12 km and 31 vertical levels (Solomos et al., 2015(Solomos et al., , 2018. Initial and boundary conditions for the atmospheric fields and the sea surface temperature are taken from the National Centers for Environmental Prediction (NCEP) global reanalysis at 1 • × 1 • resolution. For the CALIPSO-based n INP calculations we use the track-collocated meteorological profiles from the MERRA-2 model (Modern-Era Retrospective analysis for Research and Applications, version 2) which are included in the CALIPSO V4 product .

UAV in situ measurements
Two fixed-wing UAVs, the Cruiser and the Skywalker, performed aerosol measurements up to altitudes of 2.5 km a.g.l.
(2.85 km a.s.l.). Both UAVs were used to collect INP samples onto silicon wafers using electrostatic precipitation. The Cruiser can carry a payload of up to 10 kg, and it was equipped with the multi-INP sampler PEAC (programmable electrostatic aerosol collector) (Schrod et al., 2016). Sky- Figure 2. PollyXT profiles of the total particle backscatter coefficient (purple) and particle linear depolarization ratio (green) measured between 01:00 and 02:00 UTC on 21 April 2016. The extinction coefficient as well as the number and surface concentration of particles with a dry radius larger than 250 nm related to mineral dust (orange) and nondust aerosol (black) was obtained following the methodology described in Sect. 3.2.  walker X8 (a light UAV that can carry a payload of 2 kg) was equipped with a custom-built, lightweight version of a single-sampler PEAC (Schrod et al., 2017 (Schrod et al., 2016(Schrod et al., , 2017. FRIDGE is an isostatic diffusion chamber. The typical operation of FRIDGE allows for measurements at temperatures down to −30 • C and relative humidity with respect to water (RH w ) up to water supersaturation. FRIDGE was originally designed to address the condensation and deposition freezing ice nucleation modes at water saturation and below. How-ever, because condensation already begins at subsaturation, its measurements at a RH w between 95 % and 100 % encompass ice nucleation by deposition nucleation plus condensation/immersion freezing, which cannot be distinguished by this measurement technique. Recent measurements during a large-scale intercomparison experiment with controlled laboratory settings showed that the method compares well to other INP counters for various aerosol types (DeMott et al., 2018). However, sometimes FRIDGE measurements are on the lower end of observations when compared to instruments that encompass pure immersion freezing. The INP samples collected on 5, 15 and 21 April 2016 were used for comparison with the lidar-derived n INP . The samples were analyzed at −20, −25 and −30 • C and at a RH w of 95 %, 97 %, 99 % and 101 % with respect to water, or equivalently 115 % to 135 % with respect to ice (RH ice ) (Schrod et al., 2017). Hereon, the samples analyzed at a RH w < 100 % are used as a reference for the deposition mode parameterizations, and the samples analyzed at a RH w of 101 % are used as a reference for the immersion/condensation parameterizations. The errors of the INP measurements were estimated to be ∼ 20 % considering the statistical reproducibility of an individual sample, for the samples analyzed for the experiment.
Cruiser was additionally equipped with an optical particle counter (OPC, Met One Instruments, Model 212 Profiler) that measures the aerosol particle number concentration with 1 Hz resolution in eight channels ranging from 0.15 to 5 µm in radius (Mamali et al., 2018). The inlet of the OPC was preheated to keep the relative humidity below 50 % to minimize the influence of water absorption. The Cruiser-OPC measurements on 5, 9, 15 and 22 April 2016 were used to calculate the n 250,dry profiles discussed in Sect. 4.1.
The measurements from the OPC on board the Cruiser UAV were validated at the ground, using a similar OPC and a differential mobility analyzer (DMA). The first comparison showed underestimation for the bin with radius 1.5 µm to 2.5 µm and for the last bin with radius more than 5 µm. The second comparison showed that the OPC underestimates by less than 10 % the number concentration of particles with radius between 0.15 µm and 0.5 µm (Burkart et al., 2010). Moreover, there are no data provided for particles with radius less than 0.15 µm. In order to correct for this undersampling we fit a bimodal number size distribution on the in situ data and derive a corrected n 250,dry and S dry . An example of this correction is shown in Fig. 3 for the number and surface size distributions measured at 1.2 km on 5 April 2016. For the cases discussed herein we found that the corrected n 250,dry in situ values were ∼ 20 % higher than the raw measurements.

Spaceborne cloud observations
A-Train spaceborne cloud observations are complementarily used to provide us the 3-D distribution and characteristics of the clouds formed in the presence of the calculated n INP . For the spatial distribution of the clouds formed during 21 April 2016, the true color observations from the MODIS instrument (Moderate Resolution Imaging Spectroradiometer) on board Aqua satellite are used (available from NASA at https://worldview.earthdata.nasa.gov/, last access: 8 August 2019). To get a better insight into the vertical cloud structure, we use outputs from the synergistic radar-lidar retrieval DARDAR (raDAR/liDAR; Delanoë and Hogan , 2008). The DARDAR retrieval (initiated by LAT-MOS and the University of Reading) uses collocated Cloud-Sat, CALIPSO, and MODIS measurements and provides a cloud classification product (DARDAR-MASK; Ceccaldi et al., 2013) and ice cloud retrieval products (DARDAR-Cloud; Delanoë et al., 2014) on a 60 m vertical and 1.1 km horizontal resolution (available at http://www.icare.univ-lille1. fr/projects/dardar, last access: 8 August 2019). In this work, we use the DARDAR-MASK product for cloud classification, and we utilize the DARDAR-Cloud product to derive an estimation of the ice crystal number concentration (n ice ) of the scene. With increasing maximum diameter (D max ), the ice crystals become more complex and their effective density decreases (Heymsfield et al., 2010). The DARDAR algorithm describes this relationship using a combination of in situ measurements by Brown and Francis (1995) for lowdensity aggregates (D max > 300 µm) and by Mitchell (1996) for hexagonal columns (D max < 300 µm). We derive the n ice (DARDAR-Nice) following the approach presented by Sourdeval et al. (2018) on the DARDAR-Cloud parameters of the ice water content (IWC) and the normalization factor of the modified gamma size distribution (N * 0 ). The direct propagation of uncertainties for IWC and N * 0 provided by DARDAR-Cloud gives an estimate for the relative uncertainty in n ice from about 25 % in lidar-radar conditions to 50 % in lidaronly or radar-only conditions (Sourdeval et al., 2018). This estimation accounts for instrumental errors and uncertainties associated with the a priori profiles used in DARDAR-Cloud. In cases with high homogeneous nucleation rates or dominant aggregation processes, N i can be underestimated (respective overestimated) by an additional 50 % due to deviations from the assumed particle size distribution. Due to further assumptions within DARDAR-Cloud (e.g., a fixed mass-dimensional relationship), additional uncertainties can increase the error of the retrieved n ice . In Sect. 4.3, the retrieved n ice is only used as a hint to estimate the order of magnitude of the true n ice .

Results and discussion
We present here the comparison between the UAV-OPC observations and the lidar-derived n 250 profiles (Sect. 4.1). The measurements used for this comparison corresponds to one intense dust event, where the UAV measurements were conducted under cloudy conditions (9 April), and three moderate events with dust and continental mixtures, where the UAV flights were conducted under cloud-free conditions (5, 15 and 22 April). Subsequently, we present the comparison between the UAV-INP measurements and the lidar-derived n INP during 3 d with moderate dust load conditions (Sect. 4.2). From a total of six INP samples, one sample is collected during 21 April in the presence of a pure-dust event under cloudy conditions, and the remaining five samples are collected during 5 and 15 April inside dust and continental aerosol layers under cloud-free conditions. A brief description of the aerosol conditions of the measurements used is provided herein.
On 5 April 2016, a homogeneous elevated dust layer was observed above the lidar station at 1.0-1.8 km from 00:00 to 08:00 UTC, which was later on mixed into the developing planetary boundary layer (PBL). In the next hours (until 12:00 UTC), only moderate variability was observed above the station (in the lidar backscatter coefficient and δ p curtains -not shown). The UAV samples were collected between 11:37 and 11:57 UTC at 30 km west of the lidar site with westerly winds prevailing. Constant δ p of around 0.15 between 0.5 and 2.5 km supports the qualitative homogeneity between the two observation sites during this time period.
On 9 April 2016, a thick pure-dust layer (with δ p ≈ 0.3) was observed above the lidar station, as part of a major dust event above Cyprus between 8 and 11 April 2016. The mean AOD at Nicosia was 0.83 (at 500 nm) with a corresponding mean Ångström exponent of 0.17 (at 440-870 nm). During the event, ice and water clouds were frequently formed at the top of the dust layer (mainly between 3 and 6 km). The DREAM model and backward trajectory analysis revealed that this event originated from the central Sahara, with the dust particles being advected by a southwesterly flow directly towards Cyprus, reaching the island after 1 d (Schrod et al., 2017). The UAV samples were collected between 08:12 and 08:23 UTC inside the dust layer, and these observations were compared with the lidar-derived profiles at 06:50-06:59 UTC (a closer-in-time collocation between the lidar observations and the UAV measurements is not possible due to clouds with a cloud base at 4 km later on). The OPC concentrations collected that day were the highest observed during the period of the INUIT-BACCHUS-ACTRIS experiment.
On 15 April 2016 a persistent elevated dust layer was observed above Nicosia. Backward trajectory analysis (not shown) revealed that this dust event originated from Algeria and that the dust plume was transported over Greece and Turkey before reaching Cyprus. Cruiser UAVs collected samples between 06:54 and 08:45 UTC (during the boundary layer development). At that time, a pure-dust layer (δ p ≈ 0.3) was present between 2.5 and 3.8 km height. Below 2.0 km the dust was mixed with continental spherical particles from the residual layer with δ p decreasing with height (reaching ∼ 0.1 at 0.6 km). During the 2 h flight, the scene above the station changed considerably, with a 31 % increase in the aerosol optical thickness (from 0.33 to 0.48) and 16 % decrease in the Ångström exponent (from 0.31 to 0.26). The UAV measurements that day reached heights of up to 2.2 km, thus capturing only the mixed bottom layer and the lower part of the elevated dust layer. For the comparison with the lidar-derived concentrations, only the UAV measurements inside the lower part of the elevated dust layer (1.7-2.2 km) are used.
The pure-dust event on 20 to 21 April 2016 is considered the golden case of our dataset, as it has been observed simultaneously with the PollyXT lidar, the UAVs and the A-Train satellites. Additionally, it is the only pure-dust event of our dataset where we have simultaneously good lidar observations and in situ INP measurements. Figure 4 provides an overview of the times and heights of the PollyXT measurements, along with the CALIPSO overpass and UAV measurement times, between 20 and 22 April 2016. During that period atmospheric conditions supported the transport of dust from the Saharan desert and the Arabian Peninsula to the Eastern Mediterranean (δ p = 0.28±0.03) (Floutsi, 2018). The elevated dust plume arrived over the lidar site at 4-5 km height (∼ 15:00 UTC on 20 April 2016), quickly widened to stretch from 2 to 8 km height with the top of the main plume at 5 km height, and disappeared at 18:00 UTC on 21 April. On that day, ice clouds were formed within the dust plume and were present between 02:00 and 10:45 UTC above Nicosia. As shown in the figure, UAV flights were performed inside the dust layer on 21 April 2016 (OPC measurements and INP sampling). The event was captured from the A-Train satellites at 11:01 UTC (CALIPSO overpass time). Figure 5 provides an overview of the aerosol and clouds above the area, with the MODIS true color image (upper panel) and the combined DARDAR and CALIPSO L2 feature mask (lower panel). Dust is observed above the broader region at altitudes up to 6 km, and ice clouds are formed inside the dust layer south of Cyprus at altitudes greater than 4 km (T < 0 • C). The ice clouds are detected and characterized at 1 km horizontal resolution (DARDAR-MASK product), while the dust plume is detected at 20 and 80 km horizontal resolution (CALIPSO L2 product). On 22 April 2016 a transported plume was detected between 03:00 and 10:00 UTC, at altitudes of 1 to 2 km above Cyprus. The layer consisted of a mixture of dust with pollution aerosol and is characterized by a homogeneous particle linear depolarization ratio of δ p = 0.17 ± 0.03. UAV flights (OPC and INP sampling) were performed in the mixed layer during that day between 04:32 and 05:13 UTC (Fig. 4).
All in situ samples were collected at a location about 28 km to the west of the lidar site; thus the atmospheric homogeneity of the two areas had to be considered to select suitable measurement times for the comparisons. For this analysis we used the sun-photometer measurements at Agia Marina and Nicosia, backward trajectories, model fields and MODIS measurements. This was especially necessary for the case on 21 April when clouds were formed at the top of the dust layer. During that day, the CALIPSO-derived n INP at 11:01 UTC were compared to UAV-measured ones acquired approximately 1.5 h earlier (between 08:30 and 09:40 UTC). The space and time homogeneity of the CALIPSO-derived s dry and n 250,dry profiles (acquired shortly after the end of the cloudy period) is confirmed by the respective estimates from the PollyXT measurements during 01:00 to 02:00 UTC (before the beginning of the cloud formation) as shown in Fig. 6. The different measurement times of the ground-based and spaceborne lidars are marked in Fig. 4. For the CALIPSO profiles, along-track observations ±80 km away from the lidar station are used. During that time, the dust plume declined by approximately 300 m. Nevertheless, CALIPSOand PollyXT-retrieved profiles are in agreement within their error bars within the dense dust plume. Aerosol conditions were less homogeneous above and below this layer (see Fig. 4), causing stronger differences between the profiles of the four parameters from the two instruments. The comparison between the CALIPSO-derived n INP and the UAV measurements from this case is discussed in Sect. 4.2 (see Fig. 9). Table 3. Overview of the AERONET-based parameterizations used in this study for the conversion of the measured optical aerosol properties (α d , α c ) into the microphysical properties (n 250,d,dry , S d,dry , n 250,c,dry and S c,dry ). The parameterizations were introduced in Mamouri and Ansmann (2016). In the equations, α is in per megameter (Mm −1 ), c 250 in Mm cm −3 , c s in Mm m 2 cm −3 , n 250,dry in cm −3 and S dry in m 2 cm −3 . For the values of the conversion parameters (c 250,d , c s,d , c 250,c and c s,c ) see Table 2.

Evaluation of the n 250,dry retrieval
For the assessment of the lidar-based n 250 retrieval we used the OPC measurements on 5, 9, 15 and 22 April. The profiles of n 250,dry retrieved from PollyXT observations and in situ measurements are shown in Fig. 7a. The lidar dust-only profiles (orange lines) are calculated from the dust extinction profiles and Eq. (8) ( Table 3). The remaining nondust component is considered continental with n 250,c,dry provided by Eq. (10) ( Table 3). The total n 250,dry profiles (Fig. 7a, black lines) are the summation of n 250,d,dry and n 250,c,dry . The red dots correspond to the uncorrected UAV n 250,dry measurements. The blue dots correspond to the corrected UAV n 250,dry measurements (as described in Sect. 3.3). We use only the respective height ranges at which homogeneous aerosol conditions allow for a comparison of the UAV-and lidar-derived estimates. These measurements correspond to heights above 0.5 km on 5 April, above the PBL on 9 and 15 April (> 1 and > 2 km respectively), and above the nocturnal boundary layer on 22 April (> 0.7 km). It seems that the distance has little impact on the lidar-derived and the in situ-measured n 250,dry presented in Fig. 7, with most of the in situ-derived n 250,dry being well within the error bars of the lidar retrieval when considering the contributions of both mineral dust and continental pollution. On 9 April we observed the highest differences between the lidar-derived and in situ-measured n 250,dry , which may be attributed to the ∼ 1 h time difference between the in situ sampling and the lidar retrieval (limitation due to mid-level clouds as discussed already). Nevertheless, the case is included here, as it represent the strongest dust event observed during the campaign. Overall, the values of n 250,dry varied between 1 and 50 cm −1 . Figure 8 provides a quantitative comparison of the observations presented in Fig. 7 for lidar retrievals of n 250,dry considering both mineral dust and continental pollution and the corresponding in situ measurements at the same height levels. Again, we see that the results agree well within the error bars of the lidar retrieval with R 2 = 0.98. The uncertainties of the UAV-derived n 250,dry values presented in Figs. 7 and 8 correspond to the standard deviation of the 30 s average (OPC initial resolution of 1 s). The error in the OPC data due to the assumption of the refractive index and the shape of the particles used for the derivation of the particle size distribution from the OPC measurements were not taken into account in this study. Nevertheless, it is not expected to be high because the refractive index used is characteristic for dust particles (n = 1.59). We have to keep in mind the effect of a possible inhomogeneity between the two stations. In view of all uncertainty sources, the lidar-and UAVderived n 250,dry are in good agreement. In terms of absolute values, the lidar-derived n 250,dry are slightly lower than the UAV-derived ones. We conclude that lidar measurements are capable of providing reliable spatiotemporal distributions of n 250,dry in cases with dust and continental aerosol presence with an uncertainty of 20 to 40 %.
The profiles of S dry retrieved from PollyXT observations and in situ measurements are shown in Fig. 7b. The dust-only profiles (orange lines) are calculated from the dust extinction profiles and Eq. (9) ( Table 3). The remaining nondust component is considered continental with n 250,c,dry provided by Eq. (11) ( Table 3). The total S dry profiles (Fig. 7b, black lines) are the summation of S d,dry and S c,dry . These profiles are compared to the total S dry derived from the corrected in situ number size distribution (e.g., Fig. 3b). We see that the latter agree well within the uncertainty of the lidar-derived S d,dry (orange line) but do not agree well when both mineral dust and continental pollution are considered (black line). This is mainly due to the sampling cutoff of the OPC instru-ment for particles with radius smaller than 150 nm, which are mainly composed of polluted continental particles. The effect is not seen in the corrected n 250 , since the size ranges considered there are larger than 250 nm.

Evaluation of the n INP retrieval
For the assessment of the lidar-based n INP retrieval, the UAV measurements on 5, 15 and 21 April 2016 are used. The samples of 5 and 15 April were collected under the moderately mixed dust and continental conditions shown in Fig. 7. On 5 April, the sample was collected at an altitude of 1.823 km (δ p = 0.14 ± 0.02). On 15 April two samples were collected from a 0.998 km and 1.281 km altitude (δ p = 0.15±0.03). On 21 April, the pure-dust sample was collected from a 2.55 km altitude (δ p = 0.28 ± 0.03) (Fig. 4). Analysis performed in FRIDGE chamber provided the INP concentrations for these cases. The in situ samples were analyzed at −20, −25 and −30 • C. For the deposition nucleation (Figs. 9a and 10a), the samples were analyzed at a RH w of 95 %, 97 % and 99 %, leading to three values of S ice for each temperature (1.16, 1.18 and 1.23 for −20 • C; 1.21, 1.24 and 1.26 for −25 • C; and 1.27, 1.30 and 1.33 for −30 • C). For the immersion freezing (Fig. 9b), the samples were analyzed at a RH w of 101 %, leading to S ice of 1.23, 1.29 and 1.35 for the temperatures of −20, −25 and −30 • C, respectively. For T = −20 • C, RH w = 101 % and S ice = 1.23, we refer to the freezing process as condensation freezing.
The sample of 21 April was analyzed by single-particle analysis using a scanning electron microscope, which showed that 99 % of the particles were dust and 1 % was Ca sulfates and carbonaceous particles (Schrod et al., 2017). This sample is used in order to evaluate the performance of the n INP lidar estimates in a pure-dust case, where (i) the errors originating from the first step of our methodology (separation in dust and nondust aerosol components) are small (∼ 30 %) and (ii) the uncertainties induced from the D10 and U17(soot) parameterizations are minimal. Figure 9 shows the n INP on 21 April as they were calculated from the CALIPSO lidar measurements (colored symbols) and measured from the UAV-FRIDGE samples (black triangles), (panel a) for deposition nucleation (as a function of saturation over ice) and (panel b) for condensation and immersion freezing (as a function of temperature).
Likewise, we are using all the aforementioned cases in order to evaluate the performance of the n INP lidar estimates in cases with dust and continental aerosols. Figure 10 shows scatter plots of all the lidar-estimated n INP (from PollyXT and CALIPSO) against the in situ measurements for (panel a) deposition nucleation and (panel b) condensation and immersion freezing. In Fig. 10b the ratio between the lidar-derived and the in situ n INP is provided as a function of temperature. Similar results are observed for both the pure-dust (Fig. 9) and the dust and continental cases (Fig. 10), with the lidarestimated n INP during the pure-dust event to show the best agreement with the in situ measurements.
For the n INP retrievals in the deposition mode we see that using the U17-dep in a dust case the lidar-derived concentrations are in excellent agreement with the in situ observations (well within their uncertainties), with n INP values to span over 2.5 orders of magnitude (for different ice supersaturation conditions) and the retrievals to capture the whole extent of this range (Fig. 9a). The lidar-retrieved U17-dep values in this case are dominated by the dust-related n INP (estimated from Eq. 3; Table 1), with the nondust-related n INP (estimated from Eq. 6; Table 1) being 5 orders of magnitude lower. In dust and continental cases (Fig. 10a), 97 % of all the U17-dep lidar-derived n INP are within the error bars of the in situ measurements and within a factor of 10 around the 1 : 1 line (r = 0.75). The n INP sampled with the UAVs ranged between 0.02 and 20 L −1 . Using S15 parameterization, the predicted n INP values are 3 to 5 orders of magnitude larger than the in situ measurements in both dust and dustcontinental cases (r = 0.42). An overestimation was already expected as discussed in Sect. 2 and Steinke et al. (2015), but for completeness we include these results. Figures 9b and 10b show the lidar-derived immersion/condensation INPs. U17-imm dust-related n INP values are calculated using the INP parameterization of Eq. (1) (Table 1) with the S d,dry from Eq. (9) ( Table 3). The D15 dustrelated n INP are calculated using Eq. (2) ( Table 1) with the n 250,d,dry from Eq. (8) ( Table 3). The D10 continental-related n INP are calculated using Eq. (7) ( Table 1) with the n 250,c,dry from Eq. (10) ( Table 3). The D15 + D10 values for the total (dust + continental) aerosol in the scene are the summation of the aforementioned D15 (dust-related) and D10 (continental-related) n INP calculations (See Figs. A1 and A2 in Appendix). We did not include the U17-imm soot estimates in the plot since these are quite similar to the estimated values from D10 at temperatures < −18 • C (Sect. 2; Fig. 1). Consequently, for the total INP load in the scene, the estimations provided from D15 + D10 are similar to the ones provided from D15 + U17-imm(soot). In the rest of this paper, we will discuss only the joint D15 + D10 estimates, keeping in mind that the same conclusions apply for the joint D15 + U17-imm(soot) estimates.
In Figs. 9b and 10b we see that the lidar-derived n INP using D15 for dust and D10 for continental particles are in good agreement with the in situ observations, within the respective uncertainties for the samples analyzed at −20 and −25 • C. The best n INP agreement is observed for the puredust sample analyzed under condensation freezing conditions (at −20 • C): with in situ measurements of 3.6±0.1 L −1 and lidar-derived D15 + D10 estimates of 3.8 L −1 . From them, 2.4 L −1 originated from the D15 dust contribution and 1.4 L −1 from the D10 nondust contribution (although the contribution from the nondust INP at lower temperatures was insignificant with nondust concentrations 1 order of magnitude lower than the dust ones). Using all the dust and continental cases we see that, for the samples analyzed under condensation freezing conditions, the D15 + D10 estimated n INP are no more than 2.5 times higher than the in situ measurements (Fig. 10b). Larger differences are observed at the temperatures where immersion freezing dominates over condensation as the main INP pathway, with 1.5-7 times larger values at −25 • C and 4-13 times larger values at −30 • C. Indicatively, for the pure-dust case, at T = −25 • C the in situ n INP were 12 ± 3 L −1 and the D15 + D10 lidar-derived n INP were 26 L −1 (with a negative error of 14 L −1 ). At T = −30 • C, the in situ n INP were 62 ± 14 L −1 while D15 + D10 n INP estimates were 1 order of magnitude higher (242 L −1 ). Overall, in 85 % of the analyzed cases, the D15 + D10 lidar retrievals are less than an order of magnitude higher than the UAV measurements. Regarding the U17-imm lidarderived n INP values, they are overall 1 to 3 orders of magnitude higher than the in situ ones. In particular they are 3-11, 2-80 and 2-1000 times larger than the samples analyzed at FRIDGE chamber at −20, −25 and −30 • C, respectively. Nevertheless, the in situ observations are within the uncertainty of the parameterization for all the cases. Indicatively, for the pure-dust case, the U17-imm lidar-derived n INP values are 50 L −1 at T = −20 • C. Recent comparisons of n INP derived from samples analyzed in the FRIDGE chamber usually present good linear correlations but somewhat lower values with observations derived from pure immersion paths (e.g., D15) (DeMott et al., 2018). Possible reasons for these discrepancies may be (a) deficits and inadequacies in instrumentation and measurement techniques, (b) the lacking overlap of the freezing modes, (c) inconsistencies between the inlet systems of the parameterization measurement (using cutoffs) and the in situ measurements (using no cutoff), and (d) a variation in RH w (D15: 105 %; FRIDGE: 101 %) (Schrod et al., 2017).
The error bars of the lidar-based n INP estimations in Figs. 9 and 10 are calculated using Gaussian error propagation together with the typical uncertainties provided in Table 2. In DeMott et al. (2015), a standard deviation of 2 orders of mag-  nitude is reported as the uncertainty of the D15 parameterization. In the same plots, the uncertainty of the n INP from in situ data is very low. Under most experimental conditions, the repeatability of the ice nucleation in the FRIDGE chamber dominates other uncertainties. An uncertainty of 20 % has been suggested as a useful guideline for the uncertainty of the intrinsic measurements, corresponding to the statistical reproducibility of an individual sample. However, it has also been reported that natural variability by far outweighs the intrinsic uncertainty (Schrod et al., 2016). We need to consider the full uncertainty including precision and accuracy. The DeMott et al. (2018) intercomparison of INP methods saw that at all temperatures and for various test aerosols the n INP uncertainty for immersion freezing is 1 order of mag-nitude, while for deposition condensation the uncertainty is expected to be even larger.
Our analysis suggests that the D15 + D10 (and D15 + U17-imm(soot)) immersion/condensation parameterization (applicable for the temperature range −35 to −9 • C) and the U17-dep parameterization (applicable for the temperature range −50 to −33 • C) agree well with in situ observations of n INP and can provide good n INP estimates in pure-dust and dust-continental environments. The U17-imm pure immersion parameterization provides values 1-2 orders of magnitude larger; we therefore consider the n INP estimates according to D15 + D10 as the lower boundary of possible values, with the actual values being up to 1 order of magnitude larger in the temperature regime of immersion freezing.

n INP profiles from PollyXT and CALIOP during the evolution of mixed-phase clouds in a Saharan dust event
The case study of 21 April 2016 demonstrates the feasibility of the proposed methodology to provide profiles of cloud-relevant aerosol parameters up to the cloud levels, using (ground-based and spaceborne) lidar measurements. In particular for this case, the temporarily averaged PollyXT lidar observations at 01:00-02:00 UTC and the spatially averaged CALIPSO observations at 11:01 UTC provide us the information of the n 250,dry , S dry and n INP right before and after the cloud event which was formed inside the dust layer that day between 02:00 and 10:45 UTC. The profiles of n 250,dry and S dry before (PollyXT) and after (CALIPSO) the cloud event are the ones already presented in Fig. 6. As discussed above, the dust plume declined by approximately 300 m during that period while its n INP stayed relatively constant inside its dense part. Above the main dust layer the aerosol conditions were variable, with multiple thin layers present up to 8 km altitude only before the appearance of the clouds. Specifically, a contribution of nondust continental particles is observed between 5.6 and 8 km a.g.l.
at altitudes between 6.3 and 8 km (temperatures < −21 • C). These abundances are reversed inside the main dust layer (altitudes between 4 and 5.5 km; temperatures: [−20,−6] • C) where 34 % of the total n INP is attributed to dust aerosols (0.06 L −1 ) and 66 % to nondust aerosols (0.12 L −1 ). Shortly after the period analyzed here, mixed-phase clouds are observed above Nicosia at first at altitudes between 5 and 7 km and during the rest of the cloudy period mainly above 4 km (Fig. 4). Figure 11b show the lidar-derived n INP above the station shortly after the end of the cloudy conditions. At that time, the main dust layer is observed at altitudes up to 5.5 km without additional layers above it. These observations are close to the local noon, with the air temperature above the station being increased by 2.7 • , leading to temperatures of 0 • C at 3.6 km and −15 • C at 5.4 km a.g.l. At these altitudes, a relatively constant contribution of nondust particles is present (n 250,dry = 0.4 ± 0.2 cm −3 ; Fig. 6d), which leads to a gradual increase in the nondust n INP per height from 2 × 10 −4 L −1 (4 km; −2 • C) to 10 −2 L −1 (4.4 km; −5 • C) to 0.2 L −1 (5.3 km; −12 • C) (D10 estimates). Additionally, the dust concentration per altitude is constant inside the dust layer and is decreased gradually above 4.6 km (n 250,dry = 16 cm −3 ; 4-4.6 km; Fig. 6c). The dust-related n INP per height are 8×10 −3 L −1 (4 km; −2 • C), 3×10 −3 L −1 (4.4 km; −5 • C) and 0.1 L −1 (5.3 km; −12 • C) (D10 estimates). Overall, 25 % of the total n INP is attributed to dust aerosols (D15) and 75 % to nondust aerosols (D10) at altitudes between 3.8 and 5.6 km. Taking into consideration all the aerosols, the n INP before and after the cloud development is ∼ 0.6 L −1 and ∼ 0.1 L −1 respectively at 5.3 km altitude (D15 + D10 in Fig. 6). This difference is due to the increase in the air temperature during the day and the decrease in n 250,dry and S dry . Before the cloud formation, the n INP values at [6, 7.5] km are 1 order of magnitude larger than at 5.3 km (∼ 3 L −1 ), and they are 2 orders of magnitude higher at 7.8 km than at 6 km (200 L −1 ). These results indicate that the particles in the main dust layer and the thin layers above it acted as seeding INPs for the cloud that formed in that layer, affecting also its characteristics. However, further measurements are necessary to reach a more concrete conclusion, for example, measurements of the atmosphere dynamics (e.g., from a wind lidar) and observations of the cloud evolution (e.g., from a cloud radar as in the recent study of Ansmann et al., 2019). Although these measurements are absent from our ground-based instrumentation, we utilize the DARDAR-Nice product (based on the CLOUDSAT and CALIPSO observations on 21 April 2016 - Fig. 5) as a hint for the true n INP of the scene, and we compare them with the neighboring CALIPSO n INP estimates. Figure 12 shows the DARDAR n ice estimations along the A-Train track (presented in Fig. 5), and Fig. 13 shows the n INP calculations on the same curtain using the D15 + D10 (upper panel) and U17-imm (lower panel) parametrizations. Clouds are formed on top of the dust layer at latitudes of 32, 32.8 and 34 • N. The clouds observed at 32 and 32.8 • N are coupled with an aerosol layer at their cloud top, at altitudes of 6.3 and 7.3 km and temperatures of −18 and −25 • C respectively. Figure 14 shows the n ice profiles derived in these two clouds, along with the n INP profiles estimates in their vicinity. Due to the strong INP number increase with deceasing temperature, the highest n INP concentrations are observed at the top of the upper aerosol-cloud layers. We assume that the ice crystals in these two clouds nucleate close to the cloud top (where the coldest temperatures are observed) and that afterwards the crystals grow and fall through the lower heights of the clouds formed. Moreover, we consider that no secondary ice production (SIP) processes are present in these clouds, or at least their contribution to the n ice is insignificant, as the cloud top temperatures are much lower than the temperatures where SIP has been observed (between −3 and −8 • C) (Hallett and Mossop, 1974;Field et al., 2017;Sullivan et al., 2017Sullivan et al., , 2018. We compare the n INP at cloud top height with the n ice inside the cloud, having in mind that, with our hypotheses, the n ice values can be up to the n INP values if all the INPs are activated to ice crystals. For the smaller cloud, at ∼ 32 • N, n ice between 0.8 and 8 L −1 are retrieved, and n INP between 0.3 to 2 and 4 to 20 L −1 are estimated with the D15 + D10 and the U17-imm respectively. For the cloud at ∼ 32.8 • N, n ice between 0.4 and 60 L −1 are retrieved, and n INP between 3 to 20 L −1 and 100 to 400 L −1 are estimated with the D15 + D10 and the U17-imm respectively. Overall, in these two clouds the n INP estimates in the top of the clouds have an uncertainty of 1-2 orders of magnitude in their estimates and differences of 1 order of magnitude in the retrievals between each other. Additionally the retrieved DARDAR profiles provide us only with a hint of the order of magnitude of the true n ice . Nevertheless the n ice estimates are between the estimated n INP values and within the errors of the two parameterizations. These results strengthen our conclusion that we can use the lidar-derived n INP from D15 + D10 and U17-imm to estimate a minimum and maximum boundary of the n ice in a cloud formed in their presence, when immersion is the dominant mechanism.

Summary and conclusions
We present a methodology for deriving n INP profiles from lidar measurements and a comparison with in situ UAV measurements of n INP . More specifically, seven INP parameterizations are tested to obtain lidar (ground-based and spaceborne) n INP estimates representative of mineral dust and continental, pollution or soot aerosol. We prove that a compilation of the parameterizations of DeMott et al. The cloud-relevant aerosol parameters necessary for INP estimations (n 250,dry and S dry ) are derived from lidar measurements as shown by Ansmann (2015, 2016).
The comparison between the lidar-derived concentrations of dry particles with radii larger than 250 nm with coincident UAV-OPC in situ measurements showed a good agreement with slightly lower values (32 %) for the n 250,dry derived by the lidar. This effect is less pronounced at low concentrations with a squared correlation coefficient of 0.98. For the majority of the cases, we find that in situ observations and remotesensing estimates are in agreement within their uncertainty ranges.
A further step for improving the lidar-derived INP retrievals and investigating the different parametrizations used is by conducting dedicated studies with collocated lidar measurements and additional temperature and humidity profiling in order to calculate the INP concentrations at real conditions, as well as the combination of the retrieved n INP with airborne in situ ice concentration measurements.
Our methodology is validated for cases with dust presence. Additional measurements are required in order to define the optimum INP parameterizations for nondust atmospheric conditions (e.g., continental, marine, smoke). Future experimental INP campaigns with airborne in situ observations from aircrafts (including UAVs) collocated with lidar measurements at pure marine conditions and at mixed aerosol conditions could provide an ideal setup for an indepth investigation of the potential of the lidar-based INP profiles in complex and nondust atmospheric conditions.
The results presented in this study give us confidence to proceed to the next step, which is to combine cloud-relevant lidar aerosol and wind parameters and cloud radar heightresolved observations to monitor the evolution of clouds embedded in aerosol layers. This will provide a unique opportunity to better understand aerosol-cloud interactions in the field of heterogeneous ice formation.
Moreover, the study enhances the confidence for the production of global 3-D products of n 250,dry , S dry and n INP from the CALIPSO dataset. The application of our methodology to CALIPSO measurements of more than a decade long could provide valuable insight into the global heightresolved distribution of n 250,dry and n INP related to mineral dust, as well as possibly other aerosol types. This will enable worldwide studies of aerosol-cloud interactions to combine the new product with satellite radar observations (CloudSat) and the upcoming EarthCARE (Earth Cloud Aerosol and Radiation Explorer) mission.
A challenge of a new global INP climatology will be the assessment of its underestimation at high altitudes where it is known that CALIPSO observations can miss thin layers with small concentrations. A way to investigate the effect of the satellite-undetected layers in the n 250,dry , S dry and n INP CALIPSO products is the utilization of ground-based lidar network observations such as EARLINET and PollyNet. Data availability. The satellite products used in this study are the CALIPSO 5 km aerosol profile product (Vaughan et al., 2019) publicly available at the NASA Atmospheric Science Data Center (Nasa Atmospheric Science data center, 2019); the DARDAR-Mask v1.1.4 product (Ceccaldi et al., 2013) as well as the DARDAR-Cloud v2.1.1 retrieval product (Delanoë, et al., 2014), both publicly available on the AERIS/ICARE database (ICARE data and services center, 2019); and The MODIS Corrected Reflectance (True Color) images (Gumley et al., 2010) publicly available on the NASA Worldview center (NASA Worldview snapshots application center, 2019). The in situ INP data used in this study can be accessed through the BACCHUS database of INP observations (BACCHUS Ice Nucleation DataBase, 2019). The database is accessible to members only, but membership is free. The PollyXT observations (level 0 data of measured signals and level 2 data of backscatter coefficient and depolarization profiles), the in situ OPC measurement and the WRF modeled profiles used in this study can be accessed through the REACT database (Marinou et al., 2019). All datasets created during the calculation of the lidar-based number concentrations and the correction of the in situ OPC number concentrations will be provided upon request. The Fig. A1 illustrates the general idea of the methodology followed for the INP estimations from the PollyXT measurements. The equations for the conversions of the measured optical properties into the microphysical properties are provided in Table 3. The equations for the conversions of the microphysical properties to INPs are provided in Table 1. Figure A1. Overview of the data analysis scheme followed for the PollyXT measurements in this work. In the first step, we separate desert and nondesert backscatter coefficients (β d and β nd ) by means of the particle linear depolarization ratio (δ p ). The backscatter coefficients for the nondesert aerosol are estimated to be continental aerosol mixtures β c by means of, e.g., backward (BW) trajectory analysis and Ångström exponent information. The two backscatter coefficients are then converted to aerosol type-dependent particle extinction coefficients (α i ). In the next step, the extinction coefficients are converted to aerosol type-dependent profiles of particle number concentrations (n 250,i,dry ) and particle surface area concentration (S i,dry ). In the next step, ice-nucleating particle number concentrations (n INP,i ) are estimated by applying INP parameterizations from the literature indicated by D10, D15, S15 and U17 for DeMott et al. (2010DeMott et al. ( , 2015, Steinke et al. (2015) and Ullrich et al. (2017), respectively. Finally, the INP concentrations estimated for the different aerosol types are summed in order to estimate the total n INP .

A2 Methodological diagram for the analysis of the spaceborne lidar measurements
The Fig. A2 illustrates the general idea of the methodology followed for the INP estimations from the CALIPSO measurements. The equations for the conversions of the measured optical properties into the microphysical properties are provided in Table 3. The equations for the conversions of the microphysical properties to INPs are provided in Table 1. Figure A2. Overview of the data analysis scheme applied to CALIPSO measurements. In the CALIPSO case considered in this work only dust and polluted dust aerosol types have been observed. For that reason, only these combinations are considered here.
Author contributions. VA and EM conceived the presented idea. EM performed the analysis, drafted the manuscript, and designed the figures with support from MT and AT; AN guided and supervised EM on the ice nuclei mechanism and parameterizations and encouraged EM to investigate the errors of the lidar-derived INP concentrations. AA guided and supervised EM on the methodology and parameterizations for the lidar-derived INP concentrations and encouraged EM to investigate the errors of the in situ OPC concentrations and the DARDAR n ice concentrations. JS performed the analysis of the INP samples and guided EM on the limitations of these measurements. JS conceived and planned the experiment. MP was responsible for the UAV flights and the analysis of the OPC measurements. DM and AT performed the correction of the OPC measurements. RE supported essentially the PollyXT measurements during the campaign. HB supported the derivation of the PollyXT lidar backscatter coefficients. SS derived the WRF model profiles. FE derived the n ice from DARDAR. All authors provided critical feedback and helped shape the research, analysis and manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "EARLINET aerosol profiling: contributions to atmospheric and climate research". It is not associated with a conference.