Skip to main content

Comparisons of allometric and climate-derived estimates of tree coarse root carbon stocks in forests of the United States



Refined estimation of carbon (C) stocks within forest ecosystems is a critical component of efforts to reduce greenhouse gas emissions and mitigate the effects of projected climate change through forest C management. Specifically, belowground C stocks are currently estimated in the United States’ national greenhouse gas inventory (US NGHGI) using nationally consistent species- and diameter-specific equations applied to individual trees. Recent scientific evidence has pointed to the importance of climate as a driver of belowground C stocks. This study estimates belowground C using current methods applied in the US NGHGI and describes a new approach for merging both allometric models with climate-derived predictions of belowground C stocks.


Climate-adjusted predictions were variable depending on the region and forest type of interest, but represented an increase of 368.87 Tg of belowground C across the US, or a 6.4 % increase when compared to currently-implemented NGHGI estimates. Random forests regressions indicated that aboveground biomass, stand age, and stand origin (i.e., planted versus artificial regeneration) were useful predictors of belowground C stocks. Decreases in belowground C stocks were modeled after projecting mean annual temperatures at various locations throughout the US up to year 2090.


By combining allometric equations with trends in temperature, we conclude that climate variables can be used to adjust the US NGHGI estimates of belowground C stocks. Such strategies can be used to determine the effects of future global change scenarios within a C accounting framework.


The management of forest ecosystems and their associated carbon (C) stocks has become an important global strategy for reducing greenhouse gas (GHG) emissions and possibly mitigating future effects of climate change [13]. Societal demands and trends in land use, in combination with future global change scenarios, may reduce the amount of C stored in forests and associated wood products [4]. As a result, there are substantial knowledge gaps regarding the C implications of various forest management activities, which may arise from the complex pathways of C emissions and sequestration in forest ecosystems [2, 3]. In addition, the logistical and methodological constraints associated with estimating C in certain forest ecosystem components across large areas, namely belowground pools, has hampered the development of accurate estimates, creating a need for refined modeling approaches to quantify belowground C stocks.

From an ecological perspective, the use of plant characteristics such as root to shoot ratio and root mass fraction have aided investigators in understanding belowground C stocks associated with coarse roots. For example, Poorter et al. [5] presented a mean root mass fraction of 0.21 for temperate and boreal forests and Smyth et al. [6] calculated a root to shoot ratio of 0.256 across Canada. These proportions lend insight into partitioning effects and belowground C stores. Various studies have examined the degree to which these ratios are altered under various stand and environmental conditions. Litton et al. [7] found that partitioning to belowground components increased with stand density in lodgepole pine forests (Pinus contorta Dougl. ex Loud. var. latifolia Engelm.), which may in turn be related to tree size. Root mass fractions have been found to decrease with total plant biomass consistently among angiosperm and gymnosperm forests [5, 8]. In tropical systems, 15–20 year-old plantations allocated more C belowground when compared to mature broad-leaved forests [9], highlighting the importance of accounting for management scenarios in assessments of belowground C stores. However, application of these findings to forest C accounting activities has been limited as few studies measure all components of the C budget (e.g., biomass, flux, and partitioning; [7]).

The monitoring of belowground C has incorporated a number of these ecological insights using a variety of approaches at different scales. Allometric equations designed at the individual-tree level are common for determining belowground C [e.g., 10, 11]. Throughout the United States (US), belowground biomass of coarse roots is commonly estimated using the equations of Jenkins et al. [12] as a ratio of total aboveground biomass and tree diameter at breast height (DBH; [13]). Although allometric equations tend to account for a large portion of the apparent variability associated with belowground biomass (e.g., R 2 values range from 0.77 to 0.96; Litton et al. [11]), there are a few drawbacks to this approach. First, these equations tend to rely on DBH and are not explicitly constructed to estimate belowground C [14]. Secondly, allometric equations have not historically incorporated climate information that integrates differences in ecosystem productivity and allows for evaluations of future climate change scenarios on global C cycles. Highlighting this concern, Reich et al. [8] recently compiled a global dataset and concluded that forest biomass found in roots was inversely related to mean annual temperature, suggesting that climate may act as a driver of belowground C allocation.

Globally, there has been an increased interest in recent years for refining forest carbon estimation to understand greenhouse gas emissions in support of the United Nations Framework Convention on Climate Change [15, 16]. Forest C stocks in the US are estimated using data collected by the US Forest Service, Forest Inventory and Analysis (FIA) program. In the current national greenhouse gas inventory (NGHGI; [17]), belowground stocks are estimated in two stages by first quantifying total aboveground biomass using allometric equations then estimating a ratio of coarse root to total aboveground biomass [12, 13]. As observations of belowground tree biomass and C are often limited [14], relying on allometric equations has been necessary to obtain estimates from strategic-scale forest inventories such as FIA’s. At the same time, the lack of empirical information across a diverse array of tree species in temperate forests such as those found throughout North America encourages researchers to test alternative approaches for quantifying belowground biomass and C. Exploring belowground C modeling approaches that incorporate climatic attributes may both adjust our estimates of coarse root C stocks at national scales (i.e., application in the US NGHGI) while enhancing evaluations of future climate change scenarios on forest C cycles.

The overall objective of this research is to adjust belowground C estimation procedures for reporting in the US NGHGI. Specific objectives are to (1) estimate belowground C stocks by employing individual tree- and stand-level methodologies, (2) adjust estimates of belowground C stocks by combining allometric and climate-derived approaches using current and projected climate attributes, and (3) compare alternative estimation approaches for belowground C stocks for future application in the US NGHGI.


Estimates of belowground carbon (BGC) from approaches currently employed in the US NGHGI suggest that C stocks are dependent on geographic region and forest type. Mean values of belowground carbon in the US greenhouse gas inventory (BGCNGHGI) were small in short-statured, open forests such as pinyon-juniper and woodland hardwood types (typically less than 2 Mg ha−1). Mean BGCNGHGI was largest in hemlock-Sitka spruce forests in the Pacific Northwest [40.76 ± 0.96 Mg ha−1 (mean ± SE)] and redwood forests in the Pacific Southwest (59.27 ± 7.06 Mg ha−1). For climate-derived estimates of belowground C, belowground carbon from climate-derived models (BGCClim) stock estimates were slightly smaller in magnitude compared to BGCNGHGI estimates [e.g., hemlock-Sitka spruce (33.82 ± 0.80 Mg ha−1) and redwood forests (45.64 ± 5.44 Mg ha−1)] and generally showed decreasing C at lower latitudes (Fig. 1). On average, BGCClim estimates were 0.60 Mg ha−1 greater than current BGCNGHGI models when considering all forest types (Additional file 1: Table S1).

Fig. 1
figure 1

Distribution of live-tree belowground C estimates from the model of Reich et al. [8] (BGCClim; Mg ha−1)

The adjustment factors used to align climate-derived predictions of belowground C with the US NGHGI approach ranged from 0.77 to 1.60 with little variability within a region of interest (Fig. 2; Additional file 1: Table S1). Equivalence tests used to contrast the two approaches for estimating belowground C with a null hypothesis of dissimilarity and a threshold of ±25 % were rejected for 20 out of the 78 forest types examined. Equivalence tests were generally rejected for forest types that displayed relatively low and high BGC stocks, e.g., pinyon-juniper and woodland hardwoods (low C stocks) and redwood forest types (high C stocks). Mean differences were generally largest across the Pacific Northwest (Westside), indicating climate-adjusted predictions estimated less belowground C compared to allometric-derived estimates in this region. Generally, negative mean differences were observed across most forest types, indicating that climate-adjusted models predict greater belowground C stocks (e.g., mean percent difference was −5.2 % greater across all forest types; Table 1). Compared to current NGHGI models, model differences showed greater belowground C stocks occurring in the Appalachian Mountain region and areas where northern hardwood forests are common, e.g., in the upper Midwest and northeastern US states. Conversely, areas of smaller belowground C stocks were identified across the Pacific Northwest and Southeast US (Fig. 3). This was further reflected when population estimates were scaled to the state level. The states of Oregon and Washington were predicted to display the largest negative mean difference in belowground C stocks (−10.6 and −10.7 %, respectively). Conversely, the largest mean positive difference in belowground C stocks was in the states of Kentucky, Tennessee, and Oklahoma (28.0, 26.7, and 22.6 %, respectively). This represents a total estimated increase of 368.87 Tg of belowground C across the US, or a 6.4 % increase when compared to currently implemented NGHGI models (Table 2).

Fig. 2
figure 2

Adjustment factors used to amend live-tree belowground C estimates from the current US national greenhouse gas inventory (BGCNGHGI; Smith et al. [13]) with climate-derived predictions (BGCClim; Reich et al. [8]) by US region and forest type. Adjustment factors less than and greater than one indicate less and more belowground C when climate-derived predictions are used. See Additional file 1: Table S1 for mean belowground C values and adjustment factors

Table 1 Equivalence test results (critical threshold of ±25 %) comparing live-tree belowground C estimates from the current US national greenhouse gas inventory (BGCNGHGI) [13] to adjusted estimates (BGCClimAdj)
Fig. 3
figure 3

Distribution of differences between live-tree belowground C estimates and adjusted estimates (BGCNGHGI − BGCClimAdj; Mg ha−1), with red colors indicating higher estimated C and purple colors less C

Table 2 Estimates of belowground carbon stocks (Tg) and associated sampling errors (SE; %) for current US national greenhouse gas inventory (BGCNGHGI) [13] and adjusted estimates (BGCClimAdj) by state

Results from the random forests (RF) regressions indicated that aboveground biomass, stand age, and stand origin (i.e., planted versus natural regeneration) were useful predictors of climate-adjusted models of belowground carbon (BGCClimAdj), as measured by their importance scores (Table 3). These variables accounted for 87 % of the total variation in belowground C stocks. Without employing aboveground biomass, 47 % of the total variability was accounted for, indicating that surrogates of climate (e.g., latitude and longitude) and knowledge of stand structure and management history (e.g., stand age and origin) may aid in understanding belowground C stocks.

Table 3 Summary of random forests model output and their importance scores (% IncMSE) for predicting refined estimates belowground C using basic stand structure, management, and physiographic variables

The largest differences in projected live-tree belowground C stocks under future mean annual temperature (MAT) changes were positive, indicating decreases in belowground C stocks up to year 2090 (Fig. 4). Differences in projected belowground C stocks were similar across regions, with belowground stocks displaying approximately 0.1 Mg ha−1 less C in 2030 than assuming current normal climates. In 2090, differences comparing current versus future climates were highest in the US Northeast, representing 0.50 ± 0.13 Mg ha−1 less C (assuming an RCP8.5 scenario), or a reduction of approximately 3.4 % in belowground C stocks in the region.

Fig. 4
figure 4

Regional differences in projected live-tree belowground C stocks compared to climate-sensitive estimates of belowground C derived from 30-year normal conditions (1960–1990). Error bars denote one standard deviation


Live tree belowground C estimated using allometric equations developed for individual trees and climate-specific predictions made at the forest stand level show markedly different patterns across various geographic regions and forest types in the US. Adjustments to the US’ approach to estimating live tree belowground C resulted in a C density estimate of 12 Mg ha−1 across the US and a 6.4 % increase in this forest C national stock when compared to current estimation strategies employed in the US NGHGI.

The largest differences between current and adjusted estimates of live tree belowground C were observed in the states of Oregon, Washington, and California, indicating climate-adjusted predictions estimated less belowground C compared to allometric-derived estimates in this region. Comparatively, this region contains the largest belowground C stocks in the US, quantified using model imputation strategies [18]. Similar differences but of a lesser magnitude were observed across the Coastal Plain region in the southeastern US. Such differences in these regions (i.e., the Pacific Northwest and Southeast) could point to the reliance on the allometric equations driven by aboveground metrics (i.e., tree diameter at breast height) in these carbon-dense stands. This may especially be the case for common tree species in the Pacific Northwest where the development of such allometric relationships is complicated by their relatively large above and belowground dimensions. In addition, intensive forest management regimes relying on artificial regeneration are common throughout these regions such that stand origin emerged as a useful predictor of live tree belowground C. Such a finding is analogous to recent work indicating greater partitioning to belowground components in managed plantations compared to natural-origin mature forests [9]. These patterns of partitioning may be reflective of the increased nutrient demands required to sustain the elevated levels of aboveground net primary production found in plantation systems [19], particularly on sites with lower nutrient capital such as in those found across the southeastern US. The importance of variables describing forest management regime and other forest conditions such as stand age in our random forests model output suggests that accounting for these aspects of forest origin and development is useful for understanding broad scale patterns in live tree belowground C.

The utilization of nationally consistent allometric equations [e.g., 12, 20] provides a general assessment of current C stocks. However, the appeal of employing BGCClimAdj estimates in the US NGHGI centers on utilizing standard forest inventory data collected in national forest inventories along with climate information. It is important to note that the model of Reich et al. [8] used here may be considered an empirical model. Such an empirical model may work well for data-rich countries with established national forest inventories. Other strategies, including process models, could potentially be employed to examine the carbon–climate relationship in forests. Incorporation of climate data enables evaluation of future climate change scenarios [e.g., 21] and their impacts on C monitoring efforts (e.g., United Nations Framework Convention on Climate Change submissions). As there are few studies that allow empirical observation of “true” coarse root C given the time and effort involved in destructively sampling belowground components, model validation procedures are limited for assessments at national scales. As an alternative, model-based approaches can be specified to be dynamic by incorporating future forest conditions and global change scenarios to determine their implications on C stocks and sequestration patterns. Such an approach using current climate conditions could provide general estimates of C stocks with associated uncertainty bounds for the temperate forests which occupy the US.

The majority of forest types displayed negative mean differences between current NGHGI and climate-adjusted models, indicating greater live tree belowground C stocks when using the adjusted models. The larger stocks in climate-adjusted models is partially a reflection of the ability of this framework to account for temperature-related shifts in patterns of belowground allocation within a species; a relationship held constant in current NGHGI models. In particular, a key component of the climate-adjusted model is an increased level of belowground allocation with decreasing temperatures given increasing nutrient limitation with colder temperatures [8]. These adjustments are reflected in the greater climate-adjusted estimates for northern latitude and high elevation forests (Fig. 2). Given the importance of climate variables in driving other belowground processes such as coarse root decomposition [22], the integration of these variables in models for describing other belowground stocks (e.g., dead roots) may not only aid in understanding C stock differences across large geographic scales but may also inform our understanding of other belowground processes.

The forecast of decreases in belowground C stocks up to year 2090 across all US regions and forest types is a function of projected increases in MAT and its role in the distribution of coarse root C stocks. National patterns in C stocks observed with US data indicated similar trends to the models presented in Reich et al. [8] showing that the proportion of total biomass to roots is greater in increasingly cold climates. Such assessments conducted in this analysis include a dynamic climate but assume constant aboveground biomass stocks, and a stand origin and phylogeny (e.g., conifer- or hardwood-dominated forest type). For example, the economic incentives of increasing C stocks to meet increased demand for wood for bioenergy [e.g., 23] has the potential to alter the proportion of planted compared to natural-origin stands at a national level. Although aridity was not found to influence the global distribution of coarse root biomass [8], what role might precipitation and/or its interaction with temperature and other climate variables play in determining future C sequestration patterns and stocks? Through designing national-scale models that incorporate climate parameters, the ability to quantify C stocks using relationships observed between climate and C dynamics is possible but is otherwise impracticable using allometric equations alone.

A particular concern from a physiological perspective relies in assessing carbon allocation tradeoffs in concert with changing climates. Changing allocation to belowground components can alter biomass accumulation and nutrient uptake [9], but without an assessment of additional carbon stocks from other pools (e.g., foliage and soil components), the role that future climate may play in determining overall stocks may be somewhat limited. Our approach in using current belowground C stocks from temperate forests across the US, with a range of current and projected climates, management histories, and stand structures, serves as a preliminary investigation of the role that climate may play on coarse root systems (e.g., Fig. 4). Carbon allocation to soil microbial biomass may be limited at northern sites with cooler climates [24] which would seemingly influence nutrient uptake and allocation within tree components. Process models [e.g., 25] may be well suited to examine such carbon allocation-climate tradeoffs, presenting a different approach compared to the empirical models examined here. In the interim, additional biomass data collected on all components from a range of species across sites with different climate regimes and management histories will aid in improving our understanding of carbon allocation patterns related to climate [26].

The finding that surrogates of climate (e.g., latitude and longitude) and knowledge of forest structure and management history (e.g., stand age and origin) were useful in predicting belowground C is encouraging when considering approaches to constructing a NGHGI. For example, any new strategy for estimating forest C pools in the US NGHGI requires “back casting” estimates to the 1990 reporting year [17]. Obtaining climate data from past years at large geographic scales may present less of a barrier to refinement of NGHGIs than other efforts (e.g., in situ field inventories in the 1990’s) to reduce the uncertainty associated with estimates of forest C and harvested wood products. Other nations could similarly adjust their C accounting practices using insight from this analysis to determine belowground stocks. Regardless of whether or not tree-level information from a national forest inventory is available, incorporating indicators of forest management history (e.g., proportion of planted versus natural origin stands) and variability in climate within a country (e.g., general trends in temperature), estimates of belowground C stocks could be developed that theoretically represent site and regional differences. In an era where destructive samples of belowground tree components is incomplete across the world’s forests, estimation strategies that merge the attributes of both allometric equations and C–climate relationships may refine forest C stocks estimations especially given emerging science that supports the relationships between climate, forest biomass, and C pools [1, 3, 8, 27].


Numerous findings emerged from our investigation by incorporating climate variables into the estimation of belowground C stocks. First, climate variables can be used to adjust the US NGHGI estimates of belowground C stocks. Specifically, adjustment factors were specified to amend current coarse root C stocks estimated from allometric equations by incorporating mean annual temperature at various locations across the US. Second, for the US NGHGI, incorporating mean annual temperature increased national belowground C stocks by 6.4 %. In contrast, coarse root C stocks were projected to decrease through 2090, primarily due to lower partitioning to belowground components under warmer conditions. Third, whether or not a forest was planted or from natural origin, and its stand age were influential variables in determining belowground C stocks. Future work that integrates both climate and stand origin will increase our ability to predict belowground C stocks across regions containing a mixture of management and climate regimes. Finally, as a means of refining NGHGIs, climate-adjusted models depicting belowground C stocks should be adopted to incorporate the impacts of future global change and management scenarios on C sequestration patterns and stocks.


Study area

Forests across the US range are characterized by four major ecoclimatic zones, including polar, temperate humid, arid, and tropical humid types [28, 29]. The study area investigated here included forestlands across the contiguous US, spanning approximately 24° latitude (LAT) and 58° longitude (LONG) in addition to coastal Alaska (mean coordinates 57.87°N, 138.60°W). Mean annual temperature (MAT) ranged from −3.0 to 24.9 °C and precipitation (MAP) from 18 to 420 cm [30, 31]. Nine broad geographic regions were identified across the study area, ultimately containing 78 unique forest types [13] (Additional file 1: Table S1).

Forest Inventory and Analysis data

The FIA program within the US Forest Service monitors forests by establishing permanent sample plots across the US in three phases [32]. During the inventory’s first phase, sample plot locations are established at an intensity of approximately 1 plot per 2400 ha. If the plot lies partially or wholly within a forested area, field personnel visit the site and establish a phase two (P2) inventory plot. Standard P2 inventory plots consist of four 7.32-m fixed radius subplots for a total plot area of approximately 0.07 ha where standing tree and site attributes are measured. Live trees with a DBH of at least 12.7 cm are measured on these subplots. Within each subplot a 2.07-m microplot is established where saplings with a DBH between 2.5 and 12.7 cm are measured.

All data were obtained from the publically-available FIA database (FIADB; [33];; download date 14 May 2014). If an FIA plot was remeasured at any point, only the most recent measurement was used in the analysis. Using the individual tree measurements, aboveground live-tree biomass (BIOAG; Mg ha−1) was estimated by summing the bole, stump, top (excluding foliage), saplings, and woodland tree species (primarily those from dryland forests) components of each plot. Additional condition-level information including stand age (STDAGE; years) and a binary variable depicting stand origin (NAT; 1 = natural, 0 = artificial regeneration) were subsequently analyzed for each plot. As a bioclimatic measure, The Hopkins index ([34]; HI) standardizes the onset of spring for a given region and was computed for each FIA plot sampled relative to the mean LAT (40.35), LONG (−95.84), and elevation (ELEV; 2283 feet):

$$\begin{aligned} {\rm HI} = \left( {\frac{{\rm ELEV} - 2283}{100}} \right) + 4\,\left( {{\rm LAT} - 40.35} \right) + 1.25\,\left( { - 95.84 - {\rm LONG}} \right) \hfill \\ \ \hfill \\ \end{aligned}$$

The HI variable assumes that spring is delayed by one day for each 100-foot rise in ELEV, four days for each 1° increase in LAT, and 1.25 days for each 1° increase in LONG (e.g., more westward; [34]). In total, 70,126 FIA plots were analyzed for their belowground C.

Belowground C in the US greenhouse gas inventory (BGCNGHGI)

The Intergovernmental Panel on Climate Change’s (IPCC) Good Practice Guidance considers forest C stocks associated with live aboveground, live belowground, dead wood, litter, and soil organic pools [35]. Hence, the focus of this analysis is on adjusting estimates of live belowground C, defined as all coarse living roots greater than 2 mm diameter [13]. Estimates of aboveground live C in the NGHGI are calculated using the component ratio method [36], but do not ultimately influence the analysis for belowground C here.

Belowground C for live trees is estimated in two stages using allometric equations. First, total aboveground biomass is estimated as a function of tree DBH [12]:

$${\rm Biomass} = \exp \left( {\alpha_{1,i} + \alpha_{2,i} \,\text{ln\,DBH} } \right)$$

where α1,i and α2,i are parameters for one of ten species groups [four hardwood groups (aspen/alder/cottonwood/willow, soft maple/birch, mixed hardwood, and hard maple/oak/hickory/beech), five conifer groups (cedar/larch, Douglas-fir, true fir/hemlock, pine, and spruce), and one woodland species group (juniper, oak, mesquite)]. Second, belowground root biomass is estimated as a ratio (BGRATIO) of root to total aboveground biomass [12]:

$${\rm BG}_{\rm RATIO} = \left\{ \begin{aligned} \exp \left( { - 1.6911 + \frac{0.8160}{\rm DBH}} \right),\,\,{\rm for}\,\,{\rm hardwood}\,\,{\rm species} \hfill \\ \exp \left( { - 1.5619 + \frac{0.6614}{\rm DBH}} \right),\,\,{\rm for}\,\,{\rm conifer}\,\,{\rm species} \hfill \\ \end{aligned} \right.$$

Hence, parameters indicate BGRATIO will decrease for larger DBH trees and that for a fixed DBH, BGRATIO will be larger for conifer compared to hardwood species. Belowground biomass was estimated by multiplying the values obtained from Eqs. 2 and 3, then converted to C by multiplying by 0.5, assuming 50 % of biomass is C [35]. Estimates of belowground C were scaled to the plot level and are hereby abbreviated as BGCNGHGI.

Belowground C from climate-derived models (BGCClim)

Recent investigations of the global distribution of biomass within forests have provided insight for comparing size- and species-specific predictions (i.e., allometric equations and proportional ratios) with climate-sensitive estimations of belowground C [e.g., 8]. Using the relationship observed between mean annual temperature and root mass fraction, the model used to determine a climate-derived estimate of belowground C was parameterized with global data compiled from various sources including Usoltsev [37], Luo et al. [38], Cannell [39], and from over 1000 additional forest stands including the US [8]. A total of 3043 of these stands contained measurements of belowground biomass. The motivation for the development of this model was to assess the distribution of biomass in roots along a temperature and precipitation spectrum [8]. We estimated a climate-sensitive prediction of belowground biomass (BGBClim) using the model of Reich et al. [8] (Table 4). In addition to MAT and NAT, a dummy variable indicating whether or not the FIA plot was primarily dominated by hardwoods or conifers (HDWD) and stem biomass (BIOSTEM; Mg ha−1) of live trees were used to estimate BGBClim. Thirty-year (1961–1990) climate data (i.e., MAT) were obtained by specifying LAT, LONG, and ELEV of each FIA plot location to a spline surface model developed from climate station data across forests of North America [30, 31]. We assigned the HDWD variable using the FIA forest type code [40] by separating conifer-dominated forest type codes (i.e., FORTYPCD ≤ 409) with hardwood-dominated codes (FORTYPCD ≥ 500). By incorporating measures of aboveground biomass, these allometric relationships allow one to capture the variability observed across a range of stand structures and ages. Continuous independent variables were centered prior to applying the Reich et al. [8] model to the FIA plots. When standardized for a given MAT and aboveground stem biomass, the model of Reich et al. [8] indicates conifer forests tended to have a smaller root mass fraction than hardwood forests. Values for BGBClim were converted to BGCClim by multiplying by 0.5 [35].

Table 4 Model parameters and model form Reich et al. [8] used in this analysis for estimating climate-adjusted belowground biomass (BGBClim)

Climate-adjusted models of belowground C (BGCClimAdj)

Current models of belowground C in the NGHGI (i.e., BGCNGHGI) could likely be adjusted by incorporating climate-derived estimates of belowground C (i.e., BGCClim). From a NGHGI reporting perspective, estimates of BGC would need to be made on individual plots, then “back cast” to contemporary estimates using the 1990 baseline reporting year [17]. Adjustment factors were estimated to align allometric- and climate-derived estimates:

$${\text{AdjFactor}} = \frac{{\text{BGC}}_{{\rm Clim}}}{{\text{BGC}}_{{\rm NGHGI}}}$$

where AdjFactor is the ratio of climate- to allometric-derived belowground C for a specific forest type found in a given geographic region. New climate-adjusted estimates of belowground C (BGCClimAdj) are then:

$${\text{BGC}}_{\text{ClimAdj}} = {\text{BGC}}_{\text{NGHGI}} \times {\text{AdjFactor}}$$

where BGCClimAdj was computed for each FIA plot record.

Analyzing belowground C model differences

We conducted equivalence tests comparing BGCNGHGI and BGCClimAdj model predictions using two one-sided tests [41]. Equivalence tests are commonly applied in the forest science literature and are advantageous in that they can be used in model validation by assuming a null hypothesis of dissimilarity [42]. Equivalence tests are unlike statistical goodness-of-fit approaches and instead examine dissimilarity. Dissimilarity in the equivalence test was specified using a threshold of ±25 %. This threshold allows for a moderate amount of disagreement between the various model predictions, with non-equivalence suggesting biological disparities in C stocks. Differences between BGCNGHGI and BGCClimAdj models were mapped across the US to examine geographic trends in estimates of belowground C when using each approach. We computed US state-level population estimates of belowground C (Tg) using BGCNGHGI and BGCClimAdj models and compared mean percent differences for the two estimation strategies.

Nonparametric random forests (RF; [43]) were implemented in R [44] to identify variables that were effective in describing BGCClimAdj. Recognizing that not all users may implement climate data in determining forest C stocks, seven variables from the FIADB (BIOAG, STDAGE, NAT, LAT, LONG, ELEV, and HI) were chosen for incorporation into the RF. This method involved building a set of regression trees based on bootstrapped samples of the belowground C data. We similarly fit a RF model without BIOAG to examine how belowground C can be predicted without knowledge of aboveground biomass stocks.

Belowground C in future scenarios

Current CMIP5 models [45] as described in the fifth assessment report (AR5) of the IPCC [46] were obtained using three scenarios (RCP 4.5, RCP 6.0, RCP 8.5; [30]). An ensemble of 17 AR5 model predictions was used for each RCP scenario. Provided that differences in belowground C could exist in future global change scenarios (i.e., changes in MAT at various locations in the US), climate data were obtained for the 30-year normal (1961–1990) and years 2030, 2060, and 2090. To gain insight into temperature-related patterns and their influence on belowground C, climate-sensitive estimates were made using these new MAT values in future years (while holding fixed the variables HDWD, NAT, and BIOSTEM) to assess the variation in belowground C assuming future scenarios. While BIOSTEM will fluctuate in response to trends in forest growth and mortality, this simulation was specifically designed to understand the influence of MAT in future belowground C stocks.





United States’ National Greenhouse Gas Inventory


greenhouse gas


Forest Inventory and Analysis






mean annual temperature


mean annual precipitation


aboveground live-tree biomass


stand age


natural regeneration


random forests


Hopkins index




belowground carbon in the US greenhouse gas inventory


Intergovernmental Panel on Climate Change


tree diameter at breast height


belowground root biomass ratio

BGCClim :

belowground carbon from climate-derived models

BGBClim :

belowground biomass from climate-derived models


hardwood forest type


stem biomass


forest type code

BGCClimAdj :

climate-adjusted models of belowground carbon


the ratio of climate- to allometric-derived belowground carbon


  1. Ryan MG, Harmon ME, Birdsey RA, Giardina CP, Heath LS, Houghton RA, et al. A synthesis of the science on forests and carbon for U.S. forests. Ecol Soc Am. 2010;13:1–16.

    Google Scholar 

  2. Malmsheimer RW, Bowyer JL, Fried JS, Gee E, Izlar RL, Miner RA, et al. Managing forests because carbon matters: integrating energy, products, and land management policy. J For. 2011;109:S7–50.

    Google Scholar 

  3. McKinley DC, Ryan MG, Birdsey RA, Giardina CP, Harmon ME, Heath LS, et al. A synthesis of current knowledge on forests and carbon storage in the United States. Ecol Appl. 2011;21:1902–24.

    Article  Google Scholar 

  4. Joyce LA, Running SW, Breshears DD, Dale VH, Malmsheimer RW, Sampson RN, et al. Ch. 7: Forests. Climate change impacts in the United States: the third national climate assessment. In: Melillo JM, Richmond TC, Yohe GW, editors. Washington, D.C: U.S. Global Change Research Program; 2014. p. 175–94. doi:10.7930/J0Z60KZC.

  5. Poorter H, Niklas KJ, Reich PB, Oleksyn J, Poot P, Mommer L. Biomass allocation to leaves, stems and roots: meta-analyses of interspecific variation and environmental control. New Phytol. 2012;193(1):30–50. doi:10.1111/j.1469-8137.2011.03952.x.

    Article  CAS  Google Scholar 

  6. Smyth CE, Kurz WA, Neilson ET, Stinson G. National-scale estimates of forest root biomass carbon stocks and associated carbon fluxes in Canada. Glob Biogeochem Cycles. 2013;27(4):2012GB004536. doi:10.1002/2012gb004536.

    Article  Google Scholar 

  7. Litton CM, Raich JW, Ryan MG. Carbon allocation in forest ecosystems. Glob Chang Biol. 2007;13(10):2089–109. doi:10.1111/j.1365-2486.2007.01420.x.

    Article  Google Scholar 

  8. Reich PB, Luo Y, Bradford JB, Poorter H, Perry CH, Oleksyn J. Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc Natl Acad Sci USA. 2014;111(38):13721–6.

    Article  CAS  Google Scholar 

  9. Raich JW, Clark DA, Schwendenmann L, Wood TE. Aboveground tree growth varies with belowground carbon allocation in a tropical rainforest environment. PLoS One. 2014;9(6):e100275. doi:10.1371/journal.pone.0100275.

    Article  Google Scholar 

  10. Bond-Lamberty B, Wang C, Gower ST. Aboveground and belowground biomass and sapwood area allometric equations for six boreal tree species of northern Manitoba. Can J For Res. 2002;32(8):1441–50. doi:10.1139/x02-063.

    Article  Google Scholar 

  11. Litton CM, Ryan MG, Tinker DB, Knight DH. Belowground and aboveground biomass in young postfire lodgepole pine forests of contrasting tree density. Can J For Res. 2003;33(2):351–63. doi:10.1139/X02-181.

    Article  Google Scholar 

  12. Jenkins JC, Chojnacky DC, Heath LS, Birdsey RA. National-scale biomass estimators for United States tree species. For Sci. 2003;49(1):12–35.

    Google Scholar 

  13. Smith JE, Heath LS, Hoover CM. Carbon factors and models for forest carbon estimates for the 2005–2011 National Greenhouse Gas Inventories of the United States. For Ecol Manag. 2013;307:7–19.

    Article  Google Scholar 

  14. Zianis D, Muukkonen P, Mäkipää R, Mencuccini M. Biomass and stem volume equations for tree species in Europe. Silva Fenn Monog. 2005;4:63.

    Google Scholar 

  15. Beets PN, Brandon AM, Goulding CJ, Kimberley MO, Paul TSH, Searles N. The inventory of carbon stock in New Zealand’s post-1989 planted forest for reporting under the Kyoto protocol. For Ecol Manag. 2011;262(6):1119–30. doi:10.1016/j.foreco.2011.06.012.

    Article  Google Scholar 

  16. Green C, Tobin B, O’Shea M, Farrell EP, Byrne KA. Above- and belowground biomass measurements in an unthinned stand of Sitka spruce (Picea sitchensis (Bong) Carr.). Eur J For Res. 2006;126:179–88.

    Article  Google Scholar 

  17. US EPA. Inventory of U.S. Greenhouse Gas Emissions and Sinks: 1990–2012. Chapter 7. Land Use, Land-Use Change, and Forestry. Annex 3.12. Methodology for Estimating Net Carbon Stock Changes in Forest Land Remaining Forest Lands. #430-R-14-003. Washington, DC: U.S. Environmental Protection Agency; 2014. Accessed 3 July 2014.

  18. Wilson BT, Woodall CW, Griffith D. Imputing forest carbon stock estimates from inventory plots to a nationally continuous coverage. Carbon Balance Manag. 2013;8:1.

    Article  Google Scholar 

  19. McMurtrie RE, Dewar RC. New insights into carbon allocation by trees from the hypothesis that annual wood production is maximized. New Phytol. 2013;199:981–90.

    Article  CAS  Google Scholar 

  20. Ter-Mikaelian MT, Korzukhin MD. Biomass models for sixty-five North American tree species. For Ecol Manag. 1997;97(1):1–24.

    Article  Google Scholar 

  21. Melillo JM, Richmond TC, Yohe GW. Climate change impacts in the United States: the third national climate assessment. U.S. Global Change Research Program. 2014. p. 841.

  22. Garrett LG, Kimberley MO, Oliver GR, Pearce SH, Beets PN. Decomposition of coarse woody roots and branches in managed Pinus radiata plantations in New Zealand—a time series approach. For Ecol Manag. 2012;269:116–23. doi:10.1016/j.foreco.2011.12.030.

    Article  Google Scholar 

  23. Sedjo R, Tian X. Does wood bioenergy increase carbon stocks in forests? J For. 2012;110(6):304–11.

    Google Scholar 

  24. Gorissen A, Tietema A, Joosten NN, Estiarte M, Peñuelas J, Sowerby A, et al. Climate change affects carbon allocation to the soil in shrublands. Ecosystems. 2004;7(6):650–61. doi:10.1007/s10021-004-0218-4.

    Article  CAS  Google Scholar 

  25. Running SW, Gower ST. FOREST-BGC, A general model of forest ecosystem processes for regional applications. II. Dynamic carbon allocation and nitrogen budgets. Tree Physiol. 1991;9(1–2):147–60. doi:10.1093/treephys/9.1-2.147.

    Article  CAS  Google Scholar 

  26. Weiskittel AR, MacFarlane DW, Radtke PJ, Affleck DLR, Temesgen H, Woodall CW, et al. A call to improve methods for estimating tree biomass for regional and national assessments. J For. 2015;113(4):414–24. doi:10.5849/jof.14-091.

    Google Scholar 

  27. Nemani RR, Keeling CD, Hashimoto H, Jolly WM, Piper SC, Tucker CJ, et al. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science. 2003;300(5625):1560–3. doi:10.1126/science.1082750.

    Article  CAS  Google Scholar 

  28. Smith WB, Miles PD, Perry CH, Pugh SA. Forest resources of the United States, 2007. USDA For. Serv. Gen. Tech. Rep. WO-78. 2009. p. 336.

  29. Bailey RG. Ecosystem geography. New York: Springer; 1996. p. 216.

    Book  Google Scholar 

  30. US Forest Service. Research on forest climate change: potential effects of global warming on forests and plant climate relationships in western North America and Mexico. Rocky Mountain Research Station, Moscow Laboratory. 2014. Accessed 28 Mar 2014.

  31. Rehfeldt GE. A spline model of climate for the western United States. USDA For. Serv. Gen. Tech. Rep. RMRS-165. 2006.

  32. Bechtold WA, Patterson PL, editors. Forest Inventory and Analysis national sample design and estimation procedures. USDA For. Serv. Gen. Tech. Rep. SRS-GTR-80. 2005.

  33. US Forest Service. The Forest Inventory and Analysis database: database description and user guide version 6.0 for phase 2. Washington DC; 2014. p. 624.

  34. Hopkins AD. Bioclimatics, a science of life and climate relations. USDA Misc. Publ. 280. USDA, Washington, DC; 1938. p. 188.

  35. Intergovernmental Panel on Climate Change (IPCC). Good practice guidance for land use, land use change and forestry. Institute for Global Environmental Studies (IGES). 2003. Accessed 21 May 2014.

  36. Woodall CW, Heath LS, Domke GM, Nichols MC. Methods and equations for estimating aboveground volume, biomass, and carbon for trees in the U.S. forest inventory, 2010. Gen. Tech. Rep. NRS-88. Newtown Square, PA: U.S. Department of Agriculture, Forest Service, Northern Research Station (1 CD-ROM). 2011. p. 30.

  37. Usoltsev VA. Forest biomass of northern Eurasia. Russian Academy of Sciences, Ural Branch, Botanical Garden and Ministry of Education of Russian Federation, Ural State Forest Engineering University, Yekaterienburg) (in Russian). 2001.

  38. Luo Y, Wang X, Zhang X, Booth T, Lu F. Root:shoot ratios across China’s forests: forest type and climatic effects. For Ecol Manag. 2012;269:19–25.

    Article  Google Scholar 

  39. Cannell MGR. World forest biomass and primary production data. London: Academic Press; 1982. p. 391.

    Google Scholar 

  40. Woudenberg SW, Conkling BL, O’Connell BM, LaPoint EB, Turner JA, Waddell KL. The Forest Inventory and Analysis Database: database description and users manual version 4.0 for phase 2. US Department of Agriculture, Forest Service Gen. Tech. Rep. RMRS-245. 2010. p. 339.

  41. Wellek S. Testing statistical hypotheses of equivalence. London: Chapman and Hall; 2003.

    Google Scholar 

  42. Robinson AP, Froese RE. Model validation using equivalence tests. Ecol Model. 2004;176(3–4):349–58. doi:10.1016/j.ecolmodel.2004.01.013.

    Article  Google Scholar 

  43. Breimen L. Random forests. Mach Learn. 2001;45:5–32.

    Article  Google Scholar 

  44. Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2(3):18–22.

    Google Scholar 

  45. Taylor KE, Stouffer RJ, Meehl GA. An overview of CMIP5 and the experiment design B. Am Meteorol Soc. 2012;93:485–98.

    Article  Google Scholar 

  46. Intergovernmental Panel on Climate Change (IPCC). Climate change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM, editors. Cambridge: Cambridge University Press; 2013. p. 1535.

Download references

Authors’ contributions

MBR, GMD, and CWW designed the study and performed the research; MBR and GMD analyzed data and contributed new methods; all authors discussed the results of the study and contributed to writing the manuscript. All authors read and approved the final manuscript.


We thank Brian Clough, Christopher Oswalt, and John Stanovick for comments that improved this work and Georgii Alexandrov and two anonymous reviewers for their evaluation of this paper.

Compliance with ethical guidelines

Competing interests The authors declare that they have no competing interests.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Matthew B. Russell.

Additional file

Additional file 1: Table S1.

Summary of live-tree belowground C estimates (Mg ha−1) from the current US national greenhouse gas inventory (BGCNGHGI; Smith et al. 2013), climate-derived predictions (BGCClim; Reich et al. 2014), and adjustment factors (AdjFactor) used to amend BGCNGHGI estimates by taking into account climate-derived predictions.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Russell, M.B., Domke, G.M., Woodall, C.W. et al. Comparisons of allometric and climate-derived estimates of tree coarse root carbon stocks in forests of the United States. Carbon Balance Manage 10, 20 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: