Acessibilidade / Reportar erro

Modelagem da biomassa total e da lenha por unidade de área para bracatingais nativos

Modelling of total and firewood biomass by unnit of area for native mimosa scabrella Benth. Stands

Resumos

A bracatinga (Mimosa scabrella Benth.) é uma espécie pioneira nativa da Floresta Ombrófila Mista de rápido crescimento e cultivada, predominantemente, na forma de povoamentos puros ou consorciada com culturas agrícolas pelos pequenos e médios produtores da região metropolitana de Curitiba. O objetivo deste estudo foi selecionar equações estimativas da biomassa total e da lenha de bracatingais por unidade de área nas diferentes idades, disponibilizando, assim, uma ferramenta simples e de fácil aplicação durante o processo de comercialização dos povoamentos. Os dados são provenientes de 272 parcelas temporárias com idades variando de 3 a 18 anos. Foram testados 21 modelos tradicionais: 10 aritméticos e semilogarítmicos e 11 logarítmicos. Foram também desenvolvidas equações pelo processo Stepwise a partir de uma matriz de correlação. As equações ajustadas foram comparadas pelo Coeficiente de determinação ajustado (R²aj), Erro-padrão da estimativa percentual (Syx%), teste F e Distribuição gráfica de resíduos. Duas equações para cada caso foram selecionadas, sendo uma tradicional e a outra desenvolvida por Stepwise. Para a biomassa total, foram selecionadas a equação de Clutter e a equação aritmética desenvolvida por Stepwise; e para a biomassa da lenha, a equação da variável combinada, proposta por Spurr, e a equação aritmética de Stepwise. Por fim, realizou-se teste Qui-quadrado (x²) a partir de uma amostra de 17 parcelas extraídas da base de dados, em que as quatro equações selecionadas foram consideradas válidas para a estimativa da biomassa total e da lenha dos bracatingais da região metropolitana de Curitiba.

Bracatinga; Modelos estimativos; Teste de validação


Mimosa scabrella Benth. is a rapid-growth native pioneer specie from the Ombrophylous Mix Forest biome. It is cultivated predominantly in pure stand form or intercropped with agricultural crops by small and medium producers in the metropolitan region of Curitiba. The objective this study was to select estimative models for total and firewood biomass of M. scabrella stands per unit of area at different ages. It was data were collected from 272 temporary samples plots aged 3 to 18 years old. We tested 21 traditional models: 10 arithmetic and semi-logarithmic and 11 logarithmic; it was also developed equations by the Stepwise process from a correlation matrix. The adjusted equation were compared basing on R²aj, Syx%, F value and on graphic analysis of residuals. Two equations for each case were selected, being one traditional and other developed by Stepwise. The Clutter equation and arithmetic developed by Stepwise were selected for total biomass. For firewood, were selected the combined variable equation proposed by Spurr, and the arithmetic equation developed by Stepwise. The validation of equations was performed by Chi-square test from 17 plots randomly selected from the data base, in which four equations were considered valid to estimate total and firewood biomass for M. scabrella stands in the metropolitan region of Curitiba.

Mimosa scabrella; Estimative models; Validation test


Modelagem da biomassa total e da lenha por unidade de área para bracatingais nativos

Modelling of total and firewood biomass by unnit of area for native mimosa scabrella Benth. Stands

Ronan Felipe de SouzaI; Sebastião do Amaral MachadoII; Afonso Figueiredo FilhoIII

I Programa de Pós-graduação em Engenharia Florestal, Universidade Federal do Paraná. E-mail:<ronanflorestal@gmail.com>

II Departamento de Ciências Florestais da Universidade Federal do Paraná, Paraná, Brasil. E-mail:<samachado@ufpr.br>

III Universidade Estadual do Centro-Oeste, UNICENTRO, Brasil. E-mail:<afig@floresta.ufpr.br>

RESUMO

A bracatinga (Mimosa scabrella Benth.) é uma espécie pioneira nativa da Floresta Ombrófila Mista de rápido crescimento e cultivada, predominantemente, na forma de povoamentos puros ou consorciada com culturas agrícolas pelos pequenos e médios produtores da região metropolitana de Curitiba. O objetivo deste estudo foi selecionar equações estimativas da biomassa total e da lenha de bracatingais por unidade de área nas diferentes idades, disponibilizando, assim, uma ferramenta simples e de fácil aplicação durante o processo de comercialização dos povoamentos. Os dados são provenientes de 272 parcelas temporárias com idades variando de 3 a 18 anos. Foram testados 21 modelos tradicionais: 10 aritméticos e semilogarítmicos e 11 logarítmicos. Foram também desenvolvidas equações pelo processo Stepwise a partir de uma matriz de correlação. As equações ajustadas foram comparadas pelo Coeficiente de determinação ajustado (R²aj), Erro-padrão da estimativa percentual (Syx%), teste F e Distribuição gráfica de resíduos. Duas equações para cada caso foram selecionadas, sendo uma tradicional e a outra desenvolvida por Stepwise. Para a biomassa total, foram selecionadas a equação de Clutter e a equação aritmética desenvolvida por Stepwise; e para a biomassa da lenha, a equação da variável combinada, proposta por Spurr, e a equação aritmética de Stepwise. Por fim, realizou-se teste Qui-quadrado (x²) a partir de uma amostra de 17 parcelas extraídas da base de dados, em que as quatro equações selecionadas foram consideradas válidas para a estimativa da biomassa total e da lenha dos bracatingais da região metropolitana de Curitiba.

Palavras-chave: Bracatinga; Modelos estimativos; Teste de validação.

ABSTRACT

Mimosa scabrella Benth. is a rapid-growth native pioneer specie from the Ombrophylous Mix Forest biome. It is cultivated predominantly in pure stand form or intercropped with agricultural crops by small and medium producers in the metropolitan region of Curitiba. The objective this study was to select estimative models for total and firewood biomass of M. scabrella stands per unit of area at different ages. It was data were collected from 272 temporary samples plots aged 3 to 18 years old. We tested 21 traditional models: 10 arithmetic and semi-logarithmic and 11 logarithmic; it was also developed equations by the Stepwise process from a correlation matrix. The adjusted equation were compared basing on R²aj, Syx%, F value and on graphic analysis of residuals. Two equations for each case were selected, being one traditional and other developed by Stepwise. The Clutter equation and arithmetic developed by Stepwise were selected for total biomass. For firewood, were selected the combined variable equation proposed by Spurr, and the arithmetic equation developed by Stepwise. The validation of equations was performed by Chi-square test from 17 plots randomly selected from the data base, in which four equations were considered valid to estimate total and firewood biomass for M. scabrella stands in the metropolitan region of Curitiba.

Keyworks: Mimosa scabrella; Estimative models; Validation test.

1. INTRODUÇÃO

A bracatinga (Mimosa scabrella Benth) é uma espécie pioneira, heliófila, nativa da Floresta Ombrófila Mista, com dominância de Araucaria angustifolia. Em 1909, Romário Martins iniciou uma campanha de divulgação para estabelecer povoamentos de bracatinga com fins de produção de lenha na região metropolitana de Curitiba, devido ao rápido crescimento e à simplicidade de práticas culturais necessárias para o plantio e manutenção desta espécie (CARPANEZZI et al., 1997).

Laurent et al. (1990) estimaram a área plantada de bracatinga na região metropolitana de Curitiba em 60.000 ha. Atualmente, a área plantada com essa espécie tem diminuído, embora o sistema tradicional de cultivo ainda seja atividade de grande importância econômica para pequenos e médios produtores da região. Dessa forma, estudos que permitem o avanço das técnicas silviculturais e de manejo da espécie, bem como a estimativa acurada da produção, são de grande importância, contribuindo para a ampliação do retorno financeiro aos produtores.

Além da comercialização na forma do volume de madeira em metros cúbicos (m³), os bracatingais também são muito comercializados pelo peso de matéria seca dos compartimentos passíveis de utilização na forma de energia (lenha, carvão). Assim, a estimativa acurada da biomassa total e da lenha por unidade de área é de grande valia, permitindo a quantificação dos valores referentes à madeira que será comercializada, bem como valores de biomassa seca remanescente nos povoamentos que será incorporada ao solo.

Oliveira et al. (2005) desenvolveram equações estimativas do volume por unidade de área para fragmentos de floresta natural da região de Viçosa, Minas Gerais. Em seus resultados, esses autores destacaram a alta correlação do volume com as variáveis independentes área basal por ha (G), altura total média () e altura dominante (hdom), o que resultou na seleção de equações que envolveram essas variáveis como as mais acuradas. Analisando o efeito da inclusão da altura do fuste na forma logaritmizada, esses mesmos autores encontraram resultados satisfatórios, observando valores de coeficiente de determinação superiores a 0,95 para as equações.

Machado et al. (2008) desenvolveram equações estimativas para o volume por unidade de área para bracatingais da região metropolitana de Curitiba. Os resultados revelaram alta correlação das variáveis combinadas G e Ghdom com o volume total e da lenha nas formas aritmética e logarítmica, sendo estas as variáveis incluídas na equação Stepwise, desenvolvida pelos referidos autores. Dos modelos aritméticos e semilogarítmicos selecionados da literatura, os que apresentaram melhores ajustes foram aqueles que envolveram a variável G, sendo considerados inadequados os que não compunham essa variável. As equações logarítmicas apresentaram melhores ajustes do que as aritméticas, sempre exibindo em sua estrutura as variáveis G e Ghdom na forma logaritmizada.

Objetivando selecionar equações capazes de estimar com precisão a biomassa total e da lenha por unidade de área nas diferentes idades dos bracatingais nativos e manejados de forma tradicional da região metropolitana de Curitiba, foram ajustados e comparados modelos tradicionais aritméticos, semilogarítmicos e logarítmicos. Também foram desenvolvidas e testadas equações pelo processo Stepwise a partir de uma matriz de correlação linear simples.

2. MATERIAL E MÉTODOS

2.1. Caracterização da área de estudo

No total foram incluídos 14 municípios, abrangendo uma área de influência de 731.850 ha. Machado et al. (2008) destacaram que esta área de estudo está compreendida entre Tunas do Paraná como a localidade mais ao Norte e Agudos do Sul no extremo Sul. Campo Largo é a cidade mais a Oeste e Campina Grande do Sul, mais ao Leste. A altitude varia de 850 m em Agudos do Sul até 980 m acima do nível do mar, em Bocaiuva do Sul.

A unidade fitogeográfica é caracterizada como transicional entre a Floresta Ombrófila Densa e a Floresta Ombrófila Mista, com dominância de A. angustifolia. O clima dominante da região, segundo a classificação de Köppen, é do tipo Cfb, com precipitação em todas as estações; a temperatura média do mês mais quente é inferior a 22 °C. Os solos, em sua maioria, são classificados como Cambissolos álicos e apresentam-se pobres, pouco desenvolvidos, ácidos, com teores elevados de alumínio e baixa saturação de bases (BAGGIO et al., 1995).

2.2. Base de dados

A base inicial de dados foi composta por 267 parcelas temporárias, instaladas em bracatingais nativos e manejados de forma tradicional na região metropolitana de Curitiba pela Embrapa – Florestas, Companhia de Energia Elétrica do Paraná – Copel, Rochadelli et al., na década de 1990, e, mais recentemente, por Machado et al., em 2006. Outras 40 parcelas foram instaladas por ocasião deste trabalho, totalizando 307 parcelas. A superfície amostral variou de 100 m² a 400 m² e as idades, de 3 a 18 anos. Nas unidades amostrais foram mensuradas as variáveis altura total, altura de fuste e circunferência à altura do peito de todos os indivíduos, com vara graduada e fita métrica, respectivamente.

As estimativas da biomassa (peso seco) total e da lenha de cada indivíduo dentro das parcelas foram feitas através das respectivas equações desenvolvidas por Urbano et al. (2008), as quais podem ser observadas a seguir com as respectivas estatísticas de ajuste, (1) e (2). A soma dos pesos de todos os indivíduos de cada parcela resultou na biomassa para as respectivas unidades de área. Esses valores foram extrapolados para hectare, tal como se usa nos inventários volumétricos convencionais partindo de equações de volume.

BTi = -11,8301 + 0,50917 dap²

aj = 0,9600 Syx% = 21,68 (1)

BLi = -6,6694 + 0,209204 dap² + 0,012728dap²h

aj = 0,9590 Syx% = 23,47 (2)

em que:

BTi = biomassa total individual (kg);

BLi = biomassa da lenha individual (kg);

dap = diâmetro à altura do peito (cm); e

h = altura total (m).

Foram também determinados a área basal e o número de árvores por hectare em cada unidade amostral, bem como idade, índice de sítio, altura total média, altura dominante, altura de fuste média, diâmetro médio e diâmetro quadrático de cada parcela.

A detecção de valores discrepantes foi realizada seguindo três critérios: quando foi possível afirmar que houve erro de medição; quando em uma parcela foram identificados valores negativos de biomassa total e da lenha; e quando foram observados valores muito superiores e inferiores nas diferentes classes de idade. Nesse procedimento, foram identificadas 35 parcelas discrepantes, as quais foram eliminadas da base de dados. Das 272 parcelas restantes, retiraram-se 17 parcelas aleatoriamente nas diferentes idades e classes de sítio, para validação dos modelos ajustados selecionados. A base de dados final, utilizada para o ajuste dos modelos estimativos da biomassa total e da lenha por unidade de área, foi, então, de 255 parcelas, para a qual se procedeu a uma análise estatística de cada variável, bem como ao cálculo do número de amostras necessárias em função da variância da população. A precisão estipulada foi de 10%, com probabilidade fixada em 95%.

2.3. Matriz de correlação linear simples

Uma matriz de correlação linear simples foi construída para analisar a relação das variáveis independentes com a biomassa total e da lenha, nas suas formas puras e logaritmizadas. O conjunto de variáveis independentes totalizou 85 variáveis, entre variáveis puras, inversas, combinadas, potencializadas e logaritmizadas.

2.4. Modelos tradicionais testados

Foram testados 10 modelos aritméticos e semilogarítmicos tradicionais obtidos na literatura, identificados de a1 até a10, destacando-se entre eles os modelos propostos por Spurr (1952), dos quais um foi adaptado da forma alemã para volume individual (a2) e o outro composto pela altura dominante (a3), muito utilizado devido à rapidez e ao baixo custo na coleta de dados em campo. Os modelos logarítmicos testados totalizaram 11 modelos, identificados de ln1 a ln11; destacou-se entre esses modelos aquele proposto por Clutter (ln11), devido à sua ampla utilização em estudos de crescimento e produção, representando ferramenta muito útil para o manejador florestal.

Visando à correção da discrepância logarítmica para aqueles modelos logaritmizados, foi efetuado o recálculo das estatísticas de ajuste, aplicando-se o Índice de Meyer (IM).

2.5. Desenvolvimento de equações por Stepwise

Para cada uma das variáveis, duas equações foram desenvolvidas pelo processo Stepwise: uma com a variável dependente na forma pura e outra na forma logaritmizada. Nesse procedimento, as variáveis são inseridas em ordem decrescente de correlação linear simples com a variável dependente e, quando não resultam em melhora significativa, são eliminadas. O nível de significância utilizado foi de 95% para a construção das equações por Stepwise.

2.6. Seleção de equações

A melhor equação estimativa para cada variável dependente foi selecionada comparando-se o coeficiente de determinação ajustado (R²aj), erro-padrão da estimativa percentual (Syx%), teste F e distribuição gráfica de resíduos percentuais. Para aquelas equações em que a variável dependente foi logaritmizada ou transformada, efetuou-se o cálculo do Índice de Schlaegel (IA), compatível com o valor de coeficiente de determinação das demais equações.

2.7. Validação das equações selecionadas

A validação permite ao pesquisador analisar a aplicabilidade da eqação estimativa desenvolvida e selecionada como de maior acurácia. Para que fossem consideradas válidas as equações selecionadas para cada caso, adotou-se o teste de aderência de Quiquadrado (x²), a 95% de significância. A análise da distribuição de resíduos para as 17 parcelas retiradas da base de dados foi realizada, servindo também como critério de validação.

3. RESULTADOS

3.1. Caracterização da base de dados

Os resultados obtidos de suficiência amostral indicaram que a base de dados utilizada foi eficiente para a estimativa de todas as variáveis. O número de árvores foi responsável pelo maior número de amostras necessárias, com 205 amostras. Para todas as variáveis dependentes, o erro de amostragem relativo foi inferior a 1%, o número de árvores foi responsável pelo menor erro, com 4,18%; para as demais variáveis independentes, os valores de erro de amostragem relativo foram superiores a 20%, com valor máximo para o diâmetro médio (40%).

O diâmetro médio apresentou menor valor de 5,23 cm aos 4 anos no sítio III e maior valor (21,59 cm) em um povoamento de 18 anos no sítio I, com coeficiente de variação igual a 32,44%. A altura total média teve seu valor mínimo igual a 6,91 m aos 3,5 anos no sítio III e valor máximo de 20,08 m na idade de 18 anos no sítio I, com coeficiente de variação de 19,71%, o menor entre as variáveis.

A área basal apresentou coeficiente de variação de 34,12%. O valor mínimo foi de 2,06 m²/ha, observado em um povoamento de 14 anos no sítio III, e o valor máximo foi de 24,09 m²/ha, observado no sítio II, em um povoamento de 6 anos de idade. A menor densidade de árvores (67 árvores/ha) foi observada em um povoamento com idade de 14 anos no sítio III e a maior (6.161 árvores/ha), em um povoamento de 5 anos de idade, no sítio I. O coeficiente de variação desta variável foi de 72,72%, o mais elevado entre as variáveis.

A biomassa total teve o menor valor estimado no sítio III em povoamento com 14 anos, com valor de 12.469,77 kg/ha; e maior valor (113.195,60 kg/ha) para um povoamento de 7 anos, no sítio I. O coeficiente de variação foi de 36,92% e a amplitude, igual a 100.725,83 kg/ha. A biomassa da lenha apresentou coeficiente de variação igual a 38,63%. O menor valor foi de 9.693,05 kg/ha, estimado para um povoamento de 14 anos no sítio III e maior valor no sítio I, em um povoamento de 10 anos, igual a 94.373,41 kg/ha, apresentando amplitude de 84.680,36 kg/ha.

3.2. Matriz de correlação linear simples

Na Tabela 1 é apresentada a relação das variáveis independentes que apresentaram correlação superior a 0,9 para pelo menos uma das variáveis dependentes. Foram destacadas as duas maiores correlações de cada variável.

3.3. Equações estimativas da biomassa total

Os coeficientes ajustados com as respectivas estatísticas de ajuste e precisão para as equações tradicionais estimativas da biomassa total são apresentados na Tabela 2. Em geral, os resultados das equações foram satisfatórios, com altos valores de R² ajyx e baixos valores de Sxy %.

No desenvolvimento da equação aritmética pelo processo Stepwise, a variável foi a primeira a ingressar, nesta sequência: ln, N² e . A equação final ajustada, representada pela equação (3), apresentou R²aj igual a 0,9932, Syx% igual a 3,03 e valor de F igual a 9351.

em que:

BT = biomassa total (kg/ha);

G = área basal (G/ha);

= diâmetro médio à altura do peito (cm); e

N = número de indivíduos (N/ha).

No caso da equação logarítmica desenvolvida por Stepwise, o valor R²aj para a equação final representado pela equação (4) foi igual a 0,9921, com Syx% de 3,28 e valor de F igual a 8213. O Índice de Meyer calculado para essa equação foi de 1,0011. As variáveis independentes que compuseram a equação final ajustada foram ln, N², GN e ².

em que:

lnBT = logaritmo neperiano da biomassa da lenha (kg/ha);

G = área basal (m²/ha); e

= diâmetro médio à altura do peito (cm).

Na Figura 1, pode ser observada a distribuição de resíduos para as equações a6 e ln11, selecionadas entre as aritméticas como as mais acuradas, e para as equações semilogarítmica e logarítmica, desenvolvidas por Stepwise. A distribuição de resíduos para as quatro equações foi homogênea, independente e pouco dispersa ao longo da linha estimativa, o que dificultou a comparação gráfica delas. As demais equações tradicionais testadas apresentaram superestimativas para valores inferiores a 50.000 kg/ha de biomassa total.


A equação ln11 foi selecionada como a mais acurada entre as tradicionais, pela superioridade que apresentou nas estatísticas de ajuste e precisão, sendo graficamente muito semelhante à equação a6. Entre as equações de Stepwise, graficamente também não se observou diferença, sendo a semilogarítmica considerada melhor por causa do melhor resultado encontrado para R²aj, Syx % e valor de F.

As equações selecionadas foram validadas pelo teste de aderência de Qui-quadrado, com valores de x² calculado iguais a 1,2145 e 0,0374, respectivamente para a equação ln11 e Stepwise. O valor tabelado de x² foi de 26,30. A dispersão de resíduos da biomassa total para as 17 unidades amostrais utilizadas na validação das equações apresentou-se homogênea, independente e pouco dispersa ao longo da linha estimativa de ambas as equações, atestando, assim, a aderência desses resíduos para as estimativas.

3.4. Equações estimativas da biomassa da lenha

Em geral, os resultados do ajuste dos modelos tradicionais para estimativa da biomassa da lenha foram satisfatórios. Elevados valores de R²aj e baixos valores de S yx% foram observados tanto nas equações aritméticas e semilogarítmicas quanto nas equações logarítmicas, salvo exceções, como pode ser observado na Tabela 3.

No desenvolvimento da equação aritmética pelo processo Stepwise, a variável independente que primeiro ingressou foi G, seguida por N e G². A equação final ajustada, representada pela equação (5), apresentou resultado satisfatório, com R²aj igual a 0,9965, Syx% de 2,27 e valor de F igual a 24449.

em que:

BL = biomassa da lenha (kg/ha);

G = área basal (m²/ha);

= altura total média (m); e

N = número de indivíduos (N/ha).

A equação logarítmica final, desenvolvida pelo processo Stepwise para a estimativa da biomassa da lenha, é representada pela equação 6; para essa equação, o valor de R²aj foi igual a 0,9934, Syx% de 3,12 e valor de F de 13100. O Índice de Meyer calculado foi de 1,0006. A variável combinada lnG foi a primeira variável incorporada, seguida por -1 e ln .

em que:

lnBL = logaritmo neperiano da biomassa da lenha (kg/ha);

G = área basal (m²/ha);

= altura total média (m); e

= diâmetro médio à altura do peito (cm).

Na Figura 2 pode ser observada a distribuição percentual de resíduos da biomassa da lenha para as equações a2 e ln2, selecionadas entre as equações tradicionais aritméticas e logarítmicas, e para as equações desenvolvidas por Stepwise, para ambos os casos.


Da mesma forma que na biomassa total, não foi possível identificar graficamente diferenças entre as equações selecionadas, e todas apresentaram homogeneidade, independência e baixa dispersão de resíduos ao longo da linha estimativa. Foram consideradas mais acuradas, então, aquelas equações que apresentaram melhores estatísticas de ajuste e precisão, resultando nas equações a2 e aritmética de Stepwise.

Ambas as equações foram consideradas válidas pelo teste de Qui-quadrado, com valores de x² iguais a 0,6616 e 0,0312, respectivamente para a equação a2 e Stepwise. O valor tabelado de x² foi igual a 26,30. A análise gráfica das 17 parcelas separadas da base de dados comprovou também a validade dessas equações, em que a distribuição de resíduos foi homogênea, independente e pouco dispersa ao longo da linha estimativa.

4. DISCUSSÃO

4.1.Análise da base de dados

Diante dos resultados da caracterização da base de dados, é possível afirmar que essa caracterização se apresentou sólida, representativa das diferentes combinações de idades e sítios existentes na região metropolitana de Curitiba. Destaca-se, porém, o caso da área basal, que apresentou o maior valor para um povoamento de 6 anos no sítio II e que foi justificado pelo elevado número de indivíduos e diâmetro médio constatados nas unidades amostrais dessa combinação.

A matriz de correlação linear simples proporcionou resultados importantes para o processo de seleção das equações, permitindo a detecção rápida de quais variáveis independentes foram altamente correlacionadas com as dependentes. Estudos relacionados à estimativa do volume por unidade de área têm utilizado essa ferramenta para a mesma finalidade, como pode ser observado nos trabalhos de Oliveira et al. (2005) e Machado et al. (2008).

A maior correlação com a combinação G foi observada apenas para a biomassa da lenha; para a biomassa total, a melhor relação foi com a combinação G². Esse resultado corrobora o do trabalho apresentado por Urbano (2008), que destacou a maior relação da biomassa total individual com o diâmetro à altura do peito.

4.2. Seleção de equações

Para a biomassa total, como mencionado anteriormente, as equações selecionadas foram a equação logarítmica de Clutter (ln11) e a equação aritmética desenvolvida por Stepwise.

O resultado da equação de Clutter merece destaque, por ser aplicado em estudos de crescimento e produção e pelo fato de envolver as variáveis idade, sítio e área basal, as quais são de fácil obtenção em campo. Sobre as equações desenvolvidas por Stepwise, destaca-se a menor correlação linear simples das variáveis que envolveram a altura, induzindo que estas não fossem incluídas nas equações. Esse resultado difere dos encontrados por Silva (1979), Figueiredo Filho (1982), Rosot (1989), Oliveira et al. (2005) e Machado et al. (2008), nos quais as equações de Stepwise para a estimativa do volume por unidade de área envolveram sempre variáveis combinadas à altura.

Urbano et al. (2008), quando desenvolveram equações para a biomassa total individual da bracatinga, também selecionaram e indicaram para uso a equação que envolveu apenas o diâmetro. A não inclusão da variável altura nas equações ajustadas neste trabalho é explicada pelo fato de a biomassa total envolver todos os compartimentos das árvores e não apenas o compartimento fuste, como ocorre nas estimativas do volume.

Os resultados do ajuste das equações tradicionais e do processo Stepwise para estimativa da biomassa da lenha seguiram a tendência da matriz de correlação, sendo G a primeira variável a ser incluída na equação. A equação tradicional selecionada foi a equação a2, proposta por Spurr e denominada como o modelo da variável combinada. Pelo processo Stepwise, a equação selecionada foi a aritmética.

O modelo proposto por Clutter, ajustado para a estimativa da biomassa da lenha, apresentou R²aj superior a 0,90 e Syx% de 11,64. Resultado esse que merece destaque, uma vez que possibilita a estimativa da biomassa da lenha a partir da idade e do sítio, sendo necessária apenas a determinação da área basal nas diferentes idades.

Neste trabalho, os modelos logarítmicos não apresentaram resultados superiores aos aritméticos, contrariando os resultados encontrados por Silva (1979), Figueiredo Filho (1982), Rosot (1989), Oliveira et al. (2005) e Machado et al. (2008), os quais destacaram que sempre os modelos logarítmicos apresentaram maior acurácia que os modelos aritméticos. Sobre esse aspecto, Machado et al. (2002) citaram que é recomendável testar vários modelos, identificando por meio de análises estatísticas o melhor para cada caso.

5. CONCLUSÕES

As equações estimativas por unidade de área apresentaram-se muito eficientes e acuradas para estimar a biomassa total e da lenha dos bracatingais da região metropolitana de Curitiba, nas diferentes idades.

A equação de Clutter ajustada apresentou elevada acurácia para ambas as variáveis dependentes estimadas e representa alternativa importante durante o processo de comercialização dos bracatingais, permitindo a utilização do método de Bitterlich. Esse método, por sua vez, permite a redução considerável de tempo e custos no inventário florestal.

6. AGRADECIMENTOS

À EMBRAPA-FLORESTAS, COPEL e ao pesquisador Rochadelli, pela disponibilização de dados; e ao CNPq, pelo auxílio financeiro durante a realização dos trabalhos de campo.

7. REFERÊNCIAS

Recebido em 25.08.2010 aceito para publicação em 07.11.2013.

  • BAGGIO, A. J.; CARPANEZZI, A. A.; SANMIGUEL, A. A. Equações para a estimativa de peso da biomassa aérea de bracanga (Mimosa scabrella Benth.) na idade de corte. Boletim de Pesquisa Florestal, n.30/31, p.37-49, 1995.
  • CARPANEZZI, A. A.; CARPANEZZI, O. T. B.; BAGGIO, A. J. Manejo de bracatingais. In:GALVÃO, A. P. M., (Org.). Tópicos em manejo florestal sustentável Curitiba: Embrapa/Cnpf, 1997. v. único. p.157-163.
  • FIGUEIREDO FILHO, A.; MACHADO, S. A.; ROSOT, N. C. Estimativas volumétricas porUnidade de Área em Florestas Naturais. Silvicultura Em São Paulo, v. 16, p. 491504, 1982.
  • LAURENT, J. M. E. Sistema agroflorestal da barcatinga na região metropolitana de Curitiba: contexto sócio-economico. Curitiba: FAO Projeto GOP/BRA/025/FRA, 1990. 50p.
  • MACHADO, S. A.; CONCEIÇÃO, M. B.; FIGUEIREDO FILHO, A. Modelagem do volume individual para diferentes idades e regimes de desbaste em plantações de Pinus oocarpa Guarapuava. Revista Ciências Exatas e Naturais, v.4, n.2, p.41-50, 2002.
  • MACHADO, S. A.; PROFUMO AGUIAR, L.; FIGUEIREDO FILHO, A.; KOEHLER, H. S. Modelagem do volume do povoamento para mimosa scabrella benth. na região metropolitana de Curitiba. Revista Árvore, v.32, n.3, p.465478, 2008.
  • OLIVEIRA, M. L. R. et al. Equações de volume de povoamento para fragmentos florestais naturais do Município de Viçosa, Minas Gerais. Revista Árvore, v.29, n.2, p.213-225, 2005.
  • ROSOT, M. A. D. Estudo comparativo de métodos para a avaliação volumétrica por unidade de área em um povoamento de Pinus taeda L. 1989. 163f. Dissertação (Mestrado em Engenharia Florestal) Setor de Ciências Agrárias, Universidade Federal do Paraná, Curitiba, 1989.
  • SILVA, J. A. Estimativa do volume por hectare e sua composição em povoamentos de Eucalyptus spp. nos cerrados em Minas Gerais e Mato Grosso do Sul. 1979. 85f. Dissertação (Mestrado em Engenharia Florestal) - Setor de Ciências Agrárias, Universidade Federal do Paraná, Curitiba, 1979.
  • URBANO, E. et al. Modelagem da biomassa aérea em Bracatingais Nativos da Região Metropolitana de Curitiba. Revista Floresta, v.38, n.2, p.361372, 2008.

Datas de Publicação

  • Publicação nesta coleção
    14 Abr 2014
  • Data do Fascículo
    Fev 2014

Histórico

  • Recebido
    25 Ago 2010
  • Aceito
    07 Nov 2013
Sociedade de Investigações Florestais Universidade Federal de Viçosa, CEP: 36570-900 - Viçosa - Minas Gerais - Brazil, Tel: (55 31) 3612-3959 - Viçosa - MG - Brazil
E-mail: rarvore@sif.org.br