Accessibility / Report Error

Expansão da vegetação arbustivo-arbórea em áreas abertas alto-montanas adjacentes a fragmentos florestais, no Planalto Sul Catarinense

Expansion of tree-shrub vegetation in upper-montane open areas adjacent to forest fragments, in “Planalto Sul Catarinense” region

Resumo

O objetivo deste estudo foi conhecer os padrões demográficos da regeneração natural do componente arbustivo-arbóreo e determinar se esses são influenciados pelas variáveis ambientais, em três áreas abertas em regeneração natural (diâmetro na altura do peito < 5 cm), em contato com fragmentos florestais de Floresta Ombrófila Mista Alto-Montana, no Planalto Sul Catarinense. Foram instaladas três transecções, uma por área, de 20 x 100 m, subdividas em subparcelas de 10 x 10 m, perpendiculares à borda dos fragmentos, com sua extensão cobrindo a área de campo. Amostraram-se todos os indivíduos regenerantes arbustivo-arbóreos com 20 cm ou mais de altura. Em cada subparcela, foram coletados os seguintes dados ambientais: propriedades físico-químicas dos solos, grau de compactação dos solos e cota altimétrica. Um ano após o inventário, foram obtidas as taxas de dinâmica para as três comunidades e para as populações, e os padrões encontrados foram relacionados com as variáveis ambientais. Os dados foram analisados por meio de testes de Kruskal-Wallis, testes de comparações múltiplas, teste de qui-quadrado aplicado à tabela de contingência, teste de correlação, Análise de Componentes Principais (PCA) e modelos lineares mistos. Os resultados demonstraram que, no período estudado, o processo de regeneração natural do componente arbustivo-arbóreo encontra-se em fase inicial e foi caracterizado pelo incremento na abundância de indivíduos, estabilidade na riqueza, independência espacial entre a mortalidade e o recrutamento e limitação do recrutamento em função da baixa fertilidade natural do solo e elevada disponibilidade de Al.

Palavras-chave:
Floresta Ombrófila Mista Alto-Montana; Padrões de dinâmica; Análise ambiental

Abstract

The study aimed to understand the dynamic patterns of natural regeneration of shrub-tree species and determine whether these patterns are influenced by environmental variables, in three open areas in natural regeneration (diameter at breast height < 5 cm), in contact to upper-montane Araucaria Forest fragments, located in the Southern Plateau of Santa Catarina state. Three 20x100 m transects, subdivided into sub-plots of 10x10 m, were installed in each area, perpendicular to the edge of the fragments with their length covering the grassland area. All shrub-tree regenerating individuals greater than or equal to 20 cm height were sampled. In each sub-plot, the environmental data were collected: soil physico-chemical properties, the degree of soil compaction and the altimetric quota. One year after the first inventory, the dynamics rates were obtained for the three community and populations, and the found patterns related to environmental variables. The data were analyzed through Kruskal-Wallis tests, multiple comparison tests, chi-square test, applied to a contingency table, Principal Component Analysis (PCA) and linear mixed models. The results demonstrated that during the evaluated period, the natural regeneration process of the shrub-tree component is found in an early stage and it was characterized by abundance increment, richness stability, spatial independency between mortality and recruitment, and recruitment limitation in function of low soil natural fertility and high Al availability.

Keywords:
Highlands Araucaria Forest; Dynamic patterns; Environmental analysis

Introdução

A Floresta Ombrófila Mista (FOM) Alto-Montana faz parte de um mosaico formado por florestas e áreas de vegetação aberta, constituídas por campos naturais (Estepes), considerados como uma vegetação proveniente de períodos glaciais do Quaternário, mais secos e frios (BEHLING et al., 2007BEHLING, H. et al. Late-Holocene fire history in a forest-grassland mosaic in southern Brasil: implications for conservation. Applied Vegetation Science, Hoboken, v. 10, n. 1, p. 81-90, 2007.; BEHLING; PILLAR, 2007BEHLING, H.; PILLAR, V. D. Late quaternary vegetation, biodiversity and fire dynamics on the southern brazilian highland and their implication for conservation and management of modern araucaria forest and grassland ecosys tems. Philosophical Transactions of the Royal Society B, London, v. 362, n. 1478, p. 243-251, 2007.; SILVA; ANAND, 2011SILVA, L. C. R.; ANAND, M. Mechanisms of Araucaria (Atlantic) forest expansion into southern Brazilian grasslands. Ecosystems, Nova York, v. 14, n. 8, p. 1354-1371, 2011.), e por áreas desmatadas para dar origem a atividades antrópicas. De acordo com Jeske-Pieruschka et al. (2013JESKE-PIERUSCHKA, V. et al. New insights into vegetation, climate and fire history of southern Brazil revealed by a 40,000 year environmental record from the State Park Serra do Tabuleiro. Vegetation History and Archaeobotany, New York, v. 22, n. 4, p. 299-314, 2013.), nos locais em que as condições edáficas são favoráveis, o clima atual da região sul favorece a expansão do componente florestal sobre as áreas abertas, sendo este processo limitado pelas ações antropogênicas desde o Holoceno Tardio.

Atualmente, as perturbações ocasionadas pela utilização do fogo e pelo pastejo podem ser consideradas como importantes controladores da expansão da floresta sobre as áreas abertas no sul Brasil (MÜLLER et al., 2007MÜLLER, S. C. et al. Plant functional types of woody species related to fire disturbance in forest-grassland ecotones. Plant Ecology , New York, v. 189, n. 1, p. 1-14, 2007.; OVERBECK et al., 2007OVERBECK, G. E. et al. Brazil’s neglected biome: the South Brazilian Campos. Perspectives inPlant Ecology, Evolution and Systematics, Amsterdam, v. 9, n. 2, p. 101-116, 2007.). Assim, quando se deseja observar a dinâmica da expansão florestal sem a influência do manejo, o ideal é estudar áreas com exclusão de fogo e pastejo, como as áreas avaliadas do presente estudo, que foram desapropriadas para a formação do Parque Nacional de São Joaquim e, atualmente, encontram-se protegidas. Segundo Oliveira e Pillar (2004OLIVEIRA, J. M.; PILLAR, V. D. P. Vegetation dynamics on mosaics of Campos and Araucaria forest between 1974 and 1999 in Southern Brazil. Community Ecology , Budapest, v. 5, n. 2, p. 197-202, 2004.), o processo de expansão tende a ocorrer em bordas de florestas contínuas, o que também justifica a escolha das áreas estudadas, adjacentes à fragmentos florestais.

Alguns estudos sobre a expansão da vegetação florestal têm sido realizados no sul Brasil, como o de Oliveira e Pillar (2004OLIVEIRA, J. M.; PILLAR, V. D. P. Vegetation dynamics on mosaics of Campos and Araucaria forest between 1974 and 1999 in Southern Brazil. Community Ecology , Budapest, v. 5, n. 2, p. 197-202, 2004.), que observaram maior avanço em áreas onde o campo apresentava abundância de indivíduos de Baccharis uncinella DC., que podem promover o sombreamento. No entanto, ainda existem lacunas a serem preenchidas em relação ao conhecimento dos padrões de expansão da vegetação arbórea sobre áreas abertas, especialmente em relação à influência de características edáficas e topográficas sobre as taxas demográficas da vegetação. Desta forma, este trabalho teve por objetivo conhecer a dinâmica da regeneração natural de espécies arbustivo-arbóreas em áreas de vegetação aberta alto-montanas localizadas no Parque Nacional de São Joaquim, Santa Catarina, e determinar as variáveis ambientais influentes nos padrões encontrados. A hipótese é de que a expansão da vegetação arbustivo-arbórea, sintetizada pelas suas taxas demográficas, seja influenciada pelas condições edáficas e topográficas.

Material e métodos

O estudo foi realizado em três áreas de vegetação aberta, adjacentes a fragmentos florestais, no Parque Nacional de São Joaquim. O parque tem sua sede no município de Urubici, no qual foram alocadas as unidades amostrais. O clima da região, segundo classificação de Köppen, é do tipo Cfb e, de acordo com a base de dados do Worldclim (HIJMANS et al., 2005HIJMANS, R. J. et al. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology, Chichester, v. 25, n. 15, p. 1965-1978, 2005.), a precipitação média anual é 1.753 mm e a temperatura média anual é 12,7oC. A vegetação é formada por um mosaico de campos e florestas, sendo as florestas classificadas como FOM Alto-Montana, segundo o IBGE (2012IBGE. Manual técnico da vegetação brasileira. Rio de Janeiro: Fundação Instituto Brasileiro de Geografia e Estatística, 2012. 271 p.). Assim, as áreas de floresta formam fragmentos naturais na paisagem, existindo também fragmentos antrópicos, formados pela exploração madeireira e criação de gado. Os solos predominantes são classificados como Neossolos Litólicos e Cambissolos; e o relevo varia de suave a fortemente ondulado.

Com o propósito de investigar o processo de ocupação florestal nas diferentes condições campestres comuns no Sul do Brasil, a amostragem contemplou regenerantes arbustivo-arbóreos ocorrentes tanto em campos naturais, quanto em campos antrópicos. A primeira área inventariada (Área 1), localizada a 28°05’41,5”S, 49°30’14,71”W e com altitude de 1.628 m, era de floresta e foi desmatada na década de 1960, sendo utilizada para pastagem, com utilização do fogo. Ela foi desapropriada em 2007 e atualmente está protegida. A segunda área (Área 2), localizada a 28°04’46,87”S, 49°30’51,29”W e com altitude de 1.356 m, era de floresta e foi desmatada para exploração madeireira na década de 1980 e, desde essa época, foi manejada com o uso do fogo e criação de gado. Foi desapropriada em 2008 e também está protegida. A terceira área (Área 3), localizada a 28°09’49,19”S, 49°36’47,56”W e com altitude de 1.660 m, é campo (Estepe) desde os tempos remotos, com gado desde o século XIX. Foi desapropriada e protegida em 2008.

Em julho de 2014 foram instaladas três transecções, uma em cada área, de largura de 20 m por 100 m de comprimento, subdivida em subparcelas de 10 x 10 m. Cada transecção contou com 20 subparcelas e teve 2.000 m2 de área total. As transecções foram alocadas perpendiculares à borda dos fragmentos, com sua extensão cobrindo a área aberta. Foram amostrados, plaqueteados e identificados todos os indivíduos regenerantes pertencentes a espécies arbustivo-arbóreas dentro das parcelas com 20 cm ou mais de altura. O estudo não considerou indivíduos arbóreos adultos (DAP - diâmetro a altura do peito - maior ou igual a 5 cm) porque esses eram escassos nas áreas (considerando as três áreas, são menos que 10 indivíduos). As identificações botânicas foram realizadas com o auxílio de especialistas (Dr. Marcos Eduardo Guerra Sobral e Dr. Cláudio Augusto Mondin) e literatura, de forma que foi possível ter segurança de que todos os indivíduos coletados pertenciam às formas de vida consideradas no presente estudo (arbustos e árvores), e as coletas foram incorporadas ao Herbário LUSC, da Universidade do Estado de Santa Catarina. Em julho de 2015, um ano após a primeira coleta de dados, as áreas foram revisitadas para a realização do estudo da dinâmica dos regenerantes, sendo todos os indivíduos recontados para verificação da mortalidade e do recrutamento e identificação daqueles que atingiram a altura mínima de 20 cm. As taxas de mortalidade (TxMor), recrutamento (TxRec), mudança líquida em número de indivíduos (TxML) e rotatividade em número de indivíduos (TxRot) foram calculadas para cada uma das comunidades e para as populações pelas fórmulas utilizadas por Korning e Balslev (1994KORNING, J.; BALSLEV, H. Growth and mortality of trees in Amazonian tropical rain forest in Ecuador. Journal of Vegetation Science , Hoboken, v. 5, n. 1, p. 77-86, 1994.), Sheil, Jennings e Savill (2000SHEIL, D.; JENNINGS, S.; SAVILL, P. Long-term permanent plot observations of vegetation dynamics in Budongo, a Ugandan rain forest. Journal of Tropical Ecology , Cambridge, v. 16, n. 6, p. 865-882, 2000.) e Oliveira-Filho et al. (2007OLIVEIRA-FILHO, A. T. et al. Dinâmica da comunidade e populações arbóreas da borda e interior de um remanescente florestal na Serra da Mantiqueira, Minas Gerais, em um intervalo de cinco anos (1999-2004). Revista Brasileira de Botânica, São Paulo, v. 30, n. 1, p. 149-161, 2007.): TxMor = (1-((N1-NMor)/N1)1/t)×100; TxRec = (1-(1-NRec/N2)1/t)×100; TxML = ((N2/N1)1/t-1)×100; TxRot = (TxMor+TxRec)/2; onde: N1 = número inicial de indivíduos; Nmor = número de indivíduos mortos; t = tempo; Nrec = número de indivíduos recrutas; N2 = número final de indivíduos.

As espécies amostradas em ambos os anos foram classificadas em guildas de regeneração, segundo a metodologia proposta por Swaine e Whitmore (1988SWAINE, M. D.; WHITMORE, T. C. On the definition of ecological species groups in tropical rain forest. Vegetatio, New York, v. 75, n. 1/2, p. 81-86, 1988.) e adaptada por Oliveira Filho et al. (1994bOLIVEIRA-FILHO, A. T. et al. Effects of soils and topography on the distribution of tree species in a tropical riverine forest in south-eastern Brazil. Journal of Tropical Ecology, Cambridge, v. 10, n. 4, p. 483-508, 1994b.), em: pioneira (P), clímax exigente em luz (CEL) e clímax tolerante ao sombreamento (CTS), com base nas observações a campo e consulta à literatura.

Foi realizada a caracterização ambiental das áreas pela obtenção das variáveis edáficas e topográficas. As propriedades físico-químicas dos solos foram determinadas a partir da coleta de nove amostras no perfil de 0 a 20 cm, em cada uma das 60 subparcelas. As nove coletas foram homogeneizadas, retirando-se uma amostra composta por subparcela. As análises das propriedades físico-químicas dos solos foram realizadas na Universidade do Estado de Santa Catarina (UDESC), quantificando-se: teores de areia, silte e argila, pH em H2O, cálcio, magnésio, alumínio, sódio, fósforo, potássio, capacidade de troca catiônica (CTC efetiva), soma de bases e matéria orgânica. O grau de compactação do solo foi obtido por meio de testes de resistência à penetração vertical no solo, utilizando um penetrômetro de impacto digital modelo Penetrolog PLG1020 Falker, com o cone tipo 2, com dez repetições sistemáticas em cada parcela, na profundidade de 0 a 20 cm. Dessas dez leituras, foi retirada uma média por parcela. Dessa forma, em cada parcela foi obtida uma leitura de resistência à penetração a 1 cm, 2 cm, até 20 cm de profundidade. Considerou-se, na caracterização da resistência à penetração em cada parcela, a média dessas leituras de 1 cm até 20 cm, a resistência à penetração máxima nesse perfil e a resistência à penetração somente a 20 cm de profundidade. No levantamento topográfico, três variáveis topográficas foram determinadas em cada parcela: cota, desnível máximo e declividade média. A cota foi calculada pela média das cotas dos quatro vértices da parcela, o desnível máximo correspondeu à maior distância vertical entre os vértices da parcela (OLIVEIRA FILHO et al., 1994aOLIVEIRA-FILHO, A. T. et al. Differentiation of streamside and upland vegetation in an area of montane semideciduous forest in southeastern Brazil. Flora, Jena, v. 189, n. 4, p. 287-305, 1994a.) e a declividade foi calculada pela média da declividade dos quatro lados da parcela.

As taxas demográficas (mortalidade, recrutamento, rotatividade e mudança líquida), calculadas por meio das fórmulas indicadas acima, foram comparadas entre as áreas por meio do teste de Krukal-Wallis, uma vez que os dados não apresentaram distribuição normal, seguido por um teste de comparação múltipla não paramétrico. A associação entre os padrões de dinâmica (incremento, estabilidade ou redução do tamanho populacional) das espécies e as guildas de regeneração para as espécies das áreas 1 e 2 (área 3 não testada, em função da baixa abundância de regenerantes arbustivo-arbóreos) foi verificada por meio de um teste de qui-quadrado aplicado a uma tabela de contingência. A relação entre as taxas demográficas foi analisada por meio de testes de correlação. A análise da heterogeneidade ambiental entre as parcelas de cada área foi realizada por meio de uma Análise de Componentes Principais (PCA), tendo como base as variáveis ambientais amostradas nas parcelas. A influência dos gradientes ambientais sobre as taxas demográficas foi analisada por meio de modelos lineares mistos, uma vez que as taxas demográficas das parcelas de uma mesma área tendem a ser mais semelhantes entre si, devido às diferenças de histórico de uso, estágio sucessional e de altitude. Desta forma, o eixo da PCA que explicou a maior proporção da variação dos dados ambientais foi considerado como efeito fixo, e as diferentes áreas, como efeito aleatório. A área 3 não foi incluída na modelagem, uma vez que apresentou poucas parcelas com indivíduos, o que representaria uma fonte de inconsistências na análise. Foram ajustados modelos considerando diferentes estruturas de efeito aleatório (intercepto e intercepto+coeficientes angulares). Quando as variâncias não apresentaram homogeneidade ao longo dos valores ajustados e ao longo das variáveis explicativas, foram adicionados aos modelos, funções para incorporar a estrutura da heterocedasticidade existente para cada área (ZUUR et al., 2009ZUUR, A. et al. Mixed effects models and extensions in ecology with R. Amsterdam: Springer Science & Business Media, 2009. 574 p.). Ainda, quando necessário, foi aplicada a transformação logarítmica junto às variáveis dependentes, com o propósito de estabilizar a variância e a aproximação da normalidade dos resíduos. A seleção e avaliação do modelo final se deu a partir do critério de AIC e análise da normalidade e homocedasticidade dos resíduos (ZUUR et al., 2009ZUUR, A. et al. Mixed effects models and extensions in ecology with R. Amsterdam: Springer Science & Business Media, 2009. 574 p.). As análises foram realizadas por meio da linguagem de programação estatística R (R DEVELOPMENT CORE TEAM, 2015R DEVELOPMENT CORE TEAM. R: a language and environment for statistical computing. [2015]. Disponível em: <Disponível em: http://www.R-project.org > Acesso em: 08 out. 2015.
http://www.R-project.org...
).

Resultados

No período de estudo, foram observadas mudanças em número de indivíduos e tendência de estabilidade na riqueza de espécies nas áreas avaliadas, com a Área 2 apresentando maior proporção de espécies e indivíduos pioneiros do que a Área 1 (Tabela 1). A Área 3 também apresentou elevada proporção de espécies e indivíduos iniciais de sucessão, no entanto, com reduzido número de indivíduos amostrados. Apesar das três áreas apresentarem incremento na abundância em número de indivíduos, destaca-se o maior equilíbrio entre mortalidade e recrutamento na Área 2 (área de maior mortalidade), o que resultou em uma menor mudança líquida neste local (5,17%.ano-1), quando comparado com os outros (Área 1 = 21,16%.ano-1; Área 3 = 33,30%.ano-1) (Tabela 2).

Tabela 1
Número e porcentagem (%) de espécies (Ne) e indivíduos (Ni) nas guildas de regeneração, em cada uma das três áreas estudadas no Parque Nacional de São Joaquim, SC.
Table 1
Number and percentage (%) of species (Ne) and trees (Ni) in regeneration guilds, in each of studied area, at São Joaquim National Park, SC state.
Tabela 2
Dinâmica das espécies arbustivo-arbóreas regenerantes em três áreas (A1, A2 e A3) estudadas no Parque Nacional de São Joaquim, SC.
Table 2
Dynamics of regenerating tree-shrub species in three studied areas (A1, A2 and A3) at São Joaquim National Park, SC state.

Tabela 3
Teste qui-quadrado aplicado a uma tabela de contingência indicando o número de espécies observado (valores fora dos parênteses) e esperado (dentro) para as diferentes tendências de dinâmica do componente arbustivo-arbóreo regenerante (↑ = incremento populacional, ↔ = estabilidade, ↓ = decremento), considerando as diferentes guildas de regeneração, em áreas campestres alto-montanas no Parque Nacional de São Joaquim, SC.
Table 3
Chi-square test applied to a contingency table indicating the number of observed (values outside parentheses) and expected (inside) species for different dynamic trends of regenerating a treeshrub component (↑ = populational increment, ↔ = stability, ↓ = decrement), considering different regeneration guilds, in upper-montane grassland areas at São Joaquim National Park, SC state.

As cinco espécies mais abundantes foram Baccharis uncinella (P), Solanum cassioides (P), Solanum paranense (P), Myrceugenia regnelliana (CEL) e Berberis laurina (P) na Área 1; e Baccharis uncinella (P), Campovassouria cruciata (P), Senecio brasiliensis (P), Croton cf. patrum (P) e Solanum paranense (P) na Área 2. Na Área 3, as únicas espécies arbustivo-arbóreas presentes foram Siphoneugena reitzii (CEL), Baccharis uncinella (P), Gaylussacia angustifolia (P) e Baccharis nummularia (P). Dentre as espécies mais abundantes na Área 1, todas apresentaram mudança líquida positiva, com destaque para Solanum cassioides e Berberis laurina, com aumento de 37,24% e 31,71%, respectivamente, nos tamanhos populacionais. Na Área 2, Solanum paranense e Croton cf. patrum se destacaram por apresentarem, respectivamente, os maiores valores de mudança líquida positiva (88,00%.ano-1) e negativa (-32,61%.ano-1). Na Área 3, Baccharis nummularia, que era inexistente em 2014, foi a única espécie que apresentou recrutamento, com apenas dois indivíduos.

Figura 1
Relações entre as taxas demográficas (TxMor = taxa de mortalidade; TxRec = taxa de recrutamento; TxML= taxa de mudança líquida e TxRot = taxa de rotatividade), do componente arbustivo-arbóreo regenerativo em três áreas estudadas no Parque Nacional de São Joaquim, SC. A parte inferior da diagonal apresenta os gráficos de dispersão entre os valores das taxas por subparcela e a parte superior da diagonal apresenta as correlações e as significâncias das correlações entre as taxas.
Figure 1
Relationships among demographic rates (TxMor = mortality rate; TxRec = recruitment rate; TxML= net change rate and TxRot = turnover rate), of regenerating tree-shrub component in three studied area at São Joaquim National Park, SC state. The inferior part of diagonal presents the scatterplots between sub-plots rate values and superior part of diagonal shows the correlations between rates and its significance.

Figura 2
Análise de Componentes Principais com variáveis ambientais das subparcelas de três áreas estudadas no Parque Nacional de São Joaquim, SC. Argila = teor de argila; silte = teor de silte; bases = saturação por bases; Mg = teor de magnésio; pH = valor de pH; Na = teor de sódio; P = teor de fósforo; pen.med = resistência à penetração média; pen.max = resistência à penetração máxima; pen20 = resistência à penetração a 20 cm; K= teor de potássio; CTCef = CTC efetiva; areia = teor de areia; H.Al = índice H+Al; MO = teor de matéria orgânica; Al = teor de alumínio; cota = cota altitudinal; dec.med = declividade média; Sal = saturação por alumínio.
Figure 2
Principal component analysis of environmental variables from sub-plots of three studied areas at São Joaquim National Park, SC state. Argila = clay content; silte= silt content; bases = bases saturation; Mg = magnesium content; pH = pH values; Na = sodium values; P = phosphorus content; pen.med = mean penetration resistance; pen.max = maximum penetration resistance; pen20 = penetration resistance at 20 cm deep; K= potassium content; CTCef = effective cation exchange capacity; areia = sand content; H.Al = H+Al index; MO = organic matter content; Al = aluminum content; cota = elevation; dec.med = mean declivity; Sal = aluminium saturation.

Em todas as áreas, a maioria das espécies apresentou tendência de incremento ou estabilidade populacional, porém, sem relação com as guildas de regeneração (p > 0,05) (Tabela 3). Assim, todas as guildas de regeneração apresentaram as mesmas proporções de espécies com incremento, estabilidade ou decremento populacional.

Ao se analisar as relações entre as taxas demográficas, considerando as três áreas, observa-se que não existiu relação significativa entre a mortalidade e o recrutamento (r = 0,11; p = 0,47) (Figura 1), indicando que estas taxas ocorreram de forma espacialmente dissociadas. A taxa de mudança líquida foi influenciada, principalmente, pelo recrutamento (r = 0,85; p < 0,01), o que explica o incremento de abundância nas áreas; e os locais com maior rotatividade foram aqueles com maior mortalidade (r = 0,55; p < 0,01) e recrutamento (r = 0,89; p < 0,01).A PCA (Figura 2) indicou que o Eixo 1 explicou a maior proporção da variação ambiental observada nas três áreas (41,58%), e este explicou de forma significativa apenas a taxa de recrutamento (Tabela 4). As variáveis que apresentaram maiores correlações com este eixo da PCA foram o índice de H+Al (0,32), o teor de alumínio (Al) (0,32) e a soma de bases (-0,29). Assim, observa-se que, da esquerda para a direita na ordenação produzida pela PCA (PCA1), ocorre uma redução na fertilidade do solo, de parcelas com maior soma de bases e menores valores de índice H+Al e teor de Al, para parcelas com menos bases e maiores valores de Al. Quanto maiores os valores de PCA1 (direita do gráfico), menores foram as taxas de recrutamento, indicando limitação da regeneração nesses locais, como nas parcelas da Área 3, localizadas no extremo à direita.

Tabela 4
Modelos Lineares Mistos das taxas demográficas (TxMor = taxa de mortalidade; TxRec = taxa de recrutamento; TxML = taxa de mudança líquida em número de indivíduos e TxRot = taxa de rotatividade em número de indivíduos) do componente arbustivo-arbóreo regenerante em áreas campestres alto-montanas no Parque Nacional de São Joaquim, SC, em função do gradiente definido pelo Eixo 1 da PCA ambiental.
Table 4
Mixed Linear Models of demographic rates (TxMor = mortality rate; TxRec = recruitment rate; TxML= net change rate in number of individuals and TxRot = rate of turnover in number of individuals) of the shrub-tree component in upper montane grassland areas in São Joaquim National Park, SC state, in function of the gradient defined by the environmental axis of PCA.

Discussão

Os resultados demonstraram que, no período de estudo, após sete (Áreas 2 e 3) e oito (Área 1) anos da proteção integral contra distúrbios de origem antrópica, o processo de regeneração natural do componente arbustivo-arbóreo nos campos alto-montanos foi caracterizado pelo incremento na abundância de indivíduos, estabilidade na riqueza e também estabilidade na proporção da participação das espécies nas guildas de regeneração. Também houve independência espacial entre a mortalidade e o recrutamento e limitação do recrutamento em função da baixa fertilidade natural do solo e elevada disponibilidade de Al. Apesar dos padrões em comum, observou-se que o processo de expansão florestal no campo nativo (Área 3) ocorreu de forma mais lenta do que o processo de sucessão florestal secundária nos campos antrópicos (Áreas 1 e 2), sendo este aspecto parcialmente explicado pelas diferenças nas condições edáficas.

O aumento na abundância de indivíduos se deu, principalmente, em função do incremento do tamanho populacional de espécies pioneiras, o que indica que as áreas ainda se encontram em fase inicial de sucessão, uma vez que este é um padrão desta fase sucessional (SWAINE; HALL, 1983SWAINE, M. D.; HALL, J. B. Early succession on cleared forest land in Ghana. Journal of Ecology, London, v. 71, n. 2, p. 601-627, 1983.; VAN BREUGEL; BONGERS; MARTÍNEZ-RAMOS, 2007VAN BREUGEL, M.; BONGERS, F.; MARTÍNEZ-RAMOS, M. Species dynamics during early secondary forest succession: recruitment, mortality and species turnover. Biotropica, Hoboken, v. 39, n. 5, p. 610-619, 2007.). A manutenção da riqueza (ganho de apenas uma espécie nas Áreas 1 e 3) e a estabilidade da fase sucessional reforçam a ideia de que o processo de reorganização florística pode não ocorrer de forma contínua (VAN BREUGEL; BONGERS; MARTÍNEZ-RAMOS, 2007VAN BREUGEL, M.; BONGERS, F.; MARTÍNEZ-RAMOS, M. Species dynamics during early secondary forest succession: recruitment, mortality and species turnover. Biotropica, Hoboken, v. 39, n. 5, p. 610-619, 2007.) e pode ocorrer de forma mais lenta do que a reorganização estrutural.

Apesar deste padrão geral, observou-se que, para o período considerado, a dinâmica não ocorreu de forma homogênea entre as áreas, com a Área 2 apresentando maior taxa de mortalidade. Além de supostas influências relacionadas às diferenças ambientais (e.g. altitude), estas variações no padrão de dinâmica podem estar relacionadas com diferenças no grau de desenvolvimento do estágio sucessional entre as áreas, uma vez que a maior representatividade de espécies e indivíduos pioneiros na Área 2 sugere um grau de estágio sucessional mais inicial neste local. De fato, a Área 2 foi integralmente protegida mais recentemente (2008) do que a Área 1 (2007). Considerando que espécies mais iniciais apresentam usualmente um ciclo de vida mais curto, uma maior taxa de mortalidade deste grupo tem sido relatada na literatura (VAN BREUGEL; BONGERS; MARTÍNEZ-RAMOS, 2007VAN BREUGEL, M.; BONGERS, F.; MARTÍNEZ-RAMOS, M. Species dynamics during early secondary forest succession: recruitment, mortality and species turnover. Biotropica, Hoboken, v. 39, n. 5, p. 610-619, 2007.).

A independência entre a mortalidade e o recrutamento indica que, assim como em florestas em avançado em estágio de sucessão (e.g., SALAMI et al., 2014SALAMI, B. et al. Influência de variáveis ambientais na dinâmica do componente arbóreo em um fragmento de Floresta Ombrófila Mista em Lages, SC. Scientia Forestalis , Piracicaba, v. 42, n. 102, p. 197-207, 2014.), a dinâmica do componente arbóreo não é um processo espacialmente homogêneo, sugerindo a existência de fatores ecológicos, bióticos e, ou, abióticos influenciando as taxas demográficas do componente regenerativo.

A relação significativa entre a taxa de recrutamento e a fertilidade do solo (relação positiva) e também com os teores de H+Al e Al (relação negativa), sintetizada pelo Eixo 1 da PCA ambiental, sugere que a disponibilidade de nutrientes e o excesso por Al representam fatores determinantes para o estabelecimento de florestas em áreas de campos alto-montanos no sul do Brasil, atuando como filtros ambientais. De fato, Marcon et al. (2014MARCON, A. K. et al. Variação florístico-estrutural em resposta à heterogeneidade ambiental em uma floresta nebular em Ububici, Planalto Catarinense. Scientia Forestalis, Piracicaba, v. 42, n. 103, p. 439-450, 2014.), avaliando a influência de gradientes ambientais sobre a organização de uma FOM Alto-Montana no mesmo município, observaram elevada substituição florística em função da fertilidade do solo e topografia, o que reforça a ideia destes fatores na organização de comunidades arbóreas em elevada altitude. Da mesma forma, Moraes et al. (2016MORAES, D. A. et al. Edaphic filters and the functional structure of plant assemblages in grasslands in southern Brazil. Journal of Vegetation Science , Hoboken, v. 27, n. 1, p. 100-110, 2016.), avaliando a organização de comunidades de plantas em áreas campestres no Paraná a partir de atributos funcionais, concluíram que os aspectos edáficos relacionados à fertilidade natural do solo foram determinantes. Assim, na região sul do Brasil, além do clima (JESKE-PIERUSCHKA et al., 2013JESKE-PIERUSCHKA, V. et al. New insights into vegetation, climate and fire history of southern Brazil revealed by a 40,000 year environmental record from the State Park Serra do Tabuleiro. Vegetation History and Archaeobotany, New York, v. 22, n. 4, p. 299-314, 2013.), do regime de distúrbios (OVERBECK et al., 2007OVERBECK, G. E. et al. Brazil’s neglected biome: the South Brazilian Campos. Perspectives inPlant Ecology, Evolution and Systematics, Amsterdam, v. 9, n. 2, p. 101-116, 2007.), da presença de elementos nucleadores (DUARTE et al., 2006DUARTE, L. S. et al. Role of nurse plants in Araucaria Forest expansion over grassland in south Brazil. Austral Ecology, Hoboken, v. 31, n. 4, p. 520-528, 2006.; CARLUCCI; DUARTE; PILLAR, 2011CARLUCCI, M. B.; DUARTE, L. S.; PILLAR, V. D. Nurse rocks influence forest expansion over native grassland in southern Brazil. Journal of Vegetation Science, Hoboken, v. 22, n. 1, p. 111-119, 2011.; CARLUCCI et al., 2011CARLUCCI, M. B. et al. Edge expansion of Araucaria forest over southern Brazilian grasslands relies on nurse plant effect. Community Ecology, Budapest, v. 12, n. 2, p. 196-201, 2011.; MARCILIO-SILVA et al., 2015MARCILIO-SILVA, V. et al. Nurse abundance determines plant facilitation networks of subtropical forest-grassland ecotone. Austral Ecology , Hoboken, v. 40, n. 8, p. 898-908, 2015.), da distância da borda florestal (MÜLLER et al., 2012MÜLLER, S. C. et al. Woody species patterns at forest-grassland boundaries in southern Brazil. Flora, Jena, v. 207, n. 8, p. 586-598, 2012.) e da predação de sementes (IOB; VIEIRA, 2008IOB, G.; VIEIRA, E. M. Seed predation of Araucaria angustifolia (Araucariaceae) in the Brazilian Araucaria Forest: influence of deposition site and comparative role of small and ‘large’ mammals. Plant Ecology, New York, v. 198, n. 2, p. 185-196, 2008.), o presente estudo corrobora a ideia de que o estabelecimento da vegetação arbustivo-arbórea também seja influenciado pela disponibilidade de nutrientes (SOBANSKI; MARQUES, 2014SOBANSKI, N.; MARQUES, M. C. M. Effects of soil characteristics and exotic grass cover on the forest restoration of the Atlantic Forest region. Journal for Nature Conservation, Amsterdam, v. 22, n. 3, p. 217-222, 2014.) e excesso de Al no solo. Ressalta-se que a Área 3, de campo natural, na qual foram observados poucos indivíduos arbustivo-arbóreos, destacou-se por apresentar os solos de menor fertilidade e maiores teores de Al, o que contribuiu para sustentar esta suposição.

Conclusões

Conclui-se que, após a proteção integral recente das áreas campestres contra distúrbios de origem antrópica, o processo de regeneração natural do componente arbustivo-arbóreo no período de estudo ainda se encontra em fase inicial, caracterizado pela elevada participação de espécies iniciais de sucessão, pelo incremento do número de indivíduos e estabilidade na riqueza. O padrão de dinâmica variou entre as áreas em função da fase de cada área ao longo da trajetória sucessional e da natureza do campo (natural ou antrópico), embora, durante o período de avaliação, tenha sido observada estabilidade na proporção da participação das espécies nas guildas sucessionais dentro de cada área. Ainda, o recrutamento de novos indivíduos foi limitado pela disponibilidade de nutrientes no solo e elevado teor de Al. Desta forma, os resultados encontrados sugerem que o processo de expansão florestal em campos nativos é também limitado por condições edáficas, ocorrendo, assim, de forma mais lenta do que sucessão florestal em campos antrópicos.

Agradecimentos

À CAPES, pela concessão de bolsa à primeira autora, e ao CNPq, pela concessão de bolsa de produtividade para o segundo e o terceiro autores e pelo financiamento do projeto. Ao Instituto Chico Mendes de Conservação da Biodiversidade (ICMBio). À equipe do PARNA de São Joaquim, em especial ao M.Sc. Michel Omena. Aos pesquisadores Dr. Marcos Eduardo Guerra Sobral e Dr. Cláudio Augusto Mondin, pelo auxílio na identificação do material botânico.

Referências

  • BEHLING, H.; PILLAR, V. D. Late quaternary vegetation, biodiversity and fire dynamics on the southern brazilian highland and their implication for conservation and management of modern araucaria forest and grassland ecosys tems. Philosophical Transactions of the Royal Society B, London, v. 362, n. 1478, p. 243-251, 2007.
  • BEHLING, H. et al Late-Holocene fire history in a forest-grassland mosaic in southern Brasil: implications for conservation. Applied Vegetation Science, Hoboken, v. 10, n. 1, p. 81-90, 2007.
  • CARLUCCI, M. B.; DUARTE, L. S.; PILLAR, V. D. Nurse rocks influence forest expansion over native grassland in southern Brazil. Journal of Vegetation Science, Hoboken, v. 22, n. 1, p. 111-119, 2011.
  • CARLUCCI, M. B. et al Edge expansion of Araucaria forest over southern Brazilian grasslands relies on nurse plant effect. Community Ecology, Budapest, v. 12, n. 2, p. 196-201, 2011.
  • DUARTE, L. S. et al Role of nurse plants in Araucaria Forest expansion over grassland in south Brazil. Austral Ecology, Hoboken, v. 31, n. 4, p. 520-528, 2006.
  • HIJMANS, R. J. et al Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology, Chichester, v. 25, n. 15, p. 1965-1978, 2005.
  • IBGE. Manual técnico da vegetação brasileira. Rio de Janeiro: Fundação Instituto Brasileiro de Geografia e Estatística, 2012. 271 p.
  • IOB, G.; VIEIRA, E. M. Seed predation of Araucaria angustifolia (Araucariaceae) in the Brazilian Araucaria Forest: influence of deposition site and comparative role of small and ‘large’ mammals. Plant Ecology, New York, v. 198, n. 2, p. 185-196, 2008.
  • JESKE-PIERUSCHKA, V. et al New insights into vegetation, climate and fire history of southern Brazil revealed by a 40,000 year environmental record from the State Park Serra do Tabuleiro. Vegetation History and Archaeobotany, New York, v. 22, n. 4, p. 299-314, 2013.
  • KORNING, J.; BALSLEV, H. Growth and mortality of trees in Amazonian tropical rain forest in Ecuador. Journal of Vegetation Science , Hoboken, v. 5, n. 1, p. 77-86, 1994.
  • MARCILIO-SILVA, V. et al Nurse abundance determines plant facilitation networks of subtropical forest-grassland ecotone. Austral Ecology , Hoboken, v. 40, n. 8, p. 898-908, 2015.
  • MARCON, A. K. et al Variação florístico-estrutural em resposta à heterogeneidade ambiental em uma floresta nebular em Ububici, Planalto Catarinense. Scientia Forestalis, Piracicaba, v. 42, n. 103, p. 439-450, 2014.
  • MORAES, D. A. et al Edaphic filters and the functional structure of plant assemblages in grasslands in southern Brazil. Journal of Vegetation Science , Hoboken, v. 27, n. 1, p. 100-110, 2016.
  • MÜLLER, S. C. et al Plant functional types of woody species related to fire disturbance in forest-grassland ecotones. Plant Ecology , New York, v. 189, n. 1, p. 1-14, 2007.
  • MÜLLER, S. C. et al Woody species patterns at forest-grassland boundaries in southern Brazil. Flora, Jena, v. 207, n. 8, p. 586-598, 2012.
  • OLIVEIRA, J. M.; PILLAR, V. D. P. Vegetation dynamics on mosaics of Campos and Araucaria forest between 1974 and 1999 in Southern Brazil. Community Ecology , Budapest, v. 5, n. 2, p. 197-202, 2004.
  • OLIVEIRA-FILHO, A. T. et al Differentiation of streamside and upland vegetation in an area of montane semideciduous forest in southeastern Brazil. Flora, Jena, v. 189, n. 4, p. 287-305, 1994a.
  • OLIVEIRA-FILHO, A. T. et al Dinâmica da comunidade e populações arbóreas da borda e interior de um remanescente florestal na Serra da Mantiqueira, Minas Gerais, em um intervalo de cinco anos (1999-2004). Revista Brasileira de Botânica, São Paulo, v. 30, n. 1, p. 149-161, 2007.
  • OLIVEIRA-FILHO, A. T. et al Effects of soils and topography on the distribution of tree species in a tropical riverine forest in south-eastern Brazil. Journal of Tropical Ecology, Cambridge, v. 10, n. 4, p. 483-508, 1994b.
  • OVERBECK, G. E. et al Brazil’s neglected biome: the South Brazilian Campos. Perspectives inPlant Ecology, Evolution and Systematics, Amsterdam, v. 9, n. 2, p. 101-116, 2007.
  • R DEVELOPMENT CORE TEAM. R: a language and environment for statistical computing. [2015]. Disponível em: <Disponível em: http://www.R-project.org > Acesso em: 08 out. 2015.
    » http://www.R-project.org
  • SALAMI, B. et al Influência de variáveis ambientais na dinâmica do componente arbóreo em um fragmento de Floresta Ombrófila Mista em Lages, SC. Scientia Forestalis , Piracicaba, v. 42, n. 102, p. 197-207, 2014.
  • SHEIL, D.; JENNINGS, S.; SAVILL, P. Long-term permanent plot observations of vegetation dynamics in Budongo, a Ugandan rain forest. Journal of Tropical Ecology , Cambridge, v. 16, n. 6, p. 865-882, 2000.
  • SILVA, L. C. R.; ANAND, M. Mechanisms of Araucaria (Atlantic) forest expansion into southern Brazilian grasslands. Ecosystems, Nova York, v. 14, n. 8, p. 1354-1371, 2011.
  • SOBANSKI, N.; MARQUES, M. C. M. Effects of soil characteristics and exotic grass cover on the forest restoration of the Atlantic Forest region. Journal for Nature Conservation, Amsterdam, v. 22, n. 3, p. 217-222, 2014.
  • SWAINE, M. D.; HALL, J. B. Early succession on cleared forest land in Ghana. Journal of Ecology, London, v. 71, n. 2, p. 601-627, 1983.
  • SWAINE, M. D.; WHITMORE, T. C. On the definition of ecological species groups in tropical rain forest. Vegetatio, New York, v. 75, n. 1/2, p. 81-86, 1988.
  • VAN BREUGEL, M.; BONGERS, F.; MARTÍNEZ-RAMOS, M. Species dynamics during early secondary forest succession: recruitment, mortality and species turnover. Biotropica, Hoboken, v. 39, n. 5, p. 610-619, 2007.
  • ZUUR, A. et al Mixed effects models and extensions in ecology with R. Amsterdam: Springer Science & Business Media, 2009. 574 p.

Datas de Publicação

  • Publicação nesta coleção
    30 Maio 2019
  • Data do Fascículo
    Jan-Mar 2019

Histórico

  • Recebido
    28 Fev 2017
  • Aceito
    01 Ago 2017
  • Publicado
    29 Mar 2019
Universidade Federal de Santa Maria Av. Roraima, 1.000, 97105-900 Santa Maria RS Brasil, Tel. : (55 55)3220-8444 r.37, Fax: (55 55)3220-8444 r.22 - Santa Maria - RS - Brazil
E-mail: cienciaflorestal@ufsm.br