Pesticide concentrations [average minimum (av-min), average maximum (av-max) and average of averages (av-av) values; ng/g] fresh weight to make it; ng/g of fresh weight in aquatic organisms, presented by continent and country; the number of quantified aquatic systems, pesticides and sampling year were also added (when more than one aquatic system, a range of values are presented).
Contamination by pesticides is a worldwide problem that can greatly disturb the biota, directly and/or indirectly. Nonetheless, few efforts were done so far to present review-style publications that analyse and integrate monitoring data—in a global scale—and evaluate possible environmental risks. Herein, we assessed possible environmental risks through theoretical calculations, using worldwide data published at least during the last 17 years and considering different trophic levels and the maximum average environmental concentrations (in water) observed in each continent. Furthermore, hazard quotients—using the estimated average daily intake, theoretical maximum daily intake and the maximum residue limits—were calculated to estimate the potential risks to humans through direct consumption of molluscs, crustaceans and fish. In summary, several pesticides were quantified at concentrations capable to affect low to medium trophic level species, which through the food web can affect higher trophic levels; theoretical approaches considering the environmental mixtures showed that algae and invertebrates are the most sensitive groups. Moreover, fish and crustaceans evidenced the highest body concentrations. To evaluate a potential risk through direct consumption, human health risk assessments were done, and in spite of no direct risk, some hazard quotients indicate a potential risk for developing carcinogenic effects.
- aquatic organisms
- hazard quotients
Worldwide, several studies have shown contamination with pesticides within different matrices. Together with the data shown previously in chapter “Pesticides in Worldwide Aquatic Systems: Part I”, information such as the maximum concentrations in waters and the concentration of pesticides found in the different biological matrices were used to (i) assess eventual individual pesticide risk through the comparison with the well-established EC50/LC50 for aquatic organisms, (ii) predict the environmental risk from pesticide mixtures found in each continent and (iii) assess the potential risk for human health when consuming molluscs, crustaceans and fish with the quantified concentrations.
2. Aquatic organisms
Fifty-two studies were used, where 111 different species were studied. The continent with the highest percentage of available results (quantified pesticides in different organisms) is Africa (39 species), followed by Europe (35 species), Asia (26 species) and then North and South America (nine and eight species, respectively). Here, we decided to focus on the sample type (zooplankton, molluscs, crustaceans, fish and mammals) analysed per continent and country (Table 1).
|Continent/country||Number of aquatic systems||Quantified pesticides||Sampling year||Sample type||av-min||av-max||av-av||References|
|Kenya||1||7||2011||F||na||na||0.3||[1, 6, 7]|
|Nigeria||5||1–65||2003–2014||F||19.5–3618||21.5–6355||20.6–5233||[1, 7, 8, 9]|
|China||5||7–45||2003–2013||F, Mo, C||0.6–2.5||1.8–34.5||0.8–11.8||[11, 12, 13, 14, 15]|
|Baltic Sea||1||18||2003||C, F, Mo||8.1||10.8||9.4|||
|Belgium||1||5||2001||F, Mo||1.5||7.6||4.1||[21, 22]|
|France||6||2–12||2001–2008||F, Mo||0.2–0.8||0.6–4.4||0.3–1.2||[21, 25, 26]|
|Italy||3||5–22||2002–2010||C, F, Mo||0.9–6.9||2.8–21.6||1.5–10.9||[10, 27, 28]|
|Portugal||4||1–54||2011–2013||F, Mo||4.6–7.6||27.2–72.0||0.2–18.6||[30, 31, 32, 33]|
|Spain||3||3–13||1996–2015||C, F, Mo||0.1–10.4||0.8–9.0||0.3–8.0||[34, 35, 36, 37]|
|Martinique Island||1||3||2003–2013||C, F, Mo||0.3||876.4||55.9|||
|USA||3||3–9||2004–2013||F, Mo||na||na||0.1–11.1||[43, 44, 45]|
|Brazil||3||6–42||1996–2009||F, Mo, Ma||0.1–30.1||0.1–410.7||0.1–99.8||[47, 48, 49, 50]|
The data collected between 1993 and 2016 averaged from 0.1 to 5233 ng/g (Table 1). Europe is represented by 22 aquatic systems, followed by Africa, with 13, and the rest with no more than nine aquatic systems. When considering the number of pesticides quantified, Africa has more observations (382) than Europe (327) and the other continents (between 90 and 220), which is due to the higher number of species studied in Africa.
Africa stands out with average concentrations of 132 ng/g (SD 411), followed by Europe (57 ng/g, SD 271), South America (17 ng/g, SD 40) and Asia and North America (5 ng/g, SD 20). This scattered difference between concentrations is mainly due to the average values observed in Warri River (Nigeria, Africa) and in the Danube Delta (Romania, Europe) [1, 2].
Grouping data by category, insecticides prevail in 89% of biologic analyses, leaving 11% for the herbicide and fungicide categories and presenting the same pattern on all continents (Figure 1). No data are available for Oceania and Antarctica; so, when citing herein “worldwide”, these continents will not appear.
Analysing data by matrix, the most common analysis occurred in fish (74%) and molluscs (20%). The remaining studies considered zooplankton, crustaceans, turtles and mammals (Figure 1). In total, 74% of the quantified pesticides were conducted on vertebrates and the other 26% in invertebrates (Figure 1). While for the latter, 80% of the quantifications were done using the whole animal, and for vertebrates, it is further divided; specific organs or tissues were used to quantify pesticides.
Many factors account for the strong bias towards vertebrates. Invertebrates are small, less complex and as a food resource almost entirely eatable, while the same is not applicable to vertebrates. Besides that, the study goal (i.e., food control or environmental/toxicological studies) also influences the type of the tissue/organ to study (muscle, liver, gonads or gills). For example, the bubbler tissue and fat are only applicable for aquatic mammals and turtles (Figure 2).
Results per Metazoan lineages (zooplankton, mollusc, crustacean, fish, turtles and mammals) were assessed considering the average concentrations and the number of quantifications (Figure 3). Average concentrations were ~11 ng/g for zooplankton and molluscs, ~35 ng/g for mammals and ~100 ng/g for crustaceans and fishes (Figure 3A).
Among continents, Africa presented the highest concentrations for crustaceans (142 ng/g) and fishes (253 ng/g) followed by North America (136 ng/g for crustaceans). Asia, Europe and South America included data belonging to four Metazoan lineages, with similar range of concentrations (~3 to ~76 ng/g (Figure 3B)).
3. Half effective and lethal concentrations (EC50/LC50) for aquatic organisms
It is now well established that at specific concentrations all pesticides are harmful to biota, affecting algae and plants, invertebrates and vertebrates . Databases such as Pesticides Properties DataBase (PPDB) present information on the physicochemical properties, environmental fate, human health and ecotoxicological data of all active ingredients and approved pesticides .
In order to evaluate the worst-case scenario, the maximum average concentrations measured in waters from each continent were compared against the acute and chronic concentrations for aquatic animals, documented by the PPDB (see chapter “Part I: Pesticides in Worldwide Aquatic Systems”). On a global scale, 57 pesticides were registered at maximum average concentrations above the LC50 and/or EC50 settled for algae, invertebrates and/or fishes; among continents, Europe reported the highest number of pesticides (44 of 116), followed by Asia (14 of 42), Africa, Oceania and finally South America (6 of 24).
The most critical measured environmental concentrations (MEC) were registered for dicofol, ethion (Asia), metribuzin (Europe) and diazinon (Africa) with values from 2- to 200-folds higher than EC50 or LC50 set for invertebrates and algae.
4. Predictive aquatic risk assessment of pesticide mixtures
Despite the common occurrence of pesticides, mixtures in the environment, laws, conventions and recommendations still focus on individual standard parameters. Modelling approaches, based on available ecotoxicological information, can be used to estimate the impact of mixtures in the biota, completing this lack of information .
Based on the European chemical legislation REACH, the ecological risk quotient (RQ) is determined by the equation:
PNEC is derived by selecting the most sensitive trophic level—from algae, invertebrate or fish—and applying an appropriate assessment factor (AF) [52, 54]. The AF, also denoted as safety or uncertainty factor, considers intra- and inter-laboratory variation of the data, biological variance and short-term to long-term exposures, presenting stipulated values for specific conditions [55, 56]; as an example, considering the Maximum Acceptable Concentration-Quality Standards (MAC-QS) to assess short-term effects an AF = 100 should be applied .
The RQ values, classified from <0.01 (negligible) to >1 (significant), indicate a range of potential risks for concern, but do not inform about the specific biological end point for that organism which is representing a specific trophic level [53, 57]. For this reason, a second approach, which defines the most sensitive trophic level for that environmental concentration, should be applied :
RQTU values are summed per trophic level (sum of the toxic units (RQSTU)). If both RQ(MEC/PNEC) and RQSTU are >1, additional considerations are required . Based on the two reference models—concentration addition (CA) and independent action (IA)—the RQSTU/MaxTU can be used to predict the second tier, resulting in the maximum value from which CA may display higher toxicity values than IA .
In this work, the maximum of the average measured concentration of pesticides in water samples was used to assess the potential risk per continent and on a worldwide scale (Table 2). From a total of 144 pesticides quantified in water samples, 133 were used for ecological risk assessment (Table 2); the remainders, mostly isomers and metabolites, were not integrated due to lack of information on their EC50 and LC50 concentrations set for these trophic levels (algae, invertebrate and fish). The highest number of pesticides suitable for this approach are represented by insecticides (n = 118). In general, algae was the most sensitive group to herbicides and fungicides, with 75% and 61.5% of the cases, respectively, while invertebrates showed the highest sensitivity to insecticides (66.1%) (Table 2).
Globally, the RQ(MEC/PNEC) resulted in 43% of very high-risk cases, led by insecticides; fungicides were the least worrisome category, as most of the cases presented low or negligible risks (Figure 4).
The results presented above are a consequence of the highest values measured around the world. Since Europe was the continent with more values of RQ(MEC/PNEC), these results are mostly representative for this continent (see Table 1). However, this does not mean that concentrations measured on the other continents are innocuous. As observed for the number of compounds analysed per continent, Africa presented the most disturbing scenarios (52%), followed by Asia and Europe (45%) and then Oceania and South America (24%) with RQ > 1 (Figure 5).
In order to evaluate the effect of the maximum average concentrations found per individual trophic level (RQTU), further evaluation should be done through RQSTU (Table 3).
When comparing between continents, the highest RQSTU ratios were attained in Europe, for algae (16.13) and fish (33.12), and in Asia for invertebrates (324.97); however, the last one is due to a punctual concentration observed in India for ethion . Independently of that, the invertebrate group is the most sensitive trophic level, presenting the highest RQSTU values. The same pattern is observed in the other continents except in Oceania, where the highest risk is observed for the algae (0.92) by the herbicides (Table 3).
The RQ(MEC/PNEC) and RQSTU demonstrate that one or more biotest organisms are sensitive to the concentrations presented on that continent; so, the ratio RQSTU/highest RQTU was done, applying the highest sum among trophic levels (Table 4).
|Continent||No. of compounds (toxic/total)||RQTU||∑RQSTU||∑RQSTU/RQTU|
For each of these scenarios, the maximal possible ratio RQSTU/RQTU was lower than the value given by the number of mixture of toxic components, suggesting that the possible observed toxicity is due to a low number of pesticides. As we can notice, the RQSTU/RQTU ratio is very similar among continents; however, the number of toxic compounds per total, where Africa presents a significant number (52%) when compared to the others, should be also considered (Table 4).
5. Human health risks
Dietary pesticide risks can be estimated through well-established indices and defined and used by US Environmental Protection Agency (EPA) , European Food Safety Authority (EFSA) and European Union Directives [61, 62]. Realistic predictions involve several parameters, such as pesticide residue intake (PRI), as the one reported by Food and Agriculture Organisation (FAO) :
The acceptable daily intake (ADI) estimates the amount of a substance in food that can be ingested daily over a lifetime without appreciable health risk to the consumer :
The estimated average daily intake (EADI), according to EPA, should be less than the established ADI values :
The theoretical maximum daily intake (TMDI) represents the maximum concentration of a pesticide residue (mg/kg) legally permitted in food :
When no specific MRL is published, a 0.01 mg/kg value is applied . Additionally, hazard quotients (HQs)—which measure the potential exposure for developing non-carcinogenic health effects—may be calculated using several assumptions .
EADI may be divided by the acute reference dose (ARfD, mg/kg/day) —which is derived from the no-observed-adverse-effect levels (NOAEL) and based on studies of short time exposures (1–7 days) , by ADI, for long intake periods, and by TMDI, which is advised by EFSA to calculate the potential risks of unintentional compounds, such as pollutants.
In the chapter Part I: Pesticides in Worldwide Aquatic Systems, the levels/categories of pesticides per continent/country are displayed. The maximum average concentrations shown in Part I were used here to assess human health risks. Data are summarised in Table 5.
|Africa||—||1.6E + 00 (30)||5.2E + 00 (370)|
|Asia||3.2E-01 (41)||1.0E-02 (3)||9.0E-02 (173)|
|Europe||6.0E-02 (143)||1.4E-01 (10)||4.0E + 00 (168)|
|North America||1.0E-02 (16)||1.4E-01 (1)||3.0E-02 (104)|
|South America||1.8E-01 (46)||—||1.7E-01 (32)|
The highest concentrations were observed in fish from Africa (5.2 mg/kg) and Europe (4.0 mg/kg), followed then by crustaceans in Africa (1.6 mg/kg). The highest number of cases (number of quantifications found considering all the pesticides, species and countries) was registered for fish (169 average cases), followed by molluscs (62 average cases), and finally crustaceans (11 average cases). The elevated number of fish studies is likely due to their importance as a food source.
For allowing a detailed evaluation of human health hazard, the same data is displayed by pesticide and continent and organised considering molluscs, crustaceans and fishes (Table 6). The food consumption rate and the average adult body weight were defined by continent [63, 64]. For the compounds endrin ketone and aldehyde, HCH (sigma and lambda), pretilachlor and pentachlorobenzene, a MRL of 0.01 mg/kg was adopted, since no specific data was found.
|∑Aldrin + dieldrin||1.2E+00||5.6E-04||6.0E-03||1.8E-04||1.0E-04||3.0E-03||0.09||3.19||5.62||0.19|
|∑Aldrin + dieldrin||1.0E-02||1.0E-05||6.0E-03||3.5E-04||1.0E-04||3.0E-03||0||0.03||0.1||0|
|∑Aldrin + dieldrin||1.5E-01||4.7E-05||6.0E-03||1.5E-04||1.0E-04||3.0E-03||0.01||0.31||0.47||0.02|
Focusing on the molluscs results, the MEC of 15, 52, 10 and 16 pesticides (from Asia, Europe, North America and South America, respectively) were used to calculate the HQs. Due to the low ratio values, only cases with at least one ratio value above 0.05 were presented. As we can see, none of the results proved to be harmful to human through direct consumption. In other words, none of the ratios were above 1, indicating that the calculated EADI was below the reference levels (MRL, TMDI, ADI and ARfD). The highest HQ(EADI/TMDI) was obtained for methoxychlor in Asia (0.55). For HQ(EADI/ADI), the highest ratio occurred in South America for mirex with 0.54.
Looking to the crustacean data, a total MEC of eight, two, three and one cases from Africa, Asia, Europe, North America, respectively, were analysed. The same criterion, which is the case with at least one ratio value above 0.05, was applied. High HQs for chlordane (gamma) were observed in crustaceans sampled in Africa (see Table 6). In spite of that, none of the ratios were above 1.
Twenty-four (Africa), 16 (Asia), 10 (Europe), 28 (North America) and 21 (South America) MEC cases were analysed considering the fish data. Once again, only HQ ratios with at least one case above 0.05 are shown. As we can see, none of the maximum average concentrations were above the MRL values; however, several HQ > 1 are observed in Africa, bringing potential exposure for developing carcinogenic health effects. This fact may be a result of bioaccumulation processes (where concentrations increase in higher trophic levels) and/or a higher interest in this matrix (increasing the data availability and diversity). These ratios were registered for six compounds—∑aldrin + dieldrin, endrin aldehyde, paraquat dichloride, endrin, heptachlor and heptachlor epoxide—where the most preoccupant cases (HQ > 3) are for the first three pesticides cited above.
6. Final considerations
Globally, and because of these high average concentrations, several individual pesticides were quantified at levels exceeding the established LC50 for fish and EC50 for invertebrates and algae.
In addition, the review has provided clear evidence that the biological data grouped according to Metazoan lineages reached higher concentrations for fish and crustaceans (Figure 3). It is worth noting however that the same pattern was not verified for higher trophic levels including turtles and aquatic mammals which may be due to the lack of samples. Considering that globally, many of the data displayed a wide range of concentrations, coupled with the fact that many of the larger aquatic species are migratory; there is a need to address the pesticide problem from a global perspective.
As a complement to this work, all edible species were evaluated for dietary pesticide risks, as mollusc, crustaceans and fish. No direct human health risk was observed; however, in Africa, some hazard quotients (HQ) were above one, indicating a potential exposure for developing carcinogenic health effects.
In conclusion, the potentially harmful effects of pesticides should be considered not only locally (national/governmental institutions) but also on a global scale.
This study was partially supported by the European Regional Development Fund (ERDF) through the Operational Competitiveness Programme (COMPETE) and Operational Human Potential Programme (POPH) and by national Portuguese funds, through Foundation for Science and Technology (FCT), via project UID/Multi/04423/2013 and project PTDC/MAR/70436/2006 (FCOMP-01-0124.FEDER-7382) and grant (SFRH/BD/79305/2011) and grant (SFRH/BPD/120558/2016). This work was also implemented in the framework of the Structured Program of R&D&I INNOVMAR (Innovation and Sustainability in the Management and Exploitation of Marine Resources; Reference NORTE-01-0145-FEDER-000035), namely, within the research line ECOSERVICES, supported by the Northern Regional Operational Programme (NORTE2020), through ERFD. Complementary funding was provided by ICBAS-U.Porto.