SciELO - Scientific Electronic Library Online

 
vol.20 número47Alternativas de modelagem para a estimativa da altura de árvores de eucaliptoCharacterization of a Pteridium arachnoideum (Kaulf.) Maxon population and the vegetation associated with the anthropic paramo in Southern Ecuador índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista Forestal Mesoamericana Kurú

versão On-line ISSN 2215-2504

Kurú vol.20 no.47 Cartago Jul./Dez. 2023

http://dx.doi.org/10.18845/rfmk.v20i47.6819 

Artículo

Equações de taper ajustadas para Corymbia citriodora Hill & Johnson empregando dados gerados pelo método da altura relativa

Taper equations adjusted for Corimbya citriodora Hill & Johnson using data generated by the relative height method

David Lucas Camargo Vieira Terra1 

Valdir Carlos Lima de Andrade2 

1. Engenheiro Florestal e Segurança do Trabalho, Lajoal Montagens, Rio Grande do Sul, Brasil, davidlcvt7@hotmail.com

2. Curso de Engenharia Florestal, Universidade Federal do Tocantins, Gurupi, Tocantins, Brasil, vclandrade@uft.edu.br

Resumo

Neste estudo o objetivo foi avaliar o ajuste de equações para descrever o perfil do fuste de Corymbia citriodora empregando-se dados de cubagem coletados diretamente ao longo do mesmo (cubagem real) e indiretamente por meio do método da altura relativa (cubagem hr). Para isso, os diâmetros do fuste foram obtidos nas posições: 0,2 m, 0,4 m, 0,7 m, 1,3 m, 2,70 m e, sucessivamente, de 2 m em 2 m até um diâmetro em torno de 3 cm com casca. Na avaliação da descrição do perfil do fuste, com as equações de taper geradas com dados das cubagens real e hr e do volume pela aplicação da fórmula de Huber, foram adotados os critérios estatísticos: resíduo médio, Bias, correlação linear, desvio padrão dos resíduos, exatidão obtida pelo teste de qui-quadrado e raiz quadrada do resíduo médio, além de análise da distribuição dos resíduos. Tais critérios foram obtidos em uma validação cruzada utilizando-se dois bancos de dados de cubagem obtidos pela alternância das posições de medição do diâmetro do fuste. Concluiu-se por haver evidência de se obter níveis de acurácia bem próximos no emprego das equações de taper geradas a partir dos dados de cubagem real e cubagem hr, tanto na descrição do perfil do fuste de C. citriodora, quanto na quantificação do seu volume.

Palavras-chave: Afilamento; mensuração florestal; multiprodutos madeireiros

Abstract

This study evaluated the fit of equations to describe the stem shape of the tree species Corymbia citriodora using cubing data collected directly along the stem (real cubing) and indirectly by measuring the relative height (hr cubing). For this, stem diameters were obtained at different heights: 0.2 m, 0.4 m, 0.7 m, 1.3 m, 2.70 m, and successively, every 2 m up to a diameter around 3 cm with bark. Stem shape was calculated using taper equations adding real cubing data and hr cubing, also the volume through the Huber formula. The statistical criteria used included: mean residue, Bias, linear correlation, standard deviation of the residues, accuracy obtained by the chi-square test and the square root of the residue, in addition to analysis of the distribution of residues. Such criteria were obtained through cross-validation using two cubing databases obtained by alternating the measurement positions of the stem diameter. As a conclusion, there is evidence of fine levels of accuracy when using taper equations generated from real cubing data and hr cubing for both stem profile description of C. citriodora and volume quantification.

Keywords: Tapering; Forest Measurement; multi-timber products

Introdução

A quantificação de multiprodutos em povoamentos florestais permite delinear o melhor aproveitamento possível da madeira, principalmente, com um cenário favorável à tomada de decisões frente às oscilações do mercado de produtos madeireiros. Para isso, é fundamental conhecer o perfil do fuste das árvores a fim de realizar o sortimento dos diferentes usos desejáveis. Tal exigência, demanda o emprego de uma equação de taper estimada a partir de dados de cubagem que, geralmente, é realizada em árvores-amostra abatidas. Também, em algumas situações, a cubagem pode ser feita com a árvore em pé, como em (1) onde foram cubadas 110 árvores de Anadenanthera peregrina empregando-se um Pentaprisma de Weller.

Independentemente, de abater ou não as árvores-amostra para cubagem, esta atividade acomoda um grande número de dados, como se observa em (2), que acumularam 802 diâmetros medidos ao longo do tronco de 86 árvores de pinus. Também, em (3), que mediram 283 diâmetros do tronco em 35 árvores cubadas de Corymbia citriodora. Já, como em (4), onde coletaram um total 1.695 dados de cubagem de 78 árvores de pinus. Porquanto, em (5) obtiveram 2.412 dados de cubagem de 201 árvores de P. brutia. Nessas situações caracterizadas, se pode empenhar esforços visando reduzir a massa de dados de cubagem sem prejudicar a confiabilidade das equações de taper geradas.

Destaque pode ser dado ao método da altura relativa (1), que possibilita a descrição do perfil do fuste de árvores sem o emprego de uma grande base de dados de cubagem. As medições podem ser coletadas com as árvores em pé, pois exige apenas os diâmetros do fuste situados à 0,1 m, 1,3 m, altura total (H) e em um ponto relativo a H, definido como altura relativa (hr). Aliás, o problema está apenas em hr que pode ser resolvido pelo emprego de um pentaprisma de Weller ou Relascópio de Bitterlich, porque em 0,1 m e 1,3 m a medição do diâmetro do fuste é feita sem dificuldades técnicas operacionais, e em H o mesmo tem valor igual a zero para a espécie em análise.

Se apercebe, então, que o Método da Altura Relativa (Método hr) consiste de uma metodologia promissora para constar em atividades de inventários florestais, porque permite que os diâmetros em 0,1 m, 0,3 m e em hr sejam tomados de forma simultânea às demais variáveis medidas em parcelas de inventário (6). Também, nesta metodologia, o volume do fuste pode ser obtido por meio de uma cubagem simulada com as equações de taper geradas resultando em boa acuracidade, como se observa em (7), (8).

Outra alternativa no desenvolvimento do Método hr é o ajuste de modelos de taper, ou de volumes, empregando-se dados de cubagem simulada com as equações que compõem a metodologia (1), (7). Por exemplo, em (7) se realizou o ajuste de modelos de volumes para avaliar o crescimento inicial de eucalipto e acácia em uma integração de lavoura, pecuária e floresta e também em monocultivo.

Por outro lado, em (1) foi avaliado o emprego de dados de A. peregrina no ajuste do modelo de taper de Kozak (1988), evidenciando acurácia em níveis aceitáveis com emprego tanto de medições diretas do diâmetro do fuste (cubagem de árvores abatidas) quanto de diâmetros obtidos pelo Método hr (cubagem de árvores em pé).

Dentre os trabalhos com enfoque no Método hr, feitos em espécies de eucalipto, acácia, mogno africano e pinus, não se observam trabalhos feitos com dados de Corymbia citriodora. Por isso, dada à importância desta espécie para a produção de diferentes produtos madeireiros (carvão vegetal, lenha, mourões, madeira serrada, dentre outros), se conduziu este trabalho com o objetivo de comparar o ajuste de equações de taper empregando-se dados de cubagem obtidos por meio de medições diretas ao longo do fuste e indiretamente por meio do Método hr.

Material e métodos

O estudo foi realizado utilizando-se dados coletados em um plantio seminal de Corymbia citriodora localizado próximo ao município de Gurupi, estado do Tocantins, Brasil. Este plantio, implantado no espaçamento inicial de 3m X 2 m, tinha a idade de 67 meses na época da coleta de dados que se referem à cubagem feita em 24 árvores abatidas.

Em cada árvore abatida, foram feitas medições diretas do diâmetro ao longo do fuste nas posições: 0,2 m; 0,4 m; 0,7 m; 1,3 m; 2,7 m e, sucessivamente, de 2 m em 2 m até um diâmetro em torno de 3 cm. Também, se mediu a altura total (H) e o diâmetro em hr = (H − 2) / 2 exigido para o desenvolvimento do Método hr que incluiu, ainda, os diâmetros em 0,2 m, 1,3 m e H. Com estes dados, se obteve as seguintes equações de taper (1, 2, 3, 4 e 5):

Em que: d(0,0≤hi≤H)= diâmetro obtido ao longo do fuste empregando-se geometria analítica (cm); CARI, CARII e CARIII= coeficiente angular da reta entre 0,2 a 1,3m, 1,3m a hr e hr a H, respectivamente; hi=i-ésima altura no fuste onde se deseja obter o diâmetro di (m); D=DAP (cm); H=altura total (m); hr1= (H-1.7)/1.7 (m), hr2= (H-1.4)/1.4 (m) e hr3= (H-1.1)/1.1 (m).

Empregando-se as equações de 1 a 4, nas mesmas posições em que se fez medições diretas de diâmetros do fuste com as árvores abatidas (cubagem real), se realizou uma cubagem supondo-se as 24 árvores em pé (cubagem hr). Dessa forma, se obteve as informações necessárias ao ajuste de modelos de taper empregando-se dados da cubagem real e cubagem hr. Para isso, utilizando-se o programa estatístico R For Windows 3.6.2 (9), procedeu-se ao ajuste do modelo de taper de Muhairwe (10) que tem a seguinte conformação eq.5:

Em que: di= diâmetro do tronco na i-ésima altura hi de cubagem (cm); z=hi/Hb, b= coeficientes da regressão à serem estimados; e= erro aleatório do modelo; demais já foram definidos anteriormente.

Para avaliar as equações de taper, geradas com dados de cubagem real e cubagem hr, foi adotada a validação cruzada que consistiu no uso de dois bancos de dados, sendo um reservado para gerar equações e outro para aplicação destas e vice-versa. Nessa condição, conforme feito por (3), adaptou-se as informações da cubagem das 24 árvores utilizando-se, para o primeiro banco de dados (BD1), as posições: 0,2 m; 0,7 m; 4,7 m e a cada 4 m até a última posição no fuste. Já, para o segundo banco de dados (BD2), utilizou-se as posições: 0,4 m; 2,7 m e a cada 4 m até a última posição no fuste que não foi utilizado no BD1.

Empregando BD1 e BD2, em uma validação cruzada, foi possível avaliar a predição dos diâmetros do fuste por meio das equações de taper estimadas a partir do modelo de Muhairwe empregando-se dados de cubagem real e cubagem hr. Neste caso, para avaliar os valores preditos de diâmetros do fuste, foram adotados os critérios estatísticos: resíduo médio (RM), Bias (B), correlação linear entre valores preditos e real (r), desvio padrão dos resíduos (DPR), exatidão de predição obtida pelo teste de qui-quadrado ao nível de 5% de significância (Ex) e raiz quadrada do resíduo médio (RQRM), além da análise gráfica da distribuição dos resíduos (AGR), em cujas formulações e procedimentos de análises podem ser vistos em (3) e (11).

A análise também foi feita na quantificação do volume do fuste das 24 árvores aplicando-se a fórmula de Huber nos dados da cubagem real e cubagem hr. Como a divisão em dois bancos de dados resultaria em baixo número de indivíduos, a análise do volume foi feita com uma validação com dados dependentes adotando-se os mesmos critérios estatísticos da avaliação de equações de taper.

Resultados e discussão

Inicialmente, aplicando-se o método da altura relativa nas 24 árvores cubadas de C. citriodora, se descreveu o perfil do fuste utilizando-se as mesmas posições em que se efetuou às medições diretas nestas mesmas árvores abatidas. Como mostra a Figura 1, há uma tendência de sobreposição da dispersão da cubagem hr com a cubagem real evidenciando a viabilidade de emprego desses dados no ajuste de modelos de taper, como se faz com os dados da cubagem real (3), (11), (12).

Figura 1 Dispersão de diâmetro do fuste de Corymbia citriodora empregando cubagem hr e cubagem real. 

Figure 1. Stem diameter dispersion of Corymbia citriodora using real and hr cubing.

Diante disso, utilizando-se os dados obtidos em BD1 (n=106) e BD2 (n=91) com a cubagem real (CR) e cubagem hr (Chr), se procedeu ao ajuste do modelo de Muhairwe (10) gerando-se duas equações de taper para cada banco de dados (Tabela 1). Em seguida, por meio da aplicação em BD2 das equações geradas por BD1 e vice versa, se obteve uma inexpressiva superioridade das equações geradas por CR, porque os resultados são bem próximos (Tabela 2). Exceção se faz com os critérios RM e B por não resultarem em valores muito próximos entre CR e Chr, porém, se apresentam em níveis aceitáveis de acurácia expressando erros inferiores à ± 2 %.

Tabela 1 Coeficientes estimados do modelo de Muhairwe empregando-se dados obtidos pela cubagem real e cubagem hr. 

Tabela 2 Critérios estatísticos adotados na validação cruzada das equações de taper geradas a partir do ajuste do modelo de Muhairwe empregando-se dados de cubagem real e cubagem hr. 

Quanto à correlação r acima de 0,992, se observam resultados expressando altíssimo índice entre valores reais e preditos como em (3), que resultou em r=0,994 para o modelo em estudo. Valores estes muito parecidos devido à, principalmente, se tratar da mesma espécie apesar de os dados serem coletados em plantios diferentes. Também, este resultado de r>0,992, é melhor ao r=0,942 obtido por (11), sendo o único trabalho que se procedeu igualmente a este estudo, mas utilizando-se o modelo de Kozak (13).

Pelo critério Ex, os diâmetros estão sendo preditos com erros inferiores à 12% com 95% de probabilidade, tendo uma exatidão com valores muito próximos para as equações geradas empregando-se dados de CR e Chr. Valores, também, melhores aos obtidos por (1) com Ex de 22,5% para Chr e de 23,7% para CR. Resultados estes, tanto de r=0,942 como Ex, tiveram qualidade inferior ao presente estudo, provavelmente, por se tratar de uma base de dados de cubagem que se refere à uma espécie florestal amostrada no ambiente de mata nativa (Floresta Atlântica).

Por outro lado, quanto à distribuição dos resíduos (Figura 2), se observa no uso de BD1 que os resíduos se distribuem de forma homogênea e bem equilibrada ao longo do eixo 0 cm de resíduo com ±1,5cm, tanto em relação ao diâmetro predito quanto à altura ao longo do fuste, tendo melhor comportamento para CR. Houve uma tendência em ocorrer resíduos negativos na predição de diâmetros do fuste acima de 15 cm para Chr, sendo tal tendência não detectada para os resíduos em relação de hi que foi igual para ambos CR e Chr.

No uso de dados BD2 (Figura 2), a distribuição de resíduos foi desejável para ambos dados de cubagens CR e Chr com efeito bem menos indesejável com o que foi detectado com BD1. Já, nos histogramas dos resíduos, houve uma tendência de haver distribuição normal destes para ambos dados de cubagens CR e Chr com uma melhor distribuição para CR.

Figura 2 Distribuição de residuos para o modelo de Muhairwe obtido na validação cruzada. R= resíduo, di= diâmetro do fuste, hi= altura no fuste, Hist= histograma, P= predição, Chr e CR= cubagem hr e CR. 

Figure 2. Residual distribution of the Muhairwe model obtained in the cross-validation. R= residue, di= stem diameter, hi= stem height, Hist= histogram, P= prediction, Chr and CR= cubing hr and CR.

As análises feitas com o ajuste do modelo de Muhairwe, evidenciaram haver níveis de acurácia aceitáveis e muito próximos nas equações de taper geradas, tanto com o emprego de CR quanto com Chr. Com isso, se pôde inferir que: usar dados de cubagem obtidos por Cr e Chr no ajuste de uma equação de taper, se pode gerar equações diferentes em níveis aceitáveis de acuracidade. Isso, aliado ao fato de o método da altura relativa possibilitar a coleta de dados com a árvore em pé, se torna bastante vantajoso porque a viabiliza ser feita junto com as medições nas parcelas de inventário (6), o que proporciona redução no tempo e custos na coleta de dados de cubagem.

A avaliação também foi feita aplicando-se a formula de Huber nos dados obtidos por CR e Chr resultando no volume individual das 24 árvores. Se utilizou de uma validação com dados dependentes, porque a separação em dois lotes de dados iria resultar em baixo número de árvores para ajuste e validação. Como resultado, se obteve para Chr um B=0,15% melhor que CR, sendo nos critérios estatísticos r, Ex e RQRM obtido um resultado favorável ao uso de CR com valores de 99,7%, 10,08% e 6,54%, respectivamente (Tabela 3).

Tabela 3 Critérios estatísticos utilizados para avaliar o volume a partir das equações do modelo de Muhairwe ajustado com dados de cubagem real e hr. 

Os resultados obtidos por Chr são próximos aos obtidos pelo emprego de CR (Tabela 3), proporcionando um nível de acurácia na quantificação do volume também em níveis aceitáveis, o que evidencia o mesmo resultado da análise das equações de taper. Além disso, a análise da distribuição de resíduos (Figura 3), comprova os resultados estarem diferentes de forma a não destacar um uso preferencial por CR. Se observa no histograma dos resíduos, uma amplitude entre ±15% para CR e uma amplitude de ±10% para Chr. Portanto, se pode inferir que é possível também utilizar o método da altura relativa para obter o volume individual de árvores de C. citriodora, com ressalvas de não haver a necessidade de abater árvores-amostra para cubagem.

Figura 3 Distribuição de resíduos do volume obtido pelas equações do modelo de Muhairwe empregando-se a cubagem hr (Chr) e real (CR). R= resíduo, D=diâmetro do tronco em 1,3 m (DAP), Hist= histograma. 

Figure 3. Residual distribution of the volume obtained by the Muhairwe model using the cubing hr (Chr) and real (CR). R= residual, D= stem diameter in 1.3 m (DHB), Hist= histogram.

Conclusão

De acordo com os resultados obtidos neste estudo, concluiu-se que o ajuste de um modelo de taper, empregando-se dados de cubagem simulados pelo Método hr, resulta em um nível de acurácia aceitável, inexpressivamente inferior ao obtido com o uso de dados reais de cubagem de C. citriodora, inclusive na quantificação do seu volume. Também, que o emprego do Método hr, possibilitará significativa redução de tempo e custos dos procedimentos técnico-operacionais envolvidos na coleta de dados, uma vez que exige a medição de diâmetros do fuste em apenas quatro posições sem o abate de árvores-amostra, possibilitando que sejam feitas em simultâneo com as medições das parcelas de inventário.

Agradecimentos

À Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) pela bolsa de pesquisa.

Referências

1. V.C.L. Andrade, et al., ''Perfil do fuste de Anadenanthera peregrina descrito pelo método da altura relativa”, Ciência Florestal, vol.26, pp.875-887, 2016. [ Links ]

2. D. Lee, et al., ''Estimation and Validation of Taper Equations for Three Major Coniferous Species in Gangwon and North Gyeongsang Provinces of South Korea”. Journal of Forest and Environmental Science, vol.33, no.4, pp.315-321, 2017. [ Links ]

3. V.C.L. Andrade, et al., ''Equações de perfil do tronco avaliados para Corymbia citriodora Hill & Johnson”, Agrarian, vol.14, no.51, pp.82-94, 2021. [ Links ]

4. D.H. Pérez, et al., ''Funciones de ahusamiento y volumen comercial para Pinus Patula Schltdl. Et Cham. en Zacualtipán, Hidalgo”. Revista Mexicana de Ciencias Forestales, vol.4, no.16, pp.34-45, 2018. [ Links ]

5. S. De Miguel, et al., ''Evaluating marginal and conditional predictions of taper models in the absence of calibration data”. Canadian Journal Forest. Research. no.42, pp.1383-1394, 2012. [ Links ]

6. H.G. Leite, V.C.L. Andrade, ''Uso do método da altura relativa em inventário florestal de um povoamento de Pinus”. Revista Arvore, vol.28, no.6, pp.865-873, 2004. [ Links ]

7. F.L.R. Oliveira, et al., ''Crescimento inicial de eucalipto e acácia, em diferentes arranjos de integração lavoura-pecuária-floresta”. Cerne, vol.21, no.2, pp.227-233, 2015. [ Links ]

8. R.C. Gama, V.C.L. Andrade, ''Volumetria de Khaya senegalensis pelo método da altura relativa”. Advance Forestry Science, vol.6, no.3, pp.737-742, 2019. [ Links ]

9. R development core team, ''R: a language and environment for statistical computing”. Vienna: R Foundation for Statistical Computing. 2019. [ Links ]

10. C.K. Muhairwe, ''Taper equations for Eucalyptus pilularis and Eucalyptus grandis for the north coast in New South Wales, Australia”. Forest Ecology and Management, Vol.113, no.2-3, pp.251-269, 1999. [ Links ]

11. H.A. Beltran, et al., ''Volume and taper equations for commercial stems of Nothofagus obliqua and N. alpina”. Cerne, vol.23, no.3, pp.299-309, 2017. [ Links ]

12. K.J. Farias, et al., ''Quantificação do volume total por meio de modelos volumétricos e de afilamento em Eucalyptus grandis W. Hill ex Maiden”. Revista de Ciências Agroveterinárias, vol.18, no.1, pp.33-40, 2019. [ Links ]

13. A. Kozak. ''A variable-exponent taper equation”. Canadian Journal Forest. Research, vol.18, no.11, pp.1363-1368, 1988 [ Links ]

Recebido: 02 de Novembro de 2022; Aceito: 22 de Fevereiro de 2023

Creative Commons License Este é um artigo publicado em acesso aberto sob uma licença Creative Commons