Articles | Volume 21, issue 19
Atmos. Chem. Phys., 21, 14671–14686, 2021

Special issue: Ice nucleation in the boreal atmosphere

Atmos. Chem. Phys., 21, 14671–14686, 2021

Research article 05 Oct 2021

Research article | 05 Oct 2021

Two-year statistics of columnar-ice production in stratiform clouds over Hyytiälä, Finland: environmental conditions and the relevance to secondary ice production

Two-year statistics of columnar-ice production in stratiform clouds over Hyytiälä, Finland: environmental conditions and the relevance to secondary ice production
Haoran Li1,2, Ottmar Möhler3, Tuukka Petäjä2, and Dmitri Moisseev2,4 Haoran Li et al.
  • 1State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing, China
  • 2Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, Helsinki, Finland
  • 3Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology, Karlsruhe, Germany
  • 4Finnish Meteorological Institute, Helsinki, Finland

Correspondence: Haoran Li (


Formation of ice particles in clouds at temperatures of −10C or warmer was documented by using ground-based radar observations. At these temperatures, the number concentration of ice-nucleating particles (INPs) is not only expected to be small, but this number is also highly uncertain. In addition, there are a number of studies reporting that the observed number concentration of ice particles exceeds expected INP concentrations, indicating that other ice generation mechanisms, such as secondary ice production (SIP), may play an important role in such clouds. To identify formation of ice crystals and report conditions in which they are generated, W-band cloud radar Doppler spectra observations collected at the Hyytiälä station for more than 2 years were used. Given that at these temperatures ice crystals grow mainly as columns, which have distinct linear depolarization ratio (LDR) values, the spectral LDR was utilized to identify newly formed ice particles.

It is found that in 5 %–13 % of clouds, where cloud top temperatures are −12C or warmer, production of columnar ice is detected. For colder clouds, this percentage can be as high as 33 %; 40 %–50 % of columnar-ice-producing events last less than 1 h, while 5 %–15 % can persist for more than 6 h. By comparing clouds where columnar crystals are produced and to the ones where these crystals are absent, the columnar-ice-producing clouds tend to have larger values of liquid water path and precipitation intensity. The columnar-ice-producing clouds were subdivided into three categories, using the temperature difference, ΔT, between the altitudes where columns are first detected and cloud top. The cases where ΔT is less than 2 K are typically single-layer shallow clouds where needles are produced at the cloud top. In multilayered clouds where 2 K<ΔT, columns are produced in a layer that is seeded by ice particles falling from above. This classification allows us to study potential impacts of various SIP mechanisms, such as the Hallet–Mossop process or freezing breakup, on columnar-ice production. To answer the question whether the observed ice particles are generated by SIP in the observed single-layer shallow clouds, ice particle number concentrations were retrieved and compared to several INP parameterizations. It was found that the ice number concentrations tend to be 1–3 orders of magnitude higher than the expected INP concentrations.

1 Introduction

Ice production in mixed-phase clouds is critical for their radiative (Sun and Shine1994) and microphysical (Korolev et al.2017) properties. At temperatures warmer than −38C, ice crystals form on ice-nucleating particles (INPs). In situ measurements have revealed that the number concentration of available INPs decreases with the increase in ambient temperature (e.g., DeMott et al.2010; Wilson et al.2015; DeMott et al.2016; Petters and Wright2015; Schneider et al.2020). This dependence is more or less universal but can also be affected by other factors, such as the geographic location, air-mass types and aerosol compositions (e.g., DeMott et al.2010; Niemand et al.2012; Wilson et al.2015; DeMott et al.2016; Petters and Wright2015; McCluskey et al.2018). In addition, it has been found that INP concentrations at high latitudes are generally lower than at mid-latitudes (e.g., DeMott et al.2016; Wex et al.2019). Above −10C, the typical concentrations of INPs are below 10−1 L−1 and can be as low as 10−6 L−1 (Petters and Wright2015; Kanji et al.2017). A number of studies, however, have reported that the ice number concentration in clouds with the top temperature warmer than −10C can exceed the expected concentration of INPs by several orders of magnitude (e.g., Mossop1985; Hobbs and Rangno1985; Rangno and Hobbs2001). This discrepancy implies that numerical weather prediction models that rely solely on INP parameterizations cannot realistically represent ice number concentrations in moderately to slightly supercooled clouds. As a result, the inappropriate parameterization of ice production may lead to biased estimates of surface shortwave radiation budget (Young et al.2019), among other things (e.g., Zhao et al.2021; Zhao and Liu2021).

Several mechanisms have been proposed to explain this discrepancy, such as the enhanced contact nucleation driven by the thermophoretic force during the evaporation of liquid drops (Beard1992; Hobbs and Rangno1985), pre-activated INPs from evaporated ice particles nucleated above (Roberts and Hallett1968; Fridlind et al.2007) or secondary ice production (SIP) mechanisms (see recent reviews by Field et al.2017; Korolev and Leisner2020). The SIP has been studied by a number of laboratory experiments since the 1940s (e.g., Findeisen and Findeisen1943; Dye and Hobbs1966; Wildeman et al.2017). Hallett and Mossop (1974) found that numerous ice splinters can be generated when supercooled liquid drops larger than ∼25µm are collected by large ice particles within the temperature range of −8 to −3C. This is referred to as the Hallett–-Mossop (H–M) process, the most studied and most frequently implemented SIP mechanism in numerical models (Field et al.2017) despite more parameterizations for other SIP processes becoming available (e.g., Hoarau et al.2018; Sullivan et al.2018; Zhao et al.2021). The enhanced ice number concentration can also be caused by the fragmentation of large supercooled liquid droplets (e.g., Evans and Hutchinson1963; Scott and Hobbs1977; Wildeman et al.2017). It has been found that the secondary ice production efficiency is positively correlated with the size of liquid droplets (Lauber et al.2018) and is enhanced in moist environments (Keinert et al.2020). At temperatures higher than −3C the fragmentation of drizzle is still active, as shown by field observations (Lauber et al.2021). In addition, studies using optical sensors mounted on aircrafts have reported the high concentration of ice columns within the temperature range of −10 to −3C (e.g., Koenig1963; Hobbs and Rangno1990; Rangno and Hobbs2001). Based on aircraft measurements from two field campaigns, Korolev et al. (2020) concluded that the secondary ice process is highly associated with the presence of liquid droplets and aged rimed ice in turbulent regions. Recently, Yang et al. (2020) found that the ice concentration in tropical maritime stratiform clouds characterized by the top temperature above −8C is on the order of 10−1101 L−1, which cannot be fully explained by primary ice nucleation, H–M processes or droplet freezing. However, despite the advantage of offering a direct way of interpreting ice microphysics, aircraft observations are only available from a few measurement campaigns and do not seem sufficient for a long-term assessment.

The polarimetric variables, such as differential reflectivity (Zdr), specific differential phase measurements (Kdp) and linear depolarization ratio (LDR), observed by dual-polarization radars are sensitive to the shape of hydrometeors (Bringi and Chandrasekar2001) and allow the analysis of ice particles with specific habits (e.g., Matrosov et al.2001; Hogan et al.2002; Tyynelä and Chandrasekar2014; Li et al.2018). At temperatures of −10 to −2C, the depositional growth of an ice crystal is stronger at the basal faces than at the prism faces. Hence, the formation of columnar ice is preferred (Lamb and Verlinde2011). This distinct habit can produce high Zdr and Kdp as observed by dual-polarization weather radars (Hogan et al.2002; Giangrande et al.2016; Sinclair et al.2016). For vertically pointing Ka- and W-band radars, ice columns usually produce LDR values as high as −15dB, which is distinctively higher than that of most other ice particle types (Aydin and Walsh1999; Tyynelä et al.2011). Oue et al. (2015), Li and Moisseev (2020), Luke et al. (2021), and Li et al. (2021) have shown that this strong LDR signal at the slow falling part of the radar Doppler spectrum can be used to identify columnar-ice crystals. In this study, this method is applied to long-term radar Doppler spectra observations for characterizing the production of columnar-ice particles in stratiform clouds. Similarly to Luke et al. (2021), we show that this phenomenon is not uncommon. By comparing radar-based retrievals of ice number concentrations to INP parameterizations, one of which was derived from observations collected at our measurement site (Schneider et al.2020), we show that the ice number concentrations tend to be 1–3 orders of magnitude higher than the expected INP concentration. This also supports the conclusions reached by Luke et al. (2021).

The paper is organized as follows. Section 2 introduces the data used in this study. The method for identifying columnar-ice particles from radar Doppler spectra is illustrated in Sect. 3. Statistical results are presented in Sect. 4. Section 5 compares the concentrations of columnar-ice particles and INPs in single-layer shallow clouds. Conclusions are given in Sect. 6.

2 Data

2.1 Cloud radar observations

The measurements used in this study were collected at Station for Measuring Ecosystem – Atmosphere Relations II (Hari and Kulmala2005; Petäjä et al.2016) located in Hyytiälä, southern Finland (61.845 N, 24.287 E; 150 m above mean sea level, a.m.s.l.). Since November 2017, a 94 GHz dual-polarization frequency-modulated continuous-wave Doppler cloud radar (Küchler et al.2017) (HYytiälä Doppler RAdar, HYDRA-W) has been operating at the station. The radar is pointing vertically and measures radar signal spectral moments, LDR and dual-polarization Doppler spectra; see Li and Moisseev (2020) for the example of the data. The LDR decoupling is about 30 dB, so the minimum observable LDR is about −30 dB. The radar operates using three chirps that define range resolution, Doppler unambiguous velocity and spectral resolution at three range intervals. Between 102 and 996 m, the range resolution is 25.5 m, Doppler unambiguous velocity is 10.24 m s−1, and the Doppler spectral resolution is 0.02 m s−1. Between 996 and 3577 m, these values are 25.5 m, 5.12 m s−1 and 0.02 m s−1, respectively. For ranges above 3577 m, they are 34 m, 5.12 m s−1 and 0.02 m s−1, respectively. In this study, HYDRA-W observations recorded between February 2018 and April 2020 were utilized (Moisseev2020).

To remove noise from Doppler spectra observations, the spectral lines with the signal-to-noise ratio less than 5 dB were filtered out. Since both co-polar and cross-polar observations were used, i.e., to compute spectral LDR, this filtering could result in complete removal of the cross-polar signal, the power of which is typically 15–30 dB lower than that of the co-polar signal (Bringi and Chandrasekar2001; Moisseev et al.2002). In such cases, no LDR values were computed.

The antenna diameter of HYDRA-W is 0.5 m, which translates into the Fraunhofer far-field distance of 157 m (Sekelsky2002; Falconi et al.2018) and antenna beam width of 0.56. Therefore, the lowest radar range bin, which is not affected by the near-field effect, is 179 m (the fourth range bin). Data recorded at distance in the radar far-field were used in this study, therefore limiting the lowest altitude to 179 m where observations were taken.

In addition to the active remote sensing system, HYDRA-W is capable of estimating the liquid water path (LWP) by using the 89 GHz passive microwave channel observations. The brightness temperature at this band is regularly calibrated using liquid nitrogen. The site-specific relation between the measured 89 GHz brightness temperature and LWP was derived by the radar manufacturer using radio-sounding and reanalysis data.

2.2 Model temperature and humidity profiles

To obtain information on atmospheric state during the cloud observations, forecasts of the Deutscher Wetterdienst (DWD) operational global Icosahedral Nonhydrostatic (ICON) model (Zängl et al.2015) were used. The microphysics scheme in ICON is inherited from the COSMO model (Seifert2008). The ICON model output is provided over all Aerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS) cloud profiling stations (available at, last access: 30 July 2020). The model data have hourly temporal resolution with a horizontal resolution of 13 km (Prill et al.2019; Reinert et al.2020). The height resolution decreases with the increase in altitude; for example, the height resolution is 0.16 km at an altitude of 1 km. Its atmospheric products, such as temperature, relative humidity (RH) and pressure, over Hyytiälä were interpolated into the temporal and spatial resolutions of HYDRA-W (CLU2021).

3 Methods

The mean Doppler velocity (MDV) of hydrometeors observed by a vertically pointing radar is the combination of particle terminal fall velocities and the vertical component of air motion. While Doppler velocity alone could be used to identify certain types of particles (Mosimann1995; Kneifel and Moisseev2020), there are associated limitations. These limitations include uncertainties in hydrometeor classification due to similarities of terminal fall velocities of different particles (Locatelli and Hobbs1974; Barthazy and Schefold2006; Li et al.2020), the presence of a mixture of ice particle populations within the radar volume (Zawadzki et al.2001; Kalesse et al.2019; Li and Moisseev2020) and impact of air motion on the observed MDV (Protat and Williams2011). By using radar Doppler spectra instead of MDV, contributions from different particle populations can be separated (e.g., Zawadzki et al.2001; Kalesse et al.2019; Radenz et al.2019; Li and Moisseev2020; Luke et al.2021). In radar Doppler spectral power observations, the presence of multiple populations of hydrometeors, such as the co-existence of supercooled liquid water and ice (e.g., Zawadzki et al.2001; Shupe et al.2004; Luke et al.2010; Kalesse et al.2016; Li and Moisseev2019), a mixture of different ice types (e.g., Zawadzki et al.2001; Kalesse et al.2019; Radenz et al.2019; Li and Moisseev2020), could manifest as multiple spectral peaks. Even in such cases, however, classification of these particles can be ambiguous.

To further improve the hydrometeor identification, dual-polarization radar observations can be used. For slant measurements the spectral differential reflectivity has been found to be useful (Spek et al.2008). Because hydrometeors typically do not have a preferred azimuth orientation, the differential reflectivity is not very useful for the classification purposes at vertical incidence. Nonetheless, LDR can be used to identify prolate particles (Myagkov et al.2016a, b), such as ice columns (Oue et al.2015). At an elevation angle of 90, columnar-ice particles can produce LDR signals as high as −16 to −13 dB (Oue et al.2015), while other ice particles may produce values smaller than −20 dB (e.g., Tyynelä et al.2011). Furthermore, given the relatively small fall velocities of newly produced columns, in regions where they coexist with other ice particles they usually populate the slow-falling part of the Doppler spectrum. Therefore, the slow-falling part characterized by high spectral LDR (-15 dB) in the Doppler spectrum indicates the presence of columnar-ice particles (Oue et al.2015; Radenz et al.2019).

Figure 1HYDRA-W Doppler spectral (a) power and (b) LDR at 2 February 2018 07:56:38 UTC. (c) Spectral power (black line) and LDR (red crosses) at 1.124 km as marked by the black dot dashed lines in (a) and (b) and LDR (green dots) at 0.51 km as marked by the green dot dashed lines in (a) and (b). The gray, red and blue shaded areas in (c) indicate background ice falling from above, newly formed ice columns and supercooled liquid, respectively, at 1.124 km. Negative velocity indicates downward motion.


Figure 1 shows an example of observed spectral power and LDR. Two distinct populations of ice particles can be clearly identified from Fig. 1a, and the slower-falling one corresponds to ice columns as indicated by the high spectral LDR (Fig. 1b). The observed spectral power and LDR at 1.124 km (black dot dashed lines) are shown in Fig. 1c. Three distinct peaks can be identified from the spectral power, and the slower-falling ice columns are well characterized by the spectral LDR exceeding −18 dB. In contrast, the spectral LDR of faster-falling ice is around −25 dB, which mainly depends on the cross-coupling between the polarization channels (Moisseev et al.2002) and can be much higher than the LDR signal of larger aggregates (Tyynelä et al.2011). Interestingly, supercooled liquid water also seems present, as indicated by the well-defined spectral peak at around 0 m s−1 (Zawadzki et al.2001; Shupe et al.2004; Luke et al.2010; Kalesse et al.2016; Li and Moisseev2019). It appears that this liquid layer extends from ∼0.9 to ∼1.3 km (Fig. 1a). The potential mechanisms of producing these ice columns will be discussed in more detail in the following sections.

Given the spectral characteristics of ice columns as discussed above, the following criteria were set to identify ice columns in clouds.

  • Within the slowest 1 m s−1 of the Doppler spectrum, at least three spectral bins exceed the LDR of −18 dB.

  • The temperature of the radar range bin is between −10 and 0 C.

The observed radar Doppler spectrum is not only dependent on the scattering properties of hydrometeors in the radar volume, but is also affected by the turbulent broadening (Kollias et al.2011; Tridon and Battaglia2015). For example, the air at around 0.51 km seems rather turbulent, as indicated by the spectral power (green dashed line in Fig. 1a). However, it appears that this issue does not significantly affect the results of columnar-ice detection. The noisy spectral LDR values (green dots in Fig. 1c) between 0.3 and 1 m s−1 are attributed to the low signal-to-noise ratio. Such weak impact on spectral LDR due to turbulence may be explained by the distinctively high LDR values of ice columns, which contrast with much weaker LDR signals of ice aggregates (Tyynelä et al.2011).

4 Results

By utilizing HYDRA-W Doppler spectral observations recorded between February 2018 and April 2020, statistics of environmental conditions associated with columnar-ice production were derived. All detected cloud cases within the temperature range of −10 to 0 C were analyzed. From the selected events the cases where significant inversion was detected, which could cause melting (e.g., Kumjian et al.2020), were excluded. Given the data selection criteria, no rainfall or summer cloud cases were analyzed. This was done to avoid potential problems due to radar signal attenuation in rain and melting layer (Li and Moisseev2019). In total, 175 d of observations satisfying the data selection criteria were identified and analyzed.

4.1 Temperature and RH conditions in columnar-ice-producing regions

Formation and growth of ice particles require favorable environmental conditions. These conditions were assessed by using ICON forecasts, which supplemented the radar observations. Here, we define Hcolumn_top as the highest level where ice columns are detected and Tcolumn_top as the temperature at this height. As shown in Fig. 2a, Tcolumn_top values are mostly between −8 and −3C, with the highest frequency at around −5C and a median value of −4.7C. These values are within the growth region of ice columns (Lamb and Verlinde2011). Such temperature distribution also bears a good resemblance to the results obtained from the early rime-splintering laboratory experiment (Hallett and Mossop1974) and a recent statistical study (Luke et al.2021). The statistics of humidity relative to ice (RHice) and water (RHliquid) at Hcolumn_top are shown in Fig. 2b. The median values of RHice and RHliquid are 102.6 % and 98.3 %, respectively, indicating the supply of water vapor is sufficient for growth of ice particles. This finding is, however, not surprising, since the method detects ice columns, and they are growing in this temperature regime. However, the values of RHliquid and RHice should be interpreted with caution. ICON applies a liquid saturation adjustment, limiting the liquid supersaturation to saturation. RHliquid values exceeding 100 % are attributed to numerical artifacts. RHice was calculated based on the forecasted temperature, pressure as well as RHliquid and therefore can be affected by numerical artifacts as well. Given the uncertainty of ICON forecasts, we regard the presented statistics in Fig. 2 as a “sanity” check for our method.

Figure 2Statistics of (a) temperature (Tcolumn_top) and (b) relative humidity (RH) at Hcolumn_top for all identified columnar-ice-producing cases. Hcolumn_top: the highest level where ice columns are detected.


It should be noted that although ice columns can be detected by our method, Hcolumn_top may be lower than the height where ice crystals are generated. There are two potential reasons for this. Firstly, the newly formed ice particles may be less non-spherical (Korolev et al.2020; Luke et al.2021), and in this case they will have LDR values which are much smaller than our detection threshold. Secondly, at early stages of growth the radar signal of ice crystals is rather weak and does not allow accurate detection and identification of columns (Luke et al.2021). However, the altitude difference between Hcolumn_top and the actual height where columns are generated is expected to be small and not significantly affect our results.

4.2 Properties of columnar-ice-producing clouds

There are a number of questions that are associated with formation of ice crystals in clouds at these relatively warm temperatures. Above −10C, the number concentration of INPs is expected to be small and rather uncertain (DeMott et al.2010; Kanji et al.2017). Therefore, it is important to know how often and under which conditions these ice crystals form. Because ice formation could be facilitated by ice particles falling from upper cloud layers, i.e., by SIP, the location where ice columns are forming, with respect to the cloud top, should be identified. Finally, the importance of the columnar-ice production for surface precipitation should also be assessed.

Figure 3Duration of cloud observations (bars) and occurrence of columnar-ice-producing clouds (red dotted curve) over Hyytiälä as a function of cloud top temperature (CTT). The results were calculated based on the data collected from February 2018 to April 2020.


To identify such conditions, cloud top temperature (CTT), defined as the temperature of the highest detected radar return for a given measurement time, is used. Because there are cases where several cloud layers are observed and there are gaps between these layers, typically the top of the lowest one is used. However, particles forming in the upper clouds, while not detected by the radar, may seed lower cloud layers and therefore modify their properties (Vassel et al.2019). To limit the impact of such conditions on our analysis, following Seifert et al. (2009) we have used a radar echo gap of 2 km as the threshold, which defines whether the layers are connected. Recently, Proske et al. (2020) suggested that the threshold of 2 km may overestimate the cases of cloud seeding. For this reason, we have also tried the threshold of 0.5 km to determine the cloud top, but we did not see significant changes in the results.

The statistics of the recorded cloud top temperatures are presented in Fig. 3. The figure (left) shows durations of detected cloud samples within the temperature of −10 to 0 C for a given CTT range as recorded during the observation period. Because of the focus on cold cloud cases, where the temperature in an atmospheric column does not exceed 0 C, the observed cloud cases were recorded between October and April. The observations show that low-level clouds, i.e., clouds with warmer CTT, are relatively more frequent. This resembles the cloud occurrence statistics of, e.g., Hagihara et al. (2010) and Shupe et al. (2011). It appears that deeper clouds, i.e., where CTT is below −12C, are more conducive to columnar-ice production. In these cases the frequency of columnar-ice occurrence is about 25 %–33 %. For warmer clouds the frequency is lower and is around 5 %–13 %. The average occurrence is 15 %. Interestingly, our results are comparable with a recent study by Luke et al. (2021), who found that the occurrence of columnar ice over an Arctic site is between 10 % and 25 % depending on the temperature.

Figure 4Comparison of liquid water path (LWP) for clouds with and without columnar-ice production. A cloud sample was identified if the cloud base was within the temperature range of −10 to 0 C. The box plots represent the median (horizontal strip) and 5 %–95 % quantile (whisker) ranges of the distribution.


As shown in Fig. 2b, the majority of columnar-ice production cases took place in areas around liquid saturation. Although direct observations of liquid were not available, the measurements collected by the 89 GHz passive channel in HYDRA-W allows estimation of LWP. The LWP values for the cloud cases are shown in Fig. 4. The observations show a significant amount of supercooled liquid water present in the atmospheric column. The cloud cases where ice columns were detected tend to have larger LWP values, especially where CTT values were smaller than −8C. This potentially indicates that supercooled liquid droplets may play important roles in formation of ice columns. Given the mixed-phase cloud conditions, the observed columns are most probably ice needles.

Figure 5The same as Fig. 4 but for W-band radar reflectivity at the fourth range bin (179 m above the surface).


Formation of ice crystals is an efficient precipitation process (Lamb and Verlinde2011). To evaluate the impact of columnar-ice production on precipitation, the radar equivalent reflectivity factor is used as the proxy for the precipitation intensity. For clouds where the radar echo extends to the ground, reflectivity values at the fourth radar gate, 179 m above ground level (a.g.l.), were used. As shown in Fig. 5, the reflectivity increases with decreasing CTT. This is due to the link between the cloud thickness and precipitation intensity. The columnar-ice production tends to increase the precipitation intensity. This effect is more pronounced for warmer clouds, where CTT is −12C or warmer. In warmer clouds the precipitation intensity can be enhanced by as much as 10-fold. The factor of 10 increase in precipitation rate appears from the 10 dB increase in reflectivity (Falconi et al.2018). As will be discussed in the next section, warmer clouds tend to be single-layer clouds, where the crystal formation is more directly linked to precipitation formation. Colder clouds are prone to consist of the multiple cloud layers where precipitation processes are affected by multiple processes, such as riming, aggregation, and sublimation, at various levels (e.g., Houze and Medina2005; Verlinde et al.2013; Moisseev et al.2015)

4.3 Columnar-ice production in single-layer and multilayered clouds

For all detected columnar-ice-producing cases, the distribution of CTT was analyzed. As shown in Fig. 6a, ice columns can form in clouds with a wide range of CTT. The majority of cases fall in the CTT range of −20 to 0 C, with two peaks at around −15 and −5C, respectively. The peak at around −15 agrees rather well with the high occurrence of ice columns in clouds, with a CTT of −20 to −12C (Fig. 3).

Figure 6Relative occurrence of (a) CTT and (b) ΔT for columnar-ice production cases. ΔT = Tcolumn_top−CTT.


Because processes responsible for the formation of ice particles in single-layer and multilayered clouds may be different, the classification of the cloud cases was performed. Using CTT and Tcolumn_top, we define ΔT as the temperature difference between them. The larger ΔT is, the lower the inside of the observed cloud system where the columns are formed. The relative occurrence of ΔT also shows two peaks as presented in Fig. 6b. Specifically, one peak is close to ΔT=0K, indicating that ice columns are generated close to the cloud top. The second ΔT peak is around 10 K.

Given the distinct distribution of ΔT, we have grouped the recorded clouds into the following three categories.

  • Type 1: ΔT 2 K – columnar-ice production at cloud top

  • Type 2: 2 K<ΔT 12 K – multilayered cloud

  • Type 3: 12 K<ΔT – multilayered cloud

Representative events of the above cloud types are presented below.

4.3.1 Columnar-ice production at cloud top: ΔT 2 K

This type of cloud is usually single layer, and ice columns are generated close to the cloud top. Figure 7 presents such an event on 4 November 2019. The precipitation intensity is relatively light, with the CTT at around −3C. The W-band reflectivity close to the surface increases to around 0 dB between 14:00 and 16:00 UTC. This region coincides with the enhanced LDR observations, which reaches values as high as −15 dB. Such high LDR values indicate that the dominant ice particle type during this period is columns. As shown in Fig. 7b, the cloud top is turbulent and seems to be capped by an inversion layer (Fig. 7c). The observed LWP ranges between 80 and 150 g m−2. This cloud with relatively low reflectivity persisted over Hyytiälä for about 1 d (not shown). Given the warm cloud top, the primary ice nucleation may not fully explain the significant columnar-ice production (DeMott et al.2010), as will be discussed in more detail later. Regarding the SIP, the H–M process does not seem to be active since it requires falling ice particles serving as rimers to produce ice splinters (Hallett and Mossop1974).

Figure 7The columnar-ice-producing event on 4 November 2019. HYDRA-W observations of (a) equivalent reflectivity factor, (b) mean Doppler velocity, where negative velocity indicates downward motion, and (c) LDR. Panel (d) presents the (left axis) detected columnar-ice region and (right axis) LWP observed by HYDRA-W. The lines in (c) are isotherms produced by ICON.


As shown in Fig. 6, around 22 % of columnar-ice production cases are attributed to single-layer shallow clouds. Bühl et al. (2016) also observed the prevalence of high LDR values for mixed-phase clouds, with a CTT of −5C. They speculated that these particles are formed mainly by primary ice nucleation instead of the SIP. Recently, Yang et al. (2020) reported a similar event of shallow stratiform clouds over the tropical ocean and found that neither INPs nor known SIP mechanisms can fully explain the strong production of ice particles in such clouds, with top temperatures warmer than −8C. In this study, we find that such clouds also frequently occur over Hyytiälä, and more detailed analysis will be presented in Sect. 5.

4.3.2 Columnar-ice production in multilayered clouds: 2 K<ΔT 12 K

The event that took place on 13 February 2018 is representative of the second cloud type, as defined by ΔT. As shown in Fig. 8a, the precipitation intensity during this event is higher than during the discussed single-layer shallow cloud case. The cloud top temperature of the upper cloud layer is about −15 to −12C. Before 08:00 UTC, the observed LWP is close to 100 g m−2, and the mean Doppler velocity at around 0.8 km is relatively small ( 1 m s−1), which indicates that particles are unrimed or very lightly rimed. From 08:00 to 09:30 UTC, the falling snowflakes seem to be heavily rimed, as revealed by a rather high LWP (from 200 to over 400 g m−2) and mean Doppler velocity measurements (1.2–1.8 m s−1) (Kneifel and Moisseev2020). The high LWP period coincides with the region of ice columns (Fig. 8d). During this period, the observed LDR values are enhanced but still relatively small. This is due to masking of the needle LDR signal by larger snowflakes.

Figure 8The columnar-ice-producing event recorded on 13 February 2018. HYDRA-W observations of (a) equivalent reflectivity factor, (b) mean Doppler velocity, where negative velocity indicates downward motion, and (c) LDR. Panel (d) presents the (left axis) detected columnar-ice region and (right axis) LWP observed by HYDRA-W. The lines in (a) and (c) are isotherms produced by ICON. Note that the y-axis scale in (a) is different from that in (b), (c) and (d).


This type of cloud frequently occurs over Hyytiälä (Figs. 3 and 6). Although sounding measurements were absent, this cloud type seems to very similar to the one reported by Westbrook and Illingworth (2013), namely, a layer of supercooled liquid water with the top temperature of around −15C seeding low-level stratus clouds in the boundary layer. In this event, the presence of supercooled liquid water may not be directly determined; however, the enhanced LWP values are indicative of the vigorously supercooled liquid water generation. In addition, the falling ice particles between 08:00 UTC and 09:30 UTC seem to be heavily rimed, as evident from mean Doppler velocity (Kneifel and Moisseev2020). The combination of the presence of supercooled liquid water and riming indicates that the H–M process could be taking place and could be responsible for the columnar-ice production.

4.3.3 Deep multilayered clouds: ΔT> 12 K

The third cloud type is not very different from the second one and represents the tail of the observed ΔT distribution as shown in Fig. 6. This type of cloud system is a deeper precipitating system with a CTT of −60 to −40C; see Fig. 9 for an example. The presented case took place on 2 February 2018. There are several features that are worthwhile pointing out. The mean Doppler velocity observations exhibit signatures of atmospheric waves. Between 05:00 and 06:00 UTC such waves can be observed at around 1 km altitude. At the later time, the wave signatures appear at 0.5 km a.g.l. The strongest velocity variation, observed around 06:00 UTC, seems to coincide with the LWP peak. At the same time as the waves are observed, signatures of columnar-ice production are also detected, pointing to a possible connection between the two.

Figure 9Same as Fig. 8 but for the columnar-ice-producing event observed on 2 February 2018. A zoom-in view on the wave signatures between 05:30 and 07:00 UTC is presented in (b).


Deep precipitating clouds usually have a large number of ice crystals formed at the cloud top. Given the large ice flux, manifested by the higher radar reflectivity values, in this precipitation system, it is difficult for supercooled liquid droplets to survive. The supercooled water droplets can be rapidly depleted through the Wegener–-Bergeron–-Findeisen (WBF) process (Lamb and Verlinde2011; Korolev et al.2017) as well as riming (Fukuta and Takahashi1999). Nevertheless, the atmospheric waves could generate conditions needed for forming and maintaining the presence of supercooled liquid water droplets (Korolev1995; Korolev and Field2008; Majewski and French2020). Recently, Li et al. (2021) provided radar observational evidence showing that the Kelvin–Helmholtz instability is favorable for supercooled drizzle and secondary ice production. In such cases, ice needles may be generated by the H–M process (Hogan et al.2002; Houser and Bluestein2011) or freezing breakup (Luke et al.2021; Li et al.2021).

4.4 Characteristics of different columnar-ice-producing cloud types

As discussed above, we have identified three types of cloud systems where columnar-ice particles form. To better understand how columnar-ice production is related to cloud properties, persistence of columnar ice in these clouds and the amount of LWP were considered for further analysis.

4.4.1 Persistence of columnar-ice crystals

As was demonstrated by the case studies, the columnar-ice production may persist over several hours, and therefore these particles may play a major role in determining cloud properties. To document this, we have derived statistics of the columnar-ice production persistence for each cloud case. This was done by computing the duration of a continuous columnar-ice production event. Since in some cases, for example, in the presented shallow single-layer clouds (Fig. 7) formation of columns may be intermittent in nature (as can also be seen in Luke et al.2021), similar to Shupe et al. (2011), gaps of less than 30 min were accepted. In addition, cases persisting for less than 1 min were removed. As shown in Fig. 10, 40 %–50 % of columnar-ice-producing events persist for less than 1 h. However, there is a significant fraction that could last for more than 3 or even 6 h. This hints that the production of ice columns plays an important role in defining cloud properties and should be included while considering radiative or precipitation properties of such clouds (see, for example, Fig. 9).

Figure 10Relative occurrence of the persistence of columnar-ice-producing events from February 2018 to April 2020. ΔT=Tcolumn_top-CTT.


4.4.2 LWP

Presence of liquid water droplets may be a necessary condition for the formation of the observed ice columns. Compared to clouds where no ice columns were detected, columnar-ice-producing clouds have higher LWP values (Fig. 4). For different columnar-ice-producing cloud types, the observed LWP values seem to be somewhat different. Figure 11 shows LWP occurrence for these three cloud types. In general, their distributions are similar, while identifiable differences still exist. The median value of LWP for single-layer clouds is the lowest, and the second type of cloud has the highest LWP. Interestingly, while both the second and third types of clouds are multilayered, the LWP of the second type is detectably higher than for the third one. Comparing Figs. 8 and 9, we find that the radar reflectivity above the columnar-ice layer is generally higher in the third cloud type than the second one. The INP concentrations in clouds with a cold top are expected to be higher than for warmer clouds, and hence more falling ice particles are expected for deep precipitating clouds. Therefore, we speculate that one of the reasons for the difference in LWP may be the ice number concentration, which is related to the consumption of the liquid water via the WBF process (Lamb and Verlinde2011) and riming (Fukuta and Takahashi1999).

Figure 11Relative occurrence of LWP for all columnar-ice-producing events. LWP: liquid water path. ΔT=Tcolumn_top-CTT.


5 Potential role of SIP in columnar-ice production in single-layer shallow clouds

Given the rather frequent formation of ice crystals at temperatures warmer than −10C, where expected INP concentrations are low, it is important to investigate a potential role of SIP. In multilayered clouds, identified here as cloud types II and III, it has been found that the ice formation can be enhanced by the H–M process (e.g., Grazioli et al.2015; Giangrande et al.2016; Sinclair et al.2016; Keppas et al.2017; Sullivan et al.2018; Gehring et al.2020), among other mechanisms (Korolev and Leisner2020). In single-layer clouds the mechanisms for ice multiplication are less established. However, recent studies by Lauber et al. (2018) and Keinert et al. (2020) have shown that freezing fragmentation may play such a role. To study whether these columnar-ice particles can be attributed to the SIP, the estimated ice number concentrations are compared to the expected INP concentrations. If the derived ice number concentrations exceed these of INPs, we can conclude that the SIP is potentially active in such clouds. The concentration of INPs was computed by using CTT in the following different parameterizations.

  • Fletcher (1962) parameterization based on INP measurements obtained below −10C.

  • Cooper (1986) parameterization which is not directly derived from INP measurements but the observed ice number concentrations when the impact of SIP is minimized. The temperature of measurements is between −30 and −5C.

  • DeMott et al. (2010) parameterization based on INP measurements from nine sites between −35 and −9C. In our study we have used the average INP concentration–temperature relation presented in DeMott et al. (2010).

  • Schneider et al. (2020) parameterization derived from the INP measurements obtained at Hyytiälä. The temperature range is −20 to −8C. Also from this study, the average INP concentration–temperature relation was used.

As was previously discussed, the INP parameterizations differ significantly at these temperatures (−10 to 0 C), as shown in Fig. 12. It should be noted that not all parameterizations were derived using observations at these cloud temperatures, and some of them were extrapolated beyond their validity range. The most interesting comparison is to Schneider et al. (2020), which is based on observations at Hyytiälä collected during 2018, so their observation period at least partially overlaps with ours. It should also be pointed out that INP observations were carried out at the ground, where INP concentrations are typically higher.

Radar-based retrieval of particle number concentration is rather uncertain. Because of this uncertainty, the derived number concentrations should be treated as our best estimates of the order of magnitude of the ice column number concentration. As shown in Fig. 4, the observed LWP for columnar-ice-producing clouds is significantly higher than those without ice columns. Hence, pristine ice crystals are anticipated to grow in mixed-phase conditions, and ice needles rather than solid columns are expected to form (Lamb and Verlinde2011). This limits the parameter space, where we need to search for microphysical properties of ice particles to constrain our retrieval. The retrieval is based on estimating ice water content from radar observations, following Hogan et al. (2006), as

(1) IWC = 10 0.00058 Z T + 0.0923 Z - 0.00706 T - 0.992 ( g m - 3 ) ,

where Z is the W-band radar reflectivity and T is the air temperature. Because in a single-layer shallow cloud (ΔT≤2C) ice needles are the predominant precipitation particles, radar reflectivity measured close to the ground, in the fourth radar gate, 329 m a.m.s.l. or 179 m a.g.l., is used in this retrieval. The selection of the reflectivity measured close to the ground helps to limit potential attenuation problems as well. Using IWC, the number concentration of ice needles, Nneedle, can be estimated as follows (Li et al.2021):

(2) N needle = IWC 10 3 m needle ( L - 1 ) ,

where mneedle is the mass of a characteristic ice needle. The introduced uncertainty at this step depends on the definition of the characteristic needle. Here, we use mean Doppler velocity (MDV) and the velocity–mass relation to estimate mneedle. Since MDV is reflectivity weighted, mneedle would be mainly determined by larger ice particles, and therefore the resulting Nneedle is underestimated. For the purpose of this study, this underestimation is not a major issue, because we want to test whether the observed Nneedle exceeds expected INP concentrations.

There are a number of reported ice needle properties. To take into account potential differences in ice needle properties, two relations linking velocity and mass by Kajikawa (1976) for rimed needles and Heymsfield (1972) for unrimed needles were used. For rimed ice needles, the relation between terminal fall velocity and mass was derived by Kajikawa (1976) and can be written as

(3) v needle, rimed = 1.55 ( 10 3 m needle, rimed ) 0.271 ρ ( z ) ρ ( 1024 m ) 0.4 ( m s - 1 ) ,

where mneedle, rimed is the mass of rimed ice needles, and the term ρ(z)ρ(1024m)0.4 accounts for the change in air density ρ at a given height z (Heymsfield et al.2007). In our study z is 329 m, and in Kajikawa (1976) the altitude where needles were observed is 1024 m. For unrimed needles, Heymsfield (1972) derived a relation linking terminal fall velocity and needle length, L. The terminal velocity of unrimed needles at a given height can be estimated as follows (Heymsfield1972):

(4) v needle, unrimed = ( 0.0006 + 0.2796 L - 0.0497 L 2 + 0.0041 L 3 ) ρ ( z ) ρ ( 0 m ) 0.4 ( m s - 1 ) .

The needle can be modeled as a cylinder, the mass of which is

(5) m needle, unrimed = 10 - 3 π 4 ρ ice L D 2 ( g ) ,

where ρice and D denote density of unrimed needles and the minor axis, respectively. Their parameterizations have been given by Heymsfield (1972):

(6) ρ ice = 0.6 L - 0.117 ( g cm - 3 )


(7) D = 0.1973 L 0.414 ( mm ) .

Applying the power law fit to mneedle, unrimed and vneedle, unrimed values when L[0.03,5] (mm) yields

(8) v needle, unrimed = 1.09 ( 10 3 m needle, unrimed ) 0.377 ρ ( z ) ρ ( 0 m ) 0.4 ( m s - 1 ) .

The observed MDV is affected by vertical air motion. To at least partially mitigate this issue, the observed MDV is averaged over 20 min (Protat and Williams2011; Mosimann1995; Kneifel and Moisseev2020; Silber et al.2020). While this step reduces the impact of air motion by averaging Doppler velocity over updrafts and downdrafts, the residual air motion is expected to widen the retrieved distribution of Nneedle.

The derived number concentrations of ice particles are compared to expected INP concentrations (Fig. 12). The results show that the estimated ice number concentrations for rimed needles (Kajikawa1976) are generally larger than that of unrimed needles (Heymsfield1972). Regardless of the difference between rimed and unrimed needles and the INP parameterization used, there seems to be a large fraction of cases where INP concentrations are not sufficient to explain observed Nneedle. Our results are similar to the conclusion reached by Luke et al. (2021), who used a different approach for establishing the range of ice crystal concentration from radar observations. As shown in Fig. 12, the majority of Nneedle values fall in the range of 10−2101L−1, which is similar to aircraft measurements obtained in tropical stratiform clouds (Yang et al.2020).

Figure 12Probability density function (PDF) of estimated columnar ice (left) and INP (right) concentrations in shallow single-layer clouds. INPs: ice-nucleating particles.


The significant discrepancy between INP concentrations estimated from INP parameterizations and retrieved ice number concentrations indicates that primary ice nucleation does not seem to be the only mechanism responsible for the formation of ice particles in these shallow clouds. Because the analysis was performed on shallow single-layer clouds, this discrepancy may not be explained by the H–M process (Hallett and Mossop1974) since no rimers are falling from upper clouds. So it appears that other, less studied, SIP mechanisms may play an important role in amplifying ice number concentrations in such shallow clouds. This conclusion is in line with a number of other studies. For example, Knight (2012) found that the SIP may take place at −5C in the absence of rimers, for which the cause is still under investigation. Recently, a similar finding was reported for stratiform clouds over the tropical ocean by Yang et al. (2020). They have speculated that droplet collisional freezing (Hobbs1965; Alkezweeny1969) and pre-activated INPs (Roberts and Hallett1968; Mossop1970) could be responsible for this discrepancy. Recent laboratory studies (Lauber et al.2018; Keinert et al.2020) have shown that freezing breakup may be a source of secondary ice particles. Luke et al. (2021) have also suggested that freezing breakup may be more efficient than the H–M process in nature.

6 Conclusions

This study documents formation of ice particles in clouds at the temperatures of −10C or warmer. The analysis was performed using W-band cloud radar observations collected at the University of Helsinki Hyytiälä station starting from February 2018 through April 2020. The columnar-ice particles were identified using measurements of spectral LDR. It was found that columnar-ice formation is relatively frequent in clouds at temperatures of −10C or warmer. The occurrence frequency of columnar-ice particles is 5 %–13 % in clouds, with top temperatures exceeding −12C. In colder clouds, this percentage can be as high as 33 %. The columnar-ice-producing clouds tend to have a higher LWP, potentially indicating that supercooled water droplets are important for formation of the observed ice particles. It was also observed that columnar-ice production seems to have a significant impact on the surface precipitation. This effect is especially important for warmer clouds.

Using the temperature difference, ΔT, between the altitudes where columns are first detected and cloud top, the columnar-ice-producing clouds were subdivided into three categories. The first category, where ΔT is less than or equal to 2 K, represents shallow single-layer clouds. In these clouds ice particles are forming at or close to the cloud top. The other two categories, where 2 K<ΔT≤12K and ΔT>12K, represent deeper multilayered clouds. In multilayered cloud systems, columnar-ice crystals are forming in the lower cloud layer seeded by ice particles falling from upper cloud levels. It was observed that 40 %–50 % of columnar-ice production cases persist for 1 h or less, while in some cases they can persist for over 6 h. The distributions of LWP values for the three types of columnar-ice-producing clouds are somewhat different. The median LWP value is largest (221 g m−2) in clouds where 2 K<ΔT≤12K. Such high LWPs could favor riming and cause the Hallet–Mossop process. To draw a definite conclusion, however, a more thorough study, where locations of supercooled liquid layers are identified, is needed. For the single-layer shallow clouds, number concentrations of ice columns were derived from the radar observations. It was observed that the concentration of ice particles exceeds the expected concentration of INPs for a large number of cases. This indicates that a SIP mechanism is active in these clouds. Given that in single-layer shallow clouds there are no rimers that could cause the H–M process, we advocate that another SIP process may play a role here.

Code availability

The code used for processing radar data is available from GitHub (last access: 1 October 2021) and from Zenodo (Li2021).

Data availability

W-band radar data are created by Dmitri Moisseev and are available from (Moisseev2020). ICON data are generated by the European Research Infrastructure for the observation of Aerosol, Clouds and Trace Gases (ACTRIS) and are available from the ACTRIS Data Centre using the following link: (CLU2021).

Author contributions

DM conceptualized the study. HL performed the investigation and wrote the draft. OM, PT and DM contributed to reviewing and editing this paper.

Competing interests

Some of the authors are members of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors have also no other competing interests to declare.


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

Special issue statement

This article is part of the special issue “Ice nucleation in the boreal atmosphere”. It is not associated with a conference.


We would like to thank the personnel of Hyytiälä station for their support in field observation. Matti Leskinen is thanked for his work in radar maintenance. We want to thank Anniina Korpinen, Simo Tukiainen and Axel Seifert for helpful discussions on ICON forecasts.

We thank the European Research Infrastructure for the observation of ACTRIS for providing the DWD ICON (Global) model data, which were produced by the DWD and the Finnish Meteorological Institute and are available for download from (last access: 1 October 2021). ACTRIS has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement no. 654109.

Financial support

This research has been supported by the Horizon 2020 (grant nos. ACTRIS-2654109, ACTRIS PPP 739530, ACTRIS-IMP 871115, ATMO-ACCESS 101008004, and ERA-PLANET iCUPE 689443) and the Academy of Finland (ACTRIS-NF 328616, ACTRIS-CF 329274, NanoBioMass 307537, ACRoBEAR 334792, and Center of Excellence in Atmospheric Sciences, 307331, Atmosphere and Climate Competence Center, 337549), University of Helsinki (ACTRIS-HY).

Open-access funding was provided by the Helsinki University Library.

Review statement

This paper was edited by Paul Zieger and reviewed by two anonymous referees.


Alkezweeny, A.: Freezing of supercooled water droplets due to collision, J. Appl. Meteorol., 8, 994–995, 1969. a

Aydin, K. and Walsh, T. M.: Millimeter wave scattering from spatial and planar bullet rosettes, IEEE T. Geosci. Remote Sens., 37, 1138–1150, 1999. a

Barthazy, E. and Schefold, R.: Fall velocity of snowflakes of different riming degree and crystal types, Atmos. Res., 82, 391–398, 2006. a

Beard, K. V.: Ice initiation in warm-base convective clouds: An assessment of microphysical mechanisms, Atmos. Res., 28, 125–152, 1992. a

Bringi, V. N. and Chandrasekar, V.: Polarimetric Doppler weather radar: principles and applications, Cambridge University Press, 2001. a, b

Bühl, J., Seifert, P., Myagkov, A., and Ansmann, A.: Measuring ice- and liquid-water properties in mixed-phase cloud layers at the Leipzig Cloudnet station, Atmos. Chem. Phys., 16, 10609–10620,, 2016. a

CLU: Icon-iglo-12-23, icon-iglo-36-47, gdas1 model data; 2018-02-01 to 2020-04-30; from Hyytiälä, ACTRIS Data Centre [data set], (last access: 1 October 2021), 2021. a, b

Cooper, W. A.: Ice initiation in natural clouds, Meteorol. Mon., 21, 29–32,, 1986. a

DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D., Twohy, C. H., Richardson, M., Eidhammer, T., and Rogers, D.: 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

DeMott, P. J., Hill, T. C., McCluskey, C. S., Prather, K. A., Collins, D. B., Sullivan, R. C., Ruppel, M. J., Mason, R. H., Irish, V. E., Lee, T., et al.: Sea spray aerosol as a unique source of ice nucleating particles, P. Natl. Acad. Sci. USA, 113, 5797–5803, 2016. a, b, c

Dye, J. and Hobbs, P.: Effect of carbon dioxide on the shattering of freezing water drops, Nature, 209, 464–466, 1966. a

Evans, D. and Hutchinson, W.: The electrification of freezing water droplets and of colliding ice particles, Q. J. Roy. Meteor. Soc., 89, 370–375, 1963. a

Falconi, M. T., von Lerber, A., Ori, D., Marzano, F. S., and Moisseev, D.: Snowfall retrieval at X, Ka and W bands: consistency of backscattering and microphysical properties using BAECC ground-based measurements, Atmos. Meas. Tech., 11, 3059–3079,, 2018. a, b

Field, P., Lawson, P., Brown, G., Lloyd, C., Westbrook, D., Moisseev, A., 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., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Secondary ice production: Current state of the science and recommendations for the future, Meteorol. Monogr., 58, 7–1, 2017. a, b

Findeisen, W. and Findeisen, E.: Untersuchungen uber die Eissplitterbildung an Reifschichten (Ein Beitrag zur Frage der Entstehung der Gewitterelektrizitat und zur Mikrostruktur der Cumulonimben), Meteorol. Z, 60, 145–154, 1943. a

Fletcher, N. H.: The physics of rainclouds, Cambridge University Press, 1962. a

Fridlind, A. M., Ackerman, A., McFarquhar, G., Zhang, G., Poellot, M., DeMott, P., Prenni, A., and Heymsfield, A.: Ice properties of single-layer stratocumulus during the Mixed-Phase Arctic Cloud Experiment: 2. Model results, J. Geophys. Res.-Atmos., 112, D24202,, 2007. a

Fukuta, N. and Takahashi, T.: The Growth of Atmospheric Ice Crystals: A Summary of Findings in Vertical Supercooled Cloud Tunnel Studies, J. Atmos. Sci., 56, 1963–1979,<1963:TGOAIC>2.0.CO;2, 1999. a, b

Gehring, J., Oertel, A., Vignon, É., Jullien, N., Besic, N., and Berne, A.: Microphysics and dynamics of snowfall associated with a warm conveyor belt over Korea, Atmos. Chem. Phys., 20, 7373–7392,, 2020. a

Giangrande, S. E., Toto, T., Bansemer, A., Kumjian, M. R., Mishra, S., and Ryzhkov, A. V.: Insights into riming and aggregation processes as revealed by aircraft, radar, and disdrometer observations for a 27 April 2011 widespread precipitation event, J. Geophys. Res.-Atmos., 121, 5846–5863, 2016. a, b

Grazioli, J., Lloyd, G., Panziera, L., Hoyle, C. R., Connolly, P. J., Henneberger, J., and Berne, A.: Polarimetric radar and in situ observations of riming and snowfall microphysics during CLACE 2014, Atmos. Chem. Phys., 15, 13787–13802,, 2015. a

Hagihara, Y., Okamoto, H., and Yoshida, R.: Development of a combined CloudSat-CALIPSO cloud mask to show global cloud distribution, J. Geophys. Res.-Atmos., 115, D00H33, 2010. a

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

Hari, P. and Kulmala, M.: Station for measuring ecosystem-atmosphere relations: (SMEAR II), Boreal Environ. Res., 10, 315–322, 2005. a

Heymsfield, A.: Ice crystal terminal velocities, J. Atmos. Sci., 29, 1348–1357, 1972. a, b, c, d, e

Heymsfield, A. J., Bansemer, A., and Twohy, C. H.: Refinements to ice particle mass dimensional and terminal velocity relationships for ice clouds. Part I: Temperature dependence, J. Atmos. Sci., 64, 1047–1067, 2007. a

Hoarau, T., Pinty, J.-P., and Barthe, C.: A representation of the collisional ice break-up process in the two-moment microphysics LIMA v1.0 scheme of Meso-NH, Geosci. Model Dev., 11, 4269–4289,, 2018. a

Hobbs, P.: The aggregation of ice particles in clouds and fogs at low temperatures, J. Atmos. Sci., 22, 296–300, 1965. a

Hobbs, P. V. and Rangno, A. L.: Ice particle concentrations in clouds, J. Atmos. Sci., 42, 2523–2549, 1985. a, b

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

Hogan, R., Field, P., Illingworth, A., Cotton, R., and Choularton, T.: Properties of embedded convection in warm-frontal mixed-phase cloud from aircraft and polarimetric radar, Q. J. Roy. Meteor. Soc., 128, 451–476, 2002. a, b, c

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

Houser, J. L. and Bluestein, H. B.: Polarimetric Doppler radar observations of Kelvin–Helmholtz waves in a winter storm, J. Atmos. Sci., 68, 1676–1702, 2011. a

Houze Jr., R. A. and Medina, S.: Turbulence as a mechanism for orographic precipitation enhancement, J. Atmos. Sci., 62, 3599–3623, 2005. a

Kajikawa, M.: Observation of falling motion of columnar snow crystals, J. Meteorol. Soc. Jpn. Ser. II, 54, 276–284, 1976. a, b, c, d

Kalesse, H., Vogl, T., Paduraru, C., and Luke, E.: Development and validation of a supervised machine learning radar Doppler spectra peak-finding algorithm, Atmos. Meas. Tech., 12, 4591–4617,, 2019. a, b, c

Kalesse, H., Szyrmer, W., Kneifel, S., Kollias, P., and Luke, E.: Fingerprints of a riming event on cloud radar Doppler spectra: observations and modeling, Atmos. Chem. Phys., 16, 2997–3012,, 2016. a, b

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, Meteorological Monographs, 58, 1–1, 2017. a, b

Keinert, A., Spannagel, D., Leisner, T., and Kiselev, A.: Secondary Ice Production upon Freezing of Freely Falling Drizzle Droplets, J. Atmos. Sci., 77, 2959–2967,, 2020. a, b, c

Keppas, S. C., Crosier, J., Choularton, T., and Bower, K.: Ice lollies: An ice particle generated in supercooled conveyor belts, Geophys. Res. Lett., 44, 5222–5230, 2017. a

Kneifel, S. and Moisseev, D.: Long-Term Statistics of Riming in Nonconvective Clouds Derived from Ground-Based Doppler Cloud Radar Observations, J. Atmos. Sci., 77, 3495–3508, 2020. a, b, c, d

Knight, C. A.: Ice growth from the vapor at −5C, J. Atmos. Sci., 69, 2031–2040, 2012. a

Koenig, L. R.: The glaciating behavior of small cumulonimbus clouds, J. Atmos. Sci., 20, 29–47, 1963. a

Kollias, P., Rémillard, J., Luke, E., and Szyrmer, W.: Cloud radar Doppler spectra in drizzling stratiform clouds: 1. Forward modeling and remote sensing applications, J. Geophys. Res.-Atmos., 116, D13201,, 2011. a

Korolev, A. and Field, P. R.: The effect of dynamics on mixed-phase clouds: Theoretical considerations, J. Atmos. Sci., 65, 66–86, 2008. a

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

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, 5–1, 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.: The influence of supersaturation fluctuations on droplet size spectra formation, J. Atmos. Sci., 52, 3620–3634, 1995. a

Küchler, N., Kneifel, S., Löhnert, U., Kollias, P., Czekala, H., and Rose, T.: A W-band radar–radiometer system for accurate and continuous monitoring of clouds and precipitation, J. Atmos. Ocean. Tech., 34, 2375–2392, 2017. a

Kumjian, M. R., Tobin, D. M., Oue, M., and Kollias, P.: Microphysical Insights into Ice Pellet Formation Revealed by Fully Polarimetric Ka-Band Doppler Radar, J. Appl. Meteorol. Clim., 59, 1557–1580, 2020. a

Lamb, D. and Verlinde, J.: Physics and chemistry of clouds, Cambridge University Press, 2011. a, b, c, d, e, f

Lauber, A., Kiselev, A., Pander, T., Handmann, P., and Leisner, T.: Secondary ice formation during freezing of levitated droplets, J. Atmos. Sci., 75, 2815–2826, 2018. a, b, c

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

Li, H.: Hyytiala W-band radar spectra data processing tool, Zenodo [code],, 2021. a

Li, H. and Moisseev, D.: Melting layer attenuation at Ka- and W-bands as derived from multi-frequency radar Doppler spectra observations, J. Geophys. Res.-Atmos., 124, 9520–9533,, 2019. a, b, c

Li, H. and Moisseev, D.: Two layers of melting ice particles within a single radar bright band: Interpretation and implications, Geophys. Res. Lett., 47, e2020GL087499,, 2020. a, b, c, d, e

Li, H., Moisseev, D., and von Lerber, A.: How does riming affect dual-polarization radar observations and snowflake shape?, J. Geophys. Res.-Atmos., 123, 6070–6081, 2018. a

Li, H., Tiira, J., von Lerber, A., and Moisseev, D.: Towards the connection between snow microphysics and melting layer: insights from multifrequency and dual-polarization radar observations during BAECC, Atmos. Chem. Phys., 20, 9547–9562,, 2020. a

Li, H., Korolev, A., and Moisseev, D.: Supercooled liquid water and secondary ice production in Kelvin–Helmholtz instability as revealed by radar Doppler spectra observations, Atmos. Chem Phys., 21, 13593–13608,, 2021. a, b, c, d

Locatelli, J. D. and Hobbs, P. V.: Fall speeds and masses of solid precipitation particles, J. Geophys. Re., 79, 2185–2197, 1974. a

Luke, E. P., Kollias, P., and Shupe, M. D.: Detection of supercooled liquid in mixed-phase clouds using radar Doppler spectra, J. Geophys. Res.-Atmos., 115, D19201,, 2010. a, b

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

Majewski, A. and French, J. R.: Supercooled drizzle development in response to semi-coherent vertical velocity fluctuations within an orographic-layer cloud, Atmos. Chem. Phys., 20, 5035–5054,, 2020. a

Matrosov, S. Y., Reinking, R. F., Kropfli, R. A., Martner, B. E., and Bartram, B.: On the use of radar depolarization ratios for estimating shapes of ice hydrometeors in winter clouds, J. Appl. Meteorol., 40, 479–490, 2001. a

McCluskey, C. S., Ovadnevaite, J., Rinaldi, M., Atkinson, J., Belosi, F., Ceburnis, D., Marullo, S., Hill, T. C. J., 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

Moisseev, D.: Hyytiala W-band cloud radar 2017–2019 (Version 1), [data set],, 2020. a, b

Moisseev, D. N., Unal, C. M., Russchenberg, H. W., and Ligthart, L. P.: Improved polarimetric calibration for atmospheric radars, J. Atmos. Ocean. Tech., 19, 1968–1977, 2002. a, b

Moisseev, D. N., Lautaportti, S., Tyynela, J., and Lim, S.: Dual-polarization radar signatures in snowstorms: Role of snowflake aggregation, J. Geophys. Res.-Atmos., 120, 12644–12655, 2015. a

Mosimann, L.: An improved method for determining the degree of snow crystal riming by vertical Doppler radar, Atmos. Res., 37, 305–323, 1995. a, b

Mossop, S.: Concentrations of ice crystals in clouds, B. Am. Meteorol. Soc., 51, 474–480, 1970. a

Mossop, S.: The origin and concentration of ice crystals in clouds, B. Am. Meteorol. Soc., 66, 264–273, 1985. a

Myagkov, A., Seifert, P., Bauer-Pfundstein, M., and Wandinger, U.: Cloud radar with hybrid mode towards estimation of shape and orientation of ice crystals, Atmos. Meas. Tech., 9, 469–489,, 2016a. a

Myagkov, A., Seifert, P., Wandinger, U., Bühl, J., and Engelmann, R.: Relationship between temperature and apparent shape of pristine ice crystals derived from polarimetric cloud radar observations during the ACCEPT campaign, Atmos. Meas. Tech., 9, 3739–3754,, 2016b. a

Niemand, M., Möhler, O., Vogel, B., Vogel, H., Hoose, C., Connolly, P., Klein, H., Bingemer, H., 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

Oue, M., Kumjian, M. R., Lu, Y., Verlinde, J., Aydin, K., and Clothiaux, E. E.: Linear depolarization ratios of columnar ice crystals in a deep precipitating system over the Arctic observed by zenith-pointing Ka-band Doppler radar, J. Appl. Meteorol. Clim., 54, 1060–1068, 2015. a, b, c, d

Petäjä, T., O'Connor, E. J., Moisseev, D., Sinclair, V. A., Manninen, A. J., Väänänen, R., von Lerber, A., Thornton, J. A., Nicoll, K., Petersen, W., Chandrasekar, V., Smith, J. N., Winkler, P. M., Krüger, O., Hakola, H., Timonen, H., Brus, D., Laurila, T., Asmi, E., Riekkola, M. L., Mona, L., Massoli, P., Engelmann, R., Komppula, M., Wang, J., Kuang, C., Bäck, J., Virtanen, A., Levula, J., Ritsche, M., and Hickmon, N.: BAECC: A field campaign to elucidate the impact of biogenic aerosols on clouds and climate, B. Am. Meteorol. Soc., 97, 1909–1928, 2016. a

Petters, M. and Wright, T.: Revisiting ice nucleation from precipitation samples, Geophys. Res. Lett., 42, 8758–8766, 2015. a, b, c

Prill, F., Reinert, D., Rieger, D., Zaengl, G., Schroeter, J., Foerstner, J., Werchner, S., Weimer, M., Ruhnke, R., and Vogel, B.: Working with the ICON Model–Practical exercises for NWP mode and ICON-ART, ICON Model Tutorial, Deutscher Wetterdienst, Offenbach, Germany, 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

Protat, A. and Williams, C. R.: The accuracy of radar estimates of ice terminal fall speed from vertically pointing Doppler radar measurements, J. Appl. Meteorol. Clim., 50, 2120–2138, 2011. a, b

Radenz, M., Bühl, J., Seifert, P., Griesche, H., and Engelmann, R.: peakTree: a framework for structure-preserving radar Doppler spectra analysis, Atmos. Meas. Tech., 12, 4813–4828,, 2019. a, b, c

Rangno, A. L. and Hobbs, P. V.: Ice particles in stratiform clouds in the Arctic and possible mechanisms for the production of high ice concentrations, J. Geophys. Res.-Atmos., 106, 15065–15075, 2001. a, b

Reinert, D., Prill, F., Frank, H., Denhard, M., Baldauf, M., Schraff, C., Gebhardt, C., Marsigli, C., and Zängl, G.: DWD Database Reference for the Global and Regional ICON and ICON-EPS Forecasting System, Tech. rep., Deutscher Wetterdienst, available at: (last access: 1 October 2021), 2020. a

Roberts, P. and Hallett, J.: A laboratory study of the ice nucleating properties of some mineral particulates, Q. J. Roy. Meteor. Soc., 94, 25–34, 1968. a, b

Schneider, J., Höhler, K., Heikkilä, P., Keskinen, J., Bertozzi, B., Bogert, P., Schorr, T., Umo, N. S., Vogel, F., Brasseur, Z., Wu, Y., Hakala, S., Duplissy, J., Moisseev, D., Kulmala, M., Adams, M. P., Murray, B. J., Korhonen, K., Hao, L., Thomson, E. S., Castarède, D., Leisner, T., Petäjä, T., and Möhler, O.: The seasonal cycle of ice-nucleating particles linked to the abundance of biogenic aerosol in boreal forests, Atmos. Chem. Phys., 21, 3899–3918,, 2021. a, b, c, d

Scott, B. C. and Hobbs, P. V.: A theoretical study of the evolution of mixed-phase cumulus clouds, J. Atmos. Sci., 34, 812–826, 1977. a

Seifert, A.: A revised cloud microphysical parameterization for COSMO-LME, COSMO Newsletter, 7, 25–28, 2008. a

Seifert, P., Ansmann, A., Mattis, I., Althausen, D., and Tesche, M.: Lidar-Based Profiling of the Tropospheric Cloud-Ice Distributionto Study the Seeder-Feeder Mechanism and the Role of Saharan Dust as Ice Nuclei, in: Proceedings of the 8th International Symposium on Tropospheric Profiling, edited by: Apituley, A., Russchenberg, H. W. J., and Monna, W. A. A., p. 5, Delft, the Netherlands, Organizing Committee of the 8th International Symposium on Tropospheric Profiling, 2009. a

Sekelsky, S. M.: Near-field reflectivity and antenna boresight gain corrections for millimeter-wave atmospheric radars, J. Atmos. Ocean. Tech., 19, 468–477, 2002. a

Shupe, M. D., Kollias, P., Matrosov, S. Y., and Schneider, T. L.: Deriving mixed-phase cloud properties from Doppler radar spectra, J. Atmos. Ocean. Tech., 21, 660–670, 2004. a, b

Shupe, M. D., Walden, V. P., Eloranta, E., Uttal, T., Campbell, J. R., Starkweather, S. M., and Shiobara, M.: Clouds at Arctic atmospheric observatories. Part I: Occurrence and macrophysical properties, J. Appl. Meteorol. Clim., 50, 626–644, 2011. a, b

Silber, I., Fridlind, A. M., Verlinde, J., Ackerman, A. S., Cesana, G. V., and Knopf, D. A.: The prevalence of precipitation from polar supercooled clouds, Atmos. Chem. Phys., 21, 3949–3971,, 2021. a

Sinclair, V. A., Moisseev, D., and von Lerber, A.: How dual-polarization radar observations can be used to verify model representation of secondary ice, J. Geophys. Res.-Atmos., 121, 10–954, 2016. a, b

Spek, A., Unal, C., Moisseev, D., Russchenberg, H., Chandrasekar, V., and Dufournet, Y.: A new technique to categorize and retrieve the microphysical properties of ice particles above the melting layer using radar dual-polarization spectral analysis, J. Atmos. Ocean. Tech., 25, 482–497, 2008. a

Sullivan, S. C., Barthlott, C., Crosier, J., Zhukov, I., Nenes, A., and Hoose, C.: The effect of secondary ice production parameterization on the simulation of a cold frontal rainband, Atmos. Chem. Phys., 18, 16461–16480,, 2018. a, b

Sun, 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

Tridon, F. and Battaglia, A.: Dual-frequency radar Doppler spectral retrieval of rain drop size distributions and entangled dynamics variables, J. Geophys. Res.-Atmos., 120, 5585–5601, 2015. a

Tyynelä, J. and Chandrasekar, V.: Characterizing falling snow using multifrequency dual-polarization measurements, J. Geophys. Res.-Atmos., 119, 8268–8283, 2014. a

Tyynelä, J., Leinonen, J., Moisseev, D., and Nousiainen, T.: Radar backscattering from snowflakes: Comparison of fractal, aggregate, and soft spheroid models, J. Atmos. Ocean. Tech., 28, 1365–1372, 2011. a, b, c, d

Vassel, M., Ickes, L., Maturilli, M., and Hoose, C.: Classification of Arctic multilayer clouds using radiosonde and radar data in Svalbard, Atmos. Chem. Phys., 19, 5111–5126,, 2019. a

Verlinde, J., Rambukkange, M. P., Clothiaux, E. E., McFarquhar, G. M., and Eloranta, E. W.: Arctic multilayered, mixed-phase cloud processes revealed in millimeter-wave cloud radar Doppler spectra, J. Geophys. Res.-Atmos., 118, 13–199, 2013. a

Westbrook, C. and Illingworth, A.: The formation of ice in a long-lived supercooled layer cloud, Q. J. Roy. Meteor. Soc., 139, 2209–2221, 2013. a

Wex, H., Huang, L., Zhang, W., Hung, H., Traversi, R., Becagli, S., Sheesley, R. J., Moffett, C. E., Barrett, T. E., Bossi, R., Skov, H., Hünerbein, A., Lubitz, J., Löffler, M., Linke, O., Hartmann, M., Herenz, P., and Stratmann, F.: Annual variability of ice-nucleating particle concentrations at different Arctic locations, Atmos. Chem. Phys., 19, 5293–5311,, 2019. a

Wildeman, S., Sterl, S., Sun, C., and Lohse, D.: Fast dynamics of water droplets freezing from the outside in, Phys. Rev. Lett., 118, 084101,, 2017. a, b

Wilson, T. W., Ladino, L. A., Alpert, P. A., Breckels, M. N., Brooks, I. M., Browse, J., Burrows, S. M., Carslaw, K. S., Huffman, J. A., Judd, C., Kilthau, W. P., Mason, R. H., McFiggans, G., Miller, L. A., Nájera, J. J., Polishchuk, E., Rae, S., Schiller, C. L., Si, M., Vergara Temprado, J., Whale, T. F., Wong, J. P. S., Wurl, O., Yakobi-Hancock, J. D., Abbatt, J. P. D., Aller, J. Y., Bertram, A. K., Knopf, D. A., and Murray, B. J.: A marine biogenic source of atmospheric ice-nucleating particles, Nature, 525, 234–238, 2015. a, b

Yang, J., Wang, Z., Heymsfield, A. J., DeMott, P. J., Twohy, C. H., Suski, K. J., and Toohey, D. W.: High ice concentration observed in tropical maritime stratiform mixed-phase clouds with top temperatures warmer than −8C, Atmos. Res., 233, 104719,, 2020.  a, b, c, d

Young, G., Lachlan-Cope, T., O'Shea, S., Dearden, C., Listowski, C., Bower, K., Choularton, T., and Gallagher, M.: Radiative effects of secondary ice enhancement in coastal Antarctic clouds, Geophys. Res. Lett., 46, 2312–2321, 2019. a

Zängl, G., Reinert, D., Rípodas, P., and Baldauf, M.: The ICON (ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core, Q. J. Roy. Meteor. Soc., 141, 563–579, 2015. a

Zawadzki, I., Fabry, F., and Szyrmer, W.: Observations of supercooled water and secondary ice generation by a vertically pointing X-band Doppler radar, Atmos. Res., 59, 343–359, 2001. a, b, c, d, e

Zhao, X. and Liu, X.: Global Importance of Secondary Ice Production, Geophys. Res. Lett., 48, e2021GL092581,, 2021. a

Zhao, X., Liu, X., Phillips, V. T. J., and Patade, S.: Impacts of secondary ice production on Arctic mixed-phase clouds based on ARM observations and CAM6 single-column model simulations, Atmos. Chem. Phys., 21, 5685–5703,, 2021.bibl a, b

Short summary
In natural clouds, ice-nucleating particles are expected to be rare above –10 °C. In the current paper, we found that the formation of ice columns is frequent in stratiform clouds and is associated with increased precipitation intensity and liquid water path. In single-layer shallow clouds, the production of ice columns was attributed to secondary ice production, despite the rime-splintering process not being expected to take place in such clouds.
Final-revised paper