Next Article in Journal
A Novel 2-D Geometry Reconstruction Approach for Space Debris via Interpolation-Free Operation under Low SNR Conditions
Previous Article in Journal
Impact of Stereo Camera Calibration to Object Accuracy in Multimedia Photogrammetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantification of Precipitation Using Polarimetric Radar Measurements during Several Typhoon Events in Southern China

1
College of Electronic Engineering, Chengdu University of Information Technology, Chengdu 610225, China
2
State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China
3
NOAA Physical Sciences Laboratory, Boulder, CO 80305, USA
4
Cooperative Institute for Research in the Atmosphere, Fort Collins, CO 80523, USA
5
National Center for Atmospheric Research, Boulder, CO 80301, USA
6
College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China
7
Department of Hydraulic Engineering, State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China
8
Institute of Tropical and Marine Meteorology, China Meteorological Administration, Guangzhou 510080, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(12), 2058; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12122058
Submission received: 10 May 2020 / Revised: 22 June 2020 / Accepted: 24 June 2020 / Published: 26 June 2020
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
Accurate quantitative precipitation estimation (QPE) during typhoon events is critical for flood warning and emergency management. Dual-polarization radar has proven to have better performance for QPE, compared to traditional single-polarization radar. However, polarimetric radar applications have not been extensively investigated in China, especially during extreme events such as typhoons, since the operational dual-polarization system upgrade only happened recently. This paper extends a polarimetric radar rainfall system for local applications during typhoons in southern China and conducts comprehensive studies about QPE and precipitation microphysics. Observations from S-band dual-polarization radar in Guangdong Province during three typhoon events in 2017 are examined to demonstrate the enhanced radar rainfall performance. The microphysical properties of hydrometeors during typhoon events are analyzed through raindrop size distribution (DSD) data and polarimetric radar measurements. The stratiform precipitation in typhoons presents lower mean raindrop diameter and lower raindrop concentration than that of the convection precipitation. The rainfall estimates from the adapted radar rainfall algorithm agree well with rainfall measurements from rain gauges. Using the rain gauge data as references, the maximum normalized mean bias ( N M B ) of the adapted radar rainfall algorithm is 20.27%; the normalized standard error ( N S E ) is less than 40%; and the Pearson’s correlation coefficient ( C C ) is higher than 0.92. For the three typhoon events combined, the N S E and N M B are 36.66% and -15.78%, respectively. Compared with several conventional radar rainfall algorithms, the adapted algorithm based on local rainfall microphysics has the best performance in southern China.

Graphical Abstract

1. Introduction

Tropical cyclones are extreme weather phenomena often accompanied by high winds and heavy rainfall, which may result in significant socioeconomic losses in coastal areas. Therefore, numerous studies have focused on radar-based quantitative precipitation estimation (QPE) and forecast (QPF) during Atlantic hurricanes, and Australian and Indian cyclones [1,2]. However, polarimetric radar QPE applications for typhoon events in China are relatively rare, partly because the operational weather radar network has just been upgraded with dual-polarization capabilities. In addition, there is a lack of comprehensive analysis of QPE and associated precipitation microphysics during typhoon events.
Compared to traditional single-polarization radar, the dual-polarization radar can provide polarimetric measurements, including differential reflectivity Z D R , differential phase Ψ D P , co-polar correlation coefficient ρ H V , and the specific differential phase shift K D P [3,4]. These polarimetric observations directly reflect the microphysical properties of hydrometeors, including precipitation particle shape and size distributions [5], and thus can be used to estimate the precipitation rate and amount [6,7,8,9,10,11,12,13]. In fact, dual-polarization radar has become an essential tool to study severe weather and mitigate the resultant natural disasters. In this paper, using S-band polarimetric measurements during three typhoon events, we develop an adapted rainfall system that includes data quality control (QC), data segmentation, hydrometeor classification, and precipitation quantification. The QC method used in this paper can effectively identify and mitigate the ground clutter and environmental noises, through differential phase and clutter processing. The data segmentation method is applied to reduce the impact of blockage and beam broadening. Depending on the distance of the target from the radar, data from different elevation angles are selected for QPE. The hydrometeor classification algorithm (HCA) [14], based on traditional fuzzy-logic approach and K-means clustering, is used to obtain classification results with the polarimetric parameters ( Z H , Z D R , ρ H V , K D P ) and environment temperature information T observed from a local sounding station. The classification results and polarimetric parameters determine the selection of rainfall relations, which are acquired through simulation based on local raindrop size distribution (DSD) data.
Observations from the upgraded S-band dual-polarization radars are used to demonstrate the QPE performance of the adapted algorithm, and the polarimetric observations can also indicate typhoon features such as typhoon eye, strong convective area and precipitation particle size distributions. Using DSD data and S-band polarimetric measurements during the three typhoon events, this paper also examines the microphysical properties of typhoon-induced precipitation. Overall, we hope to promote polarimetric radar applications during typhoon events in China, through this comprehensive study which combines QPE and precipitation microphysical properties observed by in situ instruments.
The rest of this paper is organized as follows. Section 2 describes the study domain, datasets, and the adapted polarimetric radar rainfall algorithm. The typhoon events, radar observations, and the microphysical properties of typhoon-induced precipitation are presented in Section 3. Section 4 provides the QPE results and quantitative evaluation based on surface rain gauge measurements. The main findings of this study and suggested future work are summarized in Section 5.

2. Datasets and Methodology

2.1. Study Domain

The study domain of this research is located in Guangdong Province in southern China, which has the largest number of typhoon landfalls in China, with an average of 2 to 3 each year [15]. The whole territory of Guangdong is between 20°13′–25°31′N and 109°39′–117°19′E, with high mountains in the north and west, and plains and low hills in the south (Figure 1). Climatologically, Guangdong is located in the East Asian monsoon region, with an average annual precipitation of 1300–2500 mm. Affected by monsoon and terrain, the annual precipitation in the southern region is generally higher than the northern part. Figure 1 shows the study domain and locations of the observational instruments used in this study, overlaid with the best tracks and landing locations of the three typhoons included in this study.

2.2. Radar and Rain Gauge Data

According to the strategic plan of the China Meteorological Administration, 140 dual-polarized radars will be deployed nationwide by the end of 2020. The upgraded China New Generation of Weather Radars (CINRADs) are expected to improve the forecasting and warning capabilities against severe weather such as heavy rainfall, hail, and typhoons. Eight weather radars in Guangdong completed the dual-polarization upgrade in 2016; this is the first batch of the dual-polarization upgrade. This study uses data from two of them, i.e., Guangzhou and Yangjiang stations. The Yangjiang radar is closer to the best tracks of Typhoons Hato and Pakhar, whereas the Guangzhou radar is closer to the best track of Typhoon Mawar. Therefore, the data from the Yangjiang radar during Typhoons Hato and Pakhar are used, and the data from the Guangzhou radar during Typhoon Mawar are used. As shown in Figure 1, the Yangjiang radar (21.845°N, 111.979°E) is deployed in Jishan, and the Guangzhou radar (23.004°N, 113.355°E) is deployed in Dazhengang. The two radars are deployed on the mountains, mainly to avoid beam blockages due to high-rise buildings. The model of the radars is CINRAD/SA (S stands for S-band, A refers to the manufacturer of Beijing Metstar Radar Company). The scan strategy is volume coverage pattern (VCP) 21, which includes nine elevation angles, i.e., 0.5°, 1.5°, 2.4°, 3.3°, 4.3°, 6°, 9.9°, 14.6°, and 19.5°. The volume scan update cycle is six minutes. The detailed system specifications of CINRAD/SA are shown in Table 1.
Radar data quality control is indispensable before quantitative applications. The building occlusion, ground clutters, and environmental noises can all affect the accuracy of the application products. In this paper, a compound QC method based on differential phase and clutter processing is used to identify and mitigate the non-meteorological echoes, including ground clutters [16,17,18,19]. Differential phase processing is mainly performed to smooth the differential phase, which includes phase unfolding, clutter suppression, range filtering, and least square fitting. Sample observations from the Yangjiang radar 1.5-degree which sweep at 2112 UTC, 23 August 2017, before and after data QC, are shown in Figure 2. It can be seen that most of the ground clutter is suppressed after QC. Figure 3 shows the scatter density plots of Z D R versus Z H before and after data QC for a single plan position indicator (PPI) scan at 2112 UTC, 23 August 2017. The black dots overlaid with the color plots are the simulated radar moments. This study evaluates the data QC performance through comparison with the simulated radar moments, including Z H and Z D R . The radar moments are simulated using the T-matrix method [20]. Therein, the raindrop shape model is assumed to follow the one described by Thurai et al. [21]. The temperature used in the simulation is 20 °C, which is obtained from nearby sounding data. The overlapping ratios of Z D R versus Z H are investigated, to qualitatively demonstrate the QC performance. Apparently, the data after QC are clearer and agree better with the simulated radar moments, demonstrating the effectiveness of the implemented QC method.
In addition, a data segmentation method is applied to further decrease the impact of blockage and ground clutter. For example, since the 0.5-degree beam from the Yangjiang radar is blocked by the terrain [22,23], this study uses a higher elevation scan to produce complete QPE. However, while the selection of higher elevation angle data could decrease the impact of blockage, it may increase the difference between radar estimates and ground gauges, due to the height differences. This issue is further compounded by the beam broadening of S-band radar at long ranges. In the segmentation, data from different elevation angles are selected for precipitation estimation, depending on the distance of the target from the radar. In particular, the elevation angles of 2.4 degrees, 1.5 degrees, and 0.5 degrees are selected within the ranges of 20–35 km, 35–60 km, and 60–100 km, respectively, from the radar.
The rainfall estimates based on the dual-polarization radars are compared with rain gauge measurements, to demonstrate the performance of radar QPE. The location of the rain gauges is shown in Figure 1. There are 198 and 809 rain gauges within 100 km of the Yangjiang and Guangzhou radars, respectively. These rain gauges record rainfall rate at hourly scales.

2.3. Adapted Algorithm

This study uses five radar QPE algorithms (Figure 4), including the adapted algorithm, to estimate rainfall during these typhoon events and cross compare the performance of different algorithms. All algorithms are implemented after the data QC and segmentation processing. The HCA, only included in the adapted algorithm, is essentially based on the methodology detailed in Bechini and Chandrasekar [14]. Figure 5a shows the flowchart of the HCA-driven adapted rainfall algorithm, which selects different rainfall relations, based on the hydrometeor classification results and polarimetric parameters. Other algorithms (Figure 5b,c) without HCA compute rainfall rate through incorporating the co-polar correlation coefficient ρ H V .
The specific rainfall relations of the adapted algorithm with and without HCA are as follows:
R ( Z H ) = 0 . 02   ×   10 0 . 07729 Z H
R ( K D P ) = 51 . 64 K D P 0 . 7226
R ( Z H , Z D R ) = 0 . 01172   ×   10 0 . 0925 Z H 10 0 . 4035 Z D R
R ( Z D R , K D P ) = 75 . 08 K D P 0 . 9191 10 0 . 1861 Z D R
where Z H (dBZ) is the horizontal reflectivity factor and Z D R (dB) is the differential reflectivity. The DSD data collected by a disdrometer in Yangjiang, Guangdong, during April to July 2014, are used to obtain these specific rainfall relations through nonlinear regression. The estimates are computed with these specific rainfall relations, and then the rainfall relation coefficients are adjusted, with an overall bias ratio computed from radar estimates and corresponding rain gauge measurements. The Yangjiang blended rainfall relations are as follows:
R ( Z H ) = 0 . 0093   ×   10 0 . 07386 Z H
R ( K D P ) = 31 . 2195 K D P 0 . 9143
R ( Z H , Z D R ) = 0 . 024   ×   10 0 . 08242 Z H 10 0 . 4913 Z D R
R ( Z D R , K D P ) = 59 . 2626 K D P 0 . 9573 10 0 . 1317 Z D R
The rainfall relations of the blended Yangjiang algorithm are derived based on the same DSD data, but these rainfall relations are not adjusted. The Yangjiang Z–R relation only contains R( Z H ) (Equation (2a)), and the WSR-88D Z–R relation is as follows:
R ( Z H ) = 0 . 017   ×   10 0 . 0714 Z H
This relation is generally applied for convective rain [24].
The Yangjiang and Guangzhou radars collect volume observations every six minutes, whereas the rain gauges record the rainfall rate at hourly time intervals. Therefore, the rainfall estimates computed by various algorithms need to be interpolated to match the rain gauge measurements. The radar rainfall rates are calculated using piecewise cubic Hermite interpolating polynomial (PCHIP) into 1-min resolution and then used to compute the hourly rainfall rate to compare with corresponding rain gauge measurements.

3. Description of the Typhoon Events

3.1. Super Typhoon Hato (2017)

Super Typhoon Hato formed in the northwest Pacific, 1200 km southeast of Shantou (20.4°N, 128.0°E) at 1800 UTC, 19 August 2017. While it moved northwest toward Guangdong Province, it reached severe tropical storm strength, with a maximum wind speed of 25 m s−1 at 0000 UTC 22 August. It then landed in Zhuhai, Guangdong Province at 0450 UTC on 23 August with severe typhoon strength (the maximum wind speed of 45 m s−1). It continued to move westward through Jiangmen, Yangjiang, Maoming, and Yunfu after the landfall. Hato weakened into a tropical storm, with a maximum wind speed of 18 m s−1 at 0000 UTC on 24 August, and further depressed into a low pressure at 0600 UTC on 24 August. Hato was the strongest typhoon that caused landfall in China in 2017. Its landfall brought widespread heavy rainfall to western Guangdong. According to the China National Meteorological Center (NMC), a total of 17 rain gauges recorded a cumulative rainfall of more than 250 mm and 243 rain gauges recorded a cumulative rainfall of 100–250 mm from 0000 UTC 23 August to 0800 UTC 24 August 2017.
Figure 6 shows sample observations of Z H , Z D R , ρ H V , Ψ D P , and K D P from the Yangjiang radar at 0430 UTC on 23 August. A rainfall accumulation map is also shown for illustration purposes, which is computed by the adapted algorithm based on 1.5-degree scans from 0000 UTC to 2300 UTC, 23 August 2017. It can be clearly seen that the center of the typhoon is 150 km east of the radar, and there is a strong convective area with reflectivity ( Z H ) of over 40 dBZ south of the typhoon center (Figure 6a). The ρ H V value of this area is less than 0.95, indicating the occurrence of rain and hail near the center of typhoon (Figure 6c). The Ψ D P and K D P near the center of the typhoon are obviously greater than the other regions (Figure 6d,e).
In addition, the rainfall accumulation computed by the adapted algorithm in most of the region within a range of 200 km of the radar exceeds 900 mm (Figure 6f). The time-series comparison of radar moments computed by DSD data with T-matrix method and measured polarimetric parameters are shown in Figure 7. It shows that precipitation in the position of the disdrometer mainly occurred from 0000 UTC to 1100 UTC and 2100 UTC to 2300 UTC 23 August 2017. The measured Z H , Z D R and K D P   are close to the values of simulated radar moments during this typhoon event, which indicates the reliability of the radar data.

3.2. Severe Tropical Storm Pakhar (2017)

Severe tropical storm Pakhar formed at 0600 UTC on 24 August 2017 in the northwestern Pacific, 310 km to the north–east of Manila (15.4°N, 123.8°E), with a central pressure of 995 hPa. After crossing the Philippines and entering the South China Sea, it further intensified into a severe tropical storm and continued to move northwestward. Pakhar caused landfall in Taishan, Guangdong Province at 0100 UTC on 27 August, with a maximum wind speed of 33 m s−1. After the landfall, it moved northwestward and passed through Jiangmen, Yangjiang, and Yunfu. Pakhar further weakened into a tropical depression in Guigang at 1200 UTC on 27 August. According to the NMC, from 1200 UTC 26 to 0000 UTC 29 August, the 24-hour average precipitation of Guangdong Province was 72.1 mm. A total of 44 rain gauges recorded a cumulative rainfall of > 250 mm, and 642 rain gauges recorded a cumulative rainfall of 100–250 mm. Figure 8 shows the sample observations of Z H , Z D R , ρ H V , Ψ D P , and K D P from the Yangjiang radar 1.5-degree sweep at 0012 UTC on 27 August. Similarly, the rainfall accumulation from 0000 UTC to 1200 UTC, 27 August 2017, is illustrated, which is based on the adapted algorithm. The typhoon center is 150 km east of the radar, and the value of Z H exceeds 40 dBZ. The value of Z D R at 75 km northeast of the radar exceeds 2 dB and the value of the co-polar correlation coefficient exceeds 0.95 (Figure 8), indicating the large size of raindrops in this area. The region of heavy precipitation is located in the north of the Yangjiang radar. However, there was no rain at the disdrometer during this typhoon event.

3.3. Severe Tropical Storm Mawar (2017)

Severe tropical storm Mawar developed as a tropical depression about 690 km to the south–east of Shenzhen (18.7°N, 119.3°E) at 0000 UTC, 31 August 2017, in the northwest Pacific. First taking a northwestward course over the north of the South China Sea, Mawar became a tropical storm, with a maximum wind speed of 25 m s−1. It then landed in Shanwei, Guangdong Province at 1330 UTC on 3 September, as a tropical storm with a maximum wind speed of 23 m s−1. After the landfall, it weakened into a low pressure at 0200 UTC on 4 September. Mawar was the third typhoon that caused landfall in Guangdong Province in the two weeks during late August 2017. According to the NMC, from 0000 UTC on 3 September to 0000 UTC on 5 September, the 24-hour average precipitation of Guangdong Province was 39.9 mm, with the largest rainfall occurring in the Pearl River Delta. Figure 9 shows the sample observations of Z H , Z D R , ρ H V , Ψ D P , and K D P from the Guangzhou radar 0.5- degree sweep at 1542 UTC, 3 September 2017. The rainfall accumulation is computed by the adapted algorithm from 1600 UTC on 3 September to 1600 UTC on 4 September, and it is based on 0.5-degree scans. Moreover, Figure 10 shows the time-series comparison of radar moments simulated by DSD data with the T-matrix method, and measured polarimetric parameters. It shows that precipitation in the position of the disdrometer mainly occurred from 2330 UTC on 3 September to 1600 UTC on 4 September, and the measured Z H , Z D R and K D P are close to the values of simulated radar moments during this typhoon event.

3.4. Microphysical Properties of Typhoon-induced Precipitation

In this section, the microphysical properties of typhoon-induced rainfall are investigated. The mass-weighted mean diameter D M and the intercept parameter N W , which can reflect the DSD characteristics [25], are calculated from the DSD data:
D M =   N ( D ) D 4 d D   N ( D ) D 3 d D ,
N W = 4 4 π ρ W ( W D M 4 ) ,
W = π 6 ρ W   N ( D ) D 3 d D ,
where D (mm) stands for equivalent diameter; N ( D ) (mm−1·m−3) is the number concentration of raindrop per unit volume per unit size interval; ρ W is the water density; W is the water content. Figure 11 shows the D M versus log 10 N W of the stratiform (reflectivity Z H at the height of 3 km is less than 40 dBZ) and convective system (reflectivity Z H at the height of 3 km is larger than 40 dBZ), which were retrieved from DSD data during typhoons Hato and Mawar. The intersection of the cross indicates average values of D M and log 10 N W , and the lengths of the horizontal line and the vertical line indicate variances of D M and log 10 N W , respectively. Figure 11 also shows the D M and log 10 N W for the convective system (including maritime and continental type) and the stratiform system from Bringi et al. 2003 [25]. The maritime (continental) type of convective system presents D M ~1.5–1.75 mm (2–2.75 mm) and log 10 N W ~4–4.5 (3–3.5), and the stratiform precipitation system has D M ~1.25–1.75 mm and log 10 N W ~3–4. The raindrop characteristics for the convective system during Typhoon Hato in this study are consistent with that in maritime convective precipitation from Bringi et al. 2003 [25], while the raindrop characteristics of stratiform precipitation during Typhoon Hato are smaller (Figure 11a). For Typhoon Mawar, convective and stratiform precipitation demonstrate smaller sizes and higher concentrations of raindrops than those in Typhoon Hato (Figure 11). Differences in microphysical processes may lead to different raindrop distributions in individual typhoons. An investigation on the raindrop characteristics of typhoons landing in China shows that convective precipitation in typhoons presents a lower raindrop diameter and a higher raindrop concentration than that of the maritime convection [26]. Our results for Typhoon Mawar are in accordance with this conclusion.
The radars observed the polarimetric parameters, including Z H , Z D R and K D P , from 0000 UTC to 2300 UTC, 23 August, 0000 UTC to 1200 UTC, 27 August, and 1600 UTC, 3 September to 1600 UTC, 4 September. These polarimetric observations can reflect the microphysical properties of hydrometeors, including precipitation particle shape and size distributions. Figure 12 shows the scatter density plots of Z D R versus Z H and log 10 [ K D P / ( 10 Z H / 10 ) ] versus Z D R for three typhoon events. The scatters are mainly concentrated in the range of -0.5 to 1 dB, where the reflectivity is less than 40 dBZ, that also indicates the stratiform precipitation. The closer Z D R is to 0 dB, the smaller the precipitation particle is. The results indicate that the precipitation is mainly from the stratiform system, and the raindrop during three typhoon events is small-to-midsize particles. Many scatters are concentrated in the range of 0 to 2 dB, where log 10 K D P / 10 Z H / 10 is less than −4. The raindrop characteristics of the convective system during the three typhoon events are midsize particles.

4. QPE Results and Verification

4.1. Sample Results

The comparison of rain gauge measurements and the estimates from different radar rainfall algorithms for Typhoons Hato, Pakhar, and Mawar is shown in Figure 13. Among five rainfall algorithms, the adapted algorithm without HCA is applied to demonstrate HCA performance, which selects different rainfall relations based on the polarimetric parameters and hydrometeor classification results. For Hato and Pakhar, the data of 198 rain gauges within 100 km of the Yangjiang radar are used for comparison with the estimates computed by different algorithms, whereas the data of 809 rain gauges within 100 km of the Guangzhou radar are used for Mawar. As shown in Figure 13a, the estimates computed by the adapted algorithm agree well with the gauge measurements from 0700 to 0900 UTC, 23 August. The gauge, 712306, is located at a range of 32.05 km from the Yangjiang radar. Figure 13b shows the comparison of the accumulated rainfall corresponding to Figure 13a. Accumulation begins at 0000 UTC 23 August and the total gauge accumulation is 58.3 mm. The total accumulation of estimates from the adapted algorithm is 54.85 mm, whereas the total accumulation of estimates from the Yangjiang blended rainfall algorithm is 36.87 mm. For clarity, this paper mainly describes one of the gauges that is the best result, and other gauges (Figure 13c–j) are used to show various QPE results. The gauges illustrated in Figure 13c–j are chosen randomly, some of which show good agreement between radar estimates and gauge measurements, while others show relatively poor agreement. The adapted algorithm has better performance for QPE during typhoon events.
To further demonstrate the performance of different algorithms, the scatter plots of rain gauge measurements versus rainfall rates from various elevation angles are shown in Figure 14. The number of 2.4-degree scatter plots is less than 0.5-degree and 1.5-degree, which is related to data segmentation. The closer the scatters are to the 1:1 line, the better the estimation effect is. Figure 14a shows that most of the scatters ( C C = 93.05%) for the QPE results of the adapted algorithm are distributed around the 1:1 line, whereas the QPE results of the Yangjiang blended rainfall algorithm, Yangjiang Z–R relation, and WSR-88D Z–R relation are usually below the 1:1 line. Although the results of the adapted algorithm without HCA are similar to those of the adapted algorithm, it is only slightly ( C C = 0.55% not as good as the adapted algorithm due to the lack of HCA. In short, the adapted algorithm has better performance than other rainfall algorithms for QPE during typhoon events.

4.2. Evaluation Metrics and Results

To quantify the performance of the adapted rainfall algorithm, a set of metrics is computed, including the normalized mean bias ( N M B ), normalized standard error ( N S E ), Pearson’s correlation coefficients ( C C ), root-mean-square error ( R M S E ), and the Heidke skill score ( H S S ) [27,28,29], which are defined as follows:
N M B = N = 1 M ( R N G N ) N = 1 M G N ,
N S E = N = 1 M | R N G N | N = 1 M G N ,
C C = N = 1 M ( R N R N ¯ ) ( G N G N ¯ ) N = 1 M ( R N R N ¯ ) 2 N = 1 M ( G N G N ¯ ) 2 ,
R M S E = N = 1 M ( R N G N ) 2 M ,
H S S = 2 ( N 1 × N 4 N 2 × N 3 ) ( N 1 + N 3 ) ( N 3 + N 4 ) + ( N 2 + N 2 ) ( N 2 + N 4 )
where R N and G N represent the radar estimates and the rain gauge measurements at time frame N , respectively. M is the total sample number. The definitions of elements N 1 , N 2 , N 3 , and N 4 are given in Table 2, and H is the hourly rainfall threshold. The value of N M B can be positive or negative. A positive value indicates overestimation, and a negative value stands for underestimation. N S E and R M S E are non-negative and N S E is not affected by rainfall. The closer N M B , N S E , and R M S E are to 0, the smaller the error is. C C indicates the correlation between the radar estimate and the “ground truth”. The closer the value is to 1, the better the estimation effect is. Positive values of H S S indicate that the estimation is better, and a perfect estimation obtains a H S S of 1.
The results of N M B , N S E , C C , R M S E , and H S S of the five algorithms for Typhoons Hato, Pakhar, and Mawar are shown in Table 3. The N S E s of the adapted algorithm for the typhoon events are 38.38%, 30.4%, and 36.61%, respectively. The R M S E of the different algorithms for Hato are higher than for other typhoon events. The results show that the adapted algorithm has better performance for QPE than other algorithms, in terms of N M B , N S E , C C , R M S E , and H S S , but underestimates rainfall during the typhoon events. In general, the performance of the WSR-88D Z–R relation is better than that of the blended Yangjiang algorithm and Yangjiang Z–R relation for the three typhoon events combined, in terms of N M B , N S E , C C , and H S S ( H = 8). The WSR-88D Z–R relation has lower N S E than the blended Yangjiang algorithm, but the blended Yangjiang algorithm has higher C C for three typhoon events than the WSR-88D Z–R relation. There is a large statistical difference between two algorithms in terms of N M B , N S E , C C , and R M S E , although the Yangjiang Z–R relation and WSR-88D Z–R relation are both Z–R rainfall relation. However, for the three typhoon events combined, the blended Yangjiang algorithm, WSR-88D Z–R relation, and Yangjiang Z–R relation significantly underestimate the rainfall rate.
The adapted algorithm has the same calculation processes as the “adapted algorithm without HCA”, except for the inclusion of HCA. For the three typhoon events combined, the R M S E s of the two algorithms are 2.74 mm and 2.78 mm, respectively, and the N S E s of the two algorithms are 36.66% and 37.83%, respectively. The H S S s ( H = 0) of the two algorithms are 68.52% and 34.05% for the three typhoon events combined, and the H S S s ( H = 8) are 78.89% and 77.64%. The evaluation results from the adapted algorithm with and without HCA are similar, except the H S S s ( H = 0). To investigate the reason, the polarimetric parameters and HCA results are analyzed. The adapted algorithm directly calculates rain rate using different rainfall relations, based on the hydrometeor classification results and polarimetric parameters, whereas other algorithms without HCA compute rainfall rate considering co-polar correlation coefficient ρ H V . The reflectivity is very small or 0 where ρ H V is less than 0.9. If ρ H V is not considered, most of rain rates would be a constant, i.e., intercept of Z–R relation. Moreover, if the ice particles are partially melted, rain gauge would measure rainfall amount, but the radar estimated rainfall would be zero, as ρ H V is less than 0.9 due to the ice contamination. Hence, using ρ H V as the threshold is the main cause of the difference. In addition, the rainfall rate where “rain-hail mixture” is present, is ignored for the algorithms without HCA. As a result, the adapted algorithm has better performance for QPE than the adapted algorithm without HCA.
Compared to the adapted algorithm without HCA, the Yangjiang blended rainfall algorithm has different rainfall relations. The R M S E of the adapted algorithm without HCA is 1.17 mm less than that of the Yangjiang blended rainfall algorithm for the three typhoon events combined. The QPE results indicate that different rainfall relations can make a big difference in the QPE performance for two algorithms. The specific rainfall relations of the adapted algorithm with and without HCA, and the Yangjiang blended rainfall algorithm are both derived from DSD data collected by a disdrometer in Yangjiang from April to July 2014, but the specific rainfall relations of the adapted algorithm with and without HCA are additionally adjusted for typhoon events. The specific rainfall relations of the adapted algorithm with and without HCA have better performance for QPE than the Yangjiang blended rainfall algorithm, due to their representation of the precipitation microphysical process. As the WSR-88D Z–R relation has been adjusted and demonstrated in the past, it better reflects the precipitation microphysical process than the Yangjiang Z–R relation. The WSR-88D Z–R relation shows better evaluation results than the Yangjiang Z–R relation for three typhoon events. The statistical difference between the WSR-88D Z–R relation and the Yangjiang Z–R relation further proves the conclusion about the effect of different rainfall relations. Due to various precipitation types in the typhoon system, different rainfall equations, i.e., R( Z H ), R( K D P ), R( Z H , Z D R ), and R( Z D R , K D P ), should be mixed and applied for typhoon events. The blended Yangjiang algorithm has additional rainfall relations compared to the Yangjiang Z–R relation, and the QPE performance of the Yangjiang algorithm is better than the Yangjiang Z–R relation.

5. Conclusions

Typhoons are often accompanied by strong convective weather such as high winds and heavy rainfall, resulting in catastrophic flash floods. Accurate QPE is one of the essential requirements for hazard mitigation and disaster management during typhoon events. In this paper, an adapted radar rainfall system combining data QC and segmentation, HCA, and quantitative rainfall estimation is detailed and applied for three typhoon events in southern China. Dual-polarization radar observations from the Yangjiang and Guangzhou radars are examined to demonstrate the performance of the rainfall estimation algorithm. These radars are among the first batch of the dual-polarization upgrade of the operational weather radar network in China. Therefore, the results of this study will facilitate the operational use of polarimetric measurements in QPE applications in China.
In addition, precipitation microphysical analyses are examined, along with the detailed QPE analysis. The size and concentration of raindrops in convective precipitation during Typhoon Hato in this study are consistent with the results in Bringi et al. [25] for “maritime-like” precipitation clusters and Wen et al. [26] for landing typhoons in China. The convective and stratiform precipitation of Typhoon Mawar shows a smaller size and a higher concentration of raindrops than that of Typhoon Hato. Assuming rainfall measurements from the ground rain gauges are the “truth”, the rainfall estimates computed by different algorithms are evaluated. The results show that the QPE performance of the adapted algorithm is better than other algorithms for all the three typhoon events, in terms of N M B , N S E , C C , R M S E , and H S S . The evaluation results from the adapted algorithm with and without HCA are similar, except the H S S s ( H = 0). Although HCA combining different polarimetric observables can slightly improve the QPE performance, HCA is inadequate to quantify the rainfall rate where “rain-hail mixture” is present. In addition, different rainfall relations offer various QPE performance for typhoon events due to their representation of the precipitation microphysical process. Blending of different rainfall equations, i.e., R( Z H ), R( K D P ), R( Z H , Z D R ), and R( Z D R , K D P ) shows better performance for QPE than a single rainfall equation. In conclusion, this work provides a comprehensive radar rainfall analysis of representing typhoon events. It is expected that the local meteorological bureaus can implement the improved radar rainfall algorithms to enhance rainfall mapping, which in turn would help to mitigate losses from hazardous weather phenomena. Future work will focus on microphysical retrievals using polarimetric radar measurements during typhoon events.

Author Contributions

Q.X. and W.Z. carried out the detailed analysis. H.C. supervised this study. Q.X. and H.C. prepared the first draft of this manuscript. W.-c.L. and L.H. reviewed and edited the manuscript. Q.X., H.C., W.-c.L., W.Z., Y.M. discussed the analysis results. X.L. provided the disdrometer data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China under Grants 41875049 and 41605022, and the National Key Research and Development Program of China under Grants 2018YFC1507504-6 and 2017YFC1501502. The participation of H.C. was supported by the NOAA Physical Sciences Laboratory in Boulder, Colorado.

Acknowledgments

The authors would like to express their gratitude to the anonymous reviewers for their comments that improved this manuscript. The radar and rain gauge data were provided by the China Meteorological Administration (CMA).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wolff, D.B.; Petersen, W.A.; Tokay, A.; Marks, D.A.; Pippitt, J.L. Assessing Dual-Polarization Radar Estimates of Extreme Rainfall during Hurricane Harvey. J. Atmos. Oceanic Technol. 2019, 36, 2501–2520. [Google Scholar] [CrossRef]
  2. Baeck, M.L.; Smith, J.A. Rainfall Estimation by the WSR-88D for Heavy Rainfall Events. Wea. Forecast. 2010, 13, 416–436. [Google Scholar] [CrossRef]
  3. Bringi, V.N.; Chandrasekar, V. Polarimetric Doppler Weather Radar: Principles and Applications; Cambridge University Press: Cambridge, UK, 2001; p. 648. [Google Scholar]
  4. Seliga, T.A.; Bringi, V.N. Potential use of radar differential reflectivity measurements at orthogonal polarizations for measuring precipitation. J. Appl. Meteor. 1976, 15, 69–76. [Google Scholar] [CrossRef] [Green Version]
  5. Wen, G.; Chen, H.; Zhang, G.; Sun, J. An inverse model for raindrop size distribution retrieval with polarimetric variables. Remote Sens. 2018, 10, 1179. [Google Scholar] [CrossRef] [Green Version]
  6. Ryzhkov, A.V.; Schuur, T.J.; Burgess, D.W.; Heinselman, P.L.; Giangrande, S.E.; Zrnić, D.S. The joint polarization experiment: Polarimetric rainfall measurements and hydrometeor classification. Bull. Amer. Meteor. Soc. 2005, 86, 809–824. [Google Scholar] [CrossRef] [Green Version]
  7. Ryzhkov, A.; Diederich, M.; Zhang, P.; Simmer, C. Potential Utilization of Specific Attenuation for Rainfall Estimation, Mitigation of Partial Beam Blockage, and Radar Networking. J. Atmos. Oceanic Technol. 2014, 31, 599–619. [Google Scholar] [CrossRef]
  8. Chandrasekar, V.; Bringi, V.N.; Balakrishnan, N.; Zrnić, D.S. Error structure of multiparameter radar and surface measurements of rainfall. J. Atmos. Ocean. Technol. 1988, 5, 783–795. [Google Scholar] [CrossRef] [Green Version]
  9. Gorgucci, E.; Scarchilli, G.; Chandrasekar, V. A robust estimator of rainfall rate using differential reflectivity. J. Atomos. Ocean. Technol. 1994, 11, 586–592. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, H.; Chandrasekar, V. The quantitative precipitation estimation system for Dallas–Fort Worth (DFW) urban remote sensing network. J. Hydrol. 2015, 531, 259–271. [Google Scholar] [CrossRef] [Green Version]
  11. Gou, Y.; Ma, Y.; Chen, H.; Yin, J. Utilization of a C-band polarimetric radar for severe rainfall event analysis in complex terrain over Eastern China. Remote Sens. 2019, 11, 22. [Google Scholar] [CrossRef] [Green Version]
  12. Cifelli, R.; Chandrasekar, V.; Lim, S.; Kennedy, P.C.; Wang, Y.; Rutledge, S.A. A new dual-polarization radar rainfall algorithm: Application in Colorado precipitation events. J. Atmos. Ocean. Technol. 2011, 28, 352–364. [Google Scholar] [CrossRef] [Green Version]
  13. Huang, H.; Zhao, K.; Chen, H.; Hu, D.; Fu, P.; Lin, Q.; Yang, Z. Improved Attenuation-based Radar Precipitation Estimation Considering the Azimuthal Variabilities of Microphysical Properties. J. Hydrometeor. 2020, (in press). [Google Scholar]
  14. Bechini, R.; Chandrasekar, V. A semisupervised robust hydrometeor classification method for dual-polarization radar applications. J. Atmos. Oceanic Technol. 2015, 32, 22–47. [Google Scholar] [CrossRef]
  15. Yu, Z.; Wang, Y.; Xu, H.; Davidson, N.; Chen, Y.; Chen, Y.; Yu, H. On the relationship between intensity and rainfall distribution in tropical cyclones making landfall over China. J. Appl. Meteor. Climatol. 2017, 56, 2883–2901. [Google Scholar] [CrossRef]
  16. Wang, Y.; Chandrasekar, V. Algorithm for Estimation of the Specific Differential Phase. J. Atmos. Ocean. Technol. 2009, 26, 2565–2578. [Google Scholar] [CrossRef]
  17. Chen, H.; Chandrasekar, V. An improved dual-polarization radar rainfall algorithm (DROPS2.0): Application in NASA IFloods field campaign. J. Hydrometeor. 2017, 18, 917–937. [Google Scholar] [CrossRef]
  18. Chen, H.; Lim, S.; Chandrasekar, V.; Jang, B.J. Urban hydrological applications of dual-polarization x-band radar: Case study in korea. J. Hydrol. Eng. 2016, 22, E5016001.1–E5016001.13. [Google Scholar] [CrossRef]
  19. Gou, Y.; Chen, H.; Zheng, J. An improved self-consistent approach to attenuation correction for C-band polarimetric radar measurements and its impact on quantitative precipitation estimation. Atmos. Res. 2019, 226, 32–48. [Google Scholar] [CrossRef]
  20. Waterman, P.C. Matrix formulation of electromagnetic scattering. Proc. IEEE. 1965, 53, 805–812. [Google Scholar] [CrossRef]
  21. Thurai, M.; Huang, G.J.; Bringi, V.N.; Randeu, W.L.; Schoenhuber, M. Drop shapes, model comparisons, and calculations of polarimetric radar parameters in rain. J. Atmos. Ocean. Technol. 2007, 24, 1019–1032. [Google Scholar] [CrossRef]
  22. Maddox, R.A.; Zhang, J.; Gourley, J.J.; Howard, K.W. Weather radar coverage over the contiguous United States. Wea. Forecast. 2002, 17, 927–934. [Google Scholar] [CrossRef]
  23. Matrosov, S.Y.; Ralph, F.M.; Neiman, P.J.; White, A.B. Quantitative assessment of operational weather radar rainfall estimates over California’s Northern Sonoma County using HMT-West data. J. Hydrometeorol. 2014, 15, 393–410. [Google Scholar] [CrossRef] [Green Version]
  24. Fulton, R.A.; Breidenbach, J.P.; Seo, D.J.; Miller, D.A.; O’Bannon, T. The WSR-88D Rainfall Algorithm. Wea. Forecast. 1998, 13, 377–395. [Google Scholar] [CrossRef]
  25. Bringi, V.N.; Chandrasekar, V.; Hubbert, J.; Gorgucci, E.; Randeu, W.L.; Schoenhuber, M. Raindrop size distribution in different climatic regimes from disdrometer and dual-polarized radar analysis. J. Atmos. Sci. 2003, 60, 354–365. [Google Scholar] [CrossRef]
  26. Wen, L.; Zhao, K.; Chen, G.; Wang, M.; Zhou, B.; Huang, H.; Hu, D.; Lee, W.C.; Hu, H. Drop size distribution characteristics of seven typhoons in China. J. Geophys. Res. Atmos. 2018, 123, 6529–6548. [Google Scholar] [CrossRef]
  27. Barnston, A.G. Correspondence among the correlation, RMSE, and Heidke forecast verification measures: Refinement of the Heidke score. Wea. Forecast. 1992, 7, 699–709. [Google Scholar] [CrossRef] [Green Version]
  28. Conner, M.D.; Petty, G.W. Validation and intercomparison of SSM/I rain-rate retrieval methods over the continental Unites States. J. Appl. Meteor. 1998, 37, 679–700. [Google Scholar] [CrossRef]
  29. Anagnostou, M.N.; Kalogiros, J.; Marzano, F.S.; Anagnostou, E.N.; Montopoli, M.; Piccioti, E. Performance evaluation of a new dual-polarization microphysical algorithm based on long-term x-band radar and disdrometer observations. J. Hydrometeorol. 2013, 14, 560–576. [Google Scholar] [CrossRef]
Figure 1. (a) The best tracks of Typhoons Hato (orange line), Pakhar (yellow line), and Mawar (magenta line); (b) Locations of rain gauges (red asterisks), S-band dual-polarization radars (white triangles) and disdrometers (white pentagram) in Guangzhou and Yangjiang overlaid on the terrain. The two blue circles denote the 100 km coverage range of the radar.
Figure 1. (a) The best tracks of Typhoons Hato (orange line), Pakhar (yellow line), and Mawar (magenta line); (b) Locations of rain gauges (red asterisks), S-band dual-polarization radars (white triangles) and disdrometers (white pentagram) in Guangzhou and Yangjiang overlaid on the terrain. The two blue circles denote the 100 km coverage range of the radar.
Remotesensing 12 02058 g001
Figure 2. Sample observations from the Yangjiang radar 1.5-degree sweep at 2112 UTC, 23 August 2017: (a) measured Z H ; (b) Z H after data quality control; (c) measured Z D R ; (d) Z D R after data quality control.
Figure 2. Sample observations from the Yangjiang radar 1.5-degree sweep at 2112 UTC, 23 August 2017: (a) measured Z H ; (b) Z H after data quality control; (c) measured Z D R ; (d) Z D R after data quality control.
Remotesensing 12 02058 g002
Figure 3. (a) Z D R versus Z H before data quality control; (b) Z D R versus Z H after data quality control. The black dots are the simulated radar moments based on raindrop size distribution data.
Figure 3. (a) Z D R versus Z H before data quality control; (b) Z D R versus Z H after data quality control. The black dots are the simulated radar moments based on raindrop size distribution data.
Remotesensing 12 02058 g003
Figure 4. Differences between the various radar rainfall algorithms applied in this paper: (a) adapted algorithm; (b) adapted algorithm without hydrometeor classification algorithm (HCA); (c) Yangjiang blended rainfall algorithm; (d) Yangjiang Z–R relation; (e) WSR-88D Z–R relation. Compared to the Yangjiang blended rainfall algorithm, the adapted algorithm has different specific rainfall relations and HCA. The adapted algorithm without HCA has the same modules as the adapted algorithm except HCA. The Yangjiang blended rainfall algorithm has more rainfall relations rather than the Yangjiang Z–R relation.
Figure 4. Differences between the various radar rainfall algorithms applied in this paper: (a) adapted algorithm; (b) adapted algorithm without hydrometeor classification algorithm (HCA); (c) Yangjiang blended rainfall algorithm; (d) Yangjiang Z–R relation; (e) WSR-88D Z–R relation. Compared to the Yangjiang blended rainfall algorithm, the adapted algorithm has different specific rainfall relations and HCA. The adapted algorithm without HCA has the same modules as the adapted algorithm except HCA. The Yangjiang blended rainfall algorithm has more rainfall relations rather than the Yangjiang Z–R relation.
Remotesensing 12 02058 g004
Figure 5. The flowcharts of the rainfall algorithms: (a) adapted rainfall algorithm selects different specific rainfall relations based on the hydrometeor classification results, as well as the polarimetric parameters (adapted from Figure 2 in Chen et al. [17]); (b) adapted rainfall algorithm without HCA and Yangjiang blended rainfall algorithm; (c) Yangjiang Z–R relation and WSR-88D Z–R relation.
Figure 5. The flowcharts of the rainfall algorithms: (a) adapted rainfall algorithm selects different specific rainfall relations based on the hydrometeor classification results, as well as the polarimetric parameters (adapted from Figure 2 in Chen et al. [17]); (b) adapted rainfall algorithm without HCA and Yangjiang blended rainfall algorithm; (c) Yangjiang Z–R relation and WSR-88D Z–R relation.
Remotesensing 12 02058 g005
Figure 6. Sample observations from the Yangjiang radar 1.5-degree sweep at 0430 UTC, 23 August 2017: (a) Z H ; (b) Z D R ; (c) ρ H V ; (d) Ψ D P ; and (e) K D P . (f) The rainfall accumulation computed by the adapted algorithm from 0000 UTC to 2300 UTC, 23 August 2017.
Figure 6. Sample observations from the Yangjiang radar 1.5-degree sweep at 0430 UTC, 23 August 2017: (a) Z H ; (b) Z D R ; (c) ρ H V ; (d) Ψ D P ; and (e) K D P . (f) The rainfall accumulation computed by the adapted algorithm from 0000 UTC to 2300 UTC, 23 August 2017.
Remotesensing 12 02058 g006
Figure 7. The time-series comparison of radar moments simulated by raindrop size distribution (DSD) data with T-matrix method and measured polarimetric parameters during Super Typhoon Hato: (a) Z H , (b) Z D R , (c) K D P .
Figure 7. The time-series comparison of radar moments simulated by raindrop size distribution (DSD) data with T-matrix method and measured polarimetric parameters during Super Typhoon Hato: (a) Z H , (b) Z D R , (c) K D P .
Remotesensing 12 02058 g007
Figure 8. Same as Figure 6, but at 0012 UTC, 27 August 2017. The rainfall accumulation is from 0000 UTC to 1200 UTC, 27 August 2017.
Figure 8. Same as Figure 6, but at 0012 UTC, 27 August 2017. The rainfall accumulation is from 0000 UTC to 1200 UTC, 27 August 2017.
Remotesensing 12 02058 g008
Figure 9. Sample observations from the Guangzhou radar 0.5-degree sweep at 1542 UTC, 3 September 2017: (a) Z H ; (b) Z D R ; (c) ρ H V ; (d) Ψ D P ; and (e) K D P . (f) The rainfall accumulation computed by the adapted algorithm from 1600 UTC on 3 September to 1600 UTC on 4 September 2017.
Figure 9. Sample observations from the Guangzhou radar 0.5-degree sweep at 1542 UTC, 3 September 2017: (a) Z H ; (b) Z D R ; (c) ρ H V ; (d) Ψ D P ; and (e) K D P . (f) The rainfall accumulation computed by the adapted algorithm from 1600 UTC on 3 September to 1600 UTC on 4 September 2017.
Remotesensing 12 02058 g009
Figure 10. The time-series comparison of radar moments simulated by DSD data with T-matrix method and polarimetric parameters during Severe Tropical Storm Mawar: (a) Z H , (b) Z D R , (c) K D P .
Figure 10. The time-series comparison of radar moments simulated by DSD data with T-matrix method and polarimetric parameters during Severe Tropical Storm Mawar: (a) Z H , (b) Z D R , (c) K D P .
Remotesensing 12 02058 g010
Figure 11. D M versus log 10 N W for (a) stratiform and (b) convective rain during Typhoons Hato (red dot) and Mawar (blue triangle). The green and orange crosses indicate means (cross centers) and standard deviations of D M and log 10 N W for Typhoons Hato and Mawar, respectively. The light blue dashed line in (a) represents stratiform rainfall in Bringi et al. [25] The black rectangles in (b) represent the maritime (the left one) and continental (the right one) types of convective rainfall.
Figure 11. D M versus log 10 N W for (a) stratiform and (b) convective rain during Typhoons Hato (red dot) and Mawar (blue triangle). The green and orange crosses indicate means (cross centers) and standard deviations of D M and log 10 N W for Typhoons Hato and Mawar, respectively. The light blue dashed line in (a) represents stratiform rainfall in Bringi et al. [25] The black rectangles in (b) represent the maritime (the left one) and continental (the right one) types of convective rainfall.
Remotesensing 12 02058 g011
Figure 12. (a) Z D R versus Z H ; (b) log 10 [ K D P / ( 10 Z H / 10 ) ] versus Z D R . The black dots are the simulated radar moments based on raindrop size distribution data.
Figure 12. (a) Z D R versus Z H ; (b) log 10 [ K D P / ( 10 Z H / 10 ) ] versus Z D R . The black dots are the simulated radar moments based on raindrop size distribution data.
Remotesensing 12 02058 g012
Figure 13. The rain gauge measurements (green line) and estimates computed by various radar algorithms at five selected gauge locations: (a), (c), (e), (g), (i) show the comparison of hourly rainfall rates; (b), (d), (f), (h), (j) show the comparison of accumulated rainfall. The rain gauge IDs are 712306, 716530, 713548, 716520, and 711166. The red dashed lines, yellow lines, blue lines, pink lines and black lines represent the rainfall estimates from the adapted algorithm, the adapted algorithm without HCA, Yangjiang blended rainfall algorithm, Yangjiang Z–R relation, and WSR-88D Z–R relation, respectively.
Figure 13. The rain gauge measurements (green line) and estimates computed by various radar algorithms at five selected gauge locations: (a), (c), (e), (g), (i) show the comparison of hourly rainfall rates; (b), (d), (f), (h), (j) show the comparison of accumulated rainfall. The rain gauge IDs are 712306, 716530, 713548, 716520, and 711166. The red dashed lines, yellow lines, blue lines, pink lines and black lines represent the rainfall estimates from the adapted algorithm, the adapted algorithm without HCA, Yangjiang blended rainfall algorithm, Yangjiang Z–R relation, and WSR-88D Z–R relation, respectively.
Remotesensing 12 02058 g013
Figure 14. The scatter plots of radar rainfall estimates versus rain gauge measurements for all the events combined: (a) the adapted algorithm, (b) the adapted algorithm without HCA, (c) Yangjiang blended rainfall algorithm, (d) Yangjiang Z–R relation, and (e) WSR-88D Z–R relation.
Figure 14. The scatter plots of radar rainfall estimates versus rain gauge measurements for all the events combined: (a) the adapted algorithm, (b) the adapted algorithm without HCA, (c) Yangjiang blended rainfall algorithm, (d) Yangjiang Z–R relation, and (e) WSR-88D Z–R relation.
Remotesensing 12 02058 g014
Table 1. System specifications of the radars at Guangzhou and Yangjiang in Guangdong Province.
Table 1. System specifications of the radars at Guangzhou and Yangjiang in Guangdong Province.
TypeRadar Name
Yangjiang RadarGuangzhou Radar
Antenna altitude (m)106179
Frequency (MHZ)27832885
PolarizationDual linearDual linear
Pulse width (ns)15701570
Peak power output (kW)650650
Cut number1111
Range gates18401840
Sampling resolution (m)250250
Beam width (degree)0.920.92
Table 2. Elements N 1 , N 2 , N 3 , and N 4 are assigned the observed event counts in each category.
Table 2. Elements N 1 , N 2 , N 3 , and N 4 are assigned the observed event counts in each category.
RadarRain Gauge
Observation   >   H mm Observation     H mm
Estimate >   H   mmN1N2
Estimate ≤   H   mmN3N4
Table 3. Evaluation results of the various algorithms for Typhoons Hato, Pakhar, and Mawar. N M B is the normalized mean bias, N S E is normalized standard error, C C is the Pearson’s correlation coefficient, R M S E is the root mean square error, and H S S is the Heidke skill score.
Table 3. Evaluation results of the various algorithms for Typhoons Hato, Pakhar, and Mawar. N M B is the normalized mean bias, N S E is normalized standard error, C C is the Pearson’s correlation coefficient, R M S E is the root mean square error, and H S S is the Heidke skill score.
EventMetricsAlgorithm
Adapted AlgorithmAdapted Algorithm without HCABlended Yangjiang AlgorithmYangjiang Z–R RelationWSR-88D Z–R Relation
HatoNMB (%)−15.46−14.95−47.97−58.21−39.00
NSE (%)38.3839.5452.5161.1648.68
CC (%)92.9992.5990.7186.8087.02
RMSE (mm)4.434.485.907.336.07
HSS (%, H = 0)63.1336.5836.7436.7036.42
HSS (%, H = 8)75.8574.5858.9850.3465.61
PakharNMB (%)−20.27−20.47−56.11−64.22−47.51
NSE (%)30.432.6356.9164.9649.13
CC (%)94.7693.3791.3993.8294.08
RMSE (mm)2.552.794.425.013.90
HSS (%, H = 0)83.0632.1932.3732.2932.22
HSS (%, H = 8)89.1487.6255.3034.2669.30
MawarNMB (%)−15.53−15.97−55.19−58.84−39.65
NSE (%)36.6137.7158.1261.5047.36
CC (%)92.9292.3691.0591.1391.72
RMSE (mm)2.372.403.473.752.88
HSS (%, H = 0)68.1333.3633.4333.3133.32
HSS (%, H = 8)78.4277.1851.4846.9563.39
All eventsNMB (%)−15.78−16.00−53.54−58.99−39.93
NSE (%)36.6637.8356.7361.6147.77
CC (%)93.0592.5091.0289.7490.24
RMSE (mm)2.742.783.954.493.56
HSS (%, H = 0)68.5234.0534.1334.0333.91
HSS (%, H = 8)78.8977.6453.7347.0064.55

Share and Cite

MDPI and ACS Style

Xia, Q.; Zhang, W.; Chen, H.; Lee, W.-C.; Han, L.; Ma, Y.; Liu, X. Quantification of Precipitation Using Polarimetric Radar Measurements during Several Typhoon Events in Southern China. Remote Sens. 2020, 12, 2058. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12122058

AMA Style

Xia Q, Zhang W, Chen H, Lee W-C, Han L, Ma Y, Liu X. Quantification of Precipitation Using Polarimetric Radar Measurements during Several Typhoon Events in Southern China. Remote Sensing. 2020; 12(12):2058. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12122058

Chicago/Turabian Style

Xia, Qiulei, Wenjuan Zhang, Haonan Chen, Wen-Chau Lee, Lei Han, Yu Ma, and Xiantong Liu. 2020. "Quantification of Precipitation Using Polarimetric Radar Measurements during Several Typhoon Events in Southern China" Remote Sensing 12, no. 12: 2058. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12122058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop