Generic Additive Allometric Models and Biomass Allocation for Two Natural Oak Species in Northeastern China

: Current models for oak species could not accurately estimate biomass in northeastern China, since they are usually restricted to Mongolian oak ( Quercus mongolica Fisch. ex Ledeb.) on local sites, and speciﬁcally, no biomass models are available for Liaodong oak ( Quercus wutaishanica Mayr). The goal of this study was, therefore, to develop generic biomass models for both oak species on a large scale and evaluate the biomass allocation patterns within tree components. A total of 159 sample trees consisting of 120 Mongolian oak and 39 Liaodong oak were harvested and measured for wood (inside bark), bark, branch and foliage biomass. To account for the belowground biomass, 53 root systems were excavated following the aboveground harvest. The share of biomass allocated to different components was assessed by calculating the ratios. An aboveground additive system of biomass models and belowground equations were ﬁtted based on predictors considering diameter ( D ), tree height ( H ), crown width ( CW ) and crown length ( CL ). Model parameters were estimated by jointly ﬁtting the total and the components’ equations using the weighted nonlinear seemingly unrelated regression method. A leave-one-out cross-validation procedure was used to evaluate the predictive ability. The results revealed that stem biomass accounts for about two-thirds of the aboveground biomass. The ratio of wood biomass holds constant and that of branches increases with increasing D , H , CW and CL , while a reverse trend was found for bark and foliage. The root-to-shoot ratio nonlinearly decreased with D , ranging from 1.06 to 0.11. Tree diameter proved to be a good predictor, especially for root biomass. Tree height is more prominent than crown size for improving stem biomass models, yet it puts negative effects on crown biomass models with non-signiﬁcant coefﬁcients. Crown width could help improve the ﬁtting results of the branch and foliage biomass models. We conclude that the selected generic biomass models for Mongolian oak and Liaodong oak will vigorously promote the accuracy of biomass estimation.


Introduction
Mongolian oak (Quercus mongolica Fisch. ex Ledeb.) and Liaodong oak (Quercus wutaishanica Mayr) are very similar in both morphological characteristics and geographical distribution. Both species are valuable and occupy the secondary broadleaved deciduous forest regions widely distributed in northeastern China [1]. Furthermore, the two oak species have strong adaptability, with characteristics of resistance to drought and impoverishment, which play an important role in biodiversity protection and sustainable forestry development. To assess forest productivity, carbon stocks and dynamics, matter and energy flows and forest ecosystem services, accurate quantification of forest biomass depends on environmental conditions [7]. Furthermore, crown width (CW) or crown area, wood density and tree age have been integrated into biomass estimations to enhance the predictive abilities of models [5,22,27,28]. However, some variables, such as wood density, need additional work, extending the project time and increasing costs [29]. Furthermore, allometric models with multiple predictors require forest inventory datasets that contain the relevant factors.
The objective of this study was to develop generic biomass models for the two oak species in northeastern China. We focused on additive systems of biomass equations for the aboveground components, but we also constructed biomass models for roots. Both the aboveground and root biomass equations were validated through the leave-one-out cross-validation method and also compared to published models to explore the potential for application. Furthermore, we assessed the biomass allocation patterns and presented the trend that the RS ratios varied with diameter.

Study Area Description
This study was conducted in the Daxing'anling mountains (47 • Figure 1). The three mountains have a distinct continental monsoon climate. In the Daxing'anling mountains, the elevation varies between 425 m and 1760 m above sea level, the mean annual temperature is −2.8 • C and the annual precipitation is around 460 mm, which is mainly distributed from May to October. In the Xiaoxing'anling mountains, the elevation ranges from 500 m to 1000 m, the mean annual temperature is from −2.0 • C to 2.0 • C and the annual precipitation is from 550 mm to 670 mm, concentrated mostly in summer. In the Changbai mountains, the elevation ranges from 800 m to 1800 m, the mean annual temperature is from −7.0 • C to 3.0 • C and the annual precipitation is from 700 mm to 1400 mm, which mainly happens from June to October. The substrate soils are predominantly brown coniferous forest soil for the three regions.

Biomass Data
A total of 159 oak trees (120 and 39 Mongolian oak and Liaodong oak trees, respectively) were selected with a D between 1.5 cm and 33.0 cm. A destructive sampling method with direct field measurements was adopted to gather the biomass data, where the stem wood (inside bark), bark, branch and foliage components were measured separately. Because of the difficulty in root excavation, 53 trees with a D spanning 1.5-32.8 cm were selected for root biomass collection and used to establish allometric equations for root biomass estimation.
For each tree, the CW (twice the crown radius, i.e., twice the distance between the center of the bole and the outer edge of the crown) of the sampled trees was obtained from the average of two diameters measured along the North/South and East/West orientation. After felling with chain saws, D (1.30 m), H and crown length (CL, the distance from the first live branch to the treetop) were immediately measured and recorded. A general positive correlation between H and D was observed, as displayed in Figure 2, in which the vertical bars indicate the sample distribution across the different diameter classes.
For each felled tree, the live crown was evenly separated into three layers, and three sample branches that were representative in size and foliage amount for each layer were selected to separate, weighing the fresh weight of the branches and foliage. Subsamples of the branches and foliage were then randomly collected and weighed, and bagged for water content determination in the laboratory. Subsequently, the stem was equally divided into ten sections and each section was weighed directly in the field. At the position corresponding to a 0.1, 0.3 and 0.7 tree height, 3-to 5-cm-thick stem discs were cut, and after removing the bark, weighed to acquire the fresh weight of the stem wood. These discs and the removed bark were taken back to the laboratory for determining the water content. For root extraction, the soil was manually and carefully removed to obtain as intact roots as possible, especially fine roots. All extracted roots were weighed and recorded. The tree stump was dug out completely at last to measure the fresh weight. Subsamples of different root classes and stumps were randomly gathered and weighed, and taken back to the laboratory for determining the water content.

Biomass Data
A total of 159 oak trees (120 and 39 Mongolian oak and Liaodong oak trees, respectively) were selected with a D between 1.5 cm and 33.0 cm. A destructive sampling method with direct field measurements was adopted to gather the biomass data, where the stem wood (inside bark), bark, branch and foliage components were measured separately. Because of the difficulty in root excavation, 53 trees with a D spanning 1.5-32.8 cm were selected for root biomass collection and used to establish allometric equations for root biomass estimation.
For each tree, the CW (twice the crown radius, i.e., twice the distance between the center of the bole and the outer edge of the crown) of the sampled trees was obtained from the average of two diameters measured along the North/South and East/West orientation. After felling with chain saws, D (1.30 m), H and crown length (CL, the distance from the first live branch to the treetop) were immediately measured and recorded. A general positive correlation between H and D was observed, as displayed in Figure 2, in which the vertical bars indicate the sample distribution across the different diameter classes.

Biomass Data
A total of 159 oak trees (120 and 39 Mongolian oak and Liaodong oak trees, respectively) were selected with a D between 1.5 cm and 33.0 cm. A destructive sampling method with direct field measurements was adopted to gather the biomass data, where the stem wood (inside bark), bark, branch and foliage components were measured separately. Because of the difficulty in root excavation, 53 trees with a D spanning 1.5-32.8 cm were selected for root biomass collection and used to establish allometric equations for root biomass estimation.
For each tree, the CW (twice the crown radius, i.e., twice the distance between the center of the bole and the outer edge of the crown) of the sampled trees was obtained from the average of two diameters measured along the North/South and East/West orientation. After felling with chain saws, D (1.30 m), H and crown length (CL, the distance from the first live branch to the treetop) were immediately measured and recorded. A general positive correlation between H and D was observed, as displayed in Figure 2, in which the vertical bars indicate the sample distribution across the different diameter classes. For each felled tree, the live crown was evenly separated into three layers, and three sample branches that were representative in size and foliage amount for each layer were selected to separate, weighing the fresh weight of the branches and foliage. Subsamples All subsamples of the stem wood, bark, branches, foliage, roots and stumps were oven-dried at 105 • C for 72 h to constant weight to estimate the water content. The dry weights of the different components were then computed by subtracting the weight of the water from the fresh weights; i.e., the dry biomass of each component was calculated by multiplying the fresh weight by the respective ratio of the dry mass to fresh mass. For each sampled tree, the sum of the dry biomass of the stem wood, bark, branches and foliage yielded the aboveground biomass. Root (belowground) biomass was regarded as the sum of the root dry biomass and stump dry biomass (Table 1).

Data Analysis
Aboveground biomass allocated to wood, bark, branches and foliage was evaluated by calculating the biomass fractions for each component; i.e., the proportion of the component's biomass to the aboveground biomass. Bivariate relationships between biomass fractions and D, H, CW and CL were further graphically examined for each component. The RS ratio (i.e., root biomass/aboveground biomass) varying with the diameter class was graphically examined.

Allometric Models Specification
Allometry relationships between the component biomass and explanatory variables are generally expressed as a power function, such that where W represents the biomass (kg), X i (i = 1, 2, . . . , n) represent the easily measured explanatory variables that reflect the tree dimensions, β i (i = 1, 2, . . . , n) are model coefficients and ε is the error term, which is assumed to be normally distributed. As biomass data are significantly heteroscedastic, some approaches should be therefore taken to eliminate heteroscedasticity before parameter estimation [30]. In the present study, the weighted nonlinear regression was adopted, and the specific weighting function expressed as w = 1/ √ D Q , with a parameter Q, was derived from the residuals of the biomass models fitted separately by the ordinary least squares method [15,31].
The number of explanatory variables used in the allometric biomass models always ranges from one single predictor (D) to several or multiple predictors (e.g., H, CW and CL) [32]. Since the tree-dimension information available from the oak datasets in this study were D, H, CW and CL, we established the basic biomass models based on D alone and tested the model fitting improvement by adding the H, CW and CL variables as additional predictors. Thus, the following seven candidate models were fitted to the data: where W represents the biomass (kg), β i are coefficients and ε is the error term.
To ensure the additivity, tree aboveground biomass and its components were simultaneously fitted using the NSUR method based on the aforementioned seven nonlinear models (Equations (2)-(8)), with D, H, CW and CL as explanatory variables. Due to the difference in the sample size of the root and aboveground biomass, root biomass models were separately developed using these equations. Following the model structure specified by Parresol [15], the additive model system constrains the components' biomass to equal the aboveground biomass, and can be written as: where W m represents components' biomass (kg); m = wd, bk, br, fo and ag represent the wood, bark, branch, foliage and aboveground biomass, respectively; f m (X m , β m ) is a nonlinear equation of the biomass component m; and ε m is an error term for biomass component m.
In the additive model system, all components and aboveground biomass models were jointly fitted adopting the weighted NSUR [15,20]. Fitting was conducted using the PROC MODEL procedure of SAS 9.3 (SAS Institute Inc., Cary, NC, USA).

Model Evaluation and Selection
The aboveground additive biomass model systems and root biomass equations were fitted to the entire dataset. The model fitting was assessed by looking at two goodnessof-fit statistics: root mean square error (RMSE) and adjusted coefficient of determination (R 2 adj ). The best equation with the most suitable predictors for each component was selected according to minimum RMSE and maximum R 2 adj from the models with significant parameter estimates. The statistical significance (p < 0.05) of the estimated coefficients was taken into consideration as well. The best additive model system was reconstructed by combining the selected best biomass equations with the available predictors for each component. Since the diameter is an imperative factor in the field measurement and other factors may be missing in some forest inventory data, an additive model system using D as the only independent variable was considered as an alternative approach for predicting biomass. Both model systems were assessed by the aforementioned criteria. The leave-oneout cross-validation method was further used to test the selected aboveground and root biomass models. Graphical analyses of the predicted vs. observed values were plotted for each component and the prediction performance was evaluated by two statistical indices: mean prediction error (MPE) and mean percent standard error (MPSE). The mathematical expressions of these statistics are as follows: where y j is the j-th observed biomass value;ŷ j is the j-th biomass value estimated by the model; y is the mean of observed biomass value; n is the total number of samples; k is the number of parameters; andŷ j,−j is the j-th value predicted by the model fitting with n − 1 remaining observations, excluding the j-th observation.

Comparison to Existing Biomass Models
The measured biomass data of Mongolian oak in this study was used to compare the relative error (RE) of the existing species-specific equations reported by these studies [8][9][10], which included D and H as predictors. These models are only valid within the range of the independent variables, and the measured data within the specific range are therefore used. Table 2 provides a summary of these existing biomass models.
whereŴ i and W i represent the predicted and observed dry weight of tree i. n is the number of measured trees used for testing the existing equations. HLJ: Heilongjiang province; JL: Jilin province.

Aboveground Biomass Partitioning
There are large differences in the fractions of the different components that account for the aboveground biomass. Aboveground biomass allocated a large proportion to stems (approximately 68%), with a partitioning proportion of 53% and 15% for stem wood and bark, respectively. By contrast, crown biomass possesses a relatively small part of the aboveground biomass, with an average value of 32%, which can be separated into branches and foliage in the proportion of 25% and 7%, respectively. The relationships between the component biomass fractions and diameter, height, crown length and crown width were further examined graphically and displayed in Figure 3. In general, the ratio of wood biomass to aboveground biomass remains constant versus these factors. With the increase in D, H, CW and CL, both bark and foliage showed a decreasing percentage, while a reverse trend was observed for the branch component.

Root-to-Shoot Ratio
Among the individual trees that were sampled for both aboveground and root biomass, the RS ratio varied from 0.11 to 1.06 and was related to diameter ( Figure 4). The RS ratio nonlinearly decreased with increasing D. The mean RS ratio of the 53 individuals was 0.37.

Biomass Models
The nonlinear trend in observed biomass of the stem wood, bark, branch, foliage, aboveground biomass and roots, as a function of D, H, CL and CW, is shown in Figure 5. The fitting results of the aboveground additive systems and root biomass equations based on the seven candidate models are summarized in Table 3. Regression coefficients of the D-based models were statistically highly significant (p < 0.001) for all components, with R 2 adj ranging from 0.754 to 0.920. When adding other variables, some models yielded nonsignificant parameter estimates (p > 0.05) and were, therefore, considered invalid. For example, although both RMSE and R 2 adj showed little improvements in several root biomass models using additional variables, the coefficients were not significant. The equation with D as the only predictor was therefore most suitable for root biomass. In turn, for the aboveground components' biomass models, using additional explanatory variables can provide both significant model parameter estimates (p < 0.05) and better fitting results with a lower RMSE and higher R 2 adj . The stem wood biomass models showed the best fitting effects when using D, H, CW and CL as predictors, which had a lower RMSE (18.548 kg) and higher R 2 adj (0.958). The bark biomass equations showed less error with D, H and CW as predictors, and the corresponding RMSE and R 2 adj were 4.466 kg and 0.946, respectively. In comparison with the improvement of the wood and bark biomass models using H as an additional predictor, the branch and foliage equations yielded larger errors as well as non-significant coefficients when H was included. Both the branch and foliage biomasses were better predicted with smaller estimation errors when using D and CW; the corresponding RMSE and R 2 adj were 19.081 kg and 0.907 for branches and 3.275 kg and 0.765 for foliage, respectively.

Aboveground Biomass Partitioning
There are large differences in the fractions of the different components that account for the aboveground biomass. Aboveground biomass allocated a large proportion to stems (approximately 68%), with a partitioning proportion of 53% and 15% for stem wood and bark, respectively. By contrast, crown biomass possesses a relatively small part of the aboveground biomass, with an average value of 32%, which can be separated into branches and foliage in the proportion of 25% and 7%, respectively. The relationships between the component biomass fractions and diameter, height, crown length and crown width were further examined graphically and displayed in Figure 3. In general, the ratio of wood biomass to aboveground biomass remains constant versus these factors. With the increase in D, H, CW and CL, both bark and foliage showed a decreasing percentage, while a reverse trend was observed for the branch component.

Root-to-Shoot Ratio
Among the individual trees that were sampled for both aboveground and root biomass, the RS ratio varied from 0.11 to 1.06 and was related to diameter ( Figure 4). The RS ratio nonlinearly decreased with increasing D. The mean RS ratio of the 53 individuals was 0.37.

Biomass Models
The nonlinear trend in observed biomass of the stem wood, bark, branch, foliage, aboveground biomass and roots, as a function of D, H, CL and CW, is shown in Figure 5.    The allometric equations considered to be appropriate and available for biomass estimation are displayed in Table 4. Since D is the most important factor in forest inventorying and a desirable predictor in tree growth models, the additive system of biomass equations using D as the single predictor was chosen for potential biomass prediction (System 1). The best additive model system of the aboveground biomass equations was obtained through assembling the selected equation for each biomass component (System 2). Less obvious changes were observed for the biomass equations in System 2 compared to that in their original systems. While in comparison with System 1, the fitting results showed an improvement for both aboveground biomass and its components, especially for stem wood, where the R 2 adj increased from 0.886 to 0.959 and RMSE decreased by 12.191 kg.

Model Cross-Validation
The predicting performance of the appropriate allometric equations, including System 1 (D-based) and System 2 (multivariable-based) for the aboveground biomass components and D-based equation for roots, was assessed through the leave-one-out crossvalidation method. The predicted biomass versus observed biomass for each component is displayed in Figure 6. Biomass equations for wood, bark and aboveground biomass showed small prediction errors (the slopes were close to 1:1 line), whereas relatively large errors were obtained for the branch, foliage and root equations. For each component, the linear fitting slopes for the scatter plot of the predicted against observed biomass were essentially the same for the two additive systems of the biomass equations. The linear fitting effects based on System 2 were better than that based on System 1 for all the aboveground components, with a higher R 2 adj and lower MPE and MPSE (Figure 7), indicating the predicted values of System 2 were closer to the observed values than System 1.

Model Cross-Validation
The predicting performance of the appropriate allometric equations, including System 1 (D-based) and System 2 (multivariable-based) for the aboveground biomass components and D-based equation for roots, was assessed through the leave-one-out crossvalidation method. The predicted biomass versus observed biomass for each component is displayed in Figure 6. Biomass equations for wood, bark and aboveground biomass showed small prediction errors (the slopes were close to 1:1 line), whereas relatively large errors were obtained for the branch, foliage and root equations. For each component, the linear fitting slopes for the scatter plot of the predicted against observed biomass were essentially the same for the two additive systems of the biomass equations. The linear fitting effects based on System 2 were better than that based on System 1 for all the aboveground components, with a higher R 2 adj and lower MPE and MPSE (Figure 7), indicating the predicted values of System 2 were closer to the observed values than System 1. Figure 6. Scatterplots of the predicted (y-axis) and the observed (x-axis) values for the wood, bark, branch, leaf, aboveground and root biomass using the leave-one-out cross-validation method. Figure 6. Scatterplots of the predicted (y-axis) and the observed (x-axis) values for the wood, bark, branch, leaf, aboveground and root biomass using the leave-one-out cross-validation method. Figure 6. Scatterplots of the predicted (y-axis) and the observed (x-axis) values for the wood, bark, branch, leaf, aboveground and root biomass using the leave-one-out cross-validation method.

Discussion
Developing individual-based biomass models and estimating forest biomass has long become a prerequisite for studying productivity and in more recent times for assessing carbon accounting and its potential in mitigating climate change [33,34]. The two additive systems of aboveground biomass equations (Table 4) and the D-based root biomass equation obtained in the present study for two natural oak species (Mongolian oak and Liaodong oak) provide an operational approach for gaining precise biomass estimation in northeastern China.

Biomass Partitioning
The share of biomass of the tree components out of the aboveground biomass is commonly observed when making a comparison among species, diameter classes and tree ages [5,35]. The results of the present study indicated that wood biomass had the largest portion of aboveground biomass, and tree stems (both wood and bark included) accounted for about 68% of the total aboveground, while the crown had relatively less aboveground biomass allocation and foliage represented the smallest portion. The aboveground biomass allocation characteristic is consistent with previous studies for oak [8,9]. The relatively steady contribution of the wood to the aboveground biomass and the increased ratio of branch biomass versus tree size suggest that oak enhances the competitiveness of aboveground resource acquisition, such as light, perhaps mainly through promoting crown area growth to outcompete neighboring competitors. The foliage biomass proportion that decreased against increasing tree size may be due to foliage being borne on the relatively new branches, implying that the production of foliage biomass per unit of woody biomass decreases [5]. The increase in woody mass fraction, including the wood and branch components, and decrease in foliage with increasing tree size obtained in the present study could be explained by the fact that woody mass usually accumulates at the expense of foliage mass [35,36].

Root-Aboveground Biomass Link
Belowground biomass is commonly missing in most biomass datasets and is often indirectly obtained through the RS ratio based on aboveground biomass [22]. The RS ratio decreased significantly with an increase in diameter, namely, tree growth; this trend is similar to other relevant studies [24,25]. The mean value of the RS ratio observed in this study (0.37) is comparatively higher than the global mean value (0.26) reported by Cairns et al. [37] and the mean value (0.29) of broad-leaved species, including oak in China [38]. Mugasha et al. [2] reported a mean value of 0.40 in miombo woodlands, which is slightly higher. Since there is a nonlinear relationship between the RS ratio and diameter, the large difference among the mean RS ratios may be attributed to the sample size in the different diameter classes. Although the mean RS ratios have frequently been accepted for root biomass estimation for purposes such as productivity quantification and carbon accounting, the nonlinear trend indicates that using a fixed RS ratio to predict belowground biomass could introduce a bias [2,39]. Therefore, if root allometric biomass equations are not available, the application of RS ratios should be done with caution.

Allometric Equations
Biomass additivity has long been recognized as a desirable and important property when developing models to predict components and whole tree biomass [40] since the additive property eliminates the logical inconsistency between the sum of the estimated values for the tree components and the estimation provided by the total tree equation [14,41]. However, many biomass models reported to date are non-additive and established using the ordinary least square method [10,13]. The merits of NSUR make it an approach of selection to fit the system of biomass equations, considering the synergetic correlations for each component and correcting the inherent errors among the estimates of each model. In this case, the explanatory variables of the components' biomass models constitute the predictors of the total tree biomass model and the additive property is thus guaranteed through setting parameter restrictions [15,19]. Furthermore, an additive system of biomass equations, when estimated by considering the contemporaneous correlations between the biomass components, has greater statistical efficiency and reduces the confidence and prediction intervals of the biomass estimates, resulting in lower variance [12,15,40].
Tree diameter, the most important index in forest inventories, is always used for developing biomass models, and tree height is often added to the models as an additional predictor [27,41]. Besides these, we also tested two other commonly suggested explanatory variables (CW and CL) in explaining the components' biomass variability. Our results in Table 3 show that the root biomass equation with D as the single explanatory variable performed best with significant model coefficients, although the goodness-of-fit was not optimum. This was consistent with Djomo and Chimi [22] in that D alone was a good predictor for the belowground biomass estimation. As confirmed in many studies, the aboveground biomass models using only D as a predictor also gave good fitting results with significant coefficients for all components [24,42,43]. Using D reduces the limitations in reliability of the H and crown dimension (CW and CL) in most forest inventories, since the height and crown dimension measurements are difficult to measure in the field due to many factors, such as a dense canopy, taller trees and obstructions that hinder seeing the top of the tree [22,44].
Despite the difficulty in accurate measurements, including H as a second predictor in aboveground biomass models is often recommended as it improves the fitting effects [45,46]. In our case, the inclusion of H improved the stem wood and bark biomass estimations but reduced the fitting efficiency for branch and foliage biomass models, in which the coefficients for the H variable were nonsignificant and not always positive (Table 3). These results are in agreement with previous relevant studies [20,[47][48][49]. In general, the crown dimension varies greatly among species, reflecting different biomass allocation strategies [45]. This study found that the inclusion of CW and/or CL slightly improved the model fits for the aboveground components, especially for stem wood and bark, contradicting the results that adding crown radius could greatly improve model performance and reduce error [50].

Heteroscedasticity Elimination
The heteroscedastic residual error is a common problem when developing nonlinear biomass equations with additive error terms [41]. Logarithmic transformation of both the response and explanatory variables is opted for in most cases to deal with the heteroscedasticity of the original arithmetic data, e.g., [13,22,51]. However, the nonlinear transformation of converting arithmetic units to logarithms alters the relationship between the response and explanatory variables and can lead to substantial bias when estimating the allometric model parameters on the original scale through logarithmic linear fitting [52,53]. A correction factor is usually selected to remove the systematic bias generated by converting the estimated logarithmic value to the original scale [54]. While correcting bias, the additive property of the biomass models is simultaneously tending to be destroyed. The weighted NSUR method with a unique weighting function for each equation is likely preferred to address heteroscedasticity during biomass model development, e.g., [30,41,55,56]. In line with previous studies, our study also showed good performance of the weighted NSUR for addressing the heteroscedasticity (Figure 8).

Performance of Existing Equations
Through the literature review, several equations have been established for estimating the components' biomass applicable to this study area [8][9][10]. The biomass measured in this study was used to compute REs from these models and to test their reliability to the present study area (Figure 9). The four models all overestimated stem biomass (including both stem wood and bark), and only model type 2 could keep the mean relative error under 10%. While for branch biomass, all of these models gave underestimation, and the smallest mean RE was obtained by model type 1. Both model type 1 and 2 provided better estimations for foliage biomass, with a mean RE less than 5%, yet model type 3 and 4 yielded relatively large underestimations. Since stems account for a large portion of the aboveground biomass, the distribution of REs for aboveground biomass estimation is similar to the stem component, namely, model type 2 gave the smallest RE (slight underestimation). For belowground biomass, model type 1 provided overestimation while the other three models gained underestimation, and model type 4 provided the closest value to the measured biomass with a mean RE of −0.37%.
Although some equations performed well, most equations provided REs larger than 15%. The small sample sizes might affect the model's predictive ability, especially for model type 3 and 4. Furthermore, biomass data for developing model type 3 were gathered in Heilongjiang province, and that for model type 4 was obtained in Jilin province, while model type 1 and 2 were based on data collected in the two provinces (Table 2). However, in the western Daxing'anling mountains located in Inner Mongolia, no trees were sampled. Since allometric biomass equations vary with forest type, size and stand condition, the application of biomass models not suitable to the specific environment may lead to making large systematic errors from the observed data [26,27,40].  estimations for foliage biomass, with a mean RE less than 5%, yet model type 3 and 4 yielded relatively large underestimations. Since stems account for a large portion of the aboveground biomass, the distribution of REs for aboveground biomass estimation is similar to the stem component, namely, model type 2 gave the smallest RE (slight underestimation). For belowground biomass, model type 1 provided overestimation while the other three models gained underestimation, and model type 4 provided the closest value to the measured biomass with a mean RE of −0.37%. Although some equations performed well, most equations provided REs larger than 15%. The small sample sizes might affect the model's predictive ability, especially for model type 3 and 4. Furthermore, biomass data for developing model type 3 were gathered in Heilongjiang province, and that for model type 4 was obtained in Jilin province, while model type 1 and 2 were based on data collected in the two provinces (Table 2).

Conclusions
Appropriate allometric biomass equations are critical in reducing the uncertainties for obtaining reliable forest biomass estimation, especially for carbon accounting. In the present study, we developed two systems of oak tree biomass equations: one based on tree D only, and another based on multiple variables, including D, H, CW and CL. Both systems guarantee the additive property between the tree components and aboveground biomass. The additive systems of equations were estimated using the NSUR method by considering the inherent errors among the equations and addressed the heteroscedasticity by making each equation have its own weighting function. Diameter alone proved to be a good predictor for all of the components' biomass estimation, especially for root biomass. The system with additional crown dimensions slightly improved the prediction of the aboveground biomass components, but adding CW to the branch and foliage biomass models could yield more improvements than stem wood and bark. The additive system of biomass equations with D, H, CW and CL as predictors can be applied to natural oak species distributed in northeastern China for carbon accounting purposes, while the additive system of equations using D as a single predictor can be used when height and crown information are not available. Both systems of equations must be used with caution when predicting the biomass of trees outside the range of the data used for model development and outside the range of the site conditions.