Theoretically exploring direct and indirect chemical effects across ecological and exposure scenarios using mechanistic fate and effects modelling

A prospective ecological risk assessment method based on exposure and ecological scenarios (ERA-EES) to determine soil ecological risks around metal mining areas

Soil heavy metal (HM) contamination around metal mining areas (MMAs) is a global concern that requires a cost-effective ecological risk assessment (ERA) method for preventive management. Traditional ERAs, comparing environmental HM concentrations with benchmarks, are labor- and cost-intensive in field investigations and chemical analyses, which challenge the management demands of numerous MMAs. In this study, a prospective ecological risk assessment method based on exposure and ecological scenario (ERA-EES) was developed to predict the eco-risk levels (low/medium/high) around MMAs prior to field sampling. Five exposure scenario indicators related to soil HM exposure and three ecological scenario indicators reflecting the soil bioreceptor response were selected and combined with the analytic hierarchy process and fuzzy comprehensive evaluation methods for ERA-EES development. Case application and performance evaluation with 67 MMAs in China demonstrated that the ERA-EES method had an overall effective and conservative performance when referring to potential ecological risk index (PERI) levels, with an accuracy of 0.87, kappa coefficient of 0.7, and low or medium eco-risk levels in PERI classified to high levels in ERA-EES. Overall, the selected scenario indicators could efficiently reflect the risk levels of soil HM pollution from mining activities. Besides, more regulatory efforts should be paid to the MMAs of nonferrous metals, underground and long-term mining and those located in southern China. This work provided a convenient and cost-effective prospective ERA method under the trend of ERA being tiered and refined, facilitating the risk management of various MMAs.

(Eco)toxicological tests for assessing impacts of chemical stress to aquatic ecosystems: Facts, challenges, and future

Monitoring of chemicals in the aquatic environment by chemical analysis alone cannot completely assess and predict the effects of chemicals on aquatic species and ecosystems. This is primarily because of the increasing number of (unknown) chemical stressors and mixture effects present in the environment. In addition, the ability of ecological indices to identify underlying stressors causing negative ecological effects is limited. Therefore, additional complementary methods are needed that can address the biological effects in a direct manner and provide a link to chemical exposure, i.e. (eco)toxicological tests. (Eco)toxicological tests are defined as test systems that expose biological components (cells, individuals, populations, communities) to (environmental mixtures of) chemicals to register biological effects. These tests measure responses at the sub-organismal (biomarkers and in vitro bioassays), whole-organismal, population, or community level. We performed a literature search to obtain a state-of-the-art overview of ecotoxicological tests available for assessing impacts of chemicals to aquatic biota and to reveal datagaps. In total, we included 509 biomarkers, 207 in vitro bioassays, 422 tests measuring biological effects at the whole-organismal level, and 78 tests at the population- community- and ecosystem-level. Tests at the whole-organismal level and biomarkers were most abundant for invertebrates and fish, whilst in vitro bioassays are mostly based on mammalian cell lines. Tests at the community- and ecosystem-level were almost missing for organisms other than microorganisms and algae. In addition, we provide an overview of the various extrapolation challenges faced in using data from these tests and suggest some forward looking perspectives. Although extrapolating the measured responses to relevant protection goals remains challenging, the combination of ecotoxicological experiments and models is key for a more comprehensive assessment of the effects of chemical stressors to aquatic ecosystems.

Estimating temporal and spatial levels of PAHs in air using rain samples and SPME analysis: Feasibility evaluation in an urban scenario

There is a growing interest in evaluating the role of concentration changes of contaminants in temporal and spatial gradients. This is often relevant for fast moving environmental phases such as air and water. In this paper, small volumes of rainwater were sampled as proxy for air concentrations of Polycyclic Aromatic Hydrocarbons (PAHs): rain was collected in three sampling sites (high traffic, restricted traffic and a low traffic zone) in Como. Solid phase micro extraction (SPME) was used for the extraction to reduce required sample volumes, allowing the acquisition of more samples in time. Rain samples highlighted a spatial and temporal variability along a traffic gradient in the Como city, especially for the most abundant PAH, e.g. phenanthrene. Air concentrations were then estimated from rain concentrations. The results show that this is a cheap and promising method, although requiring rainfall/snowfall conditions, that can be used to perform monitoring campaign of air concentrations at a higher temporal and spatial resolution than the adopted standard methods (e.g. high-volume air samplers). The results could be employed for evaluation of the exposure, emission profiles and calibration of fate models.

Ecotoxicology of persistent organic pollutants in birds

Considering the explosive growth of the list of persistent organic pollutants (POPs), the scientific community is combatting increasing challenges to protect humans and wildlife from the potentially negative consequences of POPs. Herein, we characterize the main aspects and progress in the ecotoxicology of POPs in avian species since 2000. The majority of previous efforts has revealed the global occurrence of high levels of various POPs in birds. Laboratory research and epidemiological studies imply that POPs exert a broad-spectrum of side-effects on birds by interfering with their endocrine, immune and neural system, reproduction, and development, and growth. However, inconsistent results suggest that the potential effects of POP exposure on the physiological parameters in birds are multifactorial, involving a multitude of biological processes, species-specific differences, gender, age and types of compounds. Great progress has been achieved in identifying the species-specific sensitivity to dioxin-like compounds, which is attributed to different amino acid residues in the ligand-binding domain of the aryl hydrocarbon receptor. Besides the conventional concentration additivity, several studies have suggested that different classes of POPs possibly act synergistically or antagonistically based on their concentration. However, ecotoxicology information is still recorded in a scattered and inadequate manner, including lack of enough avian species, limited number of POPs investigated, and insufficient geographical representation, and thus our understanding of the effects of POPs on birds remains rudimentary, although mechanistic understanding of their mode of action is progressing. Particularly, research on what happens to wild bird populations and their ecosystems under POP stress is still unavailable. Thus, our aim is to predict and trace the effects POPs at different biological organization levels, especially from the molecular, cellular and individual levels to the population, community and ecosystem levels because of the limited and scattered information, as mentioned above.

Cross-species extrapolation of chemical sensitivity

Ecosystems are usually populated by many species. Each of these species carries the potential to show a different sensitivity towards all of the numerous chemical compounds that can be present in their environment. Since experimentally testing all possible species-chemical combinations is impossible, the ecological risk assessment of chemicals largely depends on cross-species extrapolation approaches. This review overviews currently existing cross-species extrapolation methodologies, and discusses i) how species sensitivity could be described, ii) which predictors might be useful for explaining differences in species sensitivity, and iii) which statistical considerations are important. We argue that risk assessment can benefit most from modelling approaches when sensitivity is described based on ecologically relevant and robust effects. Additionally, specific attention should be paid to heterogeneity of the training data (e.g. exposure duration, pH, temperature), since this strongly influences the reliability of the resulting models. Regarding which predictors are useful for explaining differences in species sensitivity, we review interspecies-correlation, relatedness-based, traits-based, and genomic-based extrapolation methods, describing the amount of mechanistic information the predictors contain, the amount of input data the models require, and the extent to which the different methods provide protection for ecological entities. We develop a conceptual framework, incorporating the strengths of each of the methods described. Finally, the discussion of statistical considerations reveals that regardless of the method used, statistically significant models can be found, although the usefulness, applicability, and understanding of these models varies considerably. We therefore recommend publication of scientific code along with scientific studies to simultaneously clarify modelling choices and enable elaboration on existing work. In general, this review specifies the data requirements of different cross-species extrapolation methods, aiming to make regulators and publishers more aware that access to raw- and meta-data needs to be improved to make future cross-species extrapolation efforts successful, enabling their integration into the regulatory environment.

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.

Community composition modifies direct and indirect effects of pesticides in freshwater food webs

For environmental risk assessment, the effects of pesticides on aquatic ecosystems are often assessed based on single species tests, disregarding the potential influence of community composition. We, therefore, studied the influence of changing the horizontal (the number of species within trophic levels) and vertical composition (number of trophic levels) on the ecological effects of the herbicide linuron and the insecticide chlorpyrifos, targeting producers and herbivores, respectively. We tested how adding, to a single primary producer, 4 selected competing producer species, 0-1-4 selected herbivore species, and one selected predator species resulting in 1, 2 and 3 trophic levels, changes the effects of the two pesticides. Linuron decreased producer biovolume less (17%) when the 4 producers were added, because insensitive producers compensated for the loss of sensitive producers. However, linuron decreased producer biovolume 42% and 32% more as we increased the number of herbivore species from 0 to 4 and as we increased trophic levels from 1 to 3, respectively. The indirect negative effect of linuron on herbivore biovolume was 11% and 15% lower when more producer and herbivores were added, respectively. Adding a predator increased this indirect negative effect by 22%. Chlorpyrifos decreased herbivore biovolume about 10% less when adding multiple herbivore or producer species. However, adding a predator magnified the direct negative impact on herbivores (13%). Increasing the number of producer, herbivore species and adding trophic levels increased the indirect positive impact on producer biovolume (between 10% and 35%). Our study shows that changing horizontal composition can both increase and decrease the effects of the selected pesticides, while changing vertical composition by adding number of trophic levels always increased these effects. Therefore, single species sensitivity will not always represent a worst case estimate of ecological effects. Protecting the most sensitive species may not ensure protection of ecosystems.

Periodic density as an endpoint of customized plankton community responses to petroleum hydrocarbons: A level of toxic effect should be matched with a suitable time scale

As an endpoint of community response to contaminants, average periodic density of populations (APDP) has been introduced to model species interactions in a community with 4 planktonic species. An ecological model for the community was developed by means of interspecific relationship including competition and predation to calculate the APDP. As a case study, we reported here the ecotoxicological effects of petroleum hydrocarbons (PHC) collected from Bohai oil field on densities of two algae, Platymonas subcordiformis and Isochrysis galbana, a rotifer, Brachionus plicatilis, and of a cladocera, Penilia avirostris, in single species and a microcosm experiment. Time scales expressing toxic effect increased with increasing levels of toxic effect from molecule to community. Remarkable periodic changes in densities were found during the tests in microcosm experiment, revealing a strong species reaction. The minimum time scale characterizing toxic effect at a community level should be the common cycle of population densities of the microcosm. In addition, the cycles of plankton densities shortened in general with increasing PHC, showing an evident toxic effect on the microcosm. Using APDP as the endpoint, a threshold concentration for the modeled microcosm was calculated to be 0.404 mg-PHC L^-1. The APDP was found to be more sensitive and reliable than the standing crops of populations as the endpoint. This indicated that the APDP, an endpoint at the community level, could be quantitatively related to the endpoints at the population level, and led to the quantitative concentration-toxic effect relationship at the community level.

Potential impact of chemical stress on freshwater invertebrates: A sensitivity assessment on continental and national scale based on distribution patterns, biological traits, and relatedness

Current chemical risk assessment approaches rely on a standard suite of test species to assess toxicity to environmental species. Assessment factors are used to extrapolate from single species to communities and ecosystem effects. This approach is pragmatic, but lacks resolution in biological and environmental parameters. Novel modelling approaches can help improve the biological resolution of assessments by using mechanistic information to identify priority species and priority regions that are potentially most impacted by chemical stressors. In this study we developed predictive sensitivity models by combining species-specific information on acute chemical sensitivity (LC50 and EC50), traits, and taxonomic relatedness. These models were applied at two spatial scales to reveal spatial differences in the sensitivity of species assemblages towards two chemical modes of action (MOA): narcosis and acetylcholinesterase (AChE) inhibition. We found that on a relative scale, 46% and 33% of European species were ranked as more sensitive towards narcosis and AChE inhibition, respectively. These more sensitive species were distributed with higher occurrences in the south and north-eastern regions, reflecting known continental patterns of endemic macroinvertebrate biodiversity. We found contradicting sensitivity patterns depending on the MOA for UK scenarios, with more species displaying relative sensitivity to narcotic MOA in north and north-western regions, and more species with relative sensitivity to AChE inhibition MOA in south and south-western regions. Overall, we identified hotspots of species sensitive to chemical stressors at two spatial scales, and discuss data gaps and crucial technological advances required for the successful application of the proposed methodology to invertebrate scenarios, which remain underrepresented in global conservation priorities.

A 3-Species Aquatic Community Model for Ecological Risk Assessment Using Basic Ecotoxicity Data

A simplified ecosystem model, the Aquatic Tritrophic Ecological Risk Assessment Model (A-TERAM), for the ecological risk assessment of chemicals is presented. The A-TERAM comprises a linear grazer food chain with 3 trophic levels (the algae-Daphnia-fish system). The model simulates the seasonal patterns of abundance at each level observed in the field, and it translates the direct toxic effects of chemicals on algae or Daphnia to implications for fish via ecological interactions; thus, the A-TERAM evaluates ecological risk in terms of the annual population growth rate of fish. The model also incorporates toxicokinetics for fish. The minimum input data required for the A-TERAM are basic ecotoxicity endpoints (algal growth inhibition median effect concentration [EC50] or no-observed effect concentration, Daphnia immobility EC50, and fish acute mortality median lethal concentration); however, additional ecotoxicity data (Daphnia reproduction test, fish early life test, and fish reproduction test) are also relevant for improving simulations. Comparisons made across 496 chemicals (255 nonagricultural chemicals and 241 agrochemicals) indicated that the A-TERAM, in comparison with the conventional predicted-effect concentration/predicted-no-effect concentration method, tended to evaluate higher risk to chemicals that are highly bioaccumulative and toxic to fish by 2 orders of magnitude at the largest but lower or comparable risk to chemicals that are toxic only to algae or Daphnia. Environ Toxicol Chem 2020;39:1086-1100. © 2020 SETAC.

Plants radically change the mobility of PCBs in soil: Role of different species and soil conditions

The mobility of Polychlorinated Biphenyls (PCBs) in soil cultivated with different plant species was evaluated by means of a column experiment to investigate the specific plant influence on PCB environmental fate and the potential for leaching. The soil was collected at a National Relevance Site for remediation located in Northern Italy (SIN Brescia-Caffaro) and underwent a rhizoremediation treatment for 18 months with different plant species (Festuca arundinacea, Cucurbita pepo ssp pepo and Medicago sativa). The same but unplanted soil was also considered as control for comparison. The columns were leached with tap water and PCB concentrations were measured in the leachate after 7 days of soil/water contact. Soil previously cultivated with different plant species exhibited statistically different behavior in terms of chemical leaching among the different fractions. Total PCB bulk concentrations ranged from 24 to 219 ng/L. Leachate samples were enriched in tetra- to hepta-PCBs. While PCB concentrations in the dissolved phases varied within a factor of 2 between controls and treatments, PCB associated to particulate organic carbon (POC) differed by more than one order of magnitude. More specifically, Medicago sativa enriched the soil with POC doubling PCB leaching with respect to the other plant species and the unplanted controls.

Micropollutants in Lake Como water in the context of circular economy: A snapshot of water cycle contamination in a changing pollution scenario

In this work we evaluated the contamination of the water cycle in Como Bay by measuring 38 selected pharmaceuticals in two main wastewater treatment plant in Switzerland and in Italy, two influents (River Breggia and Cosia), lake water (epilimnion and hypolimnion), as well as potable water. The collection of comparative information on the presence and environmental fate of these substances contributes to set specific environmental quality standard (EQS). The results presented show that the contamination of the lake reflects national health policies, which deeply influence the usage of chemicals. The outcomes of this study give an overall picture of contamination in the area, showing that concentrations of the measured compounds are generally low and under the commonly adopted ecotoxicological and toxicological thresholds. Only in a few cases did the contamination appear to be noteworthy, for some of the most persistent compounds (antibiotic macrolides, diclofenac, irbesartan, carbamazepine and dihydrocarbamazepine, bezafibrate, furosemide and hydrochlorothiazide). Some concern can be also be raised for the presence of antibiotics (clarithromycin) in drinking water, although at very low levels, due to the problem of antibiotic resistance.

Modelling peak exposure of pesticides in terrestrial and aquatic ecosystems: importance of dissolved organic carbon and vertical particle movement in soil

In the present work, an existing vegetation/air/litter/soil model (SoilPlusVeg) was modified to improve organic chemical fate description in terrestrial/aquatic ecosystems accounting for horizontal and vertical particulate organic carbon (POC) transport in soil. The model was applied to simulate the fate of three pesticides (terbuthylazine, chlorpyrifos and etofenprox), characterized by increasing hydrophobicity (log K_OW from about 3 to 7), in the soil compartment and more specifically, their movement towards surface and groundwater through infiltration and runoff processes. The aim was to evaluate the role of dissolved organic carbon (DOC) and POC in the soil in influencing the peak exposure of pesticides in terrestrial/aquatic ecosystems. Simulation results showed that while terbuthylazine and chlorpyrifos dominated the free water phase (C_W-FREE) of soil, etofenprox was mainly present in soil porewater as POC associated chemical. This resulted in an increase of this highly hydrophobic chemical movement towards groundwater and surface water, up to a factor of 40. The present work highlighted the importance of DOC and POC as an enhancer of mobility in the water of poor or very little mobile chemicals. Further studies are necessary to evaluate the bioavailability change with time and parameterize this process in multimedia fate models.

Towards a general framework for the assessment of interactive effects of multiple stressors on aquatic ecosystems: Results from the Making Aquatic Ecosystems Great Again (MAEGA) workshop

A workshop was held in Wageningen, The Netherlands, in September 2017 to collate data and literature on three aquatic ecosystem types (agricultural drainage ditches, urban floodplains, and urban estuaries), and develop a general framework for the assessment of multiple stressors on the structure and functioning of these systems. An assessment framework considering multiple stressors is crucial for our understanding of ecosystem responses within a multiply stressed environment, and to inform appropriate environmental management strategies. The framework consists of two components: (i) problem identification and (ii) impact assessment. Both assessments together proceed through the following steps: 1) ecosystem selection; 2) identification of stressors and quantification of their intensity; 3) identification of receptors or sensitive groups for each stressor; 4) identification of stressor-response relationships and their potential interactions; 5) construction of an ecological model that includes relevant functional groups and endpoints; 6) prediction of impacts of multiple stressors, 7) confirmation of these predictions with experimental and monitoring data, and 8) potential adjustment of the ecological model. Steps 7 and 8 allow the assessment to be adaptive and can be repeated until a satisfactory match between model predictions and experimental and monitoring data has been obtained. This paper is the preface of the MAEGA (Making Aquatic Ecosystems Great Again) special section that includes three associated papers which are also published in this volume, which present applications of the framework for each of the three aquatic systems.

Modeling the Sensitivity of Aquatic Macroinvertebrates to Chemicals Using Traits

In this study, a trait-based macroinvertebrate sensitivity modeling tool is presented that provides two main outcomes: (1) it constructs a macroinvertebrate sensitivity ranking and, subsequently, a predictive trait model for each one of a diverse set of predefined Modes of Action (MOAs) and (2) it reveals data gaps and restrictions, helping with the direction of future research. Besides revealing taxonomic patterns of species sensitivity, we find that there was not one genus, family, or class which was most sensitive to all MOAs and that common test taxa were often not the most sensitive at all. Traits like life cycle duration and feeding mode were identified as important in explaining species sensitivity. For 71% of the species, no or incomplete trait data were available, making the lack of trait data the main obstacle in model construction. Research focus should therefore be on completing trait databases and enhancing them with finer morphological traits, focusing on the toxicodynamics of the chemical (e.g., target site distribution). Further improved sensitivity models can help with the creation of ecological scenarios by predicting the sensitivity of untested species. Through this development, our approach can help reduce animal testing and contribute toward a new predictive ecotoxicology framework.

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.

The combined and interactive effects of zinc, temperature, and phosphorus on the structure and functioning of a freshwater community

Ecotoxicological studies mainly consist of single-species experiments evaluating the effects of a single stressor. However, under natural conditions aquatic communities are exposed to a mixture of stressors. The present study aimed to identify how the toxicity of zinc (Zn) is affected by increased temperature and increased phosphorus (P) supply and how these interactions vary among species, functional groups, and community structure and function. Aquatic microcosms were subjected to 3 Zn concentrations (background, no Zn added, and 75 and 300 μg Zn/L), 2 temperatures (16-19 and 21-24 °C), and 2 different P additions (low, 0.02, and high, 0.4 mg P L^-1  wk^-1 ) for 5 wk using a full factorial design. During the study, consistent interactions between Zn and temperature were only rarely found at the species level (4%), but were frequently found at the functional group level (36%), for community structure (100%) and for community function (100%; such as dissolved organic carbon concentrations and total chlorophyll). The majority of the Zn × temperature interactions were observed at 300 μg Zn/L and generally indicated a smaller effect of Zn at higher temperature. Furthermore, no clear indication was found that high P addition by itself significantly affected the overall effects of Zn on the community at any level of organization. Interestingly, though, 90% of all the Zn × temperature interactions observed at the species, group, and community composition level were found under high P addition. Collectively, the results of our study with the model chemical Zn suggest that temperature and phosphorus loading to freshwater systems should be accounted for in risk assessment, because these factors may modify the effects of chemicals on the structure and functioning of aquatic communities, especially at higher levels of biological organization. Environ Toxicol Chem 2018;37:2413-2427. © 2018 SETAC.

Combined effects of interspecies interaction, temperature, and zinc on Daphnia longispina population dynamics

Under natural conditions, organisms can experience a variety of abiotic (e.g., temperature, pH) and biotic (e.g., species interactions) factors, which can interact with toxicant effects. By ignoring species interactions conventional ecotoxicological studies (i.e., single-species tests) oversimplify the actual field situation. We investigated whether temperature and interspecific competition affected the effects of zinc (Zn) on a Daphnia longispina population. The D. longispina populations were exposed in a full factorial design to 3 different Zn treatments (background, 29, and 110 μg Zn/L), 2 different temperature regimes (cold, 17-18 °C; warm, 21-22 °C), and 2 interspecific competition levels (no interspecific Brachionus competition = no Brachionus calyciflorus added; interspecific Brachionus competition = B. calyciflorus added). Interspecific Brachionus competition and temperature by itself had a limited effect on the Daphnia abundance but significantly interacted with the highest Zn concentration. Without Brachionus competition the D. longispina juvenile and adult abundances under warm conditions were up to 5.5 and 21 times lower, respectively, in the high Zn treatment in comparison with the Zn control, whereas under cold conditions no significant Zn effect was observed. However, with Brachionus competition the highest Zn treatment was on average 2.2 times less toxic to the D. longispina juvenile abundance at higher temperatures. Under cold conditions the highest Zn treatment affected the juvenile abundance sooner and up to 9 times more negatively when simultaneously faced with Brachionus competition. It is possible that the competition for food reduced the amount of energy that could be used by D. longispina for reproduction, and the metabolic costs increased as a result of Zn stress. The present study clearly illustrated the influence of temperature and competition on the effects of a chemical stressor. Thus, not considering such factors in ecological risk assessment may underestimate or overestimate risks in aquatic ecosystems when extrapolating data from standard single-species tests to the field. Environ Toxicol Chem 2018;37:1668-1678. © 2018 SETAC.

Predicting pesticide fate in small cultivated mountain watersheds using the DynAPlus model: Toward improved assessment of peak exposure

The use of plant protection products (PPPs) in agricultural areas implies potential chemical loadings to surface waters, which can pose a risk to aquatic ecosystems and human health. Due to the spatio-temporal variability of PPP applications and of the processes regulating their transport to surface waters, aquatic organisms are typically exposed to pulses of contaminants. In small mountain watersheds, where runoff fluxes are more rapid due to the steep slopes, such exposure peaks are particularly likely to occur. In this work, a spatially explicit, dynamic model for predicting pesticide exposure in surface waters of cultivated mountain basins (DynAPlus) has been developed. The model has been applied to a small mountain watershed (133km²) located in the Italian Eastern Alps and characterized by intensive agriculture (apple orchards) around the main river and its tributaries. DynAPlus performance was evaluated for chlorpyrifos through experimental monitoring, using samples collected during the 2011 and 2012 productive seasons. The comparison between predictions and measurements resulted in a good agreement (R²=0.49, efficiency factor 0.60), although a more accurate spatial information in the input scenario (e.g., field-specific applications, rainfall amount, soil properties) would dramatically improve model performance. A set of illustrative simulations performed for three PPPs highlighted the potential role of DynAPlus in improving exposure predictions for ecological risk assessment and pesticide management practices (e.g., for active ingredient and application rate selection), as well as for planning efficient monitoring campaigns and/or interpreting monitoring data. However, some model improvements (e.g., solid erosion and transport) and a more thorough model validation are desirable to enlarge the applicability domain.

Assessing recovery of stream insects from pesticides using a two-patch metapopulation model

Pesticides can exert lethal and sublethal effects on streams organisms. Field studies have shown that non-polluted upstream patches promote population recovery from such effects. Nevertheless, the dynamics and potential carryover effects on the upstream patch are largely unknown. We used a metapopulation model with 2 patches to simulate lethal and sublethal effects on the downstream population of an insect with one generation per year, which was structured into early and late instars aquatic life stages, and an adult terrestrial life stage. We examined the implications for the recovery time of a range of scenarios covering different pesticide effect combinations, migration and exposure types. We found that recovery time responded most strongly to a reduction in reproduction in terms of the reduction of the intrinsic growth rate during the early instar aquatic life stage. For 60 of 96 scenarios with pesticide exposure in consecutive years, no recovery occurred within one year if the intrinsic growth rate was reduced by 50% or more. Without migration between patches (32 scenarios), the polluted downstream population went extinct in 5 of these scenarios. Migration lead overall to slightly faster recovery, albeit this was scenario dependent, but also to a carryover of the pesticide effect from the polluted downstream to the non-polluted upstream patch (up to 25% reduction in the minimal population size). A sensitivity analysis revealed that recovery time was most sensitive to the parameters length of the intrinsic growth phase during early instar aquatic life stage and to migration mortality of the late instar aquatic larvae, and least sensitive to the adult emigration rate and timing of pesticide application. Our study highlights the important role of sublethal effects for population responses to pesticides and that migration buffers against effects, but also carries effects over to non-polluted patches.

Postregistration monitoring of pesticides is urgently required to protect ecosystems

Current admission policies for pesticides follow a controlled experimental tiered risk assessment approach, giving results that are difficult to extrapolate to a real-world situation. Later analyses of compounds such as DDT and neonicotinoid pesticides clearly show that the actual chemical impacts frequently affect many more components of an ecosystem than a priori suggested by risk assessment. Therefore, to manage the actual risks for ecosystems imposed by manufactured compounds, it is proposed that current admission policies for chemicals be enriched by using postregistration monitoring. Such monitoring is essential to identify unexpected direct and indirect impacts on organisms by accounting for multiple propagation routes and exposures. Implementation of postregistration monitoring could build on existing monitoring networks. This approach would tackle the current policy impasse of compartment-based regulations versus exposure-based regulations, and, more importantly, would provide a safety lock for risk assessment across compartments and more likely ensure the protection of our natural environment. Environ Toxicol Chem 2017;36:860-865. © 2017 SETAC.

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.

European environmental scenarios of chemical bioavailability in freshwater systems

In exposure prediction for environmental risk assessment, the transition to more dynamic and realistic modelling approaches and scenarios has been recently identified as a major challenge, since it would allow a more accurate prediction of bioavailable concentrations and their variations in space and time. In this work, an improved version of the multimedia model ChimERA fate, including a phytoplankton compartment and equations to calculate phytoplankton, detritus and dissolved organic matter variations in time, was developed. The model was parameterized to simulate five dynamic scenarios for shallow meso-eutrophic water bodies based on a latitudinal gradient (in Europe); such scenarios include seasonal profiles of water temperature, phytoplankton biomass, detritus, and dissolved organic matter. Model runs were performed for each scenario for 8 hydrophobic chemicals (PCB congeners), with the aim of investigating the influence of scenario characteristics and compound properties on bioavailable concentrations. The key processes were adsorption/uptake by phytoplankton and deposition to sediment of detritus-bound chemicals. The northern scenarios ("Scandinavia" and "UK") showed the highest bioavailable concentrations, with annual maximum/minimum concentration up to 25; in contrast, for example, maximum concentrations in the "Mediterranean" scenario were lower by a factor of 2 to 9 with respect to the northern ones (depending on chemical hydrophobicity), due to the generally higher biomass and carbon levels, and showed only limited seasonal variability (up to a factor of 4). These results highlight the importance of including biomass and organic carbon dynamics in both modelling approaches and scenarios for the evaluation of exposure concentrations in aquatic environments.

Towards the review of the European Union Water Framework Directive: Recommendations for more efficient assessment and management of chemical contamination in European surface water resources

Water is a vital resource for natural ecosystems and human life, and assuring a high quality of water and protecting it from chemical contamination is a major societal goal in the European Union. The Water Framework Directive (WFD) and its daughter directives are the major body of legislation for the protection and sustainable use of European freshwater resources. The practical implementation of the WFD with regard to chemical pollution has faced some challenges. In support of the upcoming WFD review in 2019 the research project SOLUTIONS and the European monitoring network NORMAN has analyzed these challenges, evaluated the state-of-the-art of the science and suggested possible solutions. We give 10 recommendations to improve monitoring and to strengthen comprehensive prioritization, to foster consistent assessment and to support solution-oriented management of surface waters. The integration of effect-based tools, the application of passive sampling for bioaccumulative chemicals and an integrated strategy for prioritization of contaminants, accounting for knowledge gaps, are seen as important approaches to advance monitoring. Including all relevant chemical contaminants in more holistic "chemical status" assessment, using effect-based trigger values to address priority mixtures of chemicals, to better consider historical burdens accumulated in sediments and to use models to fill data gaps are recommended for a consistent assessment of contamination. Solution-oriented management should apply a tiered approach in investigative monitoring to identify toxicity drivers, strengthen consistent legislative frameworks and apply solutions-oriented approaches that explore risk reduction scenarios before and along with risk assessment.

Prospective Environmental Risk Assessment for Sediment-Bound Organic Chemicals: A Proposal for Tiered Effect Assessment

A broadly accepted framework for prospective environmental risk assessment (ERA) of sediment-bound organic chemicals is currently lacking. Such a framework requires clear protection goals, evidence-based concepts that link exposure to effects and a transparent tiered-effect assessment. In this paper, we provide a tiered prospective sediment ERA procedure for organic chemicals in sediment, with a focus on the applicable European regulations and the underlying data requirements. Using the ecosystem services concept, we derived specific protection goals for ecosystem service providing units: microorganisms, benthic algae, sediment-rooted macrophytes, benthic invertebrates and benthic vertebrates. Triggers for sediment toxicity testing are discussed.We recommend a tiered approach (Tier 0 through Tier 3). Tier-0 is a cost-effective screening based on chronic water-exposure toxicity data for pelagic species and equilibrium partitioning. Tier-1 is based on spiked sediment laboratory toxicity tests with standard benthic test species and standardised test methods. If comparable chronic toxicity data for both standard and additional benthic test species are available, the Species Sensitivity Distribution (SSD) approach is a more viable Tier-2 option than the geometric mean approach. This paper includes criteria for accepting results of sediment-spiked single species toxicity tests in prospective ERA, and for the application of the SSD approach. We propose micro/mesocosm experiments with spiked sediment, to study colonisation success by benthic organisms, as a Tier-3 option. Ecological effect models can be used to supplement the experimental tiers. A strategy for unifying information from various tiers by experimental work and exposure-and effect modelling is provided.

Evolutionary toxicology: Meta-analysis of evolutionary events in response to chemical stressors

The regulatory decision-making process regarding chemical safety is most often informed by evidence based on ecotoxicity tests that consider growth, reproduction and survival as end-points, which can be quantitatively linked to short-term population outcomes. Changes in these end-points resulting from chemical exposure can cause alterations in micro-evolutionary forces (mutation, drift, selection and gene flow) that control the genetic composition of populations. With multi-generation exposures, anthropogenic contamination can lead to a population with an altered genetic composition, which may respond differently to future stressors. These evolutionary changes are rarely discussed in regulatory or risk assessment frameworks, but the growing body of literature that documents their existence suggests that these important population-level impacts should be considered. In this meta-analysis we have compared existing contamination levels of polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons (PAHs) that have been documented to be associated with evolutionary changes in resident aquatic organisms to regulatory benchmarks for these contaminants. The original intent of this project was to perform a meta-analysis on evolutionary events associated with PCB and PAH contamination. However, this effort was hindered by a lack of consistency in congener selection for "total" PCB or PAH measurements. We expanded this manuscript to include a discussion of methods used to determine PCB and PAH total contamination in addition to comparing regulatory guidelines and contamination that has caused evolutionary effects. Micro-evolutionary responses often lead populations onto unique and unpredictable trajectories. Therefore, to better understand the risk of population-wide alterations occurring, we need to improve comparisons of chemical contamination between affected locations. In this manuscript we offer several possibilities to unify chemical comparisons for PCBs and PAHs that would improve comparability among evolutionary toxicology investigations, and with regulatory guidelines. In addition, we identify studies documenting evolutionary change in the presence of PCB and PAH contamination levels below applicable regulatory benchmarks.

Relative influence of chemical and non-chemical stressors on invertebrate communities: a case study in the Danube River

A key challenge for the ecological risk assessment of chemicals has been to evaluate the relative contribution of chemical pollution to the variability observed in biological communities, as well as to identify multiple stressor groups. In this study we evaluated the toxic pressure exerted by >200 contaminants to benthic macroinvertebrates in the Danube River using the Toxic Unit approach. Furthermore, we evaluated correlations between several stressors (chemical and non-chemical) and biological indices commonly used for the ecological status assessment of aquatic ecosystems. We also performed several variation partitioning analyses to evaluate the relative contribution of contaminants and other abiotic parameters (i.e. habitat characteristics, hydromorphological alterations, water quality parameters) to the structural and biological trait variation of the invertebrate community. The results of this study show that most biological indices significantly correlate to parameters related to habitat and physico-chemical conditions, but showed limited correlation with the calculated toxic pressure. The calculated toxic pressure, however, showed little variation between sampling sites, which complicates the identification of pollution-induced effects. The results of this study show that the variation in the structure and trait composition of the invertebrate community are mainly explained by habitat and water quality parameters, whereas hydromorphological alterations play a less important role. Among the water quality parameters, physico-chemical parameters such as suspended solids, nutrients or dissolved oxygen explained a larger part of the variation in the invertebrate community as compared to metals or organic contaminants. Significant correlations exist between some physico-chemical measurements (e.g. nutrients) and some chemical classes (i.e. pharmaceuticals, chemicals related to human presence) which constitute important multiple stressor groups. This study demonstrates that, in large rivers like the Danube, the variation in the invertebrate community seems to be more related to varying habitat and physico-chemical conditions than to chemical pollution.

Risk assessment of pesticides used in rice-prawn concurrent systems in Bangladesh

The objectives of the current study were to determine the occupational health hazards posed by the application of pesticides in rice-prawn concurrent systems of south-west Bangladesh and to assess their potential risks for the aquatic ecosystems that support the culture of freshwater prawns (Macrobrachium rosenbergii). Information on pesticide use in rice-prawn farming was collected through structured interviews with 38 farm owners held between January and May of 2012. The risks of the pesticide use to human health were assessed through structured interviews. The TOXSWA model was used to calculate pesticide exposure (peak and time-weighted average concentrations) in surface waters of rice-prawn systems for different spray drift scenarios and a simple first tier risk assessment based on threshold concentrations derived from single species toxicity tests were used to assess the ecological risk in the form of risk quotients. The PERPEST model was used to refine the ecological risks when the first tier assessment indicated a possible risk. Eleven synthetic insecticides and one fungicide (sulphur) were recorded as part of this investigation. The most commonly reported pesticide was sulphur (used by 29% of the interviewed farmers), followed by thiamethoxam, chlorantraniliprole, and phenthoate (21%). A large portion of the interviewed farmers described negative health symptoms after pesticide applications, including vomiting (51%), headache (18%) and eye irritation (12%). The results of the first tier risk assessment indicated that chlorpyrifos, cypermethrin, alpha-cypermethrin, and malathion may pose a high to moderate acute and chronic risks for invertebrates and fish in all evaluated spray drift scenarios. The higher tier assessment using the PERPEST model confirmed the high risk of cypermethrin, alpha-cypermethrin, and chlorpyrifos for insects and macro- and micro-crustaceans thus indicating that these pesticides may have severe adverse consequences for the prawn production yields.

Crude oil affecting the biomass of the marine copepod Calanus finmarchicus: Comparing a simple and complex population model

In the current study differences were evaluated between a complex 3D multistage population model (SINMOD) and a simpler consumer-resource population model for estimating the effects of crude oil on the marine copepod Calanus finmarchicus. The SINTEF OSCAR model was used to simulate hypothetical oil spills in the Lofoten area in 1995, 1997, and 2001. Both population models simulated a negligible effect of crude oil on the Calanus' biomass when assuming low species sensitivity. The simple model estimated a larger effect on the biomass (up to a 100% decline) compared to the complex model (maximum decline of 60-80%) at high species sensitivity to crude oil. These differences may be related to the inclusion of copepod advection in the complex model. Our study showed that if little data is available to parameterize a model, or if computational resources are scarce, the simple model could be used for risk screening. Nevertheless, the possibility of including a dilution factor for time-varying biomass should be examined to improve the estimations of the simple model. The complex model should be used for a more in depth risk analysis, as it includes physical processes such as the drift of organisms and differentiation between developmental stages.

Developing ecological scenarios for the prospective aquatic risk assessment of pesticides

The prospective aquatic environmental risk assessment (ERA) of pesticides is generally based on the comparison of predicted environmental concentrations in edge-of-field surface waters with regulatory acceptable concentrations derived from laboratory and/or model ecosystem experiments with aquatic organisms. New improvements in mechanistic effect modeling have allowed a better characterization of the ecological risks of pesticides through the incorporation of biological trait information and landscape parameters to assess individual, population and/or community-level effects and recovery. Similarly to exposure models, ecological models require scenarios that describe the environmental context in which they are applied. In this article, we propose a conceptual framework for the development of ecological scenarios that, when merged with exposure scenarios, will constitute environmental scenarios for prospective aquatic ERA. These "unified" environmental scenarios are defined as the combination of the biotic and abiotic parameters that are required to characterize exposure, (direct and indirect) effects, and recovery of aquatic nontarget species under realistic worst-case conditions. Ideally, environmental scenarios aim to avoid a potential mismatch between the parameter values and the spatial-temporal scales currently used in aquatic exposure and effect modeling. This requires a deeper understanding of the ecological entities we intend to protect, which can be preliminarily addressed by the formulation of ecological scenarios. In this article we present a methodological approach for the development of ecological scenarios and illustrate this approach by a case-study for Dutch agricultural ditches and the example focal species Sialis lutaria. Finally, we discuss the applicability of ecological scenarios in ERA and propose research needs and recommendations for their development and integration with exposure scenarios. Integr Environ Assess Manag 2016;12:510-521. © 2015 SETAC.

Modeling Macroinvertebrate Community Dynamics in Stream Mesocosms Contaminated with a Pesticide

Modeling community dynamics of aquatic invertebrates is an important but challenging task, in particular in ecotoxicological risk assessment. Systematic parameter estimation and rigorous assessment of model uncertainty are often lacking in such applications. We applied the mechanistic food web model Streambugs to investigate the temporal development of the macroinvertebrate community in an ecotoxicological mesocosm experiment with pulsed contaminations with the insecticide thiacloprid. We used Bayesian inference to estimate parameters and their uncertainty. Approx. 85% of all experimental observations lie within the 90% uncertainty intervals indicating reasonably good fits of the calibrated model. However, a validation with independent data was not possible due to lacking data. Investigation of vital rates and limiting factors in the model yielded insights into recovery dynamics. Inclusion of the emergence process and sub-lethal effects turned out to be potentially relevant model extensions. Measurements of food source dynamics, individual body size (classes), and additional knowledge on sub-lethal effects would support more accurate modeling. This application of a process-based, ecotoxicological community model with uncertainty assessment by Bayesian inference increased our process understanding of toxicant effects in macroinvertebrate communities and helped identifying potential improvements in model structure and experimental design.

Importance of environmental and biomass dynamics in predicting chemical exposure in ecological risk assessment

In ecological risk assessment, exposure is generally modelled assuming static conditions, herewith neglecting the potential role of emission, environmental and biomass dynamics in affecting bioavailable concentrations. In order to investigate the influence of such dynamics on predicted bioavailable concentrations, the spatially-resolved dynamic model "ChimERA fate" was developed, incorporating macrophyte and particulate/dissolved organic carbon (POC/DOC) dynamics into a water-sediment system. An evaluation against three case studies revealed a satisfying model performance. Illustrative simulations then highlighted the potential spatio-temporal variability of bioavailable concentrations after a pulsed emission of four chemicals in a system composed of a pond connected to its inflow and outflow streams. Changes in macrophyte biomass and POC/DOC levels caused exposure variations which were up to a factor of 4.5 in time and even more significant (several orders of magnitude) in space, especially for highly hydrophobic chemicals. ChimERA fate thus revealed to be a useful tool to investigate such variations and to identify those environmental and ecological conditions in which risk is expected to be highest.