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

: To improve solid precipitation monitoring in the hydrology and meteorology ﬁeld, 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 D m (mass-weighted mean diameter) and minimum log 10 N w (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 Z e -S (equivalent reﬂectivity factor–precipitation intensity) relationships among di ﬀ erent types of precipitation are signiﬁcantly di ﬀ erent. If only the Z e -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. S.H.; Y.S.; Writing—review


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  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

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.

Calculation of Integral Parameters
The number concentration N(D i ) (m −3 mm −1 ), which refers to the number of precipitation particles per unit volume per unit diameter, can be expressed as follows where D i (mm) is the median of the ith size bin, t (s) is the sampling time (t = 60 s in this study), A (mm 2 ) is the sampling cross section (A = 54 mm 2 for OTT PARSIVEL disdrometer), V j (m/s) is the median of the jth velocity bin, ∆D i (mm) is the width of ith size bin, and C ij 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] where D (mm) is the precipitation particle diameter, N 0 (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 The intercept parameter N 0 , shape parameter µ and slope parameter Λ can be expressed as follows where G = M 4 3 Water 2020, 12, 920

of 16
In addition, the normalized intercept parameter N w and the mass-weighted mean diameter D m are defined as follows [35] The precipitation intensity S (mm h −1 ) (liquid equivalent) can be calculated by the following formula where water density ρ w = 1 g cm 3 , and precipitation particle mass M(D i ) (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 Z e can be expressed as follows where D i (mm) is the equivalent ice sphere diameter and K i and K w are the dielectric constants for ice and water, respectively (|K i | 2 = 0.208 and |K w | 2 = 0.93 for the usual meteorological radar wavelength) [37]. Assuming D i corresponds to D i , then we can obtain the following expression where the ice density ρ i = 0.9167 g/cm 3 . Equation (11) is substituted into Equation (10) to obtain the following Equation [1] 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. Aggregates of densely rimed radiating assemblages of dendrites or dendrites

Categories Precipitation Type
Aggregates of unrimed radiating assemblages of plates, side planes, bullets, and columns Aggregates of unrimed side planes Water 2020, 12, 920 6 of 16

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 V-D relationships between different types of precipitation are obviously different, the V-D function fitted by 1-min average data can also be used to distinguish precipitation types [11,42,43]. In this study, the V fit = aD b 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 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.
Water 2020, 12, x FOR PEER REVIEW 6 of 16

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 V-D relationships between different types of precipitation are obviously different, the V-D function fitted by 1-min average data can also be used to distinguish precipitation types [11,42,43]. In this study, the Vfit = aD b 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 32 32 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.    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.  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. 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. Water 2020, 12, x FOR PEER REVIEW 8 of 16 Figure 3. The mean PSDs of rain, graupel, wet snow and dry snow.

The log 10 N w -D m 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.

The log 10 N w -D m Distributions of Different Types of Precipitation
The normalized intercept parameter log 10 N w and the mass-weighted mean diameter D m 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 log 10 N w versus D m 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 D m distribution of the three types of solid precipitation is much more widespread than that of rain, especially for dry snow.
The average log 10 N w -D m 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 log 10 N w -D m 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 log 10 N w and D m are shown in Table 4. Unfortunately, due to the lack of previous studies on the microphysical characteristics of solid precipitation, the log 10 N w -D m of graupel, wet snow and dry snow in this study cannot be compared with the statistical results in other regions. 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.  Z e (dBZ) ME SD ME SD ME SD ME SD ME SD ME SD ME SD

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.
Water 2020, 12, x FOR PEER REVIEW 10 of 16 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.

The Z e -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 Z e 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 Z e -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 Z e -S relationship. Table 6 records these specific parameters. As illustrated in this figure, when S > 1.1 mm h −1 , for the same Z e , the precipitation intensity of dry snow is greater than that of graupel and wet snow, and greater than that of rain. The Z e -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 Z e -S relationship of rain is uniformly used for the  Table 5.

The Z e -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 Z e . For rain, the equivalent reflectivity factor Z e 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 Z e = 0.189Z if the size of the hydrometeors is expressed as the equivalent ice sphere diameter [2]. Figure 6 shows the Z e -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, Z e = 0.189Z is used to convert this relationship into a Z e -S relationship. Table 6 records these specific parameters. As illustrated in this figure, when S > 1.1 mm h −1 , for the same Z e , the precipitation intensity of dry snow is greater than that of graupel and wet snow, and greater than that of rain. The Z e -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 Z e -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.
Water 2020, 12, x FOR PEER REVIEW 11 of 16 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 Z e -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 Z e -S relationship of snow from Gunn and Marshall [58] is also very similar to that of wet snow, while the Z e -S relationship of snow from Sekon and Srivastava [18] seems to be more similar to that of rain in this study. Furthermore, the Z e -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 Z e -S relationships between graupel from Gray and Male [59] and graupel in this study. Figure 6. The scatter distributions of Ze-S (gray triangle, asterisk, cross, square) and Ze-S fitted functions (yellow, cyan, red and blue solid lines) of rain, graupel, wet snow and dry snow, respectively. In addition, the Ze-S relationships (dotted or dashed lines) from previous studies are shown in this figure. The details are shown in Table 6.

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  Table 6.  [59]) Z e = 170S 1.60 Snow (Gunn and Marshall [58]) Z e = 378S 2.00 Wet snow (Imai 1960 [56]) Z e = 397S 2.00 Dry snow (Imai 1960 [56]) Z e = 101S 2.00 Snow (Sekon and Srivastava [18]) Z e = 337S 2.21 Wet snow (Puhakka [57] ) Z e = 302S 2.00 Dry snow (Puhakka [57]) The Z e -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 Z e -S relationship of snow from Gunn and Marshall [58] is also very similar to that of wet snow, while the Z e -S relationship of snow from Sekon and Srivastava [18] seems to be more similar to that of rain in this study. Furthermore, the Z e -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 Z e -S relationships between graupel from Gray and Male [59] and graupel in this study.

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.

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, µ-Λ, log 10 N w -D m and Z e -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 Z e -S relationships of different types of precipitation. If only the Z e -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.