http://dx.doi.org/10.1590/S0102-33062004000400005
Flora e estrutura do componente arbóreo da floresta do Parque Ecológico Quedas do Rio Bonito, Lavras, MG, e correlações entre a distribuição das espécies e variáveis ambientais1
Paulo Eduardo DalanesiI, 2; Ary Teixeira de Oliveira-FilhoII; Marco Aurélio Leite FontesII
IFaculdade de Engenharia Florestal, Fundação Integrada Municipal de Ensino Superior, CEP 75830-000, Mineiros, GO, Brasil
IIDepartamento de Ciências Florestais, Universidade Federal de Lavras, CEP 37200-000, Lavras, MG, Brasil
Palavras-chave: candeal, diversidade, comunidade arbórea, floresta de altitude, solos, topografia
Key words: diversity, elfin forest, tropical montane forest, soils, topography, tree community
A vegetação primitiva da região do Alto Rio Grande, ao sul de Minas Gerais (entre 21º00'- 22º20'S e 43º50'- 45º00'W), compreendia um complexo mosaico composto de manchas de floresta, cerrado, campo de altitude e campo rupestre (Eiten 1982). Esta notável variação fisionômica da vegetação deve-se a dois motivos principais. Em primeiro lugar, a região abriga uma das áreas de transição entre os domínios do cerrado e da floresta Atlântica sensu lato. Em segundo lugar, esta transição verifica-se em pleno domínio da serra da Mantiqueira, cujo relevo acidentado promove grande variação ambiental e vegetacional, com fisionomias peculiares nas maiores altitudes. De maneira geral, as fisionomias de campo rupestre e campo de altitude estão associadas aos solos mais rasos, enquanto, nos mais profundos ocorrem cerrados ou florestas, discriminados pela fertilidade e regime de água dos solos e pela freqüência de incêndios (Oliveira-Filho et al. 1994b).
O Parque Ecológico Quedas do Rio Bonito (PEQRB), localizado ao sul do município de Lavras, constitui uma valiosa amostra da vegetação primitiva da região do Alto Rio Grande, pois seus quatro tipos fisionômicos principais - floresta, cerrado, campo rupestre e campo de altitude - encontram-se bem representados e razoavelmente preservados. De acordo com Oliveira-Filho & Fluminhan-Filho (2000), os campos de altitude e campos rupestres ocorrem principalmente nas áreas mais elevadas e com solos mais rasos; o cerrado, nas áreas mais baixas e com solos mais profundos; e as florestas concentram-se no fundo dos vales. Na transição para os campos, a fisionomia das florestas diferencia-se em uma faixa ecotonal denominada candeal, devido à predominância da candeia (Eremanthus erythropappus (DC.) MacLeish).
O Parque abriga notável diversidade de espécies de plantas em uma área relativamente pequena, o que atraiu a atenção de pesquisadores da Universidade Federal de Lavras e de outras instituições, que vêm realizando ali vários levantamentos florísticos (Oliveira-Filho & Fluminhan-Filho 2000). A floresta do PEQRB distribui-se nas encostas adjacentes ao córrego dos Vilas Boas e possui área total de 74,4ha, o que a caracteriza como o maior fragmento florestal contínuo do município de Lavras (L.T. Oliveira, dados não publicados). Sua vegetação pode ser classificada como floresta estacional semidecidual montana, de acordo com o sistema do IBGE (Veloso et al. 1991), ou como floresta semidecídua alto-montana do domínio da Floresta Atlântica, de acordo com Oliveira-Filho & Fontes (2000). Estudos quantitativos e ecológicos da vegetação desta floresta iniciaram-se com Oliveira-Filho et al. (1994a), em um segmento contínuo, e Rondon-Neto et al. (2000), em uma área em regeneração.
O presente trabalho teve como objetivo básico expandir os estudos quantitativos e ecológicos da comunidade arbórea da floresta do Parque Ecológico Quedas do Rio Bonito. Ao trecho amostrado por Oliveira-Filho et al. (1994a) foram acrescentados mais dois, aumentando a área amostrada de 0,48ha para 2,40ha. As coletas de material botânico dentro e fora da área amostral também incrementaram substancialmente o levantamento florístico da área. Os objetivos específicos do trabalho foram (a) descrever a composição da flora arbórea da floresta do Parque Ecológico Quedas do Rio Bonito; (b) descrever a estrutura comunitária e estrutura fisionômica da floresta em diferentes hábitats de solo identificados em três trechos da floresta; e (c) avaliar as correlações entre a distribuição das abundâncias das espécies e variáveis ambientais.
Área de estudo - O Parque Ecológico Quedas do Rio Bonito (PEQRB) está localizado no extremo sul do município de Lavras, Minas Gerais, junto à divisa com o município de Ingaí e entre as coordenadas geográficas de 21º19'45" - 21º20'48"S e 44º58'18" - 44º59'24"WGW, com altitudes variando de 1.000 a 1.300m (Fig. 1). O Parque conta com uma área de 209,7ha e é de propriedade da Fundação Abraham Kasinski (FAK), à qual foi doado pela municipalidade em 22/julho/1994. O PEQRB está localizado a aproximadamente 13km do centro de Lavras pela estrada Lavras—Luminárias, que o delimita ao sul. Sua área compreende a microbacia do córrego dos Vilas Boas, à montante da antiga barragem de captação de água, apoiando-se em um dos contrafortes da serra do Carrapato, que, juntamente com as vizinhas serras da Bocaina, do Faria, do Campestre, da Estância e do Pombeiro, constituem, ao mesmo tempo, uma disjunção da serra do Espinhaço e sua ligação com o complexo da Mantiqueira. As principais rochas no PEQRB são os quartzitos e micaxistos, nas áreas mais altas (superfície de erosão Gondwana), e os gnaisses graníticos leucocráticos e quartzitos, nas áreas mais baixas (superfície de erosão pós-Gondwana), e os solos predominantes são Neossolos Litólicos e Cambissolos Háplicos (Curi et al. 1990).
De acordo com a classificação de Köppen, o padrão climático da região do PEQRB enquadra-se no tipo Cwb, ou seja, temperado mesotérmico com verões brandos e suaves e estiagens de inverno. A temperatura média anual situa-se em torno de 19,3ºC, a do mês mais frio (junho) oscila entre 13 e 16ºC e a do mês mais quente (janeiro), entre 21 e 23ºC. A precipitação média anual é de 1.493mm, com chuvas mal distribuídas durante o ano e concentradas no verão, 67% do total precipitando entre novembro e fevereiro. O inverno tem cerca de quatro meses, com pequeno déficit hídrico, entre 10 e 30mm, de acordo com o balanço hídrico de Thorntwaite e Mather (Golfari 1975).
Inventário da comunidade arbórea - A amostragem foi realizada em três transecções demarcadas no interior da floresta do PEQRB, com áreas de 0,48, 0,78 e 1,14ha e separadas por uma distância de 480m, totalizando área amostral de 2,40ha (Fig. 1). As transecções, denominadas A, B e C, foram compostas por, respectivamente, 16, 26 e 38 parcelas contíguas de 10×30m, totalizando 80 parcelas (Fig. 2). As transecções foram dispostas de maneira a cortar a calha do córrego dos Vilas Boas e estender-se nas duas encostas adjacentes, resultando em um alinhamento perpendicular ao sentido predominante do curso d'água. A forma e disposição das parcelas foram escolhidas de forma a obedecer ao princípio expresso por Causton (1988), segundo o qual formas retangulares, dispostas de maneira que o maior comprimento seja perpendicular à direção do gradiente principal presumido, são mais adequadas para estimar os parâmetros da vegetação. Desta forma, o lado menor das parcelas retangulares foi disposto paralelamente ao alinhamento da transecção com o objetivo de investigar a influência do gradiente topográfico na vegetação e nas características dos solos. As transecções A e B, que são mais curtas, empregaram dois alinhamentos contíguos de parcelas e a transecção C, mais longa, apenas um alinhamento.
A transecção A corresponde ao primeiro inventário realizado na floresta do PEQRB, em 1989, cujos resultados foram publicados por Oliveira-Filho et al. (1994a). Os inventários das transecções B e C foram realizados em 2000 e 2001, respectivamente. Nos inventários de A e B foram feitas medições de circunferência a altura do peito (CAP) e do solo (CAS) em cada árvore com CAS > 15,7cm. No inventário de C, foram medidas somente as CAPs de árvores com CAP > 15,7cm (equivalente a um diâmetro de 5,0cm). O presente trabalho reúne as informações dos três inventários e as padroniza pelos critérios da transecção C. Indivíduos com caules múltiplos foram incluídos quando a raiz quadrada da soma dos quadrados das CAPs era superior ao limite estabelecido. Foram registrados, para cada indivíduo, o nome da espécie, o valor da CAP (medida com fita métrica) e a altura (estimada com o auxílio de varas graduadas). Cada árvore recebeu uma plaqueta de alumínio numerada, identificando a parcela e o indivíduo. Alguns indivíduos de espécies abundantes e bem conhecidas foram identificados no campo, enquanto os demais tiveram material botânico coletado para posterior identificação. Foram realizadas também coletas adicionais em caminhadas de varredura por todo o fragmento florestal, com o propósito de complementar o levantamento florístico. A herborização do material botânico foi feita no Herbário do Departamento de Biologia da Universidade Federal de Lavras (Herbário ESAL) e a coleção de referência foi incorporada ao seu acervo. A identificação do material botânico foi realizada com a utilização de coleções botânicas já existentes no Herbário ESAL e também por meio de consultas à literatura, a especialistas e a coleções dos Herbários do Jardim Botânico do Rio de Janeiro (RB), Instituto de Botânica de São Paulo (SP) e Universidade Estadual de Campinas (UEC). As espécies de angiospermas (todas exceto Cyatheaceae, Podocarpaceae e Araucariaceae) foram classificadas nas famílias reconhecidas pelo Angiosperm Phylogeny Group II (APG II 2003).
Levantamento das variáveis ambientais - Foi realizado levantamento topográfico das três transecções amostrais. A inclinação do terreno foi medida com o auxílio de hipsômetro de Suunto, sendo anotados os desníveis entre os vértices das parcelas e calculadas as distâncias verticais. Com essas medidas, foram traçadas curvas de nível em intervalos de 5m a partir da margem do córrego, onde a cota era conhecida. Por meio das cotas foi confeccionada uma grade de superfície, formando uma representação tridimensional das três transecções (Fig. 2), com o auxílio do programa SURFER versão 5.0, da Golden Software. A cota dos quatro vértices de cada parcela permitiu a obtenção das variáveis topográficas cota média, desnível (diferença entre as cotas máxima e mínima) e distância vertical (diferença entre cota média e cota do córrego) (Botrel et al. 2002).
Os solos de cada parcela foram classificados de acordo com o novo Sistema Brasileiro de Classificação de Solos (EMBRAPA 1999) chegando até ao nível de subgrupo, ou 4° nível categórico. As classes de solo foram utilizadas em análises comparativas da comunidade arbórea, em subamostras consideradas como 'hábitats de solo'. A variável ordinal drenagem foi obtida para cada parcela por meio da atribuição de números de ordem (ranking) às classes de drenagem definidas pela EMBRAPA (1999) em cada classe de solo amostrada. Para obter as variáveis químicas e texturais do solo, coletaram-se amostras simples de cerca de 1L de solo superficial (0-20cm de profundidade) no centro de cada parcela. As amostras foram enviadas para o Laboratório de Análise de Solos do Departamento de Ciência do Solo da Universidade Federal de Lavras para análises químicas e texturais, de acordo com os métodos descritos pela EMBRAPA (1997). As variáveis do solo obtidas foram: pH; teores de P, K, Ca, Mg e Al; saturação por bases (V); matéria orgânica e teores de areia, silte e argila.
Para avaliar o efeito borda na distribuição das espécies foi produzido um 'fator borda' para cada parcela, de acordo com o método concebido por Oliveira-Filho et al. (1997). O fator borda consiste do comprimento linear do contorno da borda medido entre dois (ou mais) pontos de intercepto com um círculo cujo raio tem origem no centro da parcela. O comprimento deste raio representaria o raio de ação máxima do efeito borda sobre a comunidade arbórea, tendo sido utilizado o mesmo valor de 100m adotado por Espírito-Santo et al. (2002) e Souza et al. (2003) em fragmentos florestais da mesma região. No presente caso, a borda da floresta foi definida no limite entre o candeal e o campo. Ainda que possa ser visto como ecótono entre floresta e campo, o candeal apresenta, na área de estudos, uma fisionomia claramente florestal e apresenta limites abruptos com o campo e transições mais graduais e imprecisas com o interior da floresta.
Descrição da estrutura da comunidade e estrutura fisionômica - Para cada espécie amostrada nas parcelas foram calculados os seguintes parâmetros fitossociológicos: densidade absoluta, densidade relativa, freqüência absoluta, freqüência relativa, dominância absoluta, dominância relativa e o valor de importância (VI), obtido por meio da soma dos três parâmetros relativos citados (Mueller-Dombois & Ellenberg 1974). Estes parâmetros foram calculados tanto para a amostra total como para subamostras definidas por hábitats de solo.
Para expressar a diversidade de espécies arbóreas na amostra total e nas subamostras definidas pelos hábitats de solo, foram calculados o índice de diversidade de Shannon (H'), o índice de eqüabilidade de Pielou (J') (Krebs 1989), a densidade média de espécies por parcelas e os estimadores 'jackknife' de primeira e segunda ordem para o número total de espécies (Heltsche & Forrester 1983; Palmer 1991). Os valores de H' obtidos para as subamostras de cada hábitat de solo foram comparados aos pares pelo teste de t de Hutcheson (Zar 1996). A densidade de espécie por parcela foi comparada entre os hábitats de solo com análise de variância seguida de um teste de Tukey-Kramer.
A densidade e a área basal por hectare, bem como a altura e o diâmetro médios, foram comparados entre os hábitats de solo com análises de variância seguidas de testes de Tukey-Kramer. Foram preparadas distribuições de densidade de árvores por classes de diâmetro e altura para os hábitats de solo. No caso dos diâmetros, foram empregados intervalos de classe com amplitudes crescentes para compensar o forte decréscimo da densidade nas classes de tamanho maiores, típico da distribuição em J-invertido. Estes intervalos permitem uma melhor representação das classes diamétricas maiores e de baixa densidade, o que é desejável em comparações gráficas (Oliveira-Filho et al. 2001; Botrel et al. 2002).
Correlações entre distribuição de espécies e variáveis ambientais - Para realizar análise global das correlações entre a distribuição das abundâncias das espécies e variáveis ambientais na área amostral, foi utilizada análise correspondência canônica (CCA) (ter Braak 1987) com o programa PC-ORD for Windows versão 4.14 (McCune & Mefford 1999). A CCA requer duas matrizes de trabalho: a matriz das espécies e a matriz de variáveis ambientais. A matriz de abundâncias das espécies foi constituída do número de indivíduos por parcela das 50 espécies que apresentaram 30 ou mais indivíduos na amostra total. De acordo com as recomendações de ter Braak (1995), os valores de abundância foram transformados pela expressão log10(a + 1) para compensar os desvios causados por alguns poucos valores muito elevados. A matriz das variáveis ambientais incluiu inicialmente todas as variáveis químicas e texturais do solo, as três variáveis topográficas, as classes de drenagem e o fator borda. Após processar uma CCA preliminar, foram eliminadas as variáveis redundantes ou que apresentaram baixa correlação com os dois primeiros eixos de ordenação. As sete variáveis utilizadas pela CCA final foram: fator borda, cota média, classe de drenagem, pH, Al, V e matéria orgânica. Foi aplicado o teste de permutação de Monte Carlo para verificar a significância das correlações entre os padrões emergentes das espécies e variáveis ambientais na CCA final. As classes de solo foram discriminadas no diagrama de ordenação das parcelas para facilitar a interpretação ambiental a posteriori dos padrões emergentes das espécies.
Composição florística - A lista de espécies arbóreas encontradas na floresta do Parque Ecológico Quedas do Rio Bonito é apresentada na Tabela 1. Todos os espécimes coletados foram identificados até o nível de espécie. No total, foram registradas 384 espécies, pertencentes a 182 gêneros e 68 famílias botânicas. Dentre elas, 171 espécies foram acrescentadas na listagem florística a partir das coletas realizadas fora das transecções. As famílias que contribuíram com maior riqueza de espécies foram Myrtaceae, com 52 espécies, Fabaceae, com 38, Melastomataceae, com 30, Lauraceae, com 26, e Asteraceae, com 15. Essas cinco famílias somaram 41,4% do total das espécies. Por outro lado, 20 famílias (28,57%) foram amostradas por uma única espécie. É justo salientar que a presente listagem resulta do trabalho desenvolvido no PEQRB ao longo de mais de 20 anos por vários pesquisadores.
Variáveis topográficas e edáficas - Foram identificadas, nas três transecções amostrais da floresta do PEQRB, quatro classes de solo (Fig. 2). A primeira classe foi constituída pelos Neossolos Litólicos Distróficos típicos, que compreendem solos com horizonte A ou O hístico com menos de 40cm de espessura e assente diretamente sobre rochas ou sobre um horizonte C. Estes solos foram amostrados em 20 parcelas, distribuídas nas três transecções. A segunda classe compreendeu os Neossolos Regolíticos Distróficos típicos, que são solos incipientes com horizonte A sobrejacente ao horizonte C, às vezes com horizonte Bi com menos de 10cm de espessura e apresentando contato lítico a uma profundidade maior que 50cm. Estes solos foram amostrados em 10 parcelas da transecção C. A terceira e quarta classes identificadas foram a dos Cambissolos Háplicos Tb Distróficos lépticos, amostrados em 24 parcelas das transecções A e B, e os Cambissolos Háplicos Tb Distróficos típicos, amostrados em 26 parcelas da transecção C. Os Cambissolos, de maneira geral, são solos constituídos por material mineral com horizonte B incipiente (Bi) subjacente a qualquer tipo de horizonte superficial. As quatro classes de solos serão chamadas, daqui em diante, de Neossolos Litólicos, Neossolos Regolíticos, Cambissolos Lépticos e Cambissolos Típicos.
As médias e desvios padrões das variáveis topográficas e edáficas nas subamostras correspondentes às quatro classes de solo encontram-se na Tabela 2, juntamente com os resultados de comparações estatísticas. Para as variáveis topográficas, não houve diferença significativa entre os solos para a cota média, mas o desnível do terreno foi mais acentuado nos Neossolos Regolíticos, decrescendo no sentido dos Cambissolos Típicos, Neossolos Litólicos e Cambissolos Lépticos. A distância vertical até o córrego foi maior nos Cambissolos Típicos e menor nos Cambissolos Lépticos. As quatro classes de solos apresentam textura franco-arenosa e têm acentuado caráter distrófico. Não houve diferença significativa entre as quatro classes de solo para proporções de areia, silte e argila. Contudo, os Neossolos Regolíticos destacamse pelos maiores teores de P, K, Ca e Mg, resultando em maior saturação por bases. Os Cambissolos Lépticos destacam-se pelo pH mais baixo e teor de Al mais alto e os Cambissolos Típicos, pelo teor mais baixo de matéria orgânica.
Estrutura e diversidade da comunidade arbórea - Na área total de 2,4ha foram amostrados 4.137 indivíduos, os quais somaram área basal de 20,5m2/ha. Em valores absolutos, a densidade foi maior nos Cambissolos Lépticos e a área basal, nos Cambissolos Típicos, porém não houve diferenças significativas tanto para densidade quanto para área basal entre as quatro classes de solos amostradas (Tab. 3). Por outro lado, houve contrastes significativos entre os solos nas comparações das médias de diâmetro e altura por parcela. O DAP médio foi significativamente mais alto nos Neossolos Regolíticos que nos Neossolos Litólicos e Cambissolos Lépticos; e não diferiu significativamente entre os Cambissolos Típicos. A altura média foi significativamente mais alta nos Neossolos Regolíticos que nos Cambissolos Típicos; as duas outras classes não diferiram significativamente. Estas diferenças ficam mais evidentes quando se comparam as distribuições da densidade de árvores por classes de diâmetro e altura (Fig. 3). A densidade de árvores na classe de 5 a 9cm diâm. foi maior no Cambissolo Léptico e menor para o Neossolo Regolítico. Para as demais classes não houve diferenças marcantes. Para as classes de altura, destacam-se os dois Cambissolos com maior densidade na classe de 5 a 10m. Nas demais classes, as diferenças são insignificantes.
Os 4.137 espécimes amostrados nas parcelas distribuíram-se em 213 espécies e 58 famílias botânicas. Este número de espécies representa 75 a 82% da riqueza total esperada pelos estimadores jackknife de primeira e segunda ordem, respectivamente (Tab. 3). No levantamento florístico, que incluiu coletas feitas fora das parcelas, o número de espécies foi de 384, superando em mais de 100 espécies o número esperado pelo estimador jackknife de segunda ordem. O índice de diversidade de Shannon (H') foi de 4,56 nats.ind.-1 e o índice de eqüabilidade de Pielou (J'), de 0,85 para a amostra total. O valor de H' foi significativamente maior nos Cambissolos Típicos que nas demais classes de solo. Como a eqüabilidade de Pielou teve valores próximos em todas as classes de solo, o maior índice de Shannon dos Cambissolos Típicos deve-se certamente à maior riqueza de espécies. A densidade de espécies por parcela também foi significativamente mais elevada nos Cambissolos Típicos do que nos Neossolos Litólicos, mas não houve diferença significativa entre estes extremos e as duas outras classes de solo.
A Tabela 1 apresenta os parâmetros fitossociológicos das espécies amostradas nas 80 parcelas. As dez espécies que apresentaram os maiores valores de VI (valor de importância) foram: Tapirira obtusa, Sclerolobium rugosum, Siphoneugena densiflora, Copaifera langsdorffii, Pera glabrata, Miconia trianae, Alchornea triplinervia, Vochysia magnifica, Croton verrucosus e Amaioua guianensis. Juntas, estas espécies representaram 26,91% do VI total. As 10 espécies que apresentaram maior densidade foram: T. obtusa, S. densiflora, C. langsdorffii, Amaioua guianensis, M. trianae, V. magnifica, Protium widgrenii, P. glabrata, Gymnanthes concolor e Psychotria vellosiana. Estas espécies, juntas, representaram 29,82% dos indivíduos registrados. As cinco espécies que apresentaram maior área basal foram S. rugosum, T. obtusa, C. langsdorffii, P. glabrata e S. densiflora que, juntas, representaram 23,06% da área basal total.
Na Tabela 4 são apresentados os valores dos parâmetros fitossociológicos das 10 espécies de maior VI para as quatro classes de solo amostradas. Algumas espécies são compartilhadas nestas posições entre as classes de solo. T. obtusa e S. densiflora são as únicas compartilhadas pelos quatro solos. S. rugosum e C. langsdorffii ocorrem nos Neossolos Litólicos e Cambissolos Lépticos, M. trianae e C. verrucosus ocorrem nos Neossolos Litólicos e Cambissolos Típicos, Coccoloba warmingii ocorre nos Neossolos Regolíticos e Cambissolos Lépticos, A. triplinervia ocorre nos Cambissolos Lépticos e Típicos e Croton floribundus e V. magnifica ocorrem nos Cambissolos Típicos e Neossolos Regolíticos. As outras 16 espécies ocorrem em só uma classe de solo.
Padrões de distribuição das espécies - Os autovalores da análise de correspondência canônica (CCA) foram baixos, 0,362 (eixo 1) e 0,224 (eixo 2), confirmando a existência de gradientes 'curtos' sensu ter Braak (1995), nos quais há mais variação nas abundâncias relativas das espécies do que na substituição de espécies. Os dois primeiros eixos da CCA explicaram em conjunto apenas 16,5% (eixo 1, 10,2% e eixo 2, 6,3%) da variância global dos dados, indicando que grande proporção da variância permaneceu sem explicação. Os baixos valores observados de variância percentual para abundância de espécies são comuns em dados de vegetação e não prejudicam a significância das relações espécie-ambiente (ter Braak 1988). De fato, a CCA produziu valores muito altos para as correlações espécie-ambiente nos eixos 1 e 2: 0,923 e 0,821, respectivamente. Além disso, os testes de permutação de Monte Carlo indicaram, para os dois eixos de ordenação, que os gradientes das abundâncias das espécies diferem significativamente de padrões aleatórios (P < 0,01 para autovalores) e que as abundâncias das espécies foram significativamente correlacionadas com as variáveis ambientais utilizadas (P < 0,01 para correlações espécies-variáveis). As correlações das variáveis ambientais com o primeiro eixo de ordenação foram, em ordem decrescente de valor absoluto, para fator borda (0,927), Al (0,675), cota (-0,636), pH (-0,560), saturação por bases (-0,549) e matéria orgânica (0,530) e classes de drenagem (0,312). A variável classes de drenagem apresentou correlação mais forte com o segundo eixo de ordenação (0,850), as demais foram todas inferiores a 0,500.
No diagrama de ordenação das parcelas (Fig. 4A), pode-se observar que as quatro classes de solo foram discriminadas em setores diferentes do diagrama, o que reforça a visão das mesmas como hábitats bem definidos e com composição de espécies particular porque as classes de solos em si não interferem no resultado da CCA. O primeiro eixo da CCA discriminou, à esquerda, parcelas de Neossolos Regolíticos e Cambissolos Típicos como apresentando cotas mais elevadas (transecção C), pH e V mais altos, teores de alumínio mais baixos e efeito borda menos pronunciado. As tendências inversas se apresentam no lado direito do diagrama, onde se concentram parcelas de Cambissolos Lépticos (Transecções A e B). O segundo eixo da CCA discriminou, na parte superior, as parcelas dos Neossolos Litólicos, associadas a uma drenagem mais forte. O diagrama de ordenação das espécies (Fig. 4B) sugere que Myrcia rostrata, Daphnopsis fasciculata, Psychotria vellosiana, Clethra scabra e Eremanthus erythropappus tendem a ser mais abundantes nas áreas com drenagem mais acentuada, distantes da margem do córrego, próximas à borda da floresta e com predominância do Neossolos Litólicos. Por outro lado, espécies como Faramea cyanea, Ixora warmingii, Salacia elliptica e Copaifera langsdorffii são mais abundantes junto à borda da floresta ou nos solos com maior teor de alumínio e menor V. Outras espécies, como Pimenta pseudocaryophyllus, Myrcia fallax e Prunus myrtifolia, tendem a ser mais abundantes quanto mais distantes estiverem da borda e/ou em solos de maior V e menor teor de alumínio.
Composição florística e diversidade de espécies Devido principalmente à sua grande heterogeneidade ambiental, as matas ripárias apresentam, em geral, riqueza de espécies elevada em relação às florestas vizinhas não associadas a cursos d'água (Oliveira-Filho et al. 1990; Durigan et al. 2000). A alta riqueza de espécies encontrada na floresta do PEQRB deve-se, provavelmente, à combinação de vários fatores. Salienta-se a forte heterogeneidade ambiental incrementada pela condição ripária de parte da floresta, pelo pronunciado gradiente edáfico e topográfico e pelo efeito borda na sua transição para as formações campestres. Deve-se considerar que o PEQRB representa uma área de transição entre florestas, cerrados e campos onde o contato entre formações fisionômicas distintas pode contribuir para o enriquecimento da flora, conforme já observado por Botrel et al. (2002) no vizinho município de Ingaí. A grande heterogeneidade da vegetação e ambiente certamente não foi bem representada pelas transecções amostrais, o que foi confirmado pelos estimadores 'jackknife' projetados a partir das mesmas, pois estes foram largamente superados pelo número de espécies registrado no levantamento florístico total.
Apesar das considerações acima, a alta riqueza de espécies da floresta do PEQRB não pode ser atribuída primariamente à intensa amostragem, pois esta registrou o maior valor de riqueza de espécies entre 20 áreas de floresta da região do Alto Rio Grande de cujos inventários Pereira (2003) extraiu subamostras médias de igual tamanho. De acordo com este autor, a alta heterogeneidade ambiental do PEQRB é a razão principal de sua alta riqueza de espécies em relação às outras 19 florestas. A esta heterogeneidade ambiental soma-se a transição com as formações abertas, que incrementa a riqueza de espécies. Para se ter uma idéia da riqueza de espécies na floresta do PEQRB, podese também comparar seu total de 384 espécies arbóreas com as 730 encontradas nas mesmas 20 florestas da região, o que significa representação de 52,6% da lista total da região contida em área relativamente pequena. O índice de diversidade de Shannon para as florestas do PEQRB (H' = 4,56 nats.ind.-1) também é o maior valor encontrado na região do Alto Rio Grande e significativamente diferente dos demais 19 valores, de acordo com Pereira (2003). O índice de equabilidade de Pielou (J' = 0,85) não se destaca da mesma maneira, mas está situado na faixa dos dez maiores valores registrados para região. Como o índice de Shannon é influenciado tanto pela riqueza de espécies quanto pela equabilidade, o alto H' deve-se mais ao componente riqueza. Da mesma forma, a heterogeneidade ambiental relativamente mais alta é provavelmente a principal responsável pelo alto índice.
Observando os valores do índice de Shannon para as quatro classes de solo, verificam-se valores altos para todas elas (H' > 4,0 nats.ind.-1), mas os Cambissolos Típicos destacam-se pelo valor particularmente maior. Este fato provavelmente vincula-se à amostragem deste hábitat de solo, que se estendeu por faixa relativamente mais extensa que os demais. Esta faixa compreendeu desde ambiente mais úmido, próximo ao córrego (ca. 1.090m de altitude), até mais seco, na alta encosta (ca. 1.200m de altitude), o que não ocorreu para as outras classes de solo. Este fato provavelmente ocasiona heterogeneidade ambiental mais pronunciada nos Cambissolos Típicos. A densidade de espécies por parcela, significativamente mais baixa nos Neossolos Litólicos, provavelmente se explica pela predominância neste hábitat de solo da fisionomia do candeal, onde a dominância da candeia (Eremanthus erythropappus) contribui para reduzir a riqueza local de espécies.
Análise do perfil da flora arbórea da floresta do PEQRB nos níveis de famílias, gêneros e espécies mostra claras características das florestas alto-montanas do Domínio da Floresta Atlântica no Sudeste do Brasil, conforme os critérios de Oliveira-Filho & Fontes (2000). Destaca-se o grande número de espécies das famílias Myrtaceae, Melastomataceae, Lauraceae, Asteraceae, Cyatheaceae e Solanaceae e dos gêneros Miconia, Ilex, Cyathea, Clusia e Roupala. Entre as espécies típicas relacionam-se Araucaria angustifolia, Podocarpus sellowii, Drymis brasiliensis, Jacaranda subalpina, Eremanthus erythropappus, Eupatorium velutinum, Clethra scabra, Ilex sapotifolia, Euplassa organensis, Meliosma sellowii, Meliosma brasiliensis, Gordonia fruticosa e Ternstroemia brasiliensis. O pinheiro-brasileiro, A. angustifolia, é encontrado em uma área de floresta secundária adjacente ao Parque, tratando-se, provavelmente, da última população nativa remanescente da região e um extremo setentrional de distribuição geográfica natural da espécie. Outro aspecto notável é a redução da contribuição de Fabaceae para a flora das florestas neotropicais à medida que a altitude aumenta (Gentry 1995). Apesar de sua situação altitudinal próxima ao limite inferior das florestas alto-montanas (1.100m) sensu Oliveira-Filho & Fontes (2000), a proporção de espécies de Fabaceae no PEQRB (9,9%) é inferior à de florestas baixo-montanas da região (12-16%) (Pereira 2003).
Estrutura da comunidade arbórea - Os valores computados por Pereira (2003) para densidade e área basal por hectare das 20 áreas de floresta levantadas na região do Alto Rio Grande variam, respectivamente, entre 969 e 2.683 ind.ha-1 e 19,8 e 43,4m2ha-1. Portanto, a floresta do PEQRB, que registrou 1.643 ind.ha-1 e 20,5m2ha-1 tem densidade intermediária (10ª posição no ranking decrescente) e área basal baixa (18ª posição no ranking decrescente). Do ponto de vista ambiental, a floresta do PEQRB destaca-se no contexto das demais pelo caráter fortemente distrófico dos solos associado a textura comparativamente mais grossa, assim como longo histórico de impacto pelo fogo (Pereira 2003). Portanto, é possível que maior escassez de recursos hídricos e minerais, principalmente nas áreas mais distantes do curso d'água, contribua para biomassa relativamente mais baixa, expressa pela área basal. Os impactos do fogo também podem contribuir para reduzir a área basal. Contudo, são necessários estudos de produtividade líquida e dinâmica da comunidade arbórea para confirmar esta hipótese. Indicação que pode fortalecer esta hipótese vem das menores dimensões das árvores nos Neossolos Litólicos e Cambissolos Lépticos, que também apresentam drenagem mais forte, fertilidade mais baixa e menor proteção contra o fogo. Também observou-se concentração de árvores mais altas e de maior diâmetro mais próximo ao córrego.
Padrões de distribuição das espécies - A análise de correspondência canônica ordenou as parcelas de acordo com as correlações entre as abundâncias das espécies e as variáveis ambientais, resultando em clara separação das quatro classes de solo identificadas previamente. Esse fato confirma a consistência da classificação dos solos e do uso das quatro classes como definidoras de hábitats. Os Neossolos Litólicos associam-se positivamente às variáveis classes de drenagem e fator borda devido à sua concentração na borda da floresta, onde a declividade é acentuada e o solo mais raso, favorecendo forte drenagem. Os Cambissolos Lépticos mostraram-se semelhantes em relação à variável borda, mas diferem quanto à drenagem, pois estão situados mais próximos ao córrego. Os Cambissolos Típicos associaram-se mais fortemente às cotas mais elevadas e de forma negativa com a borda, porque altitude é mais elevada na maioria de suas parcelas que estão também inseridas no interior da floresta, formando ambiente com menos intensidade de luz. Nos Neossolos Regolíticos, a agregação foi mais evidente no lado oposto da variável drenagem provavelmente porque em metade de suas parcelas ocorre a passagem do córrego, tornando o ambiente mais úmido e mal drenado.
O fator borda foi a variável de maior correlação na distribuição das abundâncias das espécies. As bordas são, por definição, áreas de transição entre unidades de paisagem, tendo características que dependem das escalas espacial e temporal e das forças de interação que agem sobre esta unidade (Holland 1988). São áreas de hábitat mais expostas a perturbações externas e possuem, em geral, maior diversidade, densidade de espécie e produtividade biológica (Murcia 1985; Fortin 1994). As bordas de florestas influenciam na penetração de raios solares, modificando a qualidade e quantidade de luz, a temperatura, a evapotranspiração, a incidência de ventos e a redução da umidade do ar e do solo (Murcia 1995; Kapos et al. 1997). O fragmento florestal do PEQRB, apesar de ser o maior em área do município de Lavras, apresenta forma bastante alongada, aumentando seu perímetro e a ação do efeito borda, o que pode ainda incrementar sua diversidade de espécies. Segundo van den Berg & Oliveira-Filho (1999), nas florestas ripárias que têm bordas naturais com formações abertas, a heterogeneidade ambiental é maior na faixa ecotonal da borda, estando relacionada basicamente com a umidade do solo e a luminosidade. A variável drenagem também demonstrou forte influência na distribuição da abundância das espécies, o que pode interferir na disponibilidade de água e na fertilidade química dos solos. Nas florestas ripárias, geralmente a umidade e a fertilidade diminuem com o aumento da cota, principalmente quando a declividade é maior, conforme já detectado por alguns trabalhos feitos na região (van den Berg & Oliveira-Filho 1999; Botrel et al. 2002; Espírito-Santo et al. 2002; Souza et al. 2003).
A alta diversidade da comunidade arbórea da floresta do PEQRB provavelmente se deve à alta heterogeneidade ambiental incrementada pelo relevo acidentado e pelo pronunciado efeito borda nas faixas de transição entre floresta e ambientes abertos. Esta transição condiciona a existência do candeal, que representa verdadeiro ecótono entre a floresta semidecídua e os campos de altitude e rupestre. Estes fatos salientam a importância ecológica do Parque como valiosa amostra da vegetação primitiva da região do Alto Rio Grande, a qual merece atenção especial quanto à sua preservação.
Angiosperm Phylogeny Group II. 2003. An update of the Angiosperm phylogeny group classification for the orders and families of flowering plants: APG II Annals of the Botanical Journal of the Linnean Society 141(4): 399-436. [ Links ]
Botrel, R.T.; Oliveira-Filho, A.T.; Rodrigues, L.A. & Curi, N. 2002. Influência do solo e topografia sobre as variações da composição florística e estrutura da comunidade arbóreo-arbustiva de uma floresta estacional semidecidual em Ingaí, MG. Revista Brasileira de Botânica 25(2): 195-213. [ Links ]
Causton, D.R. 1988. An introduction to vegetation analysis, principles and interpretation. London, Unwin Hyman. [ Links ]
Curi, N.; Lima, J.M.; Andrade, H. & Gualberto, V. 1990. Geomorfologia, física, química e mineralogia dos principais solos da região de Lavras (MG). Ciência e Prática 14(2): 297-307. [ Links ]
Durigan, G.; Rodrigues, R.R. & Schiavini, I. 2000. A heterogeneidade ambiental definindo a metodologia de amostragem da floresta ciliar. Pp. 159-167. In: R.R. Rodrigues & H.F. Leitão-Filho (eds.). Matas ciliares: conservação e recuperação. São Paulo, EDUSP. [ Links ]
Eiten, G. 1982. Brazilian "Savannas". Pp. 25-47. In: B.J. Huntley & B.H. Walker (eds.). Ecology of tropical savannas. Berlim, Springer-Verlag. [ Links ]
EMBRAPA 1997. Manual de métodos de análise de solo. Rio de Janeiro, Empresa Brasileira de Pesquisa Agropecuária, Centro Nacional de Pesquisa de Solos. [ Links ]
EMBRAPA 1999. Sistema brasileiro de classificação de solos. Rio de Janeiro, Empresa Brasileira de Pesquisa Agropecuária, Centro Nacional de Pesquisa de Solos. [ Links ]
Espírito-Santo, F.D.B.; Oliveira-Filho, A.T.; Machado, E.L.M.; Souza, J.S.; Fontes, M.A.L. & Marques, J.J.G.S.M. 2002. Variáveis ambientais e a distribuição de espécies arbóreas em um remanescente de Floresta Estacional Semidecídua Montana no campus da Universidade Federal de Lavras, MG. Acta Botanica Brasilica 16(3): 331-356. [ Links ]
Fortin, M.J. 1994. Edge detection algorithms for two-dimensional ecological data. Ecology 75(4): 956-965. [ Links ]
Gentry, A.H. 1995. Patterns of diversity and floristic composition in neotropical montane forests. Pp. 103126. In: S.P. Churchill; H. Balslev; E. Forero & J.L. Luteyn (eds.). Biodiversity and conservation of Neotropical Montane Forests: Proceedings of Neotropical Montane Forest Biodiversity and Conservation Symposium. New York, The New York Botanical Garden. [ Links ]
Golfari, L. 1975. Zoneamento ecológico do estado de Minas Gerais para reflorestamento. Belo Horizonte, PNUD/FAO/IBDF. [ Links ]
Heltsche, J.F. & Forrester, N.E. 1983. Estimating species richness using the jackknife procedure. Biometrics 39(1): 1-12. [ Links ]
Holland, M.M. 1988. SCOPE/MAB technical consultations on landscape boundaries; report of a SCOPE/MAB workshop on ecotones. Pp. 47-106. In: A.F. di Castri; A.J. Hansen & M.M. Holland (eds.). A new look at ecotones: emerging international projects on landscape boundaries. Paris, Biology International. [ Links ]
Kapos, V.; Camargos, J.L.C.; Dande, G. 1997. Edge related changes in environment and plant response due to forest fragmentation in Central Amazonia. Pp. 45-54. In: W.F. Laurance & R.O. Bierregaard (eds.). Tropical forest remnats: ecology, management and conservation of fragmented communities. Chicago, The University of Chicago Press. [ Links ]
McCune, B. & Mefford, M.J. 1999. Multivariate analysis of ecological data. Gleneden Beach, MjM Software. [ Links ]
Mueller-Dombois, D. & Ellenberg, H. 1974. Aims and methods of vegetation ecology. New York, Wiley and Sons. [ Links ]
Murcia, C. 1995. Edges effects in fragmented forest: Implications for conservation. Trends in Ecology and Evolution 10(2): 58-62. [ Links ]
Oliveira-Filho, A.T.; Ratter, J.A. & Shepherd, G.J. 1990. Floristic composition and community structure of a central Brazilian gallery forest. Flora 184(2): 103-117. [ Links ]
Oliveira-Filho, A.T.; Almeida, R.J.; Mello, J.M. & Gavilanes, M.L. 1994a. Estrutura fitossociológica e variáveis ambientais em um trecho da mata ciliar do córrego dos Vilas Boas, Reserva Biológica do Poço Bonito, Lavras, MG. Revista Brasileira de Botânica 17(1): 67-85. [ Links ]
Oliveira-Filho, A.T.; Vilela, E.A.; Gavilanes, M.L. & Carvalho, D.A. 1994b. Comparison of the woody flora and soils of six areas of montane semideciduous forest in southern Minas Gerais, Brazil. Edinburgh Journal of Botany 51(3): 355-389. [ Links ]
Oliveira-Filho, A.T.; Mello, J.M. & Scolforo, J.R.S. 1997. Effects of past disturbance and edges on tree community structure and dynamics within a fragment of tropical semideciduous forest in southeastern Brazil over a five year period (1987-1992). Plant Ecology 131(1): 45-66. [ Links ]
Oliveira-Filho, A.T. & Fluminhan-Filho, M. 2000. Ecologia da vegetação do Parque Florestal Quedas do Rio Bonito. Cerne 5(2): 50-63. [ Links ]
Oliveira-Filho, A.T. & Fontes, M.A.L. 2000. Patterns of floristic differentiation among Atlantic forests in southeastern Brazil, and the influence of climate. Biotropica 32(4b): 793-810. [ Links ]
Oliveira-Filho, A.T.; Curi, N.; Vilela, E.A. & Carvalho, D.A. 2001. Variation in tree community composition and structure with changes in soil properties within a fragment of semideciduous forest in southeastern Brazil. Edinburgh Journal of Botany 58(1): 139-158. [ Links ]
Palmer, M.W. 1991. Estimating species richness: the second-order jackknife reconsidered. Ecology 72(4): 1512-1513. [ Links ]
Pereira, J.A.A. Efeitos dos impactos ambientais e da heterogeneidade ambiental sobre a diversidade e estrutura da comunidade arbórea de 20 fragmentos de florestas semidecíduas da região do Alto Rio Grande, Minas Gerais. Tese de doutorado. Universidade Fedeal de Minas Gerais, Belo Horizonte. [ Links ]
Rondon-Neto, R.M.; Botelho, S.A.; Fontes, M.A.L.; Davide, A.C. & Faria, J.M.R. 2000. Estrutura e composição florística da comunidade arbustivo-arbórea de uma clareira de origem antrópica, em uma floresta estacional semidecídua montana, Lavras, MG, Brasil. Cerne 6(2): 79-94. [ Links ]
Souza, J.S.; Espírito-Santo, F.D.B.; Fontes, M.A.L.; Oliveira-Filho, A.T. & Botezelli, L. 2003. Análise das variações florísticas e estruturais da comunidade arbórea de um fragmento de Floresta Semidecídua às margens do rio Capivari, Lavras, MG. Revista Árvore 27(2): 185-206. [ Links ]
ter Braak, C.J.F. 1987. The analysis of vegetation-environment relationships by canonical correspondence analysis. Vegetatio 69(1): 69-77. [ Links ]
ter Braak, C.J.F. 1988. CANOCO - A FORTRAN program for canonical community ordination by (Partial) (Detrended) (Canonical) correspondence analysis and redundancy analysis, version 2.1 Technical report LWA-88-2, TNO, Wageningen, Institute of Applied Computer Science. [ Links ]
ter Braak, C.J.F. 1995. Ordination. Pp. 91-173. In: R.H.G. Jongman; C.J.F. ter Braak & O.F.R. van Tongeren (eds.). Data analysis in community and landscape ecology. Cambridge, Cambridge University Press. [ Links ]
van den Berg, E. & Oliveira-Filho, A.T. 1999. Spatial partitioning among tree species within an area of tropical montane gallery forest in south-eastern Brazil. Flora 194(2/3): 249-266. [ Links ]
Veloso, H.P.; Rangel Filho, A.L.R. & Lima, J.C.A. 1991. Classificação da vegetação brasileira adaptada a um sistema universal. Rio de Janeiro, Instituto Brasileiro de Geografia e Estatística. [ Links ]
Recebido em 13/08/2003. Aceito em 31/03/2004
1 Trabalho integrante do Subprojeto Estratégias para conservação e manejo da biodiversidade em fragmentos de florestas semidecíduas, executado pela EMBRAPA Recursos Genéticos, a Universidade Federal de Lavras (UFLA) e Universidade de Brasília (UnB) e integrante do Projeto de Conservação e Utilização Sustentável da Diversidade Biológica Brasileira (PROBIO), supervisionado pelo Ministério do Meio Ambiente (MMA), gerido pelo Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) e apresentado pelo governo brasileiro ao GEF/BIRD.