Resumos
Os objetivos deste estudo foram ajustar, empregar e propor modelos de distribuição diamétrica com diferentes formas da função hiperbólica para a projeção da distribuição diamétrica de povoamentos desbastados. Foram utilizados dados de um experimento sobre desbastes, instalado em três povoamentos de um híbrido Eucalyptus grandis x Eucalyptus urophylla, localizado na região Nordeste do Estado da Bahia. As parcelas permanentes foram mensuradas aos 27, 40, 50, 58, 61, 76, 87, 101 e 112, 125, 135, 147 e 160 meses, com um desbaste realizado aos 58 meses e um segundo aos 130 meses. A utilização da função com dois parâmetros apresentou melhores resultados de raiz quadrada do erro médio (RMSE), bias e bias%.
Função hiperbólica truncada; Desbaste; Multiprodutos
The objectives of this study were to adjust, to employ and to propose diameter distribution models with different hyperbolic function shapes for the projection of the diameter distribution of thinned stands. It was used data from an experiment on thinning set up in three stands of an Eucalyptus grandis × Eucalyptus urophylla hybrid, located in the northeastern Bahia state. The permanent plots were measured at 27, 40, 50, 58, 61, 76, 87, 101 and 112, 125, 135, 147 and 160 months, with a thinning carried out at 58 months of age and a second thinning conducted at 130 months of age. The two-parameter function showed the best results of root mean square error (RMSE), bias and bias%.
Truncated hyperbolic function; Thinning; Multiproducts
Modelos de distribuição diamétrica para povoamentos desbastados de eucalipto utilizando a função hiperbólica
Diametric distribution model thinned eucalypt stands for using a hyperbolic function
Daniel Henrique Breda BinotiI; Mayra Luiza Marques da Silva BinotiII; Helio Garcia LeiteIII, Silvana Lages Ribeiro GarciaIV; Jovane Pereira da CruzV
IPrograma de Pós-Graduação em Ciência Florestal pela Universidade Federal de Viçosa, UFV, MG, E-mail: danielbinoti@yahoo.com.br
IIDepartamento de Engenharia Florestal, Universidade Federal dos Vales do Jequitinhonha e Mucuri, UFVJM , MG, E-mail: mayra_ufv@yahoo.com.br IIIDepartamento de Engenharia Florestal, Universidade Federal de Viçosa, UFV, MG, E-mail: hgleite@gmail.com
IVFaculdade de Viçosa, FDV, Brasil, E-mail: garciaslr@yahoo.com.br
VBahia Pulp, E-mail: jovane_cruz@bahiapulp.com
RESUMO
Os objetivos deste estudo foram ajustar, empregar e propor modelos de distribuição diamétrica com diferentes formas da função hiperbólica para a projeção da distribuição diamétrica de povoamentos desbastados. Foram utilizados dados de um experimento sobre desbastes, instalado em três povoamentos de um híbrido Eucalyptus grandis x Eucalyptus urophylla, localizado na região Nordeste do Estado da Bahia. As parcelas permanentes foram mensuradas aos 27, 40, 50, 58, 61, 76, 87, 101 e 112, 125, 135, 147 e 160 meses, com um desbaste realizado aos 58 meses e um segundo aos 130 meses. A utilização da função com dois parâmetros apresentou melhores resultados de raiz quadrada do erro médio (RMSE), bias e bias%.
Palavras-chave: Função hiperbólica truncada, Desbaste e Multiprodutos.
ABSTRACT
The objectives of this study were to adjust, to employ and to propose diameter distribution models with different hyperbolic function shapes for the projection of the diameter distribution of thinned stands. It was used data from an experiment on thinning set up in three stands of an Eucalyptus grandis × Eucalyptus urophylla hybrid, located in the northeastern Bahia state. The permanent plots were measured at 27, 40, 50, 58, 61, 76, 87, 101 and 112, 125, 135, 147 and 160 months, with a thinning carried out at 58 months of age and a second thinning conducted at 130 months of age. The two-parameter function showed the best results of root mean square error (RMSE), bias and bias%.
Keywords: Truncated hyperbolic function, Thinning and Multiproducts.
1. INTRODUÇÃO
Os povoamentos florestais destinados à produção de multiprodutos são caracterizados pelo longo período de tempo entre os investimentos iniciais e os recebimentos futuros em função da necessidade de árvores de maiores dimensões, que na maioria das vezes exige a aplicação de desbastes e desrama artificial. Nesse sentido, a predição ou projeção do crescimento e da produção dos povoamentos florestais em nível de árvore ou de classe de diâmetro é essencial para a definição da quantidade dos diversos produtos que serão produzidos pelo povoamento (GUIMARÃES, 1994; DEMOLINARI et al., 2007).
Modelos de distribuição diamétrica (MDD) estimam o crescimento e a produção por classe de diâmetro, permitindo a análise econômica de multiprodutos e a simulação de desbastes (BURKHART et al., 1981; CLUTTER et al., 1983; NOGUEIRA et al., 2005). A eficiência desses modelos é influenciada pelas características da função de densidade probabilidade utilizada. A capacidade de correlação dos parâmetros da função com características do povoamento é fator limitante para a construção das equações que compõem o MDD.
Estudos sobre modelagem da distribuição diamétrica destacam o emprego da função hiperbólica (GUIMARÃES, 2002; CAMPOS; LEITE, 2009; LEITE et al., 2010). Essa função pode, em sua formulação completa, ser descrita como:
em que α é o parâmetro de locação, βo parâmetro de escala (β > 0), γ o parâmetro de forma ( γ > 0) e x o centro de classe de diâmetro (x > 0).
Estudos de modelagem utilizando a função hiperbólica para povoamentos submetidos a desbastes ainda são escassos. Por isso, objetivou-se neste estudo ajustar e empregar modelos de distribuição diamétrica com diferentes formas da função hiperbólica para a projeção da distribuição diamétrica desse tipo de povoamento.
2. MATERIAL E MÉTODOS
2.1. Dados
Foram utilizados dados de um experimento sobre desbastes, instalado em três povoamentos de um híbrido Eucalyptus grandis x Eucalyptus urophylla (Rio Claro), em uma empresa florestal localizada na região Nordeste do Estado da Bahia, Brasil. O experimento foi instalado em três instalações, compostos por quatro tratamentos e seis blocos, sendo dois destes em cada instalação. Cada bloco continha duas repetições, totalizando 48 parcelas, com área útil de 2.600 m2, sendo o espaçamento inicial entre as árvores de 3,0 x 3,0 m. Os tratamentos aplicados corresponderam a 20%, 35% e 50% de área basal removida em cada desbaste, eliminando-se os piores indivíduos em função do crescimento volumétrico.
As parcelas permanentes foram mensuradas aos 27, 40, 50, 58, 61, 76, 87, 101 e 112, 125, 135, 147 e 160 meses, com um desbaste realizado aos 58 meses e um segundo aos 130 meses. Em cada medição, foram mensurados os diâmetros a 1,3 de altura (dap), de todas as árvores da parcela.
2.2. Descrição das Funções
A fdp da função hiperbólica de três parâmetros (3P) pode ser descrita como:
em que α é o parâmetro de locação, βo parâmetro de escala (β > 0), γ o parâmetro de forma (γ > 0) e x o centro de classe de diâmetro (x > 0). A supressão do parâmetro de locação da função caracteriza a função hiperbólica de dois parâmetros (2P), conforme apresentado na equação 2.
O parâmetro de locação pode ser substituído pelo diâmetro mínimo do povoamento (α = dmin), sendo ajustado da seguinte forma:
A função hiperbólica pode ser truncada à direita, em função do diâmetro máximo do povoamento, sendo descrita da seguinte forma:
em que T é o diâmetro máximo da parcela. Para realização dos ajustes, os dados de distribuição diamétrica observada em cada parcela em todas as medições foram agrupados em classes com amplitude do 1,0 cm. Testou-se a aderência das funções aos dados pelo teste de Kolmogorov-Smirnorv (KS) no nível de 1% de significância (GIBBONS; SUBHABRATA, 1992).
2.3. Modelos de distribuição diamétrica
A redistribuição teórica dos diâmetros por classe foi realizada por meio da alternativa tradicional de correlacionar parâmetros da função com características do povoamento em idades atuais ou futuras, conforme feito por Nogueira et al. (2005). Consideraram-se como variáveis dependentes os parâmetros da função hiperbólica nas quatro formas ajustadas em uma idade futura e como variáveis independentes, os parâmetros dessa mesma função em uma idade atual e as características do povoamento nas idades atual e futura. Os seguintes atributos do povoamento foram utilizados: idade, diâmetro máximo, diâmetro mínimo, diâmetro quadrático e número de árvores. Os ajustes das equações foram efetuados no programa Statistica 8.0 (STATSOFT, INC 2010).
2.4. Avaliação da redistribuição teórica
A redistribuição teórica dos diâmetros foi avaliada pela capacidade do sistema de equações em projetar o diâmetro máximo, o número total de árvores e os parâmetros da função hiperbólica (α, β e γ) para qualquer idade. O coeficiente de correlação entre os valores observados e estimados de frequência por classe de diâmetro e a análise gráfica dos resíduos foram utilizados para avaliar a exatidão das estimativas de frequência por classe de diâmetro. Foram avaliadas a consistência, a capacidade de projeção e a realidade biológica do sistema de equações, projetando diferentes distribuições iniciais observadas para idades futuras e comparando com os respectivos valores verificados de frequência por classe de diâmetro. Para avaliar a exatidão das estimativas geradas pelas equações que compõem o modelo de distribuição de diâmetros foram estimadas, ainda, as seguintes estatísticas: raiz quadrada do erro médio (RMSE), bias e bias% (SIIPILEHTO, 2000). Os estimadores utilizados foram:
em que e Yi e i = valores observados e estimados da variável sob análise; e n = número de casos.
3. RESULTADOS
As frequências estimadas de cada parcela e idade foram comparadas com as observadas. Todos os ajustes apresentaram aderência aos dados pelo teste Kolmogorov-Smirnorv (p>0,01).
A análise gráfica dos resíduos utilizada para escolha das equações para compor os modelos de distribuição diamétrica é apresentada na Figura 1. As equações que compõem os modelos ajustados são apresentadas a seguir:
Função hiperbólica com três parâmetros
Função hiperbólica com dois parâmetros
Função hiperbólica com α = dmin
Função hiperbólica truncada
em que I1 e I2 são as idades atual e futura, respectivamente, em meses; γ1 e γ 2 são os parâmetros de forma função hiperbólica nas idades atual e futura; β1 e β2 são os parâmetros de escala nas idades atual e futura; α1 e α2 são os parâmetros de escala nas idades atual e futura; dmax1 e dmax2 são os diâmetros máximos nas idades atual e futura, em cm; dmin1 e dmin2 são os diâmetros mínimos nas idades atual e futura, em cm; N1 e N2 são as frequências de árvores por hectare nas idades atual e futura; e Ln indica o logaritmo neperiano.
Na Figura 2, verifica-se a evolução da distribuição diamétrica realizada em uma parcela escolhida ao acaso em três das idades mensuradas com os quatro modelos propostos, com intensidade de desbaste de 35%. Os valores de RMSE, bias e bias% dos modelos ajustados são apresentados na Tabela 1.
4. DISCUSSÃO
Este estudo é o primeiro a comparar a performance do ajuste de diferentes formas da função hiperbólica para a projeção da distribuição diamétrica de povoamentos desbastados de eucalipto. Os dados utilizados descrevem bem as variações encontradas em plantios clonais de eucalipto nas condições edafoclimáticas descritas, sendo o número de parcelas e medições utilizadas representativos em estudos de modelagem (CAMPOS; LEITE, 2009).
A aderência aos dados pela função hiperbólica foi comprovada pelo teste KS, o que demonstra a flexibilidade da função na descrição da estrutura diamétrica de povoamentos desbastados. Diversos autores utilizam esses critérios para a escolha da fdp a ser utilizada (BARRA at al., 2004; LEITE et al., 2010).
A utilização do parâmetro de locação para o ajuste das funções acarreta diminuição significativa nos coeficientes de correlação entre valores observados e estimados das equações que compõem os modelos de crescimento e produção. Em função da aplicação dos desbastes, a recuperação do parâmetro de locação fica comprometida, o que pode ser observado pelo valor de bias% obtidos nas equações desse parâmetro, sugerindo a subestimação do valor deste (Tabela 1). Contudo, essa variação pode ser explicada pelo fato de o valor do parâmetro em determinadas situações convergir para a origem em distribuições mais achatadas, ou para o valor do diâmetro mínimo em distribuições mais estreitas.
Apesar da alta dispersão dos resíduos para as funções com a utilização do parâmetro de locação, a análise do histograma de resíduos demonstra a não existência de tendenciosidade na recuperação dos parâmetros em idades futuras.
5. Conclusão
Conclui-se que a função hiperbólica nas formas testadas pode ser utilizada para a descrição e projeção da estrutura diamétrica de povoamentos de eucalipto submetidos a desbastes, podendo a não supressão do parâmetro de locação apresentar redução na eficiência das estimativas dos modelos.
6. REFERÊNCIAS
Recebido em 03.09.2010 e aceito para publicação em 02.05.2011.
- BARRA, O. S. V. et al. Proposta metodológica para o ajuste ótimo da distribuição diamétrica Weibull 3P. Floresta, v.34, n.3, p.387-393, 2004.
- BURKHART, H. E.; CAO, Q. V.; WARE, Z. D. A. A comparison of growth and yield prediction models for loblolly pine Blacksburg: Virginia Polytechnic Institute and State University, Scholl of Forestry and Wildlife Resources, 1981. 59p
- CAMPOS, J. C. C.; LEITE, H. G. Mensuração florestal: Perguntas e respostas. 3.ed. Viçosa, MG: Universidade Federal de Viçosa, 2009. 548p.
- CLUTTER, J. L. et al. Timber management: a quantitative approach. New York: John Willey & Sons, 1983. 333p.
- DEMOLINARI, R. A. et al. Crescimento de plantios clonais de eucalipto não desbastados na região de Monte Dourado (PA). Revista Árvore, v.31, n.3, p.503-512, 2007.
- GIBBONS, J. D.; SUBHABRATA, C. Nonparametric statistical inference 3.ed. New York: Marcel Dekker, 1992. 544p. (Statistics: textbook and monograph, 31).
- GUIMARÃES, D. P. Desenvolvimento de um modelo de distribuição diamétrica de passo invariante para prognose e projeção da estrutura de povoamentos de eucalipto 178f. Tese (Doutorado em Ciência Florestal) - Universidade Federal de Viçosa, Viçosa, MG, 1994.
- GUIMARÃES, D. P. Uma função hiperbólica de distribuição probabilística de alta flexibilidade Planaltina: Embrapa Cerrados, 2002. 40p.
- LEITE, H. G. et al. Avaliação do ajuste das funções Weibull e hiperbólica a dados de povoamentos de eucalipto submetidos a desbaste. Revista Árvore, v.34, n.2, p.305-311, 2010.
- NOGUEIRA, G. S. et al. Modelo de distribuição diamétrica para povoamentos de Eucalyptus sp. submetidos a desbaste. Revista Árvore, v.29, n.4, p.579-589, 2005.
- SIIPILEHTO, J. A comparison of two parameter prediction methods for stand structure in Finland. Silva Fennica, v.34, n.4, p.331-349, 2000
Datas de Publicação
-
Publicação nesta coleção
10 Out 2011 -
Data do Fascículo
Out 2011
Histórico
-
Recebido
03 Set 2010 -
Aceito
02 Maio 2011