Abstract

In this study, to get an optimum rainfall estimation using polarimetric variables observed from Bislsan radar which is the first polarimetric radar in Korea, rainfall cases for 84 hours caused by different conditions, which are Changma front and typhoon, Changma front only, and typhoon only, occurred in 2011, were analyzed. And rainfall algorithms were developed by using long period drop size distributions with six different raindrop axis ratio relations. The combination of the relations between and , , and , , and and with different rainfall intensity would be an optimum rainfall algorithm if the reference of rainfall would be defined correctly. In the case the reference is not defined adequately, the relation between and , , , and and , , can be used as a representative rainfall relation. Particularly if the qualified is not available, the relation between and , , can be used as an optimum rainfall relation in Korea.

1. Introduction

Measurement of drop size distributions (DSDs) have been extensively used to calculate both radar reflectivity and rain rate for conventional radar and there is no unique relation between horizontal reflectivity () and rain rate () (hereinafter ) in the world because DSDs could vary from storm to storm and within the storm itself [1, 2]. Calculations of polarimetric parameters such as , differential reflectivity (), differential phase shift (), cross correlation coefficients (), specific differential phase (), and specific attenuation () could be obtained using T-matrix scattering techniques derived by Waterman [3] and later developed further by Mishchenko et al. [4]. The raindrop axis is one of the parameters for calculating polarimetric variables of T-matrix simulation. The variations of raindrop axis ratio in nature are intensively related to the oscillating of raindrop and hence are also connected to the polarimetric variables. There are many researches on the investigation of raindrop shape by laboratory study [5, 6], field measurements [7, 8], and modeling [9, 10].

Many researchers noticed that radar rainfall estimation is contaminated by a number of uncertainties such as hardware calibration, partial beam filling, rain attenuation, and nonweather echoes [11, 12]. To mitigate these problems, the particle identification algorithm using polarimetric parameters for improving data quality control and rainfall estimates by the discrimination of nonmeteorological artifacts such as anomalous propagation, birds, insects, and second trip echo was developed [1315]. And improvement of quantitative precipitation estimation (QPE) accuracy is one of the major points of polarization radar [1620]. Ryzhkov et al. [21] compared the rainfall relations with different drop shape assumptions and developed rainfall algorithm using polarimetric radar for the prototype WSR-88D. Cifelli et al. [22] compared the two rainfall algorithms, CSU-HIDRO (Colorado State University-Hydrometeor IDentification of Rainfall) and JPOLE- (Joint Polarization Experiment-) like, in the high plains environment. Recently, Ryzhkov et al. [23] investigated the potential utilization of for rainfall estimation using X-band and S-band radar. They found that method yields robust estimates of rain rates and rain totals even at S band which has very small attenuation.

There have been also many researches on polarimetric radar to implement it into operational uses. Based on these theoretical and other experimental researches, many countries are replacing or modifying their radars into polarimetric radar for operational use. There are three major agencies, Ministry of National Defense (MND), Ministry of Land, Infrastructure and Transportation (MOLIT), and Korea Meteorological Administration (KMA), which operate radars to monitor and forecast severe weather and flash flood operationally in Korea. Among these agencies, MOLIT installed polarimetric radars in 2009 and 2012 for the first time in Korea. KMA has installed the S-band polarimetric radar at the most northwestern part of Korea in 2014. For successful implementation of their radars for the purpose of operational uses, many researches on rainfall estimation, hydrometeor classification, and DSDs retrieval are required. However, there are few studies on these polarimetric related issues except for getting relationships using long period disdrometer data and assessment of each relation after applying a very simple quality control for differential phase shift [24]. They found that the accuracy of rainfall estimation using obtained by DSDs of Busan area in Korea was the best one comparing with relations calculated by ones of Oklahoma in US. And the quality control and unfolding of for calculating were applied to the rainfall estimation [25]. The above two studies used only 84,574 samples of DSDs excluding winter rainfall events, two and four raindrop shapes for calculation of rainfall relations.

This paper discussed how to improve the accuracy of the rainfall estimation using all polarimetric variables with different raindrop shapes and get optimum rainfall algorithm for Korean S-band polarimetric radar. In Section 2, rain gage, DSDs and radar dataset, calculation of rainfall relations, raindrop axis ratio relations, and the method of validation are described. Section 3 provides the optimum rainfall algorithms with and without rainfall category followed by the calculation of using observed and and the validations of rainfall estimation. Finally, we provide some conclusions and derived results are summarized in Section 4.

2. Data and Methodology

2.1. Rain Gage and Radar Dataset

The rainfall data from rain gages operated by the KMA were used to evaluate the accuracy of radar rainfall. Rain gages located from 5 km to 95 km within radar coverage are included in the analysis. Figure 1 shows the location of all instruments used in this study. The circle means the radar coverage, solid rectangle is the center of Bislsan radar, plus sign shows the distributed rain gages within radar coverage, and open rectangle is the location of a POSS (Precipitation Occurrence Sensor System) disdrometer that was installed around 82 km away from radar.

Radar data were collected by Bislsan polarimetric radar that was installed and operated by MOLIT in Korea since 2009. The specifications of Bislsan polarimetric radar and the quality control algorithm of were shown by You et al. [24, 25].

2.2. Calculations of Polarimetric Variables and Validation

The relations for converting radar variables into rain rate are required to get rainfall because radar could not observe the rainfall directly. In order to calculate these relations, disdrometer data which can measure the DSDs are needed. One-minute DSDs obtained by POSS (Precipitation Occurrence Sensor System) from 2001 to 2004 were used and processed to remove the unreliable data as shown by You et al. [25].

DSDs data for calculations of relationships were 114,105 samples after quality control and removal of negative . The total number concentration with respect to the drop size and the averaged parameters of gamma distributions was shown in Figure 2(a). The slope, shape, and intercept for gamma model were 8.2, 5.5, and , respectively, and median diameter was 1.3 mm. Most of the data are distributed in a wide range with a maximum rain rate of about 199.3 mm h−1 and the average was 2.47 mm h−1 (Figure 2(b)).

Polarimetric variables were calculated using T-matrix scattering techniques derived by Waterman [3] and later developed further by Mishchenko et al. [4]. To get the variables using DSDs, six raindrop shape assumptions are used.

The first raindrop axis ratio used in this study has slightly modified the relation proposed by Pruppacher and Beard [5] and will be called DS1: where , , and are the major axis, minor axis, and equivolume diameter of raindrop in mm, respectively:

Equation (2) is for equilibrium axis ratio derived from the numerical model of Beard and Chuang [9], which is in good agreement with the results from wind tunnel measurements (hereinafter DS2). The practical shapes of raindrops in turbulent flow are expected to be different shapes from the equilibrium shapes due to the drop oscillations. Oscillating drops appear to be more spherical on average than the drops with equilibrium shapes as shown by Andsager et al. [10] in laboratory studies. They figured out that the raindrops’ shape between 1.1 and 4.4 mm is better explained by the following formula:

Bringi et al. [26] suggested using (3) for drops with sizes smaller than 4.4 mm and (2) for larger sizes (hereinafter DS3). The shape-diameter relation that combines the observations of different authors was recently proposed by Brandes et al. [27, hereinafter DS4]:

The relations of raindrop axis ratio (DS5) proposed by Beard and Kubesh [28] and Thurai et al. [29] were combined, given by

The relation of raindrop axis ratio that slightly modified the relation proposed by Goddard et al. [30] was used (hereinafter DS6):

Another parameter in the T-matrix calculations is the temperature, which is assumed to be 20°C in this study. It is also necessary to take the canting angle into consideration of the T-matrix simulation because it can account for a 6% reduction in the coefficient of the relation [31] and may cause small negative biases of the estimators [32]. The distribution of canting angles of raindrops is Gaussian with a mean of and a standard deviation of , which have been recently determined by Huang et al. [33].

2.3. Validations

Because the rainfall in Korea is mostly accompanied with Changma front and Typhoon, three rainfall cases, which are caused by Changma front and typhoon, Changma front only, and typhoon only, were used for validations (Table 1).

Figure 3 shows the surface weather chart of each case. The typhoon MAERI was located at the eastern ocean of Taiwan and the Changma front was elongated from eastern China continent to the central Japan through the southern part of Korea on 0000 LST June 25 (Figure 3(a)). The MAERI moved to the north, located at the southern west sea of Korea, and made rainfall in the Korean peninsula on 0900 LST June 26 in 2011 (Figure 3(b)). Changma front was located at the southern part of Korea and brought rainfall at the analyzed area on 1200 LST July 9 in 2011 (Figure 3(c)). The rainfall was affected by Changma front all the time during case 2. The typhoon MUIFA was located at the southwestern sea of Korea on 0000 LST August 8 in 2011 and caused rainfall at the target area (Figure 3(d)).

Figure 4 shows the time series of averaged rainfall amount observed on the ground rain gages within radar coverage from 5 km to 95 km. Averaged rainfall amount refers to that obtained by averaging the amount of rainfall observed by rain gages within the radius of the radar. There are three peaks of rainfall in case of Changma front and typhoon; the first two peaks were due to Changma front and the third one was due to the influence of the typhoon. There are three peaks of rainfall accompanied with Changma front in the second rainfall event. The third event was also caused by typhoon but was relatively short. The period of the selected rainfall was 84 hours, 29 hours for Changma front and typhoon, 46 hours for Changma front only, and 9 hours for typhoon only.

The normalized error (NE), fractional root mean square error (RMSE), and correlation coefficients (CC) of rainfall relations and 121 gages were used to investigate the performance of each rainfall relation:where is the number of the RR and RG pairs and and are the averaged rain rate of radar and gage for an hour, respectively. The above statistical variables are calculated using 1-hour rainfall amount of radar and gage at the point. The point rainfall of radar was obtained by averaging rainfall over a small area (1 km × 1°) centered on each rain gage.

3. Results

3.1. Rainfall Relations with Different Raindrop Axis Ratios
3.1.1. The Characteristics of , , and with Different Drop Shapes

The histograms of occurrence frequency for polarimetric variables, , , and , calculated by DSDs data for 4 years in Busan with different raindrop axis ratio relations, DS1, DS3, and DS6, were shown in Figure 5.

The averages and modes of were around 26.2 dBZ and 32~33 dBZ for all raindrop axis assumptions (Figures 5(a)~5(c)). It means that the reflectivity is not sensitive to the drop axis relation. In case of , there were two modes of occurrence for all cases. The averages and standard deviation of each relation were 0.63 dB and 0.91 dB for DS1, 0.61 dB and 0.8 dB for DS3, and 0.66 dB and 0.82 dB for DS6. The low and high modes of with DS1, DS3, and DS6 were 0.2~0.3 dB and 1.7~1.8 dB, 0.2~0.3 dB and 1.6~1.7 dB, and 0.0~0.1 dB and 1.7~1.8 dB, respectively. The occurrence frequencies of low mode for each raindrop axis ratio were significantly different from each other (Figures 5(d)~5(f)). The averages and standard deviation of with different raindrop shapes were 0.11 and 0.04, 0.09 and 0.03, and 0.1 and 0.03, respectively. The modes of all drop shapes were the same, but the occurrence frequencies were different (Figures 5(g)~5(i)).

3.1.2. The Statistics of Rainfall Relations with Different Raindrop Axis Ratios

Because the occurrence frequencies of and with different raindrop axis ratios were different from each other, the rainfall relations using those variables should be different with drop shape assumptions. Table 2 shows the rainfall relations obtained by using different raindrop shape assumptions. The coefficients of were not significantly different with drop shape assumptions; however those of other relations were different with each drop shape.

Table 3 shows the cross correlations (hereinafter CC) and RMSEs (root mean square errors) of rainfall relations, , , , , and , obtained by calculations using DSDs data with different raindrop shapes. The statistics of and were not significantly different with raindrop shapes. The CC and RMSE of and were the worst and the best among the other rainfall relations. The statistics of , , and were much more variable with different raindrop axis ratios than the ones of and . The RMSEs of , , and with raindrop shapes were distributed from 3.030 to 3.828 mm, 2.965 to 3.523 mm, and 3.151 to 5.412 mm, respectively. The best performance of each relation occurred at DS3 for , DS1 for , DS1 for , DS1 for , and DS3 for .

In order to calculate more accurate and , the and with the best performance were chosen. Figure 6 shows the scatter plots of rainfall obtained by DSDs and and using the best statistics among raindrop axis ratio. , , and were chosen from DS3, DS1, and DS1, respectively. Comparing with and of single raindrop axis ratio relations having the best performance, new combined relations had better RMSE and CC. The RMSEs of new relations, and , had better score as much as around 0.2 mm and 0.6 mm, respectively. Even though the CC and RMSE of and with different drop shapes were not significant, the combined relations had better performance.

3.2. Validations of Rainfall Relations with Different Raindrop Axis Ratio Relations

To investigate the performance of rainfall relations, NE (normalized error), RMSE, and CC were calculated using rainfall from each relation with six raindrop shapes and that of gage rainfall.

Table 4 summarizes the relations and the statistics such as NE, RMSE, and CC. The NEs and RMSEs of calculated by each raindrop axis were distributed from 0.52 to 0.55 and from 4.625 to 4.996, respectively. The with assumption of DS5 was the best score of RMSE in other raindrop shapes. In case of , the distribution of NEs and RMSEs was from 0.31 to 0.43 and from 3.793 to 4.602, respectively. The best RMSE score of was from DS2. The NEs and RMSEs of occurred from 0.55 to 0.65 and from 5.146 to 7.141, respectively. The performance of was the worst score and had the best score in all raindrop axis ratio relations. The performances of validation were different from that of rainfall relation calculation. It would be caused by the variations of DSDs in this study.

To compare the performance between new combined and , the statistics were also calculated. Figure 7 shows the scatter plots rainfall from rainfall relation and gage rainfall with some statistics. The NE and RMSE of two relations from single raindrop shape assumption showed better results. However, it seems that the with two-raindrop axis ratio was more close to the gage rainfall in the range of weaker than 20 mm/h and the with two drop shapes was more accurate in the rainfall of higher than 20 mm/h.

3.3. Rainfall Estimation Using Specific Attenuation
3.3.1. Calculation of Specific Attenuation

The can be calculated from the radial profile of the attenuated reflectivity and the two-way PIA (Path Integrated Attenuation) along the propagation path proposed by Meneghini and Nakamura [34]:whereIf is not dependent on range, then (8) becomeswhere

Bringi et al. [35] recommended estimating PIA using byand Testud et al. [36] used (10) and (12) to obtain radial profiles of at C-band. In this study, was calculated by the method proposed by Ryzhkov et al. [23]. The constant was set by 0.6 and was by 0.027 calculated by the ratio of to obtained from DSDs.

Figure 8 shows the scatter plot of rainfall from , , and and rainfall from DSDs and an PPIs (Plan Position Indicators) at 0.5 degree elevation angle of gage rainfall and rainfall from and at 0251 KST on the 8th of August in 2011.

The relation had much better fit to the rainfall of DSDs than that of and relation. Comparing with the distribution of rainfall obtained by and , has better spatial resolution and more homogeneous pattern than those of .

3.3.2. Validations of , , and Relations

As mentioned in the previous section, has the potential to be the best choice for estimating rainfall using polarimetric variables. To verify the accuracy, , , and relations were obtained and the accuracy of , , and , which had the best performance in the relations calculated by the six raindrop axis ratio relations, were examined by comparing with .

Figure 9 shows the scatter plot of rainfall from gage rainfall and , , , and for 84 hours. The accuracy of was the best and was the second in 4 relations. was not better than in whole range of rainfall. However, in the range of higher rainfall, it seems that was closer to the gage rainfall. The error of would be caused by the missing radial profile of along the ray. The missing has occurred, if the difference of between the starting and end gate is negative. According to the results, and would compensate for this region.

3.4. Optimum Rainfall Algorithm

To find out the optimum rainfall algorithm for Korean S-band polarimetric radar, the gage rainfall amount was categorized by three steps, 0 to 5 mm/h assigned as light rain, 5 to 30 mm/h as medium rain, and higher than 30 mm/h as high rain. In previous section, the relations have different accuracy not only for the raindrop axis ratio relation but also for rainfall amount. The sample numbers for each category were 3,322 samples, 1,980 samples, and 92 samples, respectively.

Figure 10 shows the RMSEs of , , , , and with different rainfall categories defined in three steps. All relations have different RMSEs with respect to the raindrop axis ratio relations. with DS3, with DS2, and with DS3 have the best score at the high rainfall regime, at the low rainfall regime, and at the medium rainfall regime, respectively. Even though has the best performance among other relations in total rainfall events, using different rainfall relations with different rainfall regime would be an optimum rainfall algorithm for Korean S-band polarimetric radar. A possible optimum polarimetric rainfall algorithm can be expressed by

Table 5 shows the rainfall relations of and with different raindrop shape assumptions. The coefficients of each relation with respect to the drop shapes were much different comparing with relations obtained by combining two polarimetric variables.

In case of rainfall relations combined with , has better score than . Both and with DS3 showed better result than other raindrop shape assumptions. at the low, medium, and high rainfall regimes has similar RMSE to , , and , respectively.

Equations (14) would be an optimum rainfall algorithm, but there is still a problem to solve how the rainfall categories are defined using radar and gage. Gage does not have good spatial resolution to cover radar resolution like 125 m or 1 km and each rainfall relation has its own error at certain rainfall regime. Therefore, and with DS3 can be used regardless rainfall intensity. And in case bias calibration is not available or does not have enough quality for quantitative use, with DS3 can be used as a representative rainfall estimation.

4. Summary and Conclusions

Polarimetric radars will be main tools to monitor and forecast severe weather and flash flooding within several years in Korea. To get an optimum rainfall algorithm using polarimetric variables observed from Bislsan radar which is the first polarimetric radar in Korea, rainfall cases for 84 hours caused by different conditions, which are Changma front and typhoon, Changma front only, and typhoon only, that occurred in 2011 were analyzed. And rainfall relations were obtained by using long period DSDs with six different raindrop axis ratio relations.

In the analysis of , , and occurrence frequency, there were two modes of occurrence frequency and the variations of and histograms were higher than that of with raindrop ratio relations. According to these variations, the combined relations of and using with DS3, with DS1, and with DS1 were closer to the rainfall of DSDs.

To examine the performance of each rainfall relation, NE, RMSE, CC were calculated using rainfall recorded at 121 gages within radar coverage for 84 hours. The statistics of each rainfall relation were different with raindrop shape assumptions and rainfall intensity. had the worst performance and had the best score in all raindrop shapes. was calculated by observed and using the method by Ryzhkov et al. [23] and the rainfall relations using were also calculated and analyzed.

The performance of rainfall relations were compared with three different rainfall categories to find out an optimum rainfall relation for the S-band polarimetric in Korea. , , and had the best RMSE at the light rain, medium rain, and high rainfall regimes, respectively. and showed relatively good performance in all rainfall regimes. The combination of , , and with rainfall intensity would be an optimum rainfall algorithm if the reference of rainfall would be defined correctly. Regardless of rainfall intensity, and obtained by assuming DS3 can be used as a representative rainfall relation without consideration of rainfall intensity regime. Particularly if the qualified is not available, with DS3 drop shape assumption can be used as an optimum rainfall relation in Korea.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors acknowledge providing radar data, weather chart, and AWS data for this work from the Ministry of Land, Infrastructure, Transport and Korea Meteorological Administration. The authors also acknowledge providing codes for scattering simulation from Professor V. N. Bringi at Colorado State University. This work was funded by the Korea Meteorological Industry Promotion Agency under Grant KMIPA 2015-1050.