Additive Allometric Equations to Improve Aboveground Biomass Estimation for Mongolian Pine Plantations in Mu Us Sandy Land, Inner Mongolia, China

: Afforestation is conducive to improving ecosystem service functions and ecosystem di-versity in the Mu Us Sandy Land, however, the important attribute of biomass for Mongolian pine ( Pinus sylvestris var. mongolica Litv.) plantations has yet to be accurately evaluated. This study aimed to develop additive allometric biomass equations for the species and evaluate biomass partitioning patterns within tree components. A total of 131 trees were measured for stem, branch, and leaf biomass by destructively sampling and tree climbing, with the latter as a supplement. For each biomass component, we tested three equations with the diameter at breast ( D ) alone, height ( H ) as additional, and diameter in combination with height ( D 2 H ) as predictors using the weighted least squared method. Weighted nonlinear seemingly unrelated regression was adopted to ﬁt a system of additive allometric biomass equations utilizing the selected equations. A leave-one-out cross-validation method (the jackknife procedure) was used to assess the predictive ability. The biomass partitioning pattern was evaluated by calculating the ratios. The results revealed that the diameter alone is a good predictor for branches and foliage biomass estimates, while the stem requires H included to improve estimation accuracy. Mongolian pine allocates relatively more biomass to the crown (51.4%) compared to the stem (48.6%). Branch biomass fraction increased monotonously with increasing tree size while a reverse trend was observed for foliage. In conclusion, the additive models developed in this study provide a robust biomass estimation and can be extensively used to estimate Mongolian pine forests biomass in Mu Us Sandy Land.


Introduction
Desertification is one of the most serious global environmental problems and adversely affects ecosystem services and socio-economic development [1,2]. The Mu Us Sandy Land lies on the farming-pastoral ecotone [3] and in a transition zone from the Gobi Desert to forest-steppe [4], has been identified as a crucial ecological security barrier in northern China, but is very fragile and sensitive to ecological events [3]. Therefore, it is one of the important development areas in the Three-North Forest Shelterbelt program and artificial afforestation is the key to preventing degradation and stable sand [1,5]. Because Mongolian pine (Pinus sylvestris var. mongolica Litv.) has strong adaptability to cold, drought, and soil infertility, it has been one of the principal species in the program and has been introduced to the Mu Us Sandy land since the 1970s. Over four decades after plantation, the once sandy

Tree Sampling and Biomass Measurement
Developing accurate allometric biomass equations needs high sampling efforts across tree sizes and individuals [13]. However, destructive sampling is strictly restricted in this study area. In the present study, we took advantage of the action of clearing trees for use in construction executed by the local government to collect 116 felled trees. To expand the diameter range, we also collected biomass data by climbing 15 standing trees, with a diameter and height range of 29.1 to 44.2 cm and 9.62-15.25 m, respectively. The 131 trees were randomly picked according to occurring diameter classes. Only trees that had relatively intact bole and crown were selected for biomass measurements. These sampled trees were distributed in four separate locations ( Figure 1) with stand density ranges between 250 and 1225 trees per hectare.
For each selected individual tree, dendrometric variables including diameter at breast height (1.30 m height) and total tree height (from the trunk base to the tip) were firstly measured and recorded with a tape graduated in cm and m, respectively. As displayed in Figure 2, a generally positive relationship between tree height and diameter was

Tree Sampling and Biomass Measurement
Developing accurate allometric biomass equations needs high sampling efforts across tree sizes and individuals [13]. However, destructive sampling is strictly restricted in this study area. In the present study, we took advantage of the action of clearing trees for use in construction executed by the local government to collect 116 felled trees. To expand the diameter range, we also collected biomass data by climbing 15 standing trees, with a diameter and height range of 29.1 to 44.2 cm and 9.62-15.25 m, respectively. The 131 trees were randomly picked according to occurring diameter classes. Only trees that had relatively intact bole and crown were selected for biomass measurements. These sampled trees were distributed in four separate locations ( Figure 1) with stand density ranges between 250 and 1225 trees per hectare.
For each selected individual tree, dendrometric variables including diameter at breast height (1.30 m height) and total tree height (from the trunk base to the tip) were firstly measured and recorded with a tape graduated in cm and m, respectively. As displayed in Figure 2, a generally positive relationship between tree height and diameter was found. Thereafter, a nested regression method that treats the tree crown as a deducible system [20] was used to acquire branch and foliage biomass, and its reliability in biomass estimation has been validated many times [20][21][22]. Branches directly protruding from the trunk were defined as the first order, and those protruding from the first order were defined as the second order and so on. The branch without lateral branches was defined as the branch axis, which was treated as a cone-like minor tree bole for volume calculation using the sectional method. Consequently, the volume of all branches can be calculated by summing the branch axis and their lateral branches. Likewise, the number of leaves could be easily obtained.
according to size. These branch axes were sectioned into segments of 40 cm or 20 cm, and the diameters at both ends of a section were recorded. The leaf number of branch axes was concurrently counted. Finally, regression relationships between volume and leaf number, and diameter were established for branch axes. Based on the established relationships, we calculated the total volume of branches and the amounts of needle leaves synchronously.
A section length of 2 and 1 m for large and small trees, respectively, was applied to determine the stem volume. The diameter at each section position was measured and the top section, i.e., the tip, less than a segment was treated as a cone when estimating the volume. Finally, the total stem volume was obtained by summing the volumes of all sections [23]. Representative samples of branch and foliage, as well as stem discs, were collected and taken to the laboratory for basic characteristics determination. Three thin discs containing barks were taken from 1.30 m, 0.5, and 0.7 tree height of each felled trunk. Similarly, three branch discs were taken from the lower, middle, and upper canopy. Leaf fascicles of different sizes and ages were randomly taken from different trees and branches, amounting to over 1000. Sub-samples of twigs and leaves were stored in sealed bags to maintain moisture. In the laboratory, stem and branch samples were firstly saturated and then cut into moderate blocks to measure volume by the water displacement method. All samples including foliage were oven-dried at 75 °C to a constant weight. Basic density (g•cm −3 ) was then inferred by dividing dry weight from fresh volume. The total biomass In the actual measurement, we firstly measured the basal diameters of all the first-order branches, and then measured the basal diameters of the second-order branches, and the same for other orders. It is noteworthy that only the first-order branches were measured for the standing trees. Subsequently, over 200 branch axes were randomly collected according to size. These branch axes were sectioned into segments of 40 cm or 20 cm, and the diameters at both ends of a section were recorded. The leaf number of branch axes was concurrently counted. Finally, regression relationships between volume and leaf number, and diameter were established for branch axes. Based on the established relationships, we calculated the total volume of branches and the amounts of needle leaves synchronously.
A section length of 2 and 1 m for large and small trees, respectively, was applied to determine the stem volume. The diameter at each section position was measured and the top section, i.e., the tip, less than a segment was treated as a cone when estimating the volume. Finally, the total stem volume was obtained by summing the volumes of all sections [23].
Representative samples of branch and foliage, as well as stem discs, were collected and taken to the laboratory for basic characteristics determination. Three thin discs containing barks were taken from 1.30 m, 0.5, and 0.7 tree height of each felled trunk. Similarly, three branch discs were taken from the lower, middle, and upper canopy. Leaf fascicles of different sizes and ages were randomly taken from different trees and branches, amounting to over 1000. Sub-samples of twigs and leaves were stored in sealed bags to maintain moisture. In the laboratory, stem and branch samples were firstly saturated and then cut into moderate blocks to measure volume by the water displacement method. All samples including foliage were oven-dried at 75 • C to a constant weight. Basic density (g·cm −3 ) was then inferred by dividing dry weight from fresh volume. The total biomass of stems and branches was determined by multiplying each component volume based on their specific basic density. Similarly, we acquired total foliage biomass by multiplying the leaf amount (number of needles) by the average dry weight per needle leaf. Overall, tree aboveground biomass was expressed as the sum of all components, i.e., the stem, branch, and foliage ( Table 1). The stump biomass was not included since the workload to dig it out was too great. The biomass partitioning was made by dividing the observed component biomass by aboveground biomass.

Allometric Biomass Model Establishment
Assembled data for dendrometric variables like D, H, and component biomass were used to develop biomass regression equations. The most common mathematical power function y = β 0 x β 1 1 x β 2 2 · · · x β n n was adopted for biomass estimation. We preferred utilizing the nonlinear form with an additive error structure to avoid the system bias caused by log-transformed linear models [24,25]. Since biomass data always exhibit significant heteroscedasticity in model residuals [17,24], weightings were introduced to generate constant error variances over all observations in this study. A specific weighting function defined as w = 1/D Q was applied to correct each observation before parameter estimation. The power exponent Q was computed by the residuals of independent models fitted by OLS regression.
For each component, the nonlinear power model was tested to relate biomass with dendrometric variables including D, H, and D 2 H as defined in Equations (1)-(3).
where W i represents the i-th observed component biomass, β i is the regression coefficient, and ε i is the additive error term. The independent weighted least square (WLS) method was applied to each equation of each component. The "best model", i.e., the most suitable predictor for each component was explored and determined by two fit statistics, including the adjusted coefficient of determination (R 2 adj ) and the root mean square error (RMSE). Only equations with statistically significant coefficients (p < 0.05) were eligible to be selected. The best models are those that have the highest R 2 adj as well as the lowest RMSE. To realize the additive property between tree components, a system of equations for aboveground biomass and its components (stem, branch, and foliage) were simultaneously fitted utilizing the models selected above. The model for total aboveground biomass was equated by adding the equations for each specific component together [16,18]. Therefore, we defined the system of equations as displayed in Model (4): where W st , W br , W fo , and W ag are estimates of the stem, branch, foliage, and aboveground biomass, respectively; f (X i , β i ) is the selected nonlinear equation of the i-th component; ε i is the additive error term. The procedure of WNSUR was used to fit the system of equations [7,18], using the PROC MODEL procedure in SAS 9.3 (SAS Institute Inc., Cary, NC, USA). The WNSUR method considers the correlations between regressions residuals and set constraints on regression coefficients, resulting in lower variance and additive predictions [18].

Model Evaluation and Validation
The behavior of biomass estimators was evaluated by R 2 adj and RMSE, as described above. The formulas of the three goodness-of-statistics were presented as Equations (5) and (6). The predictive performance of the additive system of biomass equations was tested via Turkey's jackknife procedure, a leave-one-out cross-validation method. Graphical analyses of predicted versus observed values of biomass were conducted and the linear fitting statistic of R 2 adj was achieved to evaluate the model goodness-of-fit. A paired t-test was utilized to check the consistency between predicted and observed values. And the crossvalidation statistics of mean absolute error (MAE), mean prediction error (MPE), and total relative error (TRE) were used to diagnose the predictive ability. In addition, the measured biomass data in this study was used to compare with previous biomass equations. These statistics were displayed as follows: where y j andŷ j represent the actual and estimated biomass, respectively; y is the average value of the actual biomass;ŷ j,−j represents the j-th predicted value by the biomass equation fitted with the n-1 remaining observations; n and p are the number of observations and model coefficients, respectively.

Allometric Biomass Models
For the underlying estimation of biomass, allometric relationships were firstly examined between three size variables (D, H, D 2 H) and biomass of stem, branch, foliage, and aboveground total ( Figure 3). Poorer correlations were found between biomass values and tree heights. The results of the biomass equations independently fitted through weighted regression are summarized in Table 2. In all, equations for total aboveground biomass and its constituents gave good fitting effects with R 2 adj greater than 0.91, and all model coefficients were significant (p < 0.05). For the stem component, the equation with D alone as its predictor already achieved a satisfactory goodness-of-fit performance statistic. Adding H as an additional predictor or combining it with D as a single compound predictor could improve the estimation with R 2 adj increased to over 0.99 and RMSE decreased by about half. Branch and foliage biomass equations, on the contrary, showed the highest R 2 adj and lowest RMSE when D was used as the sole predictor. The aboveground total was best estimated when using D and H as multiple variables, which had higher R 2 adj and lower RMSE values. as an additional predictor or combining it with D as a single compound predictor could improve the estimation with R 2 adj increased to over 0.99 and RMSE decreased by about half. Branch and foliage biomass equations, on the contrary, showed the highest R 2 adj and lowest RMSE when D was used as the sole predictor. The aboveground total was best estimated when using D and H as multiple variables, which had higher R 2 adj and lower RMSE values.  Table 2. Regression coefficients with standard errors (in parentheses) and fitting statistics through the independently fitted procedure using weighted regression. With the same weighting functions adopted in WLS regression, the additive model system based on the considered appropriate independently fitted biomass equations were presented in Table 3. The total aboveground biomass equation was obtained by combining the best equation of each component. The equations for the stem, foliage, and aboveground total showed statistics that fit a little better with higher R 2 adj and lower RMSE values compared to the independent fitting results. Simultaneously fitted allometric biomass equations for branch indicated the goodness-of-fit statistics were slightly poorer than that for independently fitted ones. Applying weightings to the additive system of equations could adequately result in homoscedasticity of the variance errors, as shown by the graphical analyses of the residuals (Figure 4).  With the same weighting functions adopted in WLS regression, the additive model system based on the considered appropriate independently fitted biomass equations were presented in Table 3. The total aboveground biomass equation was obtained by combining the best equation of each component. The equations for the stem, foliage, and aboveground total showed statistics that fit a little better with higher R 2 adj and lower RMSE values compared to the independent fitting results. Simultaneously fitted allometric biomass equations for branch indicated the goodness-of-fit statistics were slightly poorer than that for independently fitted ones. Applying weightings to the additive system of equations could adequately result in homoscedasticity of the variance errors, as shown by the graphical analyses of the residuals (Figure 4).

Model Validation by a Jackknife Procedure
The additive system of biomass equations for total aboveground and its constituents was tested by the jackknife method, and relationships between the observed and predicted biomass values all showed a close correlation with R 2 adj larger than 0.92 ( Figure 5). The linear fitting results identified the linear trend of y = x, especially for the stem component and aboveground total. The paired t-test demonstrated that the predicted and observed biomass were statistically consistent for all components (p > 0.05). The predictive ability was also supported by the cross-validation statistics which showed relatively little bias ( Table 4). The corresponding TRE values were negative, resulting in conservative estimates.

Model Validation by a Jackknife Procedure
The additive system of biomass equations for total aboveground and its constituents was tested by the jackknife method, and relationships between the observed and predicted biomass values all showed a close correlation with R 2 adj larger than 0.92 ( Figure 5). The linear fitting results identified the linear trend of y = x, especially for the stem component and aboveground total. The paired t-test demonstrated that the predicted and observed biomass were statistically consistent for all components (p > 0.05). The predictive ability was also supported by the cross-validation statistics which showed relatively little bias ( Table 4). The corresponding TRE values were negative, resulting in conservative estimates.

Aboveground Biomass Partitioning
The relative proportion of component biomass varied against diameter in different ways ( Figure 6). The results showed an increasing trend in the stem biomass fraction over small trees and a slightly declining trend over large trees. In contrast, the branch component showed a monotone increasing percentage while the foliage component had a monotone declining percentage with tree size. Aboveground biomass exhibited relatively larger biomass allocation to stems with a value of approximately 48.6% when compared with branches and foliage. The partitioning proportion was the lowest for foliage biomass (23.7%) and intermediate for branch biomass (27.7%).

Aboveground Biomass Partitioning
The relative proportion of component biomass varied against diameter in different ways ( Figure 6). The results showed an increasing trend in the stem biomass fraction over small trees and a slightly declining trend over large trees. In contrast, the branch component showed a monotone increasing percentage while the foliage component had a monotone declining percentage with tree size. Aboveground biomass exhibited relatively larger biomass allocation to stems with a value of approximately 48.6% when compared with branches and foliage. The partitioning proportion was the lowest for foliage biomass (23.7%) and intermediate for branch biomass (27.7%).

Aboveground Allometry Relationships
The best statistics of goodness-of-fit (R 2 adj and RMSE) were found for the stem biomass equations, while relatively poor statistics were observed for branch and foliage components. Because of the influences of internal and external factors, branch and foliage biomass is naturally more variable than other parts like stems [15,26]. Despite this, the equations presented in this study for predicting the biomass of stem, branch, foliage, and their aggregates provided a satisfactory accuracy level.
Allometric biomass equations with D as the sole explanatory variable are easy to build and employ, and usually provide satisfactory biomass estimates for many species and regions [27][28][29]. Our results also indicate that the regression models utilizing D could yield good fitting results for component biomass and their aggregates. Moreover, D is very easy to accurately measure in the field and is almost involved in all forest inventory data. The inclusion of H as an additional predictor improved equations for stem and aboveground biomass but led to a decline in the goodness-of-fit for branch and foliage biomass. Similar results have been found as well in previous studies [7,22]. This probably suggests that tree height may be essential for stem biomass but may not be a good predictor for crown biomass. In a way, tree heights capture the variation of volume among trees [30]. The model using single compound variables D 2 H showed less efficient estimations than the one based on the multiple variables of D and H.
Due to the absence of information about tree crown size, we developed biomass equations using D and H only. In general, crown size variables including crown length and crown diameter, have been found to further improve crown biomass estimations [11,17,31]. However, crown size measurements from the ground are a time-consuming and laborious task, and crown overlap may increase the difficulty in measurements, hence crown dimensions were rarely available in forest inventory data [24]. In practice, equations with easily measured metrics are more practical to estimate forest aboveground biomass.
As reported in published studies, the NSUR procedure results in greater statistical efficiency and lower variance by considering the contemporaneous correlations between different biomass components [18,32,33]. In this study, the simultaneous fitting of the system of biomass equations generated little improvement in terms of R 2 adj and RMSE for the stem, foliage, and aboveground total compared to the independently fitted equations. Similar results were demonstrated by previous researchers who found small differences in the fitting statistics between independent and simultaneous fitting procedures that applied to the total biomass and its constituent parts [15,34]. But actually, the system of equations fitted simultaneously ensures the additivity, which has long been considered a desirable property when modeling the biomass estimates of the total and its constituent parts [16,18]. The additivity property conforms to the biological point of view that the estimates of biomass components and total are consistent [34]. In addition, simultaneous fitting of equations usually provides narrower confidence intervals relative to the equations fitted independently, since the model for total biomass must be a combination of that for component biomass resulting in parameter restrictions [19,34]. Therefore, considering the model sensitivity, fitting efficiency, and biological consistency, the estimates are more robust by the simultaneous fitting procedure [16].
Heteroscedasticity is almost a certainty in allometric biomass model development, as already discussed in previous studies [14,17,18,35]. The fact that biomass variability around the equations increased with the size of the independent variables can be considered general and natural, resulting in the heterogeneity of variance of biomass models [16,24,36]. When developing biomass equations, the residuals should satisfy the normality, and as a result, heterogeneous variances cannot be ignored and must be dealt with to achieve stable and minimum variance. The direct application of the OLS technique to biomass data, in this case, does not remain efficient and reliable because of the unstable variance [16].
An alternative method to tackle this daunting issue is to log-transform both response and explanatory variables before model construction to comply with normality and homoscedasticity [7,37,38]. Approximately 66% of the biomass equations were analyzed using OLS log-linear regression [39]. However, the application of log-transformation may alter the relationship between predictor and response variables and may arouse substantial bias if the parameter estimation of an allometric equation on its original scale is estimated by linear fitting of logarithms [40]. A more efficient way is to model the error structure as a weighting function to correct heteroscedasticity in biomass data, as recommended by several studies [16,18,19]. The weighted nonlinear equations can achieve relatively low errors compared to OLS log-linear equations [25]. In this study, the variances were modeled as a power-law function of predictor D, stabilizing the model variances for all biomass components (Figure 4). Thereby, the WNSUR procedure is more suitable for fitting biomass equations since it can provide accurate estimates with lower uncertainty than those obtained from independently fitted unweighted equations, corresponding with previous analyses of the performance of the equations, the biological consistency, and statistical efficiency [16].

Model Comparison and Biomass Allocation
An additive system of equations was developed for estimating the aboveground biomass of Mongolian pine trees in Northeastern China [28]. The biomass measured during the present study was used to test the existing equations' reliability on this study site. Since biomass models are valid only within the range of explanatory variables, the sampled trees that fell in the specific range were employed. Although the predicted and observed values showed good linear fitting results with R 2 adj larger than 0.88, the published equations did not provide accurate biomass estimates at the tree level, especially for individual components (Figure 7). By contrast with the 1:1 line, we found these equations overestimated the biomass for stem but underestimated for branches and foliage, resulting in relatively and seemingly effective biomass estimates for the aboveground total. The different biomass sampling methods (volume conversion versus direct weighing) may be the potential causes. But the primary reason may be that the biomass allocation pattern is driven by environmental conditions [41,42].
mended by several studies [16,18,19]. The weighted nonlinear equations can achieve relatively low errors compared to OLS log-linear equations [25]. In this study, the variances were modeled as a power-law function of predictor D, stabilizing the model variances for all biomass components ( Figure 4). Thereby, the WNSUR procedure is more suitable for fitting biomass equations since it can provide accurate estimates with lower uncertainty than those obtained from independently fitted unweighted equations, corresponding with previous analyses of the performance of the equations, the biological consistency, and statistical efficiency [16].

Model Comparison and Biomass Allocation
An additive system of equations was developed for estimating the aboveground biomass of Mongolian pine trees in Northeastern China [28]. The biomass measured during the present study was used to test the existing equations' reliability on this study site. Since biomass models are valid only within the range of explanatory variables, the sampled trees that fell in the specific range were employed. Although the predicted and observed values showed good linear fitting results with R 2 adj larger than 0.88, the published equations did not provide accurate biomass estimates at the tree level, especially for individual components (Figure 7). By contrast with the 1:1 line, we found these equations overestimated the biomass for stem but underestimated for branches and foliage, resulting in relatively and seemingly effective biomass estimates for the aboveground total. The different biomass sampling methods (volume conversion versus direct weighing) may be the potential causes. But the primary reason may be that the biomass allocation pattern is driven by environmental conditions [41,42].  As displayed in Figure 6, the stem alone accounted for approximately 48.6% of the aboveground biomass on average, whereas branches accounted for 27.7% and the foliage for 23.7%. The partitioning pattern within the tree indicates that Mongolian pine in dryland allocates the greatest fraction of its resources to crown (branches and foliage included) development, then to the stem. Interestingly, the pattern of aboveground biomass allocation among components largely differed from that reported by Dong, et al. [28] who highlighted that stem biomass occupied over three-quarters of the aboveground total, while branches and foliage had only 13% and 9%, respectively. The discrepancies in partitioning pattern could be explained due to the site quality, climate, and competition among trees. Previous studies have demonstrated that many tree species develop broad crowns but short trunks in arid regions [43,44]. The hydraulic resistance is considered the most severe limitation to tree height growth [45]. The drier the site conditions the sooner upward growth will be stopped [44]. In this study area, the low precipitation makes it drier than northeastern China, Mongolian pine produces a rather shorter tree height. For example, a tree with a diameter of 38.3 cm has a height of 22.3 m in northeast China [28], but a tree with a 44.2 cm diameter only has a height of less than 14 m in Mu Us Sandy Land. Due to the low soil water availability in dryland, upper branches are limited by hydraulic resistance as well. The middle and lower lateral branches become the principal part of photosynthesis. However, the lower branches may be shaded and limited by light. In this issue, the middle branches may grow more and then develop a broad crown [43]. A spreading crown of a short tree could also reduce wind load and help to enhance wind resistance [43,44]. In addition, the relatively low stand density, i.e., the weak competition in this study, may promote crown growth as well.
As usual, the stem is always the largest biomass pool among different components [46]. Although the results in this study showed a similar allocation pattern, the stem holds a relatively small portion of aboveground biomass relative to other studies [11,46]. The proportional distribution of biomass among tree components changes with tree size, namely tree diameter has a strong effect on aboveground biomass allocation patterns. More specifically, foliage mass declined substantially while branches increased continuously as tree size increased, which is consistent with other studies [7,11,42]. This indicates that the accumulation of foliage biomass per unit of woody mass tends to decrease as trees grow larger, corresponding to the fact that woody biomass increases generally at the cost of foliage biomass [42,47]. When trees are small and at the early development stage, more carbohydrates are allocated to foliage to maintain efficient photosynthesis and promote vertical growth. As the trees grow big, lateral branches of different individuals contact each other, resulting in light competition between neighboring trees. Thus, more resources are allocated to stems allowing for primary and secondary growth to enable trees to compete with neighbors [42]. In the meantime, additional carbohydrates are invested for branch growth. Trees with increasing biomass investment in branches keep a competitive advantage that allows them to out-compete neighbors by developing broad crowns to fight for light and available water [7]. Indeed, the biomass partitioning patterns stress the plasticity of trees growing in different environmental conditions in nutrient acquisition and trait formation. Therefore, the application of site-specific and species-specific allometric biomass equations is recommended to the corresponding sites whenever possible [30].

Conclusions
This study recommends using the additive systems of equations to estimate the aboveground biomass of Mongolian pine forests in Mu Us Sandy Land. Diameter at breast height alone is a perfect predictor for branch and foliage biomass estimation, while the stem requires tree height as an additional predictor to improve estimation accuracy. Mongolian pine species in sandy land allocate relatively more biomass in the crown, as a possible mechanism to improve water absorption ability. Given the plasticity of trees under different environments, the allometric equations developed in this study can be applied to make accurate biomass estimates reliable to the corresponding sites. The applicability of these equations should restrict an independent variable to the range of the data used for model development in this study.