Acessibilidade / Reportar erro

Variação espacial da estrutura da comunidade arbórea de um fragmento de floresta semidecídua em Piedade do Rio Grande, MG, Brasil

Spatial variation of tree community structure in a fragment of tropical semideciduous forest in Piedade do Rio Grande, Minas Gerais State, Brazil

Resumos

Foi realizado um levantamento da comunidade arbórea, do relevo e solos de um fragmento de floresta semidecídua com objetivo de verificar as possíveis correlações entre variações da estrutura comunitária e composição de espécies e variações do efeito borda e do regime de água e fertilidade dos solos. O fragmento possui 30 ha, com coordenadas 21º29'03" S e 44º06'05" W, em Piedade do Rio Grande, MG. Foram realizados um levantamento plani-altimétrico e uma classificação detalhada dos solos da floresta. Indivíduos arbóreos com diâmetro à altura do peito > 5 cm foram amostrados em 30 parcelas de 400 m² , sendo 12 parcelas situadas na borda do fragmento e 18 no interior do fragmento. Amostras do solo superficial (0 a 20 cm de profundidade) foram coletadas nas parcelas para análises das propriedades químicas e texturais. Foram registrados 1.778 indivíduos pertencentes a 174 espécies e três subgrupos de solos. A densidade de árvores e a diversidade de espécies foram mais elevadas nas parcelas de borda que nas de interior, mas nenhuma diferença foi detectada para área basal. Uma análise de correspondência canônica interativa dos padrões emergentes das variáveis ambientais, variáveis espaciais e abundância das espécies indicou que estas se distribuem no fragmento sob forte influência da proximidade espacial, do efeito borda, do regime de água do solo expresso pela cota topográfica e da fertilidade química dos solos expressa pela saturação por bases.

análise multivariada; correlação ambiente-vegetação; efeito borda; estrutura da comunidade arbórea; floresta tropical semidecídua; florística; fragmentação florestal


A survey of the tree community, topography and soils of a fragment of tropical semideciduous forest was carried out with the objective of assessing alleged correlations between the variations in community structure and species composition and variations in edge effect and soil water regime and chemical fertility. The fragment (30 ha of area) is situated at 21º29'03" S and 44º06'05" W, in Piedade do Rio Grande, Minas Gerais State, SE Brazil. A topographic survey and a detailed soil classification were accomplished in the forest area. Individual trees with diameter at breast height (dbh) > 5 cm were sampled in 30 plots with 400 m², 12 of which were placed at the forest edges and 18 in the forest interior. Samples of surface soil (0-20 cm of depth) were collected from each plot for analyses of chemical and textural properties. The survey registered 1,778 individuals belonging to 174 species, and three soil subgroups. Tree density and species diversity was higher in edge than in interior plots, but no difference was found for basal area. An interactive canonical correspondence analysis (CCA) of the patterns emerging from the environmental variables, spatial variables and species abundances indicated that these are distributed in the fragment under a strong influence of spatial proximity, edge effects, soil water regime expressed by topographical elevation and soil chemical fertility expressed by the bases saturation.

edge effect; environment-vegetation relationship; forest fragmentation; multivariate analysis; tree community structure; tropical semideciduous forest


ARTIGO

Variação espacial da estrutura da comunidade arbórea de um fragmento de floresta semidecídua em Piedade do Rio Grande, MG, Brasil

Spatial variation of tree community structure in a fragment of tropical semideciduous forest in Piedade do Rio Grande, Minas Gerais State, Brazil

Warley A. C. CarvalhoI; Ary Teixeira de Oliveira FilhoI,1 1 Autor para correspondência: ary@ufla.br ; Marco A. L. FontesI; Nilton CuriII

IUniversidade Federal de Lavras, Departamento de Ciências Florestais, 37200-000 Lavras, MG

IIUniversidade Federal de Lavras, Departamento de Ciências do Solo, 37200-000 Lavras, MG

RESUMO

Foi realizado um levantamento da comunidade arbórea, do relevo e solos de um fragmento de floresta semidecídua com objetivo de verificar as possíveis correlações entre variações da estrutura comunitária e composição de espécies e variações do efeito borda e do regime de água e fertilidade dos solos. O fragmento possui 30 ha, com coordenadas 21º29'03" S e 44º06'05" W, em Piedade do Rio Grande, MG. Foram realizados um levantamento plani-altimétrico e uma classificação detalhada dos solos da floresta. Indivíduos arbóreos com diâmetro à altura do peito > 5 cm foram amostrados em 30 parcelas de 400 m2 , sendo 12 parcelas situadas na borda do fragmento e 18 no interior do fragmento. Amostras do solo superficial (0 a 20 cm de profundidade) foram coletadas nas parcelas para análises das propriedades químicas e texturais. Foram registrados 1.778 indivíduos pertencentes a 174 espécies e três subgrupos de solos. A densidade de árvores e a diversidade de espécies foram mais elevadas nas parcelas de borda que nas de interior, mas nenhuma diferença foi detectada para área basal. Uma análise de correspondência canônica interativa dos padrões emergentes das variáveis ambientais, variáveis espaciais e abundância das espécies indicou que estas se distribuem no fragmento sob forte influência da proximidade espacial, do efeito borda, do regime de água do solo expresso pela cota topográfica e da fertilidade química dos solos expressa pela saturação por bases.

Palavras-chave - análise multivariada, correlação ambiente-vegetação, efeito borda, estrutura da comunidade arbórea, floresta tropical semidecídua, florística, fragmentação florestal

ABSTRACT

A survey of the tree community, topography and soils of a fragment of tropical semideciduous forest was carried out with the objective of assessing alleged correlations between the variations in community structure and species composition and variations in edge effect and soil water regime and chemical fertility. The fragment (30 ha of area) is situated at 21º29'03" S and 44º06'05" W, in Piedade do Rio Grande, Minas Gerais State, SE Brazil. A topographic survey and a detailed soil classification were accomplished in the forest area. Individual trees with diameter at breast height (dbh) > 5 cm were sampled in 30 plots with 400 m2, 12 of which were placed at the forest edges and 18 in the forest interior. Samples of surface soil (0-20 cm of depth) were collected from each plot for analyses of chemical and textural properties. The survey registered 1,778 individuals belonging to 174 species, and three soil subgroups. Tree density and species diversity was higher in edge than in interior plots, but no difference was found for basal area. An interactive canonical correspondence analysis (CCA) of the patterns emerging from the environmental variables, spatial variables and species abundances indicated that these are distributed in the fragment under a strong influence of spatial proximity, edge effects, soil water regime expressed by topographical elevation and soil chemical fertility expressed by the bases saturation.

Key words - edge effect, environment-vegetation relationship, forest fragmentation, multivariate analysis, tree community structure, tropical semideciduous forest

Introdução

Nas últimas décadas foi crescente o avanço nos estudos descritivos de comunidades florestais tropicais em todo mundo, de imprescindível importância para a conservação da biodiversidade. No Brasil, estas iniciativas têm sido cada vez mais encorajadas devido ao processo desordenado de ocupação do solo que, nas mais diversas regiões, têm transformado extensões florestais contínuas em fragmentos isolados (Botrel et al. 2002). Sendo o quinto maior país em extensão territorial, o Brasil apresenta grande variação em clima e solo, proporcionando vegetação variada com alta diversidade de habitats. O Estado de Minas Gerais, em particular, possui a maior variedade de formações vegetais do país, explicada por sua extensão e diversidade de condições geológicas, topográficas e climáticas. Porém, apesar da crescente constatação da necessidade da conservação dos recursos naturais, a cobertura vegetal primitiva foi reduzida a fragmentos esparsos recobrindo cerca de apenas 2% do território mineiro (CETEC 1983). A região do Alto Rio Grande, sul de Minas Gerais, encontra-se com sua cobertura vegetal reduzida a fragmentos de florestas e cerrados, a maioria deles muito perturbada por fogo, pecuária extensiva ou retirada seletiva de madeira para fins múltiplos (Botrel et al. 2002). As florestas semidecíduas desta região foram particularmente alteradas e reduzidas por conseqüência de quase sempre ocuparem os solos mais férteis e úmidos, mais propícios à agropecuária (Oliveira Filho et al. 1994).

A fragmentação florestal ocorre, na maioria das vezes, devido à substituição de parte da floresta por pastagens ou atividades agrícolas. As pressões antrópicas sobre estas fisionomias têm graves conseqüências, principalmente numa das mais antigas regiões de ocupação por colonizadores europeus e seus escravos africanos, como é a região do Alto Rio Grande. Portanto, é extremamente importante ampliar o conhecimento sobre a diversidade biológica remanescente nos atuais fragmentos florestais da região e sobre a organização espacial das comunidades vegetais e suas respostas à fragmentação da paisagem. Este conhecimento poderá gerar subsídios para o manejo e a conservação dos recursos naturais no longo prazo.

O presente trabalho foi desenvolvido em uma comunidade arbórea de um fragmento de floresta semidecídua no Município de Piedade do Rio Grande, sul de Minas Gerais. O objetivo foi verificar as possíveis correlações entre as variações da estrutura, da diversidade de espécies e da distribuição espacial das mesmas com a influência da borda do fragmento e de variações do substrato, aqui compreendido como o conjunto de variáveis topográficas e edáficas que determinam o status nutricional mineral e o regime de água dos solos. Também foi avaliada a autocorrelação espacial, ou seja os padrões de distribuição vinculados à mera proximidade entre indivíduos.

Material e métodos

Foi estudado um remanescente de floresta nativa com área de 30,2 ha, localizado no Município de Piedade do Rio Grande, MG, nas coordenadas 21º29'03" S e 44º06'04" WG (figura 1). O relevo é acidentado, com declividades entre 25 e 40º e altitudes variando entre 1.050 e 1.150 m. O embasamento rochoso é o gnaisse granítico. O clima é do tipo Cwb de Köppen (mesotérmico com verões brandos e suaves e estiagem de inverno). As médias anuais de precipitação e temperatura são 1.510 mm e 17,9 ºC, respectivamente, valores estes registrados para Barbacena (Brasil 1992), cidade a cerca de 50 km da área de estudos e altitude de 1.126 m. A vegetação é classificada como floresta semidecídua alto-montana a baixo-montana (segundo Oliveira Filho & Fontes 2000) ou floresta estacional semidecidual montana (segundo Veloso et al. 1991). Os poucos indivíduos de Araucaria angustifolia não são suficientes para caracterizar nem mesmo uma transição para florestas ombrófilas mistas. A floresta encontra-se em fase madura pois, segundo levantamento histórico, não houve corte raso ou retirada seletiva de grandes árvores nos últimos 150 anos.


O procedimento amostral foi estratificado e aleatório. A área amostral foi de 1,2 ha, constituída de 30 parcelas de 400 m2 distribuídas aleatoriamente por sorteio sobre ortofotocarta 1:10.000, sendo alocadas 18 parcelas de 20 ´ 20 m no interior do fragmento (estrato Interior) e 13 parcelas de 10 ´ 40 m nas bordas (estrato Borda). A localização das parcelas do Interior foi sorteada entre as coordenadas de uma grade de pontos de 20 ´ 20 m ajustada ao acaso à ortofotocarta. A situação das parcelas de Borda foi sorteada ao longo do perímetro da borda do fragmento e seu lado maior foi justaposto à borda de forma a maximizar a representação deste ambiente. Todos os indivíduos com DAP (diâmetro à altura do peito ou a 1,30 m de altura) > 5 cm equivalente a uma circunferência à altura do peito, CAP > 15,71 cm foram marcados e medidos em circunferência e altura. Indivíduos com caules múltiplos foram medidos quando a raiz da soma dos quadrados dos perímetros era > 15,71 cm, o que implica em uma área basal equivalente (Scolforo & Mello 2006). Nestes casos, caules CAP < 1 cm foram desprezados e, quando o resultado do cálculo era duvidoso no campo, o indivíduo era medido e, depois de feitos os cálculos, preservado na planilha ou então eliminado. Não foram incluídos indivíduos mortos e lianas. O levantamento florístico resultou de amostras de material botânico coletadas nos levantamentos das parcelas e em caminhamentos pelo fragmento, em visitas mensais de janeiro de 1999 a outubro de 2000. O material testemunho encontrase depositado no herbário da Universidade Federal de Lavras (Herbário ESAL). A classificação nas famílias seguiu o sistema APG II (2003).

Para descrição da estrutura das comunidades arbóreas do Interior e da Borda, foram extraídos, para cada espécie em cada estrato amostral, os valores de número de indivíduos, área basal total e número de parcelas onde ocorre. Testes de t de Student para amostras independentes foram empregados para comparar as médias de várias características da comunidade arbórea nas parcelas de Borda e Interior. Os pressupostos de populações independentes, normais e com variâncias iguais foram deliberadamente assumidos para realizar os testes de t, pois, de acordo com Zar (1996), "o pesquisador na área biológica não pode (...) ter sempre certeza de que estes pressupostos são corretos. Felizmente, muitos estudos têm mostrado que o teste de t é robusto o bastante para suportar consideráveis desvios de suas premissas teóricas". As variáveis estruturais comparadas foram densidade, área basal, diâmetros médio e máximo e alturas média e máxima. Foram preparadas distribuições de densidade de árvores por classes de diâmetro e altura para as subamostras de Borda e Interior, sendo utilizados intervalos de classe com amplitudes crescentes para compensar o decréscimo da densidade nas classes de tamanho maiores, conforme Botrel et al. (2002). O teste de t também foi empregado para comparar a densidade de árvores em cada classe de diâmetro e altura nos dois estratos amostrais.

Foram calculados, para a amostra total e seus dois estratos, o índice de diversidade de Shannon e a equabilidade de Pielou (Brower & Zar 1984) e os estimadores jackknife (de primeira e segunda ordem) do total de espécies (Heltsche & Forrester 1983, Palmer 1991). A densidade de espécies e o índice de Shannon na Borda e Interior foram comparados pelo teste de t para amostras independentes e pelo teste de t de Hutcheson, respectivamente (Zar 1996).

O solo de cada parcela foi classificado segundo o sistema da Embrapa (1999) até o nível de subgrupo (4º nível categórico), incluindo também os agrupamentos texturais e classes de drenagem. No centro de cada parcela, foram coletadas amostras de solo (0 a 20 cm de profundidade), as quais foram levadas ao Laboratório de Análise de Solos do Departamento de Ciências do Solo da UFLA para análises químicas e texturais. Foram obtidas as variáveis: pH, P, K, Ca2+, Mg2+, Al3+, saturação por bases (V), carbono (C), matéria orgânica (MO) e teores de areia, silte e argila, seguindo o protocolo da Embrapa (1997). Foram comparadas as variáveis químicas e texturais entre os três subgrupos de solos encontrados, por análises de variância complementadas por testes de Tukey-Kramer (Zar 1996). Foi realizado um levantamento planialtimétrico do fragmento que permitiu obter duas variáveis topográficas, cota média e desnível máximo, as quais foram utilizadas como meio de avaliação indireta das variáveis hidrológicas do solo (Oliveira Filho et al. 1994). Para avaliar a influência da borda na distribuição das espécies foi produzido um "fator de borda" para cada parcela, de acordo com o método utilizado por Espírito-Santo et al. (2002) e Souza et al. (2003), onde este é definido pelo comprimento linear do contorno da borda medido entre dois ou mais pontos de intercessão com o perímetro de um círculo cujo raio de 100 m tem origem no centro da parcela em questão.

Para correlacionar as variações quantitativas das espécies arbóreas nas parcelas com variáveis ambientais e espaciais, foi adotada a análise de correspondência canônica (CCA, ter Braak 1987) executada seguindo o procedimento interativo de Borcard et al. (1992), que exige a preparação de três matrizes. A primeira delas é a matriz das espécies que, no presente caso, consistiu do número de indivíduos (n) por parcela das 46 espécies com 10 ou mais indivíduos na amostra total. Os valores foram transformados em ln (n + 1) para compensar desvios causados pela baixa freqüência dos valores mais elevados (ter Braak 1995). A segunda é a matriz de variáveis ambientais que incluiu todas as variáveis topográficas e edáficas acrescidas do fator borda. A terceira é a matriz de variáveis espaciais que consistiu das coordenadas binárias do centro geométrico das parcelas, x e y, mais sete variáveis delas derivadas – x2, x3, y2, y3, x*y, x2*y e x*y2 – seguindo um modelo polinomial proposto por Borcard et al. (1992).

O passo seguinte foi a realização de quatro CCAs relacionando a matriz de espécies com (CCA1) variáveis ambientais, (CCA2) variáveis espaciais, (CCA3) variáveis ambientais e co-variáveis espaciais e (CCA4) variáveis espaciais e co-variáveis ambientais. Para tanto, foi utilizado o programa CANOCO 4.5 (ter Braak & Smilauer 2002). Nas duas primeiras CCAs, foi utilizada a rotina de seleção progressiva de variáveis ambientais associada a testes de permutação de Monte Carlo para verificar significância das mesmas. No caso das variáveis ambientais, apenas três permaneceram – cota, fator borda e saturação por bases – e, no das variáveis espaciais, só duas foram preservadas – y2 e y3. As mesmas variáveis foram utilizadas nas duas outras CCAs. Testes de permutação Monte Carlo foram aplicados às quatro CCAs e seus resultados foram empregados para estimar a proporção da variação dos dados das espécies que é explicada distintamente pelas variáveis ambientais e espaciais, bem como a proporção explicada indistintamente pelos dois conjuntos de variáveis e a variação remanescente que não é explicável por nenhuma das variáveis em questão.

As três variáveis ambientais foram também utilizadas para analisar a distribuição das abundâncias de cada uma das mesmas 46 espécies nas parcelas. Para comparar as abundâncias entre parcelas de Borda e Interior, foi utilizado o teste de t para amostras independentes e, no caso da cota média e saturação por bases, foram calculadas correlações de Pearson (Zar 1996). Os valores de abundância foram transformados pela função logarítmica porque esta aproximou sua distribuição do modelo normal.

Resultados

O levantamento florístico geral, que incluiu coletas fora das parcelas, registrou 65 famílias, 156 gêneros e 269 espécies, 174 das quais foram exclusivas das parcelas (tabelas 1 e 3 ). A família Fabaceae destacou-se pela maior riqueza de espécies (26) e gêneros (19), representando 9,7% do total de espécies registradas. Na segunda posição, as famílias Myrtaceae e Lauraceae apresentaram ambas 21 espécies. Outras famílias expressivas foram: Melastomataceae (15 espécies), Euphorbiaceae, Meliaceae e Rubiaceae (as três com 10 espécies). Na Borda, as cinco espécies de maior área basal foram Piptadenia gonoacantha, Machaerium villosum, Croton floribundus, Bauhinia forficata e Araucaria angustifolia; as mais numerosas foram P. gonoacantha, C. floribundus, B. forficata, Mollinedia widgrenii e Rollinia dolabripetala; as mais freqüentes foram M. widgrenii, Jacaranda macrantha, R. dolabripetala, P. gonoacantha e B. forficata. No Interior, exibiram maiores áreas basais Cariniana legalis, Coussapoa microcarpa, Metrodorea stipularis, Esenbeckia leiocarpa e Virola bicuhyba; apresentaram maiores densidades Galipea jasminiflora, M. stipularis, Anaxagorea silvatica, V. bicuhyba e Prunus myrtifolia; e foram mais freqüentes M. stipularis, A. silvatica, V. bicuhyba, Rollinia sylvatica e J. macrantha.

Foram registrados 1778 indivíduos com DAP > 5 cm, somando uma área basal de 50,73 m2. Dentre os parâmetros estruturais comparados, somente área basal não diferiu estatisticamente entre borda e interior do fragmento (tabela 2). A densidade foi significativamente maior na borda, enquanto os maiores diâmetros e alturas ocorreram no interior. As distribuições de densidade por classes de diâmetro e altura (figura 2) mostram valores estatisticamente superiores na borda para as menores classes de tamanho, com diferenças pronunciadas nas classes de 0 a 10 m de altura e de 5 a 33 cm de diâmetro. O índice de Shannon (H') da amostra total foi de 4,415 nats indivíduo-1 e a equabilidade de Pielou (J') foi de 0,855. Foram encontradas diferenças significativas entre Borda e Interior, tanto para H' como para a densidade de espécies por parcelas, ambos mais elevados na Borda (tabela 3 ), que também exibiu maior J', ou seja, menor dominância ecológica, e maiores projeções da riqueza de espécies pelos estimadores "jack-knife". Observe-se, porém, que estes estimadores foram ambos superados pelo levantamento florístico.


Os subgrupos de solos encontrados foram Latossolos Vermelho-Amarelos Ácricos típicos, Argissolos Vermelho-Amarelos Distróficos típicos e Argissolos Vermelho-Amarelos Eutróficos típicos, doravante referidos apenas como Latossolo, Argissolo D e Argissolo E. A distribuição dos solos no fragmento encontra-se na figura 3. As variáveis topográficas e edáficas das parcelas correspondentes aos três subgrupos de solo são fornecidas na tabela 4. A cota média foi significativamente mais elevada no Latossolo, que é acentuadamente drenado e concentrado no alto da encosta. Os bem drenados Argissolos se distribuem na encosta. O desnível do terreno não diferiu significativamente entre os solos. O pH e os teores de K e Ca foram significativamente mais altos no Argissolo E, enquanto Mg foi mais baixo no Argissolo D. Os teores de Al diminuíram significativamente do Latossolo para o Argissolo D e deste para o Argissolo E, enquanto a saturação por bases aumentou no sentido oposto. A matéria orgânica diferiu entre os extremos, mais alta nos Latossolos e mais baixa nos Argissolos D. Os Argissolos apresentaram mais areia e menos argila que o Latossolo e o silte foi mais elevado no Argissolo E. Em síntese, a microbacia apresenta um gradiente edáfico caracterizado por status nutricional com riqueza crescente do Latossolo para o Argissolo D e deste para Argissolo E, enquanto a textura evoluiu, no mesmo sentido, de maiores concentrações de argila e menores de areia (Latossolo) para concentrações crescentes de areia (Argissolos) e silte (Argissolo E).


As quatro CCAs apresentaram autovalores muito baixos para os dois primeiros eixos, isto é, todos < 0,35 (tabela 5). Isto significa que os gradientes são curtos em todos os casos, ou seja, há pouca substituição de espécies e a maioria delas se distribui ao longo dos gradientes variando apenas sua abundância (ter Braak 1995). As CCAs ambiental e espacial explicaram proporções semelhantes da variação total das espécies, entre 21–22%, o que salienta que a proximidade espacial tem um peso semelhante ao das variáveis ambientais quando se procura explicar a distribuição da abundância das árvores no fragmento. Além disso, conforme indica a figura 4, (a) apenas 12,68% da variação total é puramente ambiental e independente do espaço (CCA3); (b) 9,24% da variação é indistintamente espacial e ambiental (CCA1 – 3 ou CCA2 – 4) , ou seja, deve-se à estrutura espacial do próprio ambiente; (c) 11,84% da variação é puramente espacial e independente do ambiente (CCA4). Em decorrência disso, pouco mais de um terço da variação total (33,77%) foi explicada por ambiente e espaço em conjunto (CCA1 + 4 ou CCA2 + 3), restando (d) 66,23% de variação das espécies que é estocástica ou não explicável pelas variáveis disponíveis. Esta expressiva variância remanescente é comum em dados de vegetação e não prejudica a significância das relações espécie-ambiente (ter Braak 1988). Com efeito, os dois primeiros eixos das quatro CCAs explicaram proporções elevadas das variações da correlação espécie-ambiente, > 50% e > 80% no primeiro e segundo eixos, respectivamente (tabela 5). As correlações de Pearson entre espécies e variáveis ambientais e espaciais também foram elevadas nos dois primeiros eixos, variando entre 0,80 e 0,93. Além disso, em todas as CCAs, os testes de permutação de Monte Carlo indicaram diferenças altamente significativas entre as correlações encontradas e as geradas pelo acaso, tanto para o primeiro eixo canônico como para o conjunto de eixos canônicos (tabela 5).


Entre as variáveis ambientais, a cota se destaca como mais fortemente correlacionada com o primeiro eixo das duas CCAs ambientais, mas também tem peso expressivo do segundo eixo da CCA1 (tabela 5). O fator de borda tem correlações mais fortes com o segundo eixo da CCA1 e ambos os eixos da CCA3. A saturação por bases (V) produziu correlações relativamente mais fracas na CCA1 (ca ± 0,45), mas alcançou valor expressivo no segundo eixo da CCA3. As variáveis espaciais, y2 e y3, produziram correlações entre ± 0,5 e ± 0,85 em todos os casos exceto o de y3 no segundo eixo da CCA2, que ficou abaixo de ± 0,4. É notório que ambas as variáveis espaciais selecionadas pela sua correlação significativa com a variância da matriz de espécies tenham sido derivadas da dimensão correspondente ao eixo norte-sul e ao sentido da declividade predominante na encosta. Isto sugere que o efeito da proximidade das plantas sobre sua própria distribuição de abundâncias foi significativo no sentido paralelo à declividade, mas não pôde ser detectado no sentido perpendicular à mesma. A redundância de y com a cota pode, provavelmente, explicar boa parte das variações ambientais espacialmente estruturadas representadas pelo setor (b) da figura 3.

Para representar graficamente a análise, foi escolhida a CCA1 para explorar apenas os efeitos das variáveis ambientais. Os diagramas da figura 5 mostram a ordenação de espécies e parcelas em separado, para maior clareza. As coordenadas de ordenação das espécies e parcelas representam o ótimo aproximado de uma distribuição unimodal das mesmas no espaço multidimensional definido pelas variáveis ambientais. Estas, por sua vez, são representadas por retas divergentes indicando a direção predominante da magnitude crescente de cada variável no espaço de ordenação (ter Braak 1988). A ordenação sugere forte separação entre parcelas de Borda e Interior. As correlações com as três variáveis ambientais são muito claras, assim como o efeito borda-interior. A CCA sugere um gradiente onde espécies como Anaxagorea silvatica, Virola bicuhyba, Metrodorea stipularis, Trichilia catigua, Galipea jasminiflora, Rollinia sylvatica, Prunus myrtifolia, Cestrum schlechtendalii, Pisonia zapallo e Diatenopteryx sorbifolia tendem a ocorrer em maior abundância no interior, enquanto espécies como Hyptidendron asperrimum, Clethra scabra, Myrcia splendens, Casearia sylvestris, Vismia brasiliensis, Machaerium villosum, Protium widgrenii, Myrsine umbellata, Piptocarpha axillaris e Croton floribundus concentram-se mais na borda.


As comparações das abundâncias das espécies entre as parcelas da Borda e Interior (tabela 6) mostram que metade das 46 espécies foi significativamente mais abundante em um dos estratos, sendo 19 na Borda e quatro no Interior. Tais espécies foram também distribuídas no diagrama da CCA nos extremos da dicotomia borda-interior (figura 5). Dentre as mesmas 46 espécies, 19 mostraram alguma correlação significativa entre suas abundâncias nas parcelas e as variáveis cota e saturação por bases (tabela 6). Para cota, houve 12 correlações positivas e três negativas e, para saturação por bases, cinco positivas e duas negativas. Tais correlações foram coerentes com a distribuição das espécies no diagrama da CCA (figura 5). Das 46 espécies analisadas, 14 (30%) não mostraram nenhum padrão de distribuição significativamente relacionado com as duas variáveis ambientais nem com a dicotomia borda-interior.

Discussão

No contexto dos padrões de variações florísticas e ambientais das florestas Atlânticas sensu lato do Sudeste do Brasil descritos por Oliveira Filho & Fontes (2000), a composição de espécies arbóreas do fragmento estudado em Piedade do Rio Grande pode ser vista como uma expressão de duas transições florísticas: a primeira delas entre florestas ombrófilas e semidecíduas e a segunda entre florestas baixo-montanas e alto-montanas. Há espécies características de florestas ombrófilas, como Coussapoa microcarpa, Sloanea guianensis, Vochysia schwackeana, Virola bicuhyba, Picrasma crenata e Marlierea excoriata; espécies comuns de florestas semidecíduas, como Machaerium villosum, Ceiba speciosa, Luehea grandiflora, Platypodium elegans e Copaifera langsdorffii; espécies características de florestas alto-montanas, como Capsicodendron dinisii, Clethra scabra, Jacaranda subalpina, Euplassa organensis e Araucaria angustifolia; e espécies típicas de florestas baixo-montanas como Protium widgrenii, Vismia brasiliensis, Piptocarpha axillaris e Calyptranthes clusiifolia. A observação é reforçada pelas principais famílias encontradas, pois o perfil florístico do fragmento se aproxima das florestas alto-montanas no expressivo número de espécies das famílias Myrtaceae, Lauraceae, Melastomataceae e Solanaceae. Contudo, a abundância de Fabaceae se opõe ao caráter alto-montano, conferido pelas famílias anteriores (Oliveira Filho & Fontes 2000). O fragmento localiza-se em uma vertente com face de exposição voltada para o sul, situação denominada no Brasil de "terreno noruega", sendo este, em comparação com a face norte, mais frio, úmido e pouco ensolarado. Tal situação, segundo Hugget (1995), favorece espécies relativamente menos exigentes de luz e mais exigentes em umidade. É possível que esta particularidade tenha favorecido a ocorrência de espécies características de florestas de maior altitude na região. Semelhante a outras florestas da região situadas em maiores altitudes, o fragmento contém ainda uma expressiva abundância de Orchidaceae (e.g. Maxillaria sp., Oeceoclades maculata (Lindley) Lindley, Corymborkis flava (Swartz) Kuntze, Oncidium spp., Polystachya spp.), Bromeliaceae (e.g. Tillandsia spp. e Billbergia spp.), briófitas e pteridófitas recobrindo os troncos das árvores.

A riqueza de espécies registrada pode ser considerada como relativamente alta no contexto de 20 levantamentos semelhantes realizados na região (Pereira et al. 2007). A diversidade de espécies expressa pelo índice de Shannon (H' = 4,415 nats indivíduo-1) corresponde ao segundo valor mais elevado dos mesmos levantamentos (Pereira et al. 2007, Dalanesi et al. 2004). Tal fato pode ser atribuído, em grande parte, ao encontro das tipologias florestais discutido acima. Além disso, deve haver também uma expressiva contribuição da considerável heterogeneidade ambiental, propiciada pelo gradiente edáfico e topográfico bem como pelo efeito borda (Metzger 1999, Murcia 1995).

A diferença altitudinal de cerca de 100 m registrada no fragmento certamente contribuiu para a heterogeneidade ambiental ao propiciar a formação de gradientes edáficos e topográficos com reflexos na estrutura e composição de espécies da comunidade arbórea. De fato, a distribuição dos três subgrupos de solo no fragmento caracterizou uma catena com significativas variações das propriedades químicas e texturais do solo. É comum que conteúdo de água do solo de uma catena cresça do topo para a base da encosta (Resende et al. 1995). Isto se aplica ao presente caso, com Latossolo situado na parte mais alta do terreno. A alta permeabilidade, peculiar à maioria dos Latossolos, é ainda incrementada pela posição topográfica elevada (Curi et al. 1993). Portanto, o déficit hídrico no período mais seco do ano é, provavelmente, mais pronunciado no Latossolo que nos Argissolos, situados mais abaixo na encosta. Nestes últimos, o horizonte B-textural em geral desacelera a percolação da água, mas também aumenta a sua disponibilidade para as plantas durante a estiagem devido à maior capacidade de armazenamento (Buckman & Brady 1969, Resende et al. 1988). Isto provavelmente explica a concentração de espécies de maior deciduidade na parte mais alta do fragmento, o que é denunciado pelo claro contraste visual do dossel da floresta entre o topo e a encosta durante a estação seca. Também ocorre um incremento em fertilidade química do Latossolo para os Argissolos. Contudo, o incremento em fertilidade do Argissolo D para o Argissolo E só ocorreu em dois setores distintos da encosta. Na calha central da microbacia, o Argissolo D se estende até o fundo do vale, talvez como resultado de uma mais pronunciada erosão e lixiviação que nas encostas vizinhas onde é substituído pelo Argissolo E.

Boa parte da distribuição das espécies mostrou-se espacialmente estruturada e sem conexão com as variáveis ambientais utilizadas, ou seja, foi influenciada pela própria proximidade entre indivíduos co-específicos, resultando em autocorrelação espacial. Aglomerados populacionais são comuns em florestas tropicais e podem se formar tanto devido à dispersão concentrada em curtas distâncias, bem como ao histórico particular de cada peça que compõe o mosaico sucessional da floresta, fenômenos estes que têm pouca ou nenhuma relação com as características do substrato (Oldeman 1990). No entanto, uma outra fração da autocorrelação espacial das espécies foi inseparável da própria variação ambiental e, além disso, concorrente com a disposição espacial do gradiente catenário (norte-sul), sugerindo que os efeitos puramente espaciais participam, de fato, menos que os ambientais na explicação da variância da matriz de espécies.

A distribuição de várias espécies no fragmento foi significativamente correlacionada com as variações da fertilidade dos solos e da elevação topográfica, esta última provavelmente refletindo também um gradiente de disponibilidade de água nos solos. A influência da umidade e fertilidade dos solos na distribuição de espécies arbóreas em florestas tropicais tem largo subsídio na literatura (e.g. Botrel et al. 2002, Clark et al. 1998, Duivervoorden & Lips 1995, Espírito-Santo et al. 2002, Newbery & Proctor 1984, Newbery et al. 1986, Oliveira Filho et al. 1994, 1997b, 2001). Na área de estudo, cotas mais elevadas do terreno corresponderam, de modo geral, a sítios com menor estoque de água e nutrientes minerais. Espécies com maior abundância nestes sítios provavelmente apresentam maior capacidade de sobreviver à menor disponibilidade destes recursos. Entre estas estão Calyptranthes clusiifolia, Casearia sylvestris, Clethra scabra, Cordia sellowiana, Cupania vernalis, Ixora brevifolia, Protium widgrenii, Myrcia splendens, Machaerium villosum e Tapirira obtusa, todas elas já correlacionadas por outros autores com sítios de melhor drenagem e, ou, baixa fertilidade química (Botrel et al. 2002, Lorenzi 1992, 1998, Oliveira Filho et al. 1994, Souza et al. 2003). Um número menor de espécies apresentou o padrão oposto, ou seja, foram significativamente mais abundantes nos sítios mais úmidos e férteis da encosta, provavelmente exigindo uma maior disponibilidade de água e, ou, nutrientes minerais. Algumas delas, como Mollinedia widgrenii, Cabralea canjerana, Bauhinia forficata, Machaerium nictitans e Trichilia catigua já haviam sido caracterizadas desta forma por outros autores (Lorenzi 1992, 1998, Oliveira Filho et al. 1994, 1997a, Souza et al. 2003).

Apesar de ter sido detectada uma relação entre a distribuição das espécies e as variações catenárias dos solos, foi o efeito borda que mais fortemente se relacionou com variações espaciais das características da comunidade arbórea, influenciando na estrutura da comunidade e na diversidade e padrão de distribuição de espécies. As bordas de fragmentos florestais são áreas de habitat mais expostas a perturbações externas e possuem, em geral, maior diversidade e densidade de espécies e maior produtividade biológica (Fortin 1992, 1994). No fragmento estudado, a borda também apresentou maior diversidade em todos os aspectos analisados (riqueza, densidade e eqüabilidade). A natureza transicional dos ambientes de borda parece ser a razão principal deste incremento da diversidade (Fortin 1992, 1994, Murcia 1995, Metzger 1999). A fragmentação florestal aumenta a proporção de áreas de borda, o que ocasiona alterações ambientais importantes. A maior penetração dos raios solares e do vento aumenta a temperatura e a evapotranspiração, reduz a umidade do ar e do solo e cria condições mais favoráveis ao stress hídrico (Kapos et al. 1997). Tudo isto contribui para a instalação de biotas de composição e estrutura distintas daquelas do interior dos fragmentos (Metzger 1999).

A dicotomia em estrutura e composição de espécies da comunidade arbórea entre a borda e o interior de fragmentos florestais também é influenciada pelo histórico do processo de formação e manutenção da borda (Murcia 1995, Oliveira Filho et al. 1997b). Como, no presente caso, as bordas foram estabelecidas entre 1948 e 1952, pode-se afirmar que, após 50 anos, uma comunidade arbórea bastante diferenciada da original já se estabeleceu em resposta às alterações ambientais. Deve-se ter em mente que as alterações ambientais nas bordas são multifatoriais, isto é, envolvem um conjunto de variáveis que incluem condições de luz, temperatura, umidade atmosférica e do solo, exposição ao vento etc. Contudo, a maior exposição à luz tem um peso muito forte. No fragmento estudado, tais alterações ambientais provavelmente explicam a maior densidade de árvores de menor porte, favorecidas pela maior penetração de luz pela borda. Inversamente, a menor densidade de pequenas árvores no interior se deve ao maior sombreamento do sub-bosque pelas grandes árvores. Espécies encontradas em maior abundância na borda como Bauhinia forficata, Casearia sylvestris, Casearia lasiophylla, Clethra scabra, Croton floribundus, Luehea divaricata, Machaerium villosum, Myrcia splendens, Myrsine umbellata, Piptadenia gonoacantha e Zanthoxylum rhoifolium, são reconhecidas como muito exigentes de luz para seu estabelecimento e crescimento (Botrel et al. 2002, Lorenzi 1992, 1998, Souza et al. 2003). Por outro lado, espécies com distribuição preferencial pelo interior, como Galipea jasminiflora, Prunus myrtifolia e Virola bicuhyba, são reconhecidas como relativamente tolerantes à sombra (Lorenzi 1992, 1998, Oliveira Filho et al. 1994, 1997a). Algumas espécies fugiram ao padrão esperado, como Mollinedia widgrenii e Calyptranthes clusiifolia, que foram significativamente mais abundantes nas bordas apesar de serem mais comuns no sub-bosque sombrio das florestas da região (Oliveira Filho et al. 1994). No entanto, foi observado no campo que a grande maioria das árvores destas espécies estavam estabelecidas no ambiente mais sombrio adjacente às bordas e não nas bordas em si. Portanto, outros fatores, tais como dispersão ou condições luminosas particulares, devem estar envolvidos na distribuição preferencial destas espécies em um hábitat tão particular.

Como é de praxe em estudos do gênero, ainda restou uma proporção das variações da matriz de abundância das espécies que não pôde ser explicada pelo substrato, pelo efeito borda e nem pela autocorrelação espacial. Esta variação remanescente é muito maior que a explicada e não deve ser meramente desprezada como 'ruído'. Ela é relevante e, para explicá-la, só podemos, infelizmente, especular sobre o papel de variáveis diferenciadoras desconhecidas ou então sobre um universo neutro onde eventos estocásticos governariam a dinâmica de indivíduos igualmente adaptáveis aos habitats em questão (Hubbell 2001). Estas especulações só evidenciam o quanto ainda está longe uma compreensão mais completa da distribuição e diversidade de espécies em florestas tropicais.

Os padrões de diferenciação interna das características da comunidade arbórea do fragmento estudado refletem considerável heterogeneidade ambiental que testemunha a riqueza e complexidade da primitiva cobertura florestal contínua na região. A fragmentação florestal trouxe a intromissão do efeito borda, que incrementou fortemente a heterogeneidade ambiental pré-existente vinculada às variações do substrato. Este fato salienta a importância do tamanho dos fragmentos florestais para promover a conservação de amostras das comunidades primitivas, tamanho este que dever ser tanto maior quanto maior for a heterogeneidade do substrato em que se desenvolve a floresta a ser conservada.

Agradecimentos – Este trabalho integrou o 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). Os autores agradecem às taxonomistas Inês Cordeiro, Maria Lúcia Kawasaki e Lúcia Rossi, pelo atencioso auxílio na identificação do material botânico.

(recebido: 21 de janeiro de 2005; aceito: 10 de maio de 2007)

  • APG. 2003. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG II. Botanical Journal of the Linnean Society 141:399-436.
  • BORCARD, D., LEGENDRE, P. & DRAPEAU, P. 1992. Partialling out the spatial component of ecological variation. Ecology 73:1042-1055.
  • 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:195-213.
  • BRASIL. 1992. Normais climatológicas (1961-1990). Departamento Nacional de Meteorologia, Brasília.
  • BROWER, J.E. & ZAR, J.H. 1984. Field and laboratory methods for general ecology. W.M.C. Brow, Dubuque.
  • BUCKMAN, H.O. & BRADY, N.C. 1969. The nature and properties of soils. Macmillan, London.
  • CETEC - Fundação Centro Tecnológico de Minas Gerais. 1983. Diagnóstico ambiental do Estado de Minas Gerais. Séries de publicações Técnicas. Belo Horizonte.
  • CLARK, D.B., CLARK, D.A. & READ, J.M. 1998. Edaphic variation and the mesoscale distribution of tree species in a neotropical rain forest. Journal of Ecology 86:101-112.
  • CURI, N., LARACH, J.O.I, KÄMPF, N., MONIZ, A.C. & FONTES, L.E.F. 1993. Vocabulário de ciência do solo. Campinas.
  • DALANESI, P.E., OLIVEIRA FILHO, A.T. & FONTES, M.A.L. 2004. Flora e estrutura do componente arbóreo da floresta do Parque Florestal Quedas do Rio Bonito, Lavras, MG, e correlações entre distribuição das espécies e variáveis ambientais. Acta Botanica Brasilica 18:737-757.
  • DUIVERVOORDEN, J.F. & LIPS, J.M. 1995. A land-ecological study of soils, vegetation, and plant diversity in Colombian Amazonia. Tropenbos series 12. 1st ed., The Tropenbos Foundation, Wageningen.
  • EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária, Centro Nacional de Pesquisa de Solos. 1997. Manual de métodos de análises de solo. 2Ş ed., Embrapa, Rio de Janeiro.
  • EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária, Centro Nacional de Pesquisa de Solos. 1999. Sistema brasileiro de classificação de solos. 1Ş ed., Embrapa, Rio de Janeiro.
  • 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:331-356.
  • FORTIN, M.J. 1992. Detection of ecotones: definition and scaling factors. Tese de doutorado, University of New York, New York.
  • FORTIN, M.J. 1994. Edge detection algorithms for two-dimensional ecological data. Ecology 75: 956-965.
  • HELTSCHE, J.F. & FORRESTER, N.E. 1983. Estimating species richness using the jackknife procedure. Biometrics 39:1-12.
  • HUBBELL, S. P. 2001. The unified neutral theory of biodiversity and biogeography. Princeton Univ. Press, New Jersey.
  • HUGGET, R.J. 1995. Geoecology, an evolutionary approach. Routledge, London.
  • KAPOS, V., CAMARGOS, J.L.C. & DANDE, G. 1997. Edge related changes in environment and plant response due to forest fragmentation in Central Amazonia. In Tropical forest remnants: ecology, management and conservation of fragmented communities (W.F. Laurance & R.O. Bierregaard, eds.). The University of Chicago Press, Chicago, p.45-54.
  • LORENZI, H. 1992. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas nativas do Brasil. v.1, Editora Plantarum, Nova Odessa.
  • LORENZI, H. 1998. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas nativas do Brasil. v.2, Editora Plantarum, Nova Odessa.
  • METZGER, J.P. 1999. Estrutura da paisagem e fragmentação: análise bibliográfica. Anais da Academia Brasileira de Ciências 71:445-463.
  • MURCIA, C. 1995. Edges effects in fragmented forest: implications for conservation. Trends in Ecology and Evolution 10:58-62.
  • NEWBERY, D.McC. & PROCTOR, J. 1984. Ecological studies in four contrasting lowland rain forests in Gunugu Mulu National Park, Sarawak. IV. Association between tree distribution and soil factors. Journal of Ecology 72:475-493.
  • NEWBERY, D.McC., GARTLAN, J.S., MCKEY, D.B. & WATERMAN, P.G. 1986. The influence of drainage and soil phosphorus on the vegetation of Douala-Edea Forest reserve, Cameroun. Vegetatio 65:149-162.
  • OLDEMAN, R. A. A. 1990. Dynamics in tropical rain forests. In Tropical forests Botanical dynamics, speciation and diversity (L.B. Holm-Nielsen, I. C. Nielsen & H. Balslev, eds.). Academic Press, London, pp. 3-21.
  • 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:793-810.
  • OLIVEIRA FILHO, A.T., VILELA, E.A., CARVALHO, D.A. & GAVILANES, M.L. 1994. Effects of soils and topography on the distribution of tree species in a tropical riverine forest in south-eastern Brazil. Journal of Tropical Ecology 10:483-508.
  • OLIVEIRA FILHO, A.T., CURI, N., VILELA, E.A. & CARVALHO, D.A. 1997a. Tree species distribution along soil catenas in a riverside semideciduous forest in Southeastern Brazil. Flora 192:47-64.
  • OLIVEIRA FILHO, A.T., MELLO, J.M. & SCOLFORO, J.R.S. 1997b . Effects of past disturbance and edges on tree community structure and dynamics within a fragment of tropical semideciduous forest in south-eastern Brazil over a five-year period (19871992). Plant Ecology 131:45-66.
  • 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 south-eastern Brazil. Edinburgh Journal of Botany 58:139-158.
  • PALMER, M.W. 1991. Estimating species richness: the second-order jackknife reconsidered. Ecology 72:1512-1513.
  • PEREIRA, J.A.A., OLIVEIRA FILHO, A.T. & LEMOS FILHO, J.P. 2007. Environmental heterogeneity and disturbance by humans control much of the tree species diversity of Atlantic montane forest fragments in SE Brazil. Biodiversity and Conservation (no prelo: artigo 9063, BIOC1721, 24 p.).
  • RESENDE, M., CURI, N. & SANTANA, D.P. 1988. Pedologia e fertilidade do solo: interações e interpretações. Editora UFLA, Lavras.
  • RESENDE, M., CURI, N., REZENDE, S.B. & CORRÊA, G.F. 1995. Pedologia: bases para distinção de ambiente. Núcleo de Estudo de Planejamento e Uso da Terra, Viçosa.
  • SCOLFORO, J.S., MELLO, J.M. 2006. Inventário Florestal, v. 1. Editora UFLA, Lavras.
  • 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:185-206.
  • TER BRAAK, C.J.F. 1987. The analysis of vegetation-environment relationship by canonical correspondence analysis. Vegetatio 69:69-77.
  • TER BRAAK, C.J.F. 1988. Canoco a fortran program for canonical community ordinatin by (Partial) (Detrended) (Canonical) correspondence analysis and redundancy analysis, version 2.1. Technical report LWA-88-2 TNO, Wageningen.
  • TER BRAAK, C.J.F. 1995. Ordination. In Data analysis in community and landscape ecology (R.H.G. Jongman, C.J.F. Ter Braak & O.F.R. Van Tongeren, eds.). Cambridge University Press, Cambrigde, p.91-173.
  • TER BRAAK, C.J.F. & P. SMILAUER. 2002. CANOCO Reference manual and CanoDraw for Windows user's guide: Software for Canonical Community Ordination (version 4.5). Microcomputer Power, Ithaca.
  • VELOSO, H.P., RANGEL FILHO, A.L.R. & LIMA, J.C. 1991. Classificação da vegetação brasileira adaptada a um sistema universal. IBGE, Rio de Janeiro.
  • ZAR, J.H. 1996. Biostatistical analysis. 3rd ed., Prentice-Hall, New Jersey.
  • 1
    Autor para correspondência:
  • Datas de Publicação

    • Publicação nesta coleção
      28 Set 2007
    • Data do Fascículo
      Jun 2007

    Histórico

    • Aceito
      10 Maio 2007
    • Recebido
      21 Jan 2005
    Sociedade Botânica de São Paulo Caixa Postal 57088, 04089-972 São Paulo SP - Brasil, Tel.: (55 11) 5584-6300 - ext. 225, Fax: (55 11) 577.3678 - São Paulo - SP - Brazil
    E-mail: brazbot@gmail.com