Carbon stocks and changes of dead organic matter in China's forests
Jianxiao Zhu1,
Huifeng Hu2,
Shengli Tao1,
Xiulian Chi1,
Peng Li1,
Lai Jiang1,
Chengjun Ji1,
Jiangling Zhu1,
Zhiyao Tang ORCID: orcid.org/0000-0003-0154-64031,
Yude Pan3,
Richard A. Birdsey4,
Xinhua He ORCID: orcid.org/0000-0002-5570-34545 &
Jingyun Fang1,2
Nature Communications volume 8, Article number: 151 (2017) Cite this article
Carbon cycle
Climate-change ecology
Forests play an important role in global carbon cycles. However, the lack of available information on carbon stocks in dead organic matter, including woody debris and litter, reduces the reliability of assessing the carbon cycles in entire forest ecosystems. Here we estimate that the national DOM carbon stock in the period of 2004–2008 is 925 ± 54 Tg, with an average density of 5.95 ± 0.35 Mg C ha−1. Over the past two decades from periods of 1984−1988 to 2004−2008, the national dead organic matter carbon stock has increased by 6.7 ± 2.2 Tg carbon per year, primarily due to increasing forest area. Temperature and precipitation increase the carbon density of woody debris, but decrease that of litter. Additionally, the woody debris increases significantly with above ground biomass and forest age. Our results can improve estimates of the carbon budget in China's forests and for better understanding of effects of climate and stand characteristics on dead organic matter distribution.
Forests cover ~4200 million hectares (Mha) of the world's land surface1, and 50–90% of the total annual carbon (C) flux of terrestrial ecosystems occurs at the interface between forests and the atmosphere2. As key components of forest ecosystems, woody debris (including fine woody debris (FWD), snags, and logs) and litter represent important C stocks and have a strong influence on the structures and C dynamics of forest ecosystems3. Compared with C stocks in live biomass or soil, the C stock in dead organic matter (DOM, including woody debris and litter) has not been included in analysis of the large-scale C budget of China's forests4,5,6 because of a lack of reliable data from a sufficient number of representative sampling sites7. Thus, comprehensive assessments of the full forest C budget and its components are relatively rare, though such analyses are critically important for understanding how forest ecosystems respond to climate change, disturbances, and land management.
In natural or unmanaged forests, the DOM C stock is determined by the balance between the input and output of DOM, which is affected by various factors. In general, above ground biomass (AGB, including leaves, branches, tree tops, or whole trees) that dies or falls to the forest floor determines the potential input rate of DOM8, 9. Other factors, such as forest productivity, canopy cover, climate, and disturbance, also influence the input and output (decomposition) rates of DOM10. In addition, forest age is an important factor because it determines the accumulation of biomass and coarse woody debris (CWD, including snags and logs) stocks in forest ecosystems at a local scale11, 12.
Few studies have focused on the relationships between the driving factors and C stocks of DOM at the regional scale, except some work has been performed at the local scale11,12,13. Only a few country-scale studies have assessed the impacts of various driving forces on DOM C stocks7, 14, 15. Characterizing the relationships between potential driving factors and DOM C stocks is thus essential for understanding the current DOM C stock and for predicting the future evolution of ecosystem C stocks under different climate change scenarios.
Forests in China range from boreal forests in the northeast to tropical rain forests in the south and southwest, and include almost all of the major forest biomes in the Northern Hemisphere16. Forests are diverse in their species composition and climate conditions, thus providing unique venues for examining the spatial variations of DOM C density (C stock per unit area) and how it is impacted by climatic and biological factors. Although several studies have estimated C stocks of live biomass in China's forests4, 17,18,19, reliable estimates of DOM C stocks are still unavailable for this large and dynamic region. Therefore, a field-based estimate of the C stock of the entire forest sector in China is urgently needed to fulfill information needs that support management and policy decisions.
In this study, we documented estimates of DOM C stocks and changes in China's forests over the past two decades. We also examined the potential biotic and abiotic drivers that determined the spatial distribution of woody debris and litter. Specified regional Random Forest (RF) models for all DOM components were used to estimate the current distribution of DOM C and its change rates over the study period. We showed that the current DOM C stock in China's forests was 925 ± 54 Tg (5.95 ± 0.35 Mg C ha−1), with 429 ± 30 Tg in woody debris and 496 ± 24 Tg in litter. Live biomass was the dominant factor shaping the spatial distribution of the woody debris C, while climate determined the regional variation of the litter C stock. The DOM C stock increased by 6.7 ± 2.2 Tg C per year over the past 20 years, primarily due to increasing forest area.
DOM C estimates
To quantify DOM C stocks in China's forests, a consistent field sampling protocol was used to investigate the C stocks of AGB, soil, woody debris, and litter from 567 field plots (20 × 20 m2) at 189 forest sites across the major forest types and climatic zones in China (Fig. 1, Supplementary Fig. 1, and Supplementary Data 1). The DOM C density in the sampled plots ranged from 2.0 to 16.2 Mg C ha−1 (1 Mg = 106 g) with an overall mean of 6.1 ± 1.9 Mg C ha−1 (Fig. 2). The C density of the different DOM components was 0.0–1.3 Mg C ha−1 for FWD, 0.0–11.7 Mg C ha−1 for CWD (0.0–6.2 and 0.0–10.0 Mg C ha−1 for snags and logs, respectively), and 0.6–7.8 Mg C ha−1 for litter. The mean C density of DOM in different forest types ranged from 5.6 ± 1.9 Mg C ha−1 in conifer forests to 6.7 ± 1.7 Mg C ha−1 in evergreen broadleaf forests (Table 1). However, the C density of each DOM component varied substantially among the different forest types. For instance, the mean C density of woody debris was significantly lower in conifer forests (2.0 ± 1.4 Mg C ha−1) than in evergreen broadleaf forests (4.1 ± 1.6 Mg C ha−1) (F = 19.3, P < 0.001), but the mean C density of litter was significantly higher in conifer forests (3.5 ± 1.3 Mg C ha−1), and broadleaf and conifer mixed forests (3.6 ± 1.4 Mg C ha−1) than in evergreen broadleaf forests (2.0 ± 0.5 Mg C ha−1) (F = 12.1, P < 0.001; Table 1). For different origins of forests, primary forests had a significantly higher DOM C density (7.3 ± 2.4 Mg C ha−1) than secondary forests (5.8 ± 1.8 Mg C ha−1) and plantations (5.0 ± 1.6 Mg C ha−1) (F = 67.1, P < 0.001; Supplementary Fig. 2). However, the litter C density did not exhibit any significant differences among these three origins of forests (F = 1.1, P = 0.34).
Distribution of 189 sampling sites in this study. The sites include 128 natural forests (65 primary forests and 63 secondary forests), and 61 plantations. For details, see Supplementary Data 1
Frequency distribution of the carbon density of each component of dead organic matter (DOM). Frequency distributions are shown for a–d the primary forests, e–h secondary forests, and i–l plantations. Components of dead organic matterd, h and l include fine woody debris (FWD; a, e and i), coarse woody debris (CWD; b, f and j), and litter c, g and k
Table 1 Area, number of sites, and carbon density of each carbon component for the different forest typesa
Compared with the C density of AGB (64.5 ± 37.6 Mg C ha−1) and soil (68.1 ± 26.3 Mg C ha−1), the C density of DOM only accounted for 4.6 ± 1.5% of the total ecosystem C density (not including below ground biomass in this study) in the sampled sites, with 2.5 ± 1.4 and 2.1 ± 1.0% in litter and woody debris, respectively (Supplementary Fig. 3). The ratio of litter to total ecosystem C density was significantly different among the different forest types (F = 17.0, P < 0.001); the highest was in conifer forests (3.1 ± 1.3%) and the lowest was in evergreen broadleaf forests (1.2 ± 0.6%). However, the ratio of woody debris to total ecosystem C density was significantly lower in conifers (1.7 ± 1.0%) and broadleaf-conifer mixed forests (2.0 ± 0.8%) compared with that in deciduous broadleaf forests (2.2 ± 1.1%) and evergreen broadleaf forests (2.6 ± 0.7%) (F = 6.5, P < 0.001; Supplementary Fig. 3). The ratio of woody debris to total ecosystem C density in primary and secondary forests was similar (2.5 ± 1.1 vs. 2.2 ± 0.9%), and it was significantly higher than that in plantations (1.6 ± 0.8%) (F = 13.6, P < 0.001; Supplementary Fig. 3).
Impacts of biotic and climatic factors on DOM C density
Both FWD and CWD C densities were significantly positively correlated with AGB, forest age, and net primary production (NPP), but litter C density was negatively correlated with these factors (Fig. 3). Similarly, climatic factors, indicated by mean annual temperature (MAT) and precipitation (MAP), had significant positive effects on the C density of FWD and CWD (P < 0.001) but negative effects on the C density of litter (P < 0.001) (Fig. 3). In addition, the soil C density was positively correlated with the C density of CWD and litter (P < 0.05), but not withFWD.
Relationships between carbon density of dead organic matter, and biotic and climatic factors. The y axis is carbon density (Mg C ha−1) of a–f fine woody debris (FWD), g–l coarse woody debris (CWD), and m–r litter. All plots showed the linear relationships between a, g and m above ground biomass carbon density, b, h and n forest age, c, i and o net primary production (NPP), d, j and p mean annual temperature (MAT), and e, k and q precipitation (MAP), and (f, l and r) soil carbon density and carbon density of each component of dead organic matter. There were no significant linear relationships between carbon density of litter and forest age, and between carbon density of FWD and carbon density of soil
DOM C stocks in different inventory periods
Because there are significant correlations between the C density of each DOM component, and climate and stand characteristics (Fig. 3), we used these biotic and abiotic variables as predictors to estimate the spatial and temporal distributions of C density for each DOM component across China's forests. On the basis of the forest inventory and climatic data for the different inventory periods, from 1984−1988 to 2004−2008, we first used a specified regional RF model for each DOM component to map the grid-cell (at a resolution of 0.083°) C density of each DOM component across China's forest area for each inventory period (Supplementary Fig. 4). Then, we estimated the corresponding C density and C stock of each DOM component for each region and the whole country using the forest area for each inventory period (Table 2, Supplementary Table 1, see Methods).
Table 2 Modeled carbon density and stock of the forest dead organic matter in China for the forest inventory periods of 1984–1988 and 2004–2008. For information about the carbon density and stock during the other inventory periods, see Supplementary Table 1
On the basis of these calculations, the current (2004−2008) national DOM C density averaged 5.95 ± 0.35 Mg C ha−1 and showed a small variation among the different regions, ranging from 5.53 ± 0.36 Mg C ha−1 in southeastern forests to 6.28 ± 0.31 Mg C ha−1 in northern forests. However, the C density among the DOM components differed significantly, with woody debris accounting for >50% of the total DOM C density in southwestern, south central, and southeastern forests, while litter accounted for >60% of the total DOM C density in northwestern, northern, and northeastern forests (Fig. 4). The current national DOM C stock was estimated at 925 ± 54 Tg C (1 Tg = 1012 g), of which 429 ± 30 Tg C was stored in woody debris and 496 ± 24 Tg C was stored in litter (Table 2). Geographically, the largest DOM C stock (249 ± 17 Tg C, 27%) was in southwestern forests, followed by northeastern (180 ± 8 Tg C, 19%), south central (175 ± 13 Tg C, 19%), northern (137 ± 7 Tg C, 15%), and southeast (123 ± 8 Tg C, 13%) forests. The smallest DOM C stock was in northwestern forests (59 ± 3 Tg C, 6%) (Fig. 4).
Carbon stocks of the woody debris and litter in the different regions of China from periods 1984–1988 to 2004–2008. The filled colors represent the model-based estimates of dead organic matter carbon density. Carbon stocks are presented as means ± 1 SD, generated by 500 resamplings. Forests in China are divided into six regions: Northwest, including the provinces of Gansu, Ningxia, Qinghai, Shaanxi, and Xinjiang; North, including the provinces of Hebei, Inner Mongol and Shanxi, and the cities of Beijing and Tianjin; Northeast, including Heilongjiang, Jilin, and Liaoning provinces; Southwest, including the provinces of Guizhou, Sichuan, Yunan and Tibet, and the city of Chongqing; South Central, including the provinces of Guangdong, Guangxi, Hainan, Henan, Hubei and Hunan; and East, including the provinces of Anhui, Fujian, Jiangsu, Jiangxi, Shandong and Zhejiang, and the city of Shanghai. It should be noted that due to the lack of data, forests in Hong Kong, Macao and Taiwan were not included in the current study
On the basis of the corresponding C density and C stock of each DOM component for each inventory period, we estimated their changes between 1984−1988 and 2004−2008 (Fig. 5). The C density of woody debris has significantly increased at a rate of 0.008 ± 0.002 Mg C ha−1 per year (~0.3% per year) over the past two decades (R 2 = 0.91, P = 0.01), whereas the C density of litter has significantly decreased at a rate of 0.003 ± 0.001 Mg C ha−1 per year (~0.1% per year; R 2 = 0.87, P = 0.02). As a result, the total DOM C density has increased by 0.005 ± 0.002 Mg C ha−1 per year (~0.1% per year; R 2 = 0.77, P < 0.05) over the two decades (Fig. 5a). Geographically, the largest increase in DOM C density (0.020 ± 0.007 Mg C ha−1 per year) occurred in the northern region with DOM C density increasing at a rate of 0.014 ± 0.004 Mg C ha−1 per year in woody debris and 0.006 ± 0.004 Mg C ha−1 per year in litter (Supplementary Fig. 5). Consistent with DOM C density, the national DOM C stock has continuously increased at a rate of 6.7 ± 2.2 Tg C per year (0.9% per year) through the past 20 years (Fig. 5b), primarily due to an increase in forest area (1.2 Mha per year, 0.9% per year; Supplementary Table 1). In total, 58% (3.9 ± 0.9 Tg C per year) of the increase in DOM C stock was attributed to woody debris and 42% (2.8 ± 1.3 Tg C per year) to litter.
Changes in dead organic matter carbon over the past two decades. Changes in dead organic matter a carbon density and b carbon stock in China's forests from 1984–1988 to 2004–2008. Carbon densities and stocks are presented as means ± 1 SD, generated by 500 resamplings. The fitted lines are from linear regression. Significant fitted lines (P < 0.05) are shown as the solid lines, and non-significant fitted lines (P > 0.05) are shown as the dotted lines
Previous studies have investigated the spatial patterns of C stocks in live biomass and soils in China's forests4, 5, 20, but few studies have estimated the C stocks of full forest ecosystem components due to the lack of information on woody debris and litter. In this study, we quantified the C stock in woody debris and litter, as well as in AGB and soil in China's forests, to reduce uncertainties when assessing the C budget for forest ecosystem5, 14. Our estimates of the C density of woody debris (2.61 ± 0.18 Mg C ha−1) and litter (3.26 ± 0.17 Mg C ha−1) are very different particularly in components to those reported in Pan et al.21 (0.6 and 7.7 Mg C ha−1). These differences may be because their analyses were simply based on the data of live biomass of China's forest inventory reports that lacked direct field data on woody debris and litter.
Compared with temperate forests in other regions, the DOM C density is relatively low in China (Fig. 6). For instance, our estimate of 5.88 ± 0.35 Mg C ha−1 was less than half that in the U.S. (15–29 Mg C ha−1) and Europe (11–16 Mg C ha−1)21,22,23. The low DOM C density may be caused by low input rates and/or high decomposition rates. Started in the later 1970s, China's reforestation and afforestation programs have led to the large area of young-growth plantations with less than 40 years old in China (~40 Mha, 26% of total forest area)24, which is more than 10 times of that in USA (3 Mha) and Europe (1 Mha)23. This large area of young-aged plantations and excessive harvest occurred before the reforestation programs have lower densities of live biomass and relatively shorter time of the accumulation of woody debris, and therefore cause a lower DOM C input17, 24, 25. In addition, it does not exclude the fact that DOM investigations in other countries tended to represent older, undisturbed forests with large woody debris stocks21 or forests that were recently affected by catastrophic disturbances26, 27, which resulted in higher estimates of the detrital C stock.
Carbon density of each forest ecosystem component in China compared with the values reported in other regions. a The pie charts show the contribution (%) of dead organic matter, live biomass and soil to total ecosystem carbon stock. b Carbon densities (Mg C ha−1) of the woody debris, litter, live biomass, and soil in China and other regions. The carbon density of Europe (geographic Europe), the United States, and the world are based on data reported in refs. 21,22,23, 28. Carbon densities in this study are presented as means ± 1 SD, generated by 500 resamplings, and those in other studies are presented as means ± 1 SD from their reports (n = 3 for China, Europe and the USA; n = 4 for the world)
Comparing to the overall mean C density of AGB and soil in China's forests (i.e. 64.5 ± 37.6 and 68.1 ± 26.3 Mg C ha−1, respectively) (Table 1 and Fig. 6), the overall proportion of DOM C (6.1 ± 1.9 Mg C ha−1) was fairly small since being only 4.6 ± 1.5% of the entire ecosystem C, which is comparable to that reported by Pan et al.21 (5.4%) for China's forests but much lower than that in other temperate forests (8–18%)21,22,23, 28. Our current estimate was also confirmed by several other studies that have investigated the proportion of DOM C at the local scale: 4.4% in a subtropical evergreen broadleaf forest29, 4.0% in temperate deciduous broadleaf forests30, and 4.4% in a boreal forest31 in China. The small proportion of DOM C stock is likely a result of younger forests in China, which have lower DOM C stock and also incline to exempt intensive disturbances (e.g., wildfire), compared to temperate forests in other regions. In addition, the ratio of woody debris to AGB was small (4.9 ± 2.7%), much lower than that reported by Harmon and Hua32 (20–25% for subtropical and temperate forests). Woody debris only accounted for 2.1 ± 1.0% of the total C in China's forests, comparing to other temperate forests (4–18%) in the Northern Hemisphere33. It suggests that using empirical ratios of the global data would greatly overestimate woody debris C stock in China's forests.
Next, we discuss how stand characteristics and climatic factor influence DOM C density. In general, the rates of input (biomass production) and output (decomposition) determine the C density of woody debris and litter in natural or unmanaged forests. However, because the C of woody debris is persistent and slow to break down, its input rates ultimately outpace the decomposition rates. We found that woody debris was positively correlated with AGB, consistent with the finding from US forests in which the spatial distribution of downed woody debris was directly determined by the amount of live biomass14. Forest age was also expected to affect the amount of woody debris at the local scale, and often a greater amount of woody debris is found in young or old growth stands than in other stands ('U- shaped' pattern)34. However, our study suggests a monotonic increase of woody debris C density with forest age at the regional scale (Fig. 3). A study of forests in the US also reported similar results14. High stocks of woody debris commonly occur in mature forests that are highly stocked with live biomass35.
Besides the significant effect of stand characteristics, MAT and MAP also influenced positively FWD and CWD C density, although a US study found a negative trend in forests7. If climatic factors drive the distribution of forest live biomass36, 37, the potential amount of biomass input ultimately determines the spatial distribution of woody debris38, 39. The factor importance analysis (see Supplementary Materials and Methods) indicated that AGB was the optimal factor for predicting the spatial distribution of C density in all woody debris components (Supplementary Fig. 6). In contrast to the woody debris distribution, MAT and MAP were negatively correlated with litter distribution, which was likely related to faster losses of litter under hotter and wetter climates, which speed decomposition and turnover time40.
How much C has accumulated in China's forest DOM over the past two decades? As shown in Fig. 5, the DOM density has increased by 0.005 ± 0.002 Mg C ha−1 per year in China's forests over the study period. The enhanced C accumulation in vegetation biomass in China's forests4, 24, 41 could increase potential DOM C inputs. Environmental changes, including increases in annual temperature and summer precipitation may contribute to the DOM C inputs because they could stimulate forest biomass C gain42, but they could also cause DOM C loss as a result of increased decomposition rate, especially of litter. In addition, large-scale reforestation and afforestation practices since the late 1970s have led to a significant increase in the forest area from 131.7 Mha in 1984–1988 to 155.6 Mha in 2004–2008, with an expansion rate of 0.9% per year. Compared with the small increased rate of DOM C density (0.1% per year), the expanding forest area (Supplementary Table 1) might be a dominant contributor to the accumulation in DOM C stocks for China's forests (Fig. 5).
According to previous estimates, the C sinks of forest biomass24 and soils20 over the past two decades were 70.9 and 67.2 Tg C per year, respectively. Together with the estimated C sink of 6.7 ± 2.2 Tg C per year in DOM in the current study, China's forests could have sequestered a total of ~ 145 Tg C per year from the 1980s to 2000s, which is roughly equivalent to 15.5% of the contemporary (1986–2006) fossil CO2 emissions in the country43.
Finally, we evaluated the uncertainty of the estimated C density of each DOM component at both the plot and pixel levels. At the plot level, the coefficient of determination (R 2) and root mean square error (RMSE) for each DOM component were calculated from the validation plot measurements (Supplementary Table 2). The averaged R 2 for woody debris and litter was 0.53 and 0.46, and the averaged RMSE for these components was 1.21 and 0.95 Mg C ha−1, respectively. At the pixel-scale, spatial distribution of the standard deviation (SD) for each DOM component was calculated using a Monte Carlo simulation (Supplementary Fig. 7).
Several uncertainties of our estimates of C stocks in woody debris, litter, AGB, and soil in China's forests must be acknowledged. First, relative to the large forest area (153 Mha) in China, the relatively small number of observed forest sites (n = 189 with 567 plots) could lead to imprecise estimates, although we put substantial effort to ensure high-quality data from field measurements representing the main forest biomes of China. Second, estimates excluded areas with logging residue piles from logging camps in China, which could result in an underestimation of national woody debris C stock. Lastly, seasonal variability could lead to uncertainty in the measurement of litter C stock. These biases in the field survey could have consequences of either an overestimation or an underestimation of DOM C density and may have introduced uncertainties in our results.
In conclusion, our study is the first national estimation of the C stock and the spatial distribution of DOM C density in China's forests. Our results indicate that the total DOM C stock in China's forests was 925 ± 54 Tg C (the area-weighted mean of 5.95 ± 0.35 Mg C ha−1), with 429 ± 30 Tg C in woody debris and 496 ± 24 Tg C in litter. The estimated DOM C stock showed a net gain of 6.7 ± 2.2 Tg C per year over the past 20 years. However, the ratio of DOM C to the overall forest C was as low as 4.6 ± 1.5%. With the ongoing national reforestation and afforestation programs in China, the forest DOM C stock is likely to continue to increase and should positively contribute to the entire forest C sink in the future.
The C stocks of all ecosystem C sectors (AGB, soil organic matter, woody debris, and litter) in China's forests were estimated from 567 sampling plots with an area of 20 × 20 m2at 189 forest sites (three replicated plots at each site). These forest sites range from 18.7 to 52.8 °N in latitude, from 81.0 to 131.2 °E in longitude, and from 30 to 3846 m above sea level, and represent almost all forest types in China, from the tropical rain forests in Hainan Island and Xishuangbanna to the boreal forests in the Greater Khingan Mountains (Fig. 1). MAT varies from −3.8 to 25.7 °C, and MAP ranges from 163 to 1850 mm across these sites, which almost covers the whole climate space in the country (Supplementary Fig. 1). Among these sites, there are 65 primary forests, 63 secondary forests, and 61 plantations (Supplementary Data 1).
Field investigation
The fieldwork was conducted between 2011 and 2016 (Supplementary Data 1). According to Harmon et al.3, DOM was divided into three categories: (1) CWD, including standing snags and fallen logs, and defined as dead wood with a diameter of ≥10 cm at the maximum section, (2) FWD, defined as dead wood with a diameter of 2–10 cm, and (3) litter (defined as dead wood with a diameter of <2 cm and leaf litter). All these DOM components were collected from each forest plot and weighed. Some dead wood could not be directly weighed (i.e., logs with a diameter of ≥25 cm diameter and snags). For those dead woods, we first numbered, marked, and measured the length and diameter of the middle and both ends for the fallen logs, and measured the height and diameter at breast height (1.3 m, DBH) for the snags. We then weighed three parts of the log (the middle and both ends, three 10–20 cm length cylinders) after oven drying at 85 °C to a constant weight44. Next, we applied the mean ratio of the oven-dried weight to the volumes of the three woody cylinders to the total volume of the wood that could not be directly weighed to estimate its biomass. We also extrapolated the ratio to other large woody debris of the same species and determined the decay status to estimate the biomass.
The degree of CWD decay was defined in the following four classes: (I) recently fallen logs with wood and bark intact; (II) initial sapwood decay, bark loss, or easily removed bark and sound heartwood; (III) bark loss, heartwood decay, and large branch loss with stubs present; and (IV) rotten heartwood45, 46. The biomass of large snags was also extrapolated using the ratio of the weight to the volume of fallen logs in class (I) or (II). In addition, the volumes of individual snags and logs were calculated according to the following Eq. (1):
$${\rm{Volume}} = \frac{{\pi {d^2}L}}{4}$$
where d is the average diameter (middle and both ends) of the logs or the DBH for the snags and L is the length of the logs or the height of the snags.
The litter survey was performed at five subplots (2 × 2 m2, a total of 15 subplots for each site), and all the collected fine litter (including fresh and semi-decomposed leaf litter), reproductive structures, fallen bark, small dead wood <2 cm diameter (at the widest part) and other dead plant material (e.g. fallen leaves) in the surface litter layer were weighed after oven drying at 65 °C to a constant weight. The C concentrations of litter and woody debris in the four decay classes are shown in Supplementary Table 3.
In addition, we measured the height and DBH of all trees with a DBH of ≥5 cm in each plot. At each site, the ages of the 10 trees with the largest DBH were determined by performing a tree-ring analysis47. The age of the fifth-largest tree (DBH) at each site was used as the representative of the forest stand age48. The AGB of each tree and each plot was calculated using allometric equations (Supplementary Table 4) in relation to biomass components (leaf, bark, branch, and stem) according to the DBH (cm) and height (m) by species and regions.
We also sampled 1692 replicated soil profiles from 564 plots (three soil profiles from each plot) to estimate soil C density at the 188 forest sites across China (note: the soil data were not available at the Shennongjia site). Soils were separately sampled at depths of 0–10, 10–20, 20–30, 30–50, 50–70, and 70–100 cm. The soil bulk density at each depth was estimated using a standard container (100 cm3, 50.5 mm in diameter and 50 mm in height). The soil gravimetric moisture was measured after oven drying for 48 h at 105 °C. The soil subsamples were air dried for 2 weeks at room temperature (~25 °C) after the removal of debris and plant materials and then ground and sieved (0.15 mm) for C analysis. The C concentrations of DOM and soil were determined using an Elemental Analyzer (2400 II CHN Elemental Analyzer; Perkin-Elmer, Boston, MA, USA), and a factor of 0.5 was used to convert the live biomass to C49, 50. Thus, we obtained the C density of AGB, soil, FWD, CWD (snags and logs), and litter for each of the investigated sites (shown in Supplementary Data 1).
Biomass and climate data
National forest inventory data were obtained to determine the forest area and timber volume for each forest type for each province. The inventory data in five periods of 1984–1988, 1989–1993, 1994–1998, 1999–2003, and 2004–2008 were used in this study51,52,53,54,55. The following previously determined biomass expansion factor (ratio of biomass to volume) was used to convert forest timber volume to biomass4, 5, 24, 50:
$${\rm{BEF = }}a + \frac{b}{x}$$
where x is the timber volume per unit area (m3 ha−1 in this study) and a and b are constants for a specific forest type. We used Eq. (2) to calculate the provincial forest biomass for each forest type for each period based on the corresponding area and volume in the forest inventory data. Next, a factor of 0.5 was used to convert the biomass to C49, 50.
The digitized 1:1,000,000 vegetation map56 was used to obtain the spatial distribution of national vegetation for 161 forest types (excluding bamboo forests). Eighteen forest types defined from the forest inventory were matched within the 161 forest types defined using the digitized vegetation map and then used for the calculations (Supplementary Table 5). The grid-cell AGB C density for each forest type in each province was the area-weighted averaged value, which was used to obtain forest biomass C density distribution during different inventory periods (Supplementary Fig. 8).
The monthly air temperature and precipitation data were obtained from 728 meteorological stations (National Meteorological Information Center of the China Meteorological Administration, http://data.cma.cn 57. We interpolated the climate data to grid cells with a resolution of 0.083° using a Kriging interpolation algorithm. We obtained MAT and MAP for each site and each period.
In addition, we documented mean net primary productivity (NPP) from a 10-year average (2000–2009) of a Moderate Resolution Imaging Spectroradiometer (MODIS) NPP product58, 59 for the analysis of input C to DOM stocks. The mean NPP is the summation of 8-day calculations of net photosynthesis minus maintenance and growth respiration and had a nominal 0.083°resolution. For MAT, MAP, and NPP data, we extracted site-specific values from dataset according to longitude and latitude of the sampling sites.
We constructed RF models60 to simulate the national grid-cell C density of woody debris (FWD, snags, and logs) and litter. The RF fitted multiple decision trees to input data using a random subset of the input variables for each tree constructed for each response variable. Four factors (MAT, MAP, AGB C density, and forest type) were used as input variables to predict the DOM C density. The data set of each DOM component was randomly divided into two parts: the training dataset (70% of all sites, 132 sites) and the validation data set (30% of all sites, 57 sites). The training data set along with the four predicted factors was used to estimate the C density of each DOM component. We then acquired the coefficient of determination (R 2) and root mean square error (RMSE) between the predicted and measured values from the validation data set in each re-sampling. The training and evaluation procedures were repeated 500 times for each DOM component. As a result, we acquired 500 results for model estimation and evaluation. The R 2 and RMSE in Supplementary Table 2 represent the average value of 500 model simulations. To evaluate the performance of the predicted factors, we explored the relative importance of each variable in a RF model by calculating the percentage increase in the mean square error (%IncMSE) (Supplementary Fig. 6). We then used a Monte Carlo approach to propagate the model errors to pixel-, regional-, and national-scale estimates. Next, we applied the 500 RF models to estimate the C density of all the DOM components for the periods of 1984–1988, 1989–1993, 1994–1998, 1999–2003, and 2004–2008 according to the historical forest types, biomass data (Supplementary Fig. 8), and the related MAT and MAP data. The means and standard deviations (SDs) of the grid-cell, regional, and national C density of each DOM component were generated by all the random models. The final predictions for each DOM component and each period were the average estimates of the 500 RF models. On the basis of the results of the RF models, the spatial distribution of the national FWD, snags, logs, CWD, woody debris, litter, and total DOM C density and their pixel-scale SDs were mapped at a resolution of 0.083°.
To capture the change rate of the DOM C density and stock for each component, we generated a simple linear regression between the C density/stock and the year (1986, 1991, 1996, 2001 and 2006). The slope of the regression and its SD were defined as the change rate and the SD for each DOM C density and stock. The observed data in this study are listed in Supplementary Data 1.
Basic information on locations, forest type, stand age, dominant tree species, and the data for C density of each forest sector of the study sites are listed in Supplementary Data 1. The input data (mean annual temperature and precipitation, forest type, and forest biomass) used can be found in the corresponding references. The remaining data that support the findings of this study can be requested from the corresponding author.
Bonan, G. B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).
ADS CAS Article PubMed Google Scholar
Beer, C. et al. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science 329, 834–838 (2010).
Harmon, M. E. et al. Ecology of coarse woody debris in temperate ecosystems. Adv. Ecol. Res. 15, 133–302 (1986).
Fang, J., Chen, A., Peng, C., Zhao, S. & Ci, L. Changes in forest biomass carbon storage in China between 1949 and 1998. Science 292, 2320–2322 (2001).
Fang, J. et al. Forest biomass carbon sinks in East Asia, with special reference to the relative contributions of forest expansion and forest growth. Glob. Change Biol. 20, 2019–2030 (2014).
Piao, S. et al. The carbon balance of terrestrial ecosystems in China. Nature 458, 1009–1013 (2009).
Woodall, C. & Liknes, G. Relationships between forest fine and coarse woody debris carbon stocks across latitudinal gradients in the United States as an indicator of climate change effects. Ecol. Indic . 8, 686–690 (2008).
Tritton, L. M. Dead Wood in the Northern Hardwood Forest Ecosystem (Yale University, 1980).
Sturtevant, B. R., Bissonette, J. A., Long, J. N. & Roberts, D. W. Coarse woody debris as a function of age, stand structure, and disturbance in boreal Newfoundland. Ecol. Appl. 7, 702–712 (1997).
Bond-Lamberty, B. et al. Disturbance legacies and climate jointly drive tree growth and mortality in an intensively studied boreal forest. Glob. Change Biol. 20, 216–227 (2014).
Hall, S., Burke, I. & Hobbs, N. Litter and dead wood dynamics in ponderosa pine forests along a 160-year chronosequence. Ecol. Appl. 16, 2344–2355 (2006).
Carmona, M. R., Armesto, J. J., Aravena, J. C. & Pérez, C. A. Coarse woody debris biomass in successional and primary temperate forests in Chiloé Island, Chile. For. Ecol. Manage. 164, 265–275 (2002).
Jönsson, M. T., Fraver, S. & Jonsson, B. G. Spatio-temporal variation of coarse woody debris input in woodland key habitats in central Sweden. Silva Fenn. 45, 957–967 (2011).
Woodall, C. W. et al. Biomass and carbon attributes of downed woody materials in forests of the United States. For. Ecol. Manage. 305, 48–59 (2013).
Kurz, W. A. et al. CBM-CFS3: a model of carbon-dynamics in forestry and land-use change implementing IPCC standards. Ecol. Model. 220, 480–504 (2009).
Fang, J. et al. Forest community survey and the structural characteristics of forests in China. Ecography 35, 1059–1071 (2012).
Fang, J. Y., Guo, Z. D., Piao, S. L. & Chen, A. P. Terrestrial vegetation carbon sinks in China, 1981–2000. Sci. China Earth Sci. 50, 1341–1350 (2007).
Wang, X., Feng, Z. & Ouyang, Z. The impact of human disturbance on vegetative carbon storage in forest ecosystems in China. For. Ecol. Manage. 148, 117–123 (2001).
Zhang, C. et al. China's forest biomass carbon sink based on seven inventories from 1973 to 2008. Clim. Change 118, 933–948 (2013).
Yang, Y. et al. Increased topsoil carbon stock across China's forests. Glob. Change Biol. 20, 2687–2696 (2014).
Pan, Y. et al. A large and persistent carbon sink in the world's forests. Science 333, 988–993 (2011).
Goodale, C. L. et al. Forest carbon sinks in the Northern Hemisphere. Ecol. Appl. 12, 891–899 (2002).
Dixon, R. K. et al. Carbon pools and flux of global forest ecosystems. Science 263, 185–190 (1994).
Guo, Z. D., Hu, H. F., Li, P., Li, N. Y. & Fang, J. Y. Spatio-temporal changes in biomass carbon sinks in China's forests from 1977 to 2008. Sci. China Life Sci. 56, 661–671 (2013).
Piao, S., Fang, J., Zhu, B. & Tan, K. Forest biomass carbon stocks in China over the past 2 decades: estimation based on integrated inventory and satellite data. J. Geophys. Res. Biogeosciences 110, G01006 (2005).
ADS Article Google Scholar
Nalder, I. A. & Wein, R. W. Long-term forest floor carbon dynamics after fire in upland boreal forests of western Canada. Glob. Biogeochem. Cy 13, 951–968 (1999).
Domke, G. M., Woodall, C. W., Walters, B. F. & Smith, J. E. From models to measurements: comparing downed dead wood carbon stock estimates in the US forest inventory. PLoS ONE 8, e59949 (2013).
ADS CAS Article PubMed PubMed Central Google Scholar
Harmon, M. E., Krankina, O. N., Yatskov, M. & Matthews, E. in Assessment Methods for Soil Carbon (eds Lai, R., Kimble, J.&Stewart, B. A.) 533–552 (CRC Press, 2001).
Chen, G. S. et al. Conversion of a natural broad-leafed evergreen forest into pure plantation forests in a subtropical area: effects on carbon storage. Ann. For. Sci. 62, 659–668 (2005).
Zhu, J. X. et al. A significant carbon sink in temperate forests in Beijing: based on 20-year field measurements in three stands. Sci. China Life Sci. 58, 1135–1141 (2015).
Hu, X. et al. Impacts of fire severity and post-fire reforestation on carbon pools in boreal larch forests in Northeast China. J. Plant Ecol. 9, 1–9 (2016).
Harmon, M. E. & Hua, C. Coarse woody debris dynamics in two old-growth ecosystems. Bioscience 41, 604–610 (1991).
Vogt, K. Carbon budgets of temperate forest ecosystems. Tree. Physiol. 9, 69–86 (1991).
Coomes, D. A. & Allen, R. B. Mortality and tree-size distributions in natural mixed-age forests. J. Ecol. 95, 27–40 (2007).
Smith, J. E. et al. Methods for calculating forest ecosystem and harvested carbon with standard estimates for forest types of the United State. Report No. NE-343 (USDA Forest Service; 2006).
Reich, P. B. et al. Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc. Natl Acad. Sci. USA 111, 13721–13726 (2014).
Keith, H., Mackey, B. G. & Lindenmayer, D. B. Re-evaluation of forest biomass carbon stocks and lessons from the world's most carbon-dense forests. Proc. Natl Acad. Sci. USA 106, 11635–11640 (2009).
Aakala, T. Coarse woody debris in late-successional Picea abies forests in northern Europe: variability in quantities and models of decay class dynamics. For. Ecol. Manage. 260, 770–779 (2010).
Kiyono, Y. et al. Carbon stock estimation by forest measurement contributing to sustainable forest management in Cambodia. Jarq-Jpn. Agr. Res. Q 44, 81–92 (2010).
Berg, B. & McClaugherty, C. Plant litter: decomposition, humus formation, carbon sequestration (Springer, 2007).
Li, P. et al. The relative contributions of forest growth and areal expansion to forest biomass carbon. Biogeosciences 13, 375–388 (2016).
Fang, J., Piao, S., He, J. & Ma, W. Increasing terrestrial vegetation activity in China, 1982–1999. Sci. China Life Sci. 47, 229–240 (2004).
Zheng, T., Zhu, J., Wang, S. & Fang, J. When will China achieve its carbon emission peak? Natl Sci. Rev. 3, 8–15 (2016).
Grove, S. J. Extent and composition of dead wood in Australian lowland tropical rainforest with different management histories. For. Ecol. Manage. 154, 35–53 (2001).
Sollins, P. Input and decay of coarse woody debris in coniferous stands in western Oregon and Washington. Can. J. For. Res. 12, 18–28 (1982).
Marra, J. L. & Edmonds, R. L. Coarse woody debris and forest floor respiration in an old-growth coniferous forest on the Olympic Peninsula, Washington, USA. Can. J. For. Res. 24, 1811–1817 (1994).
Worbes, M., Staschel, R., Roloff, A. & Junk, W. J. Tree ring analysis reveals age structure, dynamics and wood production of a natural forest stand in Cameroon. For. Ecol. Manage. 173, 105–123 (2003).
Bruelheide, H. & Schmid, B. Community assembly during secondary forest succession in a Chinese subtropical forest. Ecol. Monogr. 81, 25–41 (2011).
Myneni, R. B. et al. A large carbon sink in the woody biomass of northern forests. Proc. Natl Acad. Sci. USA 98, 14784–14789 (2001).
Fang, J., Oikawa, T., Kato, T., Mo, W. & Wang, Z. Biomass carbon accumulation by Japan's forests from 1947–1995. Glob. Biogeochem. Cy. 19, GB2004 (2005).
Chinese Ministry of Forestry. Forest Resource Statistics of China (1989–1993) (Department of Forest Resource and Management, Chinese Ministry of Forestry, Beijing, China, 1994).
Chinese Ministry of Forestry. Forest Resource Report of China-The 7th National Forest Resources Inventory (China Forestry Publishing House, Beijing, China,2010).
Editorial Board of Vegetation Map of China. Vegetation Atlas of China (in Chinese) (Science Press, Beijing, 2001).
Zhao, X. et al. Satellite-indicated long-term vegetation changes and their drivers on the Mongolian Plateau. Landscape Ecol. 30, 1599–1611 (2015).
Pan, Y. et al. Improved estimates of net primary productivity from MODIS satellite data at regional and local scales. Ecol. Appl. 16, 125–132 (2006).
Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).
Article MATH Google Scholar
Zhao, M. & Running, S. W. Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 329, 940–943 (2010).
This work is in memory of Mr. Junyu Hu who performed extensive field work for this study. We thank team members from Peking University, including XinyuXiong, Tianli Zheng, QiongCai, Xuli Zhou, Wenjing Fang, Heng Zhang and Xingxing Jiang for assistance in field sampling between 2011 and 2016. This work was jointly supported by National Natural Science Foundation of China (31621091 and 31330012), Chinese Academy of Sciences (XDA05050000), and the US Forest Service (07-JV-11242300-117).
Department of Ecology, College of Urban and Environmental Science, and Key Laboratory for Earth Surface Processes of the Ministry of Education, Peking University, Beijing, 100871, China
Jianxiao Zhu, Shengli Tao, Xiulian Chi, Peng Li, Lai Jiang, Chengjun Ji, Jiangling Zhu, Zhiyao Tang & Jingyun Fang
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
Huifeng Hu & Jingyun Fang
US Department of Agriculture Forest Service, Newtown Square, PA, 19073, USA
Yude Pan
Woods Hole Research Center, Falmouth, MA, 02540, USA
Richard A. Birdsey
Centre of Excellence for Soil Biology, College of Resources and Environment, Southwest University, Chongqing, 400715, China
Xinhua He
Jianxiao Zhu
Huifeng Hu
Shengli Tao
Xiulian Chi
Peng Li
Lai Jiang
Chengjun Ji
Jiangling Zhu
Zhiyao Tang
Jingyun Fang
J.F.: Designed the research; Jx.Z., X.C., P.L., L.J., and C.J.: Performed the field work; Jx.Z. and S.T.: Analyzed the data; and Jx.Z., H.H., X.C., Jl.Z., Z.T.,Y.P., R.A.D., X.H., and J.F.: Wrote the manuscript. All contributing authors read and agreed to the final manuscript.
Correspondence to Jingyun Fang.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Data 1
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
Zhu, J., Hu, H., Tao, S. et al. Carbon stocks and changes of dead organic matter in China's forests. Nat Commun 8, 151 (2017). https://doi.org/10.1038/s41467-017-00207-1
The three-dimension zonal pattern of soil organic carbon density in China's forests
Ji Yuan
, Yu Zhang
, Chengming You
, Rui Cao
, Bo Tan
, Han Li
, Yurui Jiang
& Wanqin Yang
CATENA (2021)
Forest management effects on downed dead wood at stand and landscape scales in a temperate forest of the central United States
Shunzhong Wang
, Benjamin O. Knapp
, Shannon Ehlers
, Bradley Graham
, Xiaoming Gao
& Shelby Timm
Forest Ecology and Management (2021)
Multi-criteria decision making for sustainability assessment of boxboard production: A life cycle perspective considering water consumption, energy consumption, GHG emissions, and internal costs
Yi Man
, Yulin Han
, Yue Liu
, Ruojue Lin
& Jingzheng Ren
Journal of Environmental Management (2020)
Increasing soil carbon stocks in eight permanent forest plots in China
, Chuankuan Wang
, Zhang Zhou
, Guoyi Zhou
, Xueyang Hu
, Lai Jiang
, Yide Li
, Guohua Liu
, Chengjun Ji
, Shuqing Zhao
, Peng Li
, Jiangling Zhu
, Zhiyao Tang
, Chengyang Zheng
, Richard A. Birdsey
, Yude Pan
& Jingyun Fang
Biogeosciences (2020)
Trends in global research in forest carbon sequestration: A bibliometric analysis
Li Huang
, Mi Zhou
, Jie Lv
& Ke Chen
Journal of Cleaner Production (2020)
Forests in the Anthropocene
Editors' Highlights
Top Articles of 2019
Nature Communications ISSN 2041-1723 (online)