Pesticide mixtures in the Swedish streams: Environmental risks, contributions of individual compounds and consequences of single-substance oriented risk mitigation
Graphical abstract
Introduction
Multiple studies have demonstrated that complex pesticide mixtures are present in surface waters globally, in the US (e.g. Gilliom, 2001, Stone et al., 2014a, Stone et al., 2014b), Europe (e.g. Moschet et al., 2014, Schreiner et al., 2016, Ccanccapa et al., 2016) and elsewhere (e.g. in South America (Hunt et al., 2016), Australia (Allinson et al., 2015) and China (Zhang et al., 2011). Empirical evidence univocally shows that the combined toxic effects of such pesticide mixtures exceed the effect of each individual compound (e.g. Faust et al., 2001, Faust et al., 2003, Knauert et al., 2009 and Porsbring et al., 2010, see also reviews by Belden et al., 2007, Verbruggen and van den Brink, 2010 and Rodney et al., 2013).
Studies have repeatedly demonstrated that Concentration Addition (CA) describes the joint toxicity of pesticide mixtures well (reviewed by Belden et al., 2007, Rodney et al., 2013). This implies that all components contribute to the overall mixture toxicity, independently of whether they are present at concentrations above or below their individual No Observed Effect Concentration (NOEC) or Environmental Quality Standard (EQS). Mixtures might therefore cause toxic effects even if all components are present at concentrations below which an individual effect is detectable (e.g. Carvalho et al., 2014, Faust et al., 2001).Taken together, the available body of evidence thus clearly shows that pesticide mixtures warrant specific consideration during environmental risk assessment, monitoring and management.
Environmental risks of pesticides and pesticide mixtures are assessed in the European regulatory system from two perspectives. First, active ingredients and whole formulated pesticide products are evaluated for their environmental hazard, exposure and risk during market authorization (EFSA, 2013), following the legal frameworks that are laid down in Regulations EC 1107/2009 on the placing of plant protection products on the market (European Parliament, 2009a) and EC 546/2011 on uniform principles (European Commission, 2011a). However, ‘coincidental’ pesticide mixtures, i.e. mixtures of active ingredients that result from farmers applying different pesticide products in close proximity to a given water body or because of sequential spraying of different pesticides on the same field, are not considered in Regulation EC 1107/2009 nor in Directive 2009/128/EC (European Parliament, 2009b). However, it has been argued that the uniform principles in Regulation 546/2011/EC require authorization of plant protection products to be based on the “proposed conditions for use” and consequently – given common agricultural practice – to consider the environmental impact of the resulting pesticide mixtures (Frische et al., 2014).
Second, the impact of mixtures of pesticides (and other hazardous chemicals) on the ecological status of an aquatic system is assessed from the perspective of the Water Framework Directive (WFD) (European Parliament, 2000). In order to be classified as having a good ecological status, a water body also needs to have good chemical status, which requires that the concentrations of each of 45 priority pollutants, which are currently listed in Directive 2013/39/EC (European Parliament, 2013), do not exceed European-wide thresholds, so-called Environmental Quality Standards (EQS). In addition, in order to track progress towards the national goal of a “non-toxic environment” (adopted in 1999), Sweden also developed national Water Quality Objective(s) (WQO) for pesticides, defined as concentrations which are not expected to cause any adverse effects in the aquatic environment (Norberg, 2004, Lindström and Kreuger, 2015). These values are similar to EQS values and serve as a tool to evaluate surface water quality based on monitoring results, but are not legally binding. WQO's are derived using a method that closely follows the REACH approach for deriving Predicted No Effect Concentrations (PNEC) values, based on single species data and assessment factors between 10 and 1000, depending on the underlying ecotoxicological endpoints (Andersson et al., 2009, Andersson and Kreuger, 2011, KEMI, 2008).
Risk assessment of chemical mixtures is routinely performed using CA (Kortenkamp et al., 2009, Bopp et al., 2015). CA has also been suggested specifically for the assessment of pesticide mixtures (EFSA, 2013) and it is the recommended approach for setting EQS values for chemical mixtures within the context of the WFD (European Commission, 2011b).
According to CA the risk quotient (RQ) of a mixture, RQCA, is defined as:where cmix is the total concentration of the mixture, ECxMix is the mixture concentration causing x% effect, while ci and ECxi denote the corresponding concentrations of substance i. The ratio ci/ECxi provides a dimensionless measure of the toxicity contribution of compound i usually termed a Toxic Unit (TU). Although the concept is rooted in the idea of the mixture components sharing the same mode of action, as well as not taking possible synergistic (or antagonistic) effects into account (Cedergreen, 2014), CA has been successfully used for the risk assessment of heterogeneous mixtures (Belden et al., 2007, Kortenkamp et al., 2009, Verbruggen and van den Brink, 2010, Rodney et al., 2013, Bopp et al., 2015). The toxicity estimates in Eq. (1) (ECxmix and ECxi) in principle refer to the same ecotoxicological endpoint recorded for the same species under identical exposure conditions. However, in practice CA is often applied in a broader setting, e.g. by using data from different algal species in order to predict the toxicity to algae in general.
In the present paper, we have applied CA in order to separately calculate the risks for algae, crustaceans and fish. The corresponding CA-based mixture RQs are termed RQAlgae, RQCrust and RQFish. Moreover, by substituting the ECxi with the WQOi and ci with the Measured Environmental Concentration (MECi) we determined ecosystem-wide RQWQO values as the sum of the individual MEC/WQO ratios, following the rational outlined by (Backhaus & Faust,
Comparing the trophic-level specific RQs with RQWQO is difficult, as the latter is calculated using assessment factors to account for the different amounts of data available for each compound, while RQAlgae, RQCrust and RQFish are calculated without using any assessment factors. In order to bridge these two approaches, we therefore also calculated a mixture RQ for the most sensitive trophic level (RQMST), defined as:
RQMST provides a measure for the risk across trophic levels, but is calculated without using any assessment factors. It thus takes an interim position and bridges the trophic-level specific RQs (RQAlgae, RQCrust or RQFish) to the ecosystem-wide RQWQO. The RQMST is conceptually identical to the point of departure index (PODI), frequently used in human toxicology (Wilkinson et al., 2000).
A RQ provides a yardstick for assessing the need to act. Values of RQWQO exceeding 1 indicate the need for either a more advanced mixture risk assessment, and/or for the implementation of risk mitigation measures. We defined the corresponding critical values for RQAlgae, RQCrust, RQFish as 0.1, 0.01 and 0.01, respectively, following the assessment strategy for individual pesticides (EFSA, 2013). Defining a critical value for RQMST is not feasible at the moment, as no strategy has been suggested yet on how an overall assessment factor should be calculated that reflects the overall uncertainty in Eq. (2). The RQMST will always be higher than any of the organism-group specific RQs (Backhaus and Faust, 2012) and, because no assessment factors are applied, lower than the RQWQO.
The ratio between the total RQ of a mixture and the maximum RQ of its components has been termed the maximum cumulative ratio (MCR, Price and Han, 2011). That is,
If all components of a mixture are contributing equally to the predicted mixture risk, the MCR equals the number of compounds in the mixture. In a mixture whose TU distribution is dominated by one compound, the MCR approaches 1. Therefore, the MCR has been suggested as a tool to assess the value of performing mixture toxicity assessments (Price and Han, 2011).
Chemical risk assessment is in general based on comparing relevant exposure estimates (measured or modeled) with hazard estimates, such as NOEC's, EC50’s and EQS values. Such estimates are straightforward to calculate on the basis of monitoring results, as long as detected environmental concentrations are quantified, either above the chemical-analytical limit of quantification (LOQ) or the limit of detection (LOD). However, sometimes when the detection is below the LOQ but still above the LOD the concentration is not quantified (only given as ‘trace’) in order to save time in the laboratory. Nevertheless, reasonable assumptions on the trace concentrations present can be made using (LOQ + LOD) / 2 as a surrogate for unquantified detections between the LOQ and the LOD, as long as these two parameters are stated in the analytical protocol.
However, the situation becomes problematic if a monitored chemical is not detected. Such a result does not prove that the compound is not present, it only shows that the concentration is somewhere between zero and the LOD. Assuming a zero concentration for all non-detects will therefore underestimate the total risk, if no additional knowledge about e.g. emission or use pattern is available.
On the other hand, assuming that all non-detected compounds are present just below their LOD – the worst-case scenario that is still compatible with the recorded values – is also unrealistic. Such an approach immediately leads to the logical inconsistency that the estimated risk becomes simply dependent on the number of compounds analyzed. The same is true for setting the concentration used for the risk assessment a priori to any other value above zero.
Parametric and non-parametric statistical methods are available for data with “less-than” values, i.e. findings of concentrations < LOD. They allow the estimation of the likely contribution of non-detects to the total RQ. In this paper we used the non-parametric Kaplan-Meier (KM) method (Helsel, 2010, Helsel, 2012, Bolks et al., 2014), because it is not possible to ensure that the distributional assumptions of parametric alternatives are fulfilled in the analyzed data. A KM-adjusted sum of RQs lies between the sum of RQs that result from substituting all non-detects with their respective LOD and the sum of RQ that results from substituting all non-detects with zero.
The KM-method ignores the potential risk contribution of a compound, if its potential RQ exceeds the maximum of the RQs that are based on a quantified concentration. For such compounds, better analytical data are required for a reliable quantification of their risk contribution.
The southern part of Sweden is an area of intense agricultural activity and pesticide residues have been systematically monitored at six sites since 2002 (Lindström et al., 2015, Lindström and Kreuger, 2015). In this paper, we applied CA-based risk assessment approaches in order to estimate and characterize the environmental risks from the detected pesticide mixtures, using RQAlgae, RQCrust, RQFish, RQMST and RQWQO. The results will then be used for a broader discussion on the impact of non-detects on component-based mixture risk assessments. Finally, we explore the consequences of a single-substance oriented risk management, i.e. assuming that risk mitigation measures ensure that all individual concentrations are below their corresponding WQO's.
In order to explore how the different possibilities to incorporate (or ignore) concentrations below the LOD influence the final mixture risk estimates, we calculated all RQs for three different exposure scenarios (Table 1). Scenario 1 and 2 assumes that non-detects are present at a concentration equal to their LOD or at zero, respectively. Scenario 3 uses the KM-adjustment for compounds present < LOD.
Section snippets
Pesticide monitoring data
As part of a continuous Swedish pesticide monitoring program the Swedish University of Agricultural Sciences publishes data on pesticide concentrations in four streams draining 8–16 km2 and two rivers draining 102–488 km2 at http://jordbruksvatten.slu.se/pesticider_start.cfm (Agricultural land, 2017, Lindström and Kreuger, 2015). The chemical monitoring data was quality-checked and is made publically available as a downloadable datafile for a broader audience via GitHub (//github.com/ThomasBackhausLab/Swedish_Pesticide_Data.git
Results and discussion
The ecotoxicological risk of the pesticide mixtures found in Swedish freshwater ecosystems was previously described by Bundschuh et al. (2014) for the timeframe from 2002 to 2011. In this paper, we analyze three additional issues: Firstly, we explore the relevance of non-detects for the overall mixture risk. Secondly, we compare the specific risks for the three main organism groups, i.e. algae, crustaceans and fish with ecosystem-wide risks. Finally, we analyze the impact of successful
Conclusions
The presented risk analysis concludes that pesticide residues frequently put aquatic ecosystems in Southern Sweden at risk. This is in line with previous studies in aquatic ecosystems elsewhere (see above). Using WQO values produced by the Swedish Chemicals Agency and the Swedish University of Agricultural Sciences as our basis we conclude that the risk posed by pesticide mixtures were unacceptably high in 73% of the analyzed samples (when using the KM adjustment for non-detects). The fact that
Acknowledgements
Funding by the European Commission (project SOLUTIONS, FP7-ENV-2013, grant agreement 603437) and the Swedish Research Council FORMAS (project IMPROVE and NICE, grant nos. 2010-2014-1026 and 210-2011-1733) is gratefully acknowledged. The Swedish national pesticides monitoring program is funded by the Swedish Environmental Protection Agency.
References (60)
- et al.
Future water quality monitoring — adapting tools to deal with mixtures of pollutants in water resource management
Sci. Total Environ.
(2015) - et al.
The solutions project: challenges and responses for present and future emerging pollutants in land and water resources management
Sci. Total Environ.
(2015) - et al.
Evaluation of pesticide monitoring strategies in agricultural streams based on the toxic-unit concept - experiences from long-term measurements
Sci. Total Environ.
(2014) - et al.
Pesticides in the Ebro River basin: occurrence and risk assessment
Environ. Pollut.
(2016) - et al.
Predicting the joint algal toxicity of multi-component s-triazine mixtures at low-effect concentrations of individual toxicants
Aquat. Toxicol.
(2001) - et al.
Insecticide concentrations in stream sediments of soy production regions of South America
Sci. Total Environ.
(2016) - et al.
Analysis of the occurrence and risk assessment of polar pesticides in the Llobregat River Basin (NE Spain)
Chemosphere
(2012) - et al.
Analysis of monitoring programmes and their suitability for ecotoxicological risk assessment in four Spanish basins
Sci. Total Environ.
(2012) - et al.
Pesticide mixtures in streams of several European countries and the USA
Sci. Total Environ.
(2016) - et al.
Assessing, mapping and validating site-specific ecotoxicological risk for pesticide mixtures: a case study for small scale hot spots in aquatic and terrestrial environments
Ecotoxicol. Environ. Saf.
(2011)
Assessing the risks of exposures to multiple chemicals with a common mechanism of toxicity: how to cumulate?
Regul. Toxicol. Pharmacol.
Database
Pesticide and trace metal occurrence and aquatic benchmark exceedances in surface waters and sediments of urban wetlands and retention ponds in Melbourne, Australia
Environ. Sci. Pollut. Res. Int.
Beräkning av temporära riktvärden för 12 växtskyddsmedel i ytvatten
Preliminära riktvärden för växtskyddsmedel i ytvatten, beräkning av riktvärden för 64 växtskyddsmedel som saknar svenskt riktvärde
Predicitve environmental risk assessment of chemical mixtures: a conceptual framework
Environ. Sci. Technol.
How well can we predict the toxicity of pesticide mixtures to aquatic life?
Integr. Environ. Assess. Manag.
Baseline Assessment of Leftcensored Environmental Data Using R. Tech Notes 10, June 2014
Scientific Methodologies for the Assessment of Combined Effects of Chemicals – A Survey and Literature Review
Mixtures of chemical pollutants at European legislation safety concentrations: how safe are they?
Toxicol. Sci.
Quantifying Synergy: A Systematic Review of Mixture Toxicity Studies within Environmental Toxicology
PLoS ONE
Review of aquatic in situ approaches for stressor and effect diagnosis
Integr. Environ. Assess. Manag.
Guidance on tiered risk assessment for plant protection products for aquatic organisms in edge-of-field surface waters
EFSA J.
Guidance Document on Aquatic Ecotoxicology in the Context of the Directive 91/414/EEC, Sanco/3268/2001 Rev.4 (Final)
Commission regulation (EU) no 546/2011 of 10 June 2011 implementing regulation (EC) no 1107/2009 of the European Parliament and of the Council as regards uniform principles for evaluation and authorisation of plant protection products
Off. J. Eur. Union
Common Implementation Strategy for the Water Framework Directive (2000/60/EC) - Guidance Document No. 27 For Deriving Environmental Quality Standards. Technical Report - 2011 - 055
Directive 79/117/EEC of 21 December 1978 prohibiting the placing on the market and use of plant protection products containing certain active substances
Off. J. L
Directive 91/414/EEC of 15 July 1991 concerning the placing of plant protection products on the market
Off. J. L
Cited by (86)
A multi-scenario risk assessment strategy applied to mixtures of chemicals of emerging concern in the River Aconcagua basin in Central Chile
2024, Science of the Total EnvironmentApplication of confocal microscopy and flow cytometry to identify physiological responses of Prorocentrum micans to the herbicide glyphosate
2024, Marine Environmental ResearchWater Framework Directive micropollutant monitoring mirrors catchment land use: Importance of agricultural and urban sources revealed
2024, Science of the Total EnvironmentAssessment of the contribution of surfactants to mixture toxicity in French surface waters
2023, Science of the Total EnvironmentEstimating the aquatic risk from exposure to up to twenty-two pesticide active ingredients in waterways discharging to the Great Barrier Reef
2023, Science of the Total Environment
- 1
Visiting address: Carl Skottsbergs gata 22 B, 413 19 Göteborg, Canada.
- 2
Visiting address: Lennart Hjelms väg 9, 756 51 Uppsala, Canada.