Six SIP mechanisms have been discussed over the past decades. They are expected to be active depending on the environmental conditions like temperature, cloud particle size distribution, and updrafts: (i) the rime-splintering or Hallett–Mossop process, (ii) collisional breakup, (iii) droplet fragmentation upon freezing, (iv) thermal-shock fragmentation, (v) fragmentation of sublimating or melting ice, and (vi) activation of INPs in transient supersaturation of freezing droplets .

The rime-splintering process describes the production of secondary ice during riming. Yet, the underlying physical mechanism of this process is not well defined . The most common explanation for SIP during riming is that droplets, which freeze on an ice crystal, build up an internal pressure and break up upon freezing. In a laboratory study, observed that rime-splintering is active at temperatures between -3 and -8 ∘C when particles grow by riming. showed that at temperatures above -3 ∘C, droplets spread out on the ice surface and do not build an ice shell, while at temperatures below -8 ∘C the ice shell of rime might be too strong to break . While several observations of secondary ice could be explained by the rime-splintering process e.g.,, secondary ice was also observed in cases where these requirements were not fulfilled e.g.,, or the SIP by rime-splintering alone was expected to be too slow to explain the observed rapid glaciation e.g.,.

Collisional breakup describes the break-up of an ice crystal due to a collision with another ice crystal. This process was shown to be active when large ice crystals (e.g., graupel or aggregates) are present . Fragments of ice crystals were observed in several field studies e.g.,. found that significant SIP occurs by collisional breakup when relatively large concentrations of rimed ice crystals are present in convective clouds. He concluded that collisional breakup cannot explain high secondary ICNCs at warm temperatures. The laboratory study of showed that a maximum SIP rate occurs at -16 ∘C when ice spheres with different growth modes collide.

Droplet fragmentation during freezing can happen when liquid water gets trapped inside a freezing droplet after an ice shell formed around the droplet and expands. Due to the lower density of ice than liquid water, an internal pressure builds up. If the pressure reaches a critical point, different processes have been observed to release the pressure: a complete breakup in mainly two halves; opening and closing of cracks; bubbles forming and bursting on the ice shell; and streams being ejected from the shell, which may contain small ice fragments . An average number of splinters being ejected per freezing droplet could not be quantified so far because many ejected splinters are expected to be too small, and the processes happen so fast that the used measurement techniques until today were not able to detect all splinters . showed that droplets freezing in moving air break up 1 order of magnitude more frequently than droplets freezing in stagnant air. Therefore, high wind speed and turbulence can be expected to increase the likelihood of droplet fragmentation. The highest fragmentation rate was observed around -15 ∘C. It decreases to higher and lower temperatures, while it seems to increase again at temperatures close to 0 ∘C . The larger the droplet, the more likely it will fragment and the more splinters it will likely produce . Many field studies showed that large droplets are present in clouds before ICNCs exceed the INPCs by orders of magnitude e.g., and often explained these observations with SIP by droplet fragmentation. Even though large droplets are rare in clouds, they might start a cascading process of splinter production when produced splinters hit other droplets, which subsequently freeze and produce more splinters .

Thermal-shock fragmentation and the activation of INPs in transient supersaturation of freezing droplets both require precipitation-sized droplets to be active, while the fragmentation of sublimating ice can be important if cloud regions are affected by entrainment of dry air . observed no sublimation breakup for columnar and plate-like crystals, while breakup during melting only took place at relative humidities below 90 %.

From the six SIP mechanisms, the rime-splintering process is the best-constrained mechanism and has been implemented in some cloud microphysics schemes e.g.,. There have been fewer attempts to include collisional breakup e.g., and droplet fragmentation e.g.,. modeled all three ice multiplication processes and showed that none of the three SIP mechanisms dominates the ICNC enhancement. While collisional breakup is important when low updrafts prevail, and high nucleation rates are present, droplet fragmentation and the rime-splintering process can also be important in INP-limited regions. However, they pointed out that the processes are not very well constrained by laboratory and in situ data.

A remarkable finding of the study by is the observation of small pristine ice crystals persisting immediately above the melting layer in clouds. They calculated a spatial correlation time τcorr during which the environmental changes (e.g., air temperature, humidity, and cloud particle number concentration) are insignificant, and the shapes of the ice crystals can still be associated with the environment they are growing in. Based on their estimation, τcorr is on the order of 60–120 s. Assuming water vapor saturation over liquid, an ice column at temperatures below -3 ∘C can reach a length between 50 and 150 µm during τcorr depending on its aspect ratio. A relative humidity close to water saturation is a valid assumption in MPCs . It can be concluded that if external processes like blowing snow or the seeder–feeder process can be excluded, the discrepancy between the INPC and the concentration of columns shorter than 50 to 150 µm in MPCs must emerge from SIP.

suggested a recirculation process through the melting layer to be a possible explanation for the observed high concentration of small pristine ice crystals. Large ice crystals which fall through the melting layer as precipitation melt into drizzle-sized droplets. If convection or turbulence causes large updrafts, these droplets can be reintroduced into the cloud. Whenever these droplets freeze by the collision with ice crystals, the droplet-fragmentation-upon-freezing process can become active and produce high numbers of secondary-ice particles.

The measurements used in this study were collected during the RACLETS campaign, which took place in February and March 2019 in the region of Davos in the Swiss Alps. The objective of this campaign was to improve the understanding of the influence of topography and aerosols on the development of clouds. For this goal, a set of instruments were deployed at different locations. Figure  shows the locations and instruments which were used for the analysis of the presented case study.

Measurements of the in-cloud properties such as the cloud particle concentration as well as their size distribution and shape were taken from the HoloGondel platform on the Gotschnabahn, which is described in more detail in Sect. . As part of the MeteoSwiss observation network, a ceilometer (Vaisala, Model CL31) was installed in Klosters (1200 m) to determine the cloud base height as well as a wind profiler (Vaisala, Model Lap3000, Finland) in Wolfgang to measure the horizontal wind field . The general wind pattern on the ground was determined by data from different MeteoSwiss stations, the Snow and Avalanche Research SLF, and the Holfuy station of the Grischa paragliding club. Temperature measurements were taken from the MeteoSwiss station in Klosters and the snow drift station installed at Gotschnagrat . For the analysis of the whole clouds, a vertically pointing cloud radar (Model Mira-36, METEK GmbH, Germany; ) was installed in Wolfgang. INPCs were measured in the valley in Wolfgang (1630 m) and on the mountain top of Weissfluhjoch (2670 m), where ambient aerosol was sampled through heated inlets using a high-flow-rate impinger (Bertin Technologies, Model Coriolis µ, France). The samples were analyzed for the INPC immediately on-site with a drop-freezing technique using DRINCZ (Drop Freezing Ice Nuclei Counter Zurich; ) of the Swiss Federal Institute of Technology in Zurich (ETH) in Wolfgang and the University of Basel's LINDA (LED-based ice nucleation detection apparatus; ) at Weissfluhjoch as described in . The detection limit of the INPC (lowest concentration measurable) was calculated according to considering the concentration obtained if only one (first) drop froze. A detection limit of 5.2×10-4 and 4.9×10-4 L-1 has been determined for Wolfgang and Weissfluhjoch, respectively.

The HoloGondel platform consists of the HOLographic Imager for Microscopic Objects 3G (HOLIMO 3G in ) – which records the concentration, size distribution, and shapes of cloud particles – and a temperature and relative-humidity sensor (HygroMet4, Rotronic) in a ventilated housing (RS24T, Rotronic) as well as a pressure sensor (Fig. ). During the RACLETS campaign, the platform was installed on one of the gondolas running at the upper section between the Gotschnagrat mountain station (2280 m) and the middle station at Gotschnaboden (1790 m), covering a horizontal distance of about 830 m. To avoid any influence by the gondola stations on the measurements, data were only used if the difference between the pressure measured on the gondola and the pressure measured at the stations was more than 1.5 hPa, corresponding to a vertical distance of more than 15 m and a total distance of about 30 m between the gondola and the stations. The gondola runs with a maximum speed of about 6 m s-1, leading to a total measurement time of about 140 s per ride. To avoid influences from the gondola and its swing arm, only measurements of uphill rides were analyzed, when the setup was in front of the gondola in the direction of travel.

HOLIMO 3G is an open-path instrument, which uses digital in-line holography. Holograms of the sample volume between the two towers (see Fig. b) are recorded, of which 13.7 cm3 was considered for the analysis, producing a 3D distribution of the cloud particles between the two towers and a 2D image of each cloud particle. A more detailed description of the measurement principle can be found in and . The pixel size is 3.1 µm, which allows us to observe cloud particles larger than 6.2 µm . The differentiation between ice and liquid is based on the shape of the particles (circular vs. non-circular). This is possible for particles larger than approximately 25 µm, depending on the shape of the ice particles .

A neural network for particles larger than 25 µm and a decision tree for particles smaller than 25 µm in their major axis were used to separate cloud droplets from ice crystals and sort out artifacts, which are falsely identified as cloud particles by the software. All ice crystals were manually confirmed after the automated classification. Therefore, the uncertainty in the concentration of ice particles can be estimated with ±5 % for ice crystals larger than about 100 µm and ±15 % for ice crystals smaller than 100 µm . All ice particles were also manually classified into the habits “plates”; “columns”; “irregular”; “aged ice” (rimed particles and aggregates); and “unidentified”, referring to ice crystals, which could not be classified because they were too small or because their orientation did not allow a decision on their habit. Plates and columns are hexagonal prisms with the diameter of the hexagonal basis a and the prism height h. Plates have a dimension of a>h, while columns have a dimension of a≤h. For the uncertainty in the different habits, the counting uncertainty (N/V; N: number of crystals; V: measurement volume) was added because of their relatively low number compared to the measurement volume. The uncertainty in cloud droplets is estimated to be ±6 % as determined for the classification with the neural network in . Again, for droplets larger than 40 µm the counting uncertainty was added due to their relatively small numbers.

On 22 February 2019, a ridge was located in front of Switzerland, which guided air masses from the north to the Davos region. A trough over Russia reached relatively far to the south supporting the rise of air masses over Germany, Austria, and eastern Switzerland, which were all covered by stratiform clouds. From the cloud radar measurements (Fig. a), it can be anticipated that the cloud top height was about 4 km in the Davos region, referring to a cloud top temperature of around -8 ∘C as measured by a radiosonde launched at 10:25 UTC from Wolfgang. Between 05:00 and 09:00 UTC, the cloud radar showed a midlevel cloud above the stratiform cloud, reaching from about 4.5 to about 6.5 km, that might have acted as a seeder cloud in the early morning.

As can be inferred from the ceilometer data, light precipitation started at around 07:00 UTC in Klosters, where the valley station of the Gotschnabahn is located (Fig. b). The cloud base was located slightly above the height of the Gotschnaboden during the precipitation period, which lasted until about 10:30 UTC, when the cloud began to dissipate. Measurements with HoloGondel were taken between 08:00 and 10:00 UTC. The rides used for the analysis are shown as gray lines in Fig. . The temperatures and relative humidities at Gotschnaboden and Gotschnagrat (Fig. c and d) are derived from measurements of the highest and lowest point of the measurements on the gondola and are therefore only available when the gondola was in operation and close to one of the stations. During the measurement period, the temperature in Klosters increased from about +1.5 to +3.5 ∘C at a relative humidity of about 80 %, while the measurements on the gondola were taken mainly between 0 and -2.5 ∘C at relative humidities above 95 %. The melting layer can be inferred from the dark band of the ceilometer data below Gotschnaboden (Fig. b), which was measured above Klosters and may slightly differ at Gotschnaboden.

The main wind direction was from north to northeast as can be inferred from the wind profiler measurements taken in Wolfgang in Fig. b. However, the ground measurements show that the wind direction was strongly influenced by the orography (Fig. a). The valley north of the Gotschnabahn forced the wind direction to the northwest as can be seen from the wind measurements on Gotschnagrat. Thus, air masses measured on Gotschnabahn were pushed up the valley coming from northwest and must have been lifted through the melting layer before reaching the measurement site.

The slope inclination in the wind direction measured at Gotschnagrat is about 7∘ (averaged over Gotschnagrat and 1 km northwest). Assuming that the measured horizontal wind of about 5 m s-1 (Fig. ) was blown up the slope without friction, the vertical wind speed reached about 0.6 m s-1. The eddy dissipation rate, which is a measure for turbulence, was calculated from the wind profiler and cloud radar measurements as described in and reached values between 20 and 90 cm2 s-3 at an altitude of 1800 to 2300 m between 08:00 and 10:00 UTC over Wolfgang. Such values were reported in stratiform clouds as well as in cumulus clouds with weak updrafts . Since no measurements of turbulence are available over Klosters, we assume that similar turbulence is present at the measurement site.

Between 08:00 and 10:00 UTC, a total volume of 43.7 L was measured by HOLIMO, distributed over nine uphill rides, each contributing between 2.6 and 8.2 L. The difference in the measurement volumes is due to different onsets and offsets of the automated recording. HOLIMO observed aged ice crystals (here defined as aggregates and rimed particles), irregular particles, plates, and very few small columns (see Fig. ). About 35 % of all ice crystals were classified as unidentified because they were too small, or the particle orientation made a decision on its habit inconclusive. As discussed in Sect. , unidentified particles as well as plates smaller than 93 µm are referred to as small plates. They are marked with dots in the histogram plot of the ice crystal size distribution in Fig. .

The concentrations of ice crystals and cloud droplets stayed relatively constant over the measurement period (see Fig. ). Taking the whole measurement period into account, a mean cloud droplet number concentration (CDNC) of 156 cm-3 ± 9 cm-3 and a mean ICNC of 6 L-1 ± 1 L-1 were measured. Between the different rides the CDNC varied between 67 and 215 cm-3 and the ICNC between 4 and 18 L-1. The most commonly observed pristine ice particle habit was plates, with a mean concentration of 0.8 L-1 ± 0.2 L-1, varying between 0.1 and 1.5 L-1 between the different rides. Note that this is a lower estimate since the software has difficulties in detecting transparent particles, and depending on the size and orientation, plates may be hard to identify. Besides the ice particle habits, a mean concentration of droplets larger than 40 µm in diameter of 0.3 L-1 ± 0.1 L-1, varying between 0 and 0.8 L-1 between the different rides, was observed. It should be mentioned that the missing observations of plates or large droplets during single rides could be a stochastic effect due to the relatively low concentrations compared to the measurement volumes of each ride; i.e., for the smallest measurement volume of 2.6 L-1, one observed particle results in a particle concentration of 0.4 L-1. In the following, we consider the whole measurement period and site as one measurement volume, where the cloud particle concentrations are constant.

The INPC was measured in the valley in Wolfgang (1630 m) and on the mountain top of Weissfluhjoch (2670 m) (see Fig.  for the geographical positions) simultaneously at 08:00 and 10:00 UTC on 22 February 2019. No INP was detected for temperatures higher than -6 ∘C at both sites. Therefore, the upper estimate of the INPC at the measurement site (T>-3 ∘C) is equal to the detection limit and on the order of 5×10-4 L-1 for both sites, which lies several orders below the measured ICNC.

Since the INPC at the measurement location lies several orders of magnitude below the ICNC, the contribution from primary-ice nucleation will be neglected. Therefore, the concentration of secondary-ice crystals can be estimated from the concentration of ice crystals that have newly formed and grown at the same location. The measurements were taken at temperatures warmer than -3 ∘C. In this temperature regime, newly formed ice crystals grow into plates . Figure  shows the shapes and the size distribution of the ice crystals divided logarithmically into nine size bins. The histogram plot shows that the majority of the classified particles smaller than 93 µm were plates (54 % if unidentified particles are excluded; Fig. ), while aged crystals can only be found above that size. Hence, we take 93 µm as a threshold to divide between newly formed ice crystals and those which sedimented from above.

About 70 % of the ice crystals smaller than 93 µm (from now on referred to as “small ice”) could not be classified into a certain habit because the resolution was too low. Yet, a transparent part observed in the middle of many of these particles (see Fig. ) supports the assumption based on the temperature regime that these ice crystals are actual plates. In the following, small ice classified as plates including the unidentified particles is referred to as “small plates”.

A plate with a maximum dimension of 93 µm has a fall velocity of about 0.05 m s-1 (using the equations given in ), which is 1 order of magnitude lower than the estimated updraft of about 0.6 m s-1. This supports the assumption that small plates are unlikely to have sedimented from above. About 4 % of the small ice was classified as columns. The few observed columns have an aspect ratio close to 1 (see Fig. ), which is typical for columns growing at temperatures just below -3 ∘C at low supersaturation , and therefore they must have sedimented from slightly higher altitudes. This is possible because columns have a higher fall velocity than plates (a column with a length of 93 µm falls with about 0.3 m s-1 using the equations given in ). About 10 % of the small ice was classified as irregular particles. These could be large secondary-ice splinters. However, the fall velocity of irregular particles is hard to assess, and it remains unclear if they have fallen from above or formed at the measurement site by SIP. Based on these assumptions, the lower estimate of the secondary-ice concentration, which has formed and grown at the measurement site, is equal to the concentration of small plates as defined above and is on the order of 2.6 L-1 ± 0.6 L-1.

The time a plate needs to grow to 93 µm is hard to assess at temperatures close to 0 ∘C because small variations in the environmental conditions, e.g., a temperature fluctuation of only 0.5 ∘C, can change the growth time on the order of several minutes. To calculate the growth time, we use the general equation of the ice particle growth by vapor diffusion given in . Using the mass–size relation of hexagonal plates from and assuming water vapor saturation over liquid at a temperature of -2 ∘C, a splinter of 5 µm needs about 9 min to grow to a size of 93 µm and about 5 min to grow to a size of 60 µm. Thus, for any given time, a plate of 93 µm is 4 min older than a plate of 60 µm if the temperature is constant at -2 ∘C, and the water vapor is saturated over water. Therefore, plates with sizes between 60 and 93 µm must have newly formed by SIP within a time span of 4 min. The same calculation can be done for a bigger size interval, e.g., plates between 39 and 93 µm, which must have newly formed by SIP during a time span of about 6 min. Taking temperature variations in the measurement volume into account, i.e., varying the temperature between -2 and -1 ∘C, plates between 60 and 93 µm were newly formed by SIP during a time span of 4 to 6 min, while plates between 39 and 93 µm were newly formed by SIP during a time span of 6 to 10 min depending on the temperature. The observed concentration of small plates with sizes between 60 and 93 µm was 0.9 L-1 ± 0.3 L-1, whereas between 39 and 93 µm the concentration of small plates was 2.1 L-1 ± 0.5 L-1 averaged over the 2 h measurement period. Thus, 0.9 L-1 ± 0.3 L-1 of secondary ice has formed within 4 to 6 min and 2.1 L-1 ± 0.5 L-1 within 6 to 10 min. Taking all named uncertainties into account, the rate of secondary-ice production during our case study was 0.24 L-1 min-1 ± 0.09 L-1 min-1.

Droplets larger than 40 µm up to a size of 380 µm in diameter were observed between 0 and -2.7 ∘C. explained the occurrence of large droplets just above the melting layer with a recirculation process (see Sect. ). Ice particles fall through the melting layer and melt into drizzle drops, which can be reintroduced into the cloud by sufficiently high updrafts. The estimated updraft in this case study is about 0.6 m s-1, which is equal to the fall speed of a 150 µm droplet . Most of the observed droplets were smaller than 150 µm and could be lifted up in the clouds by the updraft, while the remaining ones could have been brought into the cloud by local turbulences. Also note that the updraft velocity is a very rough estimate; increasing it by 0.2 m s-1 is already enough to lift 86 % instead of 79 % of the observed droplets.

Droplet fragmentation requires the presence of large droplets (>≈40 µm; e.g.,) and has no temperature constraint towards high temperatures. It is therefore a possible mechanism to explain the observed secondary ice. Collisional breakup cannot be completely ruled out, but it is not expected to produce high secondary ICNCs at warm temperatures as outlined in the introduction. Furthermore, the missing observation of fragment-like ice particles (e.g., broken-off branches) is evidence that collisional breakup was not very active. Moreover, fragments from collisional breakup are not necessarily small particles and therefore could be miscounted as aged or irregular ice in this study.

The rime-splintering process and ice fragmentation during thermal shock can both be excluded from being active because they require lower temperatures. Additionally, INP activation in transient supersaturation around freezing drops is unlikely to be active in a relatively warm environment . As MPCs are supersaturated with respect to ice, the requirements for ice fragmentation during sublimation are also not fulfilled .

Ice crystals can partially melt and create fragments, which can be lifted back into the cloud. However, the ice crystals have to be sufficiently small to be lifted, while at the same time, small ice crystals are likely to completely melt before re-entering the cloud. Furthermore, the breakup rates of this process depend on temperature and humidity and largely on the initial shape of the ice crystals . observed no sublimation breakup for columnar and plate-like crystals and breakup during melting only at relative humidities below 90 %. The shapes of ice crystals in our study are mostly solid particles that have columnar and plate-like shapes (Fig. ), and the relative humidity on the gondola never dropped below 95 % (Fig. d). Taking all these considerations into account, we assume that ice fragmentation during melting has a negligible effect on the SIP.

Therefore, we expect that most of the small plates emerged from droplet fragmentation and that the orographically induced updrafts serve as a constant supply of new droplets larger than 40 µm in diameter, which originate from melted ice crystals. A schematic of this recirculation process in mountainous regions is shown in Fig. . Because new droplets are continuously provided, we can assume that our measurements reflect a steady state such that we observe freshly produced as well as aged particles simultaneously. The combination of SIP above the melting layer and high updrafts on the windward slope will transport the secondary-ice crystals to higher altitudes, where they influence the cloud microphysics.

Here we derive a parametrization of the SIP by droplet fragmentation at temperatures close to 0 ∘C, when primary-ice nucleation can be neglected and droplets freeze only by the collision with ice crystals, which either sedimented from above or formed by SIP. Like , we assume that only droplets larger than 40 µm are likely to contribute to SIP by droplet fragmentation. To calculate the splinter generation rate of a droplet with diameter dn (gsp(dn)), the droplet freezing rate by collision (fcol(dn)) has to be multiplied by the droplet fragmentation probability during freezing (pdf(dn)) and the number of splinters per fragmenting droplet (Nsp(dn)): 1gspdn=fcoldn⋅pdfdn⋅Nspdn. To obtain the total splinter generation rate (Gsp) in a volume V, which contains Nd>40µm droplets with diameters larger than 40 µm, the sum of gsp(dn) over all droplets with diameter dn (n=1, 2, …, Nd>40µm) has to be divided by V: 2Gsp=1V∑n=1Nd>40µmgspdn. The probability that a droplet with a fall velocity of v(dn) collides with an ice crystal with an equivalent particle area diameter aj and freezes is equal to the collision efficiency E(dn, aj). It needs to be multiplied by the combined cross-section of the ice crystal and the droplet (0.25⋅π⋅(dn+aj)2) and the absolute relative velocity between those two (|v(aj)-v(dn)|) to obtain the collection kernel. To calculate the freezing rate of a droplet in volume V, the sum of the collection kernel between the droplet and each ice crystal in the same volume with an equivalent particle area diameter aj (j=1, 2, …, Nice; Nice: number of ice crystals) has to be divided by V: 3fcoldn≈1V∑j=1NiceEdn,aj⋅|vaj-vdn|⋅π⋅dn+aj24. The probability of droplet fragmentation is size-dependent . For a fragmentation to occur, the surface energy has to be overcome, which is proportional to dn2. Therefore, we expect pdf(dn)=c⋅dn2 for some constant c. The constant c can be estimated by assuming a fragmentation probability of 40 % for droplets with dn=300 µm at temperatures higher than -2.5 ∘C, i.e., pdf(dn=300µm)=0.4, as measured by . This yields a droplet fragmentation probability per freezing droplet of 4pdfdn≈4.4×106m-2⋅dn2, which differs from other studies, where the probability of droplet fragmentation follows a normal distribution centered at -15 ∘C independent of droplet size, with a probability close to 0 % for slightly sub-zero temperatures e.g.,. However, the study of showed that droplet fragmentation is also likely to happen at temperatures close to 0 ∘C.

The concentration of splinters produced per fragmentation event could not be quantified until now because many of the splinters being produced during a fragmentation event might be too small to be observed with the available measurement techniques . Here we consider the maximum number of splinters observed during a breakup event for different droplet sizes, which were summarized by from different studies, to be the best estimate available. The data points suggest a linear correlation between the number of splinters being produced per droplet Nsp(dn) and the droplet diameter dn. Applying a linear regression, Nsp can be estimated as 5Nspdn≈9×104m-1⋅dn. The splinter generation rate in volume V at temperatures close to 0 ∘C in MPCs is subsequently given by 6Gsp≈9.9×1010⋅πm-3⋅1V2∑n=1Nd>40µmdn3∑j=1NiceEdn,aj⋅|vaj-vdn|⋅dn+aj2.

Several other studies observed an increase in the concentration of small ice crystals in the presence of large supercooled drops in clouds e.g.,. Updrafts were made responsible for the origin of these large droplets in all of the studies. The measurements were taken mostly at temperatures lower than during our case study, where newly formed ice crystals grow into columns. The recirculation process through the melting layer described above (Fig. ) is therefore also expected to play a role for droplet fragmentation in higher regions of the cloud if the updrafts are strong enough to lift drizzle drops high enough until they freeze. Besides droplet fragmentation, the rime-splintering process is expected to be active in the temperature regime between -3 and -8 ∘C, which makes the assignment of the observed secondary ICNC to a specific process more difficult. However, images of deformed frozen drops in all of the abovementioned studies are strong evidence that droplet fragmentation was active and should be accounted for besides the rime-splintering process.