A Comparative Study on the Vertical Structures and Microphysical Properties of Stratiform Precipitation over South China and the Tibetan Plateau

Under different water vapor and dynamic conditions, and the influence of topographies and atmospheric environments, stratiform precipitation over South China and the Tibetan Plateau can produce different features. In this study, stratiform precipitation vertical characteristics, bright-band (BB) microstructures, and the vertical variations of the raindrop size distribution (DSD) over a low-altitude site (Longmen site, 86 m) in South China and a high-altitude site (Nagqu site, 4507 m) on the Tibetan Plateau were comprehensively investigated and compared using measurements from a Ka-band millimeter-wave cloud radar (CR), a K-band microrain radar (MRR), and a Parsivel disdrometer (disdrometer). A reliable BB identification scheme was proposed on the basis of CR variables and used for stratiform precipitation sample selection and further statistics and analysis. Results indicate that melting layers over the Longmen are much higher and slightly thicker than those over the Nagqu due to significant differences in atmospheric conditions. For stratiform precipitation, vertical air motions and radar variables over the two sites show different variation trends from cloud top to the ground. Vertical air motions are very weak in the stratiform precipitation over the Longmen, whereas updrafts are more active over the Nagqu. Above the melting layer, radar equivalent reflectivity factor Ze (mean Doppler velocity VM) gradually increases (decreases) as height decreases over the two sites, but the aggregation rate for ice particles over the Longmen can be faster. In the melting layer, Ze (VM) at the BB bottom/center over the Longmen is larger (smaller) than those over the Nagqu for the reason that melted raindrops in the melting layers over the Longmen are larger than those over the Nagqu. Below the melting layer, profiles of radar variables and DSDs show completely different behaviors over the two sites, which reflects that the collision, coalescence, evaporation, and breakup processes of raindrops are different between the two sites. Over the Longmen, collision and coalescence dominate the precipitation properties; in particular, from 2.0–2.8 km, the breakup process competes with collision–coalescence processes but later is overpowered. In contrast, due to the lower BB heights over the Nagqu, collision and coalescence dominate raindrop properties. Comparisons of raindrop spectra suggest that the concentration of small (medium-to-large) raindrops over the Nagqu is much higher (slightly lower) than that over the Longmen. Therefore, the mass-weighted mean diameter Dm (the generalized intercept parameter Nw) over the Nagqu is smaller (larger) than that over the Longmen.


Introduction
Stratiform precipitation is one of the most common rain types on the Earth and plays an important role in the atmospheric water cycle, radiation budget, and climate change [1][2][3][4][5].

Instruments and Measurements
The CR is a vertically pointing solid-state Doppler radar that operates at a frequency of 33.44 GHz, corresponding to a wavelength of 8.5 mm. It profiles the aloft clouds and precipitation at a height range from 0.12 to 15.3 km every ~9 s, and can simultaneously obtain the radar Doppler spectrum (DSP), equivalent reflectivity factor ( , dBZ), mean Doppler velocity ( , m·s −1 ), spectrum width ( , m·s −1 ), and linear depolarization ratio (LDR, dB) with a spatial resolution of 30 m. Liu et al. [30][31][32] assessed and verified that the CR measurements are reliable and capable of accurately capturing different types of nonprecipitating cloud and light precipitation.
The MRR is also a zenith radar but is designed using frequency-modulated continuous wave technology and operating at 24.23 GHz with a wavelength of 1.24 cm. It is manufactured by Meteorologische Messtechnik GmbH (METEK), Germany, and has been widely used for precipitation and DSD studies [33][34][35][36][37][38]. The temporal resolution of the MRR is set to 60 s. Considering the significant difference in the 0 °C isotherm height in the atmosphere between the two sites as a result of their different altitudes (the melting level over Nagqu site is relatively low, mostly near 1 km, whereas that over Longmen site is much higher, approximately 4 to 5 km), the maximum detectable heights of the MRR were set to 3.1 and 6.2 km with vertical resolutions of 100 and 200 m at the Nagqu and Longmen sites, respectively. The MRR measurements include attenuation-corrected DSP, , , the raindrop equivalent volume diameter ( , mm), and the number concentration per unit volume per millimeter ( ( ), m −3 ·mm −1 ), where i is the spectral bin number.
The ground disdrometer used is a second-generation Parsivel disdrometer (disdrometer), which has also been used in many previous studies [39][40][41][42][43][44]. The disdrometer is an optical-based disdrometer that transmits a 54-cm 2 (18-cm long, 3-cm wide) horizontal laser beam, and can simultaneously measure the diameter and falling velocity of passing raindrops according to the attenuation and duration of received signals. The raw raindrop counts were reserved and 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 . The sampling interval of the disdrometer is 60 s. For convenience, abbreviations of three instruments and physical quantities are used in this paper, of which the detailed information can be found in Appendix A (Table A1).

Data Processing and Quality Control
Signal attenuation of CR is mainly affected by rain rate and hydrometeors distribution in rain measurements. Signal will be seriously attenuated at a high rain rate. Besides, excessive raindrop size can cause millimeter-wave Mie scattering, and the returned signal

Instruments and Measurements
The CR is a vertically pointing solid-state Doppler radar that operates at a frequency of 33.44 GHz, corresponding to a wavelength of 8.5 mm. It profiles the aloft clouds and precipitation at a height range from 0.12 to 15.3 km every~9 s, and can simultaneously obtain the radar Doppler spectrum (DSP), equivalent reflectivity factor (Z e , dBZ), mean Doppler velocity (V M , m·s −1 ), spectrum width (σ V , m·s −1 ), and linear depolarization ratio (LDR, dB) with a spatial resolution of 30 m. Liu et al. [30][31][32] assessed and verified that the CR measurements are reliable and capable of accurately capturing different types of nonprecipitating cloud and light precipitation.
The MRR is also a zenith radar but is designed using frequency-modulated continuous wave technology and operating at 24.23 GHz with a wavelength of 1.24 cm. It is manufactured by Meteorologische Messtechnik GmbH (METEK), Germany, and has been widely used for precipitation and DSD studies [33][34][35][36][37][38]. The temporal resolution of the MRR is set to 60 s. Considering the significant difference in the 0 • C isotherm height in the atmosphere between the two sites as a result of their different altitudes (the melting level over Nagqu site is relatively low, mostly near 1 km, whereas that over Longmen site is much higher, approximately 4 to 5 km), the maximum detectable heights of the MRR were set to 3.1 and 6.2 km with vertical resolutions of 100 and 200 m at the Nagqu and Longmen sites, respectively. The MRR measurements include attenuation-corrected DSP, Z e , V M , the raindrop equivalent volume diameter (D i , mm), and the number concentration per unit volume per millimeter (N(D i ), m −3 ·mm −1 ), where i is the spectral bin number.
The ground disdrometer used is a second-generation Parsivel disdrometer (disdrometer), which has also been used in many previous studies [39][40][41][42][43][44]. The disdrometer is an optical-based disdrometer that transmits a 54-cm 2 (18-cm long, 3-cm wide) horizontal laser beam, and can simultaneously measure the diameter and falling velocity of passing raindrops according to the attenuation and duration of received signals. The raw raindrop counts were reserved and 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 . The sampling interval of the disdrometer is 60 s. For convenience, abbreviations of three instruments and physical quantities are used in this paper, of which the detailed information can be found in Appendix A (Table A1).

Data Processing and Quality Control
Signal attenuation of CR is mainly affected by rain rate and hydrometeors distribution in rain measurements. Signal will be seriously attenuated at a high rain rate. Besides, excessive raindrop size can cause millimeter-wave Mie scattering, and the returned signal can be weakened [45][46][47]. Therefore, to alleviate the attenuation of the radar's DSP and Z e , an iterative correction procedure was implemented [48,49], as shown in Equations (1)-(4), where i and n are the range gate number and spectral bin number, respectively; k i (dB·km −1 ) and τ i are the correction coefficient and two-way transmissivity, respectively; Z em and Z ec denote the radar-measured and corrected equivalent reflectivity factor, respectively; DSP m and DSP c represent the radar-measured and corrected Doppler spectra, respectively; and ∆R is the radar gate length. Coefficients α and β were set to 0.00334 and 0.73, respectively [50]. MRR's high transmission frequency increases sensitivity in detecting small droplets; however, it also attenuates severely under heavy rain. The rain attenuation correction was applied for the path-integrated rain attenuation less than ten according to a recursive algorithm developed by Kunz [51]. Therefore, the attenuation can be basically ignored for the precipitation of small to medium rain rates. In addition, Peter et al. [34] found a surprisingly good agreement between radar reflectivities of an MRR and a C-band weather radar.
A disdrometer can produce some unrealistic data due to its inherent limitations [40,52]. Therefore, the following data quality control technologies were adopted. First, the data in the first two diameter classes were eliminated because of their low signal-to-noise ratio, and any raindrops with a diameter > 8 mm were discarded because of their natural breakup in the atmosphere [52]. Second, any DSDs with a total raindrop number less than 10 or a rain rate smaller than 0.002 mm·h −1 were treated as nonrainy samples and removed [27]. Third, any raindrops possessing an abnormal diameter (falling velocity) were removed by comparing the measured result with the theoretical relationship proposed by Atlas et al. [53]. If the measured diameter (falling velocity) of raindrops was 60% of the theoretical value, the raindrops were regarded as problematic data and abandoned. On the basis of the aforementioned processing, DSD was further calculated as where i and j are the raindrop diameter and velocity class orders, respectively; N(D i ) (m −3 ·mm −1 ) and n ij represent the raindrop number concentration and the raindrop count, respectively; ∆D i (mm) refers to the diameter interval; V j is the raindrop falling velocity; A i (m 2 ) is the effective sampling area; and ∆t (60 s) denotes the sampling time. The rain rate (R, mm·h −1 ), mass-weighted mean diameter (D m , mm), and generalized intercept parameter (N w , m −3 ·mm −1 ) were calculated using Equations (6)- (8).

. Stratiform Precipitation Sample Identification
Many studies have reported that the BB signature can be detected more easily in stratiform precipitation, as observed using centimeter-or millimeter-wavelength radars [54][55][56][57]. Therefore, in the present study, a BB identification scheme was implemented to distinguish stratiform precipitation samples from the entire radar dataset. Considering that the CR has a much higher spatiotemporal resolution than the MRR and the LDR-a polarized variable that can make the recognition of BB more robust-the scheme was designed on the basis of the CR measurements. A brief flowchart is set out in Figure 2, which illustrates the scheme's five steps.
Remote Sens. 2021, 13, x FOR PEER REVIEW To exhibit the BB identification effect of the proposed scheme, Figure 3 pres stratiform precipitation cases observed at the Longmen (left panels) and the Nagq panels), respectively. Specifically, Figure 3a1,a2 give the CR's s; Figure 3b1,b2 s CR's LDRs; Figure 3c1,c2 show the MRR's s; Figure 3d1,d2 present the disdr observed rain rate. The identified BB center, top, and bottom are represented solid lines and superimposed on each radar image.
It can be seen that both BBs of the two cases are well identified under the stance of moderate rainfall on the ground, with rain rates ranging from 0.1 to ~7 BB heights over the Longmen and Nagqu sites are located at approximately 4.5 an respectively. BB signatures in terms of the rapid increase in along the meltin are observed using two radars. Moreover, for the CR, the LDR was more sensitiv particle melting process and can be a better BB indicator than . Comparisons of ages from the two radars also illustrate that the CR can measure the entire per and identify fine stratiform precipitation structures owing to its higher spatial re Non-BB

Step 3. Preliminary determination of BB bottom and top
Step 2. Preliminary determination of BB center Step 1. Exclusion of non-BB profiles Step 4. Verification of BB heights Step  Step 1. Exclusion of non-BB profiles. According to BB observational features, both the radar's Z e and V M can significantly increase when hydrometeors pass through the melting layer with phase changes from ice to liquid [58,59]. Therefore, if all range gates possessed Z e < −10 dBZ or V M > −1.0 m·s −1 , a radar profile was treated as a non-BB sample and discarded [59,60].
Step 2. Preliminary determination of the BB center. A range interval of 750 m (25 consecutive range gates) was set with a threshold representing a possible maximum BB thickness observed using the CR [61,62]. For each reserved profile after Step 1, search maximums for radar reflectivities (ZCs) and relevant heights (HZCs) within each range interval from bottom to top. Each HZC was then set as the center to find two minimum reflectivities (ZU, ZL) and corresponding heights (HZU, HZL) within the upper and lower parts of each range interval, respectively. When the criterion of (ZC − ZU) × (ZC − ZL) > 18 dBZ was met, the HZC was considered to be a potential BB center height. If multiple HZCs satisfied the criterion at the same time, only the one with the largest ZC was selected. Using the same approach, it was also possible to obtain a BB center height (HLdrC) from the radar LDR profile but with the constraints of LdrC > LdrU, LdrC > LdrL, and (LdrC − LdrU) × (LdrC − LdrL) > 20 dB [63]. Finally, if the HZC and HLdrC met the consistency test, an average for HZC and HLdrC was produced as the final BB height (HC final ). The consistency test was proposed by Brandes et al. [64] and reflected the consistency between the variations of reflectivity and polarization parameter. If |HZC -HLdrC < d|, it means that the two parameters simultaneously detected the properties change of particles in the melting layer. The d (km) was calculated using Equation (9).
Step 3. Preliminary determination of the BB bottom and top. According to vertical variations in Z e , V M , σ V , and LDR near the melting layer, four thresholds were set: ∆Z e = 3 dBZ, ∆V M = 0.6 m·s −1 , ∆σ V = 0.4 m·s −1 , and ∆LDR = 3 dB [59]. For each profile with an identified BB center, gradients for the radar observations were established upward and downward from the BB center; when at least three radar variables consecutively increased (decreased) or exceeded the thresholds within 150 m, their corresponding heights HZT/HZB, HVT/HVB, HSwT/HSwB, and HLdrT/HLdrB were calculated (where "H" is the height; "T" and "B" denote the top and bottom, respectively; and "V" and "Sw" are V M and σ V , respectively). Comparing these top (bottom) heights, if their differences were within 200 m [63], their average was determined as the final BB top (bottom) height HT final (HB final ).
Step 4. Verification of BB heights. Further comparing the difference between HC final and HT final /HB final , if HT final − HC final < 750 m and HB final − HC final < 750 m, the identification results were perceived to be robust.
Step 5. Calculation of the BB thickness and radar variables. Following the aforementioned steps, the BB thickness and the radar's Z e , V M , σ V , and LDR at the heights of the BB center, top, and bottom were calculated and reserved for further statistical analysis.
To exhibit the BB identification effect of the proposed scheme, Figure 3 presents two stratiform precipitation cases observed at the Longmen (left panels) and the Nagqu (right panels), respectively. Specifically, Figure 3a1 It can be seen that both BBs of the two cases are well identified under the circumstance of moderate rainfall on the ground, with rain rates ranging from 0.1 to~7 mm·h −1 . BB heights over the Longmen and Nagqu sites are located at approximately 4.5 and 1 km, respectively. BB signatures in terms of the rapid increase in Z e along the melting layers are observed using two radars. Moreover, for the CR, the LDR was more sensitive to the particle melting process and can be a better BB indicator than Z e . Comparisons of Z e images from the two radars also illustrate that the CR can measure the entire perspective and identify fine stratiform precipitation structures owing to its higher spatial resolution and sensitivity, whereas the Z e was attenuated by rainfall to some extent. In contrast, the MRR has much shorter detectable height ranges and rougher spatial resolution, but its signal attenuation was insignificant. Overall, the identified heights for the bottom, center, and top of the BB based on the CR measurements were consistent with the melting layer structure reflected by the MRR. Remote Sens. 2021, 13, x FOR PEER REVIEW 8 of 21

Data Matching
After BB identification using the CR, all distinguished stratiform precipitation samples were further verified using disdrometer-measured rain rate. Stratiform precipitation samples corresponding to ground rain rates being smaller than 5 mm·h −1 were regarded as the reliable results and reserved because the CR is inappropriate for detecting mediumto-heavy precipitation. Then, samples with ground rain rates smaller than 0.02 mm·h −1 were discarded as noise caused by nonprecipitation echo. After checking the continuity of samples in the time series, isolated or intermittent samples with a duration of less than 10 min were discarded to avoid some spurious and erratic results [13]. Finally, all corresponding minute data of the MRR and disdrometer were selected and combined with the CR results to form the final matched stratiform precipitation datasets. Statistics show that 41 and 47 rain events over the Longmen and Nagqu sites were collected in total. The matched stratiform precipitation datasets for the three instruments as well as the corresponding ground mean rain rate and accumulated rain amounts from the observation periods over the two sites are listed in Table 1.

Data Matching
After BB identification using the CR, all distinguished stratiform precipitation samples were further verified using disdrometer-measured rain rate. Stratiform precipitation samples corresponding to ground rain rates being smaller than 5 mm·h −1 were regarded as the reliable results and reserved because the CR is inappropriate for detecting mediumto-heavy precipitation. Then, samples with ground rain rates smaller than 0.02 mm·h −1 were discarded as noise caused by nonprecipitation echo. After checking the continuity of samples in the time series, isolated or intermittent samples with a duration of less than 10 min were discarded to avoid some spurious and erratic results [13]. Finally, all corresponding minute data of the MRR and disdrometer were selected and combined with the CR results to form the final matched stratiform precipitation datasets. Statistics show that 41 and 47 rain events over the Longmen and Nagqu sites were collected in total. The matched stratiform precipitation datasets for the three instruments as well as the corresponding ground mean rain rate and accumulated rain amounts from the observation periods over the two sites are listed in Table 1.

Results
Under different water vapor and dynamic conditions and the influence of altitudes and atmospheric environments, the stratiform precipitation over South China and the Tibetan Plateau can produce different features. Therefore, on the basis of the matched datasets, stratiform precipitation vertical characteristics, BB microstructures, and DSD vertical variations over the Longmen and Nagqu sites were analyzed and compared in detail. Figure 4 presents the averaged profiles of radar variables, including the CR's Z e , V M , σ V , LDR, and DSP and the MRR's Z e and V M , for the entire set of stratiform precipitation samples over the two sites. Note that, herein, for the V M , a positive velocity represents upward motion of the targets (away to radar antenna), and negative velocity means downward motion of targets (toward from radar antenna). However, for the DSP, the signs are opposite. As seen from Figure 4a1-a4,b1-b4, the profiles of radar variables observed over the two sites exhibit different stratiform precipitation vertical structures. Specifically, melting layers with obvious BB signatures over the Longmen and Nagqu sites are located at height ranges of approximately 4.4-5.1 km and 0.8-1.3 km, respectively. For cloud precipitation above the melting layers, Z e s (V M s) over the two sites increase (decrease) as height decreases because of ice particle aggregation. To further quantitatively denote the changes of Z e and V M for cloud ice hydrometeors, average curves of Z e and V M above the BB were further fitted by a linear equation, as printed on Figure 4a1,a2,b1,b2. Fitting results suggest that the Z e profile over the Longmen site has a larger slope than that over the Nagqu, implying a faster aggregation rate of ice particles. Interestingly, at the BB top, on average, the Longmen possesses a larger Z e . but a similar V M compared to counterparts of the Nagqu. This phenomenon can be perceived as the ice hydrometeor over the Longmen is larger; however, its density can be lower. As ice particles fall into melting layers, s/LDRs first increase and then decrease and V M s (σ V s) decrease (increase) as height decreases over the two sites. In addition, the increments (decrements) in s/σ V s/LDRs (V M s) over the Longmen are larger than those over the Nagqu, implying that the stratiform precipitation of the former forms larger raindrops at the bottom of the melting layer than that of the latter. After melting, raindrops over the Longmen have a much longer falling path before touching the ground than those over the Nagqu, with the former's Z e s/s (σ V s) gradually increasing (decreasing) again as height decreases, while the latter's Z e s have no apparent changes, and its V M s (σ V s) continue to decrease (increase) as height decrease. Nevertheless, the median profiles can better show the variations of measurements at BB heights, and the averaged profile is easily affected by the extreme values. In addition, radar observations are inevitably affected by clutter near the ground; for the profiles below BB in the Nagqu site, the accuracy may need to be further confirmed. Figure 4a5,b5 illustrate that the averaged CR spectra over the two sites clearly indicate the melting process with sharp shifts in Doppler signals. Liu et al. [65] pointed out that the influence of turbulence on DSP is negligible in the case of weak precipitation. Therefore, the negative velocities of the DSP can approximately represent the updrafts. The hydrometeor signals over the Longmen possess positive velocities, suggesting that the vertical air motions are very weak in cloud precipitation. In contrast, there exist partial hydrometeor signals over the Nagqu that possess negative velocities, revealing that the updrafts are more active. Comparing the two radars' profiles, it can be seen that the CR's Z e s are attenuated by 2.5-7.0 and 2-4.5 dB over the Longmen and Nagqu, respectively. Furthermore, attenuations above the BB are more significant than those below the BB because the signal emitted by radar is more attenuated by rain rate, temperature, liquid water path, etc. as the distance increases. . Averaged profiles of radar variables for the stratiform precipitation over the Longmen and Nagqu sites, respectively. (a1-a5) show the radar equivalent reflectivity factor ( , dBZ), mean Doppler velocity ( , m·s −1 ), spectrum width ( , m·s −1 ), linear depolarization ratio (LDR, dB), and Doppler spectra (DSP) observed over the Longmen. (b1-b5) show the same but observed over the Nagqu. The functions in (a1,a2,b1,b2) represent the fitting results of / profiles above the BB. Figure 5 shows the probability density functions (PDFs) of the heights of the BB ce ter, top, and bottom as well as the BB entire, upper-half, and lower-half thicknesses f stratiform precipitation over the two sites. Corresponding statistical parameters involvi 5th-95th percentiles, averages (AVs), standard deviations (STDs), and skewnesses (SK for these BB heights and thicknesses are listed in Table 2. Comparisons indicate that t BB heights over the Longmen are 3 to 4 km higher than those over the Nagqu. Most of t BB bottom/center/top heights (5th-95th) over the two sites are concentrated within 4. negatively skewed over the Nagqu. With regard to BB thicknesses, the entire thickness of the BB over the Longmen are slightly larger than those over the Nagqu, most of the over the two sites are within 0.33-0.69 km and 0.27-0.63 km, with AVs of 0.51 and 0.44 k respectively. Previous studies perceived that the BB thickness can be affected by the me ing rate and the fall velocity of particles [16,58]. As shown in Figure 4, at the BB top, t ice particles over the Nagqu are smaller than those over the Longmen; thus, the melti process at the former site also can be faster and the BB is thinner. The BB entire thickness in these two regions are consistent with the results observed by other scholars, who studies show that the entire thickness is generally not more than 0.75 km when measur using the CR [59,61]. Moreover, the lower-and upper-half thicknesses over the Longm possess similar statistical parameters except in terms of SKs, demonstrating that the BB near-symmetrical. However, the upper-half thicknesses are slightly larger than those the lower-half thicknesses over the Nagqu. The BB thicknesses possess near-zero STD over the two sites, suggesting they are near-stable.   Figure 5 shows the probability density functions (PDFs) of the heights of the BB center, top, and bottom as well as the BB entire, upper-half, and lower-half thicknesses for stratiform precipitation over the two sites. Corresponding statistical parameters involving 5th-95th percentiles, averages (AVs), standard deviations (STDs), and skewnesses (SKs) for these BB heights and thicknesses are listed in Table 2. Comparisons indicate that the BB heights over the Longmen are 3 to 4 km higher than those over the Nagqu. The BB heights are positively skewed over the Longmen, but negatively skewed over the Nagqu. With regard to BB thicknesses, the entire thicknesses of the BB over the Longmen are slightly larger than those over the Nagqu, most of these over the two sites are within 0.33-0.69 km and 0.27-0.63 km, with AVs of 0.51 and 0.44 km, respectively. Previous studies perceived that the BB thickness can be affected by the melting rate and the fall velocity of particles [16,58]. As shown in Figure 4, at the BB top, the ice particles over the Nagqu are smaller than those over the Longmen; thus, the melting process at the former site also can be faster and the BB is thinner. The BB entire thicknesses in these two regions are consistent with the results observed by other scholars, whose studies show that the entire thickness is generally not more than 0.75 km when measured using the CR [59,61]. Moreover, the lower-and upper-half thicknesses over the Longmen possess similar statistical parameters except in terms of SKs, demonstrating that the BB is near-symmetrical. However, the upper-half thicknesses are slightly larger than those of the lower-half thicknesses over the Nagqu. The BB thicknesses possess near-zero STDs over the two sites, suggesting they are near-stable. Remote Sens. 2021, 13, x FOR PEER REVIEW 11 of 21  To further ascertain discrepancies in radar variables at different BB positions over the two sites, Figure 6 shows PDFs of radar , , , and LDR measured at the BB bottom, center, and top. Corresponding statistics of medians (MDs), AVs, STDs, and SKs and discrepancies between the center and the top (the BB upper part ∆ , ∆ = − ) and between the center and the bottom (the BB lower part ∆ , ∆ = − ) of these radar variables are also calculated and summarized in Table 3. It can be seen that as ice particles pass through the top, center, and bottom of the BB, radar variables change significantly as a result of the melting process. However, variable values at each BB position and deviations in the same part of the BB over the two sites are different. On average, s ( s) at the BB bottom/center/top over the Longmen are 1.2/0.8/2.1 dB (1.6/1.0/0.5 m·s −1 ) larger (smaller) than those over the Nagqu, indicating that the ice hydrometeors that formed above the BB and corresponding melted raindrops below the BB over the former site are larger than those over the latter site. LDRs at the three  To further ascertain discrepancies in radar variables at different BB positions over the two sites, Figure 6 shows PDFs of radar Z e , V M , σ V , and LDR measured at the BB bottom, center, and top. Corresponding statistics of medians (MDs), AVs, STDs, and SKs and discrepancies between the center and the top (the BB upper part ∆up, ∆up = top − center) and between the center and the bottom (the BB lower part ∆down, ∆down = center − down) of these radar variables are also calculated and summarized in Table 3.

Stratiform Precipitation Vertical Characteristics and BB Microstructures
is found that their PDFs shapes are similar.
As listed in Table 3, the ∆ s and ∆ s of the radar variables over both sites exhibit as |∆ | > |∆ | , |∆ | < |∆ | , |∆ | < |∆ | , and |∆ | < |∆ |, suggesting that is more sensitive to the melted ice hydrometeors attached with water film. The changes of the dielectric constant and the collision and coalescence of hydrometeors are the reasons for the significant increase of in the BB upper-half. Comparisons also illustrate that both |∆ | and |∆ | (|∆ | and |∆ |) over the Longmen are smaller (larger) than those over the Nagqu.     It can be seen that as ice particles pass through the top, center, and bottom of the BB, radar variables change significantly as a result of the melting process. However, variable values at each BB position and deviations in the same part of the BB over the two sites are different. On average, Z e s (V M s) at the BB bottom/center/top over the Longmen are 1.2/0.8/2.1 dB (1.6/1.0/0.5 m·s −1 ) larger (smaller) than those over the Nagqu, indicating that the ice hydrometeors that formed above the BB and corresponding melted raindrops below the BB over the former site are larger than those over the latter site. LDRs at the three BB positions over the Longmen are also 2.6/3.0/3.8 dB larger than those over the Nagqu, which may be due to the lower condensation height and some supercooled droplets contained in the melting layer over the Longmen [66]. Furthermore, σ V s at the bottom and center of BB over the Longmen are slightly larger than those over the Nagqu, whereas at the top of the BB, they are somewhat smaller. Comparing the Z e s/V M s of CR and MRR, it is found that their PDFs shapes are similar.
As listed in Table 3, the ∆ups and ∆downs of the radar variables over both sites exhibit as |∆upZ e | > |∆downZ e |, |∆upV M | < |∆downV M |, |∆upσ V | < |∆downσ V |, and |∆upLDR| < |∆downLDR|, suggesting that Z e is more sensitive to the melted ice hydrometeors attached with water film. The changes of the dielectric constant and the collision and coalescence of hydrometeors are the reasons for the significant increase of Z e in the BB upper-half. Comparisons also illustrate that both |∆upZ e | and |∆downZ e | (|∆upV M | and |∆downV M |) over the Longmen are smaller (larger) than those over the Nagqu.

Vertical Variations of Stratiform Precipitation Raindrop Spectra
Owing to significant differences in melting layer heights between the low-and highaltitude sites in the summer, falling raindrops can undergo diverse physical processes and possess different microphysical properties. Therefore, herein, stratiform precipitation DSDs observed using the MRR over the two sites were further compared. Figures 7 and 8 show average DSDs versus heights and contoured frequency by altitude diagrams (CFADs) of the D m and the logarithm of N w (log 10 N w ) below melting layers for all stratiform precipitation samples, respectively. In contrast, over the Nagqu, the vertical variations of DSDs are not as significant as over the Longmen. From the bottom of the BB to the ground, the concentration of small raindrops (D < 1.0 mm) sharply decreases, the concentration of medium raindrops (1.0 < D < 2.0 mm) remains near constant, and the concentration of large raindrops (D > 2.0 mm) increases. Especially, the larger the diameter, the greater the increase in concentration. It indicates that raindrops collision and coalescence processes are dominant. Figure 8b,d show that the tends to decrease gradually, whereas the increases from approximately 0.5-0.7 mm at BB bottom to 0.6-0.8 mm at 0.3 km. Moreover, all diameter bins have a decreasing trend near the ground, and the increment (decrement) of ( ) slightly increases, which is possibly attributable to the lower humidity and stronger evaporation at a low level over the Tibetan Plateau [29,68].
DSDs over the two sites clearly show that the concentration of raindrops with D < 1.4 mm over the Nagqu is much higher than that over the Longmen, with concentrations over the two sites ranging from 10 1 to 10 4.5 m −3 ·mm −1 and from 10 1 to 10 3 m −3 ·mm −1 , respectively. This finding is similar to the conclusion of Sumesh et al. [16] that small drops have a significant concentration in the high-altitude regions. On the contrary, amounts of medium-to-large raindrops (D > 2 mm) over the Longmen are greater than those over the Nagqu. Moreover, the over the Longmen is larger than that over the Nagqu, as found in the previous study [29]. However, the shows the opposite behavior, with the over the Longmen being much smaller than that over the Nagqu.

Discussion
In this section, we compare the DSD results of MRR with disdrometer. Besides, our results on microphysical characteristics of stratiform precipitation are compared with those obtained in other areas in China.

Comparison of MRR with Disdrometer
To examine the reliability of the MRR data, Figure 9 sets out the average DSDs measured using the MRR and disdrometer. Note that according to the actual spatial resolutions of the MRR deployed at the two sites, the lowest four range gates (800, 600, 400, and 200 m over the Longmen and 400, 300, 200, and 100 m over the Nagqu) were selected for comparison.
Compared with those of the disdrometer, the MRR's concentrations of small raindrops (D ≤ 0.4 mm over the Longmen and D ≤ 0.5 mm over the Nagqu) are much higher, which can be attributed to the inherent limitation of the disdrometer in underestimating small raindrops [40,69]. For raindrops within 0.4-4 mm over the Longmen and raindrops within 0.5-2.5 mm over the Nagqu, the radar results are consistent with the disdrometer results, but concentrations are underestimated to some extent. Conversely, for raindrops larger than 4 mm over the Longmen and raindrops larger than 2.5 mm over the Nagqu, radar can slightly overestimate concentrations. These discrepancies between the MRR and disdrometer have also been reported by Konwar et al., Jash et al., and Wen et al. [9,13,70]. The mismatch between observations of two instruments is because of the laser inhomogeneity problem of disdrometer, different measuring principles of disdrometer and MRR, different sampling volumes and heights, and the presence of updraftdowndraft or turbulence in the atmosphere. As can be seen, DSDs show clear but different vertical variations over the two sites. Over the Longmen, from the bottom of the BB up to 2.0 km, the concentration of small raindrops (D < 0.5 mm) gradually decreases, whereas the concentration of medium-tolarge raindrops (D > 1.6 mm) increases as height decreases, suggesting that collision and coalescence processes dominate at this height range. From 2.0 to 2.8 km, the concentration of raindrops (D < 4 mm) increases, especially for small raindrops (D < 0.6 mm). From the averaged DSDs, the concentration of large raindrops (D > 4 mm) decreases (averaged DSDs from 2.0 km to 2.8 km are not shown). It indicates that the breakup process of the large raindrops and the collision-coalescence processes of medium raindrops coexist. From 2.8 to 3.4 km, the concentration of small raindrops (D < 1.0 mm) decreases, whereas the concentration of larger raindrops (D > 1.0 mm) increases, suggesting that collision and coalescence processes dominate at this height range. Below 3.4 km, the concentration of small raindrops (D < 0.6 mm) decreases, whereas the concentration of medium-to-large raindrops (D > 1.0 mm) quickly increases, suggesting that collision and coalescence processes dominate at this height range. These results are consistent with previous studies [19,37,67], which concluded that significant growth in the collision and coalescence processes occurs only when small raindrops (D < 0.6 mm) are collected by large raindrops. As for CFADs of D m and log 10 N w over the Longmen, Figure 8a,c show the increased D m and decreased log 10 N w when raindrops fall from the bottom of the BB. In particular, from 2.0 to 2.8 km, D m (log 10 N w ) slightly increases (decrease), which represents the collision and coalescence processes overwhelming the breakup process. Near the ground, log 10 N w increases suddenly from~2.5 at 3.4 km to~2.8 at ground level, which is associated with the suddenly increased concentration of medium raindrops.
In contrast, over the Nagqu, the vertical variations of DSDs are not as significant as over the Longmen. From the bottom of the BB to the ground, the concentration of small raindrops (D < 1.0 mm) sharply decreases, the concentration of medium raindrops (1.0 < D < 2.0 mm) remains near constant, and the concentration of large raindrops (D > 2.0 mm) increases. Especially, the larger the diameter, the greater the increase in concentration. It indicates that raindrops collision and coalescence processes are dominant. Figure 8b,d show that the log 10 N w tends to decrease gradually, whereas the D m increases from approximately 0.5-0.7 mm at BB bottom to 0.6-0.8 mm at 0.3 km. Moreover, all diameter bins have a decreasing trend near the ground, and the increment (decrement) of D m (N w ) slightly increases, which is possibly attributable to the lower humidity and stronger evaporation at a low level over the Tibetan Plateau [29,68].
DSDs over the two sites clearly show that the concentration of raindrops with D < 1.4 mm over the Nagqu is much higher than that over the Longmen, with concentrations over the two sites ranging from 10 1 to 10 4.5 m −3 ·mm −1 and from 10 1 to 10 3 m −3 ·mm −1 , respectively. This finding is similar to the conclusion of Sumesh et al. [16] that small drops have a significant concentration in the high-altitude regions. On the contrary, amounts of medium-to-large raindrops (D > 2 mm) over the Longmen are greater than those over the Nagqu. Moreover, the D m over the Longmen is larger than that over the Nagqu, as found in the previous study [29]. However, the log 10 N w shows the opposite behavior, with the log 10 N w over the Longmen being much smaller than that over the Nagqu.

Discussion
In this section, we compare the DSD results of MRR with disdrometer. Besides, our results on microphysical characteristics of stratiform precipitation are compared with those obtained in other areas in China.

Comparison of MRR with Disdrometer
To examine the reliability of the MRR data, Figure 9 sets out the average DSDs measured using the MRR and disdrometer. Note that according to the actual spatial resolutions of the MRR deployed at the two sites, the lowest four range gates (800, 600, 400, and 200 m over the Longmen and 400, 300, 200, and 100 m over the Nagqu) were selected for comparison.   [70] and from the monsoon frontal rainfa Central China obtained by Zhou et al. [19], respectively.

Comparison of Stratiform Precipitation Microphysical Characteristics
Comparisons clearly show that the Longmen generally has larger s but s than the Nagqu at each height, and s over the low-and high-altitud  Compared with those of the disdrometer, the MRR's concentrations of small raindrops (D ≤ 0.4 mm over the Longmen and D ≤ 0.5 mm over the Nagqu) are much higher, which can be attributed to the inherent limitation of the disdrometer in underestimating small raindrops [40,69]. For raindrops within 0.4-4 mm over the Longmen and raindrops within 0.5-2.5 mm over the Nagqu, the radar results are consistent with the disdrometer results, but concentrations are underestimated to some extent. Conversely, for raindrops larger than 4 mm over the Longmen and raindrops larger than 2.5 mm over the Nagqu, radar can slightly overestimate concentrations. These discrepancies between the MRR and disdrometer have also been reported by Konwar et al.,Jash et al.,and Wen et al. [9,13,70]. The mismatch between observations of two instruments is because of the laser inhomogeneity problem of disdrometer, different measuring principles of disdrometer and MRR, different sampling volumes and heights, and the presence of updraft-downdraft or turbulence in the atmosphere. Figure 10 shows scatterplots of log 10 N w -D m pairs at different heights observed from the MRR and disdrometer. The average values of log 10 N w and D m are marked on each subgraph, as well as the two previous results for summer precipitation in other regions of China. The cyan pluses and green crosses represent the results from the Meiyu season over East China reported by Wen et al. [70] and from the monsoon frontal rainfalls over Central China obtained by Zhou et al. [19], respectively.

Comparison of Stratiform Precipitation Microphysical Characteristics
Remote Sens. 2021, 13, x FOR PEER REVIEW are attributed to the overpowered collision and coalescence processes over th process [71]. The vertical structures over South China have some similarities w over Central China over East China. When raindrops fall, collision and coales cesses dominate, but there is competition between breakup and collision-coales cesses from 2.0-2.8 km. Over the Tibetan Plateau, the collision and coalescence dominate raindrop properties.
Overall, the microphysical characteristics of stratiform precipitation are a many factors, such as weather systems, water vapor conditions, topography, et lision, coalescence, and evaporation processes govern the DSDs' vertical and evolutions, which make precipitation properties quite diverse.

Conclusions
Using measurements from a CR, an MRR, and a disdrometer, stratiform pre samples were selected using a reliable BB identification scheme. After data stratiform precipitation vertical characteristics, BB microstructures, and DSD ve iations over a low-altitude site (Longmen) in South China and a high-altitude sit on the Tibetan Plateau were comprehensively investigated and compared. The clusions of this study can be summarized as follows.
The Longmen site is located at low altitude, which leads to a much highe that over the Nagqu. Due to the obvious distinctions in the atmospheric environ the two sites, radar variables and vertical air motions in stratiform precipitation Comparisons clearly show that the Longmen generally has larger D m s but smaller log 10 N w s than the Nagqu at each height, and D m s over the low-and high-altitude sites mainly distribute within 1.2-1.4 and 0.7-1.0 mm, corresponding to log 10 N w s within 2.6-3.2 and 3.8-4.6 m −3 ·mm −1 , respectively. In the study of Zhou et al. [19], at 600 and 400 m, the average log 10 N w -D m pairs for Central China are located between those over our two sites, with smaller D m s than the Longmen and larger than the Nagqu; log 10 N w s are larger than those over the Longmen but smaller than those over the Nagqu. At 200 m and at ground level, log 10 N w s are larger than those over the Longmen and Nagqu sites. For results obtained by Wen et al. [70], the log 10 N w . for East China is slightly smaller than that of the Nagqu and larger than that of the Longmen at 200 m. On the ground, log 10 N w s are larger than that over our two sites. As a result of the high liquid water content in the atmosphere over South China, the collision and coalescence processes make D m s larger than those over the other aforementioned regions and, correspondingly, log 10 N w s can be smaller. Although summer precipitation in the three regions is influenced by the East Asian monsoon, the statistical results of SP are also different. Altitude, water vapor, atmospheric conditions, etc. may be the main factors that explain the difference, although diverse measuring instruments and quality control methods also play a part. However, the Tibetan Plateau is very different from the other three regions. Due to the short falling path of raindrops, collision and coalescence processes are insufficient, leading to a relatively higher concentration of small raindrops reaching the ground compared with the low altitude.
Microphysical processes of stratiform precipitation below the melting layer over three regions (North China, Central China, and East China) in China are compared with those over South China and the Tibetan Plateau. Over North China, both large and small raindrops have high concentrations at the bottom of BB. With decreasing height, the coalescence process of small raindrops and the breakup of large raindrops reach a balance [8]. The microphysical processes over Central China are different. Below BB, the collision and coalescence processes slightly overpower the breakup process, resulting in slightly decreasing Z e with decreasing height [19]. Over East China, Z e and liquid water content stay quasi-constant below BB. D m first decreases and then increases, while N w is the opposite. It indicates that large raindrops formed by the melting of ice particles above the zero-degree layer breakup into small raindrops. As height decreases, raindrop properties are attributed to the overpowered collision and coalescence processes over the breakup process [71]. The vertical structures over South China have some similarities with those over Central China over East China. When raindrops fall, collision and coalescence processes dominate, but there is competition between breakup and collision-coalescence processes from 2.0-2.8 km. Over the Tibetan Plateau, the collision and coalescence processes dominate raindrop properties.
Overall, the microphysical characteristics of stratiform precipitation are affected by many factors, such as weather systems, water vapor conditions, topography, etc. The collision, coalescence, and evaporation processes govern the DSDs' vertical and temporal evolutions, which make precipitation properties quite diverse.

Conclusions
Using measurements from a CR, an MRR, and a disdrometer, stratiform precipitation samples were selected using a reliable BB identification scheme. After data matching, stratiform precipitation vertical characteristics, BB microstructures, and DSD vertical variations over a low-altitude site (Longmen) in South China and a high-altitude site (Nagqu) on the Tibetan Plateau were comprehensively investigated and compared. The main conclusions of this study can be summarized as follows.
The Longmen site is located at low altitude, which leads to a much higher BB than that over the Nagqu. Due to the obvious distinctions in the atmospheric environment over the two sites, radar variables and vertical air motions in stratiform precipitation show diverse variation trends from cloud top to the ground. Above the melting layer, Z e s (V M s) gradually increase (decrease) as height decreases over the two sites, but the aggregation rate for ice particles can be faster over the Longmen site. At the BB top, larger ice particles with smaller densities have been formed over the Longmen because of similar V M but larger Z e . In the melting layer, Z e s (V M s) at the bottom/center of the BB over the Longmen are larger (smaller) than those over the Nagqu due to the larger ice hydrometeors melted into larger raindrops. Therefore, over the Longmen, hydrometeors melt and fall more slowly in the melting layer, forming a thicker BB thickness. Moreover, the radar variables' deviations in the same part of the BB over the two sites are also different. The ∆ups and ∆downs of the radar variables over both sites exhibit as |∆upZ e | > |∆downZ e |, |∆upV M | < |∆downV M |, |∆upσ V | < |∆downσ V |, and |∆upLDR| < |∆downLDR|, suggesting that Z e is more sensitive to the melting ice hydrometeors attached with water film. Below the bottom of the BB, radar variables have completely different behaviors. Z e s/V M s (σ V s) gradually increase (decrease) over the Longmen as height decreases, whereas over the Nagqu, there are no apparent change in Z e s, and its V M s (σ V s) continue to decrease (increase) as height decreases, which reflects that the collision, coalescence, evaporation, and breakup processes of raindrops are different between two sites. In addition, vertical air motions are very weak in the stratiform precipitation over the Longmen, whereas updrafts are more active over the Nagqu.
Falling raindrops possess diverse microphysical properties over two different altitude sites. Over the Longmen, stratiform precipitation is dominated by raindrop collision and coalescence processes; in particular, from 2.0-2.8 km, the breakup process competes with collision-coalescence processes but later is overpowered. In contrast, over the Nagqu, the dominant processes are collision and coalescence processes. Comparisons of raindrop spectra suggest that the concentration of small raindrops over the Nagqu is much higher than that over the Longmen, whereas there is a greater concentration of medium-to-large raindrops over the former than over the latter. Therefore, D m (log 10 N w ) over the Nagqu is smaller (larger) than that over the Longmen.
The MRR-measured results are consistent with the disdrometer results for medium raindrops over the two sites, but concentrations are underestimated to some extent. However, for large raindrops, radar can slightly overestimate concentrations. The discrepancies between the two instruments may be caused by differences in measuring principles, sampling heights, equipment limitations, and the natural variations of raindrops. The averaged log 10 N w -D m pairs over the Longmen and Nagqu sites are compared with two previous results for summer precipitation over Central and East China, which also exhibit diverse behaviors. For the Longmen, effective collision and coalescence processes lead to D m s larger than those over the other aforementioned regions-correspondingly, log 10 N w s can be smaller. For the Nagqu, because of the short falling path of raindrops, collision and coalescence processes are insufficient, leading to a relatively high concentration of small raindrops falling to the ground. Thus, it may be concluded that the altitude, water vapor, and atmospheric conditions are the main explanations for the difference between two sites with different altitudes.