Next Article in Journal
Process-Driven and Flow-Based Processing of Industrial Sensor Data
Next Article in Special Issue
NVIS Multicarrier Modulations for Remote-Sensor Applications
Previous Article in Journal
A Paper-Based Near-Infrared Optical Biosensor for Quantitative Detection of Protease Activity Using Peptide-Encapsulated SWCNTs
Previous Article in Special Issue
ANN-Based Estimation of Low-Latitude Monthly Ocean Latent Heat Flux by Ensemble Satellite and Reanalysis Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving Radar Rainfall Estimations with Scaled Raindrop Size Spectra in Mei-Yu Frontal Rainstorms

1
College of Meteorology and Oceanography, National University of Defense Technology, Nanjing 211101, China
2
Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(18), 5257; https://doi.org/10.3390/s20185257
Submission received: 12 August 2020 / Revised: 2 September 2020 / Accepted: 10 September 2020 / Published: 14 September 2020
(This article belongs to the Special Issue Remote Sensor Based Geoscience Applications)

Abstract

:
Hydrological calibration of raw weather radar rainfall estimation relies on in situ rainfall measurements. Raindrop size distribution (DSD) was collected during three typical Mei-Yu rainstorms in July 2014 using three particle size velocity (Parsivel) DSD sensors along the Mei-Yu front in Nanjing, Chuzhou, and the western Pacific, respectively. To improve the radar precipitation estimation in different parts of the Mei-Yu front, a scaling method was adopted to formulate the DSD model and further derive the ZR relations. The results suggest a distinct variation of DSDs in different parts of the Mei-Yu front. Compared with statistical radar Z–ARb relations obtained by mathematical fitting techniques, the use of a DSD model fitting based on a scaling law formulation theoretically shows a significant improvement in both stratiform (33.9%) and convective (2.8%) rainfall estimations of the Mei-Yu frontal system, which indicates that using a scaling law can better reflect the DSD variations in different parts of the Mei-Yu front. Polarimetric radar has indisputable advantages with multiparameter detection ability. Several dual-polarization radar estimators are also established by DSD sensor data, and the R(ZH, ZDR) estimator is proven to be more accurate than traditional Z–R relations in Mei-Yu frontal rainfall, with potential applications for operational C-band polarimetric radar.

1. Introduction

Changes in the spatial and temporal patterns of climate variables associated with global warming will have an influence on regional- and catchment-scale hydrological processes [1]. According to the Global Climate Observing System (GCOS), precipitation is getting more severe with drastic changes in most of the cities worldwide [2]. For instance, Hoerling et al. [3] found an increase in both frequency and intensity of heavy rain since 1979 in the northeastern United States. Cui et al. [4] revealed the interannual variability of heavy rain in central and southern China due to the variation of large-scale environmental conditions. Investigating the response of regional rainfall events, especially extreme precipitation, will have significant implications on climate prediction. Over the last few decades, weather radars have been principally used to collect rainfall variability because of their good areal coverage as well as high-resolution measurements both in time and space [5].
Various radar quantitative precipitation estimation (QPE) algorithms, including ZR relations (Z = ARb, where Z is the radar reflectivity factor and R is the rain rate), as well as polarimetric radar algorithms highly depend on surface-based DSD measurements [6,7,8,9,10,11]. Both Z and R are functions of raindrop size distributions, which can be directly derived from given DSD samples. Earlier studies tend to obtain the parameters [A, b] in ZR relations via different mathematical fitting techniques, but the statistical results sometimes cause severe deviations [12,13]. Other scholars establish Z–R relations by relating coefficients [A, b] to DSD fitting models based on a scaling law formulation [14,15]. Unlike other DSD models (e.g., exponential distributions [16] and gamma distributions [17]), the scaling law allows the ZR relation to be established without any DSD shapes imposed a priori [18]. More importantly, it provides valuable information about the intrinsic microphysical properties of the drop size distribution and its relations with radar parameters [A, b], which is likely to help improve radar rainfall estimations.
By using the scaling theory, Hazenberg et al. [19] studied the DSD characteristics and radar retrieval variability in mesoscale convective systems (MCSs), squall lines, and midlatitude cyclones over the Mediterranean region. However, no research has specifically examined the use of scaling law formulas in the Mei-Yu frontal system, a quasi-stationary front stretching over thousands of kilometers, extending eastward from southwestern China to the western Pacific through the Japanese archipelago [20,21], which is one of the most significant rainfall systems affecting the hydrological cycle in the East Asian monsoon region. Rainstorms associated with Mei-Yu (also called “Baiu” in Japan and “Changma” in South Korea) are a crucial summertime water resource and sometimes result in severe natural disasters in East Asian countries [22,23,24,25,26,27]. Meanwhile, this frontal rainstorm also intimidates the security of naval ships on the sea. Therefore, there is an urgent need to improve the performance of radar rainfall estimations within this weather system for flood monitoring, urban waterlogging prediction, and ocean storm monitoring. Based on the scaling method, one specific goal herein is to improve radar rainfall retrievals in different parts of the Mei-Yu frontal system by using in situ DSD sensor measurement data collected during three Mei-Yu rainstorms in the summer of 2014.
Following this introduction, a brief description of data and methods is presented in Section 2, the DSD properties among different parts of the Mei-Yu front are analyzed in Section 3, and radar rainfall retrievals based on DSD sensor measurements are further discussed in Section 4. Finally, Section 5 presents a summary and conclusions.

2. Data and Methods

2.1. Observational Sites and Instruments

The data collected for the analyses consisted of three-time series of 1-min DSD data that were measured by Parsivel2 DSD sensors at Nanjing (NJ, 118.5° E, 32.0° N), Chuzhou (CZ, 118.2° E, 32.3° N), and the western Pacific (WP, 148.0° E, 32.1° N). DSD data in WP were measured by an onboard Parsivel2 DSD sensor during a marine survey over the western Pacific. Three typical Mei-Yu front rain events were simultaneously detected by the above DSD sensors on the 2nd, 4th, and 12th of July, 2014, respectively when the ship sailed across the area shown in Figure 1. Meanwhile, the satellite image from an infrared-channel sensor on the Japanese geostationary meteorological satellite Himawari-7 clearly presents a large-scale cloud belt covering southwestern China all the way to the central Pacific (Figure 2a). The rainfall is caused by a quasi-stationary front and the presence of strong low pressure as presented on the surface weather chart (Figure 2b) provided by the Japan Meteorological Agency (JMA).
The Parsivel2 DSD sensor [28], manufactured by OTT Hydromet, Germany, is a 2nd generation optical sensor with a near-infrared (650 nm) rectangular single-beam that is 30 mm wide, 180 mm long (for a sensor horizontal area of 54 cm2), and 10 mm high. It archives equivalent drop diameters sorted unevenly into 32 diameter classes from 0 to 26 mm and 32 fall speed classes from 0 to 22.4 m s−1. The time resolution can be selected and was set to 1 min in the present study. Nevertheless, some measurement errors, such as strong winds, margin fallers, and splashing effects are found to influence the reliable utilization of DSD sensors. Following Wu et al. [29], we implemented a procedure on data quality-control by eliminating particles with diameters larger than 8 mm or falling speeds 50% faster or slower than the empirical raindrop speed-diameter relationship proposed by Gunn and Kinzer [30]. Apart from that, 1-min samples with drop numbers less than 10 or a rain rate less than 0.1 mm h−1 were also excluded [31]. In order to overcome the measurement error from ship movements and tilts, the onboard Parsivel2 sensor was adjusted and fixed to stay perpendicular to the direction of ship pitch to best ensure the beam orientation parallel to the horizontal plane. Finally, there are overall 2611 1-min effective DSD samples for NJ, 2306 1-min effective DSD samples for CZ, and 653 samples for WP (see Table 1). The samples in NJ (CZ) are plenitudinous with the max rain rate reaching 88.22 (71.90) mm h−1 on 12 July (4th). Although the onboard DSD sensor observations in WP are over a relatively short period of time, the max rain rate obtained in WP could reach 48.64 mm h−1 on 12 July, and the max reflectivity factor is up to 48.22 dBZ on 2 July.

2.2. Scaling of Raindrop Size Distribution

The raindrop size distribution was computed from the Parsivel2 DSD sensor counts as below:
N ( D i ) = 1 S e f f ( D i ) T Δ D i j = 1 L n i j V j
where N(Di) (mm−1 m−3) is the number concentration of raindrops per unit volume per unit size interval for raindrop diameter Di (mm); L is the bin number of Parsivel2 measurements and is known to be 32; nij is the number of raindrops within size bin i and velocity bin j; T is the sampling time (set to 60 s in this research), and Vj (m s−1) is the fall-speed for velocity bin j. Seff(D) (mm2) is the effective sampling area [32] that considers partially detected drops across Parsivel’s laser sheet and is computed based on beam width and length:
S e f f ( D i ) = 180 × ( 30 D i 2 )
The integral rainfall parameters, including the radar reflectivity factor Z (mm6 m−3), rain rate R (mm h−1), and total concentration of raindrops Nt (mm−3) are directly derived from measured DSDs and fall velocity as below:
Z = i = 1 L N ( D i ) D i 6 Δ D i
R = 6 π × 10 4 i = 1 L j = 1 L V j D i 3 N ( D i ) Δ D i
N t = i = 1 L N ( D i ) Δ D i
To better characterize the measured DSD samples, Sempere Torres et al. [14,15] proposed a more general formulation for DSD parameterization, that is, to express the DSD as a scaling law (SL)
N ( D , Ψ ) = Ψ α g ( D Ψ β )
where Ψ is a reference variable, commonly taken as the rain rate R due to its considerable variation in different rainfall systems (e.g., for deep convective system, the Rmax could reach up to 100 mm h−1), and this contributes to comparing the DSD properties over a wide range. Thus, Equation (6) can also be expressed as
N ( D , R ) = R α g ( D R β )
where α and β are scaling parameters and g(x) is the scaled raindrop size distribution with x = D/Rβ. Introducing Equation (7) into (4) yields the so-called self-consistency-constraint [19]:
6 π × 10 4 c 0 x 3 + d g ( x ) d x = 1
where the terminal fall-speed is assumed to follow the power-law relationship of Atlas and Ulbrich [33]
V ( D ) = c D d
with c = 3.778 and d = 0.67. Different functional shapes have been proposed for g(x). Herein both the exponential (µ = 0) and gamma distribution are employed, which can be expressed as
g ( x ) = κ x μ exp ( Λ x )
then ĸ can be obtained by introducing Equation (10) into (8)
κ = Λ 4.0 + d + μ 6 π × 10 4 c Γ ( 4.0 + d + μ )

2.3. Establishment of Z–R Relations

The ZR relationship can be established via the above SL equations as well as Z = ARb, the coefficients A and b thereby can be calculated by
A = κ Γ ( 7.0 + μ ) Λ 7.0 + μ , b = α + 7 β
where b is related to α and β, and A is dependent on g(x). Note that the SL method can naturally lead to power-law relationships between Z and R, without any functional shape imposed as a priori for g(x) [18].
The SL parameters α and β, as well as the parameters [µ, Λ] of g(x) for a given functional shape can be obtained by using the moment estimation method [17,34,35]. For a given DSD, the nth-order moment (Mn) is defined as
M n = 0 N ( D ) D n d D = 0 D n R α g ( D / R β ) d D = θ n R γ n
θ n = 0 x n g ( x ) d x , γ n = α + ( n + 1 ) β
Hence, to obtain the SL parameters α and β, γn should be first obtained by fitting a power law between different integral moments Mn and the rain rate R Equation (13), then α and β can be obtained according to Equation (14) that γn should follow a linear relationship with the moment order (n + 1).
In addition, the parameters [µ, Λ] of g(x) can be calculated with the second, fourth, and sixth truncated moments (T246) [35]. And the mass-weighted mean diameter Dm (mm) is calculated by the ratio of the fourth to the third moment of the DSD
η = M 4 2 M 2 M 6
μ = ( 7 11 η ) [ ( 7 11 η ) 2 4 ( η 1 ) ( 30 η 12 ) ] 1 / 2 2 ( η 1 )
Λ = [ ( 4 + μ ) ( 3 + μ ) M 2 M 4 ] 1 / 2
D m = i = 1 32 N ( D i ) D i 4 Δ D i i = 1 32 N ( D i ) D i 3 Δ D i = M 4 M 3

2.4. Variables of Polarimetric Radar

In addition to conventional weather radar measurements, a polarimetric radar is capable of measuring the differential reflectivity ZDR (dB) and the specific differential phase Kdp (deg km−1) between two orthogonally polarized radar returns, which helps to significantly improve the radar data quality by reducing the impact of DSD variability on the accuracy of rainfall estimations [36]. These polarimetric variables are computed as below:
Z h , v = 4 λ 4 π 4 | K w | 2 D min D max | f h h , v v ( D ) | 2 N ( D ) d D
Z D R = 10 log 10 ( Z h Z v )
K d p = 10 3 180 π λ Re { D min D max [ f h ( D ) f v ( D ) ] N ( D ) d D }
where λ is the radar wavelength; Kw is the dielectric constant factor of water and fhh,vv(D) (fh,v(D)) is the backscattering (forward-scattering) amplitude at the horizontal or vertical polarizations. In this research, the assumed C-band (50 mm) dual-polarization radar variables are derived from the Parsivel2 DSD sensor data by the use of T-matrix scattering approach for non-spherical particles [37]. The mean and standard deviation of the canting angle is set to 0° and 7° respectively with a Gaussian model. The raindrop temperature is assumed to be 20 °C for summer season and the relation of axis ratio r and its equivalent diameter D (mm) based on Brandes et al. [38] is used as below:
r = 0.9951 + 0.02510 D 0.03644 D 2 + 0.005303 D 3 0.0002492 D 4

2.5. Assessment Statistics

To evaluate the overall quality of radar QPE algorithm and eliminate the influence from different rainfall totals over three sites in the meantime, the normalized absolute error (NAE), the normalized mean bias (NB), and the coefficient of determination (r2) are examined in each site as defined below:
NAE = j = 1 N | R radar R j | j = 1 N R j
NB = j = 1 N ( R radar R j ) j = 1 N R j
r 2 = 1 j = 1 N ( R radar R j ) 2 j = 1 N ( R j R j ¯ ) 2
where N is the number of effective samples in rain event, Rj and Rradar represent the observed rain rate by DSD sensors and the radar-retrieved rain rate for each sample, respectively. Note that the specific expression of Rradar is not given here, since the radar-retrieved results are from either conventional radar or polarimetric radar (e.g., Rradar = (Zj/A)1/b for conventional radar).

2.6. Classification Scheme of Rain Types

Based on DSD sensor measurements, the DSD samples were further classified into two fundamental rain types (convective rain and stratiform rain). Several rainfall classification schemes have been well developed [34,39,40]. Nevertheless, the results of Tokay and Short [34] were obtained from rainfall clusters in the tropics, which is inappropriate in our Mei-Yu study because of the regional variability. The categorization scheme of Testud et al. [39] was based on the variability of rain rate (R) with time. For a time series of the rain rate {Ri}, the sample k was classified as stratiform only if the R values of ten adjacent values from Rk−5 to Rk + 5 were all less than 10 mm h−1, otherwise, it was assumed to be convective. And the classification scheme of Bringi et al. [40] was based on the standard deviation value (σR) of the rain rate R. A threshold of σR = 1.5 mm h−1 was used for the classification of stratiform and convective rain. In this research, two classification schemes were adopted together to separate total samples into convective and stratiform clusters. More specifically, for ten consecutive 1-min samples, if the R values of ten adjacent values were all less than 10 mm h−1 and the standard deviation σR was less than 1.5 mm h−1, then the sample was defined as stratiform, otherwise, it was classified as convective. As a result, NJ, CZ, and WP consisted of 77% (23%), 90% (10%), and 41% (59%) of stratiform (convective) rainfall DSD samples, respectively. The results indicate the presence of strong convective rain over the Pacific Ocean.

3. DSD Properties from Parsivel2

3.1. Analyses of Drop Size Spectra

In Figure 3 and Figure 4 the analyses of the raindrop size spectra are presented. Figure 3a presents the values of the exponents γn calculated by the least squares method between different integral moments Mn and the rain rate R (Equation (13), for 0 ≤ n ≤ 6). The fit was performed over the central moments (1 ≤ n ≤ 5) in order to minimize instrumental effects that might occur in measuring the smaller and larger raindrops by Parsivel DSD sensors. Except for the lowest-order moments, the overall fit in panels (a) and (b) of Figure 3 was very good. Based on the estimated values of γn, the coefficients α and β can be further obtained via Equation (14), and the results are also shown in each panel of Figure 3. For stratiform situations, it can be found that the sum of two coefficients was positive in CZ, negative in WP, and equal to zero in NJ, which suggests an uneven distribution of rainwater in different parts of the Mei-Yu front. Note that the coefficient α had a negative value under all circumstances in the Mei-Yu front. This was different from the results of squall line (α = 0.24) and MCS (α = 0.22) as reported by Hazenberg et al. [19], which suggests that the DSD characteristics in the Mei-Yu front were unique and worth our further research.
To develop deeper insight of the Mei-Yu microphysics, the scaled raindrop spectra along with two different functional shapes are presented in Figure 4. It can be seen that the normalized DSD data exhibited an exponential shape, which was well fitted by both exponential and gamma distributions, as shown by the high values of the coefficient of determination r2 calculated by Equation (25) and used for measuring the fitting efficiency of the relations. Comparing two rain-types, a wider normalized spectrum can be noticed in stratiform rain clusters compared to convective rain clusters, which indicates completely different DSD characteristics of the Mei-Yu front unlike other strong convective systems such as squall lines. Comparing the three regions from different parts of the Mei-Yu front, the max particle concentration of convective rain was largest in WP and smallest in CZ, and the spectra of stratiform rain were widest in NJ and narrowest in CZ. This reflects the distinct DSD variability in different parts of the Mei-Yu front. The SL method was able to provide information about the intrinsic microphysical properties of DSD. It was reported by Uijlenhoet et al. [41,42] that precipitation may reach an equilibrium state, where the collision-coalescence and collision-breakup processes of raindrops are nearly in balance, causing all characteristic raindrop sizes to remain constant when the obtained SL coefficients (α, β) satisfy α = 1, β = 0. From Figure 3d one can notice that CZ had the smallest β value (0.25) for stratiform rain as compared to NJ (0.27) and WP (0.30), which corresponds to its more uniform DSD and narrower width of spectra in Figure 4d. A different rainfall control mechanism, on the other hand, can be found in the stratiform rain of NJ. Note that the sum of coefficients α and β was equal to zero in NJ as presented in Figure 3b, which results in non-homogeneous precipitation [43] and distinct variations of drop size. This is also called a “size-control” mechanism.

3.2. Statistics of DSD Parameters

To further study the rainfall microphysics in different parts of Mei-Yu front, Figure 5 shows the evolution of Dm and log10(Nt) versus rain rate for DSD samples observed in three parts. The parameters Dm (blue rectangle) and log10(Nt) (red circle) are presented in six classes of rain rate with an interval of 10 mm h−1 and the ranges of both parameters are indicated by the whisker of each symbol in Figure 5. Additionally given in Table 2 is the number of samples in each rain rate class. Remarkably, both Dm and log10(Nt) increased with an increase in rain rate at the first three classes (R < 30 mm h−1). However, at the last three classes where R > 30 mm h−1, the parameters seemed to become stable with Dm values staying approximately at 2.16 (2.12) mm in NJ (CZ) and 2.08 mm in WP, while the log10(Nt) values still increased slowly with an increasing rain rate. Note that the ranges of both Dm and log10(Nt) become smaller and smaller with rain rate increases. At higher rain rate classes, the DSDs may reach an equilibrium state where coalescence and breakup processes of raindrops are in close balance [44]. Under equilibrium conditions, Dm remained nearly constant, while any increase in rain rate was mainly attributed to an increase in the concentration of raindrops [8], which was also called a “number-control” mechanism. The unique characteristics of Dm and log10(Nt) in our study were consistent with this theory but further demonstrate that the rainfall in the Mei-Yu front is unevenly distributed with larger raindrops and more abundant particles in its continental part.

3.3. µ–Λ Relations

It has been widely reported that µ and Λ show rather similar behaviors and a very strong correlation [35,45]. Studies have also shown that radar rainfall estimations can be improved after adjusting the µΛ relationship to ground DSD sensor observations [46]. Over the past few years, it was found that µΛ relations vary with respect to rain types, climate regimes, and terrains [35,46,47,48]. Hence, the µΛ relation needs to be adapted to the Mei-Yu frontal system.
Following the method of data-processing in Zhang et al. [35], the DSDs with R > 5 mm h−1 and N > 1000 were fitted via the moment method to obtain a second-degree polynomial µΛ relation in three parts of the Mei-Yu front. The relation for the NJ is given by
Λ = 0.0156μ2 + 0.636μ + 1.533
For CZ the relation is given by:
Λ = 0.0072μ2 + 0.688μ + 1.426
For WP the relation is given by:
Λ = 0.0244μ2 + 0.608μ + 1.351
For the entire Mei-Yu front is given by:
Λ = 0.0080μ2 + 0.741μ + 1.432
The fitting lines and their corresponding scatter plots are presented in Figure 6. Given the same Λ, it’s notable that the parameter µ of the WP is less than that of NJ (CZ). The μΛ relation can also be formulated as ΛDm = 4 + μ [35,45]. As presented in Figure 6, more samples of precipitation in NJ (gray circle) and CZ (gray cross) lie in the larger Dm region. Such characterization is consistent with the results presented in Figure 5.
Compared with the statistical results in other regions, presented together in Figure 6, the results that we derived in three parts of the Mei-Yu front were right between those obtained by Zhang et al. [35] and Chen et al. [45]. This could be attributed in part to DSD sensor uncertainty. The Parsivel2 sensors used in our study generally underestimate small drops [31,37,49,50], resulting in a larger µ value than 2DVD (two-dimensional video disdrometers) observations in Zhang et al. [35]. Using the same type of Parsivel instruments, the μΛ relation obtained by Chen et al. [45] overestimated the large drops and underestimated the small drops. This could be the reason that only a small area of Mei-Yu is considered in Chen et al. [45]. Thus, it is crucial to obtain the unique µΛ relationship for different parts of the Mei-Yu front.

4. QPE of Radar

4.1. Z–R Relations

It is reported that a single, unique Z–R relation does not exist, and it strongly depends on the geographical location and physical conditions of rainfall [51]. For instance, the National Weather Service’s (NWS) Weather Surveillance radar-1988 Doppler (WSR-88D) precipitation processing subsystem recommends Z = 250R1.2 [6] for tropical systems, and the Next Generation Weather radar (NEXRAD) utilizes Z = 3 00R1.40 for mid-latitudes [52]. It was also reported by Rosenfeld and Ulbrich [53] that there are significant differences in the Z–R relations between maritime precipitation and continental precipitation. However, owing to the particularity of geographical location, in situ observations of maritime precipitation are quite difficult. The Mei-Yu front extends from eastern China to WP. There are few comparisons to demonstrate the variability of Z–R relations between its continental and maritime parts. Therefore, combined with unique DSD observations, the study of Z-–R relations of precipitation in different parts of the Mei-Yu front is required to improve the quality of radar precipitation estimations in a specific region.
To compare the scaling law method with the other fitting methods, Figure 7 presents the results of Z–R fittings using the least squares method, which is implemented by deriving the coefficient b from the fit and set the coefficient A to force the total rainfall to agree (i.e., A = (ΣZ1/bR)b), because the fitting provides greater uncertainty in estimating A than it does b [54]. Additionally presented in this plot are scatter plots of rain rate R versus radar reflectivity factor Z as measured by the DSD sensor. In addition, Z = 368R1.21 for convective precipitation observed with Parsivel sensor during Mei-Yu [45] and Z = 193.73R1.54 for stratiform rain observed with 2DVD sensor during Mei-Yu [55], are also provided in the corresponding colors for comparisons with previous research in terms of different rain types. According to the results, the Z–R relations of precipitation within different parts of the Mei-Yu front were significantly different. Comparing the precipitation in NJ (CZ) and WP, coefficient A and exponent b in Z–ARb relations had an inverse relationship. A larger b value can be found for total precipitation in WP, which can be explained via Equation (12) since WP has a larger β (see Figure 3e,f).
Compared with Z = 368R1.21 for convective rain observed by Chen et al. [45], the relation overestimated the convective precipitation in three parts when rain rate was below 10 mm h−1, whereas it underestimated the convective precipitation when the rain rate was above 10 mm h−1. Moreover, the relation was quite similar to our Z = 489R1.17 relation for stratiform precipitation in WP, which further reflects the distinct microphysical variability within different parts of Mei-Yu precipitation. Compared with Z = 193.73R1.54 for stratiform rain observed by Wen et al. [55], the relation was very close but slightly lower to our Z = 308R1.50 (Z = 275R1.52) relationship for stratiform rain in NJ (CZ), which might be related to the different method of rain categorization used and the recognition of shallow rain type in Wen et al. [55]. Due to the lack of radar observations, shallow rainfall was recognized as stratiform rain in our surface DSD sensor-based categorization schemes. Though the influence from rain type categorization herein was not very obvious, it is worth further improvement, and we leave it for future research. Based on the above analyses, as well as previous studies [56,57,58], it is suggested that DSD variability related to the precipitation microphysics in different types of rain and terrain are major sources of the diversity of Z–R relations.
Besides the standard Z = 300R1.4 relation (STD) and the relations using the least-squares fitting method (LS), the optimal ZR relations for the normalized exponential (EXP) and gamma (GAM) distributions in each site were also obtained via Equation (12). Their values are shown in Table 3. The fitting results were similar for the EXP, GAM, and LS methods and varied mainly with respect to the value of coefficient A. Compared to the LS method, both coefficients A and b were larger in EXP and GAM methods. The result of GAM, however, was closer to that of the LS method in comparison to the EXP method.
For the assessment of radar rainfall retrievals, the statistics calculated via Equations (23) and (24) are shown in Table 4. Comparing the four estimation methods at the same site, the NAE in Table 4 for the GAM method was smallest under both stratiform and convective situations, while the NAE for the STD method was rather poor as compared to the other three techniques. In comparison with statistical LS estimations, we calculated the difference of NAE errors between LS and GAM method in three sites for both stratiform rain (8.7% for NJ, 9.9% for CZ, 15.3% for WP) and convective rain (1.2% for NJ, 1.2% for CZ, 0.4% for WP), and found that the scaled gamma spectra method showed an impressive improvement in both stratiform (33.9%) and convective (2.8%) rainfall estimations of total Mei-Yu frontal rainfall. This further indicated that using scaling law fitting can well reflect the variation of DSD in different parts of the Mei-Yu front. The EXP method had a better performance in stratiform rain than convective rain, which indicated that the EXP method was very sensitive to rain rates. For scaled DSD samples, the GAM method yielded smaller NAE and NB for both rain types in comparison to the EXP method, which could be due to its more representative description of raindrop size spectra. On the other hand, the NB in Table 4 for EXP and GAM methods were both negative, while the opposite results were noticed for LS and STD methods. This indicated that the normalized spectra shapes tended to underestimate the Mei-Yu precipitation, whereas the LS method and STD relation usually overestimated the Mei-Yu rainfall.

4.2. Polarimetric Radar Applications

One of the sources of uncertainty in radar precipitation estimates, including Z–R relations, resulted from variable raindrop size distributions. It is expected that uncertainties in radar QPE due to variability in precipitation microphysics and target identification can be reduced by the use of dual-polarization techniques, which transmit both vertically and horizontally polarized electromagnetic waves. Because larger raindrops are not spherical, their backscattering cross-sections are different for the two polarizations; hence, the returned power and the Doppler shift within the two channels vary. As a result, the returned polarimetric signals, such as ZDR calculated by Equation (20) and Kdp calculated by Equation (21), yield valuable information regarding hydrometeor size, shape, orientation, and microphysical phase [59]. Diverse rainfall estimators based on dual-polarization radar data have previously been reported, such as R(ZH, ZDR), R(Kdp) [37,46,59,60]. Therefore, the application of polarimetric radar as new rainfall estimators may help reduce the uncertainty due to DSD variability among different parts of the Mei-Yu front.
In this study, C-band polarimetric radar parameters ZH (mm6 m−3), ZDR (dB), and Kdp (deg km−1) were calculated following Zhang et al. [37] for the Brandes drop shape assumption using Parsivel2 sensor data observed during the Mei-Yu season (see Section 2.4). R(ZH, ZDR) and R(Kdp) are usually obtained in the form of R = aZHbZDRc and R = aKdpb respectively. By using the least-squares method, the best-fit results derived from DSD for the R(ZH, ZDR) relationship is
R = 0.0023ZH0.926ZDR−0.97
for the R(Kdp) relationship, the equation is fitted as
R = 19.2Kdp0.70
The comparison results and statistics are shown in Table 5. Compared to the other two, R(ZH, ZDR) performed the best in the Mei-Yu front with an NB of −3.8% and an NAE of 9.2% for NJ, an NB of −4.9% and an NAE of 10.7% for CZ, and an NB of 25% and an NAE of 25.5% for WP. The results indicated that polarimetric estimation was more accurate than the Z–R relationship. We recommend in this work to use R(ZH, ZDR) for estimations of Mei-Yu rainfall. Nevertheless, Chen et al. [10] derived three different rainfall estimators, R(ZH), R(ZH, ZDR), and R(Kdp), from two-year 2DVD observations of Mei-Yu systems and found that R(Kdp) was the most accurate among the three estimators. This is because the results of Chen et al. [10] were obtained only by using DSD data collected in eastern China, while the oceanic part of the entire Mei-Yu front was not considered, which further implies the microphysical variability in different parts of the Mei-Yu front is of great significance.

5. Summary and Conclusions

In the work, we derived scaling law DSD models in different parts of the Mei-Yu front by using the DSD samples measured from Parsivel2 sensors, and DSD-based relations were further derived to improve the accuracy of quantitative radar precipitation estimations. The major conclusions can be drawn as below:
(1)
The spectral width of normalized DSD in NJ stratiform rain was widest among different parts of the Mei-Yu front, resulting in a “size-control” drop size distribution. The max particle concentration of convective rain was largest in WP, indicating a strong oceanic convective rainstorm. The average Dm value in heavy rain was larger in NJ (2.16 mm) and CZ (2.12 mm) than in WP (2.08 mm) of the Mei-Yu front. Given the same Λ, the parameter µ of WP was less than that of NJ (CZ).
(2)
The Z–R relations were estimated by using four methods (STD, LS, EXP, and GAM). The scaled spectra shapes (EXP and GAM) tended to underestimate the Mei-Yu precipitation, whereas the statistical LS method and STD relation usually overestimated the Mei-Yu rainfall. The EXP method had a better performance in stratiform rain than convective rain, which indicates that it is very sensitive to rain rates. The GAM method performed best compared to the other three methods in Mei-Yu rainfall estimations. In comparison with LS estimations, the GAM method showed a considerable improvement in both stratiform (33.9%) and convective (2.8%) rainfall estimations of the Mei-Yu front.
(3)
Several polarimetric variables are calculated by the use of DSD sensor data, such as ZH,V and ZDR. Furthermore, we derived empirical R(ZH, ZDR) and R(Kdp) estimators to improve the polarimetric radar rainfall estimation in the Mei-Yu front and found that the R(ZH, ZDR) estimator demonstrated the most impressive improvement.
Notably, the primary focus herein is on the improvement of Z–R relations in Mei-Yu frontal precipitation with scaled DSD models. Rainfall estimation for polarimetric radar is only briefly discussed, so more research is required to better understand the benefits of calculating polarimetric radar variables via scaled DSD. In addition, the instrumental limitations of Parsivel sensors should be further discussed. The rainfall estimations based on Parsivel2 DSD sensors should be validated by weather radar observations. Meanwhile, some issues might occur from research to operation such as the difference of height between Parsivel and radar. We leave that for future research with effective radar data collected.

Author Contributions

Conceptualization, Z.W.; Data curation, H.Z. and H.L.; Formal analysis, Z.W.; Funding acquisition, L.Z.; Investigation, H.Z.; Methodology, Z.W.; Project administration, L.Z.; Resources, L.Z.; Software, H.Z.; Supervision, Z.W.; Validation, H.Z., Z.W. and L.Z.; Visualization, Y.X.; Writing—original draft, Z.W.; Writing—review & editing, H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The research was financially supported by the National Nature Science Foundation of China (41975066) and the Beijige Open Research Fund for Nanjing Joint Center of Atmospheric Research (NJCAR2018ZD03).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chiew, F.H.S. Estimation of rainfall elasticity of streamflow in Australia. Hydrol. Sci. J. 2006, 51, 613–625. [Google Scholar] [CrossRef]
  2. Guhathakurta, P.; Sreejith, O.P.; Menon, P.A. Impact of Climate Change on Extreme Rainfall Events and Flood Risk in India. J. Earth Syst. Sci. 2011, 120, 359. [Google Scholar] [CrossRef]
  3. Hoerling, M.; Eischeid, J.; Perlwitz, J.; Quan, X.-W.; Wolter, K.; Cheng, L. Characterizing recent trends in U.S. Heavy precipitation. J. Clim. 2016, 29, 2313–2332. [Google Scholar] [CrossRef]
  4. Cui, W.; Dong, X.; Xi, B.; Liu, M. Cloud and precipitation properties of MCSs along the Meiyu frontal zone in central and southern China and their associated large-scale environments. J. Geophys. Res. Atmos. 2020, 125, e2019JD031601. [Google Scholar] [CrossRef]
  5. Hsu, S.-Y.; Chen, T.-B.; Du, W.-C.; Wu, J.-H.; Chen, S.-C. Integrate Weather radar and Monitoring Devices for Urban Flooding Surveillance. Sensors 2019, 19, 825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rosenfeld, D.; Wolff, D.B.; Atlas, D. General probability-matched relations between radar reflectivity and rain rate. J. Appl. Meteorol. 1993, 32, 50–72. [Google Scholar] [CrossRef] [Green Version]
  7. Radhakrishna, B.; Rao, T.N. Differences in cyclonic raindrop size distribution from southwest to northeast monsoon season and from that of noncyclonic rain. J. Geophys. Res. 2010, 115, D16205. [Google Scholar] [CrossRef]
  8. Bringi, V.N.; Chandrasekar, V. Polarimetric Doppler Weather Radar: Principles and Applications; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  9. You, C.H.; Kang, M.Y.; Lee, D.I.; Uyeda, H. Rainfall estimation by S-band polarimetric radar in Korea. Part I: Preprocessing and preliminary results. Meteorol. Appl. 2014, 21, 975–983. [Google Scholar] [CrossRef]
  10. Chen, G.; Zhao, K.; Zhang, G.; Huang, H.; Liu, S.; Wen, L.; Yang, Z.; Yang, Z.; Xu, L.; Zhu, W. Improving Polarimetric C-Band Radar Rainfall Estimation with Two-dimensional Video Disdrometer Observations in Eastern China. J. Hydrometeorol. 2017, 18, 1375–1391. [Google Scholar] [CrossRef]
  11. Gou, Y.; Ma, Y.; Chen, H.; Yin, J. Utilization of a C-band Polarimetric Radar for Severe Rainfall Event Analysis in Complex Terrain over Eastern China. Remote Sens. 2019, 11, 22. [Google Scholar] [CrossRef] [Green Version]
  12. Steiner, M.; Smith, J.A. Reflflectivity, rain rate, and kinetic energy flux relationships based on raindrop spectra. J. Appl. Meteorol. 2000, 39, 1923–1940. [Google Scholar] [CrossRef]
  13. Chapon, B.; Delrieu, G.; Gosset, M.; Boudevillain, B. Variability of rain drop size distribution and its effect on the Z-R relationship: A case study for intense Mediterranean rainfall. Atmos. Res. 2008, 87, 52–65. [Google Scholar] [CrossRef]
  14. Sempere Torres, D.; Porrà, J.M.; Creutin, J.D. A general formulation for raindrop size distribution. J. Appl. Meteorol. 1994, 33, 1494–1502. [Google Scholar] [CrossRef]
  15. Sempere Torres, D.; Porrà, J.M.; Creutin, J.D. Experimental evidence of a general description for raindrop size distribution properties. J. Geophys. Res. 1998, 103, 1785–1797. [Google Scholar] [CrossRef]
  16. Marshall, J.S.; Palmer, W.M.K. The distribution of raindrops with size. J. Meteorol. 1948, 5, 165–166. [Google Scholar] [CrossRef]
  17. Ulbrich, C.W. Natural variations in the analytical form of the raindrop size distribution. J. Appl. Meteorol. 1983, 22, 1764–1775. [Google Scholar] [CrossRef] [Green Version]
  18. Uijlenhoet, R. Parameterization of Rainfall Microstructure for Radar Meteorology and Hydrology. Ph.D. Thesis, Wageningen University, Wageningen, The Netherlands, 1999. [Google Scholar]
  19. Hazenberg, P.; Yu, N.; Boudevillain, B.; Delrieu, G.; Uijlenhoet, R. Scaling of raindrop size distributions and classification of radar reflectivity–rain rate relations in intense Mediterranean precipitation. J. Hydrol. 2011, 402, 179–192. [Google Scholar] [CrossRef]
  20. Sampe, T.; Xie, S.P. Large-Scale Dynamics of the Meiyu-Baiu Rainband: Environmental Forcing by the Westerly Jet. J. Clim. 2010, 23, 113–134. [Google Scholar] [CrossRef] [Green Version]
  21. Kuwano-Yoshida, A.; Taguchi, B.; Xie, S.P. Baiu Rainband Termination in Atmospheric and Coupled Atmosphere–Ocean Models. J. Clim. 2013, 26, 10111–10124. [Google Scholar] [CrossRef] [Green Version]
  22. Ogura, Y.; Asai, T.; Dohi, K. A case study of a heavy precipitation event along the Baiu front in northern Kyushu, 23 July 1982: Nagasaki heavy rainfall. J. Meteorol. Soc. Jpn. 1985, 63, 883–900. [Google Scholar] [CrossRef] [Green Version]
  23. Ninomiya, K. Large- and meso-α-scale characteristics of Mei-Yu/Baiu front associated with intense rainfalls in 1–10 July 1991. J. Meteorol. Soc. Jpn. 2000, 78, 141–157. [Google Scholar] [CrossRef] [Green Version]
  24. Ninomiya, K.; Shibagaki, Y. Cloud system families in the Mei-Yu-Baiu front observed during 1–10 July 1991. J. Meteorol. Soc. Jpn. 2003, 81, 193–209. [Google Scholar] [CrossRef] [Green Version]
  25. Yoshizaki, M.; Kato, T.; Tanaka, Y.; Takayama, H.; Shoji, Y.; Seko, H. Analytical and numerical study of the 26 June 1998 orographic front observed in western Kyushu, Japan. J. Meteorol. Soc. Jpn. 2000, 78, 835–856. [Google Scholar] [CrossRef] [Green Version]
  26. Kato, T.; Aranami, K. Formation factors of 2004 Niigata-Fukushima and Fukui heavy rainfalls and problems in the predictions using a cloud-resolving model. SOLA 2005, 1, 1–4. [Google Scholar] [CrossRef] [Green Version]
  27. Kato, T. Structure of the band-shaped precipitation system inducing the heavy rainfall observed over northern Kyushu, Japan on 29 June 1999. J. Meteorol. Soc. Jpn. 2006, 84, 129–153. [Google Scholar] [CrossRef] [Green Version]
  28. Löffler-Mang, M.; Joss, J. An optical distrometer for measuring size and velocity of hydrometeors. J. Atmos. Ocean. Technol. 2000, 17, 130–139. [Google Scholar] [CrossRef]
  29. Wu, Z.; Zhang, Y.; Zhang, L.; Lei, H.; Xie, Y.; Wen, L.; Yang, J. Characteristics of summer season raindrop size distribution in three typical regions of western Pacific. J. Geophys. Res. Atmos. 2019, 124, 4054–4073. [Google Scholar] [CrossRef] [Green Version]
  30. Gunn, R.; Kinzer, G.D. The terminal velocity of fall for region droplets in stagnant air. J. Meteorol. 1949, 6, 243–248. [Google Scholar] [CrossRef] [Green Version]
  31. Tokay, A.; Petersen, W.A.; Gatlin, P.; Wingo, M. Comparison of raindrop size distribution measurements by collocated disdrometers. J. Atmos. Ocean. Technol. 2013, 30, 1672–1690. [Google Scholar] [CrossRef]
  32. Tokay, A.; Wolff, D.B.; Petersen, W.A. Evaluation of the New Version of the Laser-Optical Disdrometer, OTT Parsivel2. J. Atmos. Ocean. Technol. 2014, 31, 1276–1288. [Google Scholar] [CrossRef]
  33. Atlas, D.; Ulbrich, C.W. Path- and area-integrated rainfall measurement by microwave attenuation in the 1–3 cm band. J. Appl. Meteorol. 1977, 16, 1322–1331. [Google Scholar] [CrossRef] [Green Version]
  34. Tokay, A.; Short, D.A. Evidence from tropical raindrop spectra of the origin of rain from stratiform versus convective clouds. J. Appl. Meteorol. 1996, 35, 355–371. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, G.; Vivekanandan, J.; Brandes, E.A.; Meneghini, R.; Kozu, T. The shape-slope relation in observed gamma raindrop size distributions: Statistical error or useful information? J. Atmos. Ocean. Technol. 2003, 20, 1106–1119. [Google Scholar] [CrossRef] [Green Version]
  36. 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]
  37. Zhang, G.; Vivekanandan, J.; Brandes, E.A. A method for estimating rain rate and drop size distribution from polarimetric radar measurements. IEEE Trans. Geosci. Remote Sens. 2001, 39, 830–841. [Google Scholar] [CrossRef] [Green Version]
  38. 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]
  39. Bringi, V.N.; Chandrasekar, V.; Hubbert, J.; Gorgucci, E.; Randeu, W.L.; Schoenhuber, M. Raindrop size distribution in different climatic regimes from disdrometer and dual-polarized radar analysis. J. Atmos. Sci. 2003, 60, 354–365. [Google Scholar] [CrossRef]
  40. Testud, J.; Oury, S.; Amayenc, P.; Black, R.A. The concept of “normalized” distributions to describe raindrop spectra: A tool for cloud physics and cloud remote sensing. J. Appl. Meteorol. 2001, 40, 1118–1140. [Google Scholar] [CrossRef]
  41. Uijlenhoet, R.; Steiner, M.; Smith, J.A. Variability of rain drop size distributions in a squall line and implications for radar rainfall estimation. J. Hydrometeorol. 2003, 4, 43–61. [Google Scholar] [CrossRef]
  42. Uijlenhoet, R.; Smith, J.A.; Steiner, M. The microphysical structure of extreme precipitation as inferred from ground-based raindrop spectra. J. Atmos. Sci. 2003, 60, 1220–1238. [Google Scholar] [CrossRef]
  43. Jameson, A.R.; Kostinski, A.B. When is rain steady? J. Appl. Meteorol. 2002, 41, 83–90. [Google Scholar] [CrossRef]
  44. Hu, Z.; Srivastava, R. Evolution of raindrop size distribution by coalescence, breakup, and evaporation: Theory and observations. J. Atmos. Sci. 1995, 52, 1761–1783. [Google Scholar] [CrossRef] [Green Version]
  45. Chen, B.; Yang, J.; Pu, J. Statistical Characteristics of Raindrop Size Distribution in the Mei-Yu Season Observed in Eastern China. J. Meteorol. Soc. Jpn. 2013, 91, 215–227. [Google Scholar] [CrossRef] [Green Version]
  46. Cao, Q.; Zhang, G.; Brandes, E.; Schuur, T.; Ryzhkov, A.; Ikeda, K. Analysis of video disdrometer and polarimetric radar data to characterize rain microphysics in Oklahoma. J. Appl. Meteorol. Climatol. 2008, 47, 2238–2255. [Google Scholar] [CrossRef] [Green Version]
  47. Vivekanandan, J.; Zhang, G.; Brandes, E. Polarimetric radar estimators based on a constrained gamma drop size distribution model. J. Appl. Meteorol. 2004, 43, 217–230. [Google Scholar] [CrossRef] [Green Version]
  48. Atlas, D.; Ulbrich, C.W. Drop size spectra and integral remote sensing parameters in the transition from convective to stratiform rain. Geophys. Res. Lett. 2006, 33, L16803. [Google Scholar] [CrossRef]
  49. Thurai, M.; Petersen, W.A.; Tokay, A. Drop size distribution comparisons between Parsivel and 2D video disdrometers. Adv. Geosci. 2011, 30, 3–9. [Google Scholar] [CrossRef] [Green Version]
  50. Thurai, M.; William, C.R.; Bring, V.N. Examining the correlations between drop size distribution parameters using data from two side-by-side 2D-video disdrometers. Atmos. Res. 2014, 144, 95–110. [Google Scholar] [CrossRef]
  51. Campos, E.; Zawadzki, I. Instrumental uncertainties in ZR relations. J. Appl. Meteorol. 2000, 39, 1088–1102. [Google Scholar] [CrossRef]
  52. Fulton, R.; Breidenbach, J.P.; Seo, D.J.; Miller, D.; O’Bannon, T. The WSR-88D rainfall algorithm. Weather Forecast. 1997, 13, 377–395. [Google Scholar] [CrossRef]
  53. Rosenfeld, D.; Ulbrich, C.W. Cloud Microphysical Properties, Processes, and Rainfall Estimation Opportunities. Meteorol. Monogr. 2003, 30, 237–258. [Google Scholar] [CrossRef]
  54. Seed, A.W.; Nicol, J.; Austin, G.L.; Stow, C.D.; Bradley, S.G. The impact of radar and raingauge sampling errors when calibrating a weather radar. Meteorol. Appl. 1996, 3, 43–52. [Google Scholar] [CrossRef]
  55. Wen, L.; Zhao, K.; Zhang, G.; Xue, M.; Zhou, B.; Liu, S.; Chen, X. Statistical characteristics of raindrop size distributions observed in East China during the Asian summer monsoon season using 2-D video disdrometer and Micro Rain radar data. J. Geophys. Res. Atmos. 2016, 121, 2265–2282. [Google Scholar] [CrossRef] [Green Version]
  56. Brandes, E.A.; Vivekanandan, J.; Wilson, J.W. A comparison of radar reflectivity estimates of rainfall from collocated radars. J. Atmos. Ocean. Technol. 1999, 16, 1264–1272. [Google Scholar] [CrossRef]
  57. Steiner, M.; Smith, J.A.; Uijlenhoet, R. A microphysical interpretation of radar reflectivity–rain rate relationships. J. Atmos. Sci. 2004, 61, 1114–1131. [Google Scholar] [CrossRef]
  58. Lee, C.K.; Zawadzki, I. Variability of drop size distributions: Time-scale dependence of the variability and its effects on rain estimation. J. Appl. Meteorol. 2005, 44, 241–255. [Google Scholar] [CrossRef]
  59. Vasiloff, S.V.; Seo, D.-J.; Howard, K.W.; Zhang, J.; Kitzmiller, D.H.; Mullusky, M.G.; Krajewski, W.F.; Brandes, E.A.; Rabin, R.M.; Berkowitz, D.S.; et al. Improving QPE and Very Short Term QPF. Bull. Am. Meteorol. Soc. 2007, 88, 1899–1911. [Google Scholar] [CrossRef] [Green Version]
  60. Cao, Q.; Zhang, G.; Brandes, E.; Schuur, T. Polarimetric radar rain estimation through retrieval of drop size distribution using a Bayesian approach. J. Appl. Meteorol. Climatol. 2010, 49, 973–990. [Google Scholar] [CrossRef]
Figure 1. Locations of three Parsivel2 raindrop size distribution (DSD) sensors employed during the Mei-Yu frontal rain events. The western Pacific (WP) site is located in the southeast ocean area of Japan Sea, while Nanjing (NJ) and Chuzhou (CZ) are located in the East China continent.
Figure 1. Locations of three Parsivel2 raindrop size distribution (DSD) sensors employed during the Mei-Yu frontal rain events. The western Pacific (WP) site is located in the southeast ocean area of Japan Sea, while Nanjing (NJ) and Chuzhou (CZ) are located in the East China continent.
Sensors 20 05257 g001
Figure 2. Himawari satellite image and Japan Meteorological Agency (JMA) surface weather chart at (a) 00 UTC, 2 July (b) 00 UTC, 4 July and (c) 00 UTC, 12 July 2014.
Figure 2. Himawari satellite image and Japan Meteorological Agency (JMA) surface weather chart at (a) 00 UTC, 2 July (b) 00 UTC, 4 July and (c) 00 UTC, 12 July 2014.
Sensors 20 05257 g002
Figure 3. Scaling analysis of the drop size spectra for the Mei-Yu events in July 2014. The exponent γn versus the moment order n + 1 are plotted in three different parts of the Mei-Yu front with two different rain types (left panels (a), (c), and (e) for convective rain, right panels (b), (d), and (f) for stratiform rain). Additionally shown in each panel are the scaling exponents (α, β).
Figure 3. Scaling analysis of the drop size spectra for the Mei-Yu events in July 2014. The exponent γn versus the moment order n + 1 are plotted in three different parts of the Mei-Yu front with two different rain types (left panels (a), (c), and (e) for convective rain, right panels (b), (d), and (f) for stratiform rain). Additionally shown in each panel are the scaling exponents (α, β).
Sensors 20 05257 g003
Figure 4. By using the moment estimation method, the exponential (dashed) and gamma (solid) shapes in three different parts of the Mei-Yu front are obtained and presented, together with the scaled drop size spectra. The coefficients of the distributions [µ, Λ] and the coefficient of determination r2 are presented as well. The left panels (a), (c), and (e) represent convective rain, while the right panels (b), (d), and (f) represent stratiform rain.
Figure 4. By using the moment estimation method, the exponential (dashed) and gamma (solid) shapes in three different parts of the Mei-Yu front are obtained and presented, together with the scaled drop size spectra. The coefficients of the distributions [µ, Λ] and the coefficient of determination r2 are presented as well. The left panels (a), (c), and (e) represent convective rain, while the right panels (b), (d), and (f) represent stratiform rain.
Sensors 20 05257 g004
Figure 5. Evolution of Dm (mm) and log10(Nt) (m−3) versus rain rate (mm h−1) for DSD samples observed in (a) NJ, (b) CZ, and (c) WP. The red circles represent parameter log10(Nt), while the blue rectangles represent parameter Dm. Symbols are plotted at the mean for ordinate within each rain rate class. The color bar at right side represents six classes of rain rate with an interval of 10 mm h−1. The maximum and minimum of two parameters are given in the whiskers of each symbol.
Figure 5. Evolution of Dm (mm) and log10(Nt) (m−3) versus rain rate (mm h−1) for DSD samples observed in (a) NJ, (b) CZ, and (c) WP. The red circles represent parameter log10(Nt), while the blue rectangles represent parameter Dm. Symbols are plotted at the mean for ordinate within each rain rate class. The color bar at right side represents six classes of rain rate with an interval of 10 mm h−1. The maximum and minimum of two parameters are given in the whiskers of each symbol.
Sensors 20 05257 g005
Figure 6. µ–Λ scatterplots and fitting curves based on Parsivel2 DSD sensor observations in three parts of the Mei-Yu front. The gray lines correspond to the relationship ΛDm = 4 + μ given the values of Dm = 1.0, 1.5, 2.0, and 3.0 mm. The gray circles (crosses) represent precipitation in NJ (CZ), and the gray triangles represent precipitation in WP. The dashed line and dotted line represent the empirical μΛ relations reported by Zhang et al. [35] and Chen et al. [45], respectively.
Figure 6. µ–Λ scatterplots and fitting curves based on Parsivel2 DSD sensor observations in three parts of the Mei-Yu front. The gray lines correspond to the relationship ΛDm = 4 + μ given the values of Dm = 1.0, 1.5, 2.0, and 3.0 mm. The gray circles (crosses) represent precipitation in NJ (CZ), and the gray triangles represent precipitation in WP. The dashed line and dotted line represent the empirical μΛ relations reported by Zhang et al. [35] and Chen et al. [45], respectively.
Sensors 20 05257 g006
Figure 7. Z-R scatterplots and fitting curves for Mei-Yu precipitation in (a) NJ, (b) CZ, and (c) WP. The gray circles indicate convective precipitation samples, the gray dots indicate stratiform precipitation samples, the solid red line indicates the convective precipitation fitting curve, the solid blue line indicates the stratiform precipitation curve, the solid black line indicates the total precipitation curve, and the two dashed lines indicate the empirical relation Z = 368R1.21 for convective precipitation suggested by Chen et al. [45] and Z = 193.73R1.54 for stratiform rain obtained by Wen et al. [55].
Figure 7. Z-R scatterplots and fitting curves for Mei-Yu precipitation in (a) NJ, (b) CZ, and (c) WP. The gray circles indicate convective precipitation samples, the gray dots indicate stratiform precipitation samples, the solid red line indicates the convective precipitation fitting curve, the solid blue line indicates the stratiform precipitation curve, the solid black line indicates the total precipitation curve, and the two dashed lines indicate the empirical relation Z = 368R1.21 for convective precipitation suggested by Chen et al. [45] and Z = 193.73R1.54 for stratiform rain obtained by Wen et al. [55].
Sensors 20 05257 g007
Table 1. Mei-Yu frontal precipitation events detected in three different parts of the Mei-Yu front.
Table 1. Mei-Yu frontal precipitation events detected in three different parts of the Mei-Yu front.
DateRegions1-min Samples (min)Accumulated Precipitation (mm)Max Rain Rate (mm h−1)Max Reflectivity Factor (dBZ)
2 July 2014NJ82834.4136.2644.42
CZ6109.7426.4741.95
WP4548.605.7448.22
4 July 2014NJ123895.1986.8353.52
CZ1112107.3671.9049.31
WP1045.9110.5416.19
12 July 2014NJ54542.7388.2256.94
CZ58421.1822.3348.26
WP9510.7948.6452.76
Table 2. Sample distribution of precipitation in different rain rate classes in three parts of Mei-Yu front.
Table 2. Sample distribution of precipitation in different rain rate classes in three parts of Mei-Yu front.
RegionsRain Rate Classes (mm h−1)
0–1010–2020–3030–4040–5050–60
NJ203028466332519
CZ207313545231513
WP468813613119
Table 3. Parameter values of the radar Z = ARb relations based on three different rainfall estimation methods.
Table 3. Parameter values of the radar Z = ARb relations based on three different rainfall estimation methods.
Rain TypeRegionsLSEXPGAM
AbAbAb
convectiveNJ1281.622381.811161.81
CZ1151.627291.671431.67
WP1421.763861.801681.80
stratiformNJ3081.5015571.657481.65
CZ2751.5212041.683651.68
WP4891.178561.745761.74
Table 4. NB and NAE errors (%) of four rainfall estimation methods are obtained via Equations (23) and (24).
Table 4. NB and NAE errors (%) of four rainfall estimation methods are obtained via Equations (23) and (24).
Rain TypeRegionsNAENB
STDLSEXPGAMSTDLSEXPGAM
convectiveNJ40.119.130.317.916.60−29.8−5.4
CZ18.314.336.313.114.20−34.3−7.9
WP48.715.937.615.533.80−35.7−11.4
stratiformNJ52.433.631.524.943.80−27.3−18.4
CZ43.325.724.615.823.40−21.9−13.2
WP44.830.647.815.347.80−33.5−15
Table 5. NB and NAE errors (%) of Z–R (LS and GAM method), R(ZH, ZDR), and R(Kdp) for precipitation in the three parts of the Mei-Yu front.
Table 5. NB and NAE errors (%) of Z–R (LS and GAM method), R(ZH, ZDR), and R(Kdp) for precipitation in the three parts of the Mei-Yu front.
ParametersRegionsZR (LS)ZR (GAM)R(ZH, ZDR)R(Kdp)
NAENJ27.625.39.217.9
CZ28.523.910.721.1
WP39.835.425.534.4
NBNJ0−11.7−3.8−9.7
CZ0−13.5−4.9−11.3
WP0−25.42532.2

Share and Cite

MDPI and ACS Style

Zheng, H.; Wu, Z.; Zhang, L.; Xie, Y.; Lei, H. Improving Radar Rainfall Estimations with Scaled Raindrop Size Spectra in Mei-Yu Frontal Rainstorms. Sensors 2020, 20, 5257. https://doi.org/10.3390/s20185257

AMA Style

Zheng H, Wu Z, Zhang L, Xie Y, Lei H. Improving Radar Rainfall Estimations with Scaled Raindrop Size Spectra in Mei-Yu Frontal Rainstorms. Sensors. 2020; 20(18):5257. https://doi.org/10.3390/s20185257

Chicago/Turabian Style

Zheng, Hepeng, Zuhang Wu, Lifeng Zhang, Yanqiong Xie, and Hengchi Lei. 2020. "Improving Radar Rainfall Estimations with Scaled Raindrop Size Spectra in Mei-Yu Frontal Rainstorms" Sensors 20, no. 18: 5257. https://doi.org/10.3390/s20185257

APA Style

Zheng, H., Wu, Z., Zhang, L., Xie, Y., & Lei, H. (2020). Improving Radar Rainfall Estimations with Scaled Raindrop Size Spectra in Mei-Yu Frontal Rainstorms. Sensors, 20(18), 5257. https://doi.org/10.3390/s20185257

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