Evapotranspiration and Its Partitioning in Alpine Meadow of Three-River Source Region on the Qinghai-Tibetan Plateau

: The Qinghai-Tibetan Plateau (QTP) is generally considered to be the water source region for its surrounding lowlands. However, there have only been a few studies that have focused on quantifying alpine meadow evapotranspiration ( ET ) and its partitioning, which are important components of water balance. This paper used the Shuttleworth–Wallace (S–W) model to quantify soil evaporation ( E ) and plant transpiration ( T ) in a degraded alpine meadow (34 ◦ 24 (cid:48) N, 100 ◦ 24 (cid:48) E, 3963 m a.s.l) located at the QTP from September 2006 to December 2008. The results showed that the annual ET estimated by the S–W model ( ET SW ) was 511.5 mm (2007) and 499.8 mm (2008), while E estimated by the model ( E SW ) was 306.0 mm and 281.7 mm for 2007 and 2008, respectively, which was 49% and 29% higher than plant transpiration ( T SW ). Model analysis showed that ET , E , and T were mainly dominated by net radiation ( R n ), while leaf area index ( LAI ) and soil water content at a 5 cm depth ( SWC 5cm ) were the most important factors inﬂuencing ET partitioning. The study results suggest that meadow degradation may increase water loss through increasing E , and reduce the water conservation capability of the alpine meadow ecosystem.


Introduction
Hydrological processes of terrestrial ecosystems play an important role in interactions between the different spheres of the Earth (hydrosphere, biosphere, atmosphere, and geosphere), which mainly includes precipitation, evapotranspiration (ET), surface runoff, and drainage. Among them, ET is the main component of water loss from terrestrial ecosystems to the atmosphere [1,2]. ET is controlled by many environmental and biological factors; in turn, it affects not only plant growth and development but also the microclimate of plant communities. In addition, ET plays a major role in regional and global climate change [3] because it links closely to the latent heat energy, carbon, and water cycles in terrestrial ecosystems. Therefore, there has been great interest in the study of ET to better understand the links between ET and other Earth system processes [4].
ET is the combination of transpiration from vegetation (T) and evaporation from the soil surface (E), where ET partitioning is a subject of ongoing research due to the complexity of surface energy balance processes and measurements [5][6][7], and is very important in predicting the responses of ecosystem water balance to climate and vegetation coverage changes [8]. The most prevalent approach for measuring ecosystem ET is eddy

Study Site Description
The study was conducted in a degraded meadow (34 • 24 N, 100 • 24 E, 3963 m a.s.l) in Guoluo Prefecture in Qinghai Province, China, which is located at the Three-River Source Region (TRSR) of the Qinghai-Tibetan Plateau (QTP). The local climate is a typical plateau continental climate with long cold winters and short cool summers. Based on the data from 1995 to 2004, the monthly mean air temperature ranges from −12.3 • C (January) to 10.1 • C (July), with the annual mean temperature between −1.4-0.7 • C; the annual precipitation ranges between 381 and 551 mm with the mean value of 500 mm, of which 80% fell in the growing season from May to September; and the mean annual sunshine time is above 2500 h, with the annual total solar radiation ranging from 6238 to 6299 MJ·m −2 . The degraded meadow is dominated by Aconitum pendulum, Ligularia virgaurea, Pedicularis kansuensis, Oxytropis ochrantha, Ajania tenuifolia, Polygonum sibiricum, Euphorbia fischeriana, and Morina chinensis, with a mean vegetation height of less than 5 cm and a maximum canopy cover of 55% during the growing season. The soil of the study site is classified as Humic Cambisols [25].

Observation Method
The open-path eddy covariance system was installed in a flat degraded meadow at 3 m above the ground, with a fetch of more than 300 m from all directions. A threedimensional sonic anemometer (CSAT3, CSI, Logan, UT, USA) was used to measure turbulence. Variation of water vapor density was measured with the open-path CO 2 /H 2 O analyzer set at 10 Hz (Li-7500, Li-Cor, Lincoln, NE, USA). All the instruments were mounted on an observation tower of 3 m above the ground. Meanwhile, a micro-meteorological system was used to measure environmental variables including wind direction and velocity (014A and 034A-L, CSI, Logan, UT, USA), net radiation (CNR-1, Kipp&Zonen, Delft, The Netherlands), soil heat flux (HFT-3, CSI, Logan, UT, USA), air temperature and humidity (HMP45C, CSI, Logan, UT, USA), soil temperature at different depths (105T, CSI, Logan, UT, USA), precipitation (TE525MM, CSI, Logan, UT, USA), soil water content at different depths (TDR, CS615, CSI, Logan, UT, USA), and other related data. All data were recorded by using dataloggers (CR5000 and CR23X, CSI, Logan, UT, USA) at 15-min intervals. The study period was from 16 September 2006, to 31 December 2008. Data gaps were filled by linear interpolation using the preceding and following data when the gap was in the nighttime, daytime gaps were filled by the relationship between solar radiation and measured H or LE [18,34].
During the growing season, leaf area index (LAI) was determined using a leaf area meter (LI-3100, Li-Cor) about once a month, where fresh leaves were cut for five quadrats of 0.25 m × 0.25 m and the averaged LAI for the five quadrats was used in this study. The seasonal variations of LAI in 2007 and 2008 are shown in Figure 1. suensis, Oxytropis ochrantha, Ajania tenuifolia, Polygonum sibiricum, Euphorbia fisch Morina chinensis, with a mean vegetation height of less than 5 cm and a maxim cover of 55% during the growing season. The soil of the study site is classified Cambisols [25].

Observation Method
The open-path eddy covariance system was installed in a flat degraded m m above the ground, with a fetch of more than 300 m from all directions. A th sional sonic anemometer (CSAT3, CSI, Logan, UT, USA) was used to measure Variation of water vapor density was measured with the open-path CO2/H2O a at 10 Hz (Li-7500, Li-Cor, Lincoln, NE, USA). All the instruments were mou observation tower of 3 m above the ground. Meanwhile, a micro-meteorolog was used to measure environmental variables including wind direction an (014A and 034A-L, CSI, Logan, UT, USA), net radiation (CNR-1, Kipp&Zonen Netherlands), soil heat flux (HFT-3, CSI, Logan, UT, USA), air temperature an (HMP45C, CSI, Logan, UT, USA), soil temperature at different depths (105T, C UT, USA), precipitation (TE525MM, CSI, Logan, UT, USA), soil water content depths (TDR, CS615, CSI, Logan, UT, USA), and other related data. All data we by using dataloggers (CR5000 and CR23X, CSI, Logan, UT, USA) at 15-min int study period was from September 16, 2006, to December 31, 2008. Data gaps by linear interpolation using the preceding and following data when the gap nighttime, daytime gaps were filled by the relationship between solar radiation ured H or LE [18,34].
During the growing season, leaf area index (LAI) was determined using meter (LI-3100, Li-Cor) about once a month, where fresh leaves were cut for fiv of 0.25 m × 0.25 m and the averaged LAI for the five quadrats was used in this seasonal variations of LAI in 2007 and 2008 are shown in Figure 1.

Modeling
The S-W model was used to estimate E and T in this study, and the fo follows:

Modeling
The S-W model was used to estimate E and T in this study, and the formula is as follows: where ET SW , E SW , and T SW are the calculated ET, E, and T by the S-W model and PM s and PM c are terms used to describe E and T, respectively. C s and C c are the soil surface resistance coefficient and canopy resistance coefficient, respectively. PM s and PM c are calculated as follows: PM c = ∆R + (aeC p D − ∆r ac R s )/(r aa + r ac ) ∆ + fl(1 + (r sc /(r aa + r ac ))) where ∆ is the slope of the saturation vapor pressure-temperature curve (kPa· • C −1 ); ρ is the air density (kg·m −3 ); C p is the specific heat at constant pressure (J·kg −1 ·K −1 ); D is the vapor pressure deficit (kPa); and γ is the psychrometric constant. R and R s represent the available energy input above the canopy and the soil surface (W·m −2 ). The specific equations for the S-W model can be found in [13,35]. We followed the methods reported by [16,36] to calculate the soil surface resistance (r ss ), and an empirical equation was found as follows: In Equation (4), θ and θ s are the soil water content (m 3 ·m −3 ) and saturated soil water content (m 3 ·m −3 ), respectively.

Model Evaluation
In this study, the statistical analysis included linear regression, root mean square error (RMSE), and mean absolute error (MSE). RMSE and MSE are calculated as follows: where E i is the value estimated by the S-W model; O i is the observed value; and n is the number of E i and/or O i .

Variation of LAI and Environmental Variables
The leaf area index (LAI) in both years started to increase from May and reached its annual maximum in July, then decreased rapidly because plants began to senesce in September ( Figure 1). The LAI in 2008 was higher than that in 2007, with the maximum values of 1.20 (2008) and 0.96 (2007) m 2 ·m −2 , respectively.
There were no significant differences in annual variation for each environmental variable for the study period (Figure 2), and the corresponding statistical values in the growing season (May-September) and non-growing season are listed in Table 1. The annual maximum value of daily net radiation (R n ) appeared in June with 21.23 and 22.07 MJ·m −2 ·d −1 for 2007 and 2008, respectively, and reached the minimum value in winter. The annual R n in 2007 and 2008 was 2952.11 and 2939.99 MJ·m −2 , respectively, and more than 60% was received during the growing season. The annual variation of soil heat flux (G) followed the same trend of R n with higher and lower values recorded during the growing and non-growing seasons (Figure 2b), respectively, but it fluctuated within a relatively narrow range from −3.06 to 2.27 MJ·m −2 ·d −1 .
The annual variations of air temperature (T a ) and 5 cm soil temperature (T s5cm ) were strongly dependent on R n , however, T s5cm was obviously higher than T a throughout the whole year although T a followed the same trend as T s5cm (Figure 2c). The annual mean of T a was 0.2 • C and −0.6· • C, while the value of T s5cm was 3.9· • C and 3.1· • C for 2007 and 2008, respectively.    Year Growing Phase R n (MJ·m −2 ·d −1 ) Precipitation occurred mainly during the growing season (Figure 2d), which accounted for 89% and 87% of annual precipitation in 2007 and 2008, respectively. The mean daily precipitation was 2.9 mm·d −1 and 2.7 mm·d −1 during the growing season in 2007 and 2008, respectively, while it was only 0.3 mm·d −1 in the non-growing season for both years. The highest monthly precipitation appeared in June for 2007 and 2008, with the value of 147.5 mm and 101.7 mm, respectively. There was a significant variation in SWC 5cm during the growing season ( Figure 2d), which was strongly influenced by precipitation events and ET. This increased rapidly after the occurrence of rainfall and decreased due to E and T when no precipitation occurred.

Annual Variation of ET
Annual change of modeled evapotranspiration (ET SW ) approximately followed the same trend as R n ( Figure 3) and showed a large day to day variation during the growing season. ET SW started to increase from March and reached its highest annual value around July, and then decreased to the lowest value around January. The majority of ET SW occurred during the growing season, which was 409.9 mm and 395.3 mm in 2007 and 2008, accounting for 80% and 79% of annual ET SW , respectively (Table 2).

Evapotranspiration Partitioning
There was an obvious difference between seasonal variations in modeled evaporation (E SW ) and transpiration (T SW ) during the study period ( Figure 4). E SW increased rapidly from early March, reached its maximum value of 3.1 mm (2007) and 2.4 mm (2008) in June, then started to decrease, whereas a relatively lower value was observed in July and August of the growing season. After that, E SW began to increase again with a second peak appearing around October, and then decreased rapidly from late October. T SW in both 2007 and 2008 increased from late April when plants started to grow, with the maximum value in July and/or August when E SW had a relatively lower level. Then, T SW decreased to its minimum value in late October. The annual amount of E SW and T SW accounted for about 60% and 40% of ET SW in 2007, while the values were 56% and 44% in 2008, respectively. During the growing season, the amount of E SW and T SW in 2007 (2008) accounted for 53% (48%) and 47% (52%) of ET SW , respectively, indicating that soil evaporation was higher than plant transpiration, even though the vegetation was fully developed (for more details see Table 2).  To further explore the influence of ESW and TSW on the ETSW, we analyzed the monthly dynamics of ESW/ETSW of the degraded meadow throughout the study period ( Figure 5). During the period from November to April of the next year, ETSW was accounted for by the ESW (i.e., ESW/ETSW = 1.0) because the growth of plants stopped. ESW/ETSW gradually decreased from May with the plant growth and reached its minimum value of 0.35 (in August 2007 and in July 2008), then ESW/ETSW increased again until November. The average monthly ESW/ETSW was 0.50 during the growing season. To further explore the influence of E SW and T SW on the ET SW , we analyzed the monthly dynamics of E SW /ET SW of the degraded meadow throughout the study period ( Figure 5). During the period from November to April of the next year, ET SW was accounted for by the E SW (i.e., E SW /ET SW = 1.0) because the growth of plants stopped. E SW /ET SW gradually decreased from May with the plant growth and reached its minimum value of 0.35 (in August 2007 and in July 2008), then E SW /ET SW increased again until November. The average monthly E SW /ET SW was 0.50 during the growing season.
To further explore the influence of ESW and TSW on the ETSW, we analyzed the monthly dynamics of ESW/ETSW of the degraded meadow throughout the study period ( Figure 5). During the period from November to April of the next year, ETSW was accounted for by the ESW (i.e., ESW/ETSW = 1.0) because the growth of plants stopped. ESW/ETSW gradually decreased from May with the plant growth and reached its minimum value of 0.35 (in August 2007 and in July 2008), then ESW/ETSW increased again until November. The average monthly ESW/ETSW was 0.50 during the growing season.

Diurnal Variation of ET
Diurnal variations of E SW , T SW , and ET SW of clear days in January and July (used to represent winter and summer extreme conditions, respectively) for 2007 and 2008 are shown in Figure 6, where a clear day is defined as that on which the daily transmissivity was greater than 0.7 [18]. All three variables showed the same variation pattern for the two years of 2007 and 2008. In January, ET SW began to increase around 09:00 and peaked between 13:00 and 15:00 with the maximum average value of about 0.02 mm·h −1 , then began to decrease and fell to nearly zero at around 19:00 for both years. In July, however, ET SW began to increase around 07:00 and reached the maximum at about 14:00, then decreased to about zero at 21:00. Although the ET SW showed a similar pattern for July and January, the maximum average value of the former with about 0.6 mm·h −1 , which was much higher than that of the latter. In addition, it was found that daily T SW was higher than E SW in July in both 2007 and 2008 (especially in 2008). In July, daily E SW with the value of 1.8 mm in 2007 was higher than that of 1.5 mm in 2008, and the daily ratio of E SW to T SW (E SW /T SW ) was 0.80 in 2007, which was higher than that of 0.58 in 2008, while the daily ratio of E SW to ET SW (E SW /ET SW ) was 0. 44

Discussion
Our simulation results indicated that E accounts for a major part of ET d tation degradation. For environmental (biotic and abiotic) factors that may a its partitioning, and for the model validation, our research showed the follow

Discussion
Our simulation results indicated that E accounts for a major part of ET due to vegetation degradation. For environmental (biotic and abiotic) factors that may affect ET and its partitioning, and for the model validation, our research showed the following results: 1.
E/ET in our research site was more sensitive to change in LAI. E/ET decreased rapidly with the increase of LAI (paragraph 1 in Section 4.1); 2.
Grassland ecosystems with lower LAI and/or vegetation coverage may lose more water through ET (paragraph 2 in Section 4.1).

3.
Net. radiation had little effect on ET partitioning, but had a great influence on ET, E, and T (paragraph 2 in Section 4.2).

4.
Air temperature had a greater effect on T than on E (paragraph 3 in Section 4.2).

5.
Soil water content at a 5 cm depth affected both ET and ET partitioning in this degraded meadow, especially for the E (paragraph 4 in Section 4.2). 6.
Vapor pressure deficit had little effect on both ET and ET partitioning (paragraph 5 in Section 4.2). 7.
Leaf area index is an important factor influencing ET partitioning (paragraph 6 in Section 4.2). 8.
The model results had good agreement with the ET observed by the eddy covariance system (paragraph 2 in Section 4.3).

Effects of Vegetation on Evapotranspiration Partitioning
ET is mainly dependent on vegetation, meteorological conditions, and soil water [37], and the partitioning of ET into E and T is strongly influenced by changes in vegetation characteristics during the growing season [38,39]. Leaf area index (LAI) is often used to quantify terrestrial ecosystem ET as well as ET partitioning; an increase in leaf area will initially increase ET when the soil water content is high, and this response will weaken at high LAI [40,41]. In this degraded meadow, the variation pattern of E is quite different from that of T during the growing season ( Figure 4). Seasonal variation of T followed the same trend of LAI, with the higher values recorded around August at higher LAI (Figure 1), however, the highest E occurred around June. The result is consistent with many literature reports [39,42,43]. The analysis of the relationship between LAI and E/ET during the growing season is illustrated in Figure 7. The E/ET data were LAI-bin averaged because this data compilation helped to reduce or offset the errors associated with the measurements [34]. The LAI gaps were linearly interpolated to daily intervals [44]. It was found that there was a significant negative correlation between E/ET and LAI for 2007 and 2008 (Figure 7) (i.e., the contribution of soil evaporation to evapotranspiration decreased linearly with the increase in LAI), which is consistent with other alpine meadow ecosystems reported by [17]. A previous study also reported a similar negative relationship in multiple ecosystems (e.g., forests, crops, wetlands, shrubs, and grasses) [45]. Several studies have reported that E/ET initially decreased rapidly with an increase in LAI at the low vegetative cover (low LAI), while the response of E/ET to LAI decreased gradually with the increasing LAI, and finally approached a constant value [14,17,43,45]. In the present study, however, the LAI was very low, even in the peak growing season, with the maximum value of 1.20 m 2 ·m −2 (July 2008) due to the meadow degradation, therefore, E/ET was more sensitive to change in LAI. Based on the diurnal variation of ET partitioning in the peak growing season of July (Figure 6), the daily E/ET in 2007 was obviously higher than that in 2008 due to the relatively high LAI in 2008 compared with 2007. In addition, E in the growing season of 2007 was higher than that of 2008 whereas the opposite result was obtained for T (Table 2), and the regression line between E/ET and LAI for 2007 was above that for 2008, indicating that the contribution of E to ET in 2007 was higher than that in 2008, which may be due to the lower LAI for 2007 compared with 2008 that resulted in the increase in ET of this alpine meadow (Table 2).
To further investigate the relationship between vegetation and ET partitioning, we made a comparison between our results and some of the previously published studies on grassland ecosystems ( Table 3). All of these studies reported a negative relationship between vegetation conditions and E/ET, which is consistent with our research. However, it is worth noting that worse vegetation conditions corresponded to higher ET/P (the ratio of evapotranspiration to precipitation) ( Table 3). We suggest that this is due to the degradation of vegetation thus allowing more energy to reach the soil surface, which leads to increased E and ET. Furthermore, Gu et al. [9] found a curve relationship between the aboveground biomass and ET at an alpine meadow ecosystem. That is, ET increased gradually with the increase in aboveground biomass at first, but decreased thereafter despite the biomass still increasing. To further investigate the relationship between vegetation and ET partitioning, we made a comparison between our results and some of the previously published studies on grassland ecosystems ( Table 3). All of these studies reported a negative relationship between vegetation conditions and E/ET, which is consistent with our research. However, it is worth noting that worse vegetation conditions corresponded to higher ET/P (the ratio of evapotranspiration to precipitation) ( Table 3). We suggest that this is due to the degradation of vegetation thus allowing more energy to reach the soil surface, which leads to increased E and ET. Furthermore, Gu et al. [9] found a curve relationship between the aboveground biomass and ET at an alpine meadow ecosystem. That is, ET increased gradually with the increase in aboveground biomass at first, but decreased thereafter despite the biomass still increasing.

Effects of Environmental Factors on Evapotranspiration Partitioning
Except for the vegetation LAI, solar radiation, temperature, soil moisture, and air humidity will also have an impact on the partitioning of evapotranspiration [9,10,50,51]. Consequently, to comprehensively understand the control of environmental factors on the balance between evaporation and transpiration, net radiation (R n ), air temperature (T a ), 5 cm soil water content (SWC 5cm ), vapor pressure deficit (D) as well as the leaf area index (LAI) were chosen to analyze the effects of the above factors on the ET and its partitioning (E/ET and T/ET). In this study, we referred to the method by [42], where each dependent variable was multiplied by 0.5 and 2.0 in the model, respectively, then the model was rerun to see how much the output value changed (Figure 8 and Table 4). Here, multiplying by 0.5 and 2.0 is defined as "low level" and "high level", respectively, and "standard" is the observed value. Solar radiation is the most important source of energy for most biological and meteorological processes, and ET, E, and T are dependent mainly on the solar energy available to vaporize the water [50]. It was found that R n had little effect on ET partitioning (E/ET and T/ET) (Figure 8a), but had a great influence on ET, E, and T (Table 4), which is consistent with previous research results [50]. ET, E, and T were decreased about 66% compared with the standard when R n was at a low level and increased about 132% when R n was doubled ( Table 4). The results suggest that the change of R n strongly influences ET, E, and T, but there was almost no influence on E/ET and T/ET in this degraded meadow ecosystem. In this study E/ET, the ratio of soil evaporation to evapotranspiration; T/ET, the ratio of plant transpiration to evapotranspiration; ET/P, the ratio of evapotranspiration to precipitation.  Table 4. Effects of net radiation (R n ), air temperature (T a ), 5 cm soil water content (SWC 5cm ), vapor pressure deficit (D) and leaf area index (LAI) on ET, E, and T in the growing season of 2007 and 2008 for the degraded meadow.

Input Variables
Percentage of Variation −50% +100% Temperature is one of the major factors affecting the rate of ET, E, and T, and temperature-based models are widely used to estimate ET [51,52]. Our results showed that E/ET was slightly higher (or lower) than the standard when T a was at a low level (or high level), while an opposite change was found for T/ET (Figure 8b). However, T a had a positive relationship with ET, E, and T (Table 4), and the effect of T a on T was greater than ET and E (Table 4). In this study, the variation of T almost followed the same trend of LAI and T a with a higher value in about July (Figures 1, 2c and 4), indicating that T increased with the increase in LAI and T a , but E decreased rapidly with the increasing LAI, therefore the response of T to T a was more sensitive compared with the E.
Numerous studies have shown that ecosystem ET is closely related to the soil water content [10,25,47]. Our results indicated that SWC 5cm affected both ET and its partitioning in this degraded meadow, especially for the E (Figure 8d, Table 4). It was observed that E/ET and T/ET showed the opposite change trend when the SWC 5cm was multiplied by 0.5 and 2 (Figure 8d), respectively, in which increasing SWC 5cm significantly increased E/ET and decreased T/ET, and the converse was also true (Figure 8d). Soil water content is an important factor controlling soil surface resistance, and increasing SWC 5cm can reduce bare soil surface resistance to evaporation, and at the same time, increase the supply of soil moisture, resulting in an increase of E and E/ET. The previous study pointed out that transpiration will increase rapidly with the increase in soil water content when water supply is limited [53], which is inconsistent with our results. However, transpiration is strongly dependent not only on the soil water content, but also on meteorological and vegetation conditions. In this alpine meadow, most of the root system was distributed within the 0 to 10 cm surface layer, and the soil maintained a relatively high-water content throughout the growing season due to the abundant precipitation (Figure 2d), while a downward trend of SWC 5cm was observed in the peak growing season of July-August due to the high ET. Therefore, under the condition that other observed variables are included in the model, our results showed that T and T/ET decreased when only increasing SWC 5cm , and the possible reason is that T is predominantly controlled by R n and LAI, and at the same time, affected by the SWC 5cm . The model result is consistent with the actual change in transpiration, and a similar relationship was also observed between the T and SWC 5cm in another study by using micro-lysimeter experiments in an alpine meadow of the TRSR (article in print).
Vapor pressure deficit (D) is one of the principal weather variables affecting ET because D affects the evaporation demand of the atmosphere and canopy conductance [9]. Our results showed that D had a very small effect on both ET and its partitioning (Figure 8c, Table 4), perhaps because the value of measured D was very low and varied within a very narrow range from 0.03 to 1.64 kPa in this degraded meadow (Figure 2e), which was significantly lower than many other grassland values with the maximum D ranging from about 2 to 5 kPa [4]. Thus, when D was multiplied by 0.5 and 2.0, there was almost no change in ET as well as its partitioning.
LAI is one of the important parameters describing vegetation characteristics, which is widely adopted in ET partitioning [6,39]. The simulation results showed that E/ET will increase or decrease significantly compared with the "standard" level when LAI is multiplied by 0.5 or 2.0, while in contrast, an opposite trend was observed for T/ET with increasing or decreasing LAI (Figure 8e) and it was also noted that there was the same trend for E and T (Table 4). Our results are consistent with earlier studies that showed the effects of LAI on ET partitioning [41,45]. However, the effect of LAI on ET was relatively small because the E pattern was almost the opposite to that of T during the growing season (Figure 4), so increasing E may be offset by the decreasing T. Usually, increasing LAI can increase vegetation cover and lead to a decrease in bare soil surface area, then decreases E and/or increases T, and the reverse is also true.
Overall, LAI and SWC 5cm are the main important factors influencing E/ET and/or T/ET. E and T were primarily controlled by R n , LAI, and SWC 5cm , while the effect of T a on T was relatively large compared with E. D had little effect on both ET and ET partitioning.

Validation of the Shuttleworth-Wallace Model
The eddy covariance (EC) system was conducted in our flat degraded meadow. The WPL density correction was applied to water vapor flux [54], and the energy balance ratio (EBR) was calculated using the following equation [55]: where H, LE, and G are the sensible, latent, and soil heat fluxes. In this study, the term (LE + H), measured by the EC method, seemed to be underestimated since the average value of EBR was 0.79 in the study period, which fell in the median region of reported energy closures, which ranged from 0.55 to 0.99 [55]. LE was converted to ET (mm) by assuming a value for a conversion factor of 2450 J/g. In order to verify the performance of the S-W model over the alpine meadow, we compared the ET estimated by the S-W model (ET SW ) to that measured by the EC method (ET eddy ) ( Figure 9). Overall, there was a good agreement between ET SW and ET eddy in the study period, while the ET SW was underestimated compared to the ET eddy from December to April of the next year ( Figure 9). Gong et al. [56] also pointed out that the S-W model overestimated and/or underestimated ET at different growth stages in comparison with the results measured by the lysimeter. Therefore, we performed some statistical analyses between ET SW and ET eddy in the growing and non-growing seasons (Table 5). It was found that the model performance in the growing season was better than that in the non-growing season, and the model overestimation of ET occurred in the growing season, while the underestimation appeared in the non-growing season. The relationships between ET SW and ET eddy at different growth stages for 2006-2008 were summarized through statistical analyses ( Figure 10 and Table 5). The linear regression slopes (k) ranged from 1.04 to 1.06 and the values of R 2 (square of the correlation coefficient) were over 0.91. RMSE and MAE varied in the range of 0.3-0.6 and 0.2-0.5 mm·d −1 , respectively. The possible reason for this overestimation of ET might be due to the lack of energy balance closure of the eddy covariance method (EBR = 0.79), and thus the ET was underestimated by the EC method. Chen et al. [37] reported that the S-W model overestimated ET by comparing it with the measured data. Wei et al. [57] indicated that the S-W model overestimated ET by 5% when compared to the measured data, which is similar to our simulation results. On the whole, the ET was well estimated by the S-W model in this degraded meadow on the Qinghai-Tibetan Plateau.   Table 5. Statistical analyses of evapotranspiration estimated by the S-W model (ETSW) and measured by the eddy covariance (ETeddy) in the different periods.

Conclusions
We estimated evapotranspiration (ET) and its partitioning with the S-W model in a degraded alpine meadow in the TRSR, and compared the results with data obtained from eddy covariance. The validation confirmed the good performance of the S-W model for the prediction of ET and its partitioning in this study. Net radiation is the most important factor influencing ET while leaf area index (LAI) is a key factor affecting ET partitioning. Due to the vegetation degradation at our research site, the contribution of soil evaporation (E) accounted for the main part of ET, and ET was higher with the lower LAI. Our results suggest that the water lost by ET from the meadow ecosystem increased with the increasing intensity of vegetation degradation, that is, meadow degradation would increase water loss through increasing E, leading to reducing the water conservation capacity of the alpine meadow ecosystem in TRSR.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy and ethical concerns. All the data used in this study were observed by the authors of the Nankai University and Northwest Institute of Plateau Biology, Chinese Academy of Sciences. The observation equipment was purchased with funds from research projects. Data for this research are not publicly available due to a data confidentiality agreement. Researchers who are related to these projects can obtain these data, and we are not able to provide the original dataset to others at the moment.