Projeção da dinâmica da floresta natural de Terra-firme, região de Manaus-AM, com o uso da cadeia de transição probabilística de Markov

Dynamics of the Terra-firme primary forest in Manaus-AM region using the Markov probabilistic transition matrix

Resumos

To combine protection and utilization of forest resources in the tropics, the understanding of forest dynamics is essential. It is also important in the definition of strategies for rehabilitation of degraded areas. In Forestry, forest dynamics could be translated as the understanding of recruitment, mortality and biomass increment rates over time. For this study, these rates were estimated based on measurements carried out in 2000 and 2004 over two transects measuring 20 by 2500 m (5 hectares) each, in Manaus region. This paper deals with forest dynamics of a pristine forest based on the probabilistic transition matrix (the first-order Markov Chain) approach. The main objective is to report 4-year (2000 to 2004) changes in the forest structure. Diameter distribution and tree mortality will be projected ahead to 2008 (t+2), based upon a 4-year period of observations completed in 2004 (t+1) and its immediate past in 2000 (t). In terms of fresh aboveground biomass, this site accumulated 8.34 t.ha-1.ano-1. The chi ² test has shown no statistical difference (p = 0.01) between observed diameter frequency and the expected projected by Markov Chain. This result indicates that the Markov Chain approach is a reliable tool to project the forest dynamics on a short-term basis. In 2008, the total number of individuals will have a decrease of 2.7%, and the mortality rate will 15% higher than in 2004.

Biomass; Mortality; Recruitment; Increment

CIÊNCIAS FLORESTAIS

Projeção da dinâmica da floresta natural de Terra-firme, região de Manaus-AM, com o uso da cadeia de transição probabilística de Markov

Dynamics of the Terra-firme primary forest in Manaus-AM region using the Markov probabilistic transition matrix

Liliane Martins Teixeira; Jeffrey Q. Chambers; Anabel Rodrigues e Silva; Adriano José Nogueira Lima; Vilany Matilla Colares Carneiro; Joaquim dos Santos; Niro Higuchi^I

^IInstituto Nacional de Pesquisas da Amazônia - Av. André Araújo, 2936 - Aleixo. CEP 69060-001 - Manaus/AM - CP 478. e-mail: [e-mail]

RESUMO

Os fatores que envolvem os processos da dinâmica da floresta influenciam a sua biodiversidade e, portanto, a qualidade da floresta. A definição de estratégias que envolve a proteção e o uso adequado da floresta manejada e a recuperação de áreas já degradadas tornam-se possível com o estudo da estrutura e dinâmica da floresta primária por meio de informações como a mortalidade, o recrutamento e a permanência das árvores no sistema florestal. Este trabalho teve como objetivo avaliar a dinâmica de uma floresta não perturbada e fazer projeções da dinâmica florestal usando a matriz de transição probabilística (Cadeia de Markov). As taxas de recrutamento, mortalidade e incremento foram determinadas a partir de inventários florestais realizados em dois transectos, nos sentidos Norte-Sul e Leste-Oeste (20 x 2500 m cada, totalizando 10 ha), localizados no km 50 da BR 174, na estrada vicinal ZF-2, Manaus/AM, nos anos de 2000 e 2004. A floresta acumulou 8,34 t.ha^-1.ano^-1 de biomassa fresca acima do solo. De acordo com projeção para 2008, o número total de árvores diminuirá em 2,67% (de 5987 indivíduos (2004) para 5827 (2008)) e a mortalidade será 15% maior (de 264 (2004) para 311 (2008)). O teste Qui-quadrado mostrou que não há diferença significativa (1% de probabilidade) entre as informações coletadas e projetadas. Esses resultados permitem concluir que a Cadeia de Markov é um eficiente instrumento para projetar a dinâmica da floresta natural, contribuindo para o planejamento em curto prazo das atividades que utilizam os recursos florestais.

Palavras-chave: Biomassa, Recrutamento, Mortalidade, Incremento.

ABSTRACT

To combine protection and utilization of forest resources in the tropics, the understanding of forest dynamics is essential. It is also important in the definition of strategies for rehabilitation of degraded areas. In Forestry, forest dynamics could be translated as the understanding of recruitment, mortality and biomass increment rates over time. For this study, these rates were estimated based on measurements carried out in 2000 and 2004 over two transects measuring 20 by 2500 m (5 hectares) each, in Manaus region. This paper deals with forest dynamics of a pristine forest based on the probabilistic transition matrix (the first-order Markov Chain) approach. The main objective is to report 4-year (2000 to 2004) changes in the forest structure. Diameter distribution and tree mortality will be projected ahead to 2008 (t+2), based upon a 4-year period of observations completed in 2004 (t+1) and its immediate past in 2000 (t). In terms of fresh aboveground biomass, this site accumulated 8.34 t.ha^-1.ano^-1. The c ² test has shown no statistical difference (p = 0.01) between observed diameter frequency and the expected projected by Markov Chain. This result indicates that the Markov Chain approach is a reliable tool to project the forest dynamics on a short-term basis. In 2008, the total number of individuals will have a decrease of 2.7%, and the mortality rate will 15% higher than in 2004.

Keywords: Biomass, Mortality, Recruitment, Increment

INTRODUÇÃO

Estudar a dinâmica da floresta tropical úmida amazônica, manejada ou não, é um grande desafio para os florestais. Os modelos clássicos de produção florestal foram desenvolvidos para florestas temperadas e têm como principais variáveis, o índice de sítio e idade da árvore ou do povoamento (Sullivan & Clutter, 1972; Ferguson & Leech, 1978; Alder, 1980; Smith, 1983; Clutter et al., 1983). Essas duas variáveis são limitantes para o desenvolvimento de modelos de produção para as florestas da Amazônia porque são praticamente indisponíveis para o setor florestal, num curto prazo. Apesar de inúmeras tentativas, por meio da dendrocronologia ou da datação com 14C, a determinação das idades das inúmeras espécies que ocorrem numa determinada área, continua sendo um grande obstáculo para a ciência florestal.

Sem a idade da árvore ou do povoamento ou com muita dificuldade para obtê-la, a alternativa é prognosticar a dinâmica da floresta com o uso de parcelas permanentes. Na Amazônia, entretanto, as parcelas instaladas e devidamente monitoradas são poucas, mal distribuídas e recentes (as mais antigas estão na Flona de Tapajós, desde 1978). Considerando que as idades de árvores com DAP > 50 cm, na região de Manaus, podem variar de 200 a 1400 anos, segundo Chambers et al. (1998), 20-30 anos de observações podem parecer insuficientes para descrever, com confiança, a dinâmica de uma floresta da Amazônia.

Apesar de todas as dificuldades, aproximadamente um milhão de hectares de floresta amazônica são manejados anualmente, para produção madeireira sob algum tipo de manejo em regime de rendimento sustentável. A alternativa pode estar no uso de modelos de curto prazo que dependem exclusivamente da situação imediatamente anterior a atual, tendo como objetivo a projeção apenas para uma situação imediatamente posterior. Dentre os vários modelos disponíveis, o que melhor se ajusta às características das florestas da Amazônia, é a cadeia de Markov.

CADEIA DE MARKOV

A cadeia de Markov de primeira ordem é um processo estocástico, no qual as probabilidades de transição durante o intervalo de tempo (t e t+1) dependem apenas do estado do indivíduo no tempo t ou no conhecimento do passado imediato no tempo t+1 e não em qualquer outro estado prévio (Horn, 1975; Chiang, 1980 e Bruner & Moser, 1973). Shugart (1984) enfatiza que a natureza "invariável em tempo" de cada uma das probabilidades de transição é uma importante característica da cadeia de Markov, tendo muita afinidade com o comportamento dos ecossistemas florestais. De acordo com Bierzychudek (1982), um modelo de matriz de transição é um modelo classificado em tamanho ou uma forma da matriz de Leslie. A única exigência deste modelo é divisibilidade da população em grupo de estados e que existam probabilidades de movimento de um estado para outro, com o passar do tempo (Enright & Ogden, 1979).

Shugart & West (1981) apontam que a importância do entendimento dos ecossistemas florestais não é baseada nas idades, mas sim nas mudanças conhecidas no presente. Os modelos determinísticos consistindo de uma simples função matemática (linear, polinomial ou exponencial) não demonstraram ainda que sejam comprovadamente adequados, quando séries de tempo são envolvidas (Morrison, 1976).

Segundo Enright & Ogden (1979), nas florestas tropicais, o atributo tamanho pode ser mais importante do que a idade. Uma razão para isso é que o tamanho pode ser mais ecologicamente informativo do que a idade, quando esta é difícil de ser obtida com precisão. Além disso, segundo ainda os mesmos autores, a divisão de ciclos de vida em estágios de desenvolvimento pode permitir a predição do comportamento futuro mais precisamente do que a divisão em puras classes de idade. Usher (1966) usou o atributo tamanho no lugar da idade para desenvolver um modelo para o manejo de recursos renováveis. Ele afirma que um organismo que está na i-ésima classe no tempo t, pode permanecer na mesma classe, mudar para a classe seguinte (mais de uma classe também) ou morrer, no tempo t+1.

Os modelos que usam matriz de transição são apropriados para análise de muitos problemas biológicos, principalmente em estudos relacionados com a dinâmica da floresta (Enright & Ogden, 1979). Esses modelos têm sido usados intensivamente em estudos de dinâmica de populações de plantas ou animais em várias regiões do mundo. Alguns exemplos são: a demografia do jack-in-the-pulpit em Nova York (Bierzychudek, 1982); dinâmica florestal de uma população de Araucaria numa floresta tropical úmida de Papua Nova Guinea e Nothofagus em floresta montana temperada da Nova Zelândia (Enright & Ogden, 1979); sucessão de térmitas em Gana (Usher, 1979); sucessão florestal na Nova Jersey (Horn, 1975); aplicação da Cadeia de Markov em estudos de dinâmica florestal em florestas tropicais (Acevedo, 1981) e a aplicação de Markov para predizer o desenvolvimento de um povoamento florestal (Usher, 1966; Bruner & Moser, 1973; Peden et al., 1973; Buogiorno & Michie, 1980).

Alder (1980) também descreve a matriz de transição como uma possível ferramenta para análise de dados de crescimento e incremento de povoamentos multianos de florestas tropicais mistas. Na região de Manaus, Higuchi (1987) usou Markov para estudar a dinâmica da floresta das parcelas testemunhas do projeto de manejo florestal (Projeto Bionte) e Rocha (2001) nos transectos do projeto Jacaranda. A maioria dos trabalhos citados anteriormente inclui revisões razoáveis da teoria do método de Markov.

O entendimento do comportamento das taxas de crescimento, recrutamento e mortalidade, em condições naturais ou sob manejo, são importantes para definir quanto tempo uma espécie, ou um grupo de espécies demora em atingir uma determinada dimensão. Essa relação pode definir o ciclo do corte de uma floresta natural, para prescrições adequadas de tratamentos silviculturais e indicar se a floresta está respondendo ou não aos tratos silviculturais bem como, indicar se a floresta funciona como sumidouro ou fonte de dióxido de carbono. Este trabalho tem como objetivo avaliar a dinâmica de uma floresta de terra-firme não explorada na região de Manaus, AM, com o uso da cadeira probabilística de transição.

MATERIAL E MÉTODOS

O estudo foi desenvolvido na Estação Experimental de Silvicultura Tropical do Instituto Nacional de Pesquisas da Amazônia (EEST/INPA), distante aproximadamente 90 km de Manaus - AM, localizada entre os km 21 e 24 da margem esquerda da estrada vicinal ZF-2 do Distrito Agropecuário da Superintendência da Zona Franca de Manaus (SUFRAMA). As coordenadas geográficas da área estão entre 02º37' e 02º38' de latitude sul e 60º 09' e 60º 11' de longitude oeste, aproximadamente (Radam-Brasil, 1978).

De acordo com a classificação de Köppen, o clima é do tipo "Amw", caracterizado como tropical chuvoso, apresentando temperatura elevada, com variação térmica diurna maior do que a variação térmica anual, com precipitação elevada e predominantemente forte. A temperatura média anual é de 26,7ºC, apresentando valores médios para as máximas e mínimas de 31,2 ºC e 23,5 ºC, respectivamente. As chuvas apresentam um total anual da ordem de 2.101mm, distribuídos em duas épocas distintas durante o ano, uma chuvosa que ocorre entre novembro e maio e a outra que se estende até outubro. A umidade relativa apresenta uma média anual de 84 %, variando de 77-88% (Falesi, 1971).

O solo do tipo Latossolo Amarelo Distrófico e Álico de textura argilosa predomina na região (Radam-Brasil, 1978 & Chauvel 1982). O relevo é levemente ondulado e a maioria das ondulações é formada por pequenos platôs, que variam de 500 a 1000 m de altura. A diferença de nível entre as calhas dos igarapés e a superfície dos platôs varia de 70 a 80 m (Santos, 1996). Ferraz et al. (1998) observaram que a topografia dos transectos apresenta-se disposta em seqüências de platôs, encostas e baixios. Os solos dos platôs apresentaram textura argilosa; nas encostas, variaram de argilo-arenosos, próximos aos platôs à areno-argilosos, próximos aos baixios e, os solos nas áreas de baixio, apresentaram textura arenosa. Os solos nas áreas dos transectos podem ser classificados em três tipos: Latossolo Amarelo nos platôs; Podzólicos Vermelho-Amarelo nas encostas e Arenossolos hidromórficos nos baixios.

A vegetação da área é uma amostra representativa de uma floresta tropical úmida densa de terra-firme amazônica. Radam-Brasil (1978) denominou a área de sub-região dos baixos platôs da Amazônia e, com base na geomorfologia, em levantamentos florísticos e em inventários florestais, classificou a mesma em macro-ambiente de relevo tabular, no qual a cobertura florestal densa, raramente com estrato superior uniforme, é freqüentemente alterada por manchas de floresta aberta, onde os estratos arbustivo e herbáceo são compostos por regeneração natural das espécies arbóreas, palmeiras e plantas não-vasculares. A área apresenta vegetação diversificada com 737 espécies, pertencentes a 238 gêneros de 59 famílias. Dentre as famílias mais ricas em espécies estão Sapotaceae (69 espécies) e Chrysobalanaceae e Lauraceae (44 spp.). As espécies que mais se destacaram foram Oenocarpus bataua Mart., Eschweilera wachenheimii (Benoist) Sandwith e Eperua glabiflora (Ducke) R.S. Cowan (Carneiro, 2004).

Os dados foram coletados em duas toposseqüências (transectos) dispostas no sentido Norte-Sul (Figura 1) e Leste-Oeste (Figura 2) do "Projeto Jacaranda", desenvolvido pela CPST-INPA. Cada transecto possui cinco hectares (20 x 2500 m), totalizando dez hectares de floresta natural. O acompanhamento da dinâmica da floresta foi iniciado com a realização do inventário de todos os indivíduos arbóreos, onde se realizou a identificação por nome vulgar e medição do diâmetro à altura do peito (DAP) para todos os indivíduos com DAP ³10 cm.

Para este estudo duas coletas foram necessárias: a primeira foi realizada em 2000, onde se mediu o DAP de cada indivíduo, nos dois transectos. Em 2004, ocorreu a segunda medição de todos os indivíduos já registrados (crescimento em diâmetro), além da inclusão de novos indivíduos que em 2004 passaram a possuir DAP ³10 cm, entrando para a primeira classe de DAP (Recrutamento) e do registro e identificação das árvores que morreram no período entre 2000 e 2004 (Mortalidade).

A hipótese teste neste trabalho é de que a floresta tropical úmida não explorada, na região de Manaus, pode ser avaliada, com previsões futuras, com o uso da cadeia de transição probabilística de Markov. A distribuição diamétrica da floresta foi utilizada como variável principal, sendo aplicado o teste Qui-quadrado, para avaliar a hipótese.

O crescimento da floresta foi analisado por meio dos seguintes parâmetros:

a) Área Basal, medida em 2000 e 2004, para cada transecto, pela equação:

G = (p.DAP²)/4

b) Volume Comercial com Casca, estimado pela equação desenvolvida por Higuchi et al. (1997):

Ln V = - 7, 335 + 2, 121 Ln DAP [R² = 0,95 e S_yx = 0,27]

c) Biomassa Arbórea Fresca Acima do Nível do Solo, estimada pelas equações ajustadas por Higuchi et al. (1998):

Ln P = - 1,754 + 2,665 Ln DAP [R² = 0,92 e S_yx = 43] ( 5 cm £DAP < 20 cm)

Ln P = - 0,151 + 2,17 Ln DAP [R² = 0,90 e S_yx = 2035] (DAP ³20 cm).

d) Mortalidade, obtida pela razão entre o número de indivíduos mortos no período observado, pela abundância de indivíduos arbóreos referente ao ano de 2000:

TM = TM/(NT - TR)

e) Recrutamento: foram consideradas como recrutas todas as árvores que na segunda medição (2004) passaram a apresentar DAP ³10 cm. A taxa de recrutamento foi obtida pela razão entre o número de indivíduos arbóreos ingressantes no período observado e a abundância referente ao número de indivíduos de 2004:

TR = TR/(NT - TM)

G = área basal em m².ha^-1

DAP = diâmetro à altura do peito em cm.

Ln = logaritmo natural

P = peso da biomassa fresca acima do nível do solo, em kg.

TM = Número de indivíduos mortos durante o período observado.

TR = Número de indivíduos recrutados durante o período observado.

NT = Número total de indivíduos no período observado.

VOL = volume em m³

Leste-Oeste = LO

Norte-Sul = NS

RESULTADOS E DISCUSSÃO

ESTOQUES DE VOLUME, BIOMASSA ARBÓREA E ÁREA BASAL NOS TRANSECTOS LO E NS

Na Tabela 1 são apresentados os valores de área basal e volume nos dois transectos (LO e NS) nos anos de 2000 e 2004. A área basal apresentou um valor médio de 51,18 m².ha^-1 ± 3,67 m².ha^-1 (IC 95%) e volume médio de 647,94 m³.ha^-1 ± 47,08 m³.ha^-1 (IC 95 %) no ano 2000. E para o ano de 2004, os valores foram de 54,74 m².ha^-1 ± 0,99 m².ha^-1 (IC 95%) para área basal e de 692,77 m³.ha^-1 ± 12,17 m³.ha^-1 (IC 95 %) para volume. Os dois transectos tiveram incremento em termos de área basal e volume de 0,44 m².ha^-1.ano^-1 e 5,60 m³.ha^-1.ano^-1, respectivamente.

Na Tabela 2 são apresentados os estoques e as diferenças de estoques de biomassa vegetal arbórea acima do nível do solo nos dois transectos no período de 2000 a 2004. A estimativa total de biomassa (LO e NS) apresentou um saldo positivo, ou seja, houve acúmulo de biomassa durante o período estudado de 66,71 t.ha^-1, correspondendo em média a 8,34 t.ha^-1.ano^-1, indicando que a floresta retirou carbono e fixou nas árvores, convertendo-o em tecidos lenhosos, folhas etc. Observando esses valores em média os dois transectos tiveram em 2000 e 2004, 485,85 t ± 71,82 t (IC 95%) e 519,21 t ± 17,30 t (IC 95%), respectivamente.

TAXAS DE MORTALIDADE E RECRUTAMENTO

As causas mais comuns da morte de árvores em florestas tropicais são: idade, estado fitossanitário, estresse hídrico e vento. Em períodos de secas prolongadas (ex. em anos que ocorre o fenômeno El Niño), a mortalidade tende aumentar, porém a maioria das árvores permanece em pé e se decompõem com o passar do tempo. Em períodos de fortes chuvas, as árvores mais frágeis são derrubadas por ação do vento, ocasionando em conseqüência a queda de outras árvores mesmo saudáveis (Rocha et al., 2003).

No período observado (2000-2004), a taxa de recrutamento foi de 3,29% ± 0,36% (IC 95 %) e a taxa de mortalidade de 2,26% ± 0,22% (IC 95%) para os dois transectos. As taxas de recrutamento (1,46% no N-S e 1,83% no L-O) foram superiores às taxas de mortalidade (1,02 % e 1,24% para N-S e L-O, respectivamente). A média anual da taxa de recrutamento foi de 1,65%. Esse valor foi superior às taxas obtidas por Higuchi et al., (1997) e Rocha (2001) (1,1% e 0,90%, respectivamente), na mesma área de estudo, e a taxas obtidas em outros sítios, como na Costa Rica, 1,8 % (Lieberman & Lieberman, 1987), na Malásia, 1,4 % (Manokaran & Kochummen, 1987) e na Flona Tapajós, 1,4 % (Carvalho, 1992), porém, foi inferior ao encontrado por Vasconcelos (2003) de 2,3 % no estado do Acre.

Considerando os dois transectos juntos, a taxa de mortalidade média anual foi de 1,13 %, sendo superior às taxas obtidas por Higuchi et al., (1997) e Rocha (2001) (1,03 % e 0,86 %, respectivamente) na mesma área de estudo. Porém foi inferior quando comparado a Lieberman & Lieberman (1987) com 1,8 % na Costa Rica, Peralta et al., (1987) com 1,8 % a 2,25 % na Venezuela, a Carey et al., (1994) com 0,5 a 3,3 % na Malásia, a Carvalho (1992) com 1,3% na Flona Tapajós e a Vasconcelos (2003), com 2,4 % no Acre. Durante o período observado, 264 árvores morreram. (média de 132 ± 21,56 árvores por transecto). Em termos gerais, a mortalidade foi superior em 17 % quando comparada com a projeção para 2004. A projeção para mortalidade no período 2000-2004 corresponde, aproximadamente, a 83 % do levantamento realizado em 2004.

PROJEÇÃO DA DISTRIBUIÇÃO DIAMÉTRICA USANDO A CADEIA DE MARKOV

Os diâmetros em cada ocasião foram agrupados em classes com intervalos de 5 cm, onde a primeira classe foi de 10 cm £DAP < 15 cm e a última classe diamétrica representada por árvores com DAP ³75 cm.

As matrizes de transição foram montadas com base na distribuição diamétrica das ocasiões (2000 e 2004). A transição de um estado para outro é obtido quando uma árvore que foi medida em 2000 passa uma ou mais classes seguintes ou morrer, ou ainda permanecer na mesma classe de diâmetro ou estado.

A matriz de transição foi estruturada considerando 17 estados, dos quais o primeiro estado trata-se da classe R (Recrutamento); do 2° ao 14° estados correspondem às 14 classes diamétricas, a partir do DAP ³10 cm; e o 15° estado (ultima classe diamétrica) que corresponde aos indivíduos com DAP ³75 cm, podendo variar de 75 cm até o maior diâmetro medido, sendo, por isso, considerada uma classe aberta. O 16° estado refere-se ao número de árvores, pertencentes à classe de DAP ³75 cm, que na segunda ocasião (2004) apresentaram crescimento em diâmetro maior que o intervalo de 5 cm, suficiente para a mudança de uma próxima classe, isto é, se a próxima classe existisse. Por exemplo, um indivíduo que em 2000 apresentou DAP de 99 cm e em 2004 teve o DAP medido em 101 cm, considerando o intervalo de 5 cm, este indivíduo passaria para a próxima classe. E por último, o 17° estado ou classe M (mortalidade), composto pelas árvores que morreram durante o intervalo de tempo, ou seja, pelas árvores que estiveram presentes na primeira ocasião (2000) e ausentes na segunda ocasião (2004).

A Tabela 3 apresenta a transição de um estado para outro, ocorrida durante o período 2000 - 2004 usando as observações de ambos transectos. Verificou-se que o estado R obteve 396 árvores, ou seja, 396 indivíduos ingressaram no banco de dados, pois passaram a possuir o DAP ³10 cm. O estado M apresentou 264 indivíduos, ou seja, 264 indivíduos saíram do sistema por causas naturais (fatores ambientais que causaram a sua morte ou senescência). Em geral, as árvores têm maior probabilidade de permanecer nas mesmas classes diamétrica, principalmente considerando um período curto (quatro anos) e a taxa de crescimento anual observada. Por exemplo, a probabilidade associada à permanência de árvores da classe 10 < 15 cm no período de 2000-2004 foi de 86,2%. A mortalidade no período observado na mesma classe foi 4,2% para o mesmo período. Para essa classe, observa-se que não há um padrão definido, ou seja, as árvores com diâmetros menores têm a mesma probabilidade de morrer que as árvores de diâmetros maiores. De acordo com Rocha (2001), esse comportamento estocástico é típico da floresta tropical úmida, dificultando os trabalhos de modelagem da dinâmica da floresta, principalmente quando se usam modelos determinísticos.

PROJEÇÃO PARA O ANO DE 2008

O resultado final da Cadeia de Markov (Tabela 4) é a projeção da distribuição do número de árvores dos diferentes estados (recrutamento, mortalidade e as classes diamétricas), para um período seguinte de mesmo intervalo de tempo. A projeção para o período seguinte (2004 - 2008) mostrou que a probabilidade dos indivíduos da 1ª classe permanecer com DAP £15 cm é de 74,39%. Porém, 16,24% passarão para a 2ª classe, 1,19 % passarão para a 3ª classe e 0,02 % passarão para a 4ª classe.

A projeção para um período seguinte imediato (2008), com base na transição probabilística do ano de 2004 e a ocasião imediatamente anterior (2000), indica que a taxa de mortalidade será de 250 árvores, bem próxima à observada em 2004 (264 árvores). Para 2008, o número total de árvores será de 5854, sem incluir a taxa de recrutamento, mas se repetida essa taxa (396 recrutas em 2004), o número total de árvores em 2008 será de 6250, 5,6% maior que o número observado em 2004 (5897).

Para avaliar a qualidade das projeções feitas para o ano de 2008 e confirmar o uso da Cadeia de Markov como ferramenta de avaliação para a projeção da dinâmica da floresta tropical e do comportamento futuro deste povoamento florestal, foram comparadas as freqüências projetadas e observadas nos transectos pelo teste Qui-quadrado (c²) ao nível de 1% e 5% de significância. De acordo com a Tabela 4, o teste Qui-quadrado (c²) mostra que não existe diferença significativa entre as freqüências observadas e projetadas na distribuição do número de árvores em cada classe diamétrica para o ano de 2008 (c²cal = 17,94). Este resultado confirma a hipótese apresentada neste estudo, indicando que em florestas que apresentam as mesmas características similares às da área amostrada, a Cadeia de Markov pode ser utilizada como instrumento para projetar a dinâmica de povoamentos de florestas naturais, pois contribui para o planejamento em curto prazo de atividades que envolvam a utilização de recursos florestais.

CONCLUSÕES

O crescimento médio em área basal e volume foram de 0,44 m².ha^-1.ano^-1 e de 5,60 m³.ha^-1.ano^-1, respectivamente, para uma área de dez hectare de floresta natural de terra-firme no período de 2000-2004. A floresta acumulou em biomassa 8,34 t.ha^-1.ano^-1, indicando que em condições naturais, a floresta seqüestra carbono, fixando-o para a manutenção, produção e crescimento de novos tecidos. Os valores médios para as taxas de mortalidade e recrutamento foram de 1,13% e 1,65 %, respectivamente.

O teste Qui-quadrado não mostrou diferença significativa entre as freqüências observadas e projetadas para a distribuição do número de árvores em cada classe diamétrica, indicando a Cadeia de Markov como instrumento para a projeção de freqüência da distribuição diamétrica, do recrutamento e da mortalidade de florestas similares às da região de Manaus, reproduzindo o processo dinâmico das árvores por meio da variável classe diamétrica. A projeção da dinâmica da floresta pode contribuir para o manejo florestal em florestas tropicais na determinação do ciclo de corte ou ainda mostrar as tendências que a floresta apresentará num futuro próximo, para que o técnico ou empreendedor florestal saiba o momento e o volume ideais para manejar a floresta tropical.

BIBLIOGRAFIA CITADA

Acevedo, M.F. 1981. On Horn's Markovian Model of Forest Dynamics with Particular Reference to Tropical Forests. Theoretical Population Biology, 19: 230-250.

Alder, D. 1980. Forest Volume Estimation and Yield Prediction. V.2 - Yield Prediction. FAO Forestry Paper 22/2. 194pp.

Bierzychudek, P. 1982. The Demography of Jack-in-the-pulpit, a Forest Perennial that Changes Sex. Ecol. Monographs, 52(4):333-351.

Bruner, H.D.; Moser Jr, J.W. 1973. A Markov Chain Approach to the Prediction of Diameter Distributions in Uneven-aged Forest Stands. Can.J.For.Res, 3:409-417.

Buogiorno, J.; Michie, B.C. 1980. A Matrix Model of Unever-aged Forest Management. Forest Science, 26:609-625.

Carey, E.V.; Brown, S.; Gillespie, A.J.R.; Lugo, A.E. 1994. Tree Mortality in Mature Lowland Tropical Moist and Tropical Lower Montane Moist Forest of Venezuela. Biotropica, 26(3): 255-264.

Carneiro, V.M.C. 2004. Composição florística e análise estrutural da floresta primária de terra-firme na bacia do rio Cuieiras, Manaus-AM. Dissertação de Mestrado, Instituto Nacional de Pesquisas da Amazônia/ Universidade Federal do Amazonas, Manaus, Amazonas. 67pp.

Carvalho, J.O.P. 1992. Struture and dynamics of a logger over Brazilian Amazonian rainforest. Doctor of Philosophy's Thesis. University of Oxford. Oxford. 146pp.

Chambers, J.Q.; Higuchi, N; Schimel, J.P. 1998. Ancient Trees in Amazonia. Nature, 391:135-136.

Chauvel, A. 1982. Os latossolos amarelos, álicos, argilosos, dentro dos ecossistemas das bacias experimentais do INPA e da região vizinha. Acta Amazonica, 12(3): 47-60.

Chiang, C.L. 1980. An Introduction to Stochastic Processes and their Applications. Robert E. Krieger Publ. Co., Huntington, New York. 380pp.

Clutter, J.L.; Fortson, J.C.; Pienaar, L.V.; Brister, G.H.; Bailey, R.L. 1983. Timber Management: A Quantitative Approach. John Wiley and Sons, Inc. New York. 333pp.

Enright, N.; Ogden, J. 1979. Applications of Transition Matrix Models in Forest Dynamics: Araucaria in Papua New Guinea and Nothofagus in New Zealand. Australian. J. of Ecology, 4:3-23.

Falesi, I.C. 1971. Solos do Distrito Agropecuário da Suframa. Bol. Instituto de Pesquisa e Experimentação Agropecuária da Amazônia Ocidental. Série: Solos. Manaus, Amazonas.1(1). 99pp.

Ferguson, I.S.; Leech, J.W. 1978. Generalized Least Squares Estimation of Yield Functions. Forest Science, 24: 27-42.

Ferraz, J.; Ohta, S.; Salles, P.C. 1998. Distribuição dos Solos ao Longo de Dois Transectos em Floresta Primária ao Norte de Manaus (AM). In: Higuchi, N; Campos, M.A.A.; Sampaio, P.T.B.; Santos, J. (Eds). Pesquisas Florestais para Conservação da Floresta e Reabilitação de Áreas Degradadas da Amazônia.111-143.

Higuchi, N. 1987. Short-term growth of an undisturbed tropical moist forest in the Brasilian Amazon. Doctor of Philosophy's Thesis, Michigan State University, Michigan. 129p.

Higuchi, N. 1994. Utilização e Manejo dos Recursos Madeireiros das Florestas Tropicais Úmidas. Acta Amazonica, 24(3-4): 275-288.

Higuchi, N.; Santos, J.; Ribeiro, R.J.; Freitas, J.V. de; Vieira, G.; Cöic, A.; Minette, L.J. 1997. Crescimento e Incremento de uma Floresta Amazônica de Terra-Firme Manejada Experimentalmente. In: Relatório Final do Projeto Bionte. Biomassa e Nutrientes Florestais - MCT/INPA. Manaus. 89-131.

Higuchi, N.; Santos, J.; Ribeiro, R.J.; Minette, L.J.; Biot, Y. 1998. Biomassa da parte aérea da vegetação da Floresta Tropical Úmida de Terra-Firme da Amazônia Brasileira. Acta Amazonica, 28 (2): 153-166.

Horn, H.S. 1975. Markovian Properties of Forest Succession. In: M. Cody e J. Diamond, Eds. Ecology and Evolution of Communities, Harvard University Press. Cambridge, Mass. 196-211.

Lieberman, D.; Lieberman, M. 1987. Forest tree growth and dynamics at La Selva, Costa Rica (1969 - 1982). J. of Trop. Eco, 3:347-358.

Manokaran, N.; Kochummen, K.M. 1987. Recruitment, grow and mortality of trees in a lowland dipyterocarp forest in Peninsular Malaysia. J. of Trop. Eco, 3:315-330.

Morrison, D.F. 1976. Multivariate Statistical Methods. McGraw-Hill Inc.. 415pp.

Peralta, R.; Hartshorn, G.S.; Lieberman, D.; Lieberman, M. 1987. Reseña de estudios a largo plazo sobre composición florística y dinámica del bosque tropical en La Selva, Costa Rica. In: Clark, D.A.; Dirzo, R.; Fetcher, N. (Eds). Ecologia y ecofisiología de plantas en los bosques mesoamericanos. Revista de Biologia Tropical, 35 (1). 23-40.

Peden, L.M., Williams, J.S.; Frayer, W.E. 1973. A Markov Model for Stand Projection. Forest Science, 19: 303-314.

Radam-Brasil. 1978. Programa de Integração Nacional. Levantamento de Recursos Naturais. V. 18 (Manaus) – RADAM (projeto) DNPM, Ministério das Minas e Energia. Brasil. 626pp.

Rocha, R.M. 2001. Taxas de Recrutamento e Mortalidade da Floresta de Terra Firme da Bacia do Rio Cuieiras na Região de Manaus-AM. Dissertação de Mestrado. Instituto Nacional de Pesquisas da Amazônia/Universidade Federal do Amazonas, Manaus, Amazonas. 49pp.

Rocha, R.M.; Higuchi, N.; Santos, J.; Nakamura, S.; Silva, R.P.; Pinto, A.C.M.; Tribuzy, E.S. 2003. Taxas de recrutamento e mortalidade e mudanças de estoques de fitomassa da floresta primária na região de Manaus-Am. In: Higuchi, N; Santos, J. dos; Sampaio, P. B. T.; Marenco, R. A.; Ferraz, J.; Sales, P.C.; Saito, M.; Matsumoto, S. (Eds). Pesquisas Florestais na Amazônia Central. Projeto Jacaranda Fase II. 43-54.

Santos, J. 1996. Análise de modelos de regressão para estimar a fitomassa da floresta tropical úmida de terra-firme da Amazônia brasileira. Tese de Doutorado, Universidade Federal de Viçosa. Viçosa, Minas Gerais. 120pp.

Shugart, H.H.; West, D.C. 1981. Long-term Dynamics of Forest Ecosystems. Am. Scientist., 69:647-652.

Shugart, H.H. 1984. A Theory of Forest Dynamics: The Ecological Forest Succession Models. Springer-Verlag Inc. New York. 278p.

Smith, V.G. 1983. Compatible Basal Area Growth and Yield Models Consistent with Forest Growth Theory. Forest Science, 29:279-288.

Sullivan, A.D.; Clutter, J.L. 1972. A Simultaneous Growth and Yield Model for Loblolly Pine. Forest Science, 18:76-86.

Usher, M.B. 1966. A Matrix Approach to the Management of Renewable Resources, with Special Reference to Selection Forests. J. of Applied Ecology. 3(2): 355-367.

Usher, M.B. 1979. Markovian Approaches to Ecological Succession. J. of Animal Ecology, 48: 413-426.

Vasconcelos, S.S. 2003. Dinâmica de uma floresta explorada seletivamente no Projeto de Colonização Pedro Peixoto na Amazônia Ocidental. Dissertação de Mestrado. Instituto Nacional de Pesquisas da Amazônia/Universidade Federal do Amazonas, Manaus, Amazonas. 71pp.

Recebido em 09/08/2006

Aceito em 08/02/2007

Datas de Publicação

• Publicação nesta coleção
  31 Jan 2008
• Data do Fascículo
  2007

Histórico

• Aceito
  08 Fev 2007
• Recebido
  09 Ago 2006