Articles | Volume 23, issue 4
Research article
20 Feb 2023
Research article |  | 20 Feb 2023

Heavy snowfall event over the Swiss Alps: did wind shear impact secondary ice production?

Zane Dedekind, Jacopo Grazioli, Philip H. Austin, and Ulrike Lohmann

The change in wind direction and speed with height, referred to as vertical wind shear, causes enhanced turbulence in the atmosphere. As a result, there are enhanced interactions between ice particles that break up during collisions in clouds which could cause heavy snowfall. For example, intense dual-polarization Doppler signatures in conjunction with strong vertical wind shear were observed by an X-band weather radar during a wintertime high-intensity precipitation event over the Swiss Alps. An enhancement of differential phase shift (Kdp>1 km−1) around −15C suggested that a large population of oblate ice particles was present in the atmosphere. Here, we show that ice–graupel collisions are a likely origin of this population, probably enhanced by turbulence. We perform sensitivity simulations that include ice–graupel collisions of a cold frontal passage to investigate whether these simulations can capture the event better and whether the vertical wind shear had an impact on the secondary ice production (SIP) rate. The simulations are conducted with the Consortium for Small-scale Modeling (COSMO), at a 1 km horizontal grid spacing in the Davos region in Switzerland. The rime-splintering simulations could not reproduce the high ice crystal number concentrations, produced too large ice particles and therefore overestimated the radar reflectivity. The collisional-breakup simulations reproduced both the measured horizontal reflectivity and the ground-based observations of hydrometeor number concentration more accurately (∼20 L−1). During 14:30–15:45 UTC the vertical wind shear strengthened by 60 % within the region favorable for SIP. Calculation of the mutual information between the SIP rate and vertical wind shear and updraft velocity suggests that the SIP rate is best predicted by the vertical wind shear rather than the updraft velocity. The ice–graupel simulations were insensitive to the parameters in the model that control the size threshold for the conversion from ice to graupel and snow to graupel.

1 Introduction

In clouds, ice particles play an important role for Earth’s radiation budget and precipitation formation. Precipitation originates predominantly from mixed-phase clouds (MPCs) in the midlatitudes, especially over continental regions (Mülmenstädt et al.2015; Heymsfield et al.2015, 2020). The formation of ice particles, therefore, needs to be described adequately if any attempt is made to understand the evolution of MPCs and ice clouds.

Ice formation can occur through primary and secondary ice production (SIP) processes. Primary ice production includes homogeneous freezing of supercooled liquid water at temperatures (T) <-38C, while heterogeneous ice nucleation of supercooled liquid water dominates at warmer subzero temperatures (T>-38C). After the first formation of ice particles, secondary ice processes may occur. In a narrow temperature range, -3T-8C, rime splintering (Hallett and Mossop1974) can occur when supercooled cloud droplets collide with ice particles, freeze from the outside in and shatter as a result of internal pressure buildup. Rime splintering has been studied extensively in models but has been shown to be inadequate to capture SIP in wintertime orographic MPCs (Henneberg et al.2017; Dedekind et al.2021; Georgakaki et al.2022), producing ice number concentrations that are orders of magnitude less than observed. Ice–ice collisions have been more widely used in models in the last decade (Yano and Phillips2011; Phillips et al.2017; Sullivan et al.2018; Hoarau et al.2018; Sotiropoulou et al.2020; Zhao et al.2021) since they were first studied in the laboratory about 4 decades ago (Vardiman1978; Takahashi et al.1995). SIP as a result of ice–ice collisions was shown to contribute significantly to the ice crystal number concentrations and thereby explain the discrepancy between models and observations in the Arctic (Sotiropoulou et al.2020; Zhao et al.2021), Antarctic (Sotiropoulou et al.2021) and midlatitudes (Sullivan et al.2018; Dedekind et al.2021; Georgakaki et al.2022). The enhancement of smaller ice particles triggers an increase in the combined growth rates (reduced riming due to the smaller ice crystals but strong enhancement of deposition) of up to 33 %, resulting in larger latent heat release and stronger updraft velocities (Dedekind et al.2021). When ice–ice collisions occur in wintertime orographic MPCs, the general tendency is for riming to decrease. Hence, the depositional growth rate dominates the growth rates of ice particles. Due to the stronger updrafts, ice particles are lofted to higher regions within the cloud reducing the local precipitation rates.

The impact of turbulence associated with baroclinic waves on cloud water and precipitation formation is well known (Baumgartner and Reichel1975; Houze and Medina2005; Medina and Houze2015). Updrafts on the scale of ∼10 km from baroclinic waves have properties of shear-induced turbulence, and it is these small cells of enhanced updraft and turbulence that drive orographic precipitation (Medina and Houze2015). In regions associated with mountainous terrain, strong shear layers at low levels approaching a barrier were emphasized by Houze and Medina (2005) to set up turbulence which in turn aids in precipitation growth (by accretion) on the windward side of a mountain. Medina et al. (2005) showed in idealized simulations that a shear layer can develop as a response to flow over the terrain, by which they concluded that this mechanism, in actual topography, caused turbulent overturning which enhanced precipitation formation. In their simulation, the precipitation formation was linked to enhanced accretion (see also Medina and Houze2015). The probability of interactions between cloud hydrometeors, whether through riming and/or aggregation, increases with increasing turbulence and aids in the rapid formation of precipitation regardless of whether the turbulence is associated with orographic flow regimes or with warm conveyor belts (Houze and Medina2005; Gehring et al.2020). These interactions are not limited to the accretional growth of cloud hydrometeors but also include the fracturing of ice particles in ice–ice collisions, enhancing SIP. Dedekind et al. (2021) hypothesized that ice–graupel collisions are sensitive to the rate at which graupel forms, which is a function of the ice particle size and the riming rate. In the Seifert and Beheng (2006) two-moment (2M) cloud microphysics scheme used in this study, ice crystals or snow particles undergoing riming can only be converted to graupel once they reach a size of 200 µm.

Remote sensing from weather radars is used to study snowfall microphysics and hydrometeor habit (e.g., shape, phase or hydrometeor type). Although radar observations do not provide direct information on SIP, a few studies leveraged the Doppler and/or dual-polarization capabilities of weather radars to identify the occurrence of SIP and to speculate, case by case, on the possible mechanisms behind its origin. Two approaches can be found in the literature. Zawadzki et al. (2001), Oue et al. (2015) and Luke et al. (2021) exploited Doppler spectra collected by vertically pointing radars to identify the appearance of secondary populations of particles at given altitudes or temperature levels. Other approaches (Hogan et al.2002; Andrić et al.2013; Sinclair et al.2016; Kumjian and Lombardo2017) focused on the interpretation of the signature of dual-polarization variables and their respective evolution over the vertical column of precipitation. This second approach, also used in this study, leverages the fact that dual-polarization variables are complementary and impacted differently by changes in number, shape, size and density of hydrometeors. Additional information (in situ data, models or a combination of more radars) is typically needed to increase confidence in the retrievals collected.

In this paper, we propose that the vertical wind shear associated with a cold front passage over the mountains of eastern Switzerland enhanced the formation of small and numerous oblate ice particles through ice–ice collisions. The ice–ice collisions explain the peculiar signatures in the data collected by a Doppler dual-polarization radar deployed in the region. We address the following questions:

  • Can these radar signatures be attributed to high ice crystal number concentrations linked to SIP other than rime splintering?

  • By including ice–graupel collisions in the model, can we simulate the high ice crystal number concentrations that were observed?

  • Was there a correlation between the vertical wind shear and SIP?

  • How sensitive are SIP rates to the conversion rate of ice particles to graupel?

2 Methods

2.1 Weather radar and two-dimensional video disdrometer (2DVD)

The principle of dual-polarization for weather radars relies on radars transmitting pulsed horizontally and vertically polarized waves (Field et al.2016). The waves interact with precipitation, and by looking at the differences in power and phase of the echoes in each polarization, information about the orientation, size and number concentration (and phase) of the hydrometeors being sampled can be retrieved. Horizontal (vertical) reflectivity ZH (ZV) and differential reflectivity ZDR are variables that exploit the power intensity of the echoes. ZH [dBz] increases as particles get larger, denser and/or more numerous. ZDR [dB] is the difference ZHZV and can be used to distinguish oblate particles (ZDR>0) from prolate ones (ZDR<0), while it has near-zero values for spherical particles. In an environment where preferentially oriented anisotropic ice particles are dominant, ZDR deviations from near-zero values are frequently observed (Bader et al.1987; Kumjian et al.2014). When ice particles form aggregates and become larger and less oblate, ZDR decreases while ZH increases (Schneebeli et al.2013; Kumjian et al.2014; Grazioli et al.2015a). The backscattered power is different for horizontal and vertical polarizations in the presence of anisotropic particles, as is the propagation speed of the waves.

The rate of change in phase shift between the horizontal and vertical polarized echoes is expressed by the specific differential phase shift Kdp [ km−1]. This variable is complementary and not redundant; it is in fact not affected by the absolute calibration of a radar and is less affected than ZDR by the eventual presence of large isotropic particles within the sampling volume. For instance, local Kdp enhancements during snowfall have been documented (Schneebeli et al.2013; Bechini et al.2013) and in some cases been associated with SIP (e.g., Andrić et al.2013; Grazioli et al.2015a; Sinclair et al.2016). Grazioli et al. (2015a) suggested in a case study that an increase of Kdp can be due to very large number concentrations of rimed anisotropic ice crystals resulting from ice–ice collisions. A recent study (von Terzi et al.2022) suggested that the Kdp enhancement was due to a combination of secondary ice production and an appropriate temperature range T-15C (where growth of planar crystals by vapor deposition, dendrites in particular, is maximized). Dendrites have very low densities, favor aggregation (hence the increase of ZH below Kdp peaks) and can easily fracture on impact with other ice particles.

An X-band dual-polarization mobile Doppler weather radar (MXPol) of the École Polytechnique Fédérale de Lausanne Environmental Remote Sensing Laboratory (EPFL-LTE) was set up at 2133 m above mean sea level (a.m.s.l.) on a ski slope overseeing the valley of Davos (Schneebeli et al.2013) from the southern side as shown in Fig. 1. Its exact location was 46.789 N, 9.843 E (see Sect. 2.1). MXPol is well suited for deployment in complex Alpine terrain or remote locations (e.g., Schneebeli et al.2013; Grazioli et al.2015a, 2017) and was operated from September 2009 to July 2011. The radar was routinely scanning over the valley of Davos in a sequence including pseudo-horizontal scans (fixed elevation and variable azimuth) and 2D vertical cross sections (fixed azimuth scans with elevation ranging from 0 to 90, known as range height indicator or RHI scans). One RHI scan in particular, used as a data source of this study, was conducted every 5 min towards the NE, at an azimuth of 22. Only observations collected at elevation angles below 40 are used, in order to limit the effect of elevation dependencies on the polarimetric variables (Ryzhkov et al.2005). MXPol provides single (ZH) and dual-polarization (ZDR, Kdp, ρHV) measurements as well as Doppler data which have proven useful in several snowfall microphysics studies (e.g., Schneebeli et al.2013; Grazioli et al.2015a; Kumjian and Lombardo2017; Oue et al.2021). Additionally, retrieval algorithms adapted to polarimetric data allow for estimated properties such as hydrometeor type (Grazioli et al.2015b, as used in this work) or, under given assumptions, microphysical quantities such as ice crystal number concentration Nt, median volume diameter Dm or ice water content (IWC). The hydrometeor classification method discriminates between three ice-phase dominant hydrometeor types: individual ice crystals, aggregates and rimed particles.

The microphysical quantities can be estimated from a combination of ZH, ZDR, Kdp and the radar wavelength following  Murphy et al. (2020):


In these equations, IWC is expressed in [g m−3], Nt in [L−1] and Dm in [mm]. λ is the radar wavelength in [mm], Zdr=100.1ZDR is the differential reflectivity in linear units and Zdp=100.1ZH-100.1ZV is the reflectivity difference in linear units [mm6 m−3]. More details about the derivations of these equation can be found in Ryzhkov and Zrnic (2019) and Murphy et al. (2020). The main assumptions are the following:

  • The equations are derived assuming to be in the Rayleigh regime, which may not be fulfilled for the X-band wavelength and large hydrometeors.

  • The density and the size of the hydrometeors are assumed to be inversely proportional, which is not fulfilled for hail.

The retrievals have shown to be most reliable at T<-10C, for low riming degrees and in regions where the Kdp and ZDR signals are not close to 0. As recognized by Murphy et al. (2020), the errors may be large, and in situ validation efforts are needed to refine these techniques. As a final caveat, the equations developed theoretically are in practice very sensitive to the accuracy of the polarimetric variables, which can be very noisy. Kdp in particular is an estimated variable affected by mean errors on the order of 30 % (Grazioli et al.2014a).

An additional ground-based source of information for this event is provided by a two-dimensional video disdrometer, 2DVD (for more information about this instrument at this location, see Grazioli et al.2014b), which was deployed on the opposite side of the Davos valley with respect to MXPol (46.830 N, 9.810 E; 2543 m a.m.s.l.). The 2DVD measures the size and fall velocity of hydrometeors captured within its measurement area of 11 cm × 11 cm. The 2DVD is used in this study as ground reference to quantify the number concentration of snowflakes (larger than 0.2 mm in terms of maximum dimension, according to the sensitivity of the instrument itself) at a temporal resolution of 5 min.

2.2 Model setup

2.2.1 Spatial and temporal resolution

The Consortium for Small-scale Modeling (COSMO; Baldauf et al.2011) non-hydrostatic model, version 5.4.1b, was used for this case study. COSMO has been used to study wintertime (Lohmann et al.2016; Henneberg et al.2017; Dedekind et al.2021) and summertime (Dedekind2021; Eirund et al.2021) orographic MPCs in the Swiss Alps. The model domain roughly covers a region of 500 km × 600 km (44.5 to 49.5 N and 4 to 13 E) at a horizontal grid spacing of 1.1 km × 1.1 km (Fig. 1). A height-based hybrid smoothed level vertical coordinate system (Schär et al.2002) with 80 levels is used and stretched from the surface to 22 km. For this study, we simulate the cold front passage between 11:00 and 18:00 UTC and analyze the results between 13:00 and 18:00 UTC on 26 March 2010. COSMO is forced with initial and hourly boundary conditions from reanalysis data at a horizontal resolution of 7 km × 7 km, supplied by MeteoSwiss. The model time step is 4 s with an output frequency every 15 min.

Figure 1Overview of the model orography and the instrument location setup. The parallelogram (dashed black lines) is the domain of the flow-oriented vertical cross-section analysis in Sect. 3.1 following the direction of the dual-polarized Doppler radar MXPol (red dot located at 46.789 N, 9.843 E) data. The blue box is the domain used for analysis in Sect. 3.2 and 3.3. The red triangle is the location of the ground-based video disdrometer.


Simulations were conducted including several SIP processes, which consisted of ice–graupel collisions (as thoroughly discussed in Sect. 2.2.2 below) and a control simulation, referred to as the rime-splintering (RS) simulation, where only rime splintering was active. For each of these setups, five ensemble simulations are conducted by perturbing the initial temperature conditions at each grid point through the model domain with unbiased Gaussian noise at a zero mean and a standard deviation of 0.01 C (Selz and Craig2015; Keil et al.2019). To account for the uncertainty associated with simulating atmospheric processes in mountainous terrain, three cross sections were interpolated from the model output (only the outer two cross sections are shown in Fig. 1). Each cross section, of which one cross section cuts more or less through the location of the radar, is separated by 1 km which is similar to the model resolution. The direction of each of the cross sections is similar to the direction of the generated RHI cross section from the weather radar. The three cross sections are then averaged and compared to the radar data. To generate the Hovmöller diagrams, we further took the mean along the length of the cross section for both the simulations and the radar data.

2.2.2 Cloud microphysics scheme

We use a detailed two-moment bulk cloud microphysics scheme within COSMO with six hydrometeor categories, including cloud droplets, rain, ice, snow, graupel and hail (Seifert and Beheng2006). The 2M scheme has been used extensively to study the evolution, lifetime, persistence and aerosol–cloud interactions of MPCs (Seifert et al.2006; Baldauf et al.2011; Lohmann et al.2016; Possner et al.2017; Henneberg2017; Glassmeier and Lohmann2018; Sullivan et al.2018; Eirund et al.2019, 2021). We refer to ice particles as any combination of the hail, graupel, snow or ice categories. Cloud droplet activation is based on an empirical activation spectrum which depends on the cloud-base vertical velocity and the prescribed number concentration of cloud condensation nuclei (Seifert and Beheng2006). The application is appropriate in atmospheric models with a horizontal grid size and time resolution of Δx≤1 km and Δt<10 s, respectively. The warm-phase autoconversion process from Seifert and Beheng (2001) was updated with the collision efficiencies from Pinsky et al. (2001) and also takes into account the decrease in terminal velocity associated with an increase in air density. A better approximation of the collision rate between hydrometeors was also introduced by Seifert and Beheng (2006), which makes use of the Wisner approximation (Wisner et al.1972).

Primary production of ice occurs via homogeneous and heterogeneous nucleation pathways. Homogeneous freezing of cloud droplets, parameterized from the homogeneous freezing rates of Cotton and Field (2002), is calculated for 0>T-50C. At −38C most cloud droplets will freeze given the enhanced homogeneous nucleation rates at colder temperatures. As a lower bound, the homogeneous freezing of all cloud droplets occurs at T=-50C. The homogeneous nucleation of solution droplets, typically associated with cirrus cloud formation, follows Kärcher et al. (2006). Here, the number density and size of nucleated ice crystals are determined by the vertical wind speed, temperature and pre-existing cloud ice. Heterogeneous nucleation is empirically derived, which depends on the chemical composition and surface area of multiple species of aerosols, namely, organics, soot and dust (Phillips et al.2008).

Table 1Sensitivity settings for the collisional-breakup (BR) parameterization. The conversion rate (conv) is the size (in µm) at which rimed ice crystals or snowflakes are converted to graupel, α is the scale factor, FBR is the fragments generated, and γBR is the decay rate of fragment number at warmer temperatures. When γBR=5, as used in the Takahashi parameterization (Sullivan et al.2018), then T is included in the simulation name. The bold numbers represent the sensitivity simulations pertaining to the different conversion sizes of 300, 400 and 500 µm.

Download Print Version | Download XLSX

Secondary ice production through rime splintering, which is widely used in numerical weather prediction models, is the only process that is included in the standard version of COSMO (Blyth and Latham1997; Ovtchinnikov and Kogan2000; Phillips et al.2006; Milbrandt and Morrison2016; Phillips et al.2017). In COSMO, rime splintering occurs at -3T-8C (Hallett and Mossop1974) when supercooled droplets and rain drops (Dc,r25µm) collide with ice hydrometeors (Di,s,g100µm) (e.g., Seifert and Beheng2006). A default value of 350 fragments per milligram of rime is used in the rime-splintering parameterization. Another SIP process, collisional breakup, was added to COSMO and tested in several studies (Sullivan et al.2018; Dedekind et al.2021). Generally, collisional breakup refers to the collision of any two frozen hydrometeors with different densities in which the collisional kinetic energy is sufficient that the collisional impact causes shattering. Here, collisional breakup is when either ice or snow particles collide with graupel and fracture. This can increase the number of ice particles at temperatures warmer than −21C. Almost no empirical constraint – apart from Takahashi et al. (1995), who used collisions between hail-sized particles – exists for the efficiency of any form of collisional breakup. The collisional kinetic energy and the density between hydrometeors are important parameters for this efficiency. In this study of the heavy snowfall event during which high Kdp values were recorded, we use the parameterizations for ice–graupel collisional breakup (BR) from Dedekind et al. (2021, BR28 and BR2.8T) and Sotiropoulou et al. (2021, BR-Sot) in COSMO in different forms:

(4)BR=FBRα(T-252)1.2exp-(T-252)/γBRfor BR28(FBR,α,γBR)=(280,10,2.5),(5)BR=FBRα(T-252)1.2exp-(T-252)/γBRfor BR2.8T(FBR,α,γBR)=(280,100,5),(6)BR=FBR(T-252)1.2exp-(T-252/γBR)DD0for BR-Sot(FBR,D0,γBR)=(280,0.02,5),

where BR is the number of fragments generated per collision, α is the scale factor, FBR is the leading coefficient, T is the temperature in kelvin, γBR is the decay rate of the fragment number at warmer temperatures, D is the diameter of particles undergoing fracturing and D0 is the diameter of the hail particles used in Takahashi et al. (1995) (Table 1). Because of the inconsistency between the hail particles and their corresponding fall velocity used in Takahashi et al. (1995), which is described in more detail in Dedekind et al. (2021), all the parameterizations (Eqs. 4, 5 and 6) have scaling factors. Equations (4) and (5) were applied in Dedekind et al. (2021) for the BR28 and BR2.8T simulations, respectively. Equation (4) is scaled by α=10 and has a slower decay rate of fragment number at warmer temperatures represented by γBR=2.5, and Eq. (5) is scaled by 100 while using the same decay rate of fragment numbers of γBR=5 as used in Sullivan et al. (2018), which was derived from Takahashi et al. (1995) (Table 1). Equation (6), for the BR-Sot simulation, was applied in Sotiropoulou et al. (2021). They used a scaling parameter, D/D0, that was applied to the breakup parameterization from Sullivan et al. (2018) where D0=0.02 m.

Figure 2Hovmöller diagrams of wind speed and vertical velocity for panels (a) Doppler radar, (b) RS, (c) BR28, (d) BR2.8T and (e) BR-Sot between 13:00 and 17:30 UTC. The filled contours denote wind blowing towards (blue and green) and away (white and brown) from the radar, respectively. The pink line is the −21C isotherm. At warmer temperatures collisional breakup occurs. Hatching denotes the region where the air layer is dynamically unstable, determined by a bulk Richardson number of less than 0.25.


Similar to Dedekind (2021), the ice crystal number concentration (ICNC) in COSMO is limited to 2000 L−1. Furthermore, Dedekind (2021) concluded that the conversion rate from ice crystals or snow to graupel, which is a function of the riming rate of ice crystals or snow with raindrops, may contribute to enhanced collisional breakup (Seifert et al.2006). In Eq. (70) of Seifert and Beheng (2006), they specify that ice and snow crystals can only be converted to graupel once they reach Di,s500µm. However, in the current version of the 2M scheme (as used in this study), ice and snow crystals are converted to graupel already once they exceed Di,s200µm. Therefore, earlier graupel formation is promoted in the current version, which should lead to enhanced SIP though ice–graupel collisions. To test the model's sensitivity to these different thresholds for graupel formation, we set up sensitivity studies with graupel formation at Di,s300, 400 and 500 µm to understand how the conversion rate impacts SIP processes. To accomplish this, we change the ice category conversion size requirement, Di,s, during riming from 200 µm (BR2.8T) to 300 µm (BR2.8T_300), 400 µm (BR2.8T_400) or lastly 500 µm (BR2.8T_500) (Table 1). The results of these sensitivity simulations can be found in Appendix C.

To investigate the impact of vertical wind shear and updraft on SIP, the probability density functions (PDFs) for the variables from the collisional-breakup simulations are analyzed. Furthermore, the joint PDFs are calculated along with the mutual information (MI, Shannon and Weaver1949) score, which quantifies the strengths of dependencies between the SIP rate and cloud properties (e.g., Dawe and Austin2013). For this purpose, a 10 km × 10 km region was selected and masked by the levels in which SIP occurred (T>-21C, blue box in Fig. 1) from 15:15 to 16:30 UTC. This resulted in the 16 121 data points for which an expression from Hacine-Gharbi et al. (2013) was used for finding the optimal number of bins (17 bins in our case) to estimate the MI for continuous random variables.

3 Results

3.1 The case study

A synoptic system passed over Switzerland on 26 March 2010, during which we analyzed the evolution of a cold front from 13:00 to 17:30 UTC with intense precipitation from 15:00 to 17:30 UTC. Furthermore, the cold front was associated with a surface temperature drop of ∼7C (Fig. S1 in the Supplement), a southwesterly wind flow at higher altitudes, vertical wind shear closer to the surface below 4 km a.m.s.l. and the development of peculiar polarimetric radar signatures. In particular, Kdp reached values around 1.5 km−1 at certain height levels, and towards the end of the event it was exceeding 2 km−1 (Fig. 3a). A statistical analysis of Kdp in snowfall conducted with this radar and in this location over a long observation period (Schneebeli et al.2013) showed that the 80th percentile of Kdp at every height level is lower than 0.5 km−1. Considering that the distribution of Kdp is very skewed, values above 1 km−1 in snow can be considered unusually large.

Figure 3Hovmöller diagrams of the (a) spectral differential phase (Kdp), the (b) differential reflectivity (ZDR), (c) the horizontal reflectivity (ZH) and (d) the hydrometeor class categories derived from the Doppler radar between 13:00 and 17:30 UTC. The gray and black lines in panels (a) and (b) are where both Kdp and ZDR are larger than 0.5 and 1, respectively.


Wind shear, observed by the dual-polarization Doppler radar, was visible between 2 and 5 km a.m.s.l. (Fig. 2a). The wind velocity at lower altitudes shifted from southerly to northerly, which was captured in all simulations, albeit not as prominent as the observations (Figs. 2b–e). The wind shear and the associated updrafts may have contributed to an enhanced SIP rate between 3 and 5 km in the collisional-breakup simulations (Fig. 2c–e). Here, the bulk Richardson number, which is a ratio of the buoyant energy to shear-kinetic energy, is determined to assess the dynamic stability of the air layer. An air layer becomes turbulent if the Richardson number is less than the critical Richardson number of 0.25 (e.g., Stull2016). Figure 2b–e shows where the air layer was turbulent with enhanced interactions between ice particles that could have caused enhanced ice–graupel collisions. Regions of enhanced updraft, where hydrometeors can grow to larger sizes, were mostly seen immediately above the turbulent layer.

The subzero temperature in the region of enhanced Kdp was warmer than −21C, which is in the temperature range favorable for ice–graupel collisions (Takahashi et al.1995). We thus hypothesize that the in situ cloud conditions together with the vertical wind shear could have triggered higher secondary ice production rates that can be reflected in radar measurements, as Kdp is a possible indicator of high number concentrations of oblate hydrometeors in the radar sampling volume (Kennedy and Rutledge2011; Bechini et al.2013; Grazioli et al.2015a; von Terzi et al.2022).

Figure 4Hovmöller diagrams of the simulated reflectively for panels (a) RS, (b) BR28, (c) BR2.8T and (d) BR-Sot between 13:00 and 17:30 UTC. The hatched area is defined as the MPC where the cloud droplet mass concentration and ice mass concentration are greater than 10 and 0.1 mg m−3, respectively. The pink line is the homogeneous freezing line at −38C, and the shaded gray area is the cloud area fraction.


Figure 5Vertical profiles of the (a) horizontal reflectivity (ZH), (b) specific differential phase (Kdp) and differential reflectivity (ZDR). The solid lines are the mean, with shaded areas and error bars showing the 10th and 90th percentiles for the model simulations and Doppler radar, respectively, at 15:30 UTC.


3.2 Simulated vs. observed radar reflectivity

3.2.1 Model and Doppler radar comparison

Horizontal reflectivity ZH is used to compare the model to the observations throughout the cloud and to analyze the impact of secondary ice production on the simulated radar reflectivity. During the early afternoon, the median of ZH remained mostly below 20 dBz. At around 15:15 UTC larger ice hydrometeors were present (either as a result of enhanced aggregation or depositional growth) between 4 and 6 km a.m.s.l., which then started to sediment (Fig. 3c and d). A peak in ZH at 3 km a.m.s.l. was observed in the fall streaks when the cloud droplets rimed onto the sedimenting ice hydrometeors. The RS simulation overestimated ZH by at least 8 dBz throughout the vertical profile at 15:30 UTC, while all the collisional-breakup simulations captured ZH more accurately, especially between 3 and 5 km a.m.s.l. (Figs. 4 and 5a). The radar-derived hydrometeor classification showed that much of the ice hydrometeor growth occurred through aggregation and riming. At 15:30 UTC, very high median Kdp>1 km−1 and ZDR>1 dB were observed (Fig. 5b). The vertical evolution of Kdp and ZDR is similar, with a peak observed about 4 km a.m.s.l., which is 1 km above the peak in ZH (Fig. 5b). The large and colocated values of ZDR and Kdp suggest that a large population of oblate and rimed particles, without a significant presence of large isotropic hydrometeors, was present.

Later, during the event (Fig. 3) the peak of ZDR was more often above the peak of Kdp, suggesting that the population of particles in the areas of enhanced Kdp also included larger isotropic aggregates. The occurrence of peaks in polarimetric variables at certain heights above ground (Kdp in particular) has been observed during intense snowfall events (e.g., Kennedy and Rutledge2011; Schneebeli et al.2013; Grazioli et al.2015a). The Kdp enhancement in particular has often been observed near the −15C isotherm and has been interpreted as the signature of enhanced dendritic growth (Kennedy and Rutledge2011; Bechini et al.2013) in combination with secondary ice production (von Terzi et al.2022). Dendrites are prone to aggregation; therefore, the Kdp peak disappears (and ZH increases) as particles approach the ground level.

Figure 6Vertical profiles for the (a) ice number concentration (NICE), (b) secondary ice production (SIP) and (c) ice water content (IWC). The solid lines are the mean, with shaded areas and error bars showing the 10th and 90th percentiles for the model simulations and Doppler radar, respectively, at 15:30 UTC. The green triangle is the 2DVD surface observations for hydrometeors D>0.2 mm.


Figure 7Particle size distribution for panels (a) ice, (b) snow, (c) graupel and (d) raindrops for all the simulations at 15:30 UTC.


Figure 8Hovmöller diagrams of graupel and rain mixing ratio for panels (a) RS, (b) BR28, (c) BR2.8T and (d) BR-Sot between 13:00 and 17:30 UTC. The hatched area is defined as the MPC where the cloud droplet mass concentration and ice mass concentration are greater than 10 and 0.1 mg m−3, respectively. The pink line is the homogeneous freezing line at −38C, and the shaded gray area is the cloud area fraction.


3.2.2 Modeled and measured cloud properties

Throughout the vertical profile below 6 km at 15:30 UTC, the ice crystal number concentration was at least an order of magnitude larger than expected from the RS simulation with a SIP rate in excess of 20 L−1 s−1 (Fig. 6). In both Figs. 6 and S2 the observed ice crystal number concentration recorded by the disdrometer (∼20 L−1) was remarkably well represented at the surface by the BR2.8T and BR-Sot simulations (similar results are shown in Dedekind et al.2021). The ice crystal and snow number concentrations were orders of magnitudes larger for Di<0.4 mm and Ds<0.8 mm, respectively, compared to the RS simulation (Fig. 7a and b).

During the late stage of the snowfall event, at 17:00 UTC, the replenishment of graupel diminished rapidly (Figs. 8 and S2c), causing a substantial reduction in the SIP rate (Fig. S2). Less collisional breakup allowed the ice and snow crystals to grow to larger sizes, Di1.2 mm and Ds3.3 mm, respectively, primarily through deposition and/or aggregation (Fig. S3a and b). A lower ZDR is consistent with less anisotropic particles produced by aggregation and/or riming. However, the enhanced concentration of oblate particles (increase in Kdp) was in contrast to the simulations, showing a reduction in cloud content as the cloud began to dissipate earlier than in the observations. None of the simulations were able to describe the high ice particle formation event that was most likely triggered through ice–ice collisions of aggregates of dendrites given the favorable temperature range.

3.2.3 Microphysical explanations

ZH was significantly overestimated by the RS simulation between 13:00 and 17:30 UTC, which most likely was a result of the following chain of events. (1) Insufficient droplets of size 25 µm (Fig. 7d), within the narrow temperature range (-3T-8C), led to a limitation in ice particle growth by riming and therefore limited rime splintering (Fig. 6b). (2) Because rime splintering was not that active, typical for wintertime MPCs (e.g., Henneberg et al.2017; Dedekind et al.2021), the ice and snow crystals grew mainly by depositional growth and aggregation. (3) The ice and snow crystal size distributions widened substantially (Figs. 7a, b and S3a, b). Both of these categories had number concentrations of less than 100 L−1 with particle diameters of up to 0.8 and 5.1 mm, respectively, at 15:30 UTC. (4) The larger ice and snow crystal diameters resulted in enhanced ZH. These observations are consistent with other times during the day which showed even larger-sized ice and snow crystals of 0.9 and 5.2 mm, respectively (Figs. S4a and S5a), as well as higher rain mass mixing ratios (e.g., Fig. 8a). There were single grid points where snowflakes even reached diameters of 13 to 17 mm during the latter part of the day (not shown here). Additionally, excessive size sorting in the model most likely contributed to the overestimation in ZH. Size sorting typically occurs within the sedimentation parameterization of 2M schemes in regions of vertical wind shear or updraft cores (Milbrandt and McTaggart-Cowan2010; Kumjian and Ryzhkov2012). All these factors contributed to the RS simulation overestimating ZH compared to the observations (Fig. 5a).

When collisional breakup was allowed to occur in the BR28, BR2.8T and BR-Sot simulations, the ice and snow crystals did not have time to grow as large compared to the ice and snow crystals in the RS simulation. The smaller ice particles were associated with a reduced ZH which compared better to the observations than the RS simulations.

Comparing the collisional-breakup simulations showed that the BR28 simulation still generated 8 times more ice particles than the BR2.8T and BR-Sot simulations at 4 km a.m.s.l. at T-15C (Fig. 6b). At temperatures of −10C (3 km a.m.s.l.), the SIP rate decreased rapidly from 100 L−1 s−1 to almost 0 L−1 s−1 at the surface. The lower SIP rates (less ice–graupel collisions) compared to the BR2.8T and BR-Sot simulations had several implications: (1) ice crystals and snow particles had more time to grow to larger sizes as seen in the wider particle size distributions (Fig. 7); (2) the number of ice hydrometeors was an order of magnitude below (worst in the collisional-breakup simulations) the observed ground-based video disdrometer observations of ∼20 L−1 (for particles less than 0.2 mm) at 15:30 and 17:00 UTC (Figs. 6a and S2c); and (3) interestingly, the simulated ZH compared slightly better with the radar observations, although the ice hydrometeors were overestimated (Figs. 5a, 6a, S2c and S2d). Because the BR2.8T and BR-Sot simulations showed similar results, which were better than the BR28 simulation, only the BR2.8T simulation is used in the next section to link the simulated wind shear to SIP.

3.3 Linking simulated wind shear to SIP

In this section, the dependence of the SIP rate on wind patterns is examined over the region (blue box) depicted in Fig. 1 during three time periods: 13:00–14:15 UTC (early, Fig. S6), 14:30–15:45 UTC (middle, Fig. 9) and 16:00–17:15 UTC (late, Fig. S7). The early period was categorized with strong wind shear and a turbulent layer below 3 km (Fig. 2a–e). During the middle period, the turbulent layer extended to 4 km during which graupel increased in the MPC. The late period was categorized by less wind shear, causing the dissipation of the turbulent layer. The cloud entered a glaciated state during this time. We analyzed the middle period as it was the most important period, according to the model, in terms of SIP and is therefore shown in Fig. 9. Regions in which SIP did not occur (e.g., T<-21C) were masked out for this analysis.

Figure 9Probability density functions of different variables (P(x)) from the BR2.8T simulation from all model levels (top row) and at each model level (bottom row) for (a, e) secondary ice production (SIP) rate, (b, f) wind shear, (c, g) U-wind shear and (d, h) V-wind shear between 14:30–15:45 UTC. The solid and dashed white lines are the horizontal 50th percentile and the 25th and 75th percentiles, respectively, of each variable over the 10 km×10 km domain.


Figure 10Joint probability density function multiplied by bin area (P(x,yxΔy) for the model output of SIP versus (a) U wind, (b) V wind, (c) updraft, (d) U-wind shear, (e) V-wind shear and (f) wind shear for the BR2.8T simulation. White lines are the 50th percentile as a function of the x-axis variable.


Figure 11Normalized median values of the PDFs of the variables in Table 2 for the three different time periods for the BR2.8T simulation.


A strong shift between the early and middle period in the V-wind median and interquartile range occurred from 20.5 to −0.7 m s−1 and 6.6 to 21.3 m s−1, respectively, compared to the U wind (Table 2). In fact, the U wind had a small variability between the early and middle periods (Table 2). As the afternoon progressed, the median of the strongest V-wind shear extended from near the surface at 14:30 to 5 km a.m.s.l. at 16:30 UTC (Figs. 2 and 9). Updraft cells developed above this level of strong vertical wind shear between 10 and 20 m s−1 km−1. Our observation is consistent with Houze and Medina (2005) and Medina and Houze (2015), who also showed that updraft cells occurred at times and locations where the shear was strongest (>10 m s−1 km−1). During the middle time period, the variability in the V-wind shear was the largest with an interquartile range of 16.3 m s−1 km−1, which coincided with the highest SIP rate (Fig. 9d and Table 2). The joint probability density functions (PDFs) (P(SIP rate, V-wind shear)) illustrate that the correlation median between the V-wind shear and SIP rate peaked at 9 m s−1 km−1 and 80 L−1 s−1 (Fig. 10e). This peak coincided with the region where the wind shifted from southwesterly to northerly (along the valley) between the early and middle periods. This shift was a result of the change in the V-wind speed from negative to positive at 2.9 km a.m.s.l. The joint PDF between the V wind and SIP (Fig. 10b) was highest at this altitude. Here, an environment of stronger meridional wind shear tends to coincide with the environment for high secondary ice production (Fig. 10f and Table 2).

The contribution of the updraft velocity to SIP is not as clear. The increase (early to middle period) and decrease (middle to late period) of the median of the SIP rates poorly co-varied with the U-wind shear or updraft compared to the V-wind shear (Fig. 11).

Table 2The 25th, 50th and 75th percentiles and the interquartile range (IQR) between the 10th and 90th percentiles of the vertical profiles for the BR2.8T simulations.

Download Print Version | Download XLSX

To further highlight this role, we calculated the mutual information shared between different sets of variables. The mutual information (MI: I(X;Y)) between variables X and Y was further analyzed for non-linear relationships (Shannon and Weaver1949), where X [SIP rate] and Y [U wind, V wind, updraft, U-wind shear, V-wind shear, wind shear]. I(X;Y)) of 0 bits means no information is shared between X and Y; therefore, Y cannot be inferred from X. (Further information about MI can be found in Appendix A and also in Dawe and Austin (2013).) For instance, during the last period the relationship of I(SIP rate; wind shear) weakened drastically to 0.021 bits to below the level of significance (Table 3) compared to earlier periods. This was expected as the diminishing cloud liquid water caused a reduction in the riming rates. Lower riming rates limited graupel formation which in turn reduced ice–graupel collisions.

The significantly higher MI score between I(SIP rate; wind shear) during the early and middle periods was a result of the strengthening of the northerly valley winds during the early afternoon hours when the predominant wind aloft was southwesterly (generating the dominant V-wind shear). The development of the northerly winds could have been a result of low-level blocking that generated the shear layer (Medina et al.2005). The sharp change in the wind speed and direction enhanced the turbulent overturning and therefore promoted the riming of ice crystals and snow, leading to the formation of graupel which in turn enhanced the SIP rates. The higher MI score between I(SIP rate; wind shear) implies there is a stronger relationship between wind shear and SIP than the other variables analyzed here.

Table 3Mutual information between SIP rate and wind properties of the vertical profiles for the BR2.8T simulations. The significance level is calculated by taking the maximum of 1000 Monte Carlo simulations of mutual information between a random permutation of SIP rates and each variable.

Download Print Version | Download XLSX

4 Conclusions

A cold front passage on 26 March 2010, over the Swiss Alps, associated with strong vertical wind shear and intense polarimetric signatures was observed with a dual-polarization Doppler weather radar deployed at Davos. This study investigates the role of vertical wind shear on the rate of SIP by making simulations of wintertime orographic MPCs with a non-hydrostatic, limited-area model, COSMO, which has a two-moment bulk microphysics scheme with six hydrometeor categories, as well as two additional parameterizations for ice–graupel collisions (e.g., Sotiropoulou et al.2020; Dedekind2021) based on Takahashi et al. (1995). To conclude, our main findings can be summarized as follows:

  • Large values of Kdp>1 km−1 suggest that a large population of oblate particles was present, most likely originating from ice–ice collisions, at 4 km a.m.s.l. This level coincided with the −15C isotherm. At −15C dendritic growth is very fast, causing low-density dendrites to fracture and aggregate. At this time, ZDR was also positive, indicating that large isotropic particles were less present. At lower altitudes, ZH increased while Kdp (and ZDR) decreased, suggesting aggregation and/or riming were occurring.

  • The rime-splintering simulations overestimated ZH throughout the vertical profile and underestimated the disdrometers’ number concentration of hydrometeors at the surface. Both shortcomings could be explained by omission of ice–graupel collisions.

  • The breakup simulations (BR28, BR2.8T and BR-Sot) caused narrower ice crystal and snow distributions (enhanced number concentrations of smaller ice particles), resulting in a better representation of ZH. The enhanced number concentrations of ice particles meant that these simulations, in particular BR2.8T and BR-Sot, captured the disdrometer observations of ∼20 L−1 (considering the 0.2 mm observation limit) at 15:30 and 17:00 UTC.

  • During the middle period, 14:30–15:45 UTC, the V-wind shear increased by 60 %, causing conditions favorable for accretion and leading to enhanced graupel formation and SIP in the region favorable for SIP.

  • Another time period with high Kdp but low ZDR was observed at 17:00 UTC, which was not captured by the breakup simulations as the graupel mixing ratio was depleted. The breakup parameterization does not include ice–ice collisions and relies only on graupel as the colliding species. At this time, the radar signatures suggested that collisions of dendrites caused the formation of small oblate particles (increasing Kdp) but also the formation of a few, larger, isotropic aggregates (decreasing ZDR).

  • The mutual information between the SIP rate and other variables like vertical wind shear and updraft velocity suggested that the SIP rate is best predicted by the overall wind shear.

  • The sensitivity of the ice–graupel simulations to the conversion rate size restriction was measured using the Kullback–Leibler divergence. The sensitivity simulations were not sensitive to the conversion rate size thresholds.

Turbulent overturning, whether it is associated with baroclinic waves (Gehring et al.2020) or low-level blocking (Medina et al.2005; Houze and Medina2005; Medina and Houze2015), has been shown to play an important role in accreting hydrometeors to form precipitation. Here, we considered that the interactions of ice hydrometeors can lead to ice–graupel collisions, causing enhanced small ice fragments, as opposed to only growing larger through aggregation. These smaller fragments fall slower against updraft and may decrease local precipitation rates, enhancing precipitation downstream of the flow (Dedekind et al.2021). Wind shear plays a significant role in ice–graupel collisions and may even be more important when all ice–ice collisions are considered in more physically robust collisional-breakup parameterizations (Yano and Phillips2011; Phillips et al.2017). By only considering ice–graupel collisions, we are limited to mainly investigating collisional breakup in MPCs where riming can occur to form graupel. In the case where a cloud becomes glaciated and graupel cannot form through riming, our parameterization will not be able to simulate SIP, which may still prove to be very important.

Appendix A: Derivation of the simulated radar reflectivity

Here we briefly show the calculation of the simulated ZH (Appendix B, Seifert2002). Calculating ZH from the two-moment cloud microphysics scheme would not be possible without approximations and assumptions. The following relationship for the radar reflectivity of drops (Zw), using the Rayleigh approximation for the cross section of drops (Eq. A1), results in

(A1) η w = π 5 K w 2 λ R 4 0 D 6 f w ( D ) d D ,

where D is the particle diameter, λR is the wavelength of radar radiation, ηw is the volumetric liquid water content, fw(D) is the number density distribution function for liquid water and Kw2=0.93 is the dielectric constant of liquid water. The reflectivity factor for cloud water is given by the following:

(A2) Z ̃ w = λ R 4 π 5 K w 2 η w = 0 D 6 f w ( D ) d D = 6 π ρ w 2 0 x 2 f w ( x ) d x = 6 π ρ w 2 Z w ,

where ρw is the water density. Because of the backscatter behavior for the mass-equivalent diameter with regards to ρw and ρi (ice density), the same applies to graupel, which is described as a spherical ice particle

(A3) η g = π 5 K i 2 λ R 4 6 π ρ i 2 0 x 2 f g ( x ) d x 1 ,

where x is the particle mass, fg(x) is the number density distribution function for graupel and Ki2=0.176 is the dielectric constant of ice. The radar reflectivity factor for ice particles (e.g., graupel) is given by the following:

(A4) Z ̃ g = λ R 4 π 5 K w 2 η g = K i 2 K w 2 6 π ρ i 2 0 x 2 f g ( x ) d x = K i 2 K w 2 6 π ρ i 2 Z g .

For melting ice particles, however, Kw must be used instead of Ki. In our study, the surface and in situ cloud temperatures were below 0 C. More information on the reflectivity calculations for melting ice particles can be found in Seifert (2002, Appendix B). Finally, the radar reflectivity factor is given by

(A5) dBZ = 10 ln 10 ln Z radar mm 6 m - 3 ,

where Zradar is the sum of the reflectivity calculated for each individual cloud particle category (e.g., cloud drops, raindrops, ice crystals, snow crystals, graupel and hail):

(A6) Z radar = 6 π ρ w 2 Z c + Z r + ρ w 2 ρ i 2 | K i 2 | | K w 2 | Z ic + Z s + Z g .
Appendix B: Mutual Information

The entropy H of the variable x's probability density function P(x) is defined by Shannon and Weaver (1949) to be

(B1) H = - P ( x ) ln ( P ( x ) ) d x ,

where x is the information content of a single measurement of P(x)=-lnP(x). The entropy is a measure of the amount of information that is required to represent the PDF. From here, both the Kullback–Leibler divergence and the mutual information can be calculated.

The Kullback–Leibler divergence, also known as the relative entropy, measures the distance between two probability distributions, P(x) and Q(x), over a discrete random variable X. The Kullback–Leibler divergence is defined as follows:

(B2) D KL ( P Q ) = P ( x ) ln P ( x ) Q ( x ) d x .

The mutual information (MI) is a measure of the mutual dependence between two random variables X and Y (e.g., the entropy of X subtracted from the entropy of X conditioned on Y):

(B3) I ( X ; Y ) = H ( X ) - H ( X | Y ) .

MI describes, therefore, how different the joint distribution of the pair (X, Y) is from the distribution of X and Y. Combining Eqs. (B1) and (B3) yields

(B4) I ( X ; Y ) = [ P ( x ) ln ( P ( x ) ) - P ( x , y ) ln ( P ( x | y ) ) ] d x d y ,

and because P(x|y)=P(x,y)/P(y), Eq. (B5) can be reduced to

(B5) I ( X ; Y ) = - P ( x , y ) ln P ( x , y ) P ( x ) P ( y ) d x d y .

The range of the MI is described as follows:

(B6) MI = 0 ,  if  P ( x , y ) = P ( x ) P ( y ) ( X  and  Y  are completely independent ) , H ( X ) ,  if  P ( x , y ) = P ( x ) = P ( y ) ( X  and  Y  are perfectly correlated ) .
Appendix C: SIP sensitivity to conversion rates

Figure C1Particle size distribution over the number concentrations for panels (a) ice, (b) snow, (c) graupel, (d) cloud droplets and (e) raindrops for all the sensitivity simulations at 15:30 UTC. BR2.8T, BR2.8T_300, BR2.8T_400 and BR2.8T_500 represents the size restriction of 200, 300, 400 and 500 µm, respectively, before ice crystals and snow can be converted to graupel.


In this section the sensitivity of SIP to the rate of graupel formation, which is dependent on ice or snow crystals being larger than a given size when riming occurs, is analyzed. Figure C1 shows the particle size distribution (PSD) for the sensitivity studies during which the size restrictions are modified, which could slow the conversion process of the ice crystals and snow particles to graupel. The PSD over the cross section at 15:30 UTC showed little difference in the ice crystal number concentrations where we expected higher ice crystal number concentration for BR2.8T and consequently higher snow number concentrations due to enhanced aggregation (Fig. C1a and b). The largest differences from the BR2.8T_300 simulation were in the form of enhanced snow (for diameters 0.14<Ds<0.42 mm) and graupel number concentrations (for diameters 1.2<Dg<2.2 mm). However, at 15:30 UTC there is no clear signal beyond model variability, showing that the slower conversion rates to graupel affect the simulations (Fig. C1b and c). We compared the probability distributions of the total number of ice hydrometeors (NISG, i.e., ice crystals, snow and graupel), SIP rate, ice crystal and graupel number concentrations of the BR2.8T_500, BR2.8T_400 and BR2.8T_300 simulations, respectively, to that of the reference simulation, BR2.8T, between 14:15 and 15:45 UTC (Fig. 8b). We chose this time period as the largest graupel concentrations were present. The Kullback–Leibler divergence (DKL(P ∥ Q)), which measures how one probability distribution P is different from a second probability distribution Q, shows little information loss between variables in Table C1, except for graupel. A value of 0 bits means that the probability distributions are the same (e.g., no information loss). The largest DKL(BR2.8T_300  BR2.8T) was 3.33 bits for the graupel distribution and was reflected in larger differences in the SIP rate of 1.07 bits and ice number concentration of 0.05 bits (Table C1 and Fig. C1c). If the SIP rate was sensitive to the conversion rate, it is expected that the information loss would be the greatest between BR2.8T_500 and BR2.8T and not between BR2.8T_300 and BR2.8T (e.g., DKL(BR2.8T_500  BR2.8T) >DKL(BR2.8T_400  BR2.8T) >DKL(BR2.8T_300  BR2.8T)). This result leads us to conclude that the different conversion rate size thresholds from ice crystals and snow to graupel, used in the paper in conjunction with the collisional-breakup parameterization, were not significant for SIP.

Table C1The Kullback–Leibler divergence, DKL(P ∥ Q), between two probability distributions P and Q from 14:14 to 15:45 UTC. P (BR2.8T_300, BR2.8T_400 and BR2.8T_500) is the measured probability distribution against the reference probability distribution Q (BR2.8T). Each distribution consist of ∼10 800 grid points separated into 24 bins over the cross section in Fig. 1.

Download Print Version | Download XLSX

Code and data availability

The COSMO model output, radar and 2DVD datasets used for our analysis are available at (Dedekind et al.2022a), and the software to analyze the data can be found at (Dedekind et al.2022b).


The supplement related to this article is available online at:

Author contributions

ZD conducted the simulations, and JG post-processed the radar data. ZD and JG were the main authors of the paper and interpreted the data. PHA and UL contributed to the interpretation and analysis techniques. ZD and JG contributed to the writing of the study.

Competing interests

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


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


We would like to thank Alexis Berne for his valuable feedback on the interpretation of the radar data. All simulations were performed with the Consortium for Small-scale Modeling (COSMO) model. The simulations were performed and are stored at the Swiss National Supercomputing Center (CSCS) under project s1009. We sincerely thank the reviewers for their constructive feedback. Their suggestions and comments considerably improved the quality of the article.

Financial support

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

Review statement

This paper was edited by Timothy Garrett and reviewed by two anonymous referees.


Andrić, J., Kumjian, M. R., Zrnić, D. S., Straka, J. M., and Melnikov, V. M.: Polarimetric Signatures above the Melting Layer in Winter Storms: An Observational and Modeling Study, J. Appl. Meteorol. Climatol., 52, 682–700,, 2013. a, b

Bader, M. J., Clough, S. A., and Cox, G. P.: Aircraft and dual polarization radar observations of hydrometeors in light stratiform precipitation, Q. J. Roy. Meteor. Soc., 113, 491–515,, 1987. a

Baldauf, M., Seifert, A., Förstner, J., Majewski, D., Raschendorfer, M., and Reinhardt, T.: Operational Convective-Scale Numerical Weather Prediction with the COSMO Model: Description and Sensitivities, Mon. Weather Rev., 139, 3887–3905,, 2011. a, b

Baumgartner, A. and Reichel, E.: The World Water Balance: Mean Annual Global, Continental and Maritime Precipitation, Evaporation and Run-off, Elsevier Scientific Publishing Company, 179 pp., ISBN 978-0-444-99858-3, 1975. a

Bechini, R., Baldini, L., and Chandrasekar, V.: Polarimetric Radar Observations in the Ice Region of Precipitating Clouds at C-Band and X-Band Radar Frequencies, J. Appl. Meteorol. Climatol., 52, 1147–1169,, 2013. a, b, c

Blyth, A. M. and Latham, J.: A multi-thermal model of cumulus glaciation via the hallett-Mossop process, Q. J. Roy. Meteor. Soc., 123, 1185–1198,, 1997. a

Cotton, R. J. and Field, P. R.: Ice nucleation characteristics of an isolated wave cloud, Q. J. Roy. Meteor. Soc., 128, 2417–2437,, 2002. a

Dawe, J. T. and Austin, P. H.: Direct entrainment and detrainment rate distributions of individual shallow cumulus clouds in an LES, Atmos. Chem. Phys., 13, 7795–7811,, 2013. a, b

Dedekind, Z.: The Impact of the Ice Phase on Orographic Mixed-phase Clouds and Surface Precipitation in the Swiss Alps, Doctoral Thesis, ETH Zurich,, 2021. a, b, c, d

Dedekind, Z., Lauber, A., Ferrachat, S., and Lohmann, U.: Sensitivity of precipitation formation to secondary ice production in winter orographic mixed-phase clouds, Atmos. Chem. Phys., 21, 15115–15134,, 2021. a, b, c, d, e, f, g, h, i, j, k, l

Dedekind, Z., Grazioli, J., Austin, P. H., and Lohmann, U.: Dataset for Heavy snowfall event over the Swiss Alps: Did wind shear impact secondary ice production? (Version 1), Zenodo [data set],, 2022a. a

Dedekind, Z., Grazioli, J., Austin, P. H., and Lohmann, U.: Software for Heavy snowfall event over the Swiss Alps: Did wind shear impact secondary ice production? (Version 1), Zenodo [code],, 2022b. a

Eirund, G. K., Possner, A., and Lohmann, U.: Response of Arctic mixed-phase clouds to aerosol perturbations under different surface forcings, Atmos. Chem. Phys., 19, 9847–9864,, 2019. a

Eirund, G. K., Drossaart van Dusseldorp, S., Brem, B. T., Dedekind, Z., Karrer, Y., Stoll, M., and Lohmann, U.: Aerosol–cloud–precipitation interactions during a Saharan dust event – A summertime case-study from the Alps, Q. J. Roy. Meteor. Soc.,, 2021. a, b

Field, P. R., Lawson, R. P., Brown, P. R. A., Lloyd, G., Westbrook, C., Moisseev, D., Miltenberger, A., Nenes, A., Blyth, A., Choularton, T., Connolly, P., Buehl, J., Crosier, J., Cui, Z., Dearden, C., DeMott, P., Flossmann, A., Heymsfield, A., Huang, Y., Kalesse, H., Kanji, Z. A., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Secondary Ice Production: Current State of the Science and Recommendations for the Future, Meteorological Monographs, 58, 7.1–7.20,, 2016. a

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

Georgakaki, P., Sotiropoulou, G., Vignon, É., Billault-Roux, A.-C., Berne, A., and Nenes, A.: Secondary ice production processes in wintertime alpine mixed-phase clouds, Atmos. Chem. Phys., 22, 1965–1988,, 2022. a, b

Glassmeier, F. and Lohmann, U.: Precipitation Susceptibility and Aerosol Buffering of Warm- and Mixed-Phase Orographic Clouds in Idealized Simulations, J. Atmos. Sci., 75, 1173–1194,, 2018. a

Grazioli, J., Schneebeli, M., and Berne, A.: Accuracy of Phase-Based Algorithms for the Estimation of the Specific Differential Phase Shift Using Simulated Polarimetric Weather Radar Data, IEEE Geosci. Remote Sens. Lett., 11, 763–767,, 2014a. a

Grazioli, J., Tuia, D., Monhart, S., Schneebeli, M., Raupach, T., and Berne, A.: Hydrometeor classification from two-dimensional video disdrometer data, Atmos. Meas. Tech., 7, 2869–2882,, 2014b. a

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

Grazioli, J., Tuia, D., and Berne, A.: Hydrometeor classification from polarimetric radar measurements: a clustering approach, Atmos. Meas. Tech., 8, 149–170,, 2015b. a

Grazioli, J., Genthon, C., Boudevillain, B., Duran-Alarcon, C., Del Guasta, M., Madeleine, J.-B., and Berne, A.: Measurements of precipitation in Dumont d'Urville, Adélie Land, East Antarctica, The Cryosphere, 11, 1797–1811,, 2017. a

Hacine-Gharbi, A., Deriche, M., Ravier, P., Harba, R., and Mohamadi, T.: A new histogram-based estimation technique of entropy and mutual information using mean squared error minimization, Comput. Electr. Eng., 39, 918–933,, 2013. a

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

Henneberg, O.: Orographic Mixed-phase Clouds in the Swiss Alps – Occurrence, Persistence and Sensitivity, Doctoral Thesis, ETH Zurich,, 2017. a

Henneberg, O., Henneberger, J., and Lohmann, U.: Formation and Development of Orographic Mixed-Phase Clouds, J. Atmos. Sci., 74, 3703–3724,, 2017. a, b, c

Heymsfield, A. J., Bansemer, A., Poellot, M. R., and Wood, N.: Observations of Ice Microphysics through the Melting Layer, J. Atmos. Sci., 72, 2902–2928,, 2015. a

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

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

Hogan, R. J., Field, P. R., Illingworth, A. J., Cotton, R. J., and Choularton, T. W.: 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

Houze, R. A. and Medina, S.: Turbulence as a Mechanism for Orographic Precipitation Enhancement, J. Atmos. Sci., 62, 3599–3623,, 2005. a, b, c, d, e

Kärcher, B., Hendricks, J., and Lohmann, U.: Physically based parameterization of cirrus cloud formation for use in global atmospheric models, J. Geophys. Res.-Atmos., 111, 1–11,, 2006. a

Keil, C., Baur, F., Bachmann, K., Rasp, S., Schneider, L., and Barthlott, C.: Relative contribution of soil moisture, boundary-layer and microphysical perturbations on convective predictability in different weather regimes, Q. J. Roy. Meteor. Soc., 145, 3102–3115,, 2019. a

Kennedy, P. C. and Rutledge, S. A.: S-Band Dual-Polarization Radar Observations of Winter Storms, J. Appl. Meteorol. Climatol., 50, 844–858, 2011. a, b, c

Kumjian, M. R. and Lombardo, K. A.: Insights into the Evolving Microphysical and Kinematic Structure of Northeastern U.S. Winter Storms from Dual-Polarization Doppler Radar, Mon. Weather Rev., 145, 1033–1061,, 2017. a, b

Kumjian, M. R. and Ryzhkov, A. V.: The Impact of Size Sorting on the Polarimetric Radar Variables, J. Atmos. Sci., 69, 2042–2060,, 2012. a

Kumjian, M. R., Rutledge, S. A., Rasmussen, R. M., Kennedy, P. C., and Dixon, M.: High-Resolution Polarimetric Radar Observations of Snow-Generating Cells, J. Appl. Meteorol. Climatol., 53, 1636–1658,, 2014. a, b

Lohmann, U., Henneberger, J., Henneberg, O., Fugal, J. P., Bühl, J., and Kanji, Z. A.: Persistence of orographic mixed-phase clouds, Geophys. Res. Lett., 43, 10512–10519,, 2016. a, 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. USA, 118, e2021387118,, 2021. a

Medina, S. and Houze, R. A.: Small-Scale Precipitation Elements in Midlatitude Cyclones Crossing the California Sierra Nevada, Mon. Weather Rev., 143, 2842–2870,, 2015. a, b, c, d, e

Medina, S., Smull, B. F., Houze, R. A., and Steiner, M.: Cross-Barrier Flow during Orographic Precipitation Events: Results from MAP and IMPROVE, J. Atmos. Sci., 62, 3580–3598,, 2005. a, b, c

Milbrandt, J. A. and McTaggart-Cowan, R.: Sedimentation-Induced Errors in Bulk Microphysics Schemes, J. Atmos. Sci., 67, 3931–3948,, 2010. a

Milbrandt, J. A. and Morrison, H.: Parameterization of Cloud Microphysics Based on the Prediction of Bulk Ice Particle Properties. Part III: Introduction of Multiple Free Categories, J. Atmos. Sci., 73, 975–995,, 2016. a

Mülmenstädt, J., Sourdeval, O., Delanoë, J., and Quaas, J.: Frequency of occurrence of rain from liquid-, mixed-, and ice-phase clouds derived from A-Train satellite retrievals, Geophys. Res. Lett., 42, 6502–6509,, 2015. a

Murphy, A. M., Ryzhkov, A., and Zhang, P.: Columnar Vertical Profile (CVP) Methodology for Validating Polarimetric Radar Retrievals in Ice Using In Situ Aircraft Measurements, J. Atmos. Ocean. Tech., 37, 1623–1642,, 2020. a, b, c

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. Climatol., 54, 1060–1068,, 2015. a

Oue, M., Kollias, P., Matrosov, S. Y., Battaglia, A., and Ryzhkov, A. V.: Analysis of the microphysical properties of snowfall using scanning polarimetric and vertically pointing multi-frequency Doppler radars, Atmos. Meas. Tech., 14, 4893–4913,, 2021. a

Ovtchinnikov, M. and Kogan, Y. L.: An Investigation of Ice Production Mechanisms in Small Cumuliform Clouds Using a 3D Model with Explicit Microphysics. Part I: Model Description, J. Atmos. Sci., 57, 2989–3003,<2989:AIOIPM>2.0.CO;2, 2000. a

Phillips, V. T. J., Blyth, A. M., Brown, P. R. A., Choularton, T. W., and Latham, J.: The glaciation of a cumulus cloud over New Mexico, Q. J. Roy. Meteor. Soc., 127, 1513–1534,, 2006. a

Phillips, V. T. J., DeMott, P. J., and Andronache, C.: An Empirical Parameterization of Heterogeneous Ice Nucleation for Multiple Chemical Species of Aerosol, J. Atmos. Sci., 65, 2757–2783,, 2008. a

Phillips, V. T. J., Yano, J.-I., Formenton, M., Ilotoviz, E., Kanawade, V., Kudzotsa, I., Sun, J., Bansemer, A., Detwiler, A. G., Khain, A., and Tessendorf, S. A.: Ice Multiplication by Breakup in Ice–Ice Collisions. Part II: Numerical Simulations, J. Atmos. Sci., 74, 2789–2811,, 2017. a, b, c

Pinsky, M., Khain, A., and Shapiro, M.: Collision Efficiency of Drops in a Wide Range of Reynolds Numbers: Effects of Pressure on Spectrum Evolution, J. Atmos. Sci., 58, 742–764,<0742:CEODIA>2.0.CO;2, 2001. a

Possner, A., Ekman, A. M. L., and Lohmann, U.: Cloud response and feedback processes in stratiform mixed-phase clouds perturbed by ship exhaust, Geophys. Res. Lett., 44, 1964–1972,, 2017. a

Ryzhkov, A. V. and Zrnic, D. S.: Radar Polarimetry for Weather Observations, Springer Atmospheric Sciences, Springer International Publishing, Cham,, 2019. a

Ryzhkov, A. V., Giangrande, S. E., Melnikov, V. M., and Schuur, T. J.: Calibration Issues of Dual-Polarization Radar Measurements, J. Atmos. Ocean. Tech., 22, 1138–1155,, 2005. a

Schär, C., Leuenberger, D., Fuhrer, O., Lüthi, D., and Girard, C.: A New Terrain-Following Vertical Coordinate Formulation for Atmospheric Prediction Models, Mon. Weather Rev., 130, 2459–2480,<2459:ANTFVC>2.0.CO;2, 2002. a

Schneebeli, M., Dawes, N., Lehning, M., and Berne, A.: High-Resolution Vertical Profiles of X-Band Polarimetric Radar Observables during Snowfall in the Swiss Alps, J. Appl. Meteorol. Climatol., 52, 378–394,, 2013. a, b, c, d, e, f, g

Seifert, A.: Parametrisierung wolkenmikrophysikalischer Prozesse und Simulation konvektiver Mischwolken, Doctoral Thesis, Universitat Karlsruhe, Forschungszentrum Karlsruhe, (last access: 22 September 2020), 2002. a, b

Seifert, A. and Beheng, K. D.: A double-moment parameterization for simulating autoconversion, accretion and selfcollection, Atmos. Res., 59-60, 265–281,, 2001. a

Seifert, A. and Beheng, K. D.: two-moment cloud microphysics parameterization for mixed-phase clouds. Part 2: Maritime vs. continental deep convective storms, Meteorol. Atmos. Phys., 92, 67–82,, 2006. a, b, c, d, e, f

Seifert, A., Khain, A., Pokrovsky, A., and Beheng, K. D.: A comparison of spectral bin and two-moment bulk mixed-phase cloud microphysics, Atmos. Res., 80, 46–66,, 2006. a, b

Selz, T. and Craig, G. C.: Upscale Error Growth in a High-Resolution Simulation of a Summertime Weather Event over Europe, Mon. Weather Rev., 143, 813–827,, 2015. a

Shannon, C. E. and Weaver, W.: The Mathematical Theory of Communication, University of Illinois Press, (last access: 31 March 2022), 1949. a, b, c

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–10,970,, 2016. a, b

Sotiropoulou, G., Sullivan, S., Savre, J., Lloyd, G., Lachlan-Cope, T., Ekman, A. M. L., and Nenes, A.: The impact of secondary ice production on Arctic stratocumulus, Atmos. Chem. Phys., 20, 1301–1316,, 2020. a, b, c

Sotiropoulou, G., Vignon, É., Young, G., Morrison, H., O'Shea, S. J., Lachlan-Cope, T., Berne, A., and Nenes, A.: Secondary ice production in summer clouds over the Antarctic coast: an underappreciated process in atmospheric models, Atmos. Chem. Phys., 21, 755–771,, 2021. a, b, c

Stull, R.: Practical Meteorology: An Algebra-based Survey of Atmospheric Science, AVP International, University of British Columbia, google-Books-ID: xP2sDAEACAAJ, 2016. 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, c, d, e, f, g

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

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

von Terzi, L., Dias Neto, J., Ori, D., Myagkov, A., and Kneifel, S.: Ice microphysical processes in the dendritic growth layer: a statistical analysis combining multi-frequency and polarimetric Doppler cloud radar observations, Atmos. Chem. Phys., 22, 11795–11821,, 2022. a, b, c

Wisner, C., Orville, H. D., and Myers, C.: A Numerical Model of a Hail-Bearing Cloud, J. Atmos. Sci., 29, 1160–1181,<1160:ANMOAH>2.0.CO;2, 1972.  a

Yano, J.-I. and Phillips, V. T. J.: Ice–Ice Collisions: An Ice Multiplication Process in Atmospheric Clouds, J. Atmos. Sci., 68, 322–333,, 2011. a, b

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–60, 343–359,, 2001. 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. a, b

Short summary
Simulations allowing ice particles to collide with one another producing more ice particles represented surface observations of ice particles accurately. An increase in ice particles formed through collisions was related to sharp changes in the wind direction and speed with height. Changes in wind speed and direction can therefore cause more enhanced collisions between ice particles and alter how fast and how much precipitation forms. Simulations were conducted with the atmospheric model COSMO.
Final-revised paper