Intercomparison of Cloud Vertical Structures over Four Different Sites of the Eastern Slope of the Tibetan Plateau in Summer Using Ka-Band Millimeter-Wave Radar Measurements

: The eastern slope of the Tibetan Plateau is a crucial corridor of water-vapor transport from the Tibetan Plateau to Eastern China. This is also a region with active cloud initiation, and the locally hatched cloud systems have a profound impact on the radiation budget and hydrological cycle over the downstream Sichuan Basin and the middle reach of the Yangtze River. It is noteworthy that there is a strong diversiﬁcation in the characteristics and evolution of the ESTP cloud systems due to the complex terrain. Therefore, in this study, ground-based Ka-band millimeter-wave cloud radar measurements collected at the Ganzi (GZ), Litang (LT), Daocheng (DC), and Jiulong (JL) sites of the ESTP in 2019 were analyzed to compare the vertical structures of summer nonprecipitating clouds, including cloud occurrence frequency, radar reﬂectivity factor, cloud base height, cloud top height, and cloud thickness. The occurrence frequency exhibits two peaks on the ESTP with maximum values of ~20% (2–4 km) and 15% (7–9 km), respectively. The greatest (smallest) occurrence frequency occurs in the JL (GZ). The cloud occurrence frequency of all sites increases rapidly in the afternoon, and the occurrence frequency of the DC presents larger values at 2–4 km. In contrast, the occurrence frequency in the JL shows another increase from 2000 LT to midnight at 7–11 km. Stronger radar echoes occur most frequently in the LT at 5–7 km, and hydrometeor sizes and phase states vary dramatically in mixed-phase clouds. A small number of radar echoes occur at midnight in the JL. A characteristic bimodality of the cloud base height and top height for single-layer, double-layer, and triple-layer clouds was observed. Clouds show a higher base height in the GZ and higher top height in the JL. The ESTP is dominated by thin clouds with thicknesses of 200–400 m. The cloud base height, top height, and thickness exhibit an increase in the afternoon, and higher top height occurs more frequently from midnight to the next early morning in the JL because of its mountain-valley terrain. dominate in the GZ with 38% occurrence frequency. The four-site observations show that cloud base height, top height, and thicknesses increase signiﬁcantly in the afternoon due to the convection on the ESTP in summer. Higher top heights occur more frequently at midnight and early morning in the JL due to its mountain-valley terrain.


Introduction
The Tibetan Plateau (TP), known as the third pole of the world, is the highest and largest plateau in the world [1,2]. The TP plays a significant role in regulating the atmospheric circulation and modulating Asian climate [3][4][5]. Because of its high terrain elevation, the observation of clouds over the TP and its vicinity heavily relies on passive satellite instruments, and various studies have been conducted to investigate cloud formation Currently, by using satellite data, most studies have examined the large-scale cloud characteristics over the main body of the TP. For the ground-based observations, analysis chiefly focused on the cloud vertical structures at the central TP (in Nagqu) [22][23][24], as well as convective activities and precipitation of the southern TP [8][9][10]. However, there is a lack of attention on the cloud vertical structure and diurnal statistical analysis on the eastern slope of the TP (ESTP). Clouds over the ESTP can evolve and eventually form heavy precipitation, which has a profound impact on the Sichuan Basin and the middle reach of the Yangtze River [27]. Geographically, the ESTP is a transitional area between the main body of the TP and the downstream Sichuan Basin. It is also a crucial corridor for the eastward water-vapor transport from the TP. In addition, the ESTP is near the Asian monsoon region and is dominated by a special circulation system induced by steep topography and complex terrains. Under this condition, the water-vapor distribution and internal circulation are cross-scale in the area, which makes it conducive to localized cloud formation. Owing to data limitations for this complex region, knowledge of cloud characteristics on the ESTP remains insufficient. Therefore, multisite ground-based observations with high spatiotemporal resolution on the ESTP are synthesized and analyzed in this study. These observations will help understand regional cloud characteristics on the ESTP, providing an observational basis and constrains for the potential improvement of numerical weather forecast models.
In monsoon summer during June and August in 2019, with joint support from the National Natural Science Foundation of China and the Institute of Heavy Rain, China Meteorological Administration, a field experiment on clouds and precipitation observations was conducted on the ESTP [28]. In this experiment, a Ka-band millimeter-wave cloud radar (Ka-MMCR), an X-band dual-polarization radar, a wind profile radar, and a microwave radiometer, were deployed at the meteorological station of Ganzi (GZ), Litang (LT), Daocheng (DC), and Jiulong (JL), respectively. The purpose of this experiment was to provide cloud properties over the ESTP from these ground-based observation instruments.
As shown in Figure 1, the four observation sites lie on the steep ESTP, which is the junction between the main body of the plateau and the downstream Sichuan Basin. The observation region (denoted by the red dotted rectangle) slopes from northwest to southeast in the ESTP, with high mountains and deep valleys staggered. Topographically, the GZ [31.62 • N, 100.00 • E, 3395 m above sea level (ASL)] is located in the northern part of the observation region, and is relatively close to the central TP. The GZ lies in a flat valley between the mountains, with the underlying surface of deep meadows. The LT (30.00 • N, 100.27 • E, 3951 m ASL) and DC (29.05 • N, 100.30 • E, 3727 m ASL) are located on the mountains concentrating around the southeastern edge of the TP, east of the Jinsha River [7]. The relative humidity in the LT and DC might be much higher than that in the GZ, due to the accumulation of water vapor in the mountain areas. The JL (29.00 • N, 101.50 • E, 2925 m ASL) is located in the southeastern part of the observation region. It possesses a mountain-valley terrain with a large altitude gradient and is the main source of the China southwest vortex [29]. These four-site observations with different topographic features provided a chance to investigate the characteristics of regional cloud systems and developments on the ESTP during the boreal summer in 2019.
Based on the four-site Ka-MMCR observations over the ESTP, this study statistically analyzed the vertical distribution and diurnal variation of cloud occurrence frequency, radar reflectivity factor, cloud base height (CBH), CTH, and cloud thickness. The remainder of this paper is organized as follows: In Section 2, the instruments, data, and method are described. Section 3 exhibits the statistical results of the cloud occurrence frequency, radar reflectivity factor, CBH, CTH, cloud thickness, and their diurnal variations. Discussions on the results are presented in Section 4. Conclusions are drawn in Section 5.

Radar Observations and ERA5 Reanalysis Data
To continuously obtain the aloft clouds, the zenith Ka-MMCRs manufactured by the Sun Create Electronics Co., Ltd. of China were employed at the GZ, LT, DC, and JL, respectively. All radars work at 35 GHz, with a wavelength of 8.6 mm and a beamwidth of 0.3 • . The radars can provide vertical profiles of the reflectivity factor (Z e , dBZ), mean Doppler velocity (V M , m s −1 ), spectrum width (σ V , m s −1 ), the linear depolarization ratio (LDR, dB), and Doppler spectra up to 20 km above the ground with a spatial resolution of 30 m and a temporal resolution of 10 s. Three operation modes are designed in Ka-MMCR to observe clouds at different altitudes [30]. In addition, the four Ka-MMCRs calibrations were accomplished before the field experiment. Take the radar calibration in the GZ for instance: external signal source ranging from −90 dBm to −40 dBm was input to the radar receiver, to obtain the values of Z e at the height of 1 km. Meanwhile, the values of Z e were calculated by the meteorological radar equation, based on the constants and measured parameters. Maximum difference of Z e of copolarization channel is 0.47 dB, and maximum difference of Z e of quadrature polarization channel is 0.94 dB. Intensive cloud observations were simultaneously conducted at the four observation sites from 5 June to 20 August 2019.
The atmospheric circulation patterns over the sampling region were analyzed based on ERA5 reanalysis products (https://cds.climate.copernicus.eu/cdsapp#!/search?Text= ERA5 (1 January 2022) of a 0.25 • × 0.25 • spatial resolution, hourly data on pressure levels from 1979 to the present. To reach the temporal consistency, hourly ERA5 data during the observation period were matched to the ground-based Ka-MMCR observation. The temperature and moisture profiles were supplied by the reanalysis products, to reveal the regional atmospheric conditions of the observational region.

Radar Data Quality Control
During the observation, both cloud echoes and nonmeteorological echoes were totally included in the radar reflectivity factor. In the low-and mid-latitude regions, nonmeteorological echoes observed by Ka-MMCR are mostly caused by low-level plankton, which consists of dust, insects, pollen, and other targets [31,32]. Plankton targets are mainly present within the planetary boundary layer and can affect the data usage of Ka-MMCR. Due to the existence of the nonmeteorological echoes, CBH is underestimated with uncertainties. For separation of the plankton contamination from the radar dataset, Z e and LDR probability distributions of plankton and clouds were investigated [33]. According to the method from a previous study in Nagqu, the nonmeteorological echoes contribute to radar-measured Z e with a value less than 0 dBZ and LDR with a value larger than −16 dB [34]. Herein, the algorithm was adopted to separate the nonmeteorological echoes in this study.
In addition, radar range sidelobes also should be filtered before the statistical analysis. To detect middle-and high-level clouds, the wide pulses and the pulse compression technologies are used in Ka-MMCR, which results in the range sidelobes. These artifacts prevail near regions of clouds that contain great values of reflectivity factor and Doppler velocity in the radial direction [35]. The presence of these artifacts results in apparent smearing of echo from intense cloud echoes into the data of the adjacent height gates, which will further obstruct the retrieval of cloud macroscopic parameters, such as CBH (underestimated) and the CTH (overestimated). For vertical-pointing radars, the range sidelobe clutters will probably present above the top or below the base of the strong clouds with great values of reflectivity factor. In the study, the range sidelobe clutters are removed using Moran's method [35].
To present the effect of the radar data quality control, a typical case that contains nonmeteorological echoes and range sidelobe clutters is plotted in Figure 2. The observation exhibits that clouds successively moved over the Ka-MMCR deployed in the LT on 9 July 2019. In Figure 2a, the clouds are surrounded by nonmeteorological echoes (denoted by red arrows) and are seriously impacted by the range sidelobe clutters (denoted by black arrows). The nonmeteorological echoes were concentrated on the low-level lower atmosphere, occurring at the height of 0-1 km around 0700-1400 local time (LT). These echoes were frequently mistaken for fractus near the ground, which brings an obstacle to cloud layer detection. Meanwhile, relatively more intensive clouds with greater radar echoes approaching 30 dBZ developed, and the clouds were confused by the range sidelobes (marked by black arrows) during 0900-2000 LT, which were determined by the Moran's criterion. On one hand, range sidelobes (0900-1400 LT) lying at the cloud bottom present a reflectivity value of~−25 dBZ, around 45 dB smaller than the strong echo core (20 dBZ) in clouds at the height of~4.5 km. On the other hand, range sidelobes (1400-2000 LT) connecting with the cloud top show a reflectivity value of~−20 dBZ, around 50 dB smaller than the strong echo core (30 dBZ) in clouds near the ground. Furtherly, the existing range sidelobes in between 1400-1500 LT fuzz the boundaries of adjacent cloud layers, which might cause errors in the following statistics of the CBH, CTH, and cloud thickness. As shown in Figure 2b, the scattered nonmeteorological targets were well-eliminated after the quality control, and the range sidelobe clutters were also removed.

Profile Classification
Millimeter-wavelength signals are usually used to detect the small-size hydrometeors in clouds, and the occurrence of precipitation on the ESTP could lead to strong attenuation of the Ka-MMCR signal [36,37]. Because of the presence of precipitation, the cloud top height obtained by the Ka-MMCR tends to be lower than reality [38]. Meanwhile, the bottom of precipitating clouds is close to the surface, and it is difficult to verify the cloud base height using the Ka-MMCR. In general, important parameters (e.g., cloud base height and top height) during the cloud formation and development stages will be inaccurately detected by the Ka-MMCR with collocated rainfall. Therefore, nonprecipitating clouds observed by Ka-MMCR are focused in this statistical analysis. The non-precipitating cloud profiles were identified and selected from the total radar profiles by the following steps. First, the entire radar profiles were divided into two categories, namely the cloud profiles and the clear-sky profiles. A profile was defined as the cloud profile when the radar echo includes hydrometeors; otherwise, it was treated as the clear-sky profile. Subsequently, the cloud profiles were further classified into two types-the nonprecipitating cloud profile and the precipitating cloud profile-by judging if the simultaneous ground rain rate (observed by the collocated surface disdrometers) was greater than 0.01 mm h −1 .
After the classification, different types of sample profiles from the four observation sites are summarized in Table 1. The N c is the number of cloud profiles, the N npc is the number of nonprecipitating cloud profiles and the N total is the total number of all sample profiles. Statistics show that a total of 276,480 profiles (N total ) were observed by the Ka-MMCR at each observation site from 5 June to 20 August 2019. From all the here-analyzed observations, both the largest N c (144,892) and N npc (121,934) were observed in the JL.

Cloud Occurrence Frequency
Cloud occurrence frequency refers to the height-resolved distributions of cloud occurrence measured by the Ka-MMCRs. This analysis is carried out on the nonprecipitating cloud profiles. Hence, the cloud occurrence frequency is calculated as the ratio of the number of nonprecipitating cloud profiles (N npc ) to the total number of the observed profiles (N total ). Figure 3 shows the vertical distribution of cloud occurrences at four sites. Original profiles of cloud occurrence frequency are presented in Figure 3a, and smoothed profiles are presented in Figure 3b. Four original profiles roughly show a bimodal distribution, respectively. These profiles contain small values (5%-10%) below~1.5 km, which are probably caused by a small amount of clouds in the low-level atmosphere. To unravel the overall trends of the distribution of nonprecipitating clouds with different underlying surfaces, these profiles are smoothed. As shown in Figure 3b, there are two peaks located at the height of 2-4 km (main peak) and 7-9 km (secondary peak), respectively. This finding is in good agreement with the studies in the Nagqu: the values of cloud occurrence frequency show two peaks (located at 2-4 km and 6-9 km separately) in summer [39]. Nevertheless, the four-site observations on the ESTP display significant regional differences in occurrence frequency peak values and locations. The results show that maximum values (main peak) of occurrence frequency in the DC and JL are~20%, which is significantly larger than in the LT and GZ. Conversely, the main peak location of occurrence frequency in the GZ is higher than 3.5 km, which is slightly higher than in other sites. Meanwhile, the maximum value (secondary peak) of occurrence frequency in the JL is larger than 15%, which is dramatically larger than in other three sites. The secondary peak location of occurrence frequency in the JL is higher than the counterparts of the other sites, with the height of 8.5 km. In general, the occurrence frequency in the JL is larger than in other sites, with the cloudiest conditions among the four observation sites. Meanwhile, the occurrence frequency in the GZ is the smallest, which indicates that cloud formation in this region is relatively difficult. The ESTP is the submaximum elevation-gradient region of China, and strong ascending motion occurs in this region due to the topographic forcing and the upslope mountain-valley circulation [12,40]. Among the four observation sites, the JL (2925 m, AGL) is located in the canyon, with the largest altitude gradient. Deep-cut rivers flow from north to south in this region, and the rivers play an important role in guiding water-vapor transport. Thereby, more warm and humid airflow can facilitate the cloud formations in the JL. This supports the Ka-MMCR observations that the highest cloud occurrences occur both in the low-level and high-level troposphere in the JL. In contrast, the water-vapor amount in the GZ is relatively insufficient, because the moisture airflow splits into several branches during the transport. Meanwhile, the vertical movement of airflow becomes weak due to the decreasing altitude gradient in the GZ (3395 m, AGL). Under this condition, the occurrence frequency in this region presents the smallest values. These may suggest that the local underlying surface plays an important role in controlling the occurrence of the clouds. To further understand the regional differences in cloud occurrence frequency on the ESTP, the patterns of moisture convergence at four sites will be analyzed in Section 4 in detail.
To describe the cloud diurnal variations at four sites, the cloud occurrence frequency is calculated every 60 min with 30 m vertical resolution. Figure 4a-d show the temporally normalized values of cloud occurrence frequency at different AGL heights. Figure 4e presents the hourly changes of cloud occurrence frequency, which is the summation of the occurrence frequency of all heights within 60 min. Regional differences in cloud occurrence are apparent throughout the day. In the DC, observations show two cloud layers, associated with the diurnal cycle and convective instability. In the lower layer, the maximum occurrence frequency (2%) is observed between 1200 LT and 1800 LT at an altitude extending from~0.5 km to 5 km. In the upper layer, the maximum value of occurrence frequency (1%) is located at 7-10 km during 1600 LT and 1900 LT, which probably corresponds to deep convective clouds. The occurrence frequency in the JL shows a similar diurnal distribution to the DC, whereas in the lower layer, the maximum value of occurrence frequency (2%) in the JL presents a shorter period (1300-1500 LT), and is located at 3-4 km. In the upper layer, the maximum occurrence frequency (1%) in the JL is observed from 1600 LT to the early morning. In the LT site, observations show only a lower cloud layer; the maximum value of occurrence frequency (2%) is observed between 1500 LT and 1600 LT at an altitude extending from~4 km to 5 km. In contrast, the occurrence frequency in the GZ shows unapparent variation, with the maximum value of occurrence frequency (1%) smaller than other sites. Previous studies showed that diurnal differences in the intensity of convection have a significant impact on the cloud diurnal variations [41,42]. The TP is a heat source in summer, enhances the surface heating, and further enhances the convections over the TP [43,44]. Surface heating from solar radiation causes instability at lower levels, which leads to moist convection and cloud formation [45]. Under this condition, cloud occurrence frequency sharply increased in the afternoon. Furthermore, cloud occurrence frequency of the four sites does not increase at the same time due to the regional triggering time of convective processes [46]. Figure 4e shows that the occurrence frequency of the four sites increases rapidly in the afternoon (1300-1800 LT), then decreases slowly in the morning (0000-0500 LT), and reaches a minimum value around 1030 LT. In Figure 4e, the hourly occurrence frequency of clouds in the LT site begins to increase around 1030 LT earlier than other sites, probably because the convective process in the LT is triggered firstly among the four observation sites. On the other hand, cloud occurrence frequency in the JL shows a sudden increase at midnight, which is associated with the occurrence frequency increase at 7-11 km (Figure 4d). The results are similar to the findings from Zhou et al. [47], which reported that cloud occurrence frequency of the southern TP (e.g., Motuo) remained near the maximum value (40%) from evening to midnight, then decreased rapidly in the morning. The maximum values of occurrence frequency throughout the day in the DC, JL, LT, and GZ are 78%, 70%, 55%, and 50%, respectively. Nonprecipitating clouds form more frequently on the ESTP than on the southern TP.

Radar Reflectivity Factor
The radar reflectivity factor (Z e ) measured by Ka-MMCR can represent the distribution of hydrometeor size and concentration in nonprecipitating clouds. According to the meteorological radar equation [48], the size and phase state (liquid/mixed/ice phase) of hydrometeors can directly affect the value of Z e . Particularly, hydrometeors with large diameters in clouds contribute to the great values of Z e . Great values of Z e are related to strong clouds, and narrow variation ranges of Z e represent good consistency of hydrometeor size and concentration in clouds. More importantly, the temperature is a key factor of hydrometeor phase states, and regional temperature changes have an impact on the cloud intensity. Figure 5 shows the gradient changes of mean temperature from ERA5 reanalysis data in four sites between 5 June and 20 August in 2019. The temperatures decrease linearly with increasing altitudes among four sites, then the 0 • C isotherm and the −40 • C isotherm are located in between at altitudes of~1.6-2.6 km and~7.8-8.8 km, respectively. Furthermore, studies have shown that hydrometeors exist as liquid water in clouds above 0 • C, mixed-phase clouds (including supercooled water and ice crystal) might stay from 0 • C down to −40 • C, and nearly pure ice clouds (including ice crystals) are located in the high-level atmosphere below −40 • C [49]. As shown in Figure 6, the height-resolved variations of Z e are calculated in boxplots (Figure 6a-d), and the mean value profiles of Z e are superimposed (Figure 6e). Significantly, the mean values of Z e are calculated in log scale. Boxes denote the 25th and 75th percentile positions, whereas the lines and dots inside the box show the median and mean values, respectively. To ensure that each layer has samples for statistical robustness, the data were resampled to a vertical resolution of 500 m. The boxplots of Z e for non-precipitating clouds differ remarkably among the four observation sites. Above the height of −40 • C isotherm (~9 km), the values of Z e range from −25 dBZ to −5 dBZ. With the height decrease, the values of Z e of four sites firstly exhibit a slight increase, then present a rapid decrease. Compared to the same height in other sites, the boxes cover a wider dBZ range in the GZ (Figure 6a). In between the height of 0-−40 • C isotherm (~1.7-9 km), the values of Z e increase rapidly with the altitude decrease and reach the maxima of~−10 dBZ at 5-7 km. During this process, the radar reflectivity factors present a wider range, which probably indicate that particle size and concentration in clouds varies significantly among the four sites. In between the height of~2-6 km, the values of Z e decrease with the height decrease. Meanwhile, the boxplots cover a narrower dBZ range in the GZ, LT, and JL, while the dBZ range becomes wider in the DC (Figure 6c). The values of Z e maximize (−10-−5 dBZ) below of~0 • C isotherm (~2 km), which suggests that particles in clouds began to melt in this layer. Below~1.5 km, the values of Z e start to decrease rapidly. As shown in Figure 6e, the mean values of Z e of four sites present the peak values at 5-7 km. In between the height of 5-7 km, the maximum value of mean Z e in the LT is larger than that of other sites, while the maximum value of mean Z e in the GZ is smallest. In other words, mixed-phase clouds at 5-7 km present relatively stronger in the vertical direction. Among the four observation sites, nonprecipitating clouds in the LT are the strongest. To better understand the daily distribution of radar echoes at four sites, mean Z e of log scale is calculated every 60 min with 30 m vertical resolution. Figure 7a-d present the height-resolved temporal evolutions of mean Z e at different AGL heights. Figure 7e shows the hourly changes of mean values of Z e , which are the mean values of Z e of all AGL heights within 60 min. As shown in Figure 7b, the maximum Z e (~10 dBZ) occurs in the morning (0500-0800 LT), and in the afternoon to night (1800-2300 LT) in the LT site. The radar echoes with maximum value in the afternoon are mainly located at 5-7 km, which is consistent with the result of the LT site in Figure 6e. Similarly in Figure 7a, radar echoes in the GZ exhibit maximum value (~10 dBZ) in the morning (0600-0800 LT), in the afternoon (1400-1500 LT), and at night (2200-2300 LT). In the DC, the maximum Z e (~10 dBZ) occurs in the afternoon (1700-1800 LT), with the height of 2-5 km (Figure 7c). In the JL, radar echoes with maximum value (~10 dBZ) occur in the afternoon (1600-1900 LT) at 1-2 km, and in the midnight (0000-0100 LT) at 6-7 km, respectively. Importantly, large values (0 dBZ) of Z e at four sites occur in the afternoon, which is most probably associated with the convective processes [43,44]. From the Ka-MMCR observations, strong radar echoes mainly consist of ice crystal, because these clouds are located in between the height of 0-−40 • C isotherm. These results show that hydrometeor sizes and phase states vary dramatically in mixed-phase clouds on the ESTP, and the localized convections induce more complicated microphysics processes at four sites. Generally, the large values of Z e are caused by the strengthening of mixed-phase clouds. In Figure 7e, mean values of Z e increase in the afternoon, and maximize around 1500-1700 LT, then decrease from the night to the next morning. The minimum values occur around 1030 LT, which is consistent with the valley value occurrence of cloud occurrence frequency (see Figure 4e).

Cloud Macroscopic Parameters
Cloud macroscopic parameters measured by cloud radars, such as the CBH, CTH, and cloud thickness, are widely used to analyze the cloud vertical structures in recent years [38,47,50,51]. To verify the reliability of these parameters measured by the cloud radar, previous studies have compared the CBH (CTH) from the cloud radar with the ceilometer (FY-4A) and found a relatively high consistency between the two datasets [44,46]. In addition, the four Ka-MMCRs used in this study were well-calibrated before the field experiment. In this work, CBH and CTH from the Ka-MMCRs are confirmed using Wu's method for cloud boundaries of the TP [51]. In Wu's method, the lower boundary of nonprecipitating clouds is defined as the CBH, and the upper boundary of nonprecipitating clouds is defined as the CTH, corresponding to the post-quality-control data. Herein, the cloud thickness is a geometrical thickness, which is defined as the difference between the top height and base height. In the vertical direction, cloud clusters overlapped with an interval less than 90 m are treated as the same-layered cloud; otherwise, they are different-layered clouds. Then, the thickness of each layer and the space between adjacent cloud layers are determined. From the previous analysis, clouds with a layer number of no more than three occur more frequently over the TP and the adjacent areas, with the maximum occurrence frequency larger than 95% [52,53]. Therefore, the CBH, CTH, and cloud thickness of the single-layer, double-layer, and triple-layer clouds from the radars are considered, and clouds with a layer number of more than three will be neglected in this section.
The height-resolved Ka-MMCR observations enable the determination of the spatial evolution of CBH, CTH, and cloud thickness. CBH is related to water-vapor condensation height; CTH and cloud thickness present cloud vertical extension. Significantly, the vertical distribution of these parameters can reveal regional cloud formation and development on the ESTP. The probability density function (PDF) profiles of the CBH, CTH, and cloud classification results are displayed in Figure 8. The definition of different cloud layers is also interpreted. Here, Layer-1, Layer-2, and Layer-3 present the lowest-layer clouds, middle-layer clouds, and highest-layer clouds, respectively. Because the clouds in lowlevel atmosphere have little impact on the bimodal distribution of clouds in the vertical direction at four sites (see Figure 3), the profiles of CBH and CTH are processed by fitting in Figure 8a-c,e-g. From all the here-analyzed observations, the CBH profiles of all layers present a bimodality in the vertical direction (Figure 8a-c). For the lowest-layer clouds, the main peaks of the CBH profiles are located at 1-4 km, and the secondary peaks of the CBH profiles are located at 6-8 km. Note that the main peak of the CBH of the lowestlayer clouds in the GZ is higher by 1km than other sites. For the highest-layer clouds, the main peaks of CBH profiles are located at 7-10 km; particularly, the peak in the JL is higher by 1 km than in other sites. Meanwhile, different categories of clouds can be classified by the values of CBH. According to the cloud classification method proposed by the WMO [54], the CBH results are used to categorize the clouds as high (CBH > 5 km), middle (2.5 km ≤ CBH ≤ 5 km), and low (CBH < 2.5 km) clouds, as shown in Figure 8d. In the GZ, middle clouds occur most frequently (38%), then the occurrence frequency of low clouds and high clouds are 34% and 27%, respectively. Differently, low clouds occur most frequently both in the DC (46%) and the LT (38%), then middle clouds occur second-most frequently (35% in the LT, 29% in the DC), and the occurrences of high clouds are fewest (27% in the LT, 26% in the DC). In contrast, the occurrence frequency value of low clouds is maximum (39%) in the JL, followed by the high clouds (33%) and middle clouds (28%). Consistent with the CBH distributions, the CTH profiles play a bimodal pattern in the vertical direction (Figure 8e-g). For different layer clouds, the main peaks of CTH profiles in the GZ are located at 4.5-5 km, and the peaks are significantly higher than in other sites. Furthermore, the secondary peak of the CTH profile in the JL is located at~10 km, and the peak is also apparently higher than in other sites. Clouds in the JL generally show a higher (~0.5 km) top height among the four sites, which suggests that a more drastic vertical air motion happened in this region in the 2019 summer. PDF plots of the thickness of single-layer, double-layer, and triple-layer clouds are illustrated in Figure 9. The PDF profiles of all cloud layers show a unimodal distribution. As shown in Figure 9, the maximum value of the PDF corresponds to the central thickness. For the same observation site, the single-layer clouds exhibit a minimum value of PDF, the double-layer clouds present a second minimum value of PDF, and the triple-layer clouds show a maximum value of PDF (Figure 9a). Central thickness of single-layer clouds presents a maximum value, central thickness of double-layer clouds presents a second maximum value, and central thickness of triple-layer clouds presents a minimum value. Meanwhile, the thickness of single-layer clouds extends to~8 km, the thickness of double-layer clouds extends to~6.6 km, and the thickness of triple-layer clouds extends to~6 km. Because of the increase in cloud layer number, the peak value of PDF increases, the central thickness decreases, and the cloud thickness extension also decreases. In other words, clouds become thinner when the number of cloud layers increases, which is in good agreement with the compression effect of TP terrain on the vertical cloud structures reported by previous studies [10,52]. In general, the ESTP is dominated by thin clouds with thicknesses of 200-400 m. Among the four observation sites, the maximum value of PDF profiles in the LT is larger, and the cloud thickness exhibits smaller values. This result potentially shows that clouds in the LT are more strongly compressed by the high-altitude terrain. To further investigate the diurnal variations of water-vapor condensation height and cloud vertical extension, diurnal evolutions of base height of the lowest-layer clouds (CBH-1), top height of the highest-layer clouds (CTH-3), and the thickness of CTH-3 and CBH-1 are displayed in Figure 10. Significantly, the base height of CBH-1 is closely associated with local underlying surface, and the top height of CTH-3 largely exhibit the localized cloud vertical extension. The data of CBH-1, CTH-3, and the thickness of CTH-3 and CBH-1 are calculated every 60 min with 30 m vertical resolution. Statistics show that 50% (25th-75th of the data included by the box) of the base heights range from 0.24 km to 6 km (Figure 10a,d,g,j). Cloud base heights in the GZ display a gentle decrease throughout the morning until~1500 LT, and increase rapidly to the peak at around 2200 LT. In contrast, cloud base heights in the LT roughly show a bimodality during the day, with one peak from midnight to around 1300 LT and the other peak from the afternoon to the late night. In the JL, clouds exhibit much lower base heights in the middle night and early morning, which is likely associated with the warm and humid airflow. Compared with the base heights, 50% of top heights reveal a bimodal pattern ranging from 3 km to 10 km, with the peak during early morning and afternoon (Figure 8b,e,h,k). In the GZ, one peak of top height occurs at 1900-2100 LT and the other peak occurs at 0200-0400 LT. The diurnal variation of top heights in the LT is similar to that in the GZ, with the peak values in the early morning (0300-0400 LT) and late night (2000-2300 LT). The peaks of top heights in the DC occur during 0000-0300 LT and 1500-2000 LT. In the JL, cloud top heights show a minimum value at 1200 LT, rapidly increase during 1300-1700 LT, and maximize at 2000 LT, then drop from the night to the early morning. Significantly, higher cloud top heights occur more frequently from midnight to the next early morning in the JL, which might be correlated with the collision of the up-valley wind from summer monsoon airflow with the down-valley wind [55]. Interestingly, the observations show that cloud top heights increase earlier than the base heights increase in the afternoon, and the diurnal variation of cloud top heights are more significant than those of base heights. The daily changes in the thickness are similar to the evolution of top heights, with two peaks in the afternoon and early morning, respectively (Figure 10c,f,i,l). It is most probable that the convectively unstable atmosphere might thereby thicken the clouds. These results are generally in good agreement with previous finds from radiosondes on the ESTP and other regions in China [56].

Discussion and Comparison
The ESTP region is well-known to act as a convergence zone during the summer monsoon [44]. Geographically, the ESTP is an important part of the TP and absorbs the solar radiation strongly. The large sensible heating in this region results in convergences, which form ascending movements throughout the day [57]. Therefore, local atmospheric circulation presents some differences in vertical and daily changes, which might impact the cloud properties of four-site observations on the ESTP.
For instance, the diurnal variation of cloud occurrence frequency is closely related to the diurnal cycle of the convergences. The airflow ascending movements caused by convergences are beneficial to the water-vapor condensation. Figure 11 shows the diurnal variations of moisture flux divergence at four sites. The values of moisture flux divergence are averagely calculated by the ERA5 reanalysis data every 60 min. Data interpolation is used in the vertical direction. The negative values of moisture flux divergence present convergence, and smaller negative values indicate stronger convergence. As shown in Figure 11b-d, the surrounding airflows converge (denoted by a negative contour value) near the surface in the LT, DC, and JL, from the afternoon to the evening. Strong convergences occur in the afternoon, aggravate the cloud process, and result in a great increase in occurrence frequency in these sites (see Figure 4b-d), respectively. The minimum value of moisture flux divergence (<−30 × 10 −6 g/cm 2 ·hPa·s) in the DC and JL is smaller than in the LT, which indicates that stronger convergences happen in the first two sites. From the Ka-MMCR observations, the values of occurrence frequency in the DC and JL are larger than in the LT (see Figure 4e). In contrast, unapparent diurnal variation cycles of convergence occur in the GZ, which leads to unapparent variations in cloud occurrence frequency throughout the whole day (see Figure 4a). Based on the early studies, the thermal forcing of heating increases the potential temperature of the air column, and the induced mechanical forcing results in convergence in lower levels but divergence in upper levels [58]. From the ERA5 reanalysis data in the JL, airflows diverge (denoted by a positive contour value) strongly above the convergence zone, with a maximum value of moisture flux divergence larger than 18 × 10 −6 g/cm 2 ·hPa·s. Under this condition, the air ascending motion is strongly suppressed, which may weaken the local cloud formation process. Hence, from the observations, the occurrence frequency values of the lower troposphere in the JL are generally smaller than in the DC during the afternoon (see Figure 4c-d). In general, unique cloud properties on the ESTP in summer are caused by the effects of topography, irradiation, and heat source [59]. Numerous studies show large differences in cloud vertical structures and diurnal variations between the ESTP and Eastern China. By using the Cloudsat/CALIPSO data, the distribution of hydrometeors also exhibits two peaks in summer in Eastern China [8]. The first peak is located in the upper troposphere, with the maximum value at higher altitudes (12-15 km) in Eastern China than that (7-9 km) in the ESTP (see Figure 3). The second peak is located at approximately the same height (1-3 km) in Eastern China as that (2-4 km) in the ESTP. For the distribution of dBZ for nonprecipitating clouds, the observations in Eastern China (e.g., Beijing) exhibit different vertical variations during summer [47]. Firstly, the dBZ values of clouds in Beijing are generally larger than that in the ESTP, with the maxima of~40 dBZ (Beijing). The result shows that the clouds on the ESTP (adjacent region of the TP) are weaker than in other regions, which is consistent with the previous studies [8,9]. Secondly, the dBZ profile (mean value) of mixed-phase clouds present a linear decrease from~9 km downwards to~5 km (~0 • C isotherm) in Beijing. The observations on the ESTP (see Figure 5) show that the mean dBZ values of mixed-phase clouds rapidly maximize (−5-0 dBZ) in between the height of 6-9 km, and decrease to~−9 dBZ at~2 km (~0 • C isotherm). There are likely more complicated physical processes in mixed-phase clouds on the ESTP. In addition, the cloud occurrence frequency presents a sharp increase from 1500 LT in Beijing during summer and reaches the maxima (60%) in the early morning [60]. Due to the heat source of the TP in summer, cloud occurrence frequency of the ESTP exhibits a rapid increase around 1300 LT earlier than that in Beijing. Furthermore, the daily occurrence frequency of clouds shows a maximum of 80% larger than that in Beijing, which indicates ESTP has a higher water-vapor content in the summer [59].
However, there are still some uncertainties and limitations of the results, which should be mentioned as follows: a small number of the nonprecipitating cloud profiles are possibly neglected in the statistical analysis. Based on the radar profile classification method in Section 2.3, nonprecipitating cloud profiles are determined by minutely rain rate. Therefore, six adjacent cloud profiles in one minute are classified as precipitating cloud profiles or nonprecipitating cloud profiles. Therefore, a small quantity of nonprecipitating cloud profiles are out of the statistical analysis, when short-time rainfall occurs in less than one minute. Furthermore, previous analysis of the International Satellite Cloud Climatology Project (ISCCP) data revealed that regional differences in global cloud properties are closely related to the dynamic effects of thermodynamic and microphysical processes within clouds [61]. In this paper, the present findings concentrate on the macroscopic characteristics of the clouds on the ESTP, such as cloud occurrence frequency, radar reflectivity factor, CBH, CTH, and cloud thickness. To deeply investigate the regional cloud properties on the ESTP, more microphysical characteristics of clouds will be studied in further work. Studies about precipitating clouds on the ESTP by combining multisource observations (e.g., the space-borne active Ka/Ku-band radars [62]) will also be carried out in the future.

Conclusions
To better understand the cloud physical properties in the critical region of the TP, field observations on the ESTP were launched by the Institute of Heavy Rain, China Meteorological Administration, in 2019. This study explored the nonprecipitating cloud climatological characteristics across the ESTP using continuous Ka-MMCR observations collected in 2019 in four sites. For the first time, regional cloud vertical structures and diurnal variations on the ESTP were statistically analyzed in detail. The main conclusions are as follows.
Cloud occurrence frequency obtained from the Ka-MMCR measurements at GZ, LT, DC, and JL from 5 June to 20 August 2019, were compared. In general, the distribution of hydrometeor exhibits two peaks in summer on the ESTP. The first peak is located at 7-9 km, with a maximum occurrence frequency of 15% in the JL. The second peak is located at 2-4 km, with the maximum occurrence frequency of~20% in the DC and JL. The maximum occurrence frequency throughout the day in the DC, JL, LT, and GZ are 78%, 70%, 55%, and 50%, respectively. Generally, cloud occurrence frequency of the ESTP increases rapidly in the afternoon (1300-1800 LT), and decreases slowly in the morning (0000-0500 LT), then minimizes around 1030 LT. In the DC, the occurrence frequency presents larger values (2%) from 1200 LT to 1800 LT at 2-4 km. In contrast, the occurrence frequency in the JL shows an increase from 2000 LT to midnight in the upper troposphere (7-11 km).
In the terms of the distributions of Z e , nonprecipitating clouds on the ESTP present the mean values of reflectivity factors from −25 dBZ to −7 dBZ in the vertical direction. For mixed-phase clouds, the values of Z e show a rapid increase at~9 km, and maximize (−5-0 dBZ) at 5-7 km, then sequentially decrease downwards to~3 km. Compared to other sites, clouds in the GZ are much weaker, with the maximum of mean Z e with~−5 dBZ. During the day, the values of Z e increase from the afternoon (1300 LT), and maximize in the late evening. Among the four sites, stronger radar echoes occur most frequently in the LT, with the maximum Z e (10 dBZ) at 5-7 km. A small number of strong radar echoes occur at midnight (0000-0100 LT) in the JL, with the maximum Z e of 10 dBZ.
A characteristic bimodality of the cloud base height and top height for single-layer, double-layer, and triple-layer clouds was observed on the ESTP. Importantly, the main peak (1-4 km) location of base height of the lowest layer clouds in the GZ is higher by 1 km than other sites, while the main peak (7-10 km) location of base height of the highest layer clouds in the JL is higher by 1 km than other sites. For cloud top heights, the main peak (4.5-5 km) of the GZ is located higher than in other sites. The secondary peak of top heights in the JL is located at~10 km, and is also apparently higher than in other sites. In addition, the ESTP is dominated by thin clouds with thicknesses of 200-400 m. Low clouds occur most frequently in the DC (46%), JL (39%), and LT (38%), respectively. The occurrence frequency of high clouds presents a larger value of 33% in the JL. In contrast, middle clouds dominate in the GZ with 38% occurrence frequency. The four-site observations show that cloud base height, top height, and thicknesses increase significantly in the afternoon due to the convection on the ESTP in summer. Higher top heights occur more frequently at midnight and early morning in the JL due to its mountain-valley terrain.