Observing ice particle growth along fall streaks in mixed-phase clouds using spectral polarimetric radar data

. The growth of ice crystals in presence of supercooled liquid droplets represents the most important process for precipitation formation in the mid-latitudes. However, such mixed-phase interaction processes remain rela-tively unknown, as capturing the complexity in cloud dynamics and microphysical variabilities turns to be a real observational challenge. Ground-based radar systems equipped with fully polarimetric and Doppler capabilities in high temporal and spatial resolutions such as the S-band transportable atmospheric radar (TARA) are best suited to observe mixed-phase growth processes. In this paper, measurements are taken with the TARA radar during the ACCEPT campaign (analysis of the composition of clouds with extended polarization techniques). Besides the common radar observables, the 3-D wind ﬁeld is also retrieved due to TARA unique three beam conﬁguration. The novelty of this paper is to combine all these observations with a particle evolution detection algorithm based on a new fall streak retrieval technique in order to study ice particle growth within complex precipitating mixed-phased cloud systems. In the presented cases, three different growth processes of ice crystals, plate-like crystals, and needles are detected and related to the presence of supercooled liquid water. Moreover, TARA observed signatures are assessed with co-located measurements obtained from a cloud radar and radiosondes. This paper shows that it is possible to observe ice particle growth processes within complex systems taking advantage of adequate technology and state of the art retrieval algorithms. A signiﬁcant improve-ment is made towards a conclusive interpretation of ice particle growth processes and their contribution to rain production using fall streak rearranged radar data.

Abstract. The growth of ice crystals in presence of supercooled liquid droplets represents the most important process for precipitation formation in the mid-latitudes. However, such mixed-phase interaction processes remain relatively unknown, as capturing the complexity in cloud dynamics and microphysical variabilities turns to be a real observational challenge. Ground-based radar systems equipped with fully polarimetric and Doppler capabilities in high temporal and spatial resolutions such as the S-band transportable atmospheric radar (TARA) are best suited to observe mixedphase growth processes. In this paper, measurements are taken with the TARA radar during the ACCEPT campaign (analysis of the composition of clouds with extended polarization techniques). Besides the common radar observables, the 3-D wind field is also retrieved due to TARA unique three beam configuration. The novelty of this paper is to combine all these observations with a particle evolution detection algorithm based on a new fall streak retrieval technique in order to study ice particle growth within complex precipitating mixed-phased cloud systems. In the presented cases, three different growth processes of ice crystals, plate-like crystals, and needles are detected and related to the presence of supercooled liquid water. Moreover, TARA observed signatures are assessed with co-located measurements obtained from a cloud radar and radiosondes. This paper shows that it is possible to observe ice particle growth processes within complex systems taking advantage of adequate technology and state of the art retrieval algorithms. A significant improvement is made towards a conclusive interpretation of ice parti-cle growth processes and their contribution to rain production using fall streak rearranged radar data.

Introduction
The interaction of liquid water droplets and ice crystals leads to an enhanced growth of the ice crystals while falling through layers of supercooled liquid water. In the midlatitudes, this is an important process for precipitation formation (Pruppacher and Klett, 1996;Lamb and Verlinde, 2011). Nevertheless, the implications of the microphysical processes on precipitation formation are currently not well understood. Such microphysical processes involve a variety of hydrometeor sizes, shapes, phases, all affected by cloud dynamics, making the observation particularly challenging (Kollias et al., 2007;Shupe et al., 2008). Nowadays, ground-based radar measurement techniques have the advanced capabilities to observe and study microphysical processes within mixed-phase cloud systems (Kollias et al., 2007). Analyses of ice particle growth processes within mixed-phase cloud systems are primarily based on cloud radar data. Among them, Oue et al. (2016) uses polarimetric radar (K a -, and X-band) and lidar measurements to study ice particle growth processes in the cloud. Kalesse et al. (2016) analyzes Doppler spectra (K a -band) during a riming event and compare their findings of particle growth rates from the observations to simulations of such a riming event. Data from the same measurement campaign are used by Kneifel et al.
Published by Copernicus Publications on behalf of the European Geosciences Union.
(2015), where they apply a radar triple frequency method (W-, K a -, and X-band) to distinguish between different ice hydrometeor types (aggregates, rimed particles) and their particle size distribution characteristics during snowfall events. Myagkov et al. (2016) uses polarimetric radar measurements (K a -band) to derive the shape and orientation of mainly pristine ice particles generated at liquid topped mixed phase clouds.
In the studies mentioned above, the focus is on stratiform and layered clouds (dynamically stable and homogeneous) in the Arctic region or on snowfall events that do not involve any phase transition (ice particles constantly remain in the ice phase). Studies of precipitating mixed-phase cloud systems where melting of the ice particles is involved are rare. However, to understand the whole precipitation formation process the growth of the ice particles and the associated microphysical processes above the melting layer have to be studied to be related to the precipitation intensity below the melting layer. Nevertheless, the frequency bands listed in the studies above (W-, K a -, and X-band) suffer from attenuation when it comes to observations of precipitation cloud systems, Kollias et al. (2007). To avoid attenuation in precipitation, radars that operate in lower radar frequencies (Cand S-band) can be used. In the past, low-frequency radars were mainly used to estimate precipitation (Bringi and Chandrasekar, 2001) and to classify hydrometeor types (Straka et al., 2000;Chandrasekar et al., 2013). However, in the last few years, the polarimetric capabilities of those radars were also used to study ice particle microphysics and growth processes in cloud systems during snowfall events or in the ice part of the cloud above the melting (Kennedy and Rutledge, 2011;Bechini et al., 2013;Moisseev et al., 2015). Microphysical studies of ice crystal habits using a S-band radar were first done by Bader et al. (1987). They compared differential reflectivity, Z DR , signatures from the radar with aircraft measurements of ice crystal habits of precipitating stratiform clouds. It was found that the enhanced bands of Z DR values occur where the aircraft measurements detected defined pristine ice particle shapes. Later it was found that a combination of Z DR and specific differential propagation phase, K DP , can be used to identify regions of plate-like particle growth processes and onset of aggregation (Kennedy and Rutledge, 2011;Moisseev et al., 2015). Andrić et al. (2013) modeled polarimetric signatures for a winter storm case. The related simulations showed agreements with the observations; nevertheless, the model was not able to match all the features of the observations in a single run.
Low-frequency radar observation of the ice particle growth process due to the interaction of the ice crystals with supercooled liquid droplets is also found in the literature. Field et al. (2004) compared the observed Z DR signatures (Sband) to the obtained ice to liquid particle fraction measured with an aircraft and found a good agreement between liquid water presence and Z DR signatures that were caused by defined ice crystal shapes. Similar findings have also been made by Hogan et al. (2002). Hogan et al. (2002) showed that small convective cells embedded in warm frontal mixedphase cloud systems triggered ice particle formation, growth, and multiplication related to supercooled liquid water. Chandrasekar et al. (2013) discussed spectrograms measured with a C-band radar that exhibit signatures of aggregation and riming of ice particles. This shows that precipitation radar can capture particle growth in spectral domain and their related microphysics. Keat and Westbrook (2017) used Z DR signatures and co-polar correlation coefficient measurements within mixed phase clouds to separate the signals of pristine ice crystals and aggregates. Their results are validated against aircraft and K a -band radar observations.
In this work, radar data measured during the ACCEPT campaign are presented. They are obtained with TARA, which is a fully polarimetric S-band Doppler radar profiler providing measurements in high temporal and spatial resolutions. Because of TARA's unique three beam configuration, the full 3-D wind vector per sampling volume is also retrieved for each measurement. Therefore, using this radar and its configuration, it is possible to identify and study microphysical processes of ice particles within complex mixedphase clouds. On top of that, TARA measurements are rearranged along fall streaks according to the recently developed fall streak retrieval algorithm (Pfitzenmaier et al., 2017). Fall streaks are considered to be the precipitating path of a population of particles from top to bottom of the cloud system. Therefore, the analyses of microphysical changes along the fall streaks contain information of microphysical evolution of the same particle population (Pfitzenmaier et al., 2017;Kalesse et al., 2016). The retrieval technique reconstructs the precipitation path based on the obtained TARA wind information. In a next step, Doppler spectra and polarimetric radar variables rearranged along those fall streaks are used to study the microphysical evolution of a particle population. This paper demonstrates the advantage of analyzing data along fall streak instead of vertical profiles of radar data. Rearranged spectrograms of polarimetric variables allow to understand how the ice particle size and shape distribution changes versus height. This is used to improve the understanding of the impact of supercooled liquid water on ice particle growth within a cloud system. Furthermore, the spectral signatures are compared to additional measurements performed during the ACCEPT campaign (K a -band cloud radar spectral information and radiosonde temperature profiles) to assess the results from TARA data. After introducing the observation strategy of the ACCEPT campaign, in Sect. 2, some examples of microphysical observations of ice particle growth processes and their spectral signatures are given in Sects. 3 and 4. An overview of the fall streak technique is given in Sect. 5. Section 6 finally combines the retrieved microphysical information with fall streak correction to interpret growth processes of three different cases. Discussions and conclusions are provided in Sect. 7. The aim of the ACCEPT campaign (analysis of the composition of clouds with extended polarization techniques) is achieving a better understanding of microphysical processes involved in mixed-phase clouds using high-resolution polarimetric observations. The measurements were performed from October to mid-November 2014 at the Cabauw Experimental Site for Atmospheric Research (CESAR), the Netherlands. An instrumental synergy was used during ACCEPT to detect the different phases, the variety of sizes and shapes of the involved hydrometeors. For this purpose the TARA radar (Heijnen et al., 2000) at the site was extended by several other sensors: the Leipzig Aerosol and Cloud Remote Observation System (LACROS) (Bühl et al., 2013), a second vertically pointing K a -band cloud radar Mira (Görsdorf et al., 2015), and the Raman Lidars Polly XT Baars et al., 2016) and CAELI (CESAR Water Vapor, Aerosol and Cloud Lidar) (Apituley et al., 2009). A picture of the setup is provided in Fig. 1. During special observation periods, radiosondes were launched at the site.

Atmos
Some results of the ACCEPT campaign are already published. Myagkov et al. (2016) developed a retrieval to obtain the shapes of new generated ice crystals at liquid topped single layer mixed-phase clouds and compared their results to laboratory studies. Good agreement was obtained. Pfitzenmaier et al. (2017) developed a retrieval algorithm to reconstruct fall streaks within the TARA measurements to improve the study of ice particle growth due to supercooled liquid water presence within the cloud system.

The transportable atmospheric radar -TARA
TARA is a frequency modulated continuous wave (FM-CW) S-band radar profiler (3.3 GHz) that has full polarimetric and Doppler capabilities. TARA is able to provide highresolution column measurements in the Doppler domain (0.03 m s −1 ), in range (30.0 m, max. height of 10.05 km) as well as high temporal resolution (2.56 s). Measurements are performed under a fixed antenna elevation (45 • ) and a fixed azimuth (246.5 • related to the North). Measurements at 45 • elevation are considered as an optimum to get polarimetric and Doppler spectra information related to microphysical properties of the probed medium Unal et al., 2012). Therefore, TARA does not only allow the analysis of polarimetric bulk parameters as differential reflectivity, Z DR , and co-polar correlation coefficient, ρ HV , but it is also possible to study changes in the signatures of Doppler spectra, sZ (sZ hh ), and spectral differential reflectivity, sZ DR (sZ DR being the ratio of horizontal polarized to vertical polarized spectral reflectivity, sZ hh /sZ vv ). These changes in the spectra and spectrogram (height profiles of spectra) can be related to cloud microphysical variations. Therefore, it is possible to identify separate particle modes (having different velocity ranges) in the sZ and sZ DR spectra. Furthermore, the unique three beam configuration of TARA allows for the retrieval of the full 3-D wind vector per measurement. Using the Doppler spectra information of the three beams, main beam and two offset beams, the horizontal wind velocity vd h , the vertical Doppler velocity vd V , and the wind direction φ W are retrieved with a minimal temporal resolution of 2.56 s (Unal et al., 2012). The combination of all these measurements makes it possible to study ice particle growth within such complex precipitating mixed-phase cloud systems. Table 1 gives an overview of the specifications of TARA during ACCEPT.

Instrumental synergy
In addition to the TARA radar, other co-located sensors measured in parallel, allowing full comparison of the observed data signatures for similar volumes of study. In this paper, measurements of the Mira and TARA radars are compared, and data from two radiosonde launches are used. Unlike TARA, the data of high-frequency Mira type radar systems, suffer from attenuation and non-Rayleight effects during precipitation events. Low-frequency radar systems with high spatial and temporal resolution like the TARA radar have advantages measuring precipitation cloud systems. Nevertheless, the K a -band radar is more sensitive to the smaller pristine ice particles near cloud top. Therefore, it can detect particles in far ranges where TARA, operating at a low-frequency, is not able to measure these small particles. Despite of the attenuation and other effects the detected cloud top with Mira for the discussed case is at least 0.5 km higher than the one identified with TARA. TARA and the vertical pointing Mira obtain linear depolarization ratio in the spectral domain (sL DR ) at high temporal resolution (1 s for Mira). This parameter, which is the ratio of sZ vh , cross-polar received spectral reflectivity, to sZ hh , is not affected by attenuation issues at vertical incidence. Therefore, comparisons between the observed spectral polarimetric signatures from TARA and Mira can be done. These spectral polarimetric measurements are used in combination with thermodynamical conditions (temperature and humidity) within the cloud system, measured from radiosonde launches. On the one hand, radiosondes can retrieve regions of supercooled liquid water within precipitating cloud systems when other methods or sensors are limited. On the other hand, the temperature profiles in super-cooled liquid layer and the polarimetric radar signatures can help identifying the dominant ice particle shape affecting the growth process in the mixed-phase area (Bader et al., 1987;Pruppacher and Klett, 1996;Fukuta and Takahashi, 1999;Myagkov et al., 2016).

Ice particle growth processes
Ice particles within clouds can grow through three main processes (Pruppacher and Klett, 1996;Lamb and Verlinde, 2011): -Riming of ice particles. This occurs when supercooled water droplets collide and freeze onto bigger ice crystals, see Fig. 2a. Ice particles grow into large, dense and almost spherical shaped particles. In terms of conical graupel the resulting particles get slightly prolate Atmos. Chem. Phys., 18, 7843-7862, 2018 www.atmos-chem-phys.net/18/7843/2018/ (Oue et al., 2015a). Because of the sensitivity of the reflectivity to the hydrometeor size, Z hh ∼ D 6 N (D)dD where N (D) is the hydrometeor size distribution and D is the equivolumetric diameter, the growth of particles is strongly connected to an increase of reflectivity. The differential reflectivity Z DR gives information on the shape of the measured particle population (positive = oblate, ≈ 0 dB spherical, and negative = prolate particle shape). The conceptual change in particle shapes is indicated by the dotted lines in Fig. 2. Thus, riming processes affect Z DR values by turning oblate ice particle (Z DR > 0) into near spherical (Z DR around 0) or in case of conical graupel in slightly prolate particles (Z DR < 0).
-Water vapor diffusional growth of ice particles. This occurs when water vapor diffuses towards the crystals from the gas-phase, see Fig. 2b. During that process particles grow, therefore the reflectivity increases. During the diffusional growth particles keep their characteristic shape as long as they stay in the same growth regime, (Fukuta and Takahashi, 1999;Lamb and Verlinde, 2011). Their Z DR values slightly increase during the growth process. If particles grow large enough the crystals start to aggregate which leads to a decrease of polarimetric signature (Pruppacher and Klett, 1996). In our analysis we assume that the polarimetric signal is dominated by particle growth at constant atmospheric conditions and related to the shape of the dominant type of particles.
-Aggregation of ice particles. This occurs when ice crystals collide and form bigger ice crystals, see Fig. 2c. Aggregation leads to an increase of ice particle size and a change in particle shape from pristine particles to spherical shaped ice particles. The decrease of the differential reflectivity signatures depends strongly on the aggregated pristine particle type. Moisseev et al. (2015) suggested that early aggregates of dendrites align horizontally and therefore can also contribute to a high Z DR -signature. This leads to a less strong decrease of polarimetric signatures of Z DR with increasing size of the grown aggregates. Hobbs et al. (1974) indicated that needles have the tendency to clump very fast into spherical particles if their number concentration is high enough. Therefore, needles that aggregate turn faster into spherical shaped particles than dendrites.

Signatures of ice particle growth in spectral radar observations
The main three ice particle growth processes, riming, diffusional growth and aggregation result in different signatures in the spectral radar observations, see shows the diffusional growth of ice crystals and (c) the aggregation of ice crystals. The hexagonal plates and cylindrical columns represent either plate and dendrite shape or needles and columns, respectively. Size growth is indicated with bigger particles in the sketch. In the measurements an increase in particle size is related to an increase in Z values, indicated by the arrow at the bottom. The particle shapes are indicated by the dotted ellipses (the particles are modeled as spheroids). sZ, and differential reflectivity spectrum, sZ DR , is challenging due to concurrent growth process mechanisms occurring within the same resolution volume on different group of particles. Therefore, the signatures of different processes can be overlaid and difficult to separate. It is pointed out that these sketches are meant to explain spectral signatures of a S-band slantwise profiling radar. Other radar setups may have different spectral signatures. Figure 3a sketches riming signatures in the Doppler spectrum, sZ (black line), and differential reflectivity spectrum, sZ DR (red line). A Doppler spectrum represents the Doppler velocity distribution weighted by particle backscattering. Therefore, sZ is related to the particle size distribution (negative velocities indicate particle movements towards the radar, therefore, larger particles are associated with larger negative velocities), while sZ DR is related to the particle shape distribution. In a Doppler spectrum, a separated particle mode of larger and denser ice particles is often a clear indication of a riming process (Kollias et al., 2007;Chandrasekar et al., 2013;Oue et al., 2015a). As described in Fig. 3a they form a separated particle mode in the Doppler spectrum (left side) that indicates larger particles fall velocities in comparison to already existing particle mode. The particles in the mode  on the right side are growing. Ice crystals must have a certain size to be effective to rime (Wang and Ji, 2000;Ávila et al., 2009). Due to the fast growth of already large ice crystals the reflectivity values are also quite large in the right riming mode. Because not all ice crystals fulfill this criteria and might grow due to diffusional growth and/or aggregation, in the left mode the reflectivity also increases. Well defined shapes of pristine crystals are mainly responsible of the sZ DR -signatures in cloud. Riming of particles leads to a decrease of polarimetric signatures and a decrease of the sZ DR values towards zero. In case that riming leads to conical graupel sZ DR values become negative (see dotted red line at the left particle mode). In the right mode, a clear shapesize dependence is still expected and the higher sZ DR values correspond to the smaller ice crystals that have kept their pristine ice crystal shape. However, riming can also increase the sZ DR values. Oue et al. (2015a) showed that initial riming of the branches of dendrites increases the particle bulk density which leads to increased sZ DR values. Figure 3b depicts a signature of diffusional growth and aggregation in sZ and sZ DR . The small particle mode on the right side in the Doppler spectrum indicates the diffusional growth of the ice particles. Particles that grow via water vapor diffusion keep their pristine crystal shape and, therefore, the sZ DR -values of that mode stay large. A significant feature of the diffusional growth is the separated mode in the sZ DR spectrum. Such a mode of defined particle shapes is not present during the aggregation process of particles. However, at a certain stage particles start to aggregate and the Doppler spectrum width increases. The onset of aggregation is seen in the significant decrease of values in the sZ DR spectrum, because aggregation leads to more spherical shaped particles. Figure 3c illustrates an aggregation signature of ice crystals in sZ and sZ DR . Due to the domination of the large particles in the signal the spectral reflectivity can be skewed towards bigger particles. Depending on the particle size distribution a separation of the different particle modes is not always possible in the Doppler spectrum. Therefore, the signatures in sZ DR are needed to separate the different particle modes in the spectral domain if they are present. Nevertheless, it has to be mentioned that the differential reflectivity is also influenced by the density of the measured particles. Therefore, the strong decrease of sZ DR signature due to aggregation onset can be caused also by the lower density of the aggregates, compared to the pristine crystals.
The shown signatures in the spectral domain are used to identify qualitatively particle growth within mixed phase clouds. Using bulk and spectral polarimetric measurements, several cases during a complex precipitation system will be investigated next.

Methodology
The aim of this study is to analyze the microphysical process of a particle population on its way from cloud top to the bottom of the cloud system. Therefore, radar data are rearranged along so-called fall streaks before their signatures are examined. This offers a new insight on ice growth processes occurring in complex, local, and inhomogeneous cloud conditions in precipitating cloud systems. Marshall (1953) and Browne (1952) were the first that analyzed fall streak signatures within radar measurements. They investigated the structure and shape of outstanding reflectivity patterns within radar measurements. They described a fall streak like Yuter and Houze (2003), as a manifestation of an inhomogeneity in the microphysical structure of a cloud system. To be observed, the relative size and number of precipitation particles within the fall streak need to be large such that their radar reflectivity stands out as a local maximum from the immediate background reflectivity. Nevertheless, it was already pointed out that the main shape and structure of the fall streaks is influenced by the present cloud dynamics (Marshall, 1953;Browne, 1952). Because of this dynamical influence on the fall streak shape and structure, Pfitzenmaier et al.
(2017) defines a fall streak as the path of a particle population obtained from the observation of its own motion. Taking Atmos. Chem. Phys., 18, 7843-7862, 2018 www.atmos-chem-phys.net/18/7843/2018/ into account the dynamical conditions within the cloud system the fall streak retrieval is based on the mean 3-D wind field retrieved by TARA. Therefore, a mean fall streak is retrieved which is based on radar Doppler measurements only. Using directly measured 3-D wind information (vertical Doppler velocity vd V , horizontal wind speed vd h , and wind direction φ W ) the fall streak retrieval estimates the time displacement per height to reconstruct the path of a particle population. This is done using a bottom-up approach starting at the lowest (closest to the ground) valid data point at starting time t 0 . The time displacement at the height z is in principle estimated following Eq. (1). It consists of two terms: (1) The first term is the displacement time related to the antenna elevation, t α . The second term is the displacement time due to the cloud system dynamics, t dyn . The calculation of the displacement contributions t α and t dyn is done following the steps of Pfitzenmaier et al. (2017). The initial condition for the retrieval algorithm (cloud base height and averaging time for the wind profiles) are set individually per case, applying the suggested method in Pfitzenmaier et al. (2017). So, a mean fall streak per particle population is retrieved. Individual fall speeds of different particles within the population are not represented and taken into account and effects like size sorting cannot be analyzed using this method. However, it is possible to analyze particle growth processes along the particle fall streak under the assumption of continuous generation of new particles and microphysical homogeneity in the cloud system. Therefore, the rearranged data are supposed to represent all microphysical processes of the tracked particle population. The microphysical analyses of the cases are based on fall streak corrected bulk parameter profiles of reflectivity, Z, differential reflectivity, Z DR , and the co-polar correlation coefficient ρ HV as well as on the fall streak corrected Doppler spectra, sZ, and spectral differential reflectivity, sZ DR . The growth of the measured particle population is indicated by an increase of reflectivity. Changes in the Z DR -signature give information on the shape and density changes of the measured particle population. The co-polar correlation coefficient ρ HV gives information about particle homogeneity within the measured volume. The closer ρ HV is to 1 the more homogeneous are the particle shapes of the measured particle population. Therefore, changes of those parameters are used to examine the microphysical evolution of the particle population tracked from cloud top to bottom (Pfitzenmaier et al., 2017;Kalesse et al., 2016;Oue et al., 2015a). The analysis of the fall streak corrected spectrograms are used to investigate the changes in the spectra due to microphysical changes of the tracked ice particle distribution. Spectrograms provide Doppler spectra per height bin at each time step (spectral reflectivity versus Doppler velocity and height). Therefore, it is possible to identify the different signatures of riming, aggregation, and diffusional growth of different pristine crystals by analyzing their changes in signature (Fig. 3) from cloud top to bottom.
Due to the variation of small time scale dynamics (horizontal and vertical wind at 45 • elevation) on the spectra, the signatures are only analyzed qualitatively. Nevertheless, the changes in the Doppler spectra shape (broadening, modality, and amplitude) with height are detailed enough to provide information about the present particle microphysics. To ensure a good data quality the spectra are averaged over 3 time bins (∼ 8 s) per height (one time bin before and one time bin after the rearranged spectrum are used for the averaging) as well as a 10 dB threshold above the noise level is applied to Doppler bins in the spectra. This removes artifacts at the edges of the spectra due to low SNR. However sZ DR values are still noisy. Therefore, the presented and analyzed visible signatures in the spectrograms are also checked manually for consistency over time and height.

Observations and results
Figures 4 and 5 show measurements obtained with TARA radar, on 7 November 2014, from 10:00 to 12:00 UTC. In addition, the fields of reflectivity and linear depolarization ratio from the co-located vertical profiling cloud radar Mira are presented in Fig. 4a and d, respectively. The reflectivity (Z) fields in Fig. 4a and b display a precipitating cloud system. The band of enhanced Z values around 2.2 km indicates the melting layer of this frontal system. The attenuation and the non-Rayleigh scattering effect in the precipitation is clearly visible for the 35 GHz radar Mira, comparing the Z fields of both radars and, e.g., in Case 1 and Case 3 differences in Z-field above the melting layer up to 20 dBZ are visible. The Z fields in the cloud part show a high variation. Variability in the structure of the fields of differential reflectivity (Z DR ) and linear depolarization ratio (L DR ) are also visible in Fig. 4c and d. They have a high correlation with the visible fall streak signatures in the Z-field (e.g., 10:30-11:30 UTC above 4 km). Unfortunately, the TARA obtained wind fields between 10:30 and 11:20 UTC are corrupted (due to clutter in the non polarimetric offset beam measurements that influenced the retrieved wind fields especially the vertical Doppler velocity component, see Fig. 5c). Therefore, the data analysis is restricted to cases where the wind retrieval quality is good enough to apply the fall streak retrieval technique (Pfitzenmaier et al., 2017). Figure 5 displays the Doppler spectrum width and wind fields. Besides some small fluctuations in the horizontal wind field, Fig. 5b, the measurements show horizontally homogeneous wind conditions within the cloud in general. However, in the wind direction field, Fig. 5d, a vertical shear in wind direction is visible (about 30 • from 10:00 UTC onwards). Later the wind direction shear reaches the cloud part (starting at 11:35 UTC). . Overview of the radar measurements used in this paper. Panels (b, c) show reflectivity and differential reflectivity measurements, obtained with TARA radar (S-band). Panels (a, d) show the reflectivity and linear depolarization ratio measurements from the Mira radar (K a -band). The boxes 1, 2, and 3 highlight the time frames where fall streaks at 10:12:50, 11:33:31 and 11:47:30 UTC are retrieved and their rearranged data are discussed and analyzed Therefore, enhanced values of vertical Doppler velocity and horizontal wind speed are visible as well as an increase of the Doppler spectral width can be identified. The turbulence is caused by the mixing of an approaching air-mass that is related to a cold frontal cloud system, into the present airmass. Case 1 analyses an example during the enhanced Z DR and L DR signatures on top of the lower cloud layer between 10:09 and 10:18 UTC. While Case 2 focuses on the band of enhanced Z DR between 11:28 and 11:33 UTC at around 3 km. For the last case, Case 3, the fall streak structures of Z and polarimetric variables that generated within the enhanced band of Z DR at around 5 km between 11:35 and 11:47 UTC, are discussed.
6.1 Case 1: aggregation of needles, 10:09-10:18 UTC With Case 1 the analysis of along a fall streak rearranged Sband radar data is discussed to understand the present particle growth process. Figure 6a  quickly as the polarimetric variables change towards a spherical shape dependent signature. Therefore, a strong growth of particle size from pristine shaped particles into almost spherical particles is assumed to be present. To classify the particle population and the growth process more precisely, the fall streak corrected radar data are analyzed. The particle growth is visible in the fall streak corrected TARA sZ spectrogram in Fig. 6c. The sZ spectrogram shows a mono-modal particle population where sZ maximum increases from sZ = −7 dBZ at the top (3.3 km) to sZ = 10 dBZ above the melting layer (indicated by ML). The four Doppler spectra examples in Fig. 7 (left column, black spectra) show the evolution from cloud top towards the melting layer. In these examples of Doppler spectra the particle growth is seen by an increase of spectrum peak values and the broadening of the spectrum width with decreasing height. The fall streak corrected Z profile in Fig. 8a (black line) shows a 25 dBZ increase of Z, from cloud top towards the melting layer. In the profile a slight change in the slope at 3 km is visible. While above 3 km the generation and growth of particles are very fast ( Z ≈ 20 dBZ in 0.6 km) below the increase is lower. Below, the slope shows a slight increasing linear trend towards the melting layer. Such linearity in the reflectivity profile indicates aggregation of particles (Westbrook et al., 2007). Due to the homogeneous cloud conditions the contrast of the along fall streak rearranged and vertical Z profiles (light blue profile) is not huge in the cloud. A larger difference in the data is only visible in the reflectivity of the rain pattern. Therefore, using fall streak rearranged data in Figure 6. Panel (a) shows the retrieved fall streak at 10:12:50 UTC obtained with TARA. Panels (c, e) are the fall streak corrected spectrograms (sZ and sZ DR , 45 • elevation -note that the Doppler velocity contains the radial wind). Panels (b, d) show the fall streak corrected profiles of ρ HV and Z DR in black while light blue represents the vertical Z DR profile at 10:12:50 UTC. Panel (f) shows the sL DR spectrogram of the vertical pointing Mira at 10:13:00 UTC (all data displayed in the spectrograms have SNR > 10 dB). ML points out the signature of the melting layer and N the signature related to the needle particles.
Case 1 is mainly important for connecting the cloud microphysics to the rain intensity below.
The decrease of the polarimetric radar signatures in Fig. 6e and f at N (dotted circles), points out that the ice particles lose their pristine crystal shapes and become spherical particles. Values of sZ DR decrease in the spectrogram from sZ DR ≈ 1 dB near 3.25 km to sZ DR ≈ 0.2 dB around 3 km. In the spectrograms it is also visible that particles at all sizes become spherical below 3 km at the maximum value of sZ DR ≈ 0.2 dB. In comparison to the schematic sketch of the Z DR -spectrogram in Fig. 3c sZ DR values decrease throughout the whole size range (Doppler velocity range). In Fig. 7 this is visible in the column of the spectral differential reflectivity (right column, red spectra). At 3.076 km the sZ DR signature versus the whole Doppler velocity range in the spectra is close to zero dB. This observation indicates strong aggregation of particles into spherical ice particles. The fall streak corrected Z DR -profile in Fig. 6d (black line) shows the strong decrease of polarimetric signature from Z DR ≈ 1 dB to Z DR ≈ 0.1 dB at 3 km. The high values in the ρ HV profile in Fig. 6b exhibit homogeneity of the particle shapes within the measured population. ρ HV values are constantly above 0.9975 below 3 km and therefore match the low Z DR -values and sZ DR -signatures that confirm a homogeneous spherical particle population throughout the whole Doppler velocity range.
With TARA operating wavelength it is not possible to have a direct signature related to supercooled liquid water in the Doppler spectra or bulk parameters, additional information has to be used to detect its presence. The knowledge of the presence of supercooled liquid water helps in the characterization of the pristine crystal type, followed by an even better description of the involved ice crystal growth process. Comparing the observations with the temperature and dew-point temperature profiles measured by a radiosonde at 10:18 UTC (Fig. 9) a supercooled liquid water layer at cloud top is identified when both profiles match and the relative humidity reaches 100 % (light blue shaded area). The detected temperature range of −8 to −6 • C corresponds to the growth regime of needles at that height (Pruppacher and Klett, 1996)  a second indirect validation the observed sL DR values and signatures in Fig. 6f are used. They are consistent with other observations and simulated sL DR values of needle particles (Aydin and Walsh, 1999;Matrosov et al., 2001;Oue et al., 2015b). The simulated L DR values span a range from −16 to −12.5 dB that is in good agreement with the observed sL DR values in the spectrogram that reach values up to −17 dB in Fig. 6f. Therefore, the observed signatures of sL DR and the temperature values indicate the presence of needles at cloud top. Summing up needles that generated at the cloud top grow fast into a homogeneous population of spherical particles. The fall streak corrected Z-profile, Fig. 8a, and the sZspectrogram in Fig. 6c show an increase of reflectivity values. The slopes of the Z-profile are consistent with the linear increase of reflectivity in case of aggregation of ice crystals into snowflakes. In addition the observed Doppler ve-locities within Mira radar sL DR -spectrogram of −2 m s −1 , right above the melting layer also match the expected velocities for aggregates or slightly rimed particles (Mitchell, 1996;Fukuta and Takahashi, 1999). A strong aggregation or clumping of needle particles is mentioned by Hobbs et al. (1974), Rangno and Hobbs (2001) and Hogan et al. (2002). The example of sZ DR -spectra in Fig. 7 shows that such spherical particles can be observed within a short distance below particle generation. These compact and dense particles are more efficient to produce precipitation compared to the ice crystals and snowflakes produced in the time frame before Case 1.
6.2 Case 2: generation of a second particle population, 11:28-11:33 UTC The focus of Case 2 is the enhanced Z DR band between 3.2 and 2.6 km in Fig. 4c and the related microphysical processes that can be identified and connected to its presence. Figure 4b shows that in parallel the reflectivity values increase below 3.2 km, compared to above, which indicates particle growth. Because the band of enhanced Z DR decreases (particles become more spherical) towards the melting layer a particle growth process can be assumed in Case 2.
In order to go further in the process description, rearranged fall streaks are additionally retrieved and analyzed. The retrieval is done setting the initial conditions to 2.25 km for the cloud base height and 90 s for the wind averaging. The Z DR profile, Fig. 10d (black), shows at N, a strong increase from Z DR ≈ 0.2 dB at 3.2 km to a maximum of 0.65 dB at 2.9 km which indicates the generation of new particles at that height. From N on towards the melting layer, the Z profile, Fig. 8b (black), shows a linear increase of reflectivity of 15 dBZ that indicates an ongoing aggregation process (Westbrook et al., 2007). The rearranged Z DR profile in Fig. 10d shows also signatures that are in agreement with an aggregation process. Below N, the Z DR -values decrease due to the aggregation of the new generated particles into more spherical shaped particles. It is worth noting that such signatures cannot be identified in the vertical, not re-arranged, Z and Z DR profiles (light blue profiles in the Figs. 8b and 10d). This analysis therefore demonstrates the advantage of using along fall streak rearranged radar data, as obtained from Pfitzenmaier et al. (2017). Compared to Case 1, the slope of the Zprofile in Case 2 increases much slower at a rate of 15 dBZ per 1 km (in Case 1 25 dBZ per 1 km). The different slopes of the Z-profiles might indicate different aggregation processes for the two cases.
In Fig. 10b the ρ HV -minimum of 0.956 right above the increase of Z DR at 3.2 km indicates the generation of new particles. Due to a high variation of particle shapes (new generated particles, and seeded particles) within the sampling volume the values of ρ HV lead to that minimum. Below 3.1 km, the ρ HV value increases and reaches values of > 0.99 below 2.8 km. This increase of ρ HV shows that the particle population becomes more and more homogeneous, which indicates an aggregation process.
Next, rearranged fall streak data are analyzed in the spectral domain in order to discriminate between the generated particles at N and the particles seeding from above. The rearranged sZ and sZ DR spectrograms are presented in Fig. 10c and e. At 3.1 km, a broadening is observed in the sZ spectrogram. This broadening of the spectra corresponds to the generation of a second particle mode that is visible at N. The maximum sZ values increase from sZ ≈ −11 dBZ above N to sZ ≈ 7 dBZ right above the melting layer. These observations are indicating aggregation of the newly generated crystals towards the melting layer. In Fig. 10e a bimodality is visible in the sZ DR -spectrogram at N. In comparison to the sZ spectrogram the sZ DR spectrogram contains a bimodal shape from N till about 350 m lower. In the sZ-spectrogram such separated modes cannot be identified. Because of the stable second mode in the sZ DR spectrogram, the assumption of an aggregation process below N is adjusted and a separated diffusional growth of the new generated particles can be assumed, before they aggregate with the ice particles seeded from above.
The analysis of the single sZ and sZ DR spectra in Fig. 11 confirms the hypothesis of diffusional growth. Figure 11 displays four fall streak rearranged Doppler spectra (sZ, left column, black spectra) and differential reflectivity spectra (sZ DR , right column, red line) at four different heights, 3394, 3055, 2864, and 2630 m. At N, a clear broadening in the sZ spectra is visible that is caused by the development of a second particle mode. With decreasing height, the shape of sZ is again monomodal and its values rise, see sZ at 2864, and 2630 m. At 3055 m, it is seen that the newly developed particles at N have high sZ DR -values and the sZ DR spectrum shape turns bimodal. However, in comparison to the sZ that loses its bimodal shape rather quickly the sZ DR maintains it further below. At 2864 m two particle populations can still be identified in the sZ DR spectrum. Between −8.5 and −7.25 m s −1 , high sZ DR -values indicate a large amount of pristine shaped ice crystals, while at Doppler velocities < −8.5 m s −1 low sZ DR -values refer to almost spherical crystals in that part of the spectrum. This clearly shows that from N to 2864 m the generated particles keep their shape dependence and, therefore, can be separated using the sZ DR spectrum. In parallel the sZ values increases for Doppler velocities indicating large particles. Due to this growth sZ loses a clearly visible bimodality that was present at N.
These observed signatures fit to the schematic sketches of the diffusional growth of particles, Fig. 3b. The signatures of the sZ and the sZ DR clearly show that particles generated at N grow separately from the particles seeding from above. Due to the clear size-dependent growth of the smaller particles we assume a diffusional growth of the newly generated particles. Nevertheless, at lower heights the clear separation in the Atmos. Chem. Phys., 18, 7843-7862, 2018 www.atmos-chem-phys.net/18/7843/2018/ Figure 9. Radiosonde launches from the CESAR site at 10:18 UTC (a) and after precipitation system has passed at 12:54 UTC (b). Displayed are profiles of the temperature t (red) and dew-point temperature t d (blue). The light blue bars indicate areas where supercooled liquid water was detected with the radiosondes; their corresponding temperature range on the x-axis is highlighted.
sZ DR spectra disappears. This can be explained by the merging of the two particle populations and, therefore, an aggregation process between the particles, see spectra at 2630 m. At N and at 2864 m, the sZ DR -spectrogram shows negative signatures for large ice particles. The values show that the growth at N leads to prolate particles. Super cooled liquid water droplets can lead to a partial riming of large particles or turbulence could lead to this prolate orientation. However, with decreasing height the shape of these particle populations becomes spherical or even slightly oblate again. Nevertheless, the analysis of the along-fall-streak rearranged spectral data has demonstrated that such processes cannot be identified using integrated volume data. Because using integrated moment data such size dependent spectral signatures are not longer present. The particle type of the new generated particles can be identified by combining the polarimetric measurements from TARA with additional data from the radiosonde measurements and the sL DR -signatures from the Mira radar. Bands of enhanced Z DR are an indicator for the generation and growth of pristine ice particles, mainly dendrites or hexagonal plates (Bader et al., 1987;Kennedy and Rutledge, 2011;Moisseev et al., 2015). However, the visible polarimetric signatures in the radar measurements in Case 2 are caused most likely by newly generated ice needles or columns. The temperature ranges of the radiosonde launches in Fig. 9 show values between −8 and −5 • C at 3 km for both launches. This temperature corresponds to a needle or column generation regime, as in Case 1. In addition, the sL DR -spectrogram observed with Mira shows similar signatures with sL DR -values of −17 dB that agree to the simulation signatures for needles or columns (Aydin and Walsh, 1999;Matrosov et al., 2001;Oue et al., 2015b). Because the cloud conditions are not homogeneous, the hypothesis of the presence of needles is mainly built upon the radiosonde data. The radar Mira sL DR is used to confirm this hypothesis, although both radars measure different sampling volumes. Therefore, it is expected that the enhanced Z DR signature is caused by needles or columns. An indication for needle particles can be the presence of a supercooled liquid water layer at that height. Hogan et al. (2002) showed that small-scale dynamics within frontal system clouds lead to supercooled liquid droplet formation in the cloud system. Such upward motions are visible in the vertical Doppler velocity field during Case 2 in Fig. 5, and the hypothesis of needle generation is considered. Nevertheless, in comparison to Case 1, the aggregation efficiency may be lower. In Hobbs et al. (1974), a strong correlation between the particle number concentration and the aggregation efficiency for needles is mentioned. Lower number concentration would lead to less dense aggregates -this may explain the difference between Case 1 and Case 2.
Finally, it is observed, that the new generated particles at 3.1 km and the related growth of ice crystals lead to an increase of Z in the rain pattern below the melting layer, see Fig. 10a. When the Z DR -band at 3 km is not present anymore, before and just after Case 2, Fig. 4c, the reflectivity values in rain are less, Fig. 4b. The radiosonde temperature profiles and updraft patterns in the vertical Doppler velocity field show the potentiality for supercooled liquid water layer presence around 3 km. The measured sL DR -signatures show signatures of needle or columnar shaped ice particles. Therefore, the existence of a supercooled liquid water layer is possible. Furthermore, the sZ DR -spectrogram clearly exhibits two different particle populations below this height. The decrease of polarimetric signature in parallel to the increase of Z-values suggests the aggregation of the two present particle populations during Case 2. Deeper analysis of the sZ DR spectra indicate that the mode of smaller particles first grows due to the diffusional growth before merging with the other mode. When these two modes merge, the clear separation disappears and the aggregation of the particles becomes dominant. Figure 10. Panel (a) shows the retrieved fall streak at 11:33:31 UTC on top of the Z field. Panels (c, e) are the fall streak corrected spectrograms (sZ and sZ DR , 45 • elevation -note that the Doppler velocity contains the radial wind). Panels (b, d) show the fall streak corrected profiles of ρ HV and Z DR in black, while light blue represents the vertical Z DR profile at 11:33:31 UTC. Panel (f) shows the sL DR spectrogram of the vertical pointing Mira at 11:32:00 UTC (all data displayed in the spectrograms have SNR > 10 dB). ML points out the signature of the melting layer and N the signature related to the new particle population.
6.3 Case 3: particle growth of hexagonal particles, 11:35-11:47 UTC Figure 12 shows the analysis of the retrieved fall streak from 11:47:30 UTC. The focus of Case 3 is related to the increased Z DR signature at 5 km between 11:35:00 and 11:47:00 UTC. At that height the Z-values increase and indicate a particle growth. The fall streak rearranged data are retrieved using a cloud base height of 2.25 km and a 90 s averaged wind profile. The vertical profiles and fall streak corrected data are only analyzed from cloud top till 3.0 km. This is done because at 3.0 km a strong shear in wind direction is visible in Fig. 5d). Therefore, a homogeneous cloud cannot longer be assumed and the along the fall streak rearranged data would describe a different particle population below 3.0 km. So the analysis focuses at regions above 3.0 km and no link to the increased precipitation pattern is done for that case. The case is anyhow challenging. Because the data above 3.0 km show signatures of aggregation and some indications for riming are found. Therefore, the focus of Case 3 is the identification of the most probable particle growth process. Figure 4b depicts a horizontal homogeneous Z-field during Case 3, where the main increase of reflectivity is visible between 5 and 4.5 km. In the fall streak corrected sZ spectrogram, Figs. 12c and 13, the increase of sZ can be localized at the same height range. There the maximum values of sZ increase from −5 dBZ above 5 km to sZ ≈ 6 dBZ at around 4.5 km (sZ ≈ 7.5 dBZ at 4 km). The Z-profile in Fig. 12d (black line) shows constant values above 5 km (Z ≈ 10 dBZ) and below 4 km (Z ≈ 20 dBZ). The main growth process of the tracked particle population is visible between 5 km and 4 km, where the reflectivity increases linearly of 10 dBZ, which indicates an ongoing aggregation (Westbrook et al., 2007). This signature is consistent with other observations and, therefore, the studies of the microphysics are based on fall streak rearranged TARA data.
The sZ DR spectrogram exhibits a clear shape-size dependency even below the particle growth. The sZ DR signature Atmos. Chem. Phys., 18, 7843-7862, 2018 www.atmos-chem-phys.net/18/7843/2018/ Figure 11. Along the fall streak at 11:33:31 UTC rearranged sZ and sZ DR at four different altitudes. The left panel shows Doppler spectra at 3394, 3055, 2864, and 2630 m. The right column displays the corresponding spectral differential reflectivity at the same heights. All spectra are at 45 • elevation and averaged over 3 consecutive time bins (∼ 8 s).
consists of high values at the right edge where the small particles are present. This signature is constantly present between 4.5 and 3 km. The maximum values of sZ DR are just below 5 km, sZ DR ≈ 2.25 dB, sZ DR -spectra at P in Fig. 13, and decrease with the height. Nevertheless, the shape-size dependency stays. The fall streak corrected Z DR profile in Fig. 12f exhibits a maximum of 0.75 dB between 5.0 and 4.5 km. Below Z DR values are rather constant around 0.4 dB till 3.0 km. Such a bulk volume is challenging to interpret without the spectral polarimetric measurement. The profile of ρ HV in Fig. 12b increases with the observed particle growth below 5 km (Z profile) and stays constant from 4.5 to 3 km. It also exhibits a minimum before the increase around 5 km (slightly above the maximum of Z DR ). This is caused by a large variety of particle shapes before the growth, like in Case 2. Below the growth process ρ HV shows a constant value that is 0.05 lower than in Case 2. This is caused by the shape-size relation that is present in the sZ DR spectrogram. Because this signature stays constant till 3 km further growth processes of the tracked particle population like aggregation do not occur.
The constant shape-size dependency within the sZ DR spectrogram between 4.5 and 3.0 km is the most dominant feature. Comparing these signatures to the schematic sketches in Fig. 3 we can opt for riming or aggregation of particles. In Fig. 13 the Doppler spectra do not show a separated mode. Simulations of sZ DR signature indicate that the observed signature at P in Fig. 13, 4.794 km, can be produced by a mixture of pristine hexagonal shaped ice particles and aggregated or rimed particles (Spek et al., 2008). Like in Case 1 and 2, the trend of the sZ DR versus Doppler velocity can be examined, when it is stable for several heights and times. However, caution is required for the values of sZ DR because of their large variance. The absolute values of the Doppler velocities of the Mira spectra (not shown) are below 2 m s −1 . They correspond to pristine and slightly rimed particles (Mitchell, 1996). Due to the homogeneity of the particle population below 4.5 km in terms of sZ and sZ DR we assume that the particles mainly grow at P and then seed further through the cloud system. Below 4 km the sZ DR spectrogram and spectra in Figs. 12e and 13 (right column) show that the larger particles could be prolate. Negative differential reflectivity signatures can be expected for conical graupel presence (Oue et al., 2015a). However, the observations of Oue et al. (2015a) indicated that the reflectivity of graupel particles is much higher than in our observations. Riming of particles that lead to prolate particles is also possible but the presence of conical graupel cannot be confirmed. Nevertheless, the observed vertical Doppler velocities of the Mira radar are within the expected range for slightly rimed particles. In comparison to the schematic sketches in Fig. 3a no separated particle mode is visible. However, the Doppler spectra at 4031, and 3309 m in Fig. 13 exhibit a small bump where sZ DR signatures are negative. Therefore, a small amount of rimed prolate oriented particles could be present.
The identification of the dominant particle shape at 5 km can indicate the dominant growth process. As dendrites and plates have different efficiencies of aggregation and riming (Pruppacher and Klett, 1996;Lamb and Verlinde, 2011), the temperature of radiosonde launches is used to identify the possible pristine particle growth regime. In Fig. 9 both profiles indicate a temperature range between −15 and −13 • C, but a direct identification of a supercooled liquid water layer around 5.0 km can not be done (dew-point-temperature and temperature profiles do not match at that height). Neverthe- less, the temperature range at 5.0 km agrees with the growth regime of hexagonal plates or dendrites (Pruppacher and Klett, 1996;Fukuta and Takahashi, 1999). As in Case 2, the vertical Doppler velocity measurements indicate updraft patterns in that region that increases the probability of the presence of supercooled liquid water droplets (Hogan et al., 2002). However, from the observations it is not clear which pristine particle shape is dominating the growth process at that height: hexagonal plates or dendrites. Summing up all the facts, the polarimetric measurements, the temperature range, and the updraft patterns we conclude that the presence of supercooled liquid water around P is possible. Because no separated riming mode could be identified in the spectra we do not expect riming to be the dominant particle growth process. It is more likely that the constant reflectivity values and the constant shape-size relationship along the fall streak re-arranged data below P are caused by an aggregation process. However, more measurements of the particle microphysics would be necessary to get a more detailed picture on which growth process might dominate in that case.

Discussion and conclusion
The high temporal and spatial resolution of profiling S-band radar TARA made it possible to detect and identify structures of particle generation and growth within complicated and dynamically complex mixed-phase cloud systems. Using a S-band radar system has the advantage that the attenuation of the radar signal during precipitation can be neglected compared to higher frequencies operating systems (Kollias et al., 2007). The unique three beam configuration of TARA makes it possible to retrieve the full 3-D wind field for each measurement (Unal et al., 2012). Furthermore, this dynamical information is used to retrieve particle fall-streaks within the bulk parameter fields, Pfitzenmaier et al. (2017). Therefore, the bulk parameter profiles, Z, Z DR , and ρ HV , as well as spectrograms, sZ and sZ DR , are rearranged along such retrieved fall streak before they are analyzed. Consequently, these corrected spectrograms and profiles show the consistency of the particle growth within the cloud part and relate better to the increase of the rain pattern below the melting layer. Such increased rain patterns are identified during Case 1 and 2. The analysis of the along fall streak rearranged measurements allows to localize the particle growth process. This can be done by examination of the increase and the slope of Z along the fall streak. Westbrook et al. (2007) showed that the aggregation of ice crystals leads to a constant increase of Z with decreasing height. The observations were based on averaged data of homogeneous ice clouds over a long time period. Using rearranged data along fall streak such studies can be done in less homogeneous cloud conditions. The presence of supercooled liquid water leads to a much faster growth of ice particles within a shorter amount of time (fast increase of Z within short height ranges) compared to the growth of particles by aggregation in the ice phase only (Pruppacher and Klett, 1996). That may explain the increase in the Z profile within 1 km of 25 dBZ in Case 1, 15 dBZ in Case 2, and 10 dBZ in Case 3.
The fall streak corrected sZ DR measurements of TARA are able to deliver shape information and their distribution within the tracked particle population. The changes of the observed signature is in good agreement with the discussed particle growth (increase of Z). However, stable trends of sZ DR should be examined, rather than absolute values because of the large variance of this spectral polarimetric parameter.
In Case 3 the observed polarimetric signature in the sZ DR spectrogram agrees with the simulated sZ DR signature for a mixture of hexagonal plates and rimed particles or aggregates in Spek et al. (2008). Therefore, we can opt for riming and/or aggregation of particles in that case. Aggregation is assumed to be the dominant particle growth process. Nevertheless, negative sZ DR -values are observed in the spectrogram, which may indicate rimed prolate oriented particles. However, a typical separated riming mode within the Doppler spectra sZ is not visible. The Mira vertical Doppler velocities fit to the velocity range of slightly rimed particles. Observed updraft patterns indicate the possible presence of supercooled liquid water at the region of the main particle growth, like observed by Hogan et al. (2002). Therefore, riming cannot be neglected in Case 3, however, explicit riming signatures cannot be observed in the measured spectra.
The observed signatures of sL DR and sZ DR as well as the measured temperature in Case 1 are in agreement with modeled signatures of needles and their strong tendency to clump and aggregate as reported by Rangno and Hobbs (2001). In Case 2 the observed signatures are more complex. While the temperature range can be considered the same as Case 1 the polarimetric signature differs. One reason could be the difference in aggregation efficiency of the new produced particles in these two cases. In Case 2 the needles show a lower aggregation efficiency with the ice particles seeded from above, than in Case 1 where only needles were present. The decrease of the polarimetric signature is less strong than in Case 1. Another reason might be a difference in the particle concentration of the generated particles that is less in Case 2 (sZ of the increased sZ DR mode is less that the sZ values in Case 1). Then the aggregation efficiency is less and the decrease of the sZ DR signature takes longer. Nevertheless, from the observations it is not visible which mechanism is dominant or if they depend on each other. It is only seen that a second population of particles with a high shape size dependency is present in the data. These polarimetric signatures as well as the corresponding temperature profiles indicate the presence of supercooled liquid water.
Parallel radiosonde launches and the sL DR of the vertical pointing K a radar Mira are used to compare and verify the results based on the TARA measurements. The cases presented in Sect. 6 are examples of the spectral and bulk parameter characteristics observed during the whole event. These analyzed ice particle growth processes above the melting layer show a strong correlation with increased pattern of precipitation in the reflectivity fields. Therefore, the presented data with applied fall streak correction demonstrate the advantage of using also high resolution S-band radar data to increase the understanding of ice particle growth within precipitating mixed-phase cloud systems.
Author contributions. CU designed the data processing of TARA, contributed to the quality control of the data and the result discussion. YD helped organizing and performing the ACCEPT campaign and discussing the results. HR designed the scientific structure of the project, helped in the realization of ACCEPT, and contributed to the discussion of the results. LP performed the measurements taken during the ACCEPT campaign in autumn 2014, Cabauw, the Netherlands, supported by the colleagues from TROPOS. LP performed the fall streak analysis, contributed to the result discussion, and prepared the manuscript with contributions of all coauthors.
Competing interests. The authors declare that they have no conflict of interest.