Assessment of the impact of chemical pollution on endangered migratory fish in two major rivers of France, including spawning grounds

Water pollution is a one of the most contributors to aquatic biodiversity decline. Consequently, ecological risk assessment methods have been developed to investigate the effects of existing stresses on the environment, including the toxic effects of chemicals. One of the existing approaches to quantify toxic risks is called "Potentially Affected Fraction of species" (PAF), which estimates the potential loss of species within a group of species studied. In this study, the PAF method was applied to the Garonne catchment (southwest France) due to the limited information available on the involvement of water pollution in the decline of diadromous fish populations. This approach was used to quantify the potential toxic risk associated with chemical contamination of water for fish species. The objectives were to quantify this risk (1) in the Garonne and Dordogne rivers and (2) in the spawning grounds of two endangered anadromous fish species: the allis shad and the European sturgeon during the development period of their early life stages. Environmental pollution data was provided for 21 sites within the Garonne catchment between 2007 and 2022, and toxicity data was obtained specifically from freshwater toxicity tests on fish species. Then, for each site and each year, the potential toxic risk for a single substance (ssPAF) and for a mixture of substances (msPAF) was calculated and classified as high (>5 %), moderate (>1 % and < 5 %) or low (<1 %). Potential toxic risks were mostly moderate and mainly associated with: metals > other industrial pollutants and hygiene and care products > agrochemicals. In summary, this study highlights the probable involvement of water contamination on the decline, fate and restoration of diadromous fish populations in the Garonne catchment, focusing notably on the toxic effects on early life stages, a previously understudied topic.

Lipids as biomarkers to assess the nutritional and physiological status of two diadromous fish (Anguilla anguilla and Chelon auratus) at early life stages in a temperate macrotidal estuary

Estuaries are considered as key habitats for the early life stages of fish. However, in the face of massive destruction of many estuarine intertidal areas, management and conservation measures are needed. Fish condition indicators may be used as a proxy of habitat quality and provide valuable information for management of coastal areas. In this study, the larvae of golden mullet (Chelon auratus) and European glass eels (Anguilla anguilla) were sampled in three sites of the Gironde Estuary. Different lipid classes and fatty acids were quantified: phospholipids (globally, phosphatidylethanolamine and phosphatidylcholine), triglycerides, omega-3 (particularly docosahexaenoic and eicosapentaenoic acids), omega-6 and C18:1. These biomarkers provide information on the nutritional status of the larvae as well as on prey availability and larvae diet between sites. One site significantly differed from the others as it seemed to offer abundant and better-quality prey. The very high levels of omega-3 contained in mullet larvae suggested that this site provided a high amount of diatoms. However, the mullet larvae that colonized this site also showed physiological stress that could be explained by exposure to pollutants through their prey. This work constitutes an essential baseline for developing biomarkers to assess the quality of habitats in a global change context.

Identifying the impact of toxicity on stream macroinvertebrate communities in a multi-stressor context based on national ecological and ecotoxicological monitoring databases

In situ bioassays are used to measure the harmful effects induced by mixtures of toxic chemicals in watercourses. In France, national-scale biomonitoring data are available including invertebrate surveys and in-field chemical toxicity measures with caged gammarids to assess environmental toxicity of mixtures of chemicals. The main objective of our study is to present a proof-of-concept approach identifying possible links between in-field chemical toxicity, stressors and the ecological status. We used two active biomonitoring databases comprising lethal toxicity (222 in situ measures of gammarid mortality) and sublethal toxicity (101 in situ measures of feeding inhibition). We measured the ecological status of each active biomonitoring site using the I₂M₂ metric (macroinvertebrate-based multimetric index), accounted for known stressors of nutrients and organic matter, hydromorphology and chemical toxicity. We observed a negative relationship between stressors (hydromorphology, nutrients and organic matter, and chemical toxicity) and the good ecological status. This relationship was aggravated in watercourses where toxicity indicators were degraded. We validated this hypothesis for instance with nutrients and organic matter like nitrates or hydromorphological conditions like percentage of vegetation on banks. Future international assesments concerning the role of in-field toxic pollution on the ecological status in a multi-stressor context are now possible via the current methodology.

Decadal biodiversity trends in rivers reveal recent community rearrangements

While it is recognized that biodiversity currently declines at a global scale, we still have an incomplete understanding of local biodiversity trends under global change. To address this deficiency, we examined the recent decadal trends in water quality and biodiversity (taxonomic and functional) of key river organisms (diatoms, macroinvertebrates and fish) in France. We implemented regression, RLQ and fourth-corner analyses. Our results showed that nutrient loads tended to decrease, diatom richness tended to decline and macoinvertebrate richness tended to increase. The recovery of sensitive taxa in all three groups suggested a successful outcome of water quality management in France over the past decades. Our study further revealed consistent rearrangements within river communities, with a decrease in the ratio of planktonic to benthic diatoms, and corresponding functional changes in macroinvertebrate and fish trait composition, indicative of a trophic cascade in response to changes in environmental conditions.

Quantification of multi-scale links of anthropogenic pressures with PAH and PCB bioavailable contamination in French freshwaters

Aquatic ecosystems are exposed to multiple environmental pressures including chemical contamination. Polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) are persistent organic pollutants (POPs) known as preoccupying substances for the environment. Active biomonitoring (ABM) is a surveillance method for polluted aquatic ecosystems measuring bioavailable contamination. In this work, the aim was to quantify the total links between environmental pressures and bioavailable contamination (for PAHs and PCBs) at the French national scale. Based on 245 sites experimented by ABM from 2017 to 2019, environmental pressures (anthropogenic pressures and environmental parameters) were defined (point source landfill density, point source urban density, point source industry density, point source road density, nonpoint source industry density, nonpoint source road density, nonpoint source urban density, nutrients and organic matter, slope, dams, straightness, coarse sediment, summer precipitation, hydrographic network density and watershed size) and characterized by one or a combination of measures called stressor indicators. The links between environmental pressures and bioavailable POPs contamination (ABM measure) at a large spatial scale were defined and quantified via structural equation modeling. Point source urban density, nutrients and organic matter, summer precipitation, straightness and point source industry density are correlated positively with PAH bioavailable contamination. In contrast, nonpoint source urban density, nonpoint source industry density, nonpoint source road density and watershed size are positively correlated with PCB bioavailable contamination. The dominant pressures linked to PAHs and PCBs were different, respectively local and large-scale pressures were linked to PAH bioavailable contamination, and only large-scale pressures were linked to PCB bioavailable contamination.

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.

How to quantify the links between bioavailable contamination in watercourses and pressures of anthropogenic land cover, contamination sources and hydromorphology at multiple scales?

Active biomonitoring permits the quantification of biological exposure to chemicals through measurements of bioavailable concentrations in biota and biological markers of toxicity in organisms. It enables respective comparison of the levels of contamination between sites and sampling campaigns. Caged gammarids are recently proposed as relevant probes for measuring bioavailable contamination in freshwater systems. The purpose of the present study was to develop a multi-pressure and multiscale approach, considering metallic contamination levels (from data based on active biomonitoring) as a response to pressures (combination of individual stressors). These pressures were anthropogenic land cover, industry density, wastewater treatment plant density, pressures on stream hydromorphological functioning, riverside vegetation and bioavailability factors. A dataset combining active biomonitoring and potentially related pressures was established at the French national scale, with 196 samplings from 2009 to 2016. The links between pressures and metallic contamination were defined and modelled via structural equation modeling (more specifically partial least squares - path modeling). The model enabled the understanding of the respective influences of pressures on metallic bioconcentration in caged sentinel organisms. Beyond validating the local influence of industries and wastewater treatment plants on metallic contamination, this model showed a complementary effect of driving forces of anthropogenic land cover (leading to human activities). It also quantified a significant influence of pressures on stream hydromorphological functioning, presence of vegetation and physico-chemical parameters on metal bioconcentration. This hierarchical multi-pressure approach could serve as a concept on how pressures and contamination (assessed by active biomonitoring) can be connected. Its future application will enable better understanding of environmental pressures leading to contamination in freshwater ecosystems.

Multisubstance Indicators Based on Caged Gammarus Bioaccumulation Reveal the Influence of Chemical Contamination on Stream Macroinvertebrate Abundances across France

Most anthropogenic stressors affecting freshwater systems are qualitatively known. However, the quantitative assessment of contaminant exposure and effects to aquatic communities is still difficult, limiting the understanding of consequences on aquatic ecosystem functioning and the implementation of effective management plans. Here, multisubstance indicators based on caged gammarid bioaccumulated contamination data are proposed (for metals and persistent organic pollutants, POPs) to map the bioavailable contamination level of freshwater ecosystems at a large spatial scale. We assessed the ability of these indicators to highlight the relationships between chemical exposure gradients and alteration in the abundance of macroinvertebrate populations on a data set of 218 watercourses distributed throughout France. We identified spatial regional heterogeneities in the levels of bioavailable contamination of metals (18 compounds) and POPs (43 compounds). Besides this, a degradation of Gammaridae, Ephemeridae, and Hydrobiidae densities with increasing levels of metal contamination are identified relative to Baetidae, Chironomidae, and Hydropsychidae. We show here that active biomonitoring allows the establishment of multisubstance indicators of bioavailable contamination, which reliably quantify chemical exposure gradients in freshwater ecosystems. Our ability to identify species-specific responses to chemical exposure gradients demonstrates the promising possibility to further decipher the effects of chemical contamination on macroinvertebrate assemblages through this type of indicator.

Direct and indirect effects of multiple stressors on stream invertebrates across watershed, reach and site scales: A structural equation modelling better informing on hydromorphological impacts

The purpose of our approach was to take into account the nested spatial scales driving stream functioning in the description of pressures/ecological status links by analysing the results of a hierarchical model. The development of this model has allowed us to answer the following questions: Does the consideration of the indirect links between anthropogenic pressures and stream ecological status modify the hierarchy of pressure types impacting benthic invertebrates? Do the different nested scales play different roles in the anthropogenic pressures/ecological status relationship? Does this model lead to better understanding of the specific role of hydromorphology in the evaluation of stream ecological status? To achieve that goal, we used the Partial Least Square (PLS) path modelling method to develop a structural model linking variables describing (i) land use and hydromorphological alterations at the watershed scale, (ii) hydromorphological alterations at the reach scale, (iii) nutrients-organic matter contamination levels at the site scale, and (iv) substrate characteristics at the sampling site scale, to explain variation in values of a macroinvertebrate-based multimetric index: the French I₂M₂. We have highlighted the importance of land use effects exerted on both hydromorphological and chemical characteristics of streams observed at finer scales and their subsequent indirect impact on stream ecological status. Hydromorphological alterations have an effect on the substrate mosaic structure and on the concentrations of nutrients and organic matter at site scale. This result implies that stream hydromorphology can have a major indirect effect on macroinvertebrate assemblages and that the hierarchy of impacts of anthropogenic pressures on stream ecological status generally described in the literature - often determining strategic restoration priorities - has to be re-examined. Finally, the effects of nutrients and organic matter on macroinvertebrate assemblages are lower than expected when all the indirect effects of land use and hydromorphological alterations are taken into account.

Unravelling river system impairments in stream networks with an integrated risk approach

Rivers are complex systems for which it is hard to make reliable assessments of causes and responses to impairments. We present a holistic risk-based framework for river ecosystem assessment integrating all potential intervening processes and functions. Risk approaches allow us to deal with uncertainty both in the construction of indicators for magnitude of stressors and in the inference of environmental processes and their impairment. Yet, here we go further than simply replacing uncertainty by a risk factor. We introduce a more accurate and rigorous notion of risk with a transcription of uncertainty in causal relationships in probability distributions for the magnitude of impairment and the weight of different descriptors, with an associated confidence in the diagnostic. We discuss how Bayesian belief networks and Bayesian hierarchical inference allow us to deal with this risk concept to predict impairments and potential recovery of river ecosystems. We developed a comprehensive approach for river ecosystem assessment, which offers an appealing tool to facilitate diagnosis of the likely causes of impairment and predict future conditions. The ability of the risk approaches to integrate multi-scale quantitative and qualitative descriptors in the identification of multiple stressor sources and pathways in the stream network, and their impairment of specific processes and structures is illustrated for the national-level risk analysis for hydromorphology and pesticide pollution. Not only does the risk-based framework provide a more complete picture of environmental impairments, but it also offers a comprehensive, user-friendly tool to instruct the decision process.

Assessing the ecological status in the context of the European Water Framework Directive: where do we go now?

The Water Framework Directive (WFD) is now well established as the key management imperative in river basins across Europe. However, there remain significant concerns with the way WFD is implemented and there is now a need for water managers and scientists to communicate better in order to find solutions to these concerns. To address this, a Science-Policy Interface (SPI) activity was launched in 2010 led by Directorate-General for Research and Innovation and Onema (the French national agency for water and aquatic ecosystems), which provided an interactive forum to connect scientists and WFD end-users. One major aim of the SPI activity was to establish a list of the most crucial research and development needs for enhancing WFD implementation. This paper synthesises the recommendations from this event highlighting 10 priority issues relating to ecological status. For lakes, temporary streams and transitional and coastal waters, WFD implementation still suffers from a lack of WFD-compliant bioassessment methods. For rivers, special attention is required to assess the ecological impacts of hydromorphological alterations on biological communities, notably those affecting river continuity and riparian covering. Spatial extrapolation tools are needed in order to evaluate ecological status for water bodies for which no data are available. The need for more functional bioassessment tools as complements to usual WFD-compliant tools, and to connect clearly good ecological state, biodiversity and ecosystem services when implementing WFD were also identified as crucial issues.

Correlations between dioxin-like and indicators PCBs: potential consequences for environmental studies involving fish or sediment

Among the numerous PCB congeners, most of the dioxin-like PCBs (DL-PCBs) need to be characterized by hyphenated techniques. It has been shown in several instances that these congeners are well related to the total PCB content in fish. We examined datasets collected mainly in France, on freshwater and marine fish and sediments. A statistical model linking DL- and indicator PCBs was developed for a dataset composed of freshwater fishes, and proved to predict well DL-PCBs from indicator PCBs in all other fish sets, including marine ones. Type II error rates remained low in almost all fish sets. A similar correlation was observed in sediments. Non-dioxin-like PCBs elicit various adverse effects and represent 95% of the total PCBs. A European guideline for them is needed; the correlation between DL- and indicator PCBs could help develop this standard in the future.