Characteristics of Raindrop Size Distribution on the Eastern Slope of the Tibetan Plateau in Summer

: Precipitation microphysics over the Tibetan Plateau (TP) remain insufficiently understood, due to the lack of observations and studies. This paper presents a comprehensive investigation of the raindrop size distribution (DSD) for rainfall that happened on the eastern slope of TP in summer. DSD differences between different rain types and under different rain rates are investigated. Confidential empirical relationships between the gamma shape and slope parameters, and between reflectivity and rain rate are proposed. DSD properties in this area are also compared with those in other areas. The results indicate that the stratiform and convective rains contribute to different rain duration and amount, with diverse rainfall macro ‐ and microphysical properties. The rain spectra of two rain types can become broader with higher concentrations as the rain rate increases. DSDs in this area are different to those in other areas. The stratiform DSD is narrower than that in the non ‐ plateau area. The two rain types of this area both have higher number concentrations for 0.437–1.625 mm raindrops than those of the mid ‐ TP. The relationships of shape– slope parameters and reflectivity–rain rate in this area are also different from those in other areas. The rain spectra in this area can produce a larger slope parameter under the same shape parameter than in the mid ‐ TP. The convective rain can produce a smaller rain rate under the same reflectivity. The accuracy proposed reflectivity–rain rate relationship in application to quantitative rainfall estimation is also discussed. The results show that the relationship has an excellent performance when the rain rate exceeds 1 mm h ‐ 1 . averages of 𝐷 𝑙𝑜𝑔 . The results show that the convective rain has a larger mean 𝐷 (cid:3040) and a comparable 𝑁 (cid:3050) compared with those of the stratiform rain. The averages of 𝐷 𝑙𝑜𝑔 for the stratiform and convective rains are 0.932–3.468 and 1.5–3.437, respectively. The black stars, blue squares, green circles and purple pluses are averages observed in other four climatic regimes in China. Comparisons of the five pairs of averages indicates that the 𝐷 𝑠(cid:4667) for two rain types observed in our area are both different from those observed in other areas. For the stratiform rain, the mean 𝐷 (cid:3040) of our area is the same as that of the mid ‐ TP, but smaller than those of non ‐ plateau areas. The mean 𝑁 (cid:3050) of our area can be larger than that of North China, but smaller than those of other three areas. For the convective rain, the mean 𝐷 (cid:3040) of our area is slightly larger than that of North China but smaller than other areas. The mean convective 𝑁 (cid:3050) of our area is larger than that of North China, but smaller than those of other areas.


Introduction
The Tibetan Plateau (TP) is the largest and highest plateau on the planet. It plays a vital role in the Asian climate and weather, owing to its powerful dynamic and thermal effects on the tropospheric atmosphere [1,2]. In summer, rainfall events over the TP are frequent, however, the reliability of numerical simulations over the TP is much lower than that in other plain areas, due to the lack of observations and insufficient knowledge of the microphysical processes and mechanisms of the precipitation. Thus, further observing and studying the microphysical characteristics of precipitation over the TP is of great importance.
As a vital aspect of precipitation microphysics, the raindrop size distribution (DSD) can reflect the number concentration of raindrops as a function of the diameter. The DSD can be directly measured by filter paper or disdrometer on the ground. Although the ground measurements cannot be the perfect representation of the real situation of the precipitation up in the air, the long-period DSD dataset, to some extent, can reveal the microphysical properties of different kinds of precipitation in different geographical locations [3][4][5][6][7][8][9][10][11][12][13][14][15][16]. Previous findings showed that different rain types possess different DSDs, and the convective rain commonly has wider raindrop spectra than the stratiform rain, with a larger mean diameter and a higher total number concentration [17][18][19]. The maritime-like convection also generally possesses larger diameters and smaller number concentrations than the continental-like convection [4]. DSDs vary in space and time and over a range of scales, they can be apparently different, even in the same synoptic system, as a result of diverse water vapor and dynamic conditions [20][21][22]. Additionally, the observed DSD can be used to formulate the raindrop spectrum. In recent years, the most widely used mathematical formulation of DSDs is the three-parameter gamma model proposed by Ulbrich (1983) [7]. Fitting the gamma parameters and investigating their relationships in specific climatic regions can be useful to optimize the microphysical parameterization schemes of the regional numerical models [23][24][25]. A conventional technique for single-polarization weather radar to estimate the precipitation typically relies on the power law relationship between reflectivity (Z, mm 6 m -3 ) and rain rate (R, mm h -1 ) [26,27]. The Z and R are proportional to the sixth and third power of the raindrop diameter. Thus, investigating the DSD properties and deriving accurate coefficients of the Z-R relationship is also crucial for improving radar quantitative precipitation estimation (QPE).
For the DSD study over the TP, limited works have been conducted because of the lack of observational data. Until 2013, an integrated observational experiment called the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) was launched by the China Meteorological Administration and other institutions [2]. One of the specific scientific goals of the TIPEX-III was to achieve a better understanding of microphysical processes of the cloud and precipitation over the TP [2,28]. The TIPEX-III deployed ground-based laser disdrometers for continuously observing rainfall properties and DSDs in the mid-TP (Nagqu, 29.04° N, 100.30° E, 4,507 m asl) from July to August in 2013, 2014 and 2015. Based on the measurement collected from this experiment, rainfall properties have been studied by Chang et al. (2016), Chen et al. (2017) and Wu et al. (2017) [19,29,30]. The results concluded that raindrop spectra over the TP can be wider than those in basin areas, the widest spectra occur during 1700-1800 LST, due to the strong convection. The DSD exhibits a significant day/night variation; the daytime convection has broader DSDs with fewer small drops but more large drops than the nighttime case. The raindrop number concentration of stratiform rains in the mid-TP also has a more significant decrease with the increased diameter than that in South China, but convective rains can be the opposite. For precipitation physical processes over the mid-TP, Porcù et al. (2014) found that collisional breakup takes place for relatively lower rain rates than those at low altitudes [31].
Several abovementioned works have preliminarily revealed the DSD characteristics of rainfalls in the mid-TP. However, those findings cannot represent the whole situation on the TP, because the plateau has very complex geography and land features, which can lead to different rainfall properties in different areas. For example, compared to the mid-TP, the terrain and climate features of a transitional region of interest in this study, which locates between the plateau and the downstream basin at the eastern slope of TP, are quite different. Except for the rainfall events caused by similar synoptic systems as in the mid-TP, many precipitations can also be produced by local special synoptic systems, such as southwest vortex. More importantly, the rainfall microphysical properties over the eastern slope of TP have not been well investigated. In this paper, we attempt to present a comprehensive study of the DSD for rainfall that happened in this area in summer. In detail, the observation experiment, instrument and measurements, followed by methods of data quality control (QC) and rain type classification, are described in Section 2. The DSD characteristics for different rain rates and types are studied in Section 3. Relationships between three gamma parameters and between DSD integral quantities are also proposed. In Section 4, the results are compared with other studies and the uncertainties of the used methods are discussed. The conclusions are summarized in Section 5.

Experiment, Instrument, and Measurements
To enhance the observation of cloud and precipitation over the eastern slope of TP, an observation experiment was performed by the Chinese Academy of Meteorological Sciences at Daocheng (29.04° N, 100.30° E-, 3800 m asl), China from 10 July to 31 August in 2016, as a continuation of the TIPEX-Ⅲ. A key precipitation observation instrument deployed on the site is a ground-based disdrometer. Figure 1 shows the experiment location and the appearance of the disdrometer. The disdrometer is a second-generation particle size and velocity (Parsivel) disdrometer, manufactured by OTT HydroMet, Germany. It is a laser-based disdrometer that can simultaneously detect the size and falling velocity of precipitation particles. The Parsivel transmits a parallel horizontal laser beam with a sampling area of 54 cm 2 (18-cm long, 3-cm wide), and a sampling time interval of 60 s. When particles pass through the laser beam, the received signals are attenuated and used to estimate the particle diameter; the duration can be utilized to detect the falling velocity. Raindrops larger than 1 mm are not spherical in Earth's atmosphere, so the estimation of the diameter is based on different axis ratio relationships. Specifically, for raindrops smaller than 1 mm and greater than 5 mm, the axis ratios are set to 1 and 0.7, respectively, and for raindrops between 1 and 5 mm, the axis ratio linearly varies from 1 to 0.7 [32,33]. After the diameter and falling velocity are measured, the raindrops are subdivided into 32 non-equidistant diameter classes, ranging from 0.062 to 24.5 mm and 32 non-equidistant velocity classes ranging from 0.05 to 20.8 m s −1 . Detailed information about each diameter and velocity class can be seen in Yuter et al. (2006) [34].

Data QC and Rainfall Physical Quantity Calculation
Apart from biases caused by the inherent limitations of Parsivel, several problematic data can be improved using the following QC technologies.
First, the data in the first two diameter classes were abandoned, due to their poor reliability of low signal-to-noise ratios, namely, the minimum detectable diameter of the raindrop was limited to approximately 0.312 mm [33]. Second, any raindrop with a diameter larger than 8 mm was also eliminated, since it can mostly break up in the low level of Earth's atmosphere. Third, considering the realistic instrument sensitivity, samples with a total drop number less than 10 or a rain rate smaller than 0.002 mm h −1 were also removed [19]; otherwise, they were regarded as valid rain samples. Fourth, as proposed by Yuter et al. (2006) and , any raindrop that possessed a normal falling velocity (or diameter) but with an unmatched (excessively larger or smaller) diameter (or falling velocity) was treated as unrealistic data. It can be produced when a large raindrop partially passes through or multiple raindrops parallelly pass through the laser beam; it can also be induced by a strong wind shear near the ground or splashing of raindrops on the instrument surface during heavy rainfalls [34,38]. This kind of problematic data was recognized by comparing the Parsivel-measured result with the theoretical diameter/falling velocity relationship shown in Equations (1) and (2) [60]. Therefore, any measured raindrop with a diameter greater than 1 mm was removed when it was ±60% outside of the theoretical value.
where D (mm) and (m s −1 ) are the particle diameter and terminal falling velocity, respectively; h (m) is the sampling height above sea level; and ℎ is a correction factor.
The raindrop numbers in different diameter and velocity classes before and after QC for the entire data were counted and compared, as shown in Figure 2, to illustrate the QC effect. A small part of the raindrops that possessed an abnormal diameter (velocity), and some that were scattered or excessively large were removed. These unrealistic data accounted for only 1.76% of all raindrops, but they can occupy prominent proportions in the large-diameter and small-velocity classes. For the raindrops in classes with diameters from 2.75 to 7.5 mm, the problematic data accounted for 10.3% to 68%; for the raindrops in classes with velocities from 0.15 to 2.6 m s −1 , the problematic data accounted for 34.87% to 94.59%. After QC, the DSD and its integral rainfall physical quantities, namely, including raindrop total number concentration (m −3 ), rain rate (mm h −1 ), rainwater content (g m −3 ), reflectivity (mm 6 m −3 ), mass-weight mean diameter (mm) and generalized intercept parameter (m −3 mm −1 ), can be further calculated as where i and j are the orders of the diameter and velocity classes, respectively; ∆ (mm) is the interval of the diameter class; is the raindrop number within the diameter class i and velocity class j; (mm s −1 ) is the falling velocity of class j; ∆ (60 s) is the sampling time interval; is the effective sampling area for ith diameter class; and (1 g cm −3 ) is the water density. (m −3 mm −1 ) is the number concentration of raindrops per unit volume per millimeter. Note that the of the Parsivel is associated with the raindrop size and can be estimated as 180 30 0.5 [32,33]. The commonly used gamma distribution proposed by Ulbrich (1983) was implemented to fit the measured DSD, which can be expressed as follows [7], where ( ) is the intercept parameter related to the total number of raindrops. is the shape parameter; a positive (negative) indicates that the raindrop spectrum curve bends upwards (downwards). Λ (mm −1 ) is the slope parameter; a small (large) Λ can correspond to a flat (steep) spectrum. Several methods, such as the moment method (MM), the truncated moment method (TMM), the least square method (LS) and the maximum likelihood method (ML), have been used to estimate the DSD gamma parameters. Studies reported that the maximum likelihood method generally has a better performance than the others, and can be less sensitive to the DSD truncation problem and the short integration time of instrument [9,49,53,[61][62][63]. Therefore, in this study, the maximum likelihood method was implemented to estimate the , , and Λ.

Rain Type Classification
To study the DSD characteristics for the different rain types, the total dataset was further classified into two subsets, namely, stratiform and convective rains, based on the time series of the rain rate. Several classification schemes based on disdrometer measurements have been reported for the separation of convective versus stratiform rains, such as those of Tokay et al.  [4,12,64,65]. In the current study, Bringi's scheme was used to separate the two rain types. Specifically, set a moving window containing 11 rain samples in the time series. The standard deviation of rain rate of samples was calculated. If the was no more than 1.5 mm h −1 , the rain type of the center was determined as stratiform; otherwise, it was treated as convective. According to this classification principle, the window center was moved to the next sample every time until all samples in time series were identified. Note that the first five and last five samples, which cannot be set to the window centers, were judged as the same rain type as the nearest window centers.

Results
In the following subsections, the DSD characteristics for different rain rates and types were investigated in detail. Subsequently, empirical relationships between gamma parameters and rainfall physical quantities for precipitation occurring on the eastern slope of the TP are proposed.

DSD Characteristics for Different Rain Rates
DSD properties can vary largely under different rain rates [19,66]. Therefore, in this work, the total dataset was stratified into five rain-rate regimes, namely, 5 10 ℎ (R4) and 10 ℎ (R5), to investigate the DSD characteristics for different rain rates. Figure 3 presents the accumulated rain amounts and durations for the five regimes. The results show that the ground rainfalls are dominated by light precipitation more than half of the time, with 65.15% of the total duration being produced by R1 and R2, whereas they account for only 12.51% of the total rain amount. The largest contributor to the total rain amount is R3, which accounts for 49.67% of the total rain amount and 31.14% of the total duration. For rainfalls with larger rain rates, R3 and R4 contribute only 3.71% of the total duration; nevertheless, they can produce 37.82% of the total rain amount. . Accumulated durations (blue) and rain amounts (red) for the five rain-rate regimes. The percentages indicate the contributions of each rain-rate regime to the total duration and total rain amount.
The DSDs for the five rain-rate regimes were shown in Figure 4, in which the boxplots illustrate the raindrop number concentration in the different diameter classes. The results reveal that raindrop spectra can evidently change with the rain rates. Specifically, the raindrop spectra gradually become wider and flatter as the rain rate increases, owing to the contribution of large raindrops; the maximum measured diameters for R1 to R5 are 1.375, 2.75, 3.75, 4.25 and 7.5 mm. The raindrop number concentration in each diameter class also increases with the rain rate. The spectrum shapes for the five rain-rate regimes are both unimodal and generally obey the gamma distribution. The peaks of the medians both appear in the fifth diameter class, with a mean diameter of 0.562 mm. The number concentration of the raindrops in each diameter class can vary within three or more orders of magnitude, while raindrops between the 25th and 75th percentiles basically concentrate within one order of magnitude. percentiles, respectively; the red central crossbar represents the median. The black dashed whiskers extend to the maximum and minimum of the data points, respectively, and the outliers are plotted individually using the "+" symbol.
For a convenient comparison of the DSD characteristics under different rain rates, average raindrop spectra for R1 to R5 are superimposed on the same graph in Figure 5. The corresponding integral rainfall physical quantities derived from the average spectra and the fitted gamma parameters are listed in Table 1. The comparison shows that the total number concentration , reflectivity Z, rainwater content W and mass-weighted mean diameter increase with the rain rate. The largest generalized intercept parameter is found in R3. For the gamma parameters, , and Λ decrease with an increased rain rate, corresponding to the widening and flattening of the raindrop spectra. The integral rainfall parameters derived from these averaged spectra are presented in Table 1. The circles represent the center of the Parsivel diameter classes.

DSD Characteristics for Different Rain Types
In this subsection, the DSD characteristics for stratiform and convective rains are analyzed and compared. During the observation period, there were 98 rainfall events, and 4599 DSD samples were obtained, of which 4217 (91.69%) were stratiform samples and 382 (8.31%) were convective samples, respectively. The accumulated rain amounts for the stratiform and convective samples are 66.37 and 47.45 mm, with mean rain rates of 0.94 and 7.45 mm h −1 , respectively. Indicating that the rainfalls that happened in the Daocheng site are light stratiform rains most of the time. Meanwhile, although relatively sparse and short-lived, heavier convective rains can also produce a comparable rain amount.
The statistical DSDs for the stratiform and convective rains are shown in Figure 6. A comparison of the two sets of raindrop boxplots demonstrates that the convective rain has wider and flatter spectra than the stratiform rain, as a result of the contribution of larger raindrops; the maximum raindrop diameters for the convective and stratiform rains are 5.5 and 3.75 mm, respectively. The peaks of the medians for two rain types are both in the fifth diameter class, with a mean diameter of 0.562 mm. These comparison features of the two rain types are quite similar to the comparison results of the raindrop spectra between the low-rain-rate regime and the high-rain-rate regime, as concluded in Subsection 3.1. Figure 7 presents the probability distributions of values of DSD-derived R, Z, W, , , , , and Λ to discern the differences in the rainfall physical quantities and gamma parameters of two rain types. Their corresponding 5th, 25th, 75th, 95th percentiles, medians and averages are also summarized in Table 2. It can be found that, generally, the stratiform rain has smaller values of R, Z, W, and , but larger values of , and Λ than the convective rain. The of two rain types are very close.    The average DSDs and gamma distribution curves fitted by the ML for two rain types are further presented in Figure 8. Table 3 lists the corresponding rainfall physical quantities and gamma parameters for the average spectra. The results show that the average DSD of convective rain is broader and flatter than that of the stratiform rain. The former also possess higher number concentrations than the latter at each diameter class. For comparison, fitted results using the MM are also superimposed on the figure. Comparisons indicate that the MM significantly underestimates the raindrop number concentration of small raindrops, and the biases of ML are much smaller. On the contrary, for medium raindrops, the ML overestimates the results, and the MM can perform better. The correlation coefficients between the observed and fitted DSDs of the stratiform rain for the ML and MM are 0.9426 and 0.9575, respectively, and the counterparts of the convective rain are 0.9132 and 0.9574, respectively. We also use the Kolmogorov-Smirnov (KS) test to verify the fitting methods. The results show that two methods both can pass the verification. The KS p-values of ML and MM for the stratiform rain are both 0.9123, and they are 0.3841 and 0.9123, respectively, for the convective rain. Overall, the gamma distribution can adequately describe the observed rain spectra, and the fitted results for the average DSDs of two rain types can be expressed as N D 1.872 10 D . exp 3.947 , convective rain (12)

Relationship of
The gamma parameters are not mutually independent. The relations between and /Λ can provide useful information reflecting the behavior of the rain spectra, and facilitate their application in other areas. Empirical relationships between the gamma distribution parameters have been presented by researchers. Ulbrich and Atlas (1998) found that the relationship can be expressed by an exponential equation [67]. Brandes et al. (2003) proposed a second-order polynomial Λ relationship, which can be written as [68] Λ 0.0365 0.735 1.935 (13) Many other scientists have reported similar equations, but their studies showed different coefficients; such variation can be directly related to the studied geographical locations and climatological regimes [5,9,11,12,19,35,41]. Some criteria are used before data fitting to obtain more reliable for selecting fitted samples. They considered that the fitted gamma parameters of DSDs with drop counts less than this threshold may be less representative, because of the limitation of fitting approaches [3,19,35,41]. Chen et al. (2017) set the threshold to 300 since the raindrop number concentration over the TP is relatively lower [19]. The relationships of and Λ are also fitted in this study. The results show that there is a good correlation between the and Λ with a certain quadratic relationship, however, the relation between and is unapparent. Figure 9 shows the scatterplots of the and Λ. The gray plots represent the entire data, and the black plots are samples after selection by a threshold of total raindrop number proposed by Chen et al. (2017). The black plots can distribute more concentratedly samples, and have higher representativeness than the entire data. The fitting result of the black plots exhibits a confidential Λ relationship. The correlation coefficient between the observed and fitted Λ is 0.9522, and the fitted result can pass the KS-test with a high p-value of 0.9911. The Λ relationship can be expresses as The relationship of Λ for DSDs observed in the mid-TP that were proposed by Chen et al.
(2017) is compared (the blue solid line). Our relationship is quite consistent with Chen's result, but our formula can produce a larger Λ under the same .

Relationships of and
The and are independent quantities that can synthetically represent the size and number concentration of the rain spectrum. Tables 1 and 3 preliminarily indicate that and can remarkably change with the rain rate R and show differences with respect to rain types. Therefore, relationships between and R are further discussed in this subsection. In many studies, a power-law equation, namely, / , was used to formulize the and relationships, but with different coefficients of A and B, due to the different geographical locations [19,35,69,70]. In this study, we also use the power-law formula to fit the relationships of and . Figure 10 presents the scatterplots of R versus / and fitted curves for two rain types. The and for samples can both increase along with R, while they exhibit different incremental features for two rain types. Specifically, the stratiform rain has larger increments of and at small R values. As R increases, the and of can more easily tend to be stable than those of the convective rain. This difference can imply that the raindrop formation and growth can be efficient under a small rain rate for the stratiform rain. In contrast, as the rain rate increase, physical processes can gradually tend to an equilibrium to yield raindrop spectra with small increments in both size and concentration. The convective rain can continually produce larger raindrops, along with an increase in the rain rate, indicating that the raindrop growth process can remain efficient, as a result of abundant water vapor and strong updrafts.
We use the KS test to verify the fitted relationships of and . The results show that the fitting of relationship is basically confidential, while the correlation between and R for both two rain types is unapparent. The fitted relationships for two rain types are given as versus R for the stratiform rain, (b) versus R for the convective rain, (c) versus R for the stratiform rain, (d) versus R for the convective rain. The blue and red lines are the results for the stratiform and convective rains fitted using the linear least squares method.

Relationship of
The relationship between Z and R has a substantial significance for improving radar QPE. Many studies have investigated the power-law relationship of for stratiform and convective rains using disdrometers, and proposed specific coefficients of A and b for different geographical locations [30,41,65,[71][72][73][74][75][76]. Typically, the relationship of 300 . for convective rains that was proposed by Fulton et al. (1998) is used in the Next-Generation Weather Radar of the United States [73]. The relationship of 200 . presented by Marshall and Palmer (1963) has also been commonly applied to mid-latitude areas for stratiform rains [72]. Herein, we also investigate the relationships for the different rain types on the eastern slope of the TP. Figure 11 shows the statistical results of Z versus R and the fitted relationships. The scatterplots of the samples indicate that Z and R well obey the power-law relationships for both rains. The correlation coefficients between the fitted and observed results for stratiform and convective rains are 0.8729 and 0.968. The fitted relationships of can pass the KS test, and they can be expressed as 149 R . , 107.7 R . ,

Discussion
In this section, we compare our results regarding DSD characteristics with those obtained in other areas of China. Biases of the ML method and the proposed relationship are also discussed.

Comparison of DSD Characteristics
(1) Average DSDs: average DSDs of the stratiform and convective rains observed in our area are compared with results obtained in the mid-TP, East China and South China. Figure 12 shows that both DSDs of two rain types on the eastern slope of TP are different from those in the other areas. For the stratiform rain, spectrum in our area is narrower than that in other areas. The number concentrations of small raindrops (<1.4 mm) are closest to those in South China, followed by East China, and they are apparently higher than those in the mid-TP. However, for larger raindrops, their number concentrations in our area are smaller than those in the other three regions. For the convective rain, the DSD in our area is closest to that in the mid-TP. Their differences are also significant. Convective DSDs in our area generally produce higher number concentrations of small and medium raindrops (<2 mm) than those in the mid-TP. On the contrary, the number concentrations of larger raindrops of the former are lower than those of the latter. (2) domain: Figure 13 shows the scatterplots of versus for two rain types observed in our area. The gray points represent all samples. The red crosses represent the averages of . The results show that the convective rain has a larger mean and a comparable compared with those of the stratiform rain. The averages of for the stratiform and convective rains are 0.932-3.468 and 1.5-3.437, respectively. The black stars, blue squares, green circles and purple pluses are averages observed in other four climatic regimes in China. Comparisons of the five pairs of averages indicates that the for two rain types observed in our area are both different from those observed in other areas. For the stratiform rain, the mean of our area is the same as that of the mid-TP, but smaller than those of non-plateau areas. The mean of our area can be larger than that of North China, but smaller than those of other three areas. For the convective rain, the mean of our area is slightly larger than that of North China but smaller than other areas. The mean convective of our area is larger than that of North China, but smaller than those of other areas.  relationship: for comparison, the relationships for two rain types in the midand southern TP proposed by Wu et al. (2017) and Porcù et al. (2014) are also drawn in Figure 11. Comparisons show that the stratiform relationship of our area is very close to the other results. For a 30-dBZ (10 3 mm 6 m −3 ) Z, R values for our area, the mid-TP, and the southern TP are 3.38, 3.59 and 3.85 mm h −1 . For the convective rain, the relationship of our area is different from those of the other two areas, but relatively closer to that of the mid-TP. For a 40-dBZ (10 4 mm 6 m −3 ) Z, R values for our area, the mid-TP, and the southern TP are 12.57, 15.07 and 21.65 mm h −1 . Overall, the relationship of convective rain has a smaller A, but a larger b than the stratiform rain.

Biases of ML Method and Relationship
To evaluate the performance of ML method in application to DSD fitting, the relative error (RE) of the gamma-fitted compared to the observed DSDs, and the derived integral rainfall quantities are further calculated. Figure 14a,b shows boxplots of REs of the raindrop number concentration in different diameter classes for two rain types. The RE of , , , , , and are depicted in Figure 14c,d, respectively. The results show that most of the gamma-fitted DSDs overestimate the number concentration of small raindrops at 0.312 mm diameter class for both of the two rain types; whereas, the biases for raindrops at larger diameter classes rapidly decrease. For the stratiform rain, the REs ranging from 0.437-to 1.375 mm are relatively smaller than other classes; their 25th-75th percentiles and averages of REs are both within 50%; for the larger raindrops at 1.625-to 3.25 mm classes, the underestimation becomes remarkable. For the convective rain, the overestimation at 0.312 mm class is smaller compared to the stratiform rain; the 25th-75th and averages of REs of a wider scope of diameter classes (0.437 to 2.375 mm classes) are within 50%. For the six DSD-derived quantities, gamma-fitted DSDs underestimate the , , and , but overestimate the and for both rain types. The REs of are larger than the other quantities, indicating that the gamma function fitted by ML method can produce steeper raindrop spectra than the reality. The RE of relationship-derived compared to the observed rain rates (Rs) are also calculated, to estimate the biases of the proposed relationship in the application to rainfall quantitative estimation. Figure 15 shows the boxplots of RE of the rain rate in different regimes for two rain types. Note that the R regimes for two rain types are assigned to be different because of their different value ranges. Obviously, with the increase of R, the REs of both two rain types are decreased. The biases in the first two R regimes (R 1 mm h -1 ) are much smaller than those in the posterior R regimes. For the rainfalls with rain rates greater than 1 mm h -1 , the boxplots of RE of stratiform rain possess 25th-75th percentiles and averages within -27%-18%, while -45%-30% can be found for the convective rain type. Overall, the proposed relationship can produce smaller biases for the stratiform rain than that for the convective rain, and it is more accurate for the rainfalls with a rain rate greater than 1mm h -1 .

Conclusions
The precipitation microphysics in different regions of the TP can be significantly different, because of the enormous and complex terrains and different climates. However, only limited observations and studies have reported the rainfall microphysics in the mid-TP. In this paper, we present a comprehensive study of characteristics of DSDs for rainfalls that happened on the eastern slope of TP. The main findings can be summarized as follows.
(1) As many studies revealed, DSD differences between the stratiform and convective rains and under different R regimes are also apparent on the eastern slope of TP. The stratiform rain contributes to more rain duration and amount compared with the convective rain. Whereas, it has much smaller R, Z, W, and , but larger , and Λ, and produces narrower and steeper rain spectra than the convective rain. As R increase, both DSDs of two rain types become wider, corresponding to occurrences of larger raindrops and higher concentrations.
(2) The DSD characteristics on the eastern slope of TP can be different from those in other areas. The stratiform DSD in this area is narrower than in non-plateau regions, because of the smaller rainfall intensity. The convective DSD in this area also can be narrower than that in East China, but broader than that in South China. The differences of rain spectra between this area and the mid-TP are also significant. For stratiform spectra in this area, number concentrations of raindrops (within 0.437-1.625 mm) are higher than those in the mid-TP. On the contrary, for raindrops smaller than 0.437 mm and larger than 1.635 mm, the former has lower number concentrations than the latter. The differences of the convective spectra between two areas are similar to those of the stratiform spectra.
(3) The confidential relationships of Λ and for rainfalls in this area are proposed with KS test. Comparisons of these relationships with the counterparts in other areas indicate that, rain spectra in this area can produce a lager Λ under the same than rain spectra in the mid-TP. The relationship of stratiform rain is very close to that in the mid-TP, whereas, the counterpart of convective rain can produce a smaller R under the same Z. stratiform rain is slightly smaller than the latter with the same . The former and of the convective rain are both smaller than those of the latter.
(5) Estimation of the gamma fitting using ML method indicated that the gamma fitted DSDs could overestimate (underestimate) the number concentration of small (large) raindrops. The bias of relationship in application to QPE can decrease with the increase of R. Overall, for both of the two rain types, their relationships can possess acceptable accuracy when R exceeds 1 mm h -1 .