Food web modeling of a river ecosystem for risk assessment of down-the-drain chemicals: a case study with AQUATOX

Community ecological response to polycyclic aromatic hydrocarbons in Baiyangdian Lake based on an ecological model

The dynamic response of a single population to chemicals can be represented by a Weibull function. However, it is unclear whether the overall response can still be represented in this manner when scaled up to the community level. In this study, we investigated the responses of biological communities to polycyclic aromatic hydrocarbons by using an ecological model of Baiyangdian Lake in northern China. The community dynamics process was divided into the following three stages. In the first stage, toxicity, played a dominant role and strong, medium, and weak species responses were observed according to the toxicity sensitivity. In the second stage, the dynamic process was dominated by the interaction strength with three alternative dynamic pathways comprising of direct response, no response, or inverse response. In the third stage, the toxicity was again dominant, and the biomasses of all species decreased to extinction. The toxicological dynamics were far more complex at the community level than those at the single species level and they were also influenced by the interaction strength as well as toxicity. The toxicological dynamic process in the community was constantly driven by the competing effects of these two forces. In addition to the total biomass, the interaction strength was identified as a suitable community-level signal because it exhibited good indicator properties regarding ecosystem steady-state transitions. However, we found that food web stability indicators were not suitable for use as community-level signals because they were not sensitive to changes in the ecosystem state. Some ecological management suggestions have been proposed, including medium to long-term monitoring, and reduction of external pollution loads and bioindicators. The results obtained in this study increase our understanding of how chemicals interfere with community dynamics, and the interaction strength and total biomass were identified as useful holistic indicators.

Optimal Method for Biomass Estimation in a Cladoceran Species, Daphnia Magna (Straus, 1820): Evaluating Length–Weight Regression Equations and Deriving Estimation Equations Using Body Length, Width and Lateral Area

Assessing the biomass of zooplankton compensates for the difference between number of individuals and the accumulated body weight of the community, which helps assess aquatic ecosystem food web functions. Daphnia are crustaceans that play an intermediate role in biological interactions within food webs. The morphology and body specification of Daphnia differ during growth; hence, it is essential to apply species-specific equations to estimate biomass. We evaluated the length–weight regression equations used previously to estimate Daphnia magna biomass and conducted regression analyses using various body specifications and biomass measurements taken directly using devices such as a microbalance and microscopic camera. Biomass estimated using an equation from the Environmental Protection Agency was significantly different from the direct measurement: average biomass was lower, indicating that the equation possibly underestimated actual biomass. The biomass of D. magna had a higher multiple R2 value when length was compared with width and area, and a linear regression equation was the most suitable equation for biomass estimation. Because body specifications and biomass are affected by various environmental factors, the development of accurate species-specific biomass estimation equations will contribute to obtaining fundamental data with which the biological responses of zooplankton to aquatic ecosystem changes can be assessed.

Prioritizing the Effects of Emerging Contaminants on Estuarine Production under Global Warming Scenarios

Due to non-linear interactions, the effects of contaminant mixtures on aquatic ecosystems are difficult to assess, especially under temperature rise that will likely exacerbate the complexity of the responses. Yet, under the current climatic crisis, assessing the effects of water contaminants and temperature is paramount to understanding the biological impacts of mixtures of stressors on aquatic ecosystems. Here, we use an ecosystem model followed by global sensitivity analysis (GSA) to prioritize the effects of four single emerging contaminants (ECs) and their mixture, combined with two temperature rise scenarios, on the biomass production of a NE Atlantic estuary. Scenarios ran for 10 years with a time-step of 0.1 days. The results indicate that macroinvertebrate biomass was significantly explained by the effect of each single EC and by their mixture but not by temperature. Globally, the most adverse effects were induced by two ECs and by the mixture of the four ECs, although the sensitivity of macroinvertebrates to the tested scenarios differed. Overall, the present approach is useful to prioritize the effects of stressors and assess the sensitivity of the different trophic groups within food webs, which may be of relevance to support decision making linked to the sustainable management of estuaries and other aquatic systems.

Assessing chemical risk within an ecosystem services framework: Implementation and added value

An ecosystem services (ES) approach to chemical risk assessment has many potential advantages, but there are also substantial challenges regarding its implementation. We report the findings of a multi-stakeholder workshop that evaluated the feasibility of adopting an ES approach to chemical risk assessment using currently available tools and data. Also evaluated is the added value such an approach would bring to environmental decision making. The aim was to build consensus across disparate stakeholders and to co-produce a common understanding of the regulatory benefits and feasibility of implementing an ES approach in European chemicals regulation. Workshop discussions were informed by proof of concept studies and resulted in the development of a novel tiered framework for assessing chemical risk to ES delivery. There was consensus on the substantial added value of adopting an ES-based approach for regulatory decision making. Ecosystem services provide a common currency and a 'unifying approach' across environmental compartments, stressors and regulatory frameworks. The ES approach informs prioritisation of risk and remedial action and aids risk communication and risk management. It facilitates a more holistic assessment, enables ES trade-offs to be compared across alternative interventions, and supports comparative risk assessments and a socio-economic analysis of management options and decisions. Key to realising this added value is a shift away from using a single threshold value to categorise risk, towards a consideration of the exposure-effect distribution for individual ES of interest. Also required is the development of an integrated systems-level approach across regulatory frameworks and agreement on specific protection goals and scenarios for framing environmental risk assessments. The need to further develop tools for extrapolating toxicity data to service providers and ES delivery, including logic chains and ecological production functions, was highlighted. Also agreed was the need for methods and metrics for ES valuation to be used in assessing trade-offs.

How Do Indirect Effects of Contaminants Inform Ecotoxicology? A Review

Indirect effects in ecotoxicology are defined as chemical- or pollutant-induced alterations in the density or behavior of sensitive species that have cascading effects on tolerant species in natural systems. As a result, species interaction networks (e.g., interactions associated with predation or competition) may be altered in such a way as to bring about large changes in populations and/or communities that may further cascade to disrupt ecosystem function and services. Field studies and experimental outcomes as well as models indicate that indirect effects are most likely to occur in communities in which the strength of interactions and the sensitivity to contaminants differ markedly among species, and that indirect effects will vary over space and time as species composition, trophic structure, and environmental factors vary. However, knowledge of indirect effects is essential to improve understanding of the potential for chemical harm in natural systems. For example, indirect effects may confound laboratory-based ecological risk assessment by enhancing, masking, or spuriously indicating the direct effect of chemical contaminants. Progress to better anticipate and interpret the significance of indirect effects will be made as monitoring programs and long-term ecological research are conducted that facilitate critical experimental field and mesocosm investigations, and as chemical transport and fate models, individual-based direct effects models, and ecosystem/food web models continue to be improved and become better integrated.

Coupling of the AQUATOX and EFDC Models for Ecological Impact Assessment of Chemical Spill Scenarios in the Jeonju River, Korea

In this study, an ecological impact was assessed for the short-term leak scenario through the AQUATOX-EFDC model, which combines the proven ecological model AQUATOX with the hydrodynamic model EFDC. A case study of the coupled AQUATOX-EFDC model was conducted for 30-30,000 kg toluene leak scenarios in the Jeonju River in South Korea. A 21-day scenario simulation was conducted, and the impact of the toluene spill accident was evaluated by comparing the biomass between the control simulation and the perturbed simulation. As a result of the simulation, it was found that in the scenario in which 3000 kg of toluene was leaked for a day, a substantial change was expected in the range of 0-640 m from the accident site. Additionally, for a 30,000 kg leak, a substantial change was expected in the range of 0-2300 m from the accident site, and the greatest damage was observed for the fish species group, the top predators. As a result, the AQUATOX-EFDC simulation showed a significant ecological impact, and the proposed model will be helpful to understand the ecological impact and establish the management strategy for the ecological risk of the chemical spill.

The bioaccumulation, elimination, and trophic transfer of BDE-47 in the aquatic food chain of Chlorella pyrenoidosa-Daphnia magna

As a persistent organic pollutant, 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) has been widely detected in aquatic environments. However, studies on the fate and transfer of BDE-47 in the aquatic food chain remain scarce. In this study, we investigated the bioaccumulation and elimination of BDE-47 in Chlorella pyrenoidosa, as well as the trophic transfer and biomagnification of BDE-47 in the "C. pyrenoidosa-Daphnia magna" food chain, using C-14 radioactive tracer technology. After 96 h of BDE-47 exposure, the algae accumulated 88.98% ± 0.59% of the initial radioactivity from the medium, and 36.09% ± 9.22% of the accumulated residues in the algae occurred in the form of bound residues. During 96 h of elimination, only 13% ± 0.50% of accumulated radioactivity in the algae was released into the medium. After 24 h of exposure, D. magna accumulated 35.99% ± 2.55% of the initial radioactivity via water filtration from the medium, and 31.35% ± 1.92% of the accumulated radioactivity in D. magna occurred as bound residues. However, D. magna accumulated 66.89% ± 2.37% of the accumulated radioactivity in the algae via food uptake from the contaminated algae, with a high portion of radioactivity observed as bound residues (83.40% ± 0.97% of accumulated radioactivity in D. magna). This indicated a reduction in the environmental risk of BDE-47. There was obvious biomagnification in the food chain between C. pyrenoidosa and D. magna (biomagnification factors, BMFs>1), resulting in environmental hazard transfer in the aquatic food chain. However, no metabolite was found during the exposure experiment, and further studies should be carried out to investigate the intrinsic mechanisms of the trophic transfer of BDE-47, especially in multilevel food chains. Therefore, this study elucidated the effect of dietary uptake on the bioaccumulation of BDE-47 in D. magna and provided new insight for future analysis regarding the bioaccumulation and biomagnification of organic pollutants in the food chain.

Priorities and opportunities in the application of the ecosystem services concept in risk assessment for chemicals in the environment

The ecosystem services approach has gained broad interest in regulatory and policy circles for use in ecological risk assessment. Whilst identifying several challenges, scientific experts from European regulatory authorities, the chemical industry and academia considered the approach applicable to all chemical sectors and potentially contributing to greater ecological relevance for setting and assessing environmental protection goals compared to current European regulatory frameworks for chemicals. These challenges were addressed in workshops to develop a common understanding across stakeholders on how the ecosystem services concept might be used in chemical risk assessment and what would need to be done to implement it. This paper describes the consensus outcome of those discussions. Knowledge gaps and research needs were identified and prioritised, exploring the use of novel approaches from ecology, ecotoxicology and ecological modelling. Where applicable, distinction is made between prospective and retrospective ecological risk assessment. For prospective risk assessment the development of environmental scenarios accounting for chemical exposure and ecological conditions was designated as a top priority. For retrospective risk assessment the top priority research need was development of reference conditions for key ecosystem services and guidance for their derivation. Both prospective and retrospective risk assessment would benefit from guidance on the taxa and measurement endpoints relevant to specific ecosystem services and from improved understanding of the relationships between measurement endpoints from standard toxicity tests and the ecosystem services of interest (i.e. assessment endpoints). The development of mechanistic models, which could serve as ecological production functions, was identified as a priority. A conceptual framework for future chemical risk assessment based on an ecosystem services approach is presented.

How do we take the pulse of an aquatic ecosystem? Current and historical approaches to measuring ecosystem integrity

Global environmental monitoring has indicated that the structure and function of some aquatic ecosystems has been significantly altered by human activities. There are many potential causes for these changes; however, one major concern is the increasing release of anthropogenic contaminants into aquatic environments. Although toxicological responses of individual organisms are typically well characterized, few studies have focused on characterizing toxicity at the ecosystem level. In fact, because of their scale and complexity, changes in ecosystem integrity are rarely considered in assessments of risks to ecosystems. This work attempts to move the conversation forward by defining integrity of ecosystems, reviewing current and historical approaches to measuring ecosystem integrity status (e.g., structural and functional measurements), and highlighting methods that could significantly contribute to the field of ecosystem toxicology (e.g., keystone species, environmental energetics, ecotoxicological modeling, and adverse outcome pathways [AOPs]). Through a critical analysis of current and historical methodologies, the present study offers a comprehensive, conceptual framework for the assessment of risks of contaminant exposure for whole ecosystems and proposes steps to facilitate better diagnoses of the integrity of aquatic systems. Environ Toxicol Chem 2019;38:289-301. © 2018 SETAC.

Simulation-based assessment of the impact of fertiliser and herbicide application on freshwater ecosystems at the Three Gorges Reservoir in China

Dams have profound impacts on river ecosystems, amongst them inundation of land, altered dynamics of the water body or uprising reservoir backwaters influencing tributary or upstream river sections. Along the outstandingly ecologically important Yangtze River in China, the Three Gorges Reservoir (TGR) is the largest project, covering an area of 1080 km². From the beginning, the dam-project came in for criticism on increasing environmental risks due to sub-merging former industrial and urban areas. We simulated dynamics of biotic and abiotic components of the TGR ecosystem (trophic guilds of aquatic organisms, hydrodynamics, nutrients), as well as the behaviour of the herbicidal substance propanil and its metabolites 3,4-Dichloroaniline (DCA) and 3,3',4,4'-tetrachloroazoxybenzene (TCAB). A modelling environment, provided by the AQUATOX software, was adapted to the specific situation at a tributary reach to the Yangtze river 'Daning River'. As the simulated food web contained several interconnected trophic levels, a significant biomagnification of metabolites was demonstrated by our simulation studies. In particular, newly emerging stagnant downstream sections of tributaries exhibited high probabilities due to accumulating pesticides from upstream sources. The common problem of algal blooms in the TGR-region was addressed by dose-response simulation experiments with essential nutrients. Impacts on structure and abundance of populations of aquatic organisms were shown. However, even high nutrient loads resulted in only slight changes of densities of organisms of all trophic levels. Nevertheless, the probabilities for large-scale algal blooms affecting drinking water quality were considered low because of high flow velocities and discharge rates towards the Yangtze River. We see high potential of simulation-based assessments that provide information for risk managers dealing with whole catchment areas. They are put in the position to differentiate the magnitude of impacts of various factors and decide about the most effective remediation measures.

Deriving predicted no-effect concentrations (PNECs) for emerging contaminants in the river Po, Italy, using three approaches: Assessment factor, species sensitivity distribution and AQUATOX ecosystem modelling

Over the past decades, per- and polyfluoroalkyl substances (PFASs) found in environmental matrices worldwide have raised concerns due to their toxicity, ubiquity and persistence. A widespread pollution of groundwater and surface waters caused by PFASs in Northern Italy has been recently discovered, becoming a major environmental issue, also because the exact risk for humans and nature posed by this contamination is unclear. Here, the Po River in Northern Italy was selected as a study area to assess the ecological risk posed by perfluoroalkyl acids (PFAAs), a class of PFASs, considering the noticeable concentration of various PFAAs detected in the Po waters over the past years. Moreover, the Po has a large environmental and socio-economic importance: it is the largest Italian river and drains a densely inhabited, intensely cultivated and heavily industrialized watershed. Predicted no-effect concentrations (PNECs) were derived using two regulated methodologies, assessment factors (AFs) and species sensitivity distribution (SSD), which rely on published ecotoxicological laboratory tests. Results were compared to those of a novel methodology using the mechanistic ecosystem model AQUATOX to compute PNECs in an ecologically-sound manner, i.e. considering physical, chemical, biological and ecological processes in the river. The model was used to quantify how the biomasses of the modelled taxa in the river food web deviated from natural conditions due to varying inputs of the chemicals. PNEC for each chemical was defined as the lowest chemical concentration causing a non-negligible yearly biomass loss for a simulated taxon with respect to a control simulation. The investigated PFAAs were Perfluorooctanoic acid (PFOA) and Perfluorooctanesulfonic acid (PFOS) as long-chained compounds, and Perfluorobutanoic acid (PFBA) and Perfluorobutanesulfonic acid (PFBS) as short-chained homologues. Two emerging contaminants, Linear Alkylbenzene Sulfonate (LAS) and triclosan, were also studied to assess the performance of the three methodologies for chemicals whose ecotoxicology and environmental fate are well-studied. The most precautionary approach was the use of AFs generally followed by SSD and then AQUATOX, except for PFOS, for which AQUATOX yielded a much lower PNEC compared to the other approaches since, unlike the other two methodologies, it explicitly simulates sublethal toxicity and indirect ecological effects. Our findings highlight that neglecting the role of ecological processes when extrapolating from laboratory tests to ecosystems can result in under-protective threshold concentrations for chemicals. Ecosystem models can complement existing laboratory-based methodologies, and the use of multiple methods for deriving PNECs can help to clarify uncertainty in ecological risk estimates.

Responses of Aquatic Plants to Eutrophication in Rivers: A Revised Conceptual Model

Compared to research on eutrophication in lakes, there has been significantly less work carried out on rivers despite the importance of the topic. However, over the last decade, there has been a surge of interest in the response of aquatic plants to eutrophication in rivers. This is an area of applied research and the work has been driven by the widespread nature of the impacts and the significant opportunities for system remediation. A conceptual model has been put forward to describe how aquatic plants respond to eutrophication. Since the model was created, there have been substantial increases in our understanding of a number of the underlying processes. For example, we now know the threshold nutrient concentrations at which nutrients no longer limit algal growth. We also now know that the physical habitat template of rivers is a primary selector of aquatic plant communities. As such, nutrient enrichment impacts on aquatic plant communities are strongly influenced, both directly and indirectly, by physical habitat. A new conceptual model is proposed that incorporates these findings. The application of the model to management, system remediation, target setting, and our understanding of multi-stressor systems is discussed. We also look to the future and the potential for new numerical models to guide management.

Review: ecotoxicity of organic and organo-metallic antifouling co-biocides and implications for environmental hazard and risk assessments in aquatic ecosystems

Hazard assessments of Irgarol 1051, diuron, 2-(thiocyanomethylthio)benzothiazole (TCMTB), dichloro-octylisothiazolin (DCOIT), chlorothalonil, dichlofluanid, thiram, zinc pyrithione, copper pyrithione, triphenylborane pyridine (TPBP), capsaicin, nonivamide, tralopyril and medetomidine were performed to establish robust environmental quality standards (EQS), based on predicted no effect concentrations (PNECs). Microalgae, zooplankton, fish and amphibians were the most sensitive ecological groups to all the antifoulants evaluated, especially in the early life stages. No differences were identified between freshwater and seawater species. The use of toxicity tests with non-standard species is encouraged because they increase the datasets, allowing EQS to be derived from probabilistic-based PNECs whilst reducing uncertainties. The global ban of tributyltin (TBT) has been heralded as a major environmental success; however, substitute antifoulants may also pose risks to aquatic ecosystems. Environmental risk assessments (ERAs) have driven decision-makings for regulating antifouling products, but in many countries there is still a lack of regulation of antifouling biocides which should be addressed.

Ecological risk assessment of toxic organic pollutant and heavy metals in water and sediment from a landscape lake in Tianjin City, China

To estimate the ecological risk of toxic organic pollutant (formaldehyde) and heavy metals (mercury (Hg), arsenic (As), cadmium (Cd), and chromium (Cr)) in water and sediment from a landscape Lake in Tianjin City, an ecological risk assessment was performed. The risk quotient (RQ) method and the AQUATOX model were used to assess the ecological risk of formaldehyde in landscape water. Meanwhile, the RQ method and the potential ecological risk index method were used to assess the ecological risk of four heavy metals in water and sediment from the studied landscape lake, respectively. The results revealed that the maximum concentration of formaldehyde in landscape water was lower than the environmental quality standards of surface water in China. The maximum simulated concentrations of formaldehyde in phytoplankton and invertebrates were 3.15 and 22.91 μg/L, respectively, which were far less than its toxicity data values (1000 and 510 μg/L, respectively), suggesting that formaldehyde in landscape water was at a safe level for aquatic organisms. The RQ model indicated that the risks of phytoplankton and invertebrates were higher than that of fish posed by Hg and Cd in landscape water, and the risks from As and Cr were acceptable for all test organisms. Cd is the most important pollution factor among all heavy metals in sediment from studied landscape lake, and the pollution factor sequence of heavy metals was Hg > As > Cr > Cd. The values of risk index (RI) for four heavy metals in samples a and b were 43.48 and 72.66, which were much lower than the threshold value (150), suggesting that the ecological risk posed by heavy metals in sediment was negligible.

Ecological risk assessment of agricultural soils for the definition of soil screening values: A comparison between substance-based and matrix-based approaches

The Italian legislation on contaminated soils does not include the Ecological Risk Assessment (ERA) and this deficiency has important consequences for the sustainable management of agricultural soils. The present research compares the results of two ERA procedures applied to agriculture (i) one based on the "substance-based" approach and (ii) a second based on the "matrix-based" approach. In the former the soil screening values (SVs) for individual substances were derived according to institutional foreign guidelines. In the latter, the SVs characterizing the whole-matrix were derived originally by the authors by means of experimental activity. The results indicate that the "matrix-based" approach can be efficiently implemented in the Italian legislation for the ERA of agricultural soils. This method, if compared to the institutionalized "substance based" approach is (i) comparable in economic terms and in testing time, (ii) is site specific and assesses the real effect of the investigated soil on a battery of bioassays, (iii) accounts for phenomena that may radically modify the exposure of the organisms to the totality of contaminants and (iv) can be considered sufficiently conservative.

Toward refined environmental scenarios for ecological risk assessment of down-the-drain chemicals in freshwater environments

Current regulatory practice for chemical risk assessment suffers from the lack of realism in conventional frameworks. Despite significant advances in exposure and ecological effect modeling, the implementation of novel approaches as high-tier options for prospective regulatory risk assessment remains limited, particularly among general chemicals such as down-the-drain ingredients. While reviewing the current state of the art in environmental exposure and ecological effect modeling, we propose a scenario-based framework that enables a better integration of exposure and effect assessments in a tiered approach. Global- to catchment-scale spatially explicit exposure models can be used to identify areas of higher exposure and to generate ecologically relevant exposure information for input into effect models. Numerous examples of mechanistic ecological effect models demonstrate that it is technically feasible to extrapolate from individual-level effects to effects at higher levels of biological organization and from laboratory to environmental conditions. However, the data required to parameterize effect models that can embrace the complexity of ecosystems are large and require a targeted approach. Experimental efforts should, therefore, focus on vulnerable species and/or traits and ecological conditions of relevance. We outline key research needs to address the challenges that currently hinder the practical application of advanced model-based approaches to risk assessment of down-the-drain chemicals. Integr Environ Assess Manag 2017;13:233-248. © 2016 SETAC.

Integrated presentation of ecological risk from multiple stressors

Current environmental risk assessments (ERA) do not account explicitly for ecological factors (e.g. species composition, temperature or food availability) and multiple stressors. Assessing mixtures of chemical and ecological stressors is needed as well as accounting for variability in environmental conditions and uncertainty of data and models. Here we propose a novel probabilistic ERA framework to overcome these limitations, which focusses on visualising assessment outcomes by construct-ing and interpreting prevalence plots as a quantitative prediction of risk. Key components include environmental scenarios that integrate exposure and ecology, and ecological modelling of relevant endpoints to assess the effect of a combination of stressors. Our illustrative results demonstrate the importance of regional differences in environmental conditions and the confounding interactions of stressors. Using this framework and prevalence plots provides a risk-based approach that combines risk assessment and risk management in a meaningful way and presents a truly mechanistic alternative to the threshold approach. Even whilst research continues to improve the underlying models and data, regulators and decision makers can already use the framework and prevalence plots. The integration of multiple stressors, environmental conditions and variability makes ERA more relevant and realistic.

Removal of Antibiotics in Biological Wastewater Treatment Systems-A Critical Assessment Using the Activated Sludge Modeling Framework for Xenobiotics (ASM-X)

Many scientific studies present removal efficiencies for pharmaceuticals in laboratory-, pilot-, and full-scale wastewater treatment plants, based on observations that may be impacted by theoretical and methodological approaches used. In this Critical Review, we evaluated factors influencing observed removal efficiencies of three antibiotics (sulfamethoxazole, ciprofloxacin, tetracycline) in pilot- and full-scale biological treatment systems. Factors assessed include (i) retransformation to parent pharmaceuticals from e.g., conjugated metabolites and analogues, (ii) solid retention time (SRT), (iii) fractions sorbed onto solids, and (iv) dynamics in influent and effluent loading. A recently developed methodology was used, relying on the comparison of removal efficiency predictions (obtained with the Activated Sludge Model for Xenobiotics (ASM-X)) with representative measured data from literature. By applying this methodology, we demonstrated that (a) the elimination of sulfamethoxazole may be significantly underestimated when not considering retransformation from conjugated metabolites, depending on the type (urban or hospital) and size of upstream catchments; (b) operation at extended SRT may enhance antibiotic removal, as shown for sulfamethoxazole; (c) not accounting for fractions sorbed in influent and effluent solids may cause slight underestimation of ciprofloxacin removal efficiency. Using tetracycline as example substance, we ultimately evaluated implications of effluent dynamics and retransformation on environmental exposure and risk prediction.

Ecological relevance of biomarkers in monitoring studies of macro-invertebrates and fish in Mediterranean rivers

Mediterranean rivers are probably one of the most singular and endangered ecosystems worldwide due to the presence of many endemic species and a long history of anthropogenic impacts. Besides a conservation value per se, biodiversity is related to the services that ecosystems provide to society and the ability of these to cope with stressors, including climate change. Using macro-invertebrates and fish as sentinel organisms, this overview presents a synthesis of the state of the art in the application of biomarkers (stress and enzymatic responses, endocrine disruptors, trophic tracers, energy and bile metabolites, genotoxic indicators, histopathological and behavioural alterations, and genetic and cutting edge omic markers) to determine the causes and effects of anthropogenic stressors on the biodiversity of European Mediterranean rivers. We also discuss how a careful selection of sentinel species according to their ecological traits and the food-web structure of Mediterranean rivers could increase the ecological relevance of biomarker responses. Further, we provide suggestions to better harmonise ecological realism with experimental design in biomarker studies, including statistical analyses, which may also deliver a more comprehensible message to managers and policy makers. By keeping on the safe side the health status of populations of multiple-species in a community, we advocate to increase the resilience of fluvial ecosystems to face present and forecasted stressors. In conclusion, this review provides evidence that multi-biomarker approaches detect early signs of impairment in populations, and supports their incorporation in the standardised procedures of the Water Frame Work Directive to better appraise the status of European water bodies.