Spatial variations and mechanisms for the stability of water use efficiency in China

A clearer understanding of the stability of water use efficiency (WUE) and its driving factors contributes to improving water use efficiency and strengthening water resource management. However, the stability of WUE is unclear. Based on the EEMD method, this study analyses the spatial variations and mechanisms for the stability of WUE in China, especially in the National Forest Protection Project (NFPP) areas. It is found that the stable WUE was dominated by non-significant trends and increasing trends in China, accounting for 33.59% and 34.19%, respectively. The non-significant trend of stable WUE was mainly located in the Three-North shelterbelt program area, and the increasing trend of stable WUE was in Huaihe and Taihu, Taihang Mountains, and Pearl River shelterbelt program areas. Precipitation and soil moisture promoted the stable WUE in these project areas. The unstable WUE was dominated by positive reversals or negative reversals of WUE trends. The positive reversals of unstable WUE were mainly located in the Yellow River shelterbelt program areas, which was promoted by temperature and radiation, while the negative reversals of unstable WUE were mainly distributed in the Yangtze River and Liaohe shelterbelt program areas, which were mainly induced by saturation water vapor pressure difference (VPD). Our results highlight that some ecological restoration programs need to be improved to cope with the negative climate impact on the stability of WUE.

1 Introduction

Water use efficiency (WUE) is an objective evaluation index of water-carbon coupling for an ecosystem, which is defined as the ratio of carbon sequences (i.e., gross primary production (GPP)) to water loss (i.e., evapotranspiration (ET)) (Beer et al., 2009; Jiang et al., 2019; Ma J. et al., 2019; Liu et al., 2020). WUE is not only an important index of the coupled terrestrial carbon-water cycle, but also one of the important parameters reflecting the impact of global changes on terrestrial ecosystems (Hatfield and Dold, 2019; Yang L. et al., 2022). Uncovering the changes in the stability of the ecosystem WUE can provide important technical and policy implications for water conservation and carbon budgets (Cao et al., 2020). A large interannual variable in water use efficiency can lead to the instability of ecosystem functions and pose a serious challenge to nature-based climate solutions (Xue et al., 2015; Zhang et al., 2016; Wang et al., 2022). Thus, spatial variations and mechanisms for water use efficiency stability have become a foregrounded and topical issue for ecosystems.

The interannual variables of WUE are stable in the absence of climate and human drivers of change (Rödenbeck et al., 2018; Piao et al., 2020). However, the driving factors and their interactive effects influence GPP and ET in different ways, making the interannual variables in WUE threaten the stability of WUE (Ma J. et al., 2019). The instability in WUE considerably challenges the sustainability of the carbon-water cycle (Yu et al., 2021). The essence of WUE instability is that WUE has weak resistance to climate fluctuation (Wang et al., 2022). Both climate and ecosystems have obvious spatial differences, thus, WUE trends are nonlinear and their stability has a significant spatial variation (Piao et al., 2020). In areas with severe climate fluctuations or where WUE is sensitive to climate fluctuations, the stability of WUE is usually low (Liu et al., 2017; Ma J. et al., 2019). It provides a rare opportunity to explore the response of ecosystem functioning to climate change (Piao et al., 2020). However, there still exists uncertainty on the WUE stability response to climate change (Belmecheri et al., 2021). Thus, reducing the uncertainty is critical for an accurate future carbon-water cycle and its response to climate change (Belmecheri et al., 2021; Migliavacca et al., 2021).

Previous studies have shown that WUE change is influenced by climate change, such as CO₂, solar radiation, temperature, precipitation, saturated water vapor pressure, and soil moisture (Liu et al., 2020; Gonsamo et al., 2021; Wang M. et al., 2021). Precipitation is one of the principal indirect driving factors affecting WUE variables (Zhang et al., 2020). Precipitation changes directly affect the transpiration and evaporation of the ecosystem, and indirectly affect the carbon uptake process of plants by regulating the soil water content (Piao et al., 2020; Zhang et al., 2020). In addition, drought-induced reduction in vegetation production and WUE was offset by driving factors, such as warming climate and ecological restoration projects (Huang J. et al., 2015; Zhang et al., 2016; Ma J. et al., 2019). Numerous studies have shown that vapor pressure difference (VPD) is a key factor and has been demonstrated to have negative influences on the WUE (Beer et al., 2009; Ma J. et al., 2019). It is revealed that the increase in CO₂ concentration will increase the photosynthetic rate, while the transpiration rate will be weakened or have no significant effect, which will lead to the increase of vegetation WUE, with a significant CO₂ fertilization effect (El Masri et al., 2019; Ma J. et al., 2019; Gonsamo et al., 2021). Both photosynthesis and transpiration are affected by temperature, and show an opposite trend with increasing temperature (Huang J. et al., 2015). When the temperature was low, the photosynthetic rate increased with the increase of temperature, and then gradually weakened after it reached the maximum (Hatfield and Dold, 2019). The main reason was that the enzyme activity was significantly affected when the optimum temperature was reached (Hatfield and Dold, 2019). On the contrary, the increase in temperature will cause an increase in VPD and then increase the transpiration rate of vegetation (Cao et al., 2020). The radiation is one of the important factors of plant photosynthesis, and it also impacts plant WUE (Jiang et al., 2019; Wang et al., 2020). Nevertheless, the driving mechanisms underlying the stability of WUE remain unclear.

The National Forest Protection Project (NFPP) in China is one of the world’s largest ecological restoration projects (Tong et al., 2017; Lu et al., 2018; Huang et al., 2019). It is a major initiative by the Chinese government to mitigate some of the environmental damage caused by rapid economic development through the implementation of a rigorous and creative policy of large-scale conservation (Cai et al., 2020; Ding et al., 2021; Xu et al., 2022a). Currently, some researchers believe that the implementation of the NFPP has led to an increase in vegetation cover in some areas, thus improving the carbon and water cycle and the ecological environment (Tong et al., 2018; Huang et al., 2019; Yang Y. et al., 2022). Another group of scholars believe that the implementation of the NFPP, especially afforestation in some areas, has led to a reduction in soil moisture due to strong forest transpiration and the increased precipitation cannot compensate for the consumption of evapotranspiration, which places a burden on local water resources (Yang et al., 2013; Hai et al., 2022). Therefore, there is an urgent need to reveal the actual trend and stability of WUE and its driving factors in NFPP areas and to provide scientific advice for the implementation of ecological projects.

In this study, we investigated the spatial variations and mechanisms for the stability of water use efficiency in China. Specifically, we aimed to propose three key issues: (1) What are the nonlinear trends of WUE in China, especially in NFPP areas? (2) What is the stability of WUE in China over recent years? (3) Which driving factors are important in determining the stability of WUE in China?

2 Materials and methods

2.1 Data sources

The terrestrial gross primary production (GPP) data was downloaded from the Global Land Surface Satellite (GLASS) program (http://www.resdc.cn/) (Jiao et al., 2022). It is generated using the Bayesian algorithm ensemble of eight widely-used light-use efficiency models and has been widely used in global carbon cycle assessment (Xu et al., 2022a). The dataset is a global composite product that spans from 1982 to 2015 with a spatial resolution of 0.05°.

The terrestrial evapotranspiration (ET) data was downloaded by the National Science & Technology Infrastructure (http://www.nesdc.org.cn/), with a 0.1° spatial resolution from 1981–2015. To obtain more accurate data, ET data is simulated by a nonlinear complementary Relational model and verified with 13 vorticity covariance measurements and 10 river basin Nash Sutcliffe efficiency measurements, with a range of 0.72-0.94 (Ma N. et al., 2019).

The WUE was calculated as follows:

$\begin{array}{ll}\text{WUE}=\text{GPP}/\text{ET}& (1)\end{array}$

In this study, six driving factors were used to reveal the driving mechanism of the stability of WUE. The soil moisture (SM) and downward shortwave radiation (RAD) were downloaded by the Climatology Lab from 1981 to 2015 with a spatial resolution of 1/24°. Temperature and precipitation were downloaded from WorldCom (https://www.worldclim.org/) from 1981 to 2015 with a spatial resolution of 1 km. Vapor pressure deficit (VPD) was provided by Xu et al. (2021) (http://dx.doi.org/10.1016/j.scitotenv.2022.155086) from 1981 to 2015 with a spatial resolution of 0.1°. CO₂ was downloaded from the Emissions Database for Global Atmospheric Research (EDGAR) (https://edgar.jrc.ec.europa.eu/) from 1981 to 2015 with a spatial resolution of 0.1°.

2.2 Nonlinear method

The linear trends of NEP show an increasing trend with a constant rate of increase. However, the rate of NEP decomposed by EEMD increases with time. Therefore, the EEMD method can reveal the nonlinear trends in WUE (Pan et al., 2018; Xu et al., 2022b). Ensemble Empirical Mode Decomposition (EEMD) is an extension of the Empirical Mode Decomposition (EMD) method (Wu et al., 2007). The decomposition process of EMD is as follows:

First, the discrete extreme points of X(t) are interpolated to the entire period with a cubic spline function to obtain the maximum and minimum envelopes, and the arithmetic mean of the upper and lower envelopes (g₁(t)) is calculated.

$\begin{array}{ll}{g}_1(t)=X(t)-m_1(t)& (2)\end{array}$

Since g₁(t) is not stable, continue to repeat the above steps:

$\begin{array}{ll}{g}_{11}(t)=g_1(t)-m_{11}(t)& (3)\end{array}$

If the standard deviation (SD) is less than a given value (usually 0.2), the above iterative process is terminated:

$\begin{array}{ll}SD={\displaystyle \sum _{t=1}^N\frac{|g_{1k}(t)-g_{1(k-1)}(t){|}^2}{g^2{}_{1(k-1)}(t)}}& (4)\end{array}$

In this way, we extracted the first Intrinsic Mode Function IMF (IMF) component (c₁=g_1k(t)) from the original data:

$\begin{array}{ll}{g}_{1k}(t)=g_{1\left(k-1\right)}(t)-m_{1k}(t)& (5)\end{array}$

k is the number of iterations and the rest of the original data is:

$\begin{array}{ll}{r}_1=X(t)-c_1& (6)\end{array}$

Since r₁ still contains fluctuations of a longer period, the above iterative process is still repeated, and the r_i is as follows:

$\begin{array}{ll}{r}_i=r_{i-1}-c_i,i=2,3\cdots n& (7)\end{array}$

Which is

$\begin{array}{ll}X(t)={\displaystyle \sum _{i=1}^n{c}_i+r_n}& (8)\end{array}$

c_i is the i-th IMF component, and r_n is the residual.

Due to the phenomenon of frequency mixing in the EMD method, the EEMD method was developed (Wu and Huang, 2009). The EEMD method introduces white noise with a certain signal-to-noise ratio into the original time series for EMD decomposition, and the IMFs obtained by each decomposition are then aggregated. The EEMD method not only greatly improves the modal aliasing defect of EMD, but also avoids the instantaneous noise that the original data may carry.

Repeat equations 2-8, with different Gaussian white noise series assed to X(t), Finally, the original signal is decomposed into a series of IMF_i(c_i) components with frequencies from high to low and a residual r_n.

The EEMD trends in WUE as a specific time t is defined as the value increase in r_n since the start time, that is trend(t)=r_n(t)-r_n(1981). The trends as their changing rates can be calculated (Pan et al., 2018):

$\begin{array}{ll}Rat{e}_{trend}(t)=Trend(t)-Trend(t-1)& (9)\end{array}$

The number of Gaussian white noises was set to 100 times, and the amplitude of these noises was set to 0.2 standard deviations of the raw data after considering the take-off between the decomposition robustness and the required computing time.

To test whether the trend is significant, the EEMD decomposition of Gaussian white noise is verified based on the Monte Carlo method (Pan et al., 2018). The non-significant trend of WUE is assumed to show no trend over time. The significantly non-linear trend of the WUE is divided into the following four categories: increasing trends (Figure 1A), decreasing trends (Figure 1B), negative reversals (Figure 1C), and positive reversals (Figure 1D).

FIGURE 1

Figure 1 The nonlinear trends of WUE. (A) increasing trends; (B) decreasing trends; (C) negative reversals; (D) positive reversals.

2.3 Stability method

Based on nonlinear trends detected by EEMD, excluding the trends with insignificant changes, we computed the interannual anomalies of WUE by removing their nonlinear trends. The standard deviation of interannual anomalies of WUE was identified as the stability of WUE (Wang et al., 2022). A higher standard deviation suggested lower stability. Additionally, the stability of WUE was classified into 6 classes (nonsignificant stable, stable, relatively stable, generally stable, relatively unstable, and unstable) based on natural breaks. The natural interval method is based on natural groupings inherited from the data. When creating classification intervals, similar values are grouped most appropriately and differences between classes are maximized. Elements are divided into classes, for which their boundaries are set at locations where the differences in data values are relatively large.

2.4 Attribution analysis method

Multiple regression analysis was used to extract the impacts of the main drivers on the stability of WUE (Sun et al., 2015; Liu et al., 2018; Chu et al., 2019). We select 6 driving factors (CO₂, temperature (TEM), precipitation (PRE), soil moisture (SOIL), radiation (RAD), and vapor pressure deficit (VPD)) as the explanatory variables (Cai et al., 2015; Sun et al., 2015). Based on min-max normalization, all the driving factors were standardized in advance. The multiple regression analysis method is as follows:

$\begin{array}{ll}\begin{array}{c}SWU{E}_{pre}=aSWU{E}^{CO2}+bSWU{E}^{TEM}+cSWU{E}^{PRE}\\ +dSWU{E}^{RAD}+eSWU{E}^{SOIL}+fSWU{E}^{VPD}+g\end{array}& (10)\end{array}$

Where SWUE_pre indicates the predicted stability of WUE, SWUE^CO2, SWUE^TEM, SWUE^PRE, SWUE^RAD, SWUE^SOIL, SWUE^VPD represent SWUE variations that are driven by CO₂, TEM, PRE, RAD, SOIL, VPD, which were also standard deviations of the interannual anomalies without long-term nonlinear trends (Wang et al., 2022). a–f are the regression coefficients. The absolute value of the regression coefficient can represent the relative importance of the driving factors, and g is the regression constant. In this study, the largest regression coefficient of the multiple regression is the main driving factor on the stability of SWUE.

3 Result

3.1 The WUE trend in China and eight ecological restoration areas

Based on the linear method (Figure 2), WUE had an increasing trend from 1982 to 2015 with an average rate of 0.0135gC/kgH₂O/yr. WUE also had a nonlinear trend detected by EEMD and the increased rate was beyond the linear rate after 2005. In the National Forest Protection Project (NFPP) area (Figure 3), the WUE had increasing trends, the growth rate of WUE was the fastest in the Yellow River (at a rate of 0.0127 gC/kgH₂O/yr) and Taihang Mountains (at a rate of 0.0108 gC/kgH₂O/yr) shelterbelt program areas. It was the slowest in the Three-North shelterbelt program areas (at a rate of 0.0023 gC/kgH₂O/yr) and the Yangtze River shelterbelt program areas (at a rate of 0.0030 gC/kgH₂O/yr). There is a tendency for WUE to increase, with the rate of increase peaking around 2000-2005 and decreasing thereafter. It indicated that although ecological restoration projects has significantly improved water use efficiency, there is a risk of WUE reduction in the future.

FIGURE 2

Figure 2 The linear and nonlinear trend of WUE in the whole of China based on linear regression method and EEMD method.

FIGURE 3

Figure 3 The linear and nonlinear trend of WUE in eight NFPP areas based on linear regression method and EEMD method (A: Three-north areas; B: Yellow River areas; C: Yangtze River shelterbelt program areas; D: Pearl River areas; E: Liaohe River areas; F: Taihang Mountains areas; G: Huaihe and Taihu areas; H: coastal areas).

In China (Figure 4), 33.59% of WUE was nonsignificant, which is assumed to be no trend over time and was mainly located in the Northwest, Northeast, and Tibetan Plateau. 34.19% and 19.72% of WUE had increasing trends and positive reversals, which are mainly distributed in the North China Plain and the Pearl River Basin (Table 1). The decreasing trends and negative reversals of WUE only accounted for 1.86% and 10.64%, respectively, located in the Yangtze River basin. In eight NFPP areas, WUE was dominated by increasing trends in the Pearl River, Taihang Mountains, Huaihe, and Taihu shelterbelt program areas, accounting for 57.90%, 69.42%, and 72.25%, respectively. It had positive reversals in the Yellow River shelterbelt program areas, accounting for 44.63%. However, in the Yangtze River, the WUE was dominated by negative reversals, which accounted for 51.35%.

FIGURE 4

Figure 4 The spatial distribution of the nonlinear trend of WUE in China (Monotonical IN: increasing trends; IN to DE: decreasing trends; IN to DE: negative reversals; DE to IN: positive reversals; A: Three-north areas; B: Yellow River areas; C: Yangtze River shelterbelt program areas; D: Pearl River areas; E: Liaohe River areas; F: Taihang Mountains areas; G: Huaihe and Taihu areas; H: coastal areas).

TABLE 1

Table 1 The percentage of nonlinear trends of WUE in China and eight NFPP areas.

3.2 The spatial distribution of the WUE stability

EEMD decomposes WUE into four interannual variations and a residual (Figure 5), with periods of 2.7, 6.5, 27, and 38-year time scales. The SWUE (standard deviation of interannual variation) ranged from 0.30 to 0.67. The SWUE were relatively stable (ranging from 0.30-0.40) with small interannual fluctuations before 2000. There was a large instability of WUE after 2000. To be more specific, it increased sharply in 2000, 2001, and 2007, with the greatest volatility in 2007.

FIGURE 5

Figure 5 EEMD decomposition for averaged WUE in China (interannual anomalies (IMF1-IMF4); residual trend; stability).

The six classes of SWUE, namely non-significant stable, stable, relatively stable, generally stable, relatively unstable, and unstable in China, account for 33.59%, 15.09%, 12.88%, 19.39%, 12.14%, and 6.91%, respectively (Figure 6 and Table 2). Overall, WUE was stable in the northwest and Tibetan Plateau, where the trend was mainly nonsignificant; while it was unstable in the northeast, southwest, and Yangtze River basin, where the WUE trends were dominated by negative reversals. In NFPP areas, the WUE was stable in Three-North shelterbelt program areas and the stable and relatively stable WUE accounted for 17.57% and 12.03%, respectively. the relatively stable WUE accounted for 32.81% in Coastal shelterbelt program areas. However, in the Liaohe and Taihang Mountains shelterbelt programs areas, the percentage of the relatively unstable and unstable WUE accounted for 57.00% and 45.62%, respectively. In Taihang Mountains shelterbelt program areas, the WUE was dominated by monotonically increasing trends and positive reversals, while the WUE usually experienced negative reversals in Liaohe River areas. It is indicated that WUE with nonlinear trends (positive reversals and negative reversals) are generally unstable, which may be influenced by climatic or anthropogenic factors, leading to positive shifts or negative feedback.

FIGURE 6

Figure 6 The spatial distribution of stability in WUE in China (non sig: non-significant stability of WUE. A: Three-north areas; B: Yellow River areas; C: Yangtze River shelterbelt program areas; D: Pearl River areas; E: Liaohe River areas; F: Taihang Mountains areas; G: Huaihe and Taihu areas; H: coastal areas).

TABLE 2

Table 2 The percentage of stability of WUE in China and eight NFPP areas.

3.3 The dominant climatic factors on the SWUE

To investigate the mechanism of climate fluctuations on SWUE in different regions of China, we analyzed the spatial distribution of the main climatic factors on the SWUE in China (Figure 7). The results show that the spatial variation of SWUE was mainly influenced by VPD in southwestern and northeastern China and along the Yangtze River, where the WUE was unstable. Temperature and solar radiation were the dominant factors for interannual fluctuations of WUE in the North China Plain and central Yangtze River, where the WUE was relatively unstable. In the NFPP areas, VPD was the main driving factor for the interannual variables of WUE in the Yangtze River, Pearl River, and Liaohe River shelterbelt program areas, accounting for 38.25%, 49.38%, and 65.85%, respectively (Table 3). This indicates that the WUE was affected by atmospheric drought in these ecological project areas, which led to low SWUE. Precipitation and soil moisture were the main driving factors for the interannual variables of WUE in the Yellow River and Tibetan Plateau, where the WUE was stable. Temperature and radiation promoted interannual variables of WUE in the Taihang Mountains and Yangtze River shelterbelt program areas, where the WUE was relatively unstable. The overall contribution of CO₂ to SWUE was low in eight NFPP areas. It indicated that VPD, temperature, and radiation lead to unstable changes in WUE, while precipitation and soil moisture lead to stable changes. Therefore, differences in the dominant climate factors in different regions need to be considered when exploring the influence of climate on the SWUE in China.

FIGURE 7

Figure 7 The dominant driving factors of WUE in China (A: Three-north areas; B: Yellow River areas; C: Yangtze River shelterbelt program areas; D: Pearl River areas; E: Liaohe River areas; F: Taihang Mountains areas; G: Huaihe and Taihu areas; H: coastal areas).

TABLE 3

Table 3 The percentage of driving factors of WUE in China and eight NFPP areas.

4 Discussion

4.1 The nonlinear trend of WUE

The core meaning of WUE can be summarized as the ratio of productivity to water consumption (El Masri et al., 2019; Jiao et al., 2022). A deep understanding of the WUE trend is an important entry point for coupling the water-carbon cycle, energy conversion, resource use, and climate change issues, especially in the ecological restoration areas (Jiao et al., 2022). Basing on linear method, previous studies revealed that the WUE shows an increasing trend with a rate of 0.0025 gC kg/H₂O/yr globally (Xue et al., 2015; Yang L. et al., 2022; Zhao et al., 2022). In this study, the WUE was also dominated by increasing trends (34.19%) in China, while decreasing trends only accounted for 1.86%. Additionally, the nonlinear trends of WUE were detected by the EEMD method. 19.72% of WUE had positive reversals while 10.64% exhibited negative reversals, which is rarely explored by the linear methods (Jiao et al., 2022). Ignoring non-linear changes in WUE may lead to an overestimation of ecosystem productivity and an underestimation of water deficit (Hatfield and Dold, 2019). Thus, revealing the nonlinear trends of WUE provided a deep understanding of ecosystem functioning (Jiao et al., 2022). In eight NFPP areas, WUE was dominated by increasing trends or positive reversals. The time of turning points mainly occurred in 2000-2005, aligning well with the implementation of the NFPP (Chen et al., 2019). The restoration project has increased the area of forest and scrub at high WUE levels, and significantly reduced the area of farmland and grassland at low WUE levels, thus increased the overall WUE. This shows that the implementation of ecological restoration projects will improve the sustainability of WUE to a certain extent (Ding et al., 2021). However, in the Yangtze River shelterbelt program areas, the negative reversals of WUE accounted for 51.35%, indicating that the Yangtze River basin is facing reduced productively and drought risks, which limits crop production and triggers grassland fires secondary disasters, such as grassland fires and crop pests and diseases (Venkatappa et al., 2021).

4.2 The stability of WUE

Although the spatial distribution of the WUE trend has been determined using linear or nonlinear trends in previous studies, the stability of WUE remained unclear (Rödenbeck et al., 2018; Jiao et al., 2022). Jiao et al. (2022) obtained the WUE stability through the contribution of inter-annual variables to the secular trend, and revealed that the WUE was strongly stable in the north and Loess Plateau while interannual variables were found in the southwest. The trend itself can reflect the stability of the WUE. In this study, based on the EEMD method, the standard deviation of residual trends of WUE is used as the indicator of stability, and the spatial variations of SWUE were compared and analyzed. The change in WUE experienced large inter-annual fluctuations after 2000, especially from 2000 to 2005. Overall, compared to northwestern China and the Tibetan Plateau, the WUE in eastern China was much more unstable. A previous study showed that most of the ecological functioning has been improved in China (Jin et al., 2018). In this study, the unstable WUE experienced positive reversals or negative reversals. The traditional linear methods that ignore non-linear trends exaggerate the stability of WUE (Hao et al., 2021). Thus, there are hidden risks beneath apparent ecosystem degradation, and ecosystem improvement or degradation may be limited and exaggerated (Pan et al., 2018).

Previous studies have revealed that ecological restoration promotes the improvement of ecosystem functioning (Yang et al., 2017; Song et al., 2022). In this study, the WUE was stable with increasing trends in the Pearl River, Huaihe and Taihu, coastal shelterbelt program areas. Although WUE was unstable in the Yellow River and Taihang Mountains shelterbelt program areas, it experienced positive shifts. Previous studies have shown that the ecological restoration program contributed to a significant increase in vegetation productivity in the Yellow River and Taihang Mountains shelterbelt program areas (Xue et al., 2022). It is indicated that most ecological restoration projects in China enhance carbon sequestration thus leading to increasing trends or positive shifts in WUE (Xue et al., 2022). However, some ecological restoration may inhibit or reverse the trend in WUE (Lu et al., 2018). Our result demonstrated that WUE was unstable and may be easily altered from increasing to decreasing in the Yangtze River shelterbelt program areas. In the Yangtze River shelterbelt program area, the massive planting of trees stimulates an increase in the water demand of the forest, which in turn absorbs large amounts of soil water and stimulates increased evapotranspiration (Zhang et al., 2016). The simultaneous increase in productivity and evapotranspiration can lead to fluctuating changes in WUE, especially during the implementation phase of the ecological project from 2000-2005 (Tong et al., 2019). In addition, the Yangtze River is in the East Asian monsoon climate zone and is significantly influenced by the monsoon, with frequent drought disasters (Hong et al., 2014). Studies have shown that drought events in the Yangtze River have begun to increase and intensify in recent years, which deeply affected vegetation productivity (Horion et al., 2016). In areas with significant human influence and stable WUE, such as the southeastern coast, there is little room for further improvement of WUE stability. Therefore, achieving the sustainability of WUE may depend on reducing anthropogenic carbon emissions in these areas (Wang et al., 2022). Therefore, the spatial heterogeneity in WUE stability underscores the importance of implementing management strategies according to the local condition.

4.3 Driving mechanism of WUE

Climate change has a profound impact on the functioning of ecosystems and, as a result, on the coupled cycles of carbon and water (Liu et al., 2020). Some climatic factors influence the trend and stability of WUE, such as CO₂ concentration, VPD, temperature, precipitation, soil moisture, etc. (Cheng et al., 2017; Hatfield and Dold, 2019). In our study, the stability variation of WUE is mainly caused by precipitation and soil moisture, while instability is mainly induced by saturation water vapor pressure difference, temperature, and solar radiation in China.

Previous studies suggested that the atmospheric water demand, represented by VPD overrode other climatic factors exerting dominantly negative effects on WUE change in alpine meadow ecosystems (Cheng et al., 2017). In this study, the VPD played a dominant role in controlling the instability of the WUE trend in the Yangtze River and Liaohe River shelterbelt program areas, where the WUE experienced negative reversals. By controlling plant stomatal activity, atmospheric drought can affect carbon acquisition and water transpiration losses (Konings et al., 2017). In recent years, China has undergone an increase in atmospheric vapor pressure deficit (Lopez et al., 2021). VPD plays a dominant role in controlling ET in these areas, with higher VPD leading to a dramatic increase in ET (Liu et al., 2020). Additionally, high VPD should lead to partial stomatal closure and suppression of photosynthetic rates (Ding et al., 2018). All these effects could lead to a negative response of WUE to changes in VPD (Cao et al., 2020). This negative impact is becoming stronger in terms of the severity and extent of the effects, indicating that atmospheric drought is becoming increasingly harmful to productivity (Ding et al., 2018). Therefore, in areas where VPD has led to a decline in WUE, the impact of atmospheric drought on vegetation should be closely monitored and the negative impact of atmospheric drought on vegetation should be reduced through artificial measures, such as artificial rainfall and irrigation.

Temperature and radiation are the dominant controlling factors on the instability of WUE in the Taihu and Huaihe and coastal shelterbelt program areas. In the Taihu Lake and Huaihe River shelterbelt program area, the increase of temperature promotes the increase in WUE, while in some coastal areas, the increase of temperature suppresses the increase in WUE and shifts it from decreasing to increasing. The appropriate increase in temperature prolongs the growth period of plants, leading to a higher increase in GPP than ET, resulting in a monotonic increase in WUE in the Huaihe and Taihu shelterbelt program areas. However, in coastal areas, especially in subtropical areas, the temperature increase promotes ET much more than GPP, resulting in negative reversals and large interannual fluctuations in WUE. There is a threshold for the effect of temperature on water utilization and a too high or too low temperature can harm plant WUE (Hatfield and Dold, 2019). When the temperature is below the threshold, WUE increases with increasing temperature, while when the temperature is above the threshold, WUE shows a negative relationship with temperature (Hatfield and Dold, 2019). The main reason is that enzyme activity is significantly affected when the optimum temperature is reached (Huang M. et al., 2015). Conversely, an increase of temperature causes an increase in VPD and thus increases the transpiration rate of vegetation (Xue et al., 2015). Xue et al. (2015) found that globally, WUE tended to increase linearly with temperature in the cooler regions, reaching a maximum at 18.5°C and decreasing thereafter. Therefore, high temperatures can cause instability in WUE with negative reversals. Thus, more attention should be paid to the high temperature, especially heat waves. In addition, with the economically developed eastern coast and the rapid expansion of urbanization, vegetation productivity is dominated by decreasing or increasing to decreasing trends, leading to negative reversals and instability in WUE (Liu et al., 2023).

Precipitation and soil moisture promoted the stability of WUE in the Yellow River shelterbelt program areas. Precipitation and soil moisture in these areas are relatively low and are crucial elements determining ecosystem composition, structure, and function (Zhang et al., 2020). Interannual variation in vegetation productivity in semi-arid regions is closely related to interannual variables of drought and precipitation (Zhang et al., 2020; Yang L. et al., 2022). Moderate rainfall could offset the effect of drought and keep the stability of WUE in the Yellow River program area (Cai et al., 2020; Wang M. et al., 2021). Previous studies revealed that WUE was positively correlated with precipitation and specific humidity (Xue et al., 2015; Wang H. et al., 2021). Liu et al. (2020) suggested that increased soil moisture contributed to a positive trend in WUE in humid and high latitudes of northern China, which could also enhance carbon sequestration because of the increased water availability (Liu et al., 2020). In addition, the Yellow River basin is the most effective area for ecological restoration projects in China, such as Natural Forest protection, afforestation, economic compensation, etc., which enhance carbon accumulation and greatly contribute to the increasing trends and stability in WUE (Kou et al., 2021; Zhang et al., 2022).

5 Conclusions

The stable WUE was dominated by nonsignificant trends and increasing trends, accounting for 33.59% and 34.19%, respectively. The nonsignificant trend of stable WUE was mainly located in Three-North shelterbelt program areas, and the increasing trend of stable WUE was in Huaihe and Taihu, Taihang Mountains, and Pearl River shelterbelt program areas. Precipitation and soil moisture promoted stable WUE in these project areas. The unstable WUE was dominated by positive reversals or negative reversals of WUE trends. The positive reversals of unstable WUE were mainly located in the Yellow River shelterbelt program areas, which was promoted by temperature and radiation, while the negative reversals of unstable WUE were mainly distributed in the Yangtze River and Liaohe shelterbelt program areas, which was mainly induced by VPD.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

XX: Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. JL: Investigation, Methodology, Software, Supervision, Resources, Visualization, Writing – original draft. FJ: Investigation, Methodology, Software, Supervision, Resources, Visualization, Writing – original draft. KZ: Methodology, Supervision, Conceptualization, Investigation, Software, Writing – review & editing. YY: Data curation, Formal Analysis, Methodology, Project administration, Supervision, Validation, Writing – original draft. JQ: Writing – original draft. YZ: Data curation, Methodology, Supervision, Writing – original draft, Writing – review & editing. NL: Funding acquisition, Resources, Visualization, Writing – review & editing. CZ: Conceptualization, Funding acquisition, Resources, Visualization, Writing – review & editing.

Funding

This study was supported by the National Key R&D Program of China (2021YFB3901104). The Special Fund of the Jiangsu for Carbon Peak and Carbon Neutralization Science and Technology Innovation (BK20220021), and the Special Fund of the Chinese Central Government for Basic Scientific Research Operations in the commonweal Research Institute (GYZX210405).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Beer, C., Ciais, P., Reichstein, M., Baldocchi, D., Law, B. E., Papale, D., et al. (2009). Temporal and among-site variability of inherent water use efficiency at the ecosystem level. Global Biogeochem. Cycles 23, 1–13. doi: 10.1029/2008gb003233

CrossRef Full Text | Google Scholar

Belmecheri, S., Maxwell, R. S., Taylor, A. H., Davis, K. J., Guerrieri, R., Moore, D. J. P., et al. (2021). Precipitation alters the CO(2) effect on water-use efficiency of temperate forests. Glob. Chang. Biol. 27, 1560–1571. doi: 10.1111/gcb.15491

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, D., Ge, Q., Wang, X., Liu, B., Goudie, A. S., Hu., S. (2020). Contributions of ecological programs to vegetation restoration in arid and semiarid China. Environ. Res. Lett. 15, 1–11. doi: 10.1088/1748-9326/abbde9

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, H., Yang, X., Xu., X. (2015). Human-induced grassland degradation/restoration in the central Tibetan Plateau: The effects of ecological protection and restoration projects. Ecol. Eng. 83, 112–119. doi: 10.1016/j.ecoleng.2015.06.031

CrossRef Full Text | Google Scholar

Cao, R., Hu, Z., Jiang, Z., Yang, Y., Zhao, W., Wu, G., et al. (2020). Shifts in ecosystem water use efficiency on China's loess plateau caused by the interaction of climatic and biotic factors over 1985–2015. Agric. For. Meteorol. 291 (1), 1–9. doi: 10.1016/j.agrformet.2020.108100

CrossRef Full Text | Google Scholar

Chen, J. M., Ju, W., Ciais, P., Viovy, N., Liu, R., Liu, Y., et al. (2019). Vegetation structural change since 1981 significantly enhanced the terrestrial carbon sink. Nat. Commun. 10, 4259. doi: 10.1038/s41467-019-12257-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, L., Zhang, L., Wang, Y. P., Canadell, J. G., Chiew, F. H. S., Beringer, J., et al. (2017). Recent increases in terrestrial carbon uptake at little cost to the water cycle. Nat. Commun. 8, 110. doi: 10.1038/s41467-017-00114-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, H., Venevsky, S., Wu, C., Wang., M. (2019). NDVI-based vegetation dynamics and its response to climate changes at Amur-Heilongjiang River Basin from 1982 to 2015. Sci. Total Environ. 650, 2051–2062. doi: 10.1016/j.scitotenv.2018.09.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Z., Liu, Y., Wang, L., Chen, Y., Yu, P., Ma, M., et al. (2021). Effects and implications of ecological restoration projects on ecosystem water use efficiency in the karst region of Southwest China. Ecol. Eng. 170, 1–9. doi: 10.1016/j.ecoleng.2021.106356

CrossRef Full Text | Google Scholar

Ding, J., Yang, T., Zhao, Y., Liu, D., Wang, X., Yao, Y., et al. (2018). Increasingly important role of atmospheric aridity on Tibetan Alpine grasslands. Geophys. Res. Lett. 45, 2852–2859. doi: 10.1002/2017gl076803

CrossRef Full Text | Google Scholar

El Masri, B., Schwalm, C., Huntzinger, D. N., Mao, J., Shi, X., Peng, C., et al. (2019). Carbon and water use efficiencies: A comparative analysis of ten terrestrial ecosystem models under changing climate. Sci. Rep. 9, 14680. doi: 10.1038/s41598-019-50808-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonsamo, A., Ciais, P., Miralles, D. G., Sitch, S., Dorigo, W., Lombardozzi, D., et al. (2021). Greening drylands despite warming consistent with carbon dioxide fertilization effect. Glob. Chang. Biol. 27, 3336–3349. doi: 10.1111/gcb.15658

PubMed Abstract | CrossRef Full Text | Google Scholar

Hai, X., Li, J., Li, J., Liu, Y., Dong, L., Wang, X., et al. (2022). Variations in plant water use efficiency response to manipulated precipitation in a temperate grassland. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.881282

CrossRef Full Text | Google Scholar

Hao, X., Zhang, J., Fan, X., Hao, H., Li., Y. (2021). Quantifying soil moisture impacts on water use efficiency in terrestrial ecosystems of China. Remote Sens. 13 (21), 4257. doi: 10.3390/rs13214257

CrossRef Full Text | Google Scholar

Hatfield, J. L., Dold, C. (2019). Water-use efficiency: advances and challenges in a changing climate. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00103

CrossRef Full Text | Google Scholar

Hong, X., Guo, S., Zhou, Y., Xiong., L. (2014). Uncertainties in assessing hydrological drought using streamflow drought index for the upper Yangtze River basin. Stoch. Env. Res. Risk A 29, 1235–1247. doi: 10.1007/s00477-014-0949-5

CrossRef Full Text | Google Scholar

Horion, S., Prishchepov, A. V., Verbesselt, J., de Beurs, K., Tagesson, T., Fensholt., R. (2016). Revealing turning points in ecosystem functioning over the Northern Eurasian agricultural frontier. Glob. Chang. Biol. 22, 2801–2817. doi: 10.1111/gcb.13267

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Piao, S., Sun, Y., Ciais, P., Cheng, L., Mao, J., et al. (2015). Change in terrestrial ecosystem water-use efficiency over the last three decades. Glob. Chang. Biol. 21, 2366–2378. doi: 10.1111/gcb.12873

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, L., Wang, B., Niu, X., Gao, P., Song., Q. (2019). Changes in ecosystem services and an analysis of driving factors for China's Natural Forest Conservation Program. Ecol. Evol. 9, 3700–3716. doi: 10.1002/ece3.4925

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Yu, H., Guan, X., Wang, G., Guo., R. (2015). Accelerated dryland expansion under climate change. Nat. Climate Change 6, 166–171. doi: 10.1038/nclimate2837

CrossRef Full Text | Google Scholar

Jiang, Y., Still, C. J., Rastogi, B., Page, G. F. M., Wharton, S., Meinzer, F. C., et al. (2019). Trends and controls on water-use efficiency of an old-growth coniferous forest in the Pacific Northwest. Environ. Res. Lett. 14 (7), 1–12. doi: 10.1088/1748-9326/ab2612

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, F., Xu, X., Zhang, M., Gong, H., Liu, H., Wang., K. (2022). Contributory factors of the secular trends to changes in ecosystem water-use efficiency in China. J. Hydro.l 615, 128690. doi: 10.1016/j.jhydrol.2022.128690

CrossRef Full Text | Google Scholar

Jin, N., Ren, W., Tao, B., He, L., Ren, Q., Li, S., et al. (2018). Effects of water stress on water use efficiency of irrigated and rainfed wheat in the Loess Plateau, China. Sci. Total Environ. 642, 1–11. doi: 10.1016/j.scitotenv.2018.06.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Konings, A. G., Williams, A. P., Gentine., P. (2017). Sensitivity of grassland productivity to aridity controlled by stomatal and xylem regulation. Nat. Geosci. 10, 284–288. doi: 10.1038/ngeo2903

CrossRef Full Text | Google Scholar

Kou, P., Xu, Q., Jin, Z., Yunus, A. P., Luo, X., Liu., M. (2021). Complex anthropogenic interaction on vegetation greening in the Chinese Loess Plateau. Sci. Total Environ. 778, 146065. doi: 10.1016/j.scitotenv.2021.146065

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Bowman, K. W., Schimel, D. S., Parazoo, N. C., Jiang, Z., Lee, M., et al. (2017). Contrasting carbon cycle responses of the tropical continents to the 2015-2016 El Nino. Science 358 (6360), eaam5690. doi: 10.1126/science.aam5690

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Feng, X., Fu., B. (2020). Changes in global terrestrial ecosystem water use efficiency are closely related to soil moisture. Sci. Total Environ. 698, 134165. doi: 10.1016/j.scitotenv.2019.134165

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Wang, Z., Wang, Z., Zeng, Y., Xue, P., Zhang., M. (2023). Stability of the ecosystem gross primary productivity increasing in Chinese forestry ecological engineering area. Agricult. Ecosyst. Environ. 356, 108636. doi: 10.1016/j.agee.2023.108636

CrossRef Full Text | Google Scholar

Liu, R., Xiao, L. L., Liu, Z., Dai., J. C. (2018). Quantifying the relative impacts of climate and human activities on vegetation changes at the regional scale. Ecol. Indic. 93, 91–99. doi: 10.1016/j.ecolind.2018.04.047

CrossRef Full Text | Google Scholar

Lopez, J., Way, D. A., Sadok., W. (2021). Systemic effects of rising atmospheric vapor pressure deficit on plant physiology and productivity. Glob. Chang. Biol. 27, 1704–1720. doi: 10.1111/gcb.15548

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, F., Hu, H., Sun, W., Zhu, J., Liu, G., Zhou, W., et al. (2018). Effects of national ecological restoration projects on carbon sequestration in China from 2001 to 2010. Proc. Natl. Acad. Sci. U.S.A. 115, 4039–4044. doi: 10.1073/pnas.1700294115

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Jia, X., Zha, T., Bourque, C. P. A., Tian, Y., Bai, Y., et al. (2019). Ecosystem water use efficiency in a young plantation in Northern China and its relationship to drought. Agric. For. Meteorol. 275, 1–10. doi: 10.1016/j.agrformet.2019.05.004

CrossRef Full Text | Google Scholar

Ma, N., Szilagyi, J., Zhang, Y., Liu., W. (2019). Complementary-relationship-based modeling of terrestrial evapotranspiration across China during 1982–2012: validations and spatiotemporal analyses. J. Geophys. Res.: Atmos. 124, 4326–4351. doi: 10.1029/2018jd029850

CrossRef Full Text | Google Scholar

Migliavacca, M., Musavi, T., Mahecha, M. D., Nelson, J. A., Knauer, J., Baldocchi, D. D., et al. (2021). The three major axes of terrestrial ecosystem function. Nature 598, 468–472. doi: 10.1038/s41586-021-03939-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, N., Feng, X., Fu, B., Wang, S., Ji, F., Pan., S. (2018). Increasing global vegetation browning hidden in overall vegetation greening: Insights from time-varying trends. Remote Sens. Environ. 214, 59–72. doi: 10.1016/j.rse.2018.05.018

CrossRef Full Text | Google Scholar

Piao, S., Wang, X., Wang, K., Li, X., Bastos, A., Canadell, J. G., et al. (2020). Interannual variation of terrestrial carbon cycle: Issues and perspectives. Glob. Chang. Biol. 26, 300–318. doi: 10.1111/gcb.14884

PubMed Abstract | CrossRef Full Text | Google Scholar

Rödenbeck, C., Zaehle, S., Keeling, R., Heimann., M. (2018). How does the terrestrial carbon exchange respond to inter-annual climatic variations? A quantification based on atmospheric CO&lt;sub<2&lt;/sub< data. Biogeosciences 15, 2481–2498. doi: 10.5194/bg-15-2481-2018

CrossRef Full Text | Google Scholar

Song, W., Feng, Y., Wang., Z. (2022). Ecological restoration programs dominate vegetation greening in China. Sci. Total Environ. 848, 157729. doi: 10.1016/j.scitotenv.2022.157729

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y. L., Yang, Y. L., Zhang, L., Wang., Z. L. (2015). The relative roles of climate variations and human activities in vegetation change in North China. Phys. Chem. Earth 87-88, 67–78. doi: 10.1016/j.pce.2015.09.017

CrossRef Full Text | Google Scholar

Tong, X., Brandt, M., Yue, Y., Horion, S., Wang, K., Keersmaecker, W. D., et al. (2018). Increased vegetation growth and carbon stock in China karst via ecological engineering. Nat. Sustain. 1, 44–50. doi: 10.1038/s41893-017-0004-x

CrossRef Full Text | Google Scholar

Tong, X., Mu, Y., Zhang, J., Meng, P., Li., J. (2019). Water stress controls on carbon flux and water use efficiency in a warm-temperate mixed plantation. J. Hydrol. 571, 669–678. doi: 10.1016/j.jhydrol.2019.02.014

CrossRef Full Text | Google Scholar

Tong, X., Wang, K., Yue, Y., Brandt, M., Liu, B., Zhang, C., et al. (2017). Quantifying the effectiveness of ecological restoration projects on long-term vegetation dynamics in the karst regions of Southwest China. Int. J. Appl. Earth Observ. Geoinform. 54, 105–113. doi: 10.1016/j.jag.2016.09.013

CrossRef Full Text | Google Scholar

Venkatappa, M., Sasaki, N., Han, P., Abe., I. (2021). Impacts of droughts and floods on croplands and crop production in Southeast Asia - An application of Google Earth Engine. Sci. Total Environ. 795, 148829. doi: 10.1016/j.scitotenv.2021.148829

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Ding, Z., Wu, C., Song, L., Ma, M., Yu, P., et al. (2021). Divergent responses of ecosystem water-use efficiency to extreme seasonal droughts in Southwest China. Sci. Total Environ. 760, 143427. doi: 10.1016/j.scitotenv.2020.143427

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Li, M., Wang, J., Li, X., Wang., L. (2020). An analytical reductionist framework to separate the effects of climate change and human activities on variation in water use efficiency. Sci. Total Environ. 727, 138306. doi: 10.1016/j.scitotenv.2020.138306

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Li, X., Xiao, J., Ma., M. (2021). Evapotranspiration components and water use efficiency from desert to alpine ecosystems in drylands. Agric. For. Meteorol., 298–299, 108283. doi: 10.1016/j.agrformet.2020.108283

CrossRef Full Text | Google Scholar

Wang, K., Piao, S., He, Y., Liu, Y., He., H. (2022). Spatial variations and mechanisms for the stability of terrestrial carbon sink in China. Sci. China Earth Sci. 66, 227–236. doi: 10.1007/s11430-021-1003-5

CrossRef Full Text | Google Scholar

Wu, Z., Huang, N. E. (2009). Ensemble emprirical mode decomposition: A noise-assisted data analysis method. Adv. Adaptive Data Anal. 01, 1–40. doi: 10.1142/S1793536909000047

CrossRef Full Text | Google Scholar

Wu, Z., Huang, N. E., Long, S. R., Peng., C. K. (2007). On the trend, detrending, and variability of nonlinear and nonstationary time series. Proc. Natl. Acad. Sci. U.S.A. 104, 14889–14894. doi: 10.1073/pnas.0701020104

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Jiao, F., Liu, H., Gong, H., Zou, C., Lin, N., et al. (2022a). Persistence of increasing vegetation gross primary production under the interactions of climate change and land use changes in Northwest China. Sci. Total Environ. 834, 155086. doi: 10.1016/j.scitotenv.2022.155086

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Liu, J., Jiao, F., Zhang, K., Ye, X., Gong, H., et al. (2022b). Ecological engineering induced carbon sinks shifting from decreasing to increasing during 1981-2019 in China. Sci. Total Environ. 864, 161037. doi: 10.1016/j.scitotenv.2022.161037

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, B.-L., Guo, Q., Otto, A., Xiao, J., Tao, S., Li., L. (2015). Global patterns, trends, and drivers of water use efficiency from 2000 to 2013. Ecosphere 6 (10), 1–13. doi: 10.1890/es14-00416.1

CrossRef Full Text | Google Scholar

Xue, Y., Liang, H., Zhang, B., He., C. (2022). Vegetation restoration dominated the variation of water use efficiency in China. J. Hydro.l 612, 128257. doi: 10.1016/j.jhydrol.2022.128257

CrossRef Full Text | Google Scholar

Yang, L., Feng, Q., Wen, X., Barzegar, R., Adamowski, J. F., Zhu, M., et al. (2022). Contributions of climate, elevated atmospheric CO2 concentration and land surface changes to variation in water use efficiency in Northwest China. Catena 213, 106220. doi: 10.1016/j.catena.2022.106220

CrossRef Full Text | Google Scholar

Yang, W., Liu, W., Vina, A., Luo, J., He, G., Ouyang, Z., et al. (2013). Performance and prospects of payments for ecosystem services programs: evidence from China. J. Environ. Manage 127, 86–95. doi: 10.1016/j.jenvman.2013.04.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Shi, Y., Sun, W., Chang, J., Zhu, J., Chen, L., et al. (2022). Terrestrial carbon sinks in China and around the world and their contribution to carbon neutrality. Sci. China Life Sci. 65, 861–895. doi: 10.1007/s11427-021-2045-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, H. F., Yao, L., Wang, Y. B., Li., J. L. (2017). Relative contribution of climate change and human activities to vegetation degradation and restoration in North Xinjiang, China. Rangeland J. 39, 289–302. doi: 10.1071/rj16069

CrossRef Full Text | Google Scholar

Yu, L., Zhao, X., Gao, X., Jia, R., Yang, M., Yang, X., et al. (2021). Effect of natural factors and management practices on agricultural water use efficiency under drought: A meta-analysis of global drylands. J. Hydrol., 594, 125977. doi: 10.1016/j.jhydrol.2021.125977

CrossRef Full Text | Google Scholar

Zhang, L., Xiao, J., Zheng, Y., Li, S., Zhou., Y. (2020). Increased carbon uptake and water use efficiency in global semi-arid ecosystems. Environ. Res. Lett. 15 (3), 034022. doi: 10.1088/1748-9326/ab68ec

CrossRef Full Text | Google Scholar

Zhang, L., Xiao, J., Zhou, Y., Zheng, Y., Li, J., Xiao., H. (2016). Drought events and their effects on vegetation productivity in China. Ecosphere 7 (12), e01591. doi: 10.1002/ecs2.1591

CrossRef Full Text | Google Scholar

Zhang, C.-y., Zhao, L., Zhang, H., Chen, M.-n., Fang, R.-y., Yao, Y., et al. (2022). Spatial-temporal characteristics of carbon emissions from land use change in Yellow River Delta region, China. Ecol. Indic. 136, 108623. doi: 10.1016/j.ecolind.2022.108623

CrossRef Full Text | Google Scholar

Zhao, F., Wu, Y., Ma, S., Lei, X., Liao., W. (2022). Increased water use efficiency in China and its drivers during 2000–2016. Ecosystems 25, 1476–1492. doi: 10.1007/s10021-021-00727-4

CrossRef Full Text | Google Scholar