ABSTRACT
The objectives of this work were to propose additive equations for biomass components (stem and crown) and total biomass for black wattle (Acacia mearnsii De Wild.) and show the inconsistency of independently adjusted biomass equations. Two procedures were used to fit nonlinear equations of biomass: i) independent and ii) systems of equations. The second procedure, defined by the application of the seemingly unrelated regression model, has better biological and statistical properties to estimate allometric equations of biomass components and total biomass when compared with the independent estimation. An effective property of this procedure is the additivity, i.e., the estimates of component biomass are compatible with those of total biomass. Independent fitted adjusted equations do not consider the dependence between the biomass components, thus, besides the estimates being nonadditive, which is an undesirable property, they will result in estimates with larger variance.
Keywords:
nonlinear seemingly unrelated regression; error modeling; additivity
RESUMO
Os objetivos desse trabalho foram propor equações aditivas de biomassa dos componentes (fuste e copa) com a biomassa total para a espécie acácia negra (Acacia mearnsii De Wild.) e demonstrar a inconsistência de equações de biomassa ajustadas independentemente. Dois procedimentos foram utilizados para ajustar equações não lineares de biomassa: i) independente e ii) sistemas de equações. O segundo procedimento, definido pela aplicação do modelo de regressão aparentemente não relacionada, apresenta melhores propriedades biológicas e estatísticas para estimar equações alométricas de biomassa dos componentes e biomassa total, quando comparado com a estimação independente. Uma propriedade efetiva desse procedimento é a aditividade, isto é, as estimativas de biomassa dos componentes são compatíveis com as de biomassa total. As equações ajustadas independentes não consideram a dependência entre os componentes de biomassa, assim, além das estimativas não serem aditivas, propriedade indesejável, resultarão em estimativas com maior variância.
Palavraschave:
regressão não linear aparentemente não relacionada; modelagem do erro; aditividade
1. INTRODUCTION
Energy generation from forest biomass has the potential to reduce carbon emissions when compared with other energy sources, especially fossil fuels. Therefore, accurate biomass estimates have been fundamental for diverse applications. Information on biomass stocks is especially useful for forest managers to determine the role of forest biomass in the global carbon cycle and to manage energy production from biomass. Hence the importance of establishing equations and enabling forest biomass inventories.
The development and testing of allometric biomass models depend on the availability of data from destructive tree evaluations and is a timeconsuming and expensive task ( Chave et al., 2014 Chave J, RéjouMéchain M, Búrquez A, Chidumayo E, Colgan MS, Delitti WBC et al. Improved allometric models to estimate the aboveground biomass of tropical trees. Global Change Biology 2014; 20(10): 31773190. http://dx.doi.org/10.1111/gcb.12629. PMid:24817483.
http://dx.doi.org/10.1111/gcb.12629 ...
). To avoid this work being carried out frequently, it is possible to use existing biomass equations and to estimate biomass stocks using forest inventory information.
These equations must be established by asserting some assumptions. Genet et al. (2011) Genet A, Wernsdörfer H, Jonard M, Pretzsch H, Rauch M, Ponette Q et al. Ontogeny partly explains the apparent heterogeneity of published biomass equations for Fagus sylvatica in central Europe. Forest Ecology and Management 2011; 261(7): 11881202. http://dx.doi.org/10.1016/j.foreco.2010.12.034.
http://dx.doi.org/10.1016/j.foreco.2010...
pointed out that a good set of biomass equations must meet the following characteristics: i) consistency, to ensure the additivity of the estimates between biomass components and total biomass; ii) robustness, to compose an operating system that works correctly for wide variations in the sample population and with low sensitivity to the sampling process and to the formulated hypotheses; and iii) accuracy.
Consistency can be ensured by applying the seemingly unrelated regression for linear models and its variation for nonlinear models, according to Parresol (1999 Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. Forest Science 1999; 45: 573593. , ^{2001} Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
). Additive biomass equations are desirable because the estimates of the component biomass and total biomass are harmonic, overcoming the lack of consistency in the estimates.
Black wattle (Acacia mearnsii De Wild.) is one of the most grown forest trees in the state of Rio Grande do Sul. It is a versatile species with characteristics allowing the generation of economic, social, and, above all, environmental benefits. Growing black wattle forests and trading their products is a solid economic activity and, according to Stein & Tonietto (1997) Stein PP, Tonietto L. Black wattle silviculture in Brazil. In: Brown AG, Ko HC, editors. Black wattle and its utilization. Barton: Rural Industries Research and Development Corporation; 1997. , has brought benefits and prosperity to more than 40 municipalities and approximately 40 thousand families.
Equations for estimating black wattle biomass have already been proposed by several researchers, as we highlight throughout this text. In the present work a new approach is applied to fit the equations, using data from the main forest sites in the state of Rio Grande do Sul, from stands with different ages, and an expressive sample size in terms of number of trees. The new approach used to adjust equations was based on the following hypothesis: given that the estimators of equation systems include constraints to obtain the coefficients in their adjustments, their estimates might result in additivity of the components, when compared with the estimators obtained in a single equation, and therefore result in estimates of biomass components additive with estimates of total biomass.
Therefore, the objectives of this study were to propose additive equations for biomass components (stem and crown) and total biomass for black wattle and demonstrate the inconsistency of independently adjusted biomass equations.
2. MATERIAL AND METHODS
2.1. Sampling in stands of black wattle
The field research was carried out in commercial stands of black wattle (Acacia mearnsii de Wild.) in areas where the species is prevalent in the state of Rio Grande do Sul, in June and July 2014. It was possible to divide these areas into three regions: Cristal, Encruzilhada do Sul, and Piratini. These regions have their climate, relief type, and soil characterized in Mochiutti (2007) Mochiutti S. Produtividade e sustentabilidade de plantações de acácianegra (Acacia mearnsii De Wild.) no Rio Grande do Sul [tese]. Curitiba: Universidade Federal do Paraná; 2007. .
The stands were sampled in a sequence of ages after planting, to cover all crop rotation ( Table 1 ). Four circular plots with 10 m in diameter (78.54 m^{2}) were randomly allocated in each stand. All trees in the plots were felled and measured for the following characteristics: diameter at breast height, total height, stem and crown biomass. The research involved the allocation of 48 temporary plots and the measurement of the characteristics on 670 trees.
Group and systems of equations organized by age and maturity of the stands, ranges of independent variables (
The diameter at breast height (d) was measured with a dendrometric tape and the total height (h) with a tape. Biomass measurement was performed for the stem component (stem wood + bark) and crown (live and dead branches, leaves, flowers and fruits), as defined in Picard et al. (2012) Picard N, SaintAndré L, Henry M. Manual for building tree volume and biomass allometric equations: from field measurement to prediction. Rome: Food and Agricultural Organization of the Unites Nations, Centre de Coopération Internationale en Recherche Agronomique pour le Développement; 2012. . For each tree, these components were separated and weighed to obtain the wet biomass with a Portable Scale with 5g of accuracy. Total biomass was defined as the sum of the stem biomass and crown components.
Samples of the components were taken and immediately weighted with a digital scale (Hoyle) with 1g of accuracy. Samples of the crown, approximately 1,500 g, were taken at positions 0%, 25%, 50%, 75%, and 95% of the total crown length (distance of the first branch, regardless of being alive or dead, to the apex of this component). In each stem, 5 discs of 2 cm thick were removed at the positions: 0%, 25%, 50%, 75%, and 95% of the total height. The samples were dried in an oven with aircirculation and air renewed at 100° C to constant mass and weighed with a digital scale with 1g of accuracy. Dry biomass was calculated from the values of biomass moisture content.
2.2. Aboveground biomass equations
The equations for biomass components and total aboveground biomass were proposed according to the age of the stands, as shown in Table 1 . Equations estimated independently were named as equation group, and those estimated jointly as equation system.
2.2.1. Independently estimated biomass equations: procedure 1
The first step in the independent estimation of equations was the evaluation of the competing models, described by Behling (2016) Behling A. Modelagem da biomassa de árvores para assegurar aditividade dos seus componentes [tese]. Curitiba: Universidade Federal do Paraná; 2016. for the input variables d (in cm) and h (in m) for crown biomass (
The linear models were adjusted by ordinary least squares, and the nonlinear models by estimated generalized nonlinear least squares, with estimators defined by Greene (2008) Greene WH. Econometric analysis. Upper Saddle River: PrenticeHall; 2008. . The performance of each model was evaluated based on the adjusted coefficient of determination
Group of equations
Models
1, 2 and 3
4 and 5
White’s test (1980) was applied to test the hypothesis of homogeneity of residuals in each equation, with 95% probability. When the hypothesis was rejected, the weights were obtained through the variance structure, as defined by Harvey (1976) Harvey AC. Estimating regression models with multiplicative heteroscedasticity. Econometrica 1976; 44(3): 461465. http://dx.doi.org/10.2307/1913974.
http://dx.doi.org/10.2307/1913974 ...
and Parresol (2001) Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
, followed by a new equation fitting and White’s test application. The variancecovariance matrices of
http://dx.doi.org/10.1139/x00202 ...
.
2.2.2. Biomass equation systems: procedure 2
The biomass equation systems were developed for total biomass, crown biomass and stem components. The models selected for the components were those defined in [1] and [2]. The model for total biomass was defined as a function of the independent variables of each function of component i. Thus, the equation systems were defined according to Equations 3 and 4 .
We used the model of nonlinear seemingly unrelated regressions to solve the systems of Equations 3 and 4 , with parameter estimates made in SAS software. The White’s test was applied to test the hypothesis of homogeneity of residuals for the component equations and the equation of the total, at 95% probability. For the cases in which the hypothesis was rejected, the weights obtained in the previous section were used and the system was adjusted again using the estimator of weightednonlinear seemingly unrelated regressions (WNSUR), followed by the application of the White’s test. Note that the residual scatter plots were not included, since one of the objectives of the residuals analysis is already met by presenting the White’s test and its interpretation. In addition, given the appropriate conditions to perform the adjustment (homogeneity of residual variances), the authors believe that it is more important to present the variancecovariance matrices, since they are useful to make several inferences, which will be shown in the following sections.
The performance of each model was evaluated based on the adjusted coefficient of determination
http://dx.doi.org/10.1139/x00202 ...
. These matrices were presented because they show the correlations between the components measured in the same tree, i.e., the dependence between the biomass components. These correlations will determine the gains in efficiency owing to the application of WNSUR. In addition, these matrices are needed for applications combined with the fitted equations in order to make inferences about the confidence intervals of the estimates.
3. RESULTS
3.1. Independently estimated biomass equations: procedure 1
The null hypothesis of homogeneity of the residuals was rejected for most equations, according to the White’s test. Thus, we obtained the weights for the biomass of the components and total biomass using the structures of variance in Equations 5 to 19 . For the equations with null hypothesis accepted, the weights were considered equal to one.
Group 1
Group 2
Group 3
Group 4
Group 5
The statistics of the equations adjusted by procedure 1 and complemented with the weighting are shown in Equations 20 to 34 . The correction of heteroscedasticity was evidenced by the weighting via White’s test and residual graphic analysis. In addition, all coefficients were significant at 95% probability by application of ttest.
Group 1
Group 2
Group 3
Group 4
Group 5
The variancecovariance matrices of
3.2. Biomass equation systems: procedure 2
Heteroscedasticity was confirmed for the same previous cases by the rejection of the null hypothesis via White’s test. Therefore, the new fitting of the functions was carried out using procedure 2, with the application of the weights already obtained in procedure 1. The new equations and the resulted statistics are shown in Equations 35 to 49 . The correction of heteroscedasticity was evidenced by the weighting via White’s test and residual graphic analysis. All adjusted coefficients were significant at 95% probability by application of ttest.
System 1
System 2
System 3
System 4
System 5
The variancecovariance matrices of
4. DISCUSSION
4.1. Nonlinear biomass equations
The biomass components and total biomass were modeled by nonlinear biomass functions. The adjustment of these functions requires more demanding computational methods, although their implementation is made possible by the low operational cost currently available in computers and softwares. Therefore, the use of nonlinear equations is not a limiting factor and provides certain positive conditions that are not achieved through linear regression, although there must be attention to ensure if appropriate functions are obtained.
In this work, biomass was modeled in the original (nonlinear) form for the following reasons: the relationship between dependent and independent variables are nonlinear; model linearization causes an inherent bias in the transformation of data in relation to the original arithmetic scale, thus correction is necessary; the application of nonlinear models allows the use of the original values not implying in data transformation as in the cases of model linearization; nonlinear models are better than linear models.
Although certain positive aspects of the nonlinear approach are highlighted, it also provides some challenges in relation to the linear approach. It is often necessary to work on the convergence criteria to avoid premature stopping of the iterations and a subsequent loss of precision of the resulted coefficients. This requires a more experienced modeler, compared to the use of linear models. In the nonlinear approach, the first challenge refers to the estimates of the initial parameters that must be specified, and experience is required to define them or, where possible, by means of a consistent β vector. In this condition, the functions must be linearized and adjusted by ordinary least squares to obtain the initial values, as performed in this work. The second challenge refers to which algorithm to use in the iteration process, such as GaussNewton and Marquardt.
Failures in convergence processes reveal other challenges in the nonlinear approach. According to Parresol (2001) Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
, they can occur for several reasons and the most important ones are: partial derivative matrix and direction matrix can be singular, indicating that the model can be overparameterized; it is possible for the parameters to occur in a numerical space in which the logarithmic and square root arguments are not valid, resulting in excess of calculations. Logically, these problems can be avoided by the careful choice of models and the establishment of different initial β vectors.
There are still other difficulties, that is, the iteration method can lead to steps that do not improve the estimates. In this case, according to Parresol (2001) Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
, the length of step
4.2. Modeling of variance structure
Heteroscedasticity of residuals is common in the modeling of tree biomass, and it has been reported by several researchers: Parresol (1999 Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. Forest Science 1999; 45: 573593. , ^{2001} Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
), SaintAndré et al. (2005) SaintAndré L, M’Bou AT, Mabiala A, Mouvondy W, Jourdan C, Roupsard O et al. Agerelated equations for above – and below: ground biomass of a Eucalyptus hybrid in Congo. Forest Ecology and Management 2005; 205(13): 199214. http://dx.doi.org/10.1016/j.foreco.2004.10.006.
http://dx.doi.org/10.1016/j.foreco.2004...
, Basuki et al. (2009) Basuki TM, Van Laake PE, Skidmore AK, Hussin YA. Allometric equations for estimating the aboveground biomass in tropical lowland Dipterocarp forest. Forest Ecology and Management 2009; 257(8): 16841694. http://dx.doi.org/10.1016/j.foreco.2009.01.027.
http://dx.doi.org/10.1016/j.foreco.2009...
, Bi et al. (2010) Bi H, Long Y, Turner J, Lei Y, Snowdon P, Li Y et al. Additive prediction of aboveground biomass for Pinus radiata (D. Don) plantations. Forest Ecology and Management 2010; 259(12): 23012314. http://dx.doi.org/10.1016/j.foreco.2010.03.003.
http://dx.doi.org/10.1016/j.foreco.2010...
, Sileshi (2014) Sileshi GW. A critical review of forest biomass estimation models, common mistakes and corrective measures. Forest Ecology and Management 2014; 329: 237254. http://dx.doi.org/10.1016/j.foreco.2014.06.026.
http://dx.doi.org/10.1016/j.foreco.2014...
, Wayson et al. (2015) Wayson CA, Johnson KD, Cole JA, Olguín MI, Carrillo OI, Birdsey RA. Estimating uncertainty of allometric biomass equations with incomplete fit error information using a pseudodata approach: methods. Annals of Forest Science 2015; 72(6): 825834. http://dx.doi.org/10.1007/s1359501404367.
http://dx.doi.org/10.1007/s1359501404...
, Sanquetta et al. (2015b) Sanquetta CR, Wojciechowski J, Dalla Corte AP, Behling A, Péllico S No, Rodrigues AL et al. Comparison of data mining and allometric model in estimation of tree biomass. BMC Bioinformatics 2015b; 16(247): 19. PMid:26250142. , and Zhao et al. (2015) Zhao D, Kane M, Markewitz D, Teskey R, Clutter M. Additive tree biomass equations for midrotation loblolly pine plantations. Forest Science 2015; 61(4): 613623. http://dx.doi.org/10.5849/forsci.14193.
http://dx.doi.org/10.5849/forsci.14193...
. Cunia & Briggs (1984) Cunia T, Briggs RD. Forcing additivity of biomass tables: some empirical results. Canadian Journal of Forest Research 1984; 14(3): 376384. http://dx.doi.org/10.1139/x84067.
http://dx.doi.org/10.1139/x84067 ...
pointed out that it is natural to expect that the biomass variability of trees with smaller diameters is smaller than the biomass variability of trees with larger diameters, leading to the variance of residuals to be nonconstant for the biomass equations.
This should not be ignored when fitting equations by regression analysis. Constant variance is fundamental in regression analysis, since it directly affects the validation of several hypotheses tested, as demonstrated in Maddala (2001) Maddala GS. Introduction to econometrics. 3rd ed. New York: John Wiley & Sons; 2001. and Greene (2008) Greene WH. Econometric analysis. Upper Saddle River: PrenticeHall; 2008. . Under heteroscedasticity of residuals two approaches can be applied to establish constant variance. The first one refers to the application of logarithmic transformation, and the second one of regression weighting. In regression weighting, instead of using simple weights such as
http://dx.doi.org/10.1139/x00202 ...
).
The logarithmic transformation to correct heteroscedasticity is an approach easier and faster than applying regression weight, because it does not require the modeling of the variance structure. However, it reverts to problems already evidenced, caused by data transformation. For these reasons, it is logical to model the variance structure instead of applying logarithmic or even another transformation. Parresol (2001) Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
also recommended modeling the error variance structure rather than using the logarithmic transformation to correct the heteroscedasticity of biomass data. He demonstrated that transformed or weighted equations for stem and bark biomass resulted in similar standard errors of the coefficients and in the statistics used for evaluation of goodnessoffit. However, for crown biomass and total biomass, the weighted equations showed better results for standard error of the coefficients and prediction intervals of the estimates. The transformed equations resulted in slightly more elastic prediction intervals, but this would have a greater impact if considered the cumulative effect by applying the equations to hundreds or thousands of trees.
Using the variance structure, we established the weights and readjusted the equations. This led to the stabilization of variance, hence the weighted functions allowed the homoscedasticity of residuals. Homoscedasticity is indispensable to obtain appropriate confidence intervals and to validate the hypothesis tests. An effect of achieving constant variance is observed when the standard errors of the coefficients are smaller than those obtained under heteroscedasticity of the residuals (although these results were not presented in this manuscript), which was also observed by Parresol (1999 Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. Forest Science 1999; 45: 573593. , ^{2001} Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
).
4.3. Aboveground tree biomass equations for black wattle
The coefficients of variation for the equations of the crown component were greater than those evaluated for the equations of stem and total components, because the crown biomass (composed by biomass of leaves and branches) presents more variability on allometric relationship with the independent variables. Poorter et al. (2012) Poorter H, Niklas KJ, Reich PB, Oleksyn J, Poot P, Mommer L. Biomass allocation to leaves, stems and roots: metaanalyses of interspecific variation and environmental control. The New Phytologist 2012; 193(1): 3050. http://dx.doi.org/10.1111/j.14698137.2011.03952.x. PMid:22085245.
http://dx.doi.org/10.1111/j.14698137.2...
reported an interesting discussion about biomass allocation to leaves, stems, and roots. The authors discussed the pattern of variation of these components and how their growth is influenced by the environment, plant size, dynamic process, and competition. It is known that the evolutionary patterns of plants are influenced by environmental conditions and, therefore, this reflects on the proportions among tree components, as well as in their allometry, i.e., the allometric relationship of the crown and stem components with the independent variables should not be the same, because their biomass proportions are different when equated to the same independent variable.
Adjusted coefficients of determination and variation from the fitted equations through procedures 1 and 2 are within the range obtained in other black wattle studies. Coefficients of determination greater than 80% for black wattle biomass equations have been frequently reported in the literature, with approaches developed by several researchers: Dudley & Fownes (1992) Dudley NS, Fownes JH. Preliminary biomass equations for eight species of fastgrowing tropical trees. Journal of Tropical Forest Science 1992; 5(1): 6873. , Carbonera Pereira et al. (1997) Carbonera Pereira J, Schumacher MV, Hoppe JM, Caldeira MVW, Santos EM. Produção de biomassa em um povoamento de Acacia mearnsii De Wild. no Estado do Rio Grande do Sul. Revista Árvore 1997; 21(4): 521526. , Caldeira (1998) Caldeira MVW. Quantificação da biomassa e do conteúdo de nutrientes em diferentes procedências de Acácianegra (Acacia mearnsii De Wild) [dissertação]. Santa Maria: Universidade Federal de Santa Maria; 1998. , Barichello et al. (2005) Barichello LR, Schumacher MV, Vogel MLM. Quantificação da biomassa de Acacia mearnsii De Wild na região sul do Brasil. Ciência Florestal 2005; 15(2): 129135. http://dx.doi.org/10.5902/198050981830.
http://dx.doi.org/10.5902/198050981830 ...
, Saidelles (2005) Saidelles FLF. Determinação da biomassa e altura de amostragem para a quantificação de nutrientes em Acacia mearnsii De Wild [tese]. Santa Maria: Universidade Federal de Santa Maria; 2005. , Mochiutti (2007) Mochiutti S. Produtividade e sustentabilidade de plantações de acácianegra (Acacia mearnsii De Wild.) no Rio Grande do Sul [tese]. Curitiba: Universidade Federal do Paraná; 2007. , Caldeira et al. (2011) Caldeira MVW, Saidelles FLF, Schumacher MV, Godinho TO. Biomassa de povoamento de Acacia mearnsii De Wild., Rio Grande do Sul, Brasil. Scientia Forestalis 2011; 39(90): 133141. , and Sanquetta et al. (2014) Sanquetta CR, Behling A, Corte AP, Simon A, Pscheidt H, Ruzza MS et al. Estoques de biomassa e carbono em povoamentos de acácia negra em diferentes idades no Rio Grande do Sul. Scientia Forestalis 2014; 42(103): 361370. .
The coefficients estimated from equations of biomass components via procedures 1 and 2 led to close relationship with the groups or systems of equations, determined according to the stand age. This finding was also evidenced by Williams et al. (2003) Williams CJ, LePage BA, Vann DR, Tange T, Ikeda H, Ando M et al. Structure, allometry, and biomass of plantations Metasequoia glyptostroboides in Japan. Forest Ecology and Management 2003; 180(13): 287301. http://dx.doi.org/10.1016/S03781127(02)005674.
http://dx.doi.org/10.1016/S03781127(02...
, who observed that the allometry of the components was not the same for populations of the same species with different ages. The authors pointed out that stand age influences the relationship between biomass and the variables that express tree dimensions such as d and h, as discussed by Niklas (1995) Niklas KJ. Sizedependent allometry of tree height, diameter and trunktaper. Annals of Botany 1995; 75(3): 217227. http://dx.doi.org/10.1006/anbo.1995.1015.
http://dx.doi.org/10.1006/anbo.1995.101...
.
SaintAndré et al. (2005) SaintAndré L, M’Bou AT, Mabiala A, Mouvondy W, Jourdan C, Roupsard O et al. Agerelated equations for above – and below: ground biomass of a Eucalyptus hybrid in Congo. Forest Ecology and Management 2005; 205(13): 199214. http://dx.doi.org/10.1016/j.foreco.2004.10.006.
http://dx.doi.org/10.1016/j.foreco.2004...
also found that the parameters estimated from the biomass models varied clearly with the stand age. They argued that this is not only a result of a change in tree maturity, but rather a combined effect of the tree age and its sociological position. The authors observed that two trees with same d and h, but with different ages, show differences in their sociological position. In this case, the younger tree may be dominant, while the older tree is likely to be suppressed. Thus, both effects, age and sociological position, cause an increase in stem biomass of the oldest tree. This occurs because the wood density increases with the stand age and the form factor (ratio between stem volume and
In addition, there is the effect of the biomass components allocation, which varies with the stand age, being also a common characteristic. In this study, the proportion of stem biomass in relation to the total biomass increased as a function of age, i.e., from the very young to mature forest, while the proportion of the crown biomass in relation to the total biomass decreased.
Therefore, the development of biomass equations according to age range can benefit forest managers when evaluating biomass stocks per component.
All intercept and exponential coefficients of the equations were significant. The intercept values were higher for the groups or systems of equations adjusted for the younger stands than for the older ones. Angular coefficients behave inversely. This inverse proportion relationship between the intercept and the angular coefficients was also observed by Fehrmann & Kleinn (2006) Fehrmann L, Kleinn C. General considerations about the use of allometric equations for biomass estimation on the example of Norway spruce in central Europe. Forest Ecology and Management 2006; 236(23): 412421. http://dx.doi.org/10.1016/j.foreco.2006.09.026.
http://dx.doi.org/10.1016/j.foreco.2006...
, Zianis (2008) Zianis D. Predicting mean aboveground forest biomass and its associated variance. Forest Ecology and Management 2008; 256(6): 14001407. http://dx.doi.org/10.1016/j.foreco.2008.07.002.
http://dx.doi.org/10.1016/j.foreco.2008...
, ZapataCuartas et al. (2012) ZapataCuartas M, Sierra C, Alleman L. Probability distribution of allometric coefficients and Bayesian estimation of aboveground tree biomass. Forest Ecology and Management 2012; 277: 173179. http://dx.doi.org/10.1016/j.foreco.2012.04.030.
http://dx.doi.org/10.1016/j.foreco.2012...
, and Sileshi (2014) Sileshi GW. A critical review of forest biomass estimation models, common mistakes and corrective measures. Forest Ecology and Management 2014; 329: 237254. http://dx.doi.org/10.1016/j.foreco.2014.06.026.
http://dx.doi.org/10.1016/j.foreco.2014...
and thoroughly discussed by Pilli et al. (2006) Pilli P, Anfodillo T, Carrer M. Towards a functional and simplified allometry for estimating forest biomass. Forest Ecology and Management 2006; 237(13): 583593. http://dx.doi.org/10.1016/j.foreco.2006.10.004.
http://dx.doi.org/10.1016/j.foreco.2006...
.
The dependence of the coefficients of the biomass equations on the stand ages is fundamental to generalize the black wattle biomass equations. The generalization of the equations – when they are adjusted for data from regional, state and national forest inventories – increases in importance and aims to create equations for various species, sites and ages. Generalization facilitates the application of the equations, because they are robust and present fewer restrictions, mainly owing to regional characteristic and age range of the tree sampling trees. In this scenario, the possibilities of using generalized biomass equations in the present case are clear, especially because stratification by age is not required anymore. Although this is an important theme, it was not considered in this work. We considered age stratification previous to define the equations.
The equations presented in this work are intended to be used in biomass inventories of black wattle stands and come from the main grown sites in the state of Rio Grande do Sul, covering a wide range of stand ages. Additional information regarding the variancecovariance matrix of coefficients and residuals was included in the appendix. Therefore, these equations are useful for inventories of biomass components and total biomass, specially the equations fitted by procedure 2 for equation systems.
The application of the equations independently adjusted (procedure 1) is not appropriate because they generate biologically inconsistent results, referred here as loss of compatibility between biomass components and total biomass. Most equations developed for black wattle and other species in Brazilian forests were independently adjusted. These equations are not additive and certain precautions should be considered before using them.
Biomass estimation for forest inventory purposes should use equations adjusted via procedure 2. These equations result in additivity of the estimates between biomass components and total biomass, a desirable and fundamental characteristic in modeling the biomass of trees. Ensuring additivity in tree biomass estimates has been recommended by several authors: Kozak (1970) Kozak A. Methods of ensuring additivity of biomass components by regression analysis. Forestry Chronicle 1970; 46(5): 402404. http://dx.doi.org/10.5558/tfc464025.
http://dx.doi.org/10.5558/tfc464025 ...
, Chiyenda & Kozak (1984) Chiyenda SS, Kozak A. Additivity of component biomass regression equations when the underlying model is linear. Canadian Journal of Forest Research 1984; 14(3): 441446. http://dx.doi.org/10.1139/x84078.
http://dx.doi.org/10.1139/x84078 ...
, Cunia & Briggs (1984) Cunia T, Briggs RD. Forcing additivity of biomass tables: some empirical results. Canadian Journal of Forest Research 1984; 14(3): 376384. http://dx.doi.org/10.1139/x84067.
http://dx.doi.org/10.1139/x84067 ...
, Parresol (1999 Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. Forest Science 1999; 45: 573593. , ^{2001} Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
http://dx.doi.org/10.1139/x00202 ...
), Carvalho & Parresol (2003) Carvalho JP, Parresol BR. Additivity in tree biomass components of Pyrenean oak ( Quercus pyrenaica Willd.). Forest Ecology and Management 2003; 179(13): 269276. http://dx.doi.org/10.1016/S03781127(02)005492.
http://dx.doi.org/10.1016/S03781127(02...
, Dong et al. (2014) Dong L, Zhang L, Li F. A compatible system of biomass equations for three conifer species in Northeast, China. Forest Ecology and Management 2014; 329: 306317. http://dx.doi.org/10.1016/j.foreco.2014.05.050.
http://dx.doi.org/10.1016/j.foreco.2014...
, Bi et al. (2015) Bi H, Murphy S, Volkova L, Weston C, Fairman T, Li Y et al. Additive biomass equations based on complete weighing of sample trees for open eucalypt forest species in southeastern Australia. Forest Ecology and Management 2015; 349: 106121. http://dx.doi.org/10.1016/j.foreco.2015.03.007.
http://dx.doi.org/10.1016/j.foreco.2015...
, Sanquetta et al. (2015a) Sanquetta CR, Behling A, Corte APD, Péllico S No, Schikowski AB, Amaral M. Simultaneous estimation as alternative to independent modeling of tree biomass. Annals of Forest Science 2015a; 72(8): 10991112. http://dx.doi.org/10.1007/s1359501504972.
http://dx.doi.org/10.1007/s1359501504...
, Zhao et al. (2015) Zhao D, Kane M, Markewitz D, Teskey R, Clutter M. Additive tree biomass equations for midrotation loblolly pine plantations. Forest Science 2015; 61(4): 613623. http://dx.doi.org/10.5849/forsci.14193.
http://dx.doi.org/10.5849/forsci.14193...
, Zheng et al. (2015) Zheng C, Mason EG, Jia L, Wei S, Sun S, Duan J. A singletree additive biomass model of Quercus variabilis Blume forest in North China. Trees 2015; 29(3): 705716. http://dx.doi.org/10.1007/s0046801411481.
http://dx.doi.org/10.1007/s0046801411...
, Affleck & DiéguezAranda (2016) Affleck DLR, DiéguezAranda U. Additive nonlinear biomass equations: a likelihoodbased approach. Forest Science 2016; 62(1): 129140. http://dx.doi.org/10.5849/forsci.15126.
http://dx.doi.org/10.5849/forsci.15126...
, and Poudel & Temesgen (2016) Poudel KP, Temesgen H. Methods for estimating aboveground biomass and its components for Douglasfir and lodgepole pine trees. Canadian Journal of Forest Research 2016; 46(1): 7787. http://dx.doi.org/10.1139/cjfr20150256.
http://dx.doi.org/10.1139/cjfr2015025...
. In this study, additivity of the components to compose the total biomass was ensured using procedure 2. An equation for each biomass component was obtained, and the equation for the total biomass was composed with the same variables and coefficients of the equations fitted to the components.
5. CONCLUSIONS
Using nonlinear equations is appropriate because they are more accurate than linear equations.
WNSUR model presents better biological and statistical properties for estimating allometric equations of biomass components and total biomass than the independent estimation. Therefore, it should be used in forest inventories.
The independent fitted equations result in nonadditivity of the estimates between biomass components and total biomass, an undesirable property. Moreover, they will result in greater variance estimates than those obtained by WNSUR equations.
Performing the modeling of variance structure and adjusting weighted equations is important to ensure reliable estimates of each biomass component and total biomass with minimum variance.

FINANCIAL SUPPORT The authors thank the TANAGRO S.A company for cession of the data and CNPq for granting a scholarship to one of the authors.
REFERENCES
 Affleck DLR, DiéguezAranda U. Additive nonlinear biomass equations: a likelihoodbased approach. Forest Science 2016; 62(1): 129140. http://dx.doi.org/10.5849/forsci.15126.
» http://dx.doi.org/10.5849/forsci.15126  Barichello LR, Schumacher MV, Vogel MLM. Quantificação da biomassa de Acacia mearnsii De Wild na região sul do Brasil. Ciência Florestal 2005; 15(2): 129135. http://dx.doi.org/10.5902/198050981830.
» http://dx.doi.org/10.5902/198050981830  Basuki TM, Van Laake PE, Skidmore AK, Hussin YA. Allometric equations for estimating the aboveground biomass in tropical lowland Dipterocarp forest. Forest Ecology and Management 2009; 257(8): 16841694. http://dx.doi.org/10.1016/j.foreco.2009.01.027.
» http://dx.doi.org/10.1016/j.foreco.2009.01.027  Behling A. Modelagem da biomassa de árvores para assegurar aditividade dos seus componentes [tese]. Curitiba: Universidade Federal do Paraná; 2016.
 Bi H, Long Y, Turner J, Lei Y, Snowdon P, Li Y et al. Additive prediction of aboveground biomass for Pinus radiata (D. Don) plantations. Forest Ecology and Management 2010; 259(12): 23012314. http://dx.doi.org/10.1016/j.foreco.2010.03.003.
» http://dx.doi.org/10.1016/j.foreco.2010.03.003  Bi H, Murphy S, Volkova L, Weston C, Fairman T, Li Y et al. Additive biomass equations based on complete weighing of sample trees for open eucalypt forest species in southeastern Australia. Forest Ecology and Management 2015; 349: 106121. http://dx.doi.org/10.1016/j.foreco.2015.03.007.
» http://dx.doi.org/10.1016/j.foreco.2015.03.007  Caldeira MVW, Saidelles FLF, Schumacher MV, Godinho TO. Biomassa de povoamento de Acacia mearnsii De Wild., Rio Grande do Sul, Brasil. Scientia Forestalis 2011; 39(90): 133141.
 Caldeira MVW. Quantificação da biomassa e do conteúdo de nutrientes em diferentes procedências de Acácianegra (Acacia mearnsii De Wild) [dissertação]. Santa Maria: Universidade Federal de Santa Maria; 1998.
 Carbonera Pereira J, Schumacher MV, Hoppe JM, Caldeira MVW, Santos EM. Produção de biomassa em um povoamento de Acacia mearnsii De Wild. no Estado do Rio Grande do Sul. Revista Árvore 1997; 21(4): 521526.
 Carvalho JP, Parresol BR. Additivity in tree biomass components of Pyrenean oak ( Quercus pyrenaica Willd.). Forest Ecology and Management 2003; 179(13): 269276. http://dx.doi.org/10.1016/S03781127(02)005492.
» http://dx.doi.org/10.1016/S03781127(02)005492  Chave J, RéjouMéchain M, Búrquez A, Chidumayo E, Colgan MS, Delitti WBC et al. Improved allometric models to estimate the aboveground biomass of tropical trees. Global Change Biology 2014; 20(10): 31773190. http://dx.doi.org/10.1111/gcb.12629. PMid:24817483.
» http://dx.doi.org/10.1111/gcb.12629  Chiyenda SS, Kozak A. Additivity of component biomass regression equations when the underlying model is linear. Canadian Journal of Forest Research 1984; 14(3): 441446. http://dx.doi.org/10.1139/x84078.
» http://dx.doi.org/10.1139/x84078  Cunia T, Briggs RD. Forcing additivity of biomass tables: some empirical results. Canadian Journal of Forest Research 1984; 14(3): 376384. http://dx.doi.org/10.1139/x84067.
» http://dx.doi.org/10.1139/x84067  Dong L, Zhang L, Li F. A compatible system of biomass equations for three conifer species in Northeast, China. Forest Ecology and Management 2014; 329: 306317. http://dx.doi.org/10.1016/j.foreco.2014.05.050.
» http://dx.doi.org/10.1016/j.foreco.2014.05.050  Dudley NS, Fownes JH. Preliminary biomass equations for eight species of fastgrowing tropical trees. Journal of Tropical Forest Science 1992; 5(1): 6873.
 Fehrmann L, Kleinn C. General considerations about the use of allometric equations for biomass estimation on the example of Norway spruce in central Europe. Forest Ecology and Management 2006; 236(23): 412421. http://dx.doi.org/10.1016/j.foreco.2006.09.026.
» http://dx.doi.org/10.1016/j.foreco.2006.09.026  Genet A, Wernsdörfer H, Jonard M, Pretzsch H, Rauch M, Ponette Q et al. Ontogeny partly explains the apparent heterogeneity of published biomass equations for Fagus sylvatica in central Europe. Forest Ecology and Management 2011; 261(7): 11881202. http://dx.doi.org/10.1016/j.foreco.2010.12.034.
» http://dx.doi.org/10.1016/j.foreco.2010.12.034  Greene WH. Econometric analysis Upper Saddle River: PrenticeHall; 2008.
 Harvey AC. Estimating regression models with multiplicative heteroscedasticity. Econometrica 1976; 44(3): 461465. http://dx.doi.org/10.2307/1913974.
» http://dx.doi.org/10.2307/1913974  Kozak A. Methods of ensuring additivity of biomass components by regression analysis. Forestry Chronicle 1970; 46(5): 402404. http://dx.doi.org/10.5558/tfc464025.
» http://dx.doi.org/10.5558/tfc464025  Maddala GS. Introduction to econometrics 3rd ed. New York: John Wiley & Sons; 2001.
 Mochiutti S. Produtividade e sustentabilidade de plantações de acácianegra (Acacia mearnsii De Wild.) no Rio Grande do Sul [tese]. Curitiba: Universidade Federal do Paraná; 2007.
 Niklas KJ. Sizedependent allometry of tree height, diameter and trunktaper. Annals of Botany 1995; 75(3): 217227. http://dx.doi.org/10.1006/anbo.1995.1015.
» http://dx.doi.org/10.1006/anbo.1995.1015  Parresol BR. Modeling multiplicative error variance: an example predicting tree diameter from stump dimensions in bald cypress. Forest Science 1993; 39(4): 670679.
 Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. Forest Science 1999; 45: 573593.
 Parresol BR. Additivity of nonlinear biomass equations. Canadian Journal of Forest Research 2001; 31(5): 865878. http://dx.doi.org/10.1139/x00202
» http://dx.doi.org/10.1139/x00202  Picard N, SaintAndré L, Henry M. Manual for building tree volume and biomass allometric equations: from field measurement to prediction Rome: Food and Agricultural Organization of the Unites Nations, Centre de Coopération Internationale en Recherche Agronomique pour le Développement; 2012.
 Pilli P, Anfodillo T, Carrer M. Towards a functional and simplified allometry for estimating forest biomass. Forest Ecology and Management 2006; 237(13): 583593. http://dx.doi.org/10.1016/j.foreco.2006.10.004.
» http://dx.doi.org/10.1016/j.foreco.2006.10.004  Poorter H, Niklas KJ, Reich PB, Oleksyn J, Poot P, Mommer L. Biomass allocation to leaves, stems and roots: metaanalyses of interspecific variation and environmental control. The New Phytologist 2012; 193(1): 3050. http://dx.doi.org/10.1111/j.14698137.2011.03952.x. PMid:22085245.
» http://dx.doi.org/10.1111/j.14698137.2011.03952.x  Poudel KP, Temesgen H. Methods for estimating aboveground biomass and its components for Douglasfir and lodgepole pine trees. Canadian Journal of Forest Research 2016; 46(1): 7787. http://dx.doi.org/10.1139/cjfr20150256.
» http://dx.doi.org/10.1139/cjfr20150256  Saidelles FLF. Determinação da biomassa e altura de amostragem para a quantificação de nutrientes em Acacia mearnsii De Wild [tese]. Santa Maria: Universidade Federal de Santa Maria; 2005.
 SaintAndré L, M’Bou AT, Mabiala A, Mouvondy W, Jourdan C, Roupsard O et al. Agerelated equations for above – and below: ground biomass of a Eucalyptus hybrid in Congo. Forest Ecology and Management 2005; 205(13): 199214. http://dx.doi.org/10.1016/j.foreco.2004.10.006.
» http://dx.doi.org/10.1016/j.foreco.2004.10.006  Sanquetta CR, Behling A, Corte AP, Simon A, Pscheidt H, Ruzza MS et al. Estoques de biomassa e carbono em povoamentos de acácia negra em diferentes idades no Rio Grande do Sul. Scientia Forestalis 2014; 42(103): 361370.
 Sanquetta CR, Behling A, Corte APD, Péllico S No, Schikowski AB, Amaral M. Simultaneous estimation as alternative to independent modeling of tree biomass. Annals of Forest Science 2015a; 72(8): 10991112. http://dx.doi.org/10.1007/s1359501504972.
» http://dx.doi.org/10.1007/s1359501504972  Sanquetta CR, Wojciechowski J, Dalla Corte AP, Behling A, Péllico S No, Rodrigues AL et al. Comparison of data mining and allometric model in estimation of tree biomass. BMC Bioinformatics 2015b; 16(247): 19. PMid:26250142.
 Sileshi GW. A critical review of forest biomass estimation models, common mistakes and corrective measures. Forest Ecology and Management 2014; 329: 237254. http://dx.doi.org/10.1016/j.foreco.2014.06.026.
» http://dx.doi.org/10.1016/j.foreco.2014.06.026  Steel RGD, Torrie JH, Dickey DA. Principles and procedures of statistics: a biometrical approach 3rd ed. New York: McGrawHill; 1996.
 Stein PP, Tonietto L. Black wattle silviculture in Brazil. In: Brown AG, Ko HC, editors. Black wattle and its utilization Barton: Rural Industries Research and Development Corporation; 1997.
 Wayson CA, Johnson KD, Cole JA, Olguín MI, Carrillo OI, Birdsey RA. Estimating uncertainty of allometric biomass equations with incomplete fit error information using a pseudodata approach: methods. Annals of Forest Science 2015; 72(6): 825834. http://dx.doi.org/10.1007/s1359501404367.
» http://dx.doi.org/10.1007/s1359501404367  Williams CJ, LePage BA, Vann DR, Tange T, Ikeda H, Ando M et al. Structure, allometry, and biomass of plantations Metasequoia glyptostroboides in Japan. Forest Ecology and Management 2003; 180(13): 287301. http://dx.doi.org/10.1016/S03781127(02)005674.
» http://dx.doi.org/10.1016/S03781127(02)005674  ZapataCuartas M, Sierra C, Alleman L. Probability distribution of allometric coefficients and Bayesian estimation of aboveground tree biomass. Forest Ecology and Management 2012; 277: 173179. http://dx.doi.org/10.1016/j.foreco.2012.04.030.
» http://dx.doi.org/10.1016/j.foreco.2012.04.030  Zhao D, Kane M, Markewitz D, Teskey R, Clutter M. Additive tree biomass equations for midrotation loblolly pine plantations. Forest Science 2015; 61(4): 613623. http://dx.doi.org/10.5849/forsci.14193.
» http://dx.doi.org/10.5849/forsci.14193  Zheng C, Mason EG, Jia L, Wei S, Sun S, Duan J. A singletree additive biomass model of Quercus variabilis Blume forest in North China. Trees 2015; 29(3): 705716. http://dx.doi.org/10.1007/s0046801411481.
» http://dx.doi.org/10.1007/s0046801411481  Zianis D. Predicting mean aboveground forest biomass and its associated variance. Forest Ecology and Management 2008; 256(6): 14001407. http://dx.doi.org/10.1016/j.foreco.2008.07.002.
» http://dx.doi.org/10.1016/j.foreco.2008.07.002
Publication Dates

Publication in this collection
08 Apr 2019 
Date of issue
2019
History

Received
06 Apr 2017 
Accepted
27 Oct 2017