Classiﬁcation of Rainfall Types Using Parsivel Disdrometer and S-band Polarimetric Radar in Central Korea

: Tools to identify and classify stratiform and convective rains at various times of the 12 days from June 2015 to March 2016 in Jincheon, Korea, were developed by using a Parsivel disdrometer and S-band polarimetric (S-POL) radar data. Stratiform and convective rains were identiﬁed using three different methods (vertical proﬁle of reﬂectivity (VPR), the method proposed by Bringi et al.[1] (BR03), and a combination of the two (BR03-VPR)) by using a Parsivel disdrometer for its applications to radar as a reference. BR03-VPR exhibits a better classiﬁcation scheme than the VPR and BR03 methods. The rain types were compared using the drop size distribution (DSD) retrieved from polarimetric variables and reﬂectivity only. By using the DSD variables, a new convective/stratiform classiﬁcation line of the log-normalized droplet number concentration (log 10 N w ) − median volume diameter ( D 0 ) was derived for this area to classify the rainfall types using DSD variables retrieved from the polarimetric radar. For the radar variables, the method by Steiner et al.[2] (SHY95) was found to be the best method, with 0.00% misclassiﬁcation of the stratiform rains. For the convective rains, the DSD retrieval method performed better. However, for both stratiform and convective rains, the fuzzy method performed better than the SHY95 and DSD retrieval methods.


Introduction
Rainfall is generally classified into stratiform and convective regimes [3][4][5] based on the spatial and temporal scales and vertical velocity of the cloud system [6]. Stratiform rains have relatively weak vertical velocity fields, greater horizontal homogeneity, and lower rainfall intensity, whereas convective rains are associated with strong vertical velocity fields, low areal coverage, and high rainfall intensities [3,7]. These parameters are important for understanding cloud physics, as stratiform and convective rains are characterized by different rainfall growth mechanisms [8], and for understanding radar-rainfall measurements [2,9]. Moreover, they also play an important role in quantitative precipitation estimation from both ground-and space-based instruments [8].
Korea, in particular, is still lacking in detailed DSDs and rainfall classification information. A few previous studies have concentrated on the south of Korea, either on rainfall classification or rainfall characteristics [e.g., 6,[42][43][44]. However, a comparison between DSD parameters by rainfall types is needed to understand the overall characteristics of rainfall in the central part of Korea. There are only a few studies on the classification of rainfall types using DSD parameters retrieved from polarimetric variables. This study compares the existing classification methods and proposes a new algorithm or separation method by using a Parsivel disdrometer and S-band polarimetric (S-POL) radar at Jincheon (central region of South Korea). The DSD retrieval and fuzzy method proposed in this study can be used to classify stratiform and convective rains by using the S-POL radar. The next section briefly describes the instruments, data, and analysis employed in this study. Sections 3 discusses results of the classification. Finally, Section 4 summarizes and concludes this study.

Datasets
The rain gauge, Parsivel disdrometer, and S-POL radar data from sites in Jincheon (36.98 • N, 127.44 • E) and Yongin (37.21 • N, 127.29 • E), shown in Figure 1, were acquired from June 2015 to March 2016. The rainfall data from Jincheon were collected and used to evaluate the accuracy of the radar rainfall from a tipping bucket rain gauge, which is operated and quality-controlled by the Korea Meteorological Administration (KMA). Ten days of 2015 and two days of 2016 were selected for a comparison of DSD parameters using the Parsivel disdrometer and S-POL radar, as presented in Table  1. The total daily rainfall for the 12 selected days was 385.4 mm in 155 h at the Jincheon site, and the cases were selected as those with daily rainfall larger than the threshold of 15 mm. Details can be found in Loh et al. [45].
The DSD values obtained using the Parsivel disdrometer from the Jincheon site, located 28.62 km from the S-POL radar site, were analyzed and compared with the radar retrievals in a statistical and functional approach. Briefly, the Parsivel disdrometer is a modern, laser-based optical disdrometer that measures the size and fall velocity of all liquid and solid precipitation [17,46]. It can measure 32 size bins from 0.062 to 24.5 mm [6]. A detailed description of the disdrometer can be found in Löffler-Mang and Joss [46]. A total of 17,280 1-min DSD spectra, observed from the Parsivel disdrometer for 12 days at Jincheon, were used to determine the relationship between polarimetric variables for the retrieval of DSD parameters from S-POL variables and the Parsivel disdrometer.
The S-POL radar data used in this study are from the Yongin site, operated by the Weather Radar Center of KMA ( Figure 1). The S-POL radar data were updated every 10 minutes with 11 elevation angles (0.20 •   where v is in m s −1 , and D is the diameter of an equivalent-volume sphere, in mm, was applied to filter unreliable data. Kruger and Krajewski [48] found that the occurrence of outliers was larger as compared to Atlas et al. [47]. Therefore, a fall-velocity-based filter was developed where V measured (m s −1 ) is the observed fall velocity and V A (m s −1 ) is v from Equation (1) used to remove the outliers of the data affected by the wind-caused turbulence [49]. In addition, 1-min rainfall rates were neglected at less than 0.1 mm h −1 (to eliminate acoustic and wind noise [9,44]) and greater than 200 mm h −1 (to enhance the capability of the detection of small raindrops (< 1 mm) [44]).

DSD parameters
A three-parameter (µ, Λ, and N 0 ) normalized gamma DSD model was used in this study. The values of Λ and µ were obtained through the untruncated-moment method. The integration of most moment calculations is usually performed over the size range from zero to infinity as where n is the order and Γ(µ + n + 1) is the complete gamma function [6,25]. The parameters µ, Λ, and N 0 can be obtained from any three moments (e.g., 2 nd , 4 th , and 6 th ). The ratio of the three moments, η, is defined as From Equation (4), the values of µ, Λ, and N 0 can be computed as [ 6,17,25,50]. Then, the value of D 0 can be obtained from µ (Equation (5)) and Λ (Equation (6)) through [50].
In fact, D m is very close to D 0 , which is defined as the ratio of the fourth to the third moment of the DSD: The relation of D m with D 0 for a gamma distribution was also shown by Ulbrich [19] as below: The liquid water content (LWC; g m −3 ) was calculated as where ρ w (g m −3 ) is the water density [1,6,17]. The normalized intercept parameter N w (mm −1 m −3 ) is estimated from D m (Equation (10)) and LWC (Equation (11)) [1,6,17].

Quality control
The Z H and Z DR bias correction schemes were not applied but a differential phase shift (Φ DP ) unfolding stage and noise removal stage based on study of You et al. [51] was applied. A detailed description can be found in You et al. [51]. After the noise removal from Φ DP , a specific differential phase (K DP ) was obtained from the slope of 9 and 25 gates of quality-controlled Φ DP , with a gate size of 250 m. The lightly filtered estimate of K DP was selected when Z H ≥ 40 dBZ, whereas when Z H < 40 dBZ, the heavily filtered estimate was used for any particular range gate [51]. The calculation of the specific attenuation at horizontal polarization (A H ) was based on the method of Ryzhkov et al. [52] and You and Lee [53]. A H can be calculated from the attenuated reflectivity (Z a ) radial profile and two-way Path Integrated Attenuation (PI A) along the propagation path (r 1 ,r 2 ). Details of the A H calculation can be found in You and Lee [53].

DSD parameters
The retrieval of DSD parameters from polarimetric radar is important for improving the accuracy of rainfall estimations [23]. According to Brandes et al. [23] (BRA04), Z H (mm 6 mm −3 ), Z V (mm 6 m −3 ), Z DR (dB), and K DP (deg km −1 ) are the most important factors in quantitative rain estimation among polarization radar parameters and depend on the drop scattering amplitudes and DSD. Similarly, A H (dB km −1 ) is also important and is given in terms of N(D) (mm −1 m −3 ) and the total extinction cross-section (σ e (D)) (mm 2 ) as described by Bringi et al. [54]. The values of Z H , Z DR , K DP , and A H can be obtained using the transition-matrix (T-matrix) scattering method, which is required to obtain information such as raindrop shape, canting angle of raindrop, frequency, temperature, and DSDs to calculate radar variables. The following raindrop shape assumptions are combined with respect to the rain drop size used for the calculation of variables from the DSDs: b/a = 1.0048 + 0.500057D − 0.02628D 2 + 0.003682D 3 − 0.0001677D 4 (13) b/a = 1.012 − 0.01445D − 0.01028D 2 (14) where a, b, and D are the major axis, minor axis, and equivolume diameter of the raindrops in millimeters, respectively. Equation (13) was proposed by Beard and Chuang [55] for the equilibrium raindrop axis ratio derived from a numerical model. Later, Andsager et al. [56] found that raindrop shapes with a diameter between 1.1 and 4.4 mm in turbulent flow are better explained by Equation (14). Bringi et al. [1] proposed the combined drop shape assumption using Eqs. (13) and (14) to better represent rain drops in nature. You and Lee [53] found that combining Equation (13) for raindrops smaller than 1.1 mm and larger than 4.4 mm with Equation (14) for the raindrop diameter between 1.1 and 4.4 mm gave more accurate rainfall estimates in Korea, compared with other raindrop shape assumptions. Other variables in the T-matrix calculations include the temperature, frequency and distribution of the canting angles of rain drops. The temperature is assumed to be 20 • C, the frequency is 2.87 GHz and the distribution of their canting angles is Gaussian with a mean of 0 • and a standard deviation of 7 • , as determined by Huang et al. [57]. The T-matrix scattering method [58,59] was used to compute the polarimetric variables and relationships between D 0 and Z DR , and N w and Z H , using DSDs measured by a Parsivel disdrometer to determine the reference classification method for rainfall types using radar variables.

Methods for rainfall classification using disdrometer and S-POL radar
The rain types are classified by disdrometer and S-POL radar after identification through the BR03, vertical profile of reflectivity (VPR), and BR03-VPR methods.
2.4.1. General identification methods

The BR03 method
This scheme was proposed by BR03 and has been applied in numerous studies [e.g., 17,21,39,44]. BR03 used a standard deviation for R (σ R ) of 1.5 mm h −1 , as the threshold for the classification of convective and stratiform rains, based on data acquired by the 2DVD in Colorado, with a radar-observed bright band signature during a stratiform upslope event. For stratiform rains, the criteria were based on σ R ≤ 1.5 mm h −1 and R ≥ 0.5 mm h −1 . Otherwise, the rainfall was considered as convective rain if σ R > 1.5 mm h −1 and R ≥ 5 mm h −1 , which was adopted and conducted in every minute in this study, as shown in the dashed box (I) of Figure  2. This method is a disdrometer-based method applied to the Parsivel data. In addition, the remaining rainfall data (total samples of the Parsivel data) that belong neither to stratiform nor convective rains were categorized as unclassified rains. For the BR03 method, 1301 (382) samples of stratiform (convective) rains from the Parsivel disdrometer were selected at the Jincheon site.

The VPR method
The bright band is a radar signature of the melting layer, which is generally identified with stratiform rains [2,9,12]. Convective rains typically can be distinguished by considering the threshold in radar reflectivity. Almost all precipitation with a radar reflectivity above 40 dBZ was identified as convective rain [e.g., [2][3][4]9,15], with different conditions. For instance, Steiner et al. [2] classified convective rains based on the intensity (reflectivity of at least 40 dBZ) and sharpness of the peaks of echo intensity. However, other thresholds have also been used to identify convective rains. For example, CA08 considered the rainfall as convective rain if R was greater than 10 mm h −1 , and the reflectivity was higher than 38 dBZ. A threshold of 30 dBZ was used by Niu et al. [60] to classify convective rains that resulted from the weaker precipitation system in a semi-arid continental regime [61].
Stratiform and convective rains can also be identified through the VPR method, as shown in the box II of Figure 2. This method was conducted every 10 minutes. Stratiform rains can be distinguished by the presence of a bright band [e.g., 2,3,6,12,18], where the bright band was determined subjectively in this study, and rains with a reflectivity of 30 dBZ or higher are identified as convective [e.g., 60,62] by the VPR method. According to Zhang et al. [61], the maximum height of convective rains is between 7-8 km and the height is typically 5-6 km. Franco et al. [63] introduced the 30 dBZ echo-top altitude as a discriminatory variable for classifying convective rains and indicated that a Z H value ≥ 30 dBZ at altitudes higher than 1 km can be considered as convective rains for the low-convection scenario, although the maximum height for the high-reflectivity scenario (i.e., > 40 dBZ) would be the same height as that of the stratiform bright band. Therefore, weak convective rains with Z H of 30-35 dBZ at low echo-top heights (≤ 5 km) were also considered as convective rains in this study. The durations selected for all rainfall types at the Yongin site are summarized in Table 1.
This method is a radar-based method applied to the S-POL data. For the VPR method, 1798 (451) samples of stratiform (convective) rains were selected at Jincheon. In the current study, time series of reflectivity were analyzed to classify the different rain types at the Jincheon site ( Figure 3). The vertical profiles of reflectivity (VPRs) of radar were calculated by 11 PPI scan data. The elevation angles within one volume in 10 minutes are 0.2, 0.6, 1.1, 1.8, 2.8, 4.2, 6.2, 9.1, 13.2, 19.0, and 80.0 degrees. Considering these elevation angles apart from 80.0 degrees and the distance between the radar and Jincheon site, the averaged vertical distance between grid points is approximately 1.12 km. Figure 3 shows the time-height cross-sections of reflectivity from the S-POL radar at the Yongin site for 12 different dates. The reflectivity values of the selected convective rains (Table 1) Table 1) by virtue of the presence of the bright band ( Figure 3). Figure 3b indicates that the strong bright band was detected at a height of 4-5 km from 0900 to 1800 LST on Jul. 12, 2015. A particularly weak-to-moderate bright band was also observed at a height of 5-6 km from 0520 to 0920 LST on Aug. 16

The BR03-VPR method
A unified method, BR03-VPR, was applied in this study. The BR03-VPR combines the methods of BR03 and VPR. The VPR method may be less effective for rainfall classification owing to human error. Human error refers to the error of the individual analyst and does not correspond to the method or the procedure. Additionally, the VPR method is based on the radar data collected every 10 minutes, whereas the BR03 method is based on the Parsivel disdrometer data, which are updated every minute. First, the rainfall types were identified using the BR03 method (dashed box in Figure 2). Next, the rainfall types were distinguished by the appearance of the bright band and threshold values of radar reflectivity, as shown by the dotted box in Figure 2. This method was used as a reference to evaluate the accuracy of classifying rainfall types based on S-POL radar (especially for DSD retrieval and the fuzzy method as described in sub-section 2.4.3). A total of 1562 (1224 stratiform and 338 convective rains) 1-min DSD spectra were sampled by the BR03-VPR method in this study. The data was only selected when the resulting outputs, from previous two methods, were both identified as either stratiform or convective rains. Moreover, the DSD parameters for the Jincheon site were extracted in order to compare with the DSD retrieval and fuzzy methods.

Classification with DSDs variables measured by Parsivel disdrometer
Several rainfall classification methods using DSD parameters acquired either from disdrometer measurements or polarimetric radar retrievals have been proposed [6]. For example, BR09 found a classification method using a relationship between D 0 and log 10 N w at Darwin, Australia. CA08 suggested a separation method for Italy that used the reflectivity and R threshold from JWD and Pludix instruments with a relationship between logN 0 and Λ. An empirical classification method for Kapingamarangi Atoll demonstrating log 10 N 0 − R and Λ − R relationships has been proposed by TS96. YOU16 developed a new classification boundary condition, based on heavy rainfall from the Parsivel disdrometer over Busan, Korea after comparing the methods of BR09, CA08, and TS96. The four separation methods, including YOU16, were compared in this study to determine an appropriate separation line for both rainfall types corresponding to the S-POL radar data.

Classification by S-POL radar 1. DSD retrieval method
Significant progress has been made in retrieving DSD parameters from polarimetric radar measurements, quantitatively comparing the radar retrievals with disdrometer measurements, and developing DSD models [e.g., 1,6,13,23 -25]. Two approaches have commonly been suggested for the retrieval of a normalized gamma distribution.
The β (slope) method, proposed by Gorgucci et al. [64,65] and Bringi et al. [66] (BR02), retrieves N w and D 0 from Z H , Z DR , K DP , and β. β is given by where r is the axis ratio (b/a; a and b are the major and minor axes of the spheroid, respectively) and D is the equivolumetric spherical diameter in mm [64]. The goal of this method is to retrieve the DSD parameters from Z H , Z DR , K DP and β [23,64,67]. The advantage of the β method is that the three parameters of the normalized gamma DSD can be retrieved independently [67]. However, the errors in the modeled DSD and K DP estimations restrict the application of the method to cases with high R values [67], and thus the retrievals from Z H , Z DR , and K DP may be dependent [68].
The constrained-gamma method retrieves the DSD parameters based on Z H and Z DR , as proposed by Brandes et al. [23], Zhang et al. [69], Brandes et al. [70]. Brandes et al. [70] also suggested an empirical relationship (Λ = 0.0365µ 2 + 0.735µ + 1.935), which essentially reduces the three-parameter normalized gamma DSD to a two-parameter model. Some studies have shown that the constrained-gamma method is better than the β method for DSD retrieval [e.g., 23,24,29] despite having some uncertainties [71].

Fuzzy logic classification method
The fuzzy logic theory was first developed by Zadeh [72] to provide a scheme for handling issues due to the indefiniteness arising more from an intrinsic ambiguity than from a statistical variation [73]. The fuzzy set allows an element to partially belong to a set indicated by a membership degree of 0 to 1 through a membership function (MF). While 0 is defined as an element that is not included in the set at all, 1 is an element that belongs to the set [74].
Many studies have shown that a valid approach for hydrometeor classification is to use a fuzzy logic algorithm [e.g., 31,36,39,40,75]. Different hydrometeor classes were described by one-dimensional (F (i) (V j ) = P (i) (V j ), where V j is the j th additional radar variable) or two-dimensional (F (i) (Z H , V j ) = P (i) (Z H )P (i) Z (V j )) MFs based on the fuzzy logic methodology [76]. The one-dimensional unconditional MFs P (i) (Z H ) and P (i) Z (V j ) characterize distributions of Z H and V j , respectively, for the i th class. The MFs P (i) Z (V j ) characterize the conditional distribution of the variable V j for the i th class for a given Z H . The product of P (i) (Z H ) and P (i) Z (V j ) represents the two-dimensional MF characterizing the joint distribution of Z H and V j in the Z − V j plane for the i th class.
The maximal aggregation value is used to identify the hydrometeor class. For each class, the aggregation was defined as where M is the number of variables and W j is the weight assigned to the j th variable [76]. Schuur et al. [76] used five radar parameters (Z H , Z DR , co-polar correlation coefficient (ρ HV ), a texture parameter SD(Z H ) of the Z H field, and a texture parameter SD(Φ DP ) of the Φ DP field) for automatic classification. Furthermore, Mahale et al. [41] identified three-body scattering by the fuzzy logic classification of S-POL radar echoes following Park et al. [31]. They defined A i as where W ij is a weight between 0 and 1 assigned to the i th class and j th variable and P (i) (V j ) is a trapezoidal MF for the i th class and j th variable. Therefore, this may be also a suitable approach for rainfall classification, as it can accurately classify the hydrometeor. However, for rainfall classification, only a few studies have classified rain types through a fuzzy logic approach [e.g., 8,26].
The fuzzy logic approach in the current study is based on the observations of Z H , Z DR , K DP , and A H , which are highly sensitive to droplet shape, size, orientation, and composition [34,[77][78][79][80], for stratiform and convective rain classification. Hence, A i can be defined as from the Equation (17). The trapezoidal MF is defined by a < b < c < d (Figure 4), and is expressed as follows The criteria for the limits from the Parsivel data at the Jincheon site used in the trapezoidal MF are the 0.5 th , 20 th , 80 th , and 99.5 th percentiles, based on Mahale et al. [41] to obtain the values of the four limits. Details of the vertex values (a, b, c, and d) are shown in Table A1 and A2. Figure 5 shows a flowchart classifying stratiform and convective rains through the trapezoidal MF using S-POL radar data from the Yongin site. The radar classification with the trapezoidal MF method can be performed after obtaining the values of the four limits from the Parsivel data. First, the stratiform rains can be identified when Z H < 30.02 dBZ, whereas the convective rains can be verified if Z H > 39.96 dBZ. According to Hwang [81], the summation of inputs for the fuzzy logic method is advantageous, because other inputs can compensate, particularly for the cases that have weak characteristics in one input. Hence, this summation of inputs is applied in the current study for the intercept of Z H between 30.02 and 39.96 dBZ ( Figure   5). The stratiform rains can be identified if P

The characteristics of DSDs
The averaged droplet number concentration of each rain type during the selected period (Table 1) are shown in Figure 6. Notably, all the convective rains of the three different methods (VPR, BR03, and BR03-VPR) contain more raindrops than the other types (stratiform and unclassified rains) across all drop sizes. Furthermore, stratiform and unclassified rains were found to have a similar pattern, but with minor differences in magnitude, particularly for small drops (0.31−0.65 mm), as shown in the bottom left panel of Figure 6. The DSDs were distributed in a concave downward manner and shifted toward larger drops for more abundant forms in all the rainfall types. The peak number concentration of all convective rains occurred at 0.44 mm, whereas almost all the stratiform rains showed a peak number concentration of 0.56 mm, except for the BR03 method (0.69 mm). This behavior is probably due to the dead time problem of the disdrometer originating from the insensitivity of the instrument to small drops [16,45].
The mean R (mean D) were 1.98 mm h −1 (1.67 mm), 24.43 mm h −1 (2.34 mm), and 0.41 mm h −1 (1.41 mm) for the averaged stratiform rains of the three different methods (VPR, BR03, and BR03-VPR), averaged convective rains (same methods), and unclassified rains of the BR03 method, respectively. The maximum log 10 N w of the stratiform rains (5.10 mm −1 m −3 ) was higher than that of the convective rains (4.93 mm −1 m −3 ) in the BR03 method (Table 2). This might have occurred when large drops formed below the bright band in the stratiform rains [6]. However, Moreover, the mean of the entire dataset was less than that of the convective rains, as the stratiform rains had an abundance of small drops. Furthermore, the DSD of the convective rains (D max = 7.50 mm) was much broader than that of the stratiform rains (D max = 4.25 mm). Figure 7 shows the average DSDs of the different R classes for stratiform and convective rains. The number of rainfall categories for the stratiform rains (six classes: R ≤ 2 mm h −1 , 2 < R ≤ 4 mm  h −1 , 4 < R ≤ 6 mm h −1 , 6 < R ≤ 10 mm h −1 , 10 < R ≤ 20 mm h −1 , and 20 < R ≤ 40 mm h −1 ) was lower than that for the convective rains (seven classes: R ≤ 2 mm h −1 , 2 < R ≤ 4 mm h −1 , 4 < R ≤ 6 mm h −1 , 6 < R ≤ 10 mm h −1 , 10 < R ≤ 20 mm h −1 , 20 < R ≤ 40 mm h −1 , and R > 40 mm h −1 ). The convective rains had higher R (R > 40 mm h −1 ) than the stratiform rains. Moreover, the R values for all the stratiform rains were found to be in the lowest category (R ≤ 2 mm h −1 ), as shown in Figure 7a. However, only the VPR method existed for the convective rains (Figure 7b). It was expected that the R of the stratiform rains would be lower than that of the convective rains [e.g., 15,82]. For the VPR method, the number concentration increased with R only when 1.06 < D ≤ 3.75 mm for the stratiform rains. In contrast, the number concentration in the convective rains increased with R only when 1.06 < D ≤ 2.38 mm and D ≥ 4.75 mm. For the BR03 method, the number concentration in the convective rains increased for all diameter bins when R increased for smaller drops (D < 1 mm), whereas for medium drops (1 ≤ D < 3 mm), almost all types showed an increase in number concentration for all diameter bins with increasing R for the stratiform and convective rains. Finally, for the BR03-VPR method, the number concentration of the stratiform rains increased with R only when D ≤ 3.75 mm, whereas for the convective rains, the number concentration increased with R only when 1.06 < D ≤ 3.25 mm and D ≥ 4.75 mm.

Relationship between DSD and R
For stratiform rains, the BR03-VPR method showed that the average drop size spectra increased with increasing R across small and medium drops. Furthermore, for medium drops, the VPR method of the stratiform rains showed that average DSDs increased with increasing R. However, for the convective rains, the BR03 method showed that the average drop size spectra increased with increasing R across small and medium drops. In addition, all the methods of the convective rains indicated that the average DSDs increased with increasing R across all the diameter bins for medium drops, except for the VPR method. Figure 8 presents the size distribution of the unclassified rains. The figure demonstrates the average DSDs according to three classes (R ≤ 2, 2 < R ≤ 6, and 4 < R ≤ 6) for the unclassified rains according to the BR03 method over small, medium, and large drops. The number concentration increased with increasing R only when the diameter in the ranges D ≥ 0.69 mm. . However, the classification lines of BR09 and YOU16 failed to separate the stratiform and convective rains correctly for all three methods (VPR, BR03, and BR03-VPR). The YOU16 classification line seems to provide a better separation than BR09 (Table 3), which may be because the rainfall types are from the same weather systems. Table 3 shows the percentages of misclassification of both rainfall types based on the D 0 − log 10 N w , log 10 N 0 − Λ, R − log 10 N 0 , and Λ − R separation lines.

Comparison of rainfall classification
The BR09 classification line separated the stratiform rains better (with a misclassification of 0.00 to 1.61%) than the convective rains (with a misclassification of 23.67 to 36.14%). Hence, a new classification line was proposed for the Jincheon site based on visual examination to separate the stratiform and convective rains more accurately, log 10 N w = −1.09D 0 + 5.3.
(20) Figure 8. Same as Figure 7, but for unclassified rains. Table 3. Percentages of the misclassification of stratiform and convective rains based on the relations of D 0 − log 10 N w , log 10 N 0 − Λ, R − log 10 N 0 , and Λ − R for three different identification methods (VPR, BR03, and BR03-VPR) The Jincheon separation line (blue line in Figure 9a) could classify both rainfall types more accurately compared with BR09 and YOU16 for all three methods. Table 3 also demonstrates that the Jincheon separation line performed better than the others, with lower misclassification rates for both stratiform (1.23−6.99%) and convective (1.78−5.76%) rains. The BR03-VPR method had the lowest misclassification rates for stratiform and convective rains, at 1.23% and 1.78%, respectively among the three methods. Figure 9b shows scatter plots of log 10 N 0 versus the slope parameter (Λ), which referred to the CA08 study. In general, the CA08 and YOU16 classification lines failed to separate the stratiform and convective rains. Table 3 also shows the high degree of misclassification, especially for convective rains, at ∼53% for the VPR, BR03, and BR03-VPR methods. The proposed classification line is very similar to that presented for CA08 and YOU16, based on the visual judgement of the Jincheon data, as shown below Λ + 4.17 = 2.22 log 10 N 0 .

Relations
However, the Jincheon separation line also failed to classify both rainfall types at the Jincheon site. It is difficult to accurately derive the log 10 N 0 − Λ relationship because of the overlapping of stratiform and convective rains, as shown in Figure 9b. YOU16 indicated that it is not possible to cleanly classify the convective and stratiform rains because of the overlapping in log 10 N 0 − Λ space. Additionally, they also proposed that the log 10 N 0 − Λ domain by the CA08 classification line is unsuitable for Korea. The scatter plots of log 10 N 0 against R for the Jincheon site are shown in Figure 9c. The R values for the stratiform and convective rains were in the ranges of 0.0−27.5 and 0.0−169.8 mm h −1 , respectively, for the VPR method. For the BR03 and BR03-VPR methods, the R for the stratiform and convective rains were in the ranges of 0.5−14.6 and 5.0−169.8 mm h −1 , respectively. Moreover, the values of log 10 N w for the stratiform and convective rains were in the ranges of 2.7−33.4 and 2.8−23.8 mm −1 m −3 , respectively. The TS96 and YOU16 classification lines failed to separate the stratiform and convective rains (Figure 9c). Table 3 also shows that the YOU16 separation line could classify the stratiform rains better than the TS96 classification line. Hence, a new separation line was created for the Jincheon site to accurately classify the two rain types. The Jincheon separation line based on visual examination is shown below The newly derived separation lines (Equation (22)) were able to clearly separate the stratiform and convective rains. Table 3 also demonstrates that the Jincheon separation line performed better with lower misclassification rates, especially for stratiform rains (3.35−8.76%). The BR03-VPR method had the lowest misclassification rates of 1.18% and 3.35% for convective and stratiform rains, respectively. Figure 9d presents scatterplots of Λ against R for the Jincheon site. Again, TS96 and YOU16 failed to classify the two rain types. The YOU16 separation line classified the rainfall types better than TS96 because of the different microphysical processes of different regions, as the TS96 separation line is applied to the rainfall in the Kapingamarangi Atoll [6]. Table 3 also indicates that the TS96 and YOU16 classification lines had higher misclassification rates, especially for convective rains. Thus, the newly proposed separation line based on the visual examination of the data is Equation (23) accurately separates the stratiform and convective rains, particularly for the BR03-VPR method. This method had the lowest misclassification rates of 0.90% and 5.03% for stratiform and convective rains, respectively. The new proposed classification line could accurately separate the stratiform and convective rains for all of the relationships, except for the log 10 N 0 − Λ space. The proposed R − log 10 N 0 classification line was able to accurately classify the convective rains, whereas for the stratiform rains, the D 0 − log 10 N w classification line was suitably accurate. The DSD parameters D 0 and log 10 N w were retrieved from the S-POL radar for the classification of rainfall types as N w and D 0 can be retrieved from polarimetric variables.

DSD retrieval
Data from all 12 days were used to classify the rainfall types using the S-POL radar. YOU16 indicated that the type of rainfall event occurring over Korea can be assumed by the values of R, D 0 , log 10 However, Table 4 shows that the relationship between Z H /N w − D 0 was not as strong as the relationship between D 0 and Z DR . Further studies on the relationship between Z H /N w and D 0 are needed for a more detailed understanding.
BRA04 indicated that the enabling of correlations between DSD parameters for DSD retrieval from a pair of independent remote measurements such as Z H and Z DR , is useful for reducing the

Equations
Power-law relations Polynomial function was obtained, as shown in Figure 11. Nonetheless, D 0 was highly correlated with Z DR , yielding correlation coefficients of 0.81 (Table 4). Table 4 also shows that D 0 had a stronger correlation with Z DR as compared to BRA04.

Fuzzy logic classification
All parameters had positively skewed distributions for both rain types ( Figure 12). The skewness for convective rains was higher than that for stratiform rains, except for Z DR . The skewness was the largest for the K DP distribution for convective rains, at 2.85. This was followed by the skewness values of 2.71, 1.41, and 0.44 for the convective rains of A H , Z DR , and Z H , respectively. Similarly, for stratiform rains, the distribution of K DP had the highest skewness as compared to the others, followed by the skewness values of 1.71, 1.54, and 0.33 for the convective rains of Z DR , A H , and Z H , respectively. The Z H value (Figure 12a) showed that Z H <29.30 dBZ indicated stratiform rains, and Z H >37.42 dBZ indicated convective rains. Figure 12b also shows that stratiform rains can be identified when Z DR <0.19 dB, whereas convective rains correspond with Z DR >1.13 dB. For K DP , stratiform and convective rains can also be identified when K DP <0.034 and K DP >0.098 deg km −1 , respectively. Also, the presence of convective rains can be confirmed when A H >0.0022 dB km −1 , and rains are classified as stratiform when A H <0.0015 dB km −1 . The trapezoidal MF can identify stratiform and convective rains easily and effectively, especially for Z H . However, clearly distinguishing both rain types only by way of the intercept is difficult. Therefore, the same steps are applied again for the intercepts of all the DSD parameters based on the interval of Z H , where 30 ≤ Z H < 40, as shown in Figure 12. Figure 13 displays the trapezoidal MF and distribution of Z H , Z DR , K DP , and A H based on 30 ≤ Z H < 40 for both stratiform and convective rains, through the BR03-VPR method. Again, all the parameters had positively skewed distributions for convective and stratiform rains, as shown in Figure 13. The trapezoidal MF for Z H , K DP , and A H has a clear triangular intercept.

Comparison
Three different methods were compared in this study. The first method was given by SHY95, which is based on the analysis of the horizontal gradient of the radar reflectivity field. The method considers rains as convective if Z H ≥ 40 dBZ or greater than the fluctuating threshold, depending on the area-averaged background reflectivity for any grid point within the radar scan radius. The DSD retrieval method uses Equation (25) and (26) to calculate the retrieved N w and D 0 , respectively, as retrieved from the observed Z H and Z DR of S-POL radar. The third method of classification, the fuzzy method, has also compared in this section.
For example, during Jun. 26, 2015 from 0130 to 0200 LST (not shown), the SHY95 method failed to identify the convective rains at the Jincheon site, unlike the DSD retrieval and fuzzy methods. This may be because the classification criteria for SHY95 are not suitable to Korea. However, DSD retrieval and fuzzy methods can classify the convective rains with a similar pattern at the Jincheon site. Table 5 displays a comparison of the precision for stratiform and convective rain classification in percentages with the three different methods from the S-POL radar, for all 12 cases that were selected. For the stratiform rains, all the methods classify the rains, particularly the SHY95 method, with 85.47% for stratiform and 0.00% for convective rains. In contrast, the DSD retrieval method showed the lowest (higher) percentages, 61.45% (24.02%) for the appearance stratiform (convective) rains. The new fuzzy method for stratiform rains performed well with 79.61% of stratiform rains and 5.87% of convective rains, which is as good as the SHY95 method.
Meanwhile, the DSD retrieval method had a higher percentage (68.83%) for the convective rains as compared with the others (64.71 and 22.55% for the fuzzy and SHY95 methods, respectively). Table 5 shows that improvement in the fuzzy method is almost thrice that of the SHY95 method for convective rains classification. For the S-POL radar, the fuzzy method seems to identify better than the SHY95 and DSD retrieval methods, especially for the stratiform rains, with a missclassification rate of 5.87%. Hence, the fuzzy method can be considered suitable  Misc * represents misclassification. Error = 100−(accuracy rate of stratiform/convective rains)−Misc for classifying the rainfall types at the Jincheon site for all 12 cases that were selected. The fuzzy method performs the best classification for stratiform and convective rains in this study.

Conclusion
The aim of this study was to compare the existing classification methods by using Parsivel and S-POL data and to propose a new separation method for central South Korea after classification. Two identification methods of stratiform and convective rains were applied, namely BR03 and VPR, and a new unified method (BR03-VPR) was introduced for rainfall in Korea based on Parsivel disdrometer data. Generally, the BR09, YOU16, and TS96 classification methods could be applied in this study if the slope and/or intercept of the equations were altered. The newly developed separation lines performed well at the Jincheon site using DSD parameters, except for the relationship between log 10 N 0 and Λ. The best-fit line could not be obtained for the Jincheon site using the CA08 method, just as in the study of YOU16.
Three methods (SHY95, DSD retrieval, and fuzzy methods) were used to classify two rain types by the S-POL radar from the Yongin site. Briefly, SHY95 is based on Steiner et al. [2], wheares the DSD Figure 13. Same as Figure 12, but for intercept of Z H between 30 and 40 dBZ. retrieval method retrieves the N w and D 0 parameters according to BRA04 and BR09. The three different methods were compared for the Korean rainfall classification. For the stratiform rain classification, SHY95 was found to be the best method, with 0.00% misclassification, whereas the DSD retrieval method was the foremost method for convective rain classification, with 30.39% misclassification. The fuzzy method performs better than the SHY95 and DSD retrieval methods, with 5.87% and 34.31% misclassification for the stratiform and convective rains, respectively. Hence, it can be concluded that the new fuzzy method can classify stratiform and convective rains more accurately as compared with the DSD retrieval and SHY95 methods, even though only 12 cases were used in this study. The fuzzy method performs the best classification for stratiform and convective rains in this study.
Finally, it is suggested that the rainfall identification should be conducted using the new unified BR03-VPR method. This method is more accurate in classifying the rain types for both disdrometer and radar data. For the Parsivel disdrometer, the newly obtained classification line D 0 − log 10 N w from DSD data would be the best classification method for central South Korea. The fuzzy method was found to be more effective in classifying the rain types for S-POL radar from central South Korea. In future studies, the characteristics of DSDs for different rainfall systems and at different locations in South Korea will be examined to gain further understanding of the microphysical characteristics of the Korean rainfall system. This method can potentially be used to characterize rainfall type in any region of the world with a different coefficient, which depending on the different characteristics of the DSD. The accumulation of DSD data according to the rainfall types is very important for microphysical understanding and also radar calibration as well as Numerical Weather Forecast validation using radar data. Table A1. The vertex values of the trapezoidal MF for Z H (dBZ), Z DR (dB), K DP (deg km −1 ), and A H (dB km −1 ) as presented in Figure 12 Stratiform  Table A2. Same as A1, but for intercept of Z H between 30 and 40 dBZ as presented in Figure 13 Stratiform You et al. [6] Appendix A The four vertex values of trapezoidal MF (a, b, c, and d) were shown in Table A1 and A2 for overall and intercept of Z H between 30 and 40 dBZ, respectively as presented in Figure 12 and 13.