Determining the Contributions of Vegetation and Climate Change to Ecosystem WUE Variation over the Last Two Decades on the Loess Plateau, China

: Exploring the variations in the water use efficiency (WUE) is helpful in gaining an in ‐ depth understanding of the regional carbon and water cycles on the Chinese Loess Plateau (CLP). Here, we employed the spatial variations in the WUE and the quantitative contributions of the in ‐ fluencing factors, including the precipitation (P), temperature (Temp), vapor pressure deficit (VPD), sunshine duration (SD), and leaf area index (LAI), with the drought index varying over the last two decades. Results showed that the multiyear average WUE decreased significantly as the drought index increased for all of the vegetation types. Per ‐ pixel interannual variability of WUE trend was 0.024 gC ∙ m − 2 ∙ mm − 1 ∙ yr − 1 . As the drought index increased, the WUE initially increased and then de ‐ creased for the forests, grassland, and shrubland, and their peaks occurred at drought index values of 2.60–3.10. Among the influencing factors, the WUE was predominantly controlled by the LAI, with an impact and relative contribution of 0.014 gC ∙ m − 2 ∙ mm − 1 ∙ yr − 1 and 58.3%, respectively. The P and SD contributed the least to the trend in WUE, and impact and relative contribution of both were 0.001 gC ∙ m − 2 ∙ mm − 1 ∙ yr − 1 and 4.17%. Our study also demonstrated that the LAI was the dominant factor affecting the WUE trends for grassland and the Yan River due to the structural parameters and geographical location. In addition, the impact and relative contribution of the residual factors on the WUE trend were 0.004 gC ∙ m − 2 ∙ mm − 1 ∙ yr − 1 and 16.7%. Our findings suggested that comprehen ‐ sive effects such as micro ‐ geomorphic changes and nitrogen deposition could not be ignored except for vegetation and climate change. This study will clarify the spatial and temporal evolution of WUE and its influence mechanism.


Introduction
The ecosystem water use efficiency (WUE), which is defined as the ratio of the gross primary productivity (GPP) to the evapotranspiration (ET) [1,2], is a critical parameter that reflects the carbon assimilation through water loss [3,4]. The WUE has generally been recognized as an important integral indicator linking the carbon and water cycles of ecosystems and the atmosphere [5], which can quantify how much carbon an ecosystem gains relative to the water used per unit [3]. Therefore, it is necessary to investigate the WUE to deeply understand and simulate the ecosystem carbon-water coupling relationship [4]. In addition, it can provide useful information about the dynamic responses of ecosystems to environmental changes [6].
In the past few decades, field measurements, gas exchange, and stable isotope techniques have gradually been used to study the WUE at plant leaf and individual scales [7][8][9][10]. Field measurements and gas exchange can accurately obtain the carbon and water balance at the monitoring stations to better and more accurately study the water use efficiency at small scales [7][8][9]. Carbon isotopic composition is a useful method to determine the drought resistance, to provide insights of chemical, physical and metabolic processes involved in carbon transformation in the stressed plant that determine the water utilization and water conservation status for the forest species [10]. With the rapid development of the eddy covariance technique [11] and remote sensing technology [12,13], studies of the WUE have been extended to the ecosystem scale. Remote sensing technology has become an effective technical approach for studying the WUE of regional and global ecosystems because it has the powerful capability to make temporally and spatially continuous observations of the Earth's surface over large areas and long time periods [14]. Currently, numerous studies have been conducted on the WUE at the regional scale using remote sensing technology products, such as GPP and ET products. Sharma and Goyal et al. revealed the characteristics of the spatial and temporal variation in the WUE at the district scale using moderate-resolution imaging spectroradiometer (MODIS) net primary productivity and ET data [15,16]. Tang et al. analyzed the characteristics of the global WUE and concluded that changes in land use were the main cause of the significant decline in the WUE using National Aeronautics and Space Administration (NASA) TERRA and AQUA MODIS-based estimates of the GPP and ET [17]. Guo et al. analyzed WUE in the Beijing-Tianjin Sand Source Region and found that the annual growth rate of WUE in this area was very high and human activities were the dominant factor for the increase of WUE based on integrating remote-sensing and the eddy covariance technique observed ET data and GPP data by using a machine learning approach [18].
The Chinese Loess Plateau (CLP), which is located in the middle reaches of the Yellow River, is the second largest plateau in China. The CLP is considered to be one of the most typical ecologically fragile areas associated with harsh environmental conditions and intensified human activities in the world [19]. In this region, the precipitation is far less than that demanded for vegetation growth, and water resource shortages are an essential limiting factor of ecological and economic development [20]. Since 1999, the Chinese government has implemented the Grain for Green (GFG) Project to alleviate serious ecological problems, which has significantly increased the vegetation cover and reversed the soil erosion trend [21]. However, it has been found that depletion of deep soil water by vegetation growth through evapotranspiration has led to an increased soil moisture deficit in recent years [22,23]. In addition, a significant streamflow reduction in the CLP has also been observed [24,25]. Therefore, studying the carbon-water coupling and its driving factors is significant for assessing the water consumption and sustainability of vegetation on the CLP [19,26].
Currently, numerous studies about the long-term WUE variation have been conducted in the whole CLP [27][28][29][30][31][32]. Cao et al. have shown that the WUE in the CLP has experienced a significant increasing trend since 1999 [33]. The results of Tian et al. have found the WUE decreased with gradients from the southeast, but with sharp rises detected in northern Shaanxi, China [32]. Zheng et al. analyzed the responses of WUE variation to regional climate and vegetation coverage changes in the CLP-based correlation coefficients and indicated that vegetation coverage change is the primary cause for the increasing ecosystem annual WUE in the CLP [31]. Zheng et al. examined spatial responses of plant canopy WUE and ecosystem WUE variations to the mean annual precipitation, mean annual air temperature, and normalized difference vegetation index gradients; results showed that normalized difference vegetation index gradients were the important factors [30]. In summary, there is a significant spatial divergence in WUE and factors such as precipitation, temperature, and vegetation coverage changes all have an impact on WUE. However, the quantitative impacts of the different influencing factors, es-pecially the comprehensive effects with other factors (e.g., microtopographic changes, increasing carbon dioxide (CO2), and nitrogen deposition) affecting the WUE have not been discussed in detail.
Therefore, the objectives of this study were (1) to analyze the spatial variation patterns of the multiyear mean WUE and the variations in the WUE with drought index changes for different vegetation types in the CLP from 2000 to 2018; (2) to explore the spatial variation patterns of the partial correlations between the WUE and climatic and vegetation factors with the drought index change for different vegetation types in the CLP from 2000 to 2018; and (3) to quantitatively calculate the impacts of the climate, vegetation, and other factors on the WUE trends, and their spatial variations for different vegetation types and sub-basins in the CLP from 2000 to 2018.

Study Area
The CLP is located in northwestern China, and it covers an area of approximately 6.4 × 10 5 km 2 (100°54′-113°33′ E，33°43′-41°16′ N; Figure 1a). The CLP is a typical arid and semiarid area where water resources have become a bottleneck in regional social economic development and ecological growth [22]. The multiyear average annual precipitation and temperature are 300-600 mm and 4-14 °C on the CLP, respectively, (Figure 1b,c). Historically, climate change and intense human activities have led to severe degradation of the ecosystem of the CLP [19]. Since 1999, major ecological restoration projects have been implemented on the CLP. The ecological environment and the capacity of the soil and water conservation have been significantly improved as a result of dramatic changes in the underlying surface of the CLP over the past two decades [21,34].

Data and Preprocessing
In this study, we selected the PML-V2 product to obtain the ET and GPP data. The terrestrial water storage anomaly (TWSA) was derived from the average of the Gravity Recovery and Climate Experiment (GRACE) Jet Propulsion Laboratory (JPL) RL06_mascon and Center for Space Research (CSR) RL06_mascon products. The Global Land Surface Satellite (GLASS) leaf area index (LAI) product was selected to characterize the vegetation structure parameter. We obtained land cover data from the GLC_FCS30, which was the first global land cover product with fine classification [35](http://data.casearth.cn, accessed on 20 December 2020). Detailed information about these remote sensing products is provided in Table 1.
The field measurements of the precipitation (P), temperature (Temp), wind speed (WS), sunshine duration (SD), and relative humidity (RH) were measured at 295 meteorological stations and were obtained from the China Meteorological Data Service Center (http://data.cma.cn/, accessed on 18 June 2020). According to the shared ET and GPP flux database published by the China Flux Observation and Research Network (ChinaFLUX) (http://www.chinaflux.org/, accessed on 1 March 2020), the observations from four flux stations (Du2, Du3, HaM, and Sw2) located near and within the CLP were collected. The annual runoff measured at five major gauging stations (Tangnaihai, Qingtongxia, Toudaoguai, Longmen, and Huayuankou) located along the main river channel were obtained from the Hydrological Bureau of the Yellow River Conservancy Commission. Detailed information about these field measurements is provided in Table 2.
To match the spatial resolution of the GPP and ET products, we resampled the LAI and land cover products at a resolution of 500 m. All of the data from meteorological stations were spatially interpolated to a resolution of 500 m using the AUSPLINE software. We calculated the ET0 and VPD based on the meteorological data using the Food and Agriculture Organization (FAO) Penman Monteith method [36] and the equations proposed by Yuan et al. [37], respectively. We calculated the drought index using the ratio of the annual average potential evaporation (ET0) to the annual average precipitation (P) during 2000-2018 [38], and then divided it into 10 regions based on the histogram distribution of the drought index ( Figure 1e). All of the time-series data were temporally aggregated at the annual scale for subsequent analysis.

Validation of GPP and ET Product Accuracy
The accuracy assessment of the PML_V2 product globally and in China showed that GPP and ET were more accurate than the other currently available global products [39,43,44]. Nevertheless, to ensure the reliability of the subsequent analysis, we assessed the GPP using the measurements from the flux stations at the annual and 8-day scales. We evaluated the ET using the flux observations at the 8-day scale, assessed the regional ET at the annual scale using the regional ET calculated using the water balance principle based on the measured P, runoff, and derived changes in the terrestrial water storage from the TWSA for the four sub-regions, including Tangnaihai-Qingtongxia (Tang-Qing), Qingtongxia-Toudaoguai (Qing-Tou), Toudaoguai-Longmen (Tou-Long), and Longmen-Huanyuankou (Long-Hua). The accuracy metrics used were the determination coefficient (R 2 ), the root mean square error (RMSE), and the Nash-Sutcliffe efficiency coefficient (NSE).

Trend Slope Analysis
We assessed the trends of the WUE and influencing factors by calculating the linear slope values on an annual timescale during the study period per pixel (Equation (1)). We used the F test to examine the statistical significance of the WUE trends and considered the interannual trends to be statistically significant if the p-values were < 0.05. Additionally, the time series with missing WUE values were excluded from the temporal trend analysis in this study. The linear slope equation is as follows: where i is the individual year, varying from 1 to19; Xi is the WUE and influencing factors in the ith year; and n is 19, representing the total number of years. When slope > 0, X is increasing; when slope < 0, X is decreasing.

Partial Correlation Analysis
Numerous previous studies have shown that P, Temp, VPD, SD and LAI are the key factors influencing the photosynthesis and transpiration in an ecosystem. Therefore, we used the partial correlation coefficient method to analyze the correlations between the interannual variation in the WUE and the influencing factors. The partial correlation coefficient formula is where , …… is the nth (n ≤ k − 2) order partial correlation coefficient between xi and xj, and k is the number of the variables. Each term on the right side of Equation (2) represents the (n − 1)th order partial correlation coefficient.

Quantitative Attribution Method of the WUE Trend
We used multiple linear regression to quantitatively calculate the sensitivity of the WUE to the different influencing factors in the CLP, as follows: where , , , and are the sensitivities of WUE to P, Temp, VPD, SD, and LAI, respectively. WUE, P, Temp, VPD, SD, and LAI are the annual time-series datasets during 2000-2018.
The WUE trend was assumed to be the sum of the WUE trend driven by the climatic factors (P, Temp, VPD and SD), the LAI, and other factors (e.g., microtopographic changes, increasing CO2, and nitrogen deposition). The differential form is as follows: (4) where is the net slope of the WUE trend, and , , , , , and are the impacts of P, Temp, VPD, SD, LAI, and others on the WUE trend, respectively.
According to the basic concept of total differentiation, the impacts of the climatic factors and LAI in Equation (4) can be decomposed via partial differentiation [14,45]: (5) where , , , , and are the slopes of P, Temp, VPD, SD, and LAI, respectively; , , , , and are the sensitivity coefficients derived from Equation (3); and is calculated using the residual method.
The relative impact rate of each influencing factor on the WUE trend was calculated as follows: where .
is the relative impact rate of the influencing factor (P, Temp, VPD, SD, LAI, or other factors) on the WUE trends. Finally, the influencing factor that had the highest relative impact rate was defined as the factor that dominated the interannual WUE trend during 2000-2018 [46].  (Figure 2d). The validation results based on both the sites and the regional scale indicated that the PML_V2 remote sensing products generally performed well and were suitable for the analysis of the spatial-temporal pattern of the WUE and the impact mechanisms on the CLP.

The Spatial-Temporal Variations in the WUE during 2000-2018
Spatially, the annual average WUE decreased gradually from the southeast to the northwest (Figure 3a), with an annual average of 1.242 gC•m −2 •mm −1 during 2000-2018 on the CLP. Overall, the spatial pattern of the WUE was similar to that of the drought index ( Figure 1e). The annual mean WUE was higher in the southeast, where the drought index was lower, and was lower in the northwest, where the drought index was higher. This result suggested that the drought index was a vital influencing factor affecting the zonal distribution of the WUE. As shown in Figure 3b, the WUE of all of the vegetation types decreased significantly with increasing drought index. Among them, the WUE of the grassland decreased the fastest with increasing drought index, which indicated that it was the most sensitive to the change in the drought index. The WUE of the shrubland decreased the slowest with increasing drought index, which indicated that it was the least sensitive to changes in the drought index. We also found that the annual mean WUE of the forests was the highest, that of the cropland was the second highest, that of the grassland was the third highest, and that of the shrubland was the lowest. Figure 4a suggests that the WUE trend increased by about 87.32% on the CLP, and 47.36% exhibited a significant decrease during 2000-2018. In contrast, 16.11% of the study area exhibited a decreasing WUE trend, with a significant decrease in 2.92% of the areas in the western part of the CLP. The temporal variations in the annual mean values of the GPP, ET and WUE on the CLP significantly increased for both GPP and ET, and we found an insignificant positive WUE trend with a value of 0.024 gC•m −2 •mm −1 •yr −1 (Figure 4b).
The temporal variations in the WUE for different vegetation types on the CLP are shown in Figure 5a. It is evident that the interannual WUE variation characteristics of the different vegetation types from 2000 to 2018 were highly consistent, that is, they all exhibited significantly increasing trends (p < 0.01). In addition, the WUE trend of the cropland was 0.027 gC•m −2 •mm −1 •yr −1 and was the fastest increasing trend. The WUE trend of the shrubland was only 0.013 gC•m −2 •mm −1 •yr −1 and was the slowest increasing trend. Moreover, we calculated the WUE trends of the different vegetation types with the drought index (Figure 5b). We found that differences remained in the WUE trends with the drought index among the vegetation types, although the WUE trends for all of the vegetation types followed significantly increasing trends.

Spatial Pattern of Correlation between WUE and Influencing Factors
We performed partial correlation analyses to explore the relationships between the WUE and the influencing factors (P, Temp, VPD, SD and LAI) in each pixel, and the results revealed differences in the spatial patterns of the correlations between the different influencing factors and the WUE (Figure 6). The WUE was negatively correlated with P for 59.42% of the entire CLP, mainly in the northwestern region (Figure 6a). The area of the CLP where WUE was positively correlated with Temp was 63.56% (Figure 6b). A positive correlation between the WUE and VPD was observed in the southeastern part of the study area (65.8%), whereas a negative correlation was observed in the northwestern part of the study area (30.78%) (Figure 6c). Regarding the partial correlation between the WUE and SD, we found that the regions with a positive partial correlation were distributed mainly in the northern part of the CLP (Figure 6d). About 89.56% of the entire area exhibited a significant positive correlation between the WUE and LAI, and the high correlation coefficients values (R > 0.8) were distributed mainly in the southeastern part of the study area (Figure 6e). Figure 6f shows the spatial patterns of the most strongly correlated factors for each pixel. The results indicated that the regions where the most correlated factor accounted for the largest area of the entire CLP (more than 62.18%) were generally distributed in the central part of the study area. The areas where the Temp and VPD were the most correlated factors were 10.93% and 8.80%, respectively. The SD was the most correlated factor in 7.36% of the area. P was the most correlated factor over 7.31% of the area (smallest area). To identify the spatial divergence pattern of the WUE, we further analyzed the partial correlation coefficients of the WUE and the influencing factors (P, Temp, VPD, SD and LAI) for different vegetation types in relation to the drought index (Figure 7). We found that the partial correlation coefficients between the WUE and LAI were the highest for all of the vegetation types. The correlation between the WUE and the influencing factors for the cropland had the lowest sensitivity to the drought index (Figure 7a). The threshold effect existed in the relationship between the WUE and the LAI for the forests and grassland, which means that the correlation coefficients of both were characterized by increasing and then decreasing trends with increasing drought index, and the peak values occurred at drought index values of 2.79-3.10 (Figure 7b,c). Across the forest ecosystems, the correlations between the WUE and the meteorological factors (P, Temp, VPD and SD) were generally positive for drought indexes of 0 to 4.74, but they changed abruptly to negative correlations when the drought index was > 4.74 (Figure 7b). This result suggested that the WUE of the forests was negatively correlated with the meteorological factors under extreme drought conditions. We observed similar spatial variation characteristics of the correlations between the WUE and the influencing factors for the grassland and shrubland (Figure 7c,d). As the drought index increased, the correlation between the WUE and P changed from negative to positive and then to negative again; and that between the WUE and Temp changed from negative to positive for the grassland and shrubland. The links between the WUE and VPD for the grassland and shrubland were characterized by a trend that increased and then decreased with increasing drought indices.

Quantitative Impacts of Influencing Factors on the WUE Trend
Based on Equations (2)-(4), we quantitatively calculated the impacts and relative impact rates of the different influencing factors on the WUE trends during 2000-2018. Figure  8 shows the spatial patterns of the impacts and the relative impact rates of the different influencing factors on the WUE trend over the entire CLP. The annual growth of the WUE on the CLP during 2000-2018 was 0.024 gC•m −2 •mm −1 •yr −1 . The LAI had the largest effect on the WUE trend, and the areas displaying a significant positive contribution between the LAI and WUE accounted for 86.22% of the total area of the CLP. We observed a negative contribution in 10.12% of the total area, mainly in the southern part of the study area (Figure 8(a5)). In the entire CLP, the impact and relative impact rate of the LAI on the WUE trends were 0.014 gC•m −2 •mm −1 •yr −1 and 58.33%, respectively (Figure 8(b5)). The impacts and relative impact rates of the meteorological factors on the WUE trends are shown in Figure 8(a1-a4,b1-b4). It is evident that the negative impact of P was distributed mainly in the northwestern arid region of the CLP, which indicated that the variation in P led to the continuous decrease in the WUE trend in this region (Figure 8(a1)). In addition, the impact and relative impact rate of the P on the WUE trend were 0.001 gC•m −2 •mm −1 •yr −1 and 4.17%, respectively (Figure 8(b1)). The Temp contributed positively to the WUE trend in 63.91% of the CLP, and the contribution was relatively large in the southeastern part of the study area (Figure 8(a2)). The increase in the WUE trend caused by the variation in Temp was 0.002 gC•m −2 •mm −1 •yr −1 and its relative impact rate was 8.33% (Figure 8(b2)). The regions in which the WUE increased and decreased because of variations in the VPD accounted for 60.73% and 29.81% of the total area of the study area, respectively ( Figure  8(a3)). In addition, the impact and relative impact rate of the VPD on the WUE trends were 0.002 gC•m −2 •mm −1 •yr −1 and 8.33%, respectively (Figure 8(b3)). The SD made a positive contribution in the northern part of the CLP (Figure 8(a4)), and its impact and relative impact rate were 0.001 gC•m −2 •mm −1 •yr −1 and 4.17%, respectively (Figure 8(b4)). Overall, the integrated impact and relative impact rate of the meteorological factors on WUE trend were 0.006 gC•m −2 •mm −1 •yr −1 and 25.00%, respectively. As also can be observed from Figure 8a6, the area in which the WUE increased (82.24%) because of the other factors was much larger than that of the factors leading to a decrease in the WUE (14.22%), and its impact and relative influence rate were 0.004 gC•m −2 •mm −1 •yr −1 and 16.67%, respectively (Figure 8(a6,b6)). We calculated the impacts and relative impact rates of the different influencing factors on the WUE trend for different vegetation types (Figure 9). The result showed that except for the negative impact of P on the WUE in the shrubland, mainly in the arid zone, the other factors influencing the WUE trends for the different vegetation had positive impacts. By comparing and analyzing the relative impact rate of each influencing factor on the WUE for the different vegetation types, we found that the LAI was the dominant factor influencing the change in WUE for all of the vegetation types, and the magnitudes of its effect on the WUE trends of the vegetation types were in the following order: grassland > cropland > forest > shrubland. By comparing the relative impact rates of the P, Temp, SD and VPD of each vegetation type, we also concluded that the Temp had a positive effect and it was the dominant meteorological factor affecting the WUE trend (Figure 9b). Among the different vegetation types, the Temp contributed the most to the forest WUE, and its impact and relative influence rate were 0.004 gC•m −2 •mm −1 •yr −1 and 16.46%, respectively. Compared with the integrated impact of the meteorological factors, the combined effect of the other factors was smaller. We further estimated and analyzed the effects of the influencing factors on the WUE trends in the different sub-basins in the CLP, including the Yan River, Wuding River, Jing River, Kuye River, Beiluo River, Huangpuchuan River, Wei River, and Fen River sub-basins. We found that the factors played mostly positive roles ( Figure 10). Thus, we concluded that the LAI was the dominant factor influencing the WUE trends in each subbasin, and its relative impact rates in the different sub-basins were ranked as follows: Yan River > Wuding River > Kuye River > Huangpuchuan River > Jing River > Beiluo River > Fen River > Wei River. Among the four meteorological factors, the Temp had the greatest effect on the WUE in the Fen River, Wei River, Jing River, Beiluo River, Yan River, and Huangpuchuan River sub-basins, and its relative impact rates were 13.04%, 17.35%, 12.20%, 14.23% and 6.55%, respectively (Figure 10b). The SD was the most influential meteorological factor in the Wuding River and Kuye River sub-basins, and its relative impact rates were 5.59% and 9.07%, respectively (Figure 10b). Overall, the meteorological factors had the large integrated effect on the WUE trends in the Wei River, Beiluo River, and Fen River sub-basins; however, they had a small integrated effect on the WUE trends in the Kuye River, Yan River, and Wuding River sub-basins. In contrast, we found the other factors to have a large impact on the WUE trends in the Wei River, Beiluo River, and Jing River sub-basins (Figure 10b).

Analysis of Spatial Heterogeneity of Multiyear Average WUE and WUE Trend
From 2000 to 2018, there was a significant increase in the WUE (0.024 gC•m −2 •mm −1 •yr −1 ) across the CLP. This increase was mainly because the growth rate of the GPP was significantly larger than that of the ET (Figure 4b). Spatially, the WUE in the study area was characterized by a large variability (i.e., gradually decreasing from southeast to northwest), which also was reported by Han et al. [31]. This spatial distribution was basically consistent with that of the drought index, and it was caused by the sufficient precipitation supply and large LAI in the humid areas and strong vegetation photosynthesis, which resulted in a higher GPP. Simultaneously, the low precipitation, sparse vegetation, and high soil evapotranspiration in the arid areas made it difficult for water to infiltrate and be used for vegetation growth, resulting in a low WUE [47].
By investigating the spatial variation characteristics of the WUE for the different vegetation types, we found that the WUE trends of the vegetation ecosystems decreased with increasing drought index. This may have been caused by the relative abundant precipitation and low evapotranspiration in the areas with low drought indexes. Moreover, the vegetation grew well and the vigorous photosynthesis was able to convert more CO2 into organic matter in the atmosphere [48]. This, coupled with the higher vegetation cover, which shaded the soil, reduced the soil evaporation to some extent [49]. Taken together, most of the water consumption in this area was used for organic matter accumulation by the vegetation itself, resulting in a high WUE. That is, as the drought index increased, the precipitation and vegetation cover decreased and the soil evaporation increased. These reasons jointly contributed to the decrease in the total carbon sequestration and the gradual decrease in the WUE.
We also found that there were differences in the WUE among the different vegetation types. Specifically, the annual mean WUE of the forests was the highest, that of the cropland was the second highest, that of the grassland was the third highest, and that of the shrubland was the lowest. This was probably caused by the zonal differences in the vegetation, the hydrothermal conditions, and their distributions [32,50]. The forests were mostly composed of tall trees with high LAI and canopy cover values, intense photosynthesis in the leaves, and well-developed root systems of the trees. This could help the deep soil moisture to be absorbed and used for vegetation growth and to store organic matter, which is why it had the highest WUE [51]. The shrubland was the primary vegetation type planted in ecological restoration projects in arid and semi-arid areas, and the shrubland was located in areas with poor stand conditions. The insufficient soil water in these areas resulted in low vegetation coverage and LAI values, and the water was universally used for soil evaporation rather than vegetation photosynthesis. In addition, the vegetation of typical shrublands (e.g., sea-buckthorn and Caragana korshinskii) have small leaf sizes and low chlorophyll contents, making it difficult for remote sensing signals to truly reflect the vegetation cover of shorter shrubland, which resulted in a certain degree of underestimation of the GPP remote sensing products [52,53]. Consequently, the WUE of shrubland was relatively small. In previous studies, the phenomenon that forests could generate more organic matter than cropland and that the arbor WUE was higher than that of herbaceous plants were also observed by Xia et al. [54].
The growth rate of the WUE in the cropland exhibited a slightly decreasing trend as the drought index increased, whereas that of the other vegetation types exhibited increasing and then decreasing trends. Importantly, the cropland was influenced by anthropogenic irrigation activities, resulting in a more stable interannual variation in its GPP [55]. Irrigation activities in dry regions can still meet the cropland water consumption demand, so the sensitivity of the rate of cropland WUE to the drought index is low. The implementation of large-scale ecological conservation projects significantly contributed to the increase in the GPP in both the shrubland and grassland under different ecological conditions [56]. The transpiration of the vegetation was more sensitive to the soil moisture, but there was limited water available for vegetation evapotranspiration in the dry areas, and thus, the rate of increase of the ET slowed down. As a result, the rates of increase of the WUE trends in the shrubland and grassland increased.

Underlying Causes of the Effects of the Influencing Factors on the WUE Trend
The effect of the LAI on the change in the WUE was the largest for the entire Chinese Loess Plateau, the different vegetation types, and the different sub-basins from 2000 to 2018. We attributed this effect to the fact that the Chinese government has implemented ecological restoration projects on the CLP since 1999, such as planting trees and grass in the barren land, abandoning farmland and planting trees and grass, and prohibiting grazing on hills, which has led to a significant increase in the vegetation LAI [34]. The increase in the leaf area significantly increased the capacity of the vegetation to fix atmospheric CO2, and it also significantly increased the share of the vegetation transpiration and the trapped evaporation in the total evapotranspiration of the ecosystem, which led to a significant increase in the WUE.
Previous studies have concluded that in arid and semiarid areas, water is the most vital factor limiting the GPP [1]. An increasing trend in precipitation mitigates the water limitations and thus sustains the plant growth more efficiently, which improves the plant productivity and the WUE [57]. Additionally, since 2000, the precipitation in the CLP region has been relatively abundant, providing favorable climatic conditions for vegetation growth. Therefore, precipitation had a positive effect on the WUE in most areas. Nevertheless, we observed that the precipitation had a negative impact on the WUE in the northwestern part of the CLP, which was consistent with the results of previous studies [58]. We attributed this impact to the fact that the region was a severely arid region with scarce perennial precipitation. The evapotranspiration was high in the severely arid region, and soil moisture increased but evaporated rapidly after precipitation occurred, resulting in a decrease in the WUE as the total ET increased. It was difficult for the soil moisture to be used for vegetation photosynthesis. In addition, the prolonged drought stress altered the morphological, physiological, and biochemical traits to help vegetation adapt to drought [59]. As physiological responses, drought stress significantly decreased maximum photosynthesis rate and electron transfer rate [10]. The significant decrease of photosynthetic rate could make WUE decrease, and this may conform the potential efficiency of the cooling system. Because the shrubland was distributed mostly in the northwest, the increase in precipitation had a negative effect on the shrubland WUE. Meanwhile, drought could reduce the morphology of the plant, such as reduction of leaf size and vegetative growth, as well as reduce the plant vitality by reducing the photosynthesis and leaf water potential, and reduce the transpiration, vegetative growth, and leaf area of the plant. Because changes in the Temp, VPD, and SD affected the vegetation photosynthesis and transpiration by altering the leaf stomatal conductance, the mechanisms by which these three factors influenced the vegetation WUE were complex. In our study, the WUE was found to be positively correlated with temperature in most of the region. When the temperature increased, more water vapor evaporated from the bare soil into the atmosphere [1], and the transpiration was partly driven by the high VPD. This high evaporative and high VPDdriven transpiration water loss, which was decoupled from the vegetation productivity, led to a higher WUE [5]. Moreover, the moderate increase in the temperature also enhanced the activity of the enzymes required for leaf photosynthesis and thus the intensity of photosynthesis, which likewise promoted an increase in the WUE trends for the vegetation types [60]. It also had been shown that the increase in the SD allowed the vegetation to maximize the effective photosynthesis [61]. In general, the increases in Temp, VPD and SD all played roles in increasing the WUE trend of the vegetation ecosystem.
Note that the effects of the other factors (e.g., microtopographic changes, the ring-CO2, N deposition, and soil fertility enhancement) on the vegetation WUE were found to be non-negligible. To effectively suppress soil erosion on the CLP, soil and water conservation engineering measures such as the construction of terraces, silt dams, and horizontal steps were implemented on a large scale, leading to variations in the local microtopography. These variations could modify the precipitation redistribution process, making a relatively adequate water supply available for vegetation growth, which would result in a change in the vegetation WUE [8,31]. In addition, long-term accumulation of plant litter production and root development might effectively improve soil organic matter content and infiltration capacity [62]. It was also one of the essential factors affecting vegetation WUE. Numerous studies have shown that increased atmospheric CO2 concentration and nitrogen deposition contribute significantly to the increase in the vegetation GPP in China, thus increasing the vegetation WUE [4,63,64]. Future studies should further explore the effects of the influencing factors.

Uncertainties
First, we validated the remote sensing products of the PML_V2 using the eddy covariance-derived data from the flux towers and compared it with the data derived from the water balance equations. Because of the site data limitations, our validation process still had inevitable uncertainties, and further comparative optimization should be performed on more sites in the future. Second, the uncertainties in the attribution analysis also should be noted. On the one hand, the meteorological data were spatially interpolated from the actual data measured at the stations. On the other hand, when calculating the effect of each factor on the WUE trends, we assumed a linear response between the influencing factors and the WUE. The results of the multiple linear regression approach, the fitting coefficient (R 2 ), and the statistical significance are shown in Figure 11. The spatial average of the R 2 values was 0.6 on the CLP, and about 64.77% of the entire CLP region performed significantly (p < 0.05). About 72.37% of the entire CLP region performed significantly (p < 0.1). Actually, the meteorological factors were not completely independent, and their interactions were so complex that it may have confused their independent contributions to the WUE trends to some extent. Therefore, we need to consider combining nonlinear statistical models, remote sensing models, hydrological models, and other methods to conduct the attribution analysis and reduce the uncertainties caused by the attribution methods in the future.

Conclusions
The investigation of the spatio-temporal variation characteristics of WUE and its influencing factors was necessary for assessing the water consumption and sustainability of vegetation on the CLP. In this study, we analyzed the spatial characteristics and temporal evolution of the WUE and the contributions of the changes in the influencing factors for different vegetation types and different sub-basins in the CLP from 2000 to 2018, using the linear slope, the partial correlation coefficient method, and the total derivative method based on multisource remote sensing datasets. The main conclusions of this study are as follows: (1) The average annual WUE was 1.242 gC•m −2 •mm −1 , which decreased significantly with increasing drought index. The sensitivity of the WUE to drought index for the different vegetation types is ranked as follows: grassland > forest > cropland > shrubland. The annual growth rate of interannual WUE was 0.024 gC•m −2 •mm −1 yr −1 , with 87.32% of the study area showing an increasing trend. (2) The area in which the LAI was the most strongly correlated with the WUE accounted for 62.18% of the total area, and it was generally distributed in the central part of the study area, followed by the Temp, VPD, SD, and P. Compared with the cropland, the correlations between the WUE and the meteorological factors were more sensitive to the drought index in the forests and grassland. (3) The WUE trend was predominantly controlled by the LAI with an impact and relative impact rate of 0.014 gC•m −2 •mm −1 •yr −1 and 58.33%, respectively. The impacts of the LAI on the WUE trends for the different vegetation types and different sub-basins were ranked as follows: grassland > cropland > forest > shrubland, and Yan River > Wuding River > Kuye River > Huangpuchuan River > Jing River > Beiluo River > Fen River > Wei River, respectively. The impacts of the meteorological factors on the WUE trends were ranked as follows: Temp > VPD > SD > P. The integrated impact and relative impact rate of the meteorological factors on the WUE trends were 0.006 gC•m −2 •mm −1 •yr −1 and 25.00%, respectively. In addition, the comprehensive effects of the residual factors that were induced mainly by human activities on the WUE trends cannot be ignored in the CLP.
Funding: This work was financially supported by the Special Research Fund of the YRIHR (HKY-JBYW-2020-09), and National Natural Science Foundation of China (41701509, 51979118).

Data Availability Statement:
The data is available on request from the corresponding author.