Acessibilidade / Reportar erro

USE OF DESTRUCTIVE AND NON-DESTRUCTIVE METHODOLOGIES TO ESTIMATE STEM BIOMASS ACCUMULATION AND CARBON STOCK IN AN EUCALYPTUS FOREST

USO DE METODOLOGIAS DESTRUTIVAS E NÃO DESTRUTIVA PARA ESTIMAR ACÚMULO DE BIOMASSA NO TRONCO E ESTOQUE DE CARBONO EM UMA FLORESTA DE EUCALIPTO

ABSTRACT

Predicting wood biomass and carbon stock contents in planted forests can vary due to limitations associated with the measurement of parameters. Therefore, reducing possible errors generated over biomass and carbon stock quantification is an important step in obtaining reliable data. The study aimed to compare the use of destructive and non-destructive methodologies for predicting biomass and carbon stock in a planted Eucalyptus forest. Scaling was performed on 21 trees and 3 methodologies for carbon stock estimation were compared. For methodology 1, a control sample was harvested, sectioned, weighted in the field, and the carbon stock calculated based on these data. Methodology 2 was also destructive, as trees were harvested, scaled and the carbon stock predicted based on these data. Methodology 3 was non-destructive, as trees were scaled upright with the aid of equipment and the predicted carbon stock was based on these data. Biomass and carbon stock were compared by Test F and no statistical difference was observed. The data were separated according to diametric classes and compared by the Kolmogorov-Smirnov test, and again no significant difference was observed. Furthermore, three equations were generated based on the Schumacher & Hall model and compared by the identity test model and no differences between the methodologies were observed. Thus, both nondestructive and destructive methodologies herein evaluated were effective and showed equal results to the control sample. Moreover, the use of the non-destructive methodology reduces time and cost destined to predicting biomass and carbon stock.

Keywords:
Basic wood density; Forest plantation; Stem volume

RESUMO

A previsão de biomassa de madeira e estoque de carbono em florestas plantadas pode variar devido a limitações associadas à medição de parâmetros. Portanto, reduzir possíveis erros gerados na quantificação de biomassa e estoque de carbono é um passo importante na obtenção de dados confiáveis. O objetivo do estudo foi comparar o uso de metodologias destrutivas e não destrutivas para a previsão de biomassa e estoque de carbono em uma floresta plantada de eucalipto. As análises foram realizadas em 21 árvores e 3 metodologias para estimativa de estoque de carbono foram comparadas. Para a metodologia 1, uma amostra controle foi colhida, seccionada, pesada em campo e o estoque de carbono calculado com base nesses dados. A metodologia 2 também foi destrutiva, pois as árvores foram cortadas, cubadas e o estoque de carbono previsto com base nesses dados. A metodologia 3 foi não destrutiva, pois as árvores foram cubadas com auxílio de um equipamento e o estoque de carbono estimado foi baseado nesses dados. A biomassa e o estoque de carbono foram comparados pelo Teste F e nenhuma diferença estatística foi observada. Os dados foram separados de acordo com as classes diamétricas e comparados pelo teste de Kolmogorov-Smirnov, e novamente não foi observada diferença significativa. Além disso, três equações foram geradas com base no modelo de Schumacher & Hall e comparadas pelo modelo de teste de identidade e não foram observadas diferenças entre as metodologias. Assim, tanto as metodologias não destrutivas quanto as destrutivas aqui avaliadas foram eficazes e apresentaram resultados iguais à amostra controle. Além disso, o uso da metodologia não destrutiva reduz o tempo e o custo destinados à previsão de biomassa e estoque de carbono.

Palavras-Chave:
Densidade básica da madeira; Florestas plantadas; Volume do tronco

1. INTRODUCTION

Greenhouse gas (GHG) emissions have increased over the years and have caused an imbalance on Earth and, consequently, climate changes on a global scale (Olorunfemi et al., 2019Olorunfemi IE, Komofale AA, Fasinmirin JT, Olufayo AA. Biomass carbon stocks of different land use management in the forest vegetative zone of Nigeria. Acta Oecologia. 2019;95(1):45-56. doi: 10.1016/j.actao.2019.01.004
https://doi.org/10.1016/j.actao.2019.01....
). The principal anthropogenic sources of GHG are the burning of fossil fuels and the change in land use (IPCC, 2014Intergovernmental Panel on Climate Change — IPCC. Climate Change 2014: Mitigation of climate change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on climate change. UK and New York, USA: Cambridge University Press; 2014.). Given this situation, the development of strategies to reduce the concentration of atmospheric CO2 is a consensus. Forests are essential mitigators with a stock potential of 2-4 PgCO2e from the atmosphere (Qureshi et al., 2012Qureshi A, Pariva, Badola R, Hussain SA. A review of protocols used for assessment of carbon stock in forested landscapes. Environmental Science & Policy. 2012;16(1):81-89. doi: 10.1016/j.envsci.2011.11.001
https://doi.org/10.1016/j.envsci.2011.11...
), because they are able to store carbon as part of their biomass (Zhang et al., 2019Zhang H, Deng Q, Hui D, Wu J, Xiong X, Zhao J, et al. Recovery in soil carbon stock but reduction in carbon stabilization after 56-year forest restoration in degraded tropical lands. Forest Ecology and Management. 2019;441(1):1-8. doi: 10.1016/j. foreco.2019.03.037
https://doi.org/10.1016/j. foreco.2019.0...
).

Forests located in the tropics are in constant focus due to their high volumetric productivity and rapid growth (Achard et al., 2008Achard F, Eva HD, Mayaux P, Stibig HJ, Belward A. Improved estimates of net carbon emissions from land cover change in the tropics for the 1990s. Global Biogeochemical Cycles. 2008;18(1):1-11. doi: 10.1029/2003GB002142
https://doi.org/10.1029/2003GB002142...
). Therefore, accurate estimates of biomass production are needed to reduce uncertainties in the carbon stock potentials in those areas (Djomo et al., 2011Djomo NA, Knohl A, Gravenhorst G. Estimations of total ecosystem carbon pools distribution and carbon biomass current annual increment of a moist tropical forest. Forest Ecology and Management. 2011;261(8):1448-1459. doi: 10.1016/j.foreco.2011.01.031
https://doi.org/10.1016/j.foreco.2011.01...
). Estimates of volume, biomass, and carbon stock may have discrepancies associated with limitations in measuring parameters (Baccini et al., 2012Baccini A, Goetz SJ, Walker WS, Laporte NT, Sun M, Sulla-Menashe D, et al. Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nature Climate Change, 2012;2(1):182- 185. doi: 10.1038/nclimate1354
https://doi.org/10.1038/nclimate1354...
). Therefore, reducing eventual errors generated in the quantification is a significant step in obtaining reliable data (Stovall et al., 2017Stovall AEL, Vorster AF, Anderson RS, Evangelista PH, Shugart HH. Non-destructive aboveground biomass estimation of coniferous trees using terrestrial LiDAR. Remote Sensing of Environment. 2017;200(1):31-42. doi: 10.1016/j.rse.2017.08.013
https://doi.org/10.1016/j.rse.2017.08.01...
).

Producing precise and accurate biomass forecasts is challenging for several reasons. First, an impartial forest inventory project is required, with reliable measurements of trees’ attributes; further, requires biomass estimation models to accurately represent forest inventory data (Dutcâ et al., 2020Dutcâ I, Mather R, Loras F. Sampling trees to develop allometric biomass models: How does tree selection affect model prediction accuracy and precision? Ecological Indicators. 2020;117(1):106553. doi: 10.1016/j.ecolind.2020.106553
https://doi.org/10.1016/j.ecolind.2020.1...
). The methodologies usually used are defined as destructive, when trees inside a plot or trees previous selected from diametric classes are harvested and measured (Singh et al., 2011Singh V, Tewari A, Kushwaha SPS, Dadhwal V. Formulating allometric equations for estimating biomass and carbon stock in small diameter trees. Forest Ecology and Management. 2011;261(11):1945-1949. doi: 10.1016/j.foreco.2011.02.019
https://doi.org/10.1016/j.foreco.2011.02...
); and non-destructive, when is not necessary to cut trees (López-López et al., 2017López-López SF, Martínez-Trinidad T, Benavides-Meza H, Garcia-Nieto M, Santos-Posadas HM. Non-destructive method for above-ground biomass estimation of Fraxinus uhdei (Wenz.) Lingelsh in an urban forest. Urban Forestry & Urban Greening. 2017;24(1):62-70. doi: 10.1016/j.ufug.2017.03.025
https://doi.org/10.1016/j.ufug.2017.03.0...
). Non-destructive methodologies for estimating biomass are faster, cheaper, and avoid environmental problems resulting from tree felling (Mòntes, 2009Mòntes N. A non-destructive method to estimate biomass in arid environments: A comment on Flombaum and Sala. Journal of Arid Environments. 2009;73(1):599-601. doi: 10.1016/j.jaridenv.2008.08.003
https://doi.org/10.1016/j.jaridenv.2008....
).

Studies on forest biomass are carried out for different purposes, including knowing its energy potential, quantifying nutrient cycling (Silveira et al., 2008Silveira P, Koehler HS, Sanquetta CR, Arce JE. O estado da arte na estimativa de biomassa e carbono em formações florestais. Floresta. 2008;38(1):185-206. doi: 10.5380/rf.v38i1.11038
https://doi.org/10.5380/rf.v38i1.11038...
), monitoring tree growth (Zhao et al., 2018Zhao K, Suarez JC, Garcia M, Hu T, Wang C, Londo A. Utility of multitemporal lidar for forest and carbon monitoring: Tree growth, biomass dynamics, and carbon flux. Remote Sensing of Environment. 2018;204(1):883-897. doi: 10.1016/j.rse.2017.09.007
https://doi.org/10.1016/j.rse.2017.09.00...
), and carbon stock potential (Chieppa et al., 2020Chieppa J, Power SA, Tissue DT, Nielsen UN. Allometric estimates of aboveground biomass using cover and height are improved by increasing specificity of plant functional groups in Eastern Australian Rangelands. Rangeland Ecology & Management. 2020;73(3):375-383. doi: 10.1016/j.rama.2020.01.009
https://doi.org/10.1016/j.rama.2020.01.0...
). Destructive sampling, at the highest cost, is limited by capital, labor, logistics and bureaucracy, in the case of native forests.. Samples can be underrepresented in areas of complex topography and unfavorable climatic conditions (Picard et al., 2012Picard N, Saint-Andre L, Henry M. Manual for building tree volume and biomass allometric equations: from field measurement to prediction. Food and Agricultural Organization of the United Nations and Centre de Coopération Internationale en Recherche Agronomique pour le Développement. 2012 [cited 2021 Aug 21]. Available from: http://www.fao.org/docrep/018/i3058e/i3058e.pdf. E-ISBN: 978-92-5-107347-6
http://www.fao.org/docrep/018/i3058e/i30...
). Therefore, testing the accuracy of a faster, cheaper and simpler methodology for estimating biomass and carbon accumulation in forests becomes important.

Thus, this study aimed to compare the use of destructive and non-destructive methodologies for estimating biomass and carbon stock in a forest with a hybrid of Eucalyptus urophylla x Eucalyptus grandis. The hypothesis that motivated this study was the possibility of differences between the non-destructive methodologies concerning the destructive ones in the estimation of biomass and carbon.

2. MATERIAL AND METHODS

2.1. Characterization of the study site

The study was conducted on a charcoal-producing rural property in Lamim, Minas Gerais (20°47’08.56” S and 43°26’37.78” O), in the Zona da Mata (Figure 1). The tree component is a hybrid of Eucalyptus grandis x Eucalyptus urophylla, planted at a spacing of 3.0 m x 2.0 m. The plantation was 5 years old in the forest inventory and the silvicultural operations performed in the area were fertilizing and ant’s control. According to the Köppen’s classification, the climate of the region is Cwa, that is, subtropical with dry winter and hot and rainy summer. Precipitation occurs mainly between October and March, with averages of 1,435 mm per year. June and July present the lowest temperatures (12ºC), and January the highest temperatures (25ºC) (Sá Junior et al., 2012Sá Junior A, Carvalho LG, Da Silva FF, Alves MC. Application of the Köppen classification for climatic zoning in the state of Minas Gerais, Brazil. Theoretical and Applied Climatology. 2012;108(1): 1-7.).

Figure 1
Rural property in Lamim, MG, where the forest inventory was conducted to estimate volume, biomass, and carbon.
Figura 1
Propriedade rural em Lamim, MG, onde foi realizado o inventário florestal para estimativa de volume, biomassa e carbono.

2.2. Forest inventory and methodologies for estimating biomass and carbon stock

The forest inventory used simple random sampling, with 27 georeferenced plots of 300 m2 (20 x 15 m). All trees had the circumference at 1.30 m above the ground (cap) measured, converted to the diameter at 1.30 m above the ground (dbh), and separated into seven diametric classes with an amplitude of 2.5 cm. Three sample trees (chosen outside the sample units) were selected by diametric class, to perform the rigorous scaling through the destructive and non-destructive methods. A total of 21 trees were selected for strict scaling and used in the evaluated methodologies.

2.2.1. Methodology 1 - Destructive by Weighing -control

Methodology 1 was considered a control to compare with other methodologies, as it is the most accurate (Chave et al., 2014Chave J, Réjou-Mé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:3177-3190. doi: 10.1111/gcb.12629
https://doi.org/10.1111/gcb.12629...
). The sample trees were felled, and their stem was cut and weighed in the field. It is important to mention that the branches were not used. Wooden discs of 2.5 cm thickness at 0% (base), 25%, 50%, 75%, and 100% of tree’s commercial height, were removed and weighed immediately. The same samples were placed in a forced circulation oven with controlled temperature (100ºC for stem, and 40ºC for leaves) at the Madeira Panel and Energy Laboratory (LAPEM UFV) and weighed until dry weight stabilization.

The proportionality method was used to calculate the total dry biomass in the field, by section of the tree, after harvest, according to the following equation:

(Eq,1) W ( f ) = W W ( f ) * D W ( s ) / W W ( s )

Where: DW(f) = Field dry weight, in g; WW (f) = Field wet weight, in g; DW (s) = Sample dry weight, in g; WW (s) = Sample wet weight, in g.

Stem carbon stock was calculated using the 0.47 factor recommended for tropical forests (IPCC, 2006Intergovernmental Panel on Climate Change — IPCC. IPCC Guidelines for national greenhouse gas inventories. Kanagawa, JP: Prepared by the national greenhouse gas inventories programme. Institute for Global Environmental Strategies; 2006.).

2.2.2. Methodology 2 - Destructive with scaling

Sample trees were felled, and their diameters with bark were measured at heights of 0 m, 0.30 m, 0.70 m, 1.00 m, and 1.30 m, using a tape measure. From this height on, measurements were taken every 1.00 meter until the minimum commercial diameter of 3 cm. The volume in each of the sections was calculated using Smalian’s formula (Eq. 2).

(Eq.2) V b = ( S A 1 + S A 2 ) / 2 * L

Where: Vcc – Volume with bark, in m3; SA1 – Sectional area of the stem lower part, in m2; SA2 – Sectional area of the upper stem, in m2; L – Stem section length, in m.

At 0% (base), 25%, 50%, 75%, and 100% of tree’s commercial height, wooden discs of 2.5 cm thickness were removed. Their opposite wedges were used to determine the basic wood density according to ABNT NBR 11941 (ABNT, 2003Associação Brasileira de Normas Técnicas — ABNT. NBR 11941-02: Determinação da densidade básica em madeira. Rio de Janeiro; 2003.). The average value of the basic wood density of the opposite wedges was used to estimate each wooden disc biomass. The biomass of the stem was obtained by multiplying its volume with bark by the average basic wood density; carbon stock was calculated using the 0.47 factor, recommended for tree species (IPCC, 2006Intergovernmental Panel on Climate Change — IPCC. IPCC Guidelines for national greenhouse gas inventories. Kanagawa, JP: Prepared by the national greenhouse gas inventories programme. Institute for Global Environmental Strategies; 2006.).

2.2.3. Methodology 3 - Non-destructive with a Wheeler Pentaprism Caliper

Sample trees (still standing) had their diameters with bark at heights of 0 m, 0.30 m, 0.70 m, 1.00 m, and 1.30 m measured. From this height on, measurements were taken every 1.00 meter using a Wheeler® Pentaprism caliper until the minimum commercial diameter of 6.5 cm. The volume in each of the sections was calculated using Smalian’s formula.

Where: Vb – Volume with bark, in m3; SA1 – Sectional area of the stem lower part, in m2; SA2 – Sectional area of the upper stem, in m2; L – Stem section length, in m.

The volume of the stem tip (stem portion above the minimum commercial diameter of 6.5 cm) was calculated using the formula for the volume of a cone.

Where: Vcone – Cone volume, in m3; SA1 – Sectional area of the stem lower part, in m2; L – Stem section length, in m.

The volumes obtained using equations 3 and 4 were summed up to obtain the stem’s total volume. A wood sample was taken from each sampled tree using a manual auger at 1.30 m above the ground (dbh) to determine the basic wood density according to ABNT NBR 11941 (ABNT, 2003Associação Brasileira de Normas Técnicas — ABNT. NBR 11941-02: Determinação da densidade básica em madeira. Rio de Janeiro; 2003.). The biomass of the stem was obtained by multiplying the volume with bark by the basic wood density of each individual. The calculated stem biomass was converted into carbon stock by multiplying by 0.47 (IPCC, 2006Intergovernmental Panel on Climate Change — IPCC. IPCC Guidelines for national greenhouse gas inventories. Kanagawa, JP: Prepared by the national greenhouse gas inventories programme. Institute for Global Environmental Strategies; 2006.).

(Eq.3) V b = ( S A 1 + S A 2 ) / 2 * L

(Eq.4) Vcone = ( S A 1 * L ) / 3

2.3. Data Statistical Analysis

The results were interpreted with the analysis of variance (ANOVA) to compare the carbon stock between methodologies 1, 2 and 3. In case of statistical difference in carbon stock, the Test F would be applied for methodologies 2 and 3 in relation to methodology 1, separately. If significant differences were detected in carbon stock, the values would be compared by the Test T for paired samples, at 95% probability. A residual analysis was performed to compare the means estimated values for carbon stock per diameter classes obtained using Methodologies 2 and 3 with the ones obtained in Methodology 1 (control).

Kolmogorov-Smirnov test was used to compare the statistical significance between the carbon stock obtained using methodology 1 (control) with the others, per diametric classes, at 95% probability.

(Eq.5) D cal = Max ( F 0 ( x ) F e ( x ) )

Where: Dcal - value for 5% significance obtained; F0 (x) = cumulative frequency observed; Fe (x) = expected cumulative frequency.

Dcal value for 5% significance was obtained according to equation 6. If Dcal<Dtab:Ho is accepted (observed distribution equal to projected); if DcalDtab:Ho is rejected (observed distribution is not equal to the projected distribution).

(Eq.6) D tab = 1.35 / n

Where: Dtab = critical value at 5% significance and “n” is the number of observations.

All statistical analyzes were performed using the R software (R CORE TEAM, 2021R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.).

2.4. Model Identity

An equation based on Schumacher and Hall (1933)Schumacher FX, Hall FS. Logarithmic expression of timber-tree volume. Journal of Agricultural Research. 1933;47(9):719-734. model was adjusted for each one of the tested methodologies to estimate carbon stock using diameter at breast height (dbh) and commercial height (Ht) from the sampled trees.

(Eq.7) C = β 0 * d b h β 1 * H t β 2

Where: C – carbon stock, in Mg; βn – model parameters; dbh – diameter with bark measured at 1.30 m from the ground, in cm; Ht – commercial height of the sample trees, in m.

The verification of the model adequacy was carried out based on the analysis of the adjusted determination coefficient (R2adj), Bias (%), and RMSE (%).

A model identity test (Graybill, 1976Graybill FA. Theory and application of linear model. Belmont: Duxbury; 1976. v.1.) was used to group the carbon stock estimation models, in relation to the control, to a significance of 5%. The test consists of reducing the sum of squares, allowing to statistically verify, by the Test F, the significance of the difference between the total sums of squares of the regressions adjusted for each methodology alone (complete model) and the sum of the squares of the regression adjusted for the total data set (reduced model). The verification of the identity of models in the forest area becomes a useful tool in modeling analysis, in an attempt to reduce the number of equations without loss of precision in the estimates, in addition to reducing sampling costs and economy of operations in use cases of a common equation.

The tested hypotheses were:

  • -H0: the reduced model adjusted for the total data set obtained using methodologies 2 and 3 in relation to the control (methodology 1) does not statistically differ from the adjusted complete models.

  • -H1: H0 is rejected.

These analyses were performed using the Microsoft Excel.

3. RESULTS

The average carbon stock obtained using methodology 1 was 0.0438 ± 0.0308 MgC, a value similar to those found using methodologies 2 (0.0470 ± 0.0343 MgC) and 3 (0.0431 ± 0.0345 MgC) (Table 1).

Table 1
Values of volume (Vol, in m3), wood density with standard deviation as a function of the samples taken along the shaft (Dens, in g cm3-1), and carbon stock (Carb, MgC) for the 21 sample trees evaluated using the methodologies 1 (Destructive with weighing - control), 2 (Destructive with scaling) and 3 (Non-destructive with Pentaprism).
Tabela 1
Valores de volume (Vol, em m3), densidade da madeira com desvio padrão em função das amostras tomadas ao longo do fuste (Dens, em g cm3-1) e estoque de carbono (Carb, MgC) para as 21 árvores amostras avaliadas pelas metodologias 1 (Destrutivo com pesagem - controle), 2 (Destrutivo com cubagem) e 3 (Não destrutiva com Pentaprisma).

The comparison using ANOVA between methodologies 2 (Value F – 0.102 < Value P – 0.751) and 3 (Value F – 0.006 < Value P – 0.941) with the control (Methodology 1), showed no significant difference between the carbon stock data.

The residual analysis plots showed methodology 2 performed better than methodology 3, with a steady trend around the identity line (Figure 2). Methodology 3 had an overestimation of data in the two lower-class centers (6.25 and 8.75 cm).

Figure 2
Residual analysis plots for methodology 2 (2A) and 3 (2B).
Figura 2
Gráficos de análise de resíduos para as metodologias 2 (2A) e 3 (2B).

Kolmogorov-Smirnov test’s evaluation resulted in a non-statistical difference between the carbon stock by diametric class for the methodologies 2 (Dcalc – 0.010 < Dtab – 1.407) and 3 (Dcalc – 0.023 < Dtab – 1.407) in relation to the control – Methodology 1.

The adjustment of the equations to estimate the carbon stock using data of each one of the tested methodologies was considered adequate, with satisfactory R2adj, RMSE (%), and Bias (%) values (Table 2).

Table 2
Parameters and adjustment of models to estimate carbon stock.
Tabela 2
Parâmetros e ajuste de modelos para estimar o estoque de carbono.

The model identity test showed the same behavior for the combination of Methodologies 2 and 3, that is, a single equation can be used to estimate biomass and carbon stock.

4. DISCUSSIONS

The quantification of biomass accumulation is an essential step to understand the carbon dynamics in forests and their ecosystem services (Houghton et al., 2009Houghton RA, Hall F, Goetz SJ. Importance of biomass in the global carbon cycle. Journal of Geophysical Research. 2009;114(1):1-13. doi: 10.1029/2009JG000935
https://doi.org/10.1029/2009JG000935...
), as it is a relevant component of carbon stocks and assessment of climate changes potential mitigation (Huy et al., 2016Huy B, Kralicek K, Poudel KP, Phuong VT, Khoa P V, Hung ND, et al. Allometric equations for estimating tree aboveground biomass in evergreen broadleaf forests of Vietnam. Forest Ecology and Management. 2016;382(1):193-205. doi: 10.1016/j.foreco.2016.10.021
https://doi.org/10.1016/j.foreco.2016.10...
). Thus, reliable biomass estimations are essential to monitor forest conditions and support decision-making under forest management (Ubuy et al., 2018Ubuy MH, Eid T, Bollandsas OM, Birhane E. Aboveground biomass models for trees and shrubs of exclosures in thedrylands of Tigray, northern Ethiopia. Journal of Arid Environments. 2018;156(1):9-18. doi: 10.1016/j.jaridenv.2018.05.007
https://doi.org/10.1016/j.jaridenv.2018....
).

The generation of reliable data on the carbon stock potential of forests is relevant in the current political momentum, in which the Paris Agreement is already in force, and some countries that have ratified it, such as Brazil, have emission reduction targets in the forestry sector (Azevedo-Ramos et al., 2020Azevedo-Ramos C, Moutinho P, Arruda VLS, Stabile MCC, Alencar A, Castro I, et al. Lawless land in no man’s land: The undesignated public forests in the Brazilian Amazon. Land Use Policy. 2020;99(1):104863. doi: 10.1016/j.landusepol.2020.104863
https://doi.org/10.1016/j.landusepol.202...
). The Brazilian government estimates that by the year 2030, the area of commercial forests will be increased by 3 million hectares, with varied stock potential, which highlights the importance of validation of biomass and carbon estimation methodologies (Brasil, 2015Brasil. Documento base para subsidiar os diálogos estruturados sobre a elaboração de uma estratégia de implementação e financiamento da contribuição nacionalmente determinada do Brasil no Acordo de Paris. 2015 [cited 2021 May 14]. Available from: https://www.mma.gov.br/images/arquivo/80051/NDC/documento_base_ndc_2_2017.pdf
https://www.mma.gov.br/images/arquivo/80...
).

The equivalence between the results of the tested methodologies with the control one is evidenced by the low difference between the carbon stock numbers presented in the results. This fact can be explained by the number of sections measured in the rigorous scaling, which contributes to a reduction in the estimation error, due the increase in the control of the tree taper (Tonini et al., 2019Tonini H, Morales MM, Silva VP, Lulu J, Farias Neto Al. Effect of planting system and solar exposure on biomass allocation in the initial growth of eucalyptus. Ciencia Florestal. 2019;29(1):96-95. doi: 10.5902/1980509817808
https://doi.org/10.5902/1980509817808...
). The use of a Wheeler’s Pentaprism (methodology 3) also allowed the generation of reliable results for carbon stock estimation. The use of this device is recommended for scaling Eucalyptus trees up to 50 m height, with a good precision in the generated estimations (Avery and Burkhart, 1997Avery TE, Burkhart HE. Forest measurements. 4.ed. New York: McGrawHill; 1997. ISBN: 007113204x.). However, tree size can affect the results obtained with the Pentaprism due to the difficulty in collecting the diameter of the section in the correct position. Another important factor to be mentioned is the importance of operator training to maintain the estimates accuracy.

The residuals values found in the lowest diameter classes for methodology 3 (37.91% and 139.54%, respectively for the 6.25 cm and 8.75 cm classes), despite having great magnitude, does not negatively impact carbon stock estimation because these values represent classes with less biomass accumulation. In a diametric distribution for a eucalyptus plantation, most individuals are concentrated in the middle-center diameter classes. Therefore, the impact of residuals values in the population’s carbon stock estimation is not so important because it affects a smaller number of trees (Nogueira et al., 2005Nogueira GS, Leite HG, Campos JCC, Carvalho AF, Souza AL. Modelo de distribuição diamétrica para povoamentos de Eucalyptus sp. submetidos a desbaste. Revista Árvore. 2005;29(4):579-589. doi: 10.1590/S0100-67622005000400010
https://doi.org/10.1590/S0100-6762200500...
).

The destructive method has some negative points when compared to indirect methodologies. The time required to carry out fieldwork is longer than in indirect methods (Flombaum and Sala, 2007Flombaum P, Sala OE. A non-destructive and rapid method to estimate biomass and aboveground net primary production in arid environments. Journal of Arid Environments. 2007;69(1):352-358. doi: 10.1016/j.jaridenv.2006.09.008
https://doi.org/10.1016/j.jaridenv.2006....
). Destructive methodologies are also limited to smaller areas with a small number of trees to be felled (Lu et al., 2014Lu D, Chen Q, Wang G, Liu L, Li G, Moran E. A survey of remote sensing-based aboveground biomass estimation methods in forest ecosystems. International Journal of Digital Earth. 2014;9(1):63-105. doi: 10.1080/17538947.2014.990526
https://doi.org/10.1080/17538947.2014.99...
). Sampling errors can also be a problem in direct methodologies, with trees selected wrongly (Brown et al., 1989Brown S, Gillespie A, Lugo AE. Biomass estimation methods for tropical forests with applications to forest inventory data. Forest Science. 1989;35(4):881-902.), which would lead to tendency errors and subsequent overestimation or underestimation of biomass accumulation (Ribeiro et al., 2009Ribeiro SC, Jacovine LAG, Soares CPB, Martins SV, Lopes De Souza A, Nardelli AMB. Quantificação de biomassa e estimativa de estoque de carbono em uma floresta madura no município de Viçosa, Minas Gerais. Revista Árvore. 2009;33(5):917-926. doi: 10.1590/S0100-67622009000500014
https://doi.org/10.1590/S0100-6762200900...
).

For Eucalyptus forests, the pentaprism proved to be a reliable tool for carbon stock estimation, with no observed statistical difference for estimates for population or diameter classes. The search for nondestructive methodologies that reliably estimate the accumulation of biomass and carbon stock is the focus of the researchers due to faster service execution and lower cost of data collection (Huff et al., 2018Huff S, Poudel KP, Ritchie M, Temesgen H. Quantifying aboveground biomass for common shrubs in northeastern California using nonlinear mixed effect models. Forest Ecology and Management. 2018;424(1):154-163. doi: 10.1016/j.foreco.2018.04.043
https://doi.org/10.1016/j.foreco.2018.04...
; Kramer et al., 2018Kramer RD, Sillett SC, Van Pelt R. Quantifying aboveground components of Picea sitchensis for allometric comparisons among tall conifers in North American rainforests. Forest Ecology and Management. 2018;430(1):59-77. doi: 10.1016/j.foreco.2018.07.039
https://doi.org/10.1016/j.foreco.2018.07...
) and, despite the possible uncertainties surrounding them, the need for data from direct methodologies demonstrates the importance of these methods in research related to the topic.

5. CONCLUSIONS

There are no differences, according to the data, in the biomass and carbon stock estimation between destructive and non-destructive methodologies in an Eucalyptus forest.

The non-destructive methodology and the destructive one with rigorous scaling is effective, with statistically similar results to the reference methodology, which reduces time and cost in estimating biomass and carbon in eucalyptus forests without compromising the result.

6. REFERENCES

  • Achard F, Eva HD, Mayaux P, Stibig HJ, Belward A. Improved estimates of net carbon emissions from land cover change in the tropics for the 1990s. Global Biogeochemical Cycles. 2008;18(1):1-11. doi: 10.1029/2003GB002142
    » https://doi.org/10.1029/2003GB002142
  • Associação Brasileira de Normas Técnicas — ABNT. NBR 11941-02: Determinação da densidade básica em madeira. Rio de Janeiro; 2003.
  • Avery TE, Burkhart HE. Forest measurements. 4.ed. New York: McGrawHill; 1997. ISBN: 007113204x.
  • Azevedo-Ramos C, Moutinho P, Arruda VLS, Stabile MCC, Alencar A, Castro I, et al. Lawless land in no man’s land: The undesignated public forests in the Brazilian Amazon. Land Use Policy. 2020;99(1):104863. doi: 10.1016/j.landusepol.2020.104863
    » https://doi.org/10.1016/j.landusepol.2020.104863
  • Baccini A, Goetz SJ, Walker WS, Laporte NT, Sun M, Sulla-Menashe D, et al. Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nature Climate Change, 2012;2(1):182- 185. doi: 10.1038/nclimate1354
    » https://doi.org/10.1038/nclimate1354
  • Brasil. Documento base para subsidiar os diálogos estruturados sobre a elaboração de uma estratégia de implementação e financiamento da contribuição nacionalmente determinada do Brasil no Acordo de Paris. 2015 [cited 2021 May 14]. Available from: https://www.mma.gov.br/images/arquivo/80051/NDC/documento_base_ndc_2_2017.pdf
    » https://www.mma.gov.br/images/arquivo/80051/NDC/documento_base_ndc_2_2017.pdf
  • Brown S, Gillespie A, Lugo AE. Biomass estimation methods for tropical forests with applications to forest inventory data. Forest Science. 1989;35(4):881-902.
  • Chave J, Réjou-Mé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:3177-3190. doi: 10.1111/gcb.12629
    » https://doi.org/10.1111/gcb.12629
  • Chieppa J, Power SA, Tissue DT, Nielsen UN. Allometric estimates of aboveground biomass using cover and height are improved by increasing specificity of plant functional groups in Eastern Australian Rangelands. Rangeland Ecology & Management. 2020;73(3):375-383. doi: 10.1016/j.rama.2020.01.009
    » https://doi.org/10.1016/j.rama.2020.01.009
  • Djomo NA, Knohl A, Gravenhorst G. Estimations of total ecosystem carbon pools distribution and carbon biomass current annual increment of a moist tropical forest. Forest Ecology and Management. 2011;261(8):1448-1459. doi: 10.1016/j.foreco.2011.01.031
    » https://doi.org/10.1016/j.foreco.2011.01.031
  • Dutcâ I, Mather R, Loras F. Sampling trees to develop allometric biomass models: How does tree selection affect model prediction accuracy and precision? Ecological Indicators. 2020;117(1):106553. doi: 10.1016/j.ecolind.2020.106553
    » https://doi.org/10.1016/j.ecolind.2020.106553
  • Flombaum P, Sala OE. A non-destructive and rapid method to estimate biomass and aboveground net primary production in arid environments. Journal of Arid Environments. 2007;69(1):352-358. doi: 10.1016/j.jaridenv.2006.09.008
    » https://doi.org/10.1016/j.jaridenv.2006.09.008
  • Graybill FA. Theory and application of linear model. Belmont: Duxbury; 1976. v.1.
  • Houghton RA, Hall F, Goetz SJ. Importance of biomass in the global carbon cycle. Journal of Geophysical Research. 2009;114(1):1-13. doi: 10.1029/2009JG000935
    » https://doi.org/10.1029/2009JG000935
  • Huff S, Poudel KP, Ritchie M, Temesgen H. Quantifying aboveground biomass for common shrubs in northeastern California using nonlinear mixed effect models. Forest Ecology and Management. 2018;424(1):154-163. doi: 10.1016/j.foreco.2018.04.043
    » https://doi.org/10.1016/j.foreco.2018.04.043
  • Huy B, Kralicek K, Poudel KP, Phuong VT, Khoa P V, Hung ND, et al. Allometric equations for estimating tree aboveground biomass in evergreen broadleaf forests of Vietnam. Forest Ecology and Management. 2016;382(1):193-205. doi: 10.1016/j.foreco.2016.10.021
    » https://doi.org/10.1016/j.foreco.2016.10.021
  • Intergovernmental Panel on Climate Change — IPCC. Climate Change 2014: Mitigation of climate change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on climate change. UK and New York, USA: Cambridge University Press; 2014.
  • Intergovernmental Panel on Climate Change — IPCC. IPCC Guidelines for national greenhouse gas inventories. Kanagawa, JP: Prepared by the national greenhouse gas inventories programme. Institute for Global Environmental Strategies; 2006.
  • Kramer RD, Sillett SC, Van Pelt R. Quantifying aboveground components of Picea sitchensis for allometric comparisons among tall conifers in North American rainforests. Forest Ecology and Management. 2018;430(1):59-77. doi: 10.1016/j.foreco.2018.07.039
    » https://doi.org/10.1016/j.foreco.2018.07.039
  • López-López SF, Martínez-Trinidad T, Benavides-Meza H, Garcia-Nieto M, Santos-Posadas HM. Non-destructive method for above-ground biomass estimation of Fraxinus uhdei (Wenz.) Lingelsh in an urban forest. Urban Forestry & Urban Greening. 2017;24(1):62-70. doi: 10.1016/j.ufug.2017.03.025
    » https://doi.org/10.1016/j.ufug.2017.03.025
  • Lu D, Chen Q, Wang G, Liu L, Li G, Moran E. A survey of remote sensing-based aboveground biomass estimation methods in forest ecosystems. International Journal of Digital Earth. 2014;9(1):63-105. doi: 10.1080/17538947.2014.990526
    » https://doi.org/10.1080/17538947.2014.990526
  • Mòntes N. A non-destructive method to estimate biomass in arid environments: A comment on Flombaum and Sala. Journal of Arid Environments. 2009;73(1):599-601. doi: 10.1016/j.jaridenv.2008.08.003
    » https://doi.org/10.1016/j.jaridenv.2008.08.003
  • Nogueira GS, Leite HG, Campos JCC, Carvalho AF, Souza AL. Modelo de distribuição diamétrica para povoamentos de Eucalyptus sp. submetidos a desbaste. Revista Árvore. 2005;29(4):579-589. doi: 10.1590/S0100-67622005000400010
    » https://doi.org/10.1590/S0100-67622005000400010
  • Olorunfemi IE, Komofale AA, Fasinmirin JT, Olufayo AA. Biomass carbon stocks of different land use management in the forest vegetative zone of Nigeria. Acta Oecologia. 2019;95(1):45-56. doi: 10.1016/j.actao.2019.01.004
    » https://doi.org/10.1016/j.actao.2019.01.004
  • Picard N, Saint-Andre L, Henry M. Manual for building tree volume and biomass allometric equations: from field measurement to prediction. Food and Agricultural Organization of the United Nations and Centre de Coopération Internationale en Recherche Agronomique pour le Développement. 2012 [cited 2021 Aug 21]. Available from: http://www.fao.org/docrep/018/i3058e/i3058e.pdf E-ISBN: 978-92-5-107347-6
    » http://www.fao.org/docrep/018/i3058e/i3058e.pdf
  • Qureshi A, Pariva, Badola R, Hussain SA. A review of protocols used for assessment of carbon stock in forested landscapes. Environmental Science & Policy. 2012;16(1):81-89. doi: 10.1016/j.envsci.2011.11.001
    » https://doi.org/10.1016/j.envsci.2011.11.001
  • R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.
  • Ribeiro SC, Jacovine LAG, Soares CPB, Martins SV, Lopes De Souza A, Nardelli AMB. Quantificação de biomassa e estimativa de estoque de carbono em uma floresta madura no município de Viçosa, Minas Gerais. Revista Árvore. 2009;33(5):917-926. doi: 10.1590/S0100-67622009000500014
    » https://doi.org/10.1590/S0100-67622009000500014
  • Sá Junior A, Carvalho LG, Da Silva FF, Alves MC. Application of the Köppen classification for climatic zoning in the state of Minas Gerais, Brazil. Theoretical and Applied Climatology. 2012;108(1): 1-7.
  • Schumacher FX, Hall FS. Logarithmic expression of timber-tree volume. Journal of Agricultural Research. 1933;47(9):719-734.
  • Silveira P, Koehler HS, Sanquetta CR, Arce JE. O estado da arte na estimativa de biomassa e carbono em formações florestais. Floresta. 2008;38(1):185-206. doi: 10.5380/rf.v38i1.11038
    » https://doi.org/10.5380/rf.v38i1.11038
  • Singh V, Tewari A, Kushwaha SPS, Dadhwal V. Formulating allometric equations for estimating biomass and carbon stock in small diameter trees. Forest Ecology and Management. 2011;261(11):1945-1949. doi: 10.1016/j.foreco.2011.02.019
    » https://doi.org/10.1016/j.foreco.2011.02.019
  • Stovall AEL, Vorster AF, Anderson RS, Evangelista PH, Shugart HH. Non-destructive aboveground biomass estimation of coniferous trees using terrestrial LiDAR. Remote Sensing of Environment. 2017;200(1):31-42. doi: 10.1016/j.rse.2017.08.013
    » https://doi.org/10.1016/j.rse.2017.08.013
  • Tonini H, Morales MM, Silva VP, Lulu J, Farias Neto Al. Effect of planting system and solar exposure on biomass allocation in the initial growth of eucalyptus. Ciencia Florestal. 2019;29(1):96-95. doi: 10.5902/1980509817808
    » https://doi.org/10.5902/1980509817808
  • Ubuy MH, Eid T, Bollandsas OM, Birhane E. Aboveground biomass models for trees and shrubs of exclosures in thedrylands of Tigray, northern Ethiopia. Journal of Arid Environments. 2018;156(1):9-18. doi: 10.1016/j.jaridenv.2018.05.007
    » https://doi.org/10.1016/j.jaridenv.2018.05.007
  • Zhang H, Deng Q, Hui D, Wu J, Xiong X, Zhao J, et al. Recovery in soil carbon stock but reduction in carbon stabilization after 56-year forest restoration in degraded tropical lands. Forest Ecology and Management. 2019;441(1):1-8. doi: 10.1016/j. foreco.2019.03.037
    » https://doi.org/10.1016/j. foreco.2019.03.037
  • Zhao K, Suarez JC, Garcia M, Hu T, Wang C, Londo A. Utility of multitemporal lidar for forest and carbon monitoring: Tree growth, biomass dynamics, and carbon flux. Remote Sensing of Environment. 2018;204(1):883-897. doi: 10.1016/j.rse.2017.09.007
    » https://doi.org/10.1016/j.rse.2017.09.007

Publication Dates

  • Publication in this collection
    16 Mar 2022
  • Date of issue
    2022

History

  • Received
    08 Sept 2021
  • Accepted
    02 Feb 2022
Sociedade de Investigações Florestais Universidade Federal de Viçosa, CEP: 36570-900 - Viçosa - Minas Gerais - Brazil, Tel: (55 31) 3612-3959 - Viçosa - MG - Brazil
E-mail: rarvore@sif.org.br