Formation of ice particles in clouds at temperatures of -10∘C or warmer was documented by using ground-based radar observations. At these temperatures, the number concentration of ice-nucleating particles (INPs) is not only expected to be small, but this number is also highly uncertain. In addition, there are a number of studies reporting that the observed number concentration of ice particles exceeds expected INP concentrations, indicating that other ice generation mechanisms, such as secondary ice production (SIP), may play an important role in such clouds. To identify formation of ice crystals and report conditions in which they are generated, W-band cloud radar Doppler spectra observations collected at the Hyytiälä station for more than 2 years were used. Given that at these temperatures ice crystals grow mainly as columns, which have distinct linear depolarization ratio (LDR) values, the spectral LDR was utilized to identify newly formed ice particles.
It is found that in 5 %–13 % of clouds, where cloud top temperatures are -12∘C or warmer, production of columnar ice is detected. For colder clouds, this percentage can be as high as 33 %; 40 %–50 % of columnar-ice-producing events last less than 1 h, while 5 %–15 % can persist for more than 6 h. By comparing clouds where columnar crystals are produced and to the ones where these crystals are absent, the columnar-ice-producing clouds tend to have larger values of liquid water path and precipitation intensity. The columnar-ice-producing clouds were subdivided into three categories, using the temperature difference, ΔT, between the altitudes where columns are first detected and cloud top. The cases where ΔT is less than 2 K are typically single-layer shallow clouds where needles are produced at the cloud top. In multilayered clouds where 2 K<ΔT, columns are produced in a layer that is seeded by ice particles falling from above. This classification allows us to study potential impacts of various SIP mechanisms, such as the Hallet–Mossop process or freezing breakup, on columnar-ice production. To answer the question whether the observed ice particles are generated by SIP in the observed single-layer shallow clouds, ice particle number concentrations were retrieved and compared to several INP parameterizations. It was found that the ice number concentrations tend to be 1–3 orders of magnitude higher than the expected INP concentrations.
Introduction
Ice production in mixed-phase clouds is critical for their radiative and microphysical properties. At temperatures warmer than -38∘C, ice crystals form on ice-nucleating particles (INPs). In situ measurements have revealed that the number concentration of available INPs decreases with the increase in ambient temperature e.g.,. This dependence is more or less universal but can also be affected by other factors, such as the geographic location, air-mass types and aerosol compositions e.g.,. In addition, it has been found that INP concentrations at high latitudes are generally lower than at mid-latitudes e.g.,. Above -10∘C, the typical concentrations of INPs are below 10-1 L-1 and can be as low as 10-6 L-1. A number of studies, however, have reported that the ice number concentration in clouds with the top temperature warmer than -10∘C can exceed the expected concentration of INPs by several orders of magnitude e.g.,. This discrepancy implies that numerical weather prediction models that rely solely on INP parameterizations cannot realistically represent ice number concentrations in moderately to slightly supercooled clouds. As a result, the inappropriate parameterization of ice production may lead to biased estimates of surface shortwave radiation budget , among other things e.g.,.
Several mechanisms have been proposed to explain this discrepancy, such as the enhanced contact nucleation driven by the thermophoretic force during the evaporation of liquid drops , pre-activated INPs from evaporated ice particles nucleated above or secondary ice production (SIP) mechanisms see recent reviews by. The SIP has been studied by a number of laboratory experiments since the 1940s e.g.,. found that numerous ice splinters can be generated when supercooled liquid drops larger than ∼25µm are collected by large ice particles within the temperature range of -8 to -3∘C. This is referred to as the Hallett–-Mossop (H–M) process, the most studied and most frequently implemented SIP mechanism in numerical models despite more parameterizations for other SIP processes becoming available e.g.,. The enhanced ice number concentration can also be caused by the fragmentation of large supercooled liquid droplets e.g.,. It has been found that the secondary ice production efficiency is positively correlated with the size of liquid droplets and is enhanced in moist environments . At temperatures higher than -3∘C the fragmentation of drizzle is still active, as shown by field observations . In addition, studies using optical sensors mounted on aircrafts have reported the high concentration of ice columns within the temperature range of -10 to -3∘C e.g.,. Based on aircraft measurements from two field campaigns, concluded that the secondary ice process is highly associated with the presence of liquid droplets and aged rimed ice in turbulent regions. Recently, found that the ice concentration in tropical maritime stratiform clouds characterized by the top temperature above -8∘C is on the order of 10-1–101 L-1, which cannot be fully explained by primary ice nucleation, H–M processes or droplet freezing. However, despite the advantage of offering a direct way of interpreting ice microphysics, aircraft observations are only available from a few measurement campaigns and do not seem sufficient for a long-term assessment.
The polarimetric variables, such as differential reflectivity (Zdr), specific differential phase measurements (Kdp) and linear depolarization ratio (LDR), observed by dual-polarization radars are sensitive to the shape of hydrometeors and allow the analysis of ice particles with specific habits e.g.,. At temperatures of -10 to -2∘C, the depositional growth of an ice crystal is stronger at the basal faces than at the prism faces. Hence, the formation of columnar ice is preferred . This distinct habit can produce high Zdr and Kdp as observed by dual-polarization weather radars . For vertically pointing Ka- and W-band radars, ice columns usually produce LDR values as high as -15dB, which is distinctively higher than that of most other ice particle types . , , , and have shown that this strong LDR signal at the slow falling part of the radar Doppler spectrum can be used to identify columnar-ice crystals. In this study, this method is applied to long-term radar Doppler spectra observations for characterizing the production of columnar-ice particles in stratiform clouds. Similarly to , we show that this phenomenon is not uncommon. By comparing radar-based retrievals of ice number concentrations to INP parameterizations, one of which was derived from observations collected at our measurement site , we show that the ice number concentrations tend to be 1–3 orders of magnitude higher than the expected INP concentration. This also supports the conclusions reached by .
The paper is organized as follows. Section 2 introduces the data used in this study. The method for identifying columnar-ice particles from radar Doppler spectra is illustrated in Sect. 3. Statistical results are presented in Sect. 4. Section 5 compares the concentrations of columnar-ice particles and INPs in single-layer shallow clouds. Conclusions are given in Sect. 6.
DataCloud radar observations
The measurements used in this study were collected at Station for Measuring Ecosystem – Atmosphere Relations II located in Hyytiälä, southern Finland (61.845∘ N, 24.287∘ E; 150 m above mean sea level, a.m.s.l.). Since November 2017, a 94 GHz dual-polarization frequency-modulated continuous-wave Doppler cloud radar (HYytiälä Doppler RAdar, HYDRA-W) has been operating at the station. The radar is pointing vertically and measures radar signal spectral moments, LDR and dual-polarization Doppler spectra; see for the example of the data. The LDR decoupling is about 30 dB, so the minimum observable LDR is about -30 dB. The radar operates using three chirps that define range resolution, Doppler unambiguous velocity and spectral resolution at three range intervals. Between 102 and 996 m, the range resolution is 25.5 m, Doppler unambiguous velocity is 10.24 ms-1, and the Doppler spectral resolution is 0.02 ms-1. Between 996 and 3577 m, these values are 25.5 m, 5.12 ms-1 and 0.02 ms-1, respectively. For ranges above 3577 m, they are 34 m, 5.12 ms-1 and 0.02 ms-1, respectively. In this study, HYDRA-W observations recorded between February 2018 and April 2020 were utilized .
To remove noise from Doppler spectra observations, the spectral lines with the signal-to-noise ratio less than 5 dB were filtered out. Since both co-polar and cross-polar observations were used, i.e., to compute spectral LDR, this filtering could result in complete removal of the cross-polar signal, the power of which is typically 15–30 dB lower than that of the co-polar signal . In such cases, no LDR values were computed.
The antenna diameter of HYDRA-W is 0.5 m, which translates into the Fraunhofer far-field distance of 157 m and antenna beam width of 0.56∘. Therefore, the lowest radar range bin, which is not affected by the near-field effect, is 179 m (the fourth range bin). Data recorded at distance in the radar far-field were used in this study, therefore limiting the lowest altitude to 179 m where observations were taken.
In addition to the active remote sensing system, HYDRA-W is capable of estimating the liquid water path (LWP) by using the 89 GHz passive microwave channel observations. The brightness temperature at this band is regularly calibrated using liquid nitrogen. The site-specific relation between the measured 89 GHz brightness temperature and LWP was derived by the radar manufacturer using radio-sounding and reanalysis data.
Model temperature and humidity profiles
To obtain information on atmospheric state during the cloud observations, forecasts of the Deutscher Wetterdienst (DWD) operational global Icosahedral Nonhydrostatic (ICON) model were used. The microphysics scheme in ICON is inherited from the COSMO model . The ICON model output is provided over all Aerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS) cloud profiling stations (available at http://cloudnet.fmi.fi/, last access: 30 July 2020). The model data have hourly temporal resolution with a horizontal resolution of 13 km. The height resolution decreases with the increase in altitude; for example, the height resolution is 0.16 km at an altitude of 1 km. Its atmospheric products, such as temperature, relative humidity (RH) and pressure, over Hyytiälä were interpolated into the temporal and spatial resolutions of HYDRA-W .
Methods
The mean Doppler velocity (MDV) of hydrometeors observed by a vertically pointing radar is the combination of particle terminal fall velocities and the vertical component of air motion. While Doppler velocity alone could be used to identify certain types of particles , there are associated limitations. These limitations include uncertainties in hydrometeor classification due to similarities of terminal fall velocities of different particles , the presence of a mixture of ice particle populations within the radar volume and impact of air motion on the observed MDV .
By using radar Doppler spectra instead of MDV, contributions from different particle populations can be separated e.g.,. In radar Doppler spectral power observations, the presence of multiple populations of hydrometeors, such as the co-existence of supercooled liquid water and ice e.g.,, a mixture of different ice types e.g.,, could manifest as multiple spectral peaks. Even in such cases, however, classification of these particles can be ambiguous.
To further improve the hydrometeor identification, dual-polarization radar observations can be used. For slant measurements the spectral differential reflectivity has been found to be useful . Because hydrometeors typically do not have a preferred azimuth orientation, the differential reflectivity is not very useful for the classification purposes at vertical incidence. Nonetheless, LDR can be used to identify prolate particles , such as ice columns . At an elevation angle of 90∘, columnar-ice particles can produce LDR signals as high as -16 to -13 dB , while other ice particles may produce values smaller than -20 dB e.g.,. Furthermore, given the relatively small fall velocities of newly produced columns, in regions where they coexist with other ice particles they usually populate the slow-falling part of the Doppler spectrum. Therefore, the slow-falling part characterized by high spectral LDR (∼-15 dB) in the Doppler spectrum indicates the presence of columnar-ice particles .
HYDRA-W Doppler spectral (a) power and (b) LDR at 2 February 2018 07:56:38 UTC. (c) Spectral power (black line) and LDR (red crosses) at 1.124 km as marked by the black dot dashed lines in (a) and (b) and LDR (green dots) at 0.51 km as marked by the green dot dashed lines in (a) and (b). The gray, red and blue shaded areas in (c) indicate background ice falling from above, newly formed ice columns and supercooled liquid, respectively, at 1.124 km. Negative velocity indicates downward motion.
Figure shows an example of observed spectral power and LDR. Two distinct populations of ice particles can be clearly identified from Fig. a, and the slower-falling one corresponds to ice columns as indicated by the high spectral LDR (Fig. b). The observed spectral power and LDR at 1.124 km (black dot dashed lines) are shown in Fig. c. Three distinct peaks can be identified from the spectral power, and the slower-falling ice columns are well characterized by the spectral LDR exceeding -18 dB. In contrast, the spectral LDR of faster-falling ice is around -25 dB, which mainly depends on the cross-coupling between the polarization channels and can be much higher than the LDR signal of larger aggregates . Interestingly, supercooled liquid water also seems present, as indicated by the well-defined spectral peak at around 0 ms-1. It appears that this liquid layer extends from ∼0.9 to ∼1.3 km (Fig. a). The potential mechanisms of producing these ice columns will be discussed in more detail in the following sections.
Given the spectral characteristics of ice columns as discussed above, the following criteria were set to identify ice columns in clouds.
Within the slowest 1 ms-1 of the Doppler spectrum, at least three spectral bins exceed the LDR of -18 dB.
The temperature of the radar range bin is between -10 and 0 ∘C.
The observed radar Doppler spectrum is not only dependent on the scattering properties of hydrometeors in the radar volume, but is also affected by the turbulent broadening . For example, the air at around 0.51 km seems rather turbulent, as indicated by the spectral power (green dashed line in Fig. a). However, it appears that this issue does not significantly affect the results of columnar-ice detection. The noisy spectral LDR values (green dots in Fig. c) between 0.3 and 1 ms-1 are attributed to the low signal-to-noise ratio. Such weak impact on spectral LDR due to turbulence may be explained by the distinctively high LDR values of ice columns, which contrast with much weaker LDR signals of ice aggregates .
Results
By utilizing HYDRA-W Doppler spectral observations recorded between February 2018 and April 2020, statistics of environmental conditions associated with columnar-ice production were derived. All detected cloud cases within the temperature range of -10 to 0 ∘C were analyzed. From the selected events the cases where significant inversion was detected, which could cause melting e.g.,, were excluded. Given the data selection criteria, no rainfall or summer cloud cases were analyzed. This was done to avoid potential problems due to radar signal attenuation in rain and melting layer . In total, 175 d of observations satisfying the data selection criteria were identified and analyzed.
Temperature and RH conditions in columnar-ice-producing regions
Formation and growth of ice particles require favorable environmental conditions. These conditions were assessed by using ICON forecasts, which supplemented the radar observations. Here, we define Hcolumn_top as the highest level where ice columns are detected and Tcolumn_top as the temperature at this height. As shown in Fig. a, Tcolumn_top values are mostly between -8 and -3∘C, with the highest frequency at around -5∘C and a median value of -4.7∘C. These values are within the growth region of ice columns . Such temperature distribution also bears a good resemblance to the results obtained from the early rime-splintering laboratory experiment and a recent statistical study . The statistics of humidity relative to ice (RHice) and water (RHliquid) at Hcolumn_top are shown in Fig. b. The median values of RHice and RHliquid are 102.6 % and 98.3 %, respectively, indicating the supply of water vapor is sufficient for growth of ice particles. This finding is, however, not surprising, since the method detects ice columns, and they are growing in this temperature regime. However, the values of RHliquid and RHice should be interpreted with caution. ICON applies a liquid saturation adjustment, limiting the liquid supersaturation to saturation. RHliquid values exceeding 100 % are attributed to numerical artifacts. RHice was calculated based on the forecasted temperature, pressure as well as RHliquid and therefore can be affected by numerical artifacts as well. Given the uncertainty of ICON forecasts, we regard the presented statistics in Fig. as a “sanity” check for our method.
Statistics of (a) temperature (Tcolumn_top) and (b) relative humidity (RH) at Hcolumn_top for all identified columnar-ice-producing cases. Hcolumn_top: the highest level where ice columns are detected.
It should be noted that although ice columns can be detected by our method, Hcolumn_top may be lower than the height where ice crystals are generated. There are two potential reasons for this. Firstly, the newly formed ice particles may be less non-spherical , and in this case they will have LDR values which are much smaller than our detection threshold. Secondly, at early stages of growth the radar signal of ice crystals is rather weak and does not allow accurate detection and identification of columns . However, the altitude difference between Hcolumn_top and the actual height where columns are generated is expected to be small and not significantly affect our results.
Properties of columnar-ice-producing clouds
There are a number of questions that are associated with formation of ice crystals in clouds at these relatively warm temperatures. Above -10∘C, the number concentration of INPs is expected to be small and rather uncertain . Therefore, it is important to know how often and under which conditions these ice crystals form. Because ice formation could be facilitated by ice particles falling from upper cloud layers, i.e., by SIP, the location where ice columns are forming, with respect to the cloud top, should be identified. Finally, the importance of the columnar-ice production for surface precipitation should also be assessed.
Duration of cloud observations (bars) and occurrence of columnar-ice-producing clouds (red dotted curve) over Hyytiälä as a function of cloud top temperature (CTT). The results were calculated based on the data collected from February 2018 to April 2020.
To identify such conditions, cloud top temperature (CTT), defined as the temperature of the highest detected radar return for a given measurement time, is used. Because there are cases where several cloud layers are observed and there are gaps between these layers, typically the top of the lowest one is used. However, particles forming in the upper clouds, while not detected by the radar, may seed lower cloud layers and therefore modify their properties . To limit the impact of such conditions on our analysis, following we have used a radar echo gap of 2 km as the threshold, which defines whether the layers are connected. Recently, suggested that the threshold of 2 km may overestimate the cases of cloud seeding. For this reason, we have also tried the threshold of 0.5 km to determine the cloud top, but we did not see significant changes in the results.
The statistics of the recorded cloud top temperatures are presented in Fig. . The figure (left) shows durations of detected cloud samples within the temperature of -10 to 0 ∘C for a given CTT range as recorded during the observation period. Because of the focus on cold cloud cases, where the temperature in an atmospheric column does not exceed 0 ∘C, the observed cloud cases were recorded between October and April. The observations show that low-level clouds, i.e., clouds with warmer CTT, are relatively more frequent. This resembles the cloud occurrence statistics of, e.g., and . It appears that deeper clouds, i.e., where CTT is below -12∘C, are more conducive to columnar-ice production. In these cases the frequency of columnar-ice occurrence is about 25 %–33 %. For warmer clouds the frequency is lower and is around 5 %–13 %. The average occurrence is 15 %. Interestingly, our results are comparable with a recent study by , who found that the occurrence of columnar ice over an Arctic site is between 10 % and 25 % depending on the temperature.
Comparison of liquid water path (LWP) for clouds with and without columnar-ice production. A cloud sample was identified if the cloud base was within the temperature range of -10 to 0 ∘C. The box plots represent the median (horizontal strip) and 5 %–95 % quantile (whisker) ranges of the distribution.
As shown in Fig. b, the majority of columnar-ice production cases took place in areas around liquid saturation. Although direct observations of liquid were not available, the measurements collected by the 89 GHz passive channel in HYDRA-W allows estimation of LWP. The LWP values for the cloud cases are shown in Fig. . The observations show a significant amount of supercooled liquid water present in the atmospheric column. The cloud cases where ice columns were detected tend to have larger LWP values, especially where CTT values were smaller than -8∘C. This potentially indicates that supercooled liquid droplets may play important roles in formation of ice columns. Given the mixed-phase cloud conditions, the observed columns are most probably ice needles.
The same as Fig. 4 but for W-band radar reflectivity at the fourth range bin (179 m above the surface).
Formation of ice crystals is an efficient precipitation process . To evaluate the impact of columnar-ice production on precipitation, the radar equivalent reflectivity factor is used as the proxy for the precipitation intensity. For clouds where the radar echo extends to the ground, reflectivity values at the fourth radar gate, 179 m above ground
level (a.g.l.), were used. As shown in Fig. , the reflectivity increases with decreasing CTT. This is due to the link between the cloud thickness and precipitation intensity. The columnar-ice production tends to increase the precipitation intensity. This effect is more pronounced for warmer clouds, where CTT is -12∘C or warmer. In warmer clouds the precipitation intensity can be enhanced by as much as 10-fold. The factor of 10 increase in precipitation rate appears from the 10 dB increase in reflectivity . As will be discussed in the next section, warmer clouds tend to be single-layer clouds, where the crystal formation is more directly linked to precipitation formation. Colder clouds are prone to consist of the multiple cloud layers where precipitation processes are affected by multiple processes, such as riming, aggregation, and sublimation, at various levels e.g.,
Columnar-ice production in single-layer and multilayered clouds
For all detected columnar-ice-producing cases, the distribution of CTT was analyzed. As shown in Fig. a, ice columns can form in clouds with a wide range of CTT. The majority of cases fall in the CTT range of -20 to 0 ∘C, with two peaks at around -15 and -5∘C, respectively. The peak at around -15 agrees rather well with the high occurrence of ice columns in clouds, with a CTT of -20 to -12∘C (Fig. ).
Relative occurrence of (a) CTT and (b)ΔT for columnar-ice production cases. ΔT = Tcolumn_top-CTT.
Because processes responsible for the formation of ice particles in single-layer and multilayered clouds may be different, the classification of the cloud cases was performed. Using CTT and Tcolumn_top, we define ΔT as the temperature difference between them. The larger ΔT is, the lower the inside of the observed cloud system where the columns are formed. The relative occurrence of ΔT also shows two peaks as presented in Fig. b. Specifically, one peak is close to ΔT=0K, indicating that ice columns are generated close to the cloud top. The second ΔT peak is around 10 K.
Given the distinct distribution of ΔT, we have grouped the recorded clouds into the following three categories.
Type 1: ΔT≤ 2 K – columnar-ice production at cloud top
Type 2: 2 K<ΔT≤ 12 K – multilayered cloud
Type 3: 12 K<ΔT – multilayered cloud
Representative events of the above cloud types are presented below.
Columnar-ice production at cloud top: ΔT≤ 2 K
This type of cloud is usually single layer, and ice columns are generated close to the cloud top. Figure presents such an event on 4 November 2019. The precipitation intensity is relatively light, with the CTT at around -3∘C. The W-band reflectivity close to the surface increases to around 0 dB between 14:00 and 16:00 UTC. This region coincides with the enhanced LDR observations, which reaches values as high as -15 dB. Such high LDR values indicate that the dominant ice particle type during this period is columns. As shown in Fig. b, the cloud top is turbulent and seems to be capped by an inversion layer (Fig. c). The observed LWP ranges between 80 and 150 gm-2. This cloud with relatively low reflectivity persisted over Hyytiälä for about 1 d (not shown). Given the warm cloud top, the primary ice nucleation may not fully explain the significant columnar-ice production , as will be discussed in more detail later. Regarding the SIP, the H–M process does not seem to be active since it requires falling ice particles serving as rimers to produce ice splinters .
The columnar-ice-producing event on 4 November 2019. HYDRA-W observations of (a) equivalent reflectivity factor, (b) mean Doppler velocity, where negative velocity indicates downward motion, and (c) LDR. Panel (d) presents the (left axis) detected columnar-ice region and (right axis) LWP observed by HYDRA-W. The lines in (c) are isotherms produced by ICON.
As shown in Fig. , around 22 % of columnar-ice production cases are attributed to single-layer shallow clouds. also observed the prevalence of high LDR values for mixed-phase clouds, with a CTT of -5∘C. They speculated that these particles are formed mainly by primary ice nucleation instead of the SIP. Recently, reported a similar event of shallow stratiform clouds over the tropical ocean and found that neither INPs nor known SIP mechanisms can fully explain the strong production of ice particles in such clouds, with top temperatures warmer than -8∘C. In this study, we find that such clouds also frequently occur over Hyytiälä, and more detailed analysis will be presented in Sect. 5.
Columnar-ice production in multilayered clouds: 2 K<ΔT≤ 12 K
The event that took place on 13 February 2018 is representative of the second cloud type, as defined by ΔT. As shown in Fig. a, the precipitation intensity during this event is higher than during the discussed single-layer shallow cloud case. The cloud top temperature of the upper cloud layer is about -15 to -12∘C. Before 08:00 UTC, the observed LWP is close to 100 gm-2, and the mean Doppler velocity at around 0.8 km is relatively small (∼ 1 ms-1), which indicates that particles are unrimed or very lightly rimed. From 08:00 to 09:30 UTC, the falling snowflakes seem to be heavily rimed, as revealed by a rather high LWP (from 200 to over 400 gm-2) and mean Doppler velocity measurements (1.2–1.8 ms-1) . The high LWP period coincides with the region of ice columns (Fig. d). During this period, the observed LDR values are enhanced but still relatively small. This is due to masking of the needle LDR signal by larger snowflakes.
The columnar-ice-producing event recorded on 13 February 2018. HYDRA-W observations of (a) equivalent reflectivity factor, (b) mean Doppler velocity, where negative velocity indicates downward motion, and (c) LDR. Panel (d) presents the (left axis) detected columnar-ice region and (right axis) LWP observed by HYDRA-W. The lines in (a) and (c) are isotherms produced by ICON. Note that the y-axis scale in (a) is different from that in (b), (c) and (d).
This type of cloud frequently occurs over Hyytiälä (Figs. and ). Although sounding measurements were absent, this cloud type seems to very similar to the one reported by , namely, a layer of supercooled liquid water with the top temperature of around -15∘C seeding low-level stratus clouds in the boundary layer. In this event, the presence of supercooled liquid water may not be directly determined; however, the enhanced LWP values are indicative of the vigorously supercooled liquid water generation. In addition, the falling ice particles between 08:00 UTC and 09:30 UTC seem to be heavily rimed, as evident from mean Doppler velocity . The combination of the presence of supercooled liquid water and riming indicates that the H–M process could be taking place and could be responsible for the columnar-ice production.
$\,12\,{$\unit{{K}}$}}?>Deep multilayered clouds: ΔT> 12 K
The third cloud type is not very different from the second one and represents the tail of the observed ΔT distribution as shown in Fig. . This type of cloud system is a deeper precipitating system with a CTT of -60 to -40∘C; see Fig. for an example. The presented case took place on 2 February 2018. There are several features that are worthwhile pointing out. The mean Doppler velocity observations exhibit signatures of atmospheric waves. Between 05:00 and 06:00 UTC such waves can be observed at around 1 km altitude. At the later time, the wave signatures appear at 0.5 km a.g.l. The strongest velocity variation, observed around 06:00 UTC, seems to coincide with the LWP peak. At the same time as the waves are observed, signatures of columnar-ice production are also detected, pointing to a possible connection between the two.
Same as Fig. but for the columnar-ice-producing event observed on 2 February 2018. A zoom-in view on the wave signatures between 05:30 and 07:00 UTC is presented in (b).
Deep precipitating clouds usually have a large number of ice crystals formed at the cloud top. Given the large ice flux, manifested by the higher radar reflectivity values, in this precipitation system, it is difficult for supercooled liquid droplets to survive. The supercooled water droplets can be rapidly depleted through the Wegener–-Bergeron–-Findeisen (WBF) process as well as riming . Nevertheless, the atmospheric waves could generate conditions needed for forming and maintaining the presence of supercooled liquid water droplets . Recently, provided radar observational evidence showing that the Kelvin–Helmholtz instability is favorable for supercooled drizzle and secondary ice production. In such cases, ice needles may be generated by the H–M process or freezing breakup .
Characteristics of different columnar-ice-producing cloud types
As discussed above, we have identified three types of cloud systems where columnar-ice particles form. To better understand how columnar-ice production is related to cloud properties, persistence of columnar ice in these clouds and the amount of LWP were considered for further analysis.
Persistence of columnar-ice crystals
As was demonstrated by the case studies, the columnar-ice production may persist over several hours, and therefore these particles may play a major role in determining cloud properties. To document this, we have derived statistics of the columnar-ice production persistence for each cloud case. This was done by computing the duration of a continuous columnar-ice production event. Since in some cases, for example, in the presented shallow single-layer clouds (Fig. ) formation of columns may be intermittent in nature as can also be seen in, similar to , gaps of less than 30 min were accepted. In addition, cases persisting for less than 1 min were removed. As shown in Fig. , 40 %–50 % of columnar-ice-producing events persist for less than 1 h. However, there is a significant fraction that could last for more than 3 or even 6 h. This hints that the production of ice columns plays an important role in defining cloud properties and should be included while considering radiative or precipitation properties of such clouds (see, for example, Fig. ).
Relative occurrence of the persistence of columnar-ice-producing events from February 2018 to April 2020. ΔT=Tcolumn_top-CTT.
LWP
Presence of liquid water droplets may be a necessary condition for the formation of the observed ice columns. Compared to clouds where no ice columns were detected, columnar-ice-producing clouds have higher LWP values (Fig. ). For different columnar-ice-producing cloud types, the observed LWP values seem to be somewhat different. Figure shows LWP occurrence for these three cloud types. In general, their distributions are similar, while identifiable differences still exist. The median value of LWP for single-layer clouds is the lowest, and the second type of cloud has the highest LWP. Interestingly, while both the second and third types of clouds are multilayered, the LWP of the second type is detectably higher than for the third one. Comparing Figs. and , we find that the radar reflectivity above the columnar-ice layer is generally higher in the third cloud type than the second one. The INP concentrations in clouds with a cold top are expected to be higher than for warmer clouds, and hence more falling ice particles are expected for deep precipitating clouds. Therefore, we speculate that one of the reasons for the difference in LWP may be the ice number concentration, which is related to the consumption of the liquid water via the WBF process and riming .
Relative occurrence of LWP for all columnar-ice-producing events. LWP: liquid water path. ΔT=Tcolumn_top-CTT.
Potential role of SIP in columnar-ice production in single-layer shallow clouds
Given the rather frequent formation of ice crystals at temperatures warmer than -10∘C, where expected INP concentrations are low, it is important to investigate a potential role of SIP. In multilayered clouds, identified here as cloud types II and III, it has been found that the ice formation can be enhanced by the H–M process e.g.,, among other mechanisms . In single-layer clouds the mechanisms for ice multiplication are less established. However, recent studies by and have shown that freezing fragmentation may play such a role. To study whether these columnar-ice particles can be attributed to the SIP, the estimated ice number concentrations are compared to the expected INP concentrations. If the derived ice number concentrations exceed these of INPs, we can conclude that the SIP is potentially active in such clouds.
The concentration of INPs was computed by using CTT in the following different parameterizations.
parameterization based on INP measurements obtained below -10∘C.
parameterization which is not directly derived from INP measurements but the observed ice number concentrations when the impact of SIP is minimized. The temperature of measurements is between -30 and -5∘C.
parameterization based on INP measurements from nine sites between -35 and -9∘C. In our study we have used the average INP concentration–temperature relation presented in .
parameterization derived from the INP measurements obtained at Hyytiälä. The temperature range is -20 to -8∘C. Also from this study, the average INP concentration–temperature relation was used.
As was previously discussed, the INP parameterizations differ significantly at these temperatures (-10 to 0 ∘C), as shown in Fig. . It should be noted that not all parameterizations were derived using observations at these cloud temperatures, and some of them were extrapolated beyond their validity range. The most interesting comparison is to , which is based on observations at Hyytiälä collected during 2018, so their observation period at least partially overlaps with ours. It should also be pointed out that INP observations were carried out at the ground, where INP concentrations are typically higher.
Radar-based retrieval of particle number concentration is rather uncertain. Because of this uncertainty, the derived number concentrations should be treated as our best estimates of the order of magnitude of the ice column number concentration.
As shown in Fig. , the observed LWP for columnar-ice-producing clouds is significantly higher than those without ice columns. Hence, pristine ice crystals are anticipated to grow in mixed-phase conditions, and ice needles rather than solid columns are expected to form . This limits the parameter space, where we need to search for microphysical properties of ice particles to constrain our retrieval. The retrieval is based on estimating ice water content from radar observations, following , as
IWC=100.00058ZT+0.0923Z-0.00706T-0.992(gm-3),
where Z is the W-band radar reflectivity and T is the air temperature. Because in a single-layer shallow cloud (ΔT≤2∘C) ice needles are the predominant precipitation particles, radar reflectivity measured close to the ground, in the fourth radar gate, 329 m a.m.s.l. or 179 m a.g.l., is used in this retrieval. The selection of the reflectivity measured close to the ground helps to limit potential attenuation problems as well. Using IWC, the number concentration of ice needles, Nneedle, can be estimated as follows :
Nneedle=IWC103mneedle(L-1),
where mneedle is the mass of a characteristic ice needle. The introduced uncertainty at this step depends on the definition of the characteristic needle. Here, we use mean Doppler velocity (MDV) and the velocity–mass relation to estimate mneedle. Since MDV is reflectivity weighted, mneedle would be mainly determined by larger ice particles, and therefore the resulting Nneedle is underestimated. For the purpose of this study, this underestimation is not a major issue, because we want to test whether the observed Nneedle exceeds expected INP concentrations.
There are a number of reported ice needle properties. To take into account potential differences in ice needle properties, two relations linking velocity and mass by for rimed needles and for unrimed needles were used. For rimed ice needles, the relation between terminal fall velocity and mass was derived by and can be written as
vneedle, rimed=1.55(103mneedle, rimed)0.271⋅ρ(z)ρ(1024m)0.4(ms-1),
where mneedle, rimed is the mass of rimed ice needles, and the term ρ(z)ρ(1024m)0.4 accounts for the change in air density ρ at a given height z. In our study z is 329 m, and in the altitude where needles were observed is 1024 m. For unrimed needles, derived a relation linking terminal fall velocity and needle length, L. The terminal velocity of unrimed needles at a given height can be estimated as follows :
vneedle, unrimed=(0.0006+0.2796L-0.0497L2+0.0041L3)ρ(z)ρ(0m)0.4(ms-1).
The needle can be modeled as a cylinder, the mass of which is
mneedle, unrimed=10-3π4ρiceLD2(g),
where ρice and D denote density of unrimed needles and the minor axis, respectively. Their parameterizations have been given by :
ρice=0.6L-0.117(gcm-3)
and
D=0.1973L0.414(mm).
Applying the power law fit to mneedle, unrimed and vneedle, unrimed values when L∈[0.03,5] (mm) yields
vneedle, unrimed=1.09(103mneedle, unrimed)0.377⋅ρ(z)ρ(0m)0.4(ms-1).
The observed MDV is affected by vertical air motion. To at least partially mitigate this issue, the observed MDV is averaged over 20 min . While this step reduces the impact of air motion by averaging Doppler velocity over updrafts and downdrafts, the residual air motion is expected to widen the retrieved distribution of Nneedle.
The derived number concentrations of ice particles are compared to expected INP concentrations (Fig. ). The results show that the estimated ice number concentrations for rimed needles are generally larger than that of unrimed needles . Regardless of the difference between rimed and unrimed needles and the INP parameterization used, there seems to be a large fraction of cases where INP concentrations are not sufficient to explain observed Nneedle. Our results are similar to the conclusion reached by , who used a different approach for establishing the range of ice crystal concentration from radar observations. As shown in Fig. , the majority of Nneedle values fall in the range of 10-2–101L-1, which is similar to aircraft measurements obtained in tropical stratiform clouds .
Probability density function (PDF) of estimated columnar ice (left) and INP (right) concentrations in shallow single-layer clouds. INPs: ice-nucleating particles.
The significant discrepancy between INP concentrations estimated from INP parameterizations and retrieved ice number concentrations indicates that primary ice nucleation does not seem to be the only mechanism responsible for the formation of ice particles in these shallow clouds. Because the analysis was performed on shallow single-layer clouds, this discrepancy may not be explained by the H–M process since no rimers are falling from upper clouds. So it appears that other, less studied, SIP mechanisms may play an important role in amplifying ice number concentrations in such shallow clouds. This conclusion is in line with a number of other studies. For example, found that the SIP may take place at -5∘C in the absence of rimers, for which the cause is still under investigation. Recently, a similar finding was reported for stratiform clouds over the tropical ocean by . They have speculated that droplet collisional freezing and pre-activated INPs could be responsible for this discrepancy. Recent laboratory studies have shown that freezing breakup may be a source of secondary ice particles. have also suggested that freezing breakup may be more efficient than the H–M process in nature.
Conclusions
This study documents formation of ice particles in clouds at the temperatures of -10∘C or warmer. The analysis was performed using W-band cloud radar observations collected at the University of Helsinki Hyytiälä station starting from February 2018 through April 2020. The columnar-ice particles were identified using measurements of spectral LDR. It was found that columnar-ice formation is relatively frequent in clouds at temperatures of -10∘C or warmer. The occurrence frequency of columnar-ice particles is 5 %–13 % in clouds, with top temperatures exceeding -12∘C. In colder clouds, this percentage can be as high as 33 %. The columnar-ice-producing clouds tend to have a higher LWP, potentially indicating that supercooled water droplets are important for formation of the observed ice particles. It was also observed that columnar-ice production seems to have a significant impact on the surface precipitation. This effect is especially important for warmer clouds.
Using the temperature difference, ΔT, between the altitudes where columns are first detected and cloud top, the columnar-ice-producing clouds were subdivided into three categories. The first category, where ΔT is less than or equal to 2 K, represents shallow single-layer clouds. In these clouds ice particles are forming at or close to the cloud top. The other two categories, where 2 K<ΔT≤12K and ΔT>12K, represent deeper multilayered clouds. In multilayered cloud systems, columnar-ice crystals are forming in the lower cloud layer seeded by ice particles falling from upper cloud levels. It was observed that 40 %–50 % of columnar-ice production cases persist for 1 h or less, while in some cases they can persist for over 6 h. The distributions of LWP values for the three types of columnar-ice-producing clouds are somewhat different. The median LWP value is largest (221 gm-2) in clouds where 2 K<ΔT≤12K. Such high LWPs could favor riming and cause the Hallet–Mossop process. To draw a definite conclusion, however, a more thorough study, where locations of supercooled liquid layers are identified, is needed. For the single-layer shallow clouds, number concentrations of ice columns were derived from the radar observations. It was observed that the concentration of ice particles exceeds the expected concentration of INPs for a large number of cases. This indicates that a SIP mechanism is active in these clouds. Given that in single-layer shallow clouds there are no rimers that could cause the H–M process, we advocate that another SIP process may play a role here.
Code availability
The code used for processing radar data is available from GitHub https://github.com/HaoranLiHelsinki/MATLAB/tree/master/HYDRA_W (last access: 1 October 2021) and from Zenodo 10.5281/zenodo.5542678.
Data availability
W-band radar data are created by Dmitri Moisseev and are available from 10.23729/5aab6b78-90c9-49ab-8264-f4168528a0f3. ICON data are generated by the European Research Infrastructure for the observation of Aerosol, Clouds and Trace Gases (ACTRIS) and are available from the ACTRIS Data Centre using the following link: https://hdl.handle.net/21.12132/2.fdcadfd7b8ac4c4d.
Author contributions
DM conceptualized the study. HL performed the investigation and wrote the draft. OM, PT and DM contributed to reviewing and editing this paper.
Competing interests
Some of the authors are members of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors have also no other competing interests to declare.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Special issue statement
This article is part of the special issue “Ice nucleation in the boreal atmosphere”. It is not associated with a conference.
Acknowledgements
We would like to thank the personnel of Hyytiälä station for their support in field observation. Matti Leskinen is thanked for his work in radar maintenance. We want to thank Anniina Korpinen, Simo Tukiainen and Axel Seifert for helpful discussions on ICON forecasts.
We thank the European Research Infrastructure for the observation of ACTRIS for providing the DWD ICON (Global) model data, which were produced by the DWD and the Finnish Meteorological Institute and are available for download from https://cloudnet.fmi.fi/ (last access: 1 October 2021). ACTRIS has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement no. 654109.
Financial support
This research has been supported by the Horizon 2020 (grant nos. ACTRIS-2654109, ACTRIS PPP 739530, ACTRIS-IMP 871115, ATMO-ACCESS 101008004, and ERA-PLANET iCUPE 689443) and the Academy of Finland (ACTRIS-NF 328616, ACTRIS-CF 329274, NanoBioMass 307537, ACRoBEAR 334792, and Center of Excellence in Atmospheric Sciences, 307331, Atmosphere and Climate Competence Center, 337549), University of Helsinki (ACTRIS-HY).Open-access funding was provided by the Helsinki University Library.
Review statement
This paper was edited by Paul Zieger and reviewed by two anonymous referees.
References
Alkezweeny, A.: Freezing of supercooled water droplets due to collision,
J. Appl. Meteorol., 8, 994–995, 1969.
Aydin, K. and Walsh, T. M.: Millimeter wave scattering from spatial and planar
bullet rosettes, IEEE T. Geosci. Remote Sens., 37,
1138–1150, 1999.
Barthazy, E. and Schefold, R.: Fall velocity of snowflakes of different riming
degree and crystal types, Atmos. Res., 82, 391–398, 2006.
Beard, K. V.: Ice initiation in warm-base convective clouds: An assessment of
microphysical mechanisms, Atmos. Res., 28, 125–152, 1992.
Bringi, V. N. and Chandrasekar, V.: Polarimetric Doppler weather radar:
principles and applications, Cambridge University Press, 2001.Bühl, J., Seifert, P., Myagkov, A., and Ansmann, A.: Measuring ice- and liquid-water properties in mixed-phase cloud layers at the Leipzig Cloudnet station, Atmos. Chem. Phys., 16, 10609–10620, 10.5194/acp-16-10609-2016, 2016.CLU: Icon-iglo-12-23, icon-iglo-36-47, gdas1 model data; 2018-02-01 to 2020-04-30; from Hyytiälä, ACTRIS Data Centre [data set], https://hdl.handle.net/21.12132/2.fdcadfd7b8ac4c4d (last access: 1 October 2021), 2021.Cooper, W. A.: Ice initiation in natural clouds, Meteorol. Mon., 21, 29–32, 10.1175/0065-9401-21.43.29, 1986.
DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D.,
Twohy, C. H., Richardson, M., Eidhammer, T., and Rogers, D.: Predicting
global atmospheric ice nuclei distributions and their impacts on climate,
P. Natl. Acad. Sci. USA, 107, 11217–11222, 2010.
DeMott, P. J., Hill, T. C., McCluskey, C. S., Prather, K. A., Collins, D. B.,
Sullivan, R. C., Ruppel, M. J., Mason, R. H., Irish, V. E., Lee, T., et al.:
Sea spray aerosol as a unique source of ice nucleating particles, P. Natl. Acad. Sci. USA, 113, 5797–5803, 2016.
Dye, J. and Hobbs, P.: Effect of carbon dioxide on the shattering of freezing
water drops, Nature, 209, 464–466, 1966.
Evans, D. and Hutchinson, W.: The electrification of freezing water droplets
and of colliding ice particles, Q. J. Roy. Meteor.
Soc., 89, 370–375, 1963.Falconi, M. T., von Lerber, A., Ori, D., Marzano, F. S., and Moisseev, D.: Snowfall retrieval at X, Ka and W bands: consistency of backscattering and microphysical properties using BAECC ground-based measurements, Atmos. Meas. Tech., 11, 3059–3079, 10.5194/amt-11-3059-2018, 2018.
Field, P., Lawson, P., Brown, G., Lloyd, C., Westbrook, D., Moisseev, A., Miltenberger, A., Nenes, A., Blyth, A., Choularton, T., Connolly, P., Buehl, J., Crosier, J., Cui, Z., Dearden, C., DeMott, P., Flossmann, A., Heymsfield, A., Huang, Y., Kalesse, H., Kanji, Z., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Secondary
ice production: Current state of the science and recommendations for the
future, Meteorol. Monogr., 58, 7–1, 2017.
Findeisen, W. and Findeisen, E.: Untersuchungen uber die Eissplitterbildung an
Reifschichten (Ein Beitrag zur Frage der Entstehung der Gewitterelektrizitat
und zur Mikrostruktur der Cumulonimben), Meteorol. Z, 60, 145–154, 1943.
Fletcher, N. H.: The physics of rainclouds, Cambridge University Press, 1962.Fridlind, A. M., Ackerman, A., McFarquhar, G., Zhang, G., Poellot, M., DeMott,
P., Prenni, A., and Heymsfield, A.: Ice properties of single-layer
stratocumulus during the Mixed-Phase Arctic Cloud Experiment: 2. Model
results, J. Geophys. Res.-Atmos., 112, D24202, 10.1029/2007JD008646, 2007.Fukuta, N. and Takahashi, T.: The Growth of Atmospheric Ice Crystals: A Summary
of Findings in Vertical Supercooled Cloud Tunnel Studies, J.
Atmos. Sci., 56, 1963–1979,
10.1175/1520-0469(1999)056<1963:TGOAIC>2.0.CO;2, 1999.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, 10.5194/acp-20-7373-2020, 2020.
Giangrande, S. E., Toto, T., Bansemer, A., Kumjian, M. R., Mishra, S., and
Ryzhkov, A. V.: Insights into riming and aggregation processes as revealed by
aircraft, radar, and disdrometer observations for a 27 April 2011 widespread
precipitation event, J. Geophys. Res.-Atmos., 121,
5846–5863, 2016.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, 10.5194/acp-15-13787-2015, 2015.Hagihara, Y., Okamoto, H., and Yoshida, R.: Development of a combined
CloudSat-CALIPSO cloud mask to show global cloud distribution, J.
Geophys. Res.-Atmos., 115, D00H33
10.1029/2009JD012344, 2010.
Hallett, J. and Mossop, S.: Production of secondary ice particles during the
riming process, Nature, 249, 26–28, 1974.
Hari, P. and Kulmala, M.: Station for measuring ecosystem-atmosphere relations:
(SMEAR II), Boreal Environ. Res., 10, 315–322,
2005.
Heymsfield, A.: Ice crystal terminal velocities, J. Atmos.
Sci., 29, 1348–1357, 1972.
Heymsfield, A. J., Bansemer, A., and Twohy, C. H.: Refinements to ice particle
mass dimensional and terminal velocity relationships for ice clouds. Part I:
Temperature dependence, J. Atmos. Sci., 64, 1047–1067,
2007.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, 10.5194/gmd-11-4269-2018, 2018.
Hobbs, P.: The aggregation of ice particles in clouds and fogs at low
temperatures, J. Atmos. Sci., 22, 296–300, 1965.
Hobbs, P. V. and Rangno, A. L.: Ice particle concentrations in clouds, J. Atmos. Sci., 42, 2523–2549, 1985.
Hobbs, P. V. and Rangno, A. L.: Rapid development of high ice particle
concentrations in small polar maritime cumuliform clouds, J.
Atmos. Sci., 47, 2710–2722, 1990.
Hogan, R., Field, P., Illingworth, A., Cotton, R., and Choularton, T.:
Properties of embedded convection in warm-frontal mixed-phase cloud from
aircraft and polarimetric radar, Q. J. Roy.
Meteor. Soc., 128, 451–476, 2002.
Hogan, R. J., Mittermaier, M. P., and Illingworth, A. J.: The retrieval of ice
water content from radar reflectivity factor and temperature and its use in
evaluating a mesoscale model, J. Appl. Meteorol. Clim.,
45, 301–317, 2006.
Houser, J. L. and Bluestein, H. B.: Polarimetric Doppler radar observations of
Kelvin–Helmholtz waves in a winter storm, J. Atmos.
Sci., 68, 1676–1702, 2011.
Houze Jr., R. A. and Medina, S.: Turbulence as a mechanism for orographic
precipitation enhancement, J. Atmos. Sci., 62,
3599–3623, 2005.
Kajikawa, M.: Observation of falling motion of columnar snow crystals, J. Meteorol. Soc. Jpn. Ser. II, 54, 276–284, 1976.Kalesse, H., Vogl, T., Paduraru, C., and Luke, E.: Development and validation of a supervised machine learning radar Doppler spectra peak-finding algorithm, Atmos. Meas. Tech., 12, 4591–4617, 10.5194/amt-12-4591-2019, 2019.Kalesse, H., Szyrmer, W., Kneifel, S., Kollias, P., and Luke, E.: Fingerprints of a riming event on cloud radar Doppler spectra: observations and modeling, Atmos. Chem. Phys., 16, 2997–3012, 10.5194/acp-16-2997-2016, 2016.
Kanji, Z. A., Ladino, L. A., Wex, H., Boose, Y., Burkert-Kohn, M., Cziczo,
D. J., and Krämer, M.: Overview of ice nucleating particles,
Meteorological Monographs, 58, 1–1, 2017.Keinert, A., Spannagel, D., Leisner, T., and Kiselev, A.: Secondary Ice
Production upon Freezing of Freely Falling Drizzle Droplets, J.
Atmos. Sci., 77, 2959–2967, 10.1175/JAS-D-20-0081.1, 2020.
Keppas, S. C., Crosier, J., Choularton, T., and Bower, K.: Ice lollies: An ice
particle generated in supercooled conveyor belts, Geophys. Res.
Lett., 44, 5222–5230, 2017.
Kneifel, S. and Moisseev, D.: Long-Term Statistics of Riming in Nonconvective
Clouds Derived from Ground-Based Doppler Cloud Radar Observations, J.
Atmos. Sci., 77, 3495–3508, 2020.Knight, C. A.: Ice growth from the vapor at -5∘C, J. Atmos.
Sci., 69, 2031–2040, 2012.
Koenig, L. R.: The glaciating behavior of small cumulonimbus clouds, J.
Atmos. Sci., 20, 29–47, 1963.Kollias, P., Rémillard, J., Luke, E., and Szyrmer, W.: Cloud radar Doppler
spectra in drizzling stratiform clouds: 1. Forward modeling and remote
sensing applications, J. Geophys. Res.-Atmos., 116, D13201,
10.1029/2010jd015237,
2011.
Korolev, A. and Field, P. R.: The effect of dynamics on mixed-phase clouds:
Theoretical considerations, J. Atmos. Sci., 65, 66–86,
2008.Korolev, A. and Leisner, T.: Review of experimental studies of secondary ice production, Atmos. Chem. Phys., 20, 11767–11797, 10.5194/acp-20-11767-2020, 2020.
Korolev, A., McFarquhar, G., Field, P. R., Franklin, C., Lawson, P., Wang, Z., Williams, E., Abel, S. J., Axisa, D., Borrmann, S., Crosier, J., Fugal, J., Krämer, M., Lohmann, U., Schlenczek, O., Schnaiter, M., and Wendisch, M.: Mixed-phase
clouds: Progress and challenges, Meteorol. Monogr., 58, 5–1, 2017.Korolev, A., Heckman, I., Wolde, M., Ackerman, A. S., Fridlind, A. M., Ladino, L. A., Lawson, R. P., Milbrandt, J., and Williams, E.: A new look at the environmental conditions favorable to secondary ice production, Atmos. Chem. Phys., 20, 1391–1429, 10.5194/acp-20-1391-2020, 2020.
Korolev, A. V.: The influence of supersaturation fluctuations on droplet size
spectra formation, J. Atmos. Sci., 52, 3620–3634, 1995.
Küchler, N., Kneifel, S., Löhnert, U., Kollias, P., Czekala, H., and
Rose, T.: A W-band radar–radiometer system for accurate and continuous
monitoring of clouds and precipitation, J. Atmos. Ocean.
Tech., 34, 2375–2392, 2017.
Kumjian, M. R., Tobin, D. M., Oue, M., and Kollias, P.: Microphysical Insights
into Ice Pellet Formation Revealed by Fully Polarimetric Ka-Band Doppler
Radar, J. Appl. Meteorol. Clim., 59, 1557–1580, 2020.
Lamb, D. and Verlinde, J.: Physics and chemistry of clouds, Cambridge
University Press, 2011.
Lauber, A., Kiselev, A., Pander, T., Handmann, P., and Leisner, T.: Secondary
ice formation during freezing of levitated droplets, J.
Atmos. Sci., 75, 2815–2826, 2018.Lauber, A., Henneberger, J., Mignani, C., Ramelli, F., Pasquier, J. T., Wieder, J., Hervo, M., and Lohmann, U.: Continuous secondary-ice production initiated by updrafts through the melting layer in mountainous regions, Atmos. Chem. Phys., 21, 3855–3870, 10.5194/acp-21-3855-2021, 2021.Li, H.: Hyytiala W-band radar spectra data processing tool, Zenodo [code], 10.5281/zenodo.5542678, 2021.Li, H. and Moisseev, D.: Melting layer attenuation at Ka- and W-bands as
derived from multi-frequency radar Doppler spectra observations, J.
Geophys. Res.-Atmos., 124, 9520–9533, 10.1029/2019JD030316, 2019.Li, H. and Moisseev, D.: Two layers of melting ice particles within a single
radar bright band: Interpretation and implications, Geophys. Res.
Lett., 47, e2020GL087499, 10.1029/2020GL087499, 2020.
Li, H., Moisseev, D., and von Lerber, A.: How does riming affect
dual-polarization radar observations and snowflake shape?, J.
Geophys. Res.-Atmos., 123, 6070–6081, 2018.Li, H., Tiira, J., von Lerber, A., and Moisseev, D.: Towards the connection between snow microphysics and melting layer: insights from multifrequency and dual-polarization radar observations during BAECC, Atmos. Chem. Phys., 20, 9547–9562, 10.5194/acp-20-9547-2020, 2020.Li, H., Korolev, A., and Moisseev, D.: Supercooled liquid water and secondary ice production in Kelvin–Helmholtz instability as revealed by radar Doppler spectra observations, Atmos. Chem Phys., 21, 13593–13608, 10.5194/acp-21-13593-2021, 2021.
Locatelli, J. D. and Hobbs, P. V.: Fall speeds and masses of solid
precipitation particles, J. Geophys. Re., 79, 2185–2197,
1974.Luke, E. P., Kollias, P., and Shupe, M. D.: Detection of supercooled liquid in
mixed-phase clouds using radar Doppler spectra, J. Geophys.
Res.-Atmos., 115, D19201,
10.1029/2009JD012884, 2010.Luke, E. P., Yang, F., Kollias, P., Vogelmann, A. M., and Maahn, M.: New
insights into ice multiplication using remote-sensing observations of
slightly supercooled mixed-phase clouds in the Arctic, P.
Natl. Acad. Sci., 118, e2021387118, 10.1073/pnas.2021387118, 2021.Majewski, A. and French, J. R.: Supercooled drizzle development in response to semi-coherent vertical velocity fluctuations within an orographic-layer cloud, Atmos. Chem. Phys., 20, 5035–5054, 10.5194/acp-20-5035-2020, 2020.
Matrosov, S. Y., Reinking, R. F., Kropfli, R. A., Martner, B. E., and Bartram,
B.: On the use of radar depolarization ratios for estimating shapes of ice
hydrometeors in winter clouds, J. Appl. Meteorol., 40, 479–490,
2001.
McCluskey, C. S., Ovadnevaite, J., Rinaldi, M., Atkinson, J., Belosi, F., Ceburnis, D., Marullo, S., Hill, T. C. J., Lohmann, U., Kanji, Z. A., O'Dowd, C., Kreidenweis, S. M., and DeMott, P. J.:
Marine and terrestrial organic ice-nucleating particles in pristine marine to
continentally influenced Northeast Atlantic air masses, J.
Geophys. Res.-Atmos., 123, 6196–6212, 2018.Moisseev, D.: Hyytiala W-band cloud radar 2017–2019 (Version 1), Fairdata.fi [data set], 10.23729/5aab6b78-90c9-49ab-8264-f4168528a0f3, 2020.
Moisseev, D. N., Unal, C. M., Russchenberg, H. W., and Ligthart, L. P.:
Improved polarimetric calibration for atmospheric radars, J.
Atmos. Ocean. Tech., 19, 1968–1977, 2002.
Moisseev, D. N., Lautaportti, S., Tyynela, J., and Lim, S.: Dual-polarization
radar signatures in snowstorms: Role of snowflake aggregation, J.
Geophys. Res.-Atmos., 120, 12644–12655, 2015.
Mosimann, L.: An improved method for determining the degree of snow crystal
riming by vertical Doppler radar, Atmos. Res., 37, 305–323, 1995.
Mossop, S.: Concentrations of ice crystals in clouds, B. Am.
Meteorol. Soc., 51, 474–480, 1970.
Mossop, S.: The origin and concentration of ice crystals in clouds, B. Am. Meteorol. Soc., 66, 264–273, 1985.Myagkov, A., Seifert, P., Bauer-Pfundstein, M., and Wandinger, U.: Cloud radar with hybrid mode towards estimation of shape and orientation of ice crystals, Atmos. Meas. Tech., 9, 469–489, 10.5194/amt-9-469-2016, 2016a.Myagkov, A., Seifert, P., Wandinger, U., Bühl, J., and Engelmann, R.: Relationship between temperature and apparent shape of pristine ice crystals derived from polarimetric cloud radar observations during the ACCEPT campaign, Atmos. Meas. Tech., 9, 3739–3754, 10.5194/amt-9-3739-2016, 2016b.
Niemand, M., Möhler, O., Vogel, B., Vogel, H., Hoose, C., Connolly, P., Klein, H., Bingemer, H., Skrotzki, J., and Leisner, T.: A
particle-surface-area-based parameterization of immersion freezing on desert
dust particles, J. Atmos. Sci., 69, 3077–3092, 2012.
Oue, M., Kumjian, M. R., Lu, Y., Verlinde, J., Aydin, K., and Clothiaux, E. E.:
Linear depolarization ratios of columnar ice crystals in a deep precipitating
system over the Arctic observed by zenith-pointing Ka-band Doppler radar,
J. Appl. Meteorol. Clim., 54, 1060–1068, 2015.
Petäjä, T., O'Connor, E. J., Moisseev, D., Sinclair, V. A., Manninen, A. J., Väänänen, R., von Lerber, A., Thornton, J. A., Nicoll, K., Petersen, W., Chandrasekar, V., Smith, J. N., Winkler, P. M., Krüger, O., Hakola, H., Timonen, H., Brus, D., Laurila, T., Asmi, E., Riekkola, M. L., Mona, L., Massoli, P., Engelmann, R., Komppula, M., Wang, J., Kuang, C., Bäck, J., Virtanen, A., Levula, J., Ritsche, M., and Hickmon, N.: BAECC: A field campaign to elucidate the impact of
biogenic aerosols on clouds and climate, B. Am.
Meteorol. Soc., 97, 1909–1928, 2016.
Petters, M. and Wright, T.: Revisiting ice nucleation from precipitation
samples, Geophys. Res. Lett., 42, 8758–8766, 2015.
Prill, F., Reinert, D., Rieger, D., Zaengl, G., Schroeter, J., Foerstner, J.,
Werchner, S., Weimer, M., Ruhnke, R., and Vogel, B.: Working with the ICON
Model–Practical exercises for NWP mode and ICON-ART, ICON Model Tutorial,
Deutscher Wetterdienst, Offenbach, Germany, 2019.Proske, U., Bessenbacher, V., Dedekind, Z., Lohmann, U., and Neubauer, D.: How frequent is natural cloud seeding from ice cloud layers (<−35∘C) over Switzerland?, Atmos. Chem. Phys., 21, 5195–5216, 10.5194/acp-21-5195-2021, 2021.
Protat, A. and Williams, C. R.: The accuracy of radar estimates of ice terminal
fall speed from vertically pointing Doppler radar measurements, J.
Appl. Meteorol. Clim., 50, 2120–2138, 2011.Radenz, M., Bühl, J., Seifert, P., Griesche, H., and Engelmann, R.: peakTree: a framework for structure-preserving radar Doppler spectra analysis, Atmos. Meas. Tech., 12, 4813–4828, 10.5194/amt-12-4813-2019, 2019.
Rangno, A. L. and Hobbs, P. V.: Ice particles in stratiform clouds in the
Arctic and possible mechanisms for the production of high ice concentrations,
J. Geophys. Res.-Atmos., 106, 15065–15075, 2001.Reinert, D., Prill, F., Frank, H., Denhard, M., Baldauf, M., Schraff, C.,
Gebhardt, C., Marsigli, C., and Zängl, G.: DWD Database Reference for the
Global and Regional ICON and ICON-EPS Forecasting System, Tech.
rep., Deutscher Wetterdienst, available at: https://www.dwd.de/DWD/forschung/nwv/fepub/icon_database_main.pdf (last access: 1 October 2021),
2020.
Roberts, P. and Hallett, J.: A laboratory study of the ice nucleating
properties of some mineral particulates, Q. J. Roy.
Meteor. Soc., 94, 25–34, 1968.Schneider, J., Höhler, K., Heikkilä, P., Keskinen, J., Bertozzi, B., Bogert, P., Schorr, T., Umo, N. S., Vogel, F., Brasseur, Z., Wu, Y., Hakala, S., Duplissy, J., Moisseev, D., Kulmala, M., Adams, M. P., Murray, B. J., Korhonen, K., Hao, L., Thomson, E. S., Castarède, D., Leisner, T., Petäjä, T., and Möhler, O.: The seasonal cycle of ice-nucleating particles linked to the abundance of biogenic aerosol in boreal forests, Atmos. Chem. Phys., 21, 3899–3918, 10.5194/acp-21-3899-2021, 2021.
Scott, B. C. and Hobbs, P. V.: A theoretical study of the evolution of
mixed-phase cumulus clouds, J. Atmos. Sci., 34,
812–826, 1977.
Seifert, A.: A revised cloud microphysical parameterization for COSMO-LME,
COSMO Newsletter, 7, 25–28, 2008.
Seifert, P., Ansmann, A., Mattis, I., Althausen, D., and Tesche, M.: Lidar-Based Profiling of the Tropospheric Cloud-Ice Distributionto Study the Seeder-Feeder Mechanism and the Role of Saharan Dust as Ice Nuclei, in: Proceedings of the 8th International Symposium on Tropospheric Profiling, edited by: Apituley, A., Russchenberg, H. W. J., and Monna, W. A. A., p. 5, Delft, the Netherlands, Organizing Committee of the 8th International Symposium on Tropospheric Profiling, 2009.
Sekelsky, S. M.: Near-field reflectivity and antenna boresight gain corrections
for millimeter-wave atmospheric radars, J. Atmos. Ocean.
Tech., 19, 468–477, 2002.
Shupe, M. D., Kollias, P., Matrosov, S. Y., and Schneider, T. L.: Deriving
mixed-phase cloud properties from Doppler radar spectra, J.
Atmos. Ocean. Tech., 21, 660–670, 2004.
Shupe, M. D., Walden, V. P., Eloranta, E., Uttal, T., Campbell, J. R.,
Starkweather, S. M., and Shiobara, M.: Clouds at Arctic atmospheric
observatories. Part I: Occurrence and macrophysical properties, J.
Appl. Meteorol. Clim., 50, 626–644, 2011.Silber, I., Fridlind, A. M., Verlinde, J., Ackerman, A. S., Cesana, G. V., and Knopf, D. A.: The prevalence of precipitation from polar supercooled clouds, Atmos. Chem. Phys., 21, 3949–3971, 10.5194/acp-21-3949-2021, 2021.
Sinclair, V. A., Moisseev, D., and von Lerber, A.: How dual-polarization radar
observations can be used to verify model representation of secondary ice,
J. Geophys. Res.-Atmos., 121, 10–954, 2016.
Spek, A., Unal, C., Moisseev, D., Russchenberg, H., Chandrasekar, V., and
Dufournet, Y.: A new technique to categorize and retrieve the microphysical
properties of ice particles above the melting layer using radar
dual-polarization spectral analysis, J. Atmos. Ocean.
Tech., 25, 482–497, 2008.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, 10.5194/acp-18-16461-2018, 2018.
Sun, Z. and Shine, K. P.: Studies of the radiative properties of ice and
mixed-phase clouds, Q. J. Roy. Meteor. Soc.,
120, 111–137, 1994.
Tridon, F. and Battaglia, A.: Dual-frequency radar Doppler spectral retrieval
of rain drop size distributions and entangled dynamics variables, J.
Geophys. Res.-Atmos., 120, 5585–5601, 2015.
Tyynelä, J. and Chandrasekar, V.: Characterizing falling snow using
multifrequency dual-polarization measurements, J. Geophys.
Res.-Atmos., 119, 8268–8283, 2014.
Tyynelä, J., Leinonen, J., Moisseev, D., and Nousiainen, T.: Radar
backscattering from snowflakes: Comparison of fractal, aggregate, and soft
spheroid models, J. Atmos. Ocean. Tech., 28,
1365–1372, 2011.Vassel, M., Ickes, L., Maturilli, M., and Hoose, C.: Classification of Arctic multilayer clouds using radiosonde and radar data in Svalbard, Atmos. Chem. Phys., 19, 5111–5126, 10.5194/acp-19-5111-2019, 2019.
Verlinde, J., Rambukkange, M. P., Clothiaux, E. E., McFarquhar, G. M., and
Eloranta, E. W.: Arctic multilayered, mixed-phase cloud processes revealed in
millimeter-wave cloud radar Doppler spectra, J. Geophys. Res.-Atmos., 118, 13–199, 2013.
Westbrook, C. and Illingworth, A.: The formation of ice in a long-lived
supercooled layer cloud, Q. J. Roy. Meteor.
Soc., 139, 2209–2221, 2013.Wex, H., Huang, L., Zhang, W., Hung, H., Traversi, R., Becagli, S., Sheesley, R. J., Moffett, C. E., Barrett, T. E., Bossi, R., Skov, H., Hünerbein, A., Lubitz, J., Löffler, M., Linke, O., Hartmann, M., Herenz, P., and Stratmann, F.: Annual variability of ice-nucleating particle concentrations at different Arctic locations, Atmos. Chem. Phys., 19, 5293–5311, 10.5194/acp-19-5293-2019, 2019.Wildeman, S., Sterl, S., Sun, C., and Lohse, D.: Fast dynamics of water
droplets freezing from the outside in, Phys. Rev. Lett., 118,
084101, 10.1103/PhysRevLett.118.084101, 2017.
Wilson, T. W., Ladino, L. A., Alpert, P. A., Breckels, M. N., Brooks, I. M., Browse, J., Burrows, S. M., Carslaw, K. S., Huffman, J. A., Judd, C., Kilthau, W. P., Mason, R. H., McFiggans, G., Miller, L. A., Nájera, J. J., Polishchuk, E., Rae, S., Schiller, C. L., Si, M., Vergara Temprado, J., Whale, T. F., Wong, J. P. S., Wurl, O., Yakobi-Hancock, J. D., Abbatt, J. P. D., Aller, J. Y., Bertram, A. K., Knopf, D. A., and Murray, B. J.:
A marine biogenic source of atmospheric ice-nucleating particles, Nature,
525, 234–238, 2015.Yang, J., Wang, Z., Heymsfield, A. J., DeMott, P. J., Twohy, C. H., Suski,
K. J., and Toohey, D. W.: High ice concentration observed in tropical
maritime stratiform mixed-phase clouds with top temperatures warmer than -8∘C, Atmos. Res., 233, 104719, 10.1016/j.atmosres.2019.104719, 2020.
Young, G., Lachlan-Cope, T., O'Shea, S., Dearden, C., Listowski, C., Bower, K.,
Choularton, T., and Gallagher, M.: Radiative effects of secondary ice
enhancement in coastal Antarctic clouds, Geophys. Res. Lett., 46,
2312–2321, 2019.
Zängl, G., Reinert, D., Rípodas, P., and Baldauf, M.: The ICON
(ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M:
Description of the non-hydrostatic dynamical core, Q. J.
Roy. Meteor. Soc., 141, 563–579, 2015.
Zawadzki, I., Fabry, F., and Szyrmer, W.: Observations of supercooled water and
secondary ice generation by a vertically pointing X-band Doppler radar,
Atmos. Res., 59, 343–359, 2001.Zhao, X. and Liu, X.: Global Importance of Secondary Ice Production,
Geophys. Res. Lett., 48, e2021GL092581, 10.1029/2021GL092581, 2021.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, 10.5194/acp-21-5685-2021, 2021.bibl