SciELO - Scientific Electronic Library Online

 
vol.20 issue47Taper equations adjusted for Corimbya citriodora Hill & Johnson using data generated by the relative height method author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista Forestal Mesoamericana Kurú

On-line version ISSN 2215-2504

Kurú vol.20 n.47 Cartago Jul./Dec. 2023

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

Artículo

Alternativas de modelagem para a estimativa da altura de árvores de eucalipto

Modeling alternatives for height estimation in eucalypt tree species

Bruno Oliveira Lafetá1 

Ivelton Alves Pimenta1 

Milene Alves dos Santos1 

Heloisa Brenda1 

Xavier Rodrigues1 

Ivan da Costa Ilhéu Fontan1 

Ana Carolina Ferraro1 

Diego dos Santos Vieira2 

1. Departamento de Engenharia Florestal, Instituto Federal de Minas Gerais; São João Evangelista, Brasil; bruno.lafeta@ifmg.edu.br, iveltonalves@gmail.com, milenesantoskp@gmail.com, heloisabrenda.eng.florestal@gmail.com, ivan.fontan@ifmg.edu.br ana.ferraro@ifmg.edu.br

2. Departamento de Engenharia Florestal, Universidade Federal dos Vales do Jequitinhonha e Mucuri; Diamantina, Brasil; diegovieir4@gmail.com

Resumo

O objetivo do presente trabalho foi avaliar alternativas para a modelagem hipsométrica generalista de eucalipto via regressão linear, com a inclusão de índices de competição independentes da distância como variável preditora. Os dados utilizados foram provenientes de 34 parcelas distribuídas em quatro unidades de manejo florestal com 72 meses de idade. Calcularam-se 7 índices de competição independentes da distância. Avaliou-se o desempenho preditivo de 16 modelos hipsométricos generalistas, sendo propostos 14 modelos de dupla entrada (DAP e índice de competição). Todas as equações geradas foram biologicamente consistentes. Diante da ausência de informações da altura de árvores dominantes, a modelagem generalista pode ser melhor viabilizada com a inclusão do índice IC4, que representa a área basal de árvores vizinhas competidoras. Conclui-se que a inclusão de índices de competição em modelos hipsométricos ampliam a capacidade de generalização de equações para talhões de eucalipto com distintas capacidades produtivas. O índice de competição IC4 favorece a qualidade de ajuste e desempenho preditivo da modelagem hipsométrica nos sítios estudados.

Palavras chave: Competição; generalização; relação hipsométrica; regressão linear

Abstract

This study aimed to evaluate alternatives for general hypsometric modeling of eucalypt through linear regression, with inclusion of distance-independent competition indices as a predictor variable. The data was collected from 34 plots distributed in four 72-month-old forest management units. Seven distance-independent competition indices were calculated. The study evaluated the predictive performance of 16 general hypsometric models, and 14 double-entry models (DAP and competition index) were proposed. All equations generated were biologically consistent. Given the lack of information on the height of dominant trees, general modeling can be better employed with inclusion of the IC4 index, which represents basal area of neighboring competing trees. It is concluded that the inclusion of competition indices in hypsometric models increases generalization capacity of equations for eucalyptus stands with different productive capacities. The IC4 competition index enhances fit quality and predictive performance of the hypsometric modeling in the study sites.

Keywords: Competition; generalization; hypsometric relationship; linear regression

Introdução

O detalhamento e a qualidade de informações biométricas são essenciais para um adequado planejamento e prognóstico do crescimento e produtividade florestal. A busca por alternativas que forneçam estimativas exatas e reduzem o tempo das medições em campo são constantes, exigindo-se de mais pesquisas sobre o processamento de dados biométricos (1).

A altura total de uma árvore se refere à distância linear do nível do solo até o ápice da copa. Trata-se de uma variável biométrica amplamente utilizada em cálculos volumétricos de fustes e de seus múltiplos sortimentos, classificação produtiva de sítios, caracterização de crescimento e orçamentos de carbono (2). Por outro lado, sua medição é mais dispendiosa e demorada em comparação com aquela do diâmetro. A fim de se otimizar a mensuração em povoamentos equiâneos, é comum medir a altura de uma parte das árvores da parcela e adotar a modelagem hipsométrica para a estimativa de altura dos demais indivíduos, que tiveram somente o diâmetro medido.

A modelagem hipsométrica por sua vez, é uma técnica para o estabelecimento de relação funcional com vistas à estimativa da altura individual de árvores, processada frequentemente via regressão linear. A estimação da variável mais laboriosa em função daquela correlacionada de mais fácil e rápida medição, permite uma sensível redução nos custos e melhoria operacional de amostragem (3).

Os modelos hipsométricos podem ser classificados em duas categorias distintas de acordo com a aplicabilidade, tradicional (local) ou generalistas (geral). Modelos tradicionais são os mais simples, dependentes apenas do diâmetro e indicados para o ajustamento por talhão ou para um pequeno conjunto de talhões análogos e homogêneos em termos de capacidade produtiva. A forte correlação entre a altura e o diâmetro justifica o uso desta como variável preditora em modelos tradicionais. Contudo, diversos aspectos estão relacionados com a expressão do crescimento vegetal. É preciso a compreensão de que as árvores se encontram em contínuo estado de competição em um povoamento florestal, sendo o crescimento influenciado por aspectos genéticos, edáficos, fisiográficos, climáticos e demais mecanismos relacionados à complexidade de ecossistemas (4), (5).

A inclusão de características do povoamento aumenta a qualidade preditiva e realismo biológico dos modelos hipsométricos, tornando a equação adaptável a vários sítios (1), (6), (7), (8). O modelo generalista mais conhecido no setor florestal brasileiro é o de (9), que se baseia no uso de dados diamétricos e logaritmizados da altura dominante como variáveis preditoras. Todavia, a coleta de dados é um entrave para o uso desse modelo, que requer a identificação e medição da altura de mais árvores em campo, especificadamente, daquelas dominantes.

Na busca de alternativas mais eficientes para a estimativa da altura de árvores, adaptações em modelos hipsométricos clássicos podem ser adotadas para a melhoria do desempenho da modelagem sem intensificar a mensuração em campo. Apesar de incipiente, a inclusão de variáveis preditoras obtidas em nível de parcela tem sido pesquisada em substituição a altura dominante, como o diâmetro médio quadrático, diâmetro máximo, diâmetro dominante, média aritmética do diâmetro e outras (2), (10), (11). Outra alternativa promissora, é a inclusão de informações competitivas individuais da própria árvore como variável preditora para a estimativa de sua altura.

A competitividade entre árvores pode ser resumida na forma de índices ou expressões matemáticas que refletem o grau de supressão/dominância de cada árvore em relação àquelas competidoras (12), (13), (14), (15), (16). Os índices de competição independentes da distância é uma representação numérica do efeito cumulativo de todo o conjunto de características e interações bióticas e abióticas inerentes a uma árvore ou grupo de indivíduos ao longo do tempo, desconsiderando a posição espacial relativa de árvores dentro da unidade amostral (15), (17), (18).

Embora o uso de índices de competição seja comum em modelagens do crescimento e produção em nível de árvore individual (19), (20), investigações científicas sobre sua aplicação na modelagem hipsométrica são escassas e podem subsidiar a quantificação de recursos florestais. Mediante exposto, o objetivo do presente trabalho foi avaliar alternativas para a modelagem hipsométrica generalista de eucalipto via regressão linear, com a inclusão de índices de competição independentes da distância como variáveis preditoras.

Material e métodos

O trabalho foi realizado com informações de 34 parcelas distribuídas aleatoriamente em povoamentos clonais de Eucalyptus urophylla S. T. Blake pertencentes a uma empresa florestal, localizados na região central de Minas Gerais. O clima da região é do tipo Cwa pelo sistema internacional de classificação de Köppen (21). Foram estabelecidas quatros unidades de manejo florestal (UMF) com médias de temperatura e precipitação de: 22,7ºC e 1.259 mm para a UMF1 (87,91 ha); de 23,4ºC e 1.043 mm para a UMF2 (42,74 ha); de 23,4ºC e 1.043 mm para UMF3 (77,21 ha) e; de 23,0ºC e 1.189mm par UMF4 (85,57 ha), respectivamente.

Realizou-se o inventário florestal na idade de 72 meses, sendo lançadas: 11 parcelas na UMF1; 6 parcelas na UMF2; 10 parcelas na UMF3 e 7 parcelas na UMF4. A área útil de cada parcela e o espaçamento de plantio foram 400 m2 e 3,0 × 2,5 m, respectivamente. Foram mensurados o diâmetro com casca à altura do peito (DAP - à altura de 1,30 m do solo, cm) de todos os fustes. Em cada parcela, tomaram-se a altura total (Ht, m) das 15 primeiras árvores e de quatro árvores dominantes (Hd, m), conforme conceito de (22), equivalente a uma árvore dominante a cada 100 m2. Assumindo a idade índice de 72 meses, a capacidade produtiva de UMF1 a UMF4 foram de: S = 29,26 m; S = 29,81 m; S = 24,94 m e; S = 26,12 m, respectivamente.

Calcularam-se sete índices de competição independentes da distância conhecidos no setor florestal (Quadro 1). Os índices IC1, IC2 e IC3 representam a competição bilateral e assumem que os recursos de crescimento são disputados por todas as árvores vizinhas. Os demais índices (IC4 a IC7) assumem como árvores competidoras aquelas de diâmetro maior que a árvore objetivo e, portanto, representam a competição unilateral, ou seja, árvores maiores não são afetadas por aquelas vizinhas menores.

Quadro 1 Índices de competição calculados para as árvores de eucalipto. 

Os dados foram divididos aleatoriamente em grupos de parametrização (70,0 % das amostras) e de validação (30,0 % das amostras), mutualmente exclusivos, empregando método randômico de amostragem sem reposição (método holdout). Nessa divisão, todas as parcelas foram contempladas com árvores em ambos os grupos (Quadro 2).

Quadro 2 Descrição de dados utilizados no ajuste hipsométrico de eucalipto aos 72 meses de idade, separados em grupos de parametrização e validação. 

Ajustaram-se 16 modelos hipsométricos generalistas (Quadro 3), sendo dois modelos rotineiramente adotados no setor florestal baseados em informações de fustes dominantes e 14 propostos, de dupla entrada (DAP e índice de competição). Os modelos propostos foram estabelecidos a partir dos modelos tradicionais de (23) e de (24), que possuem saída hipsométrica logaritmizada e se diferem em termos de variável preditora; o primeiro inclui o inverso de DAP e o segundo, o logaritmo neperiano de DAP. Todos os modelos foram ajustados a todo conjunto de dados do grupo de parametrização, isto é, sem estratificação. O método dos Mínimos Quadrados Ordinários (MQO) foi adotado nas análises de regressão.

Quadro 3 Modelos hipsométricos generalistas testados para a estimativa da altura de eucalipto aos 72 meses de idade. 

A fim de se avaliar a presença de multicolinearidade dos modelos, calculou-se o Valor de Inflação da Variância (VIF). Os pontos que extrapolaram a tendência geral dos dados não foram eliminados das análises estatísticas a fim de se verificar a capacdade da modelagem em lidar com outliers ou ruídos. Para diagnóstico de efeito estatístico, empregou-se 5% de significância em todas as análises. Estas foram efetuadas com auxílio do software R versão 3.5.2 (25).

Resultados

Os coeficientes e a qualidade de ajuste dos modelos hipsométricos generalistas para eucalipto aos 72 meses de idade se encontram no Quadro 4. O valor de inflação da variância de todos os ajustes realizados foi baixo (VIF < 9). Todos os coeficientes das equações foram significativos (p ≤ 0,05), exceto para a equação (10). Esta equação foi descartada da análise gráfica subsequente. Em geral, as equações apresentaram poucos desvios (baixos valores de RQEM e MDA) e altas correlações (r ≥ 0,90 na fase de parametrização e r ≥ 0,89, na validação). Tomando-se como referência a fase de parametrização, a RQEM e MDA aumentaram em torno de 1,27 a 8,38% e de 0,35 a 8,64% na fase de validação. Notou-se pequeno aumento do coeficiente de correlação na fase de validação somente para a equação (2), de 0,06%.

Quadro 4 Coeficientes e qualidade de ajuste dos modelos hipsométricos generalistas para eucalipto aos 72 meses de idade. 

No que se refere às estatísticas de qualidade de ajuste, a equação (2) tendeu a oscilar menos entre fases de parametrização e validação. O modelo de (9) ajustado, ou modelo (1), exibiu as maiores correlações e os menores valores de RQEM e MDA nas fases de parametrização e validação.

Dentre os modelos propostos (Quadro 3), as equações (13) e (14) tiveram os menores valores de RQEM (≤ 1,15) e MDA (≤ 0,93), além da maior correlação (≥ 0,93), durante parametrização e validação. Estas duas equações se basearam no modelo de (24) e a inclusão dos índices de competição IC4 e IC5 como preditores, respectivamente. Uma característica em comum destes índices é a inclusão do atributo área basal das árvores vizinhas competidoras na expressão matemática. Ambas as equações foram selecionadas para as análises de inspeção gráfica subsequentes (Figura 1).

Figura 1 Distribuição de resíduos percentuais das equações hipsométricas generalistas obtidas para eucalipto aos 72 meses de idade. 

Figure 1. Distribution of percentage residuals of general hypsometric equations obtained for eucalyptus at 72 months of age.

Apenas para fins de visualização gráfica, foram omitidos os resíduos percentuais da equação (14) referentes aos valores de 46,55% (DAP = 3,44 e IC5 = 0,24) e 47,39% (DAP = 3,82 e IC5 = 0,23) nas fases de parametrização e validação, respectivamente. A distribuição dos resíduos percentuais das equações selecionadas tendeu à homocedasticidade, exceto na equação (14). As equações (1) e (13) exibiram resíduos mais próximos ao eixo das abcissas, especialmente, no intervalo entre -10 a 10%.

Discussão

Todas as equações hipsométricas geradas foram biologicamente consistentes, com coeficientes negativos associados ao inverso de DAP e positivos para as demais relações diretas com atributos biométricos (Quadro 4). É conveniente salientar que as estatísticas apresentadas neste quadro não retratam o desempenho dos modelos hipsométricos, tratam-se apenas, das médias das variações das estimativas totais ou indicativo da qualidade de ajuste (26). Portanto, fez-se necessária a validação e análise gráfica de resíduos.

A inclusão da altura dominante na modelagem hipsométrica tem sido exaustivamente utilizada para ampliar o potencial de generalização das equações até talhões com diferentes potenciais produtivos. O modelo de (9) ajustado apresentou a melhor qualidade de ajuste e desempenho estatístico para representar a relação hipsométrica de E. urophylla (Quadro 4 e Figura 1). Enfatiza-se que este modelo se difere do modelo de (23) pela inclusão da altura dominante em sua expressão e possui, como limitação, a necessidade de informações de árvores dominantes durante a amostragem, tornando o inventário florestal mais dispendioso e laborioso.

Por outro lado, a inclusão de índice competição como variável preditora, independentes da distância e baseados em aspectos da estrutura horizontal, influenciou o desempenho preditivo dos modelos ajustados. Nas circunstâncias cuja amostragem não tenha abrangida a identificação e/ou coleta de informação hipsométrica de árvores dominantes, índices de competição podem ser analisados para a definição de equações generalistas para a estimativa da altura. Verificou-se que a inclusão de índices de competição favorece, em algumas ocasiões, a qualidade de ajustamento, reduzindo os desvios e ampliando o desempenho preditivo da modelagem generalista. Esta característica possui grande importância prática no processamento de dados, pois facilita a aplicação da equação hipsométrica para amplas áreas de cultivo e diminui a quantidade de estratificações e ajustes hipsométricos.

Os modelos generalistas propostos se fundamentaram em características competitivas individuais das próprias árvores. Embora a competição seja um fenômeno de alta complexidade biológica (12), (14), (27) e tenha sido sintetizada na forma de índice, a precisão do ajuste desses modelos foi considerada satisfatória. As equações provenientes dos modelos (11) a (14), baseadas na expressão do modelo de (24), exibiram qualidade de ajuste um pouco melhor em relação àquela que se baseou na inclusão do diâmetro dominante como variável preditora.

Apesar de que o ajuste do modelo (10) tenha resultado em coeficiente de regressão não significativo e, em comparação aos demais ajustes, poucos desvios ((Quadro 4) não se verificou a presença de multicolinearidade, haja vista que, VIF acima de 10 é um indício claro de multicolinearidade (28). Ressalta-se que diante de eventual presença da multicolinearidade, os coeficientes da regressão deixam de refletir os efeitos inerentes de particular variável independente sobre aquela dependente, representando parcialmente seu real efeito. Também, não se observou problemas de sub ou sobreparametrização nas equações obtidas, o que levaria a discrepâncias na qualidade de ajuste entre as fases de validação e parametrização. A validação é uma importante etapa para a escolha de um modelo hipsométrico livre de vieses e resíduos heterogêneos.

O uso indiscriminado de modelos hipsométricos deve ser visto com bastante cautela, pois dificuldades relacionadas a inconsistência biológica são comuns durante o processamento estatístico. A escolha apropriada de um modelo hipsométrico permite a obtenção de estimativas não enviesadas de altura e, consequentemente, do volume. A inclusão do índice de competição IC4 como variável preditora em modelos hipsométricos mostrou que pode ser uma alternativa viável em substituição ao uso de informações biométricas de árvores dominantes, reduzindo soma de quadrado dos resíduos e a intensidade de fustes amostrados para a modelagem hipsométrica generalista em povoamento equiâneos. Este índice representou a competição unilateral, aquela que implica em concorrência especialmente por radiação luminosa (29), que é comum em povoamentos cuja competição já se encontra estabelecida, em idades próximas a colheita florestal.

A modelagem hipsométrica é uma etapa fundamental para a estimativa da altura de árvores e quantificação volumétrica em inventários florestais de povoamentos equiâneos. Os resultados obtidos fornecem subsídios importantes para o desenvolvimento de futuras pesquisas sobre a hipsometria em povoamentos florestais empregando, particularmente, modelos de regressão lineares, facilmente ajustados e interpretados. A inclusão de preditores que representam aspectos relacionados as interações bióticas e abióticas na modelagem hipsométrica, como índices de competição, mostrou potencial para a melhoria preditiva sem comprometer a amostragem em campo e o realismo biológico.

Conclusões

A inclusão de índices de competição independentes da distância em modelos hipsométricos ampliam a capacidade de generalização de equações para talhões de eucalipto com distintas capacidades produtivas. O índice de competição IC4, que representa à área basal de árvores vizinhas competidoras, favorece a qualidade de ajuste e desempenho preditivo da modelagem hipsométrica nos sítios estudados. Diante da ausência de informações da altura de árvores dominantes, a modelagem generalista pode ser melhor viabilizada com a inclusão desse índice como variável preditora no modelo de (24).

Agradecimentos

Ao Instituto Federal de Educação, Ciência e Tecnologia (IFMG) - Campus de São João Evangelista-MG por todo apoio logístico para a realização do presente trabalho.

Referências

1. D. Dantas, et al. ''Reduction of sampling intensity in forest inventories to estimate the total height of eucalyptus trees”, Bosque, vol. 41, no. 3, pp. 353-364, 2020. [ Links ]

2. A. C. Ferraz Filho, et al. ''Height-diameter models for Eucalyptus sp. plantations in Brazil”, Cerne, vol. 24, no. 1, pp. 9-17, 2018. [ Links ]

3. M. N. I. Sanquetta, et al. ''Ajuste de equações hipsométricas para a estimação da altura total de indivíduos jovens de teca”, Científica, vol. 43, no. 4, pp. 400-406, 2015. [ Links ]

4. S. Vatraz, et al. ''Índices de competição dependentes da distância do estrato arbóreo na Amazônia brasileira”, Espacios, vol. 37, no. 27, pp. 1-11, 2016. [ Links ]

5. S. Versace, et al. ''Prediction of competition indices in a Norway spruce and silver fir-dominated forest using LiDAR data”, Remote Sensing, vol. 11, pp. 18, 2019. [ Links ]

6. D. H. B. Binoti, et al. ''Estimation of height of eucalyptus trees with neuroevolution of augmenting topologies (NEAT)”, Revista Árvore, vol. 41, no. 3, e410314, 2017. [ Links ]

7. T. Schmitt e V. C. L. Andrade. '' Identidade de modelos hipsométricos para um plantio clonal”, Advances in Forestry Science, vol. 6, no. 2, pp. 583-587, 2019. [ Links ]

8. P. R. S. Fernandes, et al. ''Equações hipsométricas para Pinus sp. na estação experimental de ciências florestais de Itatinga, São Paulo”, Journal of Biotechnology and Biodiversity, vol. 9, no. 1, pp. 96-106, 2021. [ Links ]

9. .J .C. C. Campos, et al. Inventário florestal nacional, reflorestamento: Minas Gerais. Brasília, DF: IBDF, 1984. [ Links ]

10. M. Tomé, et al. ''Relação hipsométrica geral para Eucalyptus globulus Labill. em Portugal”, Silva Lusitana, vol. 15, no. 1, pp. 41-55, 2007. [ Links ]

11. S. J. Téo e T. C. Silva. ''General height-diameter equation depending on the stand variables, for Eucalyptus benthamii”, Floresta e Ambiente, vol. 27, no. 4, e20180302, 2020. [ Links ]

12. S. Hill, et al. ''Individual-tree and stand-based development following natural disturbance in a heterogeneously structure forest: a LiDAR-based approach”, Ecological Informatics, vol. 38, pp. 12-25, 2016. [ Links ]

13. E. A. Costa , et al .''Competition indices and their relationship with basal area increment of Araucaria”, Journal of Agricultural Science, vol. 10, no. 5, pp. 198-2010, 2018. [ Links ]

14. G. Quiñonez-Barraza, et al. ''Considering neighborhood effects improve individual dbh growth models for natural mixed-species forest in Mexico”, Annals of Forest Science, vol. 75, no. 78, pp. 1-11, 2018. [ Links ]

15. Y. Yang e S. Huang. ''Effects of competition and climate variables on modelling height to live crown for tree boreal tree species in Alberta, Canada”, European Journal of Forest Research, vol. 137, pp. 153-167, 2018. [ Links ]

16. M. Minatti, et al. ''Índices de competição dependentes da distância para Araucaria angustifólia em São João do Triunfo - PR”, Brazilian Journal of Biosystems Engineering, vol. 13, no. 3, pp. 187-194, 2019. [ Links ]

17. J. Tenzin, et al. ''Individual tree basal area increment models for broadleaved forests in Bhutan”, Forestry, vol. 90, pp. 367-380, 2017. [ Links ]

18. I. M. Lustosa Júnior, et al. ''Competition indexes to evaluate tree growth in a semi-deciduous seasonal forest”, Floresta e Ambiente , vol. 26, no. 4, pp. 1-12, 2019. [ Links ]

19. G. C. Vieira, et al. ''Prognoses of diameter and height of trees of eucalyptus using artificial intelligence”, Science of the Total Environment, vol. 619-620, pp. 1473-1481, 2018. [ Links ]

20. J. V. N. Carrijo, et al. ''The growth and production modeling of individual trees of Eucalyptus urophylla plantations”, Journal of Forestry Research, vol. 31, no. 5, pp. 1663-1672, 2020. [ Links ]

21. C. A. Alvarez et al. ''Köppen’s climate classification map for Brazil”, Meteorologische Zeitschrift, vol. 22, no. 6, pp. 711-728, 2013. [ Links ]

22. E. Assmann. The principles of forest yield study. New York: Pergamon, 1970. [ Links ]

23. R. O. Curtis. ''Height, diameter and height diameter age equations for second growth Douglas-fir”, Forest Science, vol. 13, no. 4, pp. 365-375, 1967. [ Links ]

24. A. Stoffels e J. V. Soest. ''The main problems in sample plots”, Bus Natlerbouth, vol. 25, pp. 190-199, 1953. [ Links ]

25. R Core Team. ''R: A language and environment for statistical computing”, Vienna: R Foundation for Statistical Computing. 2018. [ Links ]

26. R. R. Souza, et al. ''Forma de fuste de árvores de Eucalyptus em plantios com diferentes densidades iniciais”, Sciencia Forestalis, vol. 44, no. 109, pp. 33-40, 2016. [ Links ]

27. R, P. Sharma e A. Brunner. ''Modeling individual tree height growth of Norway spruce and Scots pine from national forest inventory data in Norway”, Scandinavian Journal of Forest Research, vol. 32, no. 6, pp. 501-514, 2016. [ Links ]

28. D. N. Gujarati. Econometria básica. São Paulo, SP: Pearson Makron Books, 2000. [ Links ]

29. C. T. Schons, et al. ''Índices de competição em nível de árvore individual para espécies da floresta ombrófila mista”, Scientia Forestalis, vol. 48, no. 128, e3275, 2020 [ Links ]

Recebido: 30 de Setembro de 2022; Aceito: 27 de Junho de 2023

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