Snow formation is a complex interplay of multiple microphysical growth processes, and the prevailing snow characteristics are inherently linked to local climate. However, the persistent shortage of observations for characterizing snow microphysics at a global scale continues to constrain our understanding of snow growth processes. Here, we investigate snow riming and aggregation signatures in stratiform precipitation through triple-frequency radar observations collected during coordinated field campaigns across Southern China, the Eastern United States, Western Europe, Northern Europe and Antarctica. The results suggest that the velocity-based riming estimates are generally consistent with triple-frequency observations, and the riming frequency increases with temperature. Our analysis of dual-frequency observations in these field campaigns qualitatively indicate the dendritic growth zone around -15 °C playing a key role in initiating enhanced snow size growth, and reveals a generally temperature-dependent snowflake growth characteristics. The snow over Eastern US is characterized by the most prominent riming growth, corresponding to moderate to heavy riming. Triple-frequency signatures of snowflakes over west Europe are consistent with Southern China, while the latter shows a higher degree of riming. The weakest snow growth signatures were found over west Antarctica, potentially owing to the scarcity of ice nucleating particles and available water vapor for deposition. In addition, our statistics reveal a latitudinal dependence for snowfall detection limitations with current spaceborne Ku- and Ka-band radars, and shed novel insights into future triple-frequency satellite missions as well as joint application of weather and spaceborne radars.
National Natural Science Foundation of China42475095China Meteorological Administration2024-K-01KDW2412State Key Laboratory of Severe Weather2025QZA04Natural Science Foundation of Anhui Province2408055UQ007Natural Science Foundation of Shandong ProvinceZR2025LQX001Introduction
Half of the Earth's surface precipitation events originate from snow . Formation and growth of snow are governed by multiple interacting microphysical processes, such as, nucleation, secondary ice production, deposition, riming, aggregation, which take effect on various pathways depending on atmospheric dynamics and thermodynamics . Hence, the prevailing snow characteristics (size, shape, density, falling velocity, etc) are dependent on regional atmospheric conditions, such as temperature, humidity, air motions, etc. For instance, aggregation is significantly enhanced around the melting layer and is affected by vertical air motions ; riming is favorable in shallow Arctic mixed-phase clouds at temperatures greater than -12 °C and is necessary for graupel formation . The concentration of ice nucleating particles (INPs) largely depends on temperature, in addition to other factors such as aerosol compositions and concentrations, and presents significant site-to-site variations .
In spite of the importance of regional climate to snow microphysics, explicitly representing microphysical processes in models is challenging owing to our knowledge gaps in cloud physics as well as the simplified microphysics schemes that are inherently uncertain and lack observational constraints . For example, the parameterized aggregation efficiency in many microphysics schemes monotonically increase with temperature , while an increased sticking efficiency at around -15 °C is more consistent with actual observations . Although riming has been implemented in some microphysics schemes , the interactions between riming and aggregation as evidenced in recent observations suggest the need of more delicate tunning of snow microphysics in models. Therefore, adequate observations are required for disentangling the process-level understanding of snow microphysics and improving the representation of specific physical processes in numerical models.
Major responsible microphysical processes for snow characteristics can be inferred from in-situ observations. In well-designed field campaigns, ice particles in natural clouds have been extensively recorded with in-situ probes. Since these instruments are mounted on specific platforms, such as balloons, aircrafts, or cable cars, profiling the microphysical processes taking place in clouds is complemented by the means of remote sensing. Thanks to the development of active remote sensing techniques, meteorological radars have shown promise in characterizing snow microphysics. Since the dawn of meteorological radars, snow has been identified from radar echoes, and then radars have proven to be a unique tool for observing snow . In recent decades, significant improvement has been made in quantitative characterization of snow microphysics thanks to the advances in understanding the scattering characteristics of snow as well as the implementation of Doppler, dual-polarization and multi-frequency radar techniques among others.
The basis of radar remote sensing of snow microphysics lies in the interactions between electromagnetic waves and ice particles, known as scattering. If the maximum dimension (D) of a particle is much smaller than the radar wavelength (D≪λ/10), the Rayleigh scattering is satisfied and the observed radar reflectivity is proportional to D6. At commonly-used weather radar wavelengths, e.g., S- (∼ 10 cm), C- (∼ 5 cm), and X-band (∼ 3 cm), the difference between radar reflectivities (in logarithmic units; dBZ) at two frequencies called dual-wavelength ratio (DWR) is usually 0 dB, since the radar wavelength is much larger than the hydrometeor dimensions. As the radar wavelength decreases to the magnitude of millimeters, e.g., K- (∼ 8 mm) and W- (∼ 3 mm) band, which are comparable to snow dimensions, the non-Rayleigh scattering appears, and DWR >0 dB. Therefore, the observed non-zero signatures of DWR are linked to snow sizes. For vertically-pointing radars, DWRKa,W and DWRX,Ka start exceeding 1 dB at particle maximum dimensions of about 0.75 and 2 mm, respectively . Neglecting the effect of radar signal attenuation, DWR values are dependent on scattering of targets, which is the basis of dual- and triple-frequency radar retrievals of snow microphysics.
and have demonstrated that the DWR between non-Rayleigh and Rayleigh scattering frequencies can be used to estimate particle median sizes. Scattering calculations of physically realistic snowflake shapes using discrete dipole approximation suggested that aggregation of unrimed snow leads to a “hook” feature (thick blue curve in Fig. 1c) which is different from the growth of rimed snow (thick yellow curve in Fig. 1c) . Such signatures have been evidenced in airborne and ground-based radar observations . Furthermore, the transition from “hook” to flat signatures in triple-frequency map is indicative of the prevailing snow growth shifting from aggregation to riming , which facilitates a process-level assessment of snow growth characteristics. identified unusual triple-frequency signatures of riming over Antarctica, and discussed responsible factors using model simulations. using triple-frequency observations validated model parameterizations of snow aggregation, which confirmed the need to represent an additional peak in snow sticking efficiency at -15 °C. In addition, forward simulations of spaceborne triple-frequency observations indicate that at least one frequency can be reliably employed to monitor the entire cloud-precipitation process .
(a) Overview of the triple-frequency radar field campaigns. (b) C/Ku/Ka/W-band radars in METRICs. (c) Triple-frequency signatures of physical parameters of snowflakes. D0, ρ, and μ denote median volume diameter, snow density, and the shape parameter of the Gamma size distribution, respectively. Since the average value of μ is around 0 , our interpretation of triple-frequency statistics bears the assumption of μ=0, namely the exponential distribution. The BAECC, AWARE and IMPACTS
Understanding of snow growth obtained from ground-based triple-frequency campaigns may also bring seminal insights into next-generation satellite missions . Spaceborne dual- or triple-frequency radar retrieval algorithms have been developed based on the non-Rayleigh scattering signals among different frequency bands and validated against in-situ observations . On the other hand, the small ice particles at cloud tops are good Rayleigh-scattering targets for cross-calibration among radars with different frequencies deployed on various platforms. For example, spaceborne radars, such as Tropical Rainfall Measuring Mission TRMM,, CloudSat , Global Precipitation Measurement GPM,, Fengyun-3G , and Earth Cloud Aerosol and Radiation Explorer EarthCARE,, have been proposed for cross-calibration of ground-based radars . However, there is no consensus on the ZRayleigh value that can be used to identify the cloud top Rayleigh scattering regions. have shown significant latitude-dependence of characteristic sizes of snow at cloud tops, suggesting the potential geographical dependence of snow sizes. In addition, cloud formation temperature and cloud type which changes with latitudes are relevant to snow size distributions and snow characteristic sizes . This impact seems to be detectable in radar observations. showed that W band reflectivity starts deviating from ZRayleigh (ZRayleigh,W) at about -20 dBZ using aircraft observations during European Cloud Radiation Experiment, while found it ranges from -20 to -10 dBZ in different field campaigns. Therefore, examining the geographical impact on snow microphysics also benefits the cross-calibration of multi-frequency radars.
Since 2014, a number of field campaigns hosting high-sensitivity triple-frequency radars have been carried out over regions spanning from Antarctica to tropics (Fig. 1a). The well-calibrated and -aligned long-term radar observations open a new opportunity to assess snow microphysical processes over various climatologies. In this study, we attempt to compare the triple-frequency radar observations obtained in these campaigns and assess the geographical fingerprints in snow microphysical processes.
The remainder of this paper is organized as follows. The second section overviews five triple-frequency field campaigns, followed by the methods of snow identification, radar reflectivity calibration and riming estimation. The results are presented in the fourth section. Conclusions are given in the final section.
Data
We examine stratiform precipitation (rainfall with bright band signatures and snowfall, cloud cases were removed for minimizing the impact of sublimation) with ground-based/airborne multi-frequency radars in five field campaigns. Namely, the Biogenic Aerosols-Effects on Clouds and Climate campaign BAECC,, TRIple frequency and Polarimetric radar Experiment for improving process observation of winter precipitation TRIPEx-pol,, the Atmospheric Radiation Measurements West Antarctic Radiation Experiment AWARE,, Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms IMPACTS,, and the Multi-frequency radar Experiment for TRopical Ice Clouds (METRICs). As shown in Fig. 1a, these campaigns were conducted over various latitudes, ranging from tropical (METRICs), mid-latitude (IMPACTS, TRIPEx-pol), high-latitude (BAECC) to polar (AWARE) regions. Note that the X/Ka/W-band setups were employed in IMPACTS, TRIPEx-pol, BAECC and AWARE, and C/Ka/W-band radars were used in METRICs. For simplicity, we use ZX/C to refer X- (IMPACTS, TRIPEx-pol, BAECC, and AWARE) or C-band (METRICs) reflectivity. Similarly, DWRX∕C,Ka denotes the dual-wavelength ratio between X- or C-band reflectivity and Ka-band reflectivity.
METRICs
From July to September 2022, METRICs took place at the Shenzhen Shiyan Observatory (22.65° N, 113.89° E), China. During the three-month campaign, five radars were deployed at the Shiyan observatory for synergetic observations, including vertically-pointing W-, Ka-, Ku-, and C-band profiling radars, and a L-band wind profiler . The cloud radars were operating with a distance to each other less than 5 m. About 5 km to the Observatory, a X-band phased-array radar and an operational dual-polarization S-band weather radar were in operation, providing reflectivity calibration basis for the C-band radar. The C-band reflectivity was calibrated by matching with the S-band radar reflectivity observations at the height of 0.5–1 km during rainfall events. Similar with previous triple-frequency radar setups, W-, Ka-, and C-band radars were used in this study. Their range resolutions are 30 m, and the time resolutions are 3, 25, 0.93 s, respectively. The sonding data including temperature, pressure and humidity from Hongkong Observatory which was launched two times per day and is about 40 km from Shiyan Observatory was used for gaseous attenuation and air density correction. The stratiform rainfall observations of 20 h were used in this study.
TRIPEx-pol
TRIPEx-pol was carried out at the Jülich Observatory (50.9° N, 6.4° E) from November 2018 to January 2019. During TRIPEx-pol, vertically pointing X, Ka and W band radars were installed at the same roof platform with the horizontal distances less than 20 m. The X- and Ka-band systems are pulsed radar systems manufactured by Metek GmbH, while the W-band radar is a frequency modulated continuous wave (FMCW) system manufactured by Radiometer Physics GmbH. The range resolution of these radars is 30 m and the time resolutions are 2, 2 and 3 s, respectively. The radar reflectivity was calibrated with the simulated reflectivity from drop size distributions (DSDs) measured by the PARSIVEL optical disdrometer during 21 rainfall periods . Vertical profiles of atmospheric temperature, pressure, and humidity were from the European Centre for Medium-Range Weather Forecasts-Integrated Forecasting System (ECMWF-IFS) forecasts, and we used the interpolated model products over the Jülich Observatory. In this study, the level 2 data products which have been well calibrated in total of 60 h stratiform rainfall and 18 h snowfall were used.
BAECC
The BAECC field campaign was conducted at the University of Helsinki's Hyytiälä Station (61.8° N, 24.3° E) from February to September 2014 . This experiment hosted comprehensive vertically pointing multi-frequency radars, including X/Ka-band scanning ARM cloud radar (X/Ka-SACR), Ka ARM zenith radar (KAZR) and W-band ARM cloud radar (MWACR). In this study, X-SACR, KAZR, MWAR observations were used. In addition, radiosondes were launched every six hours to obtain temperature, pressure, and humidity. The X-band reflectivity was corrected with a collocated operational C-band radar during snowfall , and surface DSDs observations were used to calibrate X-band radar reflectivity at 500 m during rainfall . In BAECC, triple-frequency radar observations of stratiform rainfall (16 h) and snowfall (12 h) were used.
AWARE
Within the framework of the Western Antarctic Radiation Experiment of Atmospheric Radiation Measurement (ARM), the second ARM mobile equipment (AMF2) was deployed at McMurdo Station (77.83° S, 166.67° E). The multi-frequency radars include KAZR, MWACR and X/Ka-SACR , which have been used in BAECC as well. During AWARE triple-frequency radars were in operation from December 2015 to January 2016 in which the surface precipitation was in the form of snowfall. A radiosonde was launched four times a day to acquire temperature, pressure, and humidity profiles. Absolute calibrations of the scanning radar systems were performed on site with a corner reflector, and a systematic comparison was conducted with nearby CloudSat measurements . The recorded snowfall was in total of 23 h.
IMPACTS
IMPACTS is a multi-year field campaign with an airborne multi-frequency radar setup. Sponsored by NASA, IMPACTS was conducted to study wintertime snowstorms focusing on East Coast cyclones during the winters of 2020–2023 . The Earth Resources 2 (ER-2) aircraft flew above snowfall systems, and carried nadir-pointing radars at frequencies of X-, Ku-, Ka-, and W-bands. The IMPACTS data has already been quality controlled including aircraft motions and attitudes (https://www.earthdata.nasa.gov/data/projects/impacts/collection, last access: December 2025), and has been employed in previous studies . To assure that the radar ray path is not excessively prolonged during aircraft turns, we have removed periods with rolling angles exceeding 0.2°. We also applied the X-band reflectivity threshold of 50 dBZ to identify surface echoes. To minimize the surface clutter contamination, radar observations within 0.5 km to surface were removed. For facilitating the triple-frequency analysis, radar observations were interpolated into temporal and range resolutions of 0.5 s and 26.25 m, respectively. The hourly High-Resolution Rapid Refresh HRRR, analysis data was used to provide temperature, pressure and humidity for gaseous attenuation and air density correction. The quality-controlled IMPACTS dataset includes 16.5 h snowfall and 3.5 h stratiform rainfall.
Methods
In our datasets, AWARE was focused on Antarctic snowfall, while METRICs recorded tropical rainfall. TRIPEx-pol, IMPACTS, and BAECC datasets include both rainfall and snowfall events. To generate the triple-frequency map, we identified snow in rainfall events and carefully considered different attenuation sources. To independently evaluate snow riming signatures from triple frequency observations, we quantified riming with a velocity-based approach and constructed characteristic DWRX∕C,Ka-DWRKa,W curves with different riming degrees.
Identification of Snow
Snowfall events can be identified with surface temperature provided the absence of temperature inversion over 0°, while it is essential to identify snow above the melting layer in stratiform rainfall events. In radar observations, the melting layer is characterized by distinct changes of polarimetric variables. The significant enhancement of linear depolarization ratio (LDR) as observed by vertically pointing radars has been widely used for melting layer detection, despite that there is slight frequency dependence (tens of meters) on the melting layer top height detection . Here, Ka-band LDR is employed for melting layer identification. For a robust melting layer detection, we follow the principle of identifying the significant changes of LDR gradients (dLDR) . The search for melting layer top starts from snow region, and stops at a reference level (Href) where the LDR is 3 dB above that in snow and dLDR values are consecutively positive from this level to 90 m below. To avoid potential impact of the partial melting, the melting layer top is determined 100 m above Href.
Attenuation Correction to Radar Reflectivity
Although radars in field campaigns had been well calibrated, attenuation correction needs to be considered for DWR analysis. For up-looking ground-based radars, we corrected wet radome attenuation, applied the relative calibration at the cloud top, and corrected the overestimated reflectivity due to gaseous and hydrometeor attenuation. For down-looking airborne radars in IMPACTS, we only need to account for gaseous and hydrometeor attenuation which leads to reflectivity underestimation.
Firstly, the wet radome attenuation in rainfall events was corrected by matching the observed reflectivity at 500 m and the simulated reflectivity from surface DSDs observations (BAECC, TRIPEx-pol). In METRICs, the collocated S-band weather radar (5 km away) data was used to calibrate the C-band radar reflectivity. Then, the X/C-band reflectivity profile is assumed to be well calibrated, since the attenuation from snow, rain and melting layer in stratiform rainfall is negligible .
Secondly, the relative calibration at Ka- and W-bands was made by identifying small ice particles at the cloud top. Following the algorithm given by , we identified the Rayleigh-scattering regions at Ka- and W-bands using C(X)-/Ka-band and Ka-/W-band pairs, respectively. Moving downwards from the cloud top, DWR remains constant until the non-Rayleigh scattering at the higher frequency starts appearing. The Ka-band reflectivity was adjusted by matching the C(X)-band reflectivity at cloud tops. Then, this approach was applied to the Ka/W-band radar pairs. The presence of strong radar signal attenuation, e.g., in hailstorms or rainstorms, the Ka- and W-band radars suffer from sensitivity losses and their signals may not reach the level where Rayleigh-scattering ice particles exist. Hence, the algorithm cannot detect a region of flat DWR at cloud tops, and such profiles will be discarded. Lastly, we need to consider the gaseous and hydrometeor attenuation. The water vapor and oxygen gaseous attenuation was corrected using temporally interpolated sounding observations as input to millimeter wave propagation model . The attenuation from snow and supercooled liquid water is negligible at Ka band, but can be significant at W-band. We corrected the snow attenuation using ASnow,W=0.0325ZW,lin, where ZW,lin is ZW in linear scale . This approach yielded a median W-band path-integrated snow attenuation of 0.3 dB in all field campaigns, comparing to the value of 0.1 dB from another parameterization approach . Therefore, we can reasonably deduce that the uncertainty of snow attenuation correction is less than 1 dB.
The attenuation from supercooled liquid water is not accounted in this study. This impact on Ka-band radar is relatively small, but can be significant at W-band. have shown that the W-band attenuation from a supercooled liquid water path (SLWP) of 300 g m−2 is on the magnitude of 2 dB, and found that the occurrence of SLWP exceeding 300 g m−2 in snowfall is below 15 % over central Finland. From the perspective of triple-frequency signatures, the supercooled liquid attenuation at W-band can lead to overestimated DWRKa,W in IMPACTS data (top-down view), which is similar to the effect of riming as inferred from radar Doppler velocity observations. For ground-based radars (down-top view), the uncorrected liquid water attenuation at W-band leads to underestimation of DWRKa,W. Nonetheless, riming as inferred from triple-frequency observations can be compared to those estimated from radar Doppler velocity which is immune to the attenuation effect as will be discussed later.
Doppler-velocity-based estimation of rime mass fraction
In correspondence with the triple-frequency signatures, process of riming leads to increased radar Doppler velocity . To independently characterize the riming signatures, we use radar Doppler velocity to quantify the rime mass fraction ,
FR=1-∫DminDmaxN(D)mur(D)dD∫DminDmaxN(D)mob(D)dD
where Dmax and Dmin are maximum and minimum particle sizes, respectively; mob(D) and mur(D) are masses of observed and unrimed snowflakes as a function of snow diameter D, respectively; N(D) is the particle size distribution. If mob(D) is smaller than mur(D), FR is assigned to 0. Because riming leads to increased ice density and fall velocity, FR can also be estimated from vertically-pointing radar Doppler velocity observations . Following the approach used by , updraft regions were removed, and a time-average window of 20 min was imposed to cancel out up-and-downward motions. FR was estimated using the fitting between FR and the mean Doppler velocity (MDV) observed by Ka-band radars as below,
FR=0.0791MDV4-0.5965MDV3+1.362MDV2-0.5525MDV-0.0514
To minimize the impact of varying air density (ρair), MDV was adjusted to the air condition of 1000 hPa and 0 °C (air density air, ρair,0) with a factor of (ρair,0ρair)0.54. ρair was derived from the temperature and relative humidity obtained from the interpolated sonde observations.
DWR<sub>X∕C,Ka</sub>-DWR<sub>Ka,W</sub> signatures of riming
In triple-frequency space, riming and aggregation can lead to diverse signatures . To compare our observations with previous scattering simulations, we followed the approach by . Firstly, we employed the backscatter cross sections of individual “realistic” dendritic snowflakes with various riming and aggregation degrees at different frequencies as simulated by , and generated DWR simulations using the exponential distribution. Although the Gamma function provides a better fit for snow size distributions, the shape parameter μ is around 0 in statistics , supporting the adequacy of using the exponential distribution in this study. In the exponential distribution N(D)=N0exp(-ΛD), where D, N0 and Λ are snow diameter,the intercept parameter and the inverse scale parameter, respectively, the sizes of snowflakes are controlled by Λ. Changing Λ from 2 to 367 mm−1, we can get different characteristic sizes of snow populations and map snow growth signatures into the triple-frequency space . Then, we quantified FR using the mass-size relation from dataset as input to Eq. (1). As shown in Fig. , the characteristic DWRX∕C,Ka-DWRKa,W signatures with FR =0, 0.53, and 0.72 are denoted by circled, squared, and star-shaped curves, respectively.
ResultsCloud top temperature and riming occurrence
Ice formation at the cloud top is essential for later snow growth processes. We have compared the cloud top temperature (CTT) observations in the five campaigns as shown in Fig. 2. The cloud top was determined using the highest radar echo as mostly detected by the Ka band radar. In METRICs, the thick rain layer leads to decent attenuation at Ka-band and the C-band radar occasionally observed the highest cloud top . In presence of multi-layer clouds, we have removed upper clouds with a radar echo gap exceeding 2 km for excluding the impact of seeding . As shown in Fig. a, clouds over tropics are characterized by the coldest tops, while cloud tops over high-latitudes are shallower and warmer, which are in line with Cloudsat observations . Note that due to the degradation of C-band radar sensitivity with range, the cloud top radar reflectivity in METRICs is of the magnitude of -10 to 0 dBZ, suggesting that the actual cloud top temperatures were even colder.
At temperatures colder than -70 °C, a colder cloud top is associated with more ice populations , and the prevalence of CTT between -70 and -60 °C in METRICs is expected to generate numerous pristine ice particles. In addition, not like the stratiform precipitation in baroclinic cyclones and fronts in mid-latitudes, the stratiform rainfall over tropics is closely associated with nearby convective cells . The ice particles in stratiform region may originate from the convection at younger and more vigorous stages .
Distributions of cloud top (a) temperature and (b) corresponding calibrated radar reflectivity in the five field campaigns. The frequency of riming as defined as FR ≥ 0.5 as a function of in-cloud temperature is presented in (c).
have shown that riming occurrence increases with temperature with seasonal variations. As shown in Fig. c, our statistics in general agree with their conclusions, except for rather limited riming signatures in AWARE. In addition, many rimed snow cases in IMPACTS have a surface temperature of -8 to -6 °C and thus less observations with the in-cloud temperature warmer than -6 °C, which may explain the decreased riming frequency warmer than -6 °C. Snow observations made in METRICs and IMPACTS show more frequent riming, while riming is less frequent in TRIPEx-pol and BAECC.
Dual-frequency signaturesDependence on reflectivity
Statistical results of DWRX∕C,Ka and DWRKa,W as a function of Z during stratiform precipitation are given in Fig. a, b. As expected, DWR generally increases with Z after certain thresholds below which the snow dimensions are so small that the Rayleigh scattering condition is met. As ZX/C increases from 0 dBZ to about 10 dBZ, the median values of DWRX∕C,Ka increase from 0 dB to around 1 dB (Fig. a). The interquartile ranges of DWRX∕C,Ka in all campaigns except for METRICs show great variability after ZX/C exceeds 10 dBZ, which is explained by the increased complexity of non-Rayleigh scattering at Ka-band as the snow size increases. In METRICs, the interquartile ranges of DWRX∕C,Ka remain stable near 0 dB for ZX/C values below 18 dBZ. In addition, median DWRX∕C,Ka values from both IMPACTS and METRICs do not exceed 1 dB for ZX/C<18 dBZ, implying the latitude dependence of DWR-ZX/C relation.
Observed DWR (circles) as a function of (a)ZX/C and (b)ZKa. The boxplots represent the median (horizontal line) and interquartile range (box) of DWR values within a reflectivity interval of 2 dB. Blue, purple, red, green and yellow dashed curves represent DWR fits for METRICs, IMPACTS, TRIPEx-pol, BAECC and AWARE, respectively. Black dot-dashed and dashed black curves represent DWR fits for snow over Oliktok Point (70.4958° N, 149.8868° W) in October 2016 and May 2017 as made by , respectively. The DWR fits for expressions of DWR =10a⋅(Z+20)b in are given in Table .
Similar dependence can also be found for DWRKa,W. As shown in Fig. b, DWRKa,W starts exceeding 1 dB at ZKa≈ 0, 2, 4, 8 and 12 dBZ in AWARE, BAECC, TRIPEx-pol, IMPACTS, and METRICs, respectively. We found that DWRKa,W statistics of BAECC and AWARE are similar to those observed at the Oliktok Point, Alaska , implying similar cloud microphysics over high-latitudes and polar regions. Comparing DWR observations from different campaigns, it appears that the radar reflectivity values over which the non-Rayleigh scattering starts appearing at Ka (Fig. a) and W (Fig. b) bands decrease with the increase of latitude. This latitude dependence appears to be associated with cloud top temperature (Fig. a). Namely, non-Rayleigh scattering signatures appear at a larger radar reflectivity threshold for a colder cloud top.
The average cloud top over low-latitudes is above 15 km, but is below 6 km at high-latitudes . Namely, the climatology of cloud top temperature increases towards polar regions. Given much colder cloud tops corresponding to more INPs and the ice transport from local convections over low-latitudes, ice number concentrations are expected to decrease with latitude. In contrast, significant increase of DWR(Ka,W) occurs below the dendritic growth zone across all campaigns (Fig. 4b), indicating that in-cloud temperature is more important than ice number concentration for snow aggregation. Hence, larger ZX/C(ZKa) thresholds for Ka- (W-) band non-Rayleigh scattering at METRICs and IMPACTS may be explained by more ice particles over low latitudes.
Fitted parameters for DWR=10a⋅(Z+20)b in Fig. 2.
Fitting ParametersabDWRX∕C,KaMETRICs-25.1815.54IMPACTS-9.3345.835TRIPEx-pol-8.6465.732BAECC-7.985.346AWARE-9.9636.684DWRKa,WMETRICs-10.446.956IMPACTS-6.3714.387TRIPEx-pol-5.3123.826BAECC-5.0513.73AWARE-4.183.291Dependence on temperature
, using combined active and passive satellite sensors, have shown that the effective diameter of cloud top ice increases with temperature, and the ice size is generally smaller over lower latitudes. However, snow growth characteristics beneath cloud tops over different latitudes are still poorly understood. As shown in Fig. 4a, DWRX∕C,Ka values (≥ 3 dB) are concentrated at temperatures exceeding -12 °C. More significant temperature influence is found for DWRKa,W (Fig. 4b), which can be explained by a pronounced difference in saturation vapor pressure between ice and liquid and a rapid aggregation favored by the dendritic features of snowflakes at around -15 °C. This implies that the unique temperature zone around -15 °C may be essential for initiating large snowflakes regardless of the regional climate.
Observed (circles) (a) DWRX∕C,Ka and (b) DWRKa,W(c)ZX/C as a function of temperature. The boxplots represent the median (horizontal line) and interquartile range (box) of the observations within a temperature interval of 4 °C.
At temperatures warmer than -12 °C, the rate of DWR increase with temperature is more pronounced, and ZX/C shows similar temperature dependence expect for AWARE. The enhanced ice growth is attributed to the thickened quasi-liquid layer on the ice particle surface, and therefore more efficient aggregation . In AWARE, the absence of ZX/C increase could be explained by the near-surface sublimation which lowers ZX/C but does not necessarily lead to decrease of DWR .
Observed DWRX∕C,Ka (circles) as a function of ZX/C at temperatures greater than -8 °C in (a) METRICs, (b) IMPACTS, (c) TRIPEx-pol, (d) BAECC, and (e) AWARE. The colored isolines represent the observation density distributions for each field campaign. In (f), the black dashed curve marks the fit for FY-3G Ku and Ka band radar observations of snow over the melting layer of low- to mid-latitude stratiform rainfall (Liu et al., 2024).
A zoom-in view into the ZX/C-DWRX∕C,Ka space for temperatures warmer than -8 °C (Fig. ) reveals an obvious positive correlation between DWRX∕C,Ka and ZX/C. For a given ZX/C at 20 dBZ, DWRX∕C,Ka in METRICs and IMPACTS is mostly 0–2 dB, comparing to 2–4 dB in TRIPEx-pol and BAECC, suggesting that a high ZX/C does not necessarily mean large-sized snowflakes. Snow observations with ZX/C≥ 20 dBZ are prevalent in METRICs and IMPACTS, but they are rather rare in TRIPEx-pol, BAECC and AWARE. Recently, using FY-3G spaceborne radar observations of snowflakes above melting layers of low- to mid-latitude stratiform rainfall showed that DWRKu,Ka stably increases with ZKu, and their fit (black dashed curve in Fig. 4) actually agrees well with the ground-based radar observations made in METRICs and IMPACTS. Given cloud top temperatures in METRICs and IMPACTS are lower than -70 to -50 °C (Fig. a), the fit in appears to be representative of low- to mid-latitude deep precipitating snowfall events.
Triple-frequency signatures
Observed DWRX∕C,Ka-DWRKa,W occurrence and FR for temperature ranges of -10–0 °C are given in Fig. . As the temperature rises from -20 to -10 °C (gray isolines) to -10 to 0° (white isolines), the increases of DWRX∕C,Ka are about 2–6 dB in METRICs, IMPACTS, TRIPEx-Pol and BAECC. In contrast, the weak snow growth in AWARE may be because of the katabatic winds which lower relative humidity and the sparse amount of INPs over McMurdo .
Observed DWRX∕C,Ka-DWRKa,W occurrence and FR within the temperature range of -10–0 °C from (a) METRICs, (b) TRIPEx-pol, (c) BAECC, and (d) AWARE. The DWR planes are overlapped by scattering model curves of rimed dendrite aggregates with effective liquid water paths of 0.1 kg m−2 (FR ≈ 0), 0.5 kg m−2 (FR ≈ 0.53), and 1 kg m−2 (FR ≈ 0.72) (L15_0.1, L15_0.5, and L15_1.0, respectively) as adapted from . The gray and white isolines represent the DWRX∕C,Ka-DWRKa,W occurrence density of 5 % within the temperature of -20 to -10 and -10 to 0 °C, respectively. The fitting curves are marked with different colors.
In METRICs, IMPACTS and TRIPEx-pol, 5 % of DWRX∕C,Ka observations between -10 and 0 °C (white isolines) extend to above 5–6 dB. Specifically, triple-frequency observations in IMPACTS generally follow the scattering model assuming the FR of 0.53–0.72 and the velocity-based FR estimates are on the order of 0.4–0.5, presenting obvious riming signatures. The observed DWRX∕C,Ka-DWRKa,W signatures in METRICs and IMPACTS follow the scattering model assuming the FR =0, while the velocity-based FR in METRICs is larger than TRIPEx-pol. Since the uncorrected supercooled liquid water attenuation does not affect the observed Doppler velocity but can lead to underestimated DWRKa,W, It is anticipated that the observed DWRKa,W in METRICs is underestimated. The highest DWRX∕C,Ka observations and least riming signatures were found for TRIPEx-pol, implying that the absence of riming is favorable for snow aggregation . To compare the triple-frequency signatures in these campaigns, power-law fits were made to the median values of DWRX∕C,Ka in each DWRX∕C,Ka interval. As shown in Fig. , the prefactor in AWARE is the smallest, and is the largest in TRIPEx-pol.
Although the velocity-based FR is qualitatively consistent with triple-frequency signatures of riming, FR should be used with caution. have shown that large uncertainties exist for FR <0.5, which is actually the reason why we quantified riming occurrence using FR >0.5 in Fig. 2c. We had proposed a dual-frequency approach for riming classification , while the triple-frequency observations are more favorable for quantitative retrieval of snow microphysics . In addition, the use of third frequency (X- or C-band) narrows the retrieval of snow size distributions. Although the exponential function for snow size distribution can be assumed in climatological analysis, the third frequency is needed to inform the parameters in Gamma function which is more consistent to observations .
Conclusion and Discussion
In this work, we examined snow microphysics over various geographies using triple-frequency radar observations from METRICs, TRIPEx-pol, BAECC, IMPACTS and AWARE. Our results suggest the promise of using long-term triple-frequency setup for understanding the climatology of ice formation, snow aggregation and riming processes. Our major conclusions are conceptualized in Fig. and summarized below,
Conceptual diagram of major conclusions in this study. Based on the statistics in Fig. 4c and current spaceborne radar sensitivities, we sketched the detectable temperature ranges above the melting layer over different regions with colored bars.
Majority of snow in BAECC and AWARE (high latitudes) originates from a cloud top temperature >-40 °C, and therefore the heterogeneous nucleation plays an essential role in snow formation. In contrast, precipitating clouds in other campaigns (low- to mid-latitudes) are much deeper, and the coldest cloud top was found in METRICs.
Following the conceptual model given by , our analysis suggests that the triple-frequency signatures of riming generally collaborate with the velocity-based approach, and major snow growth characteristics are dependent on regional climate. With the in-cloud Doppler measurement capability, long-term EarthCARE observations may be used for globally mapping riming signatures. The velocity-based FR estimates suggests increased frequency of snow riming with in-cloud temperature. The heaviest riming in IMPACTS, corroborating with the expected riming signatures in DWRX∕C,Ka-DWRKa,W space. In contrast, riming signatures in TRIPEx-pol and AWARE is rather weak.
DWR observations in these field campaigns qualitatively indicate the dendritic growth zone around -15° playing a key role in initiating enhanced snow size growth, and reveals a generally temperature-dependent snowflake growth characteristics. Specifically, TRIPEx-pol presents a typical “hook” signature with FR =0, indicative of favorable aggregation in absence of riming. The weakest snow growth is found for AWARE, potentially owing to the scarcity of ice nucleating particles and avaliable water vapor for deposition.
Our statistics are also indicative of the capability of current spaceborne radars in snow detection. Compared to the high sensitivity of EarthCARE W-band radar, the sensitivity of GPM-CO and FY-3G radars is on the magnitude of 10–20 dBZ. Namely, GPM-CO and FY-3G radars well detect low- to mid-latitude snow at temperatures warmer than -12 °C but most high-latitude snow cannot be observed (Fig. ).
While the promise of using multi-month triple-frequency radar observations is presented in this study, the analyzed field campaigns were conducted in different seasons. In addition, the datasets generated from relatively short observation periods are not adequate for examining seasonal variations. We advocate the sustained support for a triple-frequency super site, which would facilitate multi-year observations and greatly advance our climatological understanding of snow microphysics.
On the other hand, the minimum reflectivity where the Ka-band non-Rayleigh scattering appears over low- to mid-latitudes (METRICs and IMPACTS) is about 15–20 dBZ (ZX/C) which is well detectable by FY-3G Ku-(12 dBZ sensitivity) and Ka-band (8 dBZ sensitivity) radars . This suggests that snow above the melting layer if is well identified could be served as Rayleigh-scattering targets at both Ku and Ka bands for inter-calibration of ground-based S-band weather radars and spaceborne radars without considering the attenuation effects.
In the last decade, ground-based/airborne triple-frequency radar field campaigns have demonstrated the unprecedented value of triple-frequency radars in unraveling snow microphysics in various geographics. The recent success of METRICs campaign allowed us to do synergetic analysis of the snow growth climatology from polar regions to tropics, and significant geographical differences have been identified. Thanks to the high-sensitivity of well-matched ground-based/airborne triple-frequency radars, the snow growth processes with rather weak radar echoes are detectable. In contrast, spaceborne radars can also provide very unique W- (CloudSat, EarthCARE), Ka- (GPM, FY-3G) and Ku-band (GPM, FY-3G) coincidence observations , while the minimal detectable signal is on the order of 13–19 dBZ for GPM Ku/Ka-band radars , and 10–13 dBZ for FY-3G Ku/Ka-band radars . Namely, the majority of snow growth signatures (Figs. , ) as observed in the five campaigns are missed in these dataset. For the future Ku-Ka-W band spaceborne mission, assumed a minimal sensitivity of 0–5 dBZ which seems to be a good selection for making use of the non-Rayleigh scattering characteristics (Fig. ) and majority of solid precipitation at temperatures warmer than -10 °C can be detected (Fig. c).
Data availability
The used data including METRICs in this paper is available at 10.5281/zenodo.14575727. TRIPEx-pol observations can be accessed at 10.5281/zenodo.5025636. Observations made in BAECC and AWARE are available from ARM data center (https://www.arm.gov/data/, last access: December 2025). IMPACTS observations can be accessed at 10.5067/IMPACTS/DATA101.
Author contributions
HL and XS conceptualized the study. QL and HL performed the experiment and wrote the paper. YZ, WL, ZR, LL, AL, and CZ conducted the METRICs campaign and collected the radar observations. All the authors took part in the interpretation of the results and edits of the paper.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
We thank Stefan Kneifel, Leonie Von Terzi, and Maximilian Maahn for very helpful comments on this work.
Financial support
This research has been supported by Open Grants of the High Impact Weather Key Laboratory (special), China Meteorological Administration (grant no. 2024-K-01), State Key Laboratory of Severe Weather Meteorological Science and Technology (grant no. 2025QZA04), National Natural Science Foundation of China (grant no. 42475095), Science Funds of Changdao National Climatic Observatory (grant no. 2024cdkfz01), Basic Research Fund of CAMS (grant no. 2023Z008), Anhui Provincial Natural Science Foundation (grant no. 2408055UQ007), Shandong Provincial Natural Science Foundation (grant no. ZR2025LQX001), Open Grants of the China Meteorological Administration Aerosol‐Cloud and Precipitation Key Laboratory (grant no. KDW2412), and Alexander von Humboldt Foundation.
Review statement
This paper was edited by Timothy Garrett and reviewed by two anonymous referees.
References Anagnostou, E. N., Morales, C. A., and Dinku, T.: The use of TRMM precipitation radar observations in determining ground radar calibration biases, Journal of Atmospheric and Oceanic Technology, 18, 616–628, 2001.Barrett, A. I., Westbrook, C. D., Nicol, J. C., and Stein, T. H. M.: Rapid ice aggregation process revealed through triple-wavelength Doppler spectrum radar analysis, Atmospheric Chemistry and Physics, 19, 5753–5769, 10.5194/acp-19-5753-2019, 2019.Battaglia, A., Kollias, P., Dhillon, R., Roy, R., Tanelli, S., Lamer, K., Grecu, M., Lebsock, M., Watters, D., Mroz, K., Heymsfield, G., Li, L., and Furukawa, K.: Spaceborne cloud and precipitation radars: Status, challenges, and ways forward, Reviews of Geophysics, 58, e2019RG000686, 10.1029/2019RG000686, 2020.Bertrand, L., Kay, J. E., Haynes, J., and de Boer, G.: A global gridded dataset for cloud vertical structure from combined CloudSat and CALIPSO observations, Earth System Science Data, 16, 1301–1316, 10.5194/essd-16-1301-2024, 2024. Blaylock, B. K., Horel, J. D., and Liston, S. T.: Cloud archiving and data mining of High-Resolution Rapid Refresh forecast model output, Computers & Geosciences, 109, 43–50, 2017. Brandes, E. A., Ikeda, K., Zhang, G., Schönhuber, M., and Rasmussen, R. M.: A statistical and physical description of hydrometeor distributions in Colorado snowstorms using a video disdrometer, Journal of applied meteorology and climatology, 46, 634–650, 2007.Bringi, V. N. and Chandrasekar, V.: Polarimetric Doppler weather radar: principles and applications, Cambridge university press, 10.1016/S0169-8095(02)00024-8, 2001.Chase, R. J.: Improving active remote sensing retrievals of snowfall at microwave wavelengths: an emphasis on the global precipitation measurement mission’s dual-frequency precipitation radar, Ph.D. thesis, University of Illinois at Urbana-Champaign, 10.13016/M2J01X, 2021. Chase, R. J., Finlon, J. A., Borque, P., McFarquhar, G. M., Nesbitt, S. W., Tanelli, S., Sy, O. O., Durden, S. L., and Poellot, M. R.: Evaluation of Triple-Frequency Radar Retrieval of Snowfall Properties Using Coincident Airborne In Situ Observations During OLYMPEX, Geophysical Research Letters, 45, 5752–5760, 2018.Chellini, G., Gierens, R., and Kneifel, S.: Ice Aggregation in Low-Level Mixed-Phase Clouds at a High Arctic Site: Enhanced by Dendritic Growth and Absent Close to the Melting Level, Journal of Geophysical Research: Atmospheres, 127, e2022JD036860, 10.1029/2022JD036860, 2022. 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, Proceedings of the National Academy of Sciences, 107, 11217–11222, 2010.Dias Neto, J., Kneifel, S., Ori, D., Trömel, S., Handwerker, J., Bohn, B., Hermes, N., Mühlbauer, K., Lenefer, M., and Simmer, C.: The TRIple-frequency and Polarimetric radar Experiment for improving process observations of winter precipitation, Earth System Science Data, 11, 845–863, 10.5194/essd-11-845-2019, 2019.Ding, H., Li, H., and Liu, L.: Improved spectral processing for a multi-mode pulse compression Ka–Ku-band cloud radar system, Atmospheric Measurement Techniques, 15, 6181–6200, 10.5194/amt-15-6181-2022, 2022. Dolinar, E. K., Campbell, J. R., Marquis, J. W., Garnier, A. E., and Karpowicz, B. M.: Novel parameterization of ice cloud effective diameter from collocated CALIOP-IIR and CloudSat retrievals, Journal of Applied Meteorology and Climatology, 61, 891–907, 2022.Dunnavan, E. L., Carlin, J. T., Schvartzman, D., Ryzhkov, A. V., Bluestein, H., Emmerson, S., McFarquhar, G. M., Heymsfield, G. M., and Yorks, J.: High-Resolution Snowstorm Measurements and Retrievals Using Cross-Platform Multi-Frequency and Polarimetric Radars, Geophysical Research Letters, 50, e2023GL103692, 10.1029/2023GL103692, 2023.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, Atmospheric Measurement Techniques, 11, 3059–3079, 10.5194/amt-11-3059-2018, 2018. Field, P. R. and Heymsfield, A. J.: Importance of snow to global precipitation, Geophysical Research Letters, 42, 9512–9520, 2015.Fitch, K. E. and Garrett, T. J.: Graupel precipitating from thin Arctic clouds with liquid water paths less than 50 g m−2, Geophysical Research Letters, 49, e2021GL094075, 10.1029/2021GL094075, 2022. Grazioli, J., Madeleine, J.-B., Gallée, H., Forbes, R. M., Genthon, C., Krinner, G., and Berne, A.: Katabatic winds diminish precipitation contribution to the Antarctic ice mass balance, Proceedings of the National Academy of Sciences, 114, 10858–10863, 2017.Gultepe, I., Isaac, G., Hudak, D., Nissen, R., and Strapp, J. W.: Dynamical and Microphysical Characteristics of Arctic Clouds during BASE, Journal of Climate, 13, 1225–1254, 10.1175/1520-0442(2000)013<1225:DAMCOA>2.0.CO;2, 2000. Heymsfield, A., Bansemer, A., Heymsfield, G., Noone, D., Grecu, M., and Toohey, D.: Relationship of multiwavelength radar measurements to ice microphysics from the IMPACTS field program, Journal of Applied Meteorology and Climatology, 62, 289–315, 2023. 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, Journal of the Atmospheric Sciences, 64, 1047–1067, 2007. Heymsfield, A. J., Field, P., and Bansemer, A.: Exponential size distributions for snow, Journal of the Atmospheric Sciences, 65, 4017–4031, 2008.Heymsfield, A. J., Schmitt, C., and Bansemer, A.: Ice cloud particle size distributions and pressure-dependent terminal velocities from in situ observations at temperatures from 0 to -86 C, Journal of the Atmospheric Sciences, 70, 4123–4154, 2013. Heymsfield, A. J., Schmitt, C. G., Chen, C. C. J., Bansemer, A. R., Gettelman, A., Field, P. R., and Liu, C.: Contributions of the liquid and ice phases to global surface precipitation: Observations and global climate modeling, Journal of the Atmospheric Sciences, 77, 2629–2648, 2020.Hines, K. M., Bromwich, D. H., Silber, I., Russell, L. M., and Bai, L.: Predicting Frigid Mixed-Phase Clouds for Pristine Coastal Antarctica, Journal of Geophysical Research: Atmospheres, 126, e2021JD035112, 0.1029/2021JD035112, 2021. Hogan, R. J., Illingworth, A. J., and Sauvageot, H.: Measuring crystal size in cirrus using 35-and 94 GHz radars, Journal of Atmospheric and Oceanic Technology, 17, 27–37, 2000. Hooper, J. and Kippax, A.: Radar echoes from meteorological precipitation, Proceedings of the IEE-Part I: General, 97, 89–95, 1950. Hou, A. Y., Kakar, R. K., Neeck, S., Azarbarzin, A. A., Kummerow, C. D., Kojima, M., Oki, R., Nakamura, K., and Iguchi, T.: The global precipitation measurement mission, Bulletin of the American meteorological Society, 95, 701–722, 2014. Houze Jr., R. A.: Stratiform precipitation in regions of convection: A meteorological paradox?, Bulletin of the American Meteorological Society, 78, 2179–2196, 1997. Illingworth, A. J., Barker, H. W., Beljaars, A. C. M., Ceccaldi, M., Chepfer, H., Clerbaux, N., Cole, J. N. S., Delanoë, J. M. E., Doménech, C., Donovan, D. P., Fukuda, S., Hirakata, M., Hogan, R. J., Huenerbein, A., Kollias, P., Kubota, T., Nakajima, T., Nakajima, T. Y., Nishizawa, T., Ohno, Y., Okamoto, H., Oki, R., Sato, K., Satoh, M., Shephard, M. W., Velázquez-Blazquez, A., Wandinger, U., Wehr, T., and van Zadelhoff, G. J.: The EarthCARE satellite: The next step forward in global measurements of clouds, aerosols, precipitation, and radiation, Bulletin of the American Meteorological Society, 96, 1311–1332, 2015. 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.Karrer, M., Seifert, A., Ori, D., and Kneifel, S.: Improving the representation of aggregation in a two-moment microphysical scheme with statistics of multi-frequency Doppler radar observations, Atmospheric Chemistry and Physics, 21, 17133–17166, 10.5194/acp-21-17133-2021, 2021. Kneifel, S. and Moisseev, D.: Long-Term Statistics of Riming in Nonconvective Clouds Derived from Ground-Based Doppler Cloud Radar Observations, Journal of the Atmospheric Sciences, 77, 3495–3508, 2020.Kneifel, S., Kulie, M., and Bennartz, R.: A triple-frequency approach to retrieve microphysical snowfall parameters, Journal of Geophysical Research: Atmospheres, 116, 10.1029/2010JD015430, 2011. Kneifel, S., Lerber, A., Tiira, J., Moisseev, D., Kollias, P., and Leinonen, J.: Observed relations between snowfall microphysics and triple-frequency radar measurements, Journal of Geophysical Research: Atmospheres, 120, 6034–6055, 2015. Kulie, M. S., Hiley, M. J., Bennartz, R., Kneifel, S., and Tanelli, S.: Triple-frequency radar reflectivity signatures of snow: Observations and comparisons with theoretical ice particle scattering models, Journal of Applied Meteorology and Climatology, 53, 1080–1098, 2014. Kummerow, C., Barnes, W., Kozu, T., Shiue, J., and Simpson, J.: The tropical rainfall measuring mission (TRMM) sensor package, Journal of Atmospheric and Oceanic Technology, 15, 809–817, 1998. Leinonen, J. and Moisseev, D.: What do triple-frequency radar signatures reveal about aggregate snowflakes?, Journal of Geophysical Research: Atmospheres, 120, 229–239, 2015. Leinonen, J. and Szyrmer, W.: Radar signatures of snowflake riming: A modeling study, Earth and Space Science, 2, 346–358, 2015.Leinonen, J., Kneifel, S., Moisseev, D., Tyynelä, J., Tanelli, S., and Nousiainen, T.: Evidence of nonspheroidal behavior in millimeter-wavelength radar observations of snowfall, Journal of Geophysical Research: Atmospheres, 117, 10.1029/2012JD017680, 2012.Leinonen, J., Lebsock, M. D., Tanelli, S., Suzuki, K., Yashiro, H., and Miyamoto, Y.: Performance assessment of a triple-frequency spaceborne cloud–precipitation radar concept using a global cloud-resolving model, Atmospheric Measurement Techniques, 8, 3493–3517, 10.5194/amt-8-3493-2015, 2015.Leinonen, J., Lebsock, M. D., Tanelli, S., Sy, O. O., Dolan, B., Chase, R. J., Finlon, J. A., von Lerber, A., and Moisseev, D.: Retrieval of snowflake microphysical properties from multifrequency radar observations, Atmospheric Measurement Techniques, 11, 5471–5488, 10.5194/amt-11-5471-2018, 2018.Li, H. and Moisseev, D.: Melting layer attenuation at Ka- and W-bands as derived from multi-frequency radar Doppler spectra observations, Journal of Geophysical Research: Atmospheres, 124, 10.1029/2019JD030316, 2019.Li, H. and Moisseev, D.: Two layers of melting ice particles within a single radar bright band: Interpretation and implications, Geophysical Research Letters, 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?, Journal of Geophysical Research: Atmospheres, 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, Atmospheric Chemistry and Physics, 20, 9547–9562, 10.5194/acp-20-9547-2020, 2020.Li, H., Möhler, O., Petäjä, T., and Moisseev, D.: Two-year statistics of columnar-ice production in stratiform clouds over Hyytiälä, Finland: environmental conditions and the relevance to secondary ice production, Atmospheric Chemistry and Physics, 21, 14671–14686, 10.5194/acp-21-14671-2021, 2021.Li, H., Moisseev, D., Luo, Y., Liu, L., Ruan, Z., Cui, L., and Bao, X.: Assessing specific differential phase (KDP)-based quantitative precipitation estimation for the record- breaking rainfall over Zhengzhou city on 20 July 2021, Hydrology and Earth System Sciences, 27, 1033–1046, 10.5194/hess-27-1033-2023, 2023a.Li, Q.: A survey of snow growth signatures from tropics to Antarctica using triple-frequency radar observations, Zenodo [data set], 10.5281/zenodo.14575727, 2024.Li, Q., Li, H., Sun, X., Ruan, Z., Liu, L., Gao, W., Liao, C., and Ding, H.: On the quantification of thermodynamic phases of raining clouds: Insights from multi-year CloudSat and ground-based radar observations over Longmen, Southern China, Journal of Geophysical Research: Atmospheres, 130, e2024JD041824, 10.1029/2024JD041824, 2025. Li, Z., Wen, Y., Liao, L., Wolff, D., Meneghini, R., and Schuur, T.: Joint collaboration on comparing NOAA’s ground-based weather radar and NASA–JAXA’s spaceborne radar, Bulletin of the American Meteorological Society, 104, E1435–E1451, 2023b. Liebe, H. J.: An updated model for millimeter wave propagation in moist air, Radio Science, 20, 1069–1089, 1985. Lin, Y., Donner, L. J., and Colle, B. A.: Parameterization of riming intensity and its impact on ice fall speed using ARM data, Monthly Weather Review, 139, 1036–1047, 2011.Liu, B., Li, H., Liu, L., Shang, J., Yin, H., and Kuo, K.-S.: Sea Surface and Snowflakes as Natural Targets Connecting FY-3G and GPM-CO Dual-Frequency Radars, Geophysical Research Letters, 51, e2024GL110878, 10.1029/2024GL110878, 2024.Liu, B., Li, H., Liu, L., Shang, J., Kuo, K.-S., Lu, C., Yuan, M., and Jiang, B.: On the detection sensitivities of dual-frequency radars onboard FY-3G and GPM-CO, Atmospheric Research, 107935, 10.1016/j.atmosres.2025.107935, 2025. Liu, Y., Yau, M.-K., Shima, S.-i., Lu, C., and Chen, S.: Parameterization and explicit modeling of cloud microphysics: Approaches, challenges, and future directions, Advances in Atmospheric Sciences, 40, 747–790, 2023. Löhnert, U., Schween, J. H., Acquistapace, C., Ebell, K., Maahn, M., Barrera-Verdejo, M., Hirsikko, A., Bohn, B., Knaps, A., O’Connor, E., Simmer, C., Wahner, A., and Crewell, S.: JOYCE: Jülich observatory for cloud evolution, Bulletin of the American Meteorological Society, 96, 1157–1174, 2015. Lubin, D., Zhang, D., Silber, I., Scott, R. C., Kalogeras, P., Battaglia, A., Bromwich, D. H., Cadeddu, M., Eloranta, E., Fridlind, A., Frossard, A., Hines, K. M., Kneifel, S., Leaitch, W. R., Lin, W., Nicolas, J., Powers, H., Quinn, P. K., Rowe, P., Russell, L. M., Sharma, S., Verlinde, J., and Vogelmann, A. M.: AWARE: The atmospheric radiation measurement (ARM) west Antarctic radiation experiment, Bulletin of the American Meteorological Society, 101, E1069–E1091, 2020.Maherndl, N., Moser, M., Lucke, J., Mech, M., Risse, N., Schirmacher, I., and Maahn, M.: Quantifying riming from airborne data during the HALO-(AC)3 campaign, Atmospheric Measurement Techniques, 17, 1475–1495, 10.5194/amt-17-1475-2024, 2024. Masaki, T., Iguchi, T., Kanemaru, K., Furukawa, K., Yoshida, N., Kubota, T., and Oki, R.: Calibration of the dual-frequency precipitation radar onboard the Global Precipitation Measurement Core Observatory, IEEE Transactions on Geoscience and Remote Sensing, 60, 1–16, 2020. Mason, S., Chiu, C., Hogan, R., Moisseev, D., and Kneifel, S.: Retrievals of Riming and Snow Density From Vertically Pointing Doppler Radars, Journal of Geophysical Research: Atmospheres, 123, 13–807, 2018.Mason, S. L., Hogan, R. J., Westbrook, C. D., Kneifel, S., Moisseev, D., and von Terzi, L.: The importance of particle size distribution and internal structure for triple-frequency radar retrievals of the morphology of snow, Atmospheric Measurement Techniques, 12, 4993–5018, 10.5194/amt-12-4993-2019, 2019. Matrosov, S.: Possibilities of cirrus particle sizing from dual-frequency radar measurements, Journal of Geophysical Research: Atmospheres, 98, 20675–20683, 1993. Matrosov, S. Y.: A dual-wavelength radar method to measure snowfall rate, Journal of Applied Meteorology, 37, 1510–1521, 1998. Matrosov, S. Y., Maahn, M., and De Boer, G.: Observational and modeling study of ice hydrometeor radar dual-wavelength ratios, Journal of Applied Meteorology and Climatology, 58, 2005–2017, 2019. Maynard, R. H.: Radar and weather, Journal of Atmospheric Sciences, 2, 214–226, 1945.McMurdie, L. A., Heymsfield, G., Yorks, J. E., and Braun, S. A.: Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) Collection, NASA EOSDIS Global Hydrometeorology Resource Center Distributed Active Archive Center [data set], Huntsville, Alabama, U.S.A. 10.5067/IMPACTS/DATA101, 2019. McMurdie, L. A., Heymsfield, G. M., Yorks, J. E., Braun, S. A., Skofronick-Jackson, G., Rauber, R. M., Yuter, S., Colle, B., McFarquhar, G. M., Poellot, M., Novak, D. R., Lang, T. J., Kroodsma, R., McLinden, M., Oue, M., Kollias, P., Kumjian, M. R., Greybush, S. J., Heymsfield, A. J., Finlon, J. A., McDonald, V. L., and Nicholls, S.: Chasing snowstorms: The investigation of microphysics and precipitation for Atlantic coast-threatening snowstorms (IMPACTS) campaign, Bulletin of the American Meteorological Society, 103, E1243–E1269, 2022. Moisseev, D., von Lerber, A., and Tiira, J.: Quantifying the effect of riming on snowfall using ground-based observations, Journal of Geophysical Research: Atmospheres, 122, 4019–4037, 2017. Morrison, H. and Milbrandt, J. A.: Parameterization of cloud microphysics based on the prediction of bulk ice particle properties. Part I: Scheme description and idealized tests, Journal of the Atmospheric Sciences, 72, 287–311, 2015. Morrison, H., De Boer, G., Feingold, G., Harrington, J., Shupe, M. D., and Sulia, K.: Resilience of persistent Arctic mixed-phase clouds, Nature Geoscience, 5, 11–17, 2012.Morrison, H., van Lier-Walqui, M., Fridlind, A. M., Grabowski, W. W., Harrington, J. Y., Hoose, C., Korolev, A., Kumjian, M. R., Milbrandt, J. A., Pawlowska, H., Posselt, D. J., Prat, O. P., Reimel, K. J., Shima, S.-I., van Diedenhoven, B., and Xue, L.: Confronting the challenge of modeling cloud and precipitation microphysics, Journal of advances in modeling earth systems, 12, e2019MS001689, 10.1029/2019MS001689, 2020. Mosimann, L.: An improved method for determining the degree of snow crystal riming by vertical Doppler radar, Atmospheric research, 37, 305–323, 1995.Mroz, K., Battaglia, A., Nguyen, C., Heymsfield, A., Protat, A., and Wolde, M.: Triple-frequency radar retrieval of microphysical properties of snow, Atmospheric Measurement Techniques, 14, 7243–7254, 10.5194/amt-14-7243-2021, 2021. Mülmenstädt, J., Sourdeval, O., Delanoë, J., and Quaas, J.: Frequency of occurrence of rain from liquid-, mixed-, and ice-phase clouds derived from A-Train satellite retrievals, Geophysical Research Letters, 42, 6502–6509, 2015.Nguyen, C. M., Wolde, M., Battaglia, A., Nichman, L., Bliankinshtein, N., Haimov, S., Bala, K., and Schuettemeyer, D.: Coincident in situ and triple-frequency radar airborne observations in the Arctic, Atmospheric Measurement Techniques, 15, 775–795, 10.5194/amt-15-775-2022, 2022. Pang, S., Ruan, Z., Yang, L., Liu, X., Huo, Z., Li, F., and Ge, R.: Estimating raindrop size distributions and vertical air motions with spectral difference using vertically pointing radar, Journal of Atmospheric and Oceanic Technology, 38, 1697–1713, 2021. Asmi, E., Brus, D., Bäck, J., Chandrasekar, V., Engelmann, R., Hakola, H., Hickmon, N., Komppula, M., Krüger, O., Kuang, C., Laurila, T., Levula, J., Manninen, A. J., Massoli, P., Moisseev, D., Mona, L., Nicoll, K., O’Connor, E. J., Petersen, W., and Petäjä, T.: BAECC: A field campaign to elucidate the impact of biogenic aerosols on clouds and climate, Bulletin of the American Meteorological Society, 97, 1909–1928, 2016. Protat, A., Rauniyar, S., Delanoë, J., Fontaine, E., and Schwarzenboeck, A.: W-band (95 GHz) radar attenuation in tropical stratiform ice anvils, Journal of Atmospheric and Oceanic Technology, 36, 1463–1476, 2019. Ruan, Z., Ming, H., Ma, J., Ge, R., and Bian, L.: Analysis of the microphysical properties of a stratiform rain event using an L-Band profiler radar, Journal of Meteorological Research, 28, 268–280, 2014.Sassen, K. and Wang, Z.: Classifying clouds around the globe with the CloudSat radar: 1-year of results, Geophysical Research Letters, 35, 10.1029/2007GL032591, 2008. Schumacher, C. and Houze Jr., R. A.: Stratiform rain in the tropics as seen by the TRMM precipitation radar, Journal of Climate, 16, 1739–1756, 2003. Seifert, P., Ansmann, A., Mattis, I., Althausen, D., Tesche, M., Wandinger, U., Muller, D., and Pérez, C.: Lidar-based profiling of the tropospheric cloud-ice distribution to 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, HWJ, and Monna, WAA, Delft, the Netherlands, ISBN 978-90-6960-233-2, 2009.Song, J. I., Yum, S. S., Park, S.-H., Kim, K.-H., Park, K.-J., and Joo, S.-W.: Climatology of melting layer heights estimated from cloud radar observations at various locations, Journal of Geophysical Research: Atmospheres, 126, e2021JD034816, 10.1029/2021JD034816, 2021. Stein, T. H., Westbrook, C. D., and Nicol, J.: Fractal geometry of aggregate snowflakes revealed by triple-wavelength radar measurements, Geophysical Research Letters, 42, 176–183, 2015. Stephens, G. L., Vane, D. G., Boain, R. J., Mace, G. G., Sassen, K., Wang, Z., Illingworth, A. J., O'Connor, E. J., Rossow, W. B., Durden, S. L., Miller, S. D., Austin, R. T., Benedetti, A., and Mitrescu, C.: The CloudSat mission and the A-Train: A new dimension of space-based observations of clouds and precipitation, Bulletin of the American Meteorological Society, 83, 1771–1790, 2002.Tiira, J., Moisseev, D. N., von Lerber, A., Ori, D., Tokay, A., Bliven, L. F., and Petersen, W.: Ensemble mean density and its connection to other microphysical properties of falling snow as observed in Southern Finland, Atmospheric Measurement Techniques, 9, 4825–4841, 10.5194/amt-9-4825-2016, 2016.Tridon, F., Battaglia, A., Chase, R., Turk, F. J., Leinonen, J., Kneifel, S., Mroz, K., Finlon, J., Bansemer, A., Tanelli, S., Heymsfield, A., and Nesbitt, S.: The microphysics of stratiform precipitation during OLYMPEX: Compatibility between triple-frequency radar and airborne in situ observations, Journal of Geophysical Research: Atmospheres, 124, 8764–8792, 2019. Tridon, F., Battaglia, A., and Kneifel, S.: Estimating total attenuation using Rayleigh targets at cloud top: applications in multilayer and mixed-phase clouds observed by ground-based multifrequency radars, Atmospheric Measurement Techniques, 13, 5065–5085, 10.5194/amt-13-5065-2020, 2020.Tridon, F., Silber, I., Battaglia, A., Kneifel, S., Fridlind, A., Kalogeras, P., and Dhillon, R.: Highly supercooled riming and unusual triple-frequency radar signatures over McMurdo Station, Antarctica, Atmospheric Chemistry and Physics, 22, 12467–12491, 10.5194/acp-22-12467-2022, 2022.Turk, F. J., Ringerud, S. E., Camplani, A., Casella, D., Chase, R. J., Ebtehaj, A., Gong, J., Kulie, M., Liu, G., Milani, L., Panegrossi, G., Padullés, R., Rysman, J.-F., Sanò, P., Vahedizade, S., and Wood, N. B.: Applications of a CloudSat-TRMM and CloudSat-GPM satellite coincidence dataset, Remote Sensing, 13, 2264, 10.3390/rs13122264, 2021. 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, Journal of Geophysical Research: Atmospheres, 118, 13–199, 2013.von Terzi, L.: TRIPEx-pol level-2 data, Zenodo [data set], 10.5281/zenodo.5025636, 2021.von Terzi, L., Dias Neto, J., Ori, D., Myagkov, A., and Kneifel, S.: Ice microphysical processes in the dendritic growth layer: a statistical analysis combining multi-frequency and polarimetric Doppler cloud radar observations, Atmospheric Chemistry and Physics, 22, 11795–11821, 10.5194/acp-22-11795-2022, 2022.Widener, K. and Johnson, K.: W-band arm cloud radar (wacr) handbook, The US Department of Energy, Washington DC, USA, 10.2172/1019541, 2006. Yu, W., Tao, H., Jingchao, G., Kai, J., Rui, L., Wencheng, S., and Guosheng, L.: Simulation of the capability of Ku, Ka and W tri-frequency satellite-borne radar measuring the three-dimensional structure of cloud and precipitation, Chinese Science Bulletin, 64, 430–443, 2019.Zhang, P., Gu, S., Chen, L., Shang, J., Lin, M., Zhu, A., Yin, H., Wu, Q., Shou, Y., Sun, F., Xu, H., Yang, G., Wang, H., Li, L., Zhang, H. W., Chen, S., and Lu, N.: FY-3G satellite instruments and precipitation products: first report of China’s Fengyun rainfall mission in-orbit, Journal of Remote Sensing, 3, 0097, 10.34133/remotesensing.0097, 2023.Zhang, Y., Zheng, H., Zhang, L., Huang, Y., Liu, X., and Wu, Z.: Assessing the effect of riming on snow microphysics: The first observational study in East China, Journal of Geophysical Research: Atmospheres, 126, e2020JD033763, 10.1029/2020JD033763, 2021.Zhang, Y., Ouyang, X., Wang, M., Rosenfeld, D., Zhao, D., and Wu, X.: Evaluating ice phase microphysics in the simulation of a snowstorm over Northern China, Journal of Geophysical Research: Atmospheres, 129, e2023JD040221, 10.1029/2023JD040221, 2024.