Skip to main content

Variation of net primary productivity and its drivers in China’s forests during 2000–2018

Abstract

Background

Net primary productivity (NPP) in forests plays an important role in the global carbon cycle. However, it is not well known about the increase rate of China’s forest NPP, and there are different opinions about the key factors controlling the variability of forest NPP.

Methods

This paper established a statistics-based multiple regression model to estimate forest NPP, using the observed NPP, meteorological and remote sensing data in five major forest ecosystems. The fluctuation values of NPP and environment variables were extracted to identify the key variables influencing the variation of forest NPP by correlation analysis.

Results

The long-term trends and annual fluctuations of forest NPP between 2000 and 2018 were examined. The results showed a significant increase in forest NPP for all five forest ecosystems, with an average rise of 5.2 gC·m− 2·year− 1 over China. Over 90% of the forest area had an increasing NPP range of 0–161 gC·m− 2·year− 1. Forest NPP had an interannual fluctuation of 50–269 gC·m− 2·year− 1 for the five major forest ecosystems. The evergreen broadleaf forest had the largest fluctuation. The variability in forest NPP was caused mainly by variations in precipitation, then by temperature fluctuations.

Conclusions

All five forest ecosystems in China exhibited a significant increasing NPP along with annual fluctuations evidently during 2000–2018. The variations in China’s forest NPP were controlled mainly by changes in precipitation.

Background

Vegetation net primary productivity (NPP), as a biomass increment, plays an important role in the global carbon cycle. It contributes to understanding the contribution of NPP to the carbon exchanges that take place between the biosphere and atmosphere. Vegetation NPP indicates the atmospheric carbon fixed by plants except for the carbon released by respiration. Global vegetation NPP has increased due to climate change over the last few decades (Nemani et al. 2003; Chen et al. 2019). A significant increase in NPP was observed for grasslands and shrublands after the 1980s across China (Piao et al. 2001; Liu et al. 2017; Liang et al. 2015; Xu and Wang 2016; Wang et al. 2017). However, some research has found declines in global NPP over the past decade because of drought (Zhao and Running 2010; Medlyn 2011).

There are various process models for the estimation of vegetation NPP. The Carnegie-Ames-Stanford Approach (CASA) provides a useful tool to estimate NPP by combining grid meteorological data and the satellite-derived fraction of incoming photosynthetically active radiation absorbed by vegetation (Smith et al. 2008). The Carbon Exchange between Vegetation, Soil and Atmosphere (CEVSA) model was used to explore NPP in China’s terrestrial ecosystems, and confirmed a slow increase in NPP over the period 1981–1998 (Tao et al. 2003). The Boreal Ecosystem Productivity Simulator (BEPS) was used to map the spatio-temporal variations of NPP in China for the period from 2000 to 2010 (Chen et al. 2007; Liu et al. 2013). The Lund-Potsdam-Jena Dynamic Global Vegetation Model (LPJ-DGVM) was widely applied to predict NPP dynamics by combining terrestrial vegetation dynamics (Smith et al. 2008; Gao et al. 2016). Based on the eddy covariance technique, an eco-physiological process model was developed to calculate NPP directly from gross primary productivity (GPP) by the mean ratio of the NPP/GPP of different vegetation types in a local area (Wang et al. 2017). All the process models can be applied in other regions. However, traditional process models have many eco-physiological parameters, and the parameters need to be modified for different research areas. It is difficult to ensure the accuracy of the parameters. Additionally, there are uncertainties in the estimation of NPP because of accumulated errors during the complex process.

Statistics models for vegetation NPP are relatively simple but they can provide a credible output for a specific research area, although they may be difficult to apply to other regions (Adams et al. 2004). The Miami model is a classical statistics-based model that has been used widely for estimating vegetation NPP in different regions (Lieth 1973; Alexandrov and Matsunaga 2008). In China, Zhou and Zhang (1996) established a general statistical model. The model is suitable for predicting NPP for terrestrial vegetation over China, but its accuracy is not very good for estimating forest NPP. Until now, there has not been a specific statistics-based model to predict forest NPP accurately.

Both process models and statistics models need to identify the dominant factors affecting NPP, since different dominant factors might lead to different simulation results. For example, Zhao and Running (2010) reported that global NPP had declined over the past decade, but their findings were based on outcomes from models which assumed a strong temperature dependence (Medlyn 2011). Despite recognition of the roles of temperature and precipitation in controlling NPP dynamics, there are differing opinions over their relative influence (Nemani et al. 2003; Zhang et al. 2017). Some observations have revealed that NPP was strongly correlated with annual precipitation at a continental scale, but interannual variability in ANPP (above-ground NPP) was not related to variability in precipitation (Knapp and Smith 2001). This conclusion is controversial, since it was reported that the interannual variation in NPP could be explained by changes in precipitation in arid and semi-arid regions, and by changes in air temperature in other regions of China (Chen et al. 2013; Liang et al. 2015; Xu and Wang 2016). In addition, vegetation structure (e.g., canopy cover and leaf area index) can change due to natural growth, climate change and human disturbance (Chen et al. 2019). NPP dynamics are likely to be affected by changes in vegetation structure, but most studies ignore the effect of vegetation structure dynamics on NPP fluctuations. The effect of environment variables in NPP dynamics will most likely vary in different regions.

This paper established a statistics-based multiple regression model to estimate forest NPP across China. The objectives were to (1) explore the long-term trends and annual fluctuations of forest NPP for different forest types at a national scale during the period 2000–2018; and (2) identify the key environmental factors controlling the variability of forest NPP. This study helps to understand the variability of forest NPP, as well as the interactions between forest productivity and environmental factors.

Materials and methods

Building the statistics model

The paper assumed that the spatiotemporal patterns of forest NPP were controlled mainly by climatic, topographic and forest structure variables, although previous research has shown different relative dominance between factors affecting NPP. To establish the statistics multivariate regression model, we collected long-term observations of annual forest NPP, normalized difference vegetation index (NDVI, which was used to quantify vegetation structure), altitude, average precipitation and temperature between 1980 and 2010 at 1000 forest inventory sites (Fig. 1).

Table 1 The parameters of the statistics-based multiple regression model between forest net primary productivity (NPP) and environmental variables (T is annual temperature (°C), P is annual precipitation (mm), A is altitude (m), NDVI is normalized difference vegetation index (0 < NDVI ≤1))
Fig. 1
figure 1

Distribution of forest types and forest observation sites where environmental data were collected. Black dots show forest observation sites used to establish the model, and the red dots show forest observation sites used to check the model

Local topography significantly affects the spatial variation of climatic variables, and each elevation increase of 100 m in China leads to an average reduction of 25 gC·m− 2 for annual NPP (Chen et al. 2007). The topographic effect on NPP can be described by altitude, so altitude from a digital elevation model (DEM) was selected as an important variable in establishing the statistics-based multivariate regression model. The NDVI was extracted from the annual maximum NDVI of a 16-day cycle, 2000–2010 (https://lpdaac.usgs.gov).

Forest NPP was calculated using the traditional Volume Growth Rate Model from the data of forest inventories, 1980–2010 (Eq. 1) (Luo 1996).

$$ \mathrm{NPP}=\mathrm{Ps}+\mathrm{Pb}+\mathrm{Pl}+\Pr +\mathrm{Pu} $$
(1)

where Ps, Pb, Pl and Pr are the annual net biomass of stem, branch, leaf and root in the tree layer, respectively. Pu is the annual net biomass of shrubs and herbs under the forest canopy.

Quality control of forest NPP had been done. Firstly, forest NPP for each observation site was the average value of continuous forest inventory since 1980, and the average value reduced data uncertainty. Secondly, forest NPP for each observation site was checked by comparing with the values of the surrounding area and the values of the same forest types. Thirdly, forest NPP was checked by comparing with the results of some forest productivity models (Luo 1996; Luo et al. 2004).

Based on the relationship between forest NPP and environmental variables (annual maximum NDVI, annual precipitation, annual temperature, and altitude) at 1000 forest inventory sites, a statistics-based multiple regression model was established (Table 1), with a multiple regression R = 0.72, and a standard error of 4.34 (Eq. 2),

$$ \mathrm{NPP}=97.13\mathrm{NDVI}+0.022 PT+0.128P-9.136T-0.027A+333.67 $$
(2)

where NPP is forest NPP (gC·m− 2·year− 1), NDVI is normalized difference vegetation index (0 < NDVI ≤1), P is annual precipitation (mm), T is annual temperature (°C), and A is altitude (m).

Simulation of annual forest NPP over China

The statistics-based multiple regression model allowed us to simulate annual forest NPP over China. The data included climatic variables (annual temperature and annual precipitation), altitude and NDVI. Annual temperature and annual precipitation were obtained from 800 meteorological stations of the National Meteorological Information Center of the China Meteorological Administration (http://cdc.cma.gov.cn). Altitude was derived from a 90-m resolution DEM dataset (http://srtm.csi.cgiar.org/srtmdata/). NDVI with a 1-km spatial resolution was obtained from Land Processes Distributed Active Archive Center (https://lpdaac.usgs.gov). To reduce the residual noise caused by haze and clouds, the annual maximum NDVI was extracted from a series of 16-day compositions of NDVI in the growing season. All the input variables were interpolated or resampled into raster layers with a spatial resolution of 1 km × 1 km grid cells using the Kriging method with ArcGIS10.2 (http://www.esri.com/software/arcgis).

The distribution of annual forest NPP was obtained from overlay analysis between the Vegetation Map of the People’s Republic of China in 2008 and the simulated annual forest NPP. Prior to overlay analysis, all the forest types were merged into five dominant forest types, including evergreen broadleaf forest, evergreen needleleaf forest, broadleaf-needleleaf mixed forest, deciduous broadleaf forest, deciduous needleleaf forest.

Accuracy assessment

The observed forest NPP at 90 observation sites from forest inventories were used to check the accuracy of the multiple regression model. The check sites were selected following the principle of representativeness. The 90 check sites covered the typical forest types in China, and their spatial distribution could not be concentrated in one area.

Comparison between the predicted NPP and the observed NPP provided a fitting coefficient of R2 = 0.536 (P < 0.01) (Fig. 2). The results indicated that the multiple regression model had a good performance in predicting forest NPP with acceptable accuracy. Therefore, the model was capable of providing a reliable estimate of China’s forest NPP.

Fig. 2
figure 2

Comparison between the observed net primary productivity (NPP) with the predicted NPP at 90 forest observation sites over China

Results

Trend and fluctuation of forest NPP

The results revealed an average forest NPP of 840 gC·m− 2·year− 1 for all forest types in China. Evergreen forests had the highest NPP (1058 gC·m− 2·year− 1 for the evergreen broadleaf forest, 934 gC·m− 2·year− 1 for the evergreen needleleaf forest), because evergreen forest maintains continuous photosynthesis throughout the year. Deciduous forests had lower NPP values (759 gC·m− 2·year− 1 for the deciduous broadleaf forest, 590 gC·m− 2·year− 1 for the deciduous needleleaf forest) since deciduous forests have a narrower photosynthesis period.

Variations (long-term trend and interannual fluctuation) of forest NPP in China were examined. Results revealed an overall increasing trend in forest NPP, with an average increase rate of 5.2 gC·m− 2·year− 1 over China between 2000 and 2018. All the five major forest ecosystems showed significant increasing forest NPP, with an increase rate of 2.9–8.0 gC·m− 2·year− 1. The largest increase rate was found in evergreen broadleaf forest (8.0 gC·m− 2·year− 1), followed by evergreen needleleaf forest (6.5 gC·m− 2·year− 1), broadleaf-needleleaf mixed forest (6.3 gC·m− 2·year− 1), deciduous broadleaf forest (5.4 gC·m− 2·year− 1) and deciduous needleleaf forest (2.9 gC·m− 2·year− 1) (Fig. 3).

Fig. 3
figure 3

Long-term change and annual fluctuation of forest net primary productivity (NPP) in five forest ecosystems over China, 2000–2018. EBF is evergreen broadleaf forest; DBF is deciduous broadleaf forest; ENF is evergreen needleleaf forest; DNF is deciduous needleleaf forest; BNMF is broadleaf-needleleaf mixed forest

Interannual fluctuation in forest NPP displayed a fluctuation range of 50–269 gC·m− 2·year− 1 among the five major forest ecosystems. NPP in the evergreen broadleaf forest had the highest fluctuation range (269 gC·m− 2·year− 1), followed by the evergreen needleleaf forest (202 gC·m− 2·year− 1). NPP in the deciduous needleleaf forest had the lowest fluctuation range (50 gC·m− 2·year− 1).

The fluctuation rate (ratio of fluctuation range to average value) was calculated to reflect the degree of fluctuation. The results showed that the evergreen broadleaf forest had the largest fluctuation rate (25%), followed by the evergreen needleleaf forest (22%) and the broadleaf-needleleaf mixed forest (20%). The deciduous needleleaf forest had the smallest fluctuation rate (8%) (Fig. 3 and Table 2). The results indicated that, among the major forest ecosystems in China, the most unstable forest NPP occurred in the evergreen broadleaf forest, and the most stable forest NPP occurred in the deciduous needleleaf forest.

Table 2 The average annual net primary productivity (NPP) and its interannual fluctuation for five main forest types over China during the period 2000–2018

Comparing forest NPP with its fluctuation, our results revealed an interesting phenomenon: the greater the forest NPP, the greater the fluctuation of forest NPP for the five major forest ecosystems in China. For example, evergreen broadleaf forest had the largest NPP at 1058 gC·m− 2·year− 1, with the largest fluctuation range, 269 gC·m− 2·year− 1, as well as the largest fluctuation rate (25%) among the five major forest ecosystems. Deciduous needleleaf forest had the smallest NPP at 590 gC·m− 2·year− 1, with the smallest fluctuation amplitude, 50 gC·m− 2·year− 1, as well as the smallest fluctuation rate (8%) (Fig. 3 and Table 2). It is likely that this phenomenon is related to the fluctuation of dominant factors affecting NPP, but it needs further confirmation.

Spatial distribution of forest NPP trends

Figure 4 shows the rate of change of forest NPP during the period 2000–2018 at 1 km × 1 km spatial resolution. China’s forest NPP exhibited an increasing rate of 0–161 gC·m− 2·year− 1, which was statistically significant over 90% of forest area during the period 2000–2018. Forest NPP over China did not exhibit a consistent increasing trend everywhere. Approximately 10% of forest area in some regions of southwest and northern China showed a slight decreasing trend in forest NPP, with a range 0–20 gC·m− 2·year− 1 (Fig. 4). Generally, there was an overall increasing trend in forest NPP, despite a slightly decreasing trend in some forest areas in China.

Fig. 4
figure 4

Spatial distribution of trend rate (gC·m− 2·year− 1) of forest net primary productivity (NPP), 2000–2018. Positive values denote an increasing rate, negative values denote a decreasing rate

Key drivers of the variation in forest NPP

To identify the key drivers leading to the variation in forest NPP, the fluctuation values (deviation from the mean) of NPP and the environment variables were extracted. Correlation analysis showed a significant correlation between NPP and precipitation fluctuations, with a Pearson coefficient 0.977 (P < 0.01), Kendall tau-b coefficient 0.908 (P < 0.01), Spearman’s rho coefficient 0.973 (P < 0.01). A moderate correlation was found between NPP and temperature fluctuations, with Pearson coefficient 0.484 (P < 0.01) (Table 3). The results indicated that the variation in forest NPP was much more sensitive to annual precipitation fluctuation and annual precipitation fluctuation played a pivotal role in controlling forest NPP changes compared with other environment variables.

Table 3 The correlations between fluctuations (deviation from the mean) of net primary productivity (NPP) and the driving factors (annual precipitation, annual temperature and normalised difference vegetation index (NDVI))

Comparing the interannual change curves of forest NPP, annual precipitation, annual temperature, and annual maximum NDVI, more similar curves were discovered between NPP and annual precipitation. The results provide another credible proof that forest NPP fluctuation was sensitive to changes in precipitation. In contrast, temperature and NDVI made a smaller contribution to forest NPP (Fig. 5).

Fig. 5
figure 5

Interannual variation curves of annual net primary productivity (NPP) and the major environmental factors (annual precipitation, annual temperature and annual maximum normalised difference vegetation index (NDVI)) in Chinese forests, 2000–2018

We also examined the spatial patterns of the changes of forest NPP, precipitation and temperature. It was discovered that the increase in forest NPP was accompanied by an increase in precipitation and temperature for most forest areas in China. The decrease in forest NPP was accompanied by a decrease in precipitation for a few forest areas in southwest China (Figs. 4, 6 and 7).

Fig. 6
figure 6

Change trend rate (mm·year− 1) of annual precipitation in China, 2000–2018 (Positive values denote an increase, negative values denote a decrease)

Fig. 7
figure 7

Change trend rate (°C·year− 1) of annual temperature in China, 2000–2018 (Positive values denote an increase, negative values denote a decrease)

In brief, precipitation fluctuation was a key factor leading to variations in forest NPP, and temperature was a secondary factor, judging by the correlation coefficients, the change curves and the spatial patterns of change trends.

Discussions

There are different estimates of China’s forest NPP, ranging from 420.1–843 gC·m− 2·year− 1 according to previous studies (Table 4). It was difficult to judge whether one was more credible than the others because of the use of different methods. However, all (previous and our) studies conclude that forest NPP has shown an obvious increasing trend over China in recent decades (Wang et al. 2008; Mao et al. 2010; Yu et al. 2014). The increase in forest NPP is consistent with other findings that the terrestrial vegetation NPP of the northern hemisphere has increased over the past several decades (Nemani et al. 2003). Our results presented an average increase rate of 5.2 gC·m− 2·year− 1 for China’s forest NPP in the period 2000–2018, despite different increase rates for the five major forest ecosystems in China (Fig. 3).

Table 4 Forest net primary productivity (NPP) and its trends in China in recent decades according to previous studies

Previous studies have ignored the impact of forest structure on NPP, even though forest spatial structure might affect forest NPP (Fotis et al. 2018). The general statistical model of Zhou and Zhang (1996) as a typical climate productivity model also ignored the impact of forest structure, and forest NPP from the model of Zhou and Zhang was usually lower than the observed forest NPP. Our model selected vegetation structure (quantified by NDVI) as an important factor, since spatio-temporal change in forest structure was directly related to forest NPP dynamics. Our results confirmed that forest structure had significant effect on forest NPP dynamics (with Spearman’s rho coefficient 0.490, P < 0.05). Therefore, vegetation structure as an input variable in our model improved the accuracy of simulation.

There are a series of factors that may affect NPP, but they do not all have a strong influence on NPP dynamics. The dynamics of the dominant factors affecting NPP are the real causes of changes to NPP. Therefore, we extracted the fluctuation values (deviation from the mean) of NPP and environment variables prior to analysing their relationship. The relationship between these fluctuation values provided more reliable information on the factors controlling NPP dynamics.

Our results revealed that changes in precipitation played a key role in the variation in forest NPP over China, and temperature fluctuation was a secondary factor affecting forest NPP variation. The results disagreed with previous results that air temperature was the dominant climatic factor that controlled the interannual variability in NPP throughout China, except for arid and semi-arid regions (Liang et al. 2015). The disagreement might result from different analysis methods. The results of Liang et al. (2015) came from a simple correlation analysis between annual NPP and climate conditions (i.e., mean annual air temperature and annual cumulative precipitation) from 1982 to 2010. Our results came from the correlation analysis between the fluctuation values of annual forest NPP and environmental factors (annual precipitation, annual temperature, and annual maximal NDVI) from 2000 to 2018. Furthermore, our results from the correlation analysis were confirmed by the fluctuation curves and the spatial trends of forest NPP and environmental factors.

However, there were some uncertainties in the simulation results of our model. Our model paid more attention to natural factors, but ignored human factors. In fact, vegetation NPP in China benefited from the ecological conservation and restoration efforts, since a series of ecological protection policies had been implemented, such as Three-North Shelter Forest Program, Grain for Green Program, and Natural Forest Protection Program (Lv et al. 2015; Zhu and Zheng 2019). Although our model selected vegetation structure (quantified by NDVI) as an important factor, and vegetation structure could reflect the impact of human beings on vegetation NPP to some extent, it was not enough to reflect the overall human impact.

Conclusions

This paper established a statistics-based multiple regression model to estimate annual forest NPP across China. The model had a good performance in providing a reliable estimation in forest NPP. The trends and the fluctuations of forest NPP in the period 2000–2018 were examined. The results revealed a significant increasing forest NPP for all the five major forest ecosystems, with an average increase rate of 5.2 gC·m− 2·year− 1 over China. There was an increasing NPP trend rate of 0–161 gC·m− 2·year− 1 over 90% of the forest area. The interannual fluctuation range of forest NPP was 50–269 gC·m− 2·year− 1. The evergreen broadleaf forest had the largest fluctuation (269 gC·m− 2·year− 1), followed by the evergreen needleleaf forest, the broadleaf-needleleaf mixed forest, the deciduous broadleaf forest and the deciduous needleleaf forest. Our results revealed that the variability in forest NPP was caused mainly by variations in precipitation, then by temperature fluctuations. Therefore, precipitation played a key role in the variation in forest NPP.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

NPP:

Net primary productivity

ANPP:

Above-ground net primary productivity

GPP:

Gross primary productivity

CASA:

Carnegie-Ames-Stanford Approach

CEVSA:

Vegetation, Soil and Atmosphere model

BEPS:

Boreal Ecosystem Productivity Simulator

LPJ-DGVM:

Lund-Potsdam-Jena Dynamic Global Vegetation Model

DEM:

Digital elevation model

NDVI:

Normalized difference vegetation index

References

  • Adams B, White A, Lenton TM (2004) An analysis of some diverse approaches to modelling terrestrial net primary productivity. Ecol Model 177:353–391

    Article  CAS  Google Scholar 

  • Alexandrov GA, Matsunaga T (2008) Normative productivity of the global vegetation. Carbon Balance Manag 3:8

    Article  Google Scholar 

  • Chen C, Park T, Wang XH, Piao SL, Xu BD, Chaturvedi RK, Fuchs R, Brovkin V, Ciais P, Fensholt R, Tømmervik H, Bala G, Zhu CC, Nemani RR, Myneni RB (2019) China and India lead in greening of the world through land-use management. Nat Sustain 2:122–129

    Article  Google Scholar 

  • Chen XF, Chen JM, An SQ, Ju WM (2007) Effects of topography on simulated net primary productivity at landscape scale. J Environ Manag 85:585–596

    Article  CAS  Google Scholar 

  • Chen Z, Yu GR, Ge JP, Sun XM, Hirano T, Saigusa N, Wang QF, Zhu XJ, Zhang YP, Zhang JH, Yan JH, Wang HM, Zhao L, Wang YF, Shi PL, Zhao FH (2013) Temperature and precipitation control of the spatial variation of terrestrial ecosystem carbon exchange in the Asian region. Agric For Meteorol 182:266–276

    Article  Google Scholar 

  • Fotis AT, Morin TH, Fahey RT, Hardiman BS, Bohrer G, Curtis PS (2018) Forest structure in space and time: biotic and abiotic determinants of canopy complexity and their effects on net primary productivity. Agric For Meteorol 250:181–191

    Article  Google Scholar 

  • Gao QZ, Guo YQ, Xu HM, Ganjurjav H, Li Y, Wan YF, Qin XB, Ma X, Liu S (2016) Climate change and its impacts on vegetation distribution and net primary productivity of the alpine ecosystem in the Qinghai-Tibetan plateau. Sci Total Environ 554-555:34–41

    Article  CAS  Google Scholar 

  • Knapp AK, Smith MD (2001) Variation among biomes in temporal dynamics of aboveground primary production. Science 291:481–484

    Article  CAS  Google Scholar 

  • Liang W, Yang YT, Fan DM, Guan HD, Zhang T, Long D, Zhou Y, Bai D (2015) Analysis of spatial and temporal patterns of net primary production and their climate controls in China from 1982 to 2010. Agric For Meteorol 204:22–36

    Article  Google Scholar 

  • Lieth H (1973) Primary production: terrestrial ecosystems. Human Ecol 1:303–332

    Article  Google Scholar 

  • Liu G, Sun R, Xiao ZQ, Cui TX (2017) Analysis of spatial and temporal variation of net primary productivity and climate controls in China from 2001 to 2014. Acta Ecol Sin 37(15):4936–4945

    Google Scholar 

  • Liu YB, Ju WM, He HL, Wang SQ, Sun R, Zhang YD (2013) Changes of net primary productivity in China during recent 11 years detected using an ecological model driven by MODIS data. Front Earth Sci 7(1):112–127

    Article  CAS  Google Scholar 

  • Luo TX (1996) Patterns of net primary productivity for Chinese major forest types and their mathematical models. Ph.D. dissertation, Chinese Academy of Sciences, Beijing

    Google Scholar 

  • Luo TX, Pan YD, Ouyang H, Shi PL, Luo J, Yu ZL, Lu Q (2004) Leaf area index and net primary productivity along subtropical to alpine gradients in the Tibetan plateau. Glob Ecol Biogeogr 13:345–358

    Article  Google Scholar 

  • Lv YH, Zhang LW, Feng XM, Zeng Y, Fu BJ, Yao XL, Li JR, Wu BF (2015) Recent ecological transitions in China: greening, browning, and influential factors. Sci Report 5:8732

    Article  Google Scholar 

  • Mao JF, Dan L, Wang B, Dai YJ (2010) Simulation and evaluation of terrestrial ecosystem NPP with M-SDGVM over continental China. Adv Atmos Sci 27(2):427–442

    Article  Google Scholar 

  • Medlyn BE (2011) Comment on “drought-induced reduction in global terrestrial net primary production from 2000 through 2009”. Science 333:1093

    Article  CAS  Google Scholar 

  • Nemani RR, Keeling CD, Hashimoto H, Jolly WM, Piper SC, Tucker CJ, Myneni RB, Running SW (2003) Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 300:1560–1563

    Article  CAS  Google Scholar 

  • Ni J (2003) Net primary productivity in forests of China: scaling-up of national inventory data and comparison with model predictions. Forest Ecol Manag 176:485–495

    Article  Google Scholar 

  • Piao SL, Fang JY, Guo QH (2001) Terrestrial net primary production and its spatio-temporal patterns in China during 1982–1999. Acta Sci Nat Univ Pek 37(4):563–569

    Google Scholar 

  • Smith B, Knorr W, Widlowski JL, Pinty B, Gobron N (2008) Combining remote sensing data with process modelling to monitor boreal conifer forest carbon balances. For Ecol Manag 255(12):3985–3994

    Article  Google Scholar 

  • Tao B, Li K, Shao XM, Cao MK (2003) The temporal and spatial patterns of terrestrial net primary productivity in China. J Geogr Sci 13(2):163–171

    Article  Google Scholar 

  • Wang Q, Zhao P, Ren H, Kakubari Y (2008) Spatiotemporal dynamics of forest net primary production in China over the past two decades. Glob Planet Chang 61:267–274

    Article  Google Scholar 

  • Wang X, Tan K, Chen BZ, Du PJ (2017) Assessing the spatiotemporal variation and impact factors of net primary productivity in China. Sci Report 7:44415

    Article  Google Scholar 

  • Xu HJ, Wang XP (2016) Effects of altered precipitation regimes on plant productivity in the arid region of northern China. Ecol Inform 31:137–146

    Article  Google Scholar 

  • Yu C, Wang B, Liu H, Yang XS, Xiu ZZ (2014) Dynamic change of net production and mean net primary productivity of China’s forests. Forest Res 27(4):542–550

    Google Scholar 

  • Zhang ZY, Zhang RD, Cescatti A, Wohlfahrt G, Buchmann N, Zhu J, Chen GH, Moyano F, Pumpanen J, Hirano T, Takagi K, Merbold L (2017) Effect of climate warming on the annual terrestrial net ecosystem CO2 exchange globally in the boreal and temperate regions. Sci Report 7:3108

    Article  CAS  Google Scholar 

  • Zhao MS, Running SW (2010) Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 329:940–943

    Article  CAS  Google Scholar 

  • Zhou GS, Zhang XS (1996) Study on NPP of natural vegetation of China under global climate change. Acta Phyt Sin 20(1):11–19

    Google Scholar 

  • Zhu JJ, Zheng X (2019) The prospects of development of the three-north afforestation program (TNAP): on the basis of the results of the 40-year construction general assessment of the TNAP. Chin J Ecol 38:1600–1610

    Google Scholar 

  • Zhuang QL, Zhang TL, Xiao JF, Luo TX (2009) Quantification of net primary production of Chinese forest ecosystems with spatial statistical approaches. Mitig Adapt Strateg Glob Change 14:85–99

    Article  Google Scholar 

Download references

Acknowledgments

We gratefully acknowledge Leonie Seabrook, PhD, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Funding

This study was supported by the National Natural Science Fundation of China (No. 41571175, 31661143028) and the special funds for basic research and operation from the Chinese Academy of Meteorological Science (2017Y003).

Author information

Authors and Affiliations

Authors

Contributions

Guangsheng Zhou & Yuhe Ji conceived and designed the study. Yuhe Ji wrote the paper. Tianxiang Luo, Li Zhou and Xiaomin Lv revised and perfected the design of the experiments. Yakir Dan reviewed and edited the manuscript. All authors read and approved the manuscript.

Corresponding author

Correspondence to Guangsheng Zhou.

Ethics declarations

Ethics approval and consent to participate

The subject has no ethic risk.

Consent for publication

All the data and relevant processing methods, including the automated generation of figures of evaluation results, belong to the host institution, namely, Chinese Academy of Meteorological Science.

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ji, Y., Zhou, G., Luo, T. et al. Variation of net primary productivity and its drivers in China’s forests during 2000–2018. For. Ecosyst. 7, 15 (2020). https://doi.org/10.1186/s40663-020-00229-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40663-020-00229-0

Keywords