Stand carbon storage and net primary production in China’s subtropical secondary forests are predicted to increase by 2060
Carbon Balance and Management volume 17, Article number: 6 (2022)
Forest ecosystems play an important role in carbon sequestration, climate change mitigation, and achieving China's target to become carbon (C) neutral by 2060. However, changes in C storage and net primary production (NPP) in natural secondary forests stemming from tree growth and future climate change have not yet been investigated in subtropical areas in China. Here, we used data from 290 inventory plots in four secondary forests [evergreen broad-leaved forest (EBF), deciduous and evergreen broad-leaved mixed forest (DEF), deciduous broad-leaved forest (DBF), and coniferous and broad-leaved mixed forest (CDF)] at different restoration stages and run a hybrid model (TRIPLEX 1.6) to predict changes in stand carbon storage and NPP under two future climate change scenarios (RCP4.5 and RCP8.5).
The runs of the hybrid model calibrated and validated by using the data from the inventory plots suggest significant increase in the carbon storage by 2060 under the current climate conditions, and even higher increase under the RCP4.5 and RCP8.5 climate change scenarios. In contrast to the carbon storage, the simulated EBF and DEF NPP declines slightly over the period from 2014 to 2060.
The obtained results lead to conclusion that proper management of China’s subtropical secondary forests could be considered as one of the steps towards achieving China’s target to become carbon neutral by 2060.
Forests cover 4.03 billion hectares worldwide, which is about 30% of the terrestrial surface of our planet, and account for 45% of terrestrial C and 75% of terrestrial gross primary production (GPP) [1, 2]. Forests play an important role in the global C cycle [3, 4]. Vegetation NPP measures the direct production capacity of forest ecosystems and C sequestration efficiency in the terrestrial C cycle [5, 6]. Attaining net zero emissions by mid-century will require all nations to comply with the objectives and principles laid out in the UN Framework Convention on Climate Change (UNFCCC) and its Paris Agreement. Accurate estimation of forest C storage and potential capacity is essential for achieving this goal. The current rate of change in global ecosystems is unprecedented. Forests can sequester C in ecosystems , and this is a safe and affordable strategy for mitigating the effects of climate change . Recently, the Chinese government officially announced its goal to peak CO2 emissions before 2030 and achieve C neutrality in 2060. There is thus an urgent need to determine how the NPP and C sequestration capability of forests might change in the future (especially in 2030 and 2060) under climate change.
Natural forest restoration is an effective approach for storing C [7, 8]. On average, natural forests have 6 times and 40 times C storage capacity of agroforests and plantations, respectively (sequestering 12, 1.9, and 0.3 Pg C per 100 M ha by 2100, respectively) . Therefore, natural forest restoration on disturbed lands should be prioritized . The most effective places for natural forest restoration for C storage are in the tropics and subtropics. Subtropical forests cover 11% of the land surface and store significant amounts of C . There are also large areas of naturally restored young forest on bare land in mountainous areas . Some EBF (the climax vegetation) have been converted into secondary forests via anthropogenic disturbance [7, 10]. Thus, subtropical secondary forests have become one of the most common types of natural forests that have a complex floristic composition and structure . The C sequestration ability of plantations and some natural forests has been well characterized, but our understanding of the role of secondary forests in sequestering C and mitigating climate change is limited .
Previous studies have shown that the C sequestration ability of forest ecosystems is affected by the interaction of various factors, such as forest origin, forest type, stand age, geography, and soil conditions [13, 14]. The C sequestration capacity varies greatly among forest types . C storage and NPP were significantly higher in broad-leaved forest than in CDF in a subtropical area . Stand age affects forest C storage and is an important variable for predicting future C sequestration [14, 16]. The growth of forests typically changes with forest age, and C stocks tend to be higher in old, complex forests [12, 14, 17, 18]. Global climate change will affect the structure and functions of forest ecosystems [19, 20]. According to the IPPC5, surface mean temperature is projected to increase in the subtropical region, and precipitation is projected to decrease (RCPs) . The warming tendency from 2011 to 2100 is 0.06 °C/10 a for RCP2.6, 0.24 °C/10 a for RCP4.5, and 0.63 °C/10 a for RCP8.5 . Precipitation will decrease by about 10% in subtropical China over the next century. Climate factors, such as temperature and precipitation, affect forest productivity, plant physiology, and community composition through the specific climatic constraints associated with different forests [23, 24]. Consequently, there is much uncertainty in the direction and magnitude of these effects on different forests [20, 25, 26]. More work is needed to investigate the ecosystem production and C storage capacity in subtropical secondary forests with different stand age and predict changes under future climate change conditions [26, 27]. The stabilization scenario that assumed the imposition of emissions mitigation policies  (RCP4.5) and the high—greenhouse gas concentration scenario that does not include any specific climate mitigation target  (RCP8.5) will serve as representatives , such work would aid the ability of forest managers to manage forests under future climate change.
The C storage and NPP in forests can be estimated using conversion parameter estimation, remote sensing, and model simulation. The models for simulating C storage and productivity are grouped into empirical statistical models based on inventory data of stand characteristics and site conditions [31, 32], mechanistic models based on the interaction between the physiological process of forest growth and environmental factors [33,34,35], and hybrid models based on statistical experience, environmental change, and a combination of forest physiological and ecological processes . The use of inventory, satellite, and field data for initialization, parameterization, and validation in process-based hybrid models is a promising approach for regional estimation and prediction in future and climate change [37,38,39]. TRIPLEX1.6 is a hybrid model  with simple inputs and high simulation accuracy . Compared with models that have been developed with a static climate or the classical assumption of stable site conditions, the TRIPLEX1.6 model can project the impact of climate change. Therefore, the TRIPLEX1.6 model is a suitable method for simulating and predicting C storage and NPP in subtropical secondary forest under future climate change.
Hunan Province located the subtropical region of southern China and has various secondary forests . We used national forest inventory plot data across Hunan Province and the TRIPLEX1.6 model to simulate stand C storage and NPP of subtropical secondary forests and predict the effects of two future climate change scenarios (RCP4.5 and RCP8.5) on stand C storage and NPP in different forest types. Specifically, the objectives of this study were to (1) characterize differences in C storage and NPP among forest types; (2) predict stand C storage and NPP in different forest types in 2030 and 2060; and (3) quantifying the impacts of stand age and climate change (RCP4.5 and RCP8.5) on C storage and NPP secondary forests in 2030 and 2060. The results of this study can provide a reliable data in C storage of natural secondary forest to help forest managers achieving the C neutrality target in 2060.
TIPLEX1.6 model validation
We used 290 forest stands according to the proportion of different forest types to calibrate the TRIPLEX1.6 model, and the remainder stands were conducted to validate the prediction of stem density, DBH, C storage and NPP. The simulated values were compared against the observations of all secondary forests to test the prediction accuracy of the model. There were high correlations between the simulated and observed values for stem density (R2 = 0.99, p < 0.01), stand average DBH (R2 = 0.96, p < 0.01), stand C storage (R2 = 0.87, p < 0.01), and NPP (R2 = 0.94, p < 0.01) in all forest stands (Fig. 1).
Dynamics of stand C storage and NPP with stand age among secondary forests
Simulations and predictions of TRIPLEX1.6 revealed the dynamics of stand C storage and NPP among secondary forest types over 100 years (Fig. 2). Overall, the results showed that there were significant differences in stand C storage and NPP between CDF and the three other subtropical secondary forests. Stand C storage gradually increased in all forests in 100 years, whereas the stand C storage of CDF was significantly lower than that of the other three forests. The dynamics of stand NPP revealed a difference in the growth rate of C storage among secondary forests. The NPP of three broad-leaved forests peaked in approximately 25 years and then decreased gradually over time, whereas the NPP of CDF increased for 50 years and then stabilized. The NPP was significantly lower in CDF than in the other three forests.
As the forest develops, the age composition will change over time. Most forests gradually change from young to mature and over-mature. We then predicted the C storage and NPP among secondary forests across the age groups from 2014 to 2060 (Figs. 3 and 4). C storage of over-mature forest was consistently the highest in all age groups among the different forest types, and the C storage across all age groups gradually increased. The NPP of young forest in the three broad-leaved forests was the highest in all age groups and then decreased over time (the young forest in DBF could not be well fitted); the NPP of mature forest in CDF was the highest and then increased over time.
Effects of climatic change on stand C storage and NPP in secondary forests
We predicted stand C storage and the difference in NPP between current and future climate scenarios (RCP4.5 and RCP8.5) among secondary forests (Table 1). We found that future climate changes were predicted to result in increased stand C storage and NPP for all subtropical forests. Overall, the response under RCP8.5 and RCP4.5 was the same. Under the same climate scenario, the response of NPP in forests in 2030 was stronger than that in 2060. The response of stand C storage and NPP of the four forest types under the climate scenarios differed. The difference in stand C storage for CDF was significantly lower compared with that for the other three forests (p < 0.001) (Table 1). The difference in NPP for EBF and CDF was significantly lower than that for DBF and DEF (p < 0.001) (Table 1).
Prediction of the stand C storage and NPP of secondary forests in 2030 and 2060
Overall, the average stand C storage of subtropical secondary forests is expected to reach 65.74 t C ha−1 in 2030 and 113.71 t C ha−1 in 2060. To understand the role of the four subtropical secondary forests in reaching targets of having CO2 emissions peak in 2030 and achieving C neutrality in 2060, stand C storage and NPP of the four forest types under the three climate scenarios (current, RCP4.5, and RCP8.5) were predicted (Table 2). Under all climate scenarios, stand C storage was the highest in EBF and smallest in CDF. The growth rate of C storage was the fastest in DEF from 2014 to 2030 and in CDF from 2030 to 2060. NPP was highest in DEF and lowest in CDF.
C sequestration capacity of subtropical secondary forests
The C sequestration capacity of subtropical secondary forests varies at different restoration stages. Predictions of stand C storage in 2060 by the TRIPLEX1.6 model indicated that stand C storage of EBF was greater than that of other forest types. Similar results were obtained in Jiangxi and Hunan [15, 32]. Previous studies that measured the biomass of permanent forest plots in subtropical forest [42, 43] have suggested that the above-ground biomass of subtropical old-growth EBF in China is 210–230 t ha−1 . This value is consistent with our simulated values of C storage using a conversion coefficient of 0.5. These findings indicate that EBF plays an important role in forest C storage and productivity capacity in subtropical areas. The key role of these late restoration stage forests in maintaining C stocks by preventing emissions derived from deforestation, forest degradation, and future climate warming emphasizes the need for habitat protection and sustainable management of complex and heterogeneous subtropical secondary forests .
There were also differences in the rate of C production among forest types. The increasing trend of NPP in CDF with stand age indicated that there were interdependent or showed complementary environmental adaptations, leading to higher resource utilization efficiency; improved use of their individual niches might also explain why the NPP in the forest increased [45, 46]. These findings indicate that CDF have the potential to support high C productivity. However, in most secondary forests, C stocks typically increase rapidly during the initial phase of regeneration, and then decelerate over subsequent decades and even centuries as primary forest species gradually colonize the area and grow to maturity . The maximum NPP of EBF occurred earlier compared with that of other forests because more productive forests show an earlier growth peak .
The study of C sequestration in secondary forests is particularly important given that the proportion of secondary subtropical forests is projected to continue to increase because of increases in anthropogenic and environmental disturbances [7, 47]. Previous research has shown that the greater amount of C in primary forest compared with secondary forest stems in large part to the difference in the number of large trees . Hence, the natural regeneration of forests is the cheapest and technically easiest option for achieving targets of having CO2 emissions peak in 2030 and reaching C neutrality in 2060. Protecting woodland from fire and other human disturbances allows trees to return and forests to flourish; it also allows C stocks to build rapidly and reach the level of C storage of a mature forest in roughly 70 years . The protection of subtropical natural secondary forest not only involves restoring it to the state of mature forest but also altering forest structure and composition so that it has a more stable C sequestration capacity similar to evergreen broad-leaved forest.
The influence of stand age on stand C storage and NPP
Forest age was an important factor in modulating the vegetation C storage. Owing to nutrient limitation, stomatal constraint, declines in photosynthesis during the stand development, stand NPP decreased with stand age [14, 49]. Previous studies indicate that NPP generally decreases in old forests to about half or one-third of its maximum value; other researchers have found old forests to be as productive as young forest stands . In our study, the predicted NPP of old forests (100 years) decreased by 7.77–16.50% from its peak for all forest types except for CDF. The continual recruitment of younger trees and rapid leaf area replacement, the compensatory growth of remaining vegetation, and canopy structural heterogeneity may offset this decrease in NPP in the moderately disturbed forests in this study . These results indicate that subtropical secondary forests in China still have high potential for C sequestration compared with other forests, and make C storage will continue to increase with stand growth. Particularly, C storage is related to changes in the age composition . These indicate that forest managers should pay more attention to the C sequestration contribution of old-growth forests.
Responses of stand C storage and NPP in subtropical forests to climate change
We compared the stand C storage and NPP between future climate change (RCP4.5 and 8.5) and the current state in 2030 and 2060. Climate change had a positive effect on forest C storage, and this might be explained by the fact that higher water availability (precipitation) and heat (temperature) can promote vegetation productivity [20, 50]. Nevertheless, the C sequestration potential under future climate change was still limited in the short and medium-term. This response of C storage to climate change is consistent with previous studies [51,52,53]. Stand C storage and NPP significantly differed across forest types (p < 0.001). Climate change responses of forests were more sensitive in the middle of recovery, and this result is consistent with Wu et al. . This might be driven by species-specific physiological characteristics and site conditions [41, 54]. Under disturbances or environmental changes, species asynchrony, facilitation, and species interactions might be the main drivers of ecosystem stability [11, 55]. Differences in the response of forest types and species to climate change contribute to forest community composition. Hence, understanding these differences would aid the development of measures for mitigating the negative effects of climate change, such as changes in stand structural attributes .
Model evaluation, limitations, and uncertainty
We calibrated and validated the TRIPLEX1.6 model with 875 permanent forest plot records in 2014 to obtain the optimal parameters for EBF, DEF, DBF, and CDF, which are four types of subtropical natural secondary forests in southern China. Our focus in adjusting model parameters was on the initial value of soil parameters. Soil fertility is an important initial input variable . The soil fertility should be set based on the growth of the forest, the thickness of the humus layer observed in the field, and the degree of stress that the soil has experienced in the past. The value of soil fertility in same area are provided in previous studies, however, subtropical ecosystems are vulnerable to extreme climate events caused by global warming . For example, the heavy snow/ice event that occurred in 2008 in southern China severely damaged forest ecosystems and resulted in decreased soil fertility. Therefore, the effect of historical natural disasters should be considered when setting the initial value of soil fertility. High consistency between the observed and simulated data from the model was achieved via model parameterization and repetitious adjustments.
Our prediction of C storage indicates that it will accumulate gradually with forest growth in medium-term, and changes in NPP among subtropical secondary forests were consistent with those reported in previous studies [15, 57]. All stand variables had R2 values of 0.87–0.99, which indicated that they could be used to predict growth C dynamics across various stand ages. These findings of model performance are consistent with previous studies showing that the TRIPLEX1.6 model is robust and can be used to predict the growth of subtropical natural secondary forests and artificial forests, including C. lanceolata and P. massoniana in southern and southeast China with R2 values of 0.94 for total biomass C  and 0.91 for total biomass . Combined with previous simulation results of coniferous forests, the TRIPLEX1.6 model could simulate forest growth with a high level of accuracy for all forest types in subtropical China.
The high biophysical heterogeneity and large amount of young age group subtropical secondary forests in subtropical China have increased the difficulty of evaluating the role of stand C storage and NPP in subtropical forests [57, 58]. Because of data limitations, the renewal of natural forest and under story regeneration in secondary forests was not considered in this study, which might reduce the robustness of our simulation results. Subtropical secondary forest has a high C sequestration capacity. In this study, the existing area of vegetation was not considered. Accurate prediction of spatial distributions is a major challenge and an important topic that should be examined in future studies. TRIPLEX1.6 simulations do not consider several variables related to contingency factors (i.e., anthropogenic disturbances, infestations, and extreme weather), which decreases their accuracy . Despite these uncertainties, the TRIPLEX1.6 model suggests that the secondary forests in subtropical China have a high capacity for C sequestration and storage, but the effects of climate change (RCP4.5 and RCP8.5) on stand C storage and NPP differ among forest types.
We used inventory plots data at different restoration stages and hybrid model of TRIPLEX 1.6 to predict changes in stand C storage and NPP under the two climate change scenarios (RCP4.5 and RCP8.5). Stand C storage of EBF was predicted to be the largest, and the NPP of DEF was predicted to be the highest. The NPP in CBF increased gradually over time and showed strong C storage potential capacity. Old-growth forests (mature and over-mature forests) have an important role in forest C sequestration. We not only considered the impact of stand age and forest types on C storage and NPP, but also made predictions about the impact of climate change under carbon neutrality in 2060. There was a significant difference in the effects of climate change (RCP4.5 and RCP8.5) on stand C storage between CDF and the other three forest types, but these changes were still limited in short and medium-term. These results indicate that the role of the floristic composition and tree growth of existing forests should be considered by forest managers for increasing C sequestration in secondary forests.
This study was carried out in Hunan Province (24°38′–30°08′ N and 108°47′–114°15′ E, Fig. 5) in central subtropical China. This province features mountains in the east, south, and west with elevations ranging from 21.60 to 2115.20 m (Fig. 5), hills in the east-central region, and plains in the northeast region surrounding Dongting Lake. The topography is heterogeneous, with slopes between 0° and 70°. The total land area is 21.18 million ha, including 12.53 million ha of forests, 3.92 million ha of arable land, 1.20 million ha of land covered by water bodies, and 1.00 million ha of developed land.
This region has a continental and subtropical humid monsoon climate that is rich in light, heat, and water resources (Fig. 6). The mean annual rainfall is 1500 mm, with 75% falling between March and August. The mean annual temperature is 14.10 °C. Annual climate changes are large: summers are hot, and winters are short and cold. The soil types include red soil, yellow soil, and red-yellow soil, with small areas of calcareous soil, purple soil, moisture soil, and mountain meadow soil (Fig. 6), which are mainly developed from granite, limestone, and shale.
Forest types and stand age groups
The National Forest Inventory of China classified subtropical forests into four forest types according to the dominant tree species: pine forest, fir forest, deciduous forest, and evergreen broad-leaved forest . These forest types encompass most but not all wide spread tree species in Hunan. DEF are an important forest type in the restoration process that have not been examined extensively in previous studies because of the difficulty of determining the dominant tree species and stand structure. Thus, we divided subtropical secondary forest into four forest types: CDF, DBF, DEF, and EBF. Forest type identification was based on the technical regulations for continuous forest inventory: coniferous or broad-leaved species accounted for less than 70% of the species in CDF; the proportion of broad-leaved trees accounted for more than 70% of the species in broad-leaved forest; and the proportion of all broad-leaved species accounted for less than 70% of the species in DEF . Among the natural secondary forests in Hunan, CDF are mainly dominated by Pinus massoniana forests; the DBF are mainly dominated by Alniphyllum fortunei, Choerospondias axillaris, and Liquidambar formosana. Hance; and tree species in the EBF are mainly dominated by Cyclobalanopsis glauca Thunb, Litsea coreana Levl. var. sinensis, and Schima superba Gardn.
The study area is composed of subtropical natural uneven-aged forest. Therefore, five of the largest trees outside each plot were selected for stem core sampling, and the ring counts of these five tree samples were averaged to estimate stand age . According to the forestry standards of China, the age composition of four forest types can be divided into young forest, middle-aged forest, premature forest, mature forest, and over-mature forest. The age group of the two mixed forests was determined by the age group of the dominant tree species and their growth condition [An additional table file shows this in more detail (see Additional file 1)].
TRIPLEX1.6 model description
TRIPLEX is a hybrid model [An additional table file shows this in more detail (see Additional file 2)] that integrates three well-established process-based sub-models: the forest production sub-model 3-PG , the forest growth and yield sub-model TREEDYN3.0 , and the soil-C, soil-N (nitrogen) and soil–water-balance sub-model CENTURY4.0 . The forest production sub-model estimates monthly GPP (including above- and below-ground biomass) from photosynthetic active radiation (PAR), mean air temperature, vapor pressure deficit (VPD), soil water, the percentage of frost days, and the leaf area index. The forest growth and yield sub-model calculates tree growth and yield variables (height, DBH) using a function of the stem wood biomass increment . The soil C, soil N, and soil–water balance sub-model simulate soil C and N dynamics between litter and soil pools, simulate water balances and dynamics, and calculate monthly water loss through transpiration, evaporation, soil water content, and snow water content . TRIPLEX combines empirical and mechanistic components that can simulate various growth volumes in forest ecosystems using input data derived from site, soil, climate, and stand growth . The hybrid model compensates for the inability of simplified growth models to take into account soil factors, feedback between the forest ground and underground, and climate impacts. The simulation of the TRIPLEX model involves key processes and dynamics including PAR, GPP, forest growth, biomass, soil C, soil nitrogen, and soil water . After the model was optimized to version 1.6, a reasonable and balanced parameter generalization procedure that did not lead to a significant reduction of model accuracy but increased model practicability was described .
The TRIPLEX1.6 simulation requires input data such as latitude, longitude, soil texture, monthly climate records, tree species physiological variables (such as maximum tree height and diameter), tree species process mediators, stand structure, and certain initial site conditions. Simulation outputs include tree diameter, height, basal area, total volume, leaf area index, GPP, net ecosystem production (NEP), biomass, soil C, and N and water dynamics. TRIPLEX1.6 has been tested and applied to forest growth and biomass production on a regional scale at Zhejiang Province, China [39, 63] and northeastern China  and has also been used to predict Cunninghamia lanceolata and P. massoniana stand production in Hunan Province, China . Natural secondary forests in a subtropical region were examined in this paper, and the suitability of the model in the subtropical region has been demonstrated by previous studies .
Model input data
There were five primary data sources: permanent forest-plot (667 m2) records taken from the forest inventory, climate datasets, field observations, the literature, and assumptions. In Hunan Province, there were 875 permanent sample plots of natural secondary forest in 2014. The permanent forest plot database is a well-designed forest inventory program built from data acquired at 5-year intervals beginning in 1973. Field investigations in the permanent forest plot database include records of location, site conditions, and DBH of each tree; average stand height ; and land use, site class, dominant tree species, stand density, age, average DBH, average tree height, and volume . Plots included at least one of the four forest types and met the criteria for tree size (e.g., DBH ≥ 5.0 cm) and stage of stand development (with the exception of seedlings). Plots were required to meet various criteria regarding site, climate, and growth, and the initial data needed to be detailed and representative of the forest types, as well as of the study area as a whole . According to the forest plot database, stand biomass can be estimated using allometric equations  [An additional table file shows this in more detail (see Additional file 3)] for the main subtropical tree species. The conversion parameter between biomass and C storage was 0.5. Based on data collected in 2009 and 2014, we calculated average annual NPP (t ha−1 year−1) using the following formula:
where Bt1 and Bt2 are the plot biomass in year t1 (2009) and t2 (2014), respectively. Climate data for each permanent plot were interpolated from data collected from meteorological stations . The data recorded between 2005 and 2014 were used as climate inputs for each plot. The climate inputs were monthly frost days, monthly average air temperature, monthly sum of precipitation, and the monthly average atmospheric VPD . Atmospheric N deposition was set to 18 kg N ha−1 year−1 . Field observations were made at independent sampling plots in Hunan Province. The choice of reference literature and assumptions were based on indicator functions for the four forest types selected for this study [66,67,68].
Predictors for the TRIPLEX1.6 model were climate normal from the 96 automatic meteorological monitoring stations in Hunan Province for current (2005–2014) and future (2015–2060) periods. The future climate change scenario data were derived from the WCRP’s CMIP5 multi-model dataset provided by the National Climate Center of China Meteorological Administration (http://www.climatechange-data.cn). According to the fifth report of the IPCC, we obtained data for current scenarios, RCP4.5 (medium emissions scenario, assumes the imposition of emissions mitigation policies), and RCP8.5 (high emissions scenario, does not include any specific climate mitigation target) to simulate the responses of different forest types and age composition to future climate change. Because the CMIP5 GCMs have different horizontal resolutions, future climate data were interpolated on a common 1° × 1° grid using the nearest neighbor interpolation method.
Model initialization, parameterization, and simulations
From the initial modelling, we examined the most sensitive parameters through direct measurement, literature searches, default values, and statistical fitting. To simulate forest ecosystem processes and dynamics, TRIPLEX1.6 requires some initial values of stand variables describing conditions of forest stands and soils. There are three key variables (stand density, tree height, and DBH) related to initial conditions for forest growth and yield simulation . To ensure the robustness of the TRIPLEX1.6 model, most of the general and nonspecific site parameters were derived from previous studies. These include PAR parameters; the minimum, maximum, and optimum temperature for tree growth; stomata and canopy conductance; initial N for tree growth; the lignin-nitrogen ratio and lignin fraction of leaf and fine and coarse roots; and the fraction of soil water flow . Several new parameters such as wood C density, specific leaf area (SLA), mortality, the fraction of leaf, branch, wood, and coarse and fine roots were used and adjusted from default model values to better represent the forest ecosystems of subtropical China for this study [An additional table file shows this in more detail (see Additional file 4 and 5)].
TRIPLEX1.6 was calibrated by randomly selecting one-third of each forest type and validated by remainder data before simulation runs. The simulation was executed for DBH, stem density, C storage, and NPP. We simulated each stand from its respective year of regeneration to 2014, at which point all simulations across all stands within the natural subtropical forest region were summed. The same procedure was followed for all model runs by Zhou et al. .
Availability of data and materials
All data parameters about the TRIPLEX1.6 model on all of the conclusions of the manuscript rely to be presented in the additional supporting files (see Additional files 6 and 7). And the statistical parameters of sample plot data and climatic conditions and calculation method of forest vegetation carbon are also mentioned in the additional supporting files (see Additional files 8, 9, and 10).
Net primary production
Net ecosystem production
Gross primary productivity
Photosynthetic active radiation
Vapor pressure deficit
Diameter at breast height (1.3 cm)
Permanent sample plots
Evergreen broad-leaved forest
Deciduous and evergreen broad-leaved mixed forest
Deciduous broad-leaved forest
Conifer and broad-leaved mixed forest
UN framework convention on climate change
Li Y, Liu XJ, Xu WB, Bongers FJ, Bao WK, Chen B, Chen GK, Guo K, Lai JS, Lin DM, Mi XC, Tian XJ, Wang XH, Yan JH, Yang B, Zheng YR, Ma KP. Effects of diversity, climate and litter on soil organic carbon storage in subtropical forests. For Ecol Manag. 2020;476: 118479. https://doi.org/10.1016/j.foreco.2020.118479.
Pan YD, Birdsey RA, Phillips OL, Jackson RB. The structure, distribution, and biomass of the world’s forests. Annu Rev Ecol Evol Syst. 2013;44:593–622. https://doi.org/10.1146/annurev-ecolsys-110512-135914.
Bordin KM, Esquivel-Muelbert A, Bergamin RS, Klipel J, Picolotto RC, Frangipani MA, Müller SC. Climate and large-sized trees, but not diversity, drive above-ground biomass in subtropical forests. For Ecol Manag. 2021;490: 119126. https://doi.org/10.1016/j.foreco.2021.119126.
Siddiq Z, Hayyat MU, Khan AU, Mahmood R, Shahzad L, Ghaffar R, Cao KF. Models to estimate the above and below ground carbon stocks from a subtropical scrub forest of Pakistan. Glob Ecol Conserv. 2021;27: e01539. https://doi.org/10.1016/j.gecco.2021.e01539.
Bao G, Tuya A, Bayarsaikhan S, Dorjsuren A, Mandakh U, Bao Y, Vanchindorj B. Variations and climate constraints of terrestrial net primary productivity over Mongolia. Quatern Int. 2020;537:112–25. https://doi.org/10.1016/j.quaint.2019.06.017.
Zhang C, Wu S, Leng G. Possible NPP changes and risky ecosystem region identification in China during the 21st century based on BCC-CSM2. J Geog Sci. 2020;30(8):1219–32. https://doi.org/10.1007/s11442-020-1778-8.
Lewis SL, Wheeler CE. Regenerate natural forests to store carbon. Nature. 2019;568:25–8.
Wood A, Tolera M, Snell M, O’Hara P, Hailu A. Community forest management (CFM) in south-west Ethiopia: maintaining forests, biodiversity and carbon stocks to support wild coffee conservation. Glob Environ Chang. 2019;59: 101980. https://doi.org/10.1016/j.gloenvcha.2019.101980.
FAO, 2020. Global forest resources assessment. Rome: FAO; 2020.
Zhang J, Ge Y, Chang J, Jiang B, Jiang H, Peng CH, Zhu JR, Yuan WG, Qi LZ, Yu SQ. Carbon storage by ecological service forests in Zhejiang Province, subtropical China. For Ecol Manag. 2007;245(1):64–75. https://doi.org/10.1016/j.foreco.2007.03.042.
Ouyang S, Xiang WH, Gou MM, Chen L, Lei PF, Xiao WF, Deng XW, Zeng LX, Li JR, Zhang T, Peng CH, Forrester DI, Meyer C. Stability in subtropical forests: the role of tree species diversity, stand structure, environmental and socio-economic conditions. Glob Ecol Biogeogr. 2021;30:500–13. https://doi.org/10.1111/geb.13235.
Seedre M, Janda P, Trotsiuk V, Hedwall PO, Morrissey RC, Mikoláš M, Bače R, Čada V, Svoboda M. Biomass carbon accumulation patterns throughout stand development in primary uneven-aged forest driven by mixed-severity natural disturbances. For Ecol Manag. 2020;455: 117676. https://doi.org/10.1016/j.foreco.2019.117676.
Ouyang S, Xiang WH, Wang XP, Xiao WF, Chen L, Li SG, Sun H, Deng XW, Forrester DI, Zeng LX, Lei PF, Lei XD, Gou MM, Peng CH. Effects of stand age, richness and density on productivity in subtropical forests in China. J Ecol. 2019;107:2266–77. https://doi.org/10.1111/1365-2745.13194.
Xu L, Shi YJ, Fang HY, Zhou GM, Xu XJ, Zhou YF, Tao JX, Ji B, Xu J, Li C, Chen L. Vegetation carbon stocks driven by canopy density and forest age in subtropical forest ecosystems. Sci Total Environ. 2018;631–632:619–26. https://doi.org/10.1016/j.scitotenv.2018.03.080.
Wang SQ, Zhou L, Chen JM, Ju WM, Feng XF, Wu WX. Relationships between net primary productivity and stand age for several forest types and their influence on China’s carbon balance. J Environ Manag. 2011;92:1651–62. https://doi.org/10.1016/j.jenvman.2011.01.024.
Zheng ZT, Zhu WQ, Zhang YJ. Seasonally and spatially varied controls of climatic factors on net primary productivity in alpine grasslands on the Tibetan Plateau. Glob Ecol Conserv. 2020;21: e00814. https://doi.org/10.1016/j.gecco.2019.e00814.
Cheng J, Lee X, Theng BKG, Zhang L, Fang B, Li F. Biomass accumulation and carbon sequestration in an age-sequence of Zanthoxylum bungeanum plantations under the Grain for Green Program in karst regions, Guizhou province. Agric For Meteorol. 2015;203:88–95. https://doi.org/10.1016/j.agrformet.2015.01.004.
Ryan MG, Binkley D, Fownes JH. Age-related decline in forest productivity: pattern and process. In: Begon M, Fitter AH, editors. Advances in ecological research. Cambridge: Academic Press; 1997. p. 213–62. https://doi.org/10.1016/S0065-2504(08)60009-4.
Chen B, Zhang X, Tao J, Wu J, Wang J, Shi P, Yu C. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau. Agric For Meteorol. 2014;189–190:11–8. https://doi.org/10.1016/j.agrformet.2014.01.002.
Xie Y, Wang H, Lei X. Simulation of climate change and thinning effects on productivity of Larix olgensis plantations in northeast China using 3-PG mix model. J Environ Manage. 2020;261:110249–110249. https://doi.org/10.1016/j.jenvman.2020.110249.
IPCC. Climate change 2013: the physical science basis near-term climate change: projections and predictability. Geneva: IPCC; 2013.
Xu CH, Xu Y. The projection of temperature and precipitation over China under RCP scenarios using a CMIP5 multi-model ensemble. Atmos Ocean Sci Lett. 2012;5:527–33. https://doi.org/10.1080/16742834.2012.11447042.
Laflower DM, Hurteau MD, Koch GW, North MP, Hungate BA. Climate-driven changes in forest succession and the influence of management on forest carbon dynamics in the Puget Lowlands of Washington State, USA. For Ecol Manag. 2016;362:194–204. https://doi.org/10.1016/j.foreco.2015.12.015.
Wu C, Chen D, Shen J, Sun X, Zhang S. Estimating the distribution and productivity characters of Larix kaempferi in response to climate change. J Environ Manag. 2021;280:111633. https://doi.org/10.1016/j.jenvman.2020.111633.
Sherry RA, Weng ES, Arnone JA III, Johnson DW, Schimel DS, Verburg PS, Wallace LL, Luo YQ. Lagged effects of experimental warming and doubled precipitation on annual and seasonal aboveground biomass production in a tallgrass prairie. Glob Chang Biol. 2008;14:2923–36. https://doi.org/10.1111/j.1365-2486.2008.01703.x.
Teng MJ, Zeng LX, Hu WJ, Wang PC, Yan ZG, He W, Zhang Y, Huang ZL, Xiao WF. The impacts of climate changes and human activities on net primary productivity vary across an ecotone zone in Northwest China. Sci Total Environ. 2020;714:136691–136691. https://doi.org/10.1016/j.scitotenv.2020.136691.
Xu C, Liu M, An S, Chen JM, Yan P. Assessing the impact of urbanization on regional net primary productivity in Jiangyin County, China. J Environ Manag. 2007;85(3):597–606. https://doi.org/10.1016/j.jenvman.2006.08.015.
Thomson AM, Calvin KV, Smith SJ, Kyle GP, Volke A, Patel P, Delgado-Arias S, Bond-Lamberty B, Wise MA, Clarke LE, Edmonds JA. RCP4.5: a pathway for stabilization of radiative forcing by 2100. Clim Chang. 2011;109(1–2):77–94. https://doi.org/10.1007/s10584-011-0151-4.
Riahi K, Rao S, Krey V, Cho C, Chirkov V, Fischer G, Kindermann G, Nakicenovic N, Rafaj P. RCP 8.5—a scenario of comparatively high greenhouse gas emissions. Clim Chang. 2011;109(1–2):33–57. https://doi.org/10.1007/s10584-011-0149-y.
Lei XD, Yu L, Hong LX. Climate-sensitive integrated stand growth model (CS-ISGM) of Changbai larch (Larix olgensis) plantations. For Ecol Manag. 2016;376:265–75. https://doi.org/10.1016/j.foreco.2016.06.024.
Gupta R, Sharma LK. The process-based forest growth model 3-PG for use in forest management: a review. Ecol Model. 2019;397:55–73. https://doi.org/10.1016/j.ecolmodel.2019.01.007.
Xie Y, Wang H, Lei X. Application of the 3-PG model to predict growth of Larix olgensis plantations in northeastern China. For Ecol Manag. 2017;406:208–18. https://doi.org/10.1016/j.foreco.2017.10.018.
Bossel H. TREEDYN3 forest simulation model. Ecol Model. 1996;90:187–227. https://doi.org/10.1016/0304-3800(95)00139-5.
Gilmanov TG, Parton WJ, Ojima DS. Testing the ‘CENTURY’ ecosystem level model on data sets from eight grassland sites in the former USSR representing a wide climatic/soil gradient. Ecol Model. 1997;96:191–210. https://doi.org/10.1016/S0304-3800(96)00067-1.
Wang F, Mladenoff DJ, Forrester JA, Keough C, Parton WJ. Global sensitivity analysis of a modified CENTURY model for simulating impacts of harvesting fine woody biomass for bioenergy. Ecol Model. 2013;259:16–23. https://doi.org/10.1016/j.ecolmodel.2013.03.008.
Peng CH, Liu JX, Dang QL, Apps MJ, Jiang H. TRIPLEX: a generic hybrid model for predicting forest growth and carbon and nitrogen dynamics. Ecol Model. 2002;153:109–30. https://doi.org/10.1016/S0304-3800(01)00505-1.
Peng CH, Zhou XL, Zhao SQ, Wang XP, Zhu B, Piao SL, Fang JY. Quantifying the response of forest carbon balance to future climate change in Northeastern China: model validation and prediction. Global Planet Chang. 2009;66:179–94. https://doi.org/10.1016/j.gloplacha.2008.12.001.
Zhao MF, Xiang WH, Deng XW, Tian DL, Huang ZH, Zhou XL, Yu GR, He HL, Peng CH. Application of TRIPLEX model for predicting Cunninghamia lanceolata and Pinus massoniana forest stand production in Hunan Province, southern China. Ecol Model. 2013;250:58–71. https://doi.org/10.1016/j.ecolmodel.2012.10.011.
Zhang J, Chu ZY, Ge Y, Zhou XL, Jiang H, Chang J, Peng CH, Zheng JW, Jiang B, Zhu JR, Yu SQ. TRIPLEX model testing and application for predicting forest growth and biomass production in the subtropical forest zone of China’s Zhejiang Province. Ecol Model. 2008;219(3):264–75. https://doi.org/10.1016/j.ecolmodel.2008.07.016.
Zhou XL, Peng CH, Dang QL. Assessing the generality and accuracy of the TRIPLEX model using in situ data of boreal forests in central Canada. Environ Model Softw. 2004;19(1):35–46. https://doi.org/10.1016/S1364-8152(03)00108-7.
Dai E, Wu Z, Ge Q, Xi W, Wang X. Predicting the responses of forest distribution and aboveground biomass to climate change under RCP scenarios in southern China. Glob Chang Biol. 2016;22:3642–61. https://doi.org/10.1111/gcb.13307.
Lin D, Lai J, Muller-Landau HC, Mi X, Ma K. Topographic variation in aboveground biomass in a subtropical evergreen broad-leaved forest in China. PLoS ONE. 2012;7: e48244. https://doi.org/10.1371/journal.pone.0048244.
McEwan RW, Lin YC, Sun IF, Hsieh CF, Su SH, Chang LW, Michael Song GZ, Wang HH, Hwong JL, Lin KC, Yang KC, Chiang JM. Topographic and biotic regulation of aboveground carbon storage in subtropical broad-leaved forests of Taiwan. For Ecol Manag. 2011;262:1817–25. https://doi.org/10.1016/j.foreco.2011.07.028.
Xu YZ, Franklin SB, Wang QG, Shi Z, Luo YQ, Lu ZJ, Zhang JX, Qiao XJ, Jiang MX. Topographic and biotic factors determine forest biomass spatial distribution in a subtropical mountain moist forest. For Ecol Manag. 2015;357:95–103. https://doi.org/10.1016/j.foreco.2015.08.010.
Tilman D, Lehman CL, Thomson KT. Plant diversity and ecosystem productivity: theoretical considerations. Proc Natl Acad Sci. 1997;94:1857–61. https://doi.org/10.1073/pnas.94.5.1857.
Xiang WH, Fang X. Community structure and soil properties in subtropical secondary forests. Beijing: Science Press; 2018.
Ngo KM, Turner BL, Muller-Landau HC, Davies SJ, Larjavaara M, Nik H, Nik FB, Lum S. Carbon stocks in primary and secondary tropical forests in Singapore. For Ecol Manag. 2013;296:81–9. https://doi.org/10.1016/j.foreco.2013.02.004.
Van Tuyl S, Law BE, Turner DP, Gitelman AI. Variability in net primary production and carbon storage in biomass across Oregon forests—an assessment integrating data from forest inventories, intensive sites, and remote sensing. For Ecol Manag. 2005;209:273–91. https://doi.org/10.1016/j.foreco.2005.02.002.
Li Y, Bao WK, Bongers F, Chen B, Chen GK, Guo K, Jiang MX, Lai JS, Lin DM, Liu CJ, Liu XJ, Liu Y, Mi XC, Tian XJ, Wang XH, Xu WB, Yan JH, Yang B, Zheng YR, Ma KP. Drivers of tree carbon storage in subtropical forests. Sci Total Environ. 2019;654:684–93. https://doi.org/10.1016/j.scitotenv.2018.11.024.
Wang G, Guan D, Xiao L, Peart MR. Forest biomass-carbon variation affected by the climatic and topographic factors in Pearl River Delta, South China. J Environ Manag. 2019;232:781–8. https://doi.org/10.1016/j.jenvman.2018.11.130.
Ren Y, Chen SS, Wei XH, Xi WM, Luo YJ, Song XD, Zuo SD, Yang YS. Disentangling the factors that contribute to variation in forest biomass increments in the mid-subtropical forests of China. J For Res. 2016;27:919–30. https://doi.org/10.1007/s11676-016-0237-y.
Wang WJ, He HS, Thompson FR, Fraser JS, Dijak WD. Changes in forest biomass and tree species distribution under climate change in the northeastern United States. Landsc Ecol. 2017;32:1399–413. https://doi.org/10.1007/s10980-016-0429-z.
Wu Z, Dai E, Wu Z, Lin M. Assessing differences in the response of forest aboveground biomass and composition under climate change in subtropical forest transition zone. Sci Total Environ. 2020;706: 135746. https://doi.org/10.1016/j.scitotenv.2019.135746.
Bayer AD, Pugh TAM, Krause A, Arneth A. Historical and future quantification of terrestrial carbon sequestration from a Greenhouse-Gas-Value perspective. Glob Environ Chang. 2015;32:153–64. https://doi.org/10.1016/j.gloenvcha.2015.03.004.
Guo B, Zang WQ, Luo W. Spatial-temporal shifts of ecological vulnerability of Karst Mountain ecosystem-impacts of global change and anthropogenic interference. Sci Total Environ. 2020;741:140256–140256. https://doi.org/10.1016/j.scitotenv.2020.140256.
Luo DW, Huang JG, Jiang XY, Ma QQ, Liang HX, Guo XL, Zhang SK. Effect of climate and competition on radial growth of Pinus massoniana and Schima superba in China’s subtropical monsoon mixed forest. Dendrochronologia. 2017;46:24–34. https://doi.org/10.1016/j.dendro.2017.08.001.
Piao SL, Fang JY, Zhou LM, Zhu B, Tan K, Tao S. Changes in vegetation net primary productivity from 1982 to 1999 in China. Glob Biogeochem Cycles. 2005;19:GB2027.1-GB2027.16. https://doi.org/10.1029/2004GB002274.
Tang XL, Wang YP, Zhou GY, Zhang DQ, Liu S, Liu SZ, Zhang QM, Liu JX, Yan JH. Different patterns of ecosystem carbon accumulation between a young and an old-growth subtropical forest in Southern China. Plant Ecol. 2011;212:1385–95. https://doi.org/10.1007/s11258-011-9914-2.
Chen XY, Mulder J. Atmospheric deposition of nitrogen at five subtropical forested sites in South China. Sci Total Environ. 2007;378:317–30. https://doi.org/10.1016/j.scitotenv.2007.02.028.
NFGA. Technical regulations for continuous forest inventory. GB/T 38590–2020. Merrimack: NFGA; 2020.
Landsberg JJ, Waring RH. A generalised model of forest productivity using simplified concepts of radiation–use efficiency, carbon balance and partitioning. For Ecol Manag. 1997;95:209–28. https://doi.org/10.1016/S0378-1127(97)00026-1.
Parton WJ, Scurlock JMO, Ojima DS, Gilmanov TG, Scholes RJ, Schimel DS, Kirchner T, Menaut J-C, Seastedt T, Garcia Moya E, Kamnalrut A, Kinyamario JI. Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide. Glob Biogeochem Cycles. 1993;7:785–809. https://doi.org/10.1029/93GB02042.
Zhang H, Song TQ, Wang KL, Wang GX, Liao JX, Xu GH, Zeng FP. Biogeographical patterns of forest biomass allocation vary by climate, soil and forest characteristics in China. Environ Res Lett. 2015;10(4):44014. https://doi.org/10.1088/1748-9326/10/4/044014.
Lei XD, Tang MP, Lu YC, Hong LX, Tian DL. Forest inventory in China: status and challenges. Int For Rev. 2009;11:52–63.
Fang JY, Chen AP, Peng CH, Zhao SQ, Ci LJ. Changes in forest biomass carbon storage in China between 1949 and 1998. Science. 2001;292:2320–2.
Lin DM, Lai JS, Mi XC, Ren HB, Ma KP. Spatial variation in community structure of a subtropical evergreen broad-leaved forest: Implications for sampling design. Chin Sci Bull. 2013;58:1181–6. https://doi.org/10.1007/s11434-012-5586-0.
Xiang WH, Liu SH, Deng XW, Shen AH, Lei XD, Tian DL, Zhao MF, Peng CH. General allometric equations and biomass allocation of Pinus massoniana trees on a regional scale in southern China. Ecol Res. 2011;26:697–711. https://doi.org/10.1007/s11284-011-0829-0.
Zhao MF, Xiang WH, Peng CH, Tian DL. Simulating age-related changes in carbon storage and allocation in a Chinese fir plantation growing in southern China using the 3-PG model. For Ecol Manag. 2009;257(6):1520–31. https://doi.org/10.1016/j.foreco.2008.12.025.
We thank the professional English editing team from the Charlesworth Group (see https://www.cwauthors.com.cn) to improve the language of this manuscript.
This study was supported by the Key Research and Development Program of Hunan Province (2021NK2031); the National Key Research and Development Program of China (2016YFD0600202); Scientific Research Innovation Project of Hunan Province (CX20210895) and the Huitong Forest Ecological Station funded by the National Forestry and Grassland Administration of China (2021132078).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Age group classification standard of subtropical natural secondary forests in Hunan Province.
The structural model of forest growth and carbon simulation from TRIPLEX1.0 (modified from Peng et al., 2002). Rectangles represent key pools or state variables, ovals represent core simulation processes, dotted lines represent controls, and solid lines represent the flow of carbon (C), nitrogen (N), water, and the fluxes between the forest ecosystem and external environment. Two arrow cycles refer to two feed-backs.
Compatibility superposition relative growth equation (lnWi = ßi0 + ßi1 × lnD ) regression coefficients of biomass ( Wi, kg) and DBH (D, cm) for each part of the tree species (values in parentheses is standard error) and fitting parameters.
Site-specific parameters used in TRIPLEX1.6 for simulating subtropical forest ecosystems in Hunan Province, subtropical China (Zhao et al., 2013).
Stand growth and tree species-specific parameters TRIPLEX1.6 applied to simulate growth of evergreen broad-leaved forest, deciduous and evergreen broad-leaved mixed forest, deciduous broad-leaved forest, and coniferous and broad-leaved mixed forest.
Descriptive statistics for the spatial heterogeneity of stand growth by forests in all selected permanent forest plots in Hunan Province in 2014. Including stand-structure variables: stem density (stems ha-1), stand age (year), diameter at breast height (DBH, cm), height (m), NPP (t ha-1 yr-1) and stand biomass (t ha-1).
Descriptive statistics for spatial heterogeneity of climate and forest stand structure in all selected permanent forest sample plots in Hunan Province. Site variable: elevation (above mean sea level, m); climatic variables: annual average air temperature (°C), annual rainfall (mm yr-1) and the annual average vapor pressure deficit (mbar) between 2000 and 2014.
Simulation errors of TRIPLEX1.6 applied to subtropical forest ecosystems in southeastern China, comparing density (stems ha−1), DBH (cm), NPP (t ha-1 yr-1) and C storage (t C ha−1) between modeled values and forest inventory data collected from 875 forest stands.
Variation of monthly average temperature and precipitation in Hunan Province from 2000 to 2014.
Variation of monthly average temperature and precipitation in Hunan Province from 2000 to 2014.
About this article
Cite this article
Jin, J., Xiang, W., Zeng, Y. et al. Stand carbon storage and net primary production in China’s subtropical secondary forests are predicted to increase by 2060. Carbon Balance Manage 17, 6 (2022). https://doi.org/10.1186/s13021-022-00204-y
- Carbon (C) storage
- Climate change
- Forest restoration
- Net primary production (NPP)
- TRIPLEX model