View
219
Download
0
Category
Preview:
Citation preview
PROGRAMA DE PÓS-GRADUAÇÃO EM BIOLOGIA ANIMAL
FABRÍCIO BARRETO TERESA
Diversidade funcional de comunidades de peixes de riachos
Tese apresentada para obtenção do título de Doutor em Biologia Animal junto ao Programa de Pós-Graduação em Biologia Animal do Instituto de Biociências, Letras e Ciências Exatas da Universidade Estadual Paulista “Júlio de Mesquita Filho”, Campus de São José do Rio Preto. Orientador: Profª. Drª. Lilian Casatti
São José do Rio Preto - SP
2012
Teresa, Fabrício Barreto. Diversidade funcional de comunidades de peixes de riachos / Fabrício Barreto Teresa. - São José do Rio Preto : [s.n.], 2012. 101 f. : il. ; 30 cm. Orientador: Lilian Casatti Tese (doutorado) – Universidade Estadual Paulista, Instituto de Biociências, Letras e Ciências Exatas 1. Ecologia animal. 2. Ictiofauna. 3. Peixe de riacho. 4. Peixe – Estrutura de comunidades. 5. Diversidade biológica. I. Casatti, Lilian. II. Universidade Estadual Paulista, Instituto de Biociências, Letras e Ciências Exatas. III. Título.
CDU – 591.5
Ficha catalográfica elaborada pela Biblioteca do IBILCE
Campus de São José do Rio Preto - UNESP
Este trabalho foi realizado no Laboratório de Ictiologia, do Departamento de
Zoologia e Botânica do Instituto de Biociências, Letras e Ciências Exatas, UNESP de
São José do Rio Preto, com auxílio financeiro da FAPESP na forma de bolsa de
doutorado
(processo 2008/03583-9)
Fabrício Barreto Teresa
Diversidade funcional de comunidades de peixes de riachos Tese apresentada para obtenção do título de Doutor em Biologia Animal junto ao Programa de Pós-Graduação em Biologia Animal do Instituto de Biociências, Letras e Ciências Exatas da Universidade Estadual Paulista “Júlio de Mesquita Filho”, Campus de São José do Rio Preto.
Banca Examinadora
Profª. Drª. Lilian Casatti UNESP – São José do Rio Preto Orientadora Prof. Dr. Maurício Cetra Universidade Federal de São Carlos (UFSCar) Prof. Dr. Marcus Vinícius Cianciaruso Universidade Federal de Goiás (UFG) Profª. Drª. Katharina Eichbaum Esteves Instituto de Pesca Profª. Drª. Denise de Cerqueira e Rossa Feres UNESP – São José do Rio Preto
São José do Rio Preto
08/fevereiro/ 2012
AGRADECIMENTOS .
A defesa de uma tese marca o término de um ciclo que começou muito tempo antes do
ingresso na pós-graduação. Por isso, sou grato a todas as pessoas que contribuíram de
alguma forma para a minha formação pessoal e profissional.
A minha orientadora, Profa. Dra. Lilian Casatti pela confiança, amizade e por tudo que
me ensinou. Durante minha carreira espero um dia conseguir chegar ao menos próximo
da sua dedicação, seriedade, dinamismo e competência.
A Profa. Dra. Eliane Gonçalves-de-Freitas, que além de acompanhar o
desenvolvimento do meu projeto de doutorado, orientou-me durante a graduação e o
mestrado, ensinando-me a prática científica e despertando-me para o gosto pelo ensino.
A Profa Dra Maria Stela M.C. Noll e ao Dr. Mateus F. Ferrareze pelas críticas e
sugestões proferidas durante a qualificação.
A todos os professores da graduação e pós-graduação, em especial a Profa Dra Denise
de Cerqueira e Rossa Feres e ao Prof Dr. Francisco Langeani pela amizade, incentivo
e ensinamentos ao longo de toda a minha vida acadêmica.
Ao Prof Dr. Marcus Vinícius Cianciaruso por gentilmente ter me recebido em seu
laboratório na UFG e pela solicitude no ensino de metodologias de análise de
diversidade funcional.
Aos novos e antigos amigos do laboratório de Ictiologia André, Breno, Diego, Flávio,
Jaquelini, Mariela, Márcio, Paulo, Mateus, Lucas, em especial Ângelo, Camilo,
Cristiane, Gustavo, Fernanda, Fernando, Mônica, Renato e Roselene pela amizade,
momentos agradáveis que compartilhamos e pela ajuda em diferentes etapas da tese.
A André L.S. Navarro, Angélica Amaya, Angelo R. Manzoti, Bruno Luiz R. Silva,
Camilo Roa Fuentes, Cristiane P. Ferreira, Danilo G.H. Silva, Dayane Boracini, Diogo
B. Provete, Edis Belini, Felipe Nuvuoli, Fernanda Martins, Fernando R. Silva, Getúlio
M. Tanaka, Gustavo Baviera, José Luis V. Junior, Larissa, Leonardo Gedraite, Luciana
S. Ondei, Manuel A.M. Cabajal, Mariana Garcia, Mônica C. Bastos, Paulinho
(motorista), Renato M. Romero, Roselene S. C. Ferreira, Samuel V. Vieira e Thiago
Galego pela ajuda nos trabalhos de campo.
Ao Bruno Luiz Rodrigues da Silva, amigo de longa data pela inestimável ajuda nas
diferentes etapas deste trabalho, desde as viagens de reconhecimento até a última coleta,
sempre disposto e bem-humorado. Também sou grato a sua família pela amizade e por
acolher-me em sua casa no período final do doutorado
Ao Renato de Mei Romero, especialista em proporcionar apoio moral, pela amizade,
companheirismo e pela ajuda em diversas etapas da tese, assim como em outros
trabalhos.
Aos funcionários do IBILCE pela convivência e amizade nos últimos 12 anos, em
especial Silvia E. Kazama, Vitor B. Thomazine e Rosemar R. C. Brena da seção de
Pós-Graduação pela solicitude e paciência com os pedidos fora de prazo.
A Jane Dilvana Lima, pela amizade, incentivo, paciência e disponibilidade no
atendimento das minhas frequentes solicitações para revisão dos capítulos dessa tese;
ao Dr. Robert Vadas-Jr e Dr. Fábio Cop Ferreira pela avaliação crítica do capítulo I;
ao Diogo B. Provete e Thiago Gonçalves de Souza pela leitura crítica do capítulo III.
A todos os amigos futebolistas do IBILCE pelos agradáveis momentos compartilhados
nas “peladas” semanais que contribuíram para manter minha salubridade mental e
física.
A toda minha família em especial aos meus pais Maria Luiza Ribeiro Barreto e José
Teresa pelo amor, carinho e valores que alicerçam minha vida. Às minhas irmãs
Amanda e Débora, sobrinhas Anelise e Heloisa e cunhado Alexandre pelo carinho e
companheirismo, tornando adoráveis os momentos que passamos juntos.
À minha esposa Luciana de Souza Ondei pelo amor, paciência, companheirismo,
compreensão e pelos momentos felizes que compartilhamos. Como dizem nossos
amigos Celso e Lídia de Goiânia, sou um dependente fisiológico e comportamental de
RESUMO .
Neste estudo avaliaram-se os efeitos do desmatamento das zonas ripárias sobre a
composição e diversidade funcional de comunidades de peixes de riachos. Para isso, foram
elaborados três capítulos a partir de amostragens da ictiofauna em 126 segmentos de cinco
metros de extensão em seis riachos, sendo três florestados e três desmatados, localizados
na bacia do rio São José dos Dourados, sistema do Alto rio Paraná, Noroeste do Estado de
São Paulo, Brasil. No CAPÍTULO I duas características que descrevem a distribuição
longitudinal das espécies, a preferência por fluxo e por profundidade, foram avaliadas. As
espécies apresentaram diferentes padrões de preferência, mas as respostas ao fluxo foram
mais consistentes entre riachos florestados e desmatados do que em relação à
profundidade. No CAPÍTULO II os padrões de diversidade taxonômica e funcional foram
comparados entre os riachos florestados e desmatados, considerando a heterogeneidade
estrutural em uma escala de meso-hábitat e no CAPÍTULO III foi avaliado se as alterações
provocadas pelo desmatamento representam filtros ambientais que modulam a composição
funcional das comunidades. Os padrões de diversidade e composição funcional diferiram
em função das alterações ambientais nos riachos associadas ao desmatamento. As
características funcionais associadas com o uso de hábitat, ecologia trófica, tolerância à
hipóxia e tamanho influenciaram a distribuição das espécies dentro e entre riachos. Os
resultados sugerem que as mudanças na disponibilidade de hábitat e recursos alimentares,
assim como nas condições físico-químicas representam filtros ambientais importantes na
organização das comunidades. Ao relacionar o desmatamento com alterações em outros
componentes da diversidade que não somente taxonômico, os resultados obtidos
possibilitam uma avaliação mais completa sobre os impactos biológicos desse tipo de
alteração.
Palavras-chave: ictiofauna, desmatamento, zona ripária, filtros ambientais
ABSTRACT .
The effects of riparian deforestation on the functional composition and diversity of stream
fish communities were evaluated. This study was divided into three chapters whose data
were obtained by sampling fish fauna in 126 five meters long mesohabitats in forested (n =
3) and deforested (n = 3) streams, located in the São José dos Dourados river Basin, upper
Paraná river system, northwestern São Paulo State, Brazil. In CHAPTER I two traits that
describe the longitudinal distribution of the species, the preference for velocity and depth
were evaluated. The species showed different patterns of preference, but their responses
were more consistent between forested and deforested streams for velocity than for depth.
In CHAPTER II the patterns of taxonomic and functional diversity were compared
between forested and deforested streams in a mesohabitat scale and in CHAPTER III it
was evaluated if the environmental modifications caused by deforestation represent
environmental filters that modulate the functional composition of communities. Patterns of
functional diversity and composition differed according to the environmental changes in
streams related to deforestation. The distribution of species within and between streams
was influenced by traits such as habitat use, trophic ecology, tolerance to hypoxia and size.
The results suggest that changes in the habitat and food availability and physicochemical
conditions represent environmental filters acting in the organization of communities. By
linking deforestation to changes in other components of diversity rather than only
taxonomic, the results obtained here provide a more complete knowledge about the
biological impacts of the loss of riparian forest.
Keywords: ichthyofauna, deforestation, riparian zone, environmental filters
SUMÁRIO
INTRODUÇÃO GERAL................................................................................................ 14
CAPÍTULO ...................................................................................................................... 18
Abstract........................................................................................................................... 19
Introduction..................................................................................................................... 20
Methods.......................................................................................................................... 22
Study area.................................................................................................................. 22
Fish sampling and environmental descriptors .......................................................... 22
Development of habitat suitability criteria and test for transferability….………… 23
Results............................................................................................................................ 24
Discussion....................................................................................................................... 27
References....................................................................................................................... 31
CAPÍTULO 2................................................................................................................... 35
Abstract........................................................................................................................... 36
Introduction..................................................................................................................... 37
Material and Methods..................................................................................................... 39
Study area and general sampling design………………………………………….. 39
Data sampling............................................................................................................ 40
Trait and trait categories............................................................................................ 42
Functional diversity metrics...................................................................................... 44
Data analysis........................................................................................................... 45
Results......................................................................................................................... 45
Discussion.................................................................................................................... 48
References................................................................................................................... 52
Supporting information............................................................................................... 59
CAPÍTULO 3................................................................................................................. 64
Resumo........................................................................................................................ 65
Introdução.................................................................................................................... 66
Material e Métodos...................................................................................................... 68
Área de estudo........................................................................................................ 68
Amostragem da ictiofauna...................................................................................... 69
Caracterização das unidades amostrais…………………………………............ 70
Características funcionais……………………………………………………… 71
Análise dos dados................................................................................................... 73
Resultados.................................................................................................................... 75
Discussion.................................................................................................................... 81
Referências Bibliográficas........................................................................................... 85
Apêndice...................................................................................................................... 90
DISCUSSÃO GERAL................................................................................................... 93
CONCLUSÕES.............................................................................................................. 96
REFERÊNCIAS GERAIS............................................................................................ 97
14
INTRODUÇÃO GERAL .
A conservação de ecossistemas aquáticos de água doce é um desafio limitado
pela falta de conhecimento sobre a sua estrutura e funcionamento, bem como sobre as
formas pelas quais esses sistemas são afetados pelas interferências antrópicas
(BARLETTA et al., 2010). Dentre as atividades humanas que exercem impactos
negativos sobre os rios e riachos destacam os represamentos, substituição de vegetação
nativa por culturas agrícolas ou pecuária, introdução de espécies exóticas e poluição
(ALLAN; FLECKER, 1993). A supressão da vegetação nativa dentro das bacias de
drenagem e, especialmente nas margens dos rios (zonas ripárias), são uma das
alterações cujos efeitos ainda não são completamente compreendidos, principalmente
em ecossistemas tropicais (BOJSEN; BARRIGA, 2002; LORION; KENNEDY, 2009).
No Brasil, a importância do desmatamento como fonte de impactos aos rios e
riachos varia regionalmente (BARLETTA et al., 2010). Em regiões populosas e
industrializadas esse tipo de alteração é uma das principais ameaças. Um exemplo
marcante é a região Noroeste do Estado de São Paulo, uma das mais desmatadas do
país, apresentando apenas 4% de vegetação nativa remanescente (SMA/IF, 2005). A
perda da vegetação nativa atinge 75% das zonas ripárias dos rios e riachos desta região
(SILVA et al., 2007) que por sua vez apresentam baixa integridade física e biológica
(CASATTI et al., 2006, 2009a).
A perda da vegetação ripária pode interferir na biota aquática diretamente,
através das modificações nas condições abióticas como temperatura e oxigênio
dissolvido, afetando o desempenho dos organismos em diferentes etapas do seu ciclo de
vida (LOWE-MCCONNELL, 1999), assim como no fornecimento de material alóctone
(galhos, frutos, sementes, folhas e troncos) para dentro do canal, que servem como
alimento, sítio reprodutivo e refúgio para diversas espécies (FAUSCH; NORTHCOTE,
1992; SCHNEIDER; WINEMILLER, 2008). Os efeitos indiretos incluem mudanças na
disponibilidade de alimento (e.g., algas e detrito) em função do aumento da
produtividade primária devido à perda do sombreamento (BOJSEN; BARRIGA, 2002;
LORION; KENNEDY, 2009) e a redução da retenção de sedimentos e outras
substâncias provenientes do ambiente terrestre adjacente (GREGORY et al., 1991),
influenciando a qualidade da água e a composição do substrato (RABENI; SMALE,
1995). A perda da vegetação ripária também proporciona a proliferação de gramíneas
15
nos barrancos (PUSEY; ARTHINGTON, 2003, GROWNS et al., 2003, ROCHA et al.,
2009), que por sua vez constituem micro-hábitats favoráveis para a ocorrência de
algumas espécies oportunistas (ROCHA et al., 2009).
Os efeitos do desmatamento das zonas ripárias sobre os organismos e
populações traduzem-se também em nível de comunidade e são representadas por
mudanças na composição de espécies e nos padrões de diversidade (e.g., Pusey &
Arthington, 2003; CASATTI et al., 2009b). Bojsen; Barriga (2002) evidenciaram
substituição de espécies e aumento da abundância de peixes em riachos submetidos a
desmatamento das suas zonas ripárias na Amazônia equatoriana. Além desses efeitos,
Lorion; Kennedy (2009) observaram também maior riqueza de espécies em riachos
desprovidos de matas ripárias na Costa Rica, resultados corroborados por Teresa;
Casatti (2010) no sudeste do Brasil.
Além dos impactos sobre os padrões de composição e diversidade dentro de
uma perspectiva taxonômica, existem evidências de que o declínio e/ou favorecimento
das espécies em riachos submetidos ao desmatamento não é aleatório e dependem de
algumas características dos organismos que influenciam o seu desempenho nesses
ecossistemas (características funcionais). Por exemplo, Casatti et al. (2009b)
compararam a fauna de peixes entre riachos mais afetados pela Interferências antrópica,
com ausência de vegetação ripária arbórea e riachos submetidos a menor pressão
antrópica, com zona ripária arborizada. Esses dois grupos de riachos abrigaram
comunidades muito distintas, com as espécies tolerantes e generalistas dominando os
riachos do primeiro grupo e as espécies reofílicas sendo mais associadas ao segundo
grupo (CASATTI et al., 2009b). De forma congruente, Teresa; Casatti (2010)
evidenciaram que o aumento da riqueza de espécies nos riachos desmatados é
representado pelo aumento no número de espécies tolerantes à degradação física do
hábitat.
Se as características funcionais são importantes para a definição da persistência
ou declínio das espécies em ambientes antropizados, como sugerido anteriormente, a
utilização de métricas que incorporem esse componente biológico (i.e., diversidade
funcional) pode fornecer informações úteis para a compreensão da estrutura e
organização das comunidades ao longo de gradientes de degradação antrópica. A
diversidade funcional pode ser definida como a variabilidade de características
biológicas relacionadas com a função que as espécies desempenham no ecossistema
16
(DIAZ; CABIDO, 2001). A importância de se considerar esse componente da
biodiversidade é reforçada por estudos que têm revelado que a diversidade funcional
tem maior poder preditivo sobre os processos ecossistêmicos do que a diversidade
taxonômica (TILMAN et al., 1997; DIAZ; CABIDO, 2001; MOUILLOT et al., 2011).
A explicação para isso reside no fato de que a abordagem taxonômica assume que as
espécies são diferentes sem a ponderação sobre a equivalência ou complementaridade
do papel que elas desempenham no ambiente (DIAZ; CABIDO, 2001). Nesse contexto,
a avaliação dos padrões de diversidade funcional das comunidades em ambientes
submetidos às interferências antrópicas pode contribuir também para a predição de
mudanças em processos ecológicos. No caso dos riachos submetidos ao desmatamento,
os padrões de diversidade funcional poderiam indicar o grau de modificação de
processos (e.g., produtividade, sedimentação, dinâmica trófica) que são
conhecidamente influenciados por esse tipo de alteração (MURPHY; HALL, 1981;
SWEENEY et al., 2004).
Alguns estudos têm utilizados métricas de diversidade funcional dentro de uma
perspectiva aplicada com o intuito de verificar os efeitos biológicos das interferências
antrópicas (e.g., ERNST et al., 2006; EROS et al., 2009; FLYNN et al., 2009;
VILLÉGER et al., 2010; BARRAGÁN et al., 2011). Por exemplo, Ernst et al. (2006)
observaram redução da diversidade funcional de comunidades de anuros em fragmentos
florestais degradados por desmatamento, embora não tenham diagnosticado mudanças
na diversidade taxonômica. Nesse caso, a alteração ambiental resultou na substituição
de espécies funcionalmente distintas por espécies redundantes. Dentro dessa abordagem
aplicada, é evidente que a incorporação de métricas de diversidade funcional fornece
subsídios para uma avaliação mais completa sobre os impactos das alterações humanas
sobre a biodiversidade. Com relação aos riachos, embora os efeitos das alterações
antrópicas sobre a estrutura e composição taxonômica das comunidades sejam
relativamente bem compreendidos (BOJSEN; BARRIGA, 2002; PUSEY;
ARTHINGTON, 2003; LORION; KENNEDY, 2009), o mesmo não acontece com
relação à perspectiva funcional.
Assim, o objetivo geral deste estudo foi avaliar os efeitos do desmatamento das
zonas ripárias sobre a composição e diversidade funcional de comunidades de peixes de
riachos. Baseado no conhecimento de que o desmatamento influencia as comunidades
de peixes (BOJSEN; BARRIGA, 2002; LORION; KENNEDY, 2009) e nas evidências
17
de que tais Efeitos são dependentes das características funcionais exibidas pelas
espécies (CASATTI et al., 2009b; TERESA; CASATTI, 2010), hipotetizamos que as
modificações ambientais causadas pelo desmatamento nos riachos refletem em
mudanças na composição e diversidade funcional das comunidades. Para o
cumprimento do objetivo geral, este estudo foi dividido em três capítulos que
correspondem a manuscritos, cujos dados foram obtidos a partir das amostragens da
ictiofauna em seis riachos, sendo três florestados e três desmatados, localizados na
bacia do rio São José dos Dourados, sistema do Alto rio Paraná, Noroeste do Estado de
São Paulo, Brasil.
Uma das dificuldades para a utilização da abordagem funcional para alguns
táxons é o conhecimento limitado sobre as características funcionais exibida pelas
espécies (PETCHEY; GASTON, 2006), como é o caso dos peixes de riachos. Assim,
inicialmente (Capítulo I) foi avaliada a preferência por fluxo e profundidade, duas
características associadas com a distribuição das espécies ao longo do gradiente
longitudinal. Posteriormente (Capítulo II), os padrões de diversidade taxonômica e
funcional foram comparados entre os riachos florestados e desmatados, considerando a
heterogeneidade estrutural em escala de meso-hábitat. Finalmente (Capítulo III), foi
avaliado se as alterações provocadas pelo desmatamento representam filtros ambientais
e para isto foram comparados padrões de diferenciação funcional entre as comunidades
dos riachos florestados e desmatados com padrões obtidos aleatoriamente.
18
Manuscrito submetido para Ecohydrology
- CAPÍTULO I -
Development of habitat suitability criteria for some Neotropical stream
fishes and an assessment of their transferability to streams with
different conservation status
19
ABSTRACT
We assessed the preference of 10 fish species for depth and velocity conditions in
forested streams from southeastern Brazil using habitat suitability criteria (HSC curves).
We also tested whether preference patterns observed in forested streams can be
transferred to deforested streams. We used data from fish sampled in 62 five-meter
quadrats in three forested streams to construct preference curves. Astyanax altiparanae,
A. fasciatus, Knodus moenkhausii, and Piabina argentea showed a preference for deep
slow habitats, whereas Aspidoras fuscoguttatus, Characidium zebra, Cetopsorhamdia
iheringi, Pseudopimelodus pulcher, and Hypostomus nigromaculatus showed an
opposite pattern: preference for shallow fast habitats. Hypostomus ancistroides showed
a multimodal pattern of preference for depth and velocity. To evaluate whether patterns
observed in forested streams may be transferred to deforested streams, we sampled 64
five-meters quadrats in three deforested streams using the same methodology. The
preference for velocity was more consistent than for depth, as success in the
transferability criterion was 86% and 29% of species, respectively. This indicates that
velocity is a good predictor of species abundance in streams, regardless of their
condition.
Keywords: habitat preference; velocity; depth; riparian forest; fish assemblages
20
INTRODUCTION
The recognition that species are distributed non-randomly in habitats,
preferentially occupying suitable sites for feeding, reproduction, and survival (Grinnell,
1917; Hutchinson, 1957), constitutes the basis for developing predictive models of
species distribution along environmental gradients (e.g., Ahmadi-Nedushan et al.,
2006). This approach has a broad application for natural resources management, as
since it allows prediction of changes in biological communities from changes in
environmental conditions caused by human activities. However, the predictive power of
these models depends on knowledge of the response of species in relation to the
gradient of environmental conditions. This information can be obtained by developing
habitat suitability criteria (HSC; Bovee, 1982, 1986), a component for conducting
physical habitat simulation (PHABSIM), a tool used to assess the consequences of
changes in river flows on habitat availability (Bovee, 1982).
The lack of knowledge about traits exhibited by species has limited the use of
the functional approach in studies focusing in diversity patterns. In this way, the
definition of HSC also may be useful, since the preference of the species for habitat
features may be used as traits in functional diversity studies. The HSC can be assessed
by curves defining the degree of preference by a given species in relation to habitat
variables alone, such as depth, velocity, substrate and cover. Efforts to develop HSC for
multiple species have increased, owing to the degradation of water resources and
growing concern over biodiversity loss (e.g., Lamouroux et al., 1999; Vadas & Orth,
2001; Strakosh et al., 2003). However, this approach has not been applied in
Neotropical environments, limiting the use of predictive models and the ability to
quantitatively predict the biological impacts from changes in hydraulic variables of
aquatic environments. This is especially alarming because changes in the natural flow
regime are among the major anthropogenic threats for fish conservation in Neotropical
lotic environments (Barletta et al., 2010).
Several factors, such as environmental heterogeneity, food availability,
predation, competition, and habitat availability influence the patterns of species
occupancy in riverine patches (Power, 1984; Thomas & Bovee, 1993; Jackson et al.,
2001; Leal et al., 2011). The differential effects of these factors may limit the
applicability of HSC between different environments (Groshens & Orth, 1993; Freeman
et al., 1997; Vadas & Orth, 2001). In these cases, the development of specific criteria
21
for each type of stream are recommended (Bovee, 1986), which brings high costs and
operational difficulties. However, other studies have found consistency in preference
patterns across streams (Thomas & Bovee, 1993; Strakosh et al., 2003). This highlights
the importance of testing HSC validity to the site of interest (transferability) before
developing new criteria.
Our present goal is to evaluate the preference of different fish species for depth
and velocity conditions in forested streams from southeastern Brazil, based on the
development of HSC curves. In our study region, the clearing of native vegetation
affects about 75% of riparian zones (Silva et al., 2007), and most streams have some
degradation of their physical integrity that is directly or indirectly related to
deforestation (Casatti et al., 2006). Here, we also evaluate the consistency of habitat
preference patterns of species in forested streams and test HSC transferability to
deforested streams.
22
METHODS
Study area
All streams studied (second and third order) are located in the drainage of the
São José dos Dourados, Upper Paraná River basin, northwestern region of the state of
São Paulo, Brazil (Figure 1). The climate of the region is influenced by equatorial and
tropical masses, resulting in a tropical climate with dry and rainy periods, with higher
rainfall and higher temperatures from December to February and lower rainfall and
milder temperatures between June and August (IPT, 1999). The forested streams are
located in forest fragments, representing what is best preserved in the region, whereas
the deforested streams are located away from forest fragments where the riparian zone
is dominated by pasture.
Figure 1. Location of the study sites in the São José dos Dourados basin, São Paulo
State, Brazil (A and B). Forested sites (dark circles) are located within the greatest forest
fragments and deforested (gray circles) streams lack arboreal vegetation in their riparian
zone and are located in areas dominated by pasture.
Fish sampling and environmental descriptors
We collectively sampled 62 five-meter quadrats in three forested streams to
identify the preference of fish species for depth and velocity and thus define habitat
suitability criteria. In each stream we selected enough 19-24 quadrats to maximize
structural variation between quadrats to better represent hydrological variability while
23
reducing variation within each quadrat. The fish samplings were carried out during the
dry season (September to November 2009 and April to July 2010) using electrofishing
preceded by the up- and downstream isolation of each quadrat. The electrofishing
sampling consisted of successive passes along the segment. The sample ended when no
new individual was captured after one passage.
After the fish sampling, we characterized the segments in relation to depth and
velocity. To accomplish this, we divided each quadrat into four equidistant, transverse
transects, where were done measurements at five equidistant points from one margin to
another. We used a mechanical flowmeter model 2030 (General Oceanics®) to measure
water velocity at ~ 0.6 of total depth. Data from transects were lumped together to
obtain the average values for the whole quadrat.
Development of habitat suitability criteria and test for transferability
Only species that occurred in at least 10 samples in forested streams were
included in the habitat suitability criteria for depth and velocity. We adopted this
procedure because infrequent species may have been under-sampled, resulting in higher
estimated errors of preference. We classified the depth gradient into five classes (D1 = Ø
to 0.2 m; D2 = 0.21 to 0.30 m; D3 = 0.31 to 0.40 m; D4 = 0.41 to 0.50 m; D5 > 0.51 m)
and velocity into four classes (V1 = 0-0.20 m/s; V2 = 0.21 to 0.40 m/s; V3 = 0.41-0.60
m/s; V4 > 0.61 m/s). The limits for the classes were defined to represent discrete units of
the hydrological gradient while ensuring a comparable number of samples in each class.
We built HSC curves (type III, sensu Bovee, 1986) from preference values that were
calculated the abundance ratio of species in each class of depth and velocity, divided by
the proportionate number of habitat in the class in relation to all habitats sampled
(Freeman et al., 1997). With this procedure, we considered both the relative abundance
of species in the habitats and habitat availability. Posteriorly, we standardized the
preference values to vary from 0 to 1, thereby obtaining the suitability index (SI). For
the class with the highest preference value, we assigned an index of 1 and proportionally
smaller values to other classes.
We tested if the HSC developed for forested streams could be transferred to
deforested streams. To do this, we sampled fish and environmental descriptors in 64
five-meter quadrats in three deforested streams (between 16 and 24 segments in each
stream) during the same period and using the same methodology applied to forested
streams. We considered optimal those SI values for depth and velocity conditions
24
greater than 0.7, according to the criteria developed for forested streams. We choose 0.7
as a reference value because it included most of the raw values of habitat preference
higher than 1, i.e., whose relative abundance is higher than the proportional availability
of habitats.
To test the transferability of HSC from forested to deforested streams, we
considered only species that occurred in at least 10 quadrats in the deforested streams.
We classified the quadrats according to the optimum/non-optimum threshold defined
from the HSC curves from forested streams. For the criterion to be considered
transferable, the abundance of a species should be higher in “optimal” samples of
deforested streams, as defined for forested streams. Thus, we conducted a comparison
of the mean species abundance between optimum and non-optimum quadrats using a
one-tailed t-test.
RESULTS
We recorded 31 species in forested streams, of which 10 occurred in at least 10
samples (Table I). The evaluation of the HSC curves for these species for depth and
velocity demonstrated three main patterns (Figure 2). Astyanax altiparanae, A.
fasciatus, Knodus moenkhausii, and Piabina argentea had high SI values in slow-
flowing and deep segments. Another four species showed an opposite pattern, with
higher SI values associated with shallow and fast-flowing habitats. This was the case of
Characidium zebra, Cetopsorhamdia iheringi, Pseudopimelodus pulcher, and
Hypostomus nigromaculatus. Aspidoras fuscoguttatus and Hypostomus ancistroides
showed distinct patterns from the above species and had the highest SI values in
multiple classes of depth and in habitats with high and low water velocity, respectively
(Figure 2).
25
Table I. Frequency of occurrence and number of individuals (in parentheses) of species
occurring in at least 10 segments in forested streams (R1 to R3).
Family Species Code R1 R2 R3
Characidae Astyanax altiparanae Astalt 15 (97) 9 (26) 8 (30)
Astyanax fasciatus Astfas 6 (20) 14 (36) 17 (192)
Knodus moenkhausii Knomoe 8 (17) 10 (30) 8 (13)
Piabina argentea Piaarg 13 (35) 11 (34) 12 (42)
Crenuchidae Characidium zebra Chazeb 16 (31) 20 (66) 19 (89)
Callichthydae Aspidoras fuscoguttatus Aspfus 0 (0) 6 (11) 10 (23)
Heptapteridae Cetopsorhamdia iheringi Cetihe 5 (8) 4 (9) 6 (25)
Pseudopimelodus pulcher Psepul 11 (42) 9 (38) 0 (0)
Loricariidae Hypostomus ancistroides Hypanc 6 (10) 10 (23) 16 (38)
Hypostomus nigromaculatus Hypnig 18(69) 14 (47) 11 (87)
26
Figure 2. Habitat suitability criteria for depth and velocity of 10 species (codes in Table
I) in five classes of depth (D1 = Ø to 0.2 m; D2 = 0.21 to 0.30 m; D3 = 0.31 to 0.40 m;
D4 = 0.41 to 0.50 m; D5 > 0.51 m) and in four classes of velocity (V1 = > 0.0 to 0.2
m/s; V2 = 0.21 to 0.40 m/s; V3 = 0.41 to 0.60 m/s; V4 > 0.61 m/s) in forested streams.
The dotted lines indicate SI = 0.7, above which conditions were considered optimal.
Only seven species to which we built HSC curves for forested streams also
occurred in deforested streams. Cetopsorhamdia iheringi, Pseudopimelodus pulcher,
and Knodus moenkhausii did not occur in deforested streams and therefore could not be
included in the assessment of HSC transferability. In relation to depth, the mean
abundance of Astyanax altiparanae and Aspidoras fuscoguttatus was higher in optimal
habitats of deforested streams, considering the threshold established in forested streams
(t test, P < 0.01; Table II). The abundance of the other species did not differ among
27
optimal and non-optimal habitats in deforested streams (t test, P > 0.15; Table II). In
relation to water velocity, all species had higher abundances in optimal habitats of
deforested streams, except for Hypostomus ancistroides, according to the criteria
established for forested streams (t test, P < 0.02; Table II).
Table II. Transferability test of HSC from forested to deforested streams. The + sign
indicates that the HSC developed for forested streams was able to predict the abundance
of a given species in deforested streams, and the - sign indicates the opposite.
Species
Abundance
Depth Velocity
Astyanax altiparanae + +
Astyanax fasciatus - +
Piabina argentea - +
Characidium zebra - +
Aspidoras fuscoguttatus + +
Hypostomus ancistroides - -
Hypostomus nigromaculatus - +
DISCUSSION
In this study we established the preference of ten species of stream fishes for
depth and velocity conditions. In general, species showed well-defined preference
patterns for depth and velocity, reinforcing the importance of hydraulic variables as
predictors of fish distribution in lotic environments (Lamouroux et al., 1999; Vadas &
Orth, 2001; Schwartz & Herricks 2008; Leal et al., 2011). One of the patterns was the
association of Astyanax altiparanae, A. fasciatus, Knodus moenkhausii, and Piabina
argentea in deep, slow-flowing habitats. Indeed, the morphology of these species is
compatible with this habitat preference, as their compressed bodies, well-developed
caudal and pectoral fins, and lateral eyes represent morphological adaptations to live in
lentic habitats (Gatz, 1979; Oliveira et al., 2010). Behavioral patterns displayed by these
species also favor the occupation of lentic and deep habitats, as they are mid water
swimmers that feed on drift items in the water column (Casatti & Castro, 1998; Casatti
et al., 2001; Ceneviva-Bastos et al., 2010). Probably they are also favored in deep
habitats with a higher habitat volume (i.e., water column) and where slow velocity
28
requires less energy expenditure for movement and feeding (Casatti & Castro, 1998;
Bürhnheim, 2002).
An opposite pattern was observed for Characidium zebra, Cetopsorhamdia
iheringi, Pseudopimelodus pulcher, and Hypostomus nigromaculatus that showed
preference for shallow and fast-flowing habitats. These fish have morphological and
behavioral adaptations congruent with life in these conditions. For example, these
species share depressed bodies and expanded pectoral fins (Casatti et al. 2005; Casatti
& Castro, 2006), which allow short movements in fast-flowing environments (Watson
& Balon, 1984). In addition, the diet composed of periphyton (H. nigromaculatus) and
aquatic insects (remaining species) is consistent with the high availability of these items
in shallow, fast-flowing environments (Angermeier & Karr, 1983).
The conditions considered non-optimal in this study (SI < 0.7) are not
necessarily uninhabitable. Within this range, one could still discriminate between a
range of habitable or marginal conditions and another as uninhabitable (Thomas &
Bovee, 1993), which would represent the conditions in which a given species could not
live. Most species that preferred shallow, fast-flowing habitats were less plastic, with a
wider range of conditions that may be considered uninhabitable (SI values close to zero
in two classes of depth and velocity), compared with species that preferred deep, slow-
flowing habitats. The only exception to this pattern was Characidium zebra, whose non-
optimal, but occupied conditions could be considered habitable, as the SI values were
higher than 0.37. Species that preferred deep, slow-flowing habitats showed SI values
close to zero in a maximum of one class of depth and velocity. This indicates that
species that preferred shallow, fast-flowing habitats are more specialized and probably
more sensitive to changes in hydraulic conditions. In fact, the loss of rheophilic species
in streams from our study region has been attributed to changes in hydraulic conditions
and simplification of habitats (Casatti et al., 2006; Casatti et al., 2009a; Teresa and
Casatti, 2010). For example, Cetopsorhamdia iheringi and Pseudopimelodus pulcher
did not occur in the degraded streams studied by Casatti et al. (2009a), a fact that they
attributed to the loss of riffles from siltation and habitat simplification.
The habitat suitability criteria obtained for Hypostomus ancistroides and
Aspidoras fuscoguttatus were distinct from other species with respect to depth.
Normally, the HSC results in unimodal patterns (Thomas and Bovee, 1993; Vadas and
Orth, 2001; Strakosh et al., 2003), although different patterns can also be observed
(Freeman et al., 1997). The bimodal pattern observed for A. fuscoguttatus might be an
29
artifact of small sample sizes. Because this species occurred in only two out of three
forested streams, we did not consider the samples from the stream in which it was
absent to calculate the habitat availability. Consequently, only five samples fell into the
deepest depth class. Because the calculation of SI involves the relative abundance ratio
of the proportional availability of habitat, the occurrence of a few individuals (seven) in
two samples resulted in high SI values. One of the procedures that can be adopted in
this case is either group together two adjacent classes to increase sample size or not take
into account the SI value obtained (see Freeman et al., 1997). With the adoption of
these procedures, the preference pattern of A. fuscoguttatus would be similar to species
with a preference for shallow habitats, and therefore consistent with the pattern of
habitat use reported for this species (Araújo and Garutti, 2003). Although the criterium
developed here is consistent in deforested streams, this was due to the preference of this
fish for shallowest depth class, as no individual was recorded in the deepest depth class.
This strengthens our suggestion that the second peak of the HSC curve for this species
has no biological basis. The habitat use criterion for velocity indicated a preference for
fast-flowing habitats, but similar to Characidium zebra, other velocity conditions were
considered habitable (SI values > 0.32). This is consistent with the plasticity of this
species, which occurs in streams with different hydrological characteristics (Casatti et
al., 2009a,b).
The multimodal pattern of the depth preference curve of Hypostomus
ancistroides suggests that this factor has little influence in determining its abundance
patterns. On the other hand, velocity seems to be a good predictor of its abundance in
forested streams, via a more uniform relationship, indicating a preference for lentic
habitats, as also found by Uieda et al. (1997) and Casatti et al. (2005). Interestingly, this
pattern is uncommon for species of Hypostomus, which are usually associated with fast-
flowing habitats, like H. nigromaculatus in this study.
The analysis of HSC transferability from forested to deforested streams revealed
that the preference of fish species is partially consistent across environments of different
conservation status. Depth HSCs could be transferred to only two out of seven species
(29%), whereas velocity HSCs could be transferred to six out of seven species (86%).
These results indicate that the need to develop specific criteria for each condition is
dependent on the variable considered, with velocity having higher transferability than
depth. The failure to transfer HSCs may be due to different factors, such as (i) the
preferred conditions are not available, (ii) interaction with other factors in deforested
30
streams and (iii) failure of the criterion to describe optimal habitats. The first
explanation may, at least partially, justify the failure in transferring the criteria for
depth, especially for species whose optimum condition was assigned to the last depth
class (> 0.51 m) like Astyanax fasciatus. In deforested streams, only four segments (~
6%), three of which in the same stream, fell into the deepest depth class. According to
Bovee (1986), the successful transfer of HSC is dependent on the availability of
conditions in the site to which one wants to transfer the criterion.
The second explanation (interaction with other factors) may also apply if
environmental factors influenced fish differently in forested as deforested streams. For
example, the high abundance of marginal submerged vegetation in deforested streams,
as a result of the proliferation of exotic grasses, provide favorable microhabitats for
some species (Casatti et al., 2009) and may influence their response to other variables.
Such grasses can attenuate water flow along stream banks, providing slow-flowing
microhabitats even in stretches with high velocities. As a result, species with
preferences for slow-flowing habitats would be favored, notably Hypostomus
ancistroides in my study. In fact, individuals of this species were observed in the
vegetation on the banks of deforested streams, even in stretches with high water flow.
On the other hand, this species was restricted to the mid-channel in forested streams,
where submerged grasses are lacking and supposedly the influence of velocity is higher.
This could explain the failure to transfer the criterion for flow to H. ancistroides.
The failure in transferring HSCs also may be related to the fact that they do not
always describe a habitat preferred by species. For example, in this study, depth and
velocity co-varied in both forested and deforested streams (Pearson’s product-moment
correlation performed a posteriori, r = 0.57, P < 0.01; r = 0.29, P < 0.01, respectively),
so that the preference patterns observed for a variable (e.g., depth) can be influenced by
a second variable (velocity). If the pattern for depth preference of species is an artifact
of their association with velocity, then one would expect a distinct preference pattern
for this variable (i.e., failure to transfer HSCs) only if the velocity-depth relationship
was different between forested and deforested streams. In fact this happened, as a
higher proportion of segments classified into the fourth flow class (> 0.6 m/s) were
grouped into classes with higher depth (D2, D3, and D4) in the deforested streams, as
shown by chi-square test performed a posteriori (M-L Chi-square, P = 0.03). This could
explain the failure to transfer depth HSCs among streams for species that showed
31
preference for fast-flowing habitats, like Hypostomus nigromaculatus and Characidium
zebra.
The results indicated a high consistency in the response of species for velocity,
suggesting that this variable is a good predictor of species abundance and that its HSC
has great potential to be transferred to streams of different conservation status. The
scale of this study is relevant, because the transferability of criteria is usually tested
between streams (Thomas and Bovee, 1993) or between regional models built for
individual streams (Lamouroux et al., 1999). The data presented in this study should
represent part of an effort to develop criteria for Neotropic species, encompassing a
broader spectrum of conditions, as well as other variables known to be important for
fish, like substrate and cover (Vadas and Orth, 2001; Wright and Flecker, 2004; Leal et
al., 2011).
REFERENCES
Ahmadi-Nedushan B, St-Hilaire A, Bérubé M, Robichaud E, Thiémonge N, Bobée B.
2006. A review of statistical methods for the evaluation of aquatic habitat suitability
for instream flow assessment. River Research and Applications 22: 503–523.
Angermeier PL, Karr JR. 1983. Fish communities along environmental gradients in a
system of tropical streams. Environmental Biology of Fishes 9: 117–135.
Araujo RB, Garutti V. 2003. Ecology of a stream from upper Paraná river basin
inhabited by Aspidoras fuscoguttatus Nijssen & Isbrücker, 1976 (Siluriformes,
Callichthyidae). Brazilian Journal of Biology 63: 363–372.
Barletta M, Jaureguizar AJ, Baigun C, Fontoura NF, Agostinho AA, Almeida-Val
VMF, Val AL, Torres RA, Jimenes-Segura LF, Giarrizzo T, Fabré NN, Batista VS,
Lasso C, Taphorn DC, Costa MF, Chaves PT, Vieira JP, Corrêa MFM. 2010. Fish
and aquatic habitat conservation in South America: a continental overview with
emphasis on Neotropical systems. Journal of Fish Biology 76: 2118–2176.
Bovee KD. 1982. A guide to stream habitat analysis using the instream flow
incremental methodology. Instream Flow Information Paper 12. U.S. Fish and
Wildlife Service FWS?OBS-82/26: Washington, D. C.
Bovee KD. 1986. Development and evaluation of habitat suitability criteria for use in
the instream flow incremental methodology. Instream Flow Information Paper 21.
U.S. Fish and Wildlife Service Biological Report 86 (7): Washington, D. C.
32
Bürhnheim CM. 2002. Heterogeneidade de habitats: rasos x fundos em assembléias de
peixes de igarapés de terra firme na Amazônia Central. Revista Brasileira de
Zoologia 19: 889–905.
Casatti L, Castro RMC. 1998. A fish community of the São Francisco river headwater
riffles, southeastern Brazil. Ichthyological Exploration of Freshwaters 9: 229-242.
Casatti L, Castro RMC. 2006. Testing the ecomorphological hypothesis in a headwater
riffles fish assemblage of the rio São Francisco, southeastern Brazil. Neotropical
Ichthyology 4: 203–214.
Casatti L, Langeani F, Castro RMC. 2001. Peixes de riacho do Parque Estadual Morro
do Diabo, bacia do alto rio Paraná, SP. Biota Neotropica 1: 1–15.
Casatti L, Rocha FC, Pereira DC. 2005. Habitat use by two species of Hypostomus
(Pisces, Loricariidae) in southeastern brazilian streams. Biota Neotropica 5: 157–
165.
Casatti L, Langeani F, Silva AM, Castro RMC. 2006. Stream fish, water and habitat
quality in a pasture dominated basin, southeastern Brazil. Brazilian Journal of
Biology 66: 681–696.
Casatti L, Ferreira CP, Carvalho FR. 2009a. Grass-dominated stream sites exhibit low
fish species diversity and dominance by guppies: an assessment of two tropical
pasture river basins. Hydrobiologia 623: 173–189.
Casatti L, Veronezi-Júnior JL, Ferreira CP. 2009b. Dieta do cascudo Aspidoras
fuscoguttatus (Ostariophysi, Callichthyidae) em riachos com diferentes
características limnológicas e estruturais. Biota Neotropica 9: 1–9.
Ceneviva-Bastos M, Casatti L, Rossa-Feres DC. 2010. Meso and microhabitat analysis
and feeding habits of small nektonic characins (Teleostei: Characiformes) in
Neotropical streams. Zoologia 27: 191–200.
Freeman MC, Bowen ZH, Crance JH. 1997. Transferability of habitat suitability criteria
for fishes in warmwater streams. North American Journal of Fisheries Management
17: 20–31.
Gatz AJ Jr. 1979. Ecological morphology of freshwater stream fishes. Tulane Studies in
Zoology and Botany 21: 91–124.
Grinnell J. 1917. The niche-relationship of the California thrasher. Auk 34:427–33.
Groshens TP, Orth DJ. 1993. An assessment of the transferability of habitat suitability
criteria for smallmouth bass. Virginia Water Resources Research Center Bulletin
180.
33
Hutchinson MF. 1957. Concluding remarks. Cold Spring Harbour Symposium on
Quantitative Biology 22:415–427
IPT (Instituto de Pesquisas Tecnológicas do Estado de São Paulo). 1999. Diagnóstico
da situação atual dos Recursos Hídricos e estabelecimento de diretrizes técnicas
para a elaboração do Plano da Bacia Hidrográfica do Turvo/Grande. Relatório nº
40.515. Secretaria de Recursos Hídricos, Saneamento e Obras, São Paulo.
Jackson DA, Peres-Neto PR Olden JD. 2001. What controls who is where in freshwater
fish communities – the roles of biotic, abiotic and spatial factors. Canadian Journal
of Fisheries and Aquatic Sciences 58: 157–170.
Lamouroux N, Capra H, Pouilly M, Souchon Y. 1999. Fish habitat preferences in large
streams of southern France. Freshwater Biology 42: 673–687.
Leal CG, Junqueira NT, Pompeu PS. 2011. Morphology and habitat use by fishes of the
Rio das Velhas basin in southeastern Brazil. Environmental Biology of Fishes
90:143–157.
Oliveira EF, Goulart E, Breda L, Minte-Vera CV, Paiva LRS, Vismara MR. 2010.
Ecomorphological patterns of the fish assemblage in a tropical floodplain: effects of
trophic, spatial and phylogenetic structures. Neotropical Ichthyology 8: 569–586.
Power ME. 1984. Depth distribution of armored catfish: predator-induced resource
avoidance? Ecology 65: 523–528.
São Paulo. 2006. Relatório da qualidade ambiental do Estado de São Paulo 2006.
Secretaria do Estado do Meio Ambiente – MMA, São Paulo.
Schwartz JS, Herricks EE. 2008. Fish use of ecohydraulic-based mesohabitat units in a
low-gradient Illinois stream: implications for stream restoration. Aquatic
Conservation: Marine and Freshwater Ecosystems 18: 852–866.
Silva AM, Nalon MA, Kronka FJAN, Álvares CA, Camargo PB, Martinelli LA. 2007.
Historical land-cover/use in different slope and riparian buffer zones in watershed of
the State of São Paulo, Brazil. Scientia Agricola 64: 325–335.
Strakosh TR, Neumann RM, Jacobson RA. 2003. Development and assessment of
habitat suitability criteria for adult brown trout in southern New England rivers.
Ecology of Freshwater Fish 12: 265–274.
Teresa FB, Casatti L. 2010. Importância da vegetação ripária em região intensamente
desmatada no sudeste do Brasil: um estudo com peixes de riacho. Pan-American
Journal of Aquatic Sciences 5: 444–453.
34
Thomas JA, Bovee KD. 1993. Application and testing of a procedure to evaluate
transferability of habitat suitability criteria. Regulated Rivers: Research and
Management 8: 285–294.
Uieda VS, Buzzato P, Kikuchi RM. 1997. Partilha de recursos alimentares em peixes
em um riacho de serra no sudeste do Brasil. Anais da Academia Brasileira de
Ciências 69: 243–252.
Vadas RL Jr, Orth DJ. 2001. Formulation of habitat suitability models for stream fish
guilds: do the standard methods work? Transactions of the American Fisheries
Society 130: 217–235.
Watson DJ, Balon EK. 1984. Ecomorphological analysis of fish taxocenes in rainforest
streams of northern Borneo. Journal of Fish Biology 25: 371–384.
Wright JP, Flecker AS. 2004. Deforesting the riverscape: the effects of wood on fish
diversity in a Venezuelan piedmont stream. Biological Conservation 120: 439–447.
35
Manuscrito aceito para publicação na Ecology of Freshwater Fish
- CAPÍTULO II -
Influence of forest cover and mesohabitats types on functional and
taxonomic diversity of fish communities in Neotropical lowland
streams
36
ABSTRACT
In this study, we investigated how taxonomic and functional diversity of fish
communities is influenced by forest cover and mesohabitat types in Neotropical
lowland streams. We sampled fish fauna of 126 five meter-long mesohabitats using an
electrofishing unit in forested (n = 3) and deforested (n = 3) streams in the upper Paraná
River basin of southeastern Brazil. According to velocity depth, three mesohabitat types
have been considered: riffles (shallow and fast-flowing habitat), pools (deep and slow-
flowing habitat) and run (intermediate depth and velocity). Seven functional traits and
27 trait categories related to ecological, behavioral and life history aspects of fish were
considered. Our results indicate that forest cover and mesohabitat type influence
differently functional patterns of communities. While deforestation affects communities
primarily through changes in diversity (functional and taxonomic), mesohabitat types
determine changes in the functional composition. The increased diversity in deforested
mesohabitats is driven by a decrease of species turnover among habitat patches within
streams. This can be attributed to new feeding opportunities and microhabitat
availabilities in deforested streams and favors the occurrence species with a particular
set of traits, indicating a strong habitat-trait relationship.
Keywords: ichthyofauna, functional traits, riparian forest, guilds, environmental
impact.
37
INTRODUCTION
Deforestation of the riparian zone is one of the most important threats to the
conservation of freshwater fish species (Naiman & Décamps 1997; Sweeney et al.
2004). Deforestation primarily affects streams, whose physical and biological integrity
depends on native forests in riparian zones (Pinto et al. 2006). Changes resulting from
deforestation include modifications to aquatic communities and ecosystem services
(Sweeney et al. 2004). The consequences of deforestation on fish communities have
been primarily examined under a taxonomic perspective (Bojsen & Barriga 2002;
Casatti et al. 2010). These studies have shown changes in species composition and
indicators of species diversity such as the richness, diversity, and dominance of species
(Bojsen & Barriga 2002; Wright & Flecker 2004; Teresa & Casatti 2010). However,
these metrics represent a limited set of traits of biological communities because they do
not take accont the relative functional differences between species and often fail to
detect changes caused by anthropogenic interference (Rabeni & Smale 1995; Ernst et al.
2006).
Recently, there has been remarkable growth in the number of studies
incorporating aspects of the functional diversity of communities (Petchey & Gaston
2002; Hoeinghaus et al. 2007; Poff et al. 2010). This approach takes into account both
the functional traits of species and the ecological roles they play in the ecosystem, and it
provides additional information otherwise not obtainable through a taxonomic approach
(Eros et al. 2009). In addition, the functional structure of communities is considered a
better predictor of ecosystem functioning than species diversity itself (Tilman et al.
1997; Mokany et al. 2008; Moulliot et al. 2011). The study of the anthropogenic
impacts on the functional structure of communities is fundamental because it may
predict not only communities patterns but also processes operating in ecosystems
(Mayfield et al. 2010).
Studies examining the organization of fish communities based on their
functional composition have shown that habitat use is influenced by the functional
traits of a species (Poff & Allan 1995; Jones III et al. 1999; Irz et al. 2007; Eros et al.
2009; Montaña & Winemiller 2010). Thus, anthropogenic impacts affecting the
physical structure of habitats have the potential to influence functional aspects of fish
communities. In support of this idea, the addition and/or loss of species in deforested
streams appears to be non-random (Jones III et al. 1999; Goldstein & Meador 2005),
38
suggesting that deforestation limits or favors certain species according to their traits.
For example, Rabeni & Smale (1995) demonstrated that species that rely on hard
substrates for reproduction tend to disappear from environments subject to siltation.
In addition to functional composition, functional diversity is another way in
which community structure can be described. Indices that have been developed recently
deal with functional traits within a multivariate framework (e.g., Petchey & Gaston
2006; Villéger et al. 2008). This approach has also been useful in identifying functional
changes along gradients of deforestation (Ernst et al. 2006; Barragán et al. 2011), and
its use in streams is promising because little is known about how different components
of biodiversity in these environments respond to human impacts.
Streams are heterogeneous systems composed of a mosaic of patches recognized
at different scales (Frissel et al. 1986). At the mesohabitat scale, three physiognomies
are easily recognized according to their geomorphological and hydrological features:
riffles represent shallow, fast waters; pools are deep, slow waters that usually have fine
substrate; runs are intermediate between riffles and pools. These mesohabitats influence
the composition and structure of fish fauna (Bührnheim 2002; Langeani et al. 2005;
Schwartz & Herricks 2008; Rezende et al. 2010), as well as some characteristics of fish
species (e.g., diet, reproductive strategies and preference for substrate) (Berkman &
Rabeni 1987; Casatti & Castro 1998; Goldstein & Meador 2004). For example, fish
with less diversified diets predominate in riffles, while the increased variety of food
items available in pools would support communities composed of more generalist
species (Angermeier & Karr 1983; Berkman & Rabeni 1987). Thus, it is expected that
different mesohabitats and their associated communities respond differently to
environmental disturbances.
To our knowledge, no studies have investigated how deforestation influences
the functional diversity and composition of Neotropical stream fish assemblages in
different mesohabitat types. Therefore, in this study, we investigated whether
functional diversity metrics vary between forested and deforested mesohabitats and/or
among mesohabitats types. The study was conducted in forested and deforested streams
of southeastern Brazil, one of the most developed regions of the country, where
deforestation affects more than 75% of the riparian zone of streams (Silva et al. 2007).
To characterize the fish fauna, we used multiple traits including behavioral and
ecological (e.g., life-history, feeding tactic, current preference) aspects.
39
MATERIAL AND METHODS
Study area and general sampling design
All streams studied are located in the drainage of the São José dos Dourados
River in the upper Paraná River watershed of southeastern Brazil (Figure 1). The
forested streams (n = 3) are located within the greatest forest fragments (seasonal
semideciduous forest > 1.600 ha) of the basin and the deforested streams (n = 3) lack
arboreal vegetation in their riparian zone and are located in areas dominated by pasture
(Figure 1). There are no representative forest remnants associated with deforested
streams upstream or downstream of the study sites. These streams were selected from a
group of 80 previously visited streams in the region. They were chosen because they are
representative of forested and deforested streams, present mesohabitat variability,
represent second or third order streams with similar volume of water, have no evidence
of chemical pollution, and differ only in respect to structural features resulting from
differences in forest cover, as described below. Due to the scarcity of large forest
fragments and the difficulty of finding comparable deforested areas, sites were not
randomly distributed in the landscape but rather were clustered in specific areas.
However, all these streams were considered to be influenced by the same regional
species pool (e.g., Casatti et al. 2006; 2009), as the distances between watercourses do
not explain the variation in the species composition among streams distributed within
the study area (Casatti et al. 2012).
40
Fig. 1. Map of the study area, showing São Paulo state (A), São José dos Dourados basin (B)
and location of the forested (dark circles) and deforested (open circles) streams. Large forest
fragments (>100ha) present in the basin are shown.
We sampled 126 five meters-long mesohabitats (with between 16 to 24 sites in
each stream), 62 of which were in deforested streams and 64 in forested streams. The
mesohabitat were chosen in order to represent the structural variation present in each
stream. We sampled fish fauna and environmental descriptors in each mesohabitat from
September to November 2009 and from April to July 2010.
Data sampling
We isolated each mesohabitat upstream and downstream with fishing nets before
sampling. We sampled fish with successive electrofishing passes along mesohabitats,
and the capture of fish ended when individuals were no longer captured after the passage
of electrodes through the test mesohabitat. Fish collected from each mesohabitat were
placed in separate containers, fixed in 10% formalin and preserved in 70% ethanol after
72 hours. All of the fish species found were deposited at the Laboratory of Ichthyology
collection (DZSJRP), Department of Zoology and Botany, São Paulo State University,
41
São José do Rio Preto, Brazil (see Table A1, Supporting Information, for a species list
of each stream).
We collected data on the physical structure of habitat after fish sampling. To do
this, we sub-divided each mesohabitat into four transversal, equidistant transects where
we recorded the following data: width, measured at the surface and at the middle of the
water column; depth, measured at five equidistant points from one margin to another;
velocity (with a mechanical flowmeter), measured at the middle of the water column;
streambed substrate composition (i.e., the relative proportion of each substrate
component) visually estimated according to Cummins (1962) and the volume occupied
by the submerged vegetation (roots, leaves and stems from the submerged terrestrial
vegetation) on the banks of each margin, based on the measures of the height and the
width of vegetation on stream banks of each transect. We classified the streambed
substrate as either stable, which allows colonization by benthic fauna (in the
consolidated clay, gravel, pebble and rock) or unstable, with sand and unconsolidated
clay. The values of each mesohabitat unit were calculated by the mean of the values
obtained in the transects.
The mesohabitat units were classified according to their depth and velocity
characteristics. These two variables were chosen because they represent important
predictors of fish community structure in streams (Lamouroux et al. 1999; Bürhnheim
2002). Moreover, both forested and deforested streams had similar mesohabitats
considering these two variables. After performing a hierarchical agglomerative cluster
analysis with Euclidean distance by using UPGMA method, mesohabitats were
classified in three types: (i) shallow, fast-flowing habitat; (ii) deep, slow-flowing
habitat; (iii) intermediate conditions in relation to these two descriptors. These three
groups of mesohabitats were named riffles, pools, and runs, respectively. We then
pooled and averaged the structural descriptors data to represent the mesohabitats in each
stream.
A Principal Component Analysis (PCA) was used to describe the variation of
data in the five descriptors of habitat structure measured (depth, width, velocity,
percentage of stable substrates and volume of submerged marginal vegetation). The
variables were first standardized so that their means were zero and their standard
deviation equal to one. We also carried out a Multivariate Analysis of Variance
(MANOVA) to test the significance of the structural patterns exhibited by mesohabitat
types and streams with different forest cover.
42
Traits and trait categories
We adopted a multidimensional approach to characterize fish species traits. We
considered seven traits and 27 trait categories, including behavioral and ecological
aspects (Table 1). We defined traits and their categories based on the species
information recorded in literature and in this study. The species characterization
followed a binary approach, with each species being associated to one trait category for
each trait. We adopted this practice because there is little information on the variability
of traits for many species. This approach was conservative, but it was sufficient for an
analysis of the communities’ general patterns (Poff et al. 2006). The most common trait
categories expressed by each species were obtained from the literature, from data
collected in this study and based on the author’s experience (Table A2 – Supporting
Information).
43
Table 1. Description of seven functional traits and 27 trait categories used to
characterize fish species with their respective codes, including aspects of ecology and
feeding behavior, habitat use, life history and tolerance.
Traits Trait categories
Ecology and feeding behavior
Diet1 Omnivore, periphytivore, insectivore, carnivore or
detritivore
Feeding tactic,3 Water column feeding, substrate speculation, stealth
predation, grazing, nocturnal predation or digging
Habitat use2,4
Velocity preference Slow, intermediate or fast
Stable substrate preference High, moderate or low
Habitat use
Surface, nektonic, nectobenthic, margins or benthic
Life history5
Size Small (< 50 mm), intermediate (50-150 mm) or large
(> 150 mm)
Tolerance6
Tolerance to hypoxia Tolerant or intolerant 1 Diet analysis made by us 2 Underwater observations made by us 3 Keenleyside (1979), Sazima (1986), Sabino & Zuanon (1998) 4 Casatti & Castro (1998), Casatti et al. (2001), Casatti (2002) and personal observation 5 Biometry made by us 6 Kramer & Mehegan (1981), Araujo & Garutti (2003), Bozzetti & Schulz (2004) and
Casatti et al. (2006)
For additional information on habitat use and feeding tactics, we made
observations during 1,200 minutes while snorkeling along 100-meter stretches in the six
studied streams. We also performed a diet analysis by examining the stomach contents
of 10 specimens of each species (except for those with less than 10 individuals)
representative of the streams sampled. Only species with more than five individuals
were analyzed. We identified and classified food items into broad categories (detritus,
insects, other invertebrates, periphyton, fragments of higher plants and algae). When
44
assigning a trophic group to a species, we considered the most important items in their
diet. Species with less than five individuals had their diet characterized based on data
available in literature.
We evaluated the independence or redundancy of traits prior to the calculation of
functional diversity metrics. This allowed the identification of non-informative and
autocorrelated traits that should be removed from the analysis to avoid under or
overestimation of functional diversity (Petchey & Gaston 2006). We performed a matrix
correlation analysis (2Stage analysis) using the software PRIMER 6.0 (Clarke & Gorley
2006) from binary matrices for each trait (species as samples and trait categories as
variables). The results of this analysis indicated that the seven traits had low
autocorrelation (< 0.53), meaning that the previously selected traits represent different
aspects of species. Thus, we considered all seven traits in additional analyses.
Functional diversity metrics
We used functional richness, functional evenness (Villéger et al. 2008), and
functional dispersion (Laliberte & Legendre 2010) to represent functional diversity.
These indices measure different aspects of species composition in the multidimensional
functional space, and they are independent of species richness (Villéger et al. 2008) and
very flexible when using distinct trait types (Laliberté & Legendre 2010). The
functional richness represents the volume of the convex hull generated by positioning
species distribution in the n-dimensional species trait space (Villéger et al. 2008). The
functional evenness measures the regularity with which the functional space is occupied
by species weighted by their abundance (Villéger et al. 2008). The functional evenness
is high when abundance is evenly distributed among species or when the functional
distance between species is regular (Villéger et al. 2008). The functional dispersion was
first developed as a measure of beta diversity (Anderson et al. 2006), but it was recently
incorporated as a metric of functional diversity. In contrast to functional richness,
functional dispersion may account for the abundance of species within communities,
and it is not highly affected by outliers (Laliberté & Legendre 2010). This index
consists of measuring the mean distance of the multidimensional functional space from
each species to the centroid of all species, which in turn is dislocated towards the most
abundant species (Laliberté & Legendre 2010). We used the function dbFD in the
package FD of the R software (R Development Core Team 2011) to calculate the three
indices of functional diversity. This package incorporates changes that make the indices
45
proposed by Villéger et al. (2008) more flexible, making possible the use of exclusively
qualitative traits (Laliberté & Legendre 2010), as is the case in our study.
Data analysis
We used two-way ANOVA to investigate how the number of species within each
trait category varied depending on forest cover and mesohabitats types in the R software
(R Development Core Team 2011). Because part of the data did not fit assumptions of
normality and/or homogeneity of variance, we used a more conservative significance
level (p ≤ 0.01) to minimize type I error. We evaluated the effects of forest cover and
mesohabitats on functional diversity metrics and species richness with the same
analytical strategy. In order to evaluate the contribution of species exclusive to forested
or deforested streams to the general patterns of diversity, we also conducted an analysis
considering only the species shared between these groups of streams. In these cases, data
had normal distribution and homogeneous variances; therefore, we used p ≤ 0.05 to
attribute significance. We considered different mesohabitats types and streams as
independent samples. Prior to hypothesis testing, we pooled data by mesohabitat type
within streams, producing a matrix with 18 samples corresponding to the three types of
mesohabitats sampled in the six streams.
RESULTS
The first two axes of PCA accounted for 79.5% of the variation in the physical
habitat structure data, with axis 1 accounting for 54.7% of the variation and
representing water velocity, percentage of stable substrate, and depth (eigenvectors = -
0.536, - 0.492, 0.485, respectively; Figure 2). The second axis accounted for 24.8% of
the variation and represented the width and the amount of submerged marginal
vegetation (eigenvectors = 0.681, -0.569, respectively; Figure 2). The first axis
distinguished the mesohabitats with fast water velocity, less depth and larger proportion
of stable substrates associated with riffles than those with slow water velocity, high
depth and a low proportion of stable substrates associated with pools; intermediate
values along this gradient were associated with runs (Figure 2). Axis 2 differentiated
forested from deforested streams. Accordingly, deforested streams were associated to a
higher amount of submerged marginal vegetation, whereas forested streams were
associated to greater width (Figure 2). The separation of the samples in accordance to
46
mesohabitats and forest cover was corroborated by MANOVA that revealed a
significant main effect of both factors (p < 0.01).
Fig. 2. Ordination (PCA) of the three types of mesohabitats (triangles = riffles; circles =
runs; squares = pools) in forested streams (dark symbols) and deforested (open
symbols) according to the descriptors of physical habitat.
A total of 1,335 and 2,395 individuals and 31 and 33 species were caught in the
forested and deforested streams, respectively. Ten and twelve species were exclusive of
forested and deforested streams, respectively. We found that the richness of species
within the trait categories varied according to forest cover and mesohabitats, but this
response was specific for each trait category. The following traits were represented by a
higher number of species in deforested streams (forest cover effect: p < 0.01):
detritivorous and omnivorous diet, use of margins, use of surface, low preference for
stable substrates, small and large body size, preference for low and intermediate
velocity, and hypoxia tolerance. On the other hand, forested streams presented a higher
richness of fish with a high preference for stable substrate (effect of forest cover: p <
0.01). In relation to mesohabitats, there was an increase in the richness of species
showing omnivorous diet, water-column and digging feeding tactics, nektobenthic
habitat, low preference for stable substrates and preference for slow and intermediate
velocities along the riffle-pool gradient (effect of mesohabitat: p < 0.01). On the other
47
hand, some trait categories became more representative along the pool-riffle gradient,
such as an insectivorous diet, benthic habitat, substrate speculation, high preference for
stable substrates, and a preference for fast velocity (effect of mesohabitat: p < 0.01). A
carnivorous diet was exclusive of fish from deforested streams and found with great
frequency in pools (interaction “forest cover” vs. “mesohabitats”; p < 0.01). Forest
cover or mesohabitats did not influence (p > 0.04) the following trait categories:
periphytivorous diet, nektonic habitat, grazing, stealth predation, nocturnal predation,
moderate preference for stable substrate, and preference for intermediate velocity.
The taxonomic richness and functional richness was higher in deforested than in
forested mesohabitats (effect of forest cover: p = 0.01), but did not differ among
mesohabitats (effect of mesohabitat: p = 0.34, Figure 3). Neither forest cover nor
mesohabitat had an effect on functional evenness (p > 0.42, Figure 3). The functional
dispersion was higher in deforested streams, but this effect was mesohabitat-dependent
(interaction “forest cover” vs. mesohabitats: p = 0.01), being significant only in riffles
(Tukey test, p < 0.01; Figure 3). The qualitative results were the same for the three
indices when only those species shared between forested and deforested streams were
considered.
48
Fig. 3. Taxonomic richness, functional richness, evenness and dispersion in forested
(white bars) and deforested (gray bars) streams in relation to mesohabitats types. Data
presented as mean ± standard deviation. * Indicates statistically significant difference
between forested and deforested riffles and also among the remaining forested
mesohabitats (Tukey test, p ≤ 0.01).
DISCUSSION
The evaluation of the effects of anthropogenic impacts on biological
communities and ecosystems are of central importance to conservation. In this context,
one of the current challenges is to understand the response of different components of
biodiversity to environmental change in heterogeneous environments. Our results show
that forest cover has a consistent effect on the functional diversity of stream fish
assemblages across mesohabitats types, with deforested mesohabitats presenting a
higher diversity values. These results are in contrast with other studies that investigated
the effects of anthropogenic impacts on the taxonomic and functional diversity of
different animal groups. For example, Flynn et al. (2009) noted a reduction in
taxonomic and functional diversity of mammal and bird communities along a gradient
of land use. Additionally, Ernst et al. (2006) reported a decrease in the functional
49
diversity in amphibian communities from different deforestation regimes, but no effect
on taxonomic richness along the gradient of degradation. By contrast, taxonomic
richness increased as functional diversity decreased in estuarine fish communities of a
degraded lake associated with the Gulf of Mexico (Villéger et al. 2010). These patterns
are consistent with habitat loss and simplification in degraded ecosystems, which in turn
result in the loss of species and functional traits (Flynn et al. 2009) or in the
replacement of species with unique traits by functionally redundant species (Ernst et al.
2006; Villéger et al. 2010).
The increased diversity in deforested mesohabitats would be due to the addition
of functionally distinct species in these streams. However, while species were added to
the deforested streams, an equivalent number of species occurred exclusively in the
forested streams, corroborating other studies and evidencing a pattern of species
substitution along the gradient of deforestation (Bojsen & Barriga 2002; Casatti et al.
2009). When considering all of the species and excluding species absent in either
deforested or forested streams, the congruence of the patterns of functional diversity
suggest that the exclusive species influenced the functional diversity of both group of
streams in a similar way. For example, while a carnivorous diet (Hoplias malabaricus
and Synbranchus marmoratus) and the combination of a periphytivorous diet and the
use of the margins (Hisonotus francirochai) are unique traits of species occurring
exclusively in deforested streams, an insectivorous diet with speculation tactics
(Cetopsorhamdia iheringi and Pseudopimelodus pulcher) and with a large body size
(Pimelodus maculatus) are unique combinations represented by species exclusive to
forested streams. Thus, although the exclusive species may be important to increasing
the functional diversity within each group of streams, they were not responsible for the
differences between them. The increase of taxonomic and functional diversity in
deforested mesohabitats was mainly due to higher evenness in the distribution of
species within these streams. This conclusion is based on the fact that deforested
streams presented a higher mean richness per mesohabitat, but a similar number of
species to forested ones when considering the streams as a whole (31 and 33 species,
respectively). Furthermore, even considering only the species shared between both
groups of streams, deforested mesohabitats show a higher taxonomic and functional
richness. This indicates that deforested mesohabitats support a higher taxonomic and
functional diversity mainly by providing conditions that increase the frequency of
species that were formerly present.
50
The functional evenness in both groups of streams was relatively high (> 0.6)
regardless of mesohabitat, indicating that the functional space was evenly occupied.
These results also suggest that the species added in deforested mesohabitats filled
unoccupied niches, leading to a more even use of resources. In fact, when the number of
species in a community increases, local assemblages can accommodate new species by
increasing the volume of niche space (Barragán et al. 2011). The increase of the niche
space in deforested streams may be represented by rare or formerly absent opportunities
that become available, favoring the occurrence of opportunistic and/or invasive species
(Lorion & Kennedy 2009). Some of the new opportunities in deforested streams can be
represented by the greater availability of microhabitats (Casatti et al. 2009; Teresa &
Casatti 2010) and food items, such as detritus (Rocha et al. 2009), algae and periphyton
(Rounick et al. 1982; Bojsen & Barriga 2002; Lorion & Kennedy 2009), and this would
favor, for example, the occurrence of omnivorous species. One the most prominent
habitat alterations recorded was the greater abundance of marginal herbaceous
vegetation in deforested streams (see Figure 2). This vegetation is composed mainly of
exotic grasses (Brachiaria spp.) that tend to invade stream banks of deforested areas
(Casatti et al. 2009) and may provide substrate for shelter and reproduction (Growns et
al. 2003; Teresa & Romero, 2010). This would favor the occurrence of species usually
associated with submerged vegetation, like Crenicichla britskii, Gymnotus carapo,
Hoplias malabaricus, and Hysonotus francirochai. The high abundance of marginal
submerged vegetation may have contributed to the higher richness of species with a
preference for slow or intermediate velocities, as they can attenuate water velocity
along stream banks and provide slow-flowing microhabitats even in areas with fast
velocities. The vegetation can also trap detritus in their leaves and roots (Rocha et al.
2009), which are favored by detritivores (Casatti et al. 2009) like Corydoras aeneus,
Geophagus brasiliensis, Hypostomus ancistroides and Poecilia reticulata, all of which
became more frequent in the deforested streams.
Some of the species favored by submerged marginal vegetation (e.g., G. carapo,
G. brasiliensis and H. malabaricus) present also large body, which may also face a
lesser risk of predation by birds and mammals when in vegetated environments.
Vegetated habitats may also provide camouflage and cover to conceal prey (Padial et al.
2009), which could explain the higher richness of this trait in deforested mesohabitats.
However, it is important to note that the association between trait categories represented
by a few species with the factors investigated here should be taken with caution, as they
51
are more likely to autocorrelate with other trait categories. Hoplias malabaricus and
Synbranchus marmoratus can be cited as examples of carnivores that have been
associated with degraded environments (Casatti et al. 2006; Casatti et al. 2009, this
study), but other traits (e.g., hypoxia tolerance or margin habitat) rather than the diet
may be important in determining their association with deforested streams.
Furthermore, the higher richness of carnivorous fish in pools can also be associated with
the use of these sites as refuge during periods of inactivity (Power 1984).
Traits like tolerance to hypoxia, small body size and the use of surface may
favor the colonization of shallow, warm, vegetated microhabitats and are shared by
opportunistic species like Phalloceros harpagos, Serrapinnus notomelas and the exotic
Poecilia reticulata, which are dominant in the grass-dominated banks of degraded
streams (Casatti et al. 2009). The increased richness of fish tolerant to hypoxia in
deforested streams is probably related to rising temperatures and reduced oxygen levels,
conditions commonly found in silted deforested streams (Heartsill-Scalley & Aide
2003; Casatti et al. 2006). The higher number of species with a low preference for
stable substrates reflects the physical alterations commonly reported for structurally
degraded streams (Gregory et al. 1991; Walser & Bart 1999; Casatti et al. 2006) that
experienced higher fine-sediment deposition, as shown in the PCA plot (see Figure 2).
Despite the fact that mesohabitat types did not influence the diversity patterns of
the communities, they were important in determining their functional composition. The
dominance of fish species with an omnivorous diet, nektonics, and water column
feeders showing low preference for stable substrates and a preference for slow and
intermediated water velocity in pools is compatible with the low current velocity and
often rich food resources present in this mesohabitat (Bürhnheim 2002; Angermeier &
Karr 1983; Langeani et al. 2005; this study). These traits are shared by a group of
species, mostly small characins, which represent a significant proportion of the fish
communities in Neotropical streams (Castro 1999). The nektobenthic habitat and
digging tactic are traits found with greater frequency in pools and are shared by
detritivorous/omnivorous species that feed in soft bottom sites. In contrast to the
patterns observed in the pools, a higher number of species with a high preference for
stable substrates and fast velocity habitats were associated with riffles, which is
compatible with the hydraulic conditions and substrate composition present in these
mesohabitats (Casatti & Castro 1998; Schwartz & Herricks 2008). The relationship
between traits and conditions present in the riffles was more pronounced in forested
52
streams, as forested riffles had a lower functional dispersion than deforested riffles and
other mesohabitats. This result indicates that, on average, species are functionally more
similar to each other on these sites. In fact, two groups of species predominated in
forested riffles: benthic, grazers and periphytivorous fish (e.g., Hypostomus
nigromaculatus) and benthic, speculators and insectivorous fish (e.g., Pseudopimelodus
pulcher and Cetopsorhamdia iheringi). The fast water velocity and limited availability
of food items in the riffles (Angermeier & Schlosser 1989) might represent more
restrictive environmental filters, limiting the occurrence of a great diversity of traits and
resulting in highly similar assemblages. The increased functional dispersion in
deforested streams might indicate a relaxation of these filters. For example, the
constraint imposed by the fast velocity of water would be minimized in these streams by
an abundance of submerged marginal vegetation (Pusey & Arthington 2003).
In summary, our results indicate that forest cover and mesohabitat types
influence different aspects related to the functional patterns of communities. While
deforestation affects communities primarily through changes in diversity (both
functional and taxonomic), mesohabitat types mainly determine changes in the
functional composition. The increased diversity in deforested mesohabitats is driven by
a decrease in species turnover among habitat patches within streams. The new
opportunities available in deforested streams favor species with a particular set of traits
(e.g., detritivores, species inhabiting stream margins, tolerant to hypoxia), indicating a
strong habitat-trait relationship. Knowledge of how functional patterns and processes
operate on aquatic communities over a wider range of degradation contexts, and
whether the patterns observed herein are consistent with other scales, are promising
topics for future research.
REFERENCES
Anderson, M.J.,Ellingsen, K.E. & McArdle, B.H. 2006. Multivariate dispersion as a
measure of beta diversity. Ecology Letters. 9: 683-693.
Angermeier, P.L. & Karr, J.R. 1983. Fish communities along environmental gradients
in a system of tropical streams. Environmental Biology of Fishes 9: 117–135.
Angermeier, P. & Schlosser, I.J. 1989. Species-area relationship for stream fishes.
Ecology 70: 1450–1462.
53
Araujo, R.B. & Garutti, V. 2003. Ecology of a stream from upper Paraná river basin
inhabited by Aspidoras fuscoguttatus Nijssen and Isbrücker, 1976 (Siluriformes,
Callichthyidae). Brazilian Journal of Biology 63: 363–372.
Barragán, F., Moreno, C.E., Escobar, F., Halffter, G. & Navarrete, D. 2011. Negative
impacts of human land use on dung beetle functional diversity. PLoS One 6:
e17976.
Berkman, H.E. & Rabeni, C.F. 1987. Effect of siltation on stream fish communities.
Experimental Biology of Fishes 18: 285–294.
Bojsen, B.H. & Barriga, R. 2002. Effects of deforestation on fish community structure
in Ecuadorian Amazon streams. Freshwater Biology 47: 2246–2260.
Bozzetti, M. & Schulz, U.H. 2004. An index of biotic integrity based on fish
assemblages for subtropical streams in southern Brazil. Hydrobiologia 529: 133–
144.
Bürhnheim, C.M. 2002. Heterogeneidade de habitats: rasos x fundos em assembléias de
peixes de igarapés de terra firme na Amazônia Central. Revista Brasileira de
Zoologia 19: 889–905.
Casatti, L. 2002. Alimentação dos peixes em um riacho do Parque Estadual Morro do
Diabo, bacia do Alto Rio Paraná, sudeste do Brasil. Biota Neotropica 2: 1-15.
Casatti, L. 2010. Alterações no Código Florestal Brasileiro: impactos potenciais sobre a
ictiofauna. Biota Neotropica 10: 31–34.
Casatti, L. & Castro, R.M.C. 1998. A fish community of the São Francisco river
headwater riffles, southeastern Brazil. Ichthyological Exploration of Freshwaters 9:
229–242.
Casatti, L., Langeani, F. & Castro, R.M.C. 2001. Peixes de riacho do Parque Estadual
Morro do Diabo, bacia do alto rio Paraná, SP. Biota Neotropica 1: 1–15.
Casatti, L., F. Langeani, A.M. Silva & Castro, R.M.C. 2006. Stream fish, water and
habitat quality in a pasture dominated basin, southeastern Brazil. Brazilian Journal
of Biology 66: 681–696.
Casatti, L., Ferreira, C.P. & Carvalho, F.R. 2009. Grass-dominated stream sites exhibit
low fish species diversity and dominance by guppies: an assessment of two tropical
pasture river basins. Hydrobiologia 632: 273–283.
Casatti, L., Teresa, F.B., Gonçalves-Souza, T., Bessa, E., Manzotti, A.R., Gonçalves,
C.S. & Zeni, J.O. 2012. From forests to cattail: how does the riparian zone
influence stream fish? Neotropical Ichthyology. In press.
54
Castro, R.M.C. 1999. Evolução da ictiofauna de riachos sul-americanos: padrões gerais
e possíveis processos causais. In: Caramaschi, E.P., Mazzoni, R., Bizerril, C.R.S.F.,
Peres-Neto, P.R., eds. Ecologia de Peixes de Riachos: Estado Atual e Perspectivas.
Rio de Janeiro: Oecologia Brasiliensis, v. VI, pp. 139–155.
Clarke, K.R. & Gorley, R.N. 2006. PRIMER v6: User Manual/Tutorial. Plymouth:
PRIMER-E. 192 pp
Cummins, K.H. 1962. An evaluation of some techniques for the collection and analysis
of benthic samples with special emphasis on lotic waters. The American Midland
Naturalist 67: 477–504.
Ernst, R., Linsenmair, K.E. & Rödel, M.O. 2006. Diversity erosion beyond the species
level: dramatic loss of functional diversity after selective logging in two tropical
amphibian communities. Biological Conservation 133: 143–155.
Eros, T., Heino, J., Schmera, D.& Rask, M. 2009. Characterising functional trait
diversity and trait–environment relationships in fish assemblages of boreal lakes.
Freshwater Biology 54: 1788–1803.
Flynn, D.F.B., Gogol-Prokurat, M., Nogeire, T., Molinari, N., Richers, B.T., Lin, B.B.,
Simpson, N., Mayfield, M.M. & DeClerck, F. 2009. Loss of functional diversity
under land use intensification across multiple taxa. Ecology Letters 12: 22–33.
Goldstein, R.M. & Meador, M.R. 2004. Comparisons of fish species traits from small
streams to large rivers. Transactions of the American Fisheries Society 133: 971–
983.
Goldstein, R.M. & Meador, M.R. 2005. Multilevel assessment of fish species traits to
evaluate habitat degradation in streams of the Upper Midwest. North American
Journal of Fisheries Management 25: 180–194.
Gregory, S.V., Swanson, F.J., McKee, W.A. & Cummins, K.W. 1991. An ecosystem
perspective of riparian zones. BioScience 41: 540–551.
Growns, I., Gerke, P.C., Astles, K.L. & Pollard, D.A. 2003. A comparison of fish
assemblages associated with different riparian vegetation types in the Hawkesbury–
Nepean river system. Fisheries Management & Ecology 10: 209–220.
Heartsill-Scalley, T. & Aide, T. M. 2003. Riparian vegetation and stream condition in a
tropical agriculture-secondary forest mosaic. Ecological Applications 13: 225–234.
Hoeinghaus, D.J., Winemiller, K.O. & Birnbaum, J.S. 2007. Local and regional
determinants of stream fish assemblage structure: inferences based on taxonomic vs.
functional groups. Journal of Biogeography 34: 324–338.
55
Irz, P., Michonneau, F., Oberdorff, T., Whittier T.R., Lamoroux, N., Mouillot, D. &
Argillier, C. 2007. Fish community comparisons along environmental gradients in
lakes of France and north-east USA. Global Ecology and Biogeography 16: 350–
366.
Jones III, E.B.D., Helfman, G.S., Harper, J.O. & Bolstad, P.V. 1999. Effects of riparian
forest removal on fish assemblages in southern Appalachian streams. Conservation
Biology 13: 1454–1465.
Keenleyside, M.H.A. 1979. Diversity and adaptation in fish behaviour. Berlin: Springer.
208 pp
Kramer, D.L. & Mehegan, J.P. 1981. Aquatic surface respiration, an adaptative
response to hypoxia in the guppy, Poecilia reticulata (Pisces, Poeciliidae).
Environmental Biology of Fishes 6: 299–313.
Laliberté, E. & Legendre, P. 2010. A distance-based framework for measuring
functional diversity from multiple traits. Ecology 91: 299–305.
Lamouroux, N., Capra, H. Pouilly, M. & Souchon, Y. 1999. Fish habitat preferences in
large streams of southern France. Freshwater Biology 42: 673–687.
Langeani, F., L. Casatti, H.S. Gameiro, A.B. Carmo & Rossa-Feres, D.C. 2005. Riffle
and pool fish communities in a large stream of southeastern Brazil. Neotropical
Ichthyology 3: 305–311.
Lorion, C. M. & Kennedy, B. P. 2009. Riparian forest buffers mitigate the effects of
deforestation on fish assemblages in tropical headwater streams. Ecological
Applications 19: 468–479.
Mayfield, M.M., Bonser, S.P., Morgan, J.W. Aubin, I., McNamara, S. & Vesk, P.A.
2010. What does species richness tell us about functional trait diversity?
Predictions and evidence for responses of species and functional trait diversity to
land-use change. Global Ecology and Biogeography 19: 423–431.
Mokany, K., Ash, J. & Roxburgh, S. 2008. Functional identity is more important than
diversity in influencing ecosystem processes in a temperate native grassland.
Journal of Ecology 96: 884–893.
Montaña, C.G. & Winemiller, K.O. 2010. Local-scale habitat influences morphological
diversity of species assemblages of cichlid fishes in a tropical floodplain river.
Ecology of Freshwater Fish 19: 216–227.
56
Mouillot, D., Villéger, S., Scherer-Lorenzen, M. & Mason, N.W.H. 2011. Functional
structure of biological communities predicts ecosystem multifunctionality. PLoS
One 6: e17476.
Naiman, R.J. & Décamps, H. 1997. The ecology of interfaces: riparian zones. Annual
Review of Ecology and Systematics 28: 621–658.
Padial, A.A., Thomaz, S.M. & Agostinho, A.A. 2009. Effects of structural
heterogeneity provided by the floating macrophyte Eichhornia azurea on the
predation efficiency and habitat use of the small Neotropical fish Moenkhausia
sanctaefilomenae. Hydrobiologia 624: 161–170.
Petchey, O.L. & Gaston, K.J. 2002. Functional diversity (FD), species richness and
community composition. Ecology Letters 5: 402–411.
Petchey, O.L. & Gaston K.J. 2006. Functional diversity: back to basics and looking
forward. Ecology Letters 9: 741–758.
Pinto, B.C.T., Araujo, F.G. & Hughes, R.M. 2006. Effects of landscape and riparian
condition on a fish index of biotic integrity in a large southeastern Brazil river.
Hydrobiologia 556: 69–83.
Poff, N.L. & Allan, J.D. 1995. Functional organization of stream fish assemblages in
relation to hydrologic variability. Ecology 76: 606–627.
Poff, N.L., Olden, J.D., Vieira, N.K.M., Finn, D.S., Simmons, M.P. & Kondratieff, B.C.
2006. Functional trait niches of North American lotic insects: traits-based
ecological applications in light of phylogenetic relationships. Journal of the North
American Benthological Society 25: 730–755.
Poff, N.L., Pyne, M.I., Bledsoe, B.P., Cuhaciyan, C.C. & Carlisle, D.M. 2010.
Developing linkages between species traits and multiscaled environmental variation
to explore vulnerability of stream benthic communities to climate change. Journal
of the North American Benthological Society 29: 1441–1458.
Power, M.E. 1984. Depth distribution of armored catfish: predator-induced resource
avoidance? Ecology 65: 523–528.
Pusey, B.J.& Arthington, A.H. 2003. Importance of the riparian zone to the
conservation and management of freshwater fish: a review. Marine and Freshwater
Research 54: 1–16.
R Development Core Team. 2005. R: A language and environment for statistical
computing. R Foundation for Statistical Computing, Vienna, Austria.
57
Rabeni, C.F. & Smale, M.A. 1995. Effects on siltation on stream fishes and the
potential mitigating role of the buffering riparian zone. Hydrobiologia 303: 211–
219.
Rezende, C.F., Moraes, M., Manna, L.R., Leitão, R.P., Caramaschi, E.P. & Mazzoni, R.
2010. Mesohabitat indicator species in a coastal stream of the Atlantic rainforest,
Rio de Janeiro-Brazil. Revista de Biología Tropical 58: 1479–1487.
Rocha, F.C., Casatti, L., Carvalho, F.R. & Silva, A.M.S. 2009. Fish assemblages in
stream stretches occupied by cattail (Typhaceae, Angiospermae) stands in Southeast
Brazil. Neotropical Ichthyology 7: 241–250.
Rounick J.S., Winterbourn, M.J., & Lyon, G.L. 1982. Differential utilization of
allochthonous and autochthonous inputs by aquatic invertebrates in some New
Zealand streams: a stable carbon isotope study. Oikos 39: 191–198
Sabino, J. & Zuanon, J. 1998. A stream fish assemblage in Central Amazonia:
distribution, activity patterns and feeding behavior. Ichthyological Exploration of
Freshwaters 8: 201–210.
Sazima, I. 1986. Similarities in feeding behaviour between some marine and freshwater
fishes in two tropical communities. Journal of Fish Biology, 29: 53–65.
Schwartz, J.S. & Herricks, E.E. 2008. Fish use of ecohydraulic-based mesohabitat units
in a low-gradient Illinois stream: implications for stream restoration. Aquatic
Conservation: Marine and Freshwater Ecosystems 18: 852–866.
Silva, A.M., Nalon, M.A., Kronka, F.J.A.N., Álvares, C.A., Camargo, P.B. &
Martinelli, L.A. 2007. Historical land-cover/use in different slope and riparian
buffer zones in watershed of the State of São Paulo, Brazil. Scientia Agricola 64:
325–335.
Sweeney, B.W.; Bott, T.L.; Jackson, J.K.; Kaplan, L.A.; Newbold, J.D.; Standley, L.J.;
Hession, W.C.; Horwitz, R.J. & Colman, M.G. 2004. Riparian Deforestation,
Stream Narrowing, and Loss of Stream Ecosystem Services. Proceedings of the
National Academy of Sciences of the United States of America 101: 14132–14137.
Teresa, F.B. & Casatti, L. 2010. Importância da vegetação ripária em região
intensamente desmatada no sudeste do Brasil: um estudo com peixes de riacho. Pan-
American Journal of Aquatic Sciences 5: 444–453.
Teresa, F.B. & Romero, R.M. 2010. Influence of the riparian zone phytophysiognomies
on the longitudinal distribution of fishes: evidence from a Brazilian savanna stream.
Neotropical Ichthyology 8: 163–170.
58
Tilman, D., Lehman, C.L. & Thomson, K.T. 1997. Plant diversity and ecosystem
productivity: theoretical considerations. Proceedings of Natural Academy of
Sciences of the United States of America 94: 1857–1861.
Villéger, S., Mason, N.W.H. & Mouillot, D. 2008. New multidimensional functional
diversity indices for a multifaceted framework in functional ecology. Ecology 89:
2290–2301.
Villéger, S., Miranda, J.R., Hernáandez, D.F. & Mouillot, D. 2010. Contrasting changes
in taxonomic vs. functional diversity of tropical fish communities after habitat
degradation. Ecological Applications 20: 1512–1522.
Walser, C. A. & Bart, H.L. 1999. Influence of agriculture on instream habitat and fish
assemblage structure in Piedmont watersheds of the Chattahoochee River system.
Ecology of Freshwater Fish 8: 237–246.
Wright, J.P. & Flecker, A.S. 2004. Deforesting the riverscape: the effects of wood on
fish diversity in a Venezuelan piedmont stream. Biological Conservation 120: 439–
447.
59
SUPPORTING INFORMATION
Table A1. Abundance of species captured in forested (F1 to F3) and deforested (D1 to
D3) streams corresponding to at least 3% of the total abundance.
Family Species F1 F2 F3 D1 D2 D3
Anostomidae Leporinus friderici 3 4 0 0 0 0
Characidae Astyanax altiparanae 92 25 30 39 8 27
Astyanax fasciatus 17 32 190 279 0 62
Bryconamericus stramineus 0 0 0 19 0 66
Hemigrammus marginatus 2 0 0 0 0 1
Knodus moenkhausii 12 26 13 0 0 0
Oligosarcus pintoi 0 0 4 7 18 3
Piabina argentea 34 31 41 41 60 77
Planaltina britskii 1 1 0 0 0 0
Serrapinnus notomelas 0 0 0 40 5 1
Crenuchidae Characidium cf. gomesi 0 0 0 44 0 9
Characidium zebra 25 49 87 2 0 9
Parodontidae Parodon nasus 0 1 0 3 1 1
Curimatidae Cyphocharax modestus 0 0 0 0 6 0
Steindachnerina insculpta 0 1 0 2 0 5
Erythrinidae Hoplias malabaricus 0 0 0 10 7 3
Gymnotidae Gymnotus carapo 6 4 0 2 133 86
Sternopygidae Eigenmannia trilineata 9 0 0 0 0 0
Eigenmannia virescens 11 1 0 0 1 0
Sternopygus macrurus 3 0 0 0 0 0
Callichthyidae Aspidoras fuscoguttatus 0 6 23 12 40 0
Callichthys callichthys 0 0 0 3 3 0
Corydoras aeneus 0 1 0 70 51 33
Auchenipteridae Tatia neivai 4 0 0 0 0 0
Heptapteridae Cetopsorhamdia iheringi 7 9 25 0 0 0
Imparfinis mirini 5 2 0 5 0 3
Imparfinis schubarti 1 9 0 57 10 13
Pimellodela avanhandavae 9 6 0 7 0 2
60
Pimelodus maculatus 1 0 0 0 0 0
Pseudopimelodus pulcher 38 34 0 0 0 0
Rhamdia quelen 1 10 19 29 10 1
Loricariidae Hisonotus francirochai 0 0 0 88 19 6
Hypostomus ancistroides 9 17 38 172 105 118
Hypostomus nigromaculatus 61 41 83 39 1 2
Hypostomus sp1 1 0 0 0 0 0
Hypostomus sp2 0 0 0 1 1 1
Poeciliidae Phalloceros harpagos 0 0 0 2 20 5
Poecilia reticulata 0 4 0 102 90 48
Cichlidae Cichlasoma paranaense 1 0 0 0 3 0
Crenicichla britskii 2 2 0 8 30 10
Geophagus brasiliensis 0 0 0 41 0 31
Oreochromis niloticus 0 0 0 0 0 1
Synbranchidae Synbranchus marmoratus 0 0 0 2 2 2
61
Table A2. Functional traits of species with their respective categories. Diet: ins = insectivore, oni = omnivore, per = periphytivore, car =
carnivore, det = detritivore; Feeding tactic: spec = substrate speculation, col = water-column feeding, stea = stealth predation, graz= grazing,
nocpre = nocturnal predator, digg= digger; Habitat use: sur = surface, nec = nektonic, nectob = nectobenthic, mar = margins, bent = benthic;
Preference for water flow: + = fast, +/- = intermediate; – = slow; Preference for stable substrate: + = high, +/– = moderate, – = no preference;
Body size: small = < 50 mm, medium = 51 - 150 mm, large = > 151 mm; tolerance to hypoxia: tol = tolerant, into = intolerant.
Species Diet Feeding
tactic Habitat use
Velocity
preference
Preference for
stable substrate Body size
Hypoxia
tolerance
Aspidoras fuscoguttatus ins espec bent + +/– small tol
Astyanax altiparanae oni col nek – – medium into
Astyanax fasciatus oni col nek – – medium into
Bryconamericus stramineus ins col nek – – small into
Callichthys callychthys oni nocpre bent + – medium tol
Cetopsorhamdia iheringi ins espec bent + + medium into
Characidium gomesi ins stea bent + + medium into
Characidium zebra ins stea bent + + medium into
Cichlasoma paranaense oni digg nectob – – medium tol
Corydoras aeneus oni espec bent +/– – small tol
Crenicichla britskii ins stea marg +/– – medium into
Cyphocharax modestus det digg nectob – – medium into
Eigenmannia trilineata ins stea marg +/– – medium into
62
Eigenmannia virescens ins stea marg +/– – large into
Geophagus brasiliensis oni digg nectob – – large tol
Gymnotus carapo ins stea marg +/– – large tol
Hemigrammus marginatus ins col nek +/– – small tol
Hisonotus francirochai per graz marg +/– – small into
Hoplias malabaricus car stea marg – – large tol
Hypostomus ancistroides per graz bent +/- +/– medium tol
Hypostomus nigromaculatus per graz bent + + medium into
Hypostomus sp1 per graz bent + + medium into
Hypostomus sp2 per graz bent + + medium into
Imparfinis mirini ins espec bent + + medium into
Imparfinis schubarti ins espec bent + +/– medium into
Knodus moenkhausii ins col nek +/– – small tol
Leporinus friderici ins digg nectob +/– +/– medium into
Oligosarcus pintoi ins col marg – – medium into
Oreochromis niloticus oni digg nectob – – large tol
Parodon nasus per graz bent + +/– medium into
Phalloceros harpagos det col sur – – small into
Piabina argentea oni col nek +/– – small into
Planaltina britskii oni col nek +/– – small into
63
Pimellodela avanhandavae ins nocpre bent +/– +/– medium into
Pimelodus maculatus ins nocpre bent + +/– large into
Poecilia reticulata det col sur – – small tol
Pseudopimelodus pulcher ins espec bent + + medium into
Rhamdia quelen ins nocpre bent +/– +/– medium tol
Serrapinnus notomelas det col nek – – small tol
Steindachnerina insculpta det digg nectob – – medium into
Sternopygus macrurus ins stea marg +/– – large into
Synbranchus marmoratus car stea marg +/– – large tol
Tatia neivai ins col marg +/– – medium into
64
Manuscrito será submetido para Journal of Tropical Ecology
- CAPÍTULO III -
Diferenciação funcional de comunidades de peixes em riachos com
diferentes status de conservação
65
RESUMO
As alterações ambientais causadas pelo desmatamento nos riachos podem
representar filtros ambientais, restringindo a ocorrência de espécies a partir das
características funcionais que elas exibem. Testamos essa hipótese em riachos do
sudeste do Brasil, avaliando se riachos em contextos contrastantes com relação à
cobertura florestal podem ser diferenciados a partir da contribuição não aleatória de
grupos funcionais específicos. As amostragens foram realizadas em seis riachos, sendo
três florestados e três desmatados. Foram consideradas cinco características funcionais:
dieta, uso de hábitat, preferência por fluxo, tamanho e tolerância à hipóxia. A
diferenciação dos riachos, em especial entre os riachos florestados e desmatados, foi
maior do que o esperado ao acaso e correspondeu à contribuição diferencial de três
grupos funcionais. Nos riachos desmatados predominaram espécies tolerantes à hipóxia,
que apresentaram combinação variável das outras características. Nos riachos
florestados predominaram espécies de hábitos especializados, bentônicas, com
preferência por alto fluxo de água e dieta insetívora. As alterações nos riachos
provocadas pelo desmatamento relacionadas com a disponibilidade de hábitat, recursos
alimentares e condições físico-químicas (e.g., hipóxia) representam filtros ambientais
que restringem a ocorrência de grupos de espécies especializadas em detrimento do
beneficiamento de espécies tolerantes e oportunistas. Dessa forma, é possível predizer
que a intensificação desse tipo de alteração poderá resultar em modificação do pool
funcional regionalmente.
Palavras-chave: ictiofauna, características funcionais, filtros ambientais, partição da
diversidade
66
INTRODUÇÃO
De acordo com modelos baseados no nicho, a abundância e distribuição das
espécies dependem das suas características funcionais, que por sua vez, determinam a
forma pela qual as espécies são influenciadas pelo ambiente (Townsend & Hildrew,
1994; Goldstein & Meador, 2004; Blanc et al., 2007). Assim, no tempo ecológico, a
composição de espécies em um local é o resultado da seleção não aleatória de espécies a
partir de um pool regional, refletindo a ação de processos baseados na limitação da
similaridade e/ou filtros ambientais (Leibold, 1998; Cornwell et al., 2006).
Comunidades estruturadas por filtros ambientais apresentam espécies com
características mais similares do que esperado ao acaso, resultado da seleção de um
conjunto de espécies com características compatíveis com as condições ambientais
locais (Cornwell et al., 2006). Por outro lado, comunidades estruturadas por processos
de limitação da similaridade apresentam maior dissimilaridade interespecífica do que o
esperado ao acaso, como resultado da exclusão competitiva de espécies similares
(Weiher & Keddy, 1995).
A estruturação de comunidades de riachos é o resultado da ação preponderante
dos filtros ambientais (Jackson et al., 2001; Peres-Neto, 2004). Isso é suportado por
estudos que demonstram a associação das variáveis abióticas com a composição de
espécies (p.ex., Poff, 1997; Goldstein & Meador et al., 2004; 2005). Nesse contexto, as
restrições ambientais manifestam-se em diferentes escalas (Poff, 1997) e, podem ser
geradas também por interferências antrópicas. Por exemplo, Casatti et al. (2009)
compararam a fauna de peixes de dois grupos de riachos em estado de conservação
contrastantes com relação à estrutura física do hábitat e a cobertura florestal, em uma
região submetida à forte interferência antrópica no sudeste do Brasil. Esses dois grupos
de riachos abrigaram comunidades muito distintas, com as espécies tolerantes e
generalistas dominando os riachos degradados e as espécies especialistas e sensíveis
associadas aos riachos em melhor estado de conservação (Casatti et al., 2009). Esses
resultados fornecem evidência de que as restrições na ocorrência das espécies em
riachos variam em função do seu estado de conservação, sugerindo que modificações
antrópicas possam atuar como filtros ambientais (Mayfield et al., 2005).
Se, de fato, as interferências antrópicas modificam as condições nos riachos,
representando novos filtros ambientais dentro de um contexto regional, seria esperado
encontrar convergência funcional, ou seja, comunidades funcionalmente mais similares
67
entre riachos degradados do que o esperado a partir de uma seleção aleatória de espécies
em um pool regional (Cornwell et al., 2006). Nesse sentido, a identificação de grupos de
espécies que compartilham características funcionais responsáveis pela diferenciação
entre riachos em estado de conservação contrastantes seria uma estratégia para a
identificação dos filtros ambientais específicos atuando na organização das
comunidades. Essa abordagem seria especialmente útil para a compreensão dos
processos responsáveis pela estruturação das comunidades em regiões submetidas
impactos antrópicos.
Neste estudo, avaliamos se riachos com diferentes graus de cobertura florestal
na zona ripária podem ser diferenciados a partir da composição funcional das
comunidades de peixes. Realizamos amostragens da ictiofauna em riachos florestados e
desmatados do sudeste do Brasil, em uma região onde a degradação da zona ripária é
uma das principais ameaças à conservação da ictiofauna (Casatti et al., 2006a). O
desmatamento da zona ripária influencia diversos aspectos importantes para peixes em
riachos, tais como disponibilidade de micro-hábitats (Casatti et al., 2009), recursos
alimentares (Ferreira et al., 2011), condições estruturais e fisico-químicas (Murphy &
Hall, 1981). Assim, hipotetizamos que as mudanças causadas pelo desmatamento nos
riachos representam filtros ambientais, resultando na seleção não aleatória de espécies
que compartilham características funcionais associadas à tolerância e à capacidade de
explorar as oportunidades que se tornam disponíveis. Desta forma, predizemos que a
diferenciação funcional entre as comunidades dos riachos florestados e desmatados será
maior do que o padrão obtido aleatoriamente em função da contribuição diferencial de
grupos funcionais específicos.
A estratégia analítica utilizada neste estudo foi adaptada de Pavoine et al.
(2009a), que propuseram a partição da diversidade filogenética ao longo dos ramos da
filogenia com a finalidade de avaliar a importância das linhagens filogenéticas para a
diferenciação das comunidades. De acordo com os próprios autores, essa metodologia
pode ser adaptada para uma abordagem funcional, através da substituição da filogenia
por árvores funcionais. Nesse caso, é possível avaliar se grupos de espécies que
compartilham características funcionais são importantes para a diferenciação de
comunidades de ambientes em diferentes estados de conservação, como é o caso neste
estudo.
68
MATERIAL E MÉTODOS
Área de estudo
O estudo foi realizado na bacia do rio São José dos Dourados, sistema do alto rio
Paraná, sudeste do Brasil (Figura 1). Realizamos amostragens da ictiofauna em trechos
de aproximadamente 300 metros em riachos florestados e desmatados. Os riachos
florestados (n = 3) estão localizados dentro dos maiores fragmentos florestais presentes
na bacia (floresta estacional semi-decidual > 1.600 ha) e os riachos desmatados (n = 3)
não apresentaram vegetação arbórea na zona ripária, estando localizados em uma região
cuja pastagem é a matriz predominante (Figura 1). Não há fragmentos florestais à
montante e jusante dos trechos desmatados. Estes riachos foram selecionados porque
apresentaram características representativas de condições extremas dentro de um
gradiente de cobertura florestal e ao mesmo tempo são similares quanto ao volume de
água e disponibilidade de meso-hábitats (ver capítulo II). Em função da dificuldade de
se encontrar riachos comparáveis aos florestados, os pontos amostrais foram
distribuídos de forma não aleatória na paisagem, com os riachos desmatados
apresentando maior proximidade entre si do que entre os florestados. Entretanto, ambos
os grupos de riachos foram considerados pertencentes ao mesmo pool de espécies
regional (p.ex., Casatti et al. 2006b), já que a distância não explica a variação na
composição de espécies dentro desta bacia (Casatti et al., 2012).
69
Figura. 1. Mapa da área de estudo, destacando o estado de São Paulo (A), a bacia do rio São
José dos Dourados (B), onde estão localizados os riachos florestados (círculos escuros),
inseridos dentro dos maiores fragmentos florestais da bacia (manchas cinza) e os riachos
desmatados (círculos claros) (C).
Amostragem da ictiofauna
Amostramos a ictiofauna em 126 segmentos de cinco metros de extensão cada
(entre 16 e 24 segmentos por riacho), os quais foram distribuídos ao longo de trechos de
cerca de 300 metros em cada riacho. Os segmentos foram escolhidos visando
contemplar a variabilidade de meso-hábitats disponíveis em cada riacho. Antes do início
das coletas, isolamos cada segmento à montante e à jusante com redes. Para as coletas
utilizamos equipamento de pesca elétrica (gerador de corrente alternada, 220 V, 50-
60Hz, 3,4-4,1 A, 1.000 W), passando sucessivamente os eletrodos e puçás ao longo do
segmento. A captura em cada segmento foi finalizada quando nenhum indivíduo foi
capturado após uma passagem. Os animais capturados foram fixados no campo com
formaldeído (10%) e, posteriormente, transferidos para frascos com solução de etanol
(70%) para conservação. Todos os animais encontram-se depositados na coleção
ictiológica (DZSJRP) do Departamento de Zoologia e Botânica da Universidade
Estadual Paulista, São José do Rio Preto, SP, Brasil.
70
Caracterização das unidades amostrais
Após a coleta dos peixes, realizamos a caracterização dos riachos quanto a sua
estrutura física e condições físico-químicas. Essa caracterização foi realizada em cada
segmento (exceto para os parâmetros físico-químicos que foram medidos no ponto
médio de cada trecho), cujos valores foram agregados pela média para cada riacho.
Com relação à estrutura física, medimos a largura superficial, largura interna,
profundidade medida em cinco pontos eqüidistantes de uma margem a outra, fluxo
(medido com um fluxômetro mecânico General Oceanics , modelo 2030) no meio na
coluna d’água, composição do substrato (proporção relativa de cada um dos seus
componentes) (Cummins, 1962) e composição representatividade de estruturas
presentes nos barrancos em cada uma das margens. O substrato foi classificado em
estável, que inclui os tipos de substrato consolidados que permitem a colonização pela
fauna bentônica (argila consolidada, cascalho, seixo e rocha) e instável que inclui areia
e argila inconsolidada, elementos cuja dominância nos leitos dos riachos é indicador de
baixa integridade física (Casatti et al., 2006b). Com relação à vegetação marginal,
foram consideradas qualquer estrutura da vegetação ripária tais como raízes, folhas e
caules oriundos da vegetação terrestre submersos no hábitat interno. Nesse caso, o
volume ocupado pela vegetação submersa foi calculado a partir de medidas da altura e
da largura dos bancos de vegetação em cada transecto. Os parâmetros físico-químicos
avaliados foram oxigênio dissolvido com um oxímetro (LT Lutron, modelo DO 5510),
condutividade, pH, turbidez e temperatura com Horiba, modelo U-10.
Verificamos que os riachos florestados e desmatados diferiram com relação à
proporção de substrato estável e volume de vegetação marginal em contato com a água.
Os riachos florestados apresentaram maior proporção de substrato estável (t = 3,7, p <
0,01), enquanto os riachos desmatados apresentaram maior volume da vegetação
marginal em contato com a água (t = -5,9, p = 0,02). Além disso, os riachos florestados
apresentaram maiores valores de oxigênio dissolvido do que riachos desmatados
(Tabela 1).
71
Tabela 1. Média ± desvio padrão (DP) dos descritores estruturais e físico-químicos dos
riachos florestados (F1 a F3) e desmatados (D1 a D3). Diferenças significativas entre os
grupos de riachos são destacadas em negrito (teste t, p < 0,03).
F1 F2 F3 Média ± DP D1 D2 D3 Média ± DP
Profundidade (m) 0,43 0,31 0,33 0,36 ± 0,06 0,37 0,29 0,39 0,35 ± 0,05
Largura (m) 2,50 3,06 2,53 2,70 ± 0,32 1,74 2,59 2,30 2,21 ± 0,43
Fluxo (cm/s) 0,48 0,42 0,42 0,44 ± 0,03 0,44 0,52 0,44 0,47 ± 0,05
Substrato estável (%) 62,50 46,13 70,19 59,60 ± 12,29 34,57 4,89 10,55 16,60 ± 15,76
Vegetação marginal (dm3) 0,21 0,21 0,17 0,20 ± 0,02 1,06 0,89 0,67 0,87 ± 0,20
OD (mg/l) 8,70 9,80 9,60 9,40 ± 0,59 8,40 8,20 8,30 8,30 ± 0,10
Temperatura (°C) 24,90 20,40 17,00 20,77 ± 3,96 22,30 21,90 21,10 21,77 ± 0,61
Condutividade (mS/cm) 0,07 0,09 0,10 0,09 ± 0,02 0,13 0,10 0,09 0,11 ± 0,02
pH 7,61 8,29 6,10 7,33 ± 1,12 8,35 8,40 8,41 8,39 ± 0,03
Turbidez (NTU) 16 37 8 20,3 ± 14,98 7 10 56 24,3 ± 27,47
Características funcionais
De um total de sete características funcionais pré-selecionadas: dieta, tática
alimentar, uso de hábitat, preferência por fluxo, preferência por substrato estável,
tamanho e tolerância à hipóxia, mantivemos cinco para as análises neste estudo (Tabela
A1, Apêndice). Duas características: preferência por substrato estável e tática alimentar
foram removidas porque foram significativamente correlacionadas com a preferência
por fluxo e o uso de hábitat, respectivamente. O grau de redundância das características
foi avaliado por meio da análise de correlação múltipla de matrizes (MRM), utilizando-
se o coeficiente de correlação de Pearson. A preferência por fluxo e uso de hábitat foram
mantidas porque a caracterização das espécies baseou-se nos dados obtidos nos riachos
avaliados neste estudo, diferentemente das características removidas que foram baseadas
em informações oriundas da literatura.
As cinco características selecionadas correspondem a aspectos ecológicos e de
história de vida das espécies e estão diretamente relacionadas com processos
ecossistêmicos que podem ser direta ou indiretamente influenciados pelo
desmatamento: (i) disponibilidade de micro-hábitats e sedimentação (uso de hábitat e
preferência por fluxo); (ii) produtividade e enriquecimento orgânico (dieta e tamanho);
(iii) disponibilidade e uso de recursos, dinâmica trófica e ciclagem de nutrientes (dieta,
72
uso de hábitat e tamanho); (iv) qualidade da água (tolerância à hipóxia) (Murphy &
Hall, 1981; Smale & Rabeni, 1995; Holmlund & Halmer, 1999; Bunn et al., 1999;
Esteves et al., 2008; Lorion & Kennedy, 2009).
Dieta
A dieta foi avaliada a partir da análise do conteúdo estomacal das espécies
coletadas. Os estômagos de 10 exemplares de cada espécie (exceto para aquelas que
apresentaram abundância inferior a 10 indivíduos) foram removidos e abertos para
inspeção e identificação do conteúdo estomacal. Os itens alimentares foram
identificados e classificados em grandes categorias (detrito, insetos, peixes, perifíton,
protista, fragmentos de vegetais superiores e algas). A composição porcentual de cada
item em cada estômago foi estimada visualmente para a identificação do item
dominante cuja freqüência foi combinada com a porcentagem de ocorrência do item nos
estômagos de cada espécie para a obtenção do índice de importância alimentar (IAi,
Bennemann et al., 2006):
ni
iOD
1
Onde Oi é o número proporcional de estômagos em que o item i ocorreu, e D a
proporção de estômagos em que o item i foi dominante. Os valores presentes na matriz
de dieta (variável fuzzy) correspondem aos valores de IAi.
Uso de hábitat
As espécies foram classificadas quanto ao uso do espaço (bentônico,
nectobentônico, nectônico, superfície e margem), atribuindo-se o valor “1” aos locais
que são eventualmente ocupados pela espécie e “2” para os locais onde a espécie ocorre
preferencialmente (variável fuzzy). As informações sobre o uso de hábitat foram obtidas
a partir de observações subaquáticas realizadas nos seis riachos deste estudo em um
período que antecedeu a coleta da ictiofauna, totalizando 1.200 horas de observações
subaquáticas. Complementarmente, foram utilizadas informações da literatura (Casatti
& Castro, 1998; Casatti et al., 2001; Casatti et al.,2005; Ceneviva-Bastos et al., 2010).
Preferência por fluxo
73
A preferência por fluxo (variável fuzzy) foi medida a partir do cálculo do índice
de preferência (IP) calculado para cada espécie em cada uma das quatro classes de fluxo
(F1 = 0,0 a 0,25 cm/s; F2 = 0,26 a 0,45 cm/s; F3 = 0,46 a 0,65 cm/s; F4 > 0,66 cm/s)
baseado na metodologia proposta por Freeman et al. (1997). O valor final de IP foi
padronizado para variar de 0 a 1 dividindo-se os valores de IP de cada espécie em cada
classe pelo valor máximo de IP encontrado para a espécie:
Onde, nsi é o número de indivíduos da espécie na classe i, N é o número total de
indivíduos da espécie coletados, nhi é o número de unidades de hábitats amostrados na
classe i, H é o número total de unidades de hábitats e max refere-se ao valor máximo.
Como a preferência por condições de fluxo são consistentes entre riachos florestados e
desmatados (Capítulo I) o valor final de IP para as espécies que ocorreram em ambos os
grupos de riachos foi representado pela média dos valores obtidos em cada grupo.
Tamanho e tolerância à hipóxia
Os maiores indivíduos de cada espécie amostrada tiveram seu comprimento
padrão mensurados (variável quantitativa) com o auxílio de um paquímetro. A
caracterização das espécies quanto à tolerância à hipóxia (variável nominal) baseou-se
em informações da literatura. Foram consideradas tolerantes as espécies que são
comumente encontradas em trechos com baixos níveis de oxigênio (Kramer &
Mehegan, 1981; Araujo & Garutti, 2003; Bozzetti & Schulz, 2004; Casatti et al., 2006b;
Rocha et al., 2009; Casatti et al., 2012).
Análise dos dados
Utilizamos para as análises somente espécies que ocorreram em mais do que 3%
dos segmentos amostrados, visando reduzir a possível influência de estimativas
equivocadas de dieta e preferência por fluxo baseadas em poucos indivíduos. Como o
número de segmentos amostrados diferiu entre os riachos, a matriz de espécies final foi
composta pela abundância média das espécies considerando os segmentos amostrados
em cada riacho.
O dendrograma funcional foi construído a partir de uma matriz de distância
representando a similaridade entre as espécies baseada nas cinco características
funcionais selecionadas para este estudo. Essa matriz foi obtida a partir de uma
74
generalização da aplicação do coeficiente de Gower proposta por Pavoine et al. (2009b)
que possibilita a combinação de características cujos dados são de natureza estatística
diferentes (e.g., fuzzy, nominal e quantitativa). A distância entre cada par de espécies
representa a média da distância entre as duas espécies, considerando cada característica
separadamente, após a sua padronização (variação de 0 a 1). O dendrograma funcional
final foi construído utilizando o método das médias aritméticas não ponderadas
(UPGMA).
A diversidade funcional foi calculada a partir do índice de entropia quadrática de
Rao (QE). Este índice é baseado na soma das distâncias entre os pares de espécies de
uma comunidade considerando sua abundância (Botta-Dukát, 2005) Quando QE é
aplicado para medir diversidade funcional a partir de um dendrograma, a distância entre
duas espécies é igual à soma do comprimento do braço do dendrograma que conecta as
duas espécies a um ancestral comum (Pavoine et al. 2005). Adaptando a metodologia de
Pavoine et al. (2009a) para lidar com árvore funcional ao invés da árvore filogenética,
calculamos o valor de QE para cada partição do dendrograma. Nesse caso, os grupos
funcionais que se originam em cada partição equivalem aos ramos terminais e a
abundância corresponde à soma das abundâncias das espécies que compõem o grupo.
Para avaliar as diferenças funcionais entre as comunidades (diversidade funcional β) e a
contribuição dos grupos funcionais formados nas diferentes partições do dendrograma,
realizamos a partição de QE nos componentes γ, α e β, assumindo que a diversidade
funcional regional (γ) corresponde à média da diversidade de cada riacho (α) somado à
diversidade entre os riachos (β) (Lande, 1996). Os valores de β/γ (i.e., variabilidade
entre riachos) observados para cada partição do dendrograma foram comparados com
valores simulados a partir da aleatorização dos ramos terminais, de forma que a conexão
entre a estrutura funcional do dendrograma e a abundância dos ramos terminais fosse
permutada (999 permutações) (Pavoine et al., 2009a). Valores de β/γ maiores do que o
esperado ao acaso indicam agregação funcional, ou seja, maior semelhança funcional
dentro das comunidades do que entre comunidades, resultado esperado em comunidades
organizadas por filtros ambientais (Pavoine & Bonsall, 2011).
Complementarmente, aplicamos a análise de coordenadas principais dupla
(DPCoA) (Pavoine et al., 2004) para identificar os grupos funcionais que foram mais
importantes na diferenciação funcional entre os riachos florestados e desmatados. A
DPCoA possibilita a ordenação das espécies, cujas posições relativas são definidas
pelas suas dissimilaridades funcionais, assim como das amostras cujas posições
75
relativas refletem as suas diferenças baseadas no QE (Pavoine et al., 2004). Todas as
análises foram realizadas no software R (R Development Core Team 2009), utilizando
os pacotes “ade4” (Dray & Dufour, 2007) e as funções “aptree”, “abg.aptree” e
“rtest.aptree” fornecidas por Pavoine et al. (2009a).
RESULTADOS
Foram amostradas 43 espécies (Tabela 2) das quais oito foram removidas das
análises por apresentarem baixa ocorrência, resultando em uma matriz final composta
por 35 espécies, que representaram 99,5% da abundância total. Essas espécies foram
agrupadas de acordo com suas semelhanças funcionais em um dendrograma composto
por 34 partições (Figura 2). A diferenciação dos riachos (β/γ) baseado na composição
funcional das comunidades foi maior do que o esperado ao acaso nas partições #30 e
#33 (p < 0,02). Essas partições separaram as espécies em três e seis grupos funcionais,
respectivamente (Figura 3). O primeiro eixo da DPCoA resumiu a maior parte da
variação referente às diferenças na composição funcional das comunidades (72,6%) e
permitiu destacar os grupos funcionais foram mais importantes para a diferenciação dos
riachos florestados e desmatados (Figura 4). Os três grupos originados na partição #33
foram importantes para a dissimilaridade funcional entre os riachos florestados e
desmatados, sendo que os grupos 1 e 2 foram mais associados aos riachos desmatados,
enquanto as espécies do grupo 3 foram mais abundantes nos riachos florestados
(Figuras 4). As espécies dos grupos 1 e 2 compartilham a tolerância à hipóxia e
complementarmente, as espécies do grupo 1 podem ser caracterizadas ainda pelo grande
porte, dieta carnívora/insetívora e por utilizarem as margens como hábitat. As espécies
do grupo 2 podem ser caracterizadas pelo pequeno porte, uso da superfície com dieta
detritívora (Poecilia reticulata e Phalloceros harpagos) e pelo hábitat bentônico com
preferência por fluxo de água intermediário e dieta insetívora (demais espécies) (Figura
3). O grupo 3 foi composto por espécies intolerantes à hipóxia e com variações nas
outras características (Figura 3).
76
Tabela 2. Abundância média das espécies (com seus respectivos códigos) que
ocorreram em mais do que 3% dos segmentos nos riachos florestados (F1 a F3) e
desmatados (D1 a D3), bacia do rio São José dos Dourados, Alto Paraná.
Família Espécies Códigos F1 F2 F3 D1 D2 D3
Anostomidae Leporinus friderici Lepfri 0.2 0.2 0.0 0.0 0.0 0.0
Characidae Astyanax altiparanae Astalt 4.0 1.2 1.6 1.6 0.4 1.7
Astyanax fasciatus Astfas 0.8 1.7 10.1 11.6 0.0 3.9
Bryconamericus stramineus Brystr 0.0 0.0 0.0 0.8 0.0 4.2
Knodus moenkhausii Knomoe 0.7 1.4 0.7 0.0 0.0 0.0
Oligosarcus pintoi Olipin 0.0 0.0 0.2 0.3 0.8 0.2
Piabina argentea Piaarg 1.5 1.6 2.2 1.7 2.7 4.9
Serrapinnus notomelas Sernot 0.0 0.0 0.0 1.7 0.2 0.1
Crenuchidae Characidium cf. gomesi Chagom 0.0 0.0 0.0 1.8 0.0 0.6
Characidium zebra Chazeb 1.3 3.1 4.7 0.1 0.0 0.6
Parodontidae Parodon nasus Parnas 0.0 0.0 0.0 0.1 0.0 0.1
Curimatidae Steindachnerina insculpta Steins 0.0 0.0 0.0 0.1 0.0 0.3
Erythrinidae Hoplias malabaricus Hopmal 0.0 0.0 0.0 0.4 0.3 0.2
Gymnotidae Gymnotus carapo Gymcar 0.3 0.2 0.0 0.1 6.1 5.4
Sternopygidae Eigenmannia trilineata Eigtri 0.4 0.0 0.0 0.0 0.0 0.0
Eigenmannia virescens Eigvir 0.5 0.0 0.0 0.0 0.0 0.0
Sternopygus macrurus Stemac 0.2 0.0 0.0 0.0 0.0 0.0
Callichthyidae Aspidoras fuscoguttatus Aspfus 0.0 0.5 1.2 0.5 1.9 0.0
Callichthys callichthys Calcal 0.0 0.0 0.0 0.1 0.1 0.0
Corydoras aeneus Coraen 0.0 0.0 0.0 2.9 2.3 2.1
Auchenipteridae Tatia neivai Tatnei 0.3 0.0 0.0 0.0 0.0 0.0
Heptapteridae Cetopsorhamdia iheringi Cetihe 0.3 0.4 1.3 0.0 0.0 0.0
Imparfinis mirini Impmir 0.2 0.1 0.0 0.2 0.0 0.2
Imparfinis schubarti Impsch 0.0 0.7 0.0 2.4 0.5 0.8
Pimellodela avanhandavae Pimava 0.5 0.3 0.0 0.3 0.0 0.1
Pseudopimelodus pulcher Psepul 1.8 1.8 0.0 0.0 0.0 0.0
Rhamdia quelen Rhaque 0.0 0.5 1.0 1.2 0.5 0.1
Loricariidae Hisonotus francirochai Hisfra 0.0 0.0 0.0 3.7 0.9 0.4
Hypostomus ancistroides Hypanc 0.4 1.1 2.0 7.2 5.0 7.7
77
Hypostomus nigromaculatus Hypnig 2.9 2.2 4.6 1.6 0.0 0.1
Poeciliidae Phalloceros harpagos Phahar 0.0 0.0 0.0 0.1 0.9 0.3
Poecilia reticulata Poeret 0.0 0.2 0.0 4.3 4.1 3.0
Cichlidae Crenicichla britskii Crebri 0.1 0.1 0.0 0.3 1.5 0.6
Geophagus brasiliensis Geobra 0.0 0.0 0.0 1.7 0.0 1.9
Synbranchidae Synbranchus marmoratus Synmar 0.0 0.0 0.0 0.1 0.1 0.1
78
Figura 2. Dendrograma de similaridade funcional das 35 espécies de peixes incluídas nas análises. A abundância média das espécies dentro dos
riachos florestados (F1-F3) e desmatados (D1-D3) são representadas pelo tamanho dos quadrados pretos. Os nomes das espécies encontram-se
abreviados de acordo com a tabela 2.
79
Três dos seis grupos funcionais que se originaram na partição #30 (grupos 1, 2 e
3 ) foram importantes na diferenciação funcional entre os riachos florestados e
desmatados (Figura 4). Os grupos 1 e 2 foram os mesmos originados na partição #33. O
grupo 3 foi composto por espécies de hábitos especializados, bentônicas, com
preferência por alto fluxo de água e dieta insetívora (Figura 3). Os grupos 1 e 2 foram
mais associados aos riachos desmatados e o grupo 3 concentrou-se nos riachos
florestados (Figura 4). O grupo 4, formado por espécies que vivem associadas às
margens de dieta insetívora, e o “grupo” 5, representado por uma única espécie
(Steindacherina insculpta), foram pouco abundantes e conseqüentemente tiveram
pequena contribuição para a diferenciação funcional dos riachos (Figura 2 e 4). As
espécies do grupo 6 que compartilham o hábito nectônico, dieta generalista/insetívora
com preferência por ambientes lênticos tiveram menor importância na diferenciação
funcional dos riachos por serem homogeneamente distribuídas entre as amostras (Figura
2 e 4).
80
Figura 3. Identificação dos grupos funcionais definidos pelas partições #30 e #33 com suas respectivas tendências de mudança na
representatividade em relação ao desmatamento. As características dominantes dos grupos funcionais que foram importantes para a diferenciação
dos riachos de acordo com a DPCoA (Figura 4) também são mostradas. Os nomes das espécies encontram-se abreviados de acordo com a tabela
2.
81
Figura 4. Representação gráfica da análise de coordenadas principais dupla,
descrevendo as diferenças funcionais entre os riachos florestados (F1-F3) e desmatados
(D1-D3). Os círculos pretos correspondem às espécies cujas posições relativas refletem
suas dissimilaridades funcionais. As espécies estão conectadas de acordo com os grupos
funcionais originados na partição #33 (a) e #30 (b). A identidade das espécies e suas
respectivas características funcionais são mostradas na figura 3.
DISCUSSÃO
Neste estudo, mostramos que os riachos submetidos ao desmatamento das suas
zonas ripárias diferem de riachos florestados com relação à composição de espécies
como resultado da restrição à ocorrência de grupos de espécies baseado nas suas
características funcionais. A adaptação da estratégia analítica desenvolvida por Pavoine
et al. (2009a) para lidar com uma árvore funcional em substituição à filogenia, permitiu
diferenciar os riachos desmatados e florestados a partir da estrutura funcional das
comunidades. A partir da decomposição da diversidade ao longo da árvore funcional
constatamos maiores valores de diversidade funcional beta do que o esperado ao acaso,
indicando que as comunidades dentro de cada riacho apresentam dominância de um ou
poucos grupos funcionais em relação ao esperado a partir do pool regional (alta β e
baixa α) (Pavoine et al., 2005). Os maiores valores de diversidade beta ocorreram nas
partições basais do dendrograma de similaridade funcional (#30 e #33) e foram
conseqüência da contribuição diferencial dos grupos funcionais 1, 2 e 3 entre os riachos,
sugerindo a ocorrência de convergência funcional dentro das comunidades, o que é
compatível com a influência de filtros ambientais (Cornwell et al., 2006; Mouillot et al.,
82
2007). As espécies de comunidades estruturadas por filtros ambientais tendem a ser
mais similares entre si devido às restrições do ambiente que favorece aquelas que
compartilham características compatíveis com as condições locais. A diferenciação
funcional foi mais acentuada entre os riachos florestados e desmatados, refletindo as
diferenças nas condições ambientais que, por sua vez, representam filtros ambientais
importantes para a estruturação das comunidades.
O beneficiamento de grupos funcionais em riachos desmatados, representado
neste estudo pelos grupos 1 e 2, também foi descrito por Bojsen & Barriga (2002) e
Lorion & Kennedy (2009). Estes autores mostraram um aumento da abundância e
biomassa de peixes herbívoros e detritívoros em riachos desmatados como
conseqüência da redução do sombreamento e aumento da produtividade primária. Por
outro lado, Casatti et al. (2006a,b) e Teresa & Casatti (2010) atribuíram o declínio de
espécies reofílicas à simplificação do hábitat e ao desmatamento, o que no presente
estudo corresponde ao declínio verificado no grupo 3 (#30). Esses padrões são
coerentes com a generalização de que a degradação ambiental resulta na ampliação da
distribuição de espécies oportunistas e generalistas em substituição às espécies sensíveis
e especialistas (Rahel, 2002; Devictor et al., 2008).
A associação dos grupos funcionais tolerantes à hipóxia (grupos 1 e 2) aos
riachos desmatados pode estar relacionada com as alterações em alguns processos
ecológicos. A perda da cobertura ripária e, consequentemente, o aumento da incidência
de radiação solar proporciona o aumento da produtividade primária em riachos (Murphy
& Hall, 1981; Gücker et al., 2009; Lorion & Kennedy, 2009). Isso foi evidenciado neste
estudo de maneira indireta, pela grande abundância de macrófitas e algas nos riachos
desmatados. O aumento da produtividade, associado à disponibilidade de nutrientes,
como é típico de áreas rurais (Brodie & Mitchell, 2005), pode resultar no aumento das
taxas de decomposição e, consequentemente redução da disponibilidade de oxigênio
dissolvido. Nas margens, esse efeito pode ser potencializado, uma vez que a vegetação
herbácea abundante nos barrancos desses riachos (Casatti et al., 2009; Rocha et al.,
2009) estende-se para o hábitat interno, contribuindo para a atenuação do fluxo de água
e interceptação de detritos, condições propícias para a aceleração dos processo de
decomposição (i.e., consumo de oxigênio) (Bunn et al., 1999). Assim, a suposta maior
disponibilidade de micro-hábitats com baixa disponibilidade de oxigênio nos riachos
desmatados pode estar associado com o favorecimento das espécies tolerantes à hipóxia
dos grupos 1 e 2 e declínio das espécies intolerantes do grupo 3 de ambas partições.
83
Além da tolerância à hipóxia, as condições encontradas nos riachos desmatados
podem filtrar outras características. Por exemplo, as espécies do grupo 1 caracterizam-
se pelo grande porte (> 180 mm) e inclui as únicas espécies piscívoras registradas neste
estudo (Hoplias malabaricus e Synbranchus marmoratus). O favorecimento de
predadores em riachos desmatados pode estar associado ao aumento da produtividade
primária que aumentaria a energia disponível no sistema, possibilitando a manutenção
de peixes maiores de níveis tróficos superiores (Brocksen et al., 1968; Murphy & Hall,
1981). Além disso, conforme discutido no capítulo II, a vegetação marginal abundante
pode conferir proteção para peixes de grande porte contra predadores terrestres.
O grupo funcional 2 inclui espécies que apresentam características variáveis
quanto ao uso de hábitat, preferência por fluxo e dieta e podem ser divididas em dois
subgrupos. Em um deles encontram-se Poecilia reticulata e Phalloceros harpagos que
apresentam em comum o uso da superfície, tamanho pequeno e detrito como um
importante item alimentar. Essas características favorecem a exploração de micro-
hábitats rasos e vegetados que são abundantes em riachos desmatados (Casatti et al.,
2009). As espécies pertencentes ao outro subgrupo dentro do grupo 2 são insetívoras,
bentônicas e generalistas quanto à preferência por fluxo. Essas espécies podem ser
consideradas equivalentes funcionais tolerantes à hipóxia do grupo 3 (#30), uma vez
que compartilham dieta, tamanho e uso de hábitat. Nesse sentido, o declínio do grupo 3
e favorecimento do grupo 2 nos riachos desmatados segue a tendência constatada para
peixes (Rahel, 2002; Casatti et al., 2009; Casatti et al., 2012) e outros grupos animais
(Warren et al., 2001; Devictor et al., 2008), de substituição de espécies especialistas por
generalistas tolerantes em ambientes antropicamente degradados. A redundância
funcional das espécies bentônicas do grupo 2 e as do grupo 3 (#30) quanto à dieta e uso
de hábitat e a sua complementaridade quanto à tolerância à hipóxia podem contribuir
para a manutenção de processos relacionados ao uso de recursos nos riachos
desmatados, mesmo com a mudança na composição de espécies. Isso acontece quando
as espécies compartilham características que exercem algum efeito sobre o ecossistema
(e.g., dieta baseada em insetos bentônicos), mas diferem na resposta às variações nas
condições do ambiente (e.g., tolerância à hipóxia) (Diaz & Cabido, 2001). Embora seja
plausível supor que esses grupos contribuem para processos como ciclagem de
nutrientes, dinâmica trófica e produtividade secundária (Taylor et al., 2006; Small et al.,
2011) nos dois grupos de riachos, não é possível inferir a magnitude e equivalência (ou
divergência) das suas contribuições sem a utilização de abordagens experimentais.
84
Três grupos funcionais não foram importantes para a diferenciação funcional
dos riachos florestados e desmatados, por conterem espécies raras e pouco abundantes
(grupos 4 e 5) e por não apresentarem ocorrência influenciada pela cobertura florestal
(grupo 6). As espécies do grupo 6 compartilham o hábito nectônico, dieta
insetívora/onívora e preferência por meso-hábitats com baixo fluxo de água (Capítulo I
e II) e são representadas, na sua maioria por pequenos caracídeos. Trata-se de um grupo
dominante em riachos Neotropicais (Castro, 1999), tendo sido bem representados nos
riachos deste estudo, exceto em D2, provavelmente em função de particularidades
presentes neste riacho como maior fluxo de água e menor disponibilidade de espaço na
coluna d’água livre de macrófitas no meio do canal. Essas espécies são conhecidas pela
sua plasticidade trófica e oportunismo, sendo capazes de explorar com sucesso os
recursos disponíveis mesmo em ambientes degradados (Melo, 2004; Ceneviva-Bastos et
al., 2010; Ferreira et al., 2012), o que pode explicar a sua baixa sensibilidade ao
desmatamento.
A abordagem utilizada neste estudo propiciou a identificação de combinações de
características necessárias para a formação de grupos de espécies que apresentam
respostas similares às variações no grau de cobertura florestal. Nesse sentido, os
resultados proporcionam suporte para a idéia de que a classificação de espécies em
grupos funcionais representa uma alternativa viável para a predição da integridade de
ambientes lóticos (Welcomme et al., 2006). Adicionalmente, a incorporação de aspectos
funcionais das comunidades torna essa abordagem útil também para a predição dos
efeitos das alterações antrópicas sobre processos e serviços ecossistêmicos (Dainly et
al., 2006).
Os resultados corroboram a hipótese e indicam que as variações na composição
funcional das comunidades dependem do estado de conservação dos riachos com
relação à cobertura florestal das suas zonas ripárias. Especificamente, a influência das
variações na cobertura florestal sobre a disponibilidade de hábitat e recursos
alimentares, assim como nas condições físico-químicas (e.g., hipóxia) parecem ser
importantes na organização das comunidades. Ao demonstrar que as modificações
ambientais causadas pelo desmatamento resultam na restrição de grupos de espécies que
compartilham características funcionais específicas, antecipamos que intensificação
dessa alteração poderá resultar além da homogeneização da fauna como sugerido por
Casatti et al. (2009), na limitação do pool funcional.
85
REFERÊNCIAS BIBLIOGRÁFICAS
Araujo, R.B. & Garutti, V. 2003. Ecology of a stream from upper Paraná river basin
inhabited by Aspidoras fuscoguttatus Nijssen and Isbrücker, 1976 (Siluriformes,
Callichthyidae). Brazilian Journal of Biology, 63: 363–372.
Bennemann, S.T., Casatti, L. & Oliveira, D.C. 2006. Alimentação de peixes: proposta
para análise de itens registrados em conteúdos gástricos. Biota Neotropica, 6(2): 1-
8.
Blanck, A., Tedesco, P.A. & Lamouroux, N. 2007. Relationships between life-history
strategies of European freshwater fish species and their habitat preferences
Freshwater Biology, 52: 843–859.
Botta-Dukát, Z. 2005. Rao’s quadratic entropy as a measure of functional diversity
based on multiple traits. Journal of Vegetation Science, 16: 533–540.
Bozzetti, M. & Schulz, U.H. 2004. An index of biotic integrity based on fish
assemblages for subtropical streams in southern Brazil. Hydrobiologia, 529: 133–
144.
Brocksen, R.W., Davis, G.E. & Warren, C.E. 1968. Competition, Food consumption,
and production of sculpins and trout in laboratory stream communities. The Journal
of Wildlife Management, 32(1): 51–75.
Brodie, J.E. & Mitchell, A.W. 2005. Nutrients in Australian tropical rivers: changes
with agricultural development and implications for receiving environments Marine
and Freshwater Research, 56(3): 279–302.
Bunn, S.E., Davies, P.M. & Mosisch, T.D. 1999. Ecosystem measures of river health
and their response to riparian and catchment degradation. Freshwater Biology, 41:
333–345.
Casatti, L. & Castro, R.M.C. 1998. A fish community of the São Francisco river
headwater riffles, southeastern Brazil. Ichthyological Exploration of Freshwaters,
9: 229–242.
Casatti, L., Langeani, F. & Castro, R.M.C. 2001. Peixes de riacho do Parque Estadual
Morro do Diabo, bacia do alto rio Paraná, SP. Biota Neotropica, 1: 1–15.
Casatti, L., Rocha, F.C. & Pereira, D.C. 2005. Habitat use by two species of
Hypostomus (Pisces, Loricariidae) in southeastern brazilian streams. Biota
Neotropica, 5: 157–165.
86
Casatti, L., Langeani, F. & Ferreira, C.P. 2006a. Effects of physical habitat degradation
on the stream fish assemblage structure in a Pasture Region. Environmental
Management, 38:974–982.
Casatti, L., Langeanim F., Silva, A.M. & Castro, R.M.C. 2006b. Stream fish, water and
habitat quality in a pasture dominated basin, southeastern Brazil. Brazilian Journal
of Biology, 66: 681–696.
Casatti, L., Ferreira, C.P. & Carvalho, F.R. 2009. Grass-dominated stream sites exhibit
low fish species diversity and dominance by guppies: an assessment of two tropical
pasture river basins. Hydrobiologia, 623: 173–189.
Casatti, L., Teresa, F.B., Gonçalves-Souza, T., Bessa, E., Manzotti, A.R., Gonçalves,
C.S. & Zeni, J.O. 2012. From forests to cattail: how does the riparian zone
influence stream fish? Neotropical Ichthyology. In press.
Castro, R.M.C. 1999. Evolução da ictiofauna de riachos sul-americanos: padrões gerais
e possíveis processos causais. In: Caramaschi, E.P., Mazzoni, R., Bizerril,
C.R.S.F., Peres-Neto, P.R., eds. Ecologia de Peixes de Riachos: Estado Atual e
Perspectivas. Rio de Janeiro: Oecologia Brasiliensis, v. VI, pp. 139–155.
Ceneviva-Bastos, M., Casatti, L. & Rossa-Feres, D.C. 2010. Meso and microhabitat
analysis and feeding habits of small nektonic characins (Teleostei: Characiformes)
in Neotropical streams. Zoologia, 27: 191–200.
Cornwell, W.K., Schwilk, D.W. & Ackerly, D.D. 2006. A trait-based test for habitat
filtering: convex hull volume. Ecology, 87(6): 1465–1471.
Costa, P.F. & Schulz, U.H. 2010. The fish community as an indicator of biotic integrity
of the streams in the Sinos River basin, Brazil. Brazilian Journal of Biology, 70(4):
1195–1205.
Cummins, K.H. 1962. An evaluation of some techniques for the collection and analysis
of benthic samples with special emphasis on lotic waters. The American Midland
Naturalist, 67: 477–504.
Daily, G., Alexander, S., Ehrlich, P.R., Gouler, L., Lubchenco, J., Matson, P.A.,
Mooney, H.A., Postel, S., Schneider, S.H., Tilman, D. & Woodwell, G.M.1997.
Ecosystem services: benefits supplied to human societies by natural ecosystems.
Issues in Ecology, 2: 1–18.Devictor, V., Julliard, R., Jiguet, J.C.F., Lee, A. &
Couvet, D. 2008. Functional biotic homogenization of bird communities in
disturbed landscapes. Global Ecology and Biogeography, 17: 252–261.
87
Diaz, S. & Cabido, M. 2001. Vive la différence: plant functional diversity matters to
ecosystem processes. Trends in Ecology and Evolution, 16(8): 646–655.
Dray, S. & Dufour, A.B. 2007. The ade4 package: implementing the duality diagram for
ecologists. Journal of Statistical Software, 22: 1–20.
Ernst, R., Linsenmair, K.E. & Rödel, M.O. 2006. Diversity erosion beyond the species
level: dramatic loss of functional diversity after selective logging in two tropical
amphibian communities. Biological Conservation, 133: 143–155.
Esteves, K.E., Lobo, A.V.P. & Faria, M.D.R. 2008. Trophic structure of a fish
community along environmental gradients of a subtropical river (Paraitinga River,
Upper Tietê River Basin, Brazil). Hydrobiologia, 598: 373–387.
Ferreira, A., Paula, F.R., Ferraz, S.F.B., Gerhard, P., Kashiwaqui, E.A.L., Cyrino, J.E.P.
& Martinelli, L.A. 2012. Riparian coverage affects diets of characids in neotropical
streams. Ecology of Freshwater Fish, 21(1): 12–22.
Freeman, M.C., Bowen, Z.H. & Crance, J.H. 1997. Transferability of habitat suitability
criteria for fishes in warmwater streams. North American Journal of Fisheries
Management, 17: 20–31.
Goldstein, R.M. & Meador, M.R. 2004. Comparisons of fish species traits from small
streams to large rivers. Transactions of the American Fisheries Society, 133: 971–
983.
Gücker, B., Boëchat, I.G. & Giani, A. 2009. Impacts of agricultural land use on
ecosystem structure and whole-stream metabolism of tropical Cerrado streams.
Freshwater Biology, 54: 2069–2085.
Holmlund, C.M. & Hammer, M. 1999. Ecosystem services generated by fish
populations. Ecological Economics, 29: 253–268.
Kawaguchi, Y. & Nakano, S. 2001. Contribution of terrestrial invertebrates to the
annual resource budget for salmonids in forest and grassland reaches of a
headwater stream. Freshwater Biology, 46: 303–316.
Kramer, D.L. & Mehegan, J.P. 1981. Aquatic surface respiration, an adaptative
response to hypoxia in the guppy, Poecilia reticulata (Pisces, Poeciliidae).
Environmental Biology of Fishes, 6: 299–313.
Lande, R. 1996. Statistics and partitioning of species diversity, and similarity among
multiple communities. Oikos, 76: 5–13.
Leibold, M. 1998. Similarity and local coexistence of species in regional biotas.
Evolution & Ecology, 12(1): 95–100.
88
Lorion, C.M. & Kennedy, B. P. 2009. Riparian forest buffers mitigate the effects of
deforestation on fish assemblages in tropical headwater streams. Ecological
Applications, 19: 468–479.
Melo, C.E., Machado, F.A. & Pinto-Silva, V. 2004. Feeding habits of fish from a stream
in the savanna of Central Brazil, Araguaia Basin. Neotropical Ichthyology, 2: 37–
44.
Mouillot, D., Dumay, O. & Tomasini, J.A. 2007. Limiting similarity, niche filtering and
functional diversity in coastal lagoon fish communities. Estuarine, Coastal and
Shelf Science, 71: 443–456.
Mouillot, D., Villéger, S., Scherer-Lorenzen, M. & Mason, N.W.H. 2011. Functional
structure of biological communities predicts ecosystem multifunctionality. PLoS
One, 6: e17476.
Murphy, M.L., & Hall, J.D. 1981. Varied effects of clear-cut logging on predators and
their habitat in small streams of the Cascade Mountains, Oregon. Canadian Journal
of Fisheries and Aquatic Sciences, 38: 137–145.
Pavoine, S. & Bonsall, M.B. 2011. Measuring biodiversity to explain community
assembly: a unified approach. Biological Reviews, 86: 702-812.
Pavoine, S., Dufour, A.B. & Chessel, D. 2004. From dissimilarities among species to
dissimilarities among communities: a double principal coordinate analysis. Journal
of Theoretical Biology, 228: 523–537.
Pavoine, S., Ollier, S. & Pontier, D. 2005. Measuring diversity from dissimilarities with
Rao’s quadratic entropy: are any dissimilarity indices suitable? Theoretical
Population Biology, 67: 231–239.
Pavoine, S., Love, M. & Bonsall, M. B. 2009a. Hierarchical partitioning of evolutionary
and ecological patterns in the organization of phylogenetically structured species
assemblages: application to rockfish (genus: Sebastes) in the Southern California
Bight. Ecology Letters, 12: 898–908.
Pavoine, S., Vallet, J., Dufour, A.-B., Gachet, S. & Daniel, H. 2009b. On the challenge
of treating various types of variables: application for improving the measurement
of functional diversity. Oikos, 118: 391–402.
Petchey, O.L. & Gaston K.J. 2006. Functional diversity: back to basics and looking
forward. Ecology Letters, 9: 741–758.
89
Poff, N.L. 1997. Landscape filters and species traits: towards mechanistic understanding
and prediction in stream ecology. Journal of the North American Benthological
Society, 16(2): 391–409.
R Development Core Team. 2009. R: A language and environment for statistical
computing. R Foundation for Statistical Computing, Vienna, Austria.
Rocha, F.C., Casatti, L., Carvalho, F.R. & Silva, A.M.S. 2009. Fish assemblages in
stream stretches occupied by cattail (Typhaceae, Angiospermae) stands in
Southeast Brazil. Neotropical Ichthyology, 7: 241–250.
Tilman, D., Lehman, C.L. & Thomson, K.T. 1997. Plant diversity and ecosystem
productivity: theoretical considerations. Proceedings of Natural Academy of
Sciences of the United States of America, 94: 1857–1861.
Villéger, S., Miranda, J.R., Hernández, D.F. & Mouillot, D. 2010. Contrasting changes
in taxonomic vs. functional diversity of tropical fish communities after habitat
degradation. Ecological Applications, 20: 1512–1522.
Warren, M.S., Hill, J.K., Thomas, J.A., Asher, J., Fox, R., Huntley,B., Roy, D.B.,
Telfer, M. G., Jeffcoate, S., Harding, P., Jeffcoate, G., Willis, S.G., Greatorex-
Davies, J.N., Moss D. & Thomas, C.D. 2001 Rapid responses of British butterflies
to opposing forces of climate and habitat change. Nature, 414: 65–69.
90
APÊNDICE
Tabela A1. Características funcionais das espécies que ocorreram em mais do que 3% dos segmentos. Dieta: Pei = peixes, Alg = algas, Per =
perifíton, Pro = Protista, Det = detrito, InvAu = invertebrados autóctones; InvAl = invertebrados alóctones; Preferência por fluxo: F1 = 0,0 a 0,25
cm/s; F2 = 0,26 a 0,45 cm/s; F3 = 0,46 a 0,65 cm/s; F4 > 0,66 cm/s; Uso de habitat: Ben = bentônico, Nec = nectônico, Necb = nectobentônico,
Mar = margem, Sup = superfície; Tam = tamanho; Tolerância à hipóxia*: tol = tolerantes, int = intolerantes.
Espécies Pei Alg Per Pro Veg Det InvAu InvAl F1 F2 F3 F4 Ben Nec Necb Mar Sup Tam Tol
Aspfus 0.00 0.00 0.00 0.01 0.00 0.02 0.97 0.00 0.19 0.47 0.39 1.00 2 0 0 0 0 37.45 tol
Astalt 0.00 0.07 0.00 0.00 0.48 0.00 0.02 0.42 0.82 0.76 0.24 0.04 0 2 1 1 0 84.80 int
Astfas 0.01 0.00 0.00 0.00 0.63 0.00 0.08 0.28 1.00 0.28 0.14 0.08 0 2 1 1 0 81.30 int
Brystr 0.00 0.00 0.00 0.00 0.00 0.00 0.54 0.46 1.00 0.25 0.12 0.05 0 2 0 0 1 48.10 int
Calcal 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00 0.50 1.00 0.00 1.00 2 0 0 1 0 66.60 tol
Cetihe 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.11 1.00 2 0 0 0 0 61.60 int
Chagom 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.07 0.35 1.00 2 0 0 0 0 55.40 int
Chazeb 0.00 0.00 0.00 0.00 0.01 0.00 0.99 0.00 0.39 0.37 0.44 1.00 2 0 0 0 0 62.45 int
Coraen 0.00 0.00 0.00 0.33 0.00 0.44 0.20 0.00 0.96 1.00 0.63 0.13 2 0 0 1 0 47.00 tol
Crebri 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 1.00 0.54 0.80 0.67 0 0 1 2 0 110.70 int
Eigtri 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.60 1.00 0.00 0 0 0 2 0 123.00 int
Eigvir 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.86 1.00 0.36 0 0 0 2 0 158.50 int
Geobra 0.09 0.00 0.00 0.36 0.22 0.00 0.29 0.04 1.00 0.63 0.09 0.05 0 0 2 1 0 180.00 tol
Gymcar 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.63 0.78 1.00 0.23 0 0 0 2 0 258.00 tol
91
Hisfra 0.00 0.00 0.82 0.00 0.00 0.18 0.00 0.00 0.86 1.00 1.00 0.72 0 0 0 2 0 39.40 int
Hopmal 0.95 0.00 0.00 0.00 0.00 0.00 0.05 0.00 0.93 1.00 0.10 0.07 0 0 1 2 0 198.00 tol
Hypanc 0.00 0.00 0.77 0.00 0.15 0.08 0.00 0.00 0.95 0.71 0.65 0.69 2 0 0 1 0 91.10 tol
Hypnig 0.00 0.00 0.98 0.00 0.02 0.00 0.00 0.00 0.30 0.14 0.43 1.00 2 0 0 0 0 81.80 int
Impmir 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.13 0.27 0.33 1.00 2 0 0 1 0 83.90 int
Impsch 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.03 0.38 0.27 1.00 2 0 0 1 0 54.25 int
Knomoe 0.00 0.00 0.00 0.00 0.00 0.00 0.98 0.02 1.00 0.79 0.33 0.36 0 2 1 0 1 35.80 int
Lepfri 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 1.00 0.14 0.33 0.18 0 1 2 0 0 62.50 int
Olipin 0.00 0.00 0.00 0.00 0.00 0.00 0.97 0.00 1.00 0.49 0.30 0.00 0 1 0 2 0 52.00 int
Parnas 0.00 0.03 0.10 0.00 0.07 0.50 0.30 0.00 0.00 0.00 1.00 0.50 2 0 0 0 0 90.35 int
Phahar 0.00 0.00 0.00 0.00 0.00 0.56 0.44 0.00 0.22 1.00 0.80 0.36 0 0 0 1 2 32.10 tol
Piaarg 0.00 0.00 0.00 0.00 0.04 0.00 0.92 0.04 1.00 0.78 0.79 0.32 0 1 2 0 0 60.65 int
Pimava 0.02 0.00 0.00 0.00 0.00 0.00 0.91 0.00 0.90 0.31 0.74 0.15 2 0 0 0 0 90.400 int
Poeret 0.00 0.00 0.00 0.00 0.00 0.82 0.18 0.00 0.27 0.50 0.20 0.52 0 0 0 1 2 30.000 tol
Psepul 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.03 0.01 0.41 1.00 2 0 0 0 0 66.100 int
Rhaque 0.00 0.00 0.00 0.00 0.02 0.00 0.98 0.00 0.71 0.91 0.46 0.64 2 0 0 1 0 139.600 tol
Sernot 0.00 0.15 0.00 0.00 0.15 0.61 0.10 0.00 1.00 0.10 0.02 0.00 0 2 0 0 1 32.300 int
Steins 0.00 0.07 0.00 0.07 0.07 0.74 0.05 0.00 1.00 0.00 0.00 0.00 0 0 2 0 0 93.750 int
Stemac 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 1.00 0.00 0.00 0 0 0 2 0 202.000 int
Synmar 1.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.87 0.65 0 0 0 2 0 260.000 tol
Tatnei 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 1.00 0.75 0.33 0.33 0 0 0 1 2 57.600 int
92
* Segundo Kramer & Mehegan (1981), Araujo & Garutti (2003), Bozzetti & Schulz (2004), Casatti et al. (2006b), Rocha et al. (2009) e Casatti et
al. (2012).
93
DISCUSSÃO GERAL .
Existe certo consenso de que o desmatamento das zonas ripárias é uma das
principais ameaças à conservação da ictiofauna (NAIMAN; DÉCAMPS, 1997;
BARLETTA et al., 2010). Isso foi reforçado nos estudos apresentados aqui, onde foi
possível constatar que esse tipo de alteração afeta, além da diversidade taxonômica
(Capítulo II), outros componentes da biodiversidade, como a composição e diversidade
funcional (Capítulos II e III). As modificações que o desmatamento da zona ripária
provoca nos riachos representam filtros ambientais que limitam a ocorrência de grupos
de espécies especializadas, em especial as reofílicas (Capítulo III), enquanto favorece
outras espécies que são capazes de explorar as novas oportunidades que se tornam
disponíveis (Capítulo II e III).
As métricas utilizadas nos capítulos II e III capturam aspectos diferentes da
diversidade funcional. A dispersão funcional (FDis) é uma medida de diversidade
funcional que é ponderada pela abundância das espécies, sendo portanto, menos
sensível a presença ou ausência de características diferenciadas exibidas por espécies
raras do que outros índices, como por exemplo a riqueza funcional (FRic). FDis é um
índice que apresenta alta correlação com o índice de entropia quadrática de Rao
(LALIBERTÉ; LEGENDRE, 2010), o qual foi utilizado no capítulo III e, embora seja
calculado de forma distinta, também considera a abundância das espécies. Quando os
riachos florestados e desmatados foram comparados na escala de meso-hábitats
constatou-se que, apesar dos corredores e poços desmatados abrigarem maior número
de espécies funcionalmente distintas (maior FRic), a estrutura quantitativa (medida pelo
FDis) das suas comunidades foi similar à dos florestados. Isso, provavelmente, deve-se
a dominância de espécies nectônicas, limnofílicas e de dieta insetívora/generalista
(grupo 6 do capítulo III) nesses meso-hábitats, cuja ocorrência é pouco influenciada
pela cobertura florestal (Capítulo III).
Somente as corredeiras florestadas apresentaram estrutura quantitativa
diferenciada, com menores valores de FDis, indicando maior redundância funcional, o
que foi atribuído a maior restrição ambiental presente nesse meso-hábitat (Capítulo II),
representada pelo alto fluxo de água e a baixa variabilidade na disponibilidade de
alimento (ANGERMEIER; KARR, 1983). As espécies que ocupam as corredeiras
(grupo 3 da partição #30 do Capítulo III) apresentam especializações para lidar com
essas restrições (CASATTI; CASTRO, 1998), sendo altamente sensíveis à
94
descaracterização desses meso-hábitats. Embora as corredeiras estejam presentes nos
riachos desmatados (Capítulo II) sua distribuição é limitada, além de apresentarem
estrutura física com maior grau de alteração (e.g., sedimentação, atenuação do fluxo
pelas macrófitas marginais). Isso provavelmente explica o declínio desse grupo de
espécies quando os riachos desmatados foram considerados como um todo (Capítulo
III). Esses resultados sugerem que o comprometimento da estrutura física do hábitat
causado pelo desmatamento é uma importante fonte de alteração da estrutura funcional
das comunidades de riachos.
Os resultados obtidos no capítulo II, em especial o aumento da riqueza de
espécies e de FRic nos riachos desmatados reforçam outros estudos com um enfoque
taxonômico (LORION; KENNEDY, 2009; TERESA; CASATTI, 2010) e proporcionam
suporte à idéia de que maiores valores de diversidade não são necessariamente um
indicativo de integridade biológica. Uma explicação que é usualmente evocada quando
maiores valores de diversidade são encontrados em ambientes submetidos às
interferências antrópicas é a Hipótese do Distúrbio Intermediário (HDI) (CONNELL,
1978). Essa hipótese baseia-se na idéia de que comunidades submetidas a distúrbios
muito freqüentes são representadas por espécies tolerantes e oportunistas, enquanto que
em níveis baixos de distúrbio as espécies mais especializadas, com maior habilidade
competitiva prosperam sobre as primeiras. Em níveis intermediários de distúrbios as
comunidades são compostas por uma combinação de espécies de ambos os extremos e,
consequentemente, são mais diversas (CONNELL, 1978). Portanto, para a HDI ser uma
explicação para o padrão observado no capítulo II seria necessário assumir que os
riachos desmatados estudados representam um estágio intermediário de degradação e
que os maiores valores de diversidade devem-se à co-ocorrência de espécies de ambos
os extremos de degradação. De fato, apesar dos riachos desmatados apresentarem
elementos indicadores de degradação como menor proporção de substrato estável e
dominância de macrófitas no hábitat interno (CASATTI et al., 2009b), eles ainda
possuem variabilidade de meso-hábitats (poços, corredores e corredeiras) e substrato
estável em proporção significativa, que são elementos associados a ambientes em
melhor estado de conservação (CASATTI, 2010). Entretanto, para avaliar se a HDI
pode explicar os padrões observados seria necessário contemplar um espectro maior de
condições de desmatamento, incluindo também riachos em estado de degradação
extrema (sensu CASATTI, 2010). Dessa forma, seria possível avaliar se o padrão de
diversidade obtido no capítulo II representa um pico seguido por um declínio em
95
riachos mais degradados (compatível com a HDI), ou se representa uma tendência
linear de aumento da diversidade.
Os padrões de diversidade observados também podem ser interpretados dentro
de uma perspectiva de funcionamento dos ecossistemas. A retirada da vegetação ripária
e aumento da radiação solar, combinado à menor retenção de nutrientes provenientes do
ambiente terrestre (Brodie & Mitchell, 2005) podem criar condições propícias para o
aumento da produtividade (Bunn et al., 1999). Considerando que a diversidade
funcional correlaciona-se bem com o funcionamento dos ecossistemas (Tilman et al.,
1997), esse aumento no metabolismo geral dos riachos desmatados pode estar associado
ao aumento da diversidade observado nestes ambientes (e.g., FRic). Além da
produtividade, outros processos ecossistêmicos podem estar associados com as
mudanças na composição funcional provocadas pelo desmatamento. Por exemplo, a
retirada da vegetação ripária aumenta a sedimentação, modificando a composição do
substrato e reduzindo o fluxo (Rabeni & Smale, 1995; Casatti et al., 2006), o que por
sua vez pode resultar no favorecimento das espécies bentônicas menos especializadas
em detrimento de espécies reofílicas de hábitos mais especializados (Capítulo III).
A problemática envolvendo o desmatamento é de extrema relevância,
principalmente ao considerar que as eminentes mudanças na legislação ambiental
brasileira (revisão do Código Florestal: Projeto de Lei nº 1.876/1999) prevêem redução
da área florestada ao redor dos rios e riachos. Diversos estudos têm alertado para os
impactos desse tipo de alteração sobre diversos aspectos (Toledo et al., 2010; Tundisi &
Matsumura-Tundisi, 2010), inclusive sobre a conservação da ictiofauna (Casatti, 2010;
Magalhães et al., 2011). Os resultados obtidos nos estudos apresentados aqui permitem
antecipar que a intensificação do desmatamento, que em algumas regiões já atinge a
maioria dos cursos d’água (Casatti et al., 2006; Silva et al., 2007), pode resultar em
mudanças no pool funcional dentro de um contexto regional. Isso porque o
desmatamento provoca homogeneização funcional da ictiofauna dentro dos riachos
(Capítulo II) e supostamente entre riachos (Casatti et al., 2009b), por meio da limitação
da ocorrência de espécies funcionalmente distintas daquelas que são favorecidos neste
contexto de degradação (Capítulo III).
96
CONCLUSÕES .
Os resultados obtidos nos capítulos I, II e III complementam o conhecimento
acumulado a partir da utilização da abordagem taxonômica sobre os impactos do
desmatamento em comunidades de peixes tropicais e fornecem subsídios para uma
avaliação mais completa dos efeitos desse tipo de alteração sobre a biota aquática. Os
padrões de diversidade e composição funcional diferiram em função das alterações
provocadas nos riachos pelo desmatamento, corroborando a hipótese geral. As
características funcionais associadas com o uso de hábitat, ecologia trófica, tolerância à
hipóxia e tamanho influenciaram a distribuição das espécies dentro (entre meso-
hábitats) e entre riachos em contextos de degradação contrastantes. Nesse sentido, as
modificações na estrutura funcional das comunidades causadas pelo desmatamento
parecem ser mediadas pelos efeitos desse tipo de alteração sobre a disponibilidade de
hábitat e recursos alimentares, assim como nas condições físico-químicas (e.g.,
hipóxia), que por sua vez representam filtros ambientais importantes na organização das
comunidades.
97
REFERÊNCIAS GERAIS .
ALLAN, J. D.; FLECKER, A. S. Biodiversity conservation in running waters.
BioScience, v. 43, n. 1, p. 32–43, 1993.
ANGERMEIER, P. L.; KARR, J. R. Fish communities along environmental gradients
in a system of tropical streams. Environmental Biology of Fishes, v. 9, p. 117–135,
1983.
BARLETTA, M. et al. Fish and aquatic habitat conservation in South America: a
continental overview with emphasis on Neotropical systems. Journal of Fish Biology,
v. 76, p. 2118–2176, 2010.
Barragán, F. et al. Negative impacts of human land use on dung beetle functional
diversity. PLoS One, v. 6, p. e17976, 2011.
BOJSEN, B. H.; BARRIGA, R. Effects of deforestation on fish community structure in
Ecuadorian Amazon streams. Freshwater Biology, v. 47, p. 2246–2260, 2002.
BRODIE, J. E;& MITCHELL, A. W. Nutrients in Australian tropical rivers: changes
with agricultural development and implications for receiving environments. Marine
and Freshwater Research, v. 56, n. 3, p. 279–302, 2005.
BUNN, S. E.; DAVIES, P. M; MOSISCH, T. D. Ecosystem measures of river health
and their response to riparian and catchment degradation. Freshwater Biology, v. 41, p.
333–345, 1999.
CASATTI, L. Alterações no Código Florestal Brasileiro: impactos potenciais sobre a
ictiofauna. Biota Neotropica, v. 10, n. 4, p. 31-34, 2010.
CASATTI, L.; CASTRO, R. M. C. A fish community of the São Francisco river
headwater riffles, southeastern Brazil. Ichthyological Exploration of Freshwaters, v.
9, p. 229–242, 1998.
CASATTI, L. et al. Stream fish, water and habitat quality in a pasture dominated basin,
southeastern Brazil. Brazilian Journal of Biology, v. 66, p. 681–696, 2006.
CASATTI, L.; FERREIRA, C. P.; LANGEANI, F. A fish-based biotic integrity index
for assessment of lowland streams in southeastern Brazil. Hydrobiologia, v. 623, n. 1,
p. 173–189, 2009a.
CASATTI, L.; FERREIRA, C. P.; CARVALHO, F. R. Grass-dominated stream sites
exhibit low fish species diversity and dominance by guppies: an assessment of two
tropical pasture river basins. Hydrobiologia, v. 632, p. 273–283, 2009b.
98
DIAZ, S.; CABIDO, M. Vive la différence: plant functional diversity matters to
ecosystem processes. Trends in Ecology and Evolution, v. 16, n. 8, p. 646-655, 2001.
CONNELL, J. H. Diversity in tropical rain forest and coral reefs. Science, v. 199, p.
1302-1310, 1978.
ERNST, R.; LINSENMAIR, K. E.; RÖDEL, M. O. Diversity erosion beyond the
species level: dramatic loss of functional diversity after selective logging in two tropical
amphibian communities. Biological Conservation, v. 133, p. 143–155, 2006.
EROS, T. et al. Characterizing functional trait diversity and trait–environment
relationships in fish assemblages of boreal lakes. Freshwater Biology, v. 54, p. 1788–
1803, 2009.
FAUSCH, K. D.; NORTHCOTE, T. G. Large woody debris and salmonid habitat in a
small coastal British Columbia stream. Canadian Journal of Aquatic Sciences, v. 49,
n. 4, p. 682-693, 1992.
FLYNN, D. F. B. et al. Loss of functional diversity under land use intensification across
multiple taxa. Ecology Letters, v. 12, p. 22–33, 2009.
GREGORY, S. V. An ecosystem perspective of riparian zones. BioScience, v. 41, p.
540–551, 1991.
GROWNS, I. et al. A comparison of fish assemblages associated with different riparian
vegetation types in the Hawkesbury–Nepean river system. Fisheries Management &
Ecology, v. 10, p. 209–220, 2003.
LALIBERTÉ, E.; LEGENDRE, P. A distance-based framework for measuring
functional diversity from multiple traits. Ecology, v. 91, p. 299–305, 2010.
LORION, C. M. & KENNEDY, B. P. Riparian forest buffers mitigate the effects of
deforestation on fish assemblages in tropical headwater streams. Ecological
Applications, v. 19, p. 468–479, 2009.
LOWE-MCCONNELL, R. H. Estudos ecológicos de comunidades de peixes
tropicais. São Paulo: Edusp, 1999.
MAGALHÃES, A. L. B.; CASATTI, L.; VITULE, J. R. S. Alterações no código
florestal brasileiro favorecerão espécies não-nativas de peixes de água doce. Natureza
& Conservação, v. 9, n. 1, p. 121-124, 2011.
MOUILLOT, D. et al. Functional structure of biological communities predicts
ecosystem multifunctionality. PLoS One, v. 6, p. e17476, 2011.
99
MURPHY, M. L.; HALL, J. D. Varied effects of clear-cut logging on predators andtheir
habitat in small streams of the Cascade Mountains, Oregon. Canadian Journal of
Fisheries and Aquatic Sciences, v. 38, p. 137-145, 1981.
NAIMAN, R. J.; DÉCAMPS, H. The ecology of interfaces: riparian zones. Annual
Review of Ecology and Systematics, v. 28, p. 621–658, 1997.
PETCHEY, O. L.; GASTON, K. J. Functional diversity: back to basics and looking
forward. Ecology Letters, v. 9, p. 741–758, 2006.
PUSEY, B. J.; ARTHINGTON, A. H. Importance of the riparian zone to the
conservation and management of freshwater fish: a review. Marine and Freshwater
Research, v. 54, p 1–16, 2003.
RABENI, C. F.; SMALE, M. A. Effects on siltation on stream fishes and the potential
mitigating role of the buffering riparian zone. Hydrobiologia, v. 303, p. 211–219, 1995.
ROCHA, F. C. et al. Fish assemblages in stream stretches occupied by cattail
(Typhaceae, Angiospermae) stands in Southeast Brazil. Neotropical Ichthyology, v. 7
p. 241–250, 2009.
SCHNEIDER, K. N.; WINEMILLER, K. O. Structural complexity of woody debris
patches influences fish and macroinvertebrate species richness in a temperate
floodplain-river system. Hydrobiologia, v. 610, p. 235–244, 2008.
SILVA, A. M. Historical land-cover/use in different slope and riparian buffer zones in
watershed of the State of São Paulo, Brazil. Scientia Agricola, v. 64, p. 325–335, 2007.
SMA/IF (Secretaria do Meio Ambiente / Instituto Florestal). Inventário florestal da
vegetação natural do Estado de São Paulo. São Paulo: Imprensa Oficial do Estado de
São Paulo, 2005.
SWEENEY, B. W. et al. Riparian deforestation, stream narrowing, and loss of stream
ecosystem services. Proceedings of the National Academy of Sciences of the United
States of America, v. 101, p. 14132–14137, 2004.
TERESA, F. B.; CASATTI, L. Importância da vegetação ripária em região
intensamente desmatada no sudeste do Brasil: um estudo com peixes de riacho. Pan-
American Journal of Aquatic Sciences, v. 5, p. 444–453, 2010.
TILMAN, D.; LEHMAN, C. L.; THOMSON, K. T. Plant diversity and ecosystem
productivity: theoretical considerations. Proceedings of Natural Academy of Sciences
of the United States of America, v. 94, p. 1857–1861, 1997.
TOLEDO, L. F. et al. A revisão do Código Florestal Brasileiro: impactos negativos para
a conservação dos anfíbios. Biota Neotropica, v. 10, n. 4, p. 35–38, 2010.
100
TUNDISI, J. G.; MATSUMURA-TUNDISI, T.M. Impactos potenciais das alterações
do Código Florestal nos recursos hídricos. Biota Neotropica, v. 10, n. 4, p. 67–76,
2010.
VILLÉGER, S. et al. Contrasting changes in taxonomic vs. functional diversity of
tropical fish communities after habitat degradation. Ecological Applications, v. 20, p.
1512–1522, 2010.
Recommended