Change Trend and Restoration Potential of Vegetation Net Primary Productivity in China over the Past 20 Years

: As an important vegetation parameter and ecological index, vegetation net primary productivity (NPP) can intuitively reﬂect changes in the ecological environment and the level of the carbon budget. However, the change trend of NPP and its recovery potential in China over the past 20 years remain unclear. Here, we used trend analysis, multiple regression analysis and residual analysis methods to analyse the change trend in the NPP of China’s terrestrial ecosystems from 2000 to 2019, as well as the climax background, restoration status and restoration potential of the NPP of forest, grassland and desert ecosystems. The results showed that (1) the change in vegetation NPP in China from 2000 to 2019 showed a continuous upward trend, with a change slope of 2.39 gC/m 2 /a 2 , and the area with a positive slope of change accounted for 68.10% of the country’s land area. The contribution rates of meteorological conditions and human activities to vegetation NPP changes were 85.41% and 14.59%, respectively. (2) The results obtained by the regression analysis method of meteorological conditions based on nature reserves could reﬂect the zonal climax vegetation status to a large extent, and the obtained values had a smooth transition within each ecogeographical division and between each ecogeographical division, which truly reﬂected the law of gradual change in climate, vegetation and natural conditions. The annual total NPP of the climax background vegetation in China’s forest, grassland and desert ecosystems was approximately 2.76 ± 0.28 PgC, and the annual total NPP of the three ecosystems was 1.90 ± 0.2 PgC, 0.80 ± 0.07 PgC and 0.009 ± 0.0005 PgC, respectively. (3) The annual total vegetation NPP of the restoration status of China’s forest, grassland and desert ecosystems was 2.24 PgC, and the annual total vegetation NPP of the three was 1.54 PgC, 0.65 PgC and 0.007 PgC, respectively. Beneﬁting from the effective implementation of climate warming and humidiﬁcation and ecological engineering, the agro-pastoral zone, the Loess Plateau, the eastern Sichuan Basin and the Greater Khingan Range had the most signiﬁcant increases in the past 20 years. (4) The annual total vegetation NPP of China’s forest, grassland and desert ecosystem restoration potential was approximately 0.52 ± 0.28 PgC, which accounted for approximately 19.05% of the annual total NPP of the climax background vegetation. The annual total vegetation NPP of forest, grassland and desert ecosystems restoration status was 0.36 ± 0.2 PgC, 0.16 ± 0.07 PgC and 0.002 ± 0.0005 PgC, respectively; the restoration potential accounted for 18.80%, 9.67% and 23.95% of the climax background vegetation NPP, respectively. The deployment of ecological projects should fully consider the restrictive climate conditions for decision makers and ecological scholars, and the beneﬁts and costs of the projects should be considered comprehensively.


Introduction
In recent decades, climate change has attracted widespread attention [1,2]. Under the combined influence of natural and anthropogenic factors, the processes, structures and functions of terrestrial ecosystems have undergone significant changes. These changes will have an impact on the net primary productivity (NPP) of terrestrial vegetation, which is comprehensively manifested [3]. As an important vegetation parameter and ecological indicator, NPP can directly reflect the changes in the ecological environment and the level of the carbon budget [1,4]. Thus, determining the carbon sink or source potential of ecosystems under future climate projections is of critical importance [5].
A large number of studies since the 1990s reported that the annual total NPP of vegetation in China's terrestrial ecosystem was approximately 3.64 ± 0.57 PgC [6][7][8][9][10][11][12][13][14][15], and its spatial pattern decreased from southeast to northwest due to the influence of meteorological conditions and topography. Meteorological conditions such as precipitation and temperature were the main factors that determined the change trend of NPP [16,17]. Since the beginning of the 21st century, benefiting from climate warming and humidification, and the effective implementation of ecological engineering [18][19][20], the NPP in areas such as the Qinghai-Tibet Plateau [21][22][23][24][25] and the Loess Plateau [16,26] has significantly increased. However, due to differences in research methods, data sources and research periods, the above studies still lack a clear understanding of the spatiotemporal pattern and variation trend of vegetation NPP in China over the past 20 years.
Furthermore, the research on NPP restoration potential is still relatively weak due to the difficulty in obtaining the ecological parameters of zonal climax vegetation [27]. For example, some scholars have evaluated the status and potential of carbon sequestration in grassland [28], wetlands [29] and forests [30] in China on the basis of positioning experiments, literature or sample survey data. However, the research on the carbon sequestration potential of a single vegetation type at a point scale is difficult to extend to a regional scale. Zhang et al. [31] calculated the productivity potential of grassland based on precipitation, temperature and evapotranspiration, but the method focused on meteorological elements, and ignored the differences in zonal climax vegetation. Zhao et al. [26] used a similar habitat method on the Loess Plateau to predict the future vegetation coverage in each district. However, the classification standard of terrain and meteorological conditions in this method were greatly affected by subjective factors, and the restoration potential of each district adopted the same degree of vegetation coverage, ignoring the spatial continuity of water and heat conditions. Most areas of the Loess Plateau are subject to human disturbance, and the vegetation coverage in each district obtained by geostatistical methods was not the ideal climax vegetation coverage. Therefore, the current theoretical and technical system for research on the NPP restoration potential of large-scale vegetation is still imperfect, and it is difficult to support the differentiated implementation of ecological engineering. It was urgent to make new attempts in theory and technical methods. For a long time, the Chinese government has built a large number of nature reserves [32,33] to protect representative natural ecosystems, the purpose of which is to reduce human interference and protect the authenticity, integrity and stability of the ecosystem in the region and to restore the degraded ecosystem to an ecosystem with species composition, diversity and community structures that are close to the levels of zonal community. The nature reserves provide an ideal reference standard for obtaining the zonal ecosystem climax background.
What is the spatiotemporal pattern and changing trend of vegetation NPP in China in the past 20 years? How can the zonal climax background of vegetation be selected with the help of nature reserves, and how can the zonal climax background of China's vegetation NPP be evaluated at the national scale? On this basis, how large is the restoration status and restoration potential of vegetation NPP in China, and what are the spatial differences? To address these questions, we used the NPP data of MODIS (moderate-resolution imaging spectroradiometer) to analyse the change trend of NPP and its influencing factors. Additionally, the vegetation status of nature reserves was regarded as the zonal ecosystem climax background and used to evaluate the NPP restoration potential of China's forest, grassland and desert ecosystem vegetation. This approach has important theoretical and practical significance for research on the carbon cycle and carbon budget of terrestrial ecosystems under the background of global climate change and carbon neutrality and provides a scientific basis for the sustainable use of terrestrial ecosystems and the implementation of regional differentiation of ecological engineering.

Study Area
In this study, we focused on the dynamics of vegetation NPP in China. The land area of China is approximately 9,600,000 km 2 . The landscape of China varies significantly across its width, where the diverse landscape types and different hydrothermal conditions in the latitudinal, longitudinal and vertical terrains have formed a complex natural and geographical environment [34,35]. Geographically, the northwestern part of China is located in the hinterland of the Eurasian continent, the southeastern part of China faces the Pacific and the Qinghai-Tibet Plateau in southwestern China has some of the highest terrain on Earth. The distribution characteristics of major land covers in China showed significant spatial differences ( Figure 1). The Greater and Lesser Khingan Range in the northeast, the Changbai Mountains, the Qinling Mountains in the middle and the southeast, central and southwest mountainous and hilly areas are covered by forestland. The Northeast Plain, the North China Plain, the Jianghuai River basin and the Sichuan basin are covered by cultivated land. Eastern Inner Mongolia to the Qinghai-Tibet Plateau and northern Xinjiang are covered by grassland. Most of the northwestern area is covered by desert [36][37][38].
Remote Sens. 2022, 14, x FOR PEER REVIEW 3 of 28 ecosystem climax background and used to evaluate the NPP restoration potential of China's forest, grassland and desert ecosystem vegetation. This approach has important theoretical and practical significance for research on the carbon cycle and carbon budget of terrestrial ecosystems under the background of global climate change and carbon neutrality and provides a scientific basis for the sustainable use of terrestrial ecosystems and the implementation of regional differentiation of ecological engineering.

Study Area
In this study, we focused on the dynamics of vegetation NPP in China. The land area of China is approximately 9,600,000 km 2 . The landscape of China varies significantly across its width, where the diverse landscape types and different hydrothermal conditions in the latitudinal, longitudinal and vertical terrains have formed a complex natural and geographical environment [34,35]. Geographically, the northwestern part of China is located in the hinterland of the Eurasian continent, the southeastern part of China faces the Pacific and the Qinghai-Tibet Plateau in southwestern China has some of the highest terrain on Earth. The distribution characteristics of major land covers in China showed significant spatial differences ( Figure 1). The Greater and Lesser Khingan Range in the northeast, the Changbai Mountains, the Qinling Mountains in the middle and the southeast, central and southwest mountainous and hilly areas are covered by forestland. The Northeast Plain, the North China Plain, the Jianghuai River basin and the Sichuan basin are covered by cultivated land. Eastern Inner Mongolia to the Qinghai-Tibet Plateau and northern Xinjiang are covered by grassland. Most of the northwestern area is covered by desert [36][37][38].

Meteorological Data
The meteorological data came from the "China Surface Climate Data Daily Value Data Set (V3.0)" of the China Meteorological Data Network (http://data.cma.cn, accessed    [39], which contains daily temperature and precipitation data from 699 weather stations in China since January 1951 ( Figure 2). First, this study used the  [40][41][42][43] based on thin-slice spline theory to spatially interpolate the temperature and precipitation of the site and obtain the daily temperature and precipitation data at a resolution of 1 km across the country from 2000 to 2019. Second, the 0.25 • × 0.25 • daily precipitation grid data provided by the China Meteorological Data Network were used to control the quality of the interpolated precipitation data, and finally, the revised daily precipitation data with 1 km resolution across the country from 2000 to 2019 were obtained. Finally, the total annual precipitation (MTP) at a 1 km resolution was obtained by summing the corrected daily precipitation data at a resolution of 1 km. The average annual temperature (MAT) with a resolution of 1 km was obtained by averaging the interpolated daily temperature data with a resolution of 1 km.  [39], which contains daily temperature and precipitation data from 699 weather stations in China since January 1951 ( Figure 2). First, this study used the ANUSPLIN method [40][41][42][43] based on thin-slice spline theory to spatially interpolate the temperature and precipitation of the site and obtain the daily temperature and precipitation data at a resolution of 1 km across the country from 2000 to 2019. Second, the 0.25° × 0.25° daily precipitation grid data provided by the China Meteorological Data Network were used to control the quality of the interpolated precipitation data, and finally, the revised daily precipitation data with 1 km resolution across the country from 2000 to 2019 were obtained. Finally, the total annual precipitation (MTP) at a 1 km resolution was obtained by summing the corrected daily precipitation data at a resolution of 1 km. The average annual temperature (MAT) with a resolution of 1 km was obtained by averaging the interpolated daily temperature data with a resolution of 1 km.
(a) (b) Figure 2. The spatial distribution of China's total annual precipitation (a) and annual average temperature (b) in 2019.

The NPP Data of MODIS
The NPP data were derived from the MOD17A3HGF v006 product (https://lpdaac.usgs.gov/products/mod17a3hgfv006/, accessed on 24 December 2021) [44], which provides data on the annual value of vegetation NPP after 2000; the unit is kgC/m 2 /a, the numerical scaling factor is 0.0001, the time resolution is 1 year and the spatial resolution is 500 m. The data product level is L4 [45], and it has been processed by atmospheric correction, radiation correction, geometric correction and cloud removal. We used MRT tools and Python language to perform subset extraction, format conversion, projection transformation, image mosaicking, research area cropping and other processing on MOD17A3HGF v006 products and finally obtained the NPP data of MODIS for China from 2000 to 2019, which is called MODIS-NPP in this article.

NPP Verification Data
There are two types of NPP data used for MODIS-NPP accuracy verification ( Figure  3): (1) Grassland quadrat biomass data set. The data set includes 120 sample data collected by Fan [46,47] in the Three-River Headwaters Region (TRHR) from 2004 to 2006 and 48 grassland sample plot data collected in Inner Mongolia in 2010. The biomass data of grassland quadrats has been were converted into NPP data. The quadrat biomass data set was obtained by randomly setting five 1 m × 1 m small sample frames (alpine meadow 0.5 m × 0.5 m) during the period of maximum cultivation and growth, and collecting the plant sample in the square by the method of full rotation cutting. Information was collected on

The NPP Data of MODIS
The NPP data were derived from the MOD17A3HGF v006 product (https://lpdaac. usgs.gov/products/mod17a3hgfv006/, accessed on 24 December 2021) [44], which provides data on the annual value of vegetation NPP after 2000; the unit is kgC/m 2 /a, the numerical scaling factor is 0.0001, the time resolution is 1 year and the spatial resolution is 500 m. The data product level is L4 [45], and it has been processed by atmospheric correction, radiation correction, geometric correction and cloud removal. We used MRT tools and Python language to perform subset extraction, format conversion, projection transformation, image mosaicking, research area cropping and other processing on MOD17A3HGF v006 products and finally obtained the NPP data of MODIS for China from 2000 to 2019, which is called MODIS-NPP in this article.

NPP Verification Data
There are two types of NPP data used for MODIS-NPP accuracy verification ( Figure 3): (1) Grassland quadrat biomass data set. The data set includes 120 sample data collected by Fan [46,47] in the Three-River Headwaters Region (TRHR) from 2004 to 2006 and 48 grassland sample plot data collected in Inner Mongolia in 2010. The biomass data of grassland quadrats has been were converted into NPP data. The quadrat biomass data set was obtained by randomly setting five 1 m × 1 m small sample frames (alpine meadow 0.5 m × 0.5 m) during the period of maximum cultivation and growth, and collecting the plant sample in the square by the method of full rotation cutting. Information was collected on the latitude and longitude, topography, records, climate conditions, community characteristics, coverage of dominant species and other information. After drying for 48 h at 65 • C in a constant temperature drying oven, the dry weight (DM/m 2 ) of the plant Remote Sens. 2022, 14, 1634 5 of 26 sample was measured and recorded as the above-ground biomass of the sample. The underground biomass was sampled from the 10-30 cm deep soil layer in the sample plots, and then washed in clean water using a 0.3 mm mesh screen to separate the roots and other underground biomass from the soil. To minimise human interference, sample plots were usually selected on grasslands that were enclosed and were not grazed or mowed in the same year. (2) Reference data set for carbon cycle research of China's typical forest ecosystems [48]. This data set provides long-term carbon cycle benchmark observation data products of 10 typical forest ecosystems in China from 2005 to 2015, based on the CERN long-term dynamic monitoring database. These monitoring data were collected in keeping with CERN's protocols of observation and quality control [49,50]. There were occasional missing data in time-continuous meteorological observations; therefore, the data were processed by standardised gap filling [51]. The data set represents 80% of China's forest areas, including basic observation data sets such as biology, soil, atmosphere and moisture, as well as vegetation, soil carbon pools, productivity, respiration and carbon sinks.
Remote Sens. 2022, 14, x FOR PEER REVIEW 5 of 28 the latitude and longitude, topography, records, climate conditions, community characteristics, coverage of dominant species and other information. After drying for 48 h at 65 °C in a constant temperature drying oven, the dry weight (DM/m 2 ) of the plant sample was measured and recorded as the above-ground biomass of the sample. The underground biomass was sampled from the 10-30 cm deep soil layer in the sample plots, and then washed in clean water using a 0.3 mm mesh screen to separate the roots and other underground biomass from the soil. To minimise human interference, sample plots were usually selected on grasslands that were enclosed and were not grazed or mowed in the same year. (2) Reference data set for carbon cycle research of China's typical forest ecosystems [48]. This data set provides long-term carbon cycle benchmark observation data products of 10 typical forest ecosystems in China from 2005 to 2015, based on the CERN long-term dynamic monitoring database. These monitoring data were collected in keeping with CERN's protocols of observation and quality control [49,50]. There were occasional missing data in time-continuous meteorological observations; therefore, the data were processed by standardised gap filling [51]. The data set represents 80% of China's forest areas, including basic observation data sets such as biology, soil, atmosphere and moisture, as well as vegetation, soil carbon pools, productivity, respiration and carbon sinks.

Ecosystem Macro Structure Data
The 2015 ecosystem macro structure data ( Figure 3) came from the research results of Liu's research group at the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences [37]. The spatial distribution data of China's terrestrial ecosystem types were formed by the identification, research and classification of various ecosystem types based on the 1:100,000 scale land use/land cover data (LUCC) obtained after remote sensing interpretation [52,53]. The spatial resolution of these data is 1 km, and data are included for seven types of ecosystems, including farmland, forest, grassland, water and wetland, desert, settlement and other. The farmland ecosystem

Ecosystem Macro Structure Data
The 2015 ecosystem macro structure data ( Figure 3) came from the research results of Liu's research group at the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences [37]. The spatial distribution data of China's terrestrial ecosystem types were formed by the identification, research and classification of various ecosystem types based on the 1:100,000 scale land use/land cover data (LUCC) obtained after remote sensing interpretation [52,53]. The spatial resolution of these data is 1 km, and data are included for seven types of ecosystems, including farmland, forest, grassland, water and wetland, desert, settlement and other. The farmland ecosystem mainly includes paddy and dry fields. The forest ecosystem mainly includes forestland, shrubland, sparse woodland and other woodlands. The grassland ecosystem includes high-coverage grassland, medium-coverage grassland and low-coverage grassland. The water and wetland ecosystem includes canals, lakes, reservoir ponds, permanent glacier Remote Sens. 2022, 14, 1634 6 of 26 snow, tidal flats and beaches. The desert ecosystem includes sandy land, Gobi, saline land, alpine desert and tundra. The settlement ecosystem includes urban land, rural settlement and other construction land. The other ecosystems mainly include bare soil and bare rock texture.

Nature Reserve Data
China's nature reserves refer to areas where typical ecosystems, concentrated distribution areas of rare and endangered wild animals and plants, and natural relics of special significance are protected. The "2013-2016 National Nature Reserve Spatial Distribution Data Set" used in this study was derived from the special scientific and technological basic work "Data Integration and Standardization of Basic Science and Technology Work (2013FY110900)". The data set contains 1134 national, provincial or county-level nature reserves and 10 protection types (Table 1). It should be noted that many wild animal and wild plant nature reserves are distributed with forest, grassland and desert ecosystems, so these three ecosystems are effectively protected here. This research selected 741 national and provincial nature reserves (Figure 4), with a total area of 1,302,200 km 2 , and then used LUCC data to eliminate water bodies, glaciers, snow and artificial surfaces in nature reserves. The remaining part was used to simulate the vegetation NPP of the climax background of the forest, grassland and desert ecosystems in China. This study believes that the abovementioned nature reserves have undergone years of strict protection with little human disturbance, and the authenticity, integrity and stability of the forest, grassland and desert ecosystems in their areas have approached their zonal climax backgrounds.

Ecological Geographic Area Data
This study adopted the framework plan of China's ecogeographical zonal system formulated by Zheng Du et al. [54], that is, the three-level classification of China's ecogeographical zones, which includes 10 first-level zones, 21 second-level zones and 49 third-level zones. In the process of subregional regression analysis of climax background NPP and meteorological elements, two conditions needed to be met: (1) A certain number of nature reserves must be included in each ecogeographical zone; and (2) a more refined ecogeographical division was used to make the fitting analysis more accurate. Based on the above two conditions, we selected the 21 second-level zones of China's ecogeographical divisions as the ecogeographical zones of this study (Table 2). We adopted the principle of proximity or similar climate conditions to merge the zones in which nature reserves were not included into other zones.
There were 15 secondary zones after adjustment ( Figure 4). The adjusted zones and their adjustment bases were as follows: (1) The former IA and the IIA were merged because the area of the original IA zone was not large, the climax background in the IA zone was less distributed and the climate was not much different from the original IIA zone, so the two zones were merged. (2) Combining the original IIIA zone with the original IIIB zone, the original IIIA zone was a low mountain and hilly deciduous broad-leaved forest and an artificial forest area in Jiaodong-Liaodong. We merged the two zones because the original IIIA zone had less climax background distribution, and it was adjacent to the original IIIB zone. (3) The original HID zone and the HIID zone were merged. The original HID zone was small and surrounded by the HIID zone, so the HID zone and the HIID zone were merged. (4) The original VIIA zone and the VIIIA zone were merged with the VIA zone.
The areas of the original VIIA and VIIIA zones were small, and these three zones had less climax background distribution but had relatively small differences in climate, so we merged them into one zone. (5) The original IXA zone was an equatorial tropical humid area, mainly distributed in the Spratly Islands. This zone was not the key area of this study, so we did not include it.

Trend and Significance Test of NPP
The trend of NPP for many years was calculated using Sen's slope [55]. Sen's slo avoids the loss of time series data and the influence of data distribution on the analy results, and it can eliminate the interference of outliers on the time series [56]. It has sup rior performance in the trend analysis of time series. The calculation formula of Sen's slo is as follows: where S is Sen's slope, and and are the sequence values at times j and i, respective If S > 0, the time series data show an upwards trend; otherwise, they show a downwar trend.
Sen's slope is obtained by calculating the median value of the sequence. It can redu the noise interference very well, but it cannot realise the significance test of the sequen trend by itself [57]. The Mann-Kendall (MK) trend test method is a nonparametric te method recommended by the World Meteorological Organization (WMO) and has be widely used [58][59][60][61]. The advantage of using the MK test is that the sample does not ne to obey a certain distribution, and it is not sensitive to missing values and outliers; thu the MK method was introduced to test the significance of the trend of long-term NPP da [59,62]. The statistical test method is as follows: For the time series , = 1, 2, 3, … i, …j, …, n, define the standardised test statis :

Trend and Significance Test of NPP
The trend of NPP for many years was calculated using Sen's slope [55]. Sen's slope avoids the loss of time series data and the influence of data distribution on the analysis results, and it can eliminate the interference of outliers on the time series [56]. It has superior performance in the trend analysis of time series. The calculation formula of Sen's slope is as follows: where S is Sen's slope, and x j and x i are the sequence values at times j and i, respectively. If S > 0, the time series data show an upwards trend; otherwise, they show a downwards trend. Sen's slope is obtained by calculating the median value of the sequence. It can reduce the noise interference very well, but it cannot realise the significance test of the sequence trend by itself [57]. The Mann-Kendall (MK) trend test method is a nonparametric test method recommended by the World Meteorological Organization (WMO) and has been widely used [58][59][60][61]. The advantage of using the MK test is that the sample does not need to obey a certain distribution, and it is not sensitive to missing values and outliers; thus, the MK method was introduced to test the significance of the trend of long-term NPP data [59,62]. The statistical test method is as follows: For the time series X k , k = 1, 2, 3, . . . i, . . . j, . . . , n, define the standardised test statistic Z: where, x j and x i are the sequence values at time j and i, respectively. n is the number of data. When n ≥ 8, the test statistic S is approximately normally distributed, and its mean and variance are as follows: Given the significance level α, if |Z| > Z 1−α/2 , the hypothesis indicating that the current time series data does not have a trend is rejected, and there is an obvious trend change; Z 1−α/2 is the value corresponding to the standard normal function distribution table at the significance level α. This research considered the significance levels α = 0.1 and α = 0.05 for the significance test. When |Z| is greater than 1.65 or 1.96, the change trend of the data series passes the significance test with a confidence of 90% or 95%, respectively. Additionally, combined with the positive and negative Sen trends, the significance test results were defined as five types: Significant decrease, relatively significant decrease, no significant change, relatively significant increase and significant increase (Table 3).

Climax Background Vegetation NPP Simulation of Forest, Grassland and Desert Ecosystems
To protect representative natural ecosystems, the Chinese government has built a large number of nature reserves with the purpose of reducing human interference; protecting the authenticity, integrity and stability of the ecosystem in the zones; and restoring the Remote Sens. 2022, 14, 1634 9 of 26 degraded ecosystem climax background ecosystems in which species compositions, diversities and community structures are close to those of the zonal communities. China's national and provincial nature reserves can effectively reduce human interference and protect the authenticity, integrity and stability of the excellent ecosystems in the region. They can also restore degraded ecosystems to the climax background ecosystems in which the species compositions, diversities and community structures are close to those of the zonal communities. Therefore, the vegetation in nature reserves with forests, grasslands and deserts as the protection objects selected in this research is suitable as the climax background. The zonal climax background in this study refers to a zonal community on the surface of the earth that is compatible with meteorological conditions and can fully reflect the climate characteristics of a region. While the vegetation NPP of artificial oasis, irrigated agricultural area and artificial forestation area is relatively high, it does not meet the concept of the climax background of our research. Furthermore, to make the simulation results consistent with the actual natural environment, the simulations of sandy land, permanent snow glaciers and water bodies simulations eliminated in this research.
There are two ways to extend the climax background NPP of the local scale to the regional scale: (1) The direct assignment method. That is, the climax background value of the local scale is directly assigned to the larger-scale ecogeographic area. The advantage of this method is that it is relatively simple. The disadvantage is that each ecogeographic area uses only one value, thus ignoring the spatial differences in the climate and other factors at the cell scale, which results in large errors in the research at the pixel scale. (2) Regression analysis method for meteorological conditions. This method assumes that the various types of ecosystems in the nature reserve are less affected by human interference, and the vegetation NPP can be regarded as only affected by natural factors such as climate. It can be divided into ecological geographic zones and ecosystem types. The vegetation NPP, temperature, precipitation and other natural elements of each zone can be regression analysed at the pixel scale, and the regression relationship can be calculated for a larger area to obtain a larger area climax background NPP.
The meteorological conditions regression analysis method was used to obtain the vegetation NPP of the climax background of the forest, grass and desert ecosystems in China. The specific process mainly includes the following three stages: (1) Multiyear mean meteorological and NPP data preparation. (1) (Table 4), and each equation obtained three parameters a ij , b ij and c ij , where i and j represent the parameters of ecosystem type j in the ecogeographical zone i. (4) According to the MAT avg and MTP avg of China, the regression equations obtained by stage 3 were used to calculate the climax background vegetation NPP of China's forest, grass and desert ecosystems. NPP was calculated as follows: where NPP ij is the climax background vegetation NPP (gC/m 2 /a) of the ecosystem type j in the ecogeographical zone i; MAT avg is the average temperature of China ( • C); and MTP avg is the average precipitation of China (mm). a ij , b ij and c ij are the parameters of the linear regression equations of the ecosystem type j in the ecogeographical zone i. The obtained fitting equation, R 2 and significance values are shown in the Table 4.
For deserts in HIIA-B, forests in IIIC, forests and grasslands in IVA and grasslands in VIA, the fitting equations were not significant (p > 0.5) because the number of nature reserves for these ecosystem types in the abovementioned zones was relatively small. For this reason, this research used the F test to test the significance of the fitting equation and found that the model significance satisfies p < 0.05, so we still used the fitting equations to calculate the climax background NPP of China's forest, grass and desert ecosystems.
Theoretically, the climax background NPP as the expected target value of ecological restoration should be no worse than the data of all historical years. The maximum value of each year also objectively showed the target restoration value that the current pixel could achieve, so this pixel could be used. The maximum NPP in historical years was regarded as the climax background NPP. On this basis, the maximum value of NPP in historical years was used to correct the climax background NPP.

Restoration Potential of NPP Calculation of Forest, Grassland and Desert Ecosystems
The restoration potential is the difference between the restoration status of the NPP of the ecosystem vegetation and the NPP of the zonal climax background in the same ecosystem. The calculation formula is as follows: where NPP p is the restoration potential of vegetation NPP; NPP Top is the climax background NPP; and NPP cur is the restoration status of vegetation NPP. Considering that the vegetation NPP data in a single year are greatly affected by climate fluctuations and cannot accurately reflect the restoration status of the ecosystem, this paper used the 3-year average value of the vegetation NPP of the forest, grassland and desert ecosystems from 2017 to 2019 as the restoration status. Furthermore, to compare with the climax vegetation NPP, we eliminated the restoration status NPP data of sandy land, permanent snow glaciers and water bodies.
The word "potential" is usually regarded as a commendatory word; that is, the greater the potential is, the better the result. However, in this study, the restoration potential of vegetation is a neutral term that represents the gap between the current vegetation restoration status and the climax background of the natural zones.

Multiple Regression Residual Analysis
Residual analysis analyses the reliability, periodicity or other interference of the data through the information provided by the residual. The residual refers to the difference between the observed value and the predicted value. The multiple regression residual analysis method can be used to study the influence and relative contribution of meteorological conditions and human activities on vegetation change [63][64][65]. The main steps are as follows: (1) Based on the national MODIS-NPP (NPP obs ), MAT and MTP time series data from 2000 to 2019, with NPP obs as the dependent variable and MAT and MTP as the independent variables to fit a binary linear regression equation pixel by pixel, each equation obtained three parameters a mn , b mn and c mn . m and n represent the parameters of the mth row and nth column in the raster data. (2) Based on the MAT and MTP time series data and regression model, the vegetation NPP predicted value (NPP cc ) was calculated. (3) The difference between NPP obs and NPP cc , that is, the residual of vegetation NPP (NPP ha ), was calculated.
NPP ha kmn = NPP obs kmn − NPP cc kmn (10) where NPP cc kmn is the predicted value of vegetation NPP of the mth row and nth column in the kth year, and the unit is gC/m 2 /a. NPP obs kmn is the MODIS-NPP observation value of the mth row and nth column in the kth year, and the unit is gC/m 2 /a. a mn , b mn and c mn are model parameters of the mth row and nth column in the raster data. MAT kmn and MTP kmn are the average annual temperature and total annual precipitation of the mth row and nth column in the kth year, respectively. NPP ha kmn is the residual of the mth row and nth column in the kth year. The ratio of the trend of NPP cc kmn and NPP ha kmn to the trends of NPP obs kmn was used to characterise the contribution rates of meteorological conditions and human activities to vegetation NPP changes, respectively.

Uncertainty Estimation Approaches
In this paper, we analysed the uncertainty of the climax background vegetation NPP according to the uncertainty of the climax vegetation NPP simulation equation (Table 4). The method replaced the upper and lower limits (significance level 0.05) of the regression equation parameters a ij , b ij and c ij in Table 4 with the a ij , b ij and c ij parameters in Equation (7), respectively, and simulated the climax background vegetation NPP. The upper and lower limits of the climax vegetation NPP were obtained, that is, the uncertainty of the climax vegetation NPP.
The restoration potential is the difference between the restoration status of the NPP of the ecosystem vegetation and the zonal climax background vegetation NPP (Equation (8)). The restoration status was obtained by averaging the annual MODIS-NPP values over many years, and its uncertainty did not need to be calculated, so the uncertainty of restoration potential was determined by the uncertainty of the climax background vegetation NPP according to the principle of error conduction.

MODIS-NPP Data Accuracy Verification
The accuracy verification of vegetation NPP simulation results is one of the main difficulties faced by current regional NPP estimation model research. We used the measured NPP data to verify MODIS-NPP data [13,17]. First, according to the location information of each station (or sample), we extracted the MODIS-NPP pixel values at the same location and the same year at the station (or sample), and finally obtained 270 station-years with the measured NPP data and MODIS-NPP pairs. Then, considering the impact of the flux tower carbon footprint, this study used the MODIS-NPP average value in a 3 × 3 pixel window around the station (or sample) for comparison and verification with the measured NPP data [5,[66][67][68]. Finally, the Python statistical analysis toolkit was used to perform statistical analysis on the measured NPP data and MODIS-NPP data pairs. The results showed that the trend line of MODIS-NPP and the measured NPP was near the 1:1 line, and the coefficient of determination R 2 was 0.75. The MODIS-NPP data had high simulation accuracy and high consistency with the measured NPP data, which was suitable for NPP research on China's terrestrial ecosystem.
The MODIS-NPP value of the Xishuangbanna site ( Figure 5 in the green circle) was greater than 1400 gC/m 2 /a, which was higher than that of the other sites. The vegetation type of the Xishuangbanna site is tropical monsoon forest, the forest age is 200 year and the dominant tree species are evergreen trees, including Pometia pinnata and Terminalia myriocarpa. Because this site has better climate conditions, more green plant biomass and stronger photosynthesis, the MODIS-NPP value is higher. However, MODIS-NPP will be somewhat overestimated in grassland and forest by using a t test. The accuracy verification of vegetation NPP simulation results is one of the main difficulties faced by current regional NPP estimation model research. We used the measured NPP data to verify MODIS-NPP data [13,17]. First, according to the location information of each station (or sample), we extracted the MODIS-NPP pixel values at the same location and the same year at the station (or sample), and finally obtained 270 station-years with the measured NPP data and MODIS-NPP pairs. Then, considering the impact of the flux tower carbon footprint, this study used the MODIS-NPP average value in a 3 × 3 pixel window around the station (or sample) for comparison and verification with the measured NPP data [5,[66][67][68]. Finally, the Python statistical analysis toolkit was used to perform statistical analysis on the measured NPP data and MODIS-NPP data pairs. The results showed that the trend line of MODIS-NPP and the measured NPP was near the 1:1 line, and the coefficient of determination R 2 was 0.75. The MODIS-NPP data had high simulation accuracy and high consistency with the measured NPP data, which was suitable for NPP research on China's terrestrial ecosystem.
The MODIS-NPP value of the Xishuangbanna site ( Figure 5 in the green circle) was greater than 1400 gC/m 2 /a, which was higher than that of the other sites. The vegetation type of the Xishuangbanna site is tropical monsoon forest, the forest age is 200 year and the dominant tree species are evergreen trees, including Pometia pinnata and Terminalia myriocarpa. Because this site has better climate conditions, more green plant biomass and stronger photosynthesis, the MODIS-NPP value is higher. However, MODIS-NPP will be somewhat overestimated in grassland and forest by using a t test.

Spatial Patterns of Vegetation NPP
The spatial patterns of vegetation NPP in China had obvious spatial distribution differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 6), and this trend was related to the spatial distributions of meteorological conditions, topography, soil and vegetation [13,69]. The water and heat conditions in southern China were relatively good, forest vegetation accounted for a relatively large proportion of vegetation and the vegetation coverage was relatively high; thus the NPP in southern China was relatively high. For example, the NPP of the southern Qinghai-Tibet Plateau, Yunnan and Guangxi was greater than 500 gC/m 2 /a, while the NPP was relatively low in northwestern China because northwestern China is mainly composed of grassland and deserts, so vegetation growth is obviously restricted by precipitation. For example, vegetation NPP was less than 200 gC/m 2 /a in most areas of southwestern Qinghai, southwestern Tibet and central and western Inner Mongolia. The statistical results show that the area with vegetation NPP between 200 and 500 gC/m 2 /a accounted for the largest proportion (Table 5)

Change Trend of Vegetation NPP
The annual total NPP of vegetation in China's terrestrial ecosystems varied between 2.72 and 3.29 PgC from 2000 to 2019 (Figure 7), with a multiyear average of 3.09 PgC, which accounted for approximately 5.15% of the total global NPP of 60 PgC in the IPCC 6th Assessment Report [70]. In the past 20 years, the annual total NPP of China's terrestrial ecosystem vegetation has increased by approximately 15.86%, or 0.52 PgC. The annual total NPP of China's terrestrial ecosystem vegetation fluctuated significantly from 2000 to 2019 (p < 0.01), and the trend was 0.022 PgC/a. In the first ten years of the study period, the trend was 0.03 PgC/a, and the increasing trend was significant (p < 0.01); in the next ten years, the trend was 0.023 PgC/a, and the increasing trend was significant (p < 0.01).
The highest value appeared in 2018, and the annual total NPP was 3.29 PgC. The lowest value appeared in 2000, and the annual total NPP was 2.73 PgC. Since the reform and opening up, due to the accelerated development of industrialization, urbanization and agricultural modernization and the irrational or excessive utilization of resources, the ecosystem has been severely degraded, so the NPP in 2000 was relatively lower than that in other years. Since the beginning of the 21st century, benefiting from climate warming and humidification, and the effective implementation of ecological engineering [18][19][20], which have effectively protected and restored vegetation in some project areas, the NPP in 2018 was higher than that in other years.

Change Trend of Vegetation NPP
The annual total NPP of vegetation in China's terrestrial ecosystems varied between 2.72 and 3.29 PgC from 2000 to 2019 (Figure 7), with a multiyear average of 3.09 PgC, which accounted for approximately 5.15% of the total global NPP of 60 PgC in the IPCC 6th Assessment Report [70]. In the past 20 years, the annual total NPP of China's terrestrial ecosystem vegetation has increased by approximately 15.86%, or 0.52 PgC. The annual total NPP of China's terrestrial ecosystem vegetation fluctuated significantly from 2000 to 2019 (p < 0.01), and the trend was 0.022 PgC/a. In the first ten years of the study period, the trend was 0.03 PgC/a, and the increasing trend was significant (p < 0.01); in the next ten years, the trend was 0.023 PgC/a, and the increasing trend was significant (p < 0.01). From the perspective of the trend of vegetation NPP (Figure 8a), the average trend of vegetation NPP in China from 2000 to 2019 was 2.39 gC/m 2 /a 2 . The trend of vegetation NPP in most areas was positive, while the area with a negative trend accounted for 10.21%, and the area with no obvious change accounted for 21.69%. This result was consistent with Zhao's report that the national NPP from 2001 to 2015 appeared as a growth trend, accounting for 65.29% of the area [15]. The places with a change trend of more than 5 gC/m 2 /a 2 were mainly distributed in the agro-pastoral zone, the eastern Sichuan Basin and the Greater Khingan Range. The average increasing trend in the Loess Plateau was 7.12 gC/m 2 /a 2 . The places with a change trend of less than −5 gC/m 2 /a 2 were mainly distributed in the southern Qinghai-Tibet Plateau and Guangdong, Fujian, Yunnan and other places in the south. The decreasing trend in southern Tibet was below −10 gC/m 2 /a 2 . The NPP change trend in most areas of the Qinghai-Tibet Plateau was between 0 and 3 gC/m 2 /a 2 . The NPP in the vast area of the northwest showed a certain decreasing trend, between −3 and 0 gC/m 2 /a 2 .
The MK method was used to test the significance of the trend of vegetation NPP across the country from 2000 to 2019 on a pixel scale (Figure 8b). The spatial distribution pattern of significance was relatively consistent with the spatial distribution of change trends. (1) The area with significant increases included the Loess Plateau, which was 274.69 × 10 4 km 2 (Table 6)  The highest value appeared in 2018, and the annual total NPP was 3.29 PgC. The lowest value appeared in 2000, and the annual total NPP was 2.73 PgC. Since the reform and opening up, due to the accelerated development of industrialization, urbanization and agricultural modernization and the irrational or excessive utilization of resources, the ecosystem has been severely degraded, so the NPP in 2000 was relatively lower than that in other years. Since the beginning of the 21st century, benefiting from climate warming and humidification, and the effective implementation of ecological engineering [18][19][20], which have effectively protected and restored vegetation in some project areas, the NPP in 2018 was higher than that in other years.
From the perspective of the trend of vegetation NPP (Figure 8a), the average trend of vegetation NPP in China from 2000 to 2019 was 2.39 gC/m 2 /a 2 . The trend of vegetation NPP in most areas was positive, while the area with a negative trend accounted for 10.21%, and the area with no obvious change accounted for 21.69%. This result was consistent with Zhao's report that the national NPP from 2001 to 2015 appeared as a growth trend, accounting for 65.29% of the area [15]. The places with a change trend of more than 5 gC/m 2 /a 2 were mainly distributed in the agro-pastoral zone, the eastern Sichuan Basin and the Greater Khingan Range. The average increasing trend in the Loess Plateau was 7.12 gC/m 2 /a 2 . The places with a change trend of less than −5 gC/m 2 /a 2 were mainly distributed in the southern Qinghai-Tibet Plateau and Guangdong, Fujian, Yunnan and other places in the south. The decreasing trend in southern Tibet was below −10 gC/m 2 /a 2 . The NPP change trend in most areas of the Qinghai-Tibet Plateau was between 0 and 3 gC/m 2 /a 2 . The NPP in the vast area of the northwest showed a certain decreasing trend, between −3 and 0 gC/m 2 /a 2 . Remote Sens. 2022, 14

Climax Background NPP of Forest, Grassland and Desert Ecosystems
The climax background NPP of the forest, grassland and desert ecosystems across the country showed large spatial differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients ( Figure 9). The NPP was approximately 545.80 ± 55 gC/m 2 /a, and the annual total NPP was approximately 2.76 ± 0.28 PgC. The area of 500-1000 gC/m 2 /a accounted for the largest proportion, which was 33.52% (Table 7). Among them, tropical areas, northeastern regions and most areas of southern Tibet had the largest climax background NPP due to their extensive distribution of forests, shrubs and high-coverage grasslands and the vegetation NPP in some areas reached above 1000 gC/m 2 /a. However, the climax background NPP in most areas of the western Inner Mongolia Plateau and the northwestern Qinghai-Tibet Plateau was relatively small, generally below 200 gC/m 2 /a. According to ecosystem types, the climax background NPP per unit area of the forest ecosystem in China was approximately 849.36 ± 88.32 gC/m 2 /a, and the annual total NPP was approximately 1.90 ± 0.2 PgC. The climax background NPP per unit area of the grassland ecosystem was approximately 306.75 ± 28.14 gC/m 2 /a, the annual total NPP is approximately 0.80 ± 0.07 PgC. After eliminating the sandy and other non-vegetation areas, the NPP per unit area of the desert ecosystem was approximately 67.56 ± 4.16 gC/m 2 /a and the annual total NPP was approximately 0.009 ± 0.0005 PgC. The MK method was used to test the significance of the trend of vegetation NPP across the country from 2000 to 2019 on a pixel scale (Figure 8b). The spatial distribution pattern of significance was relatively consistent with the spatial distribution of change trends. (1) The area with significant increases included the Loess Plateau, which was 274.69 × 10 4 km 2 (Table 6), and included the agro-pastoral zone, the eastern Sichuan Basin, the Loess Plateau and the Greater Khingan Range. (2) The areas with relatively significant increases were mainly distributed around places with significant increases, such as the Qinghai-Tibet Plateau, Inner Mongolia Plateau and North China Plain. (3) The area with a significant decrease was 6.47 × 10 4 km 2 , and was found in southern Tibet. (4) The area with no significant change was 457.28 × 10 4 km 2 , in northwestern China and the western Qinghai-Tibet Plateau. (5) The area with a relatively significant increase was 166.06 × 10 4 km 2 . The area with a relatively significant decrease was 16.56 × 10 4 km 2 , scattered in southern China, including Guangdong, Fujian, Yunnan and other places.

Climax Background NPP of Forest, Grassland and Desert Ecosystems
The climax background NPP of the forest, grassland and desert ecosystems across the country showed large spatial differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 9). The NPP was approximately 545.80 ± 55 gC/m 2 /a, and the annual total NPP was approximately 2.76 ± 0.28 PgC. The area of 500-1000 gC/m 2 /a accounted for the largest proportion, which was 33.52% (Table 7). Among them, tropical areas, northeastern regions and most areas of southern Tibet had the largest climax background NPP due to their extensive distribution of forests, shrubs and high-coverage grasslands and the vegetation NPP in some areas reached above 1000 gC/m 2 /a. However, the climax background NPP in most areas of the western Inner Mongolia Plateau and the northwestern Qinghai-Tibet Plateau was relatively small, generally below 200 gC/m 2 /a.

Restoration Status and Restoration Potential NPP of Forest, Grassland and Desert Ecosystems
The spatial patterns of restoration status NPP in China had obvious spatial distribution differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 10a). The restoration status NPP of the national forest, grassland and desert ecosystems was 441.84 gC/m 2 /a, and the annual total NPP was 2.24 PgC. The area of 500-1000 gC/m 2 /a accounted for the largest proportion, which was 32.20%. The area of 0-2 gC/m 2 /a accounted for 9.95% (Table 8), which was mainly distributed in the desert ecosystem in the northwestern China.
Among them, the restoration status NPP per unit area of China's forest ecosystem was approximately 689.69 gC/m 2 /a, and the annual total NPP was approximately 1.54 PgC. The restoration status NPP per unit area of the grassland ecosystem was approximately 246.42 gC/m 2 /a, and the annual total NPP was approximately 0.65 PgC. The restoration status NPP per unit area of desert ecosystem was approximately 51.38 gC/m 2 /a, and the annual total NPP was approximately 0.007 PgC. The ecological reconstruction projects with ecological protection and vegetation restoration as the main goal have achieved  According to ecosystem types, the climax background NPP per unit area of the forest ecosystem in China was approximately 849.36 ± 88.32 gC/m 2 /a, and the annual total NPP was approximately 1.90 ± 0.2 PgC. The climax background NPP per unit area of the grassland ecosystem was approximately 306.75 ± 28.14 gC/m 2 /a, the annual total NPP is approximately 0.80 ± 0.07 PgC. After eliminating the sandy and other non-vegetation areas, the NPP per unit area of the desert ecosystem was approximately 67.56 ± 4.16 gC/m 2 /a and the annual total NPP was approximately 0.009 ± 0.0005 PgC.

Restoration Status and Restoration Potential NPP of Forest, Grassland and Desert Ecosystems
The spatial patterns of restoration status NPP in China had obvious spatial distribution differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 10a). The restoration status NPP of the national forest, grassland and desert ecosystems was 441.84 gC/m 2 /a, and the annual total NPP was 2.24 PgC. The area of 500-1000 gC/m 2 /a accounted for the largest proportion, which was 32.20%. The area of 0-2 gC/m 2 /a accounted for 9.95% (Table 8), which was mainly distributed in the desert ecosystem in the northwestern China. remarkable results since 2000 [21][22][23][24][25], and the ecosystems on the Loess Plateau and in other places have been effectively restored.
(a) (b) Figure 10. Restoration status (a) and restoration potential (b) of NPP in China's forest, grassland and desert ecosystems.
The spatial patterns of restoration potential NPP in China had obvious spatial distribution differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 10b). The restoration potential NPP of the national forest, grassland and desert ecosystem types was approximately 103.96 ± 55 gC/m 2 /a, and the annual total restoration potential NPP was approximately 0.52 ± 0.28 PgC, which accounted for approximately 19.05% of the annual total climax background NPP. The area of 50-100 gC/m 2 /a accounted for the largest proportion, which was 24.93% (Table 8).
Among them, the restoration potential NPP in China's forest ecosystem was approximately 159.67 ± 88.32 gC/m 2 /a, and the annual total vegetation NPP was approximately 0.36 ± 0.2 PgC, accounting for approximately 18.80% of the annual total climax background NPP. The restoration potential NPP in the grassland ecosystem was approximately 60.33 ± 28.14 gC/m 2 /a, and the annual total vegetation NPP was approximately 0.16 ± 0.07 PgC, accounting for approximately 19.67% of the total annual climax background NPP. The restoration potential NPP of the desert ecosystem was approximately 16.18 ± 4.16 gC/m 2 /a. The annual total vegetation NPP was approximately 0.002 ± 0.0005 PgC, accounting for approximately 23.95% of the annual total climax background NPP.   Among them, the restoration status NPP per unit area of China's forest ecosystem was approximately 689.69 gC/m 2 /a, and the annual total NPP was approximately 1.54 PgC. The restoration status NPP per unit area of the grassland ecosystem was approximately 246.42 gC/m 2 /a, and the annual total NPP was approximately 0.65 PgC. The restoration status NPP per unit area of desert ecosystem was approximately 51.38 gC/m 2 /a, and the annual total NPP was approximately 0.007 PgC. The ecological reconstruction projects with ecological protection and vegetation restoration as the main goal have achieved remarkable results since 2000 [21][22][23][24][25], and the ecosystems on the Loess Plateau and in other places have been effectively restored.
The spatial patterns of restoration potential NPP in China had obvious spatial distribution differences. Generally, there was a decreasing trend from southeast to northwest along the water and heat gradients (Figure 10b). The restoration potential NPP of the national forest, grassland and desert ecosystem types was approximately 103.96 ± 55 gC/m 2 /a, and the annual total restoration potential NPP was approximately 0.52 ± 0.28 PgC, which accounted for approximately 19.05% of the annual total climax background NPP. The area of 50-100 gC/m 2 /a accounted for the largest proportion, which was 24.93% (Table 8).
Among them, the restoration potential NPP in China's forest ecosystem was approximately 159.67 ± 88.32 gC/m 2 /a, and the annual total vegetation NPP was ap-proximately 0.36 ± 0.2 PgC, accounting for approximately 18.80% of the annual total climax background NPP. The restoration potential NPP in the grassland ecosystem was approximately 60.33 ± 28.14 gC/m 2 /a, and the annual total vegetation NPP was approximately 0.16 ± 0.07 PgC, accounting for approximately 19.67% of the total annual climax background NPP. The restoration potential NPP of the desert ecosystem was approximately 16.18 ± 4.16 gC/m 2 /a. The annual total vegetation NPP was approximately 0.002 ± 0.0005 PgC, accounting for approximately 23.95% of the annual total climax background NPP.

Methodological Approaches and Uncertainty
The results obtained by the regression analysis method of meteorological conditions based on nature reserves reflected the zonal climax vegetation status to a large extent, and the obtained values had a smooth transition within each eco-geographical division and between each eco-geographical division, which truly reflected the law of gradual change in the climate, vegetation and natural conditions. However, there was still some uncertainty associated with the NPP simulation of climax background vegetation.
(1) The MODIS NPP product, partly due to the uncertainty in the default Biome-specified Parameters Look-Up Table (BPLUT) of the MODIS photosynthesis (PSN) model [71]. Parameter optimization is a promising method that can be used to calibrate uncertain parameters of the carbon cycle model. (2) Choice of zonal climax background. In the core areas of some nature reserves, such as LiboMaolan Nature Reserve, Jinggangshan Nature Reserve and other places, due to inconvenient transportation and small population, after years of strict protection, it can be considered that the vegetation condition is close to its zonal climax background. However, vegetation in nature reserves is not always in a climax state. For example, some nature reserves do not protect good natural ecosystems, only the ecological location may be important, but the reserve could still have damaged zones. Additionally, some nature reserves still have a certain degree of human disturbance due to inadequate management and other reasons. However, in nature reserves where overmature forests exist and in grassland nature reserves that have been enclosed for many years, the vegetation NPP is not maximised due to slowed or discontinued forest growth, or a lack of appropriate human disturbance [72,73]. Therefore, some indicators, such as biomass, can be added to characterise climax vegetation in future research [74]. (3) Uncertainty of the simulation method. The regression analysis of climax vegetation status and natural factors is not only closely related to meteorological conditions such as precipitation and temperature, but also related to natural factors such as soil, vegetation and topography. Generally, adding the above natural factors for multiple regression will obtain a more accurate fitting equation, but when there are more natural factor combinations, the number of corresponding nature reserves will be fewer, which will affect the accuracy of the fitting equation. Therefore, choosing between these natural factors requires a weighing process.

Change Trend of Vegetation NPP and Its Classification
The trend of vegetation NPP in most areas of China was positive. The main factors affecting this change trend of vegetation NPP included meteorological conditions and human activities on a large scale [1,64,75]. Many studies have shown that China has experienced dramatic climate change characterised by warming and humidification in the past 20 years [76,77] (Figure 11). For example, the temperature and precipitation in the Inner Mongolia Plateau, Loess Plateau, Northeast Plain and Sichuan Basin showed an increasing trend, which promote the increase of NPP. On the Qinghai-Tibet Plateau, the temperature showed an increasing trend, which promoted the increase in NPP. Benefiting from climate warming and humidification and the effective implementation of ecological engineering [18][19][20], the NPP in areas such as the Qinghai-Tibet Plateau [21][22][23][24][25] and the Loess Plateau [16,26] has significantly increased. Satellite data showed an increasing leaf area of vegetation due to direct factors (e.g., human land-use management) and indirect factors (e.g., climate change, CO 2 fertilization, nitrogen deposition and recovery from natural disturbances). Among these, climate change and CO 2 fertilization effects seemed to be the dominant drivers [78,79]. However, the significance of change trends and impact factors were different in different regions, and the main features of the different classifications were as follows: (1) The areas with a significant increase in vegetation NPP were mainly distributed in the agro-pastoral zone, the Loess Plateau, the eastern Sichuan Basin and the Greater Khingan Range (p < 0.05). Ecological projects such as tree planting and afforestation in the abovementioned areas were the leading factors affecting the increase in NPP [80]. (2) The areas with a relatively significant decrease in vegetation NPP were mainly distributed in the Qinghai-Tibet Plateau and Inner Mongolia Plateau. There were many grassland and desert ecosystems in these places. While ecological projects such as the Returning Rangeland to Grassland and the Grain for Green Program were also implemented, the change trend was not as obvious as that in the agro-pastoral zone due to natural factors such as high cold and drought. (3) The areas with a significant decrease in vegetation NPP were mainly distributed on the southern Qinghai-Tibet Plateau. The decrease in precipitation was also an important reason affecting the decrease in NPP. Human activities such as deforestation are more frequent in southern Tibet and have caused the extensive destruction of forests and other vegetation. (4) The Guangdong, Fujian, Yunnan and other places in the south (p < 0.05), and the decreasing trend in southern Tibet was the most significant; the decreasing trend in most areas was below −10 gC/m 2 /a 2 . Due to the rapid development of cities and industries in Guangdong, southern Jiangsu, Hainan Province and other places [17], the increase in atmospheric aerosol concentration has led to a decrease in solar radiation, which has led to a decrease in the effective photosynthetic radiation of vegetation [81], which was caused by a relatively significant decrease in vegetation NPP. On the Yunnan-Guizhou Plateau, the decrease in precipitation was also an important reason for the decrease in NPP (Figure 11b). (5) The areas that failed the significance test were mainly distributed in most areas of the western Qinghai-Tibet Plateau, southern Xinjiang, western Inner Mongolia and other places. The Gobi and deserts in these places were widely distributed, the surface was bare or the vegetation coverage was low and the vegetation NPP was almost 0, so there was no obvious trend for vegetation NPP. Additionally, some places, such as the central part of the North China Plain, were mainly farmland, and the types of crops were relatively fixed, so there was no obvious trend for vegetation NPP.
We used a multiple regression residual analysis method to determine the contribution rates of climate factors and human activities to NPP. We found that the contribution rate of meteorological conditions to the change in vegetation NPP from 2000 to 2019 was 85.41%. The contribution rate of human activities to the change in vegetation NPP is 14.59%, so precipitation and temperature were the main controlling factors that determined the size and change in vegetation NPP. This result was more consistent with the research results of Hou et al. [12,16,17]. The results showed that it was reasonable to use the two meteorological conditions of precipitation and temperature to simulate the climax background NPP. Additionally, human activities had a relatively large contribution rate on the Loess Plateau, in the Three-North Shelter Forest Program Area and in the Returning Rangeland to Grassland Area, mainly through the effective implementation of large-scale afforestation, returning of grazing to grass and the balance of grass and livestock. The changes in meteorological conditions such as precipitation and temperature were largely caused by human activities [70,82]. How to clarify the impact of human activities on the climate system is the focus of research on the driving forces of changes in terrestrial ecosystems. Difficulties regarding the direction of this research remain and need to be further improved.
Yunnan-Guizhou Plateau, the decrease in precipitation was also an important reason for the decrease in NPP (Figure 11b). (5) The areas that failed the significance test were mainly distributed in most areas of the western Qinghai-Tibet Plateau, southern Xinjiang, western Inner Mongolia and other places. The Gobi and deserts in these places were widely distributed, the surface was bare or the vegetation coverage was low and the vegetation NPP was almost 0, so there was no obvious trend for vegetation NPP. Additionally, some places, such as the central part of the North China Plain, were mainly farmland, and the types of crops were relatively fixed, so there was no obvious trend for vegetation NPP.

Climax Background Vegetation NPP of Forest, Grassland and Desert Ecosystems
The spatial distribution of climax background vegetation NPP of forest, grassland and desert ecosystems in China has obvious regional characteristics. Generally, it decreases from southeast to northwest along the water and heat gradients. The maximum NPP occurs in Hainan and Yunnan, and the lowest occurs in Xinjiang and on the Qinghai-Tibet Plateau.
According to ecosystem types, the climax background NPP per unit area of the forest ecosystem in China was approximately 849.36 ± 88.32 gC/m 2 /a. For example, forest ecosystems were located in subtropical humid and temperate semi-humid monsoon climate zones due to their rich vegetation and the larger proportion forest vegetation, so their productivity was higher than that of other ecosystems.
The climax background NPP per unit area of the grassland ecosystem was approximately 306.75 ± 28.14 gC/m 2 /a. In the eastern Inner Mongolia Plateau and on the southeastern Qinghai-Tibet Plateau, grassland ecosystems were widely distributed, and the climax background vegetation NPP was mainly between 200 and 500 gC/m 2 /a.
The NPP per unit area of the desert ecosystem was approximately 67.56 ± 4.16 gC/m 2 /a. For example, the northwestern inland areas such as Gansu and Xinjiang had low productivity due to their dry climate and low annual precipitation, so the climax background vegetation NPP was lower than that of forest and grassland ecosystems.

Spatial Patterns of Vegetation NPP Restoration Potential
The spatial distribution of vegetation NPP restoration potential in China had obvious regional characteristics. Generally, it decreased from southeast to northwest along the water and heat gradients, which was similar to the spatial distribution of vegetation NPP. The vegetation NPP was relatively large in most parts of the south due to the better water and heat conditions, and the degree of human influence was also greater, so the restoration potential NPP was greater, and the restoration potential NPP values in the Yunnan-Guizhou Plateau and most of southern China were above 103.96 gC/m 2 /a. While the vegetation NPP was relatively small in arid and semiarid zones or parts of the Qinghai-Tibet Plateau, the restoration potential NPP was relatively small. The northwestern inland area was less than 50 gC/m 2 /a. In the agro-pastoral zone, the Loess Plateau, the eastern Sichuan Basin and the Greater Khingan Range areas, due to their relatively good recovery degree and local natural condition limitations, the recovery potential of most areas was less than 50 gC/m 2 /a. Therefore, the deployment of major ecological projects in the future should fully consider the limitations of climatic conditions. For areas with an annual precipitation below 300 mm and an annual average temperature below 0 • C (Figure 2), artificial ecological restoration measures should be avoided as much as possible to reduce human disturbance. The natural restoration of the ecology should be promoted.
The restoration potential is an absolute measure of the gap, but the size of the restoration potential does not indicate the difficulty of restoration. In some areas of the south where the restoration potential NPP is relatively large, it is only necessary to further strengthen the protection, restoration and construction measures of the ecosystem. With the help of excellent hydrothermal conditions, vegetation NPP can reach or approach the climax of natural vegetation within a certain period of time. In arid and semiarid regions or on the Qinghai-Tibet Plateau, which is restricted by harsh natural conditions, although the restoration potential is relatively small, ecological restoration will be a slow process and it will take a relatively long time to achieve the restoration of degraded vegetation.

Conclusions
Based on data from MODIS-NPP, meteorology and nature reserves, this paper used trend analysis, multiple regression analysis and residual analysis methods to analyse the change trend in NPP of China's terrestrial ecosystems, as well as the climax background, restoration status and restoration potential of forest, grassland and desert ecosystems from 2000 to 2019. The main conclusions are as follows: (1) The change in China's vegetation NPP showed a continuous upwards trend from 2000 to 2019. The two meteorological conditions, precipitation and temperature, contributed 85.41% to the change in vegetation NPP. It could meet the needs of the ecosystem's climax background vegetation NPP simulation in nature reserves that were less affected by human activities. (2) The results obtained by the regression analysis method of meteorological conditions based on nature reserves reflected the zonal climax vegetation status to a large extent, and the obtained values had a smooth transition within each eco-geographical division and between each eco-geographical division, which truly reflected the law of gradual changes in climate, vegetation and natural conditions. The annual total climax background NPP of China's forest, grassland and desert ecosystems was 2.76 ± 0.28 PgC, and the annual total NPP of the three ecosystems was 1.90 ± 0.2 PgC, 0.80 ± 0.07 PgC and 0.009 ± 0.0005 PgC, respectively. In future research, it is necessary to strengthen the acquisition of zonal climax background reference standards and try to select vegetation NPP in the core area of nature reserves and grassland enclosure areas that have been protected for many years. (3) Benefiting from the effective implementation of climate warming and humidification and ecological engineering, the agro-pastoral zone, the Loess Plateau, the eastern Sichuan Basin and the Greater Khingan Range had the most significant increase in the past 20 years. The total annual restoration status NPP in China's forest, grassland and desert ecosystems was 2.24 PgC, and the total annual NPP of the three ecosystems was 1.54 PgC, 0.65 PgC and 0.007 PgC, respectively. (4) The annual total restoration potential NPP of China's forest, grassland and desert ecosystems was approximately 0.52 ± 0.28 PgC, which accounted for approximately 19.05% of the total annual climax background NPP. The annual total restoration potential NPP of the forest, grassland and desert ecosystems was 0.36 ± 0.2 PgC, 0.16 ± 0.07 PgC and 0.002 ± 0.0005 PgC, respectively; additionally, the restoration potential accounted for 18.80%, 9.67% and 23.95% of the climax background NPP, respectively. The deployment of ecological projects should fully consider the restrictive climate conditions for decision makers and ecological assessment scholars, and the benefits and costs of the projects should be considered comprehensively.