Forest management required for consistent carbon sink in China’s forest plantations

Forest is the largest biomass carbon (C) pool in China, taking up a substantial amount of atmospheric carbon dioxide. Although it is well understood that planted forests (PFs) act as a large C sink, the contribution of human management to C storage enhancement remains obscure. Moreover, existing projections of forest C dynamics suffer from spatially inconsistent age and type information or neglected human management impacts. In this study, using developed PF age and type maps and data collected from 1371 forest plantation sites in China, we simulated biomass C stock change and quantified management impacts for the time period 2010–2050. Results show that future forest biomass C increment might have been overestimated by 32.5%–107.5% in former studies. We also found that age-related growth will be by far the largest contributor to PF biomass C increment from 2010 to 2050 (1.23 ± 0.002 Pg C, 1 Pg = 1015 g = 1 billion metric tons), followed by the impact of human management (0.57 ± 0.02 Pg C), while the contribution of climate is slight (0.087 ± 0.04 Pg C). Besides, an additional 0.24 ± 0.07 Pg C can be stored if current PFs are all managed by 2050, resulting in a total increase of 2.13 ± 0.05 Pg C. Forest management and age-related growth dominate the biomass C change in PFs, while the effect of climatic factors on the accumulation is minor. To achieve the ambitious goal of forest C stock enhancement by 3.5 Pg from 2020 to 2050, we advocate to improve the management of existing forests and reduce the requests for more lands for forest expansion, which helps mitigate potential conflicts with agricultural sectors. Our results highlight that appropriate planning and management are required for sustaining and enhancing biomass C sequestration in China’s PF.


Background
Global forest absorbs carbon (C) equivalent to~34% of the emission from fossil-fuel combustion and cement production, and biomass C augmentation plays a dominant role (Pan et al. 2011). For Chinese forests, approximately 80%-94% of the C sink lies in biomass, leaving only a minor fraction of sink in soil organic matter Piao et al. 2009). Recent field-based surveys also confirmed that post-afforestation soil C sink was much smaller than biomass C sink, and even soil C loss was found after afforestation (Hong et al. 2020). Compared with relatively slow and variable C stored in soil after tree-planting, biomass C accumulation in forest is more rapid and attainable than soil C enhancement. Since the Paris climate summit in 2015, it is increasingly advocated to accumulate C in terrestrial ecosystem as an efficient mitigation option (Rumpel et al. 2018), indicating that vegetation biomass is expected to contribute highly in compensating C emission. Besides, China has pledged to become C neutral before 2060 at the UN General Assembly meeting in 2020, and forest C sink enhancement plays a fundamental role to achieve this ambitious goal. Therefore, many of the former studies have estimated the forest biomass C accumulation potential in China (Xu et al. 2010;He et al. 2017;Yao et al. 2018). Nonetheless, none of the former studies examined the human management impacts in determining future biomass C sink at a large scale. Existing projections of the biomass C changes in China's forests were attributed to either age-related growth or climatic factors (Xu et al. 2010;He et al. 2017;Yao et al. 2018). However, neglecting the role of human management impact can adversely affect sustainable C accumulation as the governments would be blindly guided to pursue PF's area expansion at the expense of forest tending and nursing. Therefore, clarifying the mechanisms and drivers of forest biomass C accumulation is a prerequisite not only for sustainable C sink, but also for diverting C from the atmosphere (CO 2 ) into terrestrial ecosystems.
A novel national-wide field campaign was conducted in China, in the framework of which thousands of samples were collected during 2011-2015. Field survey was performed in all 32 provinces, autonomous regions, and municipalities of mainland China, following a consistent and rigorous protocol . Based on the substantial number of samples collected, Tang et al. (2018) analyzed spatial distributions and drivers of C pools in China's forests; however, the human management impact was neglected in the C sink projection subjected to a simple biomass-age relationship derived from pooling of all samples. This study aims to reveal forest biomass C accumulation potential as a function of human management impacts. Specifically, the directlymeasured data from the China-wide field campaign were revisited, filtered, and grouped for estimating the size, change, and potential of PF biomass C in China since 2010. The results of this study provide a new perspective of the role of human interventions to maximizing the mitigation of atmospheric CO 2 , and can enlighten planning and management practices for sustaining and enhancing biomass C sequestration in China's PFs. These results can also enlighten the path toward achieving the UN SDGs for sustainable use of terrestrial ecosystems and forests management.

Field plot data
A nationwide field campaign was conducted in China, during which site information was collected as detailed as possible, including geographical characteristics, soil properties, vegetation properties, disturbances, and human management (Yu et al. 2020a, b). Disturbances refer to impacts related to lack of management (see Additional file 1 Fig. S2 in Supporting Information), which may negatively affect forest growth, while "force majeure" disturbances were excluded (e.g. wind throw, lightning fire). This helps identifying signals of forest growth under lack of management, facilitating accurate projection of forest biomass C enhancement due to management. A total of 1371 forest plantation sites surveyed were selected for analysis, while the remaining sites were excluded due to missing information of forest age, vegetation biomass C, disturbance, or management (Additional file 1, Figs. S1 and S2). Disturbances recorded included grazing, tourism, resin tapping, and treading, amongst others (Additional file 1 Fig. S2). Human management describes management types implemented (e.g. thinning, fertilization, irrigation, pesticide application, weeding, fencing, tending, and pruning), but details of the timing and frequency of management practices were not available. Thus, we divided the sites into two categories according to the presence or absence of human management. Sites that were not managed (including unmanaged sites influenced by disturbances) were allocated to the group "unmanaged", indicating the status of PF growth in the absence of management or exposure to a historical disturbance regime. For consistency with the use of forest maps reconstructed in a previous study (Yu et al. 2020a), we further grouped the sites into 17 major forests for each category (Additional file 1 Fig. S1, Table S1).

Planted forest distribution and forest type maps
The reliability of gridded map-based biomass C projections is largely determined by the reliability of forest age and type information, which is challenging for developing reliable gridded forest maps. Essentially, but often ignored, forest age map should be used in combination with the forest type information that the age map originated from forest biomass C estimation. This is understandable as forest age and type maps derived independently provide spatial mismatch of information, which may cause large biases if used together for biomass C simulations (Yu et al. 2020a). To overcome the shortcomings in existing forest maps (e.g. inconsistent forest type and age information), we have developed a model for reconstructing spatially consistent, 1-km resolution age and type maps of PFs in China (Yu et al. 2020a), facilitating accurate PF-specific C sink estimations. More specifically, the maps were generated previously using multiple sources of data, including different land cover products (e.g. MODIS Land Cover Type, European Space Agency Climate Change Initiative Land Cover data, Global Forest Change 2000-2014 from Hansen et al. (2013)), the 6th, 7th, and 8th (1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013) national forest inventories, and the digitalized PF map produced by the State Forestry Administration of China. More details about the development of the maps can be found in Yu et al. (2020a).
The PF maps used in this study depict the distribution of forest stands only, while economic forests (for producing fruits, oil, spices, herbs, etc.) and bamboo forests were excluded. Since the PF age and type maps delineate the distribution of PF in 2005, we further updated the maps to 2010 in this study. More specifically, the new PFs from 2005 to 2010 were assigned ages (randomly from 1 to 5 years) and types in each province according to the 8th national forest inventory report. We also proportionally adjusted PF age maps to match the inventory reports for each province, which helps to reduce biases we identified previously (e.g. Fujian Province, Yu et al. 2020a). To this end, the age map grid cells are multiplied by the ratio of map-averaged forest age to inventory forest age for each province. Yao et al. (2018) proposed a series of models for forest biomass C growth simulation by taking into account the age and climate factors (i.e. mean annual air temperature and precipitation). In this study, we adopted a similar approach, but introduced elevation as an additional variable to the model to improve the biomass simulations (Additional file 1 Text S3, Table S3). Elevation is an essential variable highly correlated with chemical characteristics of both soil and atmosphere. Specifically, for each of the 17 grouped forest types, a climate-and elevationdependent exponential relationship between forest age and biomass C was used to obtain the best model fitted (Additional file 1 Table S4). The best fitted model for each PF type was identified as the model with the highest R 2 and lowest root mean square error (RMSE).

Age-biomass equations
Two age-biomass relationships were generated using unmanaged and managed sites for each PF type. The relationships were implemented in simulations assuming that the PFs were managed differently from 2010 to 2050.

Experimental design and statistical analysis
In this study, we set up simulations to distinguish and quantify the effects of age, management, and climate on PF biomass C change in China from 2010 to 2050 (Additional file 1 Text S2, Table S2). For example, age effect is derived from experiments assuming that the PFs were not subjected to management (NOman, S1), while climate was fixed at the 2010s (Additional file 1 Table S2). More simulation experiments were designed, including business-as-usual management (BAUman, S2), all forest under management (ALLman, S3), and varying climate experiments (Additional file 1 Text S2, Table S2). In total, 8400 simulation experiments were performed [3 treatments (unmanaged, managed, BAU management) × 100 case runs (100 pairs of age and type maps reconstructed) × 8 earth system model × 3 or 4 RCPs = 8400 simulations]. It should be pointed out that, instead of quantifying the impact of a specific management practice, this study highlights the estimates of biomass C accumulations at national scale by assuming the PFs were exposed to current or presumed managements and disturbances regimes. More specifically, the estimated C accumulation is derived from scenario simulations when the PFs were not managed or the national management intensity remained stable (business-as-usual scenario, BAU). Similarly, the C accumulation enhancement potential is the biomass C stock increment from applying current management intensity to the currently unmanaged PFs.

Distribution of simulated biomass C in 2010
The simulated biomass C in PF was 4250.56 ± 107.99 g C•m − 2 for mainland China (1.89 ± 0.048 Pg C in total) based on an average age of 19.71 ± 0.013 years in 2010 (Fig. 1). PF biomass C stock was higher in southern regions than in northern regions, and the highest biomass C stock was found in southeast regions (Fig. 1a). The uncertainty of biomass C is relatively high in south regions, followed by northeast regions, while it was low in central regions (Fig. 1b).

Projections of age, management, and climatic effects on forest biomass C accumulation
Age effect is represented as the PF biomass C increment from experiment simulations with absence of management and climate change (experiment S1 in Additional file 1, Table S2). Generally, the biomass C increments and uncertainties are higher in south and east regions than other regions of China (Fig. 2a). Similarly, the impacts of management on biomass C stock increment are also larger in south and east regions than central, western and north regions (experiment S2 in Additional file 1 Table S2, Fig. 2b). In comparison, climate change impact on biomass C relatively small. The climate effect on PFs for the period 2010-2050 is generally positive, and the national total impacts range between 0.038 to 0.138 Pg C (experiments S4 & S1 in Additional file 1 Table S2, Fig. 2c).

Projections of biomass C accumulation in planted forests
The projected biomass C stock increment from 2010 to 2050 is predominantly contributed from age impact (66.6%, 1.23 ± 0.002 Pg C), followed by management (31.2%, 0.57 ± 0.02 Pg C), while the contribution of climate is only slight (0.087 ± 0.04 Pg C) (Fig. 3). We also found that age, management, and climatic impacts decline quickly from 2010 to 2050. Simulation results also revealed that an additional 0.24 ± 0.07 Pg C could be stored if the unmanaged PFs were managed under current intensity, which is 13% of the C from the BAU management (Fig. 3).
Our simulations predicted that biomass C increment in China's PF would be 1.23-2.13 Pg C from 2010 to 2050, majorly depending on the management (Fig. 4). Specifically, under the BAU scenario (experiment S2, Additional file 1 Table S2), decadal biomass C increment would decrease from 0.75 ± 0.12 to 0.22 ± 0.04 Pg C from 2010 to 2050, leading to a total of 1.80 Pg C stored in biomass (Fig. 4). Assuming the PFs will remain unmanaged (experiment S1), the total biomass C increment would be 1.23 Pg C until 2050 (age effect only). However, a total of 2.04 Pg C can be accumulated if the PFs were all managed (experiment S3, age and management effects plus management potential) (Fig. 4).

Simulated biomass C stock in PF
Our estimation of the PF C stock (1.89 ± 0.048 Pg C) was 16.44%-18.03% of the total forest biomass C in China (10.48-11.49 Pg C, 5568.5-5860 g C•m − 2 ) Xu et al. 2018). This is reasonable as PF was younger and generally had a lower C density compared with natural forest (Additional file 1 Text S4). Besides, our results were derived from simulations with considerable advance over prior studies. An important improvement is that the age, type, and distribution of the PFs have been strictly validated for each major forest type in each province (Yu et al. 2020a). Former studies generally suffered from lack of spatially matched information of forest characteristics (detail types, distributions) and demography (ages). Using mismatched forest age map and forest type map would introduce biases in the projection of forest growth because 1) only grid cells with both age and type present can be simulated, leading to a loss of forest area, and 2) forest age and type information would be inconsistent.
We extrapolated biomass C increment in China's forests (6.71 Pg C in total during 2010-2050, Additional file 1 Table S5) and found that it is considerably lower than the range of 8.89-13.92 Pg C reported previously for the same period (He et al. 2017;Yao et al. 2018), albeit there was an evident overestimation in the extrapolation (Additional file 1 Text S5). This implies that former studies may have greatly overestimated the biomass C increment potential in China's forests by at least 32.5%. Overestimations in previous studies may be attributed to the method and data used. An overestimation source is treating economic forests and bamboo forests as PFs, since the two forests have been classified as PF in the national inventory. Economic forests and bamboo forests, accounting for 10.5%-15.9% of the China's forest, contributed 9.8% (186 Tg) biomass C increment from 1977 to 2008 (Guo et al. 2013). Nonetheless, the increment was almost solely attributed to area expansion, which rose from 1448 to 2579 Mha during the same period. Conversely, the average biomass C density in the two forests barely changed (1761 vs 1710 g C•m − 2 ) from 1977 to 2008, implying a relatively low biomass C increment potential from age effect alone. Thus, overestimations are expected if the two forests are treated the same using age-biomass relationships derived from samples obtained from other forests. Another source of overestimation may be sampling bias. Former studies rarely reported information about the management adopted in the sites used in developing age-biomass relationships. We found that using sites majorly sampled from intensively managed forest farm or timber forests would greatly overestimate forest C accumulation rate (Yu et al. 2020b). Due to the lack of management information, it is practically difficult to assess the samples used in previous studies.
China plans to further expand forest area by 49.49 Mha from 2021 to 2050 through afforestation to increase forest C stock by 3.5 Pg (State Forestry Administration of the People's Republic of China 2016). However, our simulations revealed that this ambitious goal may be achievable by optimizing management of existing forests, thus preventing the use of additional land and potential conflicts with agricultural sectors.

Age and management effects on biomass C accumulation
Impacts of age and management dominate the changes of the projected biomass C stock increment from 2010 to 2050 (Fig. 3). Besides, the simulations also showed that age, management, and climatic impacts decline quickly from 2010 to 2050, which is consistent with other studies (He et al. 2017;Yao et al. 2018). We predicted that forest growth saturation would result into decline in decadal biomass C accumulation from 0.77 ± 0.12 in the 2010s to 0.22 ± 0.04 Pg C in the 2040s, contributing to a total biomass C increment of 1.84 Pg C (0.046 Pg C•yr − 1 ) for the study period. Yao et al. (2018) and Xu et al. (2010) estimated that biomass C accumulation in China's forests would be around 0.14-0.16 Pg C•yr − 1 , which is approximately twice as high as the rate of 0.078 Pg C•yr − 1 reported by Hu et al. (2015). Surprisingly, the estimation of He et al. (2017) is much higher (0.34 Pg C•yr − 1 ) during the same period, being likely overestimated due to no consideration of human management (He et al. 2017). Our simulated PF C increment is about 13.5%-59.0% of the increment of total forests reported. Despite a much lower area of PF compared to natural forest, PF has contributed a relatively high proportion (40.92%) of forest C sink in China during 2004-2008(Guo et al. 2013. Due to growth saturation, the simulated PFs C sink is expected to decline quickly, which should be compensated by human intervention as suggested from our simulations. Simulation results also revealed that an additional 0.24 ± 0.07 Pg C could be stored under all-management condition till 2050, which is 13% of the C from the management under BAU (Figs. 3 and 4). Among all the PF types, P. massoniana, Populus spp., and needle-leaf and broadleaf mixed forests will be the top three that would contribute 55% of the C biomass increment if all forests are managed, although the decadal contributions of the three PFs decrease from 0.040 Pg C in 2010 to 0.014 Pg C in 2050 (Additional file 1 Table S7). Since the effect of climate on biomass C is minor (Fig. 3, Additional file 1 Table S6), we focus on age and management effects hereafter.
Chinese forests have served as a strong C sink from both area expansion and forest growth in the last few decades (Zhao et al. 2019), while area expansion has been a larger contributor to C sink than forest growth in the PF (62.2% vs 37.8%, Li et al. (2016)). Wang et al. (2020) estimated the land biosphere sink at − 1.11 ± 0.38 Pg C•yr − 1 during 2010-2016, in which the progressively forested area in the southwest China (Yunnan, Guangxi and Guizhou provinces) played a pivotal role. Our study provides a similar result; however, the spatial distribution of C sink is not fully consistent with that of Wang et al. (2020). This is because Wang et al. (2020) adopted an atmospheric inversion approach, which captures C assimilated by the land ecosystem, including signals from timber harvested and exported to other regions. Especially for southern China, timber production is an essential industry in the region. In comparison, the biomass growth approach used in this study excludes the C flux from timber production. A former study also revealed that the eight provinces in southern China (i.e. Fujian, Jiangxi, Zhejiang, Hunan, Guangdong, Guangxi, Sichuan, Yunnan), covering 20% of the land area, have contributed to approximately 51% of both area expansion and C stock increment in PFs since the 1970s (Li et al. 2016). This is consistent with our simulations that the future biomass C increment potential locates mainly in the southern region (Fig. 1). Thus, continuous contribution from area expansion requires the availability of suitable land for afforestation, especially lands in the south.
However, the perspectives about PFs duration in C sequestration are divergent in China. Optimistic projections suggested a sustained C sink in PFs for the next few decades due to the young age (Huang et al. 2012;Deng et al. 2017), but a recent study suggested that forest growth saturation, land competition for food production, and soil-water depletion challenge the longevity of C sink (Tong et al. 2020). Our simulations predicted that biomass C increment in China's PF would be 1.23-2.13 Pg C from 2010 to 2050, primarily depending on management (Figs. 3 and 4). Therefore, proper forest management may be a promising approach to enhance C storage in China's PFs. For example, our previous study revealed that biomass C accumulation in P. massoniana, C. lanceolata, and Eucalyptus sites was much faster in the presence of management compared to unmanaged sites (Yu et al. 2020b). Here, we found that an additional 0.24 Pg C can be stored if current PFs are all managed by 2050. However, despite the additional C accumulation available from managements, the economic costs may be disproportionately high, indicating that enhancement of C storage in future PFs would be challenging.

Implications for forest management targeted on biomass C accumulation
We found that the majority (65.3%) of China's PFs were managed (Additional file 1 Fig. S2, Table S1), indicating that great efforts have been made by the government. Nevertheless, the impact of management on forest C accumulation has rarely been quantified in such a large scale and high spatial coverage in China, while its contribution to forest growth is usually attributed to age effect. For example, Yao et al. (2018) estimated that age-related forest biomass C storage would be 6.69 Pg C from the 2000s to the 2040s, which actually included management impacts as the projection was derived from samples managed differently. Here, by dividing collected samples into managed and unmanaged groups, we developed biomass estimation models for each of the major PF types and quantified management impacts separately. Surprisingly, we found that management would be the second largest contributor to forest biomass C accumulation (0.57 ± 0.02 Pg C) from 2010 to 2050, with an additional 42.1% (0.24 ± 0.07 Pg C) made possible if all forests are managed. This finding also implies that endeavors are required for maintaining or even enhancing C sink from China's PFs.
Appropriate planning and management, the endeavors expected to be devoted by government and foresters, should consider practices tailored to local climate and other ecological conditions. Despite a large number of sites surveyed and used in this study, more data are required to identify optimal management practices per PF type and for each location and forest type. Our results suggest that forest management and age-related growth dominate the biomass C change in PFs, while the effect of climatic factors on the accumulation is minor. This implies that management practices and forest demographic information, which have been under-represented in process-based models, are crucial for accurate simulations. Besides, due to the assumption that future forest growth from equilibrium conditions, current processbased models may fail to capture forest regrowth signal and, therefore, underestimate forest biomass C accumulation (Yao et al. 2018). For accurate attribution of historical management practice and projection of their effects on biomass C change in China's PF, more controlled experiments are needed, and the field-derived results should be incorporated into model simulations.
To examine differences in management impacts between regions, we also quantified the sensitivity of biomass C change in the ranges of air temperature and precipitation derived from the collected sites for each PF. The sensitivity was described as biomass C change from the 2010s to the 2040s between managed and unmanaged regimes. We found that management impact in most of the PF types was sensitive to temperature, except for Pinus tabulaeformis and P. elliottii (Additional file 1 Fig. S4). For the major C increment contributor, P. massoniana, management impact on biomass C increment was higher in warmer and wetter regions (Additional file 1 Fig. S4). Regarding Populus and C. lanceolata, management impact on biomass C increment was positive and higher in warmer areas, indicating a higher C storage potential from management in warmer areas. Conversely, management impact on Pinus tabulaeformis and P. elliottii was more sensitive to precipitation, implying a higher C uptake potential from PF management in wetter regions. Thus, although management helps enhancing C uptake in PFs, its impacts will be influenced by climate change.

Uncertainties
There are some uncertainties that may affect the results of this study. First, we did not separate contributions from rising CO 2 , N deposition, and distrubances such as fire and diseases. Since the sites collected have been exposed to these atmospheric, climatic, and other environmental factors, their impacts have been reflected in PFs's growth. Thus, our simulations assume that these factors will remain the same as with the historical period. Yao et al. (2018) revealed that the effect of rising CO 2 concentration on biomass C accumulation would be 25%-47% of the age effect. Similarly, N deposition may promote forest C stock increment, albeit the effects differ between locations and forest types (Nadelhoffer et al. 1999;Thomas et al. 2010; Schulte-Uebbing and de Vries 2018). However, Yu et al. (2019a) found that atmospheric N deposition in China has been stabilized over the past decade due to joint contribution from changes in the socioeconomic structure and vigorous controls in N pollution. Thus, stabilized N deposition combined with other natural disturbances (e.g. pests, fires, droughts) are likely to offset the contribution from rising atmospheric CO 2 , although the net impact is obscure.
Second, China will continue to expand its plantation area for environmental improvement and C sequestration facilitation (Wang et al. 2007;Ouyang et al. 2016;Lu et al. 2018). According to the National Forest Management and Planning (State Forestry Administration of the People's Republic of China 2016), projected forest cover will increase to 26%, with forest C stock increasing to 13 Pg by 2050, suggesting another 49.49 Mha will be forested from 2021 to 2050. Our simulations, by adopting a fixed forest distribution map, suggest that the goal of C stock enhancement is feasible through a better management of current forests alone, without requiring more lands. Although a higher C uptake could be expected, the economic costs of expanding PF area should be examined due to challenges from the rising logistical costs, unfavorable environmental conditions, and fragmented lands (Yu et al. 2019b).
Third, the forest management potential may be lower than estimated because China's PFs are serviceorientated so that management practices and intensity implemented should be tailored to meet the foremost goal. For example, the Three-North Shelter Forest Program launched in the late 1970s aiming to alleviate soil degredation and sandstorm, while the Grain for Green Program was initiated in 2000 to protect croplands in hilly areas by reducing soil erosion and flood risk (Lu et al. 2018). Thus, the management of these forests are different from the aim of timber production, which was fast growth and high yield of forest under intensive management. Despite the differing goals, appropriate planning and management are required to reduce tree mortality while encouraging forest quality enhancement.

Conclusions
To our knowledge, this is the first study to quantify the impact of human management on biomass C accumulation at such a broad scale. Our simulation results revealed that human management is the second largest contributor to PF biomass C increment (46.34% of agerelated growth). To achieve the ambitious goal of "C neutral" before 2060 as promised at the United Nations General Assembly in 2020, China plans to further expand forest area by 49.49 Mha to increase forest C stock by 3.5 Pg from 2021 to 2050. However, our simulations suggest that the goal of C stock enhancement is feasible through better management of current forests alone, without requiring more land expansion. Thus, appropriate planning and management should be stressed for sustaining and enhancing biomass C sequestration in China's PF.