Species-specific biomass allometric models and expansion factors for indigenous and planted forests of the Mozambique highlands

Secondary Miombo woodlands and forest plantations occupy increasing areas in Mozambique, the former due to anthropogenic activities. Plantations, mainly species of Eucalyptus and Pinus, are being established on sites previously covered by secondary Miombo woodlands. This affects the evolution, cycle and spatiotemporal patterns of carbon (C) storage and stocks in forest ecosystems. The estimation of C storage, which is indispensable for formulating climate change policies on sequestrating CO2, requires tools such as biomass models and biomass conversion and expansion factors (BCEF). In Mozambique, these tools are needed for both indigenous forests and plantations. The objective of this study is to fit species-specific allometric biomass models and BCEF for exotic and indigenous tree species. To incorporate efficient inter-species variability, biomass equations were fitted using nonlinear mixed-effects models. All tree component biomass models had good predictability; however, better predictive accuracy and ability was observed for the 2-predictors biomass model with tree height as a second predictor. The majority of the variability in BCEF was explained by the variation in tree species. Miombo species had larger crown biomass per unit of stem diameter and stored larger amounts of biomass per stem volume. However, due to relatively rapid growth, larger stem diameters, heights, and stand density, the plantations stored more biomass per tree and per unit area.


Introduction
Forests sequester one third of carbon dioxide (CO 2 ) emissions from fuel and land use changes (Houghton 2007;Pan et al. 2011), store large amounts of carbon (C) in vegetation and soils, are sources of C when disturbed or destroyed by anthropogenic or natural causes, and become atmospheric C sinks during forest growth after disturbance . Land use changes are associated with ecosystem C change (Fan et al. 2016) and are the second largest cause of C losses after fossil fuel combustion (Fu et al. 2010;Lozano-García and Parras-Alcántara 2013).
Miombo is the main and most extensive forest type in southern Africa (Ribeiro et al. 2015), and in Mozambique, it comprises about two-thirds of the country´s forested land (Aquino et al. 2018). Secondary Miombo woodlands and forest plantations have been occupying increasing areas in Mozambique, the former due to anthropogenic activities. Plantations, mainly species of Eucalyptus and Pinus (Blid 2014), are being established in areas previously occupied by secondary Miombo woodlands, where shrubs and small trees predominate due to past exploitation (Magalhães 2014). This

Data acquisition
The target species with DBH ≥ 5 cm were the Miombo species B. spiciformis and J. globiflora, and plantations of E. cloeziana, P. taeda, P. elliotti, and P. patula (Table 1). A total of 120 trees, distributed among the species as shown in Table 1, were randomly selected, their diameters measured, and harvested.
Total heights (H) were recorded, and the trees divided into aboveground biomass (AGB) components, stem and crown. The stem was defined as the length of the trunk from a predefined stump height of 20 cm to a top diameter of 2.5 cm. The crown was the sum of branches, foliage, seeds, flowers, and the remaining portion of the stem, from the 2.5 cm diameter height to the tip of the tree.
The stem was dived into five equal segments and the diameter measured at the midpoint. After fresh-weighting each segment, a disc sample was removed from the top, fresh weighted, oven-dried at 105 °C to constant mass, and subsequently re-weighed. The dry mass of the whole segment was then obtained by multiplying the ratio of oven dry to fresh mass of the disc by the fresh mass of the segment. The dry mass of the stem (stem biomass), was the sum of the dry masses of the constituent segments. The volume of the stem was computed using Hohenadl's formula (Magalhães and Seifert 2015a).
For logistical reasons and to ensure the largest sample size, the foliage was not separated from the branches, as it is a time-demanding task. The crown was fresh weighed, a sample was collected, weighed in the field, and oven-dried in the laboratory. To ensure proportionality in sampling, the primary branches were divided in coarse (Ø ≥ 5 cm), and fine (Ø < 5 cm), where Ø is the diameter of the branch at the insertion point on the stem. For each size class, a primary branch with mid diameter was selected and a 3-5 cm disc removed from the insertion point. A portion of the twigs, leaves, flowers, and fruits of each primary branch was collected. The crown sample, therefore, consisted of the discs, twigs, leaves, flowers, and fruit from the selected branches. The dry mass of the crown was obtained similarly to that for each stem segment. The total aboveground biomass (AGB) was the sum of the dry masses of the stem and crown.

Data analysis
Tree component biomass models were fitted using nonlinear regression, being preferred over linear ones because biomass is a nonlinear function of stem diameter and height (Schroeder et al. 1997;Ter-Mikaelian and Korzukhin 1997;Bolte et al. 2004;Salis et al. 2006). Power or allometric functions, as an extension of nonlinear functions, were preferred over other forms of nonlinear functions because growing plants maintain the proportions between different parts, e.g., between tree biomass and predictive biometric variables such as DBH and tree height (Pilli et al. 2006). Allometric growth is observed when a constant ratio or proportion is maintained between the growth rates of different parts of the plant (Enquist 2002;Opik and Rolfe 2005). This is observed between biomass and the easily measurable variables DBH and height (Zianis and Mencuccini 2004;Fehrmann and Kleinn 2006;Pilli et al. 2006). Biomass allometric functions are represented as y = bx k , where y represents tree or component biomass, x an easily measurable tree variable, k is the ratio of growth rate between y and x, and b is the scaling coefficient. The ratio of growth rate is the scaling exponent (Zianis and Mencuccini 2004;Fehrmann and Kleinn 2006;Packard and Boardman 2008). Because both the scaling coefficient and the scaling exponent vary with species (Pilli et al. 2006) as a result of varying architecture and wood density (Ketterings et al. 2001), species-specific allometric models; are preferred over general ones. Architecture influences biomass allocation and allometry (Coll et al. 2008;Trubat et al. 2012). Therefore, biomass models were fitted using nonlinear mixed-effects (NLME) models to incorporate inter-species variability. The general form of NLME models is: where f is a nonlinear function of the parameter vector ɸ ij , a composite parameter vector incorporating both fixed and random effects, and the predictor vector x ij , M is the number of species, n i the number of trees within a species, Y ij the response vector or tree component biomass vector for the n i trees in the ith species, and ε ij the vector of multivariate normally distributed errors in species i. ɸ ij is modelled as: where β is the vector of fixed-effects parameters, δ i the vector of random-effects parameters for species i, A ij and B ij are design matrices of known constants for combining the fixed and random effects of species i.
Biomass models were fitted using allometric equations, a specific form of nonlinear equations, expressed as: Therefore, the general form of allometric mixed-effects model is expressed as: where β 0 is the fixed scaling coefficient, δ 0i the random scaling coefficient for species i, β 1 the fixed scaling exponent, δ 1i the random scaling exponent for species i, ɸ 0i the composite scaling coefficient (β 0 + δ 0i , fixed + random scaling coefficient) for species i, ɸ 1i the composite scaling exponent (β 1 + δ 1i , fixed + random scaling exponent) for species i.
The relation between response variable, the tree component biomass, and the predictors DBH and height were described by the widely used allometric biomass functions (Zianis et al. 2005;Cienciala et al. 2006;Vejpustková et al. 2015), as follows: Additionally, a model form using the root collar diameter (RCD) only as a predictor was considered to allow the estimate of harvested biomass, as after harvesting the stump dimensions are still available.
Allometric models are often fitted to logarithmically transformed data (Overman et al. 1994;Bervian et al. 2006;Cienciala et al. 2006;Blujdea et al. 2012;Makungwa et al. 2013;Paul et al. 2013aPaul et al. , b, 2019Vahedi et al. 2014;Annighöfer et al. 2016;Goussanou et al. 2016;Chen et al. 2017) to deal with heteroskedasticity, where the standard errors of a variable, over a specific time, are non-constant (Overman et al. 1994;Blujdea et al. 2012;Zapata-Cuartas et al. 2012;Vahedi et al. 2014;Goussanou et al. 2016). They facilitate model fitting by transforming the equation to a linear form (Overman et al. 1994;Blujdea et al. 2012;Zapata-Cuartas et al. 2012;Goussanou et al. 2016), allowing the use of linear least squares. In this study, logarithmic transformation of the data was avoided because the results of standard statistical tests performed on such data are often not relevant (Feng et al. 2014). Log-transformation leads to biased results (Packard and Boardman 2008). Traditional allometric method, which consists of linear fitting to logarithmically transformed data and back-transformation to power-law form, is not well suited for fitting statistical Species-specific biomass allometric models and expansion factors for indigenous and planted… models to data expressed in the arithmetic scale (Packard 2013). Castro et al. (1996) maintained that modelling raw un-transformed data gives the best results.
To address within-species heteroskedasticity, the models in Eqs. 5-8 were, instead of being log-transformed, fitted using weighted NLME models. This is justified because, quite often in regression, the error variance is functionally related to the predictors (Parresol 1999), i.e., the variability of the biomass increases with the predictors (Picard 2012).
The weight functions that describe the within-species heteroskedasticity structure were obtained by modelling the error structure of ordinary least squares (OLS), following the description by Parresol (1999Parresol ( , 2001. The squares of OLS residuals were fitted against the different combination of the predictors and it was assumed that the squares are representative of the error variance (Parresol 2001). NLME models were fitted using the NLME package (Pinheiro et al. 2019) of R software (R Core Team 2020).
Various authors have suggested that tree biomass scales against stem diameter with a universal scaling exponent of 8/3 (West et al. 1999;Enquist 2002;Pilli et al. 2006). However, several researchers have noted a lack of agreement between the scaling exponent in allometric equations and the universal scaling exponent (Chambers et al. 2001;Zianis and Mencuccini 2004;Li et al. 2005;Návar 2009). Other researchers have stressed that the universal scaling exponent is not acceptable, since the ratio of aboveground biomass to diameters in different environmental conditions is not constant (Zianis and Mencuccini 2004;Niklas 2006), and because there is a large variability in allometry across species (Zapata-Cuartas et al. 2012). Therefore, to test the allometric universal scaling law, Eq. 5 was refitted for all tree components with the scaling exponent set to 8/3. Predicted biomass values from Eq. 5 fitted without and with limitations on the scaling exponent and the observed biomass were compared using one-way ANOVA. Pairwise comparisons among biomass means were tested using Tukey HSD test.
When fitting the models, parameters, either fixed or random, with estimates not significantly different from zero at α = 0.05 were removed, and the reduced model refitted. This procedure was continued until only parameters with statistically significant estimates remained.
Tree component BCEFs were computed as the ratio of tree component biomass and stem volume (Magalhães and Mate 2018). A three-way ANOVA was carried out to test for the effects of species, DBH and total height on BCEF. Tukey HSD test was used to find species-level BCEFs statistically different from each other.
Communality analysis was carried out to quantify the variance that was unique to each predictor and the variance that was common to groups of predictors, thereby identifying which predictor accounted the most to the variability of BCEF. Commonality analysis was performed using "yhat" package (Nimon et al. 2015). The dependence of BCEF on DBH and height was analysed by testing the significance of the Pearson´s correlation coefficient.
All statistical analyses were performed at a 5% significance level.

Model evaluation and validation
Predictive accuracy determined which model(s) described the relationship between component biomass and the most used dendrometric predictors better, based on the following goodness of fit statistics: Akaike´s Information Criterion, AIC, (Akaike 1973), mean residual E , (von Gadow and Hui 1999), model precision V (Calama and Montero 2004;Correia et al. 2010a), Furnival´s index of fit, FI (Furnival 1961), and the root mean square error, RMSE (von Gadow and Hui 1999), Eqs. 9-13, respectively: where L i is the maximum likelihood of the ith model, P i the number of parameters of the ith model, Y i observed biomass, Ŷ i predicted biomass of the ith model, ε ij residual from the biomass of the jth tree of the ith model, f´(Y) is the first derivative of the transformed response variable, and MSE i the mean square error of the ith model. Square brackets indicate the geometric mean.
For ease of comparison between models, E , V, and RMSE were expressed as relative values, a percentage of observed biomass Y, which is more revealing. Model precision (V) was expressed as the standard error of the residuals.
The ideal value of E , V, and RMSE is zero (von Gadow and Hui 1999), indicating a perfect fit. Lower values of AIC and FI indicate a better fit model (Furnival 1961;Akaike 1973). Therefore, models with the lowest AIC, FI, E , V, and RMSE were judged the best.
The coefficient of determination (R 2 ) was not used to evaluate the performance of the models because it is inappropriate for demonstrating the performance or validity of nonlinear models (Spiess and Neumeyer 2010;Magalhães and Mate 2018). This is because the regression sum-ofsquares and the residual sum-of-squares do not total the sum-of-squares as in linear least squares, and thus R 2 is no longer between 0 and 100% (Magalhães and Mate 2018).
To evaluate the predictive ability, the models were validated using an independent dataset. The following statistics were used: model efficiency or EF, (Soares et al. 1995;Vanclay and Skovsgaard 1997) and Error (Calama and Montero 2004), Eqs. 14 and 15. Additionally, observed and predicted biomasses were compared using the Wilcoxon signed rank test.
EF was expressed as a percentage. The ideal values of EF and Error are 100% and 0%, respectively, indicating perfect predictive ability. Therefore, models with the largest EF, the smallest Error, and no significant difference between observed and predicted biomass, as judged by the Wilcoxon test, were considered the best.

Measured data
Component biomass plotted against DBH, stratified by species, strongly followed a power relationship for all components (Fig. 1). Tree species revealed a differing set of trees, especially for crown biomass. The power trend was strong for stems and AGB but the scatter was larger for crowns, especially for larger DBH trees. The ratio of stem biomass to AGB was larger for exotic species than for indigenous species, and amounted to 94%, 81%, 77%, and 71% for E. cloeziana, P. patula, P. taeda, and P. elliotti, respectively (Fig. 2). For the Miombo species, the largest ratio was crown biomass to AGB, and amounted to 55% and 54% for B. spiciformis and J. globiflora, respectively (Fig. 2). The share of component biomass on total aboveground biomass showed a species-specific pattern (Fig. 3). For all exotic species except for Pinus patula, the percentage of stem and crown biomass was stable over the DBH range. The percentage of stem biomass of P. patula increased slightly with diameter and a slight decrease was found for crown biomass. For the Miombo species, while the percentage of stem biomass decreased steeply with DBH, there was a steep increase for the percentage of crown biomass with DBH. For B. spiciformis, the portion of stem biomass to AGB decreased with DBH from approximately 90% to 30% and that of the crown increased from 14 to 73%. Similar results were observed for J. globiflora, with the percentage of stem biomass decreasing with DBH from 87 to 8% and the crown increasing from 13 to 92%.

Biomass allometry
Species-specific allometric biomass models with DBH, RCD (root collar diameter), and height (H) as predictors and component biomass as response variables were constructed. The weight functions and the random scaling coefficients and exponents of the models were influenced by species-specific factors and by biomass components. All fixed scaling exponents and fixed and random scaling coefficients were significant at α = 0.05 (Tables 2, 3, 4 and 5). However, the random scaling exponents were significant only for crown biomass models (Tables 2, 3, 4 and 5). All biomass allometric models had an acceptable predictive accuracy and ability (Table 6). However, stem and AGB allometric models fitted the data better than crown  (Table 6).
Using only DBH as a predictor, the models estimated tree component biomass with a satisfactory predictive accuracy, i.e., model precision (V) < 2% and mean residual (E) not statistically different from zero, and predictive ability, model efficiency (EF) > 72% and Wilcoxon test revealing that observed and predicted biomass were statistically identical. Adding height as an incorporated variable (D 2 H) or as a second variable improved the predictive accuracy and ability, however better prediction was observed when height was added as a second variable (Eq. 6, Y = 0i × D ϕ 1i × H ϕ 2i + ). Adding height as an incorporated variable (D 2 H), AIC, model precision V, and the root mean square error, RMSE, decreased up to 5%, 23%, and 26%, respectively, and the EF increased up to 15%. However, adding height as a second variable, the AIC, V, and RMSE decreased up to 10%, 28%, and 50%, respectively, and EF increased up to 28%. The allometric biomass model including height as the second predictor (Eq. 6) was the best for all components.
Satisfactory predictive accuracy and predictive ability were also attained for the models using only root collar diameter (RCD), however the models were not better as those with only DBH as a predictor (Table 6). The models fitted using the universal scaling exponent (D 8/3 ) were better than the RCD models and fitted the data similarly as the DBH models. Figure 4 shows clearly that the measured biomass, the predicted biomass from Eq. 5, and the predicted biomass using the universal scaling exponent did not differ statistically.
Using the interactive D 2 H as a predictor, an isometric scaling (b 1 = 1) was obtained for all tree components, implying that, for the DBH range considered in this study, the ratio of biomass and D 2 H was constant. A positive allometry was found when DBH or RCD were used as the sole predictor. A fixed scaling exponent (b 1 ) of ≈ 3 (range: 2.89-3.07) was found when RCD was used as the predictor. Using only DBH as a predictor, the fixed scaling exponent (b 1 ) ranged from 2.60 to 2.89. These fixed scaling exponents, obtained when biomass is scaled against DBH only, were not statistically different to the universal scaling exponent (8/3) proposed by West et al. (1999).
Crown models using only DBH as a predictor for the Miombo species had larger scaling exponents (up to 25%) than those for exotic species. The larger values for the indigenous species imply that, per unit of DBH growth, the increment of crown biomass is larger, up to 25%, compared with exotic species. Note that, for stem biomass and AGB models, the random scaling exponents were not statistically different from zero, thus the composite scaling exponents were equal to the fixed ones, denoting that a unit of DBH increase resulted in similar stem and aboveground biomass increment for all species. where δ 0i is the random scaling coefficient of the species i, (e.g., δ 0 B. spiciformis random scaling coefficient for Brachystegia spiciformis), s(δ 0 ) standard deviation of the random scaling coefficients, s(e) standard deviation of the residuals, ɸ 0i composite scaling coefficient of the species i. Non-significant random parameters are not presented

Biomass conversion and expansion factors (BCEF)
Biomass conversion and expansion factors (BCEF) are widely used to convert timber volume of forest inventories to tree biomass, and are calculated as the ratio of component biomass to stem volume. Aboveground biomass of B. spiciformis and J. globiflora were approximately 157% and 139% of the stem volume (Table 7), respectively; however, the AGB of the exotic species were only up to 69% of stem volume. This suggests that Miombo species store more than double the aboveground biomass per stem volume, i.e., the BCEF are more than double that of the exotic species. Miombo species store up to 13 times more biomass in the crown per stem volume than the introduced species.
The three-way ANOVA revealed that, for all three components, stem, crown, and shoot system, the species was a significant source of variation of the BCEF (P value = 0.0000). In addition to species, DBH and height also had significant effects on crown BCEF; height was also a significant source of aboveground BCEF variation. Tukey HSD test showed that BCEF values of the two Miombo species were statistically superior to those of all introduced species (Table 7).
From the commonality analysis, the majority (60-99%) of the regression effect was explained by variance that was unique to a species, i.e., species uniquely accounted for up to 99% of the regression effect. Tree species accounted for 59.5%, 62.5%, and 99% of the regression effect of crown, aboveground, and stem BCEF, respectively. where δ 1i and δ 2i are the random scaling exponents of the species i, (e.g., δ 1 B. Spiciformis random scaling exponent for Brachystegia spiciformis), s (δ 1 ) standard deviation of the random scaling exponents δ 1 , and s (δ 2 ) standard deviation of the random scaling exponents δ 2 . Non-significant random parameters are not presented

Parameter
Equation 5 Equation 6 Equation 7 Equation 8 Weight Tree component BCEF of all exotic or introduced species, and stem and aboveground BCEF of Miombo species were not linearly correlated to DBH nor to height (P value > 0.2). Crown BCEF of Miombo species was positively dependent on both DBH and height (Fig. 5), and Pearson´s correlation coefficients between crown BCEF and DBH and height were 0.55 (P value = 0.0002), and 0.33 (P value = 0.0349), respectively.

Biomass allocation patterns
Studies on biomass allocation patterns into different tree components are of considerable interest in forest carbon monitoring and dynamics (Fournier et al. 2003;Pajtík et al. 2008;Magalhães and Seifert 2015b). For example, with the Miombo woodlands, while the stem or bole is used in the forest industry and thus the C remains stored for some time, the branches are either left in the forest or collected and used as firewood by local communities. In either case, the stored carbon is released to the atmosphere by decomposition or burning.
The biomass allocation patterns of the exotic species in this study is in agreement with a number of studies (Cienciala et al. 2006;Correia et al. 2010a, b;Xiang et al. 2011;Vargas-Larreta et al. 2017) for Pine and Eucalyptus species; as trees grow larger, the proportion of aboveground biomass allocated to the stem increases and that allocated to the crown decreases, and/or the majority of the AGB is allocated to the stem. However, this pattern of allocation was not supported by the two Miombo species as most of the AGB was allocated to the crown, and the proportion increased with increasing tree size. An increase proportion of branch biomass, and thus crown biomass, has been reported by Pajtík et al. (2008) for Norway spruce (Picea abies (L.) H. Karst.). The proportion of AGB allocated to the stem for the Miombo species in this study, 45% and 46% for B. spiciformis and J. globiflora, respectively, is lower Table 6 Predictive accuracy and ability of the fitted models Equation 5 (Universal) = Eq. 5 fitted under the constraint of the scaling exponent being equal to the universal scaling of 8/3; i.e. assuming that tree biomass scales against stem diameter with a universal scaling exponent of 8/3 Parameter Species-specific models

Equation 5
Equation 6 Equation 7 Equation 8 Equation 5 (universal)  The different biomass allocation patterns of indigenous and exotic species could be a result of differences in stand density. The Miombo forest is composed of small trees (Table 1) and no more than 200 trees per hectare, while the plantations have upwards of 1000 trees per hectare. Xiang et al. (2011) maintained that stand density enhances the proportion of stem biomass but decreases biomass of branches. As a consequence of competition for light, exotic species have small branches and smaller quantities of leaves, resulting in smaller crowns than indigenous species which do not face limitations in light resources and growing space. Crown size has been reported to decrease with competition (van Laar and Akça 2007) as a result of natural pruning as a consequence of inter-tree competition for light and growing space and self-shading and shading by neighbouring trees.

Biomass allometry
Miombo species are the most important trees of southern Africa and of Mozambique (Ribeiro et al. 2015;Aquino et al. 2018), and Eucalyptus and Pinus species the most important reforestation species, planted mainly for commercial purposes (Blid 2014). Therefore, the accurate estimation of biomass of the Miombo species and the Eucalyptus and Pinus plantations is essential. The species-specific allometric models developed here predict biomass accurately for the exotic and indigenous species across the range of DBHs used to develop the models (Table 1) and have good applicability as tested using the independent samples. Nevertheless, it is not advisable to extrapolate beyond the data range as the results may be biased.
Crown models for the Miombo species had larger composite scaling exponents (up to 25%) compared with those of the exotics, and composite scaling coefficients were constant for all species as a result of insignificant random scaling coefficients. As a result, the Miombo species had high predicted biomass for a given diameter and height compared with the exotic species. This is consistent with the finding that: (1) Miombo species store more biomass per stem volume than exotic species; (2) Miombo species allocate more biomass to crowns than to stems as compared to exotic species (Fig. 2); and, (3) per unit of DBH growth, the increment of crown biomass of Miombo is larger than that of exotic species. However, although Miombo species had larger increments of crown biomass per unit stem diameter growth and stored larger amounts of biomass per stem volume than the exotic species, due to their rapid growth, larger stem diameters and heights (on average), and stand density, the plantations store more biomass per tree and per unit area.
For all models, AGB and stem biomass had higher predictabilities than crown biomass. This is in agreement with Jagodziński et al. (2018a, b). Branch and crown biomass components are difficult to estimate with the same accuracy as stem biomass or total aboveground biomass (Vargas-Larreta et al. 2017). The high predictability of stem biomass and AGB models is attributed to the high correlation between these component biomasses and DBH (Fig. 1).
Diameter at breast height was the main predictor for estimating tree biomass components for all species, indigenous and exotics. For all components, the inclusion of height significantly improved the predictability of aboveground biomass and stem and crown biomass. Height as a second predictor substantially improved model precision (V) of stem biomass, AGB and crown biomass by 17%, 26%, and 28%, respectively, and model efficiency (EF) was improved by 23%, 23%, and 28%, respectively. The inclusion of height led to a substantial increase in predictive accuracy and ability. This is consistent with the findings by Pajtík et al. (2008), Xiang et al. (2011) andVargas-Larreta et al. (2017). However, other researchers (Johansson 1999;Porté et al. 2002;Jenkins et al. 2003) reported that the inclusion of height does not lead to a significant increase in predictive ability.
Height is not directly measured in forest inventories and it is thus more susceptible to measuring error than DBH (Loetsch et al. 1973;Machado and Figueiredo Filho 2006;Sanquetta et al. 2006) and its measurement is not standardized (Sileshi 2014). Therefore, although with better predictability, models using height (Eqs. 6, 7) will lead to biased results because of inherent measurement errors of tree height. The inclusion of height in Eq. 4, either as a combined variable D 2 H (Eq. 7) or as a second predictor (Eq. 6), will introduce two sources of errors over and above those expected in Eq. 4 (Sileshi 2014): (1) errors due to measurement of height; and, (2) errors due to the estimation of model parameters. Moreover, 1-predictor biomass models, (DBH only models), are easier and faster to use and less expensive than 2-predictors models since the latter require height measurements which are expensive and time-consuming. In biomass modelling, height is included in the model, either as a combined variable D 2 H (Eq. 7) or as a second variable, in addition to DBH (Eq. 6) (Zianis et al. 2005;Cienciala et al. 2006;Vejpustková et al. 2015). In this study, better prediction and accuracy were obtained when height Species-specific biomass allometric models and expansion factors for indigenous and planted… was added as a second variable (Eq. 6). This contradicts Carvalho and Parresol (2003) and Bi et al. (2004) who obtained better estimates for the model with a combined predictor D 2 H. However, our results are consistent with Vahedi et al. (2014) and Vejpustková et al. (2015).

Biomass conversion and expansion factors
Larger BCEF values were found for Miombo species and lower ones for the plantation species. Slow-growing species such as B. spiciformis and J. globiflora have high wood density, whereas fast-growing species (e.g., Eucalyptus and Pinus species) have low densities (Ouédraogo et al. 2013;Yeboah et al. 2014;Ramananantoandro et al. 2016). This explains why slow-growing Miombo species (Grundy 1995;Elifuraha et al. 2008;Chiteculo and Surovy 2018) had higher biomass conversion and expansion factors compared to exotic species as these factors are directly proportional to wood density (Pajtík et al. 2008;Schepaschenko et al. 2018). In other words, BCEF = BEF × ρ, where BEF and ρ are biomass expansion factors and wood density, respectively.
BCEFs have been reported to be reversely dependent on DBH and height (Brown et al. 1989;Lehtonen et al. 2004;Dutca et al. 2010;Sanquetta et al. 2011;Magalhães and Seifert 2015a) or to be independent of tree size (Marková and Pokorný 2011;Magalhães and Seifert 2015a). In this study, BCEF values of exotic species were not dependent on size; only crown BCEF of Miombo species were dependent on DBH and height. This is because the proportion of crown biomass on AGB increased with tree size, whereas that of stem biomass decreased (Fig. 3), resulting in BEF, and consequently BCEF, to increase with increasing tree size. At the same time, the stability of the biomass conversion and expansion factors of exotic species over tree size is also related to the stability of the proportions of crown and stem biomass on AGB over tree size (Fig. 3). Dutca et al. (2010) maintained that the reverse relationship between BEF and tree size is a result of an inverse relationship between wood density and size. This suggests that the direct dependence of crown BCEF of Miombo species to tree size may be due to the direct relationship between branch wood density and size. Increasing wood density with age and tree size has been reported by various researchers (Pajtík et al. 2008;Nock et al. 2009;Henry et al. 2010;Deng et al. 2014). However, there are also studies showing that wood density does not vary with stem diameter (Ramananantoandro et al. 2016;Fajardo 2018), whereas others have reported that wood density decreases with tree age and size (Pajtík et al. 2008;Liepiņš et al. 2017). This pattern of variability in wood density explains the variability of patterns of BCEF with tree size (Fig. 3).
The majority of variation of the biomass conversion and expansion factors was attributed uniquely to species rather than to stem diameter and height, and this is consistent with Luo et al. (2014) who found that BCEF exhibited remarkable variation across forest types.

Conclusions
This study provided species-species biomass allometric models and biomass conversion and expansion factors for the Miombo species B. spiciformis and J. globiflora and for exotic species of Eucalyptus cloeziana, Pinus elliotti, Pinus patula, and Pinus taeda growing in highlands of Manica province, Mozambique. The allometric biomass models, for all tree components, had an acceptable predictive accuracy and ability. The highest predictive accuracy and ability was found with models with tree height as a second variable (Y = 0i × D ϕ 1i × H ϕ 2i + ) . The Miombo species had larger increments of crown biomass per unit of stem diameter as seen by their larger scaling exponents. Similarly, they stored larger amounts of biomass per stem volume unit, i.e., their biomass conversion and expansion factors were up to 13-fold larger than those of the exotic species. The majority of the variation of these factors (up to 99%) was attributed uniquely to species rather than to diameter or height. Although with relatively lower, yet acceptable predictive accuracy and ability, the model fitted under the constraint of the allometric universal scaling law (Y = ɸ 0i × D 8/3 + ε) showed that the universal exponent (8/3) was data acceptable for both indigenous and exotic tree species. the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.