Inland navigation and land use interact to impact European freshwater biodiversity.

Inland navigation in Europe is proposed to increase in the coming years, being promoted as a low-carbon form of transport. However, we currently lack knowledge on how this would impact biodiversity at large scales and interact with existing stressors. Here we addressed this knowledge gap by analysing fish and macroinvertebrate community time series across large European rivers comprising 19,592 observations from 4,049 sampling sites spanning the past 32 years. We found ship traffic to be associated with biodiversity declines, that is, loss of fish and macroinvertebrate taxonomic richness, diversity and trait richness. Ship traffic was also associated with increases in taxonomic evenness, which, in concert with richness decreases, was attributed to losses in rare taxa. Ship traffic was especially harmful for benthic taxa and those preferring slow flows. These effects often depended on local land use and riparian degradation. In fish, negative impacts of shipping were highest in urban and agricultural landscapes. Regarding navigation infrastructure, the negative impact of channelization on macroinvertebrates was evident only when riparian degradation was also high. Our results demonstrate the risk of increasing inland navigation on freshwater biodiversity. Integrative waterway management accounting for riparian habitats and landscape characteristics could help to mitigate these impacts.

The ecological causes of functional distinctiveness in communities.

Recent work has shown that evaluating functional trait distinctiveness, the average trait distance of a species to other species in a community offers promising insights into biodiversity dynamics and ecosystem functioning. However, the ecological mechanisms underlying the emergence and persistence of functionally distinct species are poorly understood. Here, we address the issue by considering a heterogeneous fitness landscape whereby functional dimensions encompass peaks representing trait combinations yielding positive population growth rates in a community. We identify four ecological cases contributing to the emergence and persistence of functionally distinct species. First, environmental heterogeneity or alternative phenotypic designs can drive positive population growth of functionally distinct species. Second, sink populations with negative population growth can deviate from local fitness peaks and be functionally distinct. Third, species found at the margin of the fitness landscape can persist but be functionally distinct. Fourth, biotic interactions (positive or negative) can dynamically alter the fitness landscape. We offer examples of these four cases and guidelines to distinguish between them. In addition to these deterministic processes, we explore how stochastic dispersal limitation can yield functional distinctiveness. Our framework offers a novel perspective on the relationship between fitness landscape heterogeneity and the functional composition of ecological assemblages.

A functional vulnerability framework for biodiversity conservation.

Setting appropriate conservation strategies in a multi-threat world is a challenging goal, especially because of natural complexity and budget limitations that prevent effective management of all ecosystems. Safeguarding the most threatened ecosystems requires accurate and integrative quantification of their vulnerability and their functioning, particularly the potential loss of species trait diversity which imperils their functioning. However, the magnitude of threats and associated biological responses both have high uncertainties. Additionally, a major difficulty is the recurrent lack of reference conditions for a fair and operational measurement of vulnerability. Here, we present a functional vulnerability framework that incorporates uncertainty and reference conditions into a generalizable tool. Through in silico simulations of disturbances, our framework allows us to quantify the vulnerability of communities to a wide range of threats. We demonstrate the relevance and operationality of our framework, and its global, scalable and quantitative comparability, through three case studies on marine fishes and mammals. We show that functional vulnerability has marked geographic and temporal patterns. We underline contrasting contributions of species richness and functional redundancy to the level of vulnerability among case studies, indicating that our integrative assessment can also identify the drivers of vulnerability in a world where uncertainty is omnipresent.

Direct habitat descriptors improve the understanding of the organization of fish and macroinvertebrate communities across a large catchment.

In large-scale aquatic ecological studies, direct habitat descriptors (e.g. water temperature, hydraulics in river reaches) are often approximated by coarse-grain surrogates (e.g. air temperature, discharge respectively) since they are easier to measure or model. However, as biological variability can be very strong at the habitat scale, surrogate variables may have a limited ability to capture all of this variability, which may lead to a lesser understanding of the ecological processes or patterns of interest. In this study, we aimed to compare the capacity of direct habitat descriptors vs. surrogate environmental variables to explain the organization of fish and macroinvertebrate communities across the Loire catchment in France (105 km2). For this purpose, we relied on high-resolution environmental data, extensive biological monitoring data (>1000 sampling stations) and multivariate analyses. Fish and macroinvertebrate abundance datasets were considered both separately and combined to assess the value of a cross-taxa approach. We found that fish and macroinvertebrate communities exhibited weak concordance in their organization and responded differently to the main ecological gradients. Such variations are probably due to fundamental differences in their life-history traits and mobility. Regardless of the biological group considered, direct habitat descriptors (water temperature and local hydraulic variables) consistently explained the organization of fish and macroinvertebrate communities better than surrogate descriptors (air temperature and river discharge). Furthermore, the organization of fish and macroinvertebrate communities was slightly better explained by the combination of direct or surrogate environmental variables when the two biological groups were considered together than when considered separately. Tied together, these results emphasize the importance of using a cross-taxa approach in association with high-resolution direct habitat variables to more accurately explain the organization of aquatic communities.

Invasion impacts and dynamics of a European-wide introduced species.

Globalization has led to the introduction of thousands of alien species worldwide. With growing impacts by invasive species, understanding the invasion process remains critical for predicting adverse effects and informing efficient management. Theoretically, invasion dynamics have been assumed to follow an "invasion curve" (S-shaped curve of available area invaded over time), but this dynamic has lacked empirical testing using large-scale data and neglects to consider invader abundances. We propose an "impact curve" describing the impacts generated by invasive species over time based on cumulative abundances. To test this curve's large-scale applicability, we used the data-rich New Zealand mud snail Potamopyrgus antipodarum, one of the most damaging freshwater invaders that has invaded almost all of Europe. Using long-term (1979-2020) abundance and environmental data collected across 306 European sites, we observed that P. antipodarum abundance generally increased through time, with slower population growth at higher latitudes and with lower runoff depth. Fifty-nine percent of these populations followed the impact curve, characterized by first occurrence, exponential growth, then long-term saturation. This behaviour is consistent with boom-bust dynamics, as saturation occurs due to a rapid decline in abundance over time. Across sites, we estimated that impact peaked approximately two decades after first detection, but the rate of progression along the invasion process was influenced by local abiotic conditions. The S-shaped impact curve may be common among many invasive species that undergo complex invasion dynamics. This provides a potentially unifying approach to advance understanding of large-scale invasion dynamics and could inform timely management actions to mitigate impacts on ecosystems and economies.

Does global change increase the risk of maladaptation of Atlantic salmon migration through joint modifications of river temperature and discharge?

In freshwater ecosystems, water temperature and discharge are two intrinsically associated triggers of key events in the life cycle of aquatic organisms such as the migration of diadromous fishes. However, global changes have already profoundly altered the thermal and hydrological regimes of rivers, affecting the timing of fish migration as well as the environmental conditions under which it occurs. In this study, we focused on Atlantic salmon (Salmo salar), an iconic diadromous species whose individuals migrate between marine nursery areas and continental spawning grounds. An innovative multivariate method was developed to analyse long-term datasets of daily water temperature, discharge and both salmon juvenile downstream and adult upstream migrations in three French rivers (the Bresle, Oir and Nivelle rivers). While all three rivers have gradually warmed over the last 35 years, changes in discharge have been very heterogeneous. Juveniles more frequently used warmer temperatures to migrate. Adults migrating a few weeks before spawning more frequently used warm temperatures associated with high discharges. This has already led to modifications in preferential niches of both life stages and suggests a potential mismatch between these populations' ecological preference and changes in their local environment due to global change.

Joint temporal trends in river thermal and hydrological conditions can threaten the downstream migration of the critically endangered European eel.

Climate change is modifying the hydrological and thermal regimes of rivers worldwide, threatening the triggering of organisms' key life-cycle processes. European eel (Anguilla anguilla) is a critically endangered fish species that migrates over several thousand kilometres between its rearing habitats in continental waters of Europe and North Africa and its spawning area in the Sargasso Sea. Downstream migration of adult eels occurs during periods of decreasing river water temperature associated with high discharge but changes in these environmental cues may affected eel migratory conditions. An innovative multivariate method was developed to analyse long-term datasets of daily water temperature, discharge and eel passage in two European rivers. Over the past 50 years, water temperature and discharge increased in both rivers during the downstream migration period from August to November. Silver eels preferentially migrated at temperatures between 10 and 20 °C combined with high discharge. Environmental changes have resulted in the migration of silver eels under warmer water temperatures. This example illustrates how the changes in environmental cues have led to a growing mismatch between the migratory conditions preferentially selected and those actually used, which may threaten the completion of the eel's life cycle and ultimately the persistence of this already critically endangered species.

The dimensionality and structure of species trait spaces.

Trait-based ecology aims to understand the processes that generate the overarching diversity of organismal traits and their influence on ecosystem functioning. Achieving this goal requires simplifying this complexity in synthetic axes defining a trait space and to cluster species based on their traits while identifying those with unique combinations of traits. However, so far, we know little about the dimensionality, the robustness to trait omission and the structure of these trait spaces. Here, we propose a unified framework and a synthesis across 30 trait datasets representing a broad variety of taxa, ecosystems and spatial scales to show that a common trade-off between trait space quality and operationality appears between three and six dimensions. The robustness to trait omission is generally low but highly variable among datasets. We also highlight invariant scaling relationships, whatever organismal complexity, between the number of clusters, the number of species in the dominant cluster and the number of unique species with total species richness. When species richness increases, the number of unique species saturates, whereas species tend to disproportionately pack in the richest cluster. Based on these results, we propose some rules of thumb to build species trait spaces and estimate subsequent functional diversity indices.

Rebound in functional distinctiveness following warming and reduced fishing in the North Sea.

Functionally distinct species (i.e. species with unique trait combinations in the community) can support important ecological roles and contribute disproportionately to ecosystem functioning. Yet, how functionally distinct species have responded to recent climate change and human exploitation has been widely overlooked. Here, using ecological traits and long-term fish data in the North Sea, we identified functionally distinct and functionally common species, and evaluated their spatial and temporal dynamics in relation to environmental variables and fishing pressure. Functionally distinct species were characterized by late sexual maturity, few, large offspring, and high parental care, many being sharks and skates that play critical roles in structuring food webs. Both functionally distinct and functionally common species increased in abundance as ocean temperatures warmed and fishing pressure decreased over the last three decades; however, functionally distinct species increased throughout the North Sea, but primarily in southern North Sea where fishing was historically most intense, indicating a rebound following fleet decommissioning and reduced harvesting. Yet, some of the most functionally distinct species are currently listed as threatened by the IUCN and considered highly vulnerable to fishing pressure. Alarmingly these species have not rebounded. This work highlights the relevance and potential of integrating functional distinctiveness into ecosystem management and conservation prioritization.

The future of fish-based ecological assessment of European rivers: from traditional EU Water Framework Directive compliant methods to eDNA metabarcoding-based approaches.

Most of the present EU Water Framework Directive (WFD) compliant fish-based assessment methods of European rivers are multi-metric indices computed from traditional electrofishing (TEF) samples, but this method has known shortcomings, especially in large rivers. The probability of detecting rare species remains limited, which can alter the sensitivity of the indices. In recent years, environmental (e)DNA metabarcoding techniques have progressed sufficiently to allow applications in various ecological domains as well as eDNA-based ecological assessment methods. A review of the 25 current WFD-compliant methods for river fish shows that 81% of the metrics used in these methods are expressed in richness or relative abundance and thus compatible with eDNA samples. However, more than half of the member states' methods include at least one metric related to age or size structure and would have to adapt their current fish index if reliant solely on eDNA-derived information. Most trait-based metrics expressed in richness are higher when computed from eDNA than when computed from TEF samples. Comparable values are obtained only when the TEF sampling effort increases. Depending on the species trait considered, most trait-based metrics expressed in relative abundance are significantly higher for eDNA than for TEF samples or vice versa due to over-estimation of sub-surface species or under-estimation of benthic and rare species by TEF sampling, respectively. An existing predictive fish index, adapted to make it compatible with eDNA data, delivers an ecological assessment comparable with the current approved method for 22 of the 25 sites tested. Its associated uncertainty is lower than that of current fish indices. Recommendations for the development of future fish eDNA-based indices and the associated eDNA water sampling strategy are discussed.

An innovative bivariate approach to detect joint temporal trends in environmental conditions: Application to large French rivers and diadromous fish.

Most key life-events of organisms are synchronized by complex interactions of several environmental cues to ensure optimal survival and growth of individuals and their offspring. However, global change is known to affect multiple components of ecosystems and cues at the same time. Therefore, detecting joint trends in covariate time series is a crucial challenge in global change ecology that has rarely been addressed so far. In this context, we designed an innovative combination of kernel density estimations and Mann-Kendall trend tests to detect joint temporal trends in a pair of environmental variables. This methodological framework was tested on >30 years (1976-2019) of water temperature and discharge data for 6 large French rivers (the Garonne, Dordogne, Rhône, Rhine, Loire and Vienne rivers). The implications of such trends in both temperature and discharge for diadromous species key life-cycle processes were then explored by checking if significant bivariate environmental changes occurred during seasons of upstream and downstream migration, and reproductive activities. Results were contrasted between rivers and seasons: many rivers displayed an increase in the number of days with high water temperature and low river discharge, but local discharge regulation measures could have mitigated the trend in discharge. Our findings showed that species migrating or spawning in spring were likely to be strongly impacted by the new environmental conditions in the Garonne, Loire and Rhône rivers, given the marked changes in water temperature and discharge associations detected by our new method. Conditions experienced by fall-running and spawning species have been strongly affected in all the rivers studied. This innovative methodology was implemented in a new R package, ChocR, for application to other environments and ecosystems.

The effect of lateral connectedness on the taxonomic and functional structure of fish communities in a lowland river floodplain.

In river floodplains many conservation programs focus on the main river channel as the richest in species. Lateral floodplain waterbodies, which contribute largely to functional processes in river systems, often remain overlooked and exposed to anthropogenic pressures. Although the role of hydrological connectedness between lateral waterbodies and the main river on taxonomic composition of fish communities is well understood, effects on functional community composition is much less studied. Abundance data of fish communities were gathered from 152 electrofishing sites in the main channel and lateral floodplain waterbodies of the river Lippe (Germany), over 18 years. These data were used to compare taxonomic, functional, conservation and recreational fishing aspects along the floodplain lateral connectedness gradient. Fish species richness decreased along the lateral continuum from the main river channel to isolated floodplain waterbodies. In contrast, the relative abundance of endangered and also of non-native species increased along this gradient, highlighting the ecological and conservational importance of floodplain waterbodies. Species composition in floodplain waterbodies differed across the connectedness gradient showing distinct assemblages which were not merely subsets of the main channel. The variability of life-history and feeding strategists among classes of lateral connectedness confirmed the importance of each connectivity class in contributing to the overall floodplain functional diversity. This study highlights the need of preserving fish taxonomic and functional biodiversity across the floodplain as one integrated hydrosystem. Conservation and restoration measures should therefore extend to include the whole floodplain area and the complete spectrum of differently connected floodplain waterbodies in addition to the main channel of the river.

Otolith shape analysis and daily increment validation during ontogeny of larval and juvenile European chub Squalius cephalus.

This assesses features of otoliths from laboratory-reared embryos, larvae and juvenile European chub Squalius cephalus from hatching to 180 days post-hatching (dph). We observed the development of the three pairs of otoliths (lapilli, sagittae and asterisci) and more precisely shape changes, as well as timing and deposition rate of increments of the lapilli. The lapilli and the sagittae were present at hatching, whereas the asterisci formed between 20 and 30 dph. The lapillus and sagitta shapes were round until 20 dph. From 60 dph the anterior and the posterior rostra of the sagittae were well developed, but very thin, making this otolith too fragile to manipulate for further studies of shape and validation of otolith increment deposition rate. The lapilli provided reliable age estimates for free embryos, larvae and juveniles up to 120 dph. However, caution should be taken when ageing fish older than 150 dph as an underestimation was noticeable. The regression of the number of otolith increments on age showed a slope and an intercept not significantly different from 1 and 0, respectively, which indicated that otolith growth increments were deposited on a daily basis, with the first microincrement occurring at hatching. Increment counts were consistent between three interpreters, indicating a consistent and reliable age estimate. This study validates that the otolith increment deposition rate can be used to assess hatching dates and daily growth of wild S. cephalus under 150 dph and in environments similar to the conditions used in this study.

Effects of a broad range of experimental temperatures on the population growth and body-size of five species of free-living nematodes.

The temperature-size rule postulates that the growth rates of ectotherms increase under rising temperatures, while the sizes of these organisms at maturity decrease. However, the upper temperature-tolerance range is also typically represented by a metabolic tipping point after which growth suddenly ceases. Free-living nematodes are important members of ecosystems, but little is known about their thermal tolerance. In the present study we measured the population growth rates and body-size distributions of five species of free-living bacterivorous nematodes exposed in the laboratory to a broad range of temperatures. This allowed a determination of their different thermal tolerance ranges, even of closely related species, including Plectus acuminatus (thermal optimum of 20-25 °C) and P. cf. velox (10-15 °C). With the exception of Acrobeloides nanus, which had the broadest thermal tolerance range, the population growth of the other species declined between 25 and 30 °C. Our results were consistent with the temperature-size rule, as the body-size of the tested species at maturity decreased with increasing temperature. This reduction was accompanied by a smaller number of eggs carried by mature females. Although our study was purely experimental, it suggests that heat waves or other alterations in the thermal regime affect the population dynamics and body-size structure of nematode communities in the field.

Structural and functional responses of invertebrate communities to climate change and flow regulation in alpine catchments.

Understanding and predicting how biological communities respond to climate change is critical for assessing biodiversity vulnerability and guiding conservation efforts. Glacier- and snow-fed rivers are one of the most sensitive ecosystems to climate change, and can provide early warning of wider-scale changes. These rivers are frequently used for hydropower production but there is minimal understanding of how biological communities are influenced by climate change in a context of flow regulation. This study sheds light on this issue by disentangling structural (water temperature preference, taxonomic composition, alpha, beta and gamma diversities) and functional (functional traits, diversity, richness, evenness, dispersion and redundancy) effects of climate change in interaction with flow regulation in the Alps. For this, we compared environmental and aquatic invertebrate data collected in the 1970s and 2010s in regulated and unregulated alpine catchments. We hypothesized a replacement of cold-adapted species by warming-tolerant ones, high temporal and spatial turnover in taxa and trait composition, along with reduced taxonomic and functional diversities in consequence of climate change. We expected communities in regulated rivers to respond more drastically due to additive or synergistic effects between flow regulation and climate change. We found divergent structural but convergent functional responses between free-flowing and regulated catchments. Although cold-adapted taxa decreased in both of them, greater colonization and spread of thermophilic species was found in the free-flowing one, resulting in higher spatial and temporal turnover. Since the 1970s, taxonomic diversity increased in the free flowing but decreased in the regulated catchment due to biotic homogenization. Colonization by taxa with new functional strategies (i.e. multivoltine taxa with small body size, resistance forms, aerial dispersion and reproduction by clutches) increased functional diversity but decreased functional redundancy through time. These functional changes could jeopardize the ability of aquatic communities facing intensification of ongoing climate change or new anthropogenic disturbances.

Environmental DNA reveals quantitative patterns of fish biodiversity in large rivers despite its downstream transportation.

Despite the ecological and societal importance of large rivers, fish sampling remains costly and limited to specific habitats (e.g., river banks). Using an eDNA metabarcoding approach, we regularly sampled 500 km of a large river (Rhône River). Comparisons with long-term electrofishing surveys demonstrated the ability of eDNA metabarcoding to qualitatively and quantitatively reveal fish assemblage structures (relative species abundance) but eDNA integrated a larger space than the classical sampling location. Combination of a literature review and field data showed that eDNA behaves in the water column like fine particulate organic matter. Its detection distance varied from a few km in a small stream to more than 100 km in a large river. To our knowledge, our results are the first demonstration of the capacity of eDNA metabarcoding to describe longitudinal fish assemblage patterns in a large river, and metabarcoding appears to be a reliable, cost-effective method for future monitoring.

Can Recent Global Changes Explain the Dramatic Range Contraction of an Endangered Semi-Aquatic Mammal Species in the French Pyrenees?

Species distribution models (SDMs) are the main tool to predict global change impacts on species ranges. Climate change alone is frequently considered, but in freshwater ecosystems, hydrology is a key driver of the ecology of aquatic species. At large scale, hydrology is however rarely accounted for, owing to the lack of detailed stream flow data. In this study, we developed an integrated modelling approach to simulate stream flow using the hydrological Soil and Water Assessment Tool (SWAT). Simulated stream flow was subsequently included as an input variable in SDMs along with topographic, hydrographic, climatic and land-cover descriptors. SDMs were applied to two temporally-distinct surveys of the distribution of the endangered Pyrenean desman (Galemys pyrenaicus) in the French Pyrenees: a historical one conducted from 1985 to 1992 and a current one carried out between 2011 and 2013. The model calibrated on historical data was also forecasted onto the current period to assess its ability to describe the distributional change of the Pyrenean desman that has been modelled in the recent years. First, we found that hydrological and climatic variables were the ones influencing the most the distribution of this species for both periods, emphasizing the importance of taking into account hydrology when SDMs are applied to aquatic species. Secondly, our results highlighted a strong range contraction of the Pyrenean desman in the French Pyrenees over the last 25 years. Given that this range contraction was under-estimated when the historical model was forecasted onto current conditions, this finding suggests that other drivers may be interacting with climate, hydrology and land-use changes. Our results imply major concerns for the conservation of this endemic semi-aquatic mammal since changes in climate and hydrology are expected to become more intense in the future.