News

Jun 13, 2023

Climatically promoted taxonomic homogenization of macroinvertebrates in unaffected streams varies along the river continuum

Scientific Reports volume 13, Article number: 6292 (2023) Cite this article 478 Accesses Metrics details Biotic homogenization appears to be a global consequence of anthropogenic change. However, the

Scientific Reports volume 13, Article number: 6292 (2023) Cite this article

478 Accesses

Metrics details

Biotic homogenization appears to be a global consequence of anthropogenic change. However, the underlying environmental factors contributing to homogenization are difficult to identify because their effects usually interact and confound each other. This can be the reason why there is very little evidence on the role of climate warming in homogenization. By analysing macroinvertebrate assemblages in 65 streams that were as close to natural conditions as possible, we avoided the confounding effects of common anthropogenic stressors. This approach resulted in revealing a significant effect of increased temperature (both summer and winter) on changes in macroinvertebrate compositional over the past two decades. However, homogenization was significant only at opposite ends of the river continuum (submontane brooks, low-altitude rivers). Surprisingly, species of native origin predominated overall, increasing in frequency and abundance (“winners”), while only a minority of species declined or disappeared (“losers”). We hypothesise that undisturbed conditions mitigate species declines and thus homogenization, and that the temperature increase has so far been beneficial to most native species. Although we may have only captured a transitional state due to extinction debt, this underscores the importance of maintaining ecological conditions in stream to prevent species loss due to climate change.

Although biotic homogenizations very probably occurred locally during the whole evolutionary history1, anthropogenic activities in the industrial era seem to make it a global scale phenomenon2. In such temporal changes of biotic communities, common taxa spreading to new sites (“winners”) often replace spatially restricted taxa (“losers”), which increases a taxonomic and/or functional similarity among communities over time3. Although the exact mechanism behind homogenization is not universal4, environmental changes play a key role in the whole process. It may either facilitate establishment of a “winner” in a community5,6, or directly filter “winners” and “losers” from species pools. The “winners” are typically non-native species3,7, but they can be also native species8,9. In freshwaters, various stressors have been recognized to cause biotic homogenization, i.e., damming7,10, increased productivity11, salinization12, and changes in land use13. The communities with more similar compositions and functions potentially display decreased stability and resilience to further environmental changes14,15.

Temperature is one of the key environmental variables, which affects all levels of organization in biota from cells to communities16,17. The global increase in temperature18, impacting also running water systems19, can cause a decline in cold-adapted taxa and an increase in taxa adapted to higher temperatures20,21. Thus, it potentially leads to major compositional changes and functional restructuring of local communities22,23. However, whether and to what extent such climatically induced changes contribute to biotic homogenization has been only poorly studied24,25. That is mainly because numerous anthropogenic changes, such as land use, pollution and others, may largely confound the effect of increasing temperature26,27,28.

In this study, we analyse data on benthic macroinvertebrate assemblages repeatedly collected at 65 streams along the river continuum in the three time periods (1997–2000, 2007–2008, and 2015). Two main characteristics of this dataset make it especially suitable for exploring biotic homogenization in relation to the increase of temperature. First, we included only sites in an as-close-to natural state as possible to minimize the effect of common anthropogenic stressors. Second, our network of sampling sites evenly covers a large part of the river continuum from submontane brooks to lowland rivers. In contrast, most previous reports on long-term compositional changes in running waters have focussed only on a restricted section of the river continuum23,29. As different sections of the river network vary in proportion of cold-stenotherms, responses of macroinvertebrate assemblages to long-term changes are likely to differ22,23 or can lead to different degree of biotic homogenization8 along the river continuum.

Our primary goal was to test for the existence of compositional shifts with time leading to taxonomic homogenization (H1). As the main environmental conditions, selecting “winners” and “losers”, change in the longitudinal profile of streams, we expect that the degree of homogenization would differ along the river continuum (H2). Second, based on the evidence for both species profiting and losing due to the environmental changes in streams during the last decades, we assume homogenization to be caused equally by a decline of “losers” and an increase of “winners” (H3). And most importantly, as by including only close-to-natural sites, we hypothesise that temperature should appear as the main driving factor of the compositional changes through time (H4).

We selected 65 streams of 2–7th stream Strahler order, distributed at altitudes of 150–800 m a.s.l. in all main river basins of the Czech Republic (Fig. 1). To avoid confounding effects of other anthropogenic stressors, the main criterion for selection was minimal anthropogenic influence in terms of physical conditions, including natural thermal and discharge regime, water quality and the characters of riparian zone and floodplain, all being as-close-to natural state as possible. Owing to their natural character all these sites were used for establishment of type-specific reference conditions according to requirement of the Water Framework Directive, 2000/60/EC30.

Locations of the sampling sites in the Czech Republic: low altitude rivers (LowRiv), mid-altitude rivers (MidRiv), submontane rivers (SubRiv), mid-altitude brooks (MidBro), and submontane brooks (SubBro). The map was created using Arcgis 8.3 program59.

The study sites were classified according to their size (stream order according to Strahler, distance from source), catchment area and altitude into five stream types: submontane brooks (SubBro), submontane rivers (SubRiv), mid-altitude brooks (MidBro), mid-altitude rivers (MidRiv), and low altitude rivers (LowRiv) (Table 1). Unfortunately, there are no reference stream sites representing (a) lowland brooks and large rivers (such as lowland stretches of the Morava, Elbe and Vltava Rivers), which suffered morphological degradation and poor water quality, and (b) high mountain brooks and rivers, where we lack information on the past impact of acidification30.

Sampling was done during spring (April–May) and autumn (September–October) in three periods: 1997–2000 (period 1), 2007–2008 (period 2) and 2015 (period 3), using the same standardized methodology31. Benthic macroinvertebrates were sampled using a multi-habitat sampling method, i.e., taking semi-quantitative 3-min kick samples gathered with a hand net of 500 μm mesh size32. The sampling was done proportionally in all mesohabitats (e.g. riffle, pool, macrophytes, woody debris) within a typical stream stretch of up to 100 m long. This sampling method showed a high precision when compared to other methods33. The samples were preserved in 4% formaldehyde and sorted under stereomicroscope in the laboratory.

Macroinvertebrates were identified mainly to species, except of Diptera, some Oligochaeta and Coleoptera that were identified to genera. Because of variations in some aspects of taxonomic resolution through time, some species were pooled, using operational taxonomic units. This concerned some species groups of Ephemeroptera, Plecoptera, Trichoptera and Oligochaeta, with unclear reliability of identification in the 1990s. Chironomidae were excluded from the dataset due to varying taxonomic resolution between the periods.

Water quality parameters were measured six times per year during the respective years of macroinvertebrate sampling. They included water temperature, dissolved oxygen and oxygen saturation in water, biological oxygen demand (BOD5), pH, acid neutralising capacity (ANC), electrical conductivity and nutrient concentrations (TP, PO4-P, NO3-N, NH4-N) as required for ecological status assessment according to the Water Framework Directive, 2000/60/EC34. The mean roughness of bed substrate (Phi) was estimated by the cover of individual particle sizes35. Share of riffles (Riff) was assessed based on visual estimations in the sampling stretches during the macroinvertebrate sampling.

Three climatic variables, the average air temperatures in January and July (T_Jan, T_Jul) and annual precipitation (Prec), were calculated based on values of daily average temperature and daily precipitation in the year of sampling and two consecutive years prior to sampling, using gridded data provided by the Czech Hydrometeorological Institute36.

Four variables describing land use were chosen: areas of unfavourable surfaces (Un_surf), artificial surface (Ar_surf), arable land (Ar_land) and forested area (Fore). They were computed based on the information from the CORINE Land Cover system (for details see Supplements to Methods).

Prior to statistical analyses, the biotic data from spring and autumn within a sampling period were combined for each sampling site. To remove the influence of very rare species in all analyses, we included only the species recorded at least in three samples.

To visualize the compositional change in the macroinvertebrate assemblages with time, we performed Non-Metric Multidimensional Scaling (NMDS) of all samples, using Bray–Curtis dissimilarities of ln(x + 1) transformed abundance data, and drew trajectories of the assemblage developments in the ordination space as arrows.

As a measure of beta-diversity, we used the mean distance of sites to stream type centroids, calculated as Bray–Curtis dissimilarities of ln(x + 1) transformed species abundance data. We tested the homogeneity of multivariate dispersions among the three sampling periods using one-way ANOVA, both for the whole set of sites and also separately for each stream type.

We identified the species most responsible for the compositional changes between the time periods (1–2, 2–3, and 1–3), by categorizing the relative changes in species abundance and frequency. The relative change in species abundance between two periods was calculated using ln(x + 1) transformed abundance data (for details see Supplements to Methods) and it ranged from 1 (for the species that newly occurred at sites) to -1 (for the species that completely disappeared from all sites). For example, species with the relative change in abundance of 0.25 had 25% lower mean abundance in the former than in the latter period. We categorized species according to their relative change in species abundance between the periods 1–3, considering those within the intervals of [(-1)–(-0.5)] as “losers”, and those within the intervals of [0.5–1] as “winners”. Please note that our definition of “winners” and “losers” is arbitrary and can be used only for mutual comparisons in this study. To assess the differences in species composition along the river continuum we analysed the relative change in species abundance for each of the stream types separately.

To test the effects of environmental variables, we performed distance-based Redundancy Analysis (db-RDA), using Bray–Curtis dissimilarities of ln(x + 1) transformed abundance data and forward selection procedure for environmental variables. The significance of each potentially added variable was tested by 999 permutations restricted within sites, as we were interested in the temporal within-site variability. Variables with the highest contribution to the explained variation in the species data were added sequentially in the model if the p-value of their contribution was < 0.05. We also partitioned the pure and shared effects of temperature (T_Jan, T_Jul) and other significant environmental variables in the models. Db-RDAs were done both for the whole set of sites and separately for each stream type.

All analyses and visualizations were performed in R version 3.4.437, using packages vegan 2.5–638 and ggplot239 packages.

The main gradient in taxonomic composition was associated with the stream typology based on stream size and altitude (Fig. 2), reflecting the river continuum. Individual sites showed considerable shifts in species composition over the time, heading closer to each other in the NMDS plot (Fig. 2), which might indicate an overall increase in assemblage similarity among sites. Despites this, total beta-diversity did not significantly differ between the three periods (p = 0.1742), suggesting only a trend to taxonomic homogenization among all the study sites. However, when beta-diversity was analysed for stream types separately, it significantly differed in submontane brooks (p = 0.0042) and lowland rivers (p = 0.0139) (Fig. 3).

Compositional shift of macroinvertebrate assemblages in the 65 stream sites repeatedly sampled in 1997–2000 (period 1), 2007–2008 (period 2) and 2015 (period 3). Species data were analysed by Non-Metric Multidimensional Scaling (NMDS) on Bray–Curtis dissimilarities of ln(x + 1) transformed abundances; stress = 0.153. The sites are classified into five types along a river continuum. Samples from the sampling periods are connected by an arrow for each site, heading towards the youngest time period.

Beta-diversity of macroinvertebrate assemblages in the 65 repeatedly sampled stream sites in 1997–2000 (period 1), 2007–2008 (period 2) and 2015 (period 3). Streams were divided into five types along a river continuum to be analysed separately. Beta-diversity of each group was counted as mean distance of sites to the group centroid created by principle coordinate analysis (PCoA), using Bray–Curtis dissimilarities of ln(x + 1) transformed abundance data. Homogeneity of multivariate dispersions in the three periods was tested using one-way ANOVA. n.s.—p > 0.05.

Comparing the species composition between the first (1997–2000) and the third (2015) period, there was a relatively large number of species that newly occurred or increased their mean abundance at sites by up to 50%, and we considered them as “winners” (Table 2). In total, there were 158 macroinvertebrate species classified as “winners”, most of them (ca 79%), however, only in one or two stream types (Supplementary Table S1). Only seven taxa showed as “winners” in all stream types. There were only three non-native species among “winners”: planarian Girardia tigrina, gastropod Potamopyrgus antipodarum, and isopod Proasellus coxalis, all of them increased only in a single stream type. On the other hand, only several species (from one to five species) disappeared from the sites in each stream type, all of them being rare or relatively rare (frequency of 3–5 in the first period). Abundance of some other species (from four to nine species), decreased by more than 50% in a stream type. These two groups of diminishing species were considered “losers”, summing up to 37 spp. in total (Supplementary Table S2). The species frequencies seemed to change less than species abundances among the sampling periods but they showed similar patterns (Supplementary Fig. S1).

Environmental variables explained 7.1% of the total variance for the whole dataset and 10.7–18.3% of the total variance in individual stream types (Fig. 4). Temperature variables (either July or January air temperatures, or both) were significant in the whole dataset (pure effect of 5.0%) and all stream types (pure effects of 5.5–11.6%). Pure effects of other environmental variables were lower than those of temperature: 2.7% in the whole dataset and 2.3–5.2% in individual stream types. Roughness of bed substrate and annual precipitation sums were the most significant of the other environmental variables. The effects shared between temperature and other significant variables ranged from 0 to 6.2% in individual stream types. The changes in the significant environmental variables at individual sites among the three sampling periods are given in Supplementary Fig. S2.

Total, pure and shared variation explained by temperature and other environmental variables in macroinvertebrate assemblages of the 65 stream sites compared between 1997–2000 (period 1) and 2015 (period 3). Counted as R2adj. by Distance-Based Redundancy Analysis (db-RDA) using Bray–Curtis dissimilarities of ln(x + 1) transformed abundance data, for the total dataset and separately for five stream types: low altitude rivers (LowRiv), mid-altitude rivers (MidRiv), submontane rivers (SubRiv), mid-altitude brooks (MidBro), and submontane brooks (SubBro). Significant environmental variables: mean air temperatures in January and July (T_Jan, T_Jul), annual precipitation (Prec), roughness of bed substrate (Phi), concentration of nitrates (NO3), artificial surfaces (Ar_surf), unfavourable surfaces (Un_surf), share of riffles (Riff), forested area (Fore), biochemical oxygen demand (BOD5). The models were significant at p < 0.01, permutations were restricted to sites. Temp—pure effects of temperature variables; Shared—effects shared by temperature and other environmental variables; Other—pure effects of other environmental variables.

We found that composition of macroinvertebrate assemblages in the Czech reference stream sites has undergone a substantial change over the last two decades. This agrees with other studies reporting compositional changes in stream communities from different parts of the world23,25,40. Despite the change, we found no significant taxonomic homogenization on the whole river continuum scale, disproving our H1. However, beta-diversity among the studied periods showed a decreasing tendency with time, implying that the process that may in future eventually lead to homogenization on this regional scale has probably already started. We suggest that the close-to-natural conditions at our sites and naturally high between-site variation decelerate homogenization. Since multiple stressors in running water systems have typically synergic effects16,25, we might expect more significant homogenization to be associated with non-reference state. For example, macroinvertebrates in streams from a wide range of altitudes, but affected by common anthropogenic changes, experienced taxonomic homogenization25.

When the stream types were analysed separately, homogenization was significant at both ends of the river continuum (submonate brooks, lowland rivers) but not in the middle. This confirms our H2 that the degree of homogenization differs along the river continuum. It also suggests different speed of homogenization along the river continuum, which seemed relatively slower in the middle of the river continuum. This might correspond to the rare case when no signs of taxonomic homogenization have been found in mid-altitude streams in Wales, UK, over the period of 30 years41. In contrast, both submontane brooks and lowland rivers homogenized relatively fast (within a single decade given the grain of our study), although in the former it happened a decade earlier than in the latter. We could thus confirm the recognition that biotic homogenization is not a linear process42.

Most of the species clearly increased their abundance and frequency at our sites over the studied period, while the loss of species played only a negligible role. This partly contradicted our expectation (H3) because although species loss/decrease and species gain/decrease are not necessarily reciprocal, homogenization of natural communities usually involves both these processes25,42.

In direct association with climate change, a decline of cold-adapted species has been observed in freshwater assemblages23,43. We could not confirm a significant decline in cold-stenothermal species, although the proportion of cold-stenotherms was higher among the “losers” (ca 46%) than among the whole species pool (22.8%). The decline of cold-stenotherms (mostly plecopterans and trichopterans) accounted only for a very minor compositional change even in submontane brooks and rivers, where stenothermic species naturally dominate. However, as the mountain streams are most vulnerable to climate change, we might expect some more cold-stenothermic “losers” in streams at higher elevations (> 900 m a.s.l.23,44).

“Winners” constituted the dominant proportion of species in each stream type, suggesting an immigration credit42 and increasingly favourable conditions for many macroinvertebrate species. This was quite surprising because, although “winners” have been found to prevail in some freshwater assemblages45,46, in general, homogenization due to prevailing species gain has been only rarely reported8,9. Most importantly, the vast majority of the “winners” in our study (mostly trichopterans, dipterans, clitellates, but many others) were from regional species pools, while only a very minor fraction (< 2%) was of non-native origin, their increase being limited for only one stream type. That does not accord with the predictions of establishment of non-native macroinvertebrates in stream assemblages due to climate change44 and the typical cases of homogenization due to an increase of alien species, which induce extirpations of native species4,7.

Several studies have predicted or already observed an increase of warm-stenotherms and eurytherms in stream macroinvertebrate assemblages20,23,46. However, the increase of species in our study could not be explained by known species thermal tolerances47. It is a reminder that thermal preferences, being mostly derived from the known species distributions, are probably only a very rough estimates of species temperature niches. However, analysis of a possible role of other species traits (such as substrate preference, flow velocity, and sensitivity to trophic status47) for the success of “winners” is far outside the scope of this paper.

An increase of native species in assemblages may be a sign of recovery from disturbances48, as found in some European streams following the improvement of water quality since the end of 1980s49,50. At our sites with persistently close-to-natural conditions, an increase of some species might partly reflect a recovery of ecological state in surrounding streams. The increased number of suitable sites in the landscape might support the populations of native species by higher rates of dispersal.

In our study, temperature appeared to be the most significant environmental factor explaining the between period compositional changes at sites. The significant temperature effects were found across the whole river continuum and in all individual stream types, confirming our H4. With only few exceptions45,51, the previous studies could not unequivocally assign long-term compositional changes in macroinvertebrate assemblages to climate change due to the confounding effects of changes in water quality, land use, or strong hydrological events25,49,50. As in our study, temperature was the environmental variable with the largest and most consistent temporal change, we could support the significant effect of temperature on the macroinvertebrate assemblages by a strong evidence.

Temperature affects organisms directly, modulating their metabolism and phenology16,17, and/or indirectly by changing the productivity of the whole system45. By either way, the increase in temperature allowed some previously temperature limited species to establish more abundant populations. The increase in winter temperatures, suggesting a relaxed limitation by cold extremes52, was most significant in submontane brooks and rivers, but performed across whole river continuum. The increase in warm summer temperatures was significant only in the middle and lower parts of the river continuum, having surprisingly positive effect on many species.

A long term systematic enrichment in stream macroinvertebrates has been explained by a higher primary production, leading to a higher complexity in food webs45. There are no data of primary production available for our sites, however, all the indicators of tropic status (total phosphorus, nitrates, BOD5) did not show any consistent trend that could account for the large increase in macroinvertebrates (Supplementary Table S3).

The impact of climate change on running waters include the combined effects of temperature and flow alternations53,54. Unfortunately, direct measurements of discharge were not available for our sites. However, as precipitation (Prec) decreased with time at most of our sites, we assume a decreasing trend in discharge over the studied period. Structure of sediment (Phi) was another variable indicative of change in discharge. As both Prec and Phi were the most significant explanatory variables after temperature, we suggest that the significant change in macroinvertebrate assemblages is partly also accountable to the change in discharge. In lowland rivers, the reduction in discharge favoured macroinvertebrates associated with finer sediments and those preferring slow flowing conditions, similarly as in some other European rivers50.

In small streams of mid-altitudes, the changes in discharge might have led to an occasional temporal drying during the last studied time period. Due to the humid continental climate, Czech streams are historically mostly perennial, however, there has been an increase in flow intermittence during the last few decades, most recently at an accelerating rate55,56. Drying is a strong environmental filter for aquatic fauna favouring only resistant and resilient taxa57 and promoting nested communities with higher dominance of ubiquitous taxa58, which might result in homogenization. However, there is an evidence that newly occurring intermittent conditions support a higher replacement of taxa in streams and thus increasing beta-diversity56. Occasional drying might hamper the increase of abundance and frequency in “winners”, decelerating the homogenization in mid-altitudes.

To sum up, the increase of air temperature in our study led to a large compositional change and/or, depending on the position in the river continuum, also taxonomic homogenization, all of that mainly due to a large increase in abundance and gain of species. Macroinvertebrate assemblages have been clearly undergoing massive reorganization during the last 20 years but not yet decrease and loss of species. However, we might have, in fact, observed only a temporal transitional state of the increase due to an extinction delay42. Such transitional state may be followed by a decline in abundances and loss of species with further increase of temperature after reaching a threshold51. Massive decline is probable especially if the increase of temperature is combined with strong changes in discharge conditions53,54. Nevertheless, the future development in stream assemblages remains highly unpredictable. In accordance with Van Looy et al.45, our study suggests that loss of species may be largely prevented even despite the ongoing climate change by maintaining or improving the ecological quality of streams.

Data are available from the authors upon reasonable request—contact Jindřiška Bojková ([email protected]).

Vermeij, G. J. When biotas meet: Understanding biotic interchange. Science 253, 1099–1104. https://doi.org/10.1126/science.253.5024.1099 (1991).

Article ADS CAS PubMed Google Scholar

Baiser, B., Olden, J. D., Record, S., Lockwood, J. L. & McKinney, M. L. Pattern and process of biotic homogenization in the New Pangaea. Proc. R. Soc. B Biol. Sci. 279, 4772–4777. https://doi.org/10.1098/rspb.2012.1651 (2012).

Article Google Scholar

McKinney, M. L. & Lockwood, J. L. Biotic homogenization: A few winners replacing many losers in the next mass extinction. Trends Ecol. Evol. 14, 450–453. https://doi.org/10.1016/S0169-5347(99)01679-1 (1999).

Article CAS PubMed Google Scholar

Olden, J. D. & Rooney, T. P. On defining and quantifying biotic homogenization. Glob. Ecol. Biogeogr. 15, 113–120. https://doi.org/10.1111/j.1466-822X.2006.00214.x (2006).

Article Google Scholar

Rahel, F. J. & Olden, J. D. Assessing the effects of climate change on aquatic invasive species. Conserv. Biol. 22, 521–533. https://doi.org/10.1111/j.1523-1739.2008.00950.x (2008).

Article PubMed Google Scholar

Radinger, J., Alcaraz-Hernández, J. D. & García-Berthou, E. Environmental filtering governs the spatial distribution of alien fishes in a large, human - impacted Mediterranean river. Divers. Distrib. 25(5), 701–714. https://doi.org/10.1111/ddi.12895 (2019).

Article Google Scholar

Rahel, F. J. Homogenization of freshwater faunas. Annu. Rev. Ecol. Syst. 33, 291–315 (2002).

Article Google Scholar

Buisson, L. & Grenouillet, G. Contrasted impacts of climate change on stream fish assemblages along an environmental gradient. Divers. Distrib. 15, 613–626. https://doi.org/10.1111/j.1472-4642.2009.00565.x (2009).

Article Google Scholar

Finderup Nielsen, T., Sand-Jensen, K., Dornelas, M. & Bruun, H. H. More is less: Net gain in species richness, but biotic homogenization over 140 years. Ecol. Lett. 220, 1650–1657. https://doi.org/10.1111/ele.13361 (2019).

Article Google Scholar

Pool, T. K. & Olden, J. D. Taxonomic and functional homogenization of an endemic desert fish fauna. Divers. Distrib. 18, 366–376. https://doi.org/10.1111/j.1472-4642.2011.00836.x (2012).

Article Google Scholar

Donohue, I., Jackson, A. L., Pusch, M. T. & Irvine, K. Nutrient enrichment homogenizes lake benthic assemblages at local and regional scales. Ecology 90(12), 3470–3477. https://doi.org/10.1890/09-0415.1 (2009).

Article PubMed Google Scholar

Miyazono, S., Patiño, R. & Taylor, C. M. Desertification, salinization, and biotic homogenization in a dryland river ecosystem. Sci. Total Environ. 511, 444–453. https://doi.org/10.1016/j.scitotenv.2014.12.079 (2015).

Article ADS CAS PubMed Google Scholar

Siqueira, T., Lacerda, C. G. T. & Saito, V. S. How does landscape modification induce biological homogenization in tropical stream metacommunities?. Biotropica 47, 509–516. https://doi.org/10.1111/btp.12224 (2015).

Article Google Scholar

Olden, J. D. Biotic homogenization: A new research agenda for conservation biogeography. J. Biogeogr. 33, 2027–2039. https://doi.org/10.1111/j.1365-2699.2006.01572.xKjjksdnp15 (2006).

Article Google Scholar

Clavel, J., Julliard, R. & Devictor, V. Worldwide decline of specialist species: Toward a global functional homogenization?. Front. Ecol. Environ. 9, 222–228. https://doi.org/10.1890/080216 (2011).

Article Google Scholar

Woodward, G., Perkins, D. M. & Brown, L. E. Climate change and freshwater ecosystems: Impacts across multiple levels of organization. Philos. Trans. R. Soc. B Biol. Sci. 365, 2093–2106. https://doi.org/10.1098/rstb.2010.0055 (2010).

Article Google Scholar

Seebacher, F., White, C. R. & Franklin, C. E. Physiological plasticity increases resilience of ectothermic animals to climate change. Nat. Clim. Change Nat. 5, 61–66. https://doi.org/10.1038/nclimate2457 (2015).

Article ADS Google Scholar

Sala, O. E. et al. Global biodiversity scenarios for the year 2100. Science 287, 1770–1774. https://doi.org/10.1126/science.287.5459.1770 (2000).

Article CAS PubMed Google Scholar

Arora, R., Tockner, K. & Venohr, M. Changing river temperatures in northern Germany: Trends and drivers of change. Hydrol. Process. 30, 3084–3096. https://doi.org/10.1002/hyp.10849 (2016).

Article ADS Google Scholar

Domisch, S. et al. Modelling distribution in European stream macroinvertebrates under future climates. Glob. Change Biol. 19, 752–762. https://doi.org/10.1111/gcb.12107 (2013).

Article ADS Google Scholar

Nukazawa, K., Arai, R., Kazama, S. & Takemon, Y. Projection of invertebrate populations in the headwater streams of a temperate catchment under a changing climate. Sci. Total Environ. 642, 610–618. https://doi.org/10.1016/j.scitotenv.2018.06.109 (2018).

Article ADS CAS PubMed Google Scholar

Buisson, L., Thuiller, W., Lek, S., Lim, P. & Grenouillet, G. Climate change hastens the turnover of stream fish assemblages. Glob. Change Biol. 14, 2232–2248. https://doi.org/10.1111/j.1365-2486.2008.01657.x (2008).

Article ADS Google Scholar

Haase, P. et al. Moderate warming over the past 25 years has already reorganized stream invertebrate communities. Sci. Total Environ. 658, 1531–1538. https://doi.org/10.1016/j.scitotenv.2018.12.234 (2019).

Article ADS CAS PubMed Google Scholar

Magurran, A. E., Dornelas, M., Moyes, F., Gotelli, N. J. & McGill, B. Rapid biotic homogenization of marine fish assemblages. Nat. Commun. 6, 8405. https://doi.org/10.1038/ncomms9405 (2015).

Article ADS CAS PubMed Google Scholar

Mouton, T. L., Tonkin, J. D., Stephenson, F., Verburg, P. & Floury, M. Increasing climate-driven taxonomic homogenization but functional differentiation among river macroinvertebrate assemblages. Glob. Change Biol. 26, 6904–6915. https://doi.org/10.1111/gcb.15389 (2020).

Article ADS Google Scholar

Durance, I. & Ormerod, S. J. Trends in water quality and discharge confound long-term warming effects on river macroinvertebrates. Freshw. Biol. 54, 388–405. https://doi.org/10.1111/j.1365-2427.2008.02112.x (2009).

Article CAS Google Scholar

Oliver, T. H. & Morecroft, M. D. Interactions between climate change and land use change on biodiversity: Attribution problems, risks, and opportunities. Wiley Interdiscip. Rev. Clim. Change 5, 317–335. https://doi.org/10.1002/wcc.271 (2014).

Article Google Scholar

Kuemmerlen, M. et al. An attack on two fronts: Predicting how changes in land use and climate affect the distribution of stream macroinvertebrates. Freshw. Biol. 60, 1443–1458. https://doi.org/10.1111/fwb.12580 (2015).

Article Google Scholar

Zedková, B., Rádková, V., Bojková, J., Soldán, T. & Zahrádková, S. Mayflies Ephemeroptera) as indicators of environmental changes in the past five decades: A case study from the Morava and Odra River Basins (Czech Republic). Aquat. Conserv. Mar. Freshwat. Ecosyst. 25, 622–638. https://doi.org/10.1002/aqc.2529 (2015).

Article Google Scholar

Kokeš, J. et al. The PERLA system in the Czech Republic: A multivariate approach to assess ecological status of running waters. Hydrobiologia 566, 343–354. https://doi.org/10.1007/s10750-006-0085-4 (2006).

Article Google Scholar

Czech State Norm No. 75 7701. Water quality - Methodology for sampling and treatment of macroinvertebrates from running waters using method PERLA. Úřad pro technickou normalizaci, metrologii a státní zkušebnictví, Praha (2008).

Birk, S., Strackbein, J. & Hering, D. WISER methods database. Version: March 2011. http://www.wiser.eu/results/method-database/ (2010).

Friberg, N. et al. Comparison of macroinvertebrate sampling methods in Europe. Hydrobiologia 566, 365–378. https://doi.org/10.1007/s10750-006-0083-6 (2006).

Article Google Scholar

Rosendorf, P. Metodika hodnocení všeobecných fyzikálně chemických složek ekologického stavu útvarů povrchových vod tekoucích. Výzkumný ústav vodohospodářský (T.G. Masaryka, Praha, 2011).

Furse M. T., Moss D., Wright J. F., Armitage P. D. & Gunn R. J. M. A Practical Manual to the Classification and Prediction of Macroinvertebrate Communities in Running Water in Great Britain (Freshwater Biological Association, River laboratory, 1986).

Polášek, M., Zahrádková, S., Němejcová, D., Straka, M., Bareš, M., & Opatřilová, L. Monitoring of Long-Term Changes in the Biodiversity of Running Waters at the Time of Climate Change: Proposal, Implementation and Incorporation into the ARROW Public Information System (EHP-CZ02-OV-1-018-2014) (T. G. Masaryk Water Research Institute, Brno, 2017).

R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2018).

Google Scholar

Oksanen, J., et al. vegan: Community Ecology Package. R package version 2.5–6. https://CRAN.R-project.org/package=vegan (2019).

Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).

Floury, M., Souchon, Y. & van Looy, K. Climatic and trophic processes drive long-term changes in functional diversity of freshwater invertebrate communities. Ecography 41(1), 209–218. https://doi.org/10.1111/ecog.02701 (2018).

Article Google Scholar

Larsen, S., Chase, J. M., Durance, I. & Ormerod, S. J. Lifting the veil: Richness measurements fail to detect systematic biodiversity change over three decades. Ecology 99, 1316–1326. https://doi.org/10.1002/ecy.2213 (2018).

Article PubMed Google Scholar

Hillebrand, H. et al. Biodiversity change is uncoupled from species richness trends: Consequences for conservation and monitoring. J. Appl. Ecol. 55, 169–184. https://doi.org/10.1111/1365-2664.12959 (2018).

Article Google Scholar

Jyväsjärvi, J. et al. Climate-induced warming imposes a threat to north European spring ecosystems. Glob. Change Biol. 21, 4561–4569. https://doi.org/10.1111/gcb.13067 (2015).

Article ADS Google Scholar

Domisch, S., Jähnig, S. C. & Haase, P. Climate-change winners and losers: Stream macroinvertebrates of a submontane region in Central Europe. Freshw. Biol. 56, 2009–2020. https://doi.org/10.1111/j.1365-2427.2011.02631.x (2011).

Article Google Scholar

Van Looy, K., Floury, M., Ferréol, M., Prieto-Montes, M. & Souchon, Y. Long-term changes in temperate stream invertebrate communities reveal a synchronous trophic amplification at the turn of the millennium. Sci. Total Environ. 565, 481–488. https://doi.org/10.1016/j.scitotenv.2016.04.193 (2016).

Article ADS CAS PubMed Google Scholar

Bruno, D. et al. Structural and functional responses of invertebrate communities to climate change and flow regulation in alpine catchments. Glob. Change Biol. 25, 1612–1628. https://doi.org/10.1111/gcb.14581 (2019).

Article ADS Google Scholar

Schmidt-Kloiber, A. & Hering, D. www.freshwaterecology.info—an online tool that unifies, standardises and codifies more than 20,000 European freshwater organisms and their ecological preferences. Ecol. Ind. 53, 271–282. https://doi.org/10.1016/j.ecolind.2015.02.007 (2015).

Article Google Scholar

Cardinale, B. J., Gonzalez, A., Allington, G. R. H. & Loreau, M. Is local biodiversity declining or not? A summary of the debate over analysis of species richness time trends. Biol. Cons. 219, 175–183. https://doi.org/10.1016/j.biocon.2017.12.021 (2018).

Article Google Scholar

Vaughan, I. P. & Ormerod, S. J. Large-scale, long-term trends in British river macroinvertebrates. Glob. Change Biol. 18, 2184–2194. https://doi.org/10.1111/j.1365-2486.2012.02662.x (2012).

Article ADS Google Scholar

Floury, M., Usseglio-Polatera, P., Ferreol, M., Delattre, C. & Souchon, Y. Global climate change in large European rivers: Long-term effects on macroinvertebrate communities and potential local confounding factors. Glob. Change Biol. 19, 1085–1099. https://doi.org/10.1111/gcb.12124 (2013).

Article ADS Google Scholar

Baranov, V., Jourdan, J., Pilotto, F., Wagner, R. & Haase, P. Complex and nonlinear climate-driven changes in freshwater insect communities over 42 years. Conserv. Biol. 34(5), 1241–1251. https://doi.org/10.1111/cobi.13477 (2020).

Article PubMed Google Scholar

Durance, I. & Ormerod, S. J. Climate change effects on upland stream macroinvertebrates over a 25-year period. Glob. Change Biol. 13, 942–957. https://doi.org/10.1111/j.1365-2486.2007.01340.x (2007).

Article ADS Google Scholar

Pyne, M. I. & Poff, N. L. Vulnerability of stream community composition and function to projected thermal warming and hydrologic change across ecoregions in the western United States. Glob. Change Biol. 23, 77–93. https://doi.org/10.1111/gcb.13437 (2017).

Article ADS Google Scholar

Kakouei, K. et al. Projected effects of climate change induced flow alterations on stream macroinvertebrate abundances. Ecol. Evol. 8, 3393–3409. https://doi.org/10.1002/ece3.3907 (2018).

Article PubMed PubMed Central Google Scholar

Pařil, P. et al. An unexpected source of invertebrate community recovery in intermittent streams from a humid continental climate. Freshw. Biol. 64, 1971–1983. https://doi.org/10.1111/fwb.13386 (2019).

Article Google Scholar

Crabot, J., Polášek, M., Launay, B., Pařil, P. & Datry, T. Drying in newly intermittent rivers leads to higher variability of invertebrate communities. Freshw. Biol. 66, 730–744. https://doi.org/10.1111/fwb.13673 (2020).

Article Google Scholar

Leigh, C. & Datry, T. Drying as a primary hydrological determinant of biodiversity in river systems: A broad scale analysis. Ecography 40(4), 487–499. https://doi.org/10.1111/ecog.02230 (2017).

Article Google Scholar

Datry, T. et al. Broad-scale patterns of invertebrate richness and community composition in temporary rivers: Effects of flow intermittence. Ecography 37, 94–104. https://doi.org/10.1111/j.1600-0587.2013.00287.x (2014).

Article Google Scholar

ESRI (2003) ArcGIS 8.3. Environmental Systems Research Institute, Redlands, CA, USA. http://www.esri.com. Accessed 2 Mar 2023.

Download references

We are grateful to Jiří Kokeš, Vanda Šorfová, and Michal Straka, for their help with field work and species identifications, to Pavel Sedláček and Jaroslav Sova for their help with field work, to Petr Pařil, Jan Špaček, Petr Komzák, Stanislav Větříček, Jana Holasová, and Vendula Polášková for species identifications, and to Veronika Horsáková for creating the map of sites. Our special thanks go to Světlana Zahrádková for her relentless efforts and leadership when supervising the whole PERLA and RIVERCHANGE Projects. The study was supported by the Czech Science Foundation (P505/20-17305S).

These authors contributed equally: Marie Zhai and Jindřiška Bojková.

Department of Botany and Zoology, Masaryk University, Kotlářská 2, 611 37, Brno, Czech Republic

Marie Zhai, Jindřiška Bojková, Marek Polášek, Vít Syrovátka & Michal Horsák

T. G. Masaryk Water Research Institute, p.r.i., Podbabská 2582/30, 160 00, Prague 6, Czech Republic

Denisa Němejcová & Marek Polášek

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.B., M.H. and M.Z. conceived the ideas and study design; J.B., D.N. and M.P. developed methodology and collected the data; V.S. analysed the data; M.Z. and J.B. led the writing as joint first author. All authors contributed critically to the drafts and gave final approval for publication.

Correspondence to Marie Zhai.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Zhai, M., Bojková, J., Němejcová, D. et al. Climatically promoted taxonomic homogenization of macroinvertebrates in unaffected streams varies along the river continuum. Sci Rep 13, 6292 (2023). https://doi.org/10.1038/s41598-023-32806-y

Download citation

Received: 14 December 2022

Accepted: 03 April 2023

Published: 18 April 2023

DOI: https://doi.org/10.1038/s41598-023-32806-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.