Assessment of road run-off and domestic wastewater contribution to microplastic pollution in a densely populated area (Flanders, Belgium).

Plastics are omnipresent in our daily life. Unfortunately, the produced plastics will partly end up in the environment including aquatic ecosystems. People often refer to littering or illegal waste dumping as sources of plastic emission to the environment. However, daily-life sources could also, unknowingly, contribute considerably to the total microplastic pollution in the ecosystem. Hence, there is an urgent need to study these potential sources. In this research, two common sources, i.e. domestic wastewater and road run-off from tire and road wear particles, were studied in detail to quantify the relative contribution of both domestic sources towards microplastic pollution in freshwater ecosystems in Flanders, Belgium. This assessment shows that every person (in studied area) emits on average 1145 microplastics (25-1000 µm) daily through domestic wastewater, resulting in a yearly discharge of 418,000 microplastic particles per person. The road run-off samples contained between 0.02 and 9.2 mg tire wear particles per litre per day, which corresponds to an emission of 10.8 mg tire wear particles per driven vehicle km. The gross and net emissions of both above mentioned microplastic sources were extrapolated to the whole Flanders region using an emission model. From the yearly gross microplastic pollution in the domestic wastewater, 623 kg (20%) will be discharged in the freshwater. The highest losses originated from the households that have a private drain or are not (yet) connected to an active wastewater treatment plant. In Flanders, the yearly net microplastic emission into the aquatic environment of tire wear particles is estimated to be 246 tonnes (38%), mainly from the direct run-off from the road surface. Based on the results, specific mitigation measures can be installed to reduce the emission of microplastics towards the freshwater ecosystem. Other sources should be quantified in a similar way for a more holistic strategy to counteract plastic pollution.

Risk assessment of microplastics in the ocean: Modelling approach and first conclusions.

We performed an environmental risk assessment for microplastics (<5 mm) in the marine environment by estimating the order of magnitude of the past, present and future concentrations based on global plastic production data. In 2100, from 9.6 to 48.8 particles m-3 are predicted to float around in the ocean, which is a 50-fold increase compared to the present-day concentrations. From a meta-analysis with effect data available in literature, we derived a safe concentration of 6650 buoyant particles m-3 below which adverse effects are not likely to occur. Our risk assessment (excluding the potential role of microplastics as chemical vectors) suggests that on average, no direct effects of free-floating microplastics in the marine environment are to be expected up to the year 2100. Yet, even today, the safe concentration can be exceeded in sites that are heavily polluted with buoyant microplastics. In the marine benthic compartment between 32 and 144 particles kg-1 dry sediment are predicted to be present in the beach deposition zone. Despite the scarcity of effect data, we expect adverse ecological effects along the coast as of the second half of the 21st century. From then ambient concentrations will start to outrange the safe concentration of sedimented microplastics (i.e. 540 particles kg-1 sediment). Additional ecotoxicological research in which marine species are chronically exposed to realistic environmental microplastic concentration series are urgently needed to verify our findings.

A framework for ecological risk assessment of metal mixtures in aquatic systems.

Although metal mixture toxicity has been studied relatively intensely, there is no general consensus yet on how to incorporate metal mixture toxicity into aquatic risk assessment. We combined existing data on chronic metal mixture toxicity at the species level with species sensitivity distribution (SSD)-based in silico metal mixture risk predictions at the community level for mixtures of Ni, Zn, Cu, Cd, and Pb, to develop a tiered risk assessment scheme for metal mixtures in freshwater. Generally, independent action (IA) predicts chronic metal mixture toxicity at the species level most accurately, whereas concentration addition (CA) is the most conservative model. Mixture effects are noninteractive in 69% (IA) and 44% (CA) and antagonistic in 15% (IA) and 51% (CA) of the experiments, whereas synergisms are only observed in 15% (IA) and 5% (CA) of the experiments. At low effect sizes (∼ 10% mixture effect), CA overestimates metal mixture toxicity at the species level by 1.2-fold (i.e., the mixture interaction factor [MIF]; median). Species, metal presence, or number of metals does not significantly affect the MIF. To predict metal mixture risk at the community level, bioavailability-normalization procedures were combined with CA or IA using SSD techniques in 4 different methods, which were compared using environmental monitoring data of a European river basin (the Dommel, The Netherlands). We found that the simplest method, in which CA is directly applied to the SSD (CASSD ), is also the most conservative method. The CASSD has median margins of safety (MoS) of 1.1 and 1.2 respectively for binary mixtures compared with the theoretically more consistent methods of applying CA or IA to the dose-response curve of each species individually prior to estimating the fraction of affected species (CADRC or IADRC ). The MoS increases linearly with an increasing number of metals, up to 1.4 and 1.7 for quinary mixtures (median) compared with CADRC and IADRC , respectively. When our methods were applied to a geochemical baseline database (Forum of European Geological Surveys [FOREGS]), we found that CASSD yielded a considerable number of mixture risk predictions, even when metals were at background levels (8% of the water samples). In contrast, metal mixture risks predicted with the theoretically more consistent methods (e.g., IADRC ) were very limited under natural background metal concentrations (<1% of the water samples). Based on the combined evidence of chronic mixture toxicity predictions at the species level and evidence of in silico risk predictions at the community level, a tiered risk assessment scheme for evaluating metal mixture risks is presented, with CASSD functioning as a first, simple conservative tier. The more complex, but theoretically more consistent and most accurate method, IADRC , can be used in higher tier assessments. Alternatively, the conservatism of CASSD can be accounted for deterministically by incorporating the MoS and MIF in the scheme. Finally, specific guidance is also given related to specific issues, such as how to deal with nondetect data and complex mixtures that include so-called data-poor metals. Environ Toxicol Chem 2018;37:623-642. © 2017 SETAC.

Mixture toxicity in the marine environment: Model development and evidence for synergism at environmental concentrations.

Little is known about the effect of metal mixtures on marine organisms, especially after exposure to environmentally realistic concentrations. This information is, however, required to evaluate the need to include mixtures in future environmental risk assessment procedures. We assessed the effect of copper (Cu)-Nickel (Ni) binary mixtures on Mytilus edulis larval development using a full factorial design that included environmentally relevant metal concentrations and ratios. The reproducibility of the results was assessed by repeating this experiment 5 times. The observed mixture effects were compared with the effects predicted with the concentration addition model. Deviations from the concentration addition model were estimated using a Markov chain Monte-Carlo algorithm. This enabled the accurate estimation of the deviations and their uncertainty. The results demonstrated reproducibly that the type of interaction-synergism or antagonism-mainly depended on the Ni concentration. Antagonism was observed at high Ni concentrations, whereas synergism occurred at Ni concentrations as low as 4.9 µg Ni/L. This low (and realistic) Ni concentration was 1% of the median effective concentration (EC50) of Ni or 57% of the Ni predicted-no-effect concentration (PNEC) in the European Union environmental risk assessment. It is concluded that results from mixture studies should not be extrapolated to concentrations or ratios other than those investigated and that significant mixture interactions can occur at environmentally realistic concentrations. This should be accounted for in (marine) environmental risk assessment of metals. Environ Toxicol Chem 2017;36:3471-3479. © 2017 SETAC.

Comparison of four methods for bioavailability-based risk assessment of mixtures of Cu, Zn, and Ni in freshwater.

Although chemical risk assessment is still mainly conducted on a substance-by-substance basis, organisms in the environment are typically exposed to mixtures of substances. Risk assessment procedures should therefore be adapted to fit these situations. Four mixture risk assessment methodologies were compared for risk estimations of mixtures of copper (Cu), zinc (Zn), and nickel (Ni). The results showed that use of the log-normal species sensitivity distribution (SSD) instead of the best-fit distribution and sampling species sensitivities independently for each metal instead of using interspecies correlations in metal sensitivity had little impact on risk estimates. Across 4 different monitoring datasets, between 0% and 52% of the target water samples were estimated to be at risk, but only between 0% and 15% of the target water samples were at risk because of the mixture of metals and not any single metal individually. When a natural baseline database was examined, it was estimated that 10% of the target water samples were at risk because of single metals or their mixtures; when the most conservative method was used (concentration addition [CA] applied directly to the SSD, i.e., CASSD ). However, the issue of metal mixture risk at geochemical baseline concentrations became relatively small (2% of target water samples) when a theoretically more correct method was used (CA applied to individual dose response curves, i.e., CADRC ). Finally, across the 4 monitoring datasets, the following order of conservatism for the 4 methods was shown (from most to least conservative, with ranges of median margin of safety [MoS] relative to CASSD ): CASSD > CADRC (MoS = 1.17-1.25) > IADRC (independent action (IA) applied to individual dose-response curves; MoS = 1.38-1.60) > IASSD (MoS = 1.48-1.72). Therefore, it is suggested that these 4 methods can be used in a general tiered scheme for the risk assessment of metal mixtures in a regulatory context. In this scheme, the CASSD method could serve as a first (conservative) tier to identify situations with likely no potential risk at all, regardless of the method used (the sum toxic unit expressed relative to the 5% hazardous concentration [SumTUHC5 ] < 1) and the IASSD method to identify situations of potential risk, also regardless of the method used (the multisubstance potentially affected fraction of species using the IASSD method [msPAFIA,SSD ] > 0.05). The CADRC and IADRC methods could be used for site-specific assessment for situations that fall in between (SumTUHC5 > 1 and msPAFIA,SSD < 0.05). Environ Toxicol Chem 2017;36:2123-2138. © 2017 SETAC.

The ChimERA project: coupling mechanistic exposure and effect models into an integrated platform for ecological risk assessment.

Current techniques for the ecological risk assessment of chemical substances are often criticised for their lack of environmental realism, ecological relevance and methodological accuracy. ChimERA is a 3-year project (2013-2016), funded by Cefic's Long Range Initiative (LRI) that aims to address some of these concerns by developing and testing mechanistic fate and effect models, and coupling of these models into one integrated platform for risk assessment. This paper discusses the backdrop against which this project was initiated and lists its objectives and planned methodology.

An approach to assess the regulatory relevance of microevolutionary effects in ecological risk assessment of chemicals: a case study with cadmium.

The authors suggest an approach to assess the regulatory relevance of microevolutionary effects of chemicals based on a comparison of concentrations at which microevolutionary effects have been reported in the literature and conventionally derived ecotoxicological threshold concentrations. The authors found reports of microevolutionary effects of cadmium in freshwater organisms at hardness-normalized concentrations between 0.5 µg Cd L(-1) and 6290 µg Cd L(-1) (normalized to a hardness of 50 mg CaCO3 L(-1)). These concentrations were at least 1.5 times higher than the hardness-normalized hazardous concentration for 5% of the organisms of 0.34 µg Cd L(-1). This suggests that there is no immediate need to consider microevolutionary effects of Cd in environmental risk assessments of freshwater environments. However, some other aspects should be kept in mind as well. First, microevolutionary effects have so far only been investigated at few, relatively high concentrations of Cd and not encompassing the 5% hazardous concentration. Second, different types of microevolutionary effects or investigated ecotoxicological end points may influence the conclusions of the suggested comparative approach. Finally, factors influencing the bioavailability of Cd were not commonly reported in the literature, which made normalization of concentrations at which evolutionary effects occurred impossible and affected the number of studies that could be evaluated in the suggested approach.

Assessment of marine debris on the Belgian Continental Shelf.

A comprehensive assessment of marine litter in three environmental compartments of Belgian coastal waters was performed. Abundance, weight and composition of marine debris, including microplastics, was assessed by performing beach, sea surface and seafloor monitoring campaigns during two consecutive years. Plastic items were the dominant type of macrodebris recorded: over 95% of debris present in the three sampled marine compartments were plastic. In general, concentrations of macrodebris were quite high. Especially the number of beached debris reached very high levels: on average 6429±6767 items per 100 m were recorded. Microplastic concentrations were determined to assess overall abundance in the different marine compartments of the Belgian Continental Shelf. In terms of weight, macrodebris still dominates the pollution of beaches, but in the water column and in the seafloor microplastics appear to be of higher importance: here, microplastic weight is approximately 100 times and 400 times higher, respectively, than macrodebris weight.

Brief communication: the ecosystem perspective in ecotoxicology as a way forward for the ecological risk assessment of chemicals.

One of the objectives of the European Union (EU) ecological risk assessment of chemicals (ERA) is to derive maximum environmental concentrations that are not expected to cause adverse ecological effects. To this end, related EU directives list protection goals as well as guidelines that should be used to reach these goals. It is generally accepted that the individual-level endpoints on which these guidelines are based do not correspond to the listed population- and ecosystem-level protection goals. In this article, we identify 5 research topics that are key to bridging this gap: 1) the refinement of population-level effects and recovery rates by explicitly taking into account competition and 2) predation, 3) the assessment of chemical effects on biodiversity, 4) the assessment of chemical stress on ecosystem functions and services, and 5) the quantification of the effects of chemical mixtures. In addition, we illustrate why an ecosystem perspective is needed to address these topics and to inform the risk assessment process. We propose the use of existing ecotoxicological community, food web, and ecosystem models to tackle these issues and discuss why new models are needed to predict chemical effects on biodiversity.

A re-evaluation of fifteen years of European risk assessment using effect models.

Ecological risk assessments of chemicals can be informed by a suite of effect models, including population and food web models. In the risk assessments conducted under EU regulation 793/93/EC, however, applications of such effect models are extremely scarce and toxicity-extrapolation approaches are often used instead. The objective of the present study was to re-evaluate these risk assessments using two types of effect models: species sensitivity distributions (SSDs, non-mechanistic), and food web models (mechanistic). Species sensitivity distributions significantly fitted the available toxicity data for up to 35% of the chemicals, depending on the trophic levels included and the amount of data available. Median hazardous concentrations for 5% of the species (HC5-50) estimated by the SSDs were less accurate predictors of measured community-level no observed effect concentration than food web model-derived HC5-50s, albeit data were available for seven chemicals only. For datasets with more than 10 data points, the 90% confidence interval of the estimated HC5s was narrower for the food web modeling approach than for the SSD approach. The HC5-50s predicted by the two approaches were two to five times (metals) and 10 to 100 times (organic chemicals) higher than the predicted no effect concentrations (PNECs) for the aquatic environment listed in the risk assessment reports. This suggests that the derived PNECs are protective for aquatic ecosystems.

Assessment of species specificity of moulting accelerating compounds in Lepidoptera: comparison of activity between Bombyx mori and Spodoptera littoralis by in vitro reporter and in vivo toxicity assays.

BACKGROUND: Dibenzoylhydrazine analogues have been developed successfully as a new group of insect growth regulators, called ecdysone agonists or moulting accelerating compounds. A notable feature is their high activity against lepidopteran insects, raising the question as to whether species-specific analogues can be isolated. In this study, the specificity of ecdysone agonists was addressed through a comparative analysis in two important lepidopterans, the silkworm Bombyx mori L. and the cotton leafworm Spodoptera littoralis (Boisd.). RESULTS: When collections of non-steroidal ecdysone agonists containing different mother structures (dibenzoylhydrazine, acylaminoketone, tetrahydroquinoline) were tested, in vitro reporter assays showed minor differences using cell lines derived from both species. However, when compounds with high ecdysone agonist activity were examined in toxicity assays, larvicidal activity differed considerably. Of note was the identification of three dibenzoylhydrazine analogues with > 100-fold higher activity against Bombyx than against Spodoptera larvae. CONCLUSION: The present study demonstrated that species-specific ecdysone-agonist-based insecticides can be developed, but their species specificity is not based on differences in the activation of the ecdysone receptor but rather on unidentified in vivo parameters such as permeability of the cuticle, uptake/excretion by the gut or metabolic detoxification.

Do we have to incorporate ecological interactions in the sensitivity assessment of ecosystems? An examination of a theoretical assumption underlying species sensitivity distribution models.

Species sensitivity distributions (SSDs) are statistical distributions which extrapolate single-species toxicity test results to ecosystem effects. This SSD approach assumes that ecological interactions between populations, such as grazing and competition, do not influence the sensitivity of ecosystems. The validity of this assumption in a simple freshwater pelagic ecosystem was tested using ecosystem modelling. For each of a 1000 hypothetical toxicants, a lognormal SSD was fitted to chronic single-species EC10s of the species present. As such, these distributions did not account for ecological interactions and were therefore termed 'conventional SSDs' (cSSDs). Next, sensitivity distributions that did take into account ecological interactions were constructed (eco-SSD) for the same 1000 toxicants, using an ecosystem model. For 254 of the 1000 hypothetical toxicants, mean and/or variance of the cSSD were significantly higher than mean and/or variance of the eco-SSD, as such rejecting the general validity of the tested assumption. A classification tree approach indicated that especially toxicants which directly affect phytoplankton (i.e. herbicides) may have a higher mean for cSSD than for eco-SSD. Conversely, means of eco-SSD and cSSD tend to be equal for toxicants directly affecting zooplankton and fish, e.g. insecticides. For the 254 hypothetical toxicants for which the tested assumption was false, a predicted no effect concentration (PNEC) calculated as the lowest single-species EC10 divided by an application factor of 10 was on average a factor 10 lower than the corresponding ecosystem-NOEC calculated by the ecosystem model.

Chronic toxicity of copper to five benthic invertebrates in laboratory-formulated sediment: sensitivity comparison and preliminary risk assessment.

Five benthic organisms commonly used for sediment toxicity testing were chronically (28 to 35 days) exposed to copper in standard laboratory-formulated sediment (following Organization for Economic Cooperation and Development guidelines) and lethal and sub-lethal toxicities were evaluated. Sub-lethal endpoints considered were reproduction and biomass production for Lumbriculus variegatus, growth and reproduction for Tubifex tubifex, growth and emergence for Chironomus riparius, and growth for Gammarus pulex and Hyalella azteca. Expressed on whole-sediment basis the observed lethal sensitivity ranking (from most to least sensitive) was: G. pulex>L. variegatus>H. azteca=C. riparius=T. tubifex, with median chronic lethal concentrations (LC50) between 151 and 327 mg/kg dry wt. The sub-lethal sensitivity ranking (from most to least sensitive, with the most sensitive endpoint between parentheses): C. riparius (emergence)>T. tubifex (reproduction)=L. variegatus (reproduction)>G. pulex (growth)>H. azteca (growth), with median effective concentrations (EC50) between 59.2 and 194 mg/kg dry wt. No observed effect concentrations (NOEC) or 10% effective concentrations (EC10) for the five benthic invertebrates were used to perform a preliminary risk assessment for copper in freshwater sediment by means of (a) the "assessment factor approach" or (b) the statistical extrapolation approach (species sensitivity distribution). Depending on the data (NOEC or EC10) and the methodology used, we calculated a Predicted No Effect Concentration (PNEC) for sediment between 3.3 and 47.1 mg Cu/dry wt. This range is similar to the range of natural (geochemical) background concentrations of copper in sediments in Europe, i.e. 90% of sediments have a concentration between 5 and 49 mg Cu/kg dry wt. A detailed analysis of the outcome of this preliminary exercise highlighted that multiple issues need to be explored for achieving a scientifically more sound risk assessment and for the development of robust sediment quality criteria for copper, including (i) the use of the assessment factor approach vs. the statistical extrapolation approach, (ii) the importance of bioavailability modifying factors (e.g., organic carbon, acid volatile sulfide), and (iii) the influence of prevailing geochemical (bioavailable) background concentrations on the copper sensitivity of local benthic biota.

Relevance of generic and site-specific species sensitivity distributions in the current risk assessment procedures for copper and zinc.

Species sensitivity distributions (SSD) were constructed using acute toxicity data of various cladoceran species collected in five different aquatic systems. The aim of this research was to study the relative acute cladoceran community sensitivity in different aquatic systems. Current risk assessment procedures are based upon hypothetical communities and do not take into account variation in species composition and tolerance between aquatic communities. Two metals, copper and zinc, were used as model toxicants. To establish comparative sensitivity, a standard medium (International Organization for Standardization [ISO]) was used. The generic SSD (log-normal distribution) based on toxicity data obtained in this standard medium for all species (collected at all sites) resulted in a hazardous concentrations that protects 95% of the species occurring in a (hypothetical) ecosystem (i.e., hazardous concentration protecting 95% of the species of the hypothetical ecosystem [HC5]) of 6.7 microg Cu L(-1) (90% confidence limits: 4.2-10.8) and 559 microg Zn L(-1) (375-843). This generic SSD was not significantly different from the site-specific SSDs (i.e., constructed with species only occurring at a specific site). Mean community sensitivity (the geometric mean of 48-h 50% effective concentration [EC50] values of species within a community) among sites varied within a factor of 2 (between 17.3 and 23.6 microg Cu L(-1) for Cu and between 973 and 1,808 microg Zn L(-1) for Zn), and HC5s varied within a factor of 4 for copper (between 4.5 and 17.3 microg Cu L(-1)) and 7 for zinc (between 194 and 1,341 microg Zn L(-1)). For copper, the HC50 of our generic SSD was significantly lower than the one based on literature toxicity data of cladoceran species (which were recalculated to the hardness of our standard medium). In contrast, no significant differences were observed between the generic SSD and the literature-based SSD for zinc. It is suggested that the community sensitivity of different cladoceran populations is similar among aquatic systems and is not dependent on the species composition.

Probabilistic environmental risk assessment of zinc in Dutch surface waters.

In the framework of the European Union (EU) New and Existing Chemicals Policy, a regional risk assessment for Zn according to the current technical guidance documents and a probabilistic approach, by mathematically integrating both best-fitting exposure concentrations and species-sensitivity distributions into a probabilistic risk quotient distribution using Monte Carlo analysis, was explored for The Netherlands. Zinc is an essential element, and the current probability distributions may not adequately deal with this property. The threshold Pareto distribution provided the best fit to the chronic Zn toxicity data, resulting in a predicted-no-effect concentration (PNECadd) for dissolved Zn of 34.2 microg/L, whereas use of the conventional normal distribution resulted in a PNECadd for dissolved Zn of 14.6 microg/L. The extracted exposure data resulted in a regional predicted environmental concentration (PEC) for dissolved Zn in the Dutch surface waters of 20.1 microg/L and in PECadd values for dissolved Zn of between 15.5 and 17.3 microg/L, depending on the background correction used. The conventional deterministic risk characterization identified a regional risk for Zn in the Dutch surface waters. The more comprehensive probabilistic approach used in the present study, however, identified only very limited potential risks for the Dutch region. A probabilistic median risk, that the environmental concentration is greater than the no-observed-effect concentration of a species in Dutch surface waters (0.5-0.6%), depending on the inclusion of background correction, was obtained from the best-fitting distributions. Because probabilistic approaches provide a quantifiable and improved assessment of risk and quantification of the uncertainty associated with that assessment, these techniques may be considered as a way to improve the EU risk assessment procedures for data-rich substances.