Dual-polarization radar rainfall estimation in Korea according to raindrop shapes obtained by using a 2-D video disdrometer

Abstract

Polarimetric measurements are sensitive to the sizes, concentrations, orientations, and shapes of raindrops. Thus, rainfall rates calculated from polarimetric radar are influenced by the raindrop shapes and canting. The mean raindrop shape can be obtained from long-term raindrop size distribution (DSD) observations, and the shapes of raindrops can play an important role in polarimetric rainfall algorithms based on differential reflectivity (Z_DR) and specific differential phase (K_DP). However, the mean raindrop shape is associated with the variation of the DSD, which can change depending on precipitation types and climatic regimes. Furthermore, these relationships have not been studied extensively on the Korean Peninsula. In this study, we present a method to find optimal polarimetric rainfall algorithms for the Korean Peninsula by using data provided by both a two-dimensional video disdrometer (2DVD) and the Bislsan S-band dual-polarization radar. First, a new axis-ratio relation was developed to improve radar rainfall estimations. Second, polarimetric rainfall algorithms were derived by using different axis-ratio relations. The rain gauge data were used to represent the ground truth situation, and the estimated radar-point hourly mean rain rates obtained from the different polarimetric rainfall algorithms were compared with the hourly rain rates measured by a rain gauge. The daily calibration biases of horizontal reflectivity (Z_H) and differential reflectivity (Z_DR) were calculated by comparing Z_H and Z_DR radar measurements with the same parameters simulated by the 2DVD. Overall, the derived new axis ratio was similar to the existing axis ratio except for both small particles (≤2 mm) and large particles (≥5.5 mm). The shapes of raindrops obtained by the new axis-ratio relation carried out with the 2DVD were more oblate than the shapes obtained by the existing relations. The combined polarimetric rainfall relations using Z_DR and K_DP were more efficient than the single-parameter rainfall relation for estimated 2DVD rainfall; however, the R(Z_H, Z_DR) algorithm showed the best performance for radar rainfall estimations, because the rainfall events used in the analysis consisted mainly of weak precipitation and K_DP is relatively noisy at lower rain rates (≤10 mm h¹). Some of the polarimetric rainfall algorithms can be further improved by new axis-ratio relations.

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/ doi:10.5194/amt-9-3863-2016 Author(s) 2016. CC Attribution 3.0 License.

Hae-Lim Kim1, Mi-Kyung Suk1, Hye-Sook Park2, Gyu-Won Lee3, and Jeong-Seok Ko1

1Radar Analysis Division, Weather Radar Center, Korea Meteorological Administration, Seoul, South Korea

2Satellite Analysis Division, National Meteorological Satellite Center, Korea Meteorological Administration, Seoul, South Korea

3Department of Atmospheric Sciences, Kyungpook National University, Daegu, South Korea

Correspondence to: Hae-Lim Kim ([email protected])

Received: 18 January 2016 Published in Atmos. Meas. Tech. Discuss.: 27 January 2016 Revised: 25 July 2016 Accepted: 25 July 2016 Published: 22 August 2016

Abstract. Polarimetric measurements are sensitive to the sizes, concentrations, orientations, and shapes of raindrops. Thus, rainfall rates calculated from polarimetric radar are inuenced by the raindrop shapes and canting. The mean raindrop shape can be obtained from long-term raindrop size distribution (DSD) observations, and the shapes of raindrops can play an important role in polarimetric rainfall algorithms based on differential reectivity (ZDR) and specic differential phase (KDP). However, the mean raindrop shape is associated with the variation of the DSD, which can change depending on precipitation types and climatic regimes. Furthermore, these relationships have not been studied extensively on the Korean Peninsula. In this study, we present a method to nd optimal polarimetric rainfall algorithms for the Korean Peninsula by using data provided by both a two-dimensional video disdrometer (2DVD) and the Bislsan S-band dual-polarization radar. First, a new axis-ratio relation was developed to improve radar rainfall estimations. Second, polarimetric rainfall algorithms were derived by using different axis-ratio relations. The rain gauge data were used to represent the ground truth situation, and the estimated radar-point hourly mean rain rates obtained from the different polarimetric rainfall algorithms were compared with the hourly rain rates measured by a rain gauge. The daily calibration biases of horizontal reectivity (ZH) and differential reectivity (ZDR) were calculated by comparing ZH and ZDR radar measurements with the same parameters simulated by the 2DVD. Overall, the derived new axis ratio was similar to the existing axis ratio except for both small particles ( 2 mm) and large particles ( 5.5 mm). The shapes of

Dual-polarization radar rainfall estimation in Korea according to raindrop shapes obtained by using a 2-D video disdrometer

raindrops obtained by the new axis-ratio relation carried out with the 2DVD were more oblate than the shapes obtained by the existing relations. The combined polarimetric rainfall relations using ZDR and KDP were more efcient than the single-parameter rainfall relation for estimated 2DVD rainfall; however, the R(ZH, ZDR) algorithm showed the best performance for radar rainfall estimations, because the rainfall events used in the analysis consisted mainly of weak precipitation and KDP is relatively noisy at lower rain rates ( 10 mm h1). Some of the polarimetric rainfall algorithms

can be further improved by new axis-ratio relations.

1 Introduction

Radar is a very useful monitoring tool for extreme weather forecasting, ood forecasting, and rainfall estimations because of its high spatial and temporal resolution. In particular, dual-polarization radar that provides information on the reectivity (ZH), differential reectivity (ZDR), differential phase ([Phi1]DP), specic differential phase (KDP), and cross-correlation coefcient (hv) can distinguish precipitation types by means of backscatter and the differential propagation phase of hydrometeors. Dual-polarization radar can also be used to obtain more information about the raindrop size distribution (DSD), and this in turn can help to reduce the impact of DSD variability on rainfall estimations (Cifelli et al., 2011). The DSD can be characterized by parameters such as the diameter, concentration, orientation, and shape. Disdrometers and aircraft imaging investigations have been

Published by Copernicus Publications on behalf of the European Geosciences Union.

3864 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

widely used to study the DSD variability especially as a function of the rainfall rate (Bringi et al., 2003). Polarimetric parameters are sensitive to the DSD properties. Rainfall rates estimated from polarimetric radar measurements are affected by the mean shape of raindrops and canting (Brandes et al., 2002).

The average shape of a raindrop can be inferred from its size through application of a shapesize relationship for a raindrop. Some researchers have attempted to produce the mean shape of raindrops. Keenan et al. (2001); Brandes et al. (2002); Thurai and Bringi (2005); and Marzuki et al. (2013) derived empirical relations from observational data, and Pruppacher and Beard (1970); Green (1975), and Beard and Chuang (1987) investigated the shape of raindrops falling under the inuence of gravity. These raindrop axis-ratio relations play an important role in the derivation of polarimetric radar rainfall algorithms that use ZDR and KDP (Jameson, 1983, 1985; Gorgucci et al., 2001). Polarimetric rainfall algorithm can be derived from long-term DSD measurements. However, the derivation of new empirical algorithms based on long-term DSD data is still challenging for different climatological environments.

Several different polarimetric rainfall algorithms have been developed by assuming raindrop shapes (Sachidananda and Zrni, 1987; Chandrasekar et al., 1990; Ryzhkov and Zrni, 1995; Gorgucci et al., 2001). However, radar rainfall estimations are affected by the following two different sources of errors: (i) errors in radar measurements such as from attenuation, bright bands, ground clutter, and calibration bias of ZH and ZDR (Zawadzki, 1984; Villarini and Krajewski, 2010; Sebastianelli et al., 2013) and (ii) errors in the conversion of the radar measurements into rainfall rates at the ground level. In addition, the error is also related to the gamma assumption on the retrieve of the rain rate from dual-polarization radar measurements (Adirosi et al., 2014). The DSDs variability affects the accuracy of radar-derived rainfall (Maki et al., 2005). Many earlier studies have explained the variability in the DSD by evoking both different storm types and climatic regimes, and variations of the DSD are known to be caused by climatological and physical factors (e.g., Ulbrich, 1983; Tokay and Short, 1996; Bringi et al., 2003). Variations in time and space of the DSD are associated with changes in microphysical processes such as evaporation, break-up, coalescence, and condensation. They can be also due to the variation of vertical air motion (Marzuki et al., 2013).

Radar measurements often suffer from the system biases of ZH and ZDR, and thus, accurate measurement and calibration of ZH and ZDR values are necessary to achieve accurate radar rainfall estimations (Park and Lee, 2010). An assessment of the calibration bias of polarimetric radar is possible by monitoring the hardware stability, and ZH and ZDR measurements can be corrected by using ground validation equipment such as a disdrometer. Joss et al. (1968) calibrated radar reectivity by using the vertical prole of reectivity

and disdrometer-inferred ZH. The radar reectivity was calibrated by comparisons between the radar and disdrometer reectivity to check the calibration of the WSR-88D (weather surveillance radars 1988 Doppler) at Greer, South Carolina (Ulbrich and Lee, 1999). Goddard et al. (1982) and Goddard and Cherry (1984) compared radar results with disdrometer results by using the axis-ratio relations of Pruppacher and Beard (1970) and Pruppacher and Pitter (1971), respectively; they found that radar measures of ZDR were 0.3 and 0.1 dB lower than the disdrometer estimates. In addition to disdrometers, there are various other ways to correct biases in radar data, such as by using the ZHKDP implementation of the self-consistency principle, vertically pointing measurements, and comparisons of measured data and the mean ZHZDR relationship (Kwon et al., 2015). Moreover, a variety of radar calibration methods were introduced in Atlas (2002).

In this study, we computed the mean axis-ratio relation and developed several polarimetric rainfall algorithms by using two-dimensional video disdrometer (2DVD) measurements collected from September 2011 to October 2012 in Daegu, Korea. Specically, three raindrop shape assumptions (Pruppacher and Beard, 1970; Beard and Chuang, 1987; Brandes et al., 2002) and the newly derived axis-ratio relation from 2DVD data were used to derive polarimetric rainfall retrieval algorithms. The ZH and ZDR biases of Bislsan dual-polarization radar were calibrated by comparing them with simulated ZH and ZDR values obtained from the 2DVD.

Improvements in the quantitative rainfall estimations were achieved by applying the ZH and ZDR calibration biases.

In Sect. 2, the data used in this study are described. The methodology for 2DVD data quality control, derivation of the raindrop-axis ratio from 2DVD data, and simulations of the polarimetric parameters by the T-matrix scattering method are described in Sect. 3. The results of the statistical validation of rainfall estimation are presented in Sect. 4. Finally, the conclusions are given in Sect. 5.

2 Data and instruments

2.1 Disdrometer

In this study we used the compact 2DVD version deployed on the campus of Kyungpook National University, Daegu, Korea (35.9 N, 128.5 E). The 2DVD was installed in an observation eld away from the buildings. Disdrometers are used to nd the DSD characteristics at a given location. Disdrometer data used in this study were collected by a 2DVD from September 2011 to October 2012.

The 2DVD consists of two orthogonal light sheets (referred to as A and B line-scan cameras). Line-scan cameras have single-line photo detectors. The particle shadows are detected on the photo detectors, and the particle images are recorded from two sides and at different heights when the particles are falling through the measurement area (10 cm 10 cm). The 2DVD measures drop size, fall veloc-

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3865

Table 1. Specications for dual-polarization radar in Bislsan.

Parameters Characteristics

Variables ZH, Vr, SW, ZDR, [Phi1]DP, KDP, hv

Altitude of radar antenna 1085 m

Transmitter type Klystron

Transmitter peak power 750 kWAntenna diameter 8.5 mBeam width of radar 0.95

Observation

frequency 2.785 MHz (S-band)range 150 kmgate size 125 melevations 0.5 , 0 , 0.5 , 0.8 , 1.2 , 1.6 (6 elevation)

ity, and the shape of individual particles. From these data, one can calculate the DSD and all related quantities such as the rain rate, total drop number concentration, and liquid water content. A more detailed description of the 2DVD is given in Kruger and Krajewski (2002).

2.2 Radar

The Ministry of Land, Infrastructure, and Transport (MLIT) operates the Bislsan (BSL) dual-polarization radar in Bislsan, Korea (35.7 N, 128.5 E, 1085 m). The BSL S-band radar has observation range of 150 km and a frequency of2.5 min because it is used primarily for hydrological observations and ood forecasts. The BSL S-band radar measures polarimetric variables such as ZH, ZDR, KDP, and hv in real time. Six elevation angles from 0.5 to 1.6 were available,

with a gate size resolution of 125 m and radar beam width of 0.95 . The specications of the BSL radar are shown in

Table 1.The BSL S-band dual-polarization radar was located about22.3 km (17 ) away from the 2DVD location (Fig. 1). These geographical locations were adopted to compare the two sets of observation data. We used radar data from September 2011 to October 2012. During this period, rainfall events were analyzed for ZH and ZDR calibration and for rainfall estimation.

In addition, the 0.0 -elevation plan position indicator (PPI) radar data were used to avoid effects due to beam blocking and ground echoes on the measurements, whilst being as close as possible to the 2DVD. The ZH and ZDR radar parameters were averaged over ve successive gate size resolutions and two adjacent azimuth angles, and KDP was calculated from the ltered [Phi1]DP as the slope of a least squares t.

2.3 Rain gauge

A tipping bucket rain gauge was used to validate the 2DVD rainfall estimations. The rain gauge used in this study was a RG3-M tipping bucket rain gauge from the Onset Computer Corporation. The maximum rainfall rate of the rain gauge was 127 mm (5 in) per hour, and the operating tem-

Figure 1. The locations of the Bislsan polarimetric radar and the 2DVD at the rain gauge site.

perature range was from 0 to 50 C. The bucket size of the rain gauge was 0.2 mm and the time resolution was 0.5 s. The rain gauge measurements were corrected to reduce instrumental uncertainties through eld inter-comparisons with a reference gauge. The rain gauge was installed in the same location as the 2DVD.

3 Methodology

3.1 Quality control of 2DVD data

The 2DVD observation data were useful for investigating the characteristics of rainfall. However, a number of particle outliers were measured, and these anomalous data points were due to wind turbulence, splashing, break up of drops, and

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3866 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

mismatching between camera A and B (Raupach and Berne, 2015). These results can lead to incorrect information about the particles. Therefore, before using the 2DVD data, a quality control process was needed.

Figure 2a and b show the fall velocity and oblateness distribution according to the raindrop diameter before the data quality control procedure was performed. In Fig. 2b and c, we compare the axis-ratiodiameter relation of Pruppacher and Beard (1970) with that found by the disdrometer before and after the correction, respectively. Some particles had fall velocities beyond the terminal velocity of large raindrops.In particular, the outliers appeared prominently in the small raindrop ranges. To remove these outliers, velocity-based ltering was applied to the 2DVD measurement data (Thurai and Bringi, 2005). The equations used are the following:

|Vmeasured (D) VA (D)| < 0.4VA (D)

VA (D) = 9.65 10.3exp(0.6D), (1)

where D (mm) is the drop diameter, Vmeasured (m s1) is the fall velocity as measured by the 2DVD, and VA represents the Atlas velocity formula (Atlas et al., 1973). The procedure removed about 17 % of the values that did not correspond well with the expected normal distribution (Fig. 2c). Despite the application of velocity-based ltering, signicant bias still remained in the small drop size area. This was due to instrumental limitations such as the mismatch problems with the line-scan cameras and the limited vertical resolution of the instrument. Therefore, the oblateness data corresponding to raindrop diameters smaller than 0.5 mm were removed when we calculated the new axis-ratio formula.

To analyze the reliability of the 2DVD data, we compared the rain rates calculated from the 2DVD data (Eq. 2) to col-located rain gauge measurements; the difference of accumulated rainfall represents the percent error (Eq. 3).

R = 6 104XDminDmaxD3V (D)N (D)[Delta1]D [mm h1] (2)

ARraingauge AR2DVD|ARraingauge 100[%] (3)

Here, Dmax and Dmin are the maximum and minimum diameters of the observed drops in mm, N(D) is the drop number concentration in mm1 m3, V(D) is the drop fall velocity in m s1, and D is the drop interval (D = 0.2 mm). The drop fall

velocity formula used was that derived by Atlas et al. (1973); PE is the percent error and AR is the accumulated rainfall in mm.

We analyzed the rainfall cases that occurred from September 2011 to October 2012. Figure 3 shows the accumulated rainfall computed from the 2DVD and rain gauge data for six of these cases. The six events occurred at (i) 00:00 09:00 UTC on 14 October 2011 (Fig. 3a), (ii) 14:00 23:59 UTC on 2 April 2012 (Fig. 3b), (iii) 00:0023:59 UTC on 21 April 2012 (Fig. 3c), (iv) 00:0008:00 UTC on 25

Figure 2. Distribution of fall velocity and oblateness according to drop diameter. (a) Velocity-based lter for the drop measurements. The color scale represents the drop number density (log scale). The dotted line (orange) represents the results from the velocity formula of Atlas et al. (1973). (b) Drop axis ratios for all measured drops.(c) Drop axis ratios after removing mismatched drops. The dotted line (black) represents the mean axis-ratio results from the formula given in Pruppacher and Beard (1970).

PE = |

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3867

Figure 3. Time series of accumulated rainfall measured from the rain gauge and estimated from the 2DVD: (a) 14 October 2011, (b) 2 April 2012, (c) 21 April 2012, (d) 25 April 2012. (e) 23 August 2012, and (f) 27 August 2012.

April 2012 (Fig. 3d), (v) 00:0023:59 UTC on 23 August 2012 (Fig. 3e), and (vi) 16:0023:59 UTC on 27 August 2012 (Fig. 3f). Figure 3a shows the accumulated rainfall computed from the 2DVD and rain gauge data on 14 October 2011. As shown in Figs. 3af, the percent errors from the six considered rainfall cases were 2.40, 9.78, 0.12, 4.48, 8.38 and 8.66 %, respectively. The overall distribution between the 2DVD and rain gauge results was good. In general, earlier studies have found that rainfall differences between disdrometer and rain gauge data were mostly in the range of 10 to 20 % (McFarquhar and List, 1993; Sheppard and Joe,

1994; Hagen and Yuter, 2003; Tokay et al., 2003). These differences might have been due to such issues as differences in instruments, effects due to the measurement environment, and rainfall variability. Therefore, the rainfall differences between the 2DVD and rain gauge results used in this study were limited to a maximum of 20 % error, and the 2DVD data were excluded from the analysis when the rainfall difference between the 2DVD and rain gauge results exceeded 20 %.

After the quality control process, 33 rainfall cases were selected for further investigations of the characteristics of

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3868 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Table 2. Summary of the date, type of precipitation, and accumulated rainfall comparison values between 2DVD and rain gauge measurements for the 33 rainfall events.

Date Time of Accumulated rainfall (mm) PE (%) of observation (UTC) Type 2DVD rain gauge

09/05/11 01:0015:00 C 7.07 6.86 3.16 09/09/11 15:2021:20 S 4.28 4.31 0.72 09/10/11 00:0023:59 C 19.48 17.83 9.28 09/29/11 00:0017:00 S 3.81 3.72 2.46 10/13/11 00:0023:59 S 3.93 4.11 4.59 10/14/11 00:0009:00 S 13.20 13.52 2.40 10/21/11 06:0023:59 S 53.33 65.83 18.99 04/02/12 14:0023:59 M 17.85 16.26 9.78 04/21/12 00:0023:59 S 32.17 32.13 0.12 04/25/12 00:0008:00 S 18.63 17.83 4.48 05/01/12 08:0021:00 S 4.08 3.53 15.62 05/08/12 07:0011:00 C 8.94 9.60 6.89 05/14/12 00:0015:00 S 22.16 19.79 10.88 05/28/12 06:0007:00 C 15.19 14.50 4.79 06/08/12 03:0023:59 C 17.07 19.40 11.98 06/23/12 00:0008:00 C 13.43 14.50 4.79 07/06/12 00:0018:00 C 22.43 20.18 11.13 07/12/12 17:0021:00 M 7.51 7.45 0.88 07/13/12 01:0012:00 C 25.22 25.08 0.55 07/15/12 00:0012:00 C 8.22 7.64 7.54 07/16/12 15:0023:59 M 16.40 18.42 10.97 07/21/12 09:0010:30 C 5.86 5.68 3.15 08/12/12 01:0018:00 C 20.25 18.42 9.95 08/13/12 00:0015:00 C 37.27 34.09 9.31 08/23/12 00:0023:59 M 90.88 83.86 8.38 08/24/12 00:0023:59 S 7.37 7.45 1.01 08/27/12 16:0023:59 S 13.62 12.54 8.66 08/29/12 19:0023:59 S 6.26 5.29 18.39 09/09/12 00:0023:59 S 21.56 18.65 15.61 09/16/12 00:0023:59 S 83.12 69.88 18.95 09/17/12 00:0007:00 M 63.64 58.55 8.69 10/22/12 06:0011:00 M 13.63 14.92 8.64 10/27/12 00:0009:00 S 8.77 8.20 6.94

rainfall over the Korean Peninsula. The accuracy of the 33 rainfall cases was in the range of 0.1218.99 % compared to in situ rain gauge data. The radar reectivity and 1 h rainfall rates measured by the rain gauge were used to classify the data into different precipitation types. The criterion for the reectivity described by Chang et al. (2009) was applied in the present study, and rainfall rates that had values of R > 5 mm h1 ( 5 mm h1) were considered to be

of the convective (stratiform) type. The data set consisted of 15 stratiform rainfall cases, 12 convective rainfall cases, and 6 mixed rainfall cases (total of 33 rainfall events) with 17 618 min DSD samples. The precipitation type, rainfall difference, and accumulated rainfall for the 2DVD and rain gauge results are listed in Table 2. Figure 4 shows the hourly and total accumulated rainfall for the 2DVD and rain gauge results. The overall agreement between the 2DVD and the rain gauge results was good, and the total accumulated rain-

fall recorded by the rain gauge was greater than that of the 2DVD by about 0.30 %.

3.2 Raindrop axis ratio

A very small raindrop has an approximately spherical shape that becomes oblate as its size increases. The shape of a raindrop according to the drop size can be expressed as the mean axis-ratio relation. The mean raindrop shape is associated with the measured DSD shape and diameter, which is related to the variation of DSD. The DSD variations depend on different storm types and climatic regimes (Marzuki et al, 2013), and they affects rainfall rates derived from polari-metric radar measurements of reectivity. Hence, in order to produce rainfall estimation algorithms reective of the rainfall characteristic of the Korean Peninsula, a new mean axis-ratio relation, using the 2DVD data listed in Table 2, was de-

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3869

Figure 4. One-hour (left panel) and total accumulated rainfall (right panel) from the 2DVD and rain gauge for the 33 rainfall cases.

rived as a polynomial function. The size of the diameter bin was 0.2 mm, and the oblateness data corresponding to raindrop diameters smaller than 0.5 mm were removed when we derived the new axis-ratio relation because oblateness in the small range (D < 0.5 mm) was signicantly inuenced by the calibration of the instrument (Marzuki et al., 2013). In addition, we only displayed the average values for the bins containing more than ve drops in order to have sufcient numbers of raindrops. Therefore, the measured maximum diameter could reach to about 8 mm; however, the mean axis-ratio tting was established to be within 7 mm diameter because it does not have a sufcient number of values in more than 7 mm diameter.

In order to produce the mean axis-ratio relation, various tting methods such as linear and polynomial (second-, third-, fourth-order) ts were tried. The third-order polynomial relation was deemed the most suitable for the observation data, as this, relation performed (i.e., goodness of the tting) better than others. The third-order polynomial new mean axis-ratio relation (b/a) is as shown in Eq. (4), which can be reasonably extended to 7 mm:

b/a = 0.997845 0.0208475D 0.0101085D2

+6.4332 104D3 (0.5 D 7mm), (4)

where a and b are the major and minor axis, respectively, and D is the equivalent volume diameter of the particle in mm. The statistical measures corresponding to the correlation coefcient, root-mean-square error (RMSE), and mean absolute error (MAE) were 0.989, 0.025, and 0.012, respectively.

3.3 Disdrometer-rainfall algorithms

In order to produce the polarimetric rainfall algorithms, the theoretical polarimetric variables (e.g., ZH, ZDR, and KDP)

were simulated from the 2DVD data by using the T-matrix (transition) method (Zhang et al., 2001). The polarimetric variables depend on the shape of raindrops. The polarimetric rainfall algorithms were derived by making assumptions

about different axis-ratio relations. First, we calculated the complex scattering amplitudes of raindrops at the S-band of the 10.7 cm wavelength by using the mean axis-ratio relations. Second, we calculated the scattering amplitudes for the different axis-ratio relations that were used for the production of the polarimetric variables. The dual-polarimetric variables were calculated by using Eqs. (57) (Jung et al., 2010):

Zh =

4 4 4|Kw|2 [integraldisplay]0

[vextendsingle][vextendsingle][vextendsingle]

fa ()|2 + B[vextendsingle][vextendsingle][vextendsingle]

Dmax, x

A

fb ()|2

+2C Re

fa ()fb ()

[bracketrightbig]

N (D)dD [mm6 m3], (5)

Zv =

4 4 4|Kw|2 [integraldisplay]0

Dmax, x

B

[vextendsingle][vextendsingle][vextendsingle]

fa ()|2 + A[vextendsingle][vextendsingle][vextendsingle]

fb ()|2

+2C Re

fa ()fb ()

[bracketrightbig]

N (D)dD [mm6 m3], (6)

where

A =< cos4[Phi1] >=

18(3 + 4cos2

e22 + cos4

e82)

B =< sin4[Phi1] >=

18(3 4cos2

e22 + cos4

e82)

and

C =< sin2[Phi1]cos2[Phi1] >=

18(1 cos4

e82)

KDP =

180

[integraldisplay]0

Dmax

CkRe

fa (0) fb (0)

[bracketrightbig]N

(D)dD [ km1]. (7)

Here, Ck =< cos2[Phi1] >= cos2[Phi1]e22.

fa(0) and fb(0) are complex forward-scattering amplitudes and fa() and fb() are complex backscattering amplitudes for polarization along the major and minor axes, respectively. Additionally, fa and fb are their respective conjugates,

is the mean canting angle, and is the standard deviation of the canting angle. The terms

and were assumed to be 0 and 7 , respectively (Huang et al., 2008). The

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3870 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Table 3. List of different polarimetric rainfall relations used for rainfall estimations and the mean absolute error (MAE), root-mean-square error (RMSE), and correlation coefcient for estimated rain rates vs. observations.

R(Zh) = |Zh|

Polarimetric rainfall relation Scatterplot R Re Assumptions MAE RMSE Corr.

1 0.0568 0.5876 0.96 2.40 0.93 Pruppacher and Beard (1970)

2 0.0587 0.5849 0.97 2.42 0.92 Beard and Chuang (1987)

3 0.0588 0.5851 0.97 2.42 0.92 Brandes et al. (2002)

4 0.0575 0.5870 0.97 2.41 0.92 New axis ratio (experimental t)

R(KDP) = |KDP|

Polarimetric rainfall relation Scatterplot R Re Assumptions

MAE RMSE Corr.

1 38.59 0.834 0.47 1.05 0.99 Pruppacher and Beard (1970)

2 43.66 0.765 0.67 1.43 0.97 Beard and Chuang (1987)

3 46.85 0.740 0.82 1.64 0.97 Brandes et al. (2002)

4 42.20 0.831 0.46 1.15 0.98 New axis ratio (experimental t)

R(Zh, ZDR) = Z

h 100.1 ZDR

Polarimetric rainfall relation Scatterplot R Re Assumptions

MAE RMSE Corr.

1 0.0112 0.89 4.0964 0.45 0.77 0.99 Pruppacher and Beard (1970)

2 0.0094 0.88 3.5321 0.46 0.84 0.99 Beard and Chuang (1987)

3 0.0090 0.88 3.4908 0.47 0.86 0.99 Brandes et al. (2002)

4 0.0114 0.87 3.7750 0.48 0.90 0.99 New axis ratio (experimental t)

R(KDP, ZDR) = K

DP100

Polarimetric rainfall relation Scatterplot R Re Assumptions

MAE RMSE Corr.

1 66.56 0.96 1.4041 0.29 0.45 1.00 Pruppacher and Beard (1970)

2 84.02 0.93 1.6854 0.44 0.66 1.00 Beard and Chuang (1987)

3 96.57 0.92 1.9085 0.59 0.86 0.99 Brandes et al. (2002)

4 74.80 0.97 1.5489 0.23 0.36 1.00 New axis ratio (experimental t)

maximum size of the observed particles is Dmax = 8 mm, the

radar wavelength is = 10.7 cm (S-band), the dielectric fac

tor for water is Kw = 0.93, and N(D) was calculated by us

ing the 2DVD measurements.

Polarimetric rainfall relations between R and dual-polarimetric variables were derived when the rain rate was greater than 0.1 mm h1. The derived new polarimetric rainfall relations according to different axis-ratio relations are presented in Table 3.

3.4 Calibration of the radar

The radar measurements are affected by various observational errors, such as ground echoes, beam broadening, anomalous propagation echoes, and calibration biases of radar ZH and ZDR. These errors can lead to signicant un-

certainty in precipitation estimations. Therefore, accommodation of the calibration bias of radar is necessary to improve quantitative rainfall estimations involving ZH and/or ZDR.

The radar calibration was done for light rainfall events, and the simulated theoretical ZH and ZDR variables based on the new axis-ratio relation (Eq. 4) were used.

The calibration biases of ZH and ZDR were calculated from the comparison of measured ZH and ZDR with the simulated ZH and ZDR from the 2DVD measurements. To compare polarimetric radar parameters, the cross-match point must rst be determined. This is because 2DVD data consist of point measurements, whereas radar data are measured in a sampling volume. The BSL S-band radar data were averaged over ve successive gates and two adjacent azimuth angles centered on the 2DVD location. The elevation angle of 0.0 PPI was used. The radar rainfall estimation reliabil-

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

.1 ZDR

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3871

Figure 5. Different raindrop axis-ratio relations for the oblate raindrop model. The upper right subgure illustrates the axis ratio of an oblate raindrop.

ity was assessed before and after the calibration to test the effectiveness of the calibration process.

4 Results

4.1 Comparison of raindrop axis-ratio relations

We compared the new axis-ratio experimental t with existing mean axis-ratio relations such as those of Pruppacher and Beard (1970), Beard and Chuang (1987), and Brandes et al. (2002) (see Fig. 5). Consisting in the following polynomial formulas:

b/a = 1.03 0.062D(1 D 9mm) (8) b/a = 1.0048 + 5.7 104D 2.628 102D2 +3.682 103D3 1.677 104D4(1 D 7mm) (9) b/a = 0.9951 + 0.02510D 0.03644D2 +5.303 103D3 2.492 104D4(1 D 8mm). (10) Equation (8) from Pruppacher and Beard (1970) is a linear relation derived by measurements in a wind tunnel. Equation (9) is a fourth-order polynomial formula from a numerical model. Equation (10) is a polynomial empirical relation developed by Brandes et al. (2002) that was derived by combining drop shape observations. The Pruppacher and Beard (1970) linear relation results (Eq. 8, green dash-dot line) fell below the new mean axis-ratio for 1 D < 5 mm.

The Beard and Chuang (1987) polynomial relation results (Eq. 9, black dashed line) ranged from 2.5 to 6.5 mm lower than the new mean axis-ratio relation values. The new mean axis-ratio t was more oblate than that obtained with the Brandes et al. (2002) polynomial empirical relation (Eq. 10, blue dotted line) when the raindrop sizes were greater than

5.5 mm and less than 2.5 mm. With the exception of this case, the new axis ratio yielded results similar to those from Eq. (10) for diameters ranging from 3 to 5.5 mm. The new mean axis-ratio relation was very similar to existing axis-ratio relations except for small particles ( 2 mm) and large

particles ( 5.5 mm).

The differences in the mean axis-ratio relations seemed small, and dependence of the raindrop axis ratio on climatic regimes was not clearly observed, although raindrop shapes can be inuenced by the temperature and pressure (Beard and Chuang, 1987). Differences of raindrop shape can also be caused by measurement errors, drop oscillation, event selection criterion, and the tting method used (Thurai and Bringi, 2005). A small error in the mean axis ratio can lead to significant errors in the estimated DSD and rainfall rates (Bringi and Chandrasekar, 2001). Therefore, the consideration of the accurate raindrop shapes based on long-term DSD data is necessary to improve the polarimetric rainfall estimations.

4.2 Verication of polarimetric rainfall algorithms

4.2.1 Variability of DSD in rainfall estimations

To investigate the variability of the DSD in rainfall estimations from polarimetric parameters, the rain rate Re was estimated by utilizing various combinations of polarimetric variables and then was compared with R derived from the 2DVD results (Eq. 2). The MAE, RMSE, and correlation coefcient (Corr.) were calculated according to the following equations:

MAE =

1 N

X|R Re| [mm h1] (11)

RMSE =

[parenleftbigg]

1 N

X(R Re)2[parenrightbigg]0.5 [mm h1] (12)

Corr. =

1N 1

[summationtext][bracketleftbig] R

pVar(R)Var(Re)

R

[parenrightbig] R

e

Re

[parenrightbig][bracketrightbig]

, (13)

where R is the rain rate from observed 1-min DSDs and Re is the rain rate from various polarimetric rainfall algorithms.Re can then be obtained from the 2DVD data. The N is the number of comparisons. Figure 6 shows the scatterplot of R and Re for polarimetric rainfall algorithms based on the new axis-ratio relation. The scatter indicates the effect of DSD variability on rain estimations. The comparison between rain rates from the observed R and those estimated from Re showed good overall agreement. In particular, the statistics for the scatter showed the best results (MAE = 0.23 mm h1,

RMSE = 0.36 mm h1, and Corr. = 1.00) when using the

R(KDP, ZDR) based on the new axis-ratio relation. The use of the single parameter R(ZH) resulted in an increase of the MAE (0.97 mm h1) and RMSE (2.41 mm h1), but in a decrease of the Corr. (0.93) compared to the other polarimetric rainfall relations. Other polarimetric rainfall algorithms based on existing axis-ratio relations (Eqs. 810) showed similar results. A summary of the statistics according to polarimetric rainfall algorithms and mean axis-ratio relations

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3872 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Figure 6. Scatterplot of R derived from observed DSDs (17 618 min samples) and Re estimated from combinations of polarimetric parameters. Re was then obtained from the same data set.

are presented in Table 3. The differences in performance between the dual-polarization rainfall algorithms under practical conditions can also be due to the different error structures of dual-polarization measurements.

The reectivity factor is affected by the absolute calibration error, and it requires accurate knowledge of the radar constant. The differential reectivity is a relative power measurement of any rainfall algorithm. Therefore, it can be measured without being affected by absolute calibration errors.However, the ZDR-based algorithm needs to be used in conjunction with Zh or KDP because ZDR is a relative power measurement of any rainfall algorithm. Unlike Zh and ZDR,

KDP is independent of the absolute calibration error attenuation because it is related to the phase shift of the electromagnetic wave. However, KDP is relatively noisy in light rain (a low rain rate). Thus, the advantages and disadvantages of each polarimetric variable translate into the error of rainfall algorithms (Bringi and Chandrasekar, 2001).

4.2.2 Validation of rainfall estimation algorithms

In order to evaluate radar rainfall estimations according to different rainfall relations and raindrop shapes, we compared the 1 h rain rate from the BSL S-band radar with the hourly rain rate obtained by the use of the rain gauge in Daegu, Korea. Furthermore, the rainfall estimate from the 2DVD data was included for comparison. Here, the time of observations was short and unstable rainfall events in regard to the conti-

nuity and stability of measurement data were excluded from the analysis. Therefore, Statistical validation of the radar and 2DVD rainfall estimates for the different rainfall relations were performed for 18 rainfall events among the 33 rainfall cases. The MAE and RMSE are given by

MAE =

1 N

X|RR RG| [mm h1] (14)

X(RR RG)2[parenrightbigg]0.5 [mm h1], (15)

where RR is obtained by the radar (or 2DVD) averaged 1 h rain rate (mm h1), and RG is measured by the rain gauge averaged 1 h rain rate (mm h1). The results are presented in Table 4. According to the DSD statistics, the combined polarimetric rainfall relations using ZDR and KDP were more efcient than the single-parameter rainfall relation for estimated 2DVD rainfall. In particular, when using the new axis-ratio relation, the R(KDP, ZDR) relation performed better than the others. These results can be found in Fig. 7, which shows a scatterplot of the 1 h rain rate from the rain gauge and the radar (or 2DVD) data developed by using the new mean axis-ratio relation for 18 rainfall events. The pluses represent 1 h gauge rain rates vs. radar rain rates from different rainfall relations and the squares indicate gauge vs. 2DVD rain rates. In the single rainfall relation with Zh (Fig. 7a), some amount of scatter was present (MAE = 0.96 mm h1

and RMSE = 1.24 mm h1), and the scatter decrease

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

RMSE =

[parenleftbigg]

1 N

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3873

Table 4. Mean absolute error (MAE) and root-mean-square error (RMSE) of the radar estimates of hourly rain rates for the different radar rainfall algorithms listed in Table 3.

R(ZH) = |ZH|

MAE RMSE Assumptions

RADAR 2DVD RADAR 2DVD

1 1.02 0.96 1.39 1.24 Pruppacher and Beard (1970)

2 1.03 0.96 1.39 1.25 Beard and Chuang (1987)

3 1.03 0.96 1.39 1.25 Brandes et al. (2002)

4 1.02 0.96 1.39 1.24 New axis ratio (experimental t)

R(KDP) = |KDP|

MAE RMSE Assumptions

RADAR 2DVD RADAR 2DVD

1 6.04 0.69 6.99 0.93 Pruppacher and Beard (1970)

2 7.69 0.75 8.63 0.99 Beard and Chuang (1987)

3 8.67 0.76 9.62 1.00 Brandes et al. (2002)

4 6.80 0.71 7.78 0.93 New axis ratio (experimental t)

R(ZH, ZDR) = ZH

100.1 ZDR

MAE RMSE Assumptions

RADAR 2DVD RADAR 2DVD

1 0.89 0.66 1.22 0.92 Pruppacher and Beard (1970)

2 0.85 0.69 1.18 0.96 Beard and Chuang (1987)

3 0.87 0.73 1.20 1.02 Brandes et al. (2002)

4 0.85 0.72 1.17 0.99 New axis ratio (experimental t)

R(KDP, ZDR) = KDP

100.1 ZDR

MAE RMSE Assumptions

RADAR 2DVD RADAR 2DVD

1 8.19 0.58 10.05 0.88 Pruppacher and Beard (1970)

2 10.92 0.61 12.98 0.91 Beard and Chuang (1987)

3 12.71 0.64 14.96 0.94 Brandes et al. (2002)

4 9.19 0.58 11.23 0.86 New axis ratio (experimental t)

(MAE = 0.58 mm h1 and RMSE = 0.86 mm h1) when the

KDP and ZDR parameters were used for 2DVD rainfall estimation (Fig. 7d). These results were inuenced by the variability of the DSDs, and the effect of the DSD variability was worse in rainfall estimations obtained with the R(KDP) or

R(KDP, ZDR) than those obtained with the R(Zh). However, if actual radar measurements were used, the R(ZH, ZDR) algorithm showed the best performance, while the accuracy of radar rainfall estimation was reduced when the KDP parameter was used. As shown in Fig. 8, at lower rain rates (< 10 mm h1), R(KDP) and R(KDP, ZDR) were less efcient because the measured KDP became too noisy; however,

R(KDP) and R(KDP, ZDR) relations performed better than others at higher rain rates ( 10 mm h1). Therefore, as the

rain rate increased, the uncertainty in radar estimates due to the use of KDP was reduced. Zh (instead of KDP) should be

used jointly with ZDR because R(ZH, ZDR) showed excellent performance at lower rain rates. These results show that the R(KDP, ZDR) relation will be useful for heavy rainfall, and R(Zh, ZDR) is more suited for light rainfall (Ryzhkov et al., 2005).

4.2.3 Correction of calibration bias

Radar measurements can suffer from ZH and ZDR system biases, and inaccurate calibration can lead to signicant uncertainty in radar-based quantitative precipitation estimations. In this study, ZH and ZDR calibration biases were calculated for eight rainfall events through a comparison between ZH and ZDR radar measurements and ZH and ZDR simulations obtained from the 2DVD measurements. Figure 9a and b shows time series data for the 2DVD and BSL radar ZH and

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3874 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Figure 7. Scatter plot of the 1 h rain rate from rain gauge (RG) and BSL S-band radar (or 2DVD) based on Eq. (4) for 18 rainfall cases. The pluses represents 1 h gauge rain rates vs. radar hourly rain rates from polarimetric rainfall algorithms, and squares indicate gauge and 2DVD rain rates obtained by different polarimetric rainfall algorithms.

Figure 8. Scatter plots of the 1 h rain rates from rain gauge (RG) and BSL S-band radar based on Eq. (4) for 18 rainfall cases. The pluses represents 1 h gauge rain rates vs. radar hourly rain rates from polarimetric rainfall algorithms for a light rain rate (R < 10 mm h1), and squares indicate those for a heavy rain rate (R 10 mm h

1).

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3875

Figure 9. Time series of the (a) reectivity and (b) differential reectivity obtained by 2DVD and BSL S-band radar, and scatter plots of the 2DVD estimations and radar measurements for the (c) reectivity and (d) differential reectivity. Comparison statistics including the correlation coefcient (Corr), RMSE, and bias are also presented (measurements were collected on 23 August 2012).

Figure 10. Same as Fig. 9, except that the data are for 16 September 2012.

ZDR, and Fig. 9c and d shows scatter plots of the 2DVD estimations and radar measurements for the ZH and ZDR on23 August 2012. The overall distribution of the observed ZH corresponded well with the simulated parameter; however, the BSL ZDR measurements were underestimated with respect to the simulated ZDR values. The mean bias (= bias) of

ZH and ZDR was about 0.98 dBZ and 0.10 dB, respectively. The bias of ZH and ZDR on 16 September 2012 was 0.90 dBZ and 0.13 dB, respectively (Fig. 10).

The application of absolute calibration biases is the most effective way to reduce radar rainfall errors, in particular, for rainfall estimators with both ZH and ZDR (Kwon et al., 2015). The accuracy of the radar rainfall estimations was investigated by applying the calculated ZH and ZDR bias. The

R(Zh, ZDR) algorithm based on the new mean axis-ratio relation showed the best result in terms of radar rainfall estimations. These results were evaluated by comparison with rain gauge measurements, which represent the ground truth situ-

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3876 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Table 5. Mean absolute error (MAE) and root-mean-square error (RMSE) of rainfall estimates before and after applying bias correction.

Date Calibration Radar rainfall Rain MAE RMSE bias estimation gauge (mm h1) (mm h1)

Zh ZDR Before After Obs. Before After Before After (dBZ) (dB) BC BC BC BC BC BC

1 11.10.13 0.01 0.10 20.71 18.99 16.85 0.50 0.40 0.61 0.49

2 11.10.21 1.16 0.28 48.60 48.26 49.96 0.77 0.76 0.99 0.98

3 12.04.25 1.40 0.43 19.83 18.07 17.24 1.33 1.24 2.03 1.90

4 12.05.14 2.17 0.28 12.89 15.58 19.40 0.57 0.42 0.78 0.58

5 12.08.23 0.98 0.10 72.65 81.35 83.86 0.78 0.62 1.27 0.87

6 12.09.09 0.44 0.18 18.16 17.00 14.17 0.75 0.68 1.18 1.09

7 12.09.16 0.90 0.13 44.90 53.77 55.20 0.75 0.69 0.92 0.84

8 12.10.22 1.44 0.14 18.55 15.76 14.92 2.85 2.60 3.49 3.30

Avg 1.03 0.22 1.04 0.93 1.41 1.26

*BC: bias correction.

bias correction, respectively. The bar graph is the 1 h rainfall rate measured by the rain gauge. The blue and green dashed lines are estimated accumulated rainfall obtained before and after the bias correction, respectively, and the red dashed line is the accumulated rainfall measured by the rain gauge. In comparison with the rain gauge measurement (83.86 mm), the underestimated precipitation (72.65 mm) was corrected to 81.35 mm after bias correction. When the estimated rainfall was compared to the rain gauge data, rainfall estimation was improved by about 10.38 %. Figure 11b shows the results for the 16 September 2012 event. The recorded accumulated rainfall amounts were 55.20, 44.90, and 53.77 mm for the rain gauge, radar before bias correction, and radar after bias correction results, respectively. The radar rainfall estimation was improved by about 16.07 %. For the eight considered rainfall events, the ZH and ZDR total mean biases for the BSL radar was 1.03 dBZ and 0.22 dB, respectively. Moreover, MAE changed from 1.04 to 0.93 mm h1 and RMSE decreased from 1.41 to 1.26 mm h1 after the bias correction. The bias of ZH and ZDR as well as MAE and RMSE values for each of the eight rainfall events are presented in Table 5. As shown in Table 5, both MAE and RMSE values decreased after the bias correction, which, therefore, lead to an improvement of rainfall estimation.

5 Conclusions

The purpose of this study was to nd an optimal polarimetric rainfall algorithm by using 2DVD measurements in Korea, and to improve the radar rainfall estimations by correcting the ZH and ZDR biases. First, we derived a new raindrop axis-ratio relation reecting the rainfall characteristics on the Korean Peninsula by using data from 33 rainfall events, and this was done after checking the accuracy and performing a quality control procedure for the 2DVD measurements. The derived raindrop axis-ratio relation was compared with existing

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/

Figure 11. Comparison of the 1 h rain rate (left ordinate) and accumulated rainfall (right ordinate) obtained by the BSL S-band radar and rain gauge. The R(ZH, ZDR) algorithm based on Eq. (4) was used for rainfall estimations for (a) 23 August 2012 and (b) 16 September 2012. BC represents the bias correction.

ation. Figure 11a shows the 1 h rain rate (left ordinate) and accumulated rainfall (right ordinate) estimated from the radar and rain gauge on 23 August 2012. The blue and green solid lines are estimated 1 h rainfall rates obtained before and after

H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea 3877

relations. The derived new mean axis-ratio relation was very similar to existing axis-ratio relations except for both small particles ( 2 mm) and large particles ( 5.5 mm). The de

pendence of the raindrop axis ratio on climatic regimes was not clearly observed.

The polarimetric rainfall algorithms were derived based on various assumptions about the shape of raindrops. The accuracy validation of the 1 h rainfall rate obtained through rainfall algorithms was assessed by comparing 2DVD and BSL radar data with rain gauge measurements. As a result, R(KDP, ZDR) based on the new axis-ratio relation was deemed suitable for rainfall estimations according to the DSD statistics when compared with others. This occurs because the effect of the DSD variability declined in rainfall estimations obtained with the R(KDP) or R(KDP, ZDR)

compared to those obtained with the R(Zh). However, if actual radar measurements were used, at lower rain rates (< 10 mm h1), the KDP-based algorithms displayed large statistical errors, while R(Zh, ZDR) based on the new axis-ratio relation showed the best performance. This occurs because the measured KDP parameter was relatively noisy at lower rain rates. However, the R(KDP, ZDR) algorithm based on the new axis-ratio relation performed better than the others at higher rain rates ( 10 mm h1). Therefore, in order to

produce more accurate rainfall estimations, the R(Zh, ZDR) and R(KDP, ZDR) algorithm should be classied according to rainfall intensities.

To perform radar calibration measured ZH and ZDR were compared with simulated ZH and ZDR. The calculated ZH and ZDR bias was used to reduce radar bias, and to produce more accurate rainfall estimations. After bias correction, radar rainfall estimations were closer to rain gauge measurements, meaning an improvement of the rst ones.

In this paper, different raindrops axis ratios were used to derive new polarimetric rainfall relations, and the new polarimetric rainfall algorithms were assessed with respect to their ability to produce accurate point radar rainfall estimations. The new polarimetric rainfall algorithms performed better than existing rainfall algorithms, and no large differences were observed in regard to climatic regimes. In particular, the use of the R(Zh, ZDR) algorithm lead to an improvement of radar rainfall estimations in the Korean Peninsula since in this region extreme precipitation systems such as tornados, supercells, and hail storms do not often occur.We are now working to improve the KDP quality because it is a very useful parameter for high rainfall rates estimation. A classication of rain rates based on these results will also be performed in future work. In addition, polarimetric rainfall algorithms will be developed to obtain areal rainfall estimations through long-term DSD collection efforts and several rainfall-cases-based analyses.

Acknowledgements. This research was supported by the Development and application of cross governmental dual-pol radar

harmonization (WRC-2013-A-1) project of the Weather Radar Center, Korea Meteorological Administration.

Edited by: G. VulpianiReviewed by: S. Sebastianelli and four anonymous referees

References

Adirosi, E., Gorgucci, E., Baldini, L., and Tokay, A.: Evaluation of gamma raindrop size distribution assumption through comparison of rain rates of measured and radar-equivalent gamma DSD,J. Appl. Meteorol. Clim., 53, 16181635, 2014.

Atlas, D.: Radar calibration: some simple approaches, B. Am. Meteorol. Soc., 83, 10131316, 2002.

Atlas, D., Srivastava, R. C., and Sekkon, R. S.: Doppler radar characteristics of precipitation at vertical incidence, Rev. Geophys.Space Ge, 2, 135, 1973.

Beard, K. B. and Chuang, C.: A new model for the equilibrium shape of raindrops, J. Atmos. Sci., 44, 15091524, 1987.Brandes, E. A., Zhang, G., and Vivekanandan, J.: Experiments in rainfall estimation with a polarimetric radar in a subtropical environment, J. Appl. Meteorol., 41, 674685, 2002.

Bringi, V. N. and Chandrasekar, V.: Polarimetric Doppler Weather

Radar: Principles and Applications, Cambridge University Press, New Yprk, USA, p. 636, 2001.

Bringi, V. N., Chandrasekar, V., Hubbert, J., Gorgucci, E., Randeu,W. L., and Schoenhuber, M.: Raindrop size distribution in different climatic regimes from disdrometer and dual-polarized radar analysis, J. Atmos. Sci., 60, 354365, 2003.

Chandrasekar, V., Bringi, N., Balakrishnan, and Zrni, D. S.: Error structure of multiparameter radar and surface measurements of rainfall, Part III: Specic differential phase, J. Appl. Meteorol., 32, 12881293, 1990.

Chang, W. Y., Wang, T. C., and Lin, P. L.: Characteristics of the raindrop size distribution and drop shape relation in typhoon systems in the western Pacic from the 2D video disdrometer and NCU C-band polarimetric radar, J. Atmos. Ocean. Tech., 26, 1973 1993, 2009.

Cifelli, R., Chandrasekar, V., Lim, S., Kennedy, P. C., Wang, Y., and

Rutledge, S. A.: A new dual-polarization radar rainfall algorithm: Application in Colorado precipitation events, J. Atmos. Ocean.Tech., 28, 352364, 2011.

Goddard, J. W. F. and Cherry, S. M.: The ability of dual-polarization radar (copular linear) to predict rainfall rate and microwave attenuation, Radio Sci., 19, 201208, 1984.

Goddard, J. W. F., Cherry, S. M, and Bringi, V. N.: Comparison of dual-polarized radar measurements of rain with ground-based disdrometer measurements, J. Appl. Meteorol., 21, 252 256, 1982.

Gorgucci, E., Scarchilli, G., Chandrasekar, V., and Bringi, V. N.:

Rainfall estimation from polarimetric radar measurements: Composite algorithms immune to variability in raindrop shape-size relation, J. Atmos. Ocean. Tech., 18, 17731786, 2001.

Green, A. W.: An approximation for the shapes of large rain-drops,J. Appl. Meteorol., 14, 15781583, 1975.

Hagen, M. and Yuter, S. E.: Relations between radar reectivity, liquid water content, and rainfall rate during the MAPSOP, Q. J.Roy. Meteor. Soc., 129, 477493, 2003.

www.atmos-meas-tech.net/9/3863/2016/ Atmos. Meas. Tech., 9, 38633878, 2016

3878 H.-L. Kim et al.: Dual-polarization radar rainfall estimation in Korea

Huang, G. J., Bringi, V. N., and Thurai, M.: Orientation angle distributions of drops after an 80-m fall using a 2D video disdrometer,J. Atmos. Ocean. Tech., 25, 17171723, 2008.

Jameson, A. R.: Microphysical interpretation of multiparameter radar measurements in rain, Part I: Interpretation of polarization measurements and estimation of raindrop shapes, J. Atmos. Sci., 40, 17921802, 1983.

Jameson, A. R.: Microphysical interpretation of multiparameter radar measurements in rain, Part III: Interpretation and measurement of propagation differential phase shift between orthogonal linear polarizations, J. Atmos. Sci., 42, 607614, 1985.

Joss, J., Thams, J. C., and Waldvogel, A.: The accuracy of daily rainfall measurements by radar, In: pre-prints, 13th Radar Meteorology Conf., Montreal, QC, Canada, Amer. Meteor. Soc., 448 451, 1968.

Jung, Y., Xue, M., and Zhang, G.: Simulations of polarimetric radar signatures of a supercell storm using a two-moment bulk micro-physics scheme, J. Appl. Meteorol. Clim., 49, 146163, 2010.Keenan, T. D., Carey, L. D., Zrni, D. S., and May, P. T.: Sensitivity of 5-cm wavelength polarimetric radar variables to raindrop axial ratio and drop size distribution, J. Appl. Meteorol., 40, 526545, 2001.

Kruger, A. and Krajewski, W. F.: Two-dimensional video disdrometer, J. Atmos. Sci., 19, 602617, 2002.

Kwon, S., Lee, G. W., and Kim, G.: Rainfall estimation from an operational S-band dual-polarization radar: Effect of radar calibration, J. Meteorol. Soc. Jpn., 93, 6579, 2015.

Maki, M., Park, S. G., and Bringi, V. N.: Effect of natural variations in raindrop size distributions on rain rate estimators of 3 cm wavelength polarimetric radar, J. Meteorol., Soc. Jpn., 83, 871 893, 2005.

Marzuki, M., Randeu, W. L., Kozu, T., Shimomai, T., Hashiguchi,H., and Schonhuber, M.:Raindrop axis ratios, fall velocities and size distribution over Sumatra from 2D-Video Disdrometer measurement, Atmos. Res., 119, 2337, 2013.

McFarquhar, G. M. and List, R.: The effect of curve ts for the disdrometer calibration on raindrop spectra, rainfall rate, and radar reectivity, J. Appl. Meteorol., 32, 774782, 1993.

Park, S. G. and Lee, G. W.: Calibration of radar reectivity measurements from the KMA operational radar network, Asian-Pacic J.Atmos. Sci., 46, 243259, 2010.

Pruppacher, H. and Beard, K. V.: A wind tunnel investigation of the internal circulation and shape of water drops falling at terminal velocity in air, Q. J. Roy. Meteor. Soc., 96, 247256, 1970.Pruppacher, H. R. and Pitter, R.: A semi-empirical determination of the shape of cloud and rain drops, J. Atmos. Sci, 28, 8694, 1971.

Raupach, T. H. and Berne, A.: Correction of raindrop size distributions measured by Parsivel disdrometers, using a two-dimensional video disdrometer as a reference, Atmos. Meas. Tech., 8, 343365, doi:http://dx.doi.org/10.5194/amt-8-343-2015

Web End =10.5194/amt-8-343-2015 http://dx.doi.org/10.5194/amt-8-343-2015

Web End = , 2015.

Ryzhkov, A. V. and Zrni, D. S.: Comparison of dual-polarization radar estimators of rain, J. Atmos. Ocean. Tech., 12, 249256, 1995.

Ryzhkov, A. V., Giangrande, S. E., and Schuur, T. J.: Rainfall estimation with a polarimetric prototype of WSR-88D, J. Appl. Meteorol., 44, 502515, 2005.

Sachidananda, M. and Zrni, D. S.: Rain rate estimates from differential polarization measurements, J. Atmos. Ocean. Tech., 4, 588598, 1987.

Sebastianelli, S., Russo, F., Napolitano, F., and Baldini, L.: On precipitation measurements collected by a weather radar and a rain gauge network, Nat. Hazards Earth Syst. Sci., 13, 605623, doi:http://dx.doi.org/10.5194/nhess-13-605-2013

Web End =10.5194/nhess-13-605-2013 http://dx.doi.org/10.5194/nhess-13-605-2013

Web End = , 2013.

Sheppard, B. E. and Joe, P. I.: Comparison of raindrop size distribution measurements by a Joss-Waldvogel disdrometer, a PMS 2DG spectrometer, and a POSS Doppler radar, J. Atmos. Ocean. Tech., 11, 874887, 1994.

Thurai, M. and Bringi, V. N.:Drop axis ratios from a 2D video disdrometer, J. Atmos. Ocean. Tech., 22, 966978, 2005.

Tokay, A. and Short, D. A.: Evidence from tropical raindrop spectra of the origin of rain from stratiform versus convective clouds, J. Appl. Meteorol., 35, 355371, 1996.

Tokay, A., Wolff, D. B., Wolff, K. R., and Bashor, P.: Rain gauge and disdrometer measurements during the Keys Area Microphysics Project (KAMP), J. Atmos. Ocean. Tech., 20, 14601477, 2003.

Ulbrich, C. W.: Natural variations in the analytical form of the raindrop size distribution, J. Clim. Appl. Meteorol., 22, 17641775, 1983.

Ulbrich, C. W. and Lee, L. G.: Rainfall measurement error by WSR-

88D radars due to variations in ZR law parameters and the radar constant, J. Atmos. Ocean. Tech., 16, 10171024, 1999. Villarini, G. and Krajewski, W. F.: Review of the different sources of uncertainty in single polarization radar-based estimates of rainfall, Surv. Geophys., 31, 107129, 2010.

Zawadzki, I.: Factors affecting the precision of radar measurements of rain, In: Preprints, 22nd Radar Meteorology Conference, September 1013, Zurich, Switzerland, Boston, MA, Amer. Meteor. Soc., 251256, 1984.

Zhang, G., Vivekanandan, J., and Brandes, E.: A method for estimating rain rate and drop size distribution from polarimetric radar measurements, IEEE T. Geosci. Remote, 39, 830841, 2001.

Atmos. Meas. Tech., 9, 38633878, 2016 www.atmos-meas-tech.net/9/3863/2016/