Spatial and Temporal Characteristics of Water Use Efﬁciency in Typical Ecosystems on the Loess Plateau in the Last 20 Years, with Drivers and Implications for Ecological Restoration

: The water use efﬁciency (WUE) is an essential indicator of carbon–water coupling between terrestrial ecosystems and the atmosphere, and it is an important parameter for studying ecosystem responses to global climate change. A comprehensive understanding of the water–carbon coupling process in the Loess Plateau can reﬂect the balance between the “carbon absorption” and “water consumption” in vegetation, which drives the ecosystem succession process. In recent years, scholars have gained a more comprehensive understanding of the WUE and the driving factors of the Loess Plateau. However, there is still a need to study the carbon and water coupling mechanisms of different land use types in the Loess Plateau region. In this article, based on the gross primary productivity (GPP), evapotranspiration (ET), surface cover remote sensing products, and meteorological observation data, the trend of WUE changes for different vegetation types in the Loess Plateau from 2001 to 2020 and the correlations with the Normalized Difference Vegetation Index (NDVI), precipitation, and temperature values were analyzed using the Theil–Sen median (SEN) trend analysis method and correlation coefﬁcient analysis method. The spatial distribution patterns of the changes with the drought index showed that the multi-year average WUE value of the Loess Plateau was 1.24 g C mm − 1 H 2 O, and the mean WUE values in different seasons were ranked as follows: summer > autumn > spring. The WUE growth rates of all vegetation types showed a decreasing trend with the increase in drought index, and the size of the WUE response rate for each vegetation type to drought was ranked as follows: grassland > forest > shrub > crop. The annual average WUE increase rate of the Loess Plateau was 0.02 g C mm − 1 H 2 O yr − 1 , of which 93.36% of the area showed an increasing trend. The NDVI was the dominant factor affecting the spatial and temporal variations in WUE rates in the Loess Plateau, and the correlation between the NDVI and WUE was strongest in summer. In the more arid regional ecosystems, the WUE was negatively correlated with the precipitation and temperature, but in summer the precipitation had a positive effect on the WUE. The correlation of grassland and shrub WUE rates with temperature was more sensitive to the drought index than that of the forest and crop areas, but there was also a threshold effect. Therefore, when vegetation restoration is carried out in arid and semi-arid regions, the carbon and water coupling mechanisms of different vegetation types and the reasonable allocation of regional water resources should be fully considered.


Introduction
The ecosystem water use efficiency (WUE) is a measure of ecosystem carbon-water interactions and an important indicator of ecosystem sensitivity to climate change [1]. The WUE is commonly defined as the amount of dry matter fixed per unit mass of water consumed by terrestrial ecosystems, and there are many ways to express this WUE; usually, the pathways. However, in previous studies, it has been difficult to investigate the complex spatial heterogeneity of the Loess Plateau, meaning they could not investigate the carbonwater coupling mechanisms of different vegetation types. Therefore, in this paper, we study the carbon and water coupling mechanisms of different vegetation types among various ecosystems under water resource constraints in the Loess Plateau region, with complex spatial heterogeneity and a long time series, which will help clarify the driving forces of different carbon and water cycling mechanisms in ecosystem succession and provide a reference for the sustainable carbon sequestration potential of different vegetation types.

Study Area Introduction
The Loess Plateau (34~40 • N, 103~114 • E) is the largest loess accumulation area in the world and the largest loess-covered plateau in the world, located in the central-north of China, which belongs to the northern region. It mainly includes Shanxi and Shaanxi, as well as parts of Gansu, Qinghai, Ningxia, and Henan provinces, accounting for 70% of the distribution of loess in the world, and it is the largest loess accumulation area in the world [22]. The total area of the region is 63.5 × 104 km 2 , of which 45.4 × 104 km 2 is a soil erosion area, which is the most serious soil erosion area and the most fragile ecological environment in China and even in the world ( Figure 1). The region belongs to a typical arid and semi-arid continental monsoon climate zone with an average annual temperature range of 3.6~14.3 • C. Winter and spring are influenced by dry and cold polar air masses, which are cold, dry, and windy; summer and autumn are hot and humid, with an average multi-year rainfall rate of 466 mm, mainly concentrated between June and September, accounting for about 60% of the annual precipitation [23][24][25]. The winter precipitation generally accounts for only about 5% of the total. The vegetation in the study area shows a general trend of forest to grassland overgrowth from south to north. Most of the central part is a semi-arid grassland zone. In the northwestern part, the landscape gradually evolves into a desert, with desert grassland dominating.
dominates in tropical forest zones [21]. Therefore, the CWSI plays a role in the WUE, in GPP and ET pathways. However, in previous studies, it has been difficult to invest the complex spatial heterogeneity of the Loess Plateau, meaning they could not inv gate the carbon-water coupling mechanisms of different vegetation types. Therefor this paper, we study the carbon and water coupling mechanisms of different vegeta types among various ecosystems under water resource constraints in the Loess Pla region, with complex spatial heterogeneity and a long time series, which will help cl the driving forces of different carbon and water cycling mechanisms in ecosystem su sion and provide a reference for the sustainable carbon sequestration potential of diffe vegetation types.

Study Area Introduction
The Loess Plateau (34~40°N, 103~114°E) is the largest loess accumulation area in world and the largest loess-covered plateau in the world, located in the central-nor China, which belongs to the northern region. It mainly includes Shanxi and Shaanx well as parts of Gansu, Qinghai, Ningxia, and Henan provinces, accounting for 70% o distribution of loess in the world, and it is the largest loess accumulation area in the w [22]. The total area of the region is 63.5 × 104 km 2 , of which 45.4 × 104 km 2 is a soil ero area, which is the most serious soil erosion area and the most fragile ecological env ment in China and even in the world ( Figure 1). The region belongs to a typical arid semi-arid continental monsoon climate zone with an average annual temperature r of 3.6~14.3 °C. Winter and spring are influenced by dry and cold polar air masses, w are cold, dry, and windy; summer and autumn are hot and humid, with an average m year rainfall rate of 466 mm, mainly concentrated between June and September, acco ing for about 60% of the annual precipitation [23][24][25]. The winter precipitation gene accounts for only about 5% of the total. The vegetation in the study area shows a gen trend of forest to grassland overgrowth from south to north. Most of the central par semi-arid grassland zone. In the northwestern part, the landscape gradually evolves a desert, with desert grassland dominating.

Data Sources
The carbon cycle, water cycle, vegetation succession, and meteorological data for the Loess Plateau vegetation covered in this paper were obtained from the GEE selfcontained data products. The GPP and ET data were derived from the PML_V2v0.1.7 data product (https://developers.google.com/earth-engine/datasets/catalog/CAS_IGSNRR_ PML_V2_v017, accessed on 18 July 2021) to express ecosystem carbon uptake and water depletion rates [26]. The NDVI data were derived from the MOD13A2V6 product (https:// developers.google.com/earth-engine/datasets/catalog/MODIS_006_MOD13A1, accessed on 18 July 2021) to express vegetation growth rates with an inter-data resolution of 500 m and a temporal resolution of 8 d. The land cover type data were derived from the GCL_FC30 product (http://data.casearth.cn, accessed on 18 July 2021) of the Institute of Space and Space Innovation, Chinese Academy of Sciences [8], which was derived from Landsat ground surface reflection data; a Chinese land cover product with a spatial resolution of 30 m was obtained using supervised classification. The meteorological data were obtained from the China Meteorological Science Data Center (http://data.cma.cn/, accessed on 18 July 2021). Fifty-two meteorological stations within the Loess Plateau involving the seven provinces of Shanxi, Qinghai, Ningxia, Inner Mongolia, Henan, Gansu, and Shaanxi were used, including precipitation and mean temperature data. The data used in this paper all range in time from 2000 to 2020, and the specific information is shown in Table 1. Currently, there are different definitions and understandings of vegetation WUE in academia, and the results are expressed in different ways. The calculation of the WUE using the ratio of the GPP to ET is the most classical approach [27]. The calculation method is as follows: where WUE is the water use efficiency (g C mm −1 H 2 O); GPP is the total primary productivity of the terrestrial ecosystems (g C-m −2 ); ET is the ecosystem evapotranspiration (mm-m −2 ).

Calculation of Drought Index
The crop water stress index (CWSI) has a clear physical meaning and is suitable for studying the degree of vegetation drought in semi-arid areas, which is based on the principle of surface balance and takes into account the vegetation cover rates of different substrates and meteorological elements, such as precipitation and temperature, and has been widely used in studies of agricultural drought [3]. The vegetation water deficit index (CWSI) is calculated from evapotranspiration data with the following equation: where PET indicates the potential evapotranspiration in mm and ET indicates the actual evapotranspiration in mm. The CWSI values range from 0 to 1, with larger values representing areas under greater water stress, and vice versa. The land use transfer matrix may effectively uncover and depict the structural aspects of land cover change, as well as the direction of change for various use types, and reflect the special characteristics of the process of land use type transfer change [28]. The transfer matrix equation is shown as follows: where, P ij is the area of land use type i converted to type j land use type in different periods.

Trend Analysis
A trend analysis is a method of statistically testing trends in time series [29]. In this study, a combination of the Theil-Sen median (SEN) and Mann-Kendall (MK) statistical tests was used for a further analysis of the trends in environmental indicators in the study area. The formula is: where Median() represents the median value; X j and X i represents the j-th and i-th element; if β is greater than zero, this means that the time series shows an upward trend, and vice versa for a downward trend.

Correlation Analysis
In this study, the Pearson coefficients were used to analyze the correlations of each factor, which were calculated as shown below [27,30]: where R is the Pearson coefficient; x is the influence factor; and y is the dependent variable. The magnitude of R reflects the correlation between the factors, where R > 0 indicates a positive correlation and R < 0 indicates a negative correlation.

Vegetation Cover Changes from 2001 to 2020
In this study, using the land use data for the Loess Plateau in 2001 and 2020, it can be seen from Figure 2 that the land cover types of the Loess Plateau show subtle spatial differences. The whole region shows a general trend of overgrowth of forest to grassland areas from south to north, with forest areas being concentrated in the valley plains of the Guanzhong Basin, Fen River Valley, and Luoyang Basin (Figure 2), accounting for 12.07% of the total area of the Loess Plateau in 2001 and increasing to 14.28% in 2020 (Table 2). Grassland is the most abundant type of land cover category, mainly distributed in the desert areas in the north ( Figure 2) and increasing from 49.12 to 50.24% in 2001. The arable land is mainly distributed in the loess areas of the Longdong Plateau and the plateau areas on both sides of the Weihe River in the south of the Loess Plateau ( Figure 2). Since the project of "returning crop to forest and grass areas" started, the area of arable land has been reduced from 31.73 to 28.83% ( Table 2). The area covered by shrub is the smallest, which has been stable at less than 1% for many years and shows a decreasing trend ( Table 2). land is mainly distributed in the loess areas of the Longdong Plateau and the plateau areas on both sides of the Weihe River in the south of the Loess Plateau ( Figure 2). Since the project of "returning crop to forest and grass areas" started, the area of arable land has been reduced from 31.73 to 28.83% ( Table 2). The area covered by shrub is the smallest, which has been stable at less than 1% for many years and shows a decreasing trend ( Table  2).  Although the areas of these four land cover types remained relatively stable from 2001 to 2020 (Table 2), the NDVI values showed an increasing trend (Figure 3a), and the vegetation cover of forest and scrub areas showed significant advantages compared with the other two classes ( Figure 3c). Spatially, the NDVI values showed a significant increasing trend in most areas of the Loess Plateau, especially in the central-eastern part of the plateau (Figure 3b), where grassland is the main growth area. After counting the NDVI trends for different land cover types, it can be found that grassland has the fastest growth rate, followed by crop, while the slowest growth rate is for scrub (Figure 3d).  Although the areas of these four land cover types remained relatively stable from 2001 to 2020 (Table 2), the NDVI values showed an increasing trend (Figure 3a), and the vegetation cover of forest and scrub areas showed significant advantages compared with the other two classes ( Figure 3c). Spatially, the NDVI values showed a significant increasing trend in most areas of the Loess Plateau, especially in the central-eastern part of the plateau (Figure 3b), where grassland is the main growth area. After counting the NDVI trends for different land cover types, it can be found that grassland has the fastest growth rate, followed by crop, while the slowest growth rate is for scrub (Figure 3d).

Seasonal Distribution Characteristics of WUE in the Loess Plateau
The spatial distribution of the annual average WUE rates for the Loess Plateau from 2001 to 2020 shows a clear spatial differentiation (Figure 4). The 20-year WUE average is 1.24 g C mm −1 H 2 O, which indicates that the loss of 1 mm of water through evapotranspiration in the Loess Plateau represents about 1.24 g of CO 2 fixed by vegetation. In Shanxi, southern Shaanxi, and Henan provinces, the WUE values can be above 1.5 g C mm −1 H 2 O, mainly in cultivated and forested areas; the low-value areas were mainly in the Inner Mongolia Autonomous Region, Ningxia Hui Autonomous Region, Gansu Province, and Qinghai Province. The spatial distribution of WUE in different seasons over the years showed that the distribution of the WUE rates in different seasons showed the same distribution trend as that of the WUE rates over the whole year, while the average values from high to low were: summer (1.

Seasonal Distribution Characteristics of WUE in the Loess Plateau
The spatial distribution of the annual average WUE rates for the Loess Plateau from 2001 to 2020 shows a clear spatial differentiation (Figure 4). The 20-year WUE average is 1.24 g C mm −1 H2O, which indicates that the loss of 1 mm of water through evapotranspiration in the Loess Plateau represents about 1.24 g of CO2 fixed by vegetation. In Shanxi, southern Shaanxi, and Henan provinces, the WUE values can be above 1.5 g C mm −1 H2O, mainly in cultivated and forested areas; the low-value areas were mainly in the Inner Mongolia Autonomous Region, Ningxia Hui Autonomous Region, Gansu Province, and Qinghai Province. The spatial distribution of WUE in different seasons over the years showed that the distribution of the WUE rates in different seasons showed the same distribution trend as that of the WUE rates over the whole year, while the average values from high to low were: summer (1.

Seasonal Variation Characteristics of WUE Rates in the Loess Plateau
The spatial distribution of the interannual variation of the WUE rates in the Loess Plateau from 2001 to 2020 is shown in Figure 6. The WUE of the entire Loess Plateau is increasing at the rate of 0.02 g C mm −1 H 2 O per year, and the areas with increasing trends account for 93.36% of the Loess Plateau (Figure 6a Figure 5. The data are averaged over the entire Loess Plateau region. From the figure, it can be seen that the WUE rates of the different vegetation types in descending order were forest > scrub > crop > grassland, with forest areas reaching 1.72 g C mm −1 H2O (Figure 5a). This paper also found that the WUE rates of all vegetation types decreased with the increase in CWSI, while the response rates of the WUE for each vegetation type to the increase in CWSI were ranked as grassland > crop > forest > shrub (Figure 5a). The rate of decline for the grassland WUE with the increase in CWSI was the fastest, indicating that the grassland WUE was the most sensitive to the changes in drought patterns. From the change pattern of the WUE rates with the drought index rates in different seasons, it can be seen that summer had the same response pattern as the whole-year pattern, while crop was the most sensitive to drought change responses in spring ( Figure 5b); forest was the most sensitive to drought change responses in autumn ( Figure 5d); and in both spring and autumn, grassland showed lower sensitivity to drought index changes (Figure 5b,d).  Figure 5. The data are averaged over the entire Loess Plateau region. From the figure, it can be seen that the WUE rates of the different vegetation types in descending order were forest > scrub > crop > grassland, with forest areas reaching 1.72 g C mm −1 H 2 O (Figure 5a). This paper also found that the WUE rates of all vegetation types decreased with the increase in CWSI, while the response rates of the WUE for each vegetation type to the increase in CWSI were ranked as grassland > crop > forest > shrub (Figure 5a). The rate of decline for the grassland WUE with the increase in CWSI was the fastest, indicating that the grassland WUE was the most sensitive to the changes in drought patterns. From the change pattern of the WUE rates with the drought index rates in different seasons, it can be seen that summer had the same response pattern as the whole-year pattern, while crop was the most sensitive to drought change responses in spring ( Figure 5b); forest was the most sensitive to drought change responses in autumn ( Figure 5d); and in both spring and autumn, grassland showed lower sensitivity to drought index changes (Figure 5b,d).

Seasonal Variation Characteristics of WUE Rates in the Loess Plateau
The spatial distribution of the interannual variation of the WUE rates in the Loess Plateau from 2001 to 2020 is shown in Figure 6. The WUE of the entire Loess Plateau is increasing at the rate of 0.02 g C mm −1 H2O per year, and the areas with increasing trends account for 93.36% of the Loess Plateau (Figure 6a). While the increasing areas are mainly distributed in the southern and central parts of the Loess Plateau, the decreasing areas are mainly distributed in the northwestern part of the Loess Plateau. In terms of the seasonal distribution, the WUE changes are slowest in spring, increasing at 0.01g C mm −1 H2O per year (Figure 6b), with 12.63% of the area showing a decreasing trend, especially in the western part of the Loess Plateau in Qinghai Province. The summer WUE increased at a rate of 0.03 g C mm −1 H2O per year, which was the fastest growth season of the year, and 18.94% of the regions showed an annual increase of more than 0.06 g C mm −1 H2O ( Figure  6c) but showed a more obvious decreasing trend in the southern part of the Loess Plateau. Autumn showed a similar distribution pattern as summer, with 86.25% of the area showing an increasing trend (Figure 6d). The variation curves of interannual WUE trends with the drought index rates for different vegetation types are shown in Figure 7. For the entire Loess Plateau, although the WUE rates of different vegetation types showed an increasing trend in the CWSI in the interval of 0.2-0.5, there are still large differences in the variation patterns of WUE growth rates with the drought index values for the different vegetation types. In order to consider the completeness of image information, the interval of CWSI in this paper starts from 0. However, when CWSI is less than 0.3, no artificial vegetation is subjected to water stress, so the analysis results are based on CWSI of 0.3 as the starting point. The WUE growth rate of agricultural land decreased with the increase in the drought index; however, the forest, shrub, and grassland areas all showed a pattern of increasing and then decreasing, and the peak of the WUE occurred in the CWSI range of 0.3-0.4 (Figure 7a). There were still large differences in WUE trends with drought index rates in different seasons. The rate of change of the WUE in crop areas was the largest in spring, increasing at 0.02 g C mm −1 H 2 O per year (Figure 7b), and fluctuated the least with the increase in drought index, while the rate of change for grassland areas showed an increasing trend the whole time with the increase in drought index. The summer grassland showed an increasing trend at a rate of 0.03 g C mm −1 H 2 O per year (Figure 7c) and was also the fastest growing among vegetation types. The pattern of WUE changes in autumn was similar to that for the whole year (Figure 7d). It is noteworthy that although the sensitivity of vegetation WUE to CWSI differs, there is an effect with it, which shows an increasing trend of vegetation WUE between CWSI interval of 0.2 and 0.6. The variation curves of interannual WUE trends with the drought index rates for different vegetation types are shown in Figure 7. For the entire Loess Plateau, although the WUE rates of different vegetation types showed an increasing trend in the CWSI in the interval of 0.2-0.5, there are still large differences in the variation patterns of WUE growth rates with the drought index values for the different vegetation types. In order to consider the completeness of image information, the interval of CWSI in this paper starts from 0. However, when CWSI is less than 0.3, no artificial vegetation is subjected to water stress, so the analysis results are based on CWSI of 0.3 as the starting point. The WUE growth rate of agricultural land decreased with the increase in the drought index; however, the forest, shrub, and grassland areas all showed a pattern of increasing and then decreasing, and the peak of the WUE occurred in the CWSI range of 0.3-0.4 (Figure 7a). There were still large differences in WUE trends with drought index rates in different seasons. The rate of change of the WUE in crop areas was the largest in spring, increasing at 0.02 g C mm −1 H2O per year (Figure 7b), and fluctuated the least with the increase in drought index, while the rate of change for grassland areas showed an increasing trend the whole time with the increase in drought index. The summer grassland showed an increasing trend at a rate of 0.03 g C mm −1 H2O per year (Figure 7c) and was also the fastest

Spatial Distribution Characteristics of WUE and NDVI Rates, Precipitation, and Temperature Correlations in the Loess Plateau
The correlation coefficients of the Loess Plateau with NDVI, precipitation, and temperature values are shown in Figure 8. The figure shows that there are significant differences in the geographical distribution of the correlation between different influencing factors and WUE. In terms of the relationship between the vegetation WUE and individual factors in the study area, the vegetation WUE showed a strong positive correlation with the NDVI (Figure 8a), and the correlation coefficient between the NDVI and WUE was 0.71 in 93.51% of the areas that passed the significance test (Table 3). For the correlation between the NDVI and WUE in different seasons, the NDVI was most significantly correlated with the WUE in summer, with 91.63% of the areas passing the correlation test of p < 0.05 (Figure 8g) and with a correlation coefficient of 0.46 (Table 3). Only 23.84% of the areas passed the significance test in spring compared to the other two seasons, and these areas were mainly located in the southern part of the Loess Plateau (Figure 8d). The correlation in these areas was only 0.23 (Table 3).

Spatial Distribution Characteristics of WUE and NDVI Rates, Precipitation, and Temperature Correlations in the Loess Plateau
The correlation coefficients of the Loess Plateau with NDVI, precipitation, and temperature values are shown in Figure 8. The figure shows that there are significant differences in the geographical distribution of the correlation between different influencing factors and WUE. In terms of the relationship between the vegetation WUE and individual factors in the study area, the vegetation WUE showed a strong positive correlation with the NDVI (Figure 8a), and the correlation coefficient between the NDVI and WUE was 0.71 in 93.51% of the areas that passed the significance test (Table 3). For the correlation between the NDVI and WUE in different seasons, the NDVI was most significantly correlated with the WUE in summer, with 91.63% of the areas passing the correlation test of p < 0.05 (Figure 8g) and with a correlation coefficient of 0.46 (Table 3). Only 23.84% of the areas passed the significance test in spring compared to the other two seasons, and these areas were mainly located in the southern part of the Loess Plateau (Figure 8d). The correlation in these areas was only 0.23 (Table 3). The correlation between the WUE and precipitation was mildly positive, the area passing the significance test was only 18.26% of the Loess Plateau area (Figure 8b), and the correlation coefficient was 0.28 (Table 3), but there were still some areas showing a negative correlation, which were mainly distributed in a small part of Qinghai Province in the western part of the Loess Plateau. There were great differences in the correlation between precipitation and WUE in different seasons, and the precipitation and WUE values showed a stronger negative correlation in spring and autumn (Figure 8e,k), with the difference that 29.37% more regions passed the significance test in spring than in autumn. The summer precipitation-WUE correlation was most like the distribution pattern of correlations throughout the year, with both showing positive correlations, but the correlation between the summer precipitation and WUE was more significant. Remote Sens. 2022, 14, x FOR PEER REVIEW 12 of 18   The correlation between temperature and WUE was the weakest compared with the first two, indicating that the temperature had the least influence on the WUE in the Loess Plateau. The annual correlation between the temperature and WUE was 0.36 (Table 3), and only 9.82% of the regions passed the significance test, which were mainly located in parts of Inner Mongolia in the northern part of the Loess Plateau and parts of Shaanxi Province in the south (Figure 8c). For the different seasons, the correlation between the temperature and WUE still showed more obvious differences, and spring was the strongest correlation between the temperature and WUE among the three seasons, with 0.43 (Table 3) and 32.78% of the regions passing the significance test (Figure 8f). Only the summer temperature showed a non-significant negative correlation with the WUE (Figure 8i). The regions were mainly located in Shanxi Province and parts of Inner Mongolia. The effect of the autumn temperature on the WUE was the smallest, with a correlation of only 0.16 (Table 3).

Spatial Distribution Characteristics of Spatial and Temporal Variation of the WUE in the Loess Plateau Ecosystem
The ecosystem WUE takes into account the relationship between the vegetation water consumption and dry material production, and it is an extensive physiological and ecological index used to assess the suitability of plant growth. In the water-scarce Loess Plateau region, it is important to clarify the distribution of the ecosystem WUE under different drought indices, as well as the changes and responses to the influencing factors. This paper investigates the seasonal distribution and changes of the WUE rates of different vegetation types on the Loess Plateau during 2001 to 2020. It was found that the WUE of the Loess Plateau showed obvious spatial heterogeneity in space, i.e., decreasing gradually from southeast to northwest [31,32] (Figure 4). Seasonally, the WUE of the vegetation started to increase in spring when the temperature started to pick up, the precipitation and sunshine hours kept increasing, and the surface plants entered the growth period [33]. The average temperature in summer was higher than 22 • C, which was suitable for vegetation growth, with sufficient precipitation and the longest sunshine hours, meaning vegetation growth was vigorous [34], so the vegetation WUE was higher relative to other seasons. Spatially, it was high-value areas concentrated in the forest in the south and east of Loess Plateau. In autumn, the temperature started to drop, the sunshine hours gradually shortened, the vegetation gradually withered [35], and the vegetation WUE started to retract. The autumn WUE distribution was more consistent with the annual average WUE, showing low rates in the northwest and high rates in the southeast, and the high-value areas were mainly distributed in Shanxi Province and southern Shaanxi. The southeastern part of the Loess Plateau cools down slowly in autumn, which is more favorable for vegetation growth. Figure 5 shows that the highest WUE was found in the forests of the Loess Plateau, which was due to the fact that most of the forests comprise tree vegetation with larger leaves which have higher vegetation cover, stronger photosynthesis potential, and a welldeveloped root system, making them more capable of absorbing water from the soil compared to other vegetation types [36,37]. In spring, when the temperature rises and the vegetation starts to grow, slight droughts do not have much effect on water utilization in forests, so the forest WUE changes are less influenced by the drought index. Crop is more affected by anthropogenic irrigation activities [38], but crop in arid areas can still meet the water consumption needs of crops through anthropogenic irrigation, so this land cover type is least sensitive compared to other concentrated vegetation types. Summer is the season when crop fixes the most carbon, while evaporation is less than for other small vegetation types [39], so the WUE is higher in areas with crop distribution. The grasslands have a high adaptive capacity and are also the vegetation type comprising the largest area on the Loess Plateau, but due to the short root system of the grasslands [21], precipitation has become the main limiting factor, and the grasslands are subject to more severe water stress compared to other vegetation types. The hot temperatures in summer make the grassland WUE distribution more sensitive to drought indices. Figure 9 shows that the correlations between the different vegetation types and impact factors on the Loess Plateau show different patterns of change with the drought index values. The correlation between the ecosystem WUE and NDVI with increasing drought index values (Figure 9a) has the same variation pattern as for the growth trend of the WUE (Figure 7a), and both reach the peak in the interval of the CWSI of 0.3-0.4, which indicates that the NDVI dominates the spatial pattern of the interannual variation in WUE rates in the Loess Plateau. For the different seasonal variations, the correlation between the WUE and NDVI for the ecosystems showed a trend of decreasing and then increasing with increasing drought index values in spring (Figure 9d), and the drought index values were primarily affected by grassland and brush. The link between the WUE and NDVI for each type of plant in the summer revealed a trend of rising and then declining with rising values for the drought index because the proportion of water directly used for the photosynthesis of vegetation to total evapotranspiration gradually decreased with the increasing drought index values [40,41], so the correlation between the ecosystem WUE and NDVI began to decrease. The correlation between the rainfall and ecosystem WUE in the Loess Plateau showed different degrees of negative correlations in the more arid regions in the north (Figure 7b), which was due to the fact that this region is mostly covered with barren grass. The soil water content increases after rainfall but evaporates rapidly, resulting in an increase in total evapotranspiration and difficulty in using soil water for its own photosynthesis [42], so the WUE decreases. The temperature showed a non-significant positive correlation with the ecosystem WUE, and since the temperature affects the stomatal conductance and vegetation photosynthesis mainly through evapotranspiration [43], the correlation between the WUE and temperature is stronger in forests than in several other vegetation types. The negative correlation between the temperature and ecosystem WUE in summer is due to the low precipitation and uneven water heating in the Loess Plateau. As the drought index increases and the temperature rises, this will accelerate the evapotranspiration of the vegetation and reduce the WUE, especially in the scrub, and as the drought index increases, the stronger the negative correlation will be between the temperature and the WUE (Figure 9i). Remote Sens. 2022, 14, x FOR PEER REVIEW 15 of 18

Conclusions
In this paper, we analyzed the response of the WUE changes to regional climate and vegetation cover changes during the annual period of 2001 to 2020 using the PML GPP and ET datasets. Our results showed that the multi-year mean WUE value of the Loess Plateau is 1.24 g C mm −1 H 2 O, and the spatial pattern of the WUE varies among seasons, with the mean WUE size being ranked as summer > autumn > spring in the different seasons. The WUE rates of all vegetation types showed a decreasing trend with the increase in drought index values, and it was most obvious in summer. The magnitude of the response rate of the WUE to drought for each vegetation type was ranked as: grassland > forest > scrub > crop.
The annual average WUE increase rate for the Loess Plateau was 0.02 g C mm −1 H 2 O a −1 , of which 93.36% of the area showed an increasing trend. Due to the difference in precipitation temperature in the different seasons, this led to the changes in WUE rates in different seasons still showing a big difference, and the mean WUE increase in summer was more significant. The NDVI was the dominant factor affecting the spatial and temporal variations in WUE in the Loess Plateau. The correlation between the NDVI and WUE was the strongest in summer. The change pattern of the NDVI values with the increase in drought index rates showed the same change pattern as for the growth trend of the WUE. The ecosystem WUE was negatively correlated with the precipitation temperature in drier regions, but throughout the summer, the precipitation boosted the WUE. In comparison to forest and crop regions, the link between the WUE and temperature in grassland and scrub areas was more sensitive to the values of the drought index, but there was also a threshold effect. Thus, for the arid and semi-arid Loess Plateau regions, these findings deepen our understanding of the carbon and water coupling mechanisms of different vegetation types, providing important ideas for the sustainable development of vegetation on the Loess Plateau and future ecological restoration management and governance approaches.