# The Variations of Satellite-Based Ecosystem Water Use and Carbon Use Efficiency and Their Linkages with Climate and Human Drivers in the Songnen Plain, China

Authors: Bo Li; Fang Huang; Shuai Chang; Ning Sun
Journal: Advances in Meteorology
 (2019)

---

## Abstract

Ecosystem water use efficiency (WUE) and carbon use efficiency (CUE), as two of the most important ecological indicators of an ecosystem, represent the carbon assimilation rate of unit water consumption and the capacity of transferring carbon from the atmosphere to potential carbon sinks. Revealing WUE and CUE changes and their impact factors is vital for regional carbon-water interactions and carbon budget assessment. Climate affects carbon and water processes differently. Compared to WUE, the variations in CUE in response to climate factors and human activity remain inadequately understood. In this study, ecosystem-level WUE and CUE variations in the Songnen Plain (SNP), Northeast China, during 2001–2015, were investigated using Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. The relationships between WUE, CUE, main climate factors, and human impacts were explored. The results showed that ecosystem WUE and CUE have fluctuated over time, with regional average values of 1.319 gC·kg−1H2O and 0.516, respectively. Deciduous broad-leaved forests had the highest average WUE but the lowest CUE. The multiyear average CUE of grassland ranked in first place, while the lowest WUE indicated that a lesser capacity of net productivity was generated by the use of limited water supply. WUE and CUE showed a downward trend in most areas of the SNP, indicating that the carbon sequestration capacity of the terrestrial ecosystem became weaker in the past 15 years. Annual precipitation and relative humidity had positive influences on WUE and CUE in more than 60% of the study area. The total annual sunshine duration and annual average temperature negatively affected WUE and CUE in most areas. Human activities had a positive effect on ecosystem WUE changes in the SNP but might inhibit CUE variations. Our findings aid in understanding the biological regulation mechanisms of carbon-water cycle coupling and provide a scientific basis for formulating sustainable regional development strategies and guiding water and land resources management.


---

## Body

## 1. Introduction

Water use efficiency (WUE) is defined as the carbon assimilation rate of unit water consumption [1], which links the photosynthetic production of ecosystems to evapotranspiration and reflects ecosystem carbon and water cycles and their interactions. Carbon use efficiency (CUE) is the ratio of net primary productivity (NPP) to gross primary productivity (GPP), which indicates ecosystem capacity to transfer atmospheric CO2 into biomass and carbon sequestration. CUE changes strongly affect ecosystem carbon budgets and turnover. Because WUE and CUE are related to the processes of ecosystem evapotranspiration and photosynthesis, they are regarded as the important indicators for characterizing the carbon-water coupling of ecosystems. The quantitative analysis of spatial-temporal changes of WUE and CUE and their influencing factors will aid in better understanding the effects of future climate change on the water and carbon processes of ecosystems [2].Although the lack of continuous ground observations may hinder the long-term analysis of the dynamics of vegetation development, satellite remote sensing provides useful information for investigating large-scale and long-term variabilities of ecosystem WUE and CUE. In 1999, Bastiaanssen analyzed crop WUE in India based on estimates of crop yields and evaporation (ET) using remote sensing data for the first time [3]. Mo et al. integrated a soil plant atmospheric transmission crop growth model with remote sensing data to predict crop yield, water consumption, and WUE in Hebei Province, China [4]. Zwart et al. indicated that remote sensing data could be used as input data in a water productivity model and estimated the WUE of wheat from a local to global scale [5]. Teixeira et al. calculated regional evapotranspiration with remote sensing data in the middle and lower reaches of the San Francisco River in Brazil and studied the spatial distribution of crop WUE [6]. Based on crop yields derived from remote sensing images and meteorological data, Li et al. estimated large-scale ET and WUE levels of winter wheat and summer corn in the Haihe River watershed [7]. Hu et al. found that there were significant differences in daily WUE of forest and grassland ecosystems [8]. Many attempts have been made to address the relationship between temporal-spatial variations in WUE and climatic factors using ecosystem models or satellite observations. Gang et al. evaluated the spatiotemporal dynamics of global grassland WUE from 2000 to 2013 and revealed the different responses of each grassland type to climate variations [9]. Zhang et al. applied the ecosystem model CEVSA to estimate the spatiotemporal changes of WUE in high mountainous regions in southwest China from 1954 to 2010 and analyzed its response to climate change [10]. Sun et al. posited that global WUE produced different spatial relationships with climate variables [11]. By means of continuous eddy covariance (EC) measurements and time-series MODIS data, Tang et al. investigated the seasonal variations in WUE for grassland and cropland ecosystems in China’s arid and semiarid regions and contended that the dominant climate factors were linked by temperature and precipitation [12].In many studies, CUE has been considered a constant value regardless of ecosystem type or species [13–15], but that assumption at global scale might be controversial because it ignores the influence of environmental factors [16, 17]. Piao et al. demonstrated that the CUE of different vegetation differed greatly from the south temperate zone to the tropics [18]. Nemani et al. calculated CUE at a global scale using remote sensing data [19]. Zhang et al. found that the NPP/GPP ratio exhibited a pattern that did not depend only on the main climatic characteristics (i.e., temperature and precipitation) and geographical factors (namely, latitude and altitude) [20]. Khalifa et al. revealed that WUE and CUE of the different land cover types in Ethiopia had higher magnitudes than their counterparts in Sudan, and both of them exhibited higher sensitivity to drought and variations in precipitation [21]. He et al. investigated the spatial variations in CUE from different process-based carbon cycle models in comparison with those from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data to analyze the responses of CUE to mean annual precipitation and temperature [22]. To our knowledge, the analysis of the relationships of CUE with other climatic factors, such as relative humidity, solar duration, and the comprehensive impacts of human activities, has not been well explored.The Songnen Plain (SNP) is located in the temperate semihumid to semiarid transition ecological fragile zone in Northeast China, which is highly sensitive to global change. Historically, the SNP fostered both agricultural and nomadic cultures. The combined effects of vulnerable physical basis and excessive irrational human activities have made the region suffer from a high risk of land degradation for the past century. Due to concerns over aggravating desertification, many measures for ecological and environmental protection have been taken by the Chinese government and have achieved preliminary effectiveness. Previous studies have mostly focused on the land cover and use changes and effects of agricultural activities on the environment. Little is known about the spatial-temporal patterns of ecosystem-level WUE and CUE and their response to climate change and human impacts in the SNP.The objectives of this study are to explain (1) how ecosystem WUE and CUE in the SNP varied during 2001–2015 with MODIS satellite observations and (2) how climatic factors and human activities affected variations of ecosystem WUE and CUE in the SNP. The results and conclusions complement our knowledge of coupled carbon-water cycles and carbon budgets of the temperate semiarid and semihumid transitional zone ecosystem and their underlying driving mechanisms.


## 2. Materials and Methods

### 2.1. Study Area

The SNP is located in the central part of Northeast China, with a range of 121°38′ to 128°33′E and 42°49′ to 49°12′N and a total area of 22.35 × 104 km2 (Figure 1). It is an alluvial plain across which the Songhua and Nenjiang rivers flow and is situated in the central Songliao Basin between the Xiaoxing’an Mountains, Changbai Mountains, and the Songliao Watershed. It is also known as one of the world’s three major black soil regions and soda salinized land areas. The altitude ranges from 101 to 1632 meters above sea level but is mostly below 200 meters [23]. The main landform types include those of low mountains, alluvial plains, tablelands, and intervals. The SNP belongs to a temperate continental semihumid and semiarid monsoon climate zone that is characterized by significant winds and four seasons with a hot, rainy summer and a cold, dry winter. The average temperature is from −16°C to −26°C in January and 21.2°C to 23°C in July. The annual precipitation ranges between 350 and 800 mm and gradually decreases from east to west [24]. Soils are fertile, with black soil, meadow soil, and chernozem being widely distributed in this plain. From east to west, the natural vegetation types are forest steppe, meadow steppe, and typical grassland. As a typical agropastoral transitional zone, the SNP is more sensitive to climate change. Administratively, the SNP includes 54 county-level cities in the central and western parts of Jilin Province and the southwest of Heilongjiang Province. It is among the important grain commodity bases in China.Figure 1
Location and land cover of the SNP.


### 2.2. Data and Processing

#### 2.2.1. MODIS Data

MOD17A3 yearly NPP and GPP products with 1-km spatial resolution were the basis of this study. Based on the Biome-BGC (biogeochemical cycles) model, NPP data are produced by the Numerical Terradynamic Simulation Group (NTSG) and the University of Montana (UST) (https://ladsweb.modaps.eosdis.nasa.gov/). This dataset’s accuracy and reliability have been widely verified. There are no overall bias and continuous overestimation or underestimation at global scale compared with the field measured values [25, 26]. The annual actual evapotranspiration (ETa) product (MOD16A3), with spatial resolution of 1 km and estimation based on the Penman–Monteith model, was accessed at http://ntsg.umt.edu/project/mod16. The images were a subset for the SNP in ArcGIS.The MODIS land cover type dataset (MCD12Q1) is provided by the Land Processes Distributed Active Archive Center (LP DAAC) (https://lpdaac.usgs.gov/) [27]. This product includes global land cover data from 2001 to 2016 with a spatial resolution of 500 m and adopts five different classification schemes. In this study, 2009 land cover data were used. According to the IGBP class scheme and the regional characteristics of the SNP, vegetation types were summarized as follows: deciduous coniferous forest, deciduous broad-leaved forest, mixed forest, grassland, cropland, and wetland. The data were then resampled to the 1 km resolution of MODIS products.


#### 2.2.2. Meteorological Data

Monthly meteorological data of 29 weather stations in and near the SNP from 2001 to 2015 were acquired from the National Meteorological Information Center of China (http://data.cma.cn/). Yearly average precipitation (mm), total duration of sunshine (h), average temperature (°C), and average relative humidity (%) were computed. A previous study found that Kriging interpolation had high numerical precision at different resolutions, but the results could not show the details of actual terrain and were unsuitable for the plains [28]. The inverse distance weighting (IDW) spatial interpolation technique showed better adaptability in plain areas, with relatively high accuracy of interpolation. Therefore, the gridded precipitation, sunshine hours, temperature, and humidity data with spatial resolution of 1 km were produced by the IDW interpolation method in ArcGIS in this study.


### 2.3. Calculation of WUE and CUE

Ecosystem WUE refers to the rate of carbon assimilation per unit water loss, which is an important indicator of ecological function and carbon-water cycle coupling [29]. Due to differences in research purposes and data acquisition approaches, the calculation of the “water loss” of the ecosystem differs across studies, which may cause different understandings of ecosystem WUE variability mechanisms. In this study, WUE is calculated by the following equation [30]:(1)WUE=NPPET,where NPP and ET is the net primary productivity and evapotranspiration, respectively.Ecosystem CUE describes the relationship between photosynthesis and respiration, which is an important indicator of the ability of plants to transfer carbon [31]. CUE is a key controlling factor for ecosystem carbon storage and can be used to compare carbon cycles in different ecosystems. The calculation of CUE is performed as follows:(2)CUE=NPPGPP,where GPP is the gross primary productivity, representing the ability to capture energy and carbon through photosynthetic plants and the total amount of C assimilation, and NPP reveals the energy of plants stored after losing C from GPP through autotrophic respiration [32].


### 2.4. Trend Analysis

Spatial trends of WUE and CUE were examined by applying a simple linear regression model with time as the independent variable and WUE or CUE as the dependent variables. The outputs of the trend analysis are the maps of regression slope values, as expressed by the following formula:(3)slope=n×∑i=1ni×Ai−∑i=1ni∑i=1nAin×∑i=1ni2−∑i=1ni2,where slope is the slope of the fitted regression line at each pixel; n represents year range; i is 1 for the first year, 2 for the second year, and so on; and Ai represents the WUE or CUE of the year i. A negative regression coefficient (slope < 0) indicates a decline of WUE or CUE, whereas a positive value (slope > 0) depicts an increase trend. F test was used to determine the significance of change trend.


### 2.5. Correlation Analysis

To analyze the intensity of linear correlation between WUE, CUE, and climate factors from 2001 to 2015, the correlation analysis (CA) method was used. The Spearman correlation coefficient (ρ) of each pixel was calculated by the formula as follows [33]:(4)ρ=1−6∑di2n3−n,where ρ is the correlation coefficient ranging from 0 to 1; n is the number of levels; and d is the rank difference of the two pairs of variables. The rank sum test is used to test its significance.


### 2.6. Multiple Regression Residual Analysis

The multivariate regression residuals analysis method was proposed by Evans and Geerken to discriminate land degradation driven by climate and human drivers [34, 35]. In this study, by the regression analysis of WUE, CUE, and climatic factors at the pixel scale, the values of WUE and CUE were predicted, which could be regarded as the impact of climate on WUE and CUE. Residual values were estimated by the differences between observed WUE, CUE, and their predicted values, and thereby, the effect of climate change and human activities on WUE and CUE could be distinguished. The residual can be expressed as follows:(5)ε=D1−D2,where D1 and D2 are the observed and predicted values of WUE or CUE, respectively. ε indicates the human impacts on WUE or CUE. A positive residual (ε > 0) indicates positive effects of human activities on WUE or CUE, whereas a negative value (ε < 0) depicts the negative roles of anthropogenic causes. The significance of the residual sequence variation trend was also analyzed.


## 2.1. Study Area

The SNP is located in the central part of Northeast China, with a range of 121°38′ to 128°33′E and 42°49′ to 49°12′N and a total area of 22.35 × 104 km2 (Figure 1). It is an alluvial plain across which the Songhua and Nenjiang rivers flow and is situated in the central Songliao Basin between the Xiaoxing’an Mountains, Changbai Mountains, and the Songliao Watershed. It is also known as one of the world’s three major black soil regions and soda salinized land areas. The altitude ranges from 101 to 1632 meters above sea level but is mostly below 200 meters [23]. The main landform types include those of low mountains, alluvial plains, tablelands, and intervals. The SNP belongs to a temperate continental semihumid and semiarid monsoon climate zone that is characterized by significant winds and four seasons with a hot, rainy summer and a cold, dry winter. The average temperature is from −16°C to −26°C in January and 21.2°C to 23°C in July. The annual precipitation ranges between 350 and 800 mm and gradually decreases from east to west [24]. Soils are fertile, with black soil, meadow soil, and chernozem being widely distributed in this plain. From east to west, the natural vegetation types are forest steppe, meadow steppe, and typical grassland. As a typical agropastoral transitional zone, the SNP is more sensitive to climate change. Administratively, the SNP includes 54 county-level cities in the central and western parts of Jilin Province and the southwest of Heilongjiang Province. It is among the important grain commodity bases in China.Figure 1
Location and land cover of the SNP.


## 2.2. Data and Processing

### 2.2.1. MODIS Data

MOD17A3 yearly NPP and GPP products with 1-km spatial resolution were the basis of this study. Based on the Biome-BGC (biogeochemical cycles) model, NPP data are produced by the Numerical Terradynamic Simulation Group (NTSG) and the University of Montana (UST) (https://ladsweb.modaps.eosdis.nasa.gov/). This dataset’s accuracy and reliability have been widely verified. There are no overall bias and continuous overestimation or underestimation at global scale compared with the field measured values [25, 26]. The annual actual evapotranspiration (ETa) product (MOD16A3), with spatial resolution of 1 km and estimation based on the Penman–Monteith model, was accessed at http://ntsg.umt.edu/project/mod16. The images were a subset for the SNP in ArcGIS.The MODIS land cover type dataset (MCD12Q1) is provided by the Land Processes Distributed Active Archive Center (LP DAAC) (https://lpdaac.usgs.gov/) [27]. This product includes global land cover data from 2001 to 2016 with a spatial resolution of 500 m and adopts five different classification schemes. In this study, 2009 land cover data were used. According to the IGBP class scheme and the regional characteristics of the SNP, vegetation types were summarized as follows: deciduous coniferous forest, deciduous broad-leaved forest, mixed forest, grassland, cropland, and wetland. The data were then resampled to the 1 km resolution of MODIS products.


### 2.2.2. Meteorological Data

Monthly meteorological data of 29 weather stations in and near the SNP from 2001 to 2015 were acquired from the National Meteorological Information Center of China (http://data.cma.cn/). Yearly average precipitation (mm), total duration of sunshine (h), average temperature (°C), and average relative humidity (%) were computed. A previous study found that Kriging interpolation had high numerical precision at different resolutions, but the results could not show the details of actual terrain and were unsuitable for the plains [28]. The inverse distance weighting (IDW) spatial interpolation technique showed better adaptability in plain areas, with relatively high accuracy of interpolation. Therefore, the gridded precipitation, sunshine hours, temperature, and humidity data with spatial resolution of 1 km were produced by the IDW interpolation method in ArcGIS in this study.


## 2.2.1. MODIS Data

MOD17A3 yearly NPP and GPP products with 1-km spatial resolution were the basis of this study. Based on the Biome-BGC (biogeochemical cycles) model, NPP data are produced by the Numerical Terradynamic Simulation Group (NTSG) and the University of Montana (UST) (https://ladsweb.modaps.eosdis.nasa.gov/). This dataset’s accuracy and reliability have been widely verified. There are no overall bias and continuous overestimation or underestimation at global scale compared with the field measured values [25, 26]. The annual actual evapotranspiration (ETa) product (MOD16A3), with spatial resolution of 1 km and estimation based on the Penman–Monteith model, was accessed at http://ntsg.umt.edu/project/mod16. The images were a subset for the SNP in ArcGIS.The MODIS land cover type dataset (MCD12Q1) is provided by the Land Processes Distributed Active Archive Center (LP DAAC) (https://lpdaac.usgs.gov/) [27]. This product includes global land cover data from 2001 to 2016 with a spatial resolution of 500 m and adopts five different classification schemes. In this study, 2009 land cover data were used. According to the IGBP class scheme and the regional characteristics of the SNP, vegetation types were summarized as follows: deciduous coniferous forest, deciduous broad-leaved forest, mixed forest, grassland, cropland, and wetland. The data were then resampled to the 1 km resolution of MODIS products.


## 2.2.2. Meteorological Data

Monthly meteorological data of 29 weather stations in and near the SNP from 2001 to 2015 were acquired from the National Meteorological Information Center of China (http://data.cma.cn/). Yearly average precipitation (mm), total duration of sunshine (h), average temperature (°C), and average relative humidity (%) were computed. A previous study found that Kriging interpolation had high numerical precision at different resolutions, but the results could not show the details of actual terrain and were unsuitable for the plains [28]. The inverse distance weighting (IDW) spatial interpolation technique showed better adaptability in plain areas, with relatively high accuracy of interpolation. Therefore, the gridded precipitation, sunshine hours, temperature, and humidity data with spatial resolution of 1 km were produced by the IDW interpolation method in ArcGIS in this study.


## 2.3. Calculation of WUE and CUE

Ecosystem WUE refers to the rate of carbon assimilation per unit water loss, which is an important indicator of ecological function and carbon-water cycle coupling [29]. Due to differences in research purposes and data acquisition approaches, the calculation of the “water loss” of the ecosystem differs across studies, which may cause different understandings of ecosystem WUE variability mechanisms. In this study, WUE is calculated by the following equation [30]:(1)WUE=NPPET,where NPP and ET is the net primary productivity and evapotranspiration, respectively.Ecosystem CUE describes the relationship between photosynthesis and respiration, which is an important indicator of the ability of plants to transfer carbon [31]. CUE is a key controlling factor for ecosystem carbon storage and can be used to compare carbon cycles in different ecosystems. The calculation of CUE is performed as follows:(2)CUE=NPPGPP,where GPP is the gross primary productivity, representing the ability to capture energy and carbon through photosynthetic plants and the total amount of C assimilation, and NPP reveals the energy of plants stored after losing C from GPP through autotrophic respiration [32].


## 2.4. Trend Analysis

Spatial trends of WUE and CUE were examined by applying a simple linear regression model with time as the independent variable and WUE or CUE as the dependent variables. The outputs of the trend analysis are the maps of regression slope values, as expressed by the following formula:(3)slope=n×∑i=1ni×Ai−∑i=1ni∑i=1nAin×∑i=1ni2−∑i=1ni2,where slope is the slope of the fitted regression line at each pixel; n represents year range; i is 1 for the first year, 2 for the second year, and so on; and Ai represents the WUE or CUE of the year i. A negative regression coefficient (slope < 0) indicates a decline of WUE or CUE, whereas a positive value (slope > 0) depicts an increase trend. F test was used to determine the significance of change trend.


## 2.5. Correlation Analysis

To analyze the intensity of linear correlation between WUE, CUE, and climate factors from 2001 to 2015, the correlation analysis (CA) method was used. The Spearman correlation coefficient (ρ) of each pixel was calculated by the formula as follows [33]:(4)ρ=1−6∑di2n3−n,where ρ is the correlation coefficient ranging from 0 to 1; n is the number of levels; and d is the rank difference of the two pairs of variables. The rank sum test is used to test its significance.


## 2.6. Multiple Regression Residual Analysis

The multivariate regression residuals analysis method was proposed by Evans and Geerken to discriminate land degradation driven by climate and human drivers [34, 35]. In this study, by the regression analysis of WUE, CUE, and climatic factors at the pixel scale, the values of WUE and CUE were predicted, which could be regarded as the impact of climate on WUE and CUE. Residual values were estimated by the differences between observed WUE, CUE, and their predicted values, and thereby, the effect of climate change and human activities on WUE and CUE could be distinguished. The residual can be expressed as follows:(5)ε=D1−D2,where D1 and D2 are the observed and predicted values of WUE or CUE, respectively. ε indicates the human impacts on WUE or CUE. A positive residual (ε > 0) indicates positive effects of human activities on WUE or CUE, whereas a negative value (ε < 0) depicts the negative roles of anthropogenic causes. The significance of the residual sequence variation trend was also analyzed.


## 3. Results and Discussion

### 3.1. Temporal Dynamics of WUE and CUE

#### 3.1.1. Interannual Changes

The multiyear average WUE in the SNP fluctuated with time, as shown in Figure2. Average WUE peaked in 2003 (3.881 gC·kg−1H2O) in association with higher NPP. Subsequently, the regional WUE kept decreasing and dropped to the lowest value in 2009 (0.551 gC·kg−1H2O), possibly due to high evapotranspiration in that year. WUE started to increase in 2009 and rose to 1.127 gC·kg−1H2O in 2013. It declined again after 2013, and a lower WUE value of 0.940 gC·kg−1H2O was found in 2015.Figure 2
Annual average WUE and CUE in the SNP for the period of 2001–2015.By comparison, the regional CUE was relatively stable, with an average value of 0.516. The average CUE was 0.495 in 2001 and dropped to 0.467 in 2003. It showed an upward trend for the period 2001–2005. The maximum CUE (0.552) was observed in 2005. It decreased to its lowest value (0.461) in 2008, implying the worst level of CUE. Ecosystem CUE increased from 2010 to 2013 and decreased during 2013–2015.


#### 3.1.2. WUE Variations of Different Ecosystems

During 2001–2015, the average WUE of different ecosystems followed a descending order of deciduous broadleaf forest > mixed forest > grassland > deciduous coniferous forest > cropland > wetland. The highest WUE in the SNP was found in the deciduous broad-leaved forest, with a multiyear average value of 1.328 gC·kg−1H2O. The wetland showed the lowest WUE (0.954 gC·kg−1H2O). The maximum WUE values of deciduous broad-leaved forest and mixed forest were both found in 2013 (2.019 gC·kg−1H2O and 1.669 gC·kg−1H2O, respectively), whereas the minimum values were 0.577 gC·kg−1H2O and 0.685 gC·kg−1H2O in 2001, respectively. The average WUE of mixed forest was 1.158 gC·kg−1H2O. The interannual average WUEs of grassland and cropland were 0.984 gC·kg−1H2O and 0.961 gC·kg−1H2O, respectively. In deciduous coniferous forest, the multiyear average WUE was 0.972 gC·kg−1H2O, with the highest value (1.304 gC·kg−1H2O) in 2014 and the lowest value (0.584 gC·kg−1H2O) in 2009. The interannual WUE mean of wetland was 0.954 gC·kg−1H2O, and the maximum and minimum WUEs were observed in 2005 (1.138 gC·kg−1H2O) and 2009 (0.615 gC·kg−1H2O), separately. Due to the significantly increased evapotranspiration in 2009, the average WUEs of all ecosystems were relatively low.Regarding the temporal variation of WUE in different vegetation types, the decreasing trend followed the order of deciduous coniferous forest (2.43%) > grassland (1.53%) > wetland (1.37%) > cropland (0.42%). In contrast, the increased trends of WUE were for deciduous broadleaf forest (6.98%) and mixed forest (4.53%). The increased WUE in deciduous broadleaf forest and mixed forest indicated that WUE improved and that greater net productivity might be generated by using a limited water supply. The WUE of deciduous coniferous forest, grassland, wetland, and cropland revealed the opposite trend.


#### 3.1.3. CUE Variations of Different Ecosystems

Figure3(b) illustrates the CUE changes of different ecosystems from 2001 to 2015. Grassland showed the highest CUE, with a multiyear average CUE of 0.567. The CUE of deciduous broad-leaved forest was the lowest. CUE of different ecosystems followed a descending order of grassland > wetland > deciduous coniferous forest > cropland > mixed forest > deciduous broad-leaved forest. The multiyear average CUE of deciduous coniferous forest was 0.533. The highest (0.597) and lowest (0.472) CUEs were observed in 2013 and 2003, respectively. Deciduous broad-leaved forest and mixed forest showed a similar change tendency of annual CUE, with average interannual values of 0.479 and 0.480, respectively. The highest CUE (0.560) value occurred in 2005, while their lowest values of 0.239 and 0.308 were found in 2001.Figure 3
WUE (a) and CUE (b) of different ecosystems in the SNP from 2001 to 2015.
(a)
(b)The annual average CUE values of the deciduous broad-leaved forest and the mixed forest were low before 2003, but they suddenly increased in 2004 (Figure3(b)). This might be attributed to the lower annual NPP during 2001–2003, whereas GPP showed a small difference. The average CUE of cropland was 0.510, with the peak of 0.551 in 2003 and the trough of 0.452 in 2008. The multiyear average CUE of wetland was 0.542. It was found that wetland CUE peaked in 2002 (0.801) and that the minimum value occurred in 2007 (0.472). The reason was that NPP of wetland in 2002 reached 579.2 g/m2, an increase of 27.56% from that in 2001.CUE can quantitatively depict the ability of vegetation to fix atmospheric CO2 in the surface ecosystem [34–36]. During the past 15 years, the CUE of deciduous broad-leaved forest and mixed forest showed an increasing trend, with that for deciduous broad-leaved forest rising most obviously. However, average CUE values of deciduous coniferous forest, grassland, cropland, and wetland tended to decrease in the following sequence: deciduous coniferous forest > wetland > grassland > cropland. From 2001 to 2015, deciduous coniferous forest, grassland, cropland, and wetland CUE dropped, indicating their enhanced vegetation respiration. The respiration decomposed organic matter into simple inorganic substances such as water and carbon dioxide; thereby, ecosystems might release more carbon dioxide. In contrast, the CUE of deciduous broad-leaved forest and mixed forest increased, suggesting that these ecosystems were more efficient in assimilating CO2 and reducing the CO2 in the atmosphere. Given the greater absorption of carbon dioxide, carbon sequestration was more likely.


### 3.2. Spatial Distribution and Change Trend of WUE and CUE

Figures4(a) and 4(b) show the spatial distributions of multiyear average WUE and CUE, respectively, in the SNP from 2001 to 2015. The linear trends of the WUE and CUE were also calculated at the 99% confidence level (Figures 4(c) and 4(d)). Spatially, average WUE ranged from 0.178 to 8.445 gC·kg−1H2O; higher values were mainly found in the eastern part of the SNP. The average CUE values in the southwest were greater than in other regions.Figure 4
Spatial distribution and trend of WUE and CUE in the SNP (2001–2015). (a) Annual average WUE, (b) annual average CUE, (c) WUE change trend, and (d) CUE change trend.
(a)
(b)
(c)
(d)The WUE in the SNP showed a general downward trend during the past 15 years. The area of decreased WUE occupied 99.82% of the landscape, showing a concentrated spatial pattern. Among those regions where WUE decreased, cropland and grassland constituted 81.14% and 12.06% of the total land, respectively. The pixels of increased WUE showing in a dotted pattern only accounted for 0.18% of the study area, which was mostly covered by deciduous coniferous forest, deciduous broad-leaved forest, and mixed forest in the northeastern and southwestern fringes.CUE also showed a decreasing trend. Increased CUE was mainly observed in the eastern deciduous coniferous forest, deciduous broad-leaved forest, and mixed forest area, accounting for 0.11%. Ecosystem CUE tended to decrease in 99.89% of the total area. Among those areas, 81.19% and 12.07% were occupied by cropland and grassland, respectively.


### 3.3. Direct Effects of Local Climate Factors on WUE and CUE Changes

The spatial distribution of the correlation coefficients between climatic factors, WUE, and CUE is illustrated in Figures5 and 6. In general, WUE and CUE had more obvious correlation with yearly average precipitation, annual total sunshine hours, average relative humidity, and annual average temperature (Table 1).Figure 5
Correlation coefficients between WUE and major climatic factors: (a) annual precipitation, (b) total annual sunshine duration, (c) annual average relative humidity, and (d) annual average temperature.
(a)
(b)
(c)
(d)Figure 6
Correlation coefficients between CUE and major climatic factors: (a) annual precipitation, (b) total annual sunshine duration, (c) annual average relative humidity, and (d) annual average temperature.
(a)
(b)
(c)
(d)Table 1
Area percentages of the correlation coefficients between WUE, CUE, and climate factors (significant positive correlation (r > 0, p<0.05): no significant positive correlation (r > 0, p<0.05), no significant negative correlation (r < 0, p>0.05), and significant negative correlation (r < 0, p<0.05)).
Correlation
Precipitation
Sunshine duration
Relative humidity
Temperature
WUE
Significant positive correlation
8.21
1.99
17.48
0.27
No significant positive correlation
67.18
32.11
73.28
33.62
No significant negative correlation
24.06
59.15
9.22
64.54
Significant negative correlation
0.55
6.75
0.02
1.57
CUE
Significant positive correlation
10.45
1.59
1.42
0.04
No significant positive correlation
61.82
41.45
59.54
12.98
No significant negative correlation
26.29
54.39
38.61
80.30
Significant negative correlation
1.44
2.57
0.43
6.68WUE was positively correlated with average annual precipitation in two-thirds of the total study area. Precipitation had a positive effect on WUE in more than 70% of the areas covered by mixed forest, deciduous broad-leaved forest, and cropland but was negatively correlated to WUE in 54% of the grassland areas. The duration of sunshine had an inhibitory effect on the WUE in 58.63% of the vegetation areas. The longer the total sunshine duration was, the lower WUE might be. A positive correlation between the total sunshine hours and WUE was mainly observed in the southeastern part of the SNP, but the duration of sunshine in these regions had decreased during the past 15 years. Sunshine duration was significantly and negatively correlated with WUE in the eastern area, where the total sunshine hours increased. The duration of sunshine had a negative impact on more than 60% of wetland WUE. In all, 90% of the areas in the SNP showed a positive correlation between annual average relative humidity and WUE. The positive correlation between relative humidity and WUE was mainly found in the northern part of the SNP, while there was negative correlation in the southeastern and central regions. The annual average humidity had a positive impact on approximately 98% of mixed forest WUE, and negative correlations were observed in 24% of the areas covered by deciduous coniferous forest. The annual average temperature also played a role in regulating WUE. In more than 65% of the SNP, ecosystem WUE had a negative feedback on temperature. The negative correlation areas were mainly observed in the southeast. Except for the deciduous coniferous forest, average temperature had a negative effect on the WUE of other ecosystem types.The average annual precipitation in the SNP had a positive impact on CUE in approximately 72% of the total landscape and especially in mixed forest and wetland. However, negative correlation was found in approximately 52% of areas covered by grassland. The annual total sunshine duration was negatively correlated with ecosystem CUE in the eastern region, whereas it had a positive influence on the CUE of ecosystems in the southeast region. Sunshine duration increased in approximately 85% of areas covered by mixed forest where the CUE was relatively higher. More than 65% of the wetlands showed increased sunshine, while CUE declined. The annual relative humidity was positively correlated with ecosystem CUE in the Northeast fringe and the southwestern region. By contrast, CUE values in more than 60% of the pixels had positive responses to annual relative humidity. The annual average temperature had a negative correlation with the CUE of ecosystems over most of the area (87% of the pixels).Overall, it might be concluded that annual precipitation could improve WUE and CUE of the forestland but inhibited the grasslands’ WUE and CUE. The total annual sunshine duration had a negative impact on wetland WUE and CUE changes. The annual relative humidity was not only positively correlated with the CUE of all the vegetation types in the SNP but also played a positive role in the WUE of most forestland. This indicated that the forestland WUE was improved with increased rainfall. Therefore, the net productivity generated by limited water resources would be increased. The lower wetlands WUE and CUE were potentially correlated with longer sunshine hours. The respiration of vegetation might be enhanced, and more carbon dioxide would be released. In areas with relatively lower temperatures, greater WUE and CUE of vegetation generally could be observed while vegetation respiration decreased.


### 3.4. Human Influence on WUE and CUE Changes

To analyze the impact of human activities on WUE and CUE changes, the residual trend method was adopted. Based on the results of correlation analysis, there was a weak linkage between temperature and WUE and CUE variations at an annual scale. We selected precipitation, total annual sunshine duration, and annual average relative humidity to perform a regression analysis. Pixel-scale predicted WUE and CUE were calculated by a linear regression model using WUE, CUE, and the above three climatic factor maps. The predicted WUE or CUE indicated the climatic impacts, whereas their observed values were the combined results of climatic and anthropogenic factors. Thus, the residual value between the predicted and observed values might be regarded as the human influence on WUE or CUE.According to Figure7, the regional average residual value of WUE in the SNP declined with an undulating trend during the past 15 years. It kept increasing in the first two years and started to decrease from the year 2003. During 2001–2015, the residual was positive (ε > 0), indicating that human activities played a positive role in WUE changes. The decreasing trend in WUE residual might be associated with a weakening human impact. Except for 2003, the CUE residual values were negative, implying the negative effects of anthropogenic factors on CUE changes. The declined trend of the CUE residual suggested that the effects of human activities might decrease.Figure 7
Trends of WUE and CUE residuals in the SNP from 2001 to 2015.Spatially, the regions showing positive trends of WUE and CUE residuals accounted for 0.49% and 0.4% of the total area, indicating that the impacts of human activities on ecosystem WUE and CUE might increase (Figures8(a) and 8(b)). Conversely, the areas showing negative trends of residuals occupied more than 99% of the SNP, implying the human effects on ecosystem WUE and CUE tended to weaken. The pixels showing declined residuals were mainly found in the western part of the study area, dominated by cropland, grassland, and wetland. The trends of WUE and CUE residuals passed the significance level test at p<0.1 and p<0.05, respectively.Figure 8
Spatial distribution of WUE and CUE residuals in the SNP from 2001 to 2015: (a) spatial distribution of WUE residual changes and (b) spatial distribution of CUE residual changes.
(a)
(b)


## 3.1. Temporal Dynamics of WUE and CUE

### 3.1.1. Interannual Changes

The multiyear average WUE in the SNP fluctuated with time, as shown in Figure2. Average WUE peaked in 2003 (3.881 gC·kg−1H2O) in association with higher NPP. Subsequently, the regional WUE kept decreasing and dropped to the lowest value in 2009 (0.551 gC·kg−1H2O), possibly due to high evapotranspiration in that year. WUE started to increase in 2009 and rose to 1.127 gC·kg−1H2O in 2013. It declined again after 2013, and a lower WUE value of 0.940 gC·kg−1H2O was found in 2015.Figure 2
Annual average WUE and CUE in the SNP for the period of 2001–2015.By comparison, the regional CUE was relatively stable, with an average value of 0.516. The average CUE was 0.495 in 2001 and dropped to 0.467 in 2003. It showed an upward trend for the period 2001–2005. The maximum CUE (0.552) was observed in 2005. It decreased to its lowest value (0.461) in 2008, implying the worst level of CUE. Ecosystem CUE increased from 2010 to 2013 and decreased during 2013–2015.


### 3.1.2. WUE Variations of Different Ecosystems

During 2001–2015, the average WUE of different ecosystems followed a descending order of deciduous broadleaf forest > mixed forest > grassland > deciduous coniferous forest > cropland > wetland. The highest WUE in the SNP was found in the deciduous broad-leaved forest, with a multiyear average value of 1.328 gC·kg−1H2O. The wetland showed the lowest WUE (0.954 gC·kg−1H2O). The maximum WUE values of deciduous broad-leaved forest and mixed forest were both found in 2013 (2.019 gC·kg−1H2O and 1.669 gC·kg−1H2O, respectively), whereas the minimum values were 0.577 gC·kg−1H2O and 0.685 gC·kg−1H2O in 2001, respectively. The average WUE of mixed forest was 1.158 gC·kg−1H2O. The interannual average WUEs of grassland and cropland were 0.984 gC·kg−1H2O and 0.961 gC·kg−1H2O, respectively. In deciduous coniferous forest, the multiyear average WUE was 0.972 gC·kg−1H2O, with the highest value (1.304 gC·kg−1H2O) in 2014 and the lowest value (0.584 gC·kg−1H2O) in 2009. The interannual WUE mean of wetland was 0.954 gC·kg−1H2O, and the maximum and minimum WUEs were observed in 2005 (1.138 gC·kg−1H2O) and 2009 (0.615 gC·kg−1H2O), separately. Due to the significantly increased evapotranspiration in 2009, the average WUEs of all ecosystems were relatively low.Regarding the temporal variation of WUE in different vegetation types, the decreasing trend followed the order of deciduous coniferous forest (2.43%) > grassland (1.53%) > wetland (1.37%) > cropland (0.42%). In contrast, the increased trends of WUE were for deciduous broadleaf forest (6.98%) and mixed forest (4.53%). The increased WUE in deciduous broadleaf forest and mixed forest indicated that WUE improved and that greater net productivity might be generated by using a limited water supply. The WUE of deciduous coniferous forest, grassland, wetland, and cropland revealed the opposite trend.


### 3.1.3. CUE Variations of Different Ecosystems

Figure3(b) illustrates the CUE changes of different ecosystems from 2001 to 2015. Grassland showed the highest CUE, with a multiyear average CUE of 0.567. The CUE of deciduous broad-leaved forest was the lowest. CUE of different ecosystems followed a descending order of grassland > wetland > deciduous coniferous forest > cropland > mixed forest > deciduous broad-leaved forest. The multiyear average CUE of deciduous coniferous forest was 0.533. The highest (0.597) and lowest (0.472) CUEs were observed in 2013 and 2003, respectively. Deciduous broad-leaved forest and mixed forest showed a similar change tendency of annual CUE, with average interannual values of 0.479 and 0.480, respectively. The highest CUE (0.560) value occurred in 2005, while their lowest values of 0.239 and 0.308 were found in 2001.Figure 3
WUE (a) and CUE (b) of different ecosystems in the SNP from 2001 to 2015.
(a)
(b)The annual average CUE values of the deciduous broad-leaved forest and the mixed forest were low before 2003, but they suddenly increased in 2004 (Figure3(b)). This might be attributed to the lower annual NPP during 2001–2003, whereas GPP showed a small difference. The average CUE of cropland was 0.510, with the peak of 0.551 in 2003 and the trough of 0.452 in 2008. The multiyear average CUE of wetland was 0.542. It was found that wetland CUE peaked in 2002 (0.801) and that the minimum value occurred in 2007 (0.472). The reason was that NPP of wetland in 2002 reached 579.2 g/m2, an increase of 27.56% from that in 2001.CUE can quantitatively depict the ability of vegetation to fix atmospheric CO2 in the surface ecosystem [34–36]. During the past 15 years, the CUE of deciduous broad-leaved forest and mixed forest showed an increasing trend, with that for deciduous broad-leaved forest rising most obviously. However, average CUE values of deciduous coniferous forest, grassland, cropland, and wetland tended to decrease in the following sequence: deciduous coniferous forest > wetland > grassland > cropland. From 2001 to 2015, deciduous coniferous forest, grassland, cropland, and wetland CUE dropped, indicating their enhanced vegetation respiration. The respiration decomposed organic matter into simple inorganic substances such as water and carbon dioxide; thereby, ecosystems might release more carbon dioxide. In contrast, the CUE of deciduous broad-leaved forest and mixed forest increased, suggesting that these ecosystems were more efficient in assimilating CO2 and reducing the CO2 in the atmosphere. Given the greater absorption of carbon dioxide, carbon sequestration was more likely.


## 3.1.1. Interannual Changes

The multiyear average WUE in the SNP fluctuated with time, as shown in Figure2. Average WUE peaked in 2003 (3.881 gC·kg−1H2O) in association with higher NPP. Subsequently, the regional WUE kept decreasing and dropped to the lowest value in 2009 (0.551 gC·kg−1H2O), possibly due to high evapotranspiration in that year. WUE started to increase in 2009 and rose to 1.127 gC·kg−1H2O in 2013. It declined again after 2013, and a lower WUE value of 0.940 gC·kg−1H2O was found in 2015.Figure 2
Annual average WUE and CUE in the SNP for the period of 2001–2015.By comparison, the regional CUE was relatively stable, with an average value of 0.516. The average CUE was 0.495 in 2001 and dropped to 0.467 in 2003. It showed an upward trend for the period 2001–2005. The maximum CUE (0.552) was observed in 2005. It decreased to its lowest value (0.461) in 2008, implying the worst level of CUE. Ecosystem CUE increased from 2010 to 2013 and decreased during 2013–2015.


## 3.1.2. WUE Variations of Different Ecosystems

During 2001–2015, the average WUE of different ecosystems followed a descending order of deciduous broadleaf forest > mixed forest > grassland > deciduous coniferous forest > cropland > wetland. The highest WUE in the SNP was found in the deciduous broad-leaved forest, with a multiyear average value of 1.328 gC·kg−1H2O. The wetland showed the lowest WUE (0.954 gC·kg−1H2O). The maximum WUE values of deciduous broad-leaved forest and mixed forest were both found in 2013 (2.019 gC·kg−1H2O and 1.669 gC·kg−1H2O, respectively), whereas the minimum values were 0.577 gC·kg−1H2O and 0.685 gC·kg−1H2O in 2001, respectively. The average WUE of mixed forest was 1.158 gC·kg−1H2O. The interannual average WUEs of grassland and cropland were 0.984 gC·kg−1H2O and 0.961 gC·kg−1H2O, respectively. In deciduous coniferous forest, the multiyear average WUE was 0.972 gC·kg−1H2O, with the highest value (1.304 gC·kg−1H2O) in 2014 and the lowest value (0.584 gC·kg−1H2O) in 2009. The interannual WUE mean of wetland was 0.954 gC·kg−1H2O, and the maximum and minimum WUEs were observed in 2005 (1.138 gC·kg−1H2O) and 2009 (0.615 gC·kg−1H2O), separately. Due to the significantly increased evapotranspiration in 2009, the average WUEs of all ecosystems were relatively low.Regarding the temporal variation of WUE in different vegetation types, the decreasing trend followed the order of deciduous coniferous forest (2.43%) > grassland (1.53%) > wetland (1.37%) > cropland (0.42%). In contrast, the increased trends of WUE were for deciduous broadleaf forest (6.98%) and mixed forest (4.53%). The increased WUE in deciduous broadleaf forest and mixed forest indicated that WUE improved and that greater net productivity might be generated by using a limited water supply. The WUE of deciduous coniferous forest, grassland, wetland, and cropland revealed the opposite trend.


## 3.1.3. CUE Variations of Different Ecosystems

Figure3(b) illustrates the CUE changes of different ecosystems from 2001 to 2015. Grassland showed the highest CUE, with a multiyear average CUE of 0.567. The CUE of deciduous broad-leaved forest was the lowest. CUE of different ecosystems followed a descending order of grassland > wetland > deciduous coniferous forest > cropland > mixed forest > deciduous broad-leaved forest. The multiyear average CUE of deciduous coniferous forest was 0.533. The highest (0.597) and lowest (0.472) CUEs were observed in 2013 and 2003, respectively. Deciduous broad-leaved forest and mixed forest showed a similar change tendency of annual CUE, with average interannual values of 0.479 and 0.480, respectively. The highest CUE (0.560) value occurred in 2005, while their lowest values of 0.239 and 0.308 were found in 2001.Figure 3
WUE (a) and CUE (b) of different ecosystems in the SNP from 2001 to 2015.
(a)
(b)The annual average CUE values of the deciduous broad-leaved forest and the mixed forest were low before 2003, but they suddenly increased in 2004 (Figure3(b)). This might be attributed to the lower annual NPP during 2001–2003, whereas GPP showed a small difference. The average CUE of cropland was 0.510, with the peak of 0.551 in 2003 and the trough of 0.452 in 2008. The multiyear average CUE of wetland was 0.542. It was found that wetland CUE peaked in 2002 (0.801) and that the minimum value occurred in 2007 (0.472). The reason was that NPP of wetland in 2002 reached 579.2 g/m2, an increase of 27.56% from that in 2001.CUE can quantitatively depict the ability of vegetation to fix atmospheric CO2 in the surface ecosystem [34–36]. During the past 15 years, the CUE of deciduous broad-leaved forest and mixed forest showed an increasing trend, with that for deciduous broad-leaved forest rising most obviously. However, average CUE values of deciduous coniferous forest, grassland, cropland, and wetland tended to decrease in the following sequence: deciduous coniferous forest > wetland > grassland > cropland. From 2001 to 2015, deciduous coniferous forest, grassland, cropland, and wetland CUE dropped, indicating their enhanced vegetation respiration. The respiration decomposed organic matter into simple inorganic substances such as water and carbon dioxide; thereby, ecosystems might release more carbon dioxide. In contrast, the CUE of deciduous broad-leaved forest and mixed forest increased, suggesting that these ecosystems were more efficient in assimilating CO2 and reducing the CO2 in the atmosphere. Given the greater absorption of carbon dioxide, carbon sequestration was more likely.


## 3.2. Spatial Distribution and Change Trend of WUE and CUE

Figures4(a) and 4(b) show the spatial distributions of multiyear average WUE and CUE, respectively, in the SNP from 2001 to 2015. The linear trends of the WUE and CUE were also calculated at the 99% confidence level (Figures 4(c) and 4(d)). Spatially, average WUE ranged from 0.178 to 8.445 gC·kg−1H2O; higher values were mainly found in the eastern part of the SNP. The average CUE values in the southwest were greater than in other regions.Figure 4
Spatial distribution and trend of WUE and CUE in the SNP (2001–2015). (a) Annual average WUE, (b) annual average CUE, (c) WUE change trend, and (d) CUE change trend.
(a)
(b)
(c)
(d)The WUE in the SNP showed a general downward trend during the past 15 years. The area of decreased WUE occupied 99.82% of the landscape, showing a concentrated spatial pattern. Among those regions where WUE decreased, cropland and grassland constituted 81.14% and 12.06% of the total land, respectively. The pixels of increased WUE showing in a dotted pattern only accounted for 0.18% of the study area, which was mostly covered by deciduous coniferous forest, deciduous broad-leaved forest, and mixed forest in the northeastern and southwestern fringes.CUE also showed a decreasing trend. Increased CUE was mainly observed in the eastern deciduous coniferous forest, deciduous broad-leaved forest, and mixed forest area, accounting for 0.11%. Ecosystem CUE tended to decrease in 99.89% of the total area. Among those areas, 81.19% and 12.07% were occupied by cropland and grassland, respectively.


## 3.3. Direct Effects of Local Climate Factors on WUE and CUE Changes

The spatial distribution of the correlation coefficients between climatic factors, WUE, and CUE is illustrated in Figures5 and 6. In general, WUE and CUE had more obvious correlation with yearly average precipitation, annual total sunshine hours, average relative humidity, and annual average temperature (Table 1).Figure 5
Correlation coefficients between WUE and major climatic factors: (a) annual precipitation, (b) total annual sunshine duration, (c) annual average relative humidity, and (d) annual average temperature.
(a)
(b)
(c)
(d)Figure 6
Correlation coefficients between CUE and major climatic factors: (a) annual precipitation, (b) total annual sunshine duration, (c) annual average relative humidity, and (d) annual average temperature.
(a)
(b)
(c)
(d)Table 1
Area percentages of the correlation coefficients between WUE, CUE, and climate factors (significant positive correlation (r > 0, p<0.05): no significant positive correlation (r > 0, p<0.05), no significant negative correlation (r < 0, p>0.05), and significant negative correlation (r < 0, p<0.05)).
Correlation
Precipitation
Sunshine duration
Relative humidity
Temperature
WUE
Significant positive correlation
8.21
1.99
17.48
0.27
No significant positive correlation
67.18
32.11
73.28
33.62
No significant negative correlation
24.06
59.15
9.22
64.54
Significant negative correlation
0.55
6.75
0.02
1.57
CUE
Significant positive correlation
10.45
1.59
1.42
0.04
No significant positive correlation
61.82
41.45
59.54
12.98
No significant negative correlation
26.29
54.39
38.61
80.30
Significant negative correlation
1.44
2.57
0.43
6.68WUE was positively correlated with average annual precipitation in two-thirds of the total study area. Precipitation had a positive effect on WUE in more than 70% of the areas covered by mixed forest, deciduous broad-leaved forest, and cropland but was negatively correlated to WUE in 54% of the grassland areas. The duration of sunshine had an inhibitory effect on the WUE in 58.63% of the vegetation areas. The longer the total sunshine duration was, the lower WUE might be. A positive correlation between the total sunshine hours and WUE was mainly observed in the southeastern part of the SNP, but the duration of sunshine in these regions had decreased during the past 15 years. Sunshine duration was significantly and negatively correlated with WUE in the eastern area, where the total sunshine hours increased. The duration of sunshine had a negative impact on more than 60% of wetland WUE. In all, 90% of the areas in the SNP showed a positive correlation between annual average relative humidity and WUE. The positive correlation between relative humidity and WUE was mainly found in the northern part of the SNP, while there was negative correlation in the southeastern and central regions. The annual average humidity had a positive impact on approximately 98% of mixed forest WUE, and negative correlations were observed in 24% of the areas covered by deciduous coniferous forest. The annual average temperature also played a role in regulating WUE. In more than 65% of the SNP, ecosystem WUE had a negative feedback on temperature. The negative correlation areas were mainly observed in the southeast. Except for the deciduous coniferous forest, average temperature had a negative effect on the WUE of other ecosystem types.The average annual precipitation in the SNP had a positive impact on CUE in approximately 72% of the total landscape and especially in mixed forest and wetland. However, negative correlation was found in approximately 52% of areas covered by grassland. The annual total sunshine duration was negatively correlated with ecosystem CUE in the eastern region, whereas it had a positive influence on the CUE of ecosystems in the southeast region. Sunshine duration increased in approximately 85% of areas covered by mixed forest where the CUE was relatively higher. More than 65% of the wetlands showed increased sunshine, while CUE declined. The annual relative humidity was positively correlated with ecosystem CUE in the Northeast fringe and the southwestern region. By contrast, CUE values in more than 60% of the pixels had positive responses to annual relative humidity. The annual average temperature had a negative correlation with the CUE of ecosystems over most of the area (87% of the pixels).Overall, it might be concluded that annual precipitation could improve WUE and CUE of the forestland but inhibited the grasslands’ WUE and CUE. The total annual sunshine duration had a negative impact on wetland WUE and CUE changes. The annual relative humidity was not only positively correlated with the CUE of all the vegetation types in the SNP but also played a positive role in the WUE of most forestland. This indicated that the forestland WUE was improved with increased rainfall. Therefore, the net productivity generated by limited water resources would be increased. The lower wetlands WUE and CUE were potentially correlated with longer sunshine hours. The respiration of vegetation might be enhanced, and more carbon dioxide would be released. In areas with relatively lower temperatures, greater WUE and CUE of vegetation generally could be observed while vegetation respiration decreased.


## 3.4. Human Influence on WUE and CUE Changes

To analyze the impact of human activities on WUE and CUE changes, the residual trend method was adopted. Based on the results of correlation analysis, there was a weak linkage between temperature and WUE and CUE variations at an annual scale. We selected precipitation, total annual sunshine duration, and annual average relative humidity to perform a regression analysis. Pixel-scale predicted WUE and CUE were calculated by a linear regression model using WUE, CUE, and the above three climatic factor maps. The predicted WUE or CUE indicated the climatic impacts, whereas their observed values were the combined results of climatic and anthropogenic factors. Thus, the residual value between the predicted and observed values might be regarded as the human influence on WUE or CUE.According to Figure7, the regional average residual value of WUE in the SNP declined with an undulating trend during the past 15 years. It kept increasing in the first two years and started to decrease from the year 2003. During 2001–2015, the residual was positive (ε > 0), indicating that human activities played a positive role in WUE changes. The decreasing trend in WUE residual might be associated with a weakening human impact. Except for 2003, the CUE residual values were negative, implying the negative effects of anthropogenic factors on CUE changes. The declined trend of the CUE residual suggested that the effects of human activities might decrease.Figure 7
Trends of WUE and CUE residuals in the SNP from 2001 to 2015.Spatially, the regions showing positive trends of WUE and CUE residuals accounted for 0.49% and 0.4% of the total area, indicating that the impacts of human activities on ecosystem WUE and CUE might increase (Figures8(a) and 8(b)). Conversely, the areas showing negative trends of residuals occupied more than 99% of the SNP, implying the human effects on ecosystem WUE and CUE tended to weaken. The pixels showing declined residuals were mainly found in the western part of the study area, dominated by cropland, grassland, and wetland. The trends of WUE and CUE residuals passed the significance level test at p<0.1 and p<0.05, respectively.Figure 8
Spatial distribution of WUE and CUE residuals in the SNP from 2001 to 2015: (a) spatial distribution of WUE residual changes and (b) spatial distribution of CUE residual changes.
(a)
(b)


## 4. Discussion

### 4.1. Variations of WUE and CUE

As an important indicator of terrestrial ecosystems responding to global change, WUE links carbon and hydrological cycles of vegetation ecosystems [37]. Due to the varied photosynthesis capabilities of different types of vegetation, there should exist certain differences in WUE and CUE [38]. We compared the WUE estimated by different research studies (Table 2). Lu and Li quantitatively estimated and analyzed the spatial and temporal dynamic patterns and seasonal variations of WUE for different vegetation ecosystems in western China in 2002. They found that the average annual WUE of the main vegetation ecosystems in this region showed a downward trend as follows: mountain forest > desert shrub > irrigated farmland > alpine grassland > cold desert and Gobi [42]. Qiu et al. reported that the annual average WUE of woodland was significantly higher than for other vegetation types in China. The annual average WUE of the mixed forest and deciduous broad-leaved forest was the highest, while the annual average WUE of the shrubs and grasslands was lower (<1 gC·kg−1H2O) [43]. In our current study, the annual average WUE of the main vegetation ecosystem in the SNP followed a descending order of deciduous broad-leaved forest > mixed forest > grassland > deciduous coniferous forest > cropland > wetland. The regional average WUE in the SNP was 1.073 gC·kg−1H2O for the period of 2001–2015, which was consistent with the threshold of estimated WUE in China [43]. Liu et al. found that annual average WUE of the forest or grassland in Northeast China showed less interannual variability during 2002–2013 [44]. In this temperate semihumid to semiarid transitional region of Northeast China, we found that the interannual differences in the annual average WUE were relatively large from 2001 to 2015. The average WUE of the cropland decreased by 1%, while the average WUE of mixed forest, deciduous broad-leaved forest, and deciduous coniferous forest increased by 40%, 32.9%, and 19.6%, respectively.Table 2
Comparison of WUE estimated by different researches.
Vegetation type
WUE (gC·kg−1H2O)
Data source
Wetland
—
Wei [39]
—
DLEM model [40]
0.869
Xia [41]
0.954
Our article
Cropland
0.82
Wei
0.54
DLEM model
0.908
Xia
0.961
Our article
Grassland
0.77
Wei
0.75
DLEM model
0.827
Xia
0.978
Our article
Mixed forest
0.85
Wei
0.96
DLEM model
1.077
Xia
1.158
Our article
Deciduous broad-leaved forest
1.04
Wei
0.91
DLEM model
1.018
Xia
1.328
Our article
Deciduous coniferous forest
0.71
Wei
0.96
DLEM model
0.943
Xia
0.972
Our articleCUE is sensitive to environmental conditions and climate change [45]. As the function of GPP, NPP, and respiration, the CUE of vegetation (for instance, forest) may be affected by many factors such as temperature, precipitation, photosynthetically active radiation, and soil nutrients [46]. CUE showed significant spatial variations associated with different ecosystems, geographical location, and climate [20]. Chambers et al. revealed that the CUE of temperate forests was approximately 0.50 [47]. Khalifa et al. estimated the CUE of different vegetation in the sub-Saharan area and found that the annual average CUE decreased in the following sequence: cropland > wetland > grassland > mixed forest > evergreen broad-leaved forest > deciduous broad-leaved forest [21]. Our results showed that the average annual CUE of different vegetation in the SNP followed a descending order of grassland (0.567) > wetland (0.542) > deciduous coniferous forest (0.533) > cropland (0.510) > mixed forest (0.480) > deciduous broad-leaved forest (0.479). The CUE of grassland was higher than those of the deciduous coniferous and broad-leaved forest, possibly due to the lesser investment in respiring plant tissue relative to forests; the same result was reported by Law et al. [48]. The annual average CUE of the SNP was 0.520 and showed a steady trend, which was consistent with the threshold value of the global terrestrial ecosystem CUE estimated by Zhang et al. [20]. During the past 15 years, the increase in CUE of forest ecosystems in the SNP could be a positive sign of an enhancement of their efficiency, which has great potential for carbon sequestration. However, the decreased CUE of grassland, cropland, and wetland ecosystems in the SNP indicated their diminished efficiency in assimilating CO2 and reducing the CO2 in the atmosphere for the period of 2001–2015.


### 4.2. Relationships between WUE and CUE Changes, Climate Factors, and Human Activity

Climate factors could affect the changes of carbon fluxes. There have been different conclusions about the effect of precipitation on vegetation WUE to date. In the current study, we found that the WUE of ecosystems in the SNP was significantly and positively correlated with annual average precipitation and relative humidity, while it was negatively correlated with total sunshine duration and annual average temperature. In the regions with more annual average precipitation, annual total sunshine hours tend to be shorter, possibly resulting in less evapotranspiration and higher WUE. Han and Wan noted that temperature had an inhibitory effect on WUE in the grassland region of northern China, whereas precipitation might promote WUE on the ecosystem scale [49]. The annual average relative humidity can reflect the annual average precipitation to a certain extent, and its impact mechanism on vegetation WUE and CUE is similar to that of precipitation. Jiang found that, as the drought increased, the WUE of vegetation in the Northeast China Transect (NECT) gradually increased to a certain value and then decreased, indicating that there was a specific threshold of plant WUE in terms of the drought degree [50]. However, few studies on the relationship between precipitation and WUE at a regional scale in China have been conducted. This may result in a vague understanding of whether ecosystem WUE in China shows a trend of increasing first and then decreasing with the increase of annual precipitation. It is also unclear how high the annual precipitation will be when WUE reaches its peak [51].Zhang et al. indicated that CUE tended to decrease with the increase of precipitation in the regions with annual precipitation less than 2300 mm at the global scale. Illumination is one of the driving forces of plant photosynthesis. Photosynthesis of plants can only be performed within a certain range of radiation [20]. With light intensity decreasing, the maintenance respiratory rate decreased, while the growth respiration coefficient increased, resulting in the decrease of plant CUE [52]. Marthews et al. simulated the carbon balance of six tropical forest plots in the Andes-Amazon. They indicated that forest CUE did not show a consistent trend with changing altitude and temperature. The CUE value was close to a median of 0.5, which was mainly related to the balance between growth respiration and maintenance respiration of the forest canopy [53]. Zhang et al. analyzed the global pattern of NPP to GPP ratio and noted that the ratio might increase as the temperature decreased between 20°C and 27°C [20], which is roughly similar to our finding.Human interference and management measures may also have an impact on WUE and CUE changes. According to the study of Liu and Ren, the increases of vegetation coverage and WUE in the Loess Plateau were mainly attributable to ecological restoration projects such as returning cropland to forest [54]. Hirata et al. thought that carbon sequestration of forest ecosystems was primarily sensitive to the influence of climate, given the low intensity of human disturbances [55]. On the other hand, human activities such as fertilization and irrigation may affect vegetation CUE. Ryan et al. demonstrated that biomass and respiration of plant fine roots decreased due to irrigation and fertilization, and forest CUE increased through experiments on a 20-year-old Pinus radiata plantation [56]. The increase in soil nutrient use efficiency may reduce the growth of fine roots, promote dry wood production, and increase the aboveground part of CUE, suggesting that human intervention and management measures can cause a transient increase in aboveground CUE [57]. In this study, our results revealed the opposite effects of human activities on WUE and CUE variations in the SNP. Anthropogenic causes might play positive roles in WUE changes but had negative influences on CUE. Tian et al. indicated that a total of 9155.53 km2 of grassland was converted into cropland in the SNP during 1986–2010 [58]. By using MCD12Q1 land cover data, the major vegetation type changes in the SNP from 2001 to 2015 were extracted (Figure 9). The cropland increased, the grassland decreased, and the forest remained relatively stable. It was found that the grassland, accounting for approximately 4% of the total pixels, was converted into cropland, which might provide evidence for our hypothesis. The average CUE of grasslands was higher. When the grassland was reclaimed as cropland, the regional CUE thereby decreased. Spatially, due to the high percentage of areas showing the negative residual trends of WUE and CUE, the intensity of human impact tended to weaken in the SNP during the past 15 years.Figure 9
Changes in major vegetation types in the SNP from 2001 to 2015.


### 4.3. Uncertainty

This study investigated spatial and temporal patterns of WUE and CUE in SNP, using time-series remote sensing data with a spatial resolution of 1 km. Generally, uncertainty may be involved while estimating WUE and CUE using public-domain data [59]. The relatively lower spatial resolution of the available remote sensing data may affect the accuracy, especially for high heterogeneous ecosystems. Due to the lack of ground-based measurements, the regional WUE and CUE estimates in this study could not be accurately validated. Precipitation, temperature, sunshine duration, and relative humidity were chosen as the influencing climate factors of regional ecosystem WUE and CUE variations, but solar radiation and CO2 concentration were not taken into account. The mechanism of the environmental controls of variations in ecosystem WUE and CUE is more complicated. Although a residual analysis may be used to distinguish human influence from climate change influence on WUE and CUE, the effects of other factors, such as management and technology, need to be evaluated further.


## 4.1. Variations of WUE and CUE

As an important indicator of terrestrial ecosystems responding to global change, WUE links carbon and hydrological cycles of vegetation ecosystems [37]. Due to the varied photosynthesis capabilities of different types of vegetation, there should exist certain differences in WUE and CUE [38]. We compared the WUE estimated by different research studies (Table 2). Lu and Li quantitatively estimated and analyzed the spatial and temporal dynamic patterns and seasonal variations of WUE for different vegetation ecosystems in western China in 2002. They found that the average annual WUE of the main vegetation ecosystems in this region showed a downward trend as follows: mountain forest > desert shrub > irrigated farmland > alpine grassland > cold desert and Gobi [42]. Qiu et al. reported that the annual average WUE of woodland was significantly higher than for other vegetation types in China. The annual average WUE of the mixed forest and deciduous broad-leaved forest was the highest, while the annual average WUE of the shrubs and grasslands was lower (<1 gC·kg−1H2O) [43]. In our current study, the annual average WUE of the main vegetation ecosystem in the SNP followed a descending order of deciduous broad-leaved forest > mixed forest > grassland > deciduous coniferous forest > cropland > wetland. The regional average WUE in the SNP was 1.073 gC·kg−1H2O for the period of 2001–2015, which was consistent with the threshold of estimated WUE in China [43]. Liu et al. found that annual average WUE of the forest or grassland in Northeast China showed less interannual variability during 2002–2013 [44]. In this temperate semihumid to semiarid transitional region of Northeast China, we found that the interannual differences in the annual average WUE were relatively large from 2001 to 2015. The average WUE of the cropland decreased by 1%, while the average WUE of mixed forest, deciduous broad-leaved forest, and deciduous coniferous forest increased by 40%, 32.9%, and 19.6%, respectively.Table 2
Comparison of WUE estimated by different researches.
Vegetation type
WUE (gC·kg−1H2O)
Data source
Wetland
—
Wei [39]
—
DLEM model [40]
0.869
Xia [41]
0.954
Our article
Cropland
0.82
Wei
0.54
DLEM model
0.908
Xia
0.961
Our article
Grassland
0.77
Wei
0.75
DLEM model
0.827
Xia
0.978
Our article
Mixed forest
0.85
Wei
0.96
DLEM model
1.077
Xia
1.158
Our article
Deciduous broad-leaved forest
1.04
Wei
0.91
DLEM model
1.018
Xia
1.328
Our article
Deciduous coniferous forest
0.71
Wei
0.96
DLEM model
0.943
Xia
0.972
Our articleCUE is sensitive to environmental conditions and climate change [45]. As the function of GPP, NPP, and respiration, the CUE of vegetation (for instance, forest) may be affected by many factors such as temperature, precipitation, photosynthetically active radiation, and soil nutrients [46]. CUE showed significant spatial variations associated with different ecosystems, geographical location, and climate [20]. Chambers et al. revealed that the CUE of temperate forests was approximately 0.50 [47]. Khalifa et al. estimated the CUE of different vegetation in the sub-Saharan area and found that the annual average CUE decreased in the following sequence: cropland > wetland > grassland > mixed forest > evergreen broad-leaved forest > deciduous broad-leaved forest [21]. Our results showed that the average annual CUE of different vegetation in the SNP followed a descending order of grassland (0.567) > wetland (0.542) > deciduous coniferous forest (0.533) > cropland (0.510) > mixed forest (0.480) > deciduous broad-leaved forest (0.479). The CUE of grassland was higher than those of the deciduous coniferous and broad-leaved forest, possibly due to the lesser investment in respiring plant tissue relative to forests; the same result was reported by Law et al. [48]. The annual average CUE of the SNP was 0.520 and showed a steady trend, which was consistent with the threshold value of the global terrestrial ecosystem CUE estimated by Zhang et al. [20]. During the past 15 years, the increase in CUE of forest ecosystems in the SNP could be a positive sign of an enhancement of their efficiency, which has great potential for carbon sequestration. However, the decreased CUE of grassland, cropland, and wetland ecosystems in the SNP indicated their diminished efficiency in assimilating CO2 and reducing the CO2 in the atmosphere for the period of 2001–2015.


## 4.2. Relationships between WUE and CUE Changes, Climate Factors, and Human Activity

Climate factors could affect the changes of carbon fluxes. There have been different conclusions about the effect of precipitation on vegetation WUE to date. In the current study, we found that the WUE of ecosystems in the SNP was significantly and positively correlated with annual average precipitation and relative humidity, while it was negatively correlated with total sunshine duration and annual average temperature. In the regions with more annual average precipitation, annual total sunshine hours tend to be shorter, possibly resulting in less evapotranspiration and higher WUE. Han and Wan noted that temperature had an inhibitory effect on WUE in the grassland region of northern China, whereas precipitation might promote WUE on the ecosystem scale [49]. The annual average relative humidity can reflect the annual average precipitation to a certain extent, and its impact mechanism on vegetation WUE and CUE is similar to that of precipitation. Jiang found that, as the drought increased, the WUE of vegetation in the Northeast China Transect (NECT) gradually increased to a certain value and then decreased, indicating that there was a specific threshold of plant WUE in terms of the drought degree [50]. However, few studies on the relationship between precipitation and WUE at a regional scale in China have been conducted. This may result in a vague understanding of whether ecosystem WUE in China shows a trend of increasing first and then decreasing with the increase of annual precipitation. It is also unclear how high the annual precipitation will be when WUE reaches its peak [51].Zhang et al. indicated that CUE tended to decrease with the increase of precipitation in the regions with annual precipitation less than 2300 mm at the global scale. Illumination is one of the driving forces of plant photosynthesis. Photosynthesis of plants can only be performed within a certain range of radiation [20]. With light intensity decreasing, the maintenance respiratory rate decreased, while the growth respiration coefficient increased, resulting in the decrease of plant CUE [52]. Marthews et al. simulated the carbon balance of six tropical forest plots in the Andes-Amazon. They indicated that forest CUE did not show a consistent trend with changing altitude and temperature. The CUE value was close to a median of 0.5, which was mainly related to the balance between growth respiration and maintenance respiration of the forest canopy [53]. Zhang et al. analyzed the global pattern of NPP to GPP ratio and noted that the ratio might increase as the temperature decreased between 20°C and 27°C [20], which is roughly similar to our finding.Human interference and management measures may also have an impact on WUE and CUE changes. According to the study of Liu and Ren, the increases of vegetation coverage and WUE in the Loess Plateau were mainly attributable to ecological restoration projects such as returning cropland to forest [54]. Hirata et al. thought that carbon sequestration of forest ecosystems was primarily sensitive to the influence of climate, given the low intensity of human disturbances [55]. On the other hand, human activities such as fertilization and irrigation may affect vegetation CUE. Ryan et al. demonstrated that biomass and respiration of plant fine roots decreased due to irrigation and fertilization, and forest CUE increased through experiments on a 20-year-old Pinus radiata plantation [56]. The increase in soil nutrient use efficiency may reduce the growth of fine roots, promote dry wood production, and increase the aboveground part of CUE, suggesting that human intervention and management measures can cause a transient increase in aboveground CUE [57]. In this study, our results revealed the opposite effects of human activities on WUE and CUE variations in the SNP. Anthropogenic causes might play positive roles in WUE changes but had negative influences on CUE. Tian et al. indicated that a total of 9155.53 km2 of grassland was converted into cropland in the SNP during 1986–2010 [58]. By using MCD12Q1 land cover data, the major vegetation type changes in the SNP from 2001 to 2015 were extracted (Figure 9). The cropland increased, the grassland decreased, and the forest remained relatively stable. It was found that the grassland, accounting for approximately 4% of the total pixels, was converted into cropland, which might provide evidence for our hypothesis. The average CUE of grasslands was higher. When the grassland was reclaimed as cropland, the regional CUE thereby decreased. Spatially, due to the high percentage of areas showing the negative residual trends of WUE and CUE, the intensity of human impact tended to weaken in the SNP during the past 15 years.Figure 9
Changes in major vegetation types in the SNP from 2001 to 2015.


## 4.3. Uncertainty

This study investigated spatial and temporal patterns of WUE and CUE in SNP, using time-series remote sensing data with a spatial resolution of 1 km. Generally, uncertainty may be involved while estimating WUE and CUE using public-domain data [59]. The relatively lower spatial resolution of the available remote sensing data may affect the accuracy, especially for high heterogeneous ecosystems. Due to the lack of ground-based measurements, the regional WUE and CUE estimates in this study could not be accurately validated. Precipitation, temperature, sunshine duration, and relative humidity were chosen as the influencing climate factors of regional ecosystem WUE and CUE variations, but solar radiation and CO2 concentration were not taken into account. The mechanism of the environmental controls of variations in ecosystem WUE and CUE is more complicated. Although a residual analysis may be used to distinguish human influence from climate change influence on WUE and CUE, the effects of other factors, such as management and technology, need to be evaluated further.


## 5. Conclusions

In this study, spatial and temporal variations of ecosystem WUE and CUE and their relationship with climate and human activity in the SNP were assessed. WUE and CUE differed in different terrestrial ecosystems in the SNP. The multiyear average WUE of the SNP fell in the range of 0.178–8.445 gC·kg−1H2O for the period of 2001–2015, with a mean value of 1.319 gC·kg−1H2O. From 2001 to 2015, regional multiyear average CUE values varied from 0 to 0.801. Average WUE of deciduous broad-leaved forest ranked in first place, whereas grasslands possessed the lowest WUE. In contrast, grassland CUE was the highest, while deciduous broad-leaved forest had the lowest CUE. The decreases in WUE and CUE of grassland, cropland, and wetland indicated that their capability of water resources use weakened, and their efficiency in assimilating CO2 and reducing the CO2 in the atmosphere decreased. Spatially, the ecosystems’ WUE and CUE in almost the entire study area exhibited a decreasing trend over the past 15 years. The regions with a decrease in both WUE and CUE were mainly distributed in cropland and grassland, respectively.In most areas of the SNP, ecosystem WUE and CUE were positively correlated with annual precipitation and annual average relative humidity for the period of 2001–2015. In contrast, total sunshine duration and annual average temperature had significantly negative correlations with ecosystem WUE and CUE. In the past 15 years, human impact might have facilitated the increase of ecosystem WUE while oppositely restraining the CUE in the SNP, but the extent of these effects decreased.



---

*Source: 8659138-2019-04-28.xml*