Next Article in Journal
Hydrologic and Pollutant Removal Performance of Media Layers in Bioretention
Previous Article in Journal
Seepage Characteristics of a Single Ascending Relief Well Dewatering an Overlying Aquifer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microphysical Characteristics of Winter Precipitation in Eastern China from 2014 to 2019

1
College of Meteorology and Oceanography, National University of Defense Technology, Nanjing 211101, China
2
Beijing Aerospace Control Center, Beijing 100094, China
3
Nanjing Meteorological Bureau, China Meteorological Administration, Nanjing 210019, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(3), 920; https://doi.org/10.3390/w12030920
Submission received: 18 February 2020 / Revised: 17 March 2020 / Accepted: 20 March 2020 / Published: 24 March 2020
(This article belongs to the Section Hydrology)

Abstract

:
To improve solid precipitation monitoring in the hydrology and meteorology field, 1-min precipitation data observed by the PARticle SIze VELocity (PARSIVEL) disdrometer in Nanjing, eastern China, from February 2014 to February 2019 for all days with solid precipitation, were used to study the microphysical characteristics of winter precipitation. In this study, the empirical V-D (velocity–diameter) relationships and observed surface temperature are used for matching precipitation types, and the precipitation data are divided into rain, graupel, wet snow and dry snow. The results show that dry snow and wet snow have maximum Dm (mass-weighted mean diameter) and minimum log10Nw (normalized intercept parameter), while rain shows the opposite. Additionally, the μ-Λ (shape parameter–slope parameter) curve of dry snow and wet snow is very close, and the μ value of dry snow and wet snow is higher than that of graupel and higher than that of rain for the same Λ value. Furthermore, the Ze-S (equivalent reflectivity factor–precipitation intensity) relationships among different types of precipitation are significantly different. If only the Ze-S relationship of rain is used for quantitative precipitation estimation (QPE), then, for small precipitation intensity, solid precipitation will be overestimated, while, for large precipitation intensity, it will be underestimated.

1. Introduction

The accurate monitoring of solid precipitation is of great significance for aviation safety, transportation and freezing disaster prevention, especially in the middle and high latitudes [1,2,3]. Additionally, the study of the microphysical characteristics of solid precipitation contributes to the development of weather forecasts, remote sensing of precipitation, and hydrology [1,4,5,6]. The mass-size and velocity-size relationships of solid precipitation particles determine the correctness of Numerical Weather Prediction (NWP) model parameterization, which has an impact on the prediction of precipitation location, precipitation intensity and start and end times [7,8,9]. Whether passive or active, ground-based or space-based electromagnetic propagation models are sensitive to the shape and density of solid precipitation particles, and a reasonable description of precipitation microphysical characteristics is essential to the inversion algorithm of precipitation intensity based on remote sensing [10,11,12]. The erosion effect of precipitation on soil is closely related to its type and the particle size distribution (PSD), which is very important for the flood control and water fixation abilities of ecosystems [13,14,15].
Previous studies have been conducted on the microphysical characteristics of solid precipitation. Barthazy et al. [16] calibrated the Hydrometeor Velocity and Shape Detector (HVSD) instrument and measured the fall speed, shape and PSD of snowflakes. Zhang et al. [17] observed several weather processes involving solid precipitation using a 2-dimensional video disdrometer (2DVD) and compared the theoretically calculated polarization variables with the radar-measured values to verify the effectiveness of the density correction algorithm in Oklahoma. Jia et al. [10] used a PARticle SIze VELocity (PARSIVEL) disdrometer and microscopic photography to analyze the fall velocity and PSD of rain, graupel, snow and mixed-phase precipitation particles in northern China.
With the development of measurement technology in recent years, observational experiments of solid precipitation microphysical characteristics have been carried out routinely. However, on the one hand, due to the great difference in the physical properties of different solid precipitation particles, it is still an urgent problem to invert solid precipitation by radar remote sensing [18,19,20], but on the other hand, the data collected by previous researchers were mostly limited by quantity; their statistical significance has yet be further verified, and experiments were rarely carried out in many regions of the world, such as East Asia. By distinguishing precipitation types in advance, the accuracy of the quantitative precipitation estimation (QPE) of the polarization radar can be improved [3,21]. By establishing a more accurate electromagnetic propagation model according to the specific solid precipitation type, the uncertainty of the inversion algorithm can be reduced [22]. In addition, due to differences in regional climate and other factors, conducting more local statistical experiments is also a powerful means of improving the remote sensing measurement of solid precipitation information.
To improve the application level of remote sensing for solid precipitation in eastern China, this study analyzes the microphysical characteristics of rain, graupel, dry snow and wet snow by using 6 years of winter PARSIVEL (first generation) observation data in Nanjing. The next section describes the data sources, parameter calculations, and precipitation classification methods. Section 3 presents the analysis of the differences among the microphysical characteristics of different types of precipitation. The research methods are discussed in Section 4. The final section includes the conclusion and summary.

2. Data and Method

2.1. Data Sources

In this study, winter precipitation data from February 2014 to February 2019, observed by an OTT PARSIVEL disdrometer [10,23] at a station (31.97° N, 118.81° E, 32 m above sea level) in Nanjing, eastern China, were used (as shown in Figure 1). Nanjing is located in the north subtropical monsoon climate zone, which has abundant precipitation every year from approximately mid-November into winter. The temperature difference between winter and summer is significant, and the annual minimum temperature can reach −13.1 °C. With complex weather systems and frequent natural disasters, the weather and climate characteristics of Nanjing are representative of eastern China [24]. The selected precipitation days are the days with non-liquid precipitation verified by PARSIVEL data. As shown in Table 1, a total of 29 precipitation days were identified. On these precipitation days, all 1-min precipitation data (22,386 samples) recorded by PARSIVEL were used for the statistical analysis in this study. In addition to disdrometer data, local 1-h ground temperature data for precipitation days are also used to distinguish between dry snow and wet snow (see Section 2.4).

2.2. Quality Control

The PARSIVEL disdrometer determines the size of the precipitation particles by measuring the attenuation caused by precipitation particles across its laser beam and determines the fall speed of precipitation particles by recording the time when the laser beam is blocked [25]. The output data of the instrument is a matrix of 32 by 32: 32 bins with sizes ranging from 0.062 to 24.5 mm and 32 bins with speeds ranging from 0.05 to 20.5 m/s [26,27]. Because of its good cost performance, the instrument is widely used in the QPE of radar, soil erosion and weather modification [13,28,29]. However, due to the limitations of the device performance and built-in algorithm, quality control is required to obtain the best possible measurement results. Because of the low signal-to-noise ratio, data from the first two size bins are not used [30]. For every 1-min sample, if the number of precipitation particles is less than 10, the sample is excluded [27]. In addition, the sampling area is corrected as (30 mm − 0.5Di) × 180 mm to remove the boundary effect [31]. The velocity outliers on each size bin are eliminated according to the 3σ criterion (small particles with high velocity caused by spatter, large particles with low velocity caused by overlap). Through quality control, a total of 21274 samples were obtained.

2.3. Calculation of Integral Parameters

The number concentration N(Di) (m−3 mm−1), which refers to the number of precipitation particles per unit volume per unit diameter, can be expressed as follows
N ( D i ) = j = 3 32 C i j V j A t Δ D i
where Di (mm) is the median of the ith size bin, t (s) is the sampling time (t = 60 s in this study), A (mm2) is the sampling cross section (A = 54 mm2 for OTT PARSIVEL disdrometer), Vj (m/s) is the median of the jth velocity bin, ΔDi (mm) is the width of ith size bin, and Cij is the amount in ith size bin and jth velocity bin during sixty seconds. To better describe PSD, we need to fit it with appropriate functions. The commonly used function is the three-parameter gamma function due to its universality [10,32]
N ( D ) = N 0 D μ exp ( Λ D )
where D (mm) is the precipitation particle diameter, N0 (mm−1-μ m−3) is the intercept parameter, μ is the shape parameter, and Λ (mm−1) is the slope parameter. These three parameters can be determined by the order moment method [33,34]. The nth moment of PSD can be defined as follows
M n = i = 3 32 D i n N ( D i ) Δ D i
The intercept parameter N0, shape parameter μ and slope parameter Λ can be expressed as follows
N 0 = Λ ( μ + 4 ) M 3 Γ ( μ + 4 )
μ = 11 G 8 + G ( G + 8 ) 2 ( 1 G )
Λ = ( μ + 4 ) M 3 M 4
where G = M 4 3 M 3 2 M 6 .
In addition, the normalized intercept parameter Nw and the mass-weighted mean diameter Dm are defined as follows [35]
N w = 4 4 Γ ( 4 ) M 3 5 M 4 4
D m = M 4 M 3
The precipitation intensity S (mm h−1) (liquid equivalent) can be calculated by the following formula
S = 6 × 10 - 1 ρ w A i = 3 32 M ( D i ) N ( D i ) Δ D i
where water density ρw = 1 g cm3, and precipitation particle mass M(Di) (mg) employs an empirical relationship from Locatelli and Hobbs [36] for solid precipitation (as shown in Table 2).
The equivalent reflectivity factor of ice particle Ze can be expressed as follows
Z e = | K i | 2 | K w | 2 i = 3 32 D i 6 N ( D i ) Δ D i
where D i (mm) is the equivalent ice sphere diameter and Ki and Kw are the dielectric constants for ice and water, respectively (|Ki|2 = 0.208 and |Kw|2 = 0.93 for the usual meteorological radar wavelength) [37]. Assuming Di corresponds to D i , then we can obtain the following expression
M ( D i ) = M ( D i ) = π 6 ρ i D i 3
where the ice density ρi = 0.9167 g/cm3. Equation (11) is substituted into Equation (10) to obtain the following Equation [1]
Z e = | K i | 2 | K w | 2 ( 6 π ρ i ) 2 i = 1 32 M ( D i ) 2 N ( D i ) Δ D i

2.4. Classification of Precipitation Types

Common methods for distinguishing precipitation types include using Formvar slides to collect particles and manually determine the types with a microscope [10,39,40,41] or automatically judging the precipitation particle pictures taken by combining with an image recognition algorithm [7,12]. In addition, because the VD relationships between different types of precipitation are obviously different, the VD function fitted by 1-min average data can also be used to distinguish precipitation types [11,42,43]. In this study, the Vfit = aDb relationship fitted with a 1-min sample (using the least square method) was matched with the 15 V-D relationships in Table 2, and the type with the smallest difference was considered to be the type of the 1-min sample
arg min x i = 3 32 j = 1 32 C i j ( V f i t ( D i ) V x ( D i ) ) 2
where x refers to the 15 precipitation types in Table 2. Accumulated raw particle counts by size and velocity for every precipitation type are shown in Figure 2. Through this classification method, the numbers of 1-min samples of rain, graupel and snow are 11,196, 7521 and 2557, respectively.
Furthermore, depending on the temperature near the ground, snowfall may be divided into dry snow and wet snow, whose microphysical characteristics may vary greatly [2]. According to the ground temperature data recorded on 1 h time resolution, the snow samples obtained above are further divided into dry snow and wet snow. If the temperature t ≥ 0 °C, the sample is considered to be wet snow; otherwise, it is considered to be dry snow. Finally, 2014 samples of dry snow and 545 samples of wet snow are obtained.

3. Results

3.1. The PSD of Different Types of Precipitation

Figure 3 shows the distribution of the mean particle concentrations of rain, graupel, wet snow and dry snow. The particles of the three types of solid precipitation (graupel, wet snow and dry snow) are much larger overall than those of rainfall. Raindrops larger than 6 mm in diameter have a concentration of less than 0.01 m−3 mm−1, but particles of three types of solid precipitation can exceed 13 mm in diameter at the same concentration. Among them, the size of wet snow particles is smaller than that of dry snow for the same concentration, which may be because when the temperature is above freezing, the volume of hydrometeors decreases due to the transformation from the solid to liquid phase. In addition, the concentration of wet snow is lower than that of graupel and dry snow at all diameters. The concentration of graupel is higher than that of dry snow for small hydrometeors (D < 1.625 mm) and is the opposite for large hydrometeors (D > 3.25 mm). Furthermore, wet snow and dry snow have the largest peak diameters (0.937 and 0.687 mm, respectively), followed by graupel (0.562 mm), then rain (0.437 mm), which is consistent with the conclusion of Jia et al. [10]. As shown in Table 3, the skewness and kurtosis of the PSD of the four types of precipitation are all positive. The PSD has a maximum skewness of 2.32 and a maximum kurtosis of 4.12 for rain (meaning that the proportion of hydrometeors greater than the peak diameter is the highest and hydrometeors are most concentrated at the peak diameter for rainfall), whereas the opposite is true for wet snow.

3.2. The log10Nw-Dm Distributions of Different Types of Precipitation

The normalized intercept parameter log10Nw and the mass-weighted mean diameter Dm directly reflect the concentration and size characteristics of precipitation particles, and their distribution is of great significance to the classification of precipitation types [27,28,33,44]. To study the differences in the microphysical characteristics of different types of precipitation in Nanjing, the statistical results of log10Nw versus Dm are shown in Figure 4. As seen from the figure, rain generally has the highest concentration and minimal size, while wet snow and dry snow have the lowest concentrations and maximum size. In addition, dry snow and wet snow have similar concentrations, but the particle size of dry snow is larger than that of wet snow (probably due to melting, as mentioned above). Furthermore, the Dm distribution of the three types of solid precipitation is much more widespread than that of rain, especially for dry snow.
The average log10Nw-Dm value of the rain is closer to the stratiform rain line (green dashed line) measured by Bringi et al. [28], which suggests that the winter rainfall in the Nanjing area mostly comes from the stratiform cloud system. The average log10Nw-Dm value of the three types of solid precipitation is closer to the continental convective rain area (black dashed rectangle), especially for graupel. The mean and standard deviation of specific log10Nw and Dm are shown in Table 4. Unfortunately, due to the lack of previous studies on the microphysical characteristics of solid precipitation, the log10Nw-Dm of graupel, wet snow and dry snow in this study cannot be compared with the statistical results in other regions.

3.3. The μ–Λ Relationships among Different Types of Precipitation

The three-parameter gamma function described in Equation (2) is widely used in fields such as ground-based or space-based remote sensing, numerical weather prediction and microwave communication [1,45,46,47]. However, since there are as many as three unknown parameters and usually only two data sources at most (such as dual-frequency radar on GPM [48] and ground-based radar [49]), their inversion often requires certain constraints on the parameters. It is common to assume that the shape parameter μ is a constant [50]. However, it has been proven that μ is highly variable, and it is heavily influenced by the microphysical characteristics of precipitation [27,51]. Thus, a localized μ–Λ relationship is proposed to reduce the number of required inversion parameters [29,44,52,53].
Figure 5 shows the μ–Λ scatters of the four types of precipitation and their fitting curves obtained by the least square method. As seen from the figure, for the same Λ, the μ value of dry snow and wet snow is greater than that of graupel and greater than that of rain. In addition, the μ–Λ curves of dry snow and wet snow in the interval of 5 < Λ < 25 are close to coincidental. Furthermore, excluding Tang et al. [53], the μ–Λ relationship of rain from previous studies is very close to the μ–Λ relationship obtained in this study in the small Λ value region (Λ < 15), especially that obtained by Wen et al. [44], whose dataset also comes from the winter rain in Nanjing between 2014 and 2016. The parameter values of the μ–Λ relationship of each type of precipitation in this study and previous studies are shown in Table 5. Additionally, due to the lack of research on the μ–Λ relationship of solid precipitation, the μ–Λ relationship of graupel, dry snow and wet snow obtained in this study is difficult to compare with other studies.

3.4. The Ze–S Relationships among Different Types of Precipitation

The Z–S relationship is mainly used for the QPE of weather radar [27,55]. However, instead of the reflectivity factor Z, weather radar actually measures the equivalent reflectivity factor Ze. For rain, the equivalent reflectivity factor Ze is equal to the reflectivity factor Z, because the raindrops are small and nearly spherical [37]. Nevertheless, for snow particles or graupel particles, due to their large size and complex shape, the above approximation is obviously not applicable [37]. Through calculation, it can be obtained that Ze = 0.189Z if the size of the hydrometeors is expressed as the equivalent ice sphere diameter [2].
Figure 6 shows the Ze–S relationship fitted in this study and in previous studies. It should be mentioned that most of the previous results showed a Z–S relationship. For the convenience of comparison, Ze = 0.189Z is used to convert this relationship into a Ze–S relationship. Table 6 records these specific parameters. As illustrated in this figure, when S > 1.1 mm h−1, for the same Ze, the precipitation intensity of dry snow is greater than that of graupel and wet snow, and greater than that of rain. The Ze–S relationships between graupel and wet snow are the closest, especially for the large precipitation intensity range (S > 10 mm h−1). For small precipitation intensities (S < 0.2 mm h−1), the equivalent reflectivity factor of dry snow is greater than that of wet snow, followed by graupel, and then rain for the same S. This means that if the Ze–S relationship of rain is uniformly used for the inversion of S, for heavy precipitation, the S of solid precipitation will be underestimated, while, for weak precipitation, the S of solid precipitation will be overestimated.
The Ze–S relationships of wet snow from Imal [56] and Puhakka [57] are very consistent with those of the wet snow fitted in this study. In addition, the Ze–S relationship of snow from Gunn and Marshall [58] is also very similar to that of wet snow, while the Ze–S relationship of snow from Sekon and Srivastava [18] seems to be more similar to that of rain in this study. Furthermore, the Ze–S relationships of dry snow from Imai [56] and Puhakka [57] differ significantly from those of this study, especially for heavy precipitation intensity. The same is true for the Ze–S relationships between graupel from Gray and Male [59] and graupel in this study.

4. Discussion

The principle of the precipitation classification method used in this study is that there are differences in the V–D relationships among different types of precipitation particles, which is also the method used by the PARSIVEL disdrometer to distinguish the precipitation types. Similarly, Yuter et al. [11] also studied the microphysical characteristics of coexisting rain and wet snow by distinguishing the V–D relationship differences. In addition, Jia et al. [10] analyzed the V–D relationships of raindrops, graupel particles and snowflakes and found that these particles largely depended on the fitted empirical function, indicating the reliability of using the V–D relationship to identify precipitation types. However, the fall velocity of precipitation particles fluctuates under the influence of vertical airflow. According to the differences in specific weather conditions, the empirical relationship may produce errors. In this regard, Niu et al. [32] used local air density to revise the actual fall velocity of precipitation particles to exclude the influence of airflow. Furthermore, fewer field experiments have resulted in fewer samples being available at present. For example, Locatelli and Hobbs [36] used only a few dozen datapoints to fit the V–D relationships, but they have been used in the parameterization of weather and climate models [7]. Therefore, it may be inappropriate to use the empirical relationships obtained by Locatelli and Hobbs to match the precipitation type.
The estimation method of PSD parameters also needs to be discussed. In this study, the estimation method of PSD parameters is the moment method without truncation. However, in practice, the estimation effects of truncation or no truncation, moment method, least square method or maximum likelihood method are different. Mallet and Barthes [60] carried out the simulation of different estimation methods considering noise measurements, the absence of small drops, size of the collecting area, and integration time. The results show that the estimation effect of the maximum likelihood method is better than the others. Hence, the estimation method of PSD parameters can be further improved. In addition, different time steps of classification also affect the observation results, which are 1-min in this study. Matching the type based on the V–D distribution obtained by too short a time step (meaning too few particles recorded in the time interval) will result in variability in the result, while matching the type based on the V–D distribution obtained by too long a time step (meaning too many particles recorded in the time interval) may cause the result to be wrong due to the smoothness.
The PSD data used in this study were measured by a PARSIVEL disdrometer. This instrument is widely used in the study of worldwide precipitation microphysical characteristics because of its high-cost performance and portability [29]. It was originally primarily developed as a PSD measurement of rainfall and can be used as a weather sensor to distinguish precipitation types [27]. Löffler-Mang and Joss [42] investigated the potential of a PARSIVEL disdrometer measuring the PSD of solid precipitation. By adjusting the mass–size relationship, the reflectivity factor of snow measured by PARSIVEL is in good agreement with the measured value of the C-band radar [61]. However, the greatest obstacle to using it to measure solid precipitation is the built-in algorithm based on the raindrop hypothesis. Since snowflakes tend to be asymmetrical, the PARSIVEL disdrometer may underestimate their size [10]. The velocity is also subject to error, as it is determined by the vertical size of the snowflake. In addition, the PARSIVEL disdrometer can only count in discrete velocities and size ranges, and the velocity and diameter values can only be represented by the median of each bin, resulting in a quantitative error. Moreover, due to the stickiness of snow particles, some abnormally large particles recorded also introduce errors, especially for wet snow.
There are other instruments that can be used to observe precipitation microphysical characteristics, but they also have their own shortcomings. For example, 2DVD can obtain the two-dimensional contour of precipitation particles, but it will underestimate the number of precipitation particles with diameters of less than 0.3 mm [29]; snowflake video imager (SVI) can record the image of snowflakes, but not their fall speeds [22]. To obtain more accurate measurement results, researchers are committed to developing new measurement methods and instruments, such as multiangle snow imager (MASC) [62], and precipitation microphysical characteristics sensor (PMCS) [63]. These advances will have an important impact on QPE, NWP and other fields.

5. Conclusions

In this study, the microphysical characteristics of winter precipitation in eastern China were analyzed by using the PARSIVEL disdrometer in Nanjing. The four precipitation types were classified and matched based on the V–D relationship, and their PSD, μ-Λ, log10Nw-Dm and Ze-S values were also compared.
First, the spectrum of the three types of solid precipitation (graupel, wet snow and dry snow) is much wider than that of rain. For large particles, the concentration of dry snow is greater than that of graupel and greater than that of wet snow. For small particles, the concentration of graupel is greater than that of dry snow and greater than that of wet snow. Second, the size of the dry snow particles is larger than that of the wet snow particles, but the overall concentration is the same. Graupel has a smaller size and a larger number of particles, followed by rain. Moreover, the μ–Λ curves for dry snow and wet snow are very similar in the small Λ region. In general, for the same Λ value, the μ value of wet snow is greater than that of dry snow, followed by that of graupel, and finally rain. Furthermore, there are significant differences in the Ze–S relationships of different types of precipitation. If only the Ze–S relationship of rain is used for QPE, the intensity of heavy and light solid precipitation will be underestimated and overestimated, respectively.
With the invention of instruments using new principles and the continuous improvement of the existing instruments’ performances, the observation of the microphysical properties of precipitation has become increasingly accurate. In addition to PSD, researchers have even begun to study the three-dimensional shape, riming degree and density of precipitation particles. With a more accurate description of the microscopic properties of precipitation and a deeper understanding of the precipitation evolution process, existing parameterization schemes and radiation transmission models will be further improved, especially for polarization radar.

Author Contributions

Conceptualization, K.P. and X.L.; Data curation, X.L.; Funding acquisition, X.L. and Y.W.; Investigation, H.H. and Y.W.; Methodology, K.P.; Resources, X.L.; Supervision, X.L., H.H. and S.H.; Validation, K.P., X.L., Y.S. and S.H.; Visualization, K.P. and Y.S.; Writing—original draft, K.P.; Writing—review and editing, K.P. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 41975030, 41327003, and 41475020) and Scientific Research Projects of Nanjing Meteorological Bureau, China (grant number NJ201907).

Acknowledgments

The authors would like to express their gratitude to the three anonymous reviewers for their comments that improved this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boudala, F.S.; Isaac, G.A. Parameterization of visibility in snow: Application in numerical weather prediction models. J. Geophys. Res. 2009, 114. [Google Scholar] [CrossRef] [Green Version]
  2. Rasmussen, R.; Dixon, M.; Vasiloff, S.; Hage, F.; Knight, S.; Vivekanandan, J.; Xu, M. Snow Nowcasting Using a Real-Time Correlation of Radar Reflectivity with Snow Gauge Accumulation. J. Appl. Meteorol. 2003, 42, 20–36. [Google Scholar] [CrossRef] [Green Version]
  3. Notaroš, B.; Bringi, V.; Kleinkort, C.; Kennedy, P.; Huang, G.-J.; Thurai, M.; Newman, A.; Bang, W.; Lee, G. Accurate Characterization of Winter Precipitation Using Multi-Angle Snowflake Camera, Visual Hull, Advanced Scattering Methods and Polarimetric Radar. Atmosphere 2016, 7, 81. [Google Scholar] [CrossRef] [Green Version]
  4. Hongyan, L.; Hongbin, C.; Hengchi, L.; Yuxia, W. Relationship between rain rate and radar reflectivity based on the raindrop distribution data in Beijing during 2004. Acta Meteorol. Sin. 2008, 66, 125–129. [Google Scholar] [CrossRef]
  5. Holt, A.R.; Cummings, R.J.; Upton, G.J.G.; Bradford, W.J. Rain rates, drop size information, and precipitation type, obtained from one-way differential propagation phase and attenuation along a microwave link. Radio Sci. 2008, 43, 1–18. [Google Scholar] [CrossRef]
  6. Angulo-Martínez, M.; Beguería, S.; Latorre, B.; Fernández-Raga, M. Comparison of precipitation measurements by OTT Parsivel2 and Thies LPM optical disdrometers. Hydrol. Earth Syst. Sci. 2018, 22, 2811–2837. [Google Scholar] [CrossRef] [Green Version]
  7. Garrett, T.J.; Yuter, S.E. Observed influence of riming, temperature, and turbulence on the fallspeed of solid precipitation. Geophys. Res. Lett. 2014, 41, 6515–6522. [Google Scholar] [CrossRef]
  8. Praz, C.; Roulet, Y.-A.; Berne, A. Solid hydrometeor classification and riming degree estimation from pictures collected with a Multi-Angle Snowflake Camera. Atmos. Meas. Tech. 2017, 10, 1335–1357. [Google Scholar] [CrossRef] [Green Version]
  9. Woods, C.P.; Stoelinga, M.T.; Locatelli, J.D. The IMPROVE-1 storm of 1-2 February 2001. Part III: Sensitivity of a mesoscale model simulation to the representation of snow particle types and testing of a bulk microphysical scheme with snow habit prediction. J. Atmos. Sci. 2007, 64, 3927–3948. [Google Scholar] [CrossRef]
  10. Jia, X.; Liu, Y.; Ding, D.; Ma, X.; Chen, Y.; Bi, K.; Tian, P.; Lu, C.; Quan, J. Combining disdrometer, microscopic photography, and cloud radar to study distributions of hydrometeor types, size and fall velocity. Atmos. Res. 2019, 228, 176–185. [Google Scholar] [CrossRef]
  11. Yuter, S.E.; Kingsmill, D.E.; Nance, L.B.; Martin, L.-M. Observations of Precipitation Size and Fall Speed Characteristics within Coexisting Rain and Wet Snow. J. Appl. Meteorol. Climatol. 2006, 45, 1450–1464. [Google Scholar] [CrossRef] [Green Version]
  12. Garrett, T.J.; Yuter, S.E.; Fallgatter, C.; Shkurko, K.; Rhodes, S.R.; Endries, J.L. Orientations and aspect ratios of falling snow. Geophys. Res. Lett. 2015, 42, 4617–4622. [Google Scholar] [CrossRef]
  13. Angulo-Martínez, M.; López-Vicente, M.; Vicente Serrano, S.M.; Beguería, S. Mapping rainfall erosivity at a regional scale: A comparison of interpolation methods in the ebro basin (ne spain). Hydrol. Earth Syst. Sci. 2009, 13, 1907–1920. [Google Scholar] [CrossRef] [Green Version]
  14. Nanko, K.; Moskalski, S.M.; Torres, R. Rainfall erosivity–intensity relationships for normal rainfall events and a tropical cyclone on the US southeast coast. J. Hydrol. 2016, 534, 440–450. [Google Scholar] [CrossRef] [Green Version]
  15. Vaezi, A.R.; Abbasi, M.; Bussi, G.; Keesstra, S. Modeling sediment yield in semi-arid pasture micro-catchments, NW Iran. Land Degrad. Dev. 2017, 28, 1274–1286. [Google Scholar] [CrossRef]
  16. Barthazy, E.; Göke, S.; Schefold, R.; Högi, D. An optical array instrument for shape and fall velocity measurements of hydrometeors. J. Atmos. Ocean. Technol. 2004, 21, 1400–1416. [Google Scholar] [CrossRef]
  17. Zhang, G.; Luchs, S.; Ryzhkov, A.; Xue, M.; Ryzhkova, L.; Cao, Q. Winter Precipitation Microphysics Characterized by Polarimetric Radar and Video Disdrometer Observations in Central Oklahoma. J. Appl. Meteorol. Climatol. 2011, 50, 1558–1570. [Google Scholar] [CrossRef] [Green Version]
  18. Sekhon, R.S.; Srivastava, R.C. Snow size spectra and radar reflectivity. J. Atmos. Sci. 1970, 27, 299–307. [Google Scholar] [CrossRef] [Green Version]
  19. Fujiyoshi, Y.; Endoh, T.; Yamada, T.; Tsuboki, K.; Tachibana, Y.; Wakahama, G. Determination of a Z-R relationship for snowfall using a radar and high sensitivity snow gauges. J. Appl. Meteorol. 1990, 29, 147–152. [Google Scholar] [CrossRef] [Green Version]
  20. Kubo, M.; Sougen, A.; Muramoto, K.; Fujiyoshi, Y. ZR relation for snowfall using small Doppler radar and high sensitive snow gauges. In Proceedings of the ICCAS-SICE, Fukuoka, Japan, 18–21 August 2009. [Google Scholar]
  21. Huang, G.-J.; Bringi, V.N.; Cifelli, R.; Hudak, D.; Petersen, W.A. A methodology to derive radar reflectivity–liquid equivalent snow rate relations using C-band radar and a 2D video disdrometer. J. Atmos. Ocean. Technol. 2010, 27, 637–651. [Google Scholar] [CrossRef]
  22. Newman, A.J.; Kucera, P.A.; Bliven, L.F. Presenting the Snowflake Video Imager (SVI). J. Atmos. Ocean. Technol. 2009, 26, 167–179. [Google Scholar] [CrossRef]
  23. Battaglia, A.; Rustemeier, E.; Tokay, A.; Blahak, U.; Simmer, C. PARSIVEL Snow Observations: A Critical Assessment. J. Atmos. Ocean. Technol. 2010, 27, 333–344. [Google Scholar] [CrossRef]
  24. Yin, D.; Sun, Y. Climatic Characteristics of Low-temperature Freezing and Snow Calamities in Nanjing in Recent 50 Years. J. Catastrophol. 2011, 26. [Google Scholar] [CrossRef]
  25. Ghada, W.; Buras, A.; Lüpke, M.; Schunk, C.; Menzel, A. Rain Microstructure Parameters Vary with Large-Scale Weather Conditions in Lausanne, Switzerland. Remote Sens. 2018, 10, 811. [Google Scholar] [CrossRef] [Green Version]
  26. Wu, Y.; Liu, L. Statistical characteristics of raindrop size distribution in the Tibetan Plateau and southern China. Adv. Atmos. Sci. 2017, 34, 727–736. [Google Scholar] [CrossRef]
  27. Chen, B.; Yang, J.; Pu, J. Statistical Characteristics of Raindrop Size Distribution in the Meiyu Season Observed in Eastern China. J. Meteorol. Soc. Jpn. Ser. II 2013, 91, 215–227. [Google Scholar] [CrossRef] [Green Version]
  28. Bringi, V.N.; Chandrasekar, V.; Hubbert, J. Raindrop Size Distribution in Different Climatic Regimes from Disdrometer and Dual-Polarized Radar Analysis. J. Atmos. Sci. 2003, 60, 354–365. [Google Scholar] [CrossRef]
  29. Zhang, A.; Hu, J.; Chen, S.; Hu, D.; Liang, Z.; Huang, C.; Xiao, L.; Min, C.; Li, H. Statistical Characteristics of Raindrop Size Distribution in the Monsoon Season Observed in Southern China. Remote Sens. 2019, 11, 432. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, B.; Hu, Z.; Liu, L.; Zhang, G. Raindrop Size Distribution Measurements at 4500 m on the Tibetan Plateau During TIPEX-III. J. Geophys. Res. Atmos. 2017, 122. [Google Scholar] [CrossRef]
  31. Jaffrain, J.; Berne, A. Experimental quantification of the sampling uncertainty associated with measurements from PARSIVEL disdrometers. J. Hydrometeorol. 2011, 12, 352–370. [Google Scholar] [CrossRef]
  32. Niu, S.; Jia, X.; Sang, J.; Liu, X.; Lu, C.; Liu, Y. Distributions of Raindrop Sizes and Fall Velocities in a Semiarid Plateau Climate: Convective versus Stratiform Rains. J. Appl. Meteorol. Climatol. 2010, 49, 632–645. [Google Scholar] [CrossRef]
  33. Seela, B.K.; Janapati, J.; Lin, P.-L.; Wang, P.K.; Lee, M.-T. Raindrop Size Distribution Characteristics of Summer and Winter Season Rainfall over North Taiwan. J. Geophys. Res. Atmos. 2018, 123. [Google Scholar] [CrossRef] [Green Version]
  34. Balaji Kumar, S.; Jayalakshmi, J.; Pay-Liam, L.; Krishna, R.K.; Ryuichi, S.; Pao, K.W. A Comparison Study of Summer Season Raindrop Size Distribution between Palau and Taiwan, Two Islands in Western Pacific. J. Geophys. Res. Atmos. 2017, 122. [Google Scholar] [CrossRef]
  35. Wang, M.; Zhao, K.; Xue, M.; Zhang, G.; Liu, S.; Wen, L.; Chen, G. Precipitation microphysics characteristics of a Typhoon Matmo (2014) rainband after landfall over eastern China based on polarimetric radar observations. J. Geophys. Res. Atmos. 2016, 121, 12–415. [Google Scholar] [CrossRef]
  36. Locatelli, J.D.; Hobbs, P.V. Fall speeds and masses of solid precipitation particles. J. Geophys. Res. 1974, 79, 2185–2197. [Google Scholar] [CrossRef]
  37. Smith, P.L. Equivalent radar reflectivity factors for snow and ice particles. J. Clim. Appl. Meteorol. 1984, 23, 1258–1260. [Google Scholar] [CrossRef] [Green Version]
  38. Atlas, D.; Srivastava, R.C.; Sekhon, R.S. Doppler characteristics of precipitation at vertical incidence. Rev. Geophys. Space Phys. 1973, 11, 1–35. [Google Scholar] [CrossRef]
  39. Barthazy, E.; Schefold, R. Fall velocity of snowflakes of different riming degree and crystal types. Atmos. Res. 2006, 82, 391–398. [Google Scholar] [CrossRef]
  40. Schefold, R.; Baschek, B.; Wüest, M.; Barthazy, E. Fall velocity and axial ratio of snowflakes. In Proceedings of the Second European Conference on Radar Meteorology (ERAD) in conjunction with COST 717 mid-term seminar, Netherlands, 18–22 November 2002; pp. 84–89. [Google Scholar]
  41. Schmitt, C.G.; Heymsfield, A.J. The Dimensional Characteristics of Ice Crystal Aggregates from Fractal Geometry. J. Atmos. Sci. 2010, 67, 1605–1616. [Google Scholar] [CrossRef]
  42. Löffler-Mang, M.; Joss, J. An Optical Disdrometer for Measuring Size and Velocity of Hydrometeors. J. Atmos. Ocean. Technol. 2000, 17, 130–139. [Google Scholar] [CrossRef]
  43. Gonzalez, S.; Bech, J.; Udina, M.; Codina, B.; Paci, A.; Trapero, L. Decoupling between Precipitation Processes and Mountain Wave Induced Circulations Observed with a Vertically Pointing K-Band Doppler Radar. Remote Sens. 2019, 11, 1034. [Google Scholar] [CrossRef] [Green Version]
  44. Wen, L.; Zhao, K.; Wang, M.; Zhang, G. Seasonal Variations of Observed Raindrop Size Distribution in East China. Adv. Atmos. Sci. 2019, 36, 346–362. [Google Scholar] [CrossRef]
  45. D’Adderio, L.P.; Vulpiani, G.; Porcù, F.; Tokay, A.; Meneghini, R. Comparison of GPM Core Observatory and Ground-Based Radar Retrieval of Mass-Weighted Mean Raindrop Diameter at Midlatitude. J. Hydrometeorol. 2018, 19, 1583–1598. [Google Scholar] [CrossRef]
  46. Berne, A.; Schleiss, M. Retrieval of the rain drop size distribution using telecommunication dual-polarization microwave links. In Proceedings of the 34th Conference on Radar Meteorology, Williamsburg, VA, USA, 5–9 October 2009. [Google Scholar]
  47. Chwala, C.; Kunstmann, H. Commercial microwave link networks for rainfall observation: Assessment of the current status and future challenges. Wiley Interdiscip. Rev. Water 2019, 6. [Google Scholar] [CrossRef] [Green Version]
  48. Rose, C.R.; Chandrasekar, V. A systems approach to GPM dual-frequency retrieval. IEEE Trans. Geosci. Remote Sens. 2005, 43, 1816–1826. [Google Scholar] [CrossRef]
  49. Vivekanandan, J.; Zhang, G.; Brandes, E. Polarimetric Radar Estimators Based on a Constrained Gamma Drop Size Distribution Model. J. Appl. Meteorol. 2004, 43, 217–230. [Google Scholar] [CrossRef] [Green Version]
  50. Rincon, R.F.; Lang, R.H. Microwave link dual-wavelength measurements of path-average attenuation for the estimation of drop size distributions and rainfall. IEEE Trans. Geosci. Remote Sens. 2002, 40, 760–770. [Google Scholar] [CrossRef]
  51. Jacques, T.; Stephane, O.; Blackb, R.A.; Amayencc, P.; Doud, X. The Concept of “Normalized” Distribution to Describe Raindrop Spectra: A Tool for Cloud Physics and Cloud Remote Sensing. J. Appl. Meteorl. 2001, 40, 1118–1140. [Google Scholar] [CrossRef]
  52. Cao, Q.; Zhang, G.; Brandes, E.; Schuur, T.; Ryzhkov, A.; Ikeda, K. Analysis of Video Disdrometer and Polarimetric Radar Data to Characterize Rain Microphysics in Oklahoma. J. Appl. Meteorol. Climatol. 2008, 47, 2238–2255. [Google Scholar] [CrossRef] [Green Version]
  53. Tang, Q.; Xiao, H.; Guo, C.; Feng, L. Characteristics of the raindrop size distributions and their retrieved polarimetric radar parameters in northern and southern China. Atmos. Res. 2014, 135, 59–75. [Google Scholar] [CrossRef]
  54. Zhang, G.; Vivekanandan, J.; Brandes, E.A. The Shape–Slope Relation in Observed Gamma Raindrop Size Distributions: Statistical Error or Useful Information? J. Atmos. Ocean. Technol. 2003, 20, 1106–1119. [Google Scholar] [CrossRef] [Green Version]
  55. Yu, M.; Guangheng, N.; Chandrasekar, V.; Fuqiang, T.; Haonan, C. Statistical characteristics of raindrop size distribution during rainy seasons in Beijing urban area and implications for radar rainfall estimation. Hydrol. Earth Syst. Sci. 2019, 23, 4153–4170. [Google Scholar] [CrossRef]
  56. Imai, J. Raindrop size distributions and the Z–R relationship. In Proceedings of the Eighth Weather Radar Conference, San Francisco, CA, USA, 11–14 April 1960; pp. 321–326. [Google Scholar]
  57. Puhakka. On the dependence of the Z–R relation on the temperature in snowfall. In Proceedings of the 16th Conference on Radar Meteorology, Houston, TX, USA, 22–24 April 1975; pp. 504–507.
  58. Gunn, K.L.S.; Marshall, J.S. The distribution with size of aggregate snowflakes. J. Meteorol. 1957, 15, 452–461. [Google Scholar] [CrossRef] [Green Version]
  59. Gray, D.M.; Male, D.H. Handbook of Snow: Principles, Processes, Management and Use; Pergamon Press: Oxford, UK, 1981; Volume 776, p. 776. [Google Scholar]
  60. Mallet, C.; Barthes, L. Estimation of Gamma Raindrop Size Distribution Parameters: Statistical Fluctuations and Estimation Errors. J. Atmos. Ocean. Technol. 2009, 26, 1572–1584. [Google Scholar] [CrossRef]
  61. Löffler-Mang, M.; Blahak, U. Estimation of the Equivalent Radar Reflectivity Factor from Measured Snow Size Spectra. J. Appl. Meteorol. 2001, 40, 843–849. [Google Scholar] [CrossRef]
  62. Garrett, T.J.; Fallgatter, C.; Shkurko, K.; Howlett, D. Fall speed measurement and high-resolution multi-angle photography of hydrometeors in free fall. Atmos. Meas. Tech. 2012, 5, 2625–2633. [Google Scholar] [CrossRef] [Green Version]
  63. Xichuan, L.; Taichang, G.; Lei, L. A video precipitation sensor for imaging and velocimetry of hydrometeors. Atmos. Meas. Tech. 2014, 7, 2037–2046. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The geographical position of the PARSIVEL disdrometer used in this study.
Figure 1. The geographical position of the PARSIVEL disdrometer used in this study.
Water 12 00920 g001
Figure 2. Accumulated raw particle counts by size and velocity for 15 precipitation types and their corresponding VD empirical relationship (blue lines) from Locatelli and Hobbs [36] and Atlas et al. [38].
Figure 2. Accumulated raw particle counts by size and velocity for 15 precipitation types and their corresponding VD empirical relationship (blue lines) from Locatelli and Hobbs [36] and Atlas et al. [38].
Water 12 00920 g002aWater 12 00920 g002b
Figure 3. The mean PSDs of rain, graupel, wet snow and dry snow.
Figure 3. The mean PSDs of rain, graupel, wet snow and dry snow.
Water 12 00920 g003
Figure 4. (a) The scatter distributions (gray triangle, asterisk, cross, square) and mean values (the length of the error bar is the standard deviation) of log10Nw-Dm of four types of precipitation. The dashed green line represents the fitting line for stratiform rainfall, and the dashed black rectangle represents the statistical area for continental and marine rainfall from Bringi et al. [28]; (b) The probability–density curves of Dm of four types of precipitation; and (c) The probability–density curves of log10Nw of four types of precipitation.
Figure 4. (a) The scatter distributions (gray triangle, asterisk, cross, square) and mean values (the length of the error bar is the standard deviation) of log10Nw-Dm of four types of precipitation. The dashed green line represents the fitting line for stratiform rainfall, and the dashed black rectangle represents the statistical area for continental and marine rainfall from Bringi et al. [28]; (b) The probability–density curves of Dm of four types of precipitation; and (c) The probability–density curves of log10Nw of four types of precipitation.
Water 12 00920 g004
Figure 5. The scatter distributions of μΛ (gray triangle, asterisk, cross, square) and μΛ fitted functions (yellow, cyan, wet snow and dry snow solid lines) of rain, graupel, wet snow and dry snow, respectively. The four dashed lines are derived from the μΛ relationship calculated by predecessors. The details are shown in Table 5.
Figure 5. The scatter distributions of μΛ (gray triangle, asterisk, cross, square) and μΛ fitted functions (yellow, cyan, wet snow and dry snow solid lines) of rain, graupel, wet snow and dry snow, respectively. The four dashed lines are derived from the μΛ relationship calculated by predecessors. The details are shown in Table 5.
Water 12 00920 g005
Figure 6. The scatter distributions of ZeS (gray triangle, asterisk, cross, square) and ZeS fitted functions (yellow, cyan, red and blue solid lines) of rain, graupel, wet snow and dry snow, respectively. In addition, the ZeS relationships (dotted or dashed lines) from previous studies are shown in this figure. The details are shown in Table 6.
Figure 6. The scatter distributions of ZeS (gray triangle, asterisk, cross, square) and ZeS fitted functions (yellow, cyan, red and blue solid lines) of rain, graupel, wet snow and dry snow, respectively. In addition, the ZeS relationships (dotted or dashed lines) from previous studies are shown in this figure. The details are shown in Table 6.
Water 12 00920 g006
Table 1. The precipitation days selected in this study. The corresponding ground temperature and precipitation time are also shown. The time used here is Beijing time, which is 8 h earlier than UTC.
Table 1. The precipitation days selected in this study. The corresponding ground temperature and precipitation time are also shown. The time used here is Beijing time, which is 8 h earlier than UTC.
Year DateHistorical Weather ObservationGround Temperature (°C)Precipitation Time
20140205Sleet/Small to moderate rain2~302:30~23:59
0207Sleet/Light snow−3~200:00~19:30
0212Sleet/Moderate snow0~420:30~23:59
0213Sleet/Cloudy−1~400:00~08:00
0218Heavy snow/Light to moderate snow−3~208:00~23:00
20150305Sleet/Overcast2~513:00~16:00
20160111Small rain/Sleet0~600:00~23:59
0120Sleet/Moderate to heavy snow0~419:30~23:59
0121Light snow/Light snow−1~100:00~07:30
0122Light to moderate snow/Light snow−5~209:30~16:00
0130Sleet/Overcast−1~408:30~10:00
20170208Light snow/Overcast−2~200:00~23:30
20180104Snowstorm/Moderate snow−1~100:00~06:00/15:30~23:59
0107Light rain/Sleet−1~400:30~14:00
0124Overcast/Heavy snow−1~316:30~23:59
0125Heavy snow to Snowstorm/Moderate snow−3~000:00~23:59
0126Overcast/Light snow−3~−100:00~03:00
0127Heavy snow to Snowstorm/Heavy snow−2~−106:00~23:59
0128Light snow/Overcast−7~−100:00~14:30
1208Overcast/Light snow0~2No precipitation
1230Light snow/Light snow−2~112:00~23:00
20190109Heavy snow/Sleet0~202:00~23:00
0130Light rain/Sleet1~800:00~08:30/17:30~23:59
0131Light snow/Overcast−3~400:00~08:00
0207Light rain/Moderate rain0~400:00~09:00/23:00~23:59
0208Light snow/Moderate snow−1~200:00~05:30/18:00~23:59
0209Overcast/Overcast−1~100:00~05:00/13:30~22:00
0210Light snow/Overcast−2~106:00~17:00
0221Light rain/Sleet2~700:00~02:00/16:30~23:59
0222Sleet/Overcast−1~600:00~13:00
Compared with the data recorded by the PARSIVEL disdrometer, the weather observation of the date corresponding to the shadow part in the table was incorrect. Light snow was observed for December 8, 2018, but the instrument did not record any precipitation. No nonliquid precipitation was observed for February 7, 2019, but the instrument recorded nonliquid precipitation. No precipitation was observed for February 9, 2019, but the instrument recorded a long period of precipitation.
Table 2. Empirical relationships of V (m/s)-D (mm) and M (mg)-D (mm) for different types of precipitation. The empirical relationships of graupel and snow are derived from Locatelli and Hobbs [36]. The empirical relationships of rain are derived from Atlas et al. [38]. The sample number N of each precipitation type and the total number Sum of samples for each category are also shown.
Table 2. Empirical relationships of V (m/s)-D (mm) and M (mg)-D (mm) for different types of precipitation. The empirical relationships of graupel and snow are derived from Locatelli and Hobbs [36]. The empirical relationships of rain are derived from Atlas et al. [38]. The sample number N of each precipitation type and the total number Sum of samples for each category are also shown.
CategoriesPrecipitation TypeV-DM-DNSum
RainRainV = 9.65 − 10.3exp(−0.6D)M = π/6D3.01119611196
GraupelLump graupel 1V = 1.16D0.46M = 0.042D3.01427521
Lump graupel 2V = 1.30D0.66M = 0.078D2.813
Lump graupel 3V = 1.50D0.37M = 0.140D2.77361
Conical graupelV = 1.20D0.65M = 0.073D2.64
Hexagonal graupelV = 1.10D0.57M = 0.044D2.91
SnowGraupellike snow of lump typeV = 1.10D0.28M = 0.059D2.14052557
Graupellike snow of hexagonal typeV = 0.86D0.25M = 0.021D2.494
Densely rimed dendritesV = 0.62D0.33M = 0.015D2.311
Densely rimed radiating assemblagesV = 1.10D0.12M = 0.039D2.11872
Unrimed side planesV = 0.81D0.99-1
Aggregates of unrimed radiating assemblagesV = 0.80D0.16M = 0.073D1.415
Aggregates of densely rimed radiating assemblages of dendrites or dendritesV = 0.79D0.27M = 0.037D1.918
Aggregates of unrimed radiating assemblages of plates, side planes, bullets, and columnsV = 0.69D0.41M = 0.037D1.918
Aggregates of unrimed side planesV = 0.82D0.12M = 0.040D1.4123
Table 3. Peak diameters, skewness and kurtosis of particle size distribution (PSD) for different types of precipitation.
Table 3. Peak diameters, skewness and kurtosis of particle size distribution (PSD) for different types of precipitation.
Precipitation Type Peak Diameter (mm) Skewness Kurtosis
Rain0.4372.324.12
Graupel0.5621.520.77
Wet snow0.9371.230.06
Dry snow0.6871.340.24
Table 4. The mean (ME) and standard deviation (SD) of log10N0, μ, Λ, log10Nw, Dm, S and Ze for different types of precipitation.
Table 4. The mean (ME) and standard deviation (SD) of log10N0, μ, Λ, log10Nw, Dm, S and Ze for different types of precipitation.
Precipitation Typelog10N0
(mm−1-μ m−3)
μΛ
(mm−1)
log10Nw
(mm−1 m−3)
Dm
(mm)
S
(mm h−1)
Ze
(dBZ)
MESDME SDMESDMESDME SD ME SDMESD
Rain7.864.47 7.58 6.35 15.53 12.433.410.53 0.99 0.53 0.73 2.559.58 11.01
Graupel4.222.72 2.90 4.984.536.713.300.532.581.25 2.07 2.93 23.40 13.87
Wet snow3.592.38 4.79 4.494.374.872.950.493.171.78 0.16 0.206.839.89
Dry snow2.902.64 4.10 5.383.213.592.990.633.492.00 0.27 0.409.678.89
Table 5. The μ–Λ relationships among different types of precipitation and previous studies.
Table 5. The μ–Λ relationships among different types of precipitation and previous studies.
Precipitation TypeSiteTimeμ-Λ
Rain Nanjing, East China Winter from Feb 2014 to Feb 2019μ = −0.0074Λ2 + 0.7564Λ − 1.5520
Graupelμ = −0.0073Λ2 + 0.8340Λ − 0.7664
Wet snowμ = −0.0066Λ2 + 0.8688Λ + 0.8648
Dry snowμ = −0.0106Λ2 + 0.9936Λ − 0.0820
Winter rain (Wen et al. [44]) Nanjing, East China2014 to 2016μ = −0.0110Λ2 + 0.8210Λ − 2.3700
Rain (Zhang et al. [54])Florida, America Summer of 1998μ = −0.0210Λ2 + 0.9880Λ − 2.6690
Rain (Cao et al. [52])Oklahoma, America2005 to 2007μ = −0.0201Λ2 + 0.9020Λ − 1.7180
Rain (Tang et al. [53])Beijing, North China Jul to Oct 2008Λ = 0.0075μ2 + 0.7230μ + 1.1721
Table 6. The ZeS relationship of different types of precipitation and previous studies. The previously studied ZeS relationship is converted from Ze = 0.189Z.
Table 6. The ZeS relationship of different types of precipitation and previous studies. The previously studied ZeS relationship is converted from Ze = 0.189Z.
Precipitation TypeZeS
RainZe = 311S2.34
GraupelZe = 252S2.03
Wet snowZe = 394S1.93
Dry snowZe = 145S1.33
Graupel (Gray and Male [59])Ze = 170S1.60
Snow (Gunn and Marshall [58])Ze = 378S2.00
Wet snow (Imai 1960 [56])Ze = 397S2.00
Dry snow (Imai 1960 [56])Ze = 101S2.00
Snow (Sekon and Srivastava [18])Ze = 337S2.21
Wet snow (Puhakka [57] )Ze = 302S2.00
Dry snow (Puhakka [57])Ze = 198S2.00

Share and Cite

MDPI and ACS Style

Pu, K.; Liu, X.; He, H.; Sun, Y.; Hu, S.; Wu, Y. Microphysical Characteristics of Winter Precipitation in Eastern China from 2014 to 2019. Water 2020, 12, 920. https://doi.org/10.3390/w12030920

AMA Style

Pu K, Liu X, He H, Sun Y, Hu S, Wu Y. Microphysical Characteristics of Winter Precipitation in Eastern China from 2014 to 2019. Water. 2020; 12(3):920. https://doi.org/10.3390/w12030920

Chicago/Turabian Style

Pu, Kang, Xichuan Liu, Hongbing He, Yu Sun, Shuai Hu, and Yi Wu. 2020. "Microphysical Characteristics of Winter Precipitation in Eastern China from 2014 to 2019" Water 12, no. 3: 920. https://doi.org/10.3390/w12030920

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop