Quantifying the Impacts of Anthropogenic Activities and Climate Variations on Vegetation Productivity Changes in China from 1985 to 2015

: Accurate assessment of vegetation dynamics provides important information for ecosystem management. Anthropogenic activities and climate variations are the major factors that primarily influence vegetation ecosystems. This study investigates the spatiotemporal impacts of climate factors and human activities on vegetation productivity changes in China from 1985 to 2015. Actual net primary productivity (ANPP) is used to reflect vegetation dynamics quantitatively. Climate-induced potential net primary productivity (PNPP) is used as an indicator of climate change, whereas the difference between PNPP and ANPP is considered as an indicator of human activities (HNPP). Overall, 91% of the total vegetation cover area shows declining trends for net primary productivity (NPP), while only 9% shows increasing trends before 2000 (base period). However, after 2000 (restoration period), 78.7% of the total vegetation cover area shows increasing trends, whereas 21.3% of the area shows decreasing trends. Moreover, during the base period, the quantitative contribution of climate change to NPP restoration is 0.21 grams carbon per meter square per year (gC m − 2 yr − 1 ) and to degradation is 2.41 gC m − 2 yr − 1 , while during the restoration period, climate change contributes 0.56 and 0.29 gC m − 2 yr − 1 to NPP restoration and degradation, respectively. Human activities contribute 0.36 and 0.72 gC m − 2 yr − 1 during the base period, and 0.63 and 0.31 gC m − 2 yr − 1 during the restoration period to NPP restoration and degradation, respectively. The combined effects of climate and human activities restore 0.65 and 1.11 gC m − 2 yr − 1 , and degrade 2.01 and 0.67 gC m − 2 yr − 1 during the base and restoration periods, respectively. Climate factors affect vegetation cover more than human activities, while precipitation is found to be more sensitive to NPP change than temperature. Unlike the base period, NPP per unit area increases with an increase in the human footprint pressure during the restoration period. Grassland has more variability than other vegetation classes, and the grassland changes are mainly observed in Tibet, Xinjiang, and Inner Mongolia regions. The results may help policy-makers by providing necessary guidelines for the management of forest, grassland, and agricultural activities.


Introduction
Climate fluctuations and anthropogenic activities are the two major drivers of terrestrial ecosystem changes [1,2]. Abnormal climate variations have been coupled with the negative effects of human activities over the previous few decades, resulting in a series of environmental and ecological problems, such as air pollution, biodiversity losses, soil erosion, and vegetation ecosystem degradation [3]. Therefore, quantifying and separating the influences of these factors on vegetation ecosystems by emphasizing the dynamic variations in terrestrial ecosystems is an important research direction [4][5][6]. Vegetation is a key component of terrestrial ecosystems, and plays a prominent part in regulating the carbon balance, adapting land surface conditions, and emulating regional anthropogenic activities [7,8]. An ample number of studies have already been conducted on vegetation change and their driving factors (i.e., human activities and climatic variations) in recent years to investigate the vegetation dynamics at both regional and global scales [9,10].
The impact of climatic fluctuations on vegetation ecosystems is a global concern. During the last few decades, the temperature of the globe has increased rapidly. Consequently, the frequency of rainfall events is predicted to increase at the same rate, especially in humid regions [11][12][13]. Previous studies have found that the temperature increase extends the length of the growing season, and thus supports an increase in productivity of vegetation ecosystems [14,15]. Conversely, Zhang et al. [16][17][18] claim that temperature contributes negatively to vegetation productivity. As vegetation responds differently to different climatic factors, it is essential to quantify the response of temperature and precipitation to vegetation productivity, accounting for spatiotemporal heterogeneity. Human activities also contribute strongly to vegetation dynamics and influence the vegetation ecosystems both positively and negatively. Recently, China has experienced serious environmental issues due to vegetation degradation. Those environmental problems are mainly triggered by overcultivation and overgrazing in environmentally complex and ecologically fragile regions of the country [19,20]. To restore the vegetation coverage, as well as productivity, and to minimize the human pressure on the environment, a series of vegetation-related policies were introduced by China. The Grain for Green Program launched in 1999 is one of the most prominent initiatives [21]. The major objective of all the vegetation-related projects is to increase the vegetation area, either by planting trees or by converting the bare land and wasteland into grassland [22]. Some vegetation conservation programs were also introduced to support ecological restoration. For example, the Grazing Withdrawal Program (GWP) was introduced in 2003 to minimize the grazing pressure on grassland ecosystems by employing cultivated pastures through forbidding grazing [23]. However, in some areas of China, the effects of such ecological and conservation programs have still not been quantified. Some studies indicate that human activities play a negative role in vegetation dynamics during the base period, including rapid construction, overgrazing, and deforestation. However, they contribute positively during the restoration period because of afforestation and land use planning. Therefore, both human activities, as well as climate change, can alter vegetation statistics significantly [15,24,25].
Recently, many studies have documented that the comparison of actual vegetation productivity and potential vegetation productivity determines the human appropriation of vegetation productivity. Therefore, it is assumed that the difference between potential and actual productivities of vegetation represents the anthropogenic impacts [26]. Several researchers have followed the vegetation productivity difference method to quantify and separate the contributions of human activities and climate factors. Using net primary productivity (NPP) as an indicator of vegetation, previous studies indicate an exponential increase in vegetation productivity in China during the last couple of decades [27,28]. NPP is a very important indicator of vegetation growth status and exerts a significant influence on the carbon cycle of the global biosphere [14,29]. By definition, NPP represents the amount of carbon fixed through the process of photosynthesis by plants per unit of space and time. NPP is more sensitive to anthropogenic factors and climate variations than other vegetation indicators, precisely, reflect the responses of vegetation growth to human activities and climate change [30,31]. In particular, climate-induced potential NPP (PNPP) signifies the potential climatic conditions of vegetation growth without human interruption. Actual NPP (ANPP) is controlled by both human activities and climatic factors and indicates true conditions of vegetation growth. Therefore, the difference between PNPP and ANPP (PNPP-ANPP) is used to define human-induced NPP (HNPP) [28,32]. Thus, vegetation growth variations can be determined at a pixel scale by comparing the temporal variations of PNPP, ANPP, and HNPP [33].
Previous literature is quantitatively discriminated relative contributions of climatic factors and human activities in vegetation change using the NPP-based scenario simulation method. This method is supposed to achieve reasonable accuracy, as described in several studies [34]. For example, Xu et al. [23] studied the relative contributions of anthropogenic activities and climate change to assess the grassland dynamics of the Qinghai-Tibet Plateau at different time scales and used the scenario simulation method based on three types of NPP trends. Similarly, Zhang et. al. [35] applied the same method to investigate the major factors of grassland dynamics in Xinjiang. Furthermore, Yang et. al. [36] and Gang et al. [37] applied the scenario simulation method on both regional and global scales and indicated that the NPP-based scenario simulation method can successfully be applied for quantitative evaluation of the relative contribution of anthropogenic factors and climate variations in grassland dynamics at any scale.
To provide better quality of living conditions for human well-being, vegetation restoration and conservation is very important to deal with environmental and ecological issues in a country, which are mainly triggered by climate and human factors. Although several studies have already been conducted to quantify the relative contributions of climate variations and anthropogenic factors, relatively less attention has been paid to evaluating the effectiveness of these ecological restoration programs at the national level, considering all of the green infrastructure in China. Therefore, it is important to rigorously monitor the vegetation restoration activity to investigate the efficiency of these programs. In view of the above research gap, this study investigates 31-year vegetation productivity dynamics across China, considering all the green areas in the country. The study period is divided into a base period and restoration period to evaluate the effectiveness of vegetation conservation and restoration policies. The specific objectives of this study are: (1) to evaluate the spatiotemporal change pattern and sustainability trends of vegetation growth; (2) to identify the contributions of climate change and anthropogenic activities to vegetation productivity change; (3) to classify the roles of individual climatic parameters and human footprint pressure in vegetation dynamics. The work is done for both the base period (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999) and restoration period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015). The findings of this study can provide a theoretical reference to assist policy-makers for the optimized management of green land and accomplish the sustainable use of ecological resources.

Climate Data
The China Meteorological Forcing Datasets (CMDF) product is used for climatic parameters, including accumulative annual precipitation and mean annual air temperature (http://data.tpdc.ac.cn/en/data/8028b944-daaa-4511-8769-965612652c49/) [38]. CMDF is a gridded product with a spatial resolution of 0.1 degrees and a temporal resolution of 3 hours. The air temperature data are developed by merging the Princeton forcing data and the observed data from 740 ground stations covered by the China Meteorological Department (CMD). The precipitation data are produced by combining the CMD records with the Global Land Data Assimilation System (GLDAS) precipitation data from 1979 to 1998, and the Tropical Rainfall Measuring Mission (TRMM) rainfall product (3B42) for the period ranging 1998 to 2018 [39]. These datasets have already been used in previous studies to examine the vegetation response to climate variations [40,41]. The datasets are converted to annual records (mean annual temperature and annual accumulative precipitation) and resampled to 1 km 2 spatial resolution by applying the nearest neighborhood algorithm to meet the cell size of the NPP product. Figure 1 shows the spatial patterns of mean annual records of climate data for base and restoration periods.

Human Footprint Data
Human footprints (HFP) describe the intensity of human influence on terrestrial environments. The recently published HFP (https://www.nature.com/articles/sdata201667) datasets are used to investigate the relationship between human influence and net vegetation productivity changes [42] for both base and restoration periods. This is a globally-standardized gridded dataset at a spatial resolution of 1-km 2 for the years 1993 and 2009. The values of HFP data range from 0 to 64, whereby a lower value corresponds to the low intensity of human influence and vice versa. Higher values are observed in cities. However, this study was conducted only for vegetation cover and neglects most of the city areas (non-vegetated area). Therefore, the HFP values for this study range from 0 to 36. This system has already been used to explore the impacts of human activities on eco-environmental changes, including vegetation changes [43] and biodiversity loss [44]. The 1993 HFP data are used for base period, whereas the 2009 HFP data are used for the restoration period. Although these datasets are not compatible with the study period, they can provide an overall estimation of human impacts on vegetation productivity dynamics for both periods. The spatial patterns of HFP datasets are shown in Figure 2.

Vegetation Cover Data
To obtain the vegetation area, a vegetation cover map for China was acquired from the MCD12Q1 moderate-resolution imaging spectro-radiometer (MODIS) product. As described earlier, the vegetation of China increased significantly after the implementation of the Grain for Green Program; therefore, the product was downloaded for the year 2015 to account for the maximum possible vegetation area for the analysis. MCD12Q1 encloses five classification schemes, while the MODIS land cover product for the year 2015 was reclassified according to the classification system developed by the International Geosphere-Biosphere Program (IGBP) for this study [13,45,46]. The IGBP scheme classified the land cover into 17 classes, with 11 natural vegetation and three mosaic vegetation classes. The MCD12Q1 land cover product was reclassified to obtain three major vegetation classes: mixed forestland, mixed grassland, and mixed cropland. The reclassification scheme is given in the Table 1, and the spatial distribution of the vegetation classes is shown in Figure  3.

Net Primary Productivity (NPP) Data
The monthly actual net primary productivity (ANPP) product (ChinaNPP_1985_2015) of China's terrestrial ecosystems at 1 km 2 spatial resolution for 1985-2015 was used for this study [47]. ANPP data were downloaded from http://www.geodoi.ac.cn, and then monthly images were aggregated to an annual product to obtain total annual ANPP in grams carbon per meter square per year (gC m -2 yr -1 ). These data were produced using the Carnegie-Ames-Stanford approach (CASA) model, together with the meteorological station data, making sure the same CMDF data source was used to calculate ANPP and PNPP. The ANPP data were validated by comparing them with simulated and measured data published in different articles. These data were preferred over the MODIS NPP product (MOD17A3) to consider the base period, as MODIS data are not available before the year 2000. The CASA model was applied using the following equation: where ( , ) is the NPP value of a pixel x at time t (gC m −2 ), ( , ) is the value of the photosynthetically active radiation of a pixel x at time t (MJ m −2 ), and ( , ) is the value of the light energy use efficiency of a pixel x at time t.
The absorbed photosynthetically active radiation (APAR) was calculated using the following equation [48]: where SOL(x, t) is the solar radiation value of pixel x at time t (MJ/m 2 ), FPAR(x, t) represents the incoming fraction of photosynthetically active radiation absorbed by vegetation at pixel x at time t, and the constant (0.5) is a fraction of solar radiation utilized by green vegetation. Due to limitations of long term data availability, ( , ) was calculated using two different data sources: advanced very high-resolution radiometer (AVHRR) data were used to determine ( , ) from 1985 to February 2000, and MODIS data were used to determine ( , ) from March 2000 to 2015. The calculations were made based on the following equation [47]: and ( , ) are minimum and maximum values or SR, respectively.
where ( , ) is the pixel value of NDVI at time t. The light energy use efficiency ( , ) was calculated using the following equation: where ( , ) and ( , ) are the low and high effects of temperature on light energy conversion, respectively; ( , ) is the water stress; * is the maximum possible light energy use efficiency; and an * value of 0.389 is used for AVHRR data, while 0.55 is used for MODIS data.
The ANPP product was further validated using ChinaFLUX tower data and field measurements of three forest sites.
The selected sites for the field measurements included Qilian Mountain, Changbaishan Natural Reserve, and Dunhua County. The Changbaishan Natural Reserve and Dunhua County are located in Jilin province and Qilian Mountain is in Gansu province, as shown in Figure 4a The empirical relationship between above ground biomass (W), diameter at breast height (DBH) and tree height (H) was used to calculate biomass for different trees [49]. The equation is given below: where and are the parameters of the equation. The structural parameters (H and DBH) were collected for each tree in each sample plot, where 10, 14, and 10 plots were selected, respectively, to measure the biomass [50]. The biomass of each sample plot was calculated using a species-specific formula of relative tree growth [51]. The ANPP (gC m −2 year −1 ) for each sample plot was then obtained based on the biomass difference between two successive years. A strong relationship (R 2 value of 0.84) between simulated and measured ANPP values was found, as shown in Figure 4d.
The ANPP product was also validated using ChinaFLUX tower measurements. The data of 45 sites were obtained from http://www.cnern.org.cn/data/initDRsearch?classcode=DPAPER. The selected sites had at least one year of continuous measurements. Due to the footprint issues, the simulated ANPP pixel values were extracted using a 3×3 pixel window around the flux tower. The

Methodological Framework
The Figure 5 illustrates the overall methodological framework used for this study. Three types of datasets were used for this study: NPPs, climate data, and human footprint data. The monthly ANPP product (ChinaNPP_1985_2015) of China's terrestrial ecosystems at 1 km 2 spatial resolution for 1985-2015 was used for this study (Section 2.1.4). PNPP was derived using the Thornthwait Memorial (TM) model (Section 2.2.1), whereas HNPP was taken as the difference between PNPP and ANPP. The impacts of climate and human factors were investigated based on 8 scenarios by comparing the slopes of ANPP, PNPP, and HNPP (Section 2.2.2). The relative contributions of climatic factors to NPP change were determined by calculating the pixel-based Pearson's correlation of ANPP with temperature and precipitation, and were analyzed based on 8 scenarios (Section 2.2.4).

Figure 5.
Methodological framework adopted for this study. NPP, net primary productivity.

PNPP and HNPP Calculation
The climate-influenced potential net primary productivity (PNPP) was estimated using the Thornthwaite memorial (TM) model. The TM model is a modified version of the Miami model, and it is a more accurate approach than the Miami model. It is based on the climatic data used in the Miami model and was developed using Thornthwaite's potential evaporation model [54]. The TM model employs the connection between carbon uptake and evapotranspiration to calculate pixel values of PNPP using temperature and precipitation as inputs. The basic principle of TM-modelbased PNPP estimation is expressed in the following equation: = 3000 + 25 + 0.05 where indicates total annual potential NPP (gC m −2 y −1 ), is the average annual actual evapotranspiration (mm), is the average annual potential evapotranspiration, is the total annual precipitation (mm), and is the mean annual temperature (ᵒ ) .
The saturation value of dry matter (DM) cannot exceed 3000 gDMm -2 year -1 . In principal, the TM model estimates PNPP values by employing the relationship between carbon uptake and evapotranspiration. Temperature and precipitation are used as inputs for the calculation of evapotranspiration. A thorough explanation of this model is discussed by Lieth [27,54,55] in detail.
The human-induced NPP (HNPP) is estimated by calculating the difference between PNPP and ANPP [27,32,[35][36][37]56], as expressed below The value is shown in gC m −2 yr −1 ; the negative values represent NPP gain due to human activities, whereas the positive values represent a loss in NPP due to human activities [27]. In contrast, positive values of ANPP and PNPP represent NPP gain, while negative values of ANPP and PNPP represent NPP losses in gC m −2 yr −1 . The spatial patterns of ANPP, PNPP, and HNPP for all vegetation classes are shown in Figures 6 and 7.

Slopes of Climate Variables and NPP and Scenario Development
The linear regression method was adopted to compute the spatial change patterns of vegetation productivity and climate variables. The time series of pixels were calculated to obtain the slope coefficients of a trend line using the following formula: An F-test was then applied to detect the statistical significance of spatial change patterns. The following formula was applied to test the significance of slopes calculated from the above equation: where U and Q are the regression sum of squares and residual sum of squares, respectively; and is the total number of years. The climate and human impacts on NPP change were investigated by comparing the slope of ANPP (SA), slope of PNPP (SP), and slope of HNPP (SH), as described in Table 2.  [57] and was modified in 1969 by Mandelbrots [58]. In recent studies, it was effectively applied to investigate the long term consistency of time-series data to predict future trends [59]. The following steps are involved in H value calculations: • Divide the time series ( ) ( = 1, … , ) into sub-series ( ), and for each sub-series = 1,…, τ.
• Then define the mean sequence of time series, • Compute the cumulative deviation of each , • Generate the range sequence ( ) , ( ) = max ( , ) − min ( , ) , = 1, … , • Then calculate the standard deviation sequence ( ) , • At the end, rescale the calculated range, where, H is the Hurst exponent, and can be calculated by fitting the above equation as follows: log( ) ( ) = + ( × log( )) Literature indicates that the range of Hurst exponents varies from 0 to 1, and is divided into three categories for future predictions [27,54]. If H > 0.5, the trend NPP time series would be considered to be consistent with the study period. A value H < 0.5 would indicate inconsistent behavior for NPP time series, whereas H = 0.5 would mean the time series trend of NPP is absolutely random, with no consistency.
The degree of vegetation change and its consistency were evaluated based on spatial change trends of ANPP, its significance, and Hurst values (Figure 8). The eight possible scenarios were developed for vegetation sustainability assessment. The scenarios are described in Table 3.  The spatial patterns of Pearson's correlation coefficient between climate factors and ANPP were calculated to investigate the individual response of climate variations to ANPP change by using the following equation [60]: The relationship of ANPP with precipitation (Pcor) and temperature (Tcor) was then assessed with areas of vegetation change due to climate variations to identify the dominant climatic factor. This assessment was done based on eight possible scenarios, which are expressed in Table 4.

Drivers of Vegetation Dynamics
Vegetation dynamics were investigated based on two driving forces: (1) climate variations and (2) human activities. Figure 9 indicates the spatial change patterns of three vegetation classes. Spatial patterns were developed based on eight scenarios, defined in Table 1. Figure 9 indicates the areas of NPP changes due to the individual contribution of climate variations and human activities and the combined effects of both factors, whereas the areas unaffected by these two driving forces are neglected. The overall pattern indicates an NPP decrease during the base period, and conversely increasing patterns can be seen during the restoration period. The spatial patterns of forest and agriculture classes for the base period and the restoration period are almost identical to each other; grassland shows variability in the results. To be precise, grassland restoration and degradation are observed during the base and restoration periods, respectively, driven by both climate variation and human activities. However, climate variation affected grassland ecosystems slightly more than human activities.
The results in Figure 10 illustrate that during the base period, there is no significant restoration, and only 9% of the total vegetation area reflects increasing trends. Grassland presents 8% of the total and 9% of the restored area, which is mainly driven by climate change, whereas the other two vegetation classes show only 1% restoration of the total vegetation area. In contrast, the restoration period indicates an increasing trend of over 78.7% of the total vegetation cover. Both driving forces contribute significantly to NPP restoration. However, climate variations affect a larger area than human activities for grassland and forests, while human activities affect the cropland restoration area slightly more than climate-driven NPP change. The restoration due to human activities was indeed the result of the successful implementation of ecological restoration and conservation programs. Figure 11 shows the NPP degradation patterns, where 91% of the total vegetation area during the base period shows decreasing trends, whereas 9% of the total vegetation area shows increasing trends. During the base period, the climate is found to be the dominant factor of vegetation degradation in comparison to human activities. Climate affects 52% of the total vegetation cover negatively. Human activities and the combined contribution of both factors are 11% and 28%, respectively. The dominant negative climate contribution can be seen in all three classes. However, human activities also contribute to grassland degradation. The restoration period indicates that grassland degradation is mainly driven by climate variations, however the contribution of human activities cannot be ignored. Cropland and forest show slight degradation, and human activities and climate factors contribute almost equally to NPP decrease during the restoration period.   The quantitative measures of NPP dynamics for both driving forces are given in Table 5. The combined effects of climate change and human activities are found to have more impact on quantitative changes of vegetation. For the restoration period, the combined effects of climate change and human activities contribute the highest amount of productivity per unit area in vegetation restoration (1.11 gC m −2 yr −1 ) and degradation (0.67 gC m −2 yr −1 ). Whereas, climate change contributes the lowest amount of productivity per unit area in vegetation restoration (0.56 gC m −2 yr −1 ) and degradation (0.29 gC m −2 yr −1 ). For the base period, again the combined effects of climate variations and human activities are the main contributors in vegetation productivity changes per unit area, with a combined restoration amount of 0.65 gC m -2 yr -1 and degradation amount of 2.01 gC m −2 yr −1 . However, climate variations have a minimal impact on vegetation productivity restoration (0.21 gC m −2 yr −1 ) and human activities contribute minimally (0.72 gC m −2 yr −1 ) to decreasing the vegetation productivity per unit area, as described in Table 5. It was challenging to quantitatively evaluate our results. Initially, we validated our results against land cover changes (such as plantation and deforestation) using the MODIS (MCD12Q1) product for 2001 and 2015. Then, six random pockets from the areas of significant changes after 1999 were selected and compared with high-resolution imagery of Google Earth for the same season, as shown in Figure 12. Pockets 1, 3, 4, 5, and 6 show vegetation restoration, whereas pocket 2 indicates vegetation degradation due to construction land in that area (Heilongjiang Province). Plantation is observed in pockets 1, 5, and 6 (areas of the Loess Plateau and Inner Mongolia), where our results indicate anthropogenic impact plays a dominant role (covers 89% area of these three pockets). Vegetation restoration is observed due to the combined effects of climate and human activities in pockets 3 and 4 (as some plantation patches are found in those areas). In pocket 2, our results indicate that the anthropogenic impact also plays a dominant role (covers 97% area of pocket 2). Obvious vegetation changes are observed in the six randomly selected pockets. Additionally, the pockets selected from the areas of NPP restoration due to human activities show plantation in Google Earth images, whereas NPP degradation areas indicate replacement of vegetated areas into impervious surfaces.

Vegetation Sustainability and Degree of Change
Spatial patterns of vegetation sustainability and degree of change were determined as defined in Table 3. Spatial patterns ( Figure 13) indicate that during the period from 1985 to 1999, an overall decrease is observed in all vegetation types. Individually describing the vegetation classes for the base period, it can be seen that only grassland restoration is found in Qinghai-Tibet, Xingjiang, and Inner Mongolia regions. In addition, forest and croplands do not show considerable restoration.
Conversely, strong vegetation restoration patterns are observed from 2000 to 2015. All the vegetation classes increase either consistently or inconsistently. Grassland restoration is mainly found in the middle part of China, whereas forest and cropland restorations are observed in the middle, northeastern, and southern regions of China. During the restoration period, some grassland degradation patches are observed across all of China. Similarly, forest and cropland degradation spots are also observed. However, those are not threatening under slight degradation conditions.

Individual Contributions of Climate Factors and Human Activities
The individual contributions of climate factors on vegetation changes were investigated to quantify the most sensitive parameter to vegetation changes. For this purpose, pixel to pixel Pearson's correlation between ANPP and climatic factors was calculated and climate-dominated vegetation change areas (Figure 9) were compared based on the scenarios defined in Table 4. Figures 16 and 17 indicate that precipitation is a more sensitive indicator of vegetation changes in comparison to temperature. Both factors contribute differently in different geographic locations. Another point to be noted is that approximately 40% of the base period is not affected by either climatic indicator, but NPP is changed due to climatic impact in those areas. This highlights the importance of the climatic factor, which has a significant impact on vegetation productivity changes. Forest Grassland Cropland Precipitation-dominated increase (PDI) and both factors dominating the increase (BDI) are the prominent contributors in vegetation productivity increase during the restoration period, as illustrated in Figure 16. During the base period, both precipitation and temperature contributed equally in vegetation change; however, approximately 40% of the degraded area is controlled by an unknown climatic variable.

1985-1999 2000-2015
The human influence on vegetation dynamics was both positive and negative. It can be seen that factors such as rapid construction, economic revolution, and overgrazing contribute negatively, whereas proper landscape policies and plantation through ecological restoration and conservation programs contribute positively. This study uses the human influence datasets for the years 1993 and 2009. Due to the unavailability of temporal human footprint datasets at the yearly temporal resolution, 1993 data were used to analyze the NPP changes for the base period and 2009 data were applied to assess the NPP change at different human footprint pressure levels for the restoration period. Table 6. Effect of human footprint pressure on NPP dynamics for the base and restoration periods. For this analysis, NPP dynamics due to human activities mapped in Figure 9 were incorporated with human footprint pressure data. Results for the base period (Table 6) indicate that human pressure does not have much effect on NPP restoration; however, NPP decreases in higher footprint regions. During the restoration period, higher footprint pressure led to both positive changes and negative changes ( Table 6). Figure 18 indicates NPP changes for three vegetation classes. An overall decrease and increase can be seen for all vegetation classes for the base and restoration periods, respectively. At different intensity levels of human footprints, the maximum NPP restoration is observed during the restoration period. When HFP crosses a threshold of 20, net NPP gain starts decreasing. This could be due to the highly populated areas (cities). For the base period, consistent behavior of NPP degradation is observed in the areas with HFP > 5.

Discussions
Anthropogenic activities and climate fluctuations are the two dominant factors of changes in vegetation productivity. Several studies have already been conducted to distinguish the impact of these two factors on NPP dynamics by comparing ANPP, HNPP, and PNPP. However, most of them focused on grassland dynamics [36,61,62]. The Grain for Green Program was launched in 1999 in China, and since then a strong increase in grassland productivity has been observed [53]. Therefore, it is essential to quantify the impact of ecological restoration and conservation programs in China for all types of vegetation. The scenario simulation approach coupled with NPP is quite simple and is useful for differentiating the relative contributions of anthropogenic activities and climatic factors to vegetation productivity dynamics [57]. This approach enables accurate evaluation of the spatiotemporal variations in vegetation productivity driven by human activities and climate change [63]. Therefore, an NPP product (ChinaNPP_1985_2015) for China's terrestrial ecosystems was selected as an actual NPP measure for this study [47]. Precipitation-and temperature-based Thornthwaite memorial model were applied to simulate the potential NPP of vegetation classes. The Hurst exponent was applied to assess the consistency of vegetation spatial patterns, and a pixel-based Pearson's correlation was calculated to determine the individual impact of temperature and precipitation on NPP change.
The degree and sustainability of spatiotemporal NPP variations were determined by computing pixel values of the Hurst exponent, temporal change trends (slope), and significance level (p-value) of ANPP for both base and restoration periods. The study found an overall decreasing trend for the base period and an increasing trend for the restoration period, and grassland showed greater variability in comparison to other vegetation classes. Results (Figures 6-8) suggested that only grassland class represented a considerable portion (8%) of the restored area during the base period, while 12% grassland degradation of the total vegetation area was found during the restoration period. Grassland productivity changes were mainly observed in Xinjiang, Qinghai-Tibet Plateau, and Inner Mongolia. Previous studies concluded that the grassland changes in Qinghai-Tibet plateau are mainly driven by climate change during both base and restoration periods [26,64]. Tian et al. [65] concluded that both climate and human factors contributed equally to vegetation restoration in Inner Mongolia after 1999. Similarly, Zhang et al. [35] found an overall increase in the grassland area of Xinjiang province. However, human activities (overgrazing) and climate factors contributed equally to vegetation degradation. Forests and cropland classes showed decreasing and increasing trends during base and restoration periods, respectively, in most areas of China.
As discussed earlier (see Section 3.1), vegetation productivity is primarily driven by anthropogenic activities and climate variations. NPP dynamics are generally influenced by two climatic factors: precipitation and temperature. Variations in temperature and precipitation affect the vegetation ecosystems by altering soil microbes and soil moisture, and by affecting the processes of plant respiration and photosynthesis. Consequently, these changes may control the productivity of vegetation ecosystem by regulating plant growth [2,66]. The results of this study suggest that climate change is an important factor in controlling grassland productivity (Figures 9-11 and Table 5), and several studies concluded that the contribution of climate change in vegetation dynamics is greater than human activities [36,67,68]. The outcomes of the current study partially agree with these conclusions. Our findings indicate that climate variations affect vegetation areas in China more than human factors, but contribute slightly lesser than anthropogenic factors in NPP change per unit area. Similarly, Liu et al. [69] studied grassland dynamics due to climate change and human activities from 2000 to 2013 and concluded that the relative contributions of climate and human activities for grassland dynamics of China were 41.45% and 45.22%, respectively. In the current study, during the base period, climate change and human activities contributed approximately 43% and 22% to grassland dynamics, respectively, whereas 35% of the area was affected by both factors. However, during the restoration period, climate change and human activities contributed approximately 44% and 35% to grassland dynamics, respectively. However, 21% of the area was affected by both factors. Our statistics were different from the results of Liu et al. [69]; climate variations affected a larger area but quantitatively (NPP change per unit area) contributed less than human activities. Furthermore, Zhang et al. [35] studied the grassland dynamics of Xinjiang province and concluded that climate and human factors contributed equally to NPP decrease, whereas human activities played a relatively greater role in NPP restoration. Similarly, Tian et al. [65] studied vegetation dynamics of Inner Mongolia and found that both human and climate factors played almost equal roles in vegetation restoration. The results of this study also indicated that in dessert areas, climate affected 39% of the area, human activities affected 37% of the area, and 24% of the vegetation area was affected by both factors. Therefore, the results are not very different from Tian et al. and Zhang et al. [35,65]. Similarly, Huang et al. [64] and Li et al. [26] concluded that the vegetation dynamics of Qinghai-Tibet plateau are primarily driven by climate variations, and precipitation is the leading climatic factor of vegetation growth. However, human activities had a negative impact on vegetation growth in most of the areas. On the other hand, the results partially agree with their conclusions, as the results indicate that vegetation productivity in Qinghai-Tibet plateau increased by 14.5%, and 8.8% of the total vegetation area was decreased due to human activities.
The relative contributions of climate factors suggest that precipitation is the dominant factor in controlling NPP changes (Figures 16 and 17). We also observed consistency between precipitation, ANPP, and PNPP variations in most of the areas, suggesting that vegetation productivity in China is more sensitive to precipitation than temperature, which is similar to the conclusions of Yang et al. [36]. Chen et al. [70] indicated that water is a key factor of vegetation growth, and vegetation with a high root-to-shoot ratio and low leaf area index (LAI) may be ascribed to water shortage, which severely affects vegetation productivity [71][72][73]. The results of this study also coincide with Zhou et al. [3] and Zheng et al. [74], suggesting that wetter climates are more favorable for vegetation productivity growth [3,754. During the restoration period of the current study, the pixels with declining precipitation trends were mostly located in the areas with decreased ANPP (i.e., Xinjiang, Inner Mongolia, and Tibet regions). This means that the climate became drier in the areas with increased temperature, and NPP degradation due to climate variations can be expected in such regions ( Figure 9). This observation can probably be explained as being due to warming, with a lower amount of precipitation leading to drought intensification [75]. Furthermore, some studies indicate that temperature increase can accelerate the photosynthesis process and carbon assimilation of plants. However, vegetation productivity is reduced when the temperature crosses an ecosystem threshold [73,76], whereas the increasing trend of precipitation in Inner Mongolia is considered to be a dominant climatic factor for the restoration of vegetation productivity. In a similar context, Liang et al. [29] studied the climate response to spatiotemporal dynamics of ANPP in China and concluded that vegetation ANPP in arid and semi-arid regions is mainly controlled by precipitation, which is similar to the current study. In particular, our results suggest that climate-dominated NPP restoration is either controlled only by precipitation or by the combined effects of precipitation and temperature. Additionally, temperature-dominated areas showed an overall decrease in productivity.
Several studies state that human activities play a significant role in vegetation productivity changes [2,36,77]. The successful implementation of the Grain for Green Program has resulted in considerable vegetation development due to anthropogenic activities. The Grain for Green Program was launched in 1999 to change waste, barren, and grazing land into forest and grassland in fragile regions of the country [78]. Since then, vegetation degradation has been reduced, and improving patterns of vegetation productivity are reflected in our results and in other recent studies as well [3], providing evidence of successful implementation of this project. With the vegetation restoration programs, vegetation conservation programs also supported the NPP development after 1999. For example, GWP was introduced in 2003 to minimize the grazing pressure on grassland by employing cultivated pastures and by forbidding grazing [23]. In contrast to the base period, the results (Tables 6 and 7) suggest that high human pressure leads to ANPP increase during the restoration period. This indicates that local communities play a major role in vegetation restoration and conservation after 1999, and these results are in line with Xu et. al. [34]. Moreover, Zhang et al. [35] observed that anthropogenic activities play a prominent part in grassland restoration of Xinjiang, particularly through grassland conservation programs. Similarly, Wang et al. [78] found that human activities restored NPP by 5923 gigagram carbon (GgC) in the Tibet region. Xu et al. [23] showed that after implementation of the GWP, livestock numbers decreased by 7.8% in the Qinghai-Tibet Plateau, indicating a positive aspect of decreasing the grazing pressures. Similarly, Mu et al. [79] suggested that overall vegetation productivity improved in Inner Mongolia after the successful implementation of vegetation restoration and conservation programs. However, grazing pressure has not yet been reduced effectively, which could lead to vegetation degradation in the future [79]. This may be the reason for the NPP decrease in some areas of Inner Mongolia that is indicated in the current study.
There are some methodological shortcomings and data limitations of this study, which may be improved in the future. The human and climate factors are mainly responsible for NPP dynamics and trends. However, the HNPP (PNPP-ANPP) represents not only the anthropogenic activities, but also other factors, such as forest fires, vegetation species, and vegetation disease outbreaks. The influences of some factors may play a noticeable role in influencing NPP trends in some parts of China, and require quantitative analysis in the future, although temperature and precipitation are the dominant climate controlling factors ( Figure 16) for vegetation productivity changes. Other factors may also be important for NPP variations. Some studies suggest that solar radiation also affects vegetation trends [25,59]. This should also be considered in the future, especially with grazing data, to account for the maximum possible impact factors. The applied assumptions in this study may lead to certain. However, this assessment, which considers all vegetation types in China for the base and restoration periods, may provide necessary guidelines for forest, grassland, and agricultural management activities.

Conclusions
This study was conducted to quantify the impact of the "Grain and Green Program" introduced in 1999 for the restoration of the green infrastructure of China. In response, we assessed the degree of NPP change and its consistency to investigate the future scenarios, the contributions of anthropogenic and climate factors to NPP variations, as well as the individual contributions of climate indicators and human influence. Our results show a strong improvement in NPP per unit area for all types of vegetation, including forestland, grassland, and croplands, after 1999. Based on the findings of this study, we conclude the following: • Overall, the restoration period demonstrates increasing NPP trends for 78.7% and decreasing trends for 21.3% of the total vegetation cover, whereas the base period shows noticeable degradation in 91% and slight improvement in 9% of the total vegetation area.
• The consistently and significantly increasing NPP trends cover 29.64% of the total vegetation cover during the restoration period. Conversely, consistent and significant decreasing NPP trends are observed in 17.26% of the total vegetation cover during the restoration period.
• Both climatic factors and human activities contribute significantly to vegetation productivity changes. However, in this study, climate variations influence NPP changes in vegetation areas more than human activities.
• The quantitative contribution of climate to NPP restoration is found to be 0.21 and 0.56 gC m −2 yr −1 during the base and restoration periods, respectively, whereas the quantitative contribution of climate to NPP degradation is found to be 2.41 and 0.29 gC m −2 yr −1 during the base and restoration periods, respectively.
• The quantitative contribution of human activities to NPP restoration is found to be 0.36 and 0.63 gC m −2 yr −1 during the base and restoration periods, respectively, while the quantitative contribution of human activities to NPP degradation is calculated to be 0.72 and 0.31 gC m −2 yr −1 during the base and restoration periods, respectively.
• The combined effects of climate and human activities help to restore the maximum quantity of NPP (0.65 gC m −2 yr −1 ) during the base period and during the restoration period (1.11 gC m −2 yr −1 ).
• Precipitation seems to be a dominant factor in controlling NPP changes, whereas human footprint pressure has a positive impact on NPP changes for the restoration period and a negative impact on NPP changes for the base period.
Furthermore, grassland shows more variability during both base and restoration periods than cropland and forestland. The major grassland changes are observed in the areas of Qinghai-Tibet Plateau, Xinjiang, and Inner Mongolia. Some deforestation is observed in north-eastern and southeastern parts of China during the restoration period. Our results may provide necessary guidelines and facilitate policy-makers in the management forest, grassland, and agricultural activities.
Author Contributions: Shahid Naeem and Yongqiang Zhang developed the idea for this study. Yongqiang Zhang supervised and Shahid Naeem executed the research plan. Faisal Mueen Qamar and Jing Tian were involved in various discussions, while Amir Latif and Pranesh Kumar Paul provided technical support.