Equações para estimar o estoque de carbono no fuste de árvores individuais e em plantios comerciais de eucalipto

Equations to estimate carbon stock in stems of individual trees ans stands of eucalyptus

Equações para estimar o estoque de carbono no fuste de árvores individuais e em plantios comerciais de eucalipto

Equations to estimate carbon stock in stems of individual trees ans stands of eucalyptus

Carlos Pedro Boechat Soares^I; Helio Garcia Leite^I; Eric Bastos Görgens^II

^IDepartamento de Engenharia Florestal da Universidade Federal de Viçosa (UFV). E.mail:<[e-mail]>; <[e-mail]>

^IIPrograma de Pós-Graduação em Ciência Florestal da UFV

RESUMO

O objetivo deste trabalho foi ajustar equações para estimar o carbono presente no fuste de árvores individuais de eucalipto e o estoque de carbono, por unidade de área, em diferentes condições de plantio. Para isso, foram utilizadas 532 árvores para ajustar a equação referente ao modelo alométrico de Schumacher e Hall (1933) e 95 parcelas permanentes para ajustar as equações referentes aos modelos de crescimento e produção, incluindo as seguintes variáveis independentes: idade (I), área basal (B) e índice de local (S). Após as análises, verificou-se que as equações se ajustaram bem aos dados observados, fornecendo estimativas precisas.

Palavra-chave: Eucalipto, carbono, modelos de crescimento e produção e equação alométrica.

ABSTRACT

The objective of this work was to adjust equations to estimate carbon in the stems of individual eucalyptus trees and carbon stock, by unit of area, in different planting conditions. For that, 532 trees were used to adjust the allometric model of Schumacher and Hall (1933) and 95 permanent plots to adjust equations related to growth and yield models, including the following independent variables: age (I), basal area (B) and site index (S). The analyses showed that the equations adjusted well to the observed data, providing accurate estimates.

Keywords: Eucalyptus, carbon stock, growth and yield models and allometric equation.

1. INTRODUÇÃO

A utilização de modelos para representar ou explicar uma realidade é uma estratégia muito utilizada em vários campos da ciência. Segundo Hasenauer (2000), modelos são apresentações simplificadas de uma realidade, em que o grau de detalhamento dos dados de entrada limita o detalhamento dos resultados.

Especificamente no setor florestal, os modelos de crescimento e produção volumétrica são utilizados no planejamento das atividades florestais. No entanto, o desenvolvimento e ajuste de modelos que quantifiquem o teor de carbono presente em uma floresta vêm se tornando uma necessidade, devido ao grande interesse com relação à captura de carbono por esses ecossistemas.

O interesse pela quantificação do estoque de carbono em florestas se deve à crescente demanda de alternativas que visem ao controle das emissões atmosféricas de CO₂, entre elas: a diminuição da emissão desse gás para a atmosfera e a remoção do CO₂ atmosférico por ecossistemas terrestres e pelos oceanos (SCHROEDER, 1992).

Nesse contexto, as florestas têm papel importante, pois são capazes de armazenar grandes quantidades de carbono em seus diferentes compartimentos. Segundo Grainger, citado por Schroeder (1992), somente os trópicos possuem cerca de 2 (dois) bilhões de hectares de áreas degradadas ou não florestadas, das quais 758 milhões poderiam ser reflorestadas.

A preocupação em relação à elevação da concentração de CO2 na atmosfera é legítima, uma vez que a conseqüência desse fato é o aquecimento atmosférico global, denominado efeito-estufa, o qual pode causar, segundo os especialistas, mudanças permanentes no clima, alterando o regime dos ventos, a pluviosidade e a circulação dos oceanos (SOARES e OLIVEIRA, 2002).

Embora as florestas possam desempenhar importante papel na captura do CO2 atmosférico, a atividade florestal também é responsável pelo aumento da concentração desse gás. De acordo com o uso da floresta, todo o carbono acumulado é devolvido à atmosfera, por exemplo na combustão da madeira para sua transformação em carvão ou no processo de polpação para a produção de papel.

Nesse sentido, existem muitas incertezas sobre a real contribuição para o fluxo de CO2 dos processos de desflorestamento e reflorestamento, uma vez que estimativas sobre a quantidade de carbono presente nas florestas ainda são feitas, utilizando-se métodos pouco precisos (COOPER, 1983; BROWN e LUGO, 1984; BROWN et al., 1986). Assim, este trabalho foi desenvolvido com o objetivo de ajustar equações para estimar o carbono presente no fuste de árvores individuais de eucalipto e o estoque de carbono, por unidade de área, em diferentes condições de plantio.

2. MATERIAL E MÉTODOS

Os dados para a realização deste estudo são provenientes de plantios de Eucalyptus grandis W. Hill ex Maiden, com idades entre 32 e 80 meses, pertencentes à Cenibra S/A (Celulose Nipo-Brasileira S/A), localizados na região leste do Estado de Minas Gerais.

Para o ajuste das equações de carbono para árvores individuais foram selecionadas aleatoriamente 532 árvores, distribuídas em diferentes classes de tamanho (Quadro 1). Em cada árvore-amostra, mediram-se o DAP (diâmetro com casca a 1,30 m do solo) e a altura total da árvore, bem como se realizou o procedimento de cubagem rigorosa para obtenção do volume do fuste com e sem casca, pela aplicação sucessiva da expressão de Smalian, em seções de 1 m de comprimento (HUSCH et al., 1993). Além disso, foram retirados discos de madeira e casca para determinação das densidades básicas da madeira e da casca (VITAL, 1984), que, multiplicadas pelos respectivos volumes, forneceram as estimativas de biomassa da madeira e casca do fuste.

Para obter as estimativas do carbono estocado no fuste de cada árvore-amostra (madeira e casca), multiplicou-se a estimativa da biomassa por 0,5, ou seja, considerou-se que 50% da biomassa é composta por carbono (HIGUCHI e CARVALHO JR., 1994).

As estimativas do carbono estocado nos fustes das árvores-amostra e as respectivas medidas de DAP e altura total (Ht) foram utilizadas no ajuste de equações referentes ao modelo linearizado de Schumacher e Hall (1933), cuja relação funcional é assim definida:

em que:

Ln = logaritmo neperiano;

Y = carbono total no fuste (madeira e casca) - CF, ou carbono na madeira - CM, em kg;

b_i = parâmetros do modelo, para i = 1, 2 e 3; e

e = erro aleatório.

Para avaliar a precisão das equações foram utilizados o coeficiente de determinação (R²), o erro-padrão da estimativa (S_y.x) e a análise gráfica de resíduos porcentuais (R%), obtidos por:

Utilizando dados de 95 parcelas permanentes, medidas entre 1996 e 2000, e as equações de carbono para árvores individuais previamente ajustadas, procedeu-se ao ajuste das equações referentes aos seguintes modelos de crescimento e produção (CAMPOS e LEITE, 2002):

em que:

Ln = logaritmo neperiano;

I = idade, em meses;

I₁ e I₂ = idade presente e idade futura, respectivamente, em meses;

B₁ e B₂ = área basal presente e área basal futura, respectivamente, em m²/ha;

S = índice de local, em metros; e idade índice = 60 meses;

Y = estoque de carbono no fuste (CF) ou na madeira (CM), em kg/ha;

Y₂ = estoque futuro de carbono no fuste (CF) ou na madeira (CM), em kg/ha;

b_i = parâmetros dos modelos, para i = 1, 2, ... 5; e

e = erro aleatório.

A avaliação da precisão das equações referentes aos modelos descritos anteriormente foi realizada, utilizando-se também o coeficiente de determinação (R²), o erro-padrão da estimativa (S_y.x) e a análise gráfica de resíduos porcentuais, além da determinação das idades técnicas de colheita (ITC) nas três classes de produtividade.

3. RESULTADOS E DISCUSSÃO

Analisando os coeficientes de determinação (R²) das equações ajustadas para estimar o carbono total no fuste (CF) e o carbono na madeira (CM) de árvores individuais de eucalipto, verifica-se que as equações se ajustaram bem aos dados observados. Além disso, todos os coeficientes foram significativos em nível de 95% de probabilidade (*), pelo teste "t", de Student.

Carbono Total no Fuste (CF)

Carbono na Madeira (CM)

Observando a distribuição dos resíduos das equações na Figura 1, verifica-se a exatidão das estimativas, pela amplitude dos valores dos resíduos. Quanto à homogeneidade da distribuição dos resíduos porcentuais, notou-se tendência de as equações superestimarem o carbono das árvores com DAP maiores que 30 cm.

Considerando somente os modelos de crescimento e produção, que incluem as variáveis independentes idade (I) e índice de local (S), obtiveram-se as seguintes equações para o estoque de carbono total no fuste (CF) e para o estoque de carbono na madeira (CM):

As variáveis independentes das equações ajustadas foram todas significativas a 95% de probabilidade (*), tanto para a estimação do estoque de carbono total no fuste quanto para o estoque de carbono na madeira, pelo teste "t", de Student. A inclusão da variável índice de local (S) aumentou significativamente a precisão das equações, quando comparadas com aquelas que possuíam apenas a idade (I) como variável independente, haja vista o aumento do coeficiente de determinação e a diminuição do erro-padrão da estimativa.

Os gráficos de resíduos, na Figura 2, indicam não haver tendência na estimação do estoque de carbono por hectare. No entanto, percebe-se maior exatidão das estimativas para as equações que apresentavam as variáveis independentes I e S, pela menor amplitude dos resíduos porcentuais.

Para o modelo de crescimento e produção de Clutter (1963) - modelos 5a e 5b, obtiveram-se as seguintes equações:

Pela análise dos coeficientes de determinação das equações que projetam a área basal por hectare (B₂), o estoque futuro de carbono no fuste por hectare (CF₂) e o estoque futuro de carbono na madeira por hectare (CM₂), verificou-se que as equações também se ajustaram bem a dados observados. Todas as variáveis, com exceção da variável "(1-I₁/I₂).S" na equação de área basal, foram significativas a 95% de probabilidade (*), pelo teste "t" de Student.

Os gráficos de resíduos (Figura 3), por sua vez, indicam não haver tendência na estimação dos estoques de carbono dos plantios de eucalipto, verificando-se, nessa figura, exatidão entre os valores observados e estimados, pela amplitude dos resíduos porcentuais.

De posse das equações referentes ao modelo de Clutter (1963), elaboraram-se as tabelas de produção de densidade variável (Quadros 2 e 3), considerando três classes de produtividade: S = 21,5; S = 26,5; e S = 31,5.

Com as estimativas dos Quadros 2 e 3, elaboraram-se os gráficos apresentados na Figura 4, os quais mostram as idades técnicas de colheita (ITC's) nas três classes de produtividade. De acordo com a Figura 4, verifica-se menor ITC nos melhores locais (maior índice de local - S), tanto no estoque total de carbono no fuste quanto no estoque de carbono na madeira. Cabe ressaltar que as idades técnicas de colheita foram ligeiramente menores quando se considera o estoque total de carbono no fuste (madeira + casca).

4. CONCLUSÕES

Após as análises, obtiveram-se as seguintes conclusões:

• As equações para estimar a quantidade de carbono presente no fuste de árvores individuais de eucalipto, em função do DAP (diâmetro com casca medido a 1,30 m do solo) e da altura total da árvore (Ht), apresentaram-se precisas, com coeficientes de determinação (R²) maiores que 97,9%;

• Os modelos de crescimento e produção forneceram estimativas precisas do estoque de carbono, e o modelo de Clutter (1963) apresentou coerência na determinação das idades técnicas de colheita (ITC), em razão da produtividade dos plantios.

5. REFERÊNCIAS BIBLIOGRÁFICAS

Recebido em 28.04.2004 e aceito para publicação em 10.08.2005.

• BROWN, S.; LUGO, A. E. Biomass of tropical forest: A new estimate based on forest volumes. Science, v.223, p.1290-1293, 1984.
• BROWN, S.; LUGO, A. E.; CHAPMAN, J. Biomass of tropical tree plantations and its implication for the global carbon budget. Canadian Journal of Forest Research, v.13, p.390-394, 1986.
• CAMPOS, J. C. C.; LEITE, H. G. Mensuração Florestal: perguntas e respostas. Viçosa, MG: Universidade Federal de Viçosa, 2002. 407p.
• CLUTTER, J.L. Compatible growth and yield models for lobolly pine. Forest Science, v.9, n.3, p.354-371, 1963.
• COOPER, C. F. Carbon storage in managed forest. Canadian Journal of Forest Research, v.13, n.1, p.155-165, 1983.
• HASENAUER, H. Princípios para a modelagem de ecossistemas florestais. Revista Ciência & Ambiente, n. 20, p. 53-69, 2000.
• HIGUCHI, N.; CARVALHO JR., J. A. Fitomassa, e conteúdo de carbono de espécies arbóreas da Amazônia. In: SEMINÁRIO EMISSÃO X SEQÜESTRO DE CO2 - uma nova oportunidade de negócios para o Brasil, 1994, Rio de Janeiro. Anais... Rio de Janeiro: Companhia Vale do Rio Doce, 1994. 125-153p.
• HUSCH, B.; MILLER, C. I.; BEERS, T. W. Forest mensuration 3.ed. Florida: Krieger Publishing Company, 1993. 402 p.
• SCHROEDER, P. Carbon storage potential of short rotation tropical tree plantations. Forest Ecology and Management, v.50, p.31-41, 1992.
• SCHUMACHER, F.X.; HALL, F. S. Logarithmic expression of timber-tree volume. Journal of Agricultural Research, v.47, n.9, p.719-734, 1933.
• SOARES, C. P. B.; OLIVEIRA, M. L. R. Equações para estimar a quantidade de carbono na parte aérea de árvores de eucalipto em Viçosa, Minas Gerais. Revista Árvore, v.26, n.5, p.533-539, 2002.
• VITAL, B. R. Métodos de determinação da densidade da madeira. Boletim Técnico SIF, n. 1, p. 1-21, 1984.

Datas de Publicação

• Publicação nesta coleção
  30 Nov 2005
• Data do Fascículo
  Out 2005

Histórico

• Aceito
  10 Ago 2005
• Recebido
  28 Abr 2004