Next Article in Journal
Feasibility of Bi-Temporal Airborne Laser Scanning Data in Detecting Species-Specific Individual Tree Crown Growth of Boreal Forests
Next Article in Special Issue
Impact of the Management Scale on the Technical Efficiency of Forest Vegetation Carbon Sequestration: A Case Study of State-Owned Forestry Enterprises in Northeast China
Previous Article in Journal
Flood Hazard Analysis Based on Rainfall Fusion: A Case Study in Dazhou City, China
Previous Article in Special Issue
The Different Impacts of Climate Variability and Human Activities on NPP in the Guangdong–Hong Kong–Macao Greater Bay Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantifying the Influence of Climate Change and Anthropogenic Activities on the Net Primary Productivity of China’s Grasslands

1
College of Management and Economics, Tianjin University, Tianjin 300072, China
2
The Center for Eco-Environmental Accounting, Chinese Academy of Environmental Planning, Beijing 100012, China
3
State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
4
School of Natural Resources, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
*
Authors to whom correspondence should be addressed.
Remote Sens. 2022, 14(19), 4844; https://doi.org/10.3390/rs14194844
Submission received: 15 September 2022 / Revised: 24 September 2022 / Accepted: 25 September 2022 / Published: 28 September 2022

Abstract

:
As one of China’s most common vegetation types, grasslands comprise about 27.5% of its terrestrial area and 41% of its carbon storage. Since climate change (CC) and human activities (HA) have a great effect on grasslands, quantifying the contributions of CC and HA on grassland net primary productivity (NPP) is crucial in understanding the mechanisms of grassland regional carbon balances. However, current approaches, including residual trend, biophysical model and environmental background-based methods, have limitations on different scales, especially on the national scale of China. To improve assessment accuracy, modifications to the environmental background-based method were introduced in calculating the CC and HA contributions to the actual NPP (ANPP). In this study, the grassland ANPP in national nature reserves was defined as the environmental background value (PNPP), which was only affected by CC and without HA. The pixel PNPP outside the nature reserves could be replaced by the pixel PNPP in the nature reserve with the most similar habitat in the same natural ecological geographical division. The impact of HA on grassland ANPP (HNPP) could be identified by calculating the difference between PNPP and ANPP. Finally, the contributions of CC and HA to ANPP changes were assessed by the trends of ANPP, PNPP, and HNPP. The results showed that the average grassland ANPP significantly increased from 2001 to 2020. CC contributed 71.0% to ANPP change, whereas HA contributed 29.0%. Precipitation was the main contributor to grassland growth among arid and semi-arid regions, while temperature inhibited productivity in these areas. HA was the major cause of degradation in China’s grasslands, although the effects have declined over time. The research could provide support support for government decisions. It could also provide a new and feasible research method for quantitatively evaluating grasslands and other ecosystems.

Graphical Abstract

1. Introduction

As an essential component in the ecosystem, vegetation is indispensable in maintaining climatic stability and carbon balance and is crucial in achieving carbon emissions neutrality targets [1,2,3]. Net primary productivity (NPP) is the net carbon sequestered via vegetation photosynthesis in unit time and unit area [4,5]. It characterizes the carbon sequestration capacity of vegetation and serves as a vital gauge of plant dynamics [6,7,8]. Climate change (CC) and human activities (HA) are two main influencing aspects of changes in vegetation NPP [9,10]. CC (e.g., precipitation and temperature change) significantly affects the intensity of vegetation photosynthesis and respiration [11,12], leading to an increase or decrease in vegetation NPP [7,13]. HA, such as overgrazing or land use change, may result in severe ecological degradation [9,14]. Therefore, estimating the influence of CC and HA on vegetation NPP is crucial for vegetation management.
Recent studies have investigated the influence of CC and HA on vegetation NPP [1,9,12]. Currently, the most widely used methods mainly include regression analysis methods [15,16], residual trend methods [1,17], biophysical model methods [9,10] and environmental background-based methods [18]. The regression analysis method is relatively simple and directly uses the regression model to quantify the impact of CC and HA on vegetation NPP [19]. However, this approach has difficulty identifying the respective influence of CC and HA on vegetation NPP [20]. The residual trend method indirectly evaluates the effects of anthropogenic activities on vegetation NPP by simulating the residual trend of actual NPP and potential NPP [1,21]. However, this method assumes all vegetation is affected by CC and HA, resulting in deviations from the actual situation [22,23]. For example, HA has no impact on vegetation in national nature reserves [18]. The biophysical model method employs the ThornThwaite Memorial Model and Carnegie-Ames-Stanford Approach in estimating the potential and actual NPP values and then separating the relative contributions of CC and HA in vegetation NPP [24,25]. However, numerous parameters evaluating vegetation physiology and ecology would have to be measured, thus increasing the model uncertainties [10,26]. The environmental background-based method is a new way to separate the impacts of CC and HA on vegetation NPP [18]. In this method, the vegetation NPP in nature reserves is defined as the environmental background value, which is affected only by CC and without HA interference. The impact of HA on vegetation NPP is then characterized by calculating the difference between the environmental background value and the NPP outside the reserve. This method provides a more precise approach but ignores the influence of diverse habitat conditions on vegetation NPP. Improvements and modifications to the environmental background-based approach have to be incorporated to generate more accurate estimates of the influence of CC and HA on vegetation NPP.
Grassland is a major ecosystem in China, accounting for about 27.5% of China’s total land area [27,28]. As a carbon sink, grassland ecosystems play a critical role and represent about 41% of China’s carbon storage [29,30]. However, due to climate warming, land-use changes, as well as overgrazing, many grasslands in China have suffered serious ecological degradation [10,31], causing significant decline in the regional carbon sequestration capacity [32,33]. Grassland degradation has also resulted in a number of environmental problems (e.g., soil erosion and sandstorms) [34,35]. Since the late 1990s, the Chinese government has implemented policies to restore degraded grasslands, such as the Grain-for-Green Project launched in 1999 and Regulations on Conversion of Farmland to Forests in 2016 [31,36]. Succeeding environmental policies and restoration strategies would require a more accurate estimation of the climatic and anthropogenic impacts on grassland NPP, not merely as an academic exercise, but to provide important references for grassland management [10,37].
Against this background, an improved environmental background-based method was developed in this research to identify the key influencing factors of grassland NPP change in China. The objectives of this study were as follows: (i) investigate spatiotemporal dynamics of grassland NPP from 2001 to 2020, (ii) quantify the relative contributions of CC and HA on grassland NPP, and (iii) explore the relationships between grassland NPP and major climate factors.

2. Materials and Methods

2.1. Study Area

China’s grassland ecosystem was chosen as the study area, distributed in the administrative regions of Gansu, Shaanxi, Xinjiang, Tibet, Inner Mongolia, and Hebei (Figure 1) [38,39,40]. These areas are primarily located in arid and semi-arid regions [39]. In terms of overgrazing and land-use change over the years, the ecological degradation in the study area is becoming more and more serious [41,42].

2.2. Data Sources

The main research data was composed of land use, net primary productivity, precipitation, temperature, eco-geographical zoning, and national natural reserves data. The land use data for 2000–2018, the net primary productivity data for 2001–2020 and precipitation and temperature during 2001–2020 were obtained from the European Space Agency, the United States Geological Survey and the China Meteorological Data Network, respectively. The eco-geographical zoning data and DEM data were obtained from the Resource and Environment Science and Data Center (Figure 2). The national natural reserves data was derived from China’s Ministry of Ecology and Environment (Figure 3) and is separated into three zones: core area, buffer zone and experimental area [43]. The core area is strictly protected from any human activity, but in the buffer zone, scientific research and observations are allowed. In the experimental area, visits and investigations, tourism, reproduction of rare and endangered wild animals and plants, and specific production activities are allowed [44]. All data were resampled to a spatial resolution of 1 km×1 km for unified analysis (Table 1).

2.3. Methods

2.3.1. Extracting the Spatial Range of Grassland

Due to continuous changes in the grassland ecosystem from 2001 to 2020, only the unaltered grassland areas were considered in the study. The unchanged grassland areas in 2001, 2006, 2011, 2016, and 2018 were used as the final spatial grassland range (Figure 4).

2.3.2. Calculation of NPP

Improvements to the environmental background-based method were used to calculate the NPP. Three types of NPP were analyzed: actual NPP (ANPP), potential NPP (PNPP), and human-induced NPP (HNPP). The ANPP estimates the combined influence of CC and HA. The PNPP predicts the ideal condition without human interference and is solely determined by climate conditions. The HNPP, representing the influence of HA, can then be calculated using Equation (1):
H N P P = P N P P A N P P
The ANPP was calculated by NASA, based on MOD17A3 data. For the PNPP, since all disturbances are strictly prohibited in the core areas of national nature reserves, the grasslands in the core areas could be assumed to be only affected by CC and free from HA. Therefore, the grassland PNPP in the core area was approximately equal to ANPP. In addition, if the habitats (altitude, slope, aspect, precipitation, temperature, etc.) of the two pixels were completely similar, we could assume that the PNPP of the two pixels could be approximately equal [45], see support information for hypothesis verification. Therefore, the pixel PNPP outside the core area could be replaced by the pixel PNPP in the nature reserve with the most similar habitat in the same natural ecological geographical division. Similarity pixel matching was mainly based on the principle of the similarity search algorithm of ArcGIS and realized by Matlab programming.

2.3.3. Trend Analysis

To analyze the NPP dynamics from 2001 to 2020, the Theil-Sen median slope method was used as the Equation (2) [46,47]. Mann-Kendall (M-K) was utilized to test the significance level of the trend [47].
S l o p e s e n =   Median ( Q i ) =   Median ( x j x i j i )
where xi and xj are the NPP values in the ith and jth years (j > i), respectively.

2.3.4. The Relative Contributions of CC and HA on ANPP

The contribution of CC and HA to NPP could be quantitatively evaluated by the combined slope analysis for the different NPP types [47,48]. A positive SANPP indicated grassland restoration, while a negative value denoted grassland degradation [10]. Using different SPNPP and SHNPP combinations, six scenarios were identified, characterizing the impact of CC and HA (Table 2) [20].

2.3.5. Correlation Analysis

To identify the correlations between ANPP, temperature, and precipitation, multiple correlation analysis and F-test were used [20]. Partial correlation and t-test were utilized to estimate the relationship of ANPP with temperature and rainfall [49].

3. Results

3.1. Spatiotemporal Characteristics of Grassland ANPP

3.1.1. Spatial Heterogeneity Analysis of ANPP

The average grassland ANPP from 2001 to 2020 was 148.63 g C·m−2·yr−1. Regarding the spatial distribution, the grassland ANPP value decreased considerably from the southeast to the northwest of China (Figure 5). The highest ANPP values were mainly distributed in northeastern Inner Mongolia, northeastern Hebei, Shanxi, southeastern Gansu, northeastern Sichuan, northwestern Xinjiang, and some parts of Yunnan (Figure 5).

3.1.2. Spatiotemporal Variation of ANPP

The average grassland ANPP saw a significantly increasing trend (Figure 6), indicating that the synergistic effects of CC and HA led to significant growth of grassland for the given study period. About 89.1% of grassland ANPP exhibited an increasing trend (Figure 7a). The statistically significant increasing trend accounted for 48.3% of grasslands (Figure 7b), mainly distributed in Central Inner Mongolia, Southeastern Qinghai, Northwestern Xinjiang, Northwestern Sichuan, and Northeastern Tibet. In contrast, 10.9% of grassland ANPP showed a declining trend (Figure 7a). Areas witnessing a significantly decreasing trend in ANPP only accounted for 0.6% of grasslands (Figure 7b) and were mainly found in Central and southeastern Tibet, southeastern Qinghai, and northwestern Xinjiang.

3.2. Contributions of CC and HA to ANPP

3.2.1. Changing Trends of PNPP and HNPP

To differentiate the contributions of CC and HA to ANPP, we first calculated the change trends for PNPP and HNPP from 2001 to 2020. The PNPP trend showed that CC positively affected grassland growth, accounting for 92.9% of the total area. In comparison, only 7.1% of grasslands were adversely affected by CC, mainly situated in southeastern Tibet, Northwestern Xinjiang, Northwestern Sichuan, Northeastern Yunnan (Figure 8).
The HNPP trend suggested that HA had adversely affected grassland growth (SHNPP > 0), accounting for 62.0% of the terrestrial area. In the remaining 38.0% of grasslands, anthropogenic activities positively affected grassland growth (SHNPP < 0) (Figure 9).

3.2.2. The Relative Contributions of CC and HA to ANPP

Using the calculation results for the ANPP, PNPP, and HNPP change trends, the contribution of CC and HA to ANPP changes were assessed, and the areas that were largely influenced by CC were compared to those primarily affected by HA. Of the 38.4% of the ANPP area affected by CC and HA, climate-dominated areas accounted for 47.8%, whereas human-dominated areas accounted for 13.8% (Figure 10a).
As shown in Figure 10b, in most of China’s grasslands, the influence of CC on ANPP was greater than HA. The few regions where the contributions of HA were greater than that of CC were mainly distributed in Southeast Tibet, Central and Northeastern Qinghai, Western Xinjiang, and Central Inner Mongolia. Generally, the contributions of CC to the changes of ANPP were 71.0%, whereas the effect of HA was 29.0%.

3.3. Relationships between ANPP and Climate Factors

3.3.1. Multiple Correlations

According to previous research, temperature and precipitation are the dominant climate factors influencing ANPP changes [50,51,52]. To further explore its internal influence mechanism, a multiple linear regression was applied to estimate the correlation between ANPP and annual average temperature and annual average total precipitation. The results showed that ANPP was positively correlated with temperature and precipitation. Those areas with a significant positive correlation accounted for 46.5% of the total area. In comparison, the remaining 53.5% that did not exhibit a significant positive correlation were mainly distributed in Northwestern Xinjiang, Northwestern Sichuan, Southeastern Tibet, Southeastern Qinghai, and Northern Shaanxi (Figure 11). We also further estimated the correlation between ANPP and the average temperature and total precipitation in the growing season (Figure 12), and the spatial distribution of correlation coefficient was basically consistent with Figure 11. Those areas with a significant positive correlation accounted for 40.5% of the total area. In comparison, the remaining 59.5% did not exhibit a significant positive correlation (Figure 12).

3.3.2. Partial Correlations

To spatially determine the relationship of ANPP with annual average temperature and annual average total precipitation, we mapped their partial correlation results (Figure 13 and Figure 14). Precipitation was positively correlated with ANPP for 67.0% of the grassland areas. Areas with significantly positive relationships were mainly distributed in Ningxia, Southeast Qinghai, Inner Mongolia, and Northwest Xinjiang, accounting for 30.5% (Figure 13). Temperature was positively correlated with ANPP in 73.1% of the grasslands, and regions with significant positive correlations were mainly distributed in Qinghai, Tibet, Central Gansu, and Central Ningxia, accounting for 26.9% (Figure 14). We also further estimated the partial correlation between ANPP and the average temperature and total precipitation in the growing season (Figure 14 and Figure 15). Precipitation was positively correlated with ANPP for 66.2% of the grassland areas. Areas with significantly positive relationships were mainly distributed in Ningxia, Southeastern Qinghai, Inner Mongolia, and Northwest Xinjiang, accounting for 28.3% (Figure 15). Temperature was positively correlated with ANPP in 62.0% of the grasslands, and regions with significant positive correlations were mainly distributed in Qinghai, Tibet, Central Gansu, and Central Ningxia, accounting for 38.0% (Figure 16).

4. Discussion

4.1. Effects of CC on ANPP

Climate change is a crucial determinant in vegetation dynamics [53,54,55]. Numerous studies have identified precipitation and temperature as the primary climate factors influencing ANPP changes [56,57]. Based on the multiple and partial correlation analyses between ANPP and the main climate factors, we found that the ANPP of China’s grasslands generally showed an increasing trend over the last 20 years. Precipitation was found to be conducive to grassland growth, particularly in the arid and semi-arid regions, such as Gansu, Inner Mongolia, and Xinjiang (Figure 13 and Figure 15). In contrast, our results showed that temperature hindered grassland growth in these areas (Figure 14 and Figure 16), consistent with the conclusions of previous studies [1]. Since water is the most dominant factor affecting grassland productivity in arid and semi-arid areas, temperature increase may lead to the intensification of evaporation and water reduction, leading to increased grassland degradation [31,51].
This study assessed the CC and HA influences on changes in grassland ANPP. Overall, CC positively affected ANPP in most of China’s grasslands and had a greater effect than HA (Figure 10). This finding is consistent with previous studies [1]. Our results also indicated that CC was not conducive to grassland ANPP changes in some regions, particularly in southeastern Tibet (Figure 10). This was mainly attributed to the decrease in precipitation and the increase in temperature in the area during this period, considerably hindering grassland growth [31].
In addition, we found the contribution of CC was less in grassland restoration areas, such as Central Qinghai, Central Inner Mongolia, Central Tibet and Western Xinjiang (Figure 10). Among them, the contribution of CC in Central Qinghai and Central Inner Mongolia was small, mainly because they were located in population concentration areas, and the contribution of human activities was large, and had a positive impact on grassland restoration, which was consistent with [31,32]. However, the contribution of CC in Central Tibet and other places was small, mainly due to the obvious weakening of grazing activities in the past 20 years [3]. The contribution of CC in western Xinjiang and other places was small, and inconsistent with [31,32], which might have been due to inconsistent methods.

4.2. Effects of HA on ANPP

Anthropogenic activities are also critical to grassland NPP changes [5,52,58]. Human actions, for instance, overgrazing, land use change, and excessive exploitation, can lead to significant grassland degradation [9,59]. Our results showed that HA was the main cause of degradation in China’s grasslands (Figure 10). The degraded areas were mainly located in Southeastern Tibet, Northwestern Xinjiang, and Southeastern Qinghai (Figure 10), and this was mainly due to the high grazing intensity in these areas [10,31]. This result is consistent with [31,36].
In most grassland areas in China, HA adversely affected grassland growth (SHNPP > 0), accounting for 62.0% of the terrestrial area (Figure 9). The statistically significant increasing trend accounted for 15.1% of grasslands, mainly distributed in Central Qinghai, Southwestern Xin-jiang, and Southwestern Tibet. Among them, HA in Southwestern Xinjiang and Southwestern Tibet had a great impact, mainly due to the increased grazing activities [10,31]. HA in central Qinghai had a great impact, mainly because it was located in a population concentration area.

4.3. Methods for Quantitative Assessment of the CC and HA Influence on Grassland ANPP

The quantitative and comprehensive evaluation of climate change effects and anthropogenic influences on the NPP of grasslands remains a major challenge [60,61]. Previous studies have largely analyzed the effects of CC or HA on grassland NPP separately [62,63,64,65], with only a few studies conducting comprehensive assessments. The most widely used quantitative methods (i.e., regression analysis, residual trend, biophysical model, and environmental background-based methods) all contain considerable levels of uncertainty [1,24,66]. To improve the accuracy of the grassland NPP assessment, we introduced modifications to the environmental background-based method to calculate the contributions of CC and HA to NPP.
In our method, grassland NPP in nature reserves is used as the environmental background value and was assumed to be affected only by CC and free from HA. The impacts of anthropogenic activities on vegetation NPP were then characterized by calculating the difference between the environmental background value and the NPP outside the reserve. Our results were largely consistent with [31,34], suggesting the reliability of the developed approach.
Our proposed method is superior to the regression analysis and the residual trend methods since it improves the differentiation of climatic and anthropogenic effects on grassland NPP. Our method is comparatively better than the biophysical model-based approach because it can avoid calculation uncertainties from the various model parameters. Moreover, our method provides better results compared to the traditional environmental background-based approach because the impact of eco-geographical zoning on grassland NPP can be removed. Therefore, this research could provide support for government decisions and methodological reference for the quantitative evaluation of grasslands and other ecosystems.
However, there are also some limitations to this study. First, we assumed that the grassland PNPP pixels outside the core area could be replaced by the pixel PNPP in the nature reserve with the most similar habitat in the same natural ecological geographical division. Due to the availability of data, we only judged the critical factors, such as precipitation, temperature, slope, aspect and elevation, as similar habitat conditions, ignoring the impact of biotic factors (e.g., introduced species, succession) [67]. There were certain errors in some areas where natural biotic factors interfered greatly. Second, although the method constructed in this research accurately identified the contributions of CC and HA to grassland NPP, it could not further separate the impact of different types of human activities (e.g., overgrazing or environmental pollution) on net primary productivity. Therefore, we will focus on solving the above problems in future research.

5. Conclusions

In this study, an improved environmental background-based method was developed to differentiate the influence of CC and HA on grassland NPP. Our results showed that the average grassland ANPP in China from 2001 to 2020 showed a significant increasing trend, with 89.1% of grasslands showing a rising trend, while 10.9% exhibited a decreasing trend. The findings also suggested that CC had a more considerable effect on grassland NPP changes than HA in China’s grasslands. CC contributed 71.0% to ANPP change, whereas HA contributed 29.0% to ANPP change. Precipitation was found to be conducive to grassland growth among arid and semi-arid regions, whereas temperature inhibited grassland growth in these areas. While HA was the main cause of grassland degradation, its effects on grasslands have decreased in most of China’s grasslands.

Author Contributions

Conceptualization, F.Y. and X.Z.; methodology, F.Y., X.-C.W., X.Z. and Y.Z.; software, B.P. and X.Z.; validation, B.P.; analysis, X.Z.; investigation, X.Z.; writing—original draft preparation, X.Z.; writing—review and editing, F.Y. and X.-C.W.; funding acquisition, F.Y. and X.-C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (72074156), the Major Project of the National Social Science Foundation of China (21&ZD166), Fundamental Research Funds for the Central Universities (310421102), and the State Key Laboratory of Earth Surface Processes and Resource Ecology project (370100014).

Data Availability Statement

The data presented in this study can be available on request from the corresponding author.

Acknowledgments

The authors would like to thank USGS and Resource and Environment Science and Data Center for providing the research data, and also thank Yu and Wang for patient guidance. The authors are also grateful to the editor and anonymous reviewers for their positive comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ge, W.; Deng, L.; Wang, F.; Han, J. Quantifying the contributions of human activities and climate change to vegetation net primary productivity dynamics in China from 2001 to 2016. Sci. Total Environ. 2021, 773, 145648. [Google Scholar] [CrossRef]
  2. Field, C.B.; Behrenfeld, M.J.; Randerson, J.T.; Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 1998, 281, 237–240. [Google Scholar] [CrossRef] [PubMed]
  3. Zhai, D.; Gao, X.; Li, B.; Yuan, Y.; Jiang, Y.; Liu, Y.; Li, Y.; Li, R.; Liu, W.; Xu, J. Driving climatic factors at critical plant developmental stages for Qinghai-Tibet plateau alpine grassland productivity. Remote Sens. 2022, 14, 1564. [Google Scholar] [CrossRef]
  4. Wu, C.; Hou, X.; Peng, D.; Gonsamo, A.; Xu, S. Land surface phenology of China’s temperate ecosystems over 1999–2013: Spatial–temporal patterns, interaction effects, covariation with climate and implications for productivity. Agric. For. Meteorol. 2016, 216, 177–187. [Google Scholar] [CrossRef]
  5. Zhao, W.; Wu, J.; Shen, Q.; Liu, L.; Lin, J.; Yang, J. Estimation of the net primary productivity of winter wheat based on the near-infrared radiance of vegetation. Sci. Total Environ. 2022, 838, 156090. [Google Scholar] [CrossRef] [PubMed]
  6. Potter, C.; Klooster, S.; Genovese, V. Net primary production of terrestrial ecosystems from 2000 to 2009. Clim. Chang. 2012, 115, 365–378. [Google Scholar] [CrossRef]
  7. Yu, K.; Ciais, P.; Seneviratne, S.I.; Liu, Z.; Chen, H.Y.H.; Barichivich, J.; Allen, C.D.; Yang, H.; Huang, Y.; Ballantyne, A.P. Field-based tree mortality constraint reduces estimates of model-projected forest carbon sinks. Nat. Commun. 2022, 13, 2094. [Google Scholar] [CrossRef]
  8. Restrepo, H.I.; Montes, C.R.; Bullock, B.P.; Mei, B. The effect of climate variability factors on potential net primary productivity uncertainty: An analysis with a stochastic spatial 3-PG model. Agric. For. Meteorol. 2022, 315, 108812. [Google Scholar] [CrossRef]
  9. Xie, C.; Wu, S.; Zhuang, Q.; Zhang, Z.; Hou, G.; Luo, G.; Hu, Z. Where anthropogenic activity occurs, anthropogenic activity dominates vegetation net primary productivity change. Remote Sens. 2022, 14, 1092. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Wang, Q.; Wang, Z.; Yang, Y.; Li, J. Impact of human activities and climate change on the grassland dynamics under different regime policies in the Mongolian Plateau. Sci. Total Environ. 2020, 698, 134304. [Google Scholar] [CrossRef]
  11. Hidalgo-Galvez, M.D.; Barkaoui, K.; Volaire, F.; Matias, L.; Cambrolle, J.; Fernandez-Rebollo, P.; Carbonero, M.D.; Perez-Ramos, I.M. Can trees buffer the impact of climate change on pasture production and digestibility of Mediterranean dehesas? Sci. Total Environ. 2022, 835, 155535. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, G.; Hautier, Y.; Zhang, Z.; Lu, X.; Han, X. Decoupled responses of above- and below-ground stability of productivity to nitrogen addition at the local and larger spatial scale. Glob. Chang. Biol. 2022, 28, 2711–2720. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, M.; Running, S.W. Drought-Induced Reduction in Global Terrestrial Net Primary Production from 2000 through 2009. Science 2010, 329, 940–943. [Google Scholar] [CrossRef] [PubMed]
  14. Qu, S.; Wang, L.; Lin, A.; Yu, D.; Yuan, M.; Li, C. Distinguishing the impacts of climate change and anthropogenic factors on vegetation dynamics in the Yangtze River Basin, China. Ecol. Indic. 2020, 108, 105724. [Google Scholar] [CrossRef]
  15. Du, X.; Jin, X.; Yang, X.; Yang, X.; Zhou, Y. Spatial pattern of land use change and its driving force in Jiangsu province. Int. J. Environ. Res. Public Health 2014, 11, 3215–3232. [Google Scholar] [CrossRef]
  16. Del Grosso, S. Global potential net primary production predicted from vegetation class, precipitation, and temperature. Ecology 2008, 89, 2971. [Google Scholar] [CrossRef]
  17. Liu, G.; Shao, Q.; Fan, J.; Ning, J.; Rong, K.; Huang, H.; Liu, S.; Zhang, X.; Niu, L.; Liu, J. Change Trend and Restoration Potential of Vegetation Net Primary Productivity in China over the Past 20 Years. Remote Sens. 2022, 14, 1634. [Google Scholar] [CrossRef]
  18. Ouyang, L.; Ma, H.; Wang, Z.; Wang, Z.; Yu, X. Impacts of climate change and human activities on net primary productivity of grassland in the eastern Inner Mongolia. Acta Ecol. Sin. 2020, 40, 6912–6924. [Google Scholar]
  19. Mueller, T.; Dressler, G.; Tucker, C.J.; Pinzon, J.E.; Leimgruber, P.; Dubayah, R.O.; Hurtt, G.C.; Boehning-Gaese, K.; Fagan, W.F. Human Land-Use practices lead to global Long-Term increases in photosynthetic capacity. Remote Sens. 2014, 6, 5717–5731. [Google Scholar] [CrossRef]
  20. Yin, L.; Dai, E.; Zheng, D.; Wang, Y.; Ma, L.; Tong, M. What drives the vegetation dynamics in the Hengduan Mountain region, southwest China: Climate change or human activity? Ecol. Indic. 2020, 112, 106013. [Google Scholar] [CrossRef]
  21. Qi, X.; Jia, J.; Liu, H.; Lin, Z. Relative importance of climate change and human activities for vegetation changes on China’s silk road economic belt over multiple timescales. Catena 2019, 180, 224–237. [Google Scholar] [CrossRef]
  22. Li, X.; Luo, Y.; Wu, J. Decoupling Relationship between Urbanization and Carbon Sequestration in the Pearl River Delta from 2000 to 2020. Remote Sens. 2022, 14, 526. [Google Scholar] [CrossRef]
  23. Shi, Y.; Jin, N.; Ma, X.; Wu, B.; He, Q.; Yue, C.; Yu, Q. Attribution of climate and human activities to vegetation change in China using machine learning techniques. Agric. For. Meteorol. 2020, 294, 108146. [Google Scholar] [CrossRef]
  24. Teng, M.; Zeng, L.; Hu, W.; Wang, P.; Yan, Z.; He, W.; Zhang, Y.; Huang, Z.; Xiao, W. The impacts of climate changes and human activities on net primary productivity vary across an ecotone zone in Northwest China. Sci. Total Environ. 2020, 714, 136691. [Google Scholar] [CrossRef] [PubMed]
  25. Li, C.; Wang, Y.; Wu, X.; Cao, H.; Li, W.; Wu, T. Reducing human activity promotes environmental restoration in arid and semi-arid regions: A case study in Northwest China. Sci. Total Environ. 2021, 768, 144525. [Google Scholar] [CrossRef]
  26. Li, H.; Zhang, H.; Li, Q.; Zhao, J.; Guo, X.; Ying, H.; Deng, G.; Rihan, W.; Wang, S. Vegetation productivity dynamics in response to climate change and human activities under different topography and land cover in northeast china. Remote Sens. 2021, 13, 975. [Google Scholar] [CrossRef]
  27. Liu, H.; Hou, L.; Kang, N.; Nan, Z.; Huang, J. A meta-regression analysis of the economic value of grassland ecosystem services in China. Ecol. Indic. 2022, 138, 108793. [Google Scholar] [CrossRef]
  28. Huang, L.; Ning, J.; Zhu, P.; Zheng, Y.; Zhai, J. The conservation patterns of grassland ecosystem in response to the forage-livestock balance in North China. J. Geogr. Sci. 2021, 31, 518–534. [Google Scholar] [CrossRef]
  29. Yu, G.; Li, X.; Wang, Q.; Li, S. Carbon storage and its spatial pattern of terrestrial ecosystem in china. J. Resour. Ecol. 2010, 1, 97–109. [Google Scholar]
  30. Xu, L.; Yu, G.; He, N. Increased soil organic carbon storage in Chinese terrestrial ecosystems from the 1980s to the 2010s. J. Geogr. Sci. 2019, 29, 49–66. [Google Scholar] [CrossRef]
  31. Yan, Y.; Liu, X.; Wen, Y.; Ou, J. Quantitative analysis of the contributions of climatic and human factors to grassland productivity in northern China. Ecol. Indic. 2019, 103, 542–553. [Google Scholar] [CrossRef]
  32. Ge, J.; Hou, M.; Liang, T.; Feng, Q.; Meng, X.; Liu, J.; Bao, X.; Gao, H. Spatiotemporal dynamics of grassland aboveground biomass and its driving factors in North China over the past 20 years. Sci. Total Environ. 2022, 826, 154226. [Google Scholar] [CrossRef] [PubMed]
  33. Song, Z.; Wu, Y.; Yang, Y.; Zhang, X.; Van Zwieten, L.; Bolan, N.; Li, Z.; Liu, H.; Hao, Q.; Yu, C.; et al. High potential of stable carbon sequestration in phytoliths of China’s grasslands. Glob. Chang. Biol. 2022, 28, 2736–2750. [Google Scholar] [CrossRef] [PubMed]
  34. Zhou, W.; Gang, C.; Zhou, L.; Chen, Y.; Li, J.; Ju, W.; Odeh, I. Dynamic of grassland vegetation degradation and its quantitative assessment in the northwest China. Acta Oecol. 2014, 55, 86–96. [Google Scholar] [CrossRef]
  35. Guo, Y.; Li, R.; Yang, Y.; Ma, J.; Zheng, H. Integrating future grassland degradation risk to improve the spatial targeting efficiency of payment for ecosystem services. J. Environ. Manag. 2022, 317, 115490. [Google Scholar] [CrossRef]
  36. Zhou, W.; Yang, H.; Huang, L.; Chen, C.; Lin, X.; Hu, Z.; Li, J. Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 2017, 83, 303–313. [Google Scholar] [CrossRef]
  37. Naeem, S.; Zhang, Y.; Tian, J.; Qamer, F.M.; Latif, A.; Paul, P.K. Quantifying the Impacts of Anthropogenic Activities and Climate Variations on Vegetation Productivity Changes in China from 1985 to 2015. Remote Sens. 2020, 12, 1113. [Google Scholar] [CrossRef]
  38. Wang, X.; Pan, S.; Pan, N.; Pan, P. Grassland productivity response to droughts in northern China monitored by satellite-based solar-induced chlorophyll fluorescence. Sci. Total Environ. 2022, 830, 154550. [Google Scholar] [CrossRef]
  39. Zhang, Y.; Du, R.; Chen, H.; Du, P.; Zhang, S.; Ren, W. Different characteristics of microbial diversity and special functional microbes in rainwater and topsoil before and after 2019 new coronavirus epidemic in Inner Mongolia Grassland. Sci. Total Environ. 2022, 809, 151088. [Google Scholar] [CrossRef]
  40. Lu, L.; Sun, Z.; Qimuge, E.; Ye, H.; Huang, W.; Nie, C.; Wang, K.; Zhou, Y. Using remote sensing data and Species-Environmental matching model to predict the potential distribution of grassland rodents in the northern china. Remote Sens. 2022, 14, 2168. [Google Scholar] [CrossRef]
  41. Jiang, H.; Xu, X.; Guan, M.; Wang, L.; Huang, Y.; Liu, Y. Simulation of spatiotemporal land use changes for integrated model of socioeconomic and ecological processes in china. Sustainability 2019, 11, 3627. [Google Scholar] [CrossRef]
  42. Wu, X.; Dang, X.; Meng, Z.; Fu, D.; Cong, W.; Zhao, F.; Guo, J. Mechanisms of grazing management impact on preferential water flow an infiltration patterns in a semi-arid grassland in northern China. Sci. Total Environ. 2022, 813, 152082. [Google Scholar] [CrossRef] [PubMed]
  43. Xu, J.; Sun, G.; Liu, Y. Diversity and complexity in the forms and functions of protected areas in china. J. Int. Wildl. Law Policy 2014, 17, 102–114. [Google Scholar] [CrossRef]
  44. Wandesforde-Smith, G.; Snyder, K.D.; Hart, L.A. Biodiversity conservation and protected areas in china: Science, law, and the obdurate Party-State. Soc. Sci. Electron. Publ. 2014, 17, 1476–1548. [Google Scholar] [CrossRef]
  45. Tobler, W. On the first law of geography: A reply. Ann. Assoc. Am. Geogr. 2004, 94, 304–310. [Google Scholar] [CrossRef]
  46. Wei, X.; Yang, J.; Luo, P.; Lin, L.; Lin, K.; Guan, J. Assessment of the variation and influencing factors of vegetation NPP and carbon sink capacity under different natural conditions. Ecol. Indic. 2022, 138, 108834. [Google Scholar] [CrossRef]
  47. Jiang, H.; Xu, X.; Guan, M.; Wang, L.; Huang, Y.; Jiang, Y. Determining the contributions of climate change and human activities to vegetation dynamics in agro-pastural transitional zone of northern China from 2000 to 2015. Sci. Total Environ. 2020, 718, 134871. [Google Scholar] [CrossRef]
  48. Wang, Y.; Zhang, Z.; Chen, X. The dominant driving force of forest change in the yangtze river basin, china: Climate variation or anthropogenic activities? Forests 2022, 13, 82. [Google Scholar] [CrossRef]
  49. Ma, R.; Xia, C.; Liu, Y.; Wang, Y.; Zhang, J.; Shen, X.; Lu, X.; Jiang, M. Spatiotemporal change of net primary productivity and its response to climate change in temperate grasslands of china. Front. Plant Sci. 2022, 13, 899800. [Google Scholar] [CrossRef]
  50. Ma, J.; Weng, B.; Bi, W.; Xu, D.; Xu, T.; Yan, D. Impact of climate change on the growth of typical crops in karst areas: A case study of guizhou province. Adv. Meteorol. 2019, 2019, 16. [Google Scholar] [CrossRef]
  51. Hoover, D.L.; Hajek, O.L.; Smith, M.D.; Wilkins, K.; Slette, I.J.; Knapp, A.K. Compound hydroclimatic extremes in a semi-arid grassland: Drought, deluge, and the carbon cycle. Glob. Chang. Biol. 2022, 28, 2611–2621. [Google Scholar] [CrossRef] [PubMed]
  52. Pan, Y.; Jackson, R.B.; Hollinger, D.Y.; Phillips, O.L.; Nowak, R.S.; Norby, R.J.; Oren, R.; Reich, P.B.; Luescher, A.; Mueller, K.E.; et al. Contrasting responses of woody and grassland ecosystems to increased CO2 as water supply varies. Nat. Ecol. Evol. 2022, 6, 315. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, S.; Zhang, B.; Yang, Q.; Chen, G.; Yang, B.; Lu, L.; Shen, M.; Peng, Y. Responses of net primary productivity to phenological dynamics in the Tibetan Plateau, China. Agric. For. Meteorol. 2017, 232, 235–246. [Google Scholar] [CrossRef]
  54. Zhang, M.; Yuan, N.; Lin, H.; Liu, Y.; Zhang, H. Quantitative estimation of the factors impacting spatiotemporal variation in NPP in the Dongting Lake wetlands using Landsat time series data for the last two decades. Ecol. Indic. 2022, 135, 108544. [Google Scholar] [CrossRef]
  55. Ma, S.; Li, Y.; Zhang, Y.; Wang, L.; Jiang, J.; Zhang, J. Distinguishing the relative contributions of climate and land use/cover changes to ecosystem services from a geospatial perspective. Ecol. Indic. 2022, 136, 108645. [Google Scholar] [CrossRef]
  56. Chen, T.; Tang, G.; Yuan, Y.; Guo, H.; Xu, Z.; Jiang, G.; Chen, X. Unraveling the relative impacts of climate change and human activities on grassland productivity in Central Asia over last three decades. Sci. Total Environ. 2020, 743, 140649. [Google Scholar] [CrossRef]
  57. Wang, S.; Zhang, Q.; Yue, P.; Wang, J.; Yang, J.; Wang, W.; Zhang, H.; Ren, X. Precipitation-Use efficiency and its conversion with climate types in mainland china. Remote Sens. 2022, 14, 2467. [Google Scholar] [CrossRef]
  58. Zhou, Y.; Yue, D.; Guo, J.; Chao, Z.; Meng, X. Assessing the impact of land conversion and management measures on the net primary productivity in the Bailong River Basin, in China. Catena 2021, 207, 105672. [Google Scholar] [CrossRef]
  59. Liu, H.; Cao, L.; Jia, J.; Gong, H.; Qi, X.; Xu, X. Effects of land use changes on the nonlinear trends of net primary productivity in arid and semi-arid areas, China. Land Degrad. Dev. 2021, 32, 2183–2196. [Google Scholar] [CrossRef]
  60. Yan, W.; Wang, Y.; Chaudhary, P.; Ju, P.; Zhu, Q.; Kang, X.; Chen, H.; He, Y. Effects of climate change and human activities on net primary production of wetlands on the Zoige Plateau from 1990 to 2015. Glob. Ecol. Conserv. 2022, 35, e02052. [Google Scholar] [CrossRef]
  61. Xiong, Q.; Xiao, Y.; Liang, P.; Li, L.; Zhang, L.; Li, T.; Pan, K.; Liu, C. Trends in climate change and human interventions indicate grassland productivity on the Qinghai-Tibetan Plateau from 1980 to 2015. Ecol. Indic. 2021, 129, 108010. [Google Scholar] [CrossRef]
  62. Du, J.; Fu, Q.; Fang, S.; Wu, J.; He, P.; Quan, Z. Effects of rapid urbanization on vegetation cover in the metropolises of China over the last four decades. Ecol. Indic. 2019, 107, 105458. [Google Scholar] [CrossRef]
  63. Zhao, L.; Liu, Z.; Hu, Y.; Zhou, W.; Peng, Y.; Ma, T.; Liu, L.; Li, S.; Wang, L.; Mao, X. Evaluation of reasonable stocking rate based on the relative contribution of climate change and grazing activities to the productivity of alpine grasslands in Qinghai province. Remote Sens. 2022, 14, 1455. [Google Scholar] [CrossRef]
  64. Lei, T.; Wu, J.; Li, X.; Geng, G.; Shao, C.; Zhou, H.; Wang, Q.; Liu, L. A new framework for evaluating the impacts of drought on net primary productivity of grassland. Sci. Total Environ. 2015, 536, 161–172. [Google Scholar] [CrossRef] [PubMed]
  65. Gang, C.; Zhao, W.; Zhao, T.; Zhang, Y.; Gao, X.; Wen, Z. The impacts of land conversion and management measures on the grassland net primary productivity over the Loess Plateau, Northern China. Sci. Total Environ. 2018, 645, 827–836. [Google Scholar] [CrossRef] [PubMed]
  66. Zhou, Y.; Yue, D.; Li, C.; Mu, X.; Guo, J. Identifying the spatial drivers of net primary productivity: A case study in the Bailong River Basin, China. Glob. Ecol. Conserv. 2021, 28, e01685. [Google Scholar] [CrossRef]
  67. Pan, Y.; Yu, C.; Zhang, X.; Chen, B.; Wu, J.; Tu, Y.; Miao, Y.; Luo, L. A modified framework for the regional assessment of climate and human impacts on net primary productivity. Ecol. Indic. 2016, 60, 184–191. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Remotesensing 14 04844 g001
Figure 2. Eco-geographical zoning data in China.
Figure 2. Eco-geographical zoning data in China.
Remotesensing 14 04844 g002
Figure 3. The spatial distribution of national nature reserves in China.
Figure 3. The spatial distribution of national nature reserves in China.
Remotesensing 14 04844 g003
Figure 4. The spatial scope of grasslands in China.
Figure 4. The spatial scope of grasslands in China.
Remotesensing 14 04844 g004aRemotesensing 14 04844 g004b
Figure 5. Spatial distribution of annual mean ANPP (2001–2020).
Figure 5. Spatial distribution of annual mean ANPP (2001–2020).
Remotesensing 14 04844 g005
Figure 6. Changes in the annual mean ANPP (2001–2020).
Figure 6. Changes in the annual mean ANPP (2001–2020).
Remotesensing 14 04844 g006
Figure 7. (a) Trend in ANPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Figure 7. (a) Trend in ANPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Remotesensing 14 04844 g007
Figure 8. (a) Trend in PNPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Figure 8. (a) Trend in PNPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Remotesensing 14 04844 g008
Figure 9. (a) Trend in HNPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Figure 9. (a) Trend in HNPP (2001–2020) and (b) the significance levels. (Note: K is the change trend of ANPP, and P is the significance levels).
Remotesensing 14 04844 g009
Figure 10. Spatial patterns of (a) driving factors and (b) contributions to ANPP changes.
Figure 10. Spatial patterns of (a) driving factors and (b) contributions to ANPP changes.
Remotesensing 14 04844 g010
Figure 11. Spatial distribution of (a) multiple correlation coefficient between ANPP and annual average temperature and annual average total precipitation and (b) the significance level.
Figure 11. Spatial distribution of (a) multiple correlation coefficient between ANPP and annual average temperature and annual average total precipitation and (b) the significance level.
Remotesensing 14 04844 g011
Figure 12. Spatial distribution of (a) multiple correlation coefficient between ANPP and the average temperature and total precipitation in the growing season and (b) the significance level.
Figure 12. Spatial distribution of (a) multiple correlation coefficient between ANPP and the average temperature and total precipitation in the growing season and (b) the significance level.
Remotesensing 14 04844 g012
Figure 13. Spatial patterns of (a) partial correlation coefficients between ANPP and annual average total precipitation and (b) the corresponding significance levels.
Figure 13. Spatial patterns of (a) partial correlation coefficients between ANPP and annual average total precipitation and (b) the corresponding significance levels.
Remotesensing 14 04844 g013
Figure 14. Spatial patterns of (a) partial correlation coefficients between ANPP and annual average temperature and (b) the corresponding significance levels.
Figure 14. Spatial patterns of (a) partial correlation coefficients between ANPP and annual average temperature and (b) the corresponding significance levels.
Remotesensing 14 04844 g014
Figure 15. Spatial patterns of (a) partial correlation coefficients between ANPP and total precipitation in the growing season and (b) the corresponding significance levels.
Figure 15. Spatial patterns of (a) partial correlation coefficients between ANPP and total precipitation in the growing season and (b) the corresponding significance levels.
Remotesensing 14 04844 g015
Figure 16. Spatial patterns of (a) partial correlation coefficients between ANPP and the average temperature in the growing season and (b) the corresponding significance levels.
Figure 16. Spatial patterns of (a) partial correlation coefficients between ANPP and the average temperature in the growing season and (b) the corresponding significance levels.
Remotesensing 14 04844 g016
Table 1. Data sources.
Table 1. Data sources.
DataWebsiteData TypeResolution
LULC (CCI land cover)https://www.esa-landcover-cci.org (4 January 2022)raster300 m
NPPhttps://lpdaacsvc.cr.usgs.gov (10 January 2022)raster1 km
Precipitation and temperaturehttps://data.cma.cn/ (16 January 2022)vector
Eco-geographical zoninghttps://www.resdc.cn/ (25 January 2022)vector
DEMhttps://www.resdc.cn/ (27 January 2022)raster1 km
National natural reserveshttps://www.resdc.cn/ (28 January 2022)vector
Table 2. Defined scenarios.
Table 2. Defined scenarios.
SchemeDriving FactorsContribution
Climate (%)Human (%)
SANPP > 0SPNPP > 0, SHNPP < 0Both 100   ×   | S P N P P | | S P N P P |   +   | S H N P P | 100   ×   | S H N P P | | S P N P P |   +   | S H N P P |
SPNPP < 0, SHNPP < 0Human activities0100
SPNPP > 0, SHNPP > 0Climate change1000
SANPP < 0SPNPP < 0, SHNPP > 0Both 100   ×   | S P N P P | | S P N P P |   +   | S H N P P | 100   ×   | S H N P P | | S P N P P |   +   | S H N P P |
SPNPP > 0, SHNPP > 0Human activities0100
SPNPP < 0, SHNPP < 0Climate change1000
Note: SANPP, SPNPP, and SHNPP are the slopes of ANPP, PNPP and HNPP, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, X.; Peng, B.; Zhou, Y.; Yu, F.; Wang, X.-C. Quantifying the Influence of Climate Change and Anthropogenic Activities on the Net Primary Productivity of China’s Grasslands. Remote Sens. 2022, 14, 4844. https://doi.org/10.3390/rs14194844

AMA Style

Zhou X, Peng B, Zhou Y, Yu F, Wang X-C. Quantifying the Influence of Climate Change and Anthropogenic Activities on the Net Primary Productivity of China’s Grasslands. Remote Sensing. 2022; 14(19):4844. https://doi.org/10.3390/rs14194844

Chicago/Turabian Style

Zhou, Xiafei, Binbin Peng, Ying Zhou, Fang Yu, and Xue-Chao Wang. 2022. "Quantifying the Influence of Climate Change and Anthropogenic Activities on the Net Primary Productivity of China’s Grasslands" Remote Sensing 14, no. 19: 4844. https://doi.org/10.3390/rs14194844

APA Style

Zhou, X., Peng, B., Zhou, Y., Yu, F., & Wang, X. -C. (2022). Quantifying the Influence of Climate Change and Anthropogenic Activities on the Net Primary Productivity of China’s Grasslands. Remote Sensing, 14(19), 4844. https://doi.org/10.3390/rs14194844

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop