Next Article in Journal
The Use of Deep Machine Learning for the Automated Selection of Remote Sensing Data for the Determination of Areas of Arable Land Degradation Processes Distribution
Next Article in Special Issue
Quantitative Precipitation Estimates Using Machine Learning Approaches with Operational Dual-Polarization Radar Data
Previous Article in Journal
The Influence of Camera Calibration on Nearshore Bathymetry Estimation from UAV Videos
Previous Article in Special Issue
Impact of Topography and Rainfall Intensity on the Accuracy of IMERG Precipitation Estimates in an Arid Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Synthetic Quantitative Precipitation Estimation by Integrating S- and C-Band Dual-Polarization Radars over Northern Taiwan

1
Institute for Geosciences, Section Meteorology, University of Bonn, 53121 Bonn, Germany
2
Department of Atmospheric Sciences, National Central University, Taoyuan 32001, Taiwan
3
Central Weather Bureau, Taipei 100006, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 9 November 2020 / Revised: 28 December 2020 / Accepted: 31 December 2020 / Published: 5 January 2021
(This article belongs to the Special Issue Precipitation and Water Cycle Measurements Using Remote Sensing)

Abstract

:
The key factors, namely, the radar data quality, raindrop size distribution (RSD) variability, and the data integration method, which significantly affect radar-based quantitative precipitation estimation (QPE) are investigated using the RCWF (S-band) and NCU C-POL (C-band) dual-polarization radars in northern Taiwan. The radar data quality control (QC) procedures, including the corrections of attenuation, the systematic bias, and the wet-radome effect, have large impact on the QPE accuracy. With the proper QC procedures, the values of normalized root mean square error (NRMSE) decrease about 10~40% for R(ZHH) and about 5~15% for R(KDP). The QPE error from the RSD variability is mitigated by applying seasonal coefficients derived from eight-year disdrometer data. Instead of using discrete QPEs (D-QPE) from one radar, the synthetic QPEs are derived via discretely combined QPEs (DC-QPE) from S- and C-band radars. The improvements in DC-QPE compared to D-QPE are about 1.5–7.0% and 3.5–8.5% in R(KDP) and R(KDP, ZDR), respectively. A novel algorithm, Lagrangian-evolution adjustment (LEA), is proposed to compensate D-QPE from a single radar. The LEA-QPE shows 1–4% improvements in R(KDP, ZDR) at the C-band radar, which has a larger scanning temporal gap (up to 10 min). The synthetic LEA-QPEs by combining two radars have outperformed both D-QPEs and DC-QPEs.

Graphical Abstract

1. Introduction

Accurate radar-based quantitative precipitation estimation (QPE) has been one of the longstanding goals of meteorological radar. Marshall and Palmer [1] utilized horizontal reflectivity (ZHH, mm6m−3) and a power-law relation, ZHH = aRb (Z-R) obtained from simulated radar variables based on measured raindrop size distribution (RSD), to estimate rainfall rate (R, mm h−1). However, the Z-R relation varies vastly in convective, stratiform precipitation, and different climatological regions due to the natural variability in RSD [2]. Seliga and Bringi [3] proposed dual-polarization (dual-pol) radar, which is capable of transmitting horizontal and vertical electromagnetic signals. Additional dual-pol variables from a dual-pol radar, such as the differential reflectivity (ZDR) and specific differential phase (KDP), are utilized in QPE. Consequently, the QPE has significantly been improved by the better quality of radar data and the inclusion of the RSD information [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24].
Diverse forms of power-law QPE relations were recommended using single to multiple dual-pol radar parameters (R(ZHH, ZDR), R(KDP), R(KDP, ZDR), R(ZHH, KDP), and R(ZHH, KDP, ZDR)); [6,7,8,11,12,25]. These QPE algorithms have shown pronounced improvements compared to the Z-R relation. More sophisticated QPE algorithms were explored to overcome various QPE issues. For instance, the variational-based algorithm is developed for attenuation correction and QPE simultaneously [13,19,26]. Moreover, the algorithm utilizing specific attenuation (i.e., R-A) has demonstrated its superiority in QPE as well [20,21,22,23,24].
Nevertheless, the dual-pol radar measurement uncertainties contribute to QPE error significantly [11,27]. The partial beam blockage (PBB), inadequate correction for attenuation, calibration of ZHH and ZDR, and the wet radome effect (WRE) are crucial to QPE application [11,18,19,27,28,29,30]. Furthermore, the variability of RSD also degrades the QPE performance even though dual-pol radar variables are applied. The QPE can be improved by using appropriate coefficients in the power-law relationships derived from local RSD measurements [8,27]. The RSD characteristics of Taiwan in the subtropical region are diverse among various precipitation systems. Chang et al. [31] has investigated the unique RSD of typhoon systems using a 2D-Video disdrometer. Lee et al. [32] further utilized a ten-year JWD disdrometer dataset to study the seasonal variability of RSD in northern Taiwan. The mean values of mass-weighted diameter (Dm) vary from 1.5 mm of summer to 0.97 mm of winter, and the standard deviations of Dm are about 0.34–0.49 mm. It is essential to investigate the impacts of radar data quality control (QC) procedures (e.g., attenuation effect, the WRE) and unique RSD on Taiwan’s QPE.
Taiwan’s radar network consists of three Doppler single-pol S-band radars, one dual-pol S-band radar, and eight dual-pol C-band radars [29]. These radars have two primary purposes: weather surveillance and QPE; thus, different scanning strategies with various spatiotemporal resolutions are applied accordingly [33]. Regarding the QPE, each radar has an individual scheme for operational QPE [29]. The mosaicked rainfall rate field of the whole Taiwan domain is generated from instantaneous QPE from each radar via the distance and height-weighted mean scheme. In order to further improve the accuracy of QPE from multiple radars, optimization of the QC procedures and QPE coefficients based on the radar frequencies and Taiwan RSD characteristics is needed.
Furthermore, synchronizing the QPE products from a network of multiple radars with different temporal resolutions and scanning sequences is very challenging in practice. The precipitation systems evolving vastly over one radar scan period may degrade the accuracy of accumulated rainfall estimation. Instead of deriving the accumulated rainfall by discretely integrating the individual QPE of each radar scan, a new integration method considering the movement and evolution of the precipitation system is required. The study from Ventura et al. [15] utilized the advection correction to reduce the error from low temporal resolution. Chen and Chandrasekar [17] used a Piecewise Cubic Hermite Interpolating Polynomial (PCHIP)-based interpolation methodology. Nevertheless, neither the advection correction nor the PCHIP method considers the evolution of precipitation systems.
This study aims to investigate the impact of critical factors on dual-pol QPE utilizing S- and C-band radars in northern Taiwan. The essential factors discussed here include the data QC procedures, RSD variability, and the QPE integration method. The final goal is to obtain the optimal composite rainfall estimation concurrently using two (or more in the future) different radars. Section 2 shows the radar data and related QC procedures. The RSD data, power-law QPE algorithms, and QPE integration method will be introduced in Section 3. The results of various QPEs will be examined in Section 4. Section 5 provides a brief summary, and Section 6 closes with a conclusion.

2. Quality Control Procedures for S- and C-Band Dual-Polarization Radars

RCWF (radar code of Wu-Fanshan, NEXRAD WSR-88D S-band dual-pol radar, Ryzhkov et al. [9] from the Central Weather Bureau (CWB) and NCU C-POL (National Central University C-band polarimetric radar) located in northern Taiwan are exploited in this study. Their basic information and locations are shown in Table 1 and Figure 1, respectively. RCWF was upgraded to dual-pol radar around March of 2014 and remained operating until the Soudelor typhoon damaged the radome on 8 August in 2015. NCU C-POL is a research radar upgraded from a decommission Ericsson Doppler radar and has been utilized to study the RSD characteristics of typhoon systems associated with Taiwan’s complex topography [31]. Both radars have the same signal processor (Sigmat RVP8), yet they have different configurations. As shown in Figure 1, these two radars have broad overlapping coverage in northern Taiwan. In addition, different frequencies are ideal to evaluate the performance of dual-pol QPE at different wavelengths.
Both RCWF and NCU C-POL have various data quality issues. Different distributions of PBB depending on the radar locations are displayed in Figure 1. Distinct ground clutter suppression configurations from different manufacturers induce distinct ground clutter characteristics. The degree of the attenuation effect varies according to radar wavelengths. Each radar has distinct characteristics of systematic ZHH and ZDR biases due to different maintenance and calibration procedures (e.g., transmitter and receiver). Finally, different radome materials and radar frequencies cause distinct WREs. The QC procedures for both radars shown in Figure 2 (left panel) have been applied before the QPE retrievals. These procedures include the PBB and non-meteorological signal removal, the corrections of attenuation, the systematic bias, and the WRE.

2.1. PBB Removal

Radar data with the PBB effect, mainly resulting from Taiwan’s complex topography, are first identified by a beam blockage simulation [10,18,34,35]. When the central height of the radar beams, estimated by 4/3 of the Earth’s radius model, is lower than the top of the terrain (50% of the power is blocked), the data and the following bins along the radial will be removed.

2.2. Non-Meteorological Signal Removal

Despite the Doppler spectrum clutter suppression has been applied to both RCWF and NCU C-POL radars, some residual non-meteorological signals remain after the PBB removal. The standard deviation values of the differential propagation phase shift ϕ D P of five radar bins along a radial (Std. ϕ D P ) and the correlation coefficient (CC) are subsequently used to filter out non-meteorological signals. The data with the values of Std. ϕ D P higher than 15° or with CC less than 0.85 is removed for RCWF. The corresponding threshold values of Std. ϕ D P and CC for NCU C-POL are 20° and 0.8.
The threshold values of Std. ϕ D P and CC for RCWF and NCU C-POL were obtained accordingly by examining large amounts of radar data. The values which preserve the most meteorological signal and reduce the most persistent ground clutters were selected as the threshold. The same technique has been applied to the operational Doppler dual-pol radars in Taiwan [33].

2.3. Attenuation Correction

The intrinsic ZHH ( Z HH intrinsic ) and ZDR ( Z DR intrinsic ) can be obtained after attenuation and bias corrections from the observational ZHH ( Z HH obs . ) and ZDR ( Z DR obs . ) as bellow,
Z HH intrinsic = Z HH obs . + A H + Z HH   bias   and
Z DR intrinsic = Z DR obs . + A HV + Z DR bias .
The attenuation effects of ZHH and ZDR on both RCWF and NCU C-POL are corrected by the ϕ D P -based algorithm. The horizontal attenuation (AH, dB) and differential attenuation (AHV, dB) can be estimated as follows,
A H = α Δ ϕ D P   and
A HV = β Δ ϕ D P .
The coefficients of α and β for S-band (C-band) radar are 0.0151 dB deg−1 (0.0727 dB deg−1) and 0.0025 dB deg−1 (0.0161 dB deg−1), respectively. These coefficients were obtained from the long-term RSD measurements in northern Taiwan. The RSD data and implementation of coefficient calculations will be discussed in Section 3.1 and Section 3.2, respectively.

2.4. ZHH Systematic Bias and the WRE Correction

The next step is to correct ZHH biases including the radar systematic bias (i.e., the miscalibration of transmitter and receiver) and the WRE. The self-consistency among the ZHH, ZDR, and KDP measurements of rain has been utilized to estimate the radar systematic bias [36]. The WRE is an additional bias due to the attenuation effect of the water coated over the radome while a precipitation system is located over the radar site. It has similar characteristics to the radar systematic bias. Therefore, the ZHH biases which consist of both radar systematic bias and WRE can be calculated via the self-consistency method. The radar systematic bias is obtained by averaging the ZHH biases of no-WRE events (identified by the rain gauge measurements next to the radar site). Hence, the attenuation of ZHH from the WRE can be derived by examining the excess of estimated ZHH biases from the systematic bias. On the other hand, the WRE of ZDR measurements depends on the forms in which the rain accumulates on the radome [37,38,39,40,41]. The forms could be droplets, rivulets, or non-uniform film and cause diverse WREs of ZDR. It is difficult to estimate the WRE of ZDR, and there is no practical technique either. Consequently, ZDR is also excluded from the self-consistency algorithm for the ZHH bias estimation [36] to avoid WRE-contaminated ZDR.
Here, the attenuation-corrected ZHH ( Z HH corr . , in mm6m−3) data of rain is applied to the KDP-ZHH relation to obtain the value of K D P in each radar bin as shown below,
K D P = a Z HH corr . b .
Based on the RSD measurements in Section 3.2, the coefficients a and b are 5.3 × 10 5 ( 2.5 × 10 4 ) and 0.88 (0.81) for S-band (C-band) radar. The calculated K D P of each radar bin is then integrated along a radial to derive the increments of differential phase shift ( Δ Φ D P ) as
Δ Φ D P = 2 K D P d r = 2 a Z HH corr . b d r .
dr is the radial resolution of radar data [36]. Since the phase measurement (i.e., Δ Φ D P ) is not affected by the absolute calibration of the radar system and attenuation [11,36], the ZHH bias ( Z HH bias ) including both the radar systematic and the WRE biases can be estimated via the calculated Δ Φ D P and the measured Δ Φ D P as
Z HH bias = 1 b l o g 10 ( Δ Φ D P Δ Φ D P ) .
In practice, the last valid Δ Φ D P and Δ Φ D P data along a radial below the freezing level (avoid non-liquid phase data) are used to calculate radial-wise Z HH bias . The values of Δ Φ D P is higher than Δ Φ D P when ZHH is over-calibrated (i.e., positive ZHH bias), and vice versa. Then the scan-wise representative ZHH bias is determined via the mean value of radial-wise Z HH bias within one scan for the bias correction. Only the WRE of ZHH was considered via a self-consistency technique. The WRE of ZDR was not corrected in this study.

2.5. ZDR Systematic Bias Correction

The last remaining radar measurement error in the QC processes is ZDR systematic bias. Since both radars cannot perform vertical pointing scan (birdbath scan), the ZDR systematic bias is estimated by statistical analysis. The stratiform rain consisting of small sizes of raindrops mostly has ZDR values slightly above 0 dB. This ZDR characteristic of light rain is used for ZDR bias estimation. The mean ZDR averaged from the RSD-simulated ZDR of light rain (i.e., 0.19 dB), mainly the data with ZHH values from 10 to 20 dBZ, is considered as a reference. The difference between the mean value of radar-measured ZDR in the light rain region and the reference value is thus reasonably viewed as ZDR systematic bias. For example, the mean ZDR value in the light rain region from NCU C-POL was −0.28 dB, then the ZDR bias was determined as −0.47 dB. The ZDR biases of NCU C-POL were derived case by case and the values were within the range of −0.21 to −0.48 dB in this study. RCWF is a well-calibrated radar; therefore, ZDR calibration was not applied.

3. Quantitative Precipitation Estimation from S- and C-Band Dual-Polarization Radars

The RSD data obtained from the disdrometer was used to derive the coefficients for the aforementioned QC procedures and dual-pol QPE algorithms. The procedures of processing RSD data, deriving the coefficients of QPE relations, and developing different QPE (decision-tree) products are introduced. In addition, two integration methods intending to improve QPE performance are also described in this section.

3.1. Disdrometer Data

In this research, the RSD data collected by a 2D-Video Disdrometer (2DVD, [42]) located at NCU from October of 2000 to June of 2007 is used to characterize the RSD in northern Taiwan. The 2DVD records the particle diameter (D) and the terminal velocity (Vt) of each raindrop. The data is quality controlled by the Vt–based filter technique [31,43,44]. The 6-min RSD is then calculated to ensure sufficient raindrop sampling numbers of each RSD [19,31,32,45]. Moreover, the rainfall rates of less than 1 mm h−1 are removed to eliminate inadequate RSD data set. There are a total of 14,314 quality-controlled minutely RSDs available for analysis.
Chen and Chen [46] has investigated the characteristics of Taiwan’s precipitation systems and classified them into five distinct types. They are spring rain from March to April, Mei-Yu from May to June, summer convection from July to September, winter cold front from October to February, and typhoon systems. Most of the systems can be approximated simply by month, except typhoon is case-selected. Lee et al. [32] utilized the same classification and demonstrated unique RSD characteristics of these five types of precipitation systems. The bigger mean raindrop size (i.e., mass-weighted diameter) from Mei-Yu, summer convection, and typhoon events is found compared to spring rain and winter cold front.
The NCU 2DVD data is thus classified according to Chen and Chen [46]. The rainfall intensities of 6-min RSD for each precipitation type are summarized in Figure 3. The majority of rainfall intensities are less than 15 mm h−1. The maximum rainfall rates of all types of precipitation are around or above 70 mm h−1. The Mei-Yu and typhoon season have similar maximum rainfall up to 90 mm h−1. The rainfall rates of the cold front are mostly below 5 mm h−1. These quality-controlled and classified RSD data are consequently used to derive the coefficients of various QPE relations.

3.2. QPE Coefficients

Various forms of power-law relationship using one to two dual-pol parameters were proposed to improve QPE, namely R(KDP) [6,11,25], R(ZHH, ZDR) [8,12], and R(KDP, ZDR) [7,8]. The conventional Z-R and the most common dual-pol relations used in this study are shown in the following,
R ( Z HH ) = a 1 Z HH b 1 ,
R ( Z HH , Z DR ) = a 2 Z HH b 2 Z DR c 2 ,
R ( K DP ) = a 3 K DP b 3   and
R ( K DP , Z DR ) = a 4 K DP b 4 Z DR c 4 .
Lee [27] has shown that the RSD variability is one of the significant error sources in radar-based QPE. Thus, the correspondingly precipitation-type coefficients in dual-pol QPE relations are derived for reducing QPE error from the seasonal RSD variabilities [32]. First, the S- and C-band radar variables (i.e., ZHH, ZDR, KDP, AH, AHV) are computed through the T-matrix scattering calculation [47] using the measured 6-min RSDs data set from NCU 2DVD, with the assumption of raindrop axis-ratio proposed by Brandes et al. [48] and Chang et al. [31] to non-typhoon cases and typhoon cases, respectively. Then the coefficients of Equations (8)–(11) for five types of precipitations (hereafter, seasonal coefficients) are derived by fitting the simulated radar variables and RSD-derived rainfall rates via the Levenberg–Marquardt algorithm [49] at both frequencies. Diverse coefficients due to the natural RSD variability among precipitation systems can be noticed. The coefficients based on the entire database (hereafter, all-season coefficients) are also obtained for comparison (Table 2 and Table 3). The same procedures is also utilized to derive the coefficients of the attenuation correction and self-consistency relationship in Equations (3)–(5).
The seasonal coefficients are expected to diminish the RSD variations by considering the seasonal RSD characteristics, thus outperform all-season coefficients. The benefit of optimized QPE coefficients is first examined via the simulated data. The estimated rainfall (Rest) using seasonal and all-season coefficients, respectively, are derived from the RSD-simulated radar variables, and its performance is validated against RSD-derived rainfall rates (Rt). Since radar variables and Rt are both based on the same database, namely, without radar measurement errors, the influence of the RSD variability on QPE can be investigated solely. The normalized bias (NBIAS) and normalized root mean square error (NRMSE) are calculated for validation,
NBIAS =     ( R t R est ) N R t ¯   and
NRMSE = (   ( R t R est ) 2 N ) 0.5 R t ¯ .
N is the number of RSD.
The results in Figure 4 (upper panel) show that the values of NBIAS of QPEs using seasonal coefficients are closer to 0 than using all-season coefficients at S band, while C-band QPEs have comparable results. The NRMSEs of QPEs using seasonal coefficients are consistently lower than using all-season coefficients (Figure 4 lower panel). The improvements from applying seasonal coefficients are most pronounced in R(KDP) and less evident in R(ZHH). It can be explained by the fact that the seasonal variability of Taiwan RSD is mostly owing to small to medium sizes of raindrops [32] by which KDP is dominated. Moreover, the KDP-based algorithms outperform ZHH-based algorithms at both C and S bands. R(KDP, ZDR) has the lowest NRMSE due to the additional RSD information from ZDR. The results show a good agreement with Bringi and Chandrasekar [11] and Lee [27]. Even though dual-pol radar variables have been utilized, applying seasonal coefficients in QPE relations further reduces the QPE uncertainty by diminishing the RSD variability.

3.3. Simplified Decision-Tree QPE

In practice, various dual-pol radar measurement errors cause a decline in QPE performance [27]. For instance, the noisy ZDR due to a fast-scanning radar with a low sampling rate [11,28] severely affects ZDR-based QPEs. The range-derivative KDP estimated from ϕ D P requires various filtering techniques [8,11,50,51,52,53] and introduces more uncertainty in KDP-based QPEs. In the case of light rain, KDP is too small to be estimated with sufficient accuracy. The combination of different dual-pol QPE algorithms using decision-tree logic was thus proposed [8,9,14,16,54]. The superiority of this composite method over a single dual-pol algorithm has been proved [8].
Most of the decision-tree algorithms use more than two dual-pol QPE relationships. For example, studies from [14,16,54] utilize three dual-pol QPE relationships including R(Z), R(ZHH, ZDR), and R(KDP, ZDR) for their final QPE products. In this study, however, each dual-pol QPE algorithm is combined with R(Z) only. The purpose of utilizing simplified decision-tree algorithms is to examine the performance of specific dual-pol QPE relationships solely. In Figure 2 (right panel), the decision-tree algorithms are applied to quality-controlled plan position indicator (PPI) radar data. An additional filter of ZHH less than 10 dBZ, whose contribution to accumulated rainfall is neglectable, is applied. Each algorithm has its own criteria to assign different QPE relationships (Step C). The criteria were determined theoretically based on the dual-pol radar principle [11] and empirically adopted from previous researches. The concept is to avoid using inaccurate KDP in light rain and negative-biased ZDR due to the WRE in pure rain. For clarity, the simplified decision-tree QPE combining Z-R with R(ZHH, ZDR), R(KDP) or R(KDP, ZDR) are referred to R(ZHH, ZDR), R(KDP) or R(KDP, ZDR) QPEs, hereafter. To derive the composite QPE, the PPI QPE is then interpolated into the Cartesian coordinate. The lowest available Cartesian QPE is implemented as a final composite QPE product. These procedures are similar to Kwon et al. [18].

3.4. Integration Methods for Radar-Based QPE

Conventionally, the radar-based QPE is obtained by integrating each radar scan discretely, hereafter D-QPE ( Δ T i : time difference between two scans). The D-QPE of rapidly-evolving severe precipitation systems suffers from an inadequate radar scanning rate. Some researches proposed different solutions, namely the advection correction [15] and PCHIP interpolation [17]. Nevertheless, the influence of the evolution of the precipitation systems on QPE has not been well investigated. As aforementioned the time resolutions of RCWF and NCU C-POL were not synchronized and varied between 5 to 10 min (Table 4). To overcome this issue, the D-QPEs from asynchronous radars are “discretely” combined (hereafter, DC-QPE) to increase temporal resolution. The DC-QPE then can better reveal the evolution of precipitation than D-QPE does. Figure 5 displays an example that the NCU C-POL (red circles) and RCWF (blue crosses) D-QPEs are "discretely" combined whenever it is available (marked by black triangles). Hence, the DC-QPE has more sampling numbers.
However, not every location can be covered by two radars or more. A novel algorithm, namely Lagrangian-evolution adjustment (LEA), is developed for improving D-QPE/DC-QPE in this study. The QPEs from R T i and R T i + 1 are first used to derive the advection speed and direction of precipitation movement via tracking radar echo based on the cross-correlation (TREC, [15,55]). Thus, the R T i can be forward advected to T i + 1 ( = T i + Δ T i ), and the R T i + 1 can be backward advected to T i . The evolution of the precipitation system between T i and T i + 1 is subsequently estimated via time-weighted linear interpolation. The linear interpolated QPE ( R T i + Δ t j ) at T i + Δ t j is derived as follows,
R T i + Δ t j = R T i Δ T i Δ t j Δ T i + R T i + Δ T i Δ t j Δ T i .
The interval of Δ t j is defined as one minute (Figure 5 green dots) in this study. The seamless LEA-QPEs with the considerations of the advection and evolution of the precipitation systems are implemented between T i and T i + 1 .

4. QPE Comparison between C- and S-Band Radar

Eighteen heavy rainfall events for real case studies were selected by examining the rain gauge measurements in northern Taiwan (Table 4). Most of the events are from the Mei-Yu season (East Asia rainy season, from mid-May to mid-June for Taiwan), since it contributes to most of the rainfall in Taiwan [46]. The most intensive hourly rainfall rate and highest event-accumulated rainfall are 103.5 mm h−1 and 167.0 mm on 14 June, 2015. RCWF and NCU C-POL have diverse scanning strategies and configurations. RCWF has higher azimuthal and temporal resolution compared to NCU C-POL (Table 1 and Table 4). The advantage of NCU C-POL is its location at a relatively lower altitude (156 m) than RCWF (766 m), providing closer observations to the surface and less vertical discrepancy than RCWF. QPE performances are evaluated against operational tipping bucket rain gauge measurements which are collected and maintained by CWB (black dots in Figure 1). For clarity, the S- and C-band QPEs will refer to RCWF and NCU C-POL QPEs hereafter.

4.1. The Influence of QC Procedures on C- and S-Band QPE

As discussed in previous studies [11], the KDP-based QPEs are immune to the systematic bias or attenuation effect and are expected to have better performance than ZHH-based QPEs. Nevertheless, ZHH-based QPEs are still applied in the simplified decision-tree QPE algorithms for specific purposes (e.g., avoid noisy KDP in light rain). Thus, it is essential to examine the impacts of various measurement errors from attenuation and the WRE on QPE.
To avoid WRE-contaminated ZDR error in the analysis, only R(ZHH) and R(KDP) are examined for comparison. As shown in Figure 6, both the C-band R(ZHH) and R(KDP) QPEs have higher values of NRMSE than the S-band when only the systematic bias is corrected. These results can be expected because of less attenuation and the WRE at longer wavelength radars. The values of NRMSE decrease dramatically at C band after the ensuing attenuation correction (brown areas), and the improvements are more pronounced than at S band. The results indicate that applying a proper attenuation correction can significantly mitigate the impacts of the attenuation effect on the C-band QPEs. The C- and S-band QPEs have a comparable performance after the systemic bias and attenuation corrections.
Removing the WRE further reduces QPE errors (yellow areas). The most pronounced reduction of NRMSE can be found in the C-band R(ZHH) and R(KDP) QPEs. The C-band radar is more vulnerable to the weakening of the backscattering power due to rain over the radome. The C-band QPEs show slightly better results compared to S-band after the complete data QC procedures. The ZHH correction here is also beneficial to the R(KDP) algorithms since ZHH used in light rain is corrected, and more accurate ranges of ZHH are identified for the decision-tree criterion.

4.2. C- and S-Band Seasonal and All-Season Coefficients in QPE Algorithms

The dual-pol QPEs using seasonal and all-seasonal coefficients, respectively, in real case are first examined by comparing their NBIAS (Figure 7 upper panel). The QPEs using seasonal coefficients are consistently less biased than the ones using all-season coefficients. In terms of NRMSE (lower panel), C-band QPEs of seasonal coefficients show lower errors than all-seasonal coefficients, except for the R(ZHH, ZDR) algorithm. For the S-band QPEs, the seasonal ZDR-based algorithms have slightly higher NRMSE values compared to the all-season. It is noticed that both C- and S-band R(ZHH, ZDR) algorithms using the seasonal coefficients perform much worse. The biased ZDR measurement due to the uncorrected WRE (visually examined, not shown in the paper) is postulated as the main reason.
However, the R(KDP, ZDR) algorithms are less affected by ZDR than R(ZHH, ZDR). It can be explained by the more negative (more weighted) seasonal exponents of ZDR in the R(ZHH, ZDR) relations, and the less negative (less weighted) seasonal exponents of ZDR in R(KDP, ZDR). The sophisticated seasonal R(ZHH, ZDR) algorithms are sensitive to the ZDR bias in the real case. The R(ZHH, ZDR) algorithms using seasonal coefficients at S-band shows less deterioration due to the relative insensitivity of ZDR to the WRE. Overall, R(KDP) shows the most improvements after using seasonal coefficients, consistent with the result from the theoretical verification in Section 3.2.

4.3. C- and S-Band Discrete QPE Comparison (D-QPE)

The C- and S-band D-QPEs are first compared individually at different rainfall intensities to identify the best algorithm. The NRMSE (or NBIAS) as a function of rainfall intensities are derived by gradually increasing rainfall intensity thresholds in Equations (12) and (13); namely, only the data of Rt > rainfall intensity threshold are included for the calculation. In Figure 8, the R(KDP, ZDR) algorithm shows the lowest NRMSE value over other algorithms consistently at different rainfall intensities at both C- and S-band radars. The KDP-based algorithms mostly outperform ZHH-based algorithms. The R(ZHH, ZDR) algorithm performs worst at C-band due to its higher sensitivity of ZDR to the WRE, and R(ZHH) has the highest values of NRMSE at S-band.
Subsequently, various D-QPEs from two radars are compared in Figure 9. Except for the R(ZHH, ZDR) algorithm, the C-band D-QPEs are slightly better than S-band D-QPEs, which show smaller errors in heavy rain. Better results from the C-band radar can be concluded for two reasons. First, the C-band radar is located at a relatively lower altitude and provides radar measurements closer to the surface. Second, the KDP measurement is more sensitive to rainfall rate at C band. Beware that the C-band QPEs are obtained with a longer radar scan period, otherwise the advantages of the C-band radar would be more noticeable. As the comparable QPEs from S- and C-band radar have been obtained after the complete QC procedures and seasonal coefficients are applied, the integration method is investigated in the next section.

4.4. Discretely Combined C- and S-Band QPEs (DC-QPE)

The results of DC-QPEs (Figure 10 black lines) show that both R(KDP) and R(KDP, ZDR), which are the best D-QPEs and thus chosen for further examination, have pronounced improvements. The DC-QPEs outperform both S- and C-band D-QPEs, especially in R(KDP, ZDR), with lower values of NRMSE in all rainfall intensities. The reductions of NRMSE from D-QPEs to DC-QPEs are about 1.5–7.0% in R(KDP) and about 3.5–8.5% in R(KDP, ZDR). The encouraging improvements confirm the fact that radar scanning strategy (i.e., temporal resolution) does play a crucial role in radar-based QPE. Integrating asynchronous D-QPE products can reduce the inaccuracy dramatically by increasing the radar sample number.

4.5. Lagrangian-Evolution Adjustment (LEA) QPE

Figure 11 demonstrates two C-band QPE maps derived from R(KDP, ZDR) at 06:54 UTC ( R T i ) and 07:04 UTC ( R T i + 1 ) on 19 August in 2014. The temporal resolution ( Δ T i ) in this case is 10 min. The convective core (R > 100 mm h−1) and precipitation area (R > 10 mm h−1) shifted northward noticeably within 10 min. It also clearly shows the precipitation system not only moved fast but also evolved vastly. The heavy rainfall area had expanded, and its intensity had been enhanced. An example of LEA-QPE of Δ t j = 5 min is shown in Figure 11c. The convective core (blue ⨁) is located between the cores of R T i and R T i + 1 (white and black ⨁). The movement and evolution of the precipitation structure are also noticeable. There are nine additional LEA-QPEs obtained within these 10 min using the LEA technique.
The C-band radar, which has larger values of Δ T than the S-band, shows 1–4% reductions of NRMSE after the LEA algorithm is applied (Figure 12). On the other hand, S-band LEA-QPEs have little improvement with the NRMSE values remaining similar to D-QPEs. Moreover, the C-band LEA-QPEs have consistently lower values of NRMSE than the S-band LEA-QPEs. These results indicate that the LEA algorithm has more positive feedback in the case of large Δ T . The most improvements are found in rainfall rates above 40 mm h−1 in R(KDP, ZDR). It is postulated that the evolutions of the precipitation systems can be mainly revealed by the RSD variations. Hence, the LEA algorithm can further improve ZDR-based QPEs by including the RSD evolutions from ZDR measurements.
As discussed in the previous section, the DC-QPE outperforms the D-QPE by combining the S- and C-band D-QPEs. Applying the LEA technique compensates for the disadvantage of relatively lower temporal resolution at the C-band radar. In Figure 13 (upper panel), the QPEs derived from R(KDP) show a consistent underestimation of the rainfall rate, and ones from R(KDP, ZDR) are less biased. When the LEA technique is applied to DC-QPE (LEAC-QPE), the LEAC-QPEs are slightly better than DC-QPEs for R > 0 mm h−1. The difference between LEAC-QPEs and DC-QPEs diminishes as the rainfall intensity increases. NRMSE in the lower panel of Figure 13 shows that the DC-QPEs outperform both the S- and C-band LEA-QPEs with lower NRMSE. LEAC-QPEs show nearly identical NRMSE values to DC-QPEs in both R(KDP) and R(KDP, ZDR). Overall, the LEAC-QPEs have less biased NBIAS than DC-QPEs and the lowest NRMSE as DC-QPEs.

5. Summary

The dual-pol QPEs of RCWF and NCU C-POL radars have been investigated comprehensively in this study. Three key factors influencing the QPE accuracy, namely, the radar data quality, RSD variability, and the QPE integration method, were examined. Eighteen heavy rain events over northern Taiwan were selected for QPE validation in this work.
The QC procedures, including the corrections of attenuation, the systematic bias, and the WRE for S- and C-band radars over northern Taiwan, have been documented. Different levels of ZHH corrections were applied for QPEs for 18 heavy precipitation events. The results show that attenuation and the WRE play crucial roles in QPE, even for S-band radars. The proper QC procedures significantly reduce the QPE uncertainty. The C-band QPE, in the end, performs slightly better than the S-band partially due to the lower altitude of the radar-site. The R(KDP, ZDR) algorithm outperforms other QPEs.
The all-season and seasonal coefficients in four power-law QPE relations derived from eight-year 2DVD disdrometer data were examined to investigate the impact of seasonal RSD variability on QPE. The results indicate that the seasonal coefficients pronouncedly reduce the error by diminishing the RSD variability. The R(ZHH, ZDR) algorithm has higher values of NRMSE due to the WRE of ZDR measurements.
With the proper QC procedures and seasonal coefficients, the comparable QPEs from S- and C-band radars are thus ready for further composition. Even though these QPEs were from asynchronous S- and C-band radar scanning with different temporal resolutions, two QPE integration methods have been carried out to increase the radar sample number. The D-QPEs from C- and S-band radars were first combined discretely. DC-QPE has improved the QPE accuracy significantly with the reductions of NRMSE about 1.5–7.0% in R(KDP) and about 3.5–8.5% in R(KDP, ZDR) from D-QPE to DC-QPE. However, the D-QPE/DC-QPE may still miss the motion and evolution of the precipitation systems. A newly developed algorithm, namely Lagrangian-evolution adjustment (LEA), is proposed in this study to further improve the QPE performance. The advection of precipitation systems is estimated via a tracking technique, and the evolution is derived via time-weighted linear interpolation. The LEA-QPE has shown a noticeable improvement in R(KDP, ZDR) at the C-band radar, which has larger scanning temporal gaps (up to 10 min). NRMSE reduction after the LEA algorithm is about 1–4%. Further combining the DC-QPE and LEA-QPE, namely LEAC-QPE, shows little improvement.

6. Conclusions

The individual dual-pol QPE can be significantly improved by two key factors: (1) improving the data quality by applying proper radar data quality procedures, including attenuation, the radar systematic bias and the WRE corrections, and (2) reducing the influence of the RSD variability on dual-pol QPE by applying seasonal coefficients. Consequently, the comparable QPEs from asynchronous S- and C-band radar scans can be composited using more sophisticated integration methods, such as discretely combined QPE and Lagrangian-Evolution Adjustment (LEA) proposed in this study. The synthetic LEA-QPEs by combining S- and C-band dual-pol radars have outperformed conventional QPEs.

Author Contributions

The work presented here was carried out in collaboration with all authors. J.-Y.C. performed data processing, developed algorithms, produced the manuscript; W.-Y.C. helped design the research procedures and analyze the data. P.-L.C. provided the data and participated in discussions on the analyses; J.-Y.C. and W.-Y.C. contributed to the reviewing and revising of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Ministry of Science and Technology, Taiwan: MOST109-2625-M-008-011

Acknowledgments

We thank all of the staff of the Central Weather Bureau of Taiwan for providing valuable radar and rain gauge data used in this paper.

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Marshall, J.S.; Palmer, W.M.K. The distribution of raindrops with size. J. Meteorol. 1948, 5, 165–166. [Google Scholar] [CrossRef]
  2. Radar observation of the atmosphere. L. J. Battan (The University of Chicago Press) 1973. PP X, 324; 125 figures, 21 tables. £7·15. Q. J. R. Meteorol. Soc. 1973, 99, 793. [CrossRef]
  3. Seliga, T.A.; Bringi, V.N. Potential Use of Radar Differential Reflectivity Measurements at Orthogonal Polarizations for Measuring Precipitation. J. Appl. Meteorol. 1976, 15, 69–76. [Google Scholar] [CrossRef] [Green Version]
  4. Seliga, T.A.; Bringi, V.N.; Al-Khatib, H.H. A Preliminary Study of Comparative Measurements of Rainfall Rate Using the Differential Reflectivity Radar Technique and a Raingage Network. J. Appl. Meteorol. 1981, 20, 1362–1368. [Google Scholar] [CrossRef] [Green Version]
  5. Gorgucci, E.; Chandrasekar, V.; Scarchilli, G. Radar and Surface Measurement of Rainfall during CaPE: 26 July 1991 Case Study. J. Appl. Meteorol. 1995, 34, 1570–1577. [Google Scholar] [CrossRef]
  6. Gorgucci, E.; Scarchilli, G.; Chandrasekar, V. Specific Differential Phase Estimation in the Presence of Nonuniform Rainfall Medium along the Path. J. Atmos. Ocean. Technol. 1999, 16, 1690–1697. [Google Scholar] [CrossRef]
  7. Ryzhkov, A.V.; Zrnić, D.S. Comparison of Dual-Polarization Radar Estimators of Rain. J. Atmos. Ocean. Technol. 1995, 12, 249–256. [Google Scholar] [CrossRef] [Green Version]
  8. Ryzhkov, A.V.; Giangrande, S.E.; Schuur, T.J. Rainfall Estimation with a Polarimetric Prototype of WSR-88D. J. Appl. Meteorol. 2005, 44, 502–515. [Google Scholar] [CrossRef] [Green Version]
  9. Ryzhkov, A.V.; Schuur, T.J.; Burgess, D.W.; Heinselman, P.L.; Giangrande, S.E.; Zrnic, D.S. The Joint Polarization Experiment: Polarimetric rainfall measurements and hydrometeor classification. Bull. Am. Meteorol. Soc. 2005, 86, 809–824. [Google Scholar] [CrossRef] [Green Version]
  10. Fulton, R.A.; Breidenbach, J.P.; Seo, D.; Miller, D.A.; O’Bannon, T. The WSR-88D rainfall algorithm. Weather Forecast. 1998, 13, 377–395. [Google Scholar] [CrossRef]
  11. Bringi, V.N.; Chandrasekar, V. Polarimetric Doppler Weather Radar: Principles and Application; Cambridge University Press: Cambridge, UK, 2001; 636p. [Google Scholar]
  12. Illingworth, A.J.; Thompson, R.J. The Estimation of Moderate Rain Rates with Operational Polarisation Radar. In Proceedings of the 32th Conference on Radar Meteorology, Albuquerque, NM, USA, 22–29 October 2005. [Google Scholar]
  13. Hogan, R.J. A Variational Scheme for Retrieving Rainfall Rate and Hail Reflectivity Fraction from Polarization Radar. J. Appl. Meteorol. Climatol. 2007, 46, 1544–1564. [Google Scholar] [CrossRef] [Green Version]
  14. Giangrande, S.E.; Ryzhkov, A.V. Estimation of Rainfall Based on the Results of Polarimetric Echo Classification. J. Appl. Meteorol. Climatol. 2008, 47, 2445–2462. [Google Scholar] [CrossRef]
  15. Ventura, J.F.; Kabeche, F.; Fradon, B.; Hogan, R.; Boumahmoud, A.A.; Illingworth, A.; Tabary, P. Extensive evaluation of Polarimetric Quantitative Precipitation Estimations (QPE) in ideal and less ideal conditions. In Proceedings of the 6th European Conference on Radar in Meteorology and Hydrology, Sibiu, Romania, 6–10 September 2010. [Google Scholar]
  16. 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]
  17. 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]
  18. Kwon, S.; Jung, S.H.; Lee, G. Inter-comparison of radar rainfall rate using Constant Altitude Plan Position Indicator and hybrid surface rainfall maps. J. Hydrol. 2015, 531, 234–247. [Google Scholar] [CrossRef]
  19. Chang, W.; Vivekanandan, J.; Ikeda, K.; Lin, P. Quantitative Precipitation Estimation of the Epic 2013 Colorado Flood Event: Polarization Radar-Based Variational Scheme. J. Appl. Meteorol. Climatol. 2016, 55, 1477–1495. [Google Scholar] [CrossRef]
  20. Ryzhkov, A.V.; Diederich, M.; Zhang, P.; Simmer, C. Potential utilization of specific attenuation for rainfall estimation, mitigation of partial beam blockage, and radar networking. J. Atmos. Ocean. Technol. 2014, 31, 599–619. [Google Scholar] [CrossRef]
  21. Diederich, M.; Ryzhkov, A.; Simmer, C.; Zhang, P.; Trömel, S. Use of specific attenuation for rainfall measurement at X-band radar wavelengths. Part I: Radar calibration and partial beam blockage estimation. J. Hydrometeorol. 2015, 16, 487–502. [Google Scholar] [CrossRef]
  22. Diederich, M.; Ryzhkov, A.; Simmer, C.; Zhang, P.; Trömel, S. Use of specific attenuation for rainfall measurement at X-band radar wavelengths. Part II: Rainfall estimates and comparison with rain gauges. J. Hydrometeorol. 2015, 16, 503–516. [Google Scholar] [CrossRef]
  23. Wang, Y.; Cocks, S.; Tang, L.; Ruzhkov, A.; Zhang, P.; Zhang, J.; Howard, K. A prototype quantitative precipitation estimation algorithm for operational S-Band polarimetric radar utilizing specific attenuation and specific differential phase. Part I: Algorithm description. J. Hydrometeorol. 2019, 20, 985–997. [Google Scholar] [CrossRef]
  24. Zhang, J.; Tang, L.; Cocks, S.; Zhang, P.; Ryzhkov, A.; Howard, K.; Langston, C.; Kaney, B. A Dual-polarization Radar Synthetic QPE for Operations. J. Hydrometeorol. 2020, 21, 2507–2521. [Google Scholar] [CrossRef] [Green Version]
  25. Ruzanski, E.; Chandrasekar, V. Nowcasting Rainfall Fields Derived from Specific Differential Phase. J. Appl. Meteorol. Climatol. 2012, 51, 1950–1959. [Google Scholar] [CrossRef]
  26. Chang, W.; Vivekanandan, J.; Chen Wang, T. Estimation of X-Band Polarimetric Radar Attenuation and Measurement Uncertainty Using a Variational Method. J. Appl. Meteor. Climatol. 2014, 53, 1099–1119. [Google Scholar] [CrossRef]
  27. Lee, G.W. Sources of Errors in Rainfall Measurements by Polarimetric Radar: Variability of Drop Size Distributions, Observational Noise, and Variation of Relationships between R and Polarimetric Parameters. J. Atmos. Ocean. Technol. 2006, 23, 1005–1028. [Google Scholar] [CrossRef]
  28. Melnikov, V.M. Simultaneous Transmission Mode for the Polarimetric WSR-88D; NOAA/NSSL Report. 2004; 84p. Available online: http://www.nssl.noaa.gov/publications/wsr88d_reports/SHV_statistics.pdf (accessed on 1 June 2018).
  29. Chang, P.-L.; Zhang, J.; Tang, Y.-S.; Tang, L.; Lin, P.-F.; Langston, C.; Kaney, B.; Chen, C.-R.; Howard, K. An Operational Multi-Radar Multi-Sensor QPE System in Taiwan. Bull. Am. Meteorol. Soc. 2020. [Google Scholar] [CrossRef]
  30. 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. Ocean. Technol. 2019, 36, 2501–2520. [Google Scholar] [CrossRef]
  31. Chang, W.-Y.; Wang, T.-C.-C.; Lin, P.-L. Characteristics of the raindrop size distribution and drop shape relation in typhoon systems in the western Pacific from the 2D video disdrometer and NCU C-band polarimetric radar. J. Atmos. Ocean. Technol. 2009, 26, 1973–1993. [Google Scholar] [CrossRef]
  32. Lee, M.-T.; Lin, P.-L.; Chang, W.-Y.; Seela, B.K.; Janapati, J. Microphysical characteristics and types of precipitation for different seasons over north Taiwan. J. Meteorol. Soc. Jpn. 2019, 97, 841–865. [Google Scholar] [CrossRef] [Green Version]
  33. Chang, P.-L.; Lin, P.-F.; Jou, J.-D.B.; Zhang, J. An application of reflectivity climatology in constructing radar hybrid scans over complex terrain. J. Atmos. Ocean. Technol. 2009, 26, 1315–1327. [Google Scholar] [CrossRef]
  34. O’Bannon, T. Using a ‘terrain-based’ hybrid scan to improve WSR-88D precipitation estimates. In Proceedings of the 28th Conference on Radar Meteorology, Austin, TX, USA, 7–12 September 1997; American Meteorology Society: Boston, MA, USA, 1997; pp. 506–507. [Google Scholar]
  35. Maddox, R.A.; Zhang, J.; Gourley, J.J.; Howard, K.W. Weather radar coverage over the contiguous United States. Weather Forecast. 2002, 17, 927–934. [Google Scholar] [CrossRef]
  36. Vivekanandan, J.; Zhang, G.; Ellis, S.M.; Rajopadhyaya, D.; Avery, S.K. Radar reflectivity calibration using differential propagation phase measurement. Radio Sci. 2003, 38, 14-1–14-14. [Google Scholar] [CrossRef] [Green Version]
  37. Gorgucci, E.; Bechini, R.; Baldini, L.; Cremonini, R.; Chandrasekar, V. The Influence of Antenna Radome on Weather Radar Calibration and Its Real-Time Assessment. J. Atmos. Ocean. Technol. 2013, 30, 676–689. [Google Scholar] [CrossRef]
  38. Frasier, S.J.; Kabeche, F.; Figueras i Ventura, J.; Al-Sakka, H.; Tabary, P.; Beck, J.; Bousquet, O. In-Place Estimation of Wet Radome Attenuation at X Band. J. Atmos. Ocean. Technol. 2013, 30, 917–928. [Google Scholar] [CrossRef]
  39. Salazar-Cerreño, J.L.; Chandrasekar, V.; Trabal, J.M.; Siquera, P.; Medina, R.; Knapp, E.; McLaughlin, D.J. A Drop Size Distribution (DSD)-Based Model for Evaluating the Performance of Wet Radomes for Dual-Polarized Radars. J. Atmos. Ocean. Technol. 2014, 31, 2409–2430. [Google Scholar] [CrossRef]
  40. Mancini, A.; Salazar, J.L.; Lebrón, R.M.; Cheong, B.L. A Novel Instrument for Real-Time Measurement of Attenuation of Weather Radar Radome Including Its Outer Surface. Part I: The Concept. J. Atmos. Ocean. Technol. 2018, 35, 953–973. [Google Scholar] [CrossRef]
  41. Mancini, A.; Salazar, J.L.; Lebrón, R.M.; Cheong, B.L. A Novel Instrument for Real-Time Measurement of Attenuation of Weather Radar Radome Including Its Outer Surface. Part II: Applications. J. Atmos. Ocean. Technol. 2018, 35, 975–991. [Google Scholar] [CrossRef]
  42. Schönhuber, M.; Urban, H.E.; Poiares Baptista, P.P.V.; Randeu, W.L.; Riedler, W. Weather radar versus 2D-video-disdrometer data. In Weather Radar Technology for Water Resources Management; Bragg, B., Jr., Massambani, O., Eds.; UNESCO Press: Paris, France, 1997; pp. 159–171. [Google Scholar]
  43. Kruger, A.; Krajewski, W.F. Two-dimensional video disdrometer: A description. J. Atmos. Ocean. Technol. 2002, 19, 602–617. [Google Scholar] [CrossRef]
  44. Thurai, M.; Bringi, V.N. Drop axis ratios from 2D video disdrometer. J. Atmos. Ocean. Technol. 2005, 22, 966–978. [Google Scholar] [CrossRef]
  45. Smith, P.L.; Kliche, D.V. The Bias in moment estimators for parameters of drop size distribution functions: Sampling from exponential distributions. J. Appl. Meteorol. 2005, 44, 1195–1205. [Google Scholar] [CrossRef]
  46. Chen, C.S.; Chen, Y.L. The rainfall characteristics of Taiwan. Mon. Weather Rev. 2003, 131, 1323–1341. [Google Scholar] [CrossRef]
  47. Barber, P.; Yeh, C. Scattering of electromagnetic waves by arbitrarily shaped dielectric bodies. Appl. Opt. 1975, 14, 2864–2872. [Google Scholar] [CrossRef] [PubMed]
  48. Brandes, E.A.; Zhang, G.; Vivekanandan, J. Experiments in rainfall estimation with a polarimetric radar in a subtropical environment. J. Appl. Meteorol. 2002, 41, 674–685. [Google Scholar] [CrossRef] [Green Version]
  49. Rodgers, C.D. Inverse Methods for Atmospheric Soundings: Theory and Practice; World Scientific Co. Pte. Ltd.: Singapore, 2000; 256p. [Google Scholar] [CrossRef]
  50. Hubbert, J.; Chandrasekar, V.; Bringi, V.N.; Meischner, P. Processing and interpretation of coherent dual-polarized radar measurements. J. Atmos. Ocean. Technol. 1993, 10, 155–164. [Google Scholar] [CrossRef] [Green Version]
  51. Hubbert, J.; Bringi, V.N. An iterative filtering technique for the analysis of copolar differential phase and dual- frequency radar measurements. J. Atmos. Ocean. Technol. 1995, 12, 643–648. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, Y.; Chandrasekar, V. Algorithm for estimation of the specific differential phase. J. Atmos. Ocean. Technol. 2009, 26, 2565–2578. [Google Scholar] [CrossRef]
  53. Giangrande, S.E.; McGraw, R.; Lei, L. An application of linear programming to polarimetric radar differential phase processing. J. Atmos. Ocean. Technol. 2013, 30, 1716–1729. [Google Scholar] [CrossRef]
  54. Cifelli, R.; Petersen, W.A.; Carey, L.D.; Rutledge, S.A.; da Silva Dias, M.A.F. Radar observations of the kinematic, microphysical, and precipitation characteristics of two MCSs in TRMM LBA. J. Geophys. Res. Atmos. 2002, 107, 8077. [Google Scholar] [CrossRef] [Green Version]
  55. Rinehart, R.E.; Garvey, E.T. Three-dimensional storm motion detection by conventional weather radar. Nature 1978, 273, 287–289. [Google Scholar] [CrossRef]
Figure 1. An example of (a) RCWF and (b) NCU C-POL lowest available data altitude (km) for quantitative precipitation estimation (QPE) after the partial beam blockage (PBB) and non-meteorological signal removal.
Figure 1. An example of (a) RCWF and (b) NCU C-POL lowest available data altitude (km) for quantitative precipitation estimation (QPE) after the partial beam blockage (PBB) and non-meteorological signal removal.
Remotesensing 13 00154 g001
Figure 2. The flow chart of the quality control (QC) and QPE procedures.
Figure 2. The flow chart of the quality control (QC) and QPE procedures.
Remotesensing 13 00154 g002
Figure 3. The distribution of rainfall intensities calculated from 6-min raindrop size distributions (RSDs) for different precipitation types (spring rain, Mei-Yu, convection, typhoon, and cold front).
Figure 3. The distribution of rainfall intensities calculated from 6-min raindrop size distributions (RSDs) for different precipitation types (spring rain, Mei-Yu, convection, typhoon, and cold front).
Remotesensing 13 00154 g003
Figure 4. The normalized bias (NBIAS) of estimated rainfall rates using seasonal (solid lines) and all-season (dashed lines) coefficients at (a) C-band and (b) S-band. (c) and (d) are the same but for the normalized root mean square error (NRMSE). The estimated rainfall rates are derived from the RSD-simulated radar variables and compared to RSD-derived rainfall rates for the theoretical verification.
Figure 4. The normalized bias (NBIAS) of estimated rainfall rates using seasonal (solid lines) and all-season (dashed lines) coefficients at (a) C-band and (b) S-band. (c) and (d) are the same but for the normalized root mean square error (NRMSE). The estimated rainfall rates are derived from the RSD-simulated radar variables and compared to RSD-derived rainfall rates for the theoretical verification.
Remotesensing 13 00154 g004
Figure 5. An example of time steps of available data for QPE on 14 June, 2015. Red circles represent the NCU C-POL discrete QPE (NCU C-POL D-QPE), and blue crosses for the RCWF discrete QPE (RCWF D-QPE). Black triangles are the discretely combined C- and S-band QPE (DC-QPE) and green dots are the Lagrangian-evolution adjustment QPE (LEA-QPE).
Figure 5. An example of time steps of available data for QPE on 14 June, 2015. Red circles represent the NCU C-POL discrete QPE (NCU C-POL D-QPE), and blue crosses for the RCWF discrete QPE (RCWF D-QPE). Black triangles are the discretely combined C- and S-band QPE (DC-QPE) and green dots are the Lagrangian-evolution adjustment QPE (LEA-QPE).
Remotesensing 13 00154 g005
Figure 6. The NRMSEs of the C-band (refer to NCU C-POL) and S-band (refer to RCWF hereafter) QPEs at different levels of ZHH corrections. The height of each bar indicates the NRMSE after systematic-bias correction only. The brown color represents the reductions of NRMSE after the ensuing attenuation correction, and the yellow marks the reductions after the correction of the wet radome effect (WRE). The blue color is the final NRMSE when all QC procedures, including the systematic-bias, attenuation and wet-radome corrections, are completed.
Figure 6. The NRMSEs of the C-band (refer to NCU C-POL) and S-band (refer to RCWF hereafter) QPEs at different levels of ZHH corrections. The height of each bar indicates the NRMSE after systematic-bias correction only. The brown color represents the reductions of NRMSE after the ensuing attenuation correction, and the yellow marks the reductions after the correction of the wet radome effect (WRE). The blue color is the final NRMSE when all QC procedures, including the systematic-bias, attenuation and wet-radome corrections, are completed.
Remotesensing 13 00154 g006
Figure 7. Same as Figure 4 but from real cases. Eighteen heavy precipitation events observed by NCU C-POL and RCWF are evaluated via rain gauge measurements.
Figure 7. Same as Figure 4 but from real cases. Eighteen heavy precipitation events observed by NCU C-POL and RCWF are evaluated via rain gauge measurements.
Remotesensing 13 00154 g007
Figure 8. The NRMSE of D-QPEs as a function of rainfall rates for (a) C-band and (b) S-band D-QPEs.
Figure 8. The NRMSE of D-QPEs as a function of rainfall rates for (a) C-band and (b) S-band D-QPEs.
Remotesensing 13 00154 g008
Figure 9. The NRMSE (as a function of rainfall rates) of C-/S-band D-QPEs from the (a) R(ZHH), (b) R (ZHH, ZDR) (c) R(KDP), and (d) R (KDP, ZDR) algorithm.
Figure 9. The NRMSE (as a function of rainfall rates) of C-/S-band D-QPEs from the (a) R(ZHH), (b) R (ZHH, ZDR) (c) R(KDP), and (d) R (KDP, ZDR) algorithm.
Remotesensing 13 00154 g009
Figure 10. The NRMSEs of the C-band D-QPE, S-band D-QPE and DC-QPE from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm.
Figure 10. The NRMSEs of the C-band D-QPE, S-band D-QPE and DC-QPE from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm.
Remotesensing 13 00154 g010
Figure 11. The C-band QPE maps derived from R(KDP, ZDR) at (a) 06:54 UTC ( R T i ) and (b) 07:04 UTC ( R T i + 1 ) on 19 August in 2014, and LEA-QPE at (c) 06:59 UTC as Δ t j = 5 min. The solid contour lines represent the precipitation area (R > 10 mm h−1) and ⨁ marks indicate the convective core where R > 100 mm h−1.
Figure 11. The C-band QPE maps derived from R(KDP, ZDR) at (a) 06:54 UTC ( R T i ) and (b) 07:04 UTC ( R T i + 1 ) on 19 August in 2014, and LEA-QPE at (c) 06:59 UTC as Δ t j = 5 min. The solid contour lines represent the precipitation area (R > 10 mm h−1) and ⨁ marks indicate the convective core where R > 100 mm h−1.
Remotesensing 13 00154 g011
Figure 12. The NRMSE of C-/S-band D-QPE and LEA-QPE from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm.
Figure 12. The NRMSE of C-/S-band D-QPE and LEA-QPE from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm.
Remotesensing 13 00154 g012
Figure 13. The NBIAS of C-/S-band LEA-QPE, DC-QPE and LEA DC-QPE (LEAC-QPE) derived from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm. (c) and (d) are the same but showing NRMSE.
Figure 13. The NBIAS of C-/S-band LEA-QPE, DC-QPE and LEA DC-QPE (LEAC-QPE) derived from the (a) R(KDP) and (b) R(KDP, ZDR) algorithm. (c) and (d) are the same but showing NRMSE.
Remotesensing 13 00154 g013
Table 1. The basic information for RCWF (radar code of Wu-Fanshan, NEXRAD WSR-88D S-band dual-pol radar, Ryzhkov et al. [9] from the Central Weather Bureau (CWB) and NCU C-POL) National Central University C-band polarimetric radar).
Table 1. The basic information for RCWF (radar code of Wu-Fanshan, NEXRAD WSR-88D S-band dual-pol radar, Ryzhkov et al. [9] from the Central Weather Bureau (CWB) and NCU C-POL) National Central University C-band polarimetric radar).
RCWFNCU C-POL
Wavelength (cm)10.75.3
Location (degree)121.77 N
25.07 E
212.18 N
24.97 E
Height (m)766156
Range resolution (m)250250
Beam resolution (degree)0.51
Table 2. The coefficients of four power-law rain rate relations for S-band radar. The all-season coefficients are derived from the entire database measured by a 2D-Video Disdrometer (2DVD) at National Central University (NCU). The seasonal coefficients are classified into spring, Mei-Yu, summer convection, typhoon, and winter cold front.
Table 2. The coefficients of four power-law rain rate relations for S-band radar. The all-season coefficients are derived from the entire database measured by a 2D-Video Disdrometer (2DVD) at National Central University (NCU). The seasonal coefficients are classified into spring, Mei-Yu, summer convection, typhoon, and winter cold front.
All-SeasonSpringMei-YuSummer ConvectionTyphoonWinter Cold Front
R(ZHH)a10.02790.01970.02440.04350.02820.0408
b10.66190.68740.67790.62330.66240.6173
R(ZHH,ZDR)a20.00460.00190.00180.00110.00130.0033
b20.84920.94520.95781.00170.94900.8888
c2−0.6193−0.9734−1.0434−1.1240−0.7988−0.7439
R(KDP)a347.599844.686448.051648.344864.329342.5163
b30.76050.79500.79150.77250.72780.7225
R(KDP,ZDR)a464.841161.942162.363363.363373.096460.2012
b40.98800.97820.97270.97270.94760.9486
c4−0.6921−0.6445−0.6196−0.6196−0.6039−0.5836
Table 3. Same as Table 2 but for C-band radar.
Table 3. Same as Table 2 but for C-band radar.
All-SeasonSpringMei-YuSummer ConvectionTyphoonWinter Cold Front
R(ZHH)a10.03760.02600.03160.07100.03600.0434
b10.63400.63300.65580.57610.63940.6138
R(ZHH,ZDR)a20.00350.00140.00140.00130.00100.0028
b20.88860.99220.99521.00180.98120.9199
c2−0.6575−0.9840−1.0031−1.0239−0.7714−0.7474
R(KDP)a326.234323.948025.861926.488436.167024.0925
b30.74850.78230.77840.75900.71580.7103
R(KDP,ZDR)a431.251429.845930.410629.974736.896530.3301
b40.96480.95630.95930.93810.92120.9500
c4−0.5988−0.5334−0.5418−0.5132−0.5146−0.5717
Table 4. The information of 18 heavy rainfall events investigated in this study.
Table 4. The information of 18 heavy rainfall events investigated in this study.
Date (hour)Precipitation TypeMax. Hourly Rainfall (mm h−1)/Accumulated Rainfall (mm)Scan Time (min:s)
RCWFNCU C-POL
26 April 2014 (1 h)Spring605:506:40
5 May 2014 (2 h)Mei-Yu20/335:506:40
9 May 2014 (2 h)Mei-Yu23.5/355:506:40
15 May 2014 (1 h)Mei-Yu57.55:506:40
20 May 2014 (3 h)Mei-Yu26/505:506:40
21 May 2014 (2 h)Mei-Yu38/52.55:506:40
29 May 2014 (6 h)Mei-Yu56.5/1155:506:40
5 June 2014 (6 h)Mei-Yu46/135.55:506:40
7 June 2014 (1 h)Mei-Yu355:506:40
31 July 2014 (1 h)Summer convection48.55:5010:00
13 August 2014 (1 h)Summer convection36.55:5010:00
19 August 2014 (3 h)Summer convection71.5/92.55:5010:00
24 February 2015 (4 h)Cold front29.5/58.55:5010:00
12 May 2015 (2 h)Mei-Yu41.5/50.55:5010:00
22 May 2015 (2 h)Mei-Yu19.5/285:5010:00
6 June 2015 (2 h)Mei-Yu37/45.55:5010:00
14 June 2015 (2 h)Mei-Yu103.5/1675:5010:00/12:30
23 July 2015 (2 h)Summer convection74/805:5010:00/12:30
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, J.-Y.; Chang, W.-Y.; Chang, P.-L. A Synthetic Quantitative Precipitation Estimation by Integrating S- and C-Band Dual-Polarization Radars over Northern Taiwan. Remote Sens. 2021, 13, 154. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13010154

AMA Style

Chen J-Y, Chang W-Y, Chang P-L. A Synthetic Quantitative Precipitation Estimation by Integrating S- and C-Band Dual-Polarization Radars over Northern Taiwan. Remote Sensing. 2021; 13(1):154. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13010154

Chicago/Turabian Style

Chen, Ju-Yu, Wei-Yu Chang, and Pao-Liang Chang. 2021. "A Synthetic Quantitative Precipitation Estimation by Integrating S- and C-Band Dual-Polarization Radars over Northern Taiwan" Remote Sensing 13, no. 1: 154. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13010154

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