A Method for Quantifying the Impacts of Human Activities on Net Primary Production of Grasslands in Northwest China

Accurately assessing the impact of human activities on net primary productivity (NPP) of vegetation is of great significance to the achievement of sustainable development. However, it is difficult to disentangle the effects of climate conditions and human activities on NPP, and bridging this knowledge gap largely depends on the calculation of the NPP under natural conditions. Here, we propose a method for calculating natural vegetation NPP (NNPP) based on non-human influence grids, which are obtained according to the consistent rate of climate and actual NPP (ANPP) temporal changes. We selected Northwest China as study area, and we used a light use efficiency (LUE) model to estimate ANPP and used the random forest algorithm (RF) to estimate the NNPP. The results show that NNPP is very close to ANPP, and the human activities on NPP (HNPP) based on NNPP is close to the actual situation of human activities on NPP. From 2001 to 2017, the positive HNPP accounts for 40.28% of the total grassland area, with an average value of 28.65 gC·m−2·yr−1, while the negative HNPP accounts for 59.72% of the total area, with an average value of −31.19 gC·m−2·yr−1. The grassland NPP shows an increasing trend, which is dominated by climate factors. Human activity is the dominant factor for the grassland degradation, accounting for 42.78% of the degraded area, but promoting grassland growth in 11.4% of the restored area. This study provides a new method to estimate the impacts of human activities on vegetation, and the results can be used to evaluate the effectiveness of ecological environmental governance, providing a quantitative basis for scientifically building the harmonious relationship between human and nature.


Introduction
Grassland is an important part of terrestrial ecosystem, accounting for about 25% of the total land area of the earth [1]. The grassland vegetation carbon pool accounts for 20% of the world's total [2]. Net primary productivity (NPP) of vegetation is the residual assimilation of autotrophic respiration of all organic matter produced by green plants through photosynthesis [3]. It directly reflects the quality of terrestrial ecosystems, and it is the main factor determining the role of carbon source and sink of the ecosystem [4]. Therefore, quantifying the effects of climate change and human activities on grassland is essential to develop adaptation strategies to cope with climate change and the challenges posed by human activities to ecosystems.
Human activity is an important factor affecting NPP changes [5]. At present, there are several methods to quantitatively calculate the impact of human activities on NPP (HNPP).
Human activity is an important factor affecting NPP changes [5]. At present, there are several methods to quantitatively calculate the impact of human activities on NPP (HNPP). Based on the land-use change, the NPP can be estimated [6], while this method cannot estimate the anthropogenic effects in areas where the land use unchanged. The climate-sensitive calculation method [7] obtains weights by fitting the sensitivity of NPP to meteorological factor to obtain the NPP anthropogenic influenced value, and this method requires the sensitivity of the NPP includes impacts of human activities. The most widely used method is the difference between actual NPP (ANPP) and potential NPP (PNPP) [8]. The ANPP can be estimated using the improved Carnegie-Ames-Stanford approach (CASA) model [9]. The CASA model calculates NPP based on vegetation physiological process, photosynthetic effective radiation absorbed by vegetation and light energy utilization rate. The model can be used to model long-term and large-scale NPP values [10]. The PNPP is the highest form of vegetation under field conditions, and it corresponds to the NPP of the potential vegetation, and thus PNPP is relatively high. PNPP is generally based on climatological models [11], such as Miami model [12] and Thornthwaite Memorial Model [13]. These models use meteorological factors such as temperature, precipitation, and potential evapotranspiration as input parameters, while soil, topography, and vegetation types were poorly considered [3], which strongly affect vegetation growth [14]. As a result, these methods usually lead to large errors in the simulated PNPP values [15]. Therefore, it is urgent to establish a reliable method to calculate the actual natural NPP under current climate conditions.
Theoretically, if the grids without the effects of human activities can be identified, the ANPP of these grids can be regarded as the natural vegetation NPP (NNPP). Based on this methodology, we selected the northwestern region of China as study area, which accounts for about 41.29% of the China's grassland [16]. This area is the representative of a fragile ecology, and it is vulnerable to human activities. The main goal of this research is to establish a new estimation method to calculate the HNPP, and the results can provide a new idea for quantifying impacts of human activities on vegetation NPP.

Study Area
The Northwest China ( Figure 1) (31°32′ N-49°10′ N and 73°15′ E-111°50′ E) covers an area of about 3.5 million km 2 , accounting for 36% of China's total land area. Located in deep inland, this area has a typical continental temperate climate, with a mean annual air temperature of 5.06 ℃. The precipitation is low and gradually decreases from east to west. The topography is dominated by mountains, plateaus, and basins, with an average elevation of 2267 m. The land cover types include forest, grassland, farmland, and desert.
The digital elevation model (DEM) data, population density data, and land use map were downloaded from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 23 January 2020). The vector data were obtained from the National Basic Geographic Information Center (http://www.ngcc.cn/ngcc/, accessed on 23 January 2020). The grazing data were obtained from the National Bureau of Statistics of China (http://www.stats.gov.cn/, accessed on 23 January 2020), and the livestock data includes the data of the Ningxia region from 2001 to 2017 and the Sanjiangyuan region from 2007 to 2017.
The measured NPP values were derived from field sampling (n = 27) and literatures (n = 71) [20,21]. During field sampling, the grass plot was 25 m × 25 m, and three plots of 30 cm × 30 cm were selected on the diagonal. The above-ground biomass was harvested, the roots were dug out, and the sample was dried (85 • C). After reaching a constant weight, the dry mass was weighed to obtain biomass, and the dry weight was multiplied the conversion coefficient by 0.475 to obtain the NPP [22].

Calculation of ANPP
The CASA model has been widely used to estimate ANPP. This method can estimate NPP over a large region and for long time series [23]. We used the CASA model to calculate ANPP of the grasslands as follows: where, APAR(x,t) represents the photosynthetically active radiation (MJ·m −2 ) absorbed by pixel x in month t; and ε(x,t) is the actual light-use efficiency (gC·MJ −1 ) of pixel x in month t. APAR(x,t) and ε(x,t) are calculated using Equations (2) and (3), respectively.
where, SOL(x,t) represents the total solar radiation (MJ·m −2 ) of pixel x in month t; FPAR(x,t) represents the absorption coefficient of the vegetation layer on incident photosynthetically active radiation, which is determined by the NDVI; 0.5 is the proportion of SOL intercepted by the vegetation.
where, T ε1 (x,t) and T ε2 (x,t) are the temperature stress coefficients of low and high temperature to light-use efficiency, respectively; W p (x,t) is the coefficient of water stress, reflecting the influence of the water conditions. The ε max is the maximum light-use efficiency under ideal conditions. The maximum light-use efficiency rate refers to the simulation value of the maximum light energy utilization rate of typical vegetation [24]. We used the NDVI-FPAR lookup table provided in the NASA-MOD15 algorithm to calculate the monthly FPAR. The algorithm is as follows: FPAR(x, t) = 0, NDV I ≤ 0.075 min(1.16 × NDV I(x, t) − 0.0439), NDV I > 0.075 (4) where, W p is the water stress coefficient, pre and pre max are the monthly precipitation and the maximum monthly precipitation, respectively. Wu [25] and Piao [26] suggested that there was a time lag between vegetation growth and climatic factors. Therefore, according to Luo's result [7], we set the lag of three months in the CASA model for precipitation, i.e., the precipitation pre i of this month is the average of the precipitation of this month and the precipitation of the previous three months. The algorithm was expressed as follows: LSWI is the surface water index, which is calculated as follows: where, RE f2 and RE f6 are the 2nd and 6th bands of MOD09A1 data, respectively.
where, T opt is the monthly average temperature ( • C) when the NDVI value of a certain region is the highest in a year. When the average temperature of a month is less than or equal to −10 • C, the value is 0 [9].
where, T(x,t) for the average temperature of the month; When the monthly average temperature T(x,t) in a month is 10 • C higher or 13 • C lower than T opt (x,t), the value of T ε2 (x,t) in this month is equal to half of that when the monthly average temperature T(x,t) is T opt (x,t) [9].

Calculation of NNPP
The climate factors used for estimating PNPP are precipitation and temperature in the climatological model. Therefore, the PNPP was regarded as a comprehensive indicator of multiple climatic factors, and changes in PNPP represent the climate change. Theoretically, if there were only effects of climate change, the changing patterns of vegetation NPP should be consistent with the climate, and the ANPP of the corresponding grid is NNPP. Based on this hypothesis, we extracted grid without human activities by determining whether the ANPP changing trends of the grassland vegetation were consistent with the climate change. The actual NPP with unmanned influence point is defined as NNPP, and the parameters of this grid are trained using the random forest learning algorithm, and finally the NNPP values at the grid scale were upscaled to the entire area ( Figure 2). Remote Sens. 2021, 13, x FOR PEER REVIEW 5 of 18 (1) Identify the grid without human influence Climate change is a comprehensive index and it is difficult to calculate with specific formulas. The climatological model for estimating PNPP only uses climate factors. Therefore, the PNPP estimated by the climatological model can be used as a climate index. The comparison between ANPP and climate is the comparison between ANPP and PNPP. PNPP was estimate by the Thornthwaite Memorial Model [27], and the calculation formula is as follows: where, PNPP unit is gC·m −2 ·yr −1 ; V is the average annual actual evaporation (mm); r is the annual cumulative precipitation (mm); L is the annual evaporation capacity (mm); t is the average annual temperature (°C). According to the above description, the changes of ANPP and PNPP are compared year by year. The formula is as follows: (1) Identify the grid without human influence Climate change is a comprehensive index and it is difficult to calculate with specific formulas. The climatological model for estimating PNPP only uses climate factors. Therefore, the PNPP estimated by the climatological model can be used as a climate index. The comparison between ANPP and climate is the comparison between ANPP and PNPP. PNPP was estimate by the Thornthwaite Memorial Model [27], and the calculation formula is as follows: L = 3000 + 25t + 0.05t 3 (12) where, PNPP unit is gC·m −2 ·yr −1 ; V is the average annual actual evaporation (mm); r is the annual cumulative precipitation (mm); L is the annual evaporation capacity (mm); t is the average annual temperature ( • C).
According to the above description, the changes of ANPP and PNPP are compared year by year. The formula is as follows: We defined Coincidence Ratio (CR) to examine the differences between climate change and ANPP change. If A i = T i , then CR i = 1, indicating that ANPP and PNPP showed same interannual variations, otherwise CR i = 0, indicating that they were different.
where, CR y is the Coincidence Ratio of year between climate change and ANPP change during the study period, i represents the year, CR i indicates whether the ith year coincides, and n is the total number of years in the study period. Due to the disturbance in the ecosystem, a threshold is set for the judgment. If it is greater than the threshold, it indicates that the ANPP is consistent with climate change, and this grid is an unmanned influence grid. The above judgment is based on the interannual scale and failed to detect the impact of human activities during the year (at the month scale). For example, grazing is a human activity, the grass will gradually recover, and the interannual coincidence rate may not be recognized. According to the growth pattern of grassland, the monthly NPP variation shows a single peak distribution, and its maximum value should appear in June and July. If grassland was affected by humans during the annual, its monthly NPP changes would not be consistent with this rule. Since it is impossible to determine whether the maximum NPP of each grid is in June or July, June and July were excluded in the calculation. The specific calculation method is as follows: If, ANPP ij − ANPP ij+1 ≥ 0, (j = 8, 9, 10, 11), A ij = 1, else A ij = 0 (17) where, ANPP ij is the ANPP of the j month of the ith year, i represents the year (2001~2017), j represents the month, R i represents whether the inter-monthly NPP change of the ith year follows a single peak distribution.
where, R m is the ratio of monthly NPP subject to single peak distribution during the study period, and n is the total number of years in the study period (n = 17). Theoretically, the consistency of ANPP and PNPP includes inter-annual change and intra-annual change, and the coincidence rate should be 100%. However, the 17-year period is a relatively long time, the CR s could be affected by many other factors such as extreme events and non-human disturbances. Therefore, it is necessary to manually set the thresholds of the CR y and R m values. To find the best values, we tried several values, and the results were assessed using field surveys of unmanned influenced points to validate whether there was human activity. Based on the field investigation results, we manually set the thresholds of CR y and R m as 0.88 and 0.94 respectively, that is, during the 17-year period, CR y and R m are allowed to not coincide for 2 years and 1 year respectively. The grids with the values between these thresholds, namely (CR y > 0.88) ∩ (R m > 0.94), were identified as the grids without the effects of human activities.
(2) Upscaling the NPP to regional scale This study used the random forest (RF) machine learning methods to upscale the NNPP to the entire study area. Random forest is a multifactor machine learning algorithm, which is suitable for solving classification and regression problems. Multiple random samples are obtained through multiple sampling, and corresponding decision trees are established through these samples, and the specific algorithm can be found [28]. RF has a high prediction accuracy, a good tolerance for outliers and noise, and is not prone to overfitting [29]. This method has been successfully used to model the net primary productivity [30], as well as forest biomass [31].
The random forest uses the ANPP at unmanned influence points as the dependent variable, and 10 factors, including elevation, slope, slope orientation, mean annual temperature, mean annual precipitation, mean annual solar radiation, NDVI, mean temperature in growing season, cumulative precipitation in growing season, and solar radiation in growing season, as independent variables. We defined May to October as growing season because this period is usually considered as the growing season in the Northwest China [32]. Training and parameter tuning are performed to obtain the characteristic relationship of each variable, and further simulate the grassland NNPP.

Calculation of HNPP
We compared the calculation results of HNPP based on NNPP and PNPP, which were named as HNPP_N and HNPP_P, respectively. The HNPP was calculated by the difference between ANPP and NNPP (or PNPP), which was expressed as follows: where, HNPP value > 0 suggests that the effect of human activities on the vegetation is positive; HNPP value < 0 indicates that the effect of human activities on the vegetation is negative.

Changing Trend
For each pixel, the interannual changing rate for HNPP was estimated by unitary linear regression, and the calculation formula is as follows: where n represents years (the time series was from 2001~2017, n = 17); HNPP i represents the HNPP in year i.
We use the F-test to test the significance of NPP changing trend, the formula is as follows: where, ∧ y is the regression value, y i is the value of the ith year, y is the multi-year average, and n is the number of years. The results are divided into 4 groups: significant increase (Slope > 0, p < 0.01), increase but not significant (Slope > 0, 0.01 ≤ p ≤ 0.05), and decrease but not significant (Slope < 0, 0.01 ≤ p ≤ 0.05), significant decrease (Slope < 0, p < 0.01).
In this study, the changing direction and trends of HNPP are defined based on the values of HNPP and Slope. When HNPP and slope values are positive, it indicates that the positive influence of human activities on NPP is increasing. Slope value is negative, it indicates that the positive influence is decreasing. When HNPP and Slope values are negative, it indicates that the negative influence of human activities on NPP is increasing. When the slope value is positive, it indicates that the negative influence is decreasing.

ANPP Estimation and Spatial and Temporal Distribution
The accuracy (Figure 3) of modelling actual NPP is verified using the measured data (98 sites, 71 from literatures, 27 from our field work). The correlation coefficient (R 2 = 0.63, p < 0.01) indicates that the estimate is reliable. The linear fitting scatter plot mainly distributed in low value area, which can be explained as that the most areas had low NNP values. positive influence of human activities on NPP is increasing. Slope value is negative, it indicates that the positive influence is decreasing. When HNPP and Slope values are negative, it indicates that the negative influence of human activities on NPP is increasing. When the slope value is positive, it indicates that the negative influence is decreasing.

ANPP Estimation and Spatial and Temporal Distribution
The accuracy (Figure 3) of modelling actual NPP is verified using the measured data (98 sites, 71 from literatures, 27 from our field work). The correlation coefficient (R 2 = 0.63, p < 0.01) indicates that the estimate is reliable. The linear fitting scatter plot mainly distributed in low value area, which can be explained as that the most areas had low NNP values.  Figure  5). The ANPP values were higher in southern area, which have higher precipitation in Northwest China. Other high value areas are mainly distributed in high mountainous areas. The low values are mainly located in the northern desert areas. The total area with an increasing trend was 87.08 × 10 4 km 2 , accounting for 75.77% of the total grassland area. Among them, 26.58% showed significant increase (Figure 6), mainly located in the eastern Qinghai, northern Shaanxi, and Gansu province. Moreover, nearly 24.23% of grassland indicated a decreasing trend in grassland ANPP, including only 0.60% as significant decrease, mainly scattered in most parts of northern Xinjiang.

Non-Human Affected Grids
The population density in the study area is shown as Figure 7, and the most of the unmanned influence points are distributed in areas with low population density. A total of 10,231 non-human affected grids were extracted (Figure 7). Unmanned influence points are mainly distributed in uninhabited areas such as plateaus, mountains, and deserts. Most of these places are far away from strong human activities areas such as cities and towns, located in nature reserves and ANPP low-value areas.

Non-Human Affected Grids
The population density in the study area is shown as Figure 7, and the most of the unmanned influence points are distributed in areas with low population density. A total of 10,231 non-human affected grids were extracted (Figure 7). Unmanned influence points are mainly distributed in uninhabited areas such as plateaus, mountains, and deserts. Most of these places are far away from strong human activities areas such as cities and towns, located in nature reserves and ANPP low-value areas.

Spatial and Temporal Distribution
Based on the RF machine learning algorithm, the ANPP in the non-human affected grids were upscaled to the entire area to obtain NNPP (Figure 8a). The multi-year average value was 243.60 gC·m −2 ·yr −1 and an annual total of 276.75 TgC·yr −1 . The high value of NNPP was mainly distributed in the southern Shaanxi, the southeastern Qinghai, and the northwestern Xinjiang. The average annual NNPP showed an increasing trend at a rate of 1.81 gC·m −2 ·yr −1 (Figure 9). From 2001 to 2017, 80.81% of the grassland NNPP showed an increasing trend ( Figure 10). Among them, the areas with a significant increase accounted for 34.11%, mainly in the north regions of Shaanxi, middle of Gansu and Qinghai. However, 19.19% of the grassland showed a decreasing trend in the NNPP, and only 1.86% showed a significant decreasing trend, mainly in the southern Shaanxi Province.

Spatial and Temporal Distribution
Based on the RF machine learning algorithm, the ANPP in the non-human affected grids were upscaled to the entire area to obtain NNPP (Figure 8a). The multi-year average value was 243.60 gC·m −2 ·yr −1 and an annual total of 276.75 TgC·yr −1 . The high value of NNPP was mainly distributed in the southern Shaanxi, the southeastern Qinghai, and the northwestern Xinjiang. The average annual NNPP showed an increasing trend at a rate of 1.81 gC·m −2 ·yr −1 (Figure 9). From 2001 to 2017, 80.81% of the grassland NNPP showed an increasing trend ( Figure 10). Among them, the areas with a significant increase accounted for 34.11%, mainly in the north regions of Shaanxi, middle of Gansu and Qinghai. However, 19.19% of the grassland showed a decreasing trend in the NNPP, and only 1.86% showed a significant decreasing trend, mainly in the southern Shaanxi Province.     The spatial pattern of PNPP was similar to that of NNPP, with a multi-year average of 421.96 gC·m −2 ·yr −1 ( Figure 9) and a total of 479.4 TgC·yr −1 . The spatial distribution of PNPP (Figure 8b) overall decreased from southeast to northwest.     The spatial pattern of PNPP was similar to that of NNPP, with a multi-year average of 421.96 gC·m −2 ·yr −1 ( Figure 9) and a total of 479.4 TgC·yr −1 . The spatial distribution of PNPP (Figure 8b) overall decreased from southeast to northwest.

HNPP_N and HNPP_P
The HNPP_N > 0 accounted for 40.28% of the total grassland area (Figure 11a HNPP_N (Figure 13). Among them, the areas with significant increase trends accounted for 9.69%, mainly in the southern regions of Shaanxi, Gansu and Qinghai province. However, 41.14% of grassland HNPP_N in the northwest region showed decreasing trends, and 12.64% of the regions showed significant decreasing trends, which were mainly distributed in the northern Shaanxi, Tianshan Mountains, Kunlun Mountains, and the southwestern Altai Mountains. an increasing trend in 48.96% of grassland HNPP_N (Figure 13). Among them, the areas with significant increase trends accounted for 9.69%, mainly in the southern regions of Shaanxi, Gansu and Qinghai province. However, 41.14% of grassland HNPP_N in the northwest region showed decreasing trends, and 12.64% of the regions showed significant decreasing trends, which were mainly distributed in the northern Shaanxi, Tianshan Mountains, Kunlun Mountains, and the southwestern Altai Mountains.   an increasing trend in 48.96% of grassland HNPP_N (Figure 13). Among them, the areas with significant increase trends accounted for 9.69%, mainly in the southern regions of Shaanxi, Gansu and Qinghai province. However, 41.14% of grassland HNPP_N in the northwest region showed decreasing trends, and 12.64% of the regions showed significant decreasing trends, which were mainly distributed in the northern Shaanxi, Tianshan Mountains, Kunlun Mountains, and the southwestern Altai Mountains.    (Figure 13). However, HNPP_P > 0 accounts for 11.21% of the total grassland area (Figure 12b), which mainly distributed in the Qilian, Tianshan, and Altai Mountains. The HNPP_P < 0 accounted for 88.79% of the total grassland area (Figure 11b), with an average value of −224.58 gC·m −2 ·yr −1 . Negative values are widely distributed, and the lowest values are mainly distributed in southern Shaanxi and southern Gansu. The population density in the above-mentioned areas is relatively high, and human activities are intense.

The Relative Roles of Human Activities and Climate Change in Grassland Change
The spatial distribution of grassland ANPP change, including climate variation, hu-  (Figure 13). However, HNPP_P > 0 accounts for 11.21% of the total grassland area (Figure 11b), which mainly distributed in the Qilian, Tianshan, and Altai Mountains. The HNPP_P < 0 accounted for 88.79% of the total grassland area (Figure 11b population density in the above-mentioned areas is relatively high, and human activities are intense.

The Relative Roles of Human Activities and Climate Change in Grassland Change
The spatial distribution of grassland ANPP change, including climate variation, human activities, and combined effects of climate variation and human activities during the 2001-2017 were analyzed (Table 1 and Figure 14). In order to disentangle the relative effects of human activities and climate change in the process of grassland change, six possible scenarios were defined (Table 1) and evaluated by calculating the Slope of ANPP, NNPP, and HNPP_N.  From 2001 to 2017, 84.32% of the grassland showed increasing trends, and the degraded areas only accounted for 15.68% in the Northwest China. Human activity played a minor role in grassland change (16.34%), while the improvement of grassland was mainly attribute to climate. Human activities accounted for only 11.43% of areas with increasing NPP. However, the main factor of grassland degradation was human activity, accounting for 42.73% of the degraded area. In the study area, the population are sparse (29.36 people/km 2 ). The desert is widespread in this region, and the human activities is limited to the city, oasis, and its surrounding areas. During the past 30 years, precipitation had increased slightly, and the vegetation coverage had been increasing in the region [33]. From 2001 to 2017, 84.32% of the grassland showed increasing trends, and the degraded areas only accounted for 15.68% in the Northwest China. Human activity played a minor role in grassland change (16.34%), while the improvement of grassland was mainly attribute to climate. Human activities accounted for only 11.43% of areas with increasing NPP. However, the main factor of grassland degradation was human activity, accounting for 42.73% of the degraded area. In the study area, the population are sparse (29.36 people/km 2 ). The desert is widespread in this region, and the human activities is limited to the city, oasis, and its surrounding areas. During the past 30 years, precipitation had increased slightly, and the vegetation coverage had been increasing in the region [33]. This is consistent with the conclusion that grassland restoration is mainly dominated by climate change in previous study [34]. Grassland degradation is mainly dominated by human activities [15] and is distributed in major pastoral areas, such as central Inner Mongolia, Tianshan Mountains, and Altai Mountains. Overgrazing exceeds the carrying capacity of the grasslands, and then leads to grassland degradation [35].

NNPP and PNPP
During the study period, the average values of NNPP, PNPP, and ANPP were 243.60 gC·m −2 ·yr −1 , 421.96 gC·m −2 ·yr −1 , and 236.99 gC·m −2 ·yr −1 , respectively. PNPP was much higher than NNPP and ANPP. The area with higher PNPP values than NNPP accounted for 88.51% of the study area ( Figure 15). This may be explained as that the Thornthwaite Memorial Model just considers the climatic factors and ignores factors such as vegetation type, water and temperature stress on plant growth, and photosynthesis efficiency. The areas with lower PNPP values than NNPP are mainly distributed in the mountainous areas. These areas are widely covered by snow and glaciers, which absorb heat of the mountainous areas during the melting of snow and glaciers in summer and have a great impact on the temperature. Therefore, the actual evaporation in the Thornthwaite Memorial Model may be exaggerated, resulting in a smaller PNPP. aspect, solar radiation, and solar radiation of the growing season. Therefore, NNPP estimated by random forest is more reliable than PNPP.

HNPP_N Is More Consistent with the Actual Situation Than HNPP_T
The mean values of HNPP_N and HNPP_P were −7.45 gC·m −2 ·yr −1 and −189.15 gC·m −2 ·yr− 1 , accounting for 3.14% and 79.81% of ANPP, respectively. The former showed that human activities had weak negative effects, and the latter indicated strong negative effects. The spatial pattern of HNPP also showed great differences ( Figure 16). According to the theoretical premise of this study, the NNPP should have a strong correlation with ANPP. The correlation coefficient 0.95, which is much higher than the correlation between ANPP and PNPP (R 2 = 0.72). The NNPP estimation considers the parameters that were not included in the climatological model such as NDVI, altitude, slope, aspect, solar radiation, and solar radiation of the growing season. Therefore, NNPP estimated by random forest is more reliable than PNPP.

HNPP_N Is More Consistent with the Actual Situation Than HNPP_T
The mean values of HNPP_N and HNPP_P were −7.45 gC·m −2 ·yr −1 and −189.15 gC·m −2 ·yr −1 , accounting for 3.14% and 79.81% of ANPP, respectively. The former showed that human activities had weak negative effects, and the latter indicated strong negative effects. The spatial pattern of HNPP also showed great differences ( Figure 16).

HNPP_N Is More Consistent with the Actual Situation Than HNPP_T
The mean values of HNPP_N and HNPP_P were −7.45 gC·m −2 ·yr −1 and −189.15 gC·m −2 ·yr− 1 , accounting for 3.14% and 79.81% of ANPP, respectively. The former showed that human activities had weak negative effects, and the latter indicated strong negative effects. The spatial pattern of HNPP also showed great differences ( Figure 16). The HNPP_N > 0 (24.068 gC·m −2 ·yr −1 ) and HNPP_P < 0 (−197.120 gC·m −2 ·yr −1 ) in the Three-River Source Region (Figure 11), where the main human activity is grazing. Since 2003, the "Animal Reduction Project" has been implemented in the Three-River Source Region. By 2015, the stocking pressure index has decreased by 34.3% [36]. According to the statistical data from 2007 to 2017, and the inventory showed a significant decline of The HNPP_N > 0 (24.068 gC·m −2 ·yr −1 ) and HNPP_P < 0 (−197.120 gC·m −2 ·yr −1 ) in the Three-River Source Region (Figure 11), where the main human activity is grazing. Since 2003, the "Animal Reduction Project" has been implemented in the Three-River Source Region. By 2015, the stocking pressure index has decreased by 34.3% [36]. According to the statistical data from 2007 to 2017, and the inventory showed a significant decline of 146,700 stock unit (SU) per year ( Figure 17). The livestock reduction project has achieved significant result [3]. In addition to livestock reduction measures, the Three-River Source region has implemented a series of ecological restoration measures, such as grass fences, artificial replanting of grass seeds, artificial rainfall, and rodent control. The total area of returning farmland to forest and grassland is approximately 102.07 km 2 [33], and comprehensive management of 3.49 × 10 7 km 2 was implemented for improve grassland yield, forage quality, disaster resistance, and species diversity [3]. Therefore, the human activities in this area had a positive effect on grassland, indicating that the HNPP_N results were more accurate.  Figure 17). The livestock reduction project has achieved significant result [3]. In addition to livestock reduction measures, the Three-River Source region has implemented a series of ecological restoration measures, such as grass fences, artificial replanting of grass seeds, artificial rainfall, and rodent control. The total area of returning farmland to forest and grassland is approximately 102.07 km 2 [33], and comprehensive management of 3.49 × 10 7 km 2 was implemented for improve grassland yield, forage quality, disaster resistance, and species diversity [3]. Therefore, the human activities in this area had a positive effect on grassland, indicating that the HNPP_N results were more accurate.  (Figure 18). The increasing livestock leads to an increase in the negative impact of human activities and the degradation of grassland vegetation in the region [36]. During this period, the negative impact of HNPP_N had been increasing, while HNPP_P had been decreasing, and the former was Ningxia is one of the important pastoral areas in China. From 2001 to 2017, this area showed an increasing trend of livestock numbers ( Figure 18). The increasing livestock leads to an increase in the negative impact of human activities and the degradation of grassland vegetation in the region [36]. During this period, the negative impact of HNPP_N had been increasing, while HNPP_P had been decreasing, and the former was more consistent with the actual situation ( Figure 16). Ningxia is one of the important pastoral areas in China. From 2001 to 2017, this area showed an increasing trend of livestock numbers ( Figure 18). The increasing livestock leads to an increase in the negative impact of human activities and the degradation of grassland vegetation in the region [36]. During this period, the negative impact of HNPP_N had been increasing, while HNPP_P had been decreasing, and the former was more consistent with the actual situation ( Figure 16).

Conclusions
The NNPP estimated by ANPP based on the grids without effects of human activities proposed in this study can accurately reflect the ANPP of vegetation. Using this method, we found that human activities have a weak negative impact on the NPP of grassland in Northwest China, and the grassland degradation was mainly caused by human activities, but the restoration is mainly caused by climate warming and humidification. This method

Conclusions
The NNPP estimated by ANPP based on the grids without effects of human activities proposed in this study can accurately reflect the ANPP of vegetation. Using this method, we found that human activities have a weak negative impact on the NPP of grassland in Northwest China, and the grassland degradation was mainly caused by human activities, but the restoration is mainly caused by climate warming and humidification. This method can be used to evaluate the effectiveness of ecological environment treatment and provide quantitative basis for scientifically building the strategy to cope with climate change.
This study selected grassland in Northwest China as study area. This method may be also applicable to other types of vegetation and regions, and the main differences are the indices to identify grids without the effects of human activities. For example, forest land is a perennial plant, and the interannual changing pattern of NPP is different from that of grassland. In addition, the random forest algorithm is used to upscale the NPP in this study, other approaches such as deep learning methods may be also useful tools to study the temporal and spatial changes in NPP.