Resumos

Função de afilamento para Virola surinamensis (Roll.) Warb

Tapering function for Virola surinamensis (Roll.) Warb

Helio Garcia Leite^I; João Ricardo Vasconcellos Gama^II; Jovane Pereira da Cruz^II; Agostinho Lopes de Souza^I

^IProfessor do Departamento de Engenharia Florestal, UFV, 36570-000 Viçosa-MG

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

RESUMO

Este estudo foi conduzido com o propósito de gerar equação de afilamento para Virola surinamensis (Roll.) Warb, que está entre as espécies mais exportadas pela indústria madeireira do estuário amazônico. Foram testados os modelos de Kozak (1969), Baldwin (1991), Demaerschalk (1972) e Garay (1979). Todos esses modelos estimaram o diâmetro comercial com precisão, porém os três primeiros resultaram em inconsistência ao estimar a altura comercial. O modelo de Garay foi selecionado e utilizado para gerar uma curva de cubagem para a espécie.

Palavras-chave: Várzea, manejo florestal, função de afilamento e Virola surinamensis.

ABSTRACT

This study was carried out with the aim of generating a taper equation for Virola surinamensis (Roll.) Warb, which is among the main exported species by the wood industries from the Amazonian estuary. The Kozak (1969), Baldwin (1991), Demaerschalk (1972) and Garay (1979) mathematical models were tested. All the tested models estimated the diameter with good precision, however, the first three models resulted in inconsistency while estimating commercial height. The Garay model was chosen and used to generate a cubage curve for the species.

Keywords: Floodplain, forest management, taper and Virola surinamensis.

1. INTRODUÇÃO

Na área do estuário amazônico, ao norte do Brasil, a pressão sobre os recursos da floresta tem causado redução no estoque de madeira de algumas espécies de interesse econômico. Como exemplo, tem-se a exploração da Virola surinamensis, que está entre as espécies madeireiras mais exportadas pela indústria madeireira do estuário amazônico.

Cada vez mais o mercado exige produtos diversificados e de qualidade. Como os produtos florestais não fogem à regra, faz-se necessário utilizar técnicas de quantificação e qualificação desses produtos, entre as quais se têm os modelos de afilamento que possibilitam a determinação do volume total e comercial de madeira para qualquer diâmetro ou altura especificados.

Os modelos de afilamento são relações funcionais que permitem descrever o perfil longitudinal de um tronco. São funções consideradas dinâmicas, porque possibilitam a estimativa do diâmetro (di) de dada altura (hi) e altura a qualquer diâmetro especificado, como também servem para fazer o sortimento e estudar a evolução da forma da árvore ao longo de sua vida (KOZAK et al.,1969).

A integral de um modelo de afilamento permite estimar o volume entre quaisquer segmentos ao longo do fuste – procedimento denominado sortimento – e, desse modo, maximiza-se o aproveitamento de uma árvore e, conseqüentemente, de um povoamento florestal. O estudo do sortimento de povoamentos de florestas de produção (nativas ou plantadas) é uma das ferramentas importantes no planejamento da utilização dos recursos florestais e na avaliação do manejo empregado, em razão de permitir a quantificação da matéria-prima, classificando-a de acordo com as dimensões e qualidade, para os múltiplos usos da madeira (FISCHER, 1997).

Segundo Guimarães e Leite (1992), é uma técnica já bastante difundida nos países com maior tradição florestal, onde constituí fator importante nas decisões de manejo e planejamento florestais e, portanto, na maximização do aproveitamento da madeira produzida. A utilização de modelos de afilamento é uma conseqüência natural do desenvolvimento do setor florestal no Brasil, em que o aperfeiçoamento das técnicas de inventário e a necessidade de flexibilização da informação dos estoques de madeira despertaram maior interesse no aprimoramento das informações de inventário (RIOS,1997).

No Brasil, as pesquisas que analisam o perfil do tronco de árvores e que convertem florestas em multiprodutos da madeira são mais direcionadas para Eucalyptus e Pinus. Quanto às espécies nativas, destacam-se os trabalhos de Garcia et al. (1993) e Moura (1994). Como a madeira de virola é utilizada para a fabricação de grandes, médias e pequenas peças, o sortimento possibilitará o direcionamento do estoque conforme as alternativas mais lucrativas, sendo possível maximizar financeiramente os estoques de virola existentes nas florestas de várzea do estuário amazônico.

Nesse sentido, o presente estudo teve como objetivo avaliar a acuracidade de alguns modelos de afilamento para quantificar o sortimento de Virola surinamensis.

2. MATERIAL E MÉTODOS

2.1. Área de estudo

A área de estudo foi uma floresta não explorada de várzea baixa localizada no rio Medonho, Município de Afuá, ao norte do Estado do Pará. O clima da região é quente e úmido do tipo A_mi, segundo Köeppen; a precipitação média anual, de 2.500 mm; e a temperatura média anual e a umidade relativa do ar correspondem a 26 ºC e 85%, respectivamente. O período mais seco vai de setembro a novembro e o mais chuvoso, de fevereiro a abril (SUDAM, 1984).

A região apresenta extensa rede hidrográfica, com topografia plana a suavemente ondulada. Os solos são originados de formação geológica do quaternário, e o material que os compõe é oriundo de sedimentos holocênicos, com elevada fertilidade natural e conteúdo de bases trocáveis variáveis; na classificação de Vieira (1988), são solos hidromórficos, gleizados e pouco húmicos.

2.2. Informações sobre a espécie

A virola é uma espécie da família Myristicaceae, de porte médio (até 40 m de altura e DAP < 1,0 m), fuste monopodial e copa cuneiforme. Seus ramos apresentam folhas alternas, verdes glabras, com base obtusa e ápice acuminado. A inflorescência é em forma de panículas axilares ou subaxilares, com pedicelos opostos lateralmente; possuí fascículos de 8 a 15 flores nas extremidades dos ramos; flores femininas raras, com ovário de forma ovóide e estilete curto; o estigma é emarginado, bífido e ereto. O fruto é em forma de elipse, com 14 a 16 mm de comprimento; curto estipilado, com frutificação de fevereiro a julho. O tronco apresenta ramificações regulares, verticilada, quase horizontal; a casca é espessa, esbranquiçada e de coloração castanha no interior (Van-ROOSMALEN et al., 1996).

A semente é recalcitrante, apresenta dispersão primária endozoocórica/barocórica e secundária hidrocórica. No Brasil, distribui-se desde a metade oriental da região amazônica, principalmente no estuário amazônico, Estados de Roraima, do Maranhão até Pernambuco, como também nas Güianas, na Venezuela, Colômbia e em Trinidad e Antilhas. É uma espécie que cresce nos solos inundáveis de várzea, com área de ocorrência natural, situando-se entre 0º e 5º sul (LOPES et al., 1997).

A madeira de virola é leve, com densidade em torno de 0,50 g/cm³; seu cerne varia de bege-claro a castanho-escuro; alburno bem desenvolvido, sem gosto e com cheiro distinto (OLIVEIRA, 1972). É uma das espécies mais processadas nas indústrias madeireiras do estuário amazônico, por ter boa valorização econômica no mercado internacional, sendo uma das espécies preferidas para a fabricação de compensado, embalagens, artigos de esporte, brinquedos, lápis, palitos, bobinas e carretéis, entre outros utensílios (SUDAM, 1981). No período de 1980-98, mesmo apresentando taxa decrescente de comercialização devido à política de contingenciamento, a virola foi a terceira espécie arbórea mais exportada para o mercado internacional de madeiras tropicais (ANGELO et al., 2001).

2.3. Coleta de dados e análises

Com a árvore em pé, foi medido o diâmetro à altura de 1,30 m (DAP), ou logo acima da sapopema, de 93 indivíduos com diâmetro variando de 14 a 80 cm. Posteriormente, as árvores foram derrubadas, e se realizaram as mensurações da altura de cepa, fuste (parte comercial) e copa; os diâmetros foram medidos com casca em 15 secções ao longo do fuste, nos pontos correspondentes a 1, 2, 3, 4, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90 e 100% da altura de fuste. Para a obtenção do volume individual, utilizou-se o método de Smalian. Foram testados os modelos, seguidos pelas respectivas expressões de cálculo de d, h e V:

a) Kozak et al. (1969)

b) Demaerschalk (1972)

c) Baldwin (1991)

d) Garay (1979)

em que:

em que ; d= diâmetro na altura h; altura total; Tx = 0 para d com casca e Tx = 1 para d sem casca; b_i = parâmetros de regressão, com i= 0,1,...,4; e=erro aleatório, sendo e ~ N(0,s²).

O volume, no modelo de Garay (1979), foi obtido mediante o procedimento de Smalian, utilizando-se os diâmetros estimados pela equação do d desse modelo, para alturas relativas observadas na cubagem. Em uma aplicação prática, recomenda-se adotar seções de 0,1 m para conseguir determinar o volume até um diâmetro preestabelecido. Outra alternativa consiste em estimar a altura h onde o diâmetro d ocorre e, então, estabelecer intervalos em torno de 1,0 m, para estimativa dos diâmetros comerciais e obtenção do volume via Smalian.

As equações ajustadas foram avaliadas por meio do coeficiente de determinação - R², coeficiente de correlação múltipla - e análise gráfica dos resíduos. Nos modelos não-lineares, o coeficiente de correlação foi obtido por:

em que:

= diâmetro estimado pela equação; e

Y_i = diâmetro observado.

Com exceção do modelo de Kozak, os ajustes foram feitos utilizando-se o procedimento Quase Newton, disponível no Software Statistica 5.0.

Além dos modelos de afilamento, foi ajustado um modelo hipsométrico, empregando-se a relação funcional do modelo de Chapman-Richards, da forma:

3. RESULTADOS E DISCUSSÃO

Ao analisar os coeficientes de correlação, entre valores estimados e valores observados, pode-se inferir, em princípio, que os modelos de Kozak e de Garay resultaram em maior precisão. Porém, a variável dependente desses modelos é diferente das variáveis dependentes dos demais (Quadro 1). Além disso, o modelo de Kozak é linear, enquanto os demais são não-lineares. Assim, as estatísticas apresentadas no Quadro 1 servem apenas para fins de comparação relativa com resultados de outros estudos.

Analisando as Figuras 1 e 2, verifica-se que alguns modelos resultam em estimativas inconsistentes de altura comercial, enquanto outros resultam em estimativas menos precisas de diâmetro comercial ou de volume. Isso, de certo modo, é esperado, uma vez que o modelo pode ser bom para estimar uma variável, porém deficiente para outra, ou seja, se um modelo estima d com precisão, não implica estimar h também com precisão. O fato é que se a equação estima o diâmetro ao longo do tronco com precisão e não-tendenciosidade, em geral ele é adequado para estimar volume.

Optou-se pelo modelo de Garay, uma vez que a equação resultante indicou estimativas consistentes das três variáveis implícitas no estudo (d, h e V), com precisão aceitável (Figuras 1 e 2).

Os coeficientes de correlação entre valores de d/dap observados e estimados foram relativamente baixos em comparação com coeficientes observados na literatura. Isso ocorreu devido à grande amplitude de variação observada nesses valores. Apesar disso, é importante verificar que a tendência de variação de h/Ht em relação a d/dap é bem definida (Figura 1).

As equações oriundas dos modelos de Kozak, Demaerschalk e Baldwin resultaram em estimativas inconsistentes de altura comercial (Figura 2). O único modelo que resultou em estimativas precisas e livres de viés foi o de Garay. O próprio Garay (1979) já havia demonstrado a eficiência desse modelo para diferentes espécies. Esse modelo, derivado da função Chapman-Richards, é flexível o suficiente para descrever as variações de forma em árvores de diferentes espécies e tamanhos e consistente em relação às variáveis d, h e V.

Para permitir um uso prático da equação obtida a partir do modelo selecionado, foi gerada uma relação hipsométrica utilizando o modelo de Chapman-Richards, sendo a equação resultante apresentada na Figura 3.

A partir do modelo de Garay (Quadro 1), segue-se que:

Substituindo a relação hipsométrica da Figura 3 nessas duas últimas expressões, tem-se, respectivamente:

As expressões 1 e 2 permitem estimar o diâmetro, d, em qualquer altura h, ou a altura comercial, h, para determinado limite de diâmetro comercial, d, utilizando como input apenas o dap. Aplicando essas expressões, na faixa de variação de dap observada nos dados foram obtidas as curvas apresentadas na Figura 4. Essa figura contém as variações de forma da espécie estudada e pode ser utilizada para cubar árvores que tiveram apenas o dap medido no campo. O uso desse tipo de gráfico, em ambiente Excel, permite uma melhor identificação do perfil do tronco em árvores de diferentes diâmetros (daps).

4. CONCLUSÕES

A equação:

d = dap1,86696(1+0,14998Ln(1-0,99554h ^0,013315Ht ^-0,013315))obtida a partir do ajuste do modelo de Garay (1979) é indicada para quantificar multiprodutos de árvores de Virola surinamensis. A eficiência de um modelo de afilamento deve ser avaliada com base nas variáveis d e h.

5. REFERÊNCIAS BIBLIOGRÁFICAS

Recebido em 1º.06.2004 e aceito para publicação em 10.11.2005.

Datas de Publicação

Histórico