Classificação da Capacidade Produtiva de Povoamentos de Eucalipto Utilizando
Diâmetro Dominante
Classificação da Capacidade Produtiva de Povoamentos de Eucalipto Utilizando
Diâmetro Dominante
Helio Leite*, Renato Castro*, 1Antonilmar Silva*, Carlos Júnior*, Daniel
Binoti*, Ana Flávia Castro* e Mayra Binoti*
*Engenheiro Florestal
Universidade Federal de Viçosa. Departamento de Engenharia Florestal, Av. PH
Rolfs, s/n, CEP 36570-000 Viçosa, BRASIL
1Cenibra. BR 381 Km, 172 s/n, 35195000 Belo Oriente, MG, BRASIL
Sumário
Este estudo foi conduzido para avaliar as seguintes hipóteses: H0(1): a
classificação da capacidade produtiva de povoamentos equiénios de eucalipto não
desbastados pode ser feita empregando o diâmetro das árvores identificadas como
dominantes; H0(2): índices de local definidos com base no diâmetro dominante
são mais correlacionados com a área basal e a produção em volume, do que os
índices definidos com base na altura das árvores dominantes; e H0(3): a
variável altura dominante pode ser substituída com eficiência, em modelos
hipsométricos, pela variável diâmetro dominante. Para avaliação dessas
hipóteses foram empregados dados de 523 parcelas permanentes de inventários
florestais contínuos conduzidos em povoamentos de clones de eucalipto, com
espaçamento inicial de 3,0 x 3,0m. As curvas de índices de local foram
construídas empregando o método da curva-guia assumindo uma idade índice de 70
meses, ao estimar a altura dominante (A) ou diâmetro dominante (B) em função da
idade. Em seguida, o modelo de Clutter foi ajustado com o índice de local
obtido pelas duas alternativas (A e B). A partir dos modelos ajustados pelas
duas alternativas foram também simulados diferentes combinações de idades e
pesos de desbaste. Em ambos os casos, o modelo foi ajustado com todos os dados
disponíveis (geral) e por classes de capacidade de produção (produtividade
alta, média e baixa). As hipóteses 1 e 3 foram aceites e a hipótese 2 foi
rejeitada.
Palavras-chave: Índice de local; Eucalyptus sp.; método da curva-guia
Site Quality Classification of Eucalyptus Stands Using Dominant Diameter
Abstract
The present study was conducted to evaluate the following hypotheses: H0(1):
the productive capacity classification of no-thinned even-aged eucalyptus
stands can be performed using the diameter of the trees identified as dominant;
H0(2): site index based on the dominant diameter are more correlated to the
basal area and to the volume production, compared to the site index based on
the dominant trees height; and H0(3): the dominant height variable can be
efficiently replaced in hypsometric models by the dominant diameter variable.
The data used in this study came from 523 permanent sample plots measured in
continuous forest inventories carried out in eucalyptus stands, with initial
spacing of 3,0 x 3,0m. The site index curves were constructed employing the
guide curve method assuming the index-age of 70 months to estimate the dominant
height (A) or dominant diameter (B) according to age. Then, the Clutter model
was adjusted using the site index achieved by the two alternatives (A and B).
Different thinning ages and weights combinations were also simulated from the
two alternatives adjusted models. In both cases, the model was adjusted with
all the available data (general) and per productive capacity classes (high,
average and low productivity). Hypotheses 1 and 3 were accepted and hypothesis
2 was rejected.
Key words: Site index; Eucalyptus sp.; guide curve method
Classification de la Qualité du Site d'Eucalyptus Utilisant le Diamètre
Dominante
Résumé
Cette étude a été menée pour évaluer les hypothèses suivantes: H0 (1): la
classification de la capacité productive des peuplements équiens d'eucalyptus
peut être faite en utilisant le diamètre des arbres identifiés comme dominants;
H0 (2): les index du site définis sur la base du diamètre dominant ont plus de
corrélation avec la surface terrière et le volume de production, que les index
fixés en fonction de la hauteur des arbres dominants; et H0 (3): la variable
hauteur dominante peut être remplacée efficacement dans les modèles
hypsométriques, par la variable diamètre dominant. Pour évaluer ces hypothèses,
on a utilisé les données de 523 parcelles permanentes localisées dans des
peuplements de clones d'eucalyptus avec un espacement initial de 3,0 x 3,0m
sujets à des inventaires forestiers continus. Les courbes de l'indice de
station ont été construites en utilisant la méthode de la courbe guide, en
supposant un âge-indice de 70 mois, afin d'estimer la hauteur dominante (A) ou
le diamètre dominant (B) en fonction de cet âge. Ainsi, le modèle de Clutter a
été ajusté avec l'indice local obtenu par les deux variantes (A et B). À partir
des modèles établis par les deux possibilités, différentes combinaisons d'âges
et d'éclaircies ont également été simulées. Dans les deux cas, le modèle était
établi avec toutes les données disponibles (en général)et par classes de
capacité de production (productivité élevée, moyenne et basse). Les hypothèses
1 et 3 ont été acceptées et l'hypothèse 2 a été rejetée.
Mots clés: Index du site; Eucalyptus sp.; méthode de la courbe guide
Introdução
A capacidade produtiva de um local refere-se ao potencial desse lugar para
produzir madeira ou outro tipo de produto, sob a condição ambiental existente
(clima e solo) e técnicas silviculturais empregadas. A interação entre esses
fatores expressa a qualidade do local para uma espécie ou grupo de espécies
(CAMPOS e LEITE, 2009) sendo um indicador da produtividade potencial. Na
maioria das vezes a classificação é feita por meio de um índice, definido como
a altura média das árvores dominantes, ou na média das alturas das maiores
árvores, ou ainda, na média de uma quantidade pré-definida de árvores de maior
porte, em uma idade de referência. Isto tem sido feito com base na
pressuposição de que a altura não é afetada de modo significativo pela
densidade do povoamento e devido à existência de correlação significativa com a
produção volumétrica (HUSCH et al., 2003; CARMEAN, 1975). Esses índices são
utilizados em modelagem de crescimento e produção como uma das variáveis
independentes do modelo (OLIVEIRA et al., 2009).
POKHAREL e FROESE (2009) mencionam que uma das razões para o emprego de índices
de local é a possibilidade de medição da altura com exatidão e custo mínimo.
Porém, conforme STOUT e SHUMWAY (1982), em alguns casos, o emprego da variável
altura para classificação da capacidade produtiva pode resultar em erro de não
amostragem. Ainda, em certos casos, a identificação da posição para leitura da
altura total é dificultada pela alta densidade de plantio ou condições
ambientais adversas (pouca radiação solar direta e ventos), resultando em erros
que posteriormente serão inseridos na classificação.
O principal componente do custo de medição de uma parcela permanente em
povoamentos de eucalipto é a altura total. Por isso é comum, no Brasil, medir
apenas algumas alturas em cada parcela visando o ajuste de modelos
hipsométricos, em nível de parcela, talhão ou povoamento. Nesses modelos são
incluídas, além do diâmetro, variáveis independentes como a altura dominante da
parcela, o índice de local e a idade. A alternativa mais frequente consiste em
correlacionar a altura das árvores individuais com os seus diâmetros a 1,3m e
com a altura dominante da parcela.
O diâmetro a 1,3m de altura (dap) também pode ser utilizado para classificação
da capacidade produtiva (AVERY e BURKHART, 1994). Esta variável não tem sido
utilizada por ser mais influenciada por fatores ambientais do que a altura e
ser dependente da densidade do povoamento. Apesar disso, o diâmetro pode estar
mais correlacionado com a área basal e com o volume, do que a altura dominante
média. Uma justificativa para não determinar índices de local com base no
diâmetro médio é o fato desse diâmetro ser afetado diretamente pela densidade,
especialmente em plantios submetidos a desbaste ou com grande ocorrência de
mortalidade regular. Mas, se as árvores dominantes forem identificadas e
marcadas no povoamento, haveria eficiência na classificação ao empregar o
diâmetro dessas árvores (dapd) ao invés de suas alturas?
O diâmetro de árvores individuais é medido empregando suta ou fita métrica (ou
diamétrica), enquanto a altura é obtida empregando clinômetros ou hipsômetros.
Nos dois casos podem ocorrer erros de medição (SOARES et al., 2007) e o
controle desse erro é mais fácil no caso do diâmetro. Apesar disto, até o ano
de 2011 não foram encontrados trabalhos indicando o uso do diâmetro das árvores
dominantes na construção de curvas de índices de local.
Índices de local são utilizados como inputs em modelos de crescimento e
produção (CLUTTER, 1963). Nesses modelos é assumido que o crescimento depende
da idade, do grau de ocupação da área (geralmente área basal) e da capacidade
produtiva, expressa pelo índice de local. A partir destes modelos é possível
simular o crescimento do povoamento com ou sem intervenções de desbaste
(GÖRGENS et al., 2007).
Por meio dos desbastes, pode-se aumentar a área basal individual, devido ao
aumento do espaço vital e à maior disponibilidade de luz e nutrientes,
compensando-se a retirada de árvores menores e mantendo a área basal total do
povoamento quase inalterada, obtendo, no final, árvores de maiores diâmetros
(GLUFKE et al., 1997). O modelo de CLUTTER é largamente utilizado na área
florestal para prognose da produção, e suas qualidades são amplamente estudadas
e discutidas na literatura florestal (GÖRGENS et al., 2007). O seu uso para
simulação de desbaste é possível por se tratar de um modelo de densidade
variável (CAMPOS e LEITE, 2009), assumindo uma tendência de crescimento da
floresta após o desbaste similar ao crescimento antes do desbaste (GÖRGENS et
al., 2007).
Em povoamentos com a mesma espécie, sob mesmo tipo de solo, com as mesmas
condições ambientais e um mesmo arranjo espacial, podem ocorrer grandes
variações da variável área basal entre parcelas permanentes. Na maioria das
vezes são observadas correlações (lineares ou exponenciais) entre a área basal
e o índice de local determinado com base na altura dominante. Matematicamente a
área basal está mais relacionada com o dap do que com a altura, o mesmo
acontecendo com o volume. Assim, podem ser formuladas as hipóteses: H0(1): a
classificação da capacidade produtiva de povoamentos equiénios não desbastados
pode ser feita empregando o diâmetro das árvores identificadas como dominantes
no povoamento (dapd), ao invés do uso da variável altura dominante,
tradicionalmente utilizada, e H0(2): índices de local definidos com base no
dapd são mais correlacionados com a produção, em área basal e em volume, do que
os índices definidos com base na altura das árvores dominantes (Hd).
Em muitas empresas florestais brasileiras é comum o emprego de modelos
hipsométricos que utiliza a variável independente Hd para estimar a altura de
árvores individuais em inventários florestais. Estes modelos têm sido ajustados
por estrato, sendo os estratos definidos com base em idade, regime de corte,
espaçamento e capacidade produtiva do lugar. Ao decidir pela classificação da
capacidade produtiva empregando o dap da árvore dominante, ao invés de sua
altura dominante, a variável LnHd deixa de estar disponível e poderia ser
substituída por Lndapd. Portanto surge uma terceira hipótese, H0(3): a variável
Hd pode ser substituída com eficiência, em modelos hipsométricos, pela variável
dapd.
Este estudo foi conduzido empregando dados de inventários conduzidos em
povoamentos de clones de eucalipto no Brasil para avaliar as hipóteses H0(1),
H0(2) e H0(3), sempre comparando os resultados obtidos pelo método tradicional,
que utiliza a variável Hd em relação aos resultados utilizando a variável dapd.
Material e métodos
Os dados utilizados neste estudo foram provenientes de inventários florestais
contínuos realizados entre os anos de 2001 e 2006, em povoamentos de clones de
Eucalyptus urophylla x Eucalyptus grandis, com espaçamento inicial de 3,0 x
3,0m, na região de Santa Bárbara, Estado de Minas Gerais, Brasil. A vegetação
original da área onde se encontram os povoamentos (4.162ha) era composta de
campos de matas de candeia (Eremanthus sp.), campos de altitude e áreas de
transição de cerrado. Os solos predominantes são do tipo Latossolo, com 46 a
63% de argila, densidade entre 1,06 e 1,13g/cm3, teor de matéria orgânica entre
1,73 e 2,16g/kg e soma de bases variando entre 0,15 e 0,35 molc/kg. As
temperaturas mínima, média e máxima são de 17, 22 e 27°C, com umidade relativa
média de 65% e precipitação média anual de 1440 mm. As médias de radiações
global e fotossinteticamente ativa são de 15,9 mj/m2/dia e 33,6 mmol/m2/s. Nos
4.162ha a produtividade média, aos 7 anos, varia entre 26 e 44m3.ha-1.ano-1, de
madeira para produção de celulose com diâmetro mínimo de 6 cm sem casca.
Foram utilizados dados de 523 parcelas permanentes, mensuradas em idades entre
30 e 118 meses. Nesse intervalo de tempo as árvores dominantes atingiram até 40
m de altura e 30 cm de diâmetro e a amplitude de volumes observada foi de 52,34
a 475,65m3.ha-1. A sobrevivência no período de avaliação variou de 90 a 100%.
A altura dominante foi definida como a média das alturas totais de cinco
árvores dominantes de cada parcela (Hd). Essas árvores foram identificadas por
ocasião da primeira medição da parcela, sendo selecionadas árvores sadias, com
fuste bem definido, sem bifurcação e tortuosidade, livre de danos e com a copa
situada no dossel superior, recebendo radiação solar direta. O diâmetro
dominante (dapd) foi definido como a média aritmética dos diâmetros dessas
cinco árvores dominantes. Foram mensurados também, em cada parcela e medição, o
diâmetro médio quadrático (q), a média das alturas totais e a média dos
diâmetros (dapm).
O comportamento das variáveis dap, q, dapd, Ht e Hd ao longo dos anos, foi
avaliado empregando o modelo apresentado por PIENAAR e SCHIVER (1981):
em que Y1 e Y2 são as estimativas da variável nas idades atual (I1) e futura
(I2), β0 e β1 são os parâmetros do modelo, com ε ~ N (0, σ²).
A qualidade dos ajustamentos foi avaliada com base nas estatísticas: média das
diferenças absolutas (MDA), coeficiente de correlação entre valores observados
e valores estimados da variável dependente , raiz quadrada do erro médio
(RQEM), bias, bias% e desvio padrão do erro de predição (Sb) (SIIPILEHTO,
2000), além de análises gráficas de resíduos, sendo:
em que: = valor observado e estimado da variável sob análise; = erro relativo
de predição; n = número de casos.
Curvas de índices de local foram construídas empregando o método da curva-guia
(CLUTTER et al., 1983), assumindo uma idade-índice de 70 meses e o modelo de
Schumacher (eq. 8), em que Y é a variável dependente, I é a idade em meses e as
demais como definidas anteriormente.
Este modelo tem sido empregado com eficiência na região onde foram obtidos os
dados. Foram construídas curvas de índices de local para Y = Hd e Y = dapd. A
qualidade dos ajustes foi avaliada pelo coeficiente de determinação ajustado
(R²aj), MDA, RQEM, bias, bias%, Sb e pela análise gráfica da relação entre
idade e valores de Y, observados e estimados.
Para avaliar a hipótese H0(1), cada parcela foi classificada em relação à
capacidade produtiva (alta, média ou baixa produtividade) pelas alternativas A
(índices de local definidos com base na variável explicativa Hd) e B (índices
de local com base na variável dapd). Em seguida foram estimadas as percentagens
de casos com coincidência entre essas duas classificações.
Para avaliar as tendências de crescimento dos povoamentos, em área basal e
volume, foi utilizado o modelo de equações simultâneas de CLUTTER (1963) (eq.
10). Esse modelo tem sido empregado com sucesso na região de estudo para
diferentes clones na mesma forma definida por Clutter, para Pinus taeda, em
1963, sendo:
em que: = índice de local na idade I1, pela alternativa j, em m para a
alternativa A e em cm para alternativa B; B1 = área basal inicial, em m2.ha-1;
B2 = área basal futura, em m2.ha-1; V2 = volume futuro com casca, em m3.ha-1;
β0, β1, β2, β3 = parâmetros do modelo que projeta o volume; α1 e α2 =
parâmetros do modelo que projeta a área basal; e ε = erro aleatório.
O modelo de Clutter foi ajustado pelo método dos mínimos quadrados em dois
estágios, utilizando os índices de local obtidos pelas alternativas A e B.
Dessa forma, pode-se avaliar a segunda hipótese (H0(2)), ao verificar se
haveria ganho nas estimativas de produção em área basal e em volume, ao
utilizar índices de local definidos com base no dapd, ao invés de Hd.
Em ambos os casos o modelo foi ajustado com todos os dados disponíveis (geral)
e por classes de capacidade produtiva (alta, média e baixa produtividade).
Esses ajustes foram feitos para verificar o ganho de exatidão ao construir
curvas de produção na forma estratificada. A qualidade dos ajustamentos foi
avaliada com base no coeficiente de determinação, no sinal e significância dos
coeficientes das regressões e nas estatísticas MDA, RQEM, bias, bias% e Sb além
de análises gráficas de resíduos para as variáveis área basal e volume.
Foram elaborados quadros de produção para as duas alternativas, conforme DIAS
et al. (2005); definidas as idades técnicas de corte (ITC) e avaliado o
realismo biológico do sistema de equações ajustado, conforme CAMPOS e LEITE
(2009). Foram feitas simulações de desbastes aos 60 meses, com intensidade de
40% da área basal e aos 100 meses, com intensidade de 35%, com o objetivo de
verificar a eficiência da alternativa B, em relação à tradicional alternativa
A.
O modelo hipsométrico (eq. 11) foi ajustado considerando Z = Hde Z = dapd.
Este modelo é empregado com frequência em inventários florestais conduzidos na
região de estudo. O objetivo ao ajustar este modelo foi de avaliar a eficiência
da variável dapd em substituição à variável Hd, referente à terceira hipótese
(H0(3)).
As estimativas de Ht foram avaliadas pelas estatísticas MDA, RQEM, bias, bias%
e Sb.
As análises de regressão lineares e não-lineares foram realizadas com o auxílio
do programa STATISTICA 8.0 e o ajuste do modelo de CLUTTER (eq. 10) com o
auxílio do programa EVIEWS 5.0.
Resultados e discussão
Crescimento em altura e diâmetro
As equações e as medidas de precisão dos ajustes para as variáveis Ht, Hd, dap,
q e dapd, estão relacionadas no Quadro 1. Todas as estimativas dos parâmetros
foram significativas pelo teste t (p<0,01), precisas e não tendenciosas (Quadro
1 e Figura 1). Para todas as variáveis, os valores estimados foram iguais aos
observados pelo teste F, conforme GRAYBILL (1976).
Quadro_1
' Estimativas dos coeficientes dos modelos de regressão (eq. 1) ajustados para
crescimento em alturas e diâmetros de árvores de um povoamento de clones de
Eucalyptus, e as respectivas medidas de precisão
Figura 1 ' Relação entre valores estimados e valores observados de altura média
(a), altura dominante (b), média de diâmetros (c), diâmetro médio (d) e
diâmetro dominante (e).
As tendências estimadas de crescimento das variáveis de alturas e diâmetros
(Quadro_1) são indicadas na Figura 2. Com base nos resultados, verificou-se
que, assim como a variável Hd, a variável dapd não atingiu precocemente um
valor assintótico, que indicaria a estagnação do crescimento. Por esse motivo
foram empregados modelos exponenciais para relacionar essas variáveis com a
idade. Analisando as tendências de crescimento em altura (Figura 2-a),
verificou-se que a altura dominante assumiu valores superiores às médias das
alturas para todas as idades, e ao longo do tempo esta diferença tendeu a ser
maior. CAMPOS e LEITE (2009) destacam a vantagem da variável Hd em relação à
média das alturas (Ht), por ela não ser influenciada pela competição, uma vez
que a média das alturas estagna mais cedo. Este mesmo comportamento pôde ser
verificado ao avaliar o crescimento em diâmetro dominante (Figura 2-b), que
também se mostra menos sensível à competição ao longo do tempo, ao contrário do
diâmetro médio e da média dos diâmetros.
Figura 2 ' Curvas de crescimento em altura total (Ht) e altura dominante (Hd) -
(a) e média dos diâmetros (dap), diâmetro médio (q) e diâmetro dominante (dapd)
- (b) em clones de Eucalyptus sp., sob arranjo espacial de 3,0 x 3,0 m
Classificação da capacidade produtiva
O modelo de Schumacher (eq. 8), escolhido para a classificação dos índices de
local (S), ajustou-se bem aos dados nas duas alternativas. Todas as estimativas
dos parâmetros foram significativas pelo teste t (p< 0,01). As equações e
estatísticas de ajuste obtidas com as alternativas A e B estão apresentadas no
Quadro 2.
Quadro 2 ' Estatísticas dos modelos de regressão (eq. 8) para classificação da
capacidade produtiva, utilizando as variáveis independentes Hd (A)e dapd (B)
A partir da expressão e do conceito de índice de local foi gerado o feixe de
curvas anamórficas para as alternativas A e B. Três classes de capacidade
produtiva (alta, média e baixa) foram definidas para as duas alternativas
(Figura 3). O número de parcelas igualmente classificadas pelas alternativas A
e B foi apresentado no Quadro 3.
Figura_3 ' Curvas de índice de local para clones de Eucalyptus sp., construídas
pelo método da curva-guia para as alternativas (A), a partir da variável
dependente altura dominante (Hd) e (B) pela variável diâmetro dominante (dapd)
Quadro 3 ' Número de parcelas classificadas como alta, média e baixa capacidade
produtiva, pelas alternativas (A) e (B)
O uso da variável dapd para classificação da capacidade produtiva mostrou-se
estatisticamente coerente. A exatidão das estimativas de Hd obtidas neste
estudo foi compatível com as estimativas obtidas para outros povoamentos de
clones de eucalipto, conforme aquelas obtidas por Soares et al. (2004) e
Demolinari et al. (2007). Analisando as medidas de precisão das equações,
verificou-se que os modelos das duas alternativas (A: empregando Hd e B:
empregando dapd) ajustaram-se bem aos dados. O crescimento da variável dapd em
função da idade também se comportou de forma exponencial (Figura_3), semelhante
ao crescimento em Hd, indicando a possibilidade de sua utilização para fins de
classificação da capacidade produtiva.
Ao comparar a classificação feita pelas alternativas verificou-se que em 71%
das parcelas permanentes houve coincidência em termos de classe de local (alta,
média e baixa capacidade produtiva). A hipótese H0(1) é verdadeira, visto que a
classificação da capacidade produtiva de povoamentos de clones de eucalipto,
empregando o diâmetro das árvores dominantes, foi tão eficiente quanto a
classificação utilizando a altura dominante.
Crescimento e produção florestal
As equações que projetam a área basal e o volume com casca por hectare, nas
duas alternativas, são apresentadas no Quadro 4 e as respetivas estatísticas no
Quadro 5. Conforme Campos e Leite (2009), duas características do modelo de
Clutter devem ser imediatamente verificadas: o parâmetro β1 deverá ser negativo
e o coeficiente α2, positivo. Apenas o parâmetro α2 da classe de baixa
capacidade produtiva da alternativa A não foi estatisticamente significativo, o
que indica que não houve efeito da capacidade produtiva sobre o crescimento em
área basal. Neste caso particular, o termo (1-I1.I2-1)S1 foi eliminado do
modelo. Na Figura 4 estão relacionadas às áreas basais e aos volumes estimados
em função dos observados, empregando as duas alternativas.
Quadro 4 ' Equações ajustadas pelo modelo de Clutter (eq. 10) para povoamentos
de clones de Eucalyptus sp., para todos os dados (geral) e com estratificação
por classe de produtividade (alta, média e baixa) definida pelas alternativas A
e B, e suas medidas de precisão*
Quadro 5 ' Estatística das equações ajustadas pelo modelo de Clutter (eq. 10)
para povoamentos de clones de Eucalyptus sp., para todos os dados (geral) e com
estratificação por classe de produtividade (alta, média e baixa) definidas
pelas alternativas A e B
Figura 4 ' Relação entre valores observados e estimados de área basal e volume
empregando índices de local determinados com base na altura dominante (a e b) e
no diâmetro das árvores dominantes (c e d), para classes de capacidade
produtiva Alta, Média e Baixa
A partir dos quadros de produção obtidas pelas duas alternativas, determinou-se
a produtividade média (m³.ha-1.ano-1), na idade técnica de corte (ITC), para
cada índice de local, além da produção volumétrica correspondente (Quadro 6).
Quadro 6 ' Produção volumétrica total nas idades técnicas de corte de
povoamentos de clones de Eucalyptus sp., utilizando o modelo de Clutter e
índices de local determinados com base na altura dominante (A) e no diâmetro
das árvores dominantes (B)
As idades técnicas de corte definidas pela máxima produtividade média foram
menores nos locais de maior capacidade produtiva em relação aos locais de baixa
capacidade. Conforme esperado, isto ocorreu devido à antecipação da idade de
estagnação do crescimento da floresta estar ocorrendo em idade mais jovem nos
locais mais produtivos. Esta consistência foi observada ao empregar as duas
alternativas de classificação (A e B), principalmente ao ajustar o modelo de
Clutter por estrato.
Ao empregar os índices de local das duas alternativas no modelo de projeção de
área basal e volume, pode-se observar a mesma eficiência na estimação destas
variáveis para as alternativas A e B respetivamente. Este resultado foi o mesmo
ao comparar as estimativas obtidas pelo modelo geral e pelos modelos ajustados
por estrato. A hipótese H0(2) não é verdadeira, visto que não apresentaram
resultados superiores aos encontrados com a variável Hd, mas similares.
Relação hipsométrica
O modelo hipsométrico (eq. 11) foi ajustado considerando Hd e dapd. As
estatísticas estão relacionadas no Quadro 7.
Quadro 7 ' Estatísticas dos modelos de regressão ajustados para equação
hipsométrica (eq. 11), utilizando as variáveis independentes Hd e dapd
A variável Hd tem sido utilizada em modelos hipsométricos, principalmente em
inventários florestais. O uso desta variável permite diferenciar a altura de
árvores de mesmo dap em diferentes locais, ou seja, o efeito da capacidade
produtiva sobre a altura de árvores individuais é explicado pela altura média
das árvores dominantes de cada lugar. Ao decidir pelo emprego da variável dapd
em lugar da variável Hd assume-se que as alturas das árvores dominantes não
seriam medidas nas parcelas de inventário. Considerando a necessidade de
diferenciar a altura de árvores individuais em diferentes locais, naturalmente
seria empregada a variável dapd. Ao empregar esta variável no modelo
hipsométrico, não houve perda de precisão, tornando a hipótese H0(3)
verdadeira.
Simulação de desbastes
Embora os povoamentos não tenham sido submetidos a desbaste, a eficiência da
classificação feita com base no diâmetro dominante (dapd), em relação à
classificação usual (com base em Hd), foi avaliada simulando desbastes com
remoção de 40% e 35% da área basal remanescente nas idades de 60 e 100 meses,
respetivamente, para os modelos geral e estratificados. Os modelos de
crescimento e produção ajustados, empregando índices de local definidos pelas
alternativas A e B, foram eficientes em descrever o comportamento do
crescimento do povoamento antes e após a aplicação de desbaste, não sendo
observadas diferenças de tendência de crescimento entre as alternativas.
Pequenas diferenças foram observadas em locais de maior capacidade produtiva,
ao empregar o modelo geral. Este fato não inviabiliza o emprego da alternativa
B, uma vez que os modelos de crescimento e produção em geral são ajustados por
estrato.
A Figura 5 ilustra a simulação do crescimento das duas alternativas.
Figura 5 ' Simulação de desbastes em povoamentos de clones de eucalipto pelo
modelo de Clutter em três locais, empregando um modelo geral (a,b,c) e modelos
estratificados ajustados por classe de local (d,e,f), utilizando índices de
local definidos pela altura dominante (A) e pelo diâmetro das árvores
dominantes (B)
Estes resultados indicam que o crescimento em altura dominante e diâmetro
dominante foram semelhantes e pouco influenciados pela competição. O emprego da
variável dapdem lugar de Hd pode resultar em redução de custos de inventário e
modelagem.
Apesar dos resultados obtidos neste estudo é importante lembrar que nos casos
onde a mortalidade regular for significativa poderão ocorrer inconsistências na
classificação por meio do dapd.
Conclusões
Com base nos resultados deste estudo, pode-se afirmar que a hipótese H0(1) é
verdadeira. A classificação da capacidade produtiva de povoamentos de clones de
eucalipto, empregando o diâmetro das árvores dominantes, é eficiente. Esta
inferência se baseia no fato da classificação das unidades amostrais ter sido
coerente com aquela obtida pelo método usual (empregando Hd).
Índices de local definidos com base no diâmetro dominante não são mais
correlacionados com a área basal e a produção em volume do que os índices
definidos com base na altura das árvores dominantes, visto que apresentaram
resultados similares aos encontrados com a variável Hd. A hipótese H0(2) foi
rejeitada.
O uso da variável dapd em lugar de Hd pode ser empregado em modelos
hipsométricos. Não há perda de precisão ao empregar a variável alternativa. A
hipótese H0(3) foi aceita.
