Simulated rimed aggregates generated with https://github.com/jleinonen/aggregation in setting "aggregation followed by riming". Aggregates were built from between 10 to 700 monomer crystals of columns, dendrites, needles, plates or rosettes with mean sizes of 100 or 200 micrometer. Then they were exposed to ELWP = 2.0 kg m⁻². Monomer crystals are composed of cubical elements with resolution 20 micrometer. Frozen rime droplets are also represented by 20 micrometer cubes. The data set contains folders with evolution (evol) and shape files for each monomer crystal type. For each particle one evolution and one corresponding shape file exists. The evolution (evol) file contains particle mass, rime mass, area, size, fall speed (Heymsfield&Westbrook, 2010), fall speed (Khvorostyanov&Curry, 2005) for each step during the aggregation and riming process. The corresponding shape file contains the x,y,z positions of the cubical elements that compose the particle for each step. For further documentation see readme.

Riming is a key process of precipitation formation in ice-containing clouds, but quantifying riming from observations is challenging, limiting our ability to evaluate the riming process in numerical weather models. One challenge for radar observations is that riming changes both the physical properties (mass, area cross-section) and scattering properties of ice particles. These changes need to be implemented consistently as a function of riming in radar forward operators, which are required for retrievals and model evaluation in observation space. In this study, mass–size, cross-section area–size, and backscattering cross-section relations are developed as a function of the normalized rime mass for aggregates composed of various monomer types (columns, dendrites, needles, plates, and rosettes). The proposed framework allows us to simulate scattering properties of aggregated ice particles consistently as a function of riming in retrievals and radar forward operators. The parameterizations are developed from a large data set of simulated rimed aggregates of different sizes and monomer crystal types. The backscattering cross-section parameterization (the “riming-dependent parameterization”) is evaluated for radar frequencies of 35.6 and 94.0 GHz and is based on the Self-Similar Rayleigh–Gans approximation (SSRGA), which is increasingly used to calculate microwave scattering of ice crystals and snowflakes. Compared with parameterizations from the literature that do not consider riming, the riming-dependent parameterization leads to significantly smaller biases in terms of backscattering cross-section. When using the particle masses and scattering properties of the individual particles simulated by the aggregation and riming model as a reference, the bias of our parameterization is below 1 dB when integrating over an exponential particle size distribution with sizes from 0.1–10 mm.

The dendritic growth layer (DGL), defined as the temperature region between -20 and -10 °C, plays an important role for ice depositional growth, aggregation and potentially secondary ice processes. The DGL has been found in the past to exhibit specific observational signatures in polarimetric and vertically pointing radar observations. However, consistent conclusions about their physical interpretation have often not been reached. In this study, we exploit a unique 3-months dataset of mid-latitude winter clouds observed with vertically pointing triple-frequency (X-, Ka-, W-band) and polarimetric W-band Doppler radars. In addition to standard radar moments, we also analyse the multi-wavelength and polarimetric Doppler spectra. New variables, such as the maximum of the spectral differential reflectivity (ZDR) (sZDRmax), allows us to analyse the ZDR signal of asymmetric ice particles independent of the presence of low ZDR producing aggregates. This unique dataset enables us to investigate correlations between enhanced aggregation and evolution of small ice particles in the DGL. For this, the multi-frequency observations are used to classify all profiles according to their maximum average aggregate size within the DGL. The strong correlation between aggregate class and specific differential phase shift (KDP) confirms the expected link between ice particle concentration and aggregation. Interestingly, no correlation between aggregation class and sZDRmax is visible. This indicates that aggregation is rather independent of the aspect ratio and density of ice crystals. A distinct reduction of mean Doppler velocity in the DGL is found to be strongest for cases with largest aggregate sizes. Analyses of spectral edge velocities suggest that the reduction is the combined result of the formation of new ice particles with low fall velocity and a weak updraft. It appears most likely that this updraft is the result of latent heat released by enhanced depositional growth. Clearly, the strongest correlations of aggregate class with other variables are found inside the DGL. Surprisingly, no correlation between aggregate class and concentration or aspect ratio of particles falling from above into the DGL could be found. Only a weak correlation between the mean particle size falling into the DGL and maximum aggregate size within the DGL is apparent. In addition to the correlation analysis, the dataset also allows study of the evolution of radar variables as a function of temperature. We find the ice particle concentration continuously increasing from -18 °C towards the bottom of the DGL. Aggregation increases more rapidly from -15 °C towards warmer temperatures. Surprisingly, KDP and sZDRmax are not reduced by the intensifying aggregation below -15 °C but rather reach their maximum values in the lower half of the DGL. Also below the DGL, KDP and sZDRmax remain enhanced until -4 °C. Only there, additional aggregation appears to deplete ice crystals and therefore reduce KDP and sZDRmax. The simultaneous increase of aggregation and particle concentration inside the DGL necessitates a source mechanism for new ice crystals. As primary ice nucleation is expected to decrease towards warmer temperatures, secondary ice processes are a likely explanation for the increase in ice particle concentration. Previous laboratory experiments strongly point towards ice collisional fragmentation as a possible mechanism for new particle generation. The presence of an updraft in the temperature region of maximum depositional growth might also suggest an important positive feedback mechanism between ice microphysics and dynamics which might further enhance ice particle growth in the DGL.