Climatology of Cloud Phase, Cloud Radiative Effects and Precipitation Properties over the Tibetan Plateau

: Current passive sensors fail to accurately identify cloud phase, thus largely limiting the quantiﬁcation of radiative contributions and precipitation of different cloud phases over the Tibet Plateau (TP), especially for the mixed-phase and supercooled water clouds. By combining the 4 years of (January 2007–December 2010) cloud phase (2B-CLDCLASS-LIDAR), radiative ﬂuxes (2B-FLXHR-LIDAR), and precipitation (2C-PRECIP-COLUMN) products from CloudSat, this study systematically quantiﬁes the radiative contribution of cloud phases and precipitation over the TP. Statistical results indicate that the ice cloud frequently occurs during the cold season, while mixed-phase cloud fraction is more frequent during the warm season. In addition, liquid clouds exhibit a weak seasonal variation, and the relative cloud fraction is very low, but supercooled water cloud has a larger cloud distribution (the value reaches about 0.24) than those of warm water clouds in the eastern part of the TP during the warm season. Within the atmosphere, the ice cloud has the largest radiative contribution during the cold season, the mixed-phase cloud is the second most important cloud phase for the cloud radiative contribution during the warm season, and supercooled water clouds’ contribution is particularly important during the cold season. In particular, the precipitation frequency over the TP is mainly dominated by the ice and mixed-phase clouds and is larger over the southeastern part of the TP during the warm season.


Introduction
Cloud cover plays an important role in climatic systems, having a significant effect on the radiation budget and corresponding water cycles [1][2][3][4]. It is therefore crucial to understand the main characteristics and physical processes of clouds, such as various macro-physical effects (e.g., cloud fraction (CF), cloud thickness) and microphysical properties (e.g., cloud droplet number concentration, particle size, phase), as well as the complicated dynamical [5,6] and microphysical processes [7]. However, due to an incomplete understanding of physical processes, processes related to clouds have been poorly represented in climate and weather models [8]; this lack of research has been identified as the greatest source of uncertainty in climate predictions driven by climate models [9,10].
Previous investigations on cloud properties have focused on the cloud phase, a key cloud parameter. The cloud radiative effect and precipitation properties have resulted in findings which vary due to cloud phase changes [11][12][13][14][15]. Clouds are composed entirely 2 of 21 of ice when temperatures are below −40 • C, or liquid particles when temperatures are above 0 • C (freezing temperature) [16,17]. In addition, the cloud may be composed of a combination of pure ice, liquid particles, or both (that is, mixed-phase) under this temperature range (−40-0 • C) [7,[18][19][20][21][22][23][24]. In particular, the liquid water cloud is considered a supercooled water cloud under a temperature lower than 0 • C [25]. As there are differences between liquid and ice particles in refractive indices, sizes, concentration, and shapes [26,27], clouds have a distinct radiative effect depending on their phases [28][29][30]. For example, a low-level water cloud (e.g., stratocumulus) can exert a strong negative net radiative effect at the surface [31,32]. In contrast, the thin cirrus causes a positive radiative effect of about 2 W/m 2 at the surface and atmosphere, especially in the central Pacific Tropics [33][34][35]. For the mixed-phase cloud, cloud radiative effect (CRE) is closely associated with ice-liquid partitioning [36]. In addition, observation-based studies have concluded that precipitation driven by warm water clouds is mainly distributed over tropical oceans [14]. It is therefore crucial for climate models to reasonably represent cloud phase partitioning to reduce uncertainty in climate predictions caused by the cloud phase. Furthermore, it is necessary for climate models to accurately evaluate contributions of different cloud thermodynamic phases to the radiation budget and precipitation processes. However, current climate models lack a level of detail in the simulations of mixed-phase and supercooled water clouds [37,38]. Consequently, their impacts on the climate system have not been well quantified at global or regional scales [22]. In particular, most climate models oversimplify cloud phase partitioning as a function of temperature, resulting in substantial biases in CRE models [19,39].
In this study, we focus on the Tibet Plateau (TP), also known as the "Asian water tower", an area that has experienced significant climate change over the past decades [40,41]. Across this region, significant climate warming has already been observed with glacier melting, the expansion of glacier-fed lake areas, and permafrost degradation [42,43]. Previous studies have indicated that the obvious climate change over the TP region is related to not only an increase in human activities and greenhouse gas emissions, but also to cloud radiation feedback [44][45][46]. For example, based on 6-hourly weather observations of 71 stations over the TP during 1961-2003, Duan and Wu [47] found that climate warming over the TP is possibly affected by changes in total cloud cover. Liu et al. [45] also concluded that the interaction between clouds and radiation is a primary factor in the warming of the TP. Yan et al. [48] also identified an evident seasonal variation in the cloud distribution over the TP, recording a bimodal structure from June to August, mainly due to the strong development of deep convection during the summer. The cloud-top height of deep convective clouds over the TP exhibited daily changes, leading to periodic changes in precipitation over the TP [49]. In addition, previous studies have pointed out that changes in the concentration of atmospheric aerosols and ozone may affect regional cloud and precipitation properties [50,51]. Based on model simulation, MacIntosh et al. [51] showed that changes in tropospheric and stratospheric ozone can influence the cloud cover and precipitation responses via altering the radiative forcing at the top of the atmosphere. Recent field campaigns have also been devoted to studying the cloud-aerosol-radiation interactions [52,53]. For example, the observations of aerosols above clouds and their interactions (ORACLES) project is focus on the biomass burning aerosol impact on the subtropical stratocumulus of the Southeast Atlantic Ocean [52]. Additionally, the Antarctic Circumnavigation Expedition: Pre-industrial aerosol climate effect (ACE-SPACE) project combined in situ measurements, satellite observation and models to investigate the aerosol-cloud-radiation interactions over the Southern Ocean [53]. Given a limited observation, however, few studies have focused on the cloud phase and its role in regulating the cloud radiative effect and precipitation feature over the TP region, and should be further addressed to accurately predict the future changes of TP climate.
To overcome the limitations of passive sensors in identifying cloud phase, this study uses the unique advantage of space-based lidar and radar: cloud-aerosol lidar with orthogonal polarization (CALIOP) and cloud profiling radar (CPR) [54], to evaluate the impacts of the cloud phase over the TP. In particular, we investigated how radiative contributions and precipitation frequency vary with cloud phase over the TP by combining the cloud phase (2B-CLDCLASS-LIDAR and ECMWF-AUX), radiative fluxes (2B-FLXHR-LIDAR), and precipitation (2C-PRECIP-COLUMN) products from CALIPSO and CloudSat. This paper is organized as follows: In Section 2, all the datasets, methods, and definitions of the cloud phase and CRE are briefly described. Section 3.1 provides distributions of the cloud fraction, CRE, cloud optical depth of different cloud phases over the TP. Section 3.2 describes their radiative contributions over the TP. Section 3.3 further discusses the corresponding net heating rate profiles, and Section 3.4 precipitation frequency over the TP. Lastly, the discussion and conclusions are given in Sections 4 and 5, respectively.

CALIOP Instrumentation Aboard CALIPSO
Space-borne lidar CALIOP incorporated on CALIPSO operates at two wavelengths (532 nm and 1064 nm) is sensitive to optically thin hydrometeor layers and aerosol layers [55][56][57]. As combining radar and lidar measurements may improve cloud parameter simulations (e.g., CF) in atmospheric models, it is conducive to reducing the uncertainties in weather forecasts and climate predictions [58,59].

The CPR Instrument Aboard CloudSat
As CloudSat CPR operated at 94 GHz, can penetrate the optically thick cloud layer, it is capable of observing vertical structures and measuring the vertical distribution of liquid water and ice water content within clouds [60,61]. It is worth noting that CloudSat's CPR cannot always profile optically thick clouds due to attenuation effects [62][63][64]. In addition, CloudSat products have a blind zone over the lowest atmospheric altitudes (c.750 m), leading to unexpected uncertainty in low-level clouds [65]. Further information relating satellites investigated in this study is presented in Table 1.

Satellite Products and Reanalysis Dataset
The 2B-CLDCLASS-LIDAR cloud classification product that utilizes the combination of lidar backscatter and radar reflectivity has the capability of classifying the types and phases of clouds (e.g., ice, liquid, and mixed-phase). The flag in the "CloudPhase" parameter of this product provides the information related to the cloud phases, i.e., "1" represents an ice cloud, "2" represents a mixed-phase cloud, and "3" represents a water cloud. In particular, this study further collocated the temperature profiles from the ECMWF-AUX product [66], which is an intermediate dataset that consists of the ancillary ECMWF state variables interpolated across each CloudSat CPR bin, to group the liquid phase clouds into warm water clouds and supercooled-water clouds based on whether cloud top temperatures are higher or lower than the threshold temperature (0 • C) [25]. Note that the super-cooled water cloud in the mixed-phase is not considered. Based on the cloud phase products, Matus and L'Ecuyer [30] quantified the cloud radiative effects for the ice, water, and mixed-phase clouds at a global scale. In addition to this product, the 2B-FLXHR-LIDAR product is also used in our study. By employing a broadband, two-stream, plane-parallel doubling-adding radiative transfer model [67], this product uses cloud and aerosol properties retrieved from CloudSat and CALIPSO as inputs to provide us with calculated radiative fluxes at the surface (SFC), top of the atmosphere (TOA), cloud optical depth, and atmospheric heating rate profiles [68]. Combining the CRE parameters (TOACRE and BOACRE in this product) and cloud phase information obtained from 2B-CLDCLASS-LIDAR, this study quantified the radiative contributions of different cloud phases to the total cloud radiative effect. There have been related studies on the uncertainty of cloud radiative effect and cloud heating rate, and by using collocated SW and LW flux observations from the Clouds and the Earth's Radiant Energy System (CERES), the 2B-FLXHR-LIDAR product is evaluated. Henderson et al. [69] pointed out that the global mean outgoing longwave radiation and outgoing shortwave radiation estimated from the collocated CERES observations and 2B-FLXHR-LIDAR calculations agree within 5 and 4 Wm −2 , respectively, and root-mean-square differences are 16 and 6 Wm −2 on monthly/5 • scales. The cloud radiative effect is usually defined as the radiative flux differences between clear-sky and all of the possible sky conditions at a given location (e.g., TOA) [30]. Here, it is noteworthy that the radiative contribution of each cloud phase is calculated only for single-layer clouds. Given two or more cloud layers existing simultaneously in the same location but at different atmospheric layers, it needs to be considered as a multilayered cloud (MLC) system [70,71]. Thus, this study accounted for the radiative contributions for five different cloud types: Single-layer ice clouds, single-layer warm water clouds, single-layer supercooled water clouds, single-layer mixed-phase clouds, and multilayered clouds. In addition to these products, we also use the precipitation parameter "Precip_flag" in the 2C-PRECIP-COLUMN dataset [72], to further quantify the precipitation frequency for different cloud phases over the TP. Particularly, flags 1 (rain possible), 2 (rain probable), and 3 (rain certain) in Precip_flag were considered to determine precipitation frequency in this study [73][74][75][76]. Other information is shown in Table 2.

Methodology
Here, the new version of 4 years (January 2007-December 2010) of joint products (e.g., 2C-PRECIP-COLUMN, 2B-CLDCLASS-LIDAR, and 2B-FLXHR-LIDAR) of CloudSat and CALIPSO are collected to perform our analysis. The vertical and horizontal resolutions of these products are 240 m and 1.1 km, respectively. In addition, the differences in cloud properties between warm (May to September) and cold (November to March) seasons were investigated in this study. Seasons were classified based on the work of Huang et al. [77], indicating that semi-arid areas have a rapid temperature increase during the cold season. As the TP is characterized by arid and semi-arid areas, rapid temperature increases are closely related to cold and warm seasons. The relative cloud fraction (RCF) of different cloud phases can be obtained as: where, RCF shows the relative cloud fraction of different cloud phases, and N is number of cloudy profiles. N cp are the number of cloud phase profiles. In particular, the RCF of a multilayered cloud is the ratio of the number of multilayered cloud profiles to the number of cloudy profiles. In this study, the averaged cloud radiative effect (CRE) in a given grid box or region was calculated as: where CRE i is the cloud radiative effect for the ith sample profile, and N is the number of sample profiles.
Cloud fractions under single-layer (CF s ) or multi-layer (CF m ) conditions can be obtained as: where, N and M are numbers of cloudy-sky and clear-sky profiles, respectively. N s and N m are the cloud phase samples of single-layer cloud and multi-layer cloud profiles, respectively. These two formulas are applicable to four different cloud phases (e.g., ice cloud). In particular, the cloud fraction of multilayered clouds is the ratio of the number of multilayered cloud profiles to the number of all-sky profiles. Then, the radiative contribution (Cont i ) for a given cloud type i in a given grid box or region may be derived from Equation (2): ) × 100 (5) where CRE i represents the averaged cloud radiative effect caused by a given cloud type (e.g., single-layer warm water cloud). In addition, the CH, which is the averaged net radiative heating rate in a given region, is expressed as follows: Here, HR i cloud and HR j clear are the ith and jth net radiative heating rate profiles for cloudy-sky and clear-sky conditions, respectively. N and M are numbers of cloudy-sky and clear-sky profiles, respectively. CF = N (N+M) represents the total cloud fraction in a given region. Thus, CH shows the weighted cloud net radiative heating rate profile.

The CF and CRE of Different Cloud Phase over the TP
Total cloud fraction and RCF for ice clouds, mixed-phase clouds, warm water clouds, supercooled water clouds, and multilayered clouds during the cold and warm seasons over the TP are shown in Figure 1a-c. The total cloud fraction at a given grid-box is defined as the ratio of the number of cloudy profiles to the number of all sample profiles, while the RCF for each cloud phase is calculated based on the Equation (1). Our results indicated that total clouds occurred more frequently during the warm season than during the cold season, especially in the southeastern part of the TP (Figure 1(a1-a3)). A low value (approximately 36%) of total cloud fraction over the southwestern part of the TP appeared during the cold season (see Figure 1(a1)) [78]. Ice clouds (Figure 1(b1,b3)) were recorded to have a high RCF value (even >84%) during the cold season, being recorded over most regions of the TP; these clouds had considerable seasonal variability, except for the southeastern area of the TP. Compared with ice clouds, a relatively small value (<20%) and weak seasonal variation of the RCF of warm water clouds were found over the TP (Figure 1(c1,c3)).

The CF and CRE of Different Cloud Phase over the TP
Total cloud fraction and RCF for ice clouds, mixed-phase clouds, warm water clouds, supercooled water clouds, and multilayered clouds during the cold and warm seasons over the TP are shown in Figure 1a-c. The total cloud fraction at a given grid-box is defined as the ratio of the number of cloudy profiles to the number of all sample profiles, while the RCF for each cloud phase is calculated based on the Equation (1). Our results indicated that total clouds occurred more frequently during the warm season than during the cold season, especially in the southeastern part of the TP (Figure 1(a1-a3)). A low value (approximately 36%) of total cloud fraction over the southwestern part of the TP appeared during the cold season (see Figure 1(a1)) [78]. Ice clouds (Figure 1(b1,b3)) were recorded to have a high RCF value (even >84%) during the cold season, being recorded over most regions of the TP; these clouds had considerable seasonal variability, except for the southeastern area of the TP. Compared with ice clouds, a relatively small value (<20%) and weak seasonal variation of the RCF of warm water clouds were found over the TP ( Figure  1(c1,c3)).

Figure 1.
Four-year total cloud fraction and relative cloud fraction (RCF) of different cloud phases over the Tibet Plateau (TP). The spatial resolution is a 0.5 • × 0.5 • grid (the two columns on the left represent the cold season, and the two columns on the right represent the warm season). a1 and a3 represent the total cloud fraction in the cold and warm seasons, respectively. Others represent the RCF of different cloud phases. Here, each cloud phase is a single-layer plus multiple layers; for example, the ice cloud is the total ice cloud.
Identifying supercooled water clouds or mixed-phase clouds from passive satellite sensors (e.g., MODIS) is challenging when multilayered clouds exist within the satellite observation field of view [79,80]. As CloudSat and CALIPSO have a vertically resolved ability, results indicated a weak seasonal variability in RCF of supercooled water clouds; RCF was still larger than those of warm water clouds, especially in the southeastern area of the TP during the warm season ( Figure 1(c2,c4)). Due to weak satellite observation constraints, mixed-phase clouds (a significant atmospheric component) are poorly simulated in climate models [81]. With advanced observational active sensors, our results indicate that mixed-phase clouds are widely distributed during the warm season compared with those during the cold season ( Figure 1(b2,b4)). In particular, mixed-phase clouds were mainly formulated at the southeastern area of the TP during the warm season when the RCF was approximately 50%. During the cold season, however, mixed-phase clouds recorded a relatively small RCF over the TP. Findings from our investigation also indicated that a large RCF for multilayered clouds during the warm season was distributed in the southeastern area of the TP, possibly being controlled by a strong ascending motion [71,82]. The sample of the cloud fraction of cloud phases is shown in Figure S1.
Furthermore, seasonal net cloud radiative effects ( Figure 2(a1-f4)) corresponding to all clouds, each cloud phase, and multilayered clouds at the TOA ( Figure 2(a1-c4)) and in the ATM (Figure 2(d1-f4)) over the TP were recorded. The cloud phases only refer to single-layer clouds, such as single-layer ice, warm water, mixed-phase, or supercooled water clouds. As the CRE at the SFC has a similar distribution with those at the TOA, it is not included in this study (see Figure S2). In addition, the CRE in the ATM is defined as the CRE difference between the TOA and the SFC. Total clouds at the TOA have a weak cooling effect (that is, a negative cloud radiative effect) over most regions of the TP during the cold season (from −10 Wm −2 to −30 Wm −2 ). Compared with the CRE during the cold season, the total cloud during the warm season exhibited a stronger cooling effect, especially over the southeastern part of the TP, where the CRE approximately reached −250 Wm −2 . The CRE over the north part of the TP maintained a moderate cooling effect during the warm season, ranging from −60 Wm −2 to −170 Wm −2 . Similar to the cloud fraction, the CRE of the ice cloud also showed evident seasonal variability. Ice clouds recorded a weak heating effect (around 6 Wm −2 ) over the northwestern part of the TP (Figure 2(b1)), and a moderate cooling effect in the southeastern part of the TP (<−30 Wm −2 ). However, during the warm season, a distinct cooling effect (exceed −50 Wm −2 ) appeared over most regions of the TP (see Figure 2(b3)). Although ice clouds tended to occur in the multilayered cloud system (MLS) [71], this study focused on the CRE of the single-layer cloud system. This finding may partly explain the reason that a stronger cooling effect of the ice cloud occurring during the warm season of the TP (Figure 2(b1,b3)) while the ice cloud occurred more frequently during the cold season, as all ice cloud samples (single-layer and multilayered ice clouds) were included in Figure 1(b1,b3). Moreover, the ice cloud optical depth, as well as the cloud fraction of single-layer ice cloud, also play a significant role [83,84]. As such, the ice cloud radiative effect depends largely on their optical and microphysical properties [84]. In Figure 2(b1), ice clouds have a weak heating effect over the northwestern part of the TP, while they have a moderate cooling effect in the southeastern part of the TP (about 0 -−30 Wm −2 ). Incorporating the ice cloud properties retrieved from CloudSat and CALIPSO into a radiative transfer model (termed libRadtran [83]), Hong et al. [83] also found that ice clouds have a global-averaged net warming effect (~5.1 ± 3.8 Wm −2 ) as their longwave warming effect (~21.8 ± 5.4 Wm −2 ) exceeding the shortwave cooling effect (~−16.7 ± 1.7 Wm −2 ). The net radiative effect of warming (or cooling) is still determined by the ice cloud optical depth ( Figure 3).    Four-year cloud radiative effect (CRE) (Wm −2 ) of different cloud phases over the TP. (a1-c4) represent the top of the atmosphere (TOA); (d1-f4) represent the ATM. The spatial resolution is a 0.5° × 0.5° grid (The two columns on the left represent the cold season and the two columns on the right represent the warm season. Here, ice clouds, liquid clouds, and mixed-phase clouds refer to the case of the single-layer clouds).

Figure 3.
Four-year short-wave cloud optical depth of four different cloud phases over the TP. (a1-a3) represent the cold season; (b1-b3) represent the warm season. The spatial resolution is a 0.5° × 0.5° grid (Here, ice clouds, liquid clouds refer to the case of the single-layer clouds. Part of the data for water clouds during the cold season is missing). In relation to mixed-phase clouds, previous studies have concluded that the CRE is highly dependent on ice-liquid partitioning [36]. However, due to technological limitations and complex physical processes, it is still difficult to accurately partition ice-liquid in the mixed-phase using passive sensors to provide reasonable outputs for global climate models [85,86]. Observation constraints from CloudSat and CALIPSO provide a unique opportunity to quantify the radiative effect of mixed-phase clouds at a global scale. For example, Matus and L'Ecuyer [30] found that mixed-phase clouds account for only one-tenth of the globally averaged total cloud fraction, but it exerts a global net CRE of −3.4 Wm −2 on annual average, and this effect is comparable with results of ice clouds or multilayered clouds. Over the TP, however, this study further indicated that the CRE of mixed-phase clouds play an important role more than those of other cloud phases during the warm season at the TOA. In addition, the CRE of the mixed-phase cloud showed apparent seasonal variability (Figure 2(b2,b4)), which may be caused by the variability in the frequency and optical depth varying with different seasons. During the warm season, the CRE of mixed-phase clouds reached up to −50 Wm −2 over most parts of the TP. During the cold season, however, the CRE of mixed-phase clouds became smaller with a range from −3 Wm −2 to 3 Wm −2 . Although liquid water clouds usually widely occur in tropical/subtropical oceans and dominate total CRE in the regions [30], Figure 2(c1,c2) show that warm water clouds and supercooled water clouds over the TP had a similar CRE value and distributions during the cold season except for those at the southeastern boundary of the TP where there is a stronger cooling effect (>−13 Wm −2 ) by supercooled water clouds rather than warm water clouds at the TOA. Our results also recorded a difference in CRE values and distributions over the TP for warm water and supercooled water clouds during the warm season. In particular, supercooled water clouds showed a cooling effect over the southeastern part of the TP by −30 Wm −2 , which is caused by a relatively larger occurrence frequency of the supercooled water cloud than those of the warm water cloud (see Figure 2(c2,c4)). To date, few studies have focused on the supercooled water clouds and their radiative effects. This study demonstrated that supercooled water clouds play a crucial role in regulating the radiative budget over the southeast part of the TP, especially during the warm season. For multilayered clouds, the seasonal variability in the CRE was similar to that in cloud fraction (Figure 2(a2,a4)).
CRE in the ATM represents the amount of radiative energy absorbed by the whole atmospheric column. From Figure 2(d1,d3), it is evident that total clouds during the warm season exerted a positive CRE that reached above 36 Wm −2 (Figure 2(d3)), indicating that the atmosphere is warming in the TP regions. On the contrary, the total clouds during the cold season showed a weak heating effect over most regions of the TP ( Figure  2(d1)). In relation to different cloud phases, the single-layer ice and the mixed-phase clouds exhibited prominent heating effects during the warm season. For the liquid clouds, the cloud radiative effects in the atmosphere were relatively smaller, resulting in a weak cooling (or heating) effect in the TP. Compared with warm water clouds, the cooling effect of supercooled water clouds was stronger during the cold season, especially at the southeastern boundary of the TP (reached up to −3 Wm −2 ). Meanwhile, supercooled water clouds also may slightly heat the atmosphere over most of the region of the TP during the warm season. During the cold season, the CRE of multilayered clouds approximately reached −6 Wm −2 , indicating a cooling effect in the atmosphere while reversely turning to a heating effect during the warm season (Figure 2(d4)).
Cloud optical depth results ( Figure 3) indicated that ice clouds have a larger cloud optical depth in the southwestern area of the TP during the warm season, corresponding to a stronger cooling effect at the TOA. During the cold season, the cloud optical thickness reached approximately 4 in the southeastern area of TP, corresponding to a relatively stronger cooling effect than other regions. In particular, the sample size for calculating the cloud optical depth is the sample size of each cloud phase in a single-layer. By incorporating ice cloud properties retrieved from CloudSat and CALIPSO into a radiative transfer model (termed libRadtran [83]), Hong et al. [83] showed that radiative effects induced a transition from heating to cooling (or vice versa) at the turning point of the ice cloud optical depth threshold, indicating that ice cloud radiative effect is closely related to ice cloud optical depth. Warm water and supercooled water clouds also recorded similar results, thereby sufficiently illustrating the strong relationship between cloud optical depth and cloud radiative effects. In addition, Gordon et al. [87] showed that the cloud optical depth of the relatively cold low cloud shows a positive correlation with temperature, while the cloud optical depth of the relatively warm low cloud shows a negative correlation with temperature. The results of Figure 3 show that the optical depth of supercooled liquid water clouds during the warm season is greater than that of the cold season, while the optical depth of warm water clouds during the warm season is very small. This result is similar to the above. On the whole, the optical depth of water and ice clouds during the warm season is greater than that during the cold season [88,89]. Figure 4 shows the regionally averaged cloud fractions of the total and different cloud phases over the TP. As a comparison, the statistical results are presented during the cold and warm seasons. Here, the statistical results of the cloud fraction are represented by two different conditions (see Figure 4). In addition, cloud fractions of different cloud phases under the single-layer and multi-layer conditions in Figure 4 were based on the Equations (3) and (4), respectively. Note that in the case of a multi-layer cloud profile with multiple cloud phases, the different cloud phase profiles overlap. For example, if there are both ice clouds and mixed-phase clouds in a profile, we will add 1 to both ice cloud samples and mixed phase cloud samples. Thus, in the case of multi-layer condition, the sum of cloud fraction of different cloud phases may be larger than 1.

Regional Averaged Cloud Fraction and Radiative Contribution
Our results indicate that fractions of total and multilayered clouds during the warm season were higher than those during the cold season ( Figure 4). This finding is consistent with the results of Kukulies et al. [78], indicating that the occurrence frequencies of total cloud over the TP during the warm and cold seasons were nearly 80% and 60%, respectively. Among the cloud phases investigated in this study, ice cloud dominated total clouds during all seasons for both the single-and multi-layered cases. The single-layer mixed and liquid phase clouds exhibited a larger frequency occurrence during the warm season than during the cold season. In addition, a multilayered cloud system showed a tendency to frequently occur during the warm season since multilayered clouds may be induced by stronger upward air motions [82]; single-layered ice clouds frequently occur during the cold season, and they easily overlap with other cloud types during the warm season [71]. by two different conditions (see Figure 4). In addition, cloud fractions of different cloud phases under the single-layer and multi-layer conditions in Figure 4 were based on the Equations (3) and (4), respectively. Note that in the case of a multi-layer cloud profile with multiple cloud phases, the different cloud phase profiles overlap. For example, if there are both ice clouds and mixed-phase clouds in a profile, we will add 1 to both ice cloud samples and mixed phase cloud samples. Thus, in the case of multi-layer condition, the sum of cloud fraction of different cloud phases may be larger than 1.

Figure 4.
Four-year cloud fraction of different cloud phases over the TP (The first four columns represent the ice cloud, the mixed-phase cloud, the warm water cloud, and the supercooled water cloud, respectively, and represent only a single-layer, while the last four columns represent multilayer condition. And the multilayer condition further classifies the multilayered clouds into four different cloud phases if the multilayered clouds include the given cloud phases. In addition, the Four-year cloud fraction of different cloud phases over the TP (The first four columns represent the ice cloud, the mixed-phase cloud, the warm water cloud, and the supercooled water cloud, respectively, and represent only a single-layer, while the last four columns represent multilayer condition. And the multilayer condition further classifies the multilayered clouds into four different cloud phases if the multilayered clouds include the given cloud phases. In addition, the left column has total clouds and multilayered clouds, (a) represents the cold season, and (b) represents the warm season).
The fraction of the liquid water clouds also recoded a significant increase under multilayered cloud conditions during the warm season ( Figure 4). Compared with singlelayer conditions, it is also noteworthy that the proportions of liquid cloud and mixedphase clouds are similar during all seasons under multilayered cloud conditions [78]. Furthermore, our result further demonstrates that the fraction of supercooled water clouds was higher than that of the warm water clouds over the TP.
A regional average of CRE under different cloud phases was calculated over the TOA, SFC, and ATM in the TP ( Figure 5). Total clouds at both the TOA and the SFC exhibited a cooling effect, recording a trend of increase in the CRE with a change from the cold season to the warm season (refer to Figure S3 for detailed values of the CRE). The strongest cooling effect of the total cloud occurred during the warm season (WS), approximately −116.2 Wm −2 and −94.8 Wm −2 at the SFC and TOA, respectively. As shown in Figure 2, clouds during the warm season usually have a stronger impact on the radiative budget. Indeed, as highlighted in Figure 5, CRE of total clouds during the warm season was approximately triple or quadruple than that of the cold season. In addition, CRE of the multilayered cloud during the warm season was approximately quadruple that during the cold season, except for in the ATM.
Our results indicate that differences in the CRE between the cold and warm seasons in the ATM was evident. It is also found that warm water and supercooled water clouds during the cold season have a weak cooling effect in the ATM, approximately −0.3 Wm −2 and −0.6 Wm −2 , respectively. However, clouds generally heat the ATM during the warm season. Results in Figure 5 indicate that the stronger cloud heating effect in the ATM during the warm season is mainly caused by the ice, mixed-phase, and multilayered clouds, whereas liquid clouds contribute less heating or cooling effects. Besides, the cloud cooling effects at the TOA and SFC for different cloud phases were higher during the warm season compared with that during the cold season.
A regional average of CRE under different cloud phases was calculated over the TOA, SFC, and ATM in the TP ( Figure 5). Total clouds at both the TOA and the SFC exhibited a cooling effect, recording a trend of increase in the CRE with a change from the cold season to the warm season (refer to Figure S3 for detailed values of the CRE). The strongest cooling effect of the total cloud occurred during the warm season (WS), approximately −116.2 Wm −2 and −94.8 Wm −2 at the SFC and TOA, respectively. As shown in Figure 2, clouds during the warm season usually have a stronger impact on the radiative budget. Indeed, as highlighted in Figure 5, CRE of total clouds during the warm season was approximately triple or quadruple than that of the cold season. In addition, CRE of the multilayered cloud during the warm season was approximately quadruple that during the cold season, except for in the ATM. Our results indicate that differences in the CRE between the cold and warm seasons in the ATM was evident. It is also found that warm water and supercooled water clouds during the cold season have a weak cooling effect in the ATM, approximately −0.3 Wm −2 and −0.6 Wm −2 , respectively. However, clouds generally heat the ATM during the warm season. Results in Figure 5 indicate that the stronger cloud heating effect in the ATM during the warm season is mainly caused by the ice, mixed-phase, and multilayered clouds, Results for different cloud phases indicate that mixed-phase clouds have a stronger cooling effect at the TOA and SFC, with ice clouds possibly resulting in a stronger heating effect in the ATM. Although the cooling effect for the liquid clouds was relatively weak, supercooled water clouds recorded a stronger cooling effect in the TOA and SFC than that of the warm water cloud. Overall, our results indicate a considerable seasonal variability at the SFC compared with the TOA.
Results in Figure 6 further illustrate radiation contributions associated to different cloud phases to total cloud radiative effects based on Equation (2) (refer to Figure S4 for more details). At the TOA, multilayered clouds and mixed-phase clouds almost accounted for >50% of the contribution, especially during the warm season, whose contribution even exceeds 60% (68.1% max value). Compared with ice or liquid cloud phase, the mixed-phase cloud has a larger contribution during the warm season, resulting in a more crucial role in regulating the energy balance of the TP region. In addition, the radiation contribution of ice cloud was larger during the cold season than that during the warm season, especially in the ATM where its contribution peaked at 80%. Supercooled water clouds also recorded a considerable contribution at the TOA (18.2%) and SFC (12.5%) during the cold season, with the contribution even exceeds those results from the warm water clouds. This result indicates that supercooled water clouds play an important role in affecting the energy balance of the atmosphere. whereas liquid clouds contribute less heating or cooling effects. Besides, the cloud cooling effects at the TOA and SFC for different cloud phases were higher during the warm season compared with that during the cold season.
Results for different cloud phases indicate that mixed-phase clouds have a stronger cooling effect at the TOA and SFC, with ice clouds possibly resulting in a stronger heating effect in the ATM. Although the cooling effect for the liquid clouds was relatively weak, supercooled water clouds recorded a stronger cooling effect in the TOA and SFC than that of the warm water cloud. Overall, our results indicate a considerable seasonal variability at the SFC compared with the TOA.
Results in Figure 6 further illustrate radiation contributions associated to different cloud phases to total cloud radiative effects based on Equation (2) (refer to Figure S4 for more details). At the TOA, multilayered clouds and mixed-phase clouds almost accounted for >50% of the contribution, especially during the warm season, whose contribution even exceeds 60% (68.1% max value). Compared with ice or liquid cloud phase, the mixedphase cloud has a larger contribution during the warm season, resulting in a more crucial role in regulating the energy balance of the TP region. In addition, the radiation contribution of ice cloud was larger during the cold season than that during the warm season, especially in the ATM where its contribution peaked at 80%. Supercooled water clouds also recorded a considerable contribution at the TOA (18.2%) and SFC (12.5%) during the cold season, with the contribution even exceeds those results from the warm water clouds. This result indicates that supercooled water clouds play an important role in affecting the energy balance of the atmosphere.

Heating Rate Profiles of the Cloud Phase over the TP
Previous investigations have focused on the cloud heating rate over the specific regions, cloud phases, or specific cloud types [90][91][92][93]. However, due to observation limitations, few studies have evaluated the heating rate for various cloud phases over the TP Here, ice clouds, water clouds, and mixed-phase clouds refer to the case of single-layer clouds).

Heating Rate Profiles of the Cloud Phase over the TP
Previous investigations have focused on the cloud heating rate over the specific regions, cloud phases, or specific cloud types [90][91][92][93]. However, due to observation limitations, few studies have evaluated the heating rate for various cloud phases over the TP region. Vertical profiles for the net heating rate caused by the four single-layer cloud phases over different seasons are shown in Figure 7. As heating rates were obtained from the cloudy and clear sky by weighting the cloud fraction based on Equation (6), differences in heating rates represent the impact of different cloud phases on the temperature profile. Results for ice clouds indicated a strong heating effect at low-and high-level atmospheres (e.g., approximately 2-9 km and 18-24 km), attaining 1.3 K/day during the warm season. However, it was also shown (Figure 7a) that ice clouds located between 9-18 km above the TP had a strong cooling effect. Compared with the warm season, ice clouds during the cold season more easily cooled the atmosphere, and peak heating rates were approximately recorded at 12 km, attaining −0.5 K/day. Similar vertical distribution heating rates at about 30 • N for ice clouds were recorded by Dolinar et al. [94] using the 2C-ICE product over a four year period (2007 to 2010). Results for single-layer mixed-phase clouds (Figure 7b) indicated that mixed-phase clouds recorded a relatively smaller heating rate compared with ice clouds. The relatively weaker heating effect (~1.1 K/day) corresponds to a smaller CRE within the atmosphere ( Figure 5). Similar to ice clouds, mixed-phase clouds during the warm season also recorded a cooling tendency on the middle-level atmosphere (9-16 km), having a stronger effect (up to −0.4 K/day). On the contrary, the mixed-phase cloud heating rate during the cold season recorded a weak heating effect at low altitudes. Compared with the ice clouds and mixed-phase clouds, the heating rate of the liquid clouds was smaller at the high-level atmosphere due to lower occurrence frequencies. In addition, warm water clouds at low altitudes also recorded a very weak heating effect (<6 km), especially during the cold season (heating rate <0.1 K/day). Compared with warm water clouds, supercooled water clouds recorded a stronger heating effect with altitudes lower than 6km (~0.4 K/day), and a stronger cooling effect in the mid-level atmosphere (~−0.15 K/day) during the warm season. Overall, differences in the heating rate profiles between different cloud phases during different seasons were mainly controlled by the vertical distribution of the cloud phase and their microphysical, macrophysical, and optical properties. As the important components of mixed-phase clouds, altocumulus and altostratus, for example, have different cloud radiative effects and heating rates due to differences in the dominant cloud phases [95].
Contributions of different cloud phases to cloud net heating rate at different heights are shown in Figure 8 (as well as net heating rate contribution of multilayered clouds). Results for single-layer mixed-phase clouds (Figure 7b) indicated that mixed-phase clouds recorded a relatively smaller heating rate compared with ice clouds. The relatively weaker heating effect (~1.1 K/day) corresponds to a smaller CRE within the atmosphere ( Figure 5). Similar to ice clouds, mixed-phase clouds during the warm season also recorded a cooling tendency on the middle-level atmosphere (9-16 km), having a stronger effect (up to −0.4 K/day). On the contrary, the mixed-phase cloud heating rate during the cold season recorded a weak heating effect at low altitudes. Compared with the ice clouds and mixed-phase clouds, the heating rate of the liquid clouds was smaller at the high-level atmosphere due to lower occurrence frequencies. In addition, warm water clouds at low altitudes also recorded a very weak heating effect (<6 km), especially during the cold season (heating rate <0.1 K/day). Compared with warm water clouds, supercooled water clouds recorded a stronger heating effect with altitudes lower than 6km (~0.4 K/day), and a stronger cooling effect in the mid-level atmosphere (~−0.15 K/day) during the warm season. Overall, differences in the heating rate profiles between different cloud phases during different seasons were mainly controlled by the vertical distribution of the cloud phase and their microphysical, macrophysical, and optical properties. As the important components of mixed-phase clouds, altocumulus and altostratus, for example, have different cloud radiative effects and heating rates due to differences in the dominant cloud phases [95].
Contributions of different cloud phases to cloud net heating rate at different heights are shown in Figure 8 (as well as net heating rate contribution of multilayered clouds). Here, heating rate contribution, calculated using Equation (5), is similar to that of CRE. These results clearly show that the contribution for each cloud phase varies with height and season. Overall, ice clouds almost recorded the largest contribution among all cloud phases, notably exceeding 70% during the cold season. During the warm season, MLS dominated the heating rate for the majority of heights. By using the same dataset to quantify cloud radiative effects of the MLS, LÜ et al. [96] recorded that, by averaging global overlapped cloud layers, a tendency of heating low-level atmospheres more easily than single-layer cloud systems was recorded. However, the dominant contribution of MLS at in the low-level atmosphere was not identified, possibly being related to different overlap cloud types and overlap percentages over the TP region [71]. Overall, the contribution from mixed-phase clouds was larger during the warm season than during the cold season. The contribution of supercooled water clouds exceeded those of warm water cloud, being more important during the warm season. The contribution of supercooled water clouds exceeded those of warm water cloud, being more important during the warm season.

Precipitation Frequency of Cloud Phase over the TP
Different cloud phases have a different precipitation regimes characteristic at a global scale [15]. However, due to the observational limitations, few studies have focused on the effects of different cloud phases on the precipitation over the TP. This study further researched the precipitation frequency of different cloud phases over the TP during the cold and warm seasons. Based on the definition of the cloud phase precipitation produced by the bottom cloud phase [14], this study defined the precipitation frequency as the profile of a certain bottom cloud phase with precipitation divided by the cloudy-sky profiles. Figures 9 and 10 illustrate that the ice and the mixed-phase clouds produced more precipitation, while the precipitation frequency produced by the liquid water clouds was very low, which is similar to the results of recent previous investigations [14,15]. searched the precipitation frequency of different cloud phases over the TP during the cold and warm seasons. Based on the definition of the cloud phase precipitation produced by the bottom cloud phase [14], this study defined the precipitation frequency as the profile of a certain bottom cloud phase with precipitation divided by the cloudy-sky profiles. Figures 9 and 10 illustrate that the ice and the mixed-phase clouds produced more precipitation, while the precipitation frequency produced by the liquid water clouds was very low, which is similar to the results of recent previous investigations [14,15].  During the cold season, the precipitation caused by the ice clouds mainly occurred in the northeastern (mainly concentrated in Qilian Mountains) and southeastern part of effects of different cloud phases on the precipitation over the TP. This study further r searched the precipitation frequency of different cloud phases over the TP during the col and warm seasons. Based on the definition of the cloud phase precipitation produced b the bottom cloud phase [14], this study defined the precipitation frequency as the profi of a certain bottom cloud phase with precipitation divided by the cloudy-sky profile Figures 9 and 10 illustrate that the ice and the mixed-phase clouds produced more precip itation, while the precipitation frequency produced by the liquid water clouds was ver low, which is similar to the results of recent previous investigations [14,15].  During the cold season, the precipitation caused by the ice clouds mainly occurre in the northeastern (mainly concentrated in Qilian Mountains) and southeastern part o During the cold season, the precipitation caused by the ice clouds mainly occurred in the northeastern (mainly concentrated in Qilian Mountains) and southeastern part of the TP (approximately 0.01-0.02) (Figure 9). In particular, precipitation from ice clouds also frequently occurred in the western tail of the TP. Kapnick et al. [97] showed that precipitation in Karakoram Mountain occurred mainly from December to May of the following year, which may be related to the high precipitation frequency from ice clouds in the western tail boundary of the TP. Compared with the cold season, the precipitation frequency of the ice cloud during the warm season occurred throughout the TP, while the peak precipitation (reached up to 0.12) occurred in the middle and west of the TP. The central area of the TP is mainly the Qiangtang Plateau where there are plateau lakes with the largest number of lakes in the world and the highest lake level. The Qiangtang Plateau also provides abundant water vapor conditions to trigger precipitable conditions in the summer, possibly being one reason why ice clouds in the central TP had a greater precipitation frequency during the warm season. Dong et al. [98] showed that summer precipitation in the southwestern area of the TP was mainly caused by convective storms in the Indian subcontinent, possibly being a significant factor for the distinct difference in precipitation from ice clouds between the cold and warm seasons of over the southwestern part of the TP. In particular, precipitation frequency was very high at the big bend of the Yarlung Zangbo River, possibly associated to a large volume of water vapor transported from the Bay of Bengal.
For the mixed-phase cloud (Figure 9(a2) and Figure 10(a2)), the mixed-phase cloud has a relatively smaller precipitation frequency compared with that of the ice cloud during the cold season. The large precipitation frequency of the mixed-phase cloud mainly occurred at the eastern boundary of the TP and the Yarlung Zangbo River, exceeding 0.02 during the cold season. During the warm season, the mixed-phase clouds dominantly occupied the precipitation over the eastern part of the TP, and the precipitation frequency was significantly greater than that of mixed-phase clouds during the cold season, exceeding 12%. We can see that the RCF of mixed phase clouds in Figure 1 also shows the distribution characteristics corresponding to its precipitation frequency. The warm-moist air convergence from the Indian Ocean is usually transported into the TP as a heat source during the summer, inducing deep convective cloud systems mainly in the eastern part of the TP [48,78,86]. Indeed, Yang et al. [99] also derived similar conclusions using precipitation and δ18O data. Different from the ice cloud, the precipitation frequency of the mixed-phase clouds in the west and middle part of the TP was less than the ice cloud's precipitation frequency, which is about 0.06-0.09 during the warm season. Overall, the high-cloud precipitation over the TP was mainly caused by the ice cloud and mixed-phase clouds. Besides, the precipitation mainly occurred at the southeastern part of TP during the warm season [100].
It is also notable that seasonal variability of precipitation frequency is predominantly dominated by ice and mixed-phase clouds during the cold season. While precipitation frequency over the southeastern area of the TP was mainly dominated by mixed-phase clouds during the warm season. Precipitation frequency of ice clouds is mainly distributed in the middle part of the TP during the warm season. Compared with ice and mixed-phase clouds, the frequency of liquid clouds inducing precipitation was generally smaller over the TP due to a lower frequency. In addition, warm water and supercooled water clouds at the south-central border of the TP have a very lower precipitation frequency (<0.02), especially during the warm season.

Discussion
Based on the statistic results of Figure 1, we found that clouds are more distributed during the warm season than the cold season, and this finding may be due to a plentiful supply of moisture during the warm season [101]. In addition, based on long-term cloud records from NASA Clouds and CERES, Naud et al. [102] recorded a higher cloud fraction in the northern part of the TP in the winter and a lower cloud fraction in the summer, findings which are consistent with our results. Our results also indicated that cold rain caused by ice clouds dominated precipitation frequency over the TP during the cold season [14] (see Figure 9); Choi et al. [19] highlighted that an increase in the dust frequency during the cold season (e.g., winter and spring seasons) over mid-latitude regions of the Northern Hemisphere is derived from a decrease in the supercooled water cloud fraction that effectively glaciates supercooled water clouds by lifted dust aerosols. Additionally, based on the station observation data, an increase in dust and carbon contents (black, organic, and elements) have been detected during non-monsoon periods [103,104]. Spatial and seasonal distributions of mixed-phase clouds are closely related to those of deep convective clouds [105,106]. From the results of CRE and heating rate, we can see that cloud radiative effects and heating rates of cloud phases show different distribution characteristics, which is probably linked with aerosol loading and related dynamic processes [19]. Based on the observation dataset, due to low ice nuclei aerosol (e.g., dust) loading [7,19,21], there are more distributions of super-cooled water clouds over the Southern Ocean. By incorporating CERES-CloudSat-CALIPSO-MODIS (CCCM) datasets into a radiative transfer model, Bodas-Salcedo et al. [29] found that clouds with supercooled liquid tops contribute 27%-38% to the total reflected solar radiation over the Southern Ocean, whereas most climate models were not able to reproduce the importance of the supercooled water cloud on the radiative budget over the Southern Ocean. For the TP regions, the super-cooled water cloud is more than the warm water cloud, especially in the southeastern part of the TP. In addition, The spatial and seasonal variability in the cloud fraction over the TP (Figure 4) may be induced by water vapor and aerosol loading conditions [50,102]. Besides, more ice clouds during the cold season may be linked with atmospheric temperature and microphysical properties (e.g., ice nuclei aerosol) [21,102].
Our statistic results also showed that the precipitation frequency of ice and mixedphase clouds is more than that of the water clouds (see Figures 9 and 10). Heymsfield et al. [15] also demonstrated that ice clouds produce the largest proportion of precipitation, while warm water clouds produce the least precipitation on a global scale. In addition, precipitation over the Qilian Mountains during the cold season, which are located at the northeastern boundary of the TP, may be induced by the southward movement of the westerly belt. However, the precipitation on the northwest slope in the Qilian Mountains is mainly triggered by the water vapor from the westerly belt.

Conclusions
In order to understand the global energy budget and water cycles, it is important to fully understand temporal and spatial distributions of cloud phases and their corresponding radiation and precipitation properties. However, conventional passive sensors (e.g., MODIS) are not able to provide accurate cloud phase information and vertical profiles, especially for mixed-phase and supercooled water clouds [107]. Furthermore, current climate models lack partitioning skills and they are unable to simulate mixed-phase and supercooled water clouds, resulting in considerable outputs uncertainty [36,108]. In addition, few studies have evaluated the impact of different cloud phases on the Earth's radiation budget and their precipitation frequencies at regional or global scales (e.g., TP). By using space-based lidar and radar (e.g., CALIOP and CPR) methods, we investigated the spatial distributions of various cloud phases, assessing their impacts on radiative contributions and precipitation frequencies over the TP. Although some statistical results agree well with those of previous studies, new insights are also gained. New insights are gained specifically for the TP. This paper mainly show following results: Total cloud cover recorded a decreasing trend from the southeast to the northwest, possibly due to topographical characteristics across the TP [50]. In addition, our results also indicated that ice clouds frequently occurred during the cold season over the TP and the RCF was maximal (~0.84) among the cloud phases. RCFs of mixed-phase and supercooled water clouds were higher during the warm season. Our results also found that ice clouds have a larger radiative contribution during the cold season. It should be noted that liquid clouds have an atmospheric cooling effect during the cold season over the TP. Furthermore, our results demonstrated that heating rate profiles obtained from different cloud phases varied with different seasons. For example, ice clouds recorded a tendency to heat low-and high-level atmospheres during both cold and warm seasons, while warm water clouds heated the atmosphere near the surface. The heating effect of supercooled water clouds was not negligible during the warm season, especially over the low-level atmosphere. The cloud optical depth and cloud stability have a certain influence on the cloud phase. This study also examined the impacts of the cloud phases on precipitation frequency. The results showed that the precipitation frequency over the TP mainly occurred during the warm season and the precipitation occurrence was dominated by the ice and mixed-phase clouds. Maussion et al. [109] also concluded that precipitation occurred more frequently during summer. In addition, the precipitation dominated by the ice clouds was mainly concentrated in the central part of the TP, while the precipitation dominated by the mixed-phase clouds was mainly produced in the southeastern part of the TP during the warm season.
Previous studies have shown that cloud phase feedback causes more solar radiation to be reflected into space [86,110], with recent studies mainly focusing on the radiative effects of the ice, warm, and mixed-phase clouds [30]. However, our results indicated that supercooled water clouds have an important radiative impact over the TP. In particular, cloud phase partitioning in mixed-phase clouds requires further investigation to advance cloud phase parameterization schemes. It is therefore important to examine the effects of aerosols and dynamic factors on cloud phase changes, not only to accurately quantify cloud radiative impacts of various cloud phases, but to reduce biases of climate feedbacks [7,21,81].
Supplementary Materials: The following are available online at https://www.mdpi.com/2072-429 2/13/3/363/s1, Figure S1: Four-year sample of cloud fraction of different cloud phases over the TP, Figure S2: Four-year CRE of different cloud phases over the TP, Figure S3: Four-year CRE of different cloud phases at TOA, ATM, and SFC over the TP, Figure S4: Four-year cloud radiation contribution of different cloud phases at TOA, ATM, and SFC over the TP.