# Radar-Derived Quantitative Precipitation Estimation Based on Precipitation Classification

Authors: Lili Yang; Yi Yang; Peng Liu; Lina Wang
Journal: Advances in Meteorology
 (2016)

---

## Abstract

A method for improving radar-derived quantitative precipitation estimation is proposed. Tropical vertical profiles of reflectivity (VPRs) are first determined from multiple VPRs. Upon identifying a tropical VPR, the event can be further classified as either tropical-stratiform or tropical-convective rainfall by a fuzzy logic (FL) algorithm. Based on the precipitation-type fields, the reflectivity values are converted into rainfall rate using aZ-R relationship. In order to evaluate the performance of this rainfall classification scheme, three experiments were conducted using three months of data and two study cases. In Experiment I, the Weather Surveillance Radar-1988 Doppler (WSR-88D) defaultZ-R relationship was applied. In Experiment II, the precipitation regime was separated into convective and stratiform rainfall using the FL algorithm, and correspondingZ-R relationships were used. In Experiment III, the precipitation regime was separated into convective, stratiform, and tropical rainfall, and the correspondingZ-R relationships were applied. The results show that the rainfall rates obtained from all three experiments match closely with the gauge observations, although Experiment II could solve the underestimation, when compared to Experiment I. Experiment III significantly reduced this underestimation and generated the most accurate radar estimates of rain rate among the three experiments.


---

## Body

## 1. Introduction

Owing to its significant impacts on human activities, precipitation is one of the most important factors in meteorological analysis; accordingly, its study has attracted considerable attention. Quantitative precipitation estimation (QPE) plays a very important role in generating warnings of meteorological, hydrological, and geological disasters. The accuracy of numerical weather predictions can be significantly improved when the precipitation rate is assimilated into the weather prediction models [1, 2]. Currently, precipitation rate data is obtained primarily from rain gauge networks, which can provide real-time rainfall estimates. However, the spatial distribution of rain gauges is sparse in remote regions, particularly in mountainous areas, which can result in insufficient spatial resolution for accurate mapping of the rainfall patterns. Consequently, obtaining QPE products with high spatial and temporal resolution is currently one of the main tasks in conducting hydrometeorological studies [3, 4]. Weather radars, in particular, can provide precipitation estimates with high spatial and temporal resolution over a large area [5]. Radar reflectivity,Z (mm6/m3), is related to rainfall rate,R (mm/h), through a variable power law relationship (Z-R relationship), which can be used to obtain QPE products.NumerousZ-R relationships based on raindrop size distribution observations have been proposed [6–8].Z-R relationships of the formZ =ARB can be derived from pairs ofZ andR observations, assuming that sufficient corresponding pairedZ-R data samples are available. The default WSR-88DZ-R relationship,Z = 300R1.4 [9–11], was derived based on summer convective events in Florida, USA. However, this defaultZ-R relationship cannot fully address the complicated characteristics of rain in, for example, China [12–14]. It is known that different types of precipitation correspond to different raindrop size distributions (DSDs), and precipitation types can also be classified based on the formation mechanism, duration, and internal structure. The precipitation regime can therefore be divided into different rainfall types. Convective rainfall systems are associated with strong vertical velocity fields, small areal coverage, and high rainfall intensities, and an echo that can be extended to a relatively high altitude and with a high center of mass. Conversely, stratiform rainfall systems have relatively weak vertical velocity fields, great horizontal homogeneity, and low rainfall intensity. Although the intensity of stratiform rainfall is much weaker than that in adjacent convective cells, stratiform rain typically covers a larger area and always contributes with a significant portion (40–50%) of the total rainfall, even in convective-dominant systems [15]. Tropical rain systems generally occur in lower altitudes and very humid air, with condensation occurring continuously throughout the air descent. The main feature of tropical rain systems is that their reflectivity increases monotonically as height decreases, yielding a maximum reflectivity at the lowest levels of the vertical profile of reflectivity (VPR), resulting in considerable error when conducting precipitation estimation using the defaultZ-R relationship of WSR-88D. Moreover, multiple types of precipitation may coexist within a wide range of rainfall events.Radar research scientists worldwide have undertaken extensive studies to improve the accuracy of radar QPE [16–23]. Typically, precipitation is first classified into various types, and corresponding multipleZ-R relationships are then used [19, 24–26]. Precipitation can generally be classified into either convective or stratiform regimes, which has been a particular focus of many studies [27–29]. Such basic classifications typically adopt a fixed threshold or set of boundary conditions for rainfall recognition. However, these methods are known to be sensitive to the used threshold values and because there is a considerable overlap between stratiform and convective regimes in many respects, using fixed boundary conditions or thresholds often leads to false positives. Moreover, it is often difficult to directly identify the boundary between stratiform and convective regimes [30]. Yang et al. [30] used the reference values of Steiner et al. [29] to derive a statistical relationship between four characteristic parameters and two rain types, in order to develop the probability distribution of convective and stratiform rainfall based on a fuzzy-logic (FL) algorithm, which was able to improve the accuracy of QPE [31]. Despite efforts such as these, many techniques can still underestimate rainfall to some extent. Based on the stratiform and convective precipitation classification, many researchers have suggested focusing on tropical rain as an avenue to mitigate the underestimation of QPE, and some studies have been able to reduce QPE underestimation by employing an automated radar technique to identify tropical precipitation, in order to divide precipitation into tropical rain and stratiform and convective types [32, 33]. Qi et al. [34] converted reflectivity into precipitation utilizing the tropical rainZ-R relationship instead of the default WSR-88DZ-R relationship, which improved the accuracy to some extent. Similarly, Zhang et al. [35] improved the accuracy of radar-derived QPE by using an automatic identification method to identify tropica1, stratiform, and convective precipitation types based on radar observations and model outputs. They then applied variousZ-R relationships in order to convert reflectivity into rain rates, based on precipitation types. The results showed that their method improved the accuracy of radar-derived QPE.In this study, precipitation is first divided into stratiform and convective precipitation types using the FL algorithms proposed by Yang et al. [30], and then tropical rain is identified based on the VPR characteristics [36, 37], as discussed by Xu et al. [33]. Tropical rainfall is further divided into tropical-stratiform and tropical-convective rainfall, and the respective adaptiveZ-R relationships are then applied for QPE. A three-month dataset is used to compare the newly proposed tropical rain identification scheme with the results obtained with the default WSR-88DZ-R relationship and a convective/stratiform rain identification scheme. Section 2 describes the precipitation classification method and the correspondingZ-R relationships and discusses the benefits of precipitation classification. Section 3 illustrates the evaluation results of using long-term datasets using descriptive statistics. Two case studies are also provided in Section 4 to evaluate the performance of the proposed precipitation classification method. Finally, Section 5 presents the main results of the study.


## 2. Algorithm Description

Underestimation of rainfall can occur when precipitation is divided into general stratiform and convective types. This study aims to improve the accuracy of radar QPE by building upon the convective, stratiform, and tropical precipitation technique described by Xu et al. [33]. The proposed method uses volume scan reflectivity data from a Doppler radar to recognize tropical VPRs, by first dividing rain data into stratiform and convective precipitation (using an FL algorithm) and then classifying tropical rain events as stratiform, convective, or tropical precipitation. Figure 1 shows a flow chart illustrating this process of precipitation classification.Figure 1
Flow chart of the precipitation classification method devised in this study.


### 2.1. Identification of Tropical VPRs

This study uses VPRs to identify the potential presence of tropical rain. As a first step, volume scan reflectivity data from Doppler radar are simply classified into assumed stratiform or convective types, based on the vertically integrated liquid (VIL) water content. In particular, the data are classified as belonging to the stratiform type when the VIL value is below 6.5 kg/m2 [38]; otherwise, they are classified as convective. Stratiform reflectivity observations larger than 10 dBZ from all tilts in an annular region between two predefined ranges, r1 and r2 (with the empirical values of 20 and 80 km, resp.), are used to obtain an averaged VPR [39]. The number of data samples within a 20-km radius of the radar is limited, owing to a terrain clearance rule incorporated into a quality control (QC) scheme [40], whereby any radar bin whose bottom is within 50 m of the ground is removed. Furthermore, the region of interest has to be far enough from the radar receiver, to avoid the cone of silence in the receiver’s immediate vicinity. The accuracy of precipitation classification is also limited by the decrease in vertical resolution beyond an 80-km radius.The second step combines the output of the first step with the method of bright band identification proposed by Zhang et al. [39], to determine whether a bright band is present. The specific method adopted for this step can be summarized as follows:(1)
The local maximum reflectivity (Zpeak) and its height (hpeak) in the VPR are found. The search for the local maximum in the VPR starts at 500 m above the 0°C isotherm altitude (h_0°C) obtained from the Weather Research and Forecasting (WRF) model and continues downward. A 500-m cushion is used to account for uncertainties in the model’s 0°C height.(2)
Once the local maximum is found, the heighthtop (hbottom) above (below) the maximum level at which reflectivity decreases monotonically by a given percentage (default = 10%) [39] of the maximum reflectivity (in dBZ) is determined, along with its corresponding reflectivity Ztop (Zbottom).(3)
If the following three criteria are met, a bright band is considered to exist:(1)htop-hbottom≤D0,htop-hpeak≤D1,hpeak-hbottom≤D1.The parametersD0 and D1 are assumed to be dependent on both the vertical resolution of the reflectivity observations and the radar scan strategy. Here, D0 and D1 are set to 1.5 and 1 km, respectively.(4)
If a bright band exists, its top (BBt) and bottom (BBb) heights are defined as follows:(2)BBt=htop,htop-hpeak≤Dt,hpeak+Dt,htop-hpeak>Dt,BBb=hbottom,hpeak-hbottom≤Db,hpeak-Db,hpeak-hbottom>Db.Here, Dt and Db serve as top and bottom caps of the bright band (BB), with default values of 500 and 700 m, respectively [39].To avoid excessive corrections of the bright band, the height of the local maximum reflectivity is searched again if the minimum reflectivity is lower than a preset threshold (default = 28 dBZ). This height is defined as the lowest height at which the reflectivity is larger than or equal to the threshold [38]. This parameter can sometimes lead to insufficient or excessive bright band correction, an issue that will be solved in the future using observations from a dual-polarization radar [41].The third step involves the correction of any existing bright bands [42]. If the values of BBtBBb and their corresponding reflectivities (ZBt and ZBb, resp.) are valid, BBb<hpeak, and BBt>hpeak, and then the slope (S) of reflectivity between BBt and BBb is calculated as follows:(3)S=ZBt-ZBbBBt-BBb.If the reflectivity heighthi is located between BBt and BBb, the reflectivity is reassigned as(4)Zi=ZBt-S×BBt-hi.IfBBt does not exist, the value of BBb is valid, and BBb<hpeak, and the reflectivities located above BBb (with value larger than ZBb) are reassigned to ZBb.The fourth step, which involves reflectivity observations from all tilts in an annular region between 20 and 80 km, involves recalculation of the average stratiform VPR after bright band correction.The fifth step involves determining whether a tropical VPR has been captured. Following Xu et al. [33], the identification (or diagnosis) of tropical VPR is made with a least-squares fit from the bottom of the bright band to the bottom of the VPR, in cases where a bright band exists; otherwise, the fit is made fromh_0°C to the bottom of the VPR. The slope (α) of the VPR is then calculated. The resulting VPR is considered as tropical if three out of five values of α (for the present and the four previous moments) are lower than or equal to 0; otherwise, it is regarded as a nontropical VPR.


### 2.2. Convective and Stratiform Rainfall Classification

To improve the accuracy of radar QPE, this study adopted the FL algorithm developed by Yang et al. [30] to classify convective and stratiform rainfall. Brief details of the FL algorithm are given below.First, the volume scan reflectivity data from the Doppler radar are interpolated to a Cartesian coordinate system with the same latitude and longitude range (0.01°) of the original data. Four features, denoted asF1–F4, are calculated based on the reflectivity data. These features are, respectively, (a) the reflectivity at 2 km, (b) the reflectivity standard deviation in the horizontal direction, (c) the product of the radar top height and the reflectivity value at 2 km, and (d) the vertically integrated liquid water content [30].As a second step, membership functions are used to determine the degree to which each feature belongs to each rain type, using a fuzzification process. Linear functions are used as membership functions for convection:(5)μk,Cx,a,b=0,x≤a,x-ab-a,a<x<b,1,x≥b,k=1,2,…,Nand stratification:(6)μk,Sx,a,b=1,x≤a,b-xb-a,a<x<b,0,x≥b,k=1,2,…,N,where N is the input feature parameter, subscript k is set to 1, 2, 3, or 4, C and S identify the convective and stratiform cases, respectively, x is the feature value, a is the left breaking point, and b is the right breaking point. Different feature parameters are assigned different breaking points. ForF1, the parameters are set to a=20 dBZ and b=45 dBZ; forF2, they area = 1 dBZ and b=14 dBZ; for F3, they are a = 100 km·dBZ andb = 500 km·dBZ; and, for F4, they area = 0.5 kg/km2 and b=5.0 kg/km2 [30]. From (5) and (6), it is clear that(7)μk,Sx,a,b=1-μk,Cx,a,b.Third, the weighted average values of the measurement belonging to each specific class are obtained as follows:(8)Pe=∑k=1NWk·Pk,e∑k=1NWk,where Pk,e=μk,ε(Xk), e identifies the convective or stratiform case, Wk is the weighting factor (with equal weighting assigned in this study for simplicity), N is the number of input feature parameters, and Pc and Ps indicate the weighted average values of the convective and stratiform classes, respectively. As the sum of Pc and Ps is 1, the value of Pc can be used to express classifications for convection in a probabilistic manner. Grid points for which Pc≥0.5 are classified as convective rainfall; all other grid points are classified as stratiform rainfall.


### 2.3. Identification of Tropical Rainfall Regions

To improve the accuracy of the radar QPE, tropical rainfall is considered to be present if tropical VPR is present and the grid at 1 km altitude meets the following two criteria [33]: (1) the reflectivity is greater than 25 dBZ and (2) the surface wet bulb temperature is greater than 2°C. The surface wet bulb temperature is calculated as follows:(9)E=6.11×107.5×Tdc/237.7+TdcWBc=0.00066×P×Tc+4098×E/Tdc×237.7+Tdc20.00066×P+4089×E/237.7+Tdc2,where Tdc is the dew point temperature (in °C),E is the vapor pressure from dew point (in mb), WBc is the surface wet bulb temperature (in °C), P is the pressure (in mb), and Tc is the air temperature (in °C). The above parameters are simulated using the WRF model or from other data.In this way, stratiform rainfall is further divided into stratiform and tropical-stratiform rainfall, and convective rainfall is further divided into convective and tropical-convective rainfall.


### 2.4. Radar QPE Calculation

Three sets of experiments were conducted to test the proposed method. In Experiment I, only the WSR-88D defaultZ-R relationship (Z = 300R1.4) was used. In Experiment II, as the rain type is identified by a FL algorithm, the correspondingZ-R relationships were applied to the rain rate calculation at each grid point. For each grid point, the relationship Z = 300R1.4 was applied for the convective rainfall calculation if Pc≥0.5. Otherwise, the relationship Z = 200R1.6 [6] was applied, for stratiform rainfall calculation. In Experiment III, stratiform and convective rainfall in tropical rainfall events were further divided into stratiform and tropical-stratiform rainfall and convective and tropical-convective rainfall, respectively. The relationships were the same as those used in Experiment II, except that the tropical rainfall relationship Z = 30.7R1.66 was adopted (Yadong Wang, personal communication). In an attempt to exclude the impact of nonprecipitation echoes and hail, quality control was conducted prior to precipitation classification. Additionally, reflectivity at 1 km was used to reduce the errors caused by spatial inconsistencies and the effects of bright band measurements not yet properly revised.


### 2.5. Performance Evaluation

The mean deviation error (BIAS), relative absolute error (RAE) rate, and root mean square error (RMSE) are selected as statistical indicators for quantitatively assessing the quality of radar QPE in the different experiments [43]. The mean deviation error is defined as follows:(10)BIAS=1n∑i=1nRai-Rgi,where Rai is the radar QPE of the ith station, Rgi is the observed precipitation of the ith site, and n is the total number of samples involved in the performance assessment.The radar QPE error may be positive or negative; therefore, its mean will not reflect the radar QPE performance. To objectively reflect performance, the mean absolute error rate is used and is defined as follows:(11)RAE=∑i=1nRai-Rgi∑i=1nRgi×100%.In addition, downpour and drizzle may coexist during a precipitation process; as the RAE of downpour is larger than that of drizzle, the root mean square error primarily represents the performance of the downpour component of the rainfall:(12)RMSE=∑i=1nRai-Rgi2n1/2.As the above three equations imply, a BIAS closer to 0 corresponds to more similar (in the average) QPE and observed rain types, whereas a smaller RAE corresponds to a higher radar QPE accuracy for the entire precipitation process. Similarly, a smaller RMSE indicates a more accurate radar QPE for the heavy precipitation portion of the process.


## 2.1. Identification of Tropical VPRs

This study uses VPRs to identify the potential presence of tropical rain. As a first step, volume scan reflectivity data from Doppler radar are simply classified into assumed stratiform or convective types, based on the vertically integrated liquid (VIL) water content. In particular, the data are classified as belonging to the stratiform type when the VIL value is below 6.5 kg/m2 [38]; otherwise, they are classified as convective. Stratiform reflectivity observations larger than 10 dBZ from all tilts in an annular region between two predefined ranges, r1 and r2 (with the empirical values of 20 and 80 km, resp.), are used to obtain an averaged VPR [39]. The number of data samples within a 20-km radius of the radar is limited, owing to a terrain clearance rule incorporated into a quality control (QC) scheme [40], whereby any radar bin whose bottom is within 50 m of the ground is removed. Furthermore, the region of interest has to be far enough from the radar receiver, to avoid the cone of silence in the receiver’s immediate vicinity. The accuracy of precipitation classification is also limited by the decrease in vertical resolution beyond an 80-km radius.The second step combines the output of the first step with the method of bright band identification proposed by Zhang et al. [39], to determine whether a bright band is present. The specific method adopted for this step can be summarized as follows:(1)
The local maximum reflectivity (Zpeak) and its height (hpeak) in the VPR are found. The search for the local maximum in the VPR starts at 500 m above the 0°C isotherm altitude (h_0°C) obtained from the Weather Research and Forecasting (WRF) model and continues downward. A 500-m cushion is used to account for uncertainties in the model’s 0°C height.(2)
Once the local maximum is found, the heighthtop (hbottom) above (below) the maximum level at which reflectivity decreases monotonically by a given percentage (default = 10%) [39] of the maximum reflectivity (in dBZ) is determined, along with its corresponding reflectivity Ztop (Zbottom).(3)
If the following three criteria are met, a bright band is considered to exist:(1)htop-hbottom≤D0,htop-hpeak≤D1,hpeak-hbottom≤D1.The parametersD0 and D1 are assumed to be dependent on both the vertical resolution of the reflectivity observations and the radar scan strategy. Here, D0 and D1 are set to 1.5 and 1 km, respectively.(4)
If a bright band exists, its top (BBt) and bottom (BBb) heights are defined as follows:(2)BBt=htop,htop-hpeak≤Dt,hpeak+Dt,htop-hpeak>Dt,BBb=hbottom,hpeak-hbottom≤Db,hpeak-Db,hpeak-hbottom>Db.Here, Dt and Db serve as top and bottom caps of the bright band (BB), with default values of 500 and 700 m, respectively [39].To avoid excessive corrections of the bright band, the height of the local maximum reflectivity is searched again if the minimum reflectivity is lower than a preset threshold (default = 28 dBZ). This height is defined as the lowest height at which the reflectivity is larger than or equal to the threshold [38]. This parameter can sometimes lead to insufficient or excessive bright band correction, an issue that will be solved in the future using observations from a dual-polarization radar [41].The third step involves the correction of any existing bright bands [42]. If the values of BBtBBb and their corresponding reflectivities (ZBt and ZBb, resp.) are valid, BBb<hpeak, and BBt>hpeak, and then the slope (S) of reflectivity between BBt and BBb is calculated as follows:(3)S=ZBt-ZBbBBt-BBb.If the reflectivity heighthi is located between BBt and BBb, the reflectivity is reassigned as(4)Zi=ZBt-S×BBt-hi.IfBBt does not exist, the value of BBb is valid, and BBb<hpeak, and the reflectivities located above BBb (with value larger than ZBb) are reassigned to ZBb.The fourth step, which involves reflectivity observations from all tilts in an annular region between 20 and 80 km, involves recalculation of the average stratiform VPR after bright band correction.The fifth step involves determining whether a tropical VPR has been captured. Following Xu et al. [33], the identification (or diagnosis) of tropical VPR is made with a least-squares fit from the bottom of the bright band to the bottom of the VPR, in cases where a bright band exists; otherwise, the fit is made fromh_0°C to the bottom of the VPR. The slope (α) of the VPR is then calculated. The resulting VPR is considered as tropical if three out of five values of α (for the present and the four previous moments) are lower than or equal to 0; otherwise, it is regarded as a nontropical VPR.


## 2.2. Convective and Stratiform Rainfall Classification

To improve the accuracy of radar QPE, this study adopted the FL algorithm developed by Yang et al. [30] to classify convective and stratiform rainfall. Brief details of the FL algorithm are given below.First, the volume scan reflectivity data from the Doppler radar are interpolated to a Cartesian coordinate system with the same latitude and longitude range (0.01°) of the original data. Four features, denoted asF1–F4, are calculated based on the reflectivity data. These features are, respectively, (a) the reflectivity at 2 km, (b) the reflectivity standard deviation in the horizontal direction, (c) the product of the radar top height and the reflectivity value at 2 km, and (d) the vertically integrated liquid water content [30].As a second step, membership functions are used to determine the degree to which each feature belongs to each rain type, using a fuzzification process. Linear functions are used as membership functions for convection:(5)μk,Cx,a,b=0,x≤a,x-ab-a,a<x<b,1,x≥b,k=1,2,…,Nand stratification:(6)μk,Sx,a,b=1,x≤a,b-xb-a,a<x<b,0,x≥b,k=1,2,…,N,where N is the input feature parameter, subscript k is set to 1, 2, 3, or 4, C and S identify the convective and stratiform cases, respectively, x is the feature value, a is the left breaking point, and b is the right breaking point. Different feature parameters are assigned different breaking points. ForF1, the parameters are set to a=20 dBZ and b=45 dBZ; forF2, they area = 1 dBZ and b=14 dBZ; for F3, they are a = 100 km·dBZ andb = 500 km·dBZ; and, for F4, they area = 0.5 kg/km2 and b=5.0 kg/km2 [30]. From (5) and (6), it is clear that(7)μk,Sx,a,b=1-μk,Cx,a,b.Third, the weighted average values of the measurement belonging to each specific class are obtained as follows:(8)Pe=∑k=1NWk·Pk,e∑k=1NWk,where Pk,e=μk,ε(Xk), e identifies the convective or stratiform case, Wk is the weighting factor (with equal weighting assigned in this study for simplicity), N is the number of input feature parameters, and Pc and Ps indicate the weighted average values of the convective and stratiform classes, respectively. As the sum of Pc and Ps is 1, the value of Pc can be used to express classifications for convection in a probabilistic manner. Grid points for which Pc≥0.5 are classified as convective rainfall; all other grid points are classified as stratiform rainfall.


## 2.3. Identification of Tropical Rainfall Regions

To improve the accuracy of the radar QPE, tropical rainfall is considered to be present if tropical VPR is present and the grid at 1 km altitude meets the following two criteria [33]: (1) the reflectivity is greater than 25 dBZ and (2) the surface wet bulb temperature is greater than 2°C. The surface wet bulb temperature is calculated as follows:(9)E=6.11×107.5×Tdc/237.7+TdcWBc=0.00066×P×Tc+4098×E/Tdc×237.7+Tdc20.00066×P+4089×E/237.7+Tdc2,where Tdc is the dew point temperature (in °C),E is the vapor pressure from dew point (in mb), WBc is the surface wet bulb temperature (in °C), P is the pressure (in mb), and Tc is the air temperature (in °C). The above parameters are simulated using the WRF model or from other data.In this way, stratiform rainfall is further divided into stratiform and tropical-stratiform rainfall, and convective rainfall is further divided into convective and tropical-convective rainfall.


## 2.4. Radar QPE Calculation

Three sets of experiments were conducted to test the proposed method. In Experiment I, only the WSR-88D defaultZ-R relationship (Z = 300R1.4) was used. In Experiment II, as the rain type is identified by a FL algorithm, the correspondingZ-R relationships were applied to the rain rate calculation at each grid point. For each grid point, the relationship Z = 300R1.4 was applied for the convective rainfall calculation if Pc≥0.5. Otherwise, the relationship Z = 200R1.6 [6] was applied, for stratiform rainfall calculation. In Experiment III, stratiform and convective rainfall in tropical rainfall events were further divided into stratiform and tropical-stratiform rainfall and convective and tropical-convective rainfall, respectively. The relationships were the same as those used in Experiment II, except that the tropical rainfall relationship Z = 30.7R1.66 was adopted (Yadong Wang, personal communication). In an attempt to exclude the impact of nonprecipitation echoes and hail, quality control was conducted prior to precipitation classification. Additionally, reflectivity at 1 km was used to reduce the errors caused by spatial inconsistencies and the effects of bright band measurements not yet properly revised.


## 2.5. Performance Evaluation

The mean deviation error (BIAS), relative absolute error (RAE) rate, and root mean square error (RMSE) are selected as statistical indicators for quantitatively assessing the quality of radar QPE in the different experiments [43]. The mean deviation error is defined as follows:(10)BIAS=1n∑i=1nRai-Rgi,where Rai is the radar QPE of the ith station, Rgi is the observed precipitation of the ith site, and n is the total number of samples involved in the performance assessment.The radar QPE error may be positive or negative; therefore, its mean will not reflect the radar QPE performance. To objectively reflect performance, the mean absolute error rate is used and is defined as follows:(11)RAE=∑i=1nRai-Rgi∑i=1nRgi×100%.In addition, downpour and drizzle may coexist during a precipitation process; as the RAE of downpour is larger than that of drizzle, the root mean square error primarily represents the performance of the downpour component of the rainfall:(12)RMSE=∑i=1nRai-Rgi2n1/2.As the above three equations imply, a BIAS closer to 0 corresponds to more similar (in the average) QPE and observed rain types, whereas a smaller RAE corresponds to a higher radar QPE accuracy for the entire precipitation process. Similarly, a smaller RMSE indicates a more accurate radar QPE for the heavy precipitation portion of the process.


## 3. Long-Term Data Statistics

Many previous studies have shown that extreme precipitation is likely to take place in the Yangtze River–Huaihe River valley, the Yangtze River valley, and on the southeast coast of China, with persistent, heavy rainfall events being more abundant in June and July [44–46]. Tropical rain has been found in these regions, although underestimation occurs when conventional QPE is utilized. We evaluated three months of data obtained from the Yangtze River–Huaihe River valley, in order to assess the degree to which the proposed identification process improves estimation accuracy when compared with analyses conducted using the default WSR-88DZ-R relationship and convective-stratiform precipitation segmentation rain identification scheme. Doppler radar data for this study were taken from the Hefei Doppler radar (117.258°E, 31.867°N, 165.5 m), a Chinese New Generation Radar S-band radar instrument (CINRAD WSR-98D/SA). The WSR-98D/SA is a 10-cm wavelength Doppler radar with a 1° half-power beam width. The data obtained consisted of volume scans of radar reflectivity, radial velocity, and spectrum width collected in a polar coordinate system, at increasing elevation angles. During periods of precipitation, the radar operates in a 360° azimuthal volume scan mode, with the elevation angle increasing from 0.5° to 19.5° and the temporal resolution of the data depending on the operational mode of the radar. Bin spacing is 250 m in the radial direction, with reflectivity values averaged over four bins, to increase the number of independent measurements collected for each recorded value. Accordingly, reflectivity values are recorded at 1-km intervals along the radar beam, whereas velocity parameters are recorded at 250 m intervals. Each volume scan takes approximately 6 min. The radar parameters are shown in Table 1.Table 1
Hefei CINRAD radar parameters.
Wavelength Beam width Scan mode Elevation angle Bin spacing Obtained data Time resolution 10 cm 1° half-power beam width 360° azimuthal volume 0.5°–19.5° 250 m for velocity and spectral width; 1000 m for reflectivity Reflectivity, radial velocity, spectral width Approximately 6 min


### 3.1. Data and Domain

The data from the Hefei Doppler radar and gauge of Anhui Province were collected from June to August in 2010. The horizontal domain (115.758°E–118.758°E, 30.367°N–33.367°N) in this study was a 300 × 300 km grid centered at the radar with a 1 km horizontal resolution (a 301 × 301 grid). The vertical domain consisted of 73 layers and had an altitude of 18 km (0.25-km resolution). The center of the horizontal domain is the site of Hefei radar. The domain is shown in Figure2.Figure 2
Study domain (the “+” represents the Hefei radar location).


### 3.2. Statistical Results for Three Months of Data

To compare the estimates obtained with different schemes, some statistical indicators were obtained for the June–August, 2010, period, and are presented in Table2. As shown, Experiment II performed consistently better than Experiment I, but Experiment III produced the results closest to the observed precipitation and with the greatest precision. It should be mentioned that there is a significant difference between the observed precipitation and the radar estimated rainfall when the intensity of the observed precipitation is less than 2 mm/h; therefore, rainfall with intensity below 2 mm/h was excluded.Table 2
Error analysis of the June–August, 2010, data.
Error Exp. I Exp. II Exp. III BIAS (mm) −1.299 −1.124 −0.219 RAE (%) 31.0 28.3 27.1 RMSE (mm) 2.191 2.100 1.877The statistical indicators in Table2 were obtained from all precipitation types; to further understand the particular effects of tropical rain, the statistical indicators specific of tropical rain events that occurred in the June–August, 2010, period are presented in Table 3. In other words, in contrast with Table 2, the statistical indicators in Table 3 do not consider data from grid points not experiencing tropical rainfall. The results in Table 3 match those of Table 2 in that radar QPE is noticeably improved using the FL method proposed here. This improvement becomes even more significant when in the presence of tropical rainfall events.Table 3
Error analysis of stations with tropical rain in the June–August 2010 period.
Error Exp. I Exp. II Exp. III BIAS (mm) −1.912 −1.791 0.041 RAE (%) 32.3 31.0 22.7 RMSE (mm) 3.354 3.314 2.327To further understand the various precipitation events that took place in June–August, 2010, and better describe how QPE is affected under three rainfall conditions (stratiform, convective, and a hybrid precipitation case in which convective and stratiform precipitation are present at the same time), the three-month period was divided into 12 precipitation cases (Table4). The analysis of the several events is shown in Figure 5, which presents the error curves for each case. As shown, Experiment III consistently produced the lowest error values (with the exception of case 3), particularly in cases 1, 6, and 7.Table 4
Twelve precipitation events selected for performance evaluation.
Event Precipitation duration Binary classification Tropical rain present? 1 June 07–09 Stratiform Yes 2 June 13-14 Stratiform Yes 3 June 17–20 Convective Yes 4 June 23-24 Stratiform Yes 5 June 27–30 Stratiform Yes 6 July 02–04 Convective and stratiform Yes 7 July 06–19 Convective and stratiform Yes 8 July 21-22 Convective and stratiform Yes 9 July 24–August 02 Convective and stratiform Yes 10 August 04–11 Convective and stratiform Yes 11 August 14–20 Convective and stratiform Yes 12 August 21–31 Convective and stratiform YesFigure4 shows the same type of information in Figure 3, but this time considering only tropical rainfall events. Thus, as shown, separating precipitation into convective and stratiform using the FL algorithm produces better results than simply using the defaultZ-R relationship, although separating the precipitation into convective, stratiform, and tropical rain produces results even closer to the observed precipitation, not only for the overall precipitation process but also for the heavy precipitation portion only.Figure 3
Error statistics considering all stations, for each precipitation process. Label “2cls” represents classification in two precipitation classes (Experiment II), and “3cls” represents classification in three precipitation classes (Experiment III). The units of BIAS and RMSE are mm; RAE is dimensionless.
(a) (b) (c)Figure 4
Error statistics considering only tropical rain events. Label “2cls” represents classification in two precipitation classes (Experiment II), and “3cls” represents classification in three precipitation classes (Experiment III). The units of BIAS and RMSE are mm; RAE is dimensionless.
(a) (b) (c)Figure 5
Hourly errors for all stations. Cases 1, 3, 6, and 7. The abscissas represent rainfall duration (unit: h).In cases 1, 6, and 7, Experiment III produced the best results; however, its results for case 3 are not very good. These four cases were therefore selected to further study the different results of the three sets of experiments. The corresponding errors for each station are shown in Figure5, as a function of time; cases 1, 3, 6, and 7 are represented from the top to the bottom row, respectively, and the BIAS, RAE, and RMSE values are represented in columns, from left to right, respectively. It should be noted that these results are the same ones presented in Figure 3; as shown, in cases 1, 6, and 7, Experiment III not only produced the smallest total errors but also consistently exhibited the lowest errors on an hourly basis. However, the results of Experiment III are not better than those of II or I in case 3. Further analysis revealed that case 3 represents an intermittent rainfall event that was primarily formed by scattered convective cells in which either the precipitation ranges or the rainfall durations were scattered, with each rainfall lasting less than 6 h and with limited rainfall per hour. The shortcomings in case 3 may also relate to the selected reflectivity-rainfall rate (Z-R) relationships, as these were obtained by statistics related to seasonal climate and geography. Therefore, further work will have to be performed to better capture local real-time dynamicZ-R relationships in the future.To better understand the differences between gauge measurements and QPE in terms of the amount of precipitation, scatter plots of gauge versus radar-derived 6-h accumulated precipitation estimates except for case 3, for which 4-h accumulated precipitation values are used, because in this case the rainfall lasted less than 6 h are shown in Figure6. Figure 6(a) shows the results of case 1, in which the strongest precipitation process intensity is 14 mm/h. It is seen that the precipitation produced by all three experiments corresponds well with observations, in particular the rainfall estimate from Experiment III, and that most of the results approximate the ideal line. In addition, smaller rainfall intensities correspond to bigger errors, particularly at rainfall intensities less than 2 mm/h. Figure 6(b) shows the results for case 3, in which either the precipitation ranges or rainfall durations were scattered. It is clear that the radar-retrieved accumulation precipitation estimations produced in Experiment III mostly improved the problem of underestimation, although overestimation is demonstrated in some stations. Moreover, the results from Experiment II and Experiment I are similar in that convective rainfall precipitation processes are mainly obtained from the FL method classification, a result consistent with Table 4. Figure 6(c) shows the results for case 6, in which the precipitation area is broader and the strongest intensity is 70 mm/h. The correspondence between gauge and QPE improves from Experiment I to Experiment III. Figure 6(d) shows the results for case 7, in which the strongest intensity is 60 mm/h. It is clear that the precipitation for all three experiments corresponds well with observations, especially in the case of Experiment III, which shows a significantly reduced underestimation when compared with Experiments I and II.Figure 6
Scatter plot of gauge versus radar-derived 6-h accumulated precipitation estimates (case 3 uses 4-h accumulated precipitation values). (a) Case 1. (b) Case 3. (c) Case 6. (d) Case 7. (The diagonal line represents the ideal line, corresponding to an equality between the radar and gauge estimates.)
(a) (b) (c) (d)


## 3.1. Data and Domain

The data from the Hefei Doppler radar and gauge of Anhui Province were collected from June to August in 2010. The horizontal domain (115.758°E–118.758°E, 30.367°N–33.367°N) in this study was a 300 × 300 km grid centered at the radar with a 1 km horizontal resolution (a 301 × 301 grid). The vertical domain consisted of 73 layers and had an altitude of 18 km (0.25-km resolution). The center of the horizontal domain is the site of Hefei radar. The domain is shown in Figure2.Figure 2
Study domain (the “+” represents the Hefei radar location).


## 3.2. Statistical Results for Three Months of Data

To compare the estimates obtained with different schemes, some statistical indicators were obtained for the June–August, 2010, period, and are presented in Table2. As shown, Experiment II performed consistently better than Experiment I, but Experiment III produced the results closest to the observed precipitation and with the greatest precision. It should be mentioned that there is a significant difference between the observed precipitation and the radar estimated rainfall when the intensity of the observed precipitation is less than 2 mm/h; therefore, rainfall with intensity below 2 mm/h was excluded.Table 2
Error analysis of the June–August, 2010, data.
Error Exp. I Exp. II Exp. III BIAS (mm) −1.299 −1.124 −0.219 RAE (%) 31.0 28.3 27.1 RMSE (mm) 2.191 2.100 1.877The statistical indicators in Table2 were obtained from all precipitation types; to further understand the particular effects of tropical rain, the statistical indicators specific of tropical rain events that occurred in the June–August, 2010, period are presented in Table 3. In other words, in contrast with Table 2, the statistical indicators in Table 3 do not consider data from grid points not experiencing tropical rainfall. The results in Table 3 match those of Table 2 in that radar QPE is noticeably improved using the FL method proposed here. This improvement becomes even more significant when in the presence of tropical rainfall events.Table 3
Error analysis of stations with tropical rain in the June–August 2010 period.
Error Exp. I Exp. II Exp. III BIAS (mm) −1.912 −1.791 0.041 RAE (%) 32.3 31.0 22.7 RMSE (mm) 3.354 3.314 2.327To further understand the various precipitation events that took place in June–August, 2010, and better describe how QPE is affected under three rainfall conditions (stratiform, convective, and a hybrid precipitation case in which convective and stratiform precipitation are present at the same time), the three-month period was divided into 12 precipitation cases (Table4). The analysis of the several events is shown in Figure 5, which presents the error curves for each case. As shown, Experiment III consistently produced the lowest error values (with the exception of case 3), particularly in cases 1, 6, and 7.Table 4
Twelve precipitation events selected for performance evaluation.
Event Precipitation duration Binary classification Tropical rain present? 1 June 07–09 Stratiform Yes 2 June 13-14 Stratiform Yes 3 June 17–20 Convective Yes 4 June 23-24 Stratiform Yes 5 June 27–30 Stratiform Yes 6 July 02–04 Convective and stratiform Yes 7 July 06–19 Convective and stratiform Yes 8 July 21-22 Convective and stratiform Yes 9 July 24–August 02 Convective and stratiform Yes 10 August 04–11 Convective and stratiform Yes 11 August 14–20 Convective and stratiform Yes 12 August 21–31 Convective and stratiform YesFigure4 shows the same type of information in Figure 3, but this time considering only tropical rainfall events. Thus, as shown, separating precipitation into convective and stratiform using the FL algorithm produces better results than simply using the defaultZ-R relationship, although separating the precipitation into convective, stratiform, and tropical rain produces results even closer to the observed precipitation, not only for the overall precipitation process but also for the heavy precipitation portion only.Figure 3
Error statistics considering all stations, for each precipitation process. Label “2cls” represents classification in two precipitation classes (Experiment II), and “3cls” represents classification in three precipitation classes (Experiment III). The units of BIAS and RMSE are mm; RAE is dimensionless.
(a) (b) (c)Figure 4
Error statistics considering only tropical rain events. Label “2cls” represents classification in two precipitation classes (Experiment II), and “3cls” represents classification in three precipitation classes (Experiment III). The units of BIAS and RMSE are mm; RAE is dimensionless.
(a) (b) (c)Figure 5
Hourly errors for all stations. Cases 1, 3, 6, and 7. The abscissas represent rainfall duration (unit: h).In cases 1, 6, and 7, Experiment III produced the best results; however, its results for case 3 are not very good. These four cases were therefore selected to further study the different results of the three sets of experiments. The corresponding errors for each station are shown in Figure5, as a function of time; cases 1, 3, 6, and 7 are represented from the top to the bottom row, respectively, and the BIAS, RAE, and RMSE values are represented in columns, from left to right, respectively. It should be noted that these results are the same ones presented in Figure 3; as shown, in cases 1, 6, and 7, Experiment III not only produced the smallest total errors but also consistently exhibited the lowest errors on an hourly basis. However, the results of Experiment III are not better than those of II or I in case 3. Further analysis revealed that case 3 represents an intermittent rainfall event that was primarily formed by scattered convective cells in which either the precipitation ranges or the rainfall durations were scattered, with each rainfall lasting less than 6 h and with limited rainfall per hour. The shortcomings in case 3 may also relate to the selected reflectivity-rainfall rate (Z-R) relationships, as these were obtained by statistics related to seasonal climate and geography. Therefore, further work will have to be performed to better capture local real-time dynamicZ-R relationships in the future.To better understand the differences between gauge measurements and QPE in terms of the amount of precipitation, scatter plots of gauge versus radar-derived 6-h accumulated precipitation estimates except for case 3, for which 4-h accumulated precipitation values are used, because in this case the rainfall lasted less than 6 h are shown in Figure6. Figure 6(a) shows the results of case 1, in which the strongest precipitation process intensity is 14 mm/h. It is seen that the precipitation produced by all three experiments corresponds well with observations, in particular the rainfall estimate from Experiment III, and that most of the results approximate the ideal line. In addition, smaller rainfall intensities correspond to bigger errors, particularly at rainfall intensities less than 2 mm/h. Figure 6(b) shows the results for case 3, in which either the precipitation ranges or rainfall durations were scattered. It is clear that the radar-retrieved accumulation precipitation estimations produced in Experiment III mostly improved the problem of underestimation, although overestimation is demonstrated in some stations. Moreover, the results from Experiment II and Experiment I are similar in that convective rainfall precipitation processes are mainly obtained from the FL method classification, a result consistent with Table 4. Figure 6(c) shows the results for case 6, in which the precipitation area is broader and the strongest intensity is 70 mm/h. The correspondence between gauge and QPE improves from Experiment I to Experiment III. Figure 6(d) shows the results for case 7, in which the strongest intensity is 60 mm/h. It is clear that the precipitation for all three experiments corresponds well with observations, especially in the case of Experiment III, which shows a significantly reduced underestimation when compared with Experiments I and II.Figure 6
Scatter plot of gauge versus radar-derived 6-h accumulated precipitation estimates (case 3 uses 4-h accumulated precipitation values). (a) Case 1. (b) Case 3. (c) Case 6. (d) Case 7. (The diagonal line represents the ideal line, corresponding to an equality between the radar and gauge estimates.)
(a) (b) (c) (d)


## 4. Case Studies

The cases of heavy rainfall events in Hefei (valid at 0800 UTC, July 22, 2009) and Liuzhou (valid at 1300 UTC, June 12, 2008) were analyzed, to test how easily the proposed method could be applied in different places and times. The Liuzhou Doppler radar (109.456°E, 24.357°N, 346.8 m) is of the same type as the Hefei Doppler radar.


### 4.1. Heavy Rainfall in Hefei on July 21–23, 2009

Most parts of the Yangtze River–Huaihe River valley and the southern part of the Yangtze River experienced heavy rainfall on July 21–23, 2009 [25], with local downpours of approximately 150 mm of rain. During this period, there was a wide range of precipitation, including heavy rain that was distributed primarily in the north of Anhui Province, with the precipitation intensity reaching 40 mm/h.The distributions of convective and stratiform rainfall based on the classification of the FL algorithm are shown in Figure7(a), which shows a rainfall process overall dominated by widespread stratiform rainfall embedded with convective rainfall. A stratiform cloud covering a large area appeared throughout the entire precipitation process and was relatively stable. Only a few parts of Hefei could be identified as convective rainfall regions, embedded in a wide region of stratiform rainfall. The convective rainfall regions changed frequently, and their coverage area was very limited. Thus, the range of estimated precipitation is relatively large.Figure 7
Hefei heavy rain (valid at 0800 UTC, July 22, 2009). (a) Results of the FL algorithm (the abscissas represent longitude, the vertical axis represents latitude, the “+” sign represents the radar location, and the colors of “Conv” and “Stra” represent convective and stratiform rainfall regions, resp.). (b) Distribution of the convective, stratiform, and tropical rainfall. (c) Averaged stratiform VPR (the abscissas represent the value of reflectivity in dBZ, the vertical axis represents the altitude in km, the long dashed line is the 0°C isotherm altitude, the red solid line is the VPR before bright band correction, and the blue short dashed line is the VPR after bright band correction). (d) Scatter plot of gauge versus radar-derived 1-h accumulated precipitation (the abscissas represent gauge results in mm, the vertical axis represents radar-derived 1-h accumulated precipitation in mm, and the diagonal line represents the ideal line).
(a) (b) (c) (d)Although the problem of underestimation is mitigated, to some extent, by classifying the precipitation into convective and stratiform precipitation, significant differences between observation and estimation remain; the primary objective of classifying precipitation into convective, stratiform, and tropical precipitation is the improvement of QPE accuracy. Figure7(b) shows the results obtained for the three types of cloud classification. It is clear that some convective and stratiform precipitation obtained from the FL algorithms are again identified as tropical precipitation, with the corresponding cases redefined as a stratiform, convective, or tropical mixed precipitation.As can be seen from Figure7(c), the characteristics of the precipitation event in Hefei on July 21–23, 2009, include the presence of a bright band; it was identified as convective and stratiform mixed precipitation using the FL algorithm, and the VPR slope after bright band correction was found to be negative; in other words, tropical rainfall may have been present in this case.To demonstrate the differences between the three experiments’ results, a scatter plot of gauge versus radar-derived 1-h accumulated precipitation is shown in Figure7(d). Even though all three experiments produce some degree of underestimation, the correspondence between gauge and radar-derived 1-h accumulated precipitation is improved from Experiment I (greatest underestimation) to Experiment III (least underestimation). Thus, Experiment III produces a significantly more accurate estimation of rainfall compared to the other two experiments, particularly when the rainfall intensity is high.


### 4.2. Liuzhou Rainfall Event during June 11-12, 2008

The primary difference between the Liuzhou precipitation event and the Hefei event is that the Liuzhou precipitation case took place at a different location.Liuzhou, in the Guangxi Province, received heavy rain at 0800 UTC on June 11, 2008. The 12-h rainfall accumulation from 0400 to 1600 UTC on June 12 reached 233 mm, with heavy rainfall (307 mm) also occurring in Liujiang County. The water level in Liujiang exceeded the warning level of 2.86 m until 1900 UTC on June 12. The heavy rainfall caused waterlogging in some areas and a number of road traffic interruptions until June 13, when the rain began to ease.Figure8 shows an analysis of the heavy rainfall in Liuzhou during June 11-12, 2008, valid at 1300 UTC on June 12. Comparing the 1-h rainfall accumulation from the gauges with the convective and stratiform distribution obtained from the FL algorithm (Figure 8(a)), it is clear that the entire process was based primarily on convective rainfall that did not aggregate but was found in the vicinity of the radar. Heavy rainfall reached 25 mm/h just next to the radar.Figure 8
Heavy rainfall in Liuzhou during June 11-12, 2008, as recorded at (or analyzed for) 1300 UTC on June 12: (a) results of the FL algorithm (the abscissas represent longitude, the vertical axis represents latitude, the “+” sign represents the radar location, and colors of “Conv” and “Stra” represent convective and stratiform rainfall, resp.). (b) Distribution of convective, stratiform, and tropical rainfall. (c) Averaged stratiform VPR (the abscissas represent the value of reflectivity in dBZ, the vertical axis represents the altitude in km, the long dashed line is the 0°C isotherm altitude, the red solid line is the VPR before bright band correction. (d) Scatter plot of gauge versus radar-derived 1-h accumulated precipitation (the abscissas represent gauge results in mm, the vertical axis represents radar-derived 1-h accumulated precipitation in mm, and the diagonal line represents the ideal line).
(a) (b) (c) (d)The distributions of the three types of rainfall are illustrated in Figure8(b), which shows that the Liuzhou event was a mixed stratiform, convective, and tropical precipitation case. Most of the convective and stratiform precipitation obtained from the FL algorithm is again distinguishable as tropical rain, although there is some convective precipitation near the radar, which recognized the area as being dominated by convective precipitation.As the VPR shown in Figure8(c) implies, no bright band is observed for this event. The value ofh_0°C is 5042.6 m, the VPR slope is negative, and tropical precipitation is possibly present.For comparison between the three experiments, a scatter plot of gauge versus radar-derived 1-h accumulated precipitation is provided in Figure8(d). It is clear that the precipitation of all three experiments corresponds well with observations. Experiment II performs better than Experiment I, and Experiment III performs the best of all, particularly for large rainfall intensity. Thus, the rainfall estimated by Experiment III is found to correspond most closely to the observed rainfall.


## 4.1. Heavy Rainfall in Hefei on July 21–23, 2009

Most parts of the Yangtze River–Huaihe River valley and the southern part of the Yangtze River experienced heavy rainfall on July 21–23, 2009 [25], with local downpours of approximately 150 mm of rain. During this period, there was a wide range of precipitation, including heavy rain that was distributed primarily in the north of Anhui Province, with the precipitation intensity reaching 40 mm/h.The distributions of convective and stratiform rainfall based on the classification of the FL algorithm are shown in Figure7(a), which shows a rainfall process overall dominated by widespread stratiform rainfall embedded with convective rainfall. A stratiform cloud covering a large area appeared throughout the entire precipitation process and was relatively stable. Only a few parts of Hefei could be identified as convective rainfall regions, embedded in a wide region of stratiform rainfall. The convective rainfall regions changed frequently, and their coverage area was very limited. Thus, the range of estimated precipitation is relatively large.Figure 7
Hefei heavy rain (valid at 0800 UTC, July 22, 2009). (a) Results of the FL algorithm (the abscissas represent longitude, the vertical axis represents latitude, the “+” sign represents the radar location, and the colors of “Conv” and “Stra” represent convective and stratiform rainfall regions, resp.). (b) Distribution of the convective, stratiform, and tropical rainfall. (c) Averaged stratiform VPR (the abscissas represent the value of reflectivity in dBZ, the vertical axis represents the altitude in km, the long dashed line is the 0°C isotherm altitude, the red solid line is the VPR before bright band correction, and the blue short dashed line is the VPR after bright band correction). (d) Scatter plot of gauge versus radar-derived 1-h accumulated precipitation (the abscissas represent gauge results in mm, the vertical axis represents radar-derived 1-h accumulated precipitation in mm, and the diagonal line represents the ideal line).
(a) (b) (c) (d)Although the problem of underestimation is mitigated, to some extent, by classifying the precipitation into convective and stratiform precipitation, significant differences between observation and estimation remain; the primary objective of classifying precipitation into convective, stratiform, and tropical precipitation is the improvement of QPE accuracy. Figure7(b) shows the results obtained for the three types of cloud classification. It is clear that some convective and stratiform precipitation obtained from the FL algorithms are again identified as tropical precipitation, with the corresponding cases redefined as a stratiform, convective, or tropical mixed precipitation.As can be seen from Figure7(c), the characteristics of the precipitation event in Hefei on July 21–23, 2009, include the presence of a bright band; it was identified as convective and stratiform mixed precipitation using the FL algorithm, and the VPR slope after bright band correction was found to be negative; in other words, tropical rainfall may have been present in this case.To demonstrate the differences between the three experiments’ results, a scatter plot of gauge versus radar-derived 1-h accumulated precipitation is shown in Figure7(d). Even though all three experiments produce some degree of underestimation, the correspondence between gauge and radar-derived 1-h accumulated precipitation is improved from Experiment I (greatest underestimation) to Experiment III (least underestimation). Thus, Experiment III produces a significantly more accurate estimation of rainfall compared to the other two experiments, particularly when the rainfall intensity is high.


## 4.2. Liuzhou Rainfall Event during June 11-12, 2008

The primary difference between the Liuzhou precipitation event and the Hefei event is that the Liuzhou precipitation case took place at a different location.Liuzhou, in the Guangxi Province, received heavy rain at 0800 UTC on June 11, 2008. The 12-h rainfall accumulation from 0400 to 1600 UTC on June 12 reached 233 mm, with heavy rainfall (307 mm) also occurring in Liujiang County. The water level in Liujiang exceeded the warning level of 2.86 m until 1900 UTC on June 12. The heavy rainfall caused waterlogging in some areas and a number of road traffic interruptions until June 13, when the rain began to ease.Figure8 shows an analysis of the heavy rainfall in Liuzhou during June 11-12, 2008, valid at 1300 UTC on June 12. Comparing the 1-h rainfall accumulation from the gauges with the convective and stratiform distribution obtained from the FL algorithm (Figure 8(a)), it is clear that the entire process was based primarily on convective rainfall that did not aggregate but was found in the vicinity of the radar. Heavy rainfall reached 25 mm/h just next to the radar.Figure 8
Heavy rainfall in Liuzhou during June 11-12, 2008, as recorded at (or analyzed for) 1300 UTC on June 12: (a) results of the FL algorithm (the abscissas represent longitude, the vertical axis represents latitude, the “+” sign represents the radar location, and colors of “Conv” and “Stra” represent convective and stratiform rainfall, resp.). (b) Distribution of convective, stratiform, and tropical rainfall. (c) Averaged stratiform VPR (the abscissas represent the value of reflectivity in dBZ, the vertical axis represents the altitude in km, the long dashed line is the 0°C isotherm altitude, the red solid line is the VPR before bright band correction. (d) Scatter plot of gauge versus radar-derived 1-h accumulated precipitation (the abscissas represent gauge results in mm, the vertical axis represents radar-derived 1-h accumulated precipitation in mm, and the diagonal line represents the ideal line).
(a) (b) (c) (d)The distributions of the three types of rainfall are illustrated in Figure8(b), which shows that the Liuzhou event was a mixed stratiform, convective, and tropical precipitation case. Most of the convective and stratiform precipitation obtained from the FL algorithm is again distinguishable as tropical rain, although there is some convective precipitation near the radar, which recognized the area as being dominated by convective precipitation.As the VPR shown in Figure8(c) implies, no bright band is observed for this event. The value ofh_0°C is 5042.6 m, the VPR slope is negative, and tropical precipitation is possibly present.For comparison between the three experiments, a scatter plot of gauge versus radar-derived 1-h accumulated precipitation is provided in Figure8(d). It is clear that the precipitation of all three experiments corresponds well with observations. Experiment II performs better than Experiment I, and Experiment III performs the best of all, particularly for large rainfall intensity. Thus, the rainfall estimated by Experiment III is found to correspond most closely to the observed rainfall.


## 5. Conclusions

This paper introduced a method for improving the precision of radar QPE by using an FL algorithm to divide events into stratiform and convective precipitation types and then further identifying tropical rain based on the VPR characteristics. First, the proposed method corrects any bright bands that are present. Volume scan reflectivity data from Doppler radar are used to recognize tropical VPRs, using the precipitation classification technique of Xu et al. [33]. Volume scan reflectivity data from Doppler radar are then interpolated to the Cartesian coordinate system over the same latitude and longitude range (0.01°), and the rain is divided into stratiform precipitation and convective precipitation using an FL algorithm [30]. The precipitation is then considered to be tropical rainfall if it has a tropical VPR and the grid at an altitude of 1 km meets the following criteria: (1) the reflectivity is greater than 25 dBZ and (2) the surface wet bulb temperature is greater than 2°C. The reflectivity at 1 km is used to assign an adaptiveZ-R relationship. Three sets of experiments were presented in this paper. In Experiment I, for each weather event only the default WSR-88DZ-R relationship (Z = 300R1.4) was used. In Experiment II, the rain type was identified as convective or stratiform using an FL algorithm; then, the relationships Z = 300R1.4 and Z = 200R1.6 were applied for convective and stratiform rainfall calculations, respectively. In Experiment III, stratiform and convective rainfall in tropical rainfall events were further divided into stratiform and tropical-stratiform rainfall and convective and tropical-convective rainfall, respectively. The relationships adopted were the same as those used in Experiment II, except that the tropical rainfall relationship was Z = 30.7R1.66. A long sequence (three months) of data and two additional study cases were conducted, in order to determine the impact on QPE of the method devised for this study. The conclusions of the study can be summarized as follows:(1)
Although most radar QPEs were lower than the directly observed precipitation for all the rainfall events analyzed in this study, the results of all three experiments corresponded closely to the observed precipitation. Separating the precipitation into convective and stratiform types using the FL algorithm enabled better performances compared to simply using the defaultZ-R relationship. Separating the precipitation into convective, stratiform, and tropical rain produced results that were the closest to the observed precipitation, not only during the whole precipitation processes but also individually within the heavy precipitation portions of these events.(2)
Given that a large-scale precipitation event usually contains multiple types of precipitation, large estimation errors can be introduced if a singleZ-R relationship is used for the entire process. The FL algorithm used here can skillfully separate convective precipitation from stratiform precipitation and thus improve the accuracy of radar QPE to some extent. Nevertheless, the problem of underestimation remains. This study built upon the characteristics of tropical precipitation and the convective, stratiform precipitation classification techniques using an FL algorithm; the resulting method was shown to have improved accuracy, not only during the whole precipitation processes but also individually within the heavy precipitation portions of these events.(3)
The proposed tropical rain identification scheme produced similar results when applied to the cases of heavy rainfall in Hefei (valid at 0800 UTC on July 22, 2009) and in Liuzhou (valid at 1300 UTC on June 12, 2008). In future work, the newly proposed scheme will be further tested, to evaluate its general applicability to different times and locations.



---

*Source: 2457489-2016-11-24.xml*