Microphysical and Polarimetric Radar Signatures of an Epic Flood Event in Southern China

An extremely heavy rainfall event hit Guangdong province, China, from 27 August to 1 September 2018. There were two different extreme rain regions, respectively, at the Pearl River estuary and eastern Guangdong, and a record-breaking daily precipitation of 1056.7 mm was observed at Gaotan station on 30 August. This paper utilizes a suite of observations from soundings, a gauge network, disdrometers, and polarimetric radars to gain insights to the two rainfall centers. The large-scale meteorological forcing, rainfall patterns, and microphysical processes, as well as radar-based precipitation signatures are investigated. It is concluded that a west-moving monsoon depression played a critical role in sustaining the moisture supply to the two extreme rain regions, and the combined orographic enhancement further contributed to the torrential rainfall over Gaotan station. The raindrop size distributions (DSD) observed at Zhuhai and Huidong stations, as well as the observed polarimetric radar signatures indicate that the rainfall at Doumen region was characterized by larger raindrops but a lower number concentration compared with that at Gaotan region. In addition, the dual-polarization radars are used to quantify precipitation intensity during this extreme event, providing timely information for flood warning and emergency management decision-making.


Introduction
Torrential rainfall events are one of the most severe disasters around the world [1,2]. The extreme rainfall and induced floods, landslides, debris flows gravely threaten life and property. The precipitation microphysics such as raindrop size distribution (DSD) serves as a fundamental bridge in deriving radar quantitative precipitation estimation (QPE) algorithms, which is critical for improving the accuracy of precipitation estimation and predictions [3,4]. Accurate precipitation estimates are also important input to the flash flood guidance systems for flood forecast, as well as subsequent warning operations and emergency management decision-making [5]. Therefore, a better understanding of precipitation microphysics and accurate quantitative precipitation estimation for extreme rain are important for flood warning and emergency management decision-making.

Data
The observational data used in this study include rainfall measurements from a dense gauge network, two S-band polarimetric radars, and two second-generation Particle Size and Velocity (Parsivel 2 ) disdrometers in Guangdong. The instrument locations and the two special gauge stations (Gaotan and Doumen) are shown in Figure 1. The National Centers for Environmental Prediction (NCEP) final operational model global analysis data (NCEP-FNL) available every 6 h with a resolution of 0.25 • × 0.25 • at 31 vertical levels (http://rda.ucar.edu/datasets/ds083.2//#!access) are used to resolve the synoptic condition [22], along with the sounding data collected at 00:00 and 12:00 UTC at Qingyuan (QY) and Shantou (ST) (No. 59280, and No. 59316 from the University of Wyoming: http://weather.uwyo.edu/upperair/sounding.html).

Data
The observational data used in this study include rainfall measurements from a dense gauge network, two S-band polarimetric radars, and two second-generation Particle Size and Velocity (Parsivel 2 ) disdrometers in Guangdong. The instrument locations and the two special gauge stations (Gaotan and Doumen) are shown in Figure 1. The National Centers for Environmental Prediction (NCEP) final operational model global analysis data (NCEP-FNL) available every 6 h with a resolution of 0.25° × 0.25° at 31 vertical levels (http://rda.ucar.edu/datasets/ds083.2//#!access) are used to resolve the synoptic condition [22], along with the sounding data collected at 00:00 and 12:00 UTC at Qingyuan (QY) and Shantou (ST) (No. 59280, and No. 59316 from the University of Wyoming: http://weather.uwyo.edu/upperair/sounding.html).

Raindrop Size Distribution
In order to understand the microphysics of the extreme rain, two disdrometers nearest to Doumen and Gaotan gauge stations are used to provide the DSD observation. These two disdrometers are located at Zhuhai and Huidong, collocated with rain gauges within 20 m. Both disdrometers are optical disdrometers with a 54 cm 2 horizontal sample area and are configured with 1-min sampling resolution to measure the DSD and fall velocity of raindrops [23,24]. The disdrometer performance has been assessed and improved since it was invented, and many previous studies have been conducted with this device [25][26][27][28]. In particular, the velocity and particle sizes are divided into 32

Raindrop Size Distribution
In order to understand the microphysics of the extreme rain, two disdrometers nearest to Doumen and Gaotan gauge stations are used to provide the DSD observation. These two disdrometers are located at Zhuhai and Huidong, collocated with rain gauges within 20 m. Both disdrometers are optical disdrometers with a 54 cm 2 horizontal sample area and are configured with 1-min sampling resolution to measure the DSD and fall velocity of raindrops [23,24]. The disdrometer performance has been assessed and improved since it was invented, and many previous studies have been conducted with this device [25][26][27][28]. In particular, the velocity and particle sizes are divided into 32 non-uniform Remote Sens. 2020, 12, 2772 4 of 22 bins, varying from 0.05 to 20.8 m s −1 for bin-center velocity and 0.062-24.5 mm for bin-center diameter (for detailed information in the user manual, https://www.manualslib.com/products/Ott-Parsivel2-5889584.html). The direct measurements from disdrometer are the number of raindrops at each velocity (i) and diameter (j) bin. Here, we take the bin-center value of each bin as the corresponding value. Several parameters used to describe the characteristics of DSD are calculated in the following.
The total number of raindrops can be calculated as follows: where n i,j is the number of drops at each bin. The number concentration of raindrops per unit volume for the jth diameter bin N(D j ) can be calculated as follows: where N(D j ) is in m −3 mm −1 ; A is the sampling area in m 2 ; ∆t is the sampling time interval in s; A and ∆t are, respectively, 0.0054 m 2 and 60 s in this study; ∆D j (mm) is the diameter interval from D j to D j+1 for the jth diameter bin; V i (m s −1 ) is the fall speed for the ith velocity class. Due to the measurement error, especially for larger size drops [23], the empirical terminal velocity-diameter (V -D) relation in Atlas et al. [29] is adopted in this study: The drops with velocity out the range of ±60% V(D j ) are removed from the analysis [30]. The total number concentration N t (m −3 ), the mass weighted diameter D m (mm), and normalized intercept parameter N w (m −3 mm −1 ) [14] are derived as: D m is closely related to the drop size; N t and N w are related to the number of raindrops. All these parameters are important in representing the DSD characteristics.
The integral rainfall parameters including rain rate R (mm h −1 ) and liquid water content W (g m −3 ) are calculated based on the following equations: where ρ w is the water density (1.0 g cm −3 ). Meanwhile, a series of polarimetric radar variables are simulated at S-band frequency based on the DSD measurements using the T-matrix method [31][32][33], including horizontal reflectivity Z h (mm 6 m −3 , or Z H in dBZ), differential reflectivity Z dr (dB), and specific differential phase K dp (degree km −1 ). The drop shape model used in the simulation is the one proposed by Thurai et al. [34] and temperature is 20°C. The canting angle is not taken into account (i.e., canting angle is 0) as the DSD measurements are near ground.
We also want to note that to minimize the measurement errors and improve data reliability, several quality control procedures were applied on the 1-min DSD data. First, because of the low signal-to-noise ratios, the first two diameter bins are always empty, so the data in first two bins are eliminated in the analysis [23]. Second, the 1-min sample data with total raindrop number smaller than 10 or the derived rain rate less than 0.1 mm h −1 are considered noise and removed [23]. Then, if the continuous data satisfying the above conditions last less than 5 min, they will be ignored to avoid the spurious and erratic measurements [25,35]. Additionally, threshold on simulated radar parameter (i.e., Z H < 55 dBZ) is used to further guarantee the creditability of the measured DSD data. The DSDs with the radar parameters out of the range are deleted to avoid mixed phase hydrometeors.

Radar Quantitative Precipitation Estimation
The two S-band radars are located at Guangzhou (hereafter referred as GZRD, 23.004 • N, 113.355 • E, 179 m) and Meizhou (MZRD, 24.256 • N, 115.975 • E, 423 m). Both radars are configured with a 6-min time resolution and 250-m range gate spacing, and have undergone rigorous quality control to ensure the data quality [36]. These two radars are used to monitor the evolution of the storm system and associated microphysical signatures. Moreover, these two radars are applied to estimate the rainfall to show the great potential of radar quantitative precipitation estimation.
For polarimetric radar, R(Z H ), R(K dp ), R(Z H , Z dr ), R(Z dr , K dp ) are the four relations commonly used to estimate rainfall [14]. The parameters (i.e., coefficients and exponents) in these relations are determined by the local precipitation microphysics, and are usually derived from the in situ DSD data. The representative of these parameters and how to combine different relations are the key issue in deriving radar QPE [14].
In this study, four rainfall algorithms are applied to quantify the precipitation intensity and amounts during this event. These algorithms belong to two categories: one is R(Z H ) relation, i.e., WRS-88D Z-R relationship [13] and localized Z-R relationship; another is a combination of the four rainfall relations (i.e., R(Z H ), R(K dp ), R(Z H , Z dr ), R(Z dr , K dp )), i.e., the "adapted algorithm" described by Xia et al. [36] and the localized blended relation. Both Z-R relations are commonly used for single-polarized radar. The "adapted algorithm" is derived from DSD data and adjusted based on gauge observation and it has been demonstrated during several typhoon events in Southern China, which showed great performance [36].
Based on the nonlinear least-square method with the DSD data from two stations, the localized relations are fitted as follows: R(K dp ) = 40.4615K 0.7703 dp (10) R(Z dr , K dp ) = 84.4318K 0.9377 DP 10 −0.1588Z dr (12) where Equation (9) is referred to as the localized Z-R relation, and the localized blended relation is a combination of Equations (10)-(12), using the same logic described by Xia et al. [36].

Synoptic Environment during This Epic Rainfall Event
From 27 to 31 August 2018, a monsoon depression evolved, moving west from the East China Sea. The resultant precipitation system was moving from west to east along the coast, causing heavy precipitation nearly over the whole Guangdong province, especially in the south and southeast parts of Guangdong on 29 and 30 August. The synoptic evolution of this precipitation system was analyzed using the NCEP-FNL data at 12:00 UTC 29 August and 12:00 UTC 30 August ( Figure 2). Additionally, Figure 3 shows the stratification curve and convective available potential energy (CAPE) based on the sounding data at QY and ST stations. Figure 3a illustrates the sounding curves at 12:00 UTC 29 August at QY station. The large convective available potential energy (CAPE, ≈860 J kg −1 ) and low level of free convection (LFC, ≈980 hPa) suggested that the atmospheric environment was favorable for the development of convection [37]. At 20:00 UTC 30 August, the CAPE became larger and LFC became lower, indicating that the convection had further developed. The soundings at ST station also showed a conducive environment for convection development though the CAPE decreased slightly from the 29th to 30th.  black vector arrows, vector scale at upper right), and the temperature advection (shaded, ×10 −4 °C s −1 ) at 850 hPa (c) at 12:00 UTC 29 August, and (d) at 12:00 UTC 30 August. The red (blue) shadings in panels (c) and (d) indicate warm (cold) temperature advection. The 1000-300 hPa vertically integrated water vapor transport (IVT, shaded, kg m −1 s −1 , also in black vector arrows, vector scale at upper right) (e) at 12:00 UTC 29 August, and (f) at 12:00 UTC 30 August. The dashed red rectangle indicates the area of the southeast of Guangdong province.  Figure 4 shows the spatial distribution of the cumulative rainfall in Guangdong during the period from 00:00 UTC 27 August to 00:00 UTC 2 September based on the gauge hourly precipitation It can be seen that the entire region of southern China was under the influence of a monsoon depression across the mid-and lower-troposphere on 29 August. The depression moved west, bringing strong moistures transported to Guangdong province driven by the persistent south and southwest wind (Figure 2a,b). In the coastal areas of Guangdong, warm advection from the ground and lower troposphere continued to transport warm air to inland areas (Figure 2c,d). Due to the combined evolution of depression and topography lifting effect, the precipitation system was further strengthened on 30 August. The wind and temperature changes between Guangdong and the adjacent South China Sea were more violent (Figure 2c,d). Meanwhile, the vertically integrated vapor transport (IVT) over the coastal zone was much stronger on 30 August than that on 29 August (Figure 2e,f). Figure 3a illustrates the sounding curves at 12:00 UTC 29 August at QY station. The large convective available potential energy (CAPE, ≈860 J kg −1 ) and low level of free convection (LFC, ≈980 hPa) suggested that the atmospheric environment was favorable for the development of convection [37]. At 20:00 UTC 30 August, the CAPE became larger and LFC became lower, indicating that the convection had further developed. The soundings at ST station also showed a conducive environment for convection development though the CAPE decreased slightly from the 29th to 30th. Figure 4 shows the spatial distribution of the cumulative rainfall in Guangdong during the period from 00:00 UTC 27 August to 00:00 UTC 2 September based on the gauge hourly precipitation data using Inverse Distance Weighted (IDW) interpolation method. This result was only used for visualization so as to get a general sense of the geographic distribution of rainfall. Obviously, this rainstorm affected most regions of Guangdong, with the extreme rainfall occurring at the southeast part and gradually decreased to the northwest of Guangdong. The gauge maximum cumulative rainfall was 1394.6 mm recorded at Gaotan station (23.1883 • N, 115.3044 • E), which set a new record of cumulative rainfall during a single event in Guangdong. Meanwhile, the maximum daily rainfall of 1056.7mm was also recorded at Gaotan station from 21:00 UTC 29 August to 21:00 UTC 30 August, which broke the historical daily rainfall record of 924.3 mm on 17 August 2013 [7].    Figure 5 shows the distribution of the daily rainfall that exceeded 100 mm from 28 to 31 August. It is clear that the rainfall is mainly concentrated on 29 and 30 August. There is a clear rain band along the coast on 29 August, with a gradient of accumulated rainfall decreasing from the coast towards the inland region, which can be partially attributed to the land-sea roughness contrasts [38]. On 30 August, the rain band moved a little towards northeast and inland region, with extremely high accumulated value (1041.1 mm) concentrated at Gaotan station, which could be due to the combined effect of terrain and evolved synoptic condition. To further reveal the two rainfall processes, Figure 6 shows the time series of rainfall collected at two gauges where extreme daily rainfall concentrated on these two days. One station is located at Doumen (22.1967 • N, 1153.1131 • E) in Zhuhai, representing the rainfall on 29 August, when the daily rainfall reached 458.6 mm. Another station is Gaotan station, representing the most extreme rainfall pattern both for 30 August and the whole event.

Precipitation Pattern Observed by a Gauge Network
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 23 Figure 5. Daily rainfall accumulation (in mm, a day is defined from 0000 UTC to 0000 UTC of the next day) for (a) 28 August, (b) 29 August, (c) 30 August, (d) 31 August, (e) a zoomed area for 30 August, (f) a zoomed area for 29 August. Each dot represents a gauge station with daily rainfall accumulation exceeding 100 mm.

DSDs Time Series at Two Observation Stations
As there are no disdrometers deployed at the Doumen and Gaotan stations, two nearby Parsivel 2 disdrometers were selected to analyze the raindrop size distribution from 00:00 UTC 29 August to 12:00 UTC 31 August. One disdrometer is located at Zhuhai (22.2750 °N, 113.5669 °E) near the Doumen station (45 km away); the other is located at Huidong (23.0261 °N, 114.6681 °E) near the Gaotan station (67 km away) (see Figure 1). The comparisons of hourly rainfall between the disdrometers and the gauges at Zhuhai and Huidong stations (within 20 m for each pair) are shown in Figure 7. Although the disdrometer might slightly underestimate the rainfall compared with gauge data because of the absence of some records and the overlap of drops along the laser beam, overall, they captured the rainfall pattern very well. As such, we ignored the impact on DSD observations in this study. Moreover, the rainfall pattern at Zhuhai station is similar to that of Doumen station, and the rainfall pattern of Huidong station is similar to Gaotan station. The similar patterns indicate that these two disdrometers could, respectively, represent the DSD characteristic of Doumen and Gaotan stations. Figure 8 shows the time series of the drop size distribution ( : m −3 mm −1 ) in the logarithmic scale from 00:00 UTC 29 August to 12:00 UTC 31 August, derived by the 1-min disdrometer observations at the two stations. In line with the synoptic condition shown in Figure 2, it can be seen that the Zhuhai station experienced stronger rainfall during 29 August, and it reached the maximum hourly rainfall of about 40 mm (Figure 7). The maximum raindrops greater than 7 mm were recorded from 12:00 to 14:00 UTC 29 August. After that, rainfall at Zhuhai station began to weaken with the decreasing number of raindrop concentration and size of raindrops, while that at Huidong station gradually increased. After 23:00 UTC 29 August, this long-lasting storm disappeared at Zhuhai station following with two short storms characterized with lower concentration but bigger raindrops around 04:00 UTC 30 August and around 11:00 UTC 31 August. The rain at Huidong station showed an increasing trend from 13:00 UTC 29 August, and the number of raindrop concentration and size of raindrops also show an increasing trend, and the hourly rainfall reached the highest from 18:00 to 19:00 UTC 30 August. Compared to the extreme rainfall period at Zhuhai station, the number of large The differences between the two stations are clear: Doumen station experienced heavy rainfall mainly on 29 August, while Gaotan station experienced a heavier and longer rainfall from late 29 August to early 31 August. Doumen station experienced three major rainfall stages on 29 August: first from 00:00 to 06:00 UTC, second from 06:00 to 14:00 UTC, and third from 14:00 to 23:00 UTC. The third stage was longer and stronger than the first two stages, with a maximum hourly rainfall of 70.2 mm from 1800 to 1900 UTC 29 August. After 0800 UTC 30 August, rainfall at Doumen station almost disappeared. At Gaotan station, rainfall exhibited two major stages: first from 00:00 to 20:00 UTC on 29 August, then from 20:00 UTC 29 August to 08:00 UTC 31 August. The first stage, with only one hourly rainfall exceeding 20 mm is much weaker than the second. During the second stage, the hourly rainfall showed an increasing trend, and reached a maximum of 98.6 mm from 1900 to 2000 UTC 30 August. Moreover, the cumulative rainfall of the second stage (1206.6 mm) accounted for 87% rainfall of the whole event at Gaotan station.

DSDs Time Series at Two Observation Stations
As there are no disdrometers deployed at the Doumen and Gaotan stations, two nearby Parsivel 2 disdrometers were selected to analyze the raindrop size distribution from 00:00 UTC 29 August to 12:00 UTC 31 August. One disdrometer is located at Zhuhai (22.2750 • N, 113.5669 • E) near the Doumen station (45 km away); the other is located at Huidong (23.0261 • N, 114.6681 • E) near the Gaotan station (67 km away) (see Figure 1). The comparisons of hourly rainfall between the disdrometers and the gauges at Zhuhai and Huidong stations (within 20 m for each pair) are shown in Figure 7. Although the disdrometer might slightly underestimate the rainfall compared with gauge data because of the absence of some records and the overlap of drops along the laser beam, overall, they captured the rainfall pattern very well. As such, we ignored the impact on DSD observations in this study. Moreover, the rainfall pattern at Zhuhai station is similar to that of Doumen station, and the rainfall pattern of Huidong station is similar to Gaotan station. The similar patterns indicate that these two disdrometers could, respectively, represent the DSD characteristic of Doumen and Gaotan stations.
Remote Sens. 2020, 12, x FOR PEER REVIEW 12 of 23 drops (D > 5 mm) at Huidong station is much less during the strong rainfall period which we hypothesize is due to the collision caused by sea-inland orographic effect [38][39][40].     Figure 2, it can be seen that the Zhuhai station experienced stronger rainfall during 29 August, and it reached the maximum hourly rainfall of about 40 mm (Figure 7). The maximum raindrops greater than 7 mm were recorded from 12:00 to 14:00 UTC 29 August. After that, rainfall at Zhuhai station began to weaken with the decreasing number of raindrop concentration and size of raindrops, while that at Huidong station gradually increased. After 23:00 UTC 29 August, this long-lasting storm disappeared at Zhuhai station following with two short storms characterized with lower concentration but bigger raindrops around 04:00 UTC 30 August and around 11:00 UTC 31 August. The rain at Huidong station showed an increasing trend from 13:00 UTC 29 August, and the number of raindrop concentration and size of raindrops also show an increasing trend, and the hourly rainfall reached the highest from 18:00 to 19:00 UTC 30 August. Compared to the extreme rainfall period at Zhuhai station, the number of large drops (D > 5 mm) at Huidong station is much less during the strong rainfall period which we hypothesize is due to the collision caused by sea-inland orographic effect [38][39][40]. Figure 9a shows the scatterplot of log 10 N w vs. D m of the two stations. The maritime convective clusters, continental convective clusters, and the stratiform rain line described in Bringi et al. [41] as well as the convection-stratiform separation line suggested by Thompson et al. [42] are also shown in Figure 9a. The scatter points of two stations show similar boundaries, of which there are both convective rain and stratiform rain types. Meanwhile, there are more points classified as stratiform rain by the classification suggested by Thompson et al. [42], and the averaged log 10 N w -D m pairs of two stations are close to the stratiform line described by Bringi et al. [41]. Moreover, for convective rain type, there are more points located at the continental convective cluster rather than maritime cluster for both stations. Compared with the DSD characteristics studied at Yangjiang from July to August by Tang et al. [43], the log 10 N w of both two stations are smaller no matter what rain types, suggesting a lower drop concentration of this rain system. Figure 7. Hourly rainfall measured by gauges and collocated disdrometers: (a) Zhuhai station, (b) Huidong station from 00:00 UTC 29 August to 12:00 UTC 31 August.  of the two stations. The maritime convective clusters, continental convective clusters, and the stratiform rain line described in Bringi et al. [41] as well as the convection-stratiform separation line suggested by Thompson et al. [42] are also shown in Figure 9a. The scatter points of two stations show similar boundaries, of which there are both convective rain and stratiform rain types. Meanwhile, there are more points classified as stratiform rain by the classification suggested by Thompson et al. [42], and the averaged log pairs of two stations are close to the stratiform line described by Bringi et al. [41]. Moreover, for convective rain type, there are more points located at the continental convective cluster rather than maritime cluster for both stations. Compared with the DSD characteristics studied at Yangjiang from July to August by Tang et al. [43], the log of both two stations are smaller no matter what rain types, suggesting a lower drop concentration of this rain system.

The Distribution of D m and N w
The probability density functions (PDFs) of and log with the mean values are presented in Figure 9b,c. The PDF of at Zhuhai station peaks around 1.25 mm while that at Huidong peaks around at 1.75 mm. Additionally, the PDF of log at Zhuhai station peaks around 3.25 while that at Huidong peaks around at 3.75. This leads to a higher mean value of (1.75 mm) and lower mean value of log (3.18) at Zhuhai station than those at Huidong station. This result indicates that the rain at Zhuhai station has larger drops, but a lower drop number concentration compared to Huidong station. and are represented by circles (gray for Zhuhai and red for Huidong). The two rectangles correspond to the maritime (cyan) and continental convective (orange) clusters reported by Bringi et al. [41]. The black and purple dash-dot lines represent the characteristics of stratiform rain described in Bringi et al. [41] and the convection/stratiform separation line from Thompson et al. [42], respectively. The square symbols represent the average values of various types of rain, from Tang et al. [43]. CV stands for convective rain, ST is stratiform rain. ZH is short for Zhuhai, and HD is short for Huidong. It is the same in the following figures and tables.

The DSD Spectra
The characteristics of DSD of this rain system are investigated to understand the precipitation microphysical processes. To eliminate to the effect caused by rain rate and duration, the DSD measurements are divided into six classes according to the associated rain rate ( ): C1, 0.1 1; (c) PDFs of log 10 N w at Huidong and Zhuhai stations. Both data are from 00:00 UTC 29 August to 12:00 UTC 31 August. The black dot and pink circle represent Zhuhai data and Huidong data, respectively. The mean values of log 10 N w and D m are represented by circles (gray for Zhuhai and red for Huidong). The two rectangles correspond to the maritime (cyan) and continental convective (orange) clusters reported by Bringi et al. [41]. The black and purple dash-dot lines represent the characteristics of stratiform rain described in Bringi et al. [41] and the convection/stratiform separation line from Thompson et al. [42], respectively. The square symbols represent the average values of various types of rain, from Tang et al. [43]. CV stands for convective rain, ST is stratiform rain. ZH is short for Zhuhai, and HD is short for Huidong. It is the same in the following figures and tables.
The probability density functions (PDFs) of D m and log 10 N w with the mean values are presented in Figure 9b,c. The PDF of D m at Zhuhai station peaks around 1.25 mm while that at Huidong peaks around at 1.75 mm. Additionally, the PDF of log 10 N w at Zhuhai station peaks around 3.25 while that at Huidong peaks around at 3.75. This leads to a higher mean value of D m (1.75 mm) and lower mean value of log 10 N w (3.18) at Zhuhai station than those at Huidong station. This result indicates that the rain at Zhuhai station has larger drops, but a lower drop number concentration compared to Huidong station.

The DSD Spectra
The characteristics of DSD of this rain system are investigated to understand the precipitation microphysical processes. To eliminate to the effect caused by rain rate and duration, the DSD measurements are divided into six classes according to the associated rain rate (R): C1, 0.1 ≤ R < 1; C2, 1 ≤ R < 5; C3, 5 ≤ R < 10; C4, 10 ≤ R < 25; C5, 25 ≤ R < 50 mm h −1 . Such classification is based on the rain rate distribution of this system, as well as several previous studies [25,35]. The DSD sample relative frequency, rain rate statistics and mean DSD parameters for each class are summarized in Tables 1 and 2. For lower rain rate classes (C1-C3), the rain rate of Huidong station is lower than that of Zhuhai station. For C5, the rain rates of two stations are close to each other and for the other two high rain rate classes (C4 and C6), the rain rate of Huidong station is higher than that of Zhuhai station, which caused the higher hourly rainfall at from 18:00 to 19:00 UTC 30 August. However, due to the lower relative frequency of high rain rate classes, the overall average rain rate of Huidong station is smaller than that of Zhuhai station. It can be shown from Figure 10 that the total averaged DSD of Zhuhai station has wider spread as well as lower concentration at small and medium-sized drops (D < 2 mm) but higher concentration at big and large drops (D > 2 mm) compared with that of Huidong station, which corresponds to a larger D m , more W, but a lower N t for Zhuhai station compared with those for Huidong station (Table 2). For the same rain class, the DSD of Huidong shows the same characteristics with higher concentration at smaller drops as well as lower concentration at larger drops. Additionally, the diameters where the raindrop concentrations of two stations become larger as the rain classes increased, from 1.2 to 3.8 mm, due to the higher the concentration of larger drops at higher rain rates for both stations.

Polarimetric Radar Signatures and Rainfall Analysis
As the rain patterns of Zhuhai station and Huidong station are similar to that of Doumen station and Gaotan station, the DSD characteristics of Zhuhai and Huidong areas can partly represent the microphysics of Doumen and Gaotan station, which suggests a high probability of a higher drop concentration at Gaotan station and a larger drops at Doumen station. However, the total rainfall difference between Huidong and Gaotan station is significant (≈400 mm vs. ≈1000 mm), which means the DSD of Gaotan may be different and need more observation to supplement the current study. In the following study, the dual-polarized radar observations are used to further understand the rainfall.

The Polarimetric Radar Signatures
To further understand the differences between the two rain patterns on 29 and 30 August, the polarimetric radar fields over the two extreme stations are compared. Figure 11 Table 3 summarizes the relative contributions of each size class to N t , W, and R, which are obtained through dividing the parameters calculated for the raindrops of each size class by the corresponding parameters calculated for all the size classes. For both stations, the small raindrops (D < 1 mm) have the highest concentration (51.99% and 57.92%) of the total data set and the second highest relative frequency (Td 38.08% and 44.26%), but relatively low contributions to total water content and rainfall. The contribution of small and medium-sized raindrops (D: 1-2 mm) to the total number is the greatest, accounting for 47.48% (45.99%) at Zhuhai (Huidong) station. However, for Zhuhai station, the middle-sized drops (D: 2-3 mm) accounts for the most rainfall (34.39%), while the small and medium-sized raindrops (D: 1-2 mm) account for the most for Huidong station (40.50%). Moreover, for Zhuhai station, the medium-sized and larger drops (D > 2 mm) account more for rainfall and water content compared with Huidong station, as the higher relative frequency. This result is consistent with the conclusion shown in Figures 8 and 9 that larger drops are less likely to occur at Huidong, where complex topography could exacerbate the collision.

Polarimetric Radar Signatures and Rainfall Analysis
As the rain patterns of Zhuhai station and Huidong station are similar to that of Doumen station and Gaotan station, the DSD characteristics of Zhuhai and Huidong areas can partly represent the microphysics of Doumen and Gaotan station, which suggests a high probability of a higher drop concentration at Gaotan station and a larger drops at Doumen station. However, the total rainfall difference between Huidong and Gaotan station is significant (≈400 mm vs. ≈1000 mm), which means the DSD of Gaotan may be different and need more observation to supplement the current study. In the following study, the dual-polarized radar observations are used to further understand the rainfall.

The Polarimetric Radar Signatures
To further understand the differences between the two rain patterns on 29 and 30 August, the polarimetric radar fields over the two extreme stations are compared. Figure 11 shows the time series of polarimetric variables over Doumen and Gaotan stations observed by GZRD and MZRD 0.5-degree sweeps, respectively. The size of radar pixel over Doumen (Gaotan) station is around 0.25 × 1.49 km (0.25 × 2.23 km), and the horizontal/vertical distances are 93 km/0.81 km (137 km/1.21 km) between the radar pixel and the station. Generally speaking, Figure 11 shows that the timing changes of radar variables are consistent with the rainfall changes at both stations ( Figure 6). During the extreme rainfall periods of Doumen station and Gaotan station, the mean values of Z H over both stations are similar while the mean value of Z dr over Doumen station is higher than that over Gaotan station, indicating that bigger drops are more likely to occur at Doumen station than Gaotan station, which is similar to the DSD analysis of Zhuhai and Huidong.
Remote Sens. 2020, 12, x FOR PEER REVIEW 16 of 23 between the radar pixel and the station. Generally speaking, Figure 11 shows that the timing changes of radar variables are consistent with the rainfall changes at both stations ( Figure 6). During the extreme rainfall periods of Doumen station and Gaotan station, the mean values of over both stations are similar while the mean value of over Doumen station is higher than that over Gaotan station, indicating that bigger drops are more likely to occur at Doumen station than Gaotan station, which is similar to the DSD analysis of Zhuhai and Huidong. To eliminate the random error of one radar pixel, two extreme regions were defined as circular areas with a radius of 5 km centered on Doumen station and Gaotan station, respectively, hereafter referred to Doumen rain region and Gaotan rain region. The polarimetric variables from 0.5-degree sweep over these two regions from 00:00 UTC 29 August to 12:00 UTC 31 August were analyzed. Figure 12 shows the joint frequency distributions of vs. , vs. , as well as log / vs. observed by radars over two regions as well as the simulated radar moments based on DSD data at Zhuhai station and Huidong station. The frequency at each pair of data bin (e.g., the bin width is 1 dBZ for and 0.1 dB for , 0.1 degree km −1 for , and 0.1 for log10( / )) refers to the To eliminate the random error of one radar pixel, two extreme regions were defined as circular areas with a radius of 5 km centered on Doumen station and Gaotan station, respectively, hereafter referred to Doumen rain region and Gaotan rain region. The polarimetric variables from 0.5-degree sweep over these two regions from 00:00 UTC 29 August to 12:00 UTC 31 August were analyzed. Figure 12 shows the joint frequency distributions of Z dr vs. Z H , K dp vs. Z H , as well as log 10 (K dp /Z h ) vs. Z dr observed by radars over two regions as well as the simulated radar moments based on DSD data at Zhuhai station and Huidong station. The frequency at each pair of data bin (e.g., the bin width is 1 dBZ for Z H and 0.1 dB for Z dr , 0.1 degree km −1 for K dp , and 0.1 for log 10 (K dp /Z h )) refers to the total number of radar pixels with values of Z dr vs. Z H (K dp vs. Z H or log 10 (K dp /Z h ) vs. Z dr ) falling into the bins, and the frequency distribution is further normalized by dividing by the maximum frequency among all the pairs of data bins.
Remote Sens. 2020, 12, x FOR PEER REVIEW 17 of 23 0.1 degree km −1 , increases faster as the increases at Doumen rain region which indicates that there are more relative larger drops at Doumen rain region. Then, considering the higher rain rate and smaller raindrops at Gaotan station, the number concentration must be larger. The extreme lower value of log / could confirm this, which means the concentration of small drops at Gaotan is higher than that at Doumen. It can be concluded from the polarimetric radar signatures that compared with Doumen rain region, the raindrops at Gaotan rain region have higher concentration but smaller drops, which is consistent with the DSD characteristic at Huidong station.  Frequencies for the observed Z H , Z dr at Doumen rain region peaked at approximately 28 dBZ, and 0.3 dB, respectively. While at Gaotan rain region they peaked at high values, around 35 dBZ for Z H and 0.8 dB for Z dr . The Z H maximum values with the 10% normalized frequency of two regions are similar around 52 dBZ, which are consistent with the extreme rainfall periods of two stations. Meanwhile, the Z dr maximum values with the 10% normalized frequency at Gaotan rain region is smaller (2.1 dB vs. 2.5 dB). Meanwhile, for the same Z H , the range and mean Z dr at Doumen rain region is larger than that at Gaotan station, which indicates that larger drops more likely occurred at Doumen compared with Gaotan. For both regions, the slope between Z H and K dp shows a sharp transition over sampling grids with Z H exceeding 40 dBZ, resulting in a rapid increase of K dp along with the increase of Z H . Although the peak frequencies for K dp at two rain regions are similar around 0.1 degree km −1 , K dp increases faster as the Z H increases at Doumen rain region which indicates that there are more relative larger drops at Doumen rain region. Then, considering the higher rain rate and smaller raindrops at Gaotan station, the number concentration must be larger. The extreme lower value of log 10 (K dp /Z h ) could confirm this, which means the concentration of small drops at Gaotan is higher than that at Doumen. It can be concluded from the polarimetric radar signatures that compared with Doumen rain region, the raindrops at Gaotan rain region have higher concentration but smaller drops, which is consistent with the DSD characteristic at Huidong station.

Radar-Based Quantitative Precipitation Estimation (QPE)
In this part, a set of various radar rainfall algorithms including the "adapted algorithm" [36], the localized blended algorithm, the localized relation and WSR-88D Z-R relation [13], were applied to quantify the precipitation intensity and amounts during this event. After estimating the instantaneous rainfall rates using various radar rainfall algorithms, the rainfall accumulations were computed at hourly scale. To quantify the performances of different algorithms, a set of metrics was computed, including the bias (BIAS), normalized mean bias (NMB), normalized mean absolute error (NMAE), Pearson's correlation coefficients (CC), which are defined as follows: where R N and G N represent the radar estimates of different algorithms and the rain gauge measurements at time frame N, respectively. M is the total sample number. The comparison of estimates from different rainfall algorithms and rain gauge measurements at Gaotan and Doumen station is shown in Figure 13 and the evaluation results are shown in Table 4. As shown in Figure 13, all the algorithms have similar patterns to the rain gauge measurements, while QPE results of the "adapted algorithm" is better than other algorithms with lowest BIAS, NMB, NMAE and the highest CC (Table 4). This is in line with the findings during typhoon case studies [36]. Nevertheless, we should note that all the algorithms are underestimating the rain rates and accumulations during this flood event and the results still need to be optimized. To further analyze the QPE results of various radar rainfall algorithms, the scatter plots of radar rainfall estimates versus gauge measurements at all rain gauge stations less than 100 km from the radars are shown in Figure 14 and the evaluation results are shown in Table 5. Most of the hourly rainfall rates are from 0 to 40 mm hr -1 . Again, the adapted algorithm has the best performance, while all the rainfall algorithms show underestimation compared to the gauge measurements, especially during heavy rain periods.  To further analyze the QPE results of various radar rainfall algorithms, the scatter plots of radar rainfall estimates versus gauge measurements at all rain gauge stations less than 100 km from the radars are shown in Figure 14 and the evaluation results are shown in Table 5. Most of the hourly rainfall rates are from 0 to 40 mm hr −1 . Again, the adapted algorithm has the best performance, while all the rainfall algorithms show underestimation compared to the gauge measurements, especially during heavy rain periods. Remote Sens. 2020, 12, x FOR PEER REVIEW 20 of 23 Figure 14. The scatter plots of radar estimated rainfall versus gauge measurements at gauge locations less than 100 km from the radar: (a) the "adapted algorithm" [36]; (b) localized blended rainfall algorithm; (c) localized Z-R relation; (d) WSR-88D Z-R relation [13]. The color density represents the observation sample numbers in logarithmic unit.

Discussion
Although the analysis of DSD and polarimetric radar signatures shows similar results, the relatively long distances between disdrometers and the extreme rain centers might induce some uncertainty in the representation of extreme rainfall DSD characteristics. As mentioned, the polarimetric radar signatures at the two rainfall centers are quite different, even for the same precipitation system, which is likely due to the complex falling processes as a result of orographic enhancement and dynamic cloud microphysics involved in this extreme event [44]. More Figure 14. The scatter plots of radar estimated rainfall versus gauge measurements at gauge locations less than 100 km from the radar: (a) the "adapted algorithm" [36]; (b) localized blended rainfall algorithm; (c) localized Z-R relation; (d) WSR-88D Z-R relation [13]. The color density represents the observation sample numbers in logarithmic unit.

Discussion
Although the analysis of DSD and polarimetric radar signatures shows similar results, the relatively long distances between disdrometers and the extreme rain centers might induce some uncertainty in the representation of extreme rainfall DSD characteristics. As mentioned, the polarimetric radar signatures at the two rainfall centers are quite different, even for the same precipitation system, which is likely due to the complex falling processes as a result of orographic enhancement and dynamic cloud microphysics involved in this extreme event [44]. More disdrometers and in situ measurements would be required to fully resolve the three-dimensional structure of precipitation in such complex terrains.
In addition, although the polarimetric radars could provide more insights into the two extreme rain regions, the precipitation estimated by the current radar algorithms underestimate the rainfall compared to the ground gauges. The differences of sample areas between radars and gauges may be one reason, especially when the gauges are far from the radars. In such cases, a network of short-range X-band polarimetric radars would be useful for better QPE. Additionally, the adapted algorithms have relatively better performance mainly because they are derived using DSD and gauge data collected during many storm events in Southern China [36]. However, at the same time, the DSD observed at both Huidong and Zhuhai stations during this event showed a lower number concentration compared with previous long-term studies in Yangjiang [43]. This reminds us that the rainfall algorithms should be appropriately developed based on local rainfall characteristics, which is still under investigation.

Conclusions
An epic flood event occurred in southern China from 27 August to 1 September 2018, with a maximum accumulative rainfall of 1394.6 mm recorded at Gaotan station. Multi-observations including reanalysis data, soundings, gauges, disdrometers, and polarimetric radars were used to comprehensively analyze this extreme precipitation event. Based on the reanalysis and gauge data, the synoptic environment and precipitation pattern were investigated. The DSD time series, the mass weighted diameter D m , and normalized intercept parameter N w of Huidong and Zhuhai stations were derived to achieve a better understanding of the precipitation microphysics. It was found that the depression coupled with land-sea interaction and orographic enhancement kept providing moisture to Guangdong province, leading to the extreme rainfall on 29 at Doumen station and on 30 August at Gaotan station. Although both locations are under the same synoptic system, the DSD observed at Huidong station showed a higher number concentration at smaller drop size and lower number concentration at larger drop size compared with the Zhuhai station. This is likely due to the stronger collision caused by sea-inland orographic effect near the Huidong station, which is similar to the microphysical signatures observed by polarimetric radars at Gaotan and Doumen rain regions. Future study is needed to reveal the mechanism of the sea-inland orographic effect on the precipitation microphysics. In addition, the radar based QPE results show that the rainfall algorithms have great influence on the accuracy. Extra work is required to improve local radar rainfall estimates in order to further help with the severe weather warning operations.