Prévia do material em texto
UNIVERSIDADE FEDERAL DE GOIÁS (UFG) INSTITUTO DE CIÊNCIAS BIOLÓGICAS (ICB) PROGRAMA DE PÓS-GRADUAÇÃO EM ECOLOGIA E EVOLUÇÃO MAISA CARVALHO VIEIRA Estruturação da comunidade zooplanctônica em reservatórios tropicais GOIÂNIA 2022 MAISA CARVALHO VIEIRA Estruturação da comunidade zooplanctônica em reservatórios tropicais Tese apresentada ao Programa de Pós- Graduação em Ecologia e Evolução, da Universidade Federal de Goiás (UFG), como requisito para obtenção do título de Doutora em Ecologia e Evolução. Área de concentração: Ecologia e Evolução Orientador: Prof. Dr. Luis Mauricio Bini Co-orientador: Prof. Dr. Ludgero Cardoso Galli Vieira Goiânia – GO 2022 Por acaso você chegou, fez morada em meu coração e na minha alma. O maior motivo da minha existência. Dedico à Isis. 9 AGRADECIMENTOS Fazer um doutorado não é fácil, fazer um doutorado com uma gravidez, puerpério e toda mudança que isso gera na vida de uma mulher é ainda mais difícil e, junto a isso tudo uma pandemia. Digamos que meu doutorado foi nível ‘hard” e se hoje estou aqui concluindo é graças a muitas pessoas que estiveram nos bastidores. Primeiramente, eu agradeço a Deus, por ter sido meu refúgio nos momentos de angústia e dor, por ter sido minha fortaleza nos momentos que me senti perdida e por ter me guiado à felicidade e manutenção da saúde mental. Agradeço minha mãe por ter sido minha principal rede de apoio fora de casa. Agradeço ao meu companheiro, Hugo, que além de ser meu amigo, namorado, marido, pai da minha filha, ainda é revisor dos meus textos, conselheiro, “campoman”, entre várias outras qualidades que poderia citar aqui em várias e várias linhas. Sem vocês, com certeza, eu não teria conseguido sozinha. Agradeço também ao meu pai que acreditou em mim, me ensinou enfrentar os desafios e aprender com os erros. Agradeço a Carla, que se prontificou em identificar minhas amostras quando eu descobri que não poderia por conta da gravidez. Eu nem sei como agradecer você Carlinha, você é uma amiga incrível! Agradeço também a minha amiga Rafa G., que sempre que precisei estendeu sua mão para me ajudar, sempre me ouviu, me aconselhou e me acolheu. Foi um pouco de psicóloga, “campogirl”/fotógrafa, várias risadas e fofocas. Também agradeço ao Jean, que sempre tirou minhas dúvidas em relação a estatística, você é muito bom cara! Agradeço também ao restante dos meus amigos do laboratório Lea“t”, Matheus, Fagner, Letícia, Sara, Aline, Flávio, Gisele, pelas boas conversas, risadas e cafés. Agradeço também aos amigos do Nepal que me ajudaram na amostragem de campo, Leo B., Hasley, Gustavo G., Leo G., Carol, João Paulo, Thalia, Gleicon e Pedro (motorista). Não é fácil subir e descer barco em 40 represas, ficar 10 remando o tempo todo porque não tínhamos motor, empurrar a van quando atolou (por diversas vezes) e ainda espirrar lama na cara, levar mordida de cachorro em fazenda, correr de bode, colocar a mão em água muito duvidosa, entre outras coisas que aconteceram no campo (as fotos embaixo resumem bem meu agradecimento a vocês rsrs). Fonte: Acervo pessoal, PAD-DF (2018). Agradeço ao professor Júnior, por ter cedido o laboratório AquaRiparia (UnB) para algumas análises físico-químicas da água, e a Daphne na Embrapa Cerrados, por ter me auxiliado nas análises de nitrogênio e fósforo total. Agradeço ao professor Fabrício, que foi minha inspiração para entrar na academia, um excelente professor, e se hoje estou terminando um doutorado é porque lá na minha graduação me incentivou a tal. Agradeço ao prof. Luiz Felipe Velho, que me forneceu os dados do terceiro capítulo da minha tese. Também agradeço aos professores da UnB que tiveram papeis fundamentais na minha formação, entre eles o prof. Antônio Felipe (in memoriam), com seu empolgamento pelo conhecimento científico e que faz muita falta entre nós, e o prof. Ludgero, que é uma das grandes referências para mim hoje, excelente profissional e acessível. Também agradeço aos professores da UFG pelos ensinamentos e exemplos, em especial as prof. Priscila e Jascieli, duas mulheres fortes nas quais me espelho. 11 Agradeço ao meu orientador Luis Mauricio Bini, por ter me dado a oportunidade de aprender com ele. Com certeza, foi uma das pessoas que mais me surpreendeu no doutorado, e com muita sabedoria soube conduzir diversas situações. Você com certeza faz toda diferença com seus orientandos, é compreensível, é um excelente professor, e tudo que se propõe a fazer, faz com eficiência e maestria. E agora, eu entendo a admiração do Ludgero por você. Agradeço aos apoios financeiros, à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) pela bolsa de estudo, à Fundação de Amparo à Pesquisa do Distrito Federal (FAPDF) pelo financiamento do projeto que resultou em dois capítulos dessa tese. Agradeço ao professor Adriano Melo e Fernando Lansac-Tôha por avaliarem meus relatórios de acompanhamento durante o período do doutorado. Também agradeço aos membros da banca por aceitarem o convite e contribuírem para a melhoria desse trabalho. Agradeço aos meus amigos de fora da academia, principalmente minha amiga Juliana (Xuliana), que sempre ouviu meus desabafos do dia a dia. Sua amizade é muito importante para mim (mesmo que você nem vá ver isso daqui). E não poderia me esquecer também das pessoas que me ajudaram com a Isis (sem ser minha mãe e o Hugo). À Rosana, que foi a primeira pessoa que confiei em cuidar da Isis. Foi em um momento muito delicado da pandemia, onde eu precisava me esconder no escritório para conseguir estudar, enquanto ela ficava com a Isis para mim, fingindo que eu não estava em casa (para que as coisas fluíssem melhor). À Maria Clara, que ama a Isis mais do que tudo, e sempre que precisei deixar a Isis aos seus cuidados nunca me negou ajuda, mesmo não recebendo nada por isso (você tem minha eterna gratidão). 12 À Vanessa e ao Gustavo, meus vizinhos/primos, que sempre me socorreram em reuniões ou em imprevistos de última hora, e se prontificaram em me ajudar com ela. Meu muito obrigada a todos vocês! 13 SUMÁRIO Resumo .......................................................................................................................... 14 Abstract ......................................................................................................................... 15 Introdução Geral .......................................................................................................... 16 Capítulo 1: Correlates of zooplankton metacommunity structure in small reservoirs ....................................................................................................................... 26 Abstract ........................................................................................................................... 27 Introduction .................................................................................................................... 28 Methods .......................................................................................................................... 32 Results ............................................................................................................................ 36 Discussion ....................................................................................................................... 39 References ...................................................................................................................... 45 Supplementary Material .................................................................................................56 Capítulo 2: Diversidade beta temporal de comunidades locais zooplanctônicas em pequenos reservatórios ................................................................................................. 71 Resumo ........................................................................................................................... 72 Abstract ........................................................................................................................... 73 Introdução ....................................................................................................................... 74 Metodologia .................................................................................................................... 77 Resultados ....................................................................................................................... 82 Discussão ........................................................................................................................ 94 Referências ..................................................................................................................... 98 Material Suplementar ................................................................................................... 110 Capítulo 3: Evidence that dams promote biotic differentiation of zooplankton communities in two Brazilian reservoirs .................................................................. 128 Abstract ......................................................................................................................... 130 Introduction .................................................................................................................. 131 Material and methods ................................................................................................... 134 Results .......................................................................................................................... 139 Discussion ..................................................................................................................... 146 Conclusions and caveats ............................................................................................... 150 References .................................................................................................................... 151 Supplementary Material ............................................................................................... 159 Considerações Finais .................................................................................................. 174 14 Resumo Nesta tese, investigamos os padrões e processos das metacomunidades zooplanctônicas e suas variações espaciais e temporais utilizando dados de 39 pequenos reservatórios (Capítulo 1 e 2) e dois reservatórios hidrelétricos (Capítulo 3). No primeiro capítulo, quantificamos a importância de diferentes variáveis preditoras (variáveis ambientais locais, de uso e ocupação do solo, espaciais e comunidades biológicas passadas) na estruturação de metacomunidades zooplanctônicas. Encontramos que a comunidade passada foi importante preditor e que a importância dos mecanismos variaram ao longo do tempo, sugerindo necessidade de incluir dados temporais nas análises de metacomunidades. No segundo capítulo, avaliamos a diversidade beta temporal e sua relação com a variabilidade ambiental, a área do reservatório, a vegetação remanescente, a densidade média do fitoplâncton e a quantidade de pequenos reservatórios próximos ao estudado. Nós não encontramos relações significativas com os preditores, mas encontramos altos valores de diversidade beta temporal, e de perdas de abundância de táxons entre períodos de seca e chuva e ganhos de abundância de táxons entre períodos de chuva e seca. Esses resultados sugerem a importância da variação sazonal e hidrológica sobre a diversidade beta temporal, e indicam que nosso conhecimento sobre os determinantes da dinâmica temporal dessas comunidades ainda é limitado. No terceiro capítulo avaliamos se a diversidade beta da comunidade zooplanctônica aumentaria após o barramento de dois reservatórios hidrelétricos. Nossos resultados foram consistentes com a hipótese de diferenciação biótica e indicaram a importância de manter amostragens temporais. Palavras-Chave: Metacomunidade; Diversidade beta; Composição de espécies; Pequenos reservatórios; Reservatório hidrelétrico; Zooplâncton; Temporal 15 Abstract In this thesis, we investigated the patterns and processes of zooplanktonic metacommunities and their spatial and temporal variations using data from 39 small reservoirs (Chapters 1 and 2) and two hydroelectric reservoirs (Chapter 3). In the first chapter, we quantify the importance of different predictor variables (local environmental variables, land use and occupation, spatial variables, and past biological communities) in the zooplanktonic metacommunities structuring. We found that past community was an important predictor and that the importance of mechanisms varied over time, suggesting the need to include temporal data in metacommunity analyses. In the second chapter, we evaluated the temporal beta diversity and its relationship with environmental variability, the reservoir area, the remaining vegetation, the average phytoplankton density, and the small reservoirs number close to the one studied. We did not find significant relationships with the predictors, but we did find high values of temporal beta diversity, and taxon abundance losses between rain and dry periods and taxon abundance gains between dry and rain periods. These results suggest the importance of seasonal and hydrological variation on temporal beta diversity and indicate that our knowledge about the temporal dynamics determinants of these communities is still limited. In the third chapter, we evaluated whether the beta diversity of the zooplankton community would increase after the damming of two hydroelectric reservoirs. Our results were consistent with the biotic differentiation hypothesis and indicated the importance of maintaining temporal sampling. Keywords: Metacommunity; Beta diversity; Species composition; Small reservoirs; Hydroelectric reservoir; Zooplankton; Temporal 16 Introdução Geral Um dos principais objetivos da ecologia é compreender a estruturação das comunidades biológicas (ou seja, riqueza, composição e distribuição da abundância entre espécies). Para isso, a teoria de metacomunidades prevê que a estruturação das comunidades locais é determinada por interações bióticas, filtragem ambiental e dispersão (Leibold et al., 2004, 2010; Vellend, 2010; Leibold & Chase, 2017). As variações na abundância e na composição de espécies entre os locais (diversidade beta espacial) e períodos de tempos (diversidade beta temporal) de uma mesma região (Anderson et al., 2006, 2011; Legendre & Gauthier, 2014), também são fundamentais para o arcabouço da teoria de metacomunidades (Tuomisto & Ruokolainen, 2006; Heino et al., 2015). Diferentes mecanismos podem atuar para determinar a estrutura das comunidades, entre eles, os mais estudados são baseados em processos ambientais e espaciais (e.g., Soininen et al., 2011; Zhao et al., 2017; Rocha et al., 2020; Sinclair et al., 2021). Os processos ambientais estão relacionados ao nicho ecológico (Hutchinson, 1957) e as espécies são “filtradas” pelas condições locais segundo seus requerimentos ecológicos (Mittelbach & Schemske, 2015). Em estudos empíricos, as variáveis ambientais locais (e.g., limnológicas para ambientes aquáticos) são utilizadas frequentemente para representar o processo de filtragem ambiental, a heterogeneidade ambiental e a variabilidade ambiental (e.g., Declerck et al., 2011; Hatosy et al., 2013; Padial et al., 2014; Zorzal-Almeida et al., 2017). Frequentemente,os processos espaciais são relacionados a dispersão dos organismos (por exemplo, Cottenie, 2005). Entretanto, sua inferência é complexa, pois os processos espaciais também podem depender de fatores como a extensão do gradiente ambiental (Heino et al., 2015), 17 extensão espacial da região em estudo (Cottenie 2005), tamanho do propágulo e do organismo (De Bie et al. 2012; Padial et al. 2014; Martin et al. 2021), ou outras variáveis espacialmente autocorrelacionadas que não foram mensuradas no estudo (Peres-Neto e Legendre 2010; Diniz-Filho et al. . 2012). Contudo, considerando a complexidade biológica e dos ambientes, os mecanismos relacionados aos processos ambientais e espaciais supracitados na estruturação das comunidades podem ser particionados, entre eles, em variáveis ambientais locais (e.g. variáveis limnológicas), em variáveis na escala de paisagem (e.g. uso e ocupação do solo), tamanho da área, em relação a produtividade primária e a quantidade de locais “fontes” de populações próximos ao estudado. Recentemente, uma outra dimensão é relacionada à fatores históricos como possíveis mecanismos da estruturação das comunidades locais (Andersson et al., 2014; Castillo-Escrivà et al., 2017; Oliveira et al., 2020; Ortega et al., 2021), no qual a comunidade do presente pode ser fortemente influenciada pela comunidade do passado. A força dessa relação pode depender de algumas características, como o tempo de geração de espécies, o intervalo de tempo entre as amostragens e variabilidade ambiental temporal. Ambientes aquáticos de água doce são ótimos modelos para o estudo dos padrões e processos das comunidades biológicas, pois além da sua delimitação física mais facilmente determinada (do que ambientes terrestres, por exemplo), também há uma grande variabilidade ambiental e biológica ao longo do espaço e do tempo (Reynolds, 1998; Tundisi, 2003; Leibold et al., 2004; Fernandes et al., 2013). Além disso, os ecossistemas aquáticos de água doce estão entre os ambientes naturais mais ameaçados do mundo, por exemplo, por causa dos represamentos de rios (Dudgeon et al., 2006; Reid et al., 2019). O barramento promove diferentes alterações hidrológicas (Ward, 1989) e essas alterações podem levar a grandes mudanças na biodiversidade, 18 dependendo do tamanho, disposição da barragem e modo de operação (Pereira et al., 2020; Picapedra et al., 2020). Muitos estudos já relataram que a construção de barragens altera a biodiversidade de peixes, bentos, plâncton e microrganismos aquáticos (Agostinho et al., 2008; Pelicice et al., 2015; Simões et al., 2015; Rechisky et al., 2013; Braghin et al., 2018; Schmidt et al., 2020). Dessa forma, a diversidade beta é uma ferramenta importante para avaliar tendências negativas (homogeneização biótica) e positivas (diferenciação biótica) do impacto do represamento nas comunidades biológicas. Considerando comunidades zooplanctônicas, o efeito da barragem ocorre principalmente pelas mudanças na velocidade do fluxo, sedimentação, estratificação térmica, disponibilidade de alimentos, qualidade da água, entre outros parâmetros (ver mais em Braghin et al., 2018; Schmidt et al., 2020). Essas comunidades são muito diversas e possuem diferentes funções nos ecossistemas, por exemplo, participação no fluxo de energia nos ecossistemas aquáticos e ciclagem de nutrientes (Allan, 1976; Ferdous & Muktadir, 2009). Os principais organismos que fazem parte dessa comunidade são os rotíferos, microcrustáceos (cladóceros e copépodes) e amebas testáceas. Desse modo, o objetivo geral dessa tese foi investigar os padrões e os processos da estruturação das comunidades zooplanctônicas (e de seus respectivos grupos, microcrustáceos, rotíferos e amebas testáceas) em reservatórios tropicais. Essa tese está dividida em três capítulos. No primeiro capítulo, “Correlates of zooplankton metacommunity structure in small reservoirs”, quantificamos a importância relativa de diferentes variáveis preditoras (variáveis ambientais locais, de uso e ocupação do solo, espaciais e efeito histórico biológico) na estruturação de metacomunidades zooplanctônicas (e de seus grupos) de pequenos reservatórios no Distrito Federal e em Goiás. No segundo capítulo, “Diversidade beta temporal de comunidades locais 19 zooplanctônicas em pequenos reservatórios”, calculamos a diversidade beta temporal do zooplâncton (e de seus grupos) e suas partições de perdas e ganhos de abundâncias de espécies em pequenos reservatórios no Distrito Federal e em Goiás, e avaliamos a relação com a variabilidade ambiental, a área do reservatório, a vegetação remanescente, a densidade média do fitoplâncton e a quantidade de pequenos reservatórios próximos ao estudado. No terceiro capítulo, “Evidence that dams promote biotic differentiation of zooplankton communities in two Brazilian reservoirs”, testamos a hipótese de que a diversidade beta da comunidade zooplanctônica aumentaria (diferenciação biótica) após o represamento em dois reservatórios hidrelétricos (reservatório Santo Antônio do Jari, na divisa do Amapá com o Pará, e reservatório Serra do Facão, em Goiás). No primeiro capítulo, os resultados indicaram que o efeito histórico foi importante para explicar a metacomunidade zooplanctônica, e em menor grau, as variáveis ambientais locais e espaciais. No segundo capítulo, encontramos altos valores de diversidade beta temporal, com características marcantes de perdas de abundância de táxons entre períodos de seca e chuva, e ganhos de abundância de táxons entre períodos de chuva e seca. Nossas variáveis preditoras não explicaram a variação da comunidade ao longo do tempo. No terceiro capítulo, encontramos que houve uma tendência à diferenciação biótica da comunidade zooplanctônica após o represamento, sendo mais evidente no reservatório de armazenamento (Serra do Facão). Referências Agostinnho, A. A., F. M. Pelicice, L. C. Gomes, 2008. Dams and the fish fauna of the Neotropical region: impacts and management related to diversity and fisheries. Brazilian Journal of Biology 68: 1119-1131. 20 Allan, J. D., 1976. Life History Patterns in Zooplankton. The American Naturalist 110: 165–180. Anderson, M. J., T. O. Crist, J. M. Chase, M. Vellend, B. D. Inouye, A. L. Freestone, N. J. Sanders, H. V. Cornell, L. S. Comita, K. F. Davies, S. P. Harrison, N. J. B. Kraft, J. C. Stegen, & N. G. Swenson, 2011. Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecology Letters 14: 19–28. Anderson, M. J., K. E. Ellingsen, & B. H. McArdle, 2006. Multivariate dispersion as a measure of beta diversity. Ecology Letters 9: 683–693. Andersson, M. G. I., M. Berga, E. S. Lindstrom, & S. Langenheder, 2014. The spatial structure of bacterial communities is influenced by historical environmental conditions. Ecology 95: 1134–1140. Braghin, L. S. M., B. A., Almeida, D. C. Amaral, T. F. Canella, B. C. G. Gimenez, C. C. Bonecker, 2018. Effects of dams decrease zooplankton functional β-diversity in river-associated lakes. Freshwater Biology 63: 721-730. Castillo-Escrivà, A., L. Valls, C. Rochera, A. Camacho, & F. Mesquita-Joanes, 2017. Disentangling environmental, spatial, and historical effects on ostracod communities in shallow lakes. Hydrobiologia 787: 61–72. Cottenie, K., 2005. Integrating environmental and spatial processes in ecological community dynamics. Ecology Letters 8:1175–1182. De Bie T., L. De Meester, L. Brendonck, M. B. Goddeeris, D. Ercken, H. Hampel, L. Denys, L. Vanhecke, K. Van der Gucht, J. Van Wichelen, W. Vyverman, S. A. J. Declerck, 2012. Body size and dispersal mode as key traits determining metacommunity structure of aquatic organisms. Ecology Letters 15:740–747. Declerck, S. A. J., J. S. Coronel, P. Legendre, & L. Brendonck, 2011. Scale dependency of processes structuring metacommunitiesof cladocerans in temporary pools of High- 21 Andes wetlands. Ecography 34: 296–305. Diniz-Filho J. A. F., T. Siqueira, A. A. Padial, T. F. Rangel, V. L. Landeiro, L. M. Bini, 2012. Spatial autocorrelation analysis allows disentangling the balance between neutral and niche processes in metacommunities. Oikos 121:201–210. Dudgeon, D., A. H. Arthington, M. O. Gessner, Z.-I. Kawabata, D. J. Knowler, C. Lévêque, R. J. Naiman, A.-H. Prieur-Richard, D. Soto, M. L. J. Stiassny, & C. A. Sullivan, 2006. Freshwater biodiversity: importance, threats, status and conservation challenges. Biological Reviews of the Cambridge Philosophical Society 81: 163–182. Ferdous, Z., & A. K. M. Muktadir, 2009. A Review: Potentiality of Zooplankton as Bioindicator. American Journal of Applied Sciences Science Publications 6: 1815– 1819. Fernandes, I. M., R. Henriques-Silva, J. Penha, J. Zuanon, & P. R. Peres-Neto, 2013. Spatiotemporal dynamics in a seasonal metacommunity structure is predictable: the case of floodplain-fish communities. Ecography 37: 464-475. Hatosy, S. M., J. B. H. Martiny, R. Sachdeva, J. Steele, J. A. Fuhrman, & A. C. Martiny, 2013. Beta diversity of marine bacteria depends on temporal scale. Ecology 94: 1898–1904. Heino, J., A. S. Melo, L. M. Bini, F. Altermatt, S. A. Al-Shami, D. G. Angeler, N. Bonada, C. Brand, M. Callisto, K. Cottenie, O. Dangles, D. Dudgeon, A. Encalada, E. Göthe, M. Grönroos, N. Hamada, D. Jacobsen, V. L. Landeiro, R. Ligeiro, R. T. Martins, M. L. Miserendino, C. S. Md Rawi, M. E. Rodrigues, F. de O. Roque, L. Sandin, D. Schmera, L. F. Sgarbi, J. P. Simaika, T. Siqueira, R. M. Thompson, & C. R. Townsend, 2015. A comparative analysis reveals weak relationships between ecological factors and beta diversity of stream insect metacommunities at two spatial levels. Ecology and Evolution 5: 1235–1248. 22 Hutchinson, G. E., 1957. Concluding remarks. Quantitative Biology 22: 415–427. Legendre, P., & O. Gauthier, 2014. Statistical methods for temporal and space-time analysis of community composition data. Proceedings of the Royal Society B: Biological Sciences 281: 20132728. Leibold, M. A., & J. M. Chase, 2017. Metacommunity Ecology. Princeton University Press, Princeton, NJ. Leibold, M. A., E. P. Economo, & P. Peres-Neto, 2010. Metacommunity phylogenetics: Separating the roles of environmental filters and historical biogeography. Ecology Letters 13: 1290–1299. Leibold, M. A., M. Holyoak, N. Mouquet, P. Amarasekare, J. M. Chase, M. F. Hoopes, R. D. Holt, J. B. Shurin, R. Law, D. Tilman, M. Loreau, & A. Gonzalez, 2004. The metacommunity concept: a framework for multi-scale community ecology. Ecology Letters 7: 601–613. Martin G.K., B.E. Beisner, F. J. J. Chain, M. E. Cristescu, P. A. del Giorgio, A. M. Derry, 2021. Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation. Ecology 102: e03224. Mittelbach, G. G., & D. W. Schemske, 2015. Ecological and evolutionary perspectives on community assembly. Trends in Ecology & Evolution 30: 241–247. Oliveira, F. J. M., D. P. Lima-Junior, & L. M. Bini, 2020. Current environmental conditions are weak predictors of fish community structure compared to community structure of the previous year. Aquatic Ecology 54: 729–740. Ortega, J. C. G., J. Geijer, J. Bergsten, J. Heino, J. Herrmann, F. Johansson, & L. M. Bini, 2021. Spatio-temporal variation in water beetle assemblages across temperate freshwater ecosystems. Science of The Total Environment 792: 148071. Padial, A. A., F. Ceschin, S. A. J. Declerck, L. De Meester, C. C. Bonecker, F. A. 23 Lansac-Tôha, L. Rodrigues, L. C. Rodrigues, S. Train, L. F. M. Velho, & L. M. Bini, 2014. Dispersal ability determines the role of environmental, spatial and temporal drivers of metacommunity structure. PLoS ONE 9: 1–8. Pelicice, F. M., P. S. Pompeu, A. A. Agostinho, 2015. Large reservoirs as ecological barriers to downstream movements of Neotropical migratory fish. Fish and Fisheries 16: 697-715. Pereira, H. R., L. F. Gomes, H. O. Barbosa, F. M. Pelicice, J. C. Nabout, F. B. Teresa, & L. C. G. Vieira, 2020. Research on dams and fishes: determinants, directions, and gaps in the world scientific production. Hydrobiologia 847: 579–592. Peres-Neto, P.R. & P. Legendre, 2010. Estimating and controlling for spatial structure in the study of ecological communities. Glob Ecol Biogeogr 19:174–184. Picapedra, P. H. S., C. Fernandes, J. Taborda, G. Baumgartner, & P. V. Sanches, 2020. A long-term study on zooplankton in two contrasting cascade reservoirs (Iguaçu River, Brazil): Effects of inter-annual, seasonal, and environmental factors. PeerJ 8:e8979. Rechisky, E. L., D. W. Welch, A. D. Porter, M. C. Jacobs-Scott, & P. M. Winchell, 2013. Influence of multiple dam passage on survival of juvenile Chinook salmon in the Columbia River estuary and coastal ocean. Proceedings of the National Academy of Sciences of the United States of America 110: 6883–6888. Reid, A. J., A. K. Carlson, I. F. Creed, E. J. Eliason, P. A. Gell, P. T. J. Johnson, K. A. Kidd, T. J. Maccormack, J. D. Olden, S. J. Ormerod, J. P. Smol, W. W. Taylor, K. Tockner, J. C. Vermaire, D. Dudgeon, & S. J. Cooke, 2019. Emerging threats and persistent conservation challenges for freshwater biodiversity. Biological Reviews 94: 849–873. Reynolds, C. S., 1998. The state of freshwater ecology. Freshwater Biology 39: 741– 753. 24 Rocha, B. da S., C. A. de Souza, K. B. Machado, L. C. G. Vieira, & J. C. Nabout, 2020. The relative influence of the environment, land use, and space on the functional and taxonomic structures of phytoplankton and zooplankton metacommunities in tropical reservoirs. Freshwater Science 39: 321–333. Schmidt, J., P. D. B. Andrade, A. A. Padial, 2020. Zooplankton trajectory before, during and after a hydropower dam construction. Acta Limnologica Brasiliensia 32: e18. Simões, N.R., A.H. Nunes, J.D. Dias, F.A. Lânsac-Tôha, L.F.M. Velho & C.C. Bonecker, 2015. Impact of reservoirs on zooplankton diversity and implications for the conservation of natural aquatic environments. Hydrobiologia 758: 3-17. Sinclair, J. S., S. E. Arnott, W. A. Nelson, & K. B. Brougham, 2021. Consistent and transient drivers of freshwater zooplankton communities. Journal of Biogeography 48: 811–822. Sluss, T. D., G. A. Cobbs, & J. H. Thorp, 2008. Impact of turbulence on riverine zooplankton: a mesocosm experiment. Freshwater Biology 53: 1999–2010. Soininen, J., J. J. Korhonen, J. Karhu, & A. Vetterli, 2011. Disentangling the spatial patterns in community composition of prokaryotic and eukaryotic lake plankton. Limnology and Oceanography 56: 508–520. Taura, Y. M., & I. C. Duggan, 2020. The Relative Effects of Willow Invasion, Willow Control and Hydrology on Wetland Zooplankton Assemblages. Wetlands 40: 2585– 2595. Tundisi, J. G., 2003. Água no século XXI: Enfrentando a escassez. Rima, São Carlos, SP. Tuomisto, H., & K. Ruokolainen, 2006. Analyzing or explaining beta diversity? Understanding the targets of different methods of analysis. Ecology 87: 2697–2708. Vellend, M., 2010. Conceptual Synthesis in Community Ecology. The Quarterly 25 Review of Biology 85: 183–206. Ward, J. V., 1989. The Four-Dimensional Nature of Lotic Ecosystems Journal of the North American Benthological Society 8: 2–8. Zhao, K., K. Song, Y. Pan, L. Wang, L. Da, & Q. Wang, 2017. Metacommunity structure of zooplankton in river networks: Roles of environmental and spatial factors. Ecological Indicators 73: 96–104. Zorzal-Almeida, S., L. M. Bini, & D. C. Bicudo, 2017. Beta diversity of diatoms is driven by environmental heterogeneity, spatial extent and productivity. Hydrobiologia 800: 7–16. 26 CAPÍTULO 1: Correlates of zooplankton metacommunity structure in small reservoirs27 Abstract Different mechanisms act to structure metacommunities, among them, the most studied are those related to environmental and spatial drivers. Also, most studies on metacommunities are based on data obtained in snapshot surveys, hindering the test of other potential mechanisms (e.g., historical effects). Here, we quantified the relative importance of different sets of variables that could potentially structure zooplankton communities in 39 small reservoirs in the Brazil Midwest. In addition to local environmental, spatial and land use variables, we tested the effect of earlier zooplankton composition on (later) zooplankton metacommunity structure. The relative importance of these sets of explanatory variables varied over time. Unexpectedly, land use variables were not correlated to zooplankton community. The local environmental variables that best correlated with zooplankton metacommunity structure were those related to clear/turbid water states. Earlier zooplankton community data were significant predictors (specially for microcrustaceans), indicating the role of historical effects on later zooplankton metacommunity structure. Spatial variables were also often significant. Our results highlight the importance of considering multiple surveys in studies aiming to unveil the determinants of zooplankton metacommunity structure. Keywords: Artificial shallow lakes; Species composition; Historical effects; Variance partitioning 28 Introduction Local communities encompassing a metacommunity are structured by environmental factors, biotic interactions, and dispersal (Leibold et al. 2004, 2010; Vellend 2010; Leibold and Chase 2017). This basic idea assumes that the species of a given regional pool have specific environmental requirements and, therefore, those occurring in local communities have passed through multiple environmental filters (Mittelbach and Schemske 2015). Dispersal, in turn, depending on its intensity, can increase or decrease the relationship between species abundances and environmental gradients (Heino et al., 2015; Thompson et al. 2020). Dispersal rates should be sufficient for species to colonize and develop viable populations in local communities that have suitable environmental characteristics (considering their ecological requirements). On the other hand, when dispersal rates are too high (e.g., through mass effects), the strengths of the relationships between species abundances and environmental gradients tend to be reduced because many species may occur in environmentally suboptimal local communities (Heino et al. 2015). Despite the conceptual advance promoted by metacommunity theory (e.g., Ricklefs 2008), empirical results show that our power to predict species distributions is limited (e.g., Soininen 2014, 2016). Although the low predictive power may be related to stochastic processes, this inference - based on large residual variation - assumes that all relevant environmental variables have been measured and that biotic interactions can be overlooked (Vellend et al. 2014). Thus, much of the unexplained variation may also result from the absence of relevant explanatory variables (i.e., a problem of model misspecification). Furthermore, most studies on metacommunities in aquatic environments are, in general, based on snapshot surveys (but see Eros et al. 2012; 29 Wojciechowski et al. 2017; Sinclair et al. 2021, for some exceptions) and, therefore, they do not consider how previous ecological conditions have influenced current or more recent local communities (Vass and Langenheder 2017; Langenheder and Lindström 2019). Recently, however, some metacommunity studies have considered earlier environmental (Andersson et al. 2014) and biological data (Mergeay et al. 2011; Castillo-Escrivà et al. 2017; Oliveira et al. 2020; Ortega et al. 2021) as potential explanatory factors of the variation in later local communities. The strength of the relationship between earlier and later communities may depend on several factors, such as species generation time, time interval between samplings and temporal environmental variability. For example, Jacquemin and Pyron (2011) found large differences between stream fish communities sampled in 1940 and those sampled later (1996 and 2007) and attributed these differences to anthropogenic impacts (e.g., changes in water quality and species introductions) that occurred between the years. On the other hand, aquatic metacommunities may be temporally autocorrelated (i.e., local communities at time t = 1 are correlated to those at time t + 1) even when large environmental changes (e.g., hydrological variations) occur between sampling times (Oliveira et al. 2020). Similarly, even when long periods are considered, correlations between more recent and past metacommunities can be detected (Mergeay et al. 2011; Castillo-Escrivà et al. 2017). Correlations between the same metacommunity studied over time can emerge from priority effects (Fukami 2015). These effects occur when the first species that colonize a local community develop higher abundances and, therefore, preempt resources and restrict the establishment of late-arriving species (Mergeay et al. 2011). Even after severe environmental disturbances that can reduce the abundances of early 30 colonists and that, consequently, would open windows of opportunity for other species, the priority effects can be maintained when the species of the initial communities develop dormant propagule banks (e.g., resting egg banks; Hairston 1996; De Meester et al. 2002; Mergeay et al. 2011; Castillo-Escrivà et al. 2017). Thus, early-arriving species would have advantages as they would quickly reestablish their populations from the propagule banks. Alternatively, local communities assessed successively over time may be similar because strong structuring environmental gradients are maintained or re- established between samplings. For example, strong environmental gradients across a watershed (from small headwater to higher order streams) that influence fish communities tend to be maintained over time (Oliveira et al. 2020). In zooplankton communities, assuming temporal changes in environmental gradients, short generation times (Wetzel 1975; Allan 1976; Gillooly 2000) and rapid responses to environmental changes (Jeppesen et al. 2011), low correlations between data from local communities obtained in successive samplings can be expected. Empirically, these changes can be verified, for example, with an ordination analysis that would show a low juxtaposition of the scores of the (same) sampling units during two or more sampling events, considering both environmental and zooplankton data (e.g., Retting et al. 2006; Li et al. 2019; Picapedra et al. 2020). On the other hand, many zooplankton species can form a resting egg bank (Hairston 1996; Gyllström and Hansson 2004; Havel and Shurin 2004) and their communities can be structured by priority effects despite changes in environmental gradients. In these cases, the results from ordination analyses would show high congruence. Zooplankton communities can also be influenced by variables obtained at the landscape scale, in addition to those obtained locally. With the intensification of land use (i.e., agriculture, pasture, urban areas), natural and artificial aquatic environments 31 tend to receive high nutrient loads, sediments, and pesticides (Dodson et al. 2005; Van Egeren et al. 2011). Despite their socioeconomic importance and high biodiversity (Czerniawski and Domagała 2014; Czerniawski and Kowalska-Góralska 2018), small reservoirs (i.e., reservoirs used in agriculture) are especially impacted by activities carried out in adjacent terrestrial ecosystems (Joniak and Kuczyńska-Kippen 2010; Chen et al. 2019). The small size, depth, and thepossibility to receive nutrient-rich effluents, make small reservoirs prone to eutrophication, which favors the dominance of a few tolerant zooplankton species (e.g., Le Quesne et al. 2020). Here, we quantified the relative importance of different groups of variables (local scale environment, land use, spatial variables, and earlier zooplankton community data) that could structure zooplankton communities (microcrustaceans, rotifers and testate amoebae) in small reservoirs. Zooplankton groups are considered excellent ecological indicators (Jeppesen et al. 2011). Therefore, we predicted that local environmental factors would be significantly correlated with zooplankton community structure. However, because land use variables may be strongly correlated with water quality, we expected they would also be important in structuring local zooplankton communities (Van Egeren et al. 2011; Santos et al. 2021). Finally, we expected a relationship between data of local communities that were obtained in two consecutive times (i.e., the metacommunity structure at time t would be similar to that observed at time t+1). This expectation can be justified considering, among others, the maintenance of spatial environmental gradients between two sampling periods and the role of priority effects. For example, microcrustaceans and rotifers can develop large resting egg banks that can quickly, after disturbances (e.g., alternating dry and rainy periods), determine the composition of local communities (De Meester et al. 2002; Allen et al. 2011). However, experimental studies have shown that the effect of resistance egg bank on 32 local community structures over time may be greater for microcrustaceans than for rotifers (Lopes et al. 2016). Thus, considering the different zooplanktonic groups, we expected local communities at an earlier time to be correlated with local community structures at a later time especially for microcrustaceans and, to a lesser degree, for rotifers. Methods Study area This study was carried out in small reservoirs (i.e., artificial shallow lakes) located in the Rio Preto watershed (Federal District and Goiás State; Brazil; Figure 1). The Alto Rio Preto Microbasin (MARP) belongs to the São Francisco River Basin and has an area of 3,470 km², with 2,344.045 km of drainage (Schrage and Uagoda 2017). The Rio Preto watershed is characterized by flat terrains, which facilitates the intensive use of the landscape (Figure S1). The small reservoirs are mainly used for irrigation. The climatic characteristics of the studied area are typical of the Brazilian Midwest, with two well-defined seasons: a rainy one (from October to April) and a dry one (from May to September; Bustamante et al. 2012). The sampling campaigns were carried out during dry and rainy seasons. The number of ponds analyzed varied because some were dry or inaccessible in some months. Thus, 23 small reservoirs were sampled in 05/2017 (dry), 39 in 02/2018 (rainy), 38 in 09/2018 (dry), and 38 in 02/2019 (rainy). 33 Figure 1. Small reservoirs sampled in the Rio Preto watershed (MARP), Distrito Federal and Goiás (GO), Brazil. Zooplankton At each reservoir, we collected zooplankton samples by filtering 1,000 L of water through a plankton net with a 68 µm mesh size. Subsequently, the samples were fixed in a 4% formaldehyde solution and buffered with sodium tetraborate. For quantitative and qualitative analyses, zooplankton samples were concentrated to a known volume (minimum 50 mL and maximum depending on the number of organisms and amount of sediment in the sample). We used a Hensen-Stempel pipette to remove 10 % of the total volume of each sample and, using Sedgwick-Rafter chambers, we counted each aliquot under an optical microscope (Olympus CX31 - 400x); however, samples with low densities were fully counted (Bottrell et al. 1976). Qualitative analysis 34 was performed by counting species that were not registered in the quantitative analysis. The microcrustaceans (copepods and cladocerans), rotifers and testate amoebae were identified to the lowest possible taxonomic level (usually species) using specialized taxonomic keys (Koste 1978; Vucetich 1978; Ogden and Hedley 1980; Reid 1985; Matsumura-Tundisi 1986; Paggi 1995; Velho and Lansac-Tôha 1996; El Moor-Loureiro 1997; Gomes-Sousa 2008). Densities were expressed as individuals/m3. Local environmental variables The following local environmental variables were obtained using a Horiba® water quality meter: temperature, turbidity, pH, dissolved oxygen, and conductivity. Water transparency and depth were measured with a Secchi disk. We recorded presence-absence data of aquatic macrophytes (emergent/floating leaves, submerged, and free floating) in each small reservoir. Water samples were analyzed for orthophosphate, total phosphorus, inorganic nitrogen, and total nitrogen using standard methods (APHA 2005; ASTM 2016). However, not all variables were measured in all sampling campaigns (Table S1). The dimensions of the reservoirs (perimeter (km) and area (km²)) were obtained using data from Google Earth®. We also calculated the shoreline development index (Häkanson 2004). Land use variables We delimited the watershed area upstream of each small reservoir using the Digital Elevation Model (DEM) of the study area (Advanced Land Observing Satellite, 35 Phased Array L-band Synthetic Radar, https://asf.alaska.edu/). After, we used data from the MapBiomas Project (https://mapbiomas.org) in ArcMap® to obtain the following variables: watershed area (km2), remaining vegetation (%), pasture (%), agriculture (%), agriculture/pasture mosaic (%), and urban occupation (%). Data analysis All analyzes were performed using the R software (R Core Team 2021). The relationships between community abundance data at time t and our set of explanatory matrices (local environment, land use, space, and community data at time t-1; see below) were tested using partial Redundancy Analysis (pRDA) and variance partitioning (as described in Peres-Neto et al. 2006). For these analyses, we used the functions “rda” and “varpart” of the vegan package (Oksanen et al. 2020), respectively. The explanatory variables of each set were selected, before variation partitioning, following the method described in Blanchet et al. (2008) and Bauman et al. (2018a). To do so, we used the functions “mem.sel” and “forward.sel” from the adespatial package (Dray et al. 2021). Community abundance data were log-transformed, and Hellinger standardized prior to analyses. As mentioned above, the abundance community data at time t-1were used as an explanatory matrix of the zooplankton community data at time t. Thus, the community data gathered, for example, in 05/2017, were used as explanatory variables of the community data observed in 02/2018 and so on, successively, for the other pairs of sampling dates. However, owing to the high dimensionality of the community data, we first summarized the data obtained at time t-1 using a Bray-Curtis based Principal Coordinate Analysis (PCoA). Then, the PCoA axes scores, calculated with the https://asf.alaska.edu/ https://mapbiomas.org/ 36 “cmdscale” function of the “stats” package (R Core Team 2021), were used as candidate variables (i.e., as potential explanatory variables of the zooplankton community obtained in the most recent campaign). We used the geographic coordinates of the small reservoirs to generate Moran’s Eigenvector Maps (MEM) (Bauman et al. 2018b), those with positive correlations were then used as spatial predictors. Different connectivity and weighting matrices were tested (Table S2). This procedure was used to model the spatial component in the best possible way. For that, we used the functions “listw.candidates”, “listw.select”from the “adespatial” package (Dray et al. 2021). All significance tests were based on 9999 permutations. Results We found a total of 285 zooplankton taxa during the four sampling campaigns (95 rotifers, 112 testate amoebae, and 78 microcrustaceans, including the larval and juvenile forms of copepods) (Table S3). Among rotifers, Lecane Bulla (63.31%), Keratella americana (51.08%), L. leontina (43.17%), Polyarthra vulgaris (43.17%) and L. signifera (39.57%) were the most frequent. As for microcrustaceans, immature copepods (nauplii and juveniles of Cyclopidae and Diaptomidae) were highly frequent (57.55% up to 84.89%), as well as Ilyocryptus spinifer (46.05%). Finally, Arcella vulgaris (74.1%), A. costata (63.31%), Lesquereusia spiralis (61.15%), Centropyxis spinosa (51.8%), and Netzelia corona (51 .8%) were the most frequent taxa among testate amoebae (Table S3). The numbers of taxa occurring in only one sampling campaign were 33 (34.73%), 19 (24 .35%) and 33 (29.46%) for rotifers, microcrustaceans and testate amoebae, respectively. 37 Turbidity, pH, depth, water transparency, and macrophytes (submerged and emergent) were often selected by the forward procedure. Among land use, the most often selected variables were watershed area and the percentages of agricultural and urbanized areas in the landscape. The first PCoAs axis based on past communities was selected most of the time. Spatial variables (MEMs) were not selected for zooplankton communities in 05/2017, rotifers in 09/2018, testate amoebae in 05/2017, and microcrustaceans in 05/2017. In general, MEMs representing large spatial scales were the most often selected (Table S4, S5, S6, and S7). Variation in zooplankton community structure (i.e., considering the whole community) in 05/2017 was significantly correlated with the local environmental variables only (adjusted R² = 8.7%; P = 0.001). In 02/2018, the variation in the zooplankton community was explained by spatial variables (adjusted R² = 4.5%; P = 0.018) and by the zooplankton community observed in 05/2017 (i.e., PCoA axes obtained with the community data in 05/2017; adjusted R² = 4.3%; P = 0.024). For the zooplankton community in 09/2018, the spatial and environmental fraction were significant (adjusted R² = 8.4%; P = 0.001 and R² = 0,026; P = 0.001, respectively). The spatial (adjusted R² = 3.1; P = 0.034) and past community variables (adjusted R² = 4.9%; P = 0.003) were significantly correlated with the zooplankton community data recorded in 02/2019 (Figure 2). 38 Figure 2. Relative contributions (% of variance accounted) of spatial (S), local environmental (LE), early community (EC) and land use (LU) variables in the structuring the zooplankton metacommunity of the Rio Preto watershed. Values in bold indicate that the components were significant (< 0.05%). We did not use a set of explanatory variables in variation partitioning when their variables were not previously selected. Components of Venn diagrams without values (empty) indicate values close to or <0. Variation in rotifer community structure in 05/2017 was explained by spatial (adjusted R² = 8.1%; P = 0.001) and local environmental variables (adjusted R² = 12.3%; P = 0.001). Spatial variables (adjusted R² = 13.2 %; P = 0.001) were 39 significantly correlated with the rotifer community recorded in 02/2018. No variables were significantly correlated with the variation in rotifer community structure in 09/2018. In 02/2019, only the local environmental fraction was significant (adjusted R² = 3.1%; P = 0.038; Figure 2). Testate amoebae community data in 05/2017 and 02/2018 were significantly correlated with local environmental (adjusted R² = 16.9%; P = 0.001) and spatial variables (adjusted R² = 12.3%; P = 0.001), respectively. In 09/2018, the spatial (adjusted R² = 5.8%; P = 0.001) and local environmental (adjusted R² = 4.1%; P = 0.001) fractions were significant. The spatial (adjusted R² = 4.4%; P = 0.005) and local environmental (adjusted R² = 7.3 %; P = 0.002) fractions were also significant for the testate amoebae data recorded in 02/2019 (Figure 2). The variation in the microcrustacean community structure in 05/2017 were significantly correlated with the local environmental variables (adjusted R² = 3.1 %; P = 0.042). In 02/2018, the microcrustacean data were correlated with the spatial variables (adjusted R² = 9.1%; P = 0.001) and with the PCoA axes extracted from the communities recorded in the previous sampling campaign (adjusted R² = 8.5%; P = 0.001). In 09/2018, the variation in the microcrustacean community was significantly explained only by the community recorded in 02/2018 (adjusted R² = 6.8%, P = 0.009). In 02/2019, the spatial variables (adjusted R² = 4.2%, P = 0.03) and the PCoA axes from the microcrustacean community data recorded in 09/2018 (adjusted R² = 9.8%; P = 0.001) explained significant portions of the variation among local microcrustacean communities (Figure 2). Discussion 40 We found that the relative importance of the different sets of explanatory variables varied among the sampling campaigns. In general, only land use variables were not significantly correlated with zooplankton community structure in any sampling campaign. On the other hand, spatial, local environmental variables, and the community data obtained at an earlier time were more often correlated with the community structure. Other studies on aquatic metacommunities also found that variation partitioning fractions varied through time (Eros et al. 2012; Erős et al. 2014; Ortega et al. 2021; Li et al. 2021). The temporal variation in the relative importance of explanatory variables suggests that mechanisms underlying zooplankton community structure change over time and that analyses based on snapshot surveys may overlook temporally variable mechanisms (Sinclair et al. 2021; but see Cottenie et al. 2003). In general, pH, turbidity, transparency, depth, and aquatic macrophytes were the main variables that contributed to the pure local environmental fraction. Several studies have demonstrated the influence of these variables in structuring zooplankton communities in different aquatic ecosystems (Jeppesen et al. 1997; Cottenie et al. 2003; Dejen et al. 2004; Shurin et al. 2010; Zhao et al. 2017; Sinclair et al. 2021). The importance of pH in structuring zooplankton communities, for example, has been emphasized for different types of aquatic ecosystems (Shurin et al. 2010; Gray and Arnott 2011), including water bodies used in agriculture (Le Quesne et al. 2020). However, despite the predominance of low pH (< 6.5 on average) in our study area – a condition generally associated with low zooplanktonic species richness – the total species richness was high. These results demonstrate the importance of the small reservoirs for the maintenance of zooplankton biodiversity (see also Le Quesne et al. 2020). 41 The local environmental variables that were important in structuring the zooplankton communities (i.e., turbidity, transparency, depth, and aquatic macrophytes) are generally related to feedback mechanisms that determine the alternative stable states in shallow aquatic environments (Scheffer et al. 1993). Thus, similarly to what was observed in previous studies (Cottenie et al. 2001, 2003), spatial variations in turbidity and in other variables that characterize the states of clear (with submerged aquatic macrophytes) and turbid (without submerged aquatic macrophytes) waters may have determined the spatial variations in zooplankton community structure. The frequent response of testate amoebae local communities to environmental gradients was, to some extent, an unexpected result. In general, the analysis of testate amoebae in zooplankton studies is comparatively rare. Therefore, it isnot yet possible to assess whether this group can be considered an ecological indicator as reliable as the traditionally evaluated zooplankton groups (e.g., microcrustaceans and rotifers). Results obtained in other studies, however, indicate that testate amoebae can be used to infer environmental changes (Velho et al. 2003; Mitchell et al. 2008; Arrieira et al. 2015; Qin et al. 2016; Schwind et al. 2017; Nasser et al. 2019; Wang et al. 2020). The evidence regarding the effects of land-use variables on aquatic communities, as well as their relative predictive power in comparison with local-scale variables (e.g., pH and water transparency), is ambivalent (e.g., compare results from Dodson et al. 2005; Van Egeren et al. 2011; Zorzal-Almeida et al. 2017; Firmiano et al. 2020 with those from Mantyka-Pringle et al. 2014; Machado et al. 2016; Barbosa et al. 2019; Rocha et al. 2020; Santos et al. 2021). Land-use variables were not significantly correlated with zooplankton community structure in our study. This result was unexpected considering the large variation in land use (Figure S1). Thus, at least in our 42 study, the evidence points to the primacy of local-scale variables as predictors of zooplankton structure as compared to land-use variables. The spatial fraction was significant in most of the analyzes performed in this study. Similar results were observed in other studies (Beisner et al. 2006; De Bie et al. 2012; Gomes et al. 2020; Rocha et al. 2020). However, the mechanisms associated with a significant spatial fraction cannot be easily inferred (e.g., Vellend et al. 2014; Brown et al. 2017). For example, in addition to processes related to dispersal (e.g., Cottenie 2005), the effects of spatially autocorrelated, but unmeasured, environmental variables on the structuring of metacommunities cannot be ruled out (Peres-Neto and Legendre 2010; Diniz-Filho et al. 2012). However, studies have shown that the average size of propagules can predict the relative importance of environmental and spatial variables. Particularly, the dispersal limitation related spatial processes would be less important for organisms with smaller propagule sizes (De Bie et al. 2012; Padial et al. 2014; Martin et al. 2021). However, we did not find support for this prediction since the effects of spatial variables were similar for microcrustaceans and rotifers. Considering the small spatial distances between the small reservoirs (longest distance equal to 60 km), in general, it is unlikely that the significant spatial fraction detected in this study indicates that dispersal limitation is a major factor in structuring the zooplankton metacommunity (e.g., Shurin 2000; Havel and Shurin 2004). Even acknowledging scale effects, the importance of dispersal as a structuring process of zooplankton communities is still a controversial issue considering the available literature (e.g., Gray and Arnott 2011). Metacommunity studies are predominantly based on data obtained at a single time. This widely used sampling design, however, has been criticized. For example, the mechanisms that structure metacommunities may vary over time (Brown et al. 2017; 43 Zhao et al. 2017; Sinclair et al. 2021) or local communities observed at time t may reflect past environmental conditions (Fischer et al. 2001; Andersson et al. 2014; Vass and Langenheder 2017) or be a result of priority effects (Fukami 2015). Consistent with these considerations, the zooplankton dataset obtained at time t-1 was an important predictor of the metacommunity at time t (especially for microcrustaceans). Previous studies have proposed that priority effects could be investigated by quantifying the effects of past environmental conditions on local communities (Andersson et al. 2014; Vass and Langenheder 2017). However, using data from local communities themselves obtained at an earlier time, as done in this study, seems to be a more direct alternative (see also Mergeay et al. 2011; Castillo-Escrivà et al. 2017; Oliveira et al. 2020; Ortega et al. 2021). Assuming strong environmental effects on local communities and temporal invariability of environmental gradients, the effect of earlier community data would emerge simply through the action of niche-based mechanisms (which act similarly in the two analyzed times). However, if the spatial fraction detected in this study is not related to unmeasured environmental variables, this explanation seems unlikely since the local environmental fraction was not preponderant. The priority effects, where the first colonizers dominate the composition of local communities (Fukami 2015), could also explain the fraction associated with the earlier community data. In turn, priority effects can be maximized by a storage effect, for example, by a large propagule bank (Mergeay et al. 2011). For rotifers and testate amoebae, the composition of previous communities did not influence the succeeding metacommunities. These results suggest a high turnover in composition. Although we confirmed our hypothesis that stronger correlations between local communities at earlier and later times would be observed for microcrustaceans in 44 comparison to rotifers (considering, for example, the experimental results of Lopes et al. 2016), the fact that we were unable to detect a relationship for the latter group was surprising. This is so because several studies indicate that propagule bank sizes of rotifers are much higher than those of microcrustaceans (e.g., Frisch et al. 2009 and references in Lopes et al. 2016). We speculate that the benthic-pelagic decoupling for rotifers (as inferred by the absence of relationship between community structure during consecutive sampling campaigns) may be caused by the stronger effects of aerial dispersal on local community structure (which indeed tend to be stronger for rotifers than for microcrustaceans). In conclusion, our main results indicate the importance of including past community as potential predictors of zooplankton community structure in small reservoirs. Including the earlier community data as additional predictors in multivariate canonical analyses (Mergeay et al. 2011; Castillo-Escrivà et al. 2017; Oliveira et al. 2020; Ortega et al. 2021) is a straightforward way to measure the potential impacts of historical effects. However, for this to be possible, we echo previous calls (Langenheder et al. 2012; Eros et al. 2012) for the need of considering the temporal dimension in studies aiming to evaluate the determinants metacommunity structure. Acknowledgements This study was partially financed by the Fundação de Apoio a Pesquisa do Distrito Federal (FAPDF) (proc. 6067/2013-3), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES; scholarships to MCV and JCGO; Finance Code 001) and the Brazilian Council of Research (CNPq; grants to LMB and LCGV). This work has been developed in the context of Instituto Nacional de Ciência e Tecnologia (INCT) in 45 Ecology, Evolution and Biodiversity Conservation, supported by MCTIC/CNpq (proc. 465610/2014-5) and Fundação de Amparo a Pesquisa de Goiás (FAPEG). We thanks to Hasley R. Pereira, Leonardo B. Silva, Rafaela V. Granzotti, Gustavo F. Granjeiro, João Paulo A. Motta and Thallia S. Silva for helping with sampling. We declare that our manuscript complies with the current Brazilian laws. References Allan JD (1976) Life History Patterns in Zooplankton. Am Nat 110:165–180 Allen MR, VanDyke JN, Cáceres CE (2011) Metacommunity assembly and sorting in newly formed lake communities. Ecology 92:269–275 Andersson MGI, Berga M, Lindstrom ES, Langenheder S (2014) The spatial structure of bacterial communities is influenced by historical environmental conditions. Ecology 95:1134–1140 APHA (2005) Standard methods for the examination of waterand waste water. American P, Washington Arrieira RL, Schwind LTF, Bonecker CC, Lansac-Tôha FA (2015) Use of functional diversity to assess determinant assembly processes of testate amoebae community. Aquat Ecol 49:561–571 ASTM D8083-16 (2016) ASTM (American Society for Testing and Materials). In: Water by High Temperature Catalytic Combustion and Chemiluminescence Detection, ASTM. West Conshohocken, PA Barbosa H de O, Borges PP, Dala-Corte RB, et al (2019) Relative importance of local and landscape variables on fish assemblages in streams of Brazilian savanna. Fish Manag Ecol 26:119–130 Bauman D, Drouet T, Dray S, Vleminckx J (2018a) Disentangling good from bad 46 practices in the selection of spatial or phylogenetic eigenvectors. Ecography 41: 1638-1649 Bauman D, Drouet T, Fortin M-J, Dray S (2018b) Optimizing the choice of a spatial weighting matrix in eigenvector-based methods. Ecology 99:2159–2166 Beisner BE, Peres-Neto PR, Lindström ES, et al (2006) The role of environmental and spatial processes in structuring lake communities from bacteria to fish. Ecology 87:2985–2991 Blanchet FG, Legendre P, Borcard D (2008) Forward selection of explanatory variables. Ecology 89:2623–2632 Bottrell HH, Ducan A, Gliwicz ZM, et al (1976) Review of some problems in zooplankton production studies. Nor J Zool 24:419–456 Brown BL, Sokol ER, Skelton J, Tornwall B (2017) Making sense of metacommunities: dispelling the mythology of a metacommunity typology. Oecologia 183:643–652 Bustamante MMC, Nardoto GB, Pinto AS, et al (2012) Os impactos potenciais das mudanças climáticas no funcionamento biogeoquímico dos ecossistemas de cerrado (in Portuguese). Brazilian J Biol 72:655–671 Castillo-Escrivà A, Valls L, Rochera C, et al (2017) Disentangling environmental, spatial, and historical effects on ostracod communities in shallow lakes. Hydrobiologia 787:61–72 Chen W, He B, Nover D, et al (2019) Farm ponds in southern China: Challenges and solutions for conserving a neglected wetland ecosystem. Sci Total Environ 659:1322–1334 Cottenie K (2005) Integrating environmental and spatial processes in ecological community dynamics. Ecol Lett 8:1175–1182 Cottenie K, Michels E, Nuytten N, Meester L De (2003) Zooplankton metacommunity 47 structure : regional vs . local processes in highly interconnected ponds. Ecology 84:991–1000 Cottenie K, Nuytten N, Michels E, De Meester L (2001) Zooplankton community structure and environmental conditions in a set of interconnected ponds. Hydrobiologia 442:339–350 Czerniawski R, Domagała J (2014) Small dams profoundly alter the spatial and temporal composition of zooplankton communities in running waters. Int Rev Hydrobiol 99:300–311 Czerniawski R, Kowalska-Góralska M (2018) Spatial changes in zooplankton communities in a strong human-mediated river ecosystem. PeerJ 2018:e5087 De Bie T, De Meester L, Brendonck L, et al (2012) Body size and dispersal mode as key traits determining metacommunity structure of aquatic organisms. Ecol Lett 15:740–747 De Meester L, Gómez A, Okamura B, Schwenk K (2002) The Monopolization Hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecologica 23:121–135 Dejen E, Vijverberg J, Nagelkerke LAJ, Sibbing FA (2004) Temporal and spatial distribution of microcrustacean zooplankton in relation to turbidity and other environmental factors in a large tropical lake (L. Tana, Ethiopia). Hydrobiologia 513:39–49 Diniz-Filho JAF, Siqueira T, Padial AA, et al (2012) Spatial autocorrelation analysis allows disentangling the balance between neutral and niche processes in metacommunities. Oikos 121:201–210 Dodson SI, Lillie RA, Will-Wolf S (2005) Land use, water chemistry, aquatic vegetation, and zooplankton community structure of shallow lakes. Ecol Appl 48 15:1191–1198 Dray S, Bauman D, Blanchet G, et al (2021) adespatial: Multivariate Multiscale Spatial Analysis. R package version 0.3-14. https://CRAN.R- project.org/package=adespatial El Moor-Loureiro LMA (1997) Manual de identificação de cladóceros límnicos do Brasil (in Portuguese). Ed. Univer. Brasília, Brasília Eros T, Sály P, Takács P, et al (2012) Temporal variability in the spatial and environmental determinants of functional metacommunity organization – stream fish in a human-modified landscape. Freshw Biol 57:1914–1928 Erős T, Sály P, Takács P, et al (2014) Quantifying temporal variability in the metacommunity structure of stream fishes: The influence of non-native species and environmental drivers. Hydrobiologia 722:31–43 Firmiano KR, Cañedo‐Argüelles M, Gutiérrez‐Cánovas C, et al (2020) Land use and local environment affect macroinvertebrate metacommunity organization in Neotropical stream networks. J Biogeogr 00:jbi.14020 Fischer JM, Klug JL, Ives AR, Frost TM (2001) Ecological history affects zooplankton community responses to acidification. Ecology 82:2984–3000 Frisch D, Arechederra A, Green AJ (2009) Recolonisation potential of zooplankton propagule banks in natural and agriculturally modified sections of a semiarid temporary stream (Doñana, Southwest Spain). Hydrobiologia 624:115–123 Fukami T (2015) Historical Contingency in Community Assembly: Integrating Niches, Species Pools, and Priority Effects. Annu Rev Ecol Evol Syst 46:1–23 Gillooly JF (2000) Effect of body size and temperature on generation time in zooplankton. J Plankton Res 22:241–251 Gomes-Sousa MB (2008) Guia das Tecamebas - Bacia do Rio Peruaçu - Minas Gerais: 49 subsídio para a conservação e monitoramento da Bacia do Rio São Francisco (in Portuguese). Editora UFMG, Belo Horizonte Gomes LF, Barbosa JC, de Oliveira Barbosa H, et al (2020) Environmental and spatial influences on stream zooplankton communities of the Brazilian Cerrado. Community Ecol 21:25–31 Gray DK, Arnott SE (2011) Does dispersal limitation impact the recovery of zooplankton communities damaged by a regional stressor? Ecol Appl 21:1241– 1256 Gyllström M, Hansson LA (2004) Dormancy in freshwater zooplankton: Induction, termination and the importance of benthic-pelagic coupling. Aquat Sci 66:274–295 Hairston NG (1996) Zooplankton egg banks as biotic reservoirs in changing environments. Limnol Oceanogr 41:1087–1092 Häkanson L (2004) Lakes: Form and Function. Blackburn Press, Caldwell Havel JE, Shurin JB (2004) Mechanisms, effects, and scales of dispersal in freshwater zooplankton. Limnol Oceanogr 49:1229–1238 Heino J, Melo AS, Siqueira T, et al (2015) Metacommunity organisation, spatial extent and dispersal in aquatic systems: Patterns, processes and prospects. Freshw Biol 60:845–869 Jacquemin SJ, Pyron M (2011) Fishes of Indiana streams: Current and historic assemblage structure. Hydrobiologia 665:39–50 Jeppesen E, Jensen JP, Søndergaard M, et al (1997) Top-down control in freshwater lakes: the role of nutrient state, submerged macrophytes and water depth. Hydrobiologia 342:151–164 Jeppesen E, Nõges P, Davidson TA, et al (2011) Zooplankton as indicators in lakes: A scientific-based plea for including zooplankton in the ecological quality assessment 50 of lakes according to the European Water Framework Directive (WFD). Hydrobiologia 676:279–297 Joniak T, Kuczyńska-Kippen N (2010) The chemistry of water and bottom sediments in relation to zooplankton biocenosis in small agricultural ponds. Oceanol Hydrobiol Stud 39:85–96 Koste W (1978) Rotatoria: Die Radertiere Mitteuropas. Gebruder Borntraeger, Berlin Langenheder S, Berga M, Östman Ö, Székely AJ (2012) Temporal variation of β- diversity and assembly mechanisms in a bacterial metacommunity. ISME J 6:1107–1114 Langenheder S, Lindström ES (2019) Factors influencing aquatic and terrestrial bacterial community assembly. Environ Microbiol Rep 11:306–315 Le Quesne KS, Özkundakci D, Duggan IC (2020) Life on the farm: are zooplanktoncommunities in natural ponds and constructed dams the same? Mar Freshw Res 72:330-340. Leibold MA, Chase JM (2017) Metacommunity Ecology. Princeton University Press, Princeton Leibold MA, Economo EP, Peres-Neto P (2010) Metacommunity phylogenetics: Separating the roles of environmental filters and historical biogeography. Ecol Lett 13:1290–1299 Leibold MA, Holyoak M, Mouquet N, et al (2004) The metacommunity concept: A framework for multi-scale community ecology. Ecol Lett 7:601–613 Li C, Feng W, Chen H, et al (2019) Temporal variation in zooplankton and phytoplankton community species composition and the affecting factors in Lake Taihu-a large freshwater lake in China. Environ Pollut 245:1050–1057 Li Z, Heino J, Song Z, et al (2021) Spatio-temporal variation of macroinvertebrate 51 metacommunity organization in a monsoon-climate region. J Biogeogr 48:3118– 3130 Lopes PM, Bozelli R, Bini LM, et al (2016) Contributions of airborne dispersal and dormant propagule recruitment to the assembly of rotifer and crustacean zooplankton communities in temporary ponds. Freshw Biol 61:658–669 Machado KB, Teresa FB, Vieira LCG, et al (2016) Comparing the effects of landscape and local environmental variables on taxonomic and functional composition of phytoplankton communities. J Plankton Res 38:1334–1346 Mantyka-Pringle CS, Martin TG, Moffatt DB, et al (2014) Understanding and predicting the combined effects of climate change and land-use change on freshwater macroinvertebrates and fish. J Appl Ecol 51:572–581 Martin GK, Beisner BE, Chain FJJ, et al (2021) Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation. Ecology 102:e03224 Matsumura-Tundisi T (1986) Latitudinal distribution of Calanoida copepods in freshwater aquatic systems of Brazil. Rev Bras Zool 46:527–553 Mergeay J, De Meester L, Eggermont H, Verschuren D (2011) Priority effects and species sorting in a long paleoecological record of repeated community assembly through time. Ecology 92:2267–2275 Mitchell EAD, Charman DJ, Warner BG (2008) Testate amoebae analysis in ecological and paleoecological studies of wetlands: past, present and future. Biodivers Conserv 17:2115–2137 Mittelbach GG, Schemske DW (2015) Ecological and evolutionary perspectives on community assembly. Trends Ecol Evol 30:241–247 Nasser NA, Cullen J, Patterson CW, et al (2019) Inter-annual Arcellinida (testate lobose 52 amoebae) assemblage dynamics within lacustrine environments. Limnologica 76:60–71 Ogden CG, Hedley RH (1980) An Atlas of Freshwater Testate Amoebae. Oxford University Press, London Oksanen J, Blanchet FG, Friendly M, et al (2020) Package “vegan” Community Ecology Package. R package version 2.5-7. https://CRAN.R- project.org/package=vegan Oliveira FJM, Lima-Junior DP, Bini LM (2020) Current environmental conditions are weak predictors of fish community structure compared to community structure of the previous year. Aquat Ecol. 54, 729-740 Ortega JCG, Geijer J, Bergsten J, et al (2021) Spatio-temporal variation in water beetle assemblages across temperate freshwater ecosystems. Sci Total Environ 792:148071 Padial AA, Ceschin F, Declerck SAJ, et al (2014) Dispersal ability determines the role of environmental, spatial and temporal drivers of metacommunity structure. PLoS One 9:1–8 Paggi JC (1995) Crustacea Cladocera. In: Lopretto EC, Tell G (eds) Ecosistemas de aguas continentales: metodologias para su estudio (in Spanish). Ediciones Sur, La Plata, pp 909–951 Peres-Neto PR, Legendre P (2010) Estimating and controlling for spatial structure in the study of ecological communities. Glob Ecol Biogeogr 19:174–184 Peres-Neto PR, Legendre P, Dray S, Borcard D (2006) Variation partitioning of species data matrices: Estimation and comparison of fractions. Ecology 87:2614–2625 Picapedra PHS, Fernandes C, Taborda J, Baumgartner G, Sanches PV (2020) A long- term study on zooplankton in two contrasting cascade reservoirs (Iguaçu River, 53 Brazil): effects of inter-annual, seasonal, and environmental factors. PeerJ 8:e8979 Qin Y, Payne R, Yang X, et al (2016) Testate amoebae as indicators of water quality and contamination in shallow lakes of the Middle and Lower Yangtze Plain. Environ Earth Sci 75:1–11 R Core Team (2021) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna. https://www.R-project.org/ Reid JW (1985) Chave de identificação e lista de referências bibliográficas para as espécies continentais sulamericanas de vida livre da ordem Cyclopoida (Crustacea, Copepoda). Bol Zool 9:17–143 Rettig JE, Schuman LS, McCloskey JK (2005) Seasonal Patterns of Abundance: Do Zooplankton in Small Ponds do the Same Thing Every Spring–Summer? Hydrobiologia, 556, 193-207 Ricklefs RE (2008) Disintegration of the Ecological Community. Am Nat 172:741–750 Rocha B da S, de Souza CA, Machado KB, et al (2020) The relative influence of the environment, land use, and space on the functional and taxonomic structures of phytoplankton and zooplankton metacommunities in tropical reservoirs. Freshw Sci 39:321–333 Santos EF, Abra J, Castilho-Noll MSM (2021) Does land use influence the local and regional structure of the rotifer assemblage? Hydrobiologia 848:1059–1072 Scheffer M, Hosper SH, Meijer ML, et al (1993) Alternative equilibria in shallow lakes. Trends Ecol Evol 8:275–279 Schrage TJ, Uagoda RES (2017) Distribuição espacial de depressões na bacia do Alto rio Preto (GO, DF, MG) e suas relações com controles geológicos e pedogeomorfológicos (in Portuguese). Rev Bras Geomorfol 18: 379-395 Schwind LTF, Arrieira RL, Simões NR, et al (2017) Productivity gradient affects the 54 temporal dynamics of testate amoebae in a neotropical floodplain. Ecol Indic 78:264–269 Shurin JB (2000) Dispersal limitation, invasion resistance, and the structure of pond zooplankton communities. Ecology 81:3074–3086 Shurin JB, Winder M, Adrian R, et al (2010) Environmental stability and lake zooplankton diversity - contrasting effects of chemical and thermal variability. Ecol Lett 13:453–463 Sinclair JS, Arnott SE, Nelson WA, Brougham KB (2021) Consistent and transient drivers of freshwater zooplankton communities. J Biogeogr 48:811–822 Soininen J (2014) A quantitative analysis of species sorting across organisms and ecosystems. Ecology 95:3284–3292 Soininen J (2016) Spatial structure in ecological communities - a quantitative analysis. Oikos 125:160–166 Thompson PL, Guzman LM, De Meester L, et al (2020) A process-based metacommunity framework linking local and regional scale community ecology. Ecol Lett 23:1314–1329 Van Egeren SJ, Dodson SI, Torke B, Maxted JT (2011) The relative significance of environmental and anthropogenic factors affecting zooplankton community structure in Southeast Wisconsin Till Plain lakes. Hydrobiologia 668:137–146 Vass M, Langenheder S (2017) The legacy of the past: Effects of historical processes on microbial metacommunities. Aquat Microb Ecol 79:13–19 Velho LFM, Lansac-Tôha FA (1996) Testate Amoebae (Rhizopodea-Sarcodina) from Zooplankton of the High Paraná River Floodplain, State of Mato Grosso do Sul, Brazil: II. Family Difflugidae. Stud Neotrop Fauna Environ 31:179–192 Velho LFM, Lansac-Tôha FA, Bini LM (2003) Influence of Environmental 55 Heterogeneity on the Structure of Testate Amoebae (Protozoa, Rhizopoda) Assemblages in the Plankton of the Upper Paraná River Floodplain, Brazil. Int Rev Hydrobiol 88:154–166 Vellend M (2010) Conceptual Synthesis in Community Ecology. Q Rev Biol 85:183– 206 Vellend M, Srivastava DS, Anderson KM, et al (2014) Assessing the relative importance of neutral stochasticity in ecological communities. Oikos 123:1420– 1430 Vucetich MC (1978) Estudio de tecamebianos argentinos, en