#### Abstract

Three methods of enforcing additivity of tree component biomass estimates into total tree biomass estimates for *Androstachys johnsonii* Prain were studied and compared, namely, the conventional (CON) method (a method that consists of using the same independent variables for all tree component models, and for total tree model, and the same weights to enforce additivity), seemingly unrelated regression (SUR) with parameter restriction, and nonlinear seemingly unrelated regression (NSUR) with parameter restriction. The CON method was found to be statistically superior to any other method of enforcing additivity, yielding excellent fit statistics and unbiased biomass estimates. The NSUR method ranked second best but was found to be biased. The SUR method was found to be the worst; it exhibited large bias and had a poor fit for the biomass. Therefore, we recommend that only the CON and NSUR methods should be used for further estimates, provided that their limitations are considered, that is, exclusion of contemporaneous correlations for the CON method and consideration of the significant bias of the NSUR method.

#### 1. Introduction

In the early 1960s,* Androstachys johnsonii* Prain (*A. johnsonii*) had already been reported to be almost completely restricted to Mozambique [1], presumably due to overexploitation. Five decades later, there is still a lack of studies on this species in any branch of forest science, particularly tree allometry, supporting that this species may indeed be restricted to Mozambique. Moreover, in Mozambique, biomass models have rarely been reported for any tree species, and no such models have been described for* A. johnsonii*. Isolated studies in Mozambique have focused on Miombo woodlands and some forest plantations and are often part of research theses or forestry projects; therefore, the results are not always made public. Generally, such studies have only considered the aboveground biomass and have not included a breakdown of further tree components.

*A. johnsonii* woodlands (Mecrusse) are very important. Besides being restricted to Mozambique [1], Mecrusse has important socioeconomic value to local communities, which use stakes of* A. johnsonii* in the construction of homes, shelters, and furniture and sell them for income generation. At the global scale, Mecrusse is reported to be a tipping point in regional ecological and socioeconomic development [2], hence the importance of modelling and estimating its biomass.

The estimation of aboveground biomass is important to predict the amount of carbon that is sequestered [3–5], to assess nutrient cycling and fluxes and energy wood potentials [4, 6], and to provide estimates for the different tree components [5]. These types of estimates are important for several reasons as follows: (i) stem wood biomass is an important quantity because this component is the only one used in the forest industry, and the carbon therefore remains stored for a long time and is not released into the atmosphere; (ii) in many species, branches and foliage are left in the forest and decompose, releasing CO_{2} and nutrients; (iii) in some species, especially broadleaf species, the branches are collected by members of local communities for use as firewood, which will result in release of CO_{2}; (iv) the stump and root system are left in the forest, allowing the stump to either sprout (regrow), continuing the sequestration process, or decompose along with the roots, releasing CO_{2} and nutrients; and (v) in some tree species, belowground biomass can account for more than one-third of the total biomass [7]. Hence, it is critical to estimate the biomass of all tree components as well as the total tree biomass in order to assess the global carbon balance.

However, the biomass estimates of the considered tree components often do not sum to the estimate of the total tree biomass, and a desired and logical feature of the tree component regression equations is that the predictions of the components sum to the prediction for the total tree. This feature is called additivity. Various authors, such as Goicoa et al. [5], Kozak [8], Cunia [9], Cunia and Briggs [10, 11], Jacobs and Cunia [12], Parresol [4, 13], and Carvalho and Parresol [14], have proposed and/or discussed various methods to ensure the property of additivity.

The objective of this study was to fit independent linear and nonlinear tree component and total tree biomass models and compare three methods of enforcing the property of additivity (the conventional (CON) method, seemingly unrelated regression (SUR) with parameter restriction, and nonlinear seemingly unrelated regression (NSUR) with parameter restriction) for* A. johnsonii* tree species.

The CON method consists of using the same independent variables for all tree component models, and for total tree model, and the same weights to enforce additivity [4]. The SUR and NSUR methods consist in first fitting and selecting the best linear and nonlinear models, respectively, for each tree component. The total tree model is a function of the independent variables used in each component model. Then, all models, including the total, are fitted again simultaneously using joint-generalized least squares under the restriction of the coefficients of regression which ensures the additivity property [4].

#### 2. Materials and Methods

##### 2.1. Study Area

Mecrusse is a forest type where the main species, many times the only one, in the upper canopy is* A. johnsonii* [15].

In Mozambique, Mecrusse woodlands are mainly found in Inhambane and Gaza provinces and in Massangena, Chicualacuala, Mabalane, Chigubo, Guijá, Mabote, Funhalouro, Panda, Mandlakazi, and Chibuto districts. The eastern-most Mecrusse patches, covering the last five districts, were defined as the study area. The study area had an extension of 4,502,828 ha [16], of which 226,013 ha (5%) was covered by Mecrusse woodlands.

The climate is dry and tropical throughout the study area, except in the west part of the Panda district and the southwest part of the Mandlakazi district, where the climate is humid and tropical [16–21]. The climate has two seasons: the warm or rainy season from October to March and the cool or dry season from March to September [17–21].

The mean annual temperature is generally greater than 24°C, and the mean annual precipitation varies from 400 to 950 mm [16–21]. According to FAO classification [22], the soils in the study area are mainly Ferralic Arenosols covering more than 70% of the study area [16]. Arenosols, Umbric Fluvisols, and Stagnic soils are also predominant in the northern-most part of the study area [16].

The study area is characterized by a shortage of water resources as well as precipitation; thus, of the five districts that made up the study area, only the districts of Chibuto and Mandlakazi have water resources [16–21], either from precipitation or from lakes and rivers.

##### 2.2. Data Acquisition

A total of 93 trees (2 to 6 per plot) selected across all size classes (Table 1) were destructively sampled within 23 circular plots randomly located in the study area. Diameter at breast height (DBH), total tree height (), crown height (CH), and live crown length (LCL) were measured on the felled trees. Trees were divided into the following tree components: (1) root system, (2) stem wood, (3) stem bark, and (4) crown. Tree components were sampled and the dry weights were estimated as follows.

###### 2.2.1. Root System

The stump height was predefined as being 20 cm for all trees and considered as part of the taproot, as recommended by Parresol [13] and because in larger* A. johnsonii* trees this stump height (20 cm) is affected by the root buttress; therefore, the root collar was also considered part of the taproot. The root system was divided into 3 subcomponents: fine lateral roots, coarse lateral roots, and taproot. Lateral roots with diameters at insertion point on the taproot < 5 cm were considered as fine roots and those with diameters ≥ 5 cm were considered as coarse roots.

First, the root system was partially excavated to the first node, using hoes, shovels, and picks, to expose the primary lateral roots (Figures 1(a)–1(c)). The primary lateral roots were numbered and separated from the taproot with a chainsaw (Figures 1(b) and 1(c)) and removed from the soil, one by one. This procedure was repeated in the subsequent nodes until all primary roots were removed from the taproot and the soil. Finally, the taproot was excavated and removed (Figures 1(d)–1(f)). The complete removal of the root system was relatively easy because 90% of the lateral roots of* A. johnsonii* are located in the first node, which is located close to ground level (Figures 1(a)–1(d)); the lateral roots grow horizontally to the ground level, and do not grow downwards; and because the taproots had, at most, only 4 nodes and at least 1 node (at ground level).

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

Fresh weight was obtained for the taproot, each coarse lateral root and for all fine lateral roots. A sample was taken from each subcomponent, fresh weighed, marked, packed in a bag, and taken to the laboratory for oven drying. For the taproot, the samples were two discs, one taken immediately below the ground level and another from the middle of the taproot. For the coarse lateral roots, two discs were also taken, one from the insertion point on the taproot and another from the middle of it. For fine roots the sample was 5 to 10% of the fresh weight of all fine lateral roots. Oven drying of all samples was done at 105°C to constant weight (i.e., to, approximately, 0% moisture content), hereafter, referred to as dry weight.

###### 2.2.2. Stem Wood and Stem Bark

Felled trees were scaled up to a 2.5 cm top diameter. The stem was defined as the length of the trunk from the stump to the height that corresponded to 2.5 cm diameter. The remainder (from the height corresponding to 2.5 cm diameter to the tip of the tree) was considered a fine branch. The stem was divided into sections, the first with 1.1 m length, the second with 1.7 m, and the remaining with 3 m, except the last, whose length depended on the length of the stem. Discs were removed on the bottom and top of the first section and on the top of the remaining sections; that is, discs were removed at heights of 0.2 m (stump height), 1.3 m (breast height), and 3 m, and the successive discs were removed at intervals of 3 m to the top of the stem, and their fresh weights were measured using a digital scale.

Diameters over and under bark were taken from the discs in the North-South direction (previously marked on the standing tree) with the help of a ruler. The volumes over and under the bark of the stem were obtained by summing up the volumes of each section calculated using Smalian’s formula [27, 28]. Bark volume was obtained from the difference between volume over bark and volume under bark.

The discs were dipped in drums filled with water for its saturation (3 to 4 months) and subsequent determination of the saturated volume and basic density. The saturated volume of the discs was obtained based on the water displacement method [29] using Archimedes’ principle. This procedure was done twice: before and after debarking; hence, we obtained saturated volume under and over the bark.

Wood discs and respective barks were oven dried at 105°C to constant weight. Basic density was obtained by dividing the oven dry weight of the discs (with and without bark) by the relevant saturated wood volume [27, 30]. Therefore, two distinct basic densities were calculated: (1) basic density of the discs with bark and (2) basic density of the discs without bark.

We estimated the basic density at point of geometric centroid of each section using the regression function of density over height [31]. This density value was taken as representative of each section [31].

###### 2.2.3. Crown

The crown was divided into two subcomponents: branches and foliage. Primary branches, originating from the stem, were classified in two categories: primary branches with diameters at the insertion point on the stem ≥ 2.5 cm were classified as large branches, and those with diameters < 2.5 cm were classified as fine branches.

Large branches were sampled similarly to coarse roots, and fine branches and foliage were sampled similarly to fine roots. All the leaves from each tree were collected and fresh weighed together and a sample was taken for oven drying. The subcomponents branches and foliage were not treated as separated components because in the preliminary analysis the weight of the foliage did not show significant variation with DBH,* H*, CH, and LCL, exhibiting, therefore, poor fits.

###### 2.2.4. Tree Component Dry Weights

Dry weights of coarse and fine roots, large and fine branches, and foliage were obtained from the “fresh weight/oven dry weight” ratio of the respective samples by multiplying it by the relevant subcomponent total fresh weight. Dry weights of the root system and crown were obtained by summing up the relevant subcomponents’ dry weights.

Dry weights of each stem section (with and without bark) were obtained by multiplying respective densities by relevant stem section volumes. Stem (wood + bark) and stem wood dry weights were obtained by summing up each section’s dry weight with and without bark, respectively. The dry weight of the stem bark was obtained from the difference between the dry weights of stem and stem wood. Finally, the total tree biomass was obtained by adding the component dry weights.

##### 2.3. Data Analysis

Several linear and nonlinear regression model forms were tested for each tree component and for the total tree using weighted least squares (WLS). The weight functions were obtained by iteratively finding the optimal weight that homogenises the residuals and improves other fit statistics. Independent tree component models were fitted with the statistical software package R [32] and the functions lm and nls for linear models and nonlinear models (the latter of which using the Gauss-Newton algorithm). The best linear and nonlinear biomass equations selected are given in (1) and (2), respectively. Among the tested weight functions (, , , , , and ), the best weight function was found to be , for all tree component equations (linear or nonlinear). Although, the selected weight function might not be the best one among all possible weights, it is the best approximation found. ConsiderThe CON method used the same independent variables for all tree component models and the total tree model and used the same weight functions [4], achieving additivity automatically [5]. For this method, the most frequent best linear model form in (1) among tree components was used for all other components and for total tree biomass. The most frequent model form in (1) is ; therefore, the structural system of equations for the CON method is given inThe SUR method consisted of first fitting and selecting the best linear models for each tree component. The total tree model was a function (sum) of the independent variables used in each tree component model. Then, all models, including the total, were fitted again simultaneously using joint-generalized least squares (also known as SUR) under the restriction of the coefficients of regression, which ensured additivity.

The best linear model forms were found to be for belowground, stem wood, and stem bark biomasses and for the crown biomass. Summing up the best model forms from each tree component, the model form obtained for the total tree biomass was .

However, the system of equations obtained by combining the best linear model forms per component under parameter restriction will not yield effective and precise estimates because, according to SAS Institute Inc. [33], for SUR to be effective, the models must use different regressors. This requirement is not verified, as three of the four components have identical regressors. Indeed, according to Srivastava and Giles [34], applying SUR to system of the best equations given above is of no benefit when the component equations have identical explanatory variables. Moreover, as stated by Greene [35] and Bhattacharya [36], a system of linear SUR equations with identical regressors yields ineffective estimates of coefficient vectors when compared to equation-by-equation ordinary least squares (OLS).

To eliminate the ineffectiveness caused by identical regressors, SUR was applied using second best regression equations for belowground and stem wood biomasses such that the different tree component equations could have different regressors. The resulting system of equations of biomass additivity is given in (4). However, the results of SUR using the best independent model forms are given in Tables 2 and 3, for demonstration proposes of the ineffectiveness caused by identical regressors. ConsiderNote from the equations in (4) that the intercepts of all tree component biomass models are forced (constrained, restricted) to sum to the intercept of the total tree biomass model, the coefficients of regression for the regressor in the root system and stem wood biomass models are constrained to sum to the coefficient of regression for in the total tree biomass model, and the coefficients for the regressors , , and in the root system, stem bark, and crown biomass models, respectively, are constrained to be equal to the coefficients of the same regressors in the total tree biomass model, thereby achieving additivity.

The NSUR method had the same characteristics and was performed using the same procedures as the SUR method except that the system of equations was composed of nonlinear models. For reference, please see Brandeis et al. [3], Parresol [13], Carvalho and Parresol [14], and Carvalho [37]. The system of equations (including the total tree biomass) obtained by combining the best nonlinear model forms per component under parameter restriction is given byNote that the coefficients of regression of each regressor in each tree component model are forced (constrained, restricted) to be equal to coefficients of the equivalent regressor in total tree model, allowing additivity.

The systems of equations in (4) and (5) were fitted using PROC SYSLIN and PROC MODEL in SAS software [33], respectively, using the ITSUR option. Restrictions (constraints) were imposed on the regression coefficients by using SRESTRICT and RESTRICT statements in PROC SYSLIN and PROC MODEL procedures, respectively. The start values of the parameters in PROC MODEL were obtained by fitting the logarithmized models of each component in Microsoft Excel.

##### 2.4. Model Evaluations and Comparison

The best tree component and total tree biomass equation were selected by running various possible regressions on combinations of the independent variables (DBH, , and LCL) and evaluating them using the following goodness of fit statistics: adjusted coefficient of determination (Adj.), standard deviation of the residuals () and CV of the residuals, mean relative standard error (MRSE), mean residual (MR), and graphical analysis of the residuals. The computation and interpretation of these fit statistics were previously described by Goicoa et al. [5], Gadow and Hui [38], Meyer [39], Magalhães [40], and Ruiz-Peinado et al. [41]. The best models are those with highest Adj., smallest , and CV of the residuals, MRSE, and MR and with the residual plots showing no heteroscedasticity, no dependencies or systematic discrepancies.

In addition to the goodness of fit statistics described above, the methods of enforcing additivity were compared using percent standard error of the expected value and percent standard error of the predicted value, as computed in Table 4. The smaller the percent standard error of the expected and percent standard error of the predicted values is, the better the model is in predicting the biomass.

SUR and NSUR methods were used instead of, for example, simply summing the best component biomass models (i.e., Harmonization procedure [42]), because in the latter case the total biomass is not modelled and therefore its fit statistics are unknown and because the sum of tree component models with the best fits does not guarantee good fit in the total model and might produce biased estimates for whole tree biomass [6] and, further, because SUR and NSUR, unlike the CON method, take into account the contemporaneous correlation among residuals of the component equations [4, 13, 14, 24].

Nevertheless, the standard deviation and CV of the residuals for the harmonization approach (HAR) were compared with those obtained for SUR and NSUR approaches. Since, in HAR procedure, the total tree biomass is obtained simply by summing the best component models, the standard deviation of the residuals can be computed using the variance of a sum (6) [4, 13]. Consider where and are the standard deviation of the residuals of the total tree biomass model and of the th tree component biomass model and is the covariance of th and th tree component biomass models.

The CV of the residuals is, therefore, computed aswhere is the average total tree biomass (per tree).

#### 3. Results

##### 3.1. Independent Tree Component and Total Tree Models

The fit statistics and the coefficient of regression for the best tree component and total tree models are given in Table 5 for linear and nonlinear models.

All linear and nonlinear regression equations yielded satisfactory fit statistics. The linear models presented an adjusted varying from 84.24% for stem bark and crown biomass regressions to 97.61% for total tree biomass regression; the precision, as measured by the coefficient of variation (CV) of the residuals, varied from 14.29% for total tree biomass regression to 46.34% for crown biomass regression. On the other hand, the adjusted for nonlinear models varied from 84.42% for crown biomass regression to 97.60% for total tree biomass regression, and the CV of the residuals varied from 15.18% to 46.05%. For either linear or nonlinear models, the biases, as measured by the mean residual (MR), were found to be statistically not significant using Student’s* t*-test, and relatively poor fit statistics were found for stem bark and crown biomass regressions.

##### 3.2. Forcing Additivity

In the models in (1), the most frequent best linear model form is , which was found to be the best for the root system, stem wood, stem bark, and total tree biomasses. This model form was also ranked as the second model form for crown biomass. Therefore, to enforce additivity using the CON approach, this model form was generalized for all tree components and for total tree biomasses, as can be seen from (3). Tables 6 and 7 illustrate the regression coefficients and the goodness of fit statistics, respectively, for the CON method presented in (3), the SUR method in (4), and the NSUR method in (5).

###### 3.2.1. The CON Method

The results of the CON method were the same as for equation-by-equation WLS in Table 5, except that the model for crown biomass was replaced in order to have the same regressors as the remaining tree components. Better performances were found for total tree, belowground, and stem wood biomass regressions.

The graphs of the residuals against predicted values for the CON method are presented in Figure 2 and did not show any particular trend or heteroscedasticity. The cluster of points was contained in a horizontal band, showing no particular trend, with the residuals almost evenly distributed under and over the axis of abscissas, meaning that there were not obvious model defects.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

###### 3.2.2. SUR Method

As can be verified from Table 7, the adjusted varied from 52.84% for stem bark to 86.88% for total tree biomass regression, and the CVs of the residuals varied from 39.56% for total tree biomass regression to 74.38% for stem bark. All tree components and total tree models were found to be biased, and all of these models underestimated the biomass, except for the crown biomass, which was overestimated, as was observed from the mean residual (MR). Using Student’s* t*-test, these biases (MRs) were found to be statistically significant (statistically different from zero).

The biases, model defects (under and/or overestimation, heteroscedasticity), and patterns that indicated systematic discrepancies are illustrated by the graph of the residuals in Figure 3. Analyses of the residuals for SUR did not reveal heteroscedasticity but showed that the residuals were mostly agglomerated over the axis of abscissas, meaning that the designed models predicted biomass values smaller than the observed ones, underestimating the biomass (producing positive residuals). This happened to all tree components, except for the crown biomass.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

###### 3.2.3. NSUR Method

Component models (Tables 6 and 7) showed an adjusted varying from 78.12% for stem bark to 92.76% for roots. The lowest adjusted was found for the total tree model (67.76%). The CVs of the residuals varied from 21.82% for total tree to 48.25% for the stem bark model. All tree components (except the crown) and the total tree models were biased, underestimating the biomass significantly, as shown by the observation that the MRs were significantly different from zero.

Overall, the distribution of the residuals (Figure 4) was satisfactory. Minor defects were found for crown and total tree models.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

Comparing the different methods based on the relative standard errors of the expected and predicted values for total tree biomass, computed from 11 randomly selected trees from different diameter classes (Table 8), we found that the conventional method had the smallest average standard errors of the expected and predicted total tree biomass values (2.02% and 2.20%, resp.), followed by the NSUR method (3.52% and 3.68%, resp.), and lastly the SUR method (7.72% and 7.75%, resp.). These data indicated that the CON method yielded narrower confidence and prediction intervals than the NSUR and SUR methods.

#### 4. Discussion

##### 4.1. Independent Tree Component and Total Tree Models

Linear and nonlinear models were fitted for tree component and total tree biomass estimation. The difference between the performance of the selected linear and nonlinear tree component models is negligible. However, Salis et al. [43], Ter-Mikaelian and Korzukhin [44], and Schroeder et al. [45] found nonlinear models to perform better than the linear ones.

The crown models for both linear and nonlinear models were found to be less accurate and precise than the other tree component models, as evaluated by its adjusted and CVs of the residuals, suggesting more variability. According to Pardé [46], as cited by Carvalho and Parresol [14], this is because of the variability of the internal crown structure, number of branches, and variation in wood density along the branches.

##### 4.2. Additivity

Based on all of the results and analyses, the CON method was significantly superior to the SUR and NSUR methods and showed the best fit statistics for every tree component and total tree biomass models, including the largest adjusted , smallest CV of the residuals, and no significant model bias or defects. However, although the CON method was found to be statistically superior, it should be noted that it holds only under the assumption of independence among components [4], implying that the residuals are interrelated, therefore not taking into account contemporaneous correlations.

Among the methods that consider contemporaneous correlations (i.e., SUR and NSUR), NSUR appeared superior to SUR. For all tree components, NSUR was superior to SUR, with the higher adjusted , smaller CV of the residuals, and less bias. However, the total tree model of the SUR method presented a higher adjusted when compared to the total tree model of the NSUR method. Figure 5 shows that the SUR method described the data quite poorly, whereas the CON and NSUR methods described the data satisfactorily, even for the total tree model, for which the SUR method had the higher adjusted value than the NSUR method. The predicted regression lines in Figure 5 were obtained from 11 randomly selected trees from different diameter classes (2 trees per diameter class, except the last where only 1 tree was selected due to fewer representative trees); therefore the lines are function of changing all variables (refer to Table 8), hence exhibiting waves, since other variables (, CH, and LCL) did not necessarily increase as DBH increased.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

As shown in Figure 5, the regression lines for the CON and NSUR methods followed the same trend, and the CON method described the data slightly better than the NSUR method. Additionally, for all components, the SUR regression line was the poorest fit, especially for belowground, stem wood, stem bark, and total tree biomasses.

As shown in Table 7, the variance of the NSUR system was almost five times larger than that of the SUR system; however, the covariance errors for all tree components were almost two times smaller for the NSUR method; this last observation may explain the better fit of the NSUR method as all of the components in the NSUR method had larger values and smaller CVs of the residuals than those in the SUR method.

Several authors, including Parresol [4, 13], Carvalho and Parresol [14], Goicoa et al. [5], Carvalho [37], and Návar-Cháidez et al. [42], have compared different methods of enforcing additivity. All of these authors have concluded that either SUR or NSUR achieves more efficient estimates and should be the choice for additivity. However, Parresol [4] suggests that the constraint of additivity (restriction of the parameters) may compromise the efficiency of the results, a conclusion supported by SAS Institute Inc. [33], which states that restrictions should be consistent and not redundant; that is, the data must be consistent with the restriction. In fact, the lower efficiency and precision of the SUR and NSUR estimates when compared to the CON method are associated with the imposed restriction as* t*-test results based on all the restrictions imposed on SUR and NSUR were highly significant, indicating that the data were not consistent with the restriction and that the models did not fit as well with the restriction imposed. For greater details, please see Tables 9 and 10, which are SAS outputs that test the significance of the restrictions imposed for weighted SUR and weighted NSUR, respectively.

Carvalho [37] compared methods of enforcing additivity and found that the bias (MR) for stem wood was slightly larger when the models were fitted simultaneously using SUR than when tree components were fitted separately using ordinary least squares (OLS), even though other fit statistics had improved with SUR. Similar findings were observed by Goicoa et al. [5], who found that the SUR method was highly biased as it exhibited large MR and mean relative standard error (MRSE) values.

Parresol [4, 13], Carvalho and Parresol [14], and Carvalho [37] found that multivariate procedures (SUR and/or NSUR) produce more reliable estimates than when equations are estimated independently (e.g., the CON method or independent tree component models). However, Repola [47] found no significant improvements in parameter estimates using SUR when compared to the case where the models are fitted independently.

Due to the large bias found using the SUR method, this method cannot be used for biomass estimation of any tree component. However, because the bias is far smaller than that of the SUR method, the NSUR method can be used for biomass estimation as long as the bias is considered. Moreover, while the NSUR method is not superior to the CON method in terms of bias, it does have the advantage of considering contemporaneous correlation, which the CON method does not.

The CV of the residuals for total tree biomass model obtained using the HAR procedure for the linear and nonlinear models was 72.4 and 69%, respectively, 83 and 216% larger than the CV for total tree model obtained using the SUR and NSUR procedures, respectively. This shows that, in this case, the model for total biomass obtained using SUR or NSUR procedure provides more precise results than what would be obtained by summing up the individual component models to the total (HAR procedure).

##### 4.3. Extrapolation

The models fitted in this research (separately or simultaneously) are based on a dataset of 93 trees with diameters varying from 5 to 32 cm.* A. johnsonii* trees can reach diameters at breast height (DBHs) larger than 35 cm. In a forest inventory of* A. johnsonii* tree with a minimum DBH of 10 cm, Magalhães and Soto [48] (unpublished data) found only 13 trees per ha with DBHs larger than or equal to 30 cm, corresponding to only 5% of trees per ha. In this study, using 23 plots randomly distributed in the study area and a minimum DBH of 5 cm, we found only 19 trees per ha with diameters larger than or equal to 32.5 cm, corresponding to 1.54% of the total number of trees per ha. This implied that no serious bias would be added when extrapolating the models (independent tree component or NSUR models) outside the diameter range used to fit the models since very few trees were found outside the diameter range.

The models can also be safely applicable and valid over the whole range of areas where* A. johnsonii* occurs and outside the study area. This is true because the study area covered the entire range of soil and climate variations where* A. johnsonii* occurs (despite the apparent lack of large variations). For example, besides the Chibuto, Mandlakazi, Panda, Funhalouro, and Mabote districts that comprised the study area, Mecrusse (*A. johnsonii* stands) is also found in Mabalane, Massangena, and Chicualacuala districts. However, in these latter districts, the soils were nearly identical to those of the study area, composed mainly of Ferralic Arenosols [16, 22]. Similarities were also found with regard to climate and hydrology, especially with regard to rain shortage [17–21].

##### 4.4. Effect of the Measurement Procedures on the Estimates

In this study, wood density was obtained by dividing oven dry weight (at 105°C) of the discs (with and without bark) by the relevant saturated wood volume [27, 30] (air not included). It is noteworthy to mention that different definitions of the weight and volume of the discs would potentially influence the estimates of density and therefore biomass. For example, Husch et al. [28] define density as the ratio of oven dry weight and green volume (air included). Compared to the definition adopted by us, such a definition would potentially lead to large values of wood density and consequently wood biomass, as saturated volume is the maximum volume [49] and is expected to be larger than green volume.

Moura et al. [49] found no significant differences between those densities, as according to these authors, the densities must be quite the same because volume is not expected to vary above the fibre saturation point (FSP). FSP of a wood is here defined as the maximum possible amount of water that the composite polymers of the cell wall can hold at a particular temperature and pressure [50], excluding, therefore, free and adsorbed water. Differences in wood density and biomass estimates could also be found if the discs were dried to different moisture content (e.g., 12%) or if a different drying temperature was used (e.g., 65°C).

Stem was defined as the length from the top of the stump to the height corresponding to 2.5 cm diameter. Differences among stem definitions (e.g., different stump height or different minimum top diameter, stump considered as part of the stem) would affect the biomass estimates, especially stem and root biomasses. Different estimates of root biomass could also be found if the root system was partially removed, as performed by many authors (e.g., [41, 51–54]), if the depths of excavation were predefined [41, 51, 55], if fine roots were excluded [56, 57], and if root sampling procedures were applied, for example, where only a number of roots from each root system are fully excavated, and then the information from the excavated roots is used to estimate biomass for the roots not excavated [58–60].

#### 5. Conclusions

This study showed that CON method was found to be unbiased and to fit the tree component and total tree biomass well; however, the CON method had the disadvantage of not considering contemporaneous correlations. Among methods that consider contemporaneous correlations, NSUR was far superior to SUR and fit the biomass reasonably well; however, both methods were significantly biased. The CON method can be used safely as long as its limitation is considered. The NSUR method can also be used as long as the bias is accepted and taken into account. Moreover, we recommend that the SUR method should not be used due to its bias and poor description of biomass data. Since the data sets used to build the models (both independent and simultaneous) represented many variations (all diameters, soils, and climatic ranges), the selected models can be used for extrapolation.

#### Appendix

See Figure 1 and Tables 2 and 3.

The Table 3 shows the results of SUR using the system of the best linear model forms. However, as 3 of the 5 equations in the system use the same independent variables, the SUR is not effective: it has lower adjusted , sometimes negative, larger CV of the residuals, larger system variance () and larger component covariance errors () when compared to the results in Tables 6 and 7. The most jeopardized equations are those sharing the same independent variables.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This study was funded by the Swedish International Development Cooperation Agency (SIDA). Thanks are extended to Professor Almeida Sitoe for his advices during the preparation of the field work and to the field and laboratory team involved in felling, excavation, and fresh- and dry-weighing of trees and tree samples: Albino Américo Mabjaia, Salomão dos Anjos Baptista, Amélia Saraiva Monguela, João Paulino, Alzido Macamo, Jaime, Bule, Adolfo Zunguze, Viriato Chiconele, Murrombe, Paulo Goba, Dinísio Júlio, Gerente Guarinare, Sá Nogueira Lisboa, Francisco Ussivane, Nkassa Amade, Cândida Zita, and José Alfredo Amanze. The authors would also like to thank the members of local communities, community leaders, and Madeiraarte Forest Concession for the unconditional help. Acknowledges are also due to two anonymous reviewers whose insightful comments and suggestions have improved considerably this paper.