Next Article in Journal
Advanced Treatment of Laundry Wastewater by Electro-Hybrid Ozonation–Coagulation Process: Surfactant and Microplastic Removal and Mechanism
Next Article in Special Issue
Influence of Salinity on Copper Toxicity in Paracentrotus lividus and Arbacia lixula Embryos
Previous Article in Journal
Collaborative Ecological Flow Decision Making under the Bengbu Sluice Based on Ecological-Economic Objectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ecotoxicological and Chemical Approach to Assessing Environmental Effects from Pesticide Use in Organic and Conventional Rice Paddies

1
Institute for Environmental Protection and Research, 00144 Rome, Italy
2
Research and Science Department, Italian Space Agency, 00133 Rome, Italy
3
Regional Environmental Protection Agency (ARPA), 10135 Torino, Italy
*
Author to whom correspondence should be addressed.
Water 2022, 14(24), 4136; https://doi.org/10.3390/w14244136
Submission received: 18 November 2022 / Revised: 14 December 2022 / Accepted: 15 December 2022 / Published: 19 December 2022
(This article belongs to the Special Issue Ecotoxicological Risk in Aquatic Environments)

Abstract

:
Despite laws and directives for the regulation and restriction of pesticides in farming, the large use of Plant-Protection Products (PPPs) in paddy fields is a relevant worldwide cause of environmental contamination. The aim of this work is to evaluate the environmental impact due to the use of PPPs by using an integrated approach based on chemical analyses and ecotoxicological hazard assessment, supported by statistical tools, in order to overcome the issues related to traditional tabular evaluation. Samples of soil and water of seven conventional and organic paddies located in Northern Italy were examined for two years. The results evidenced a direct relationship between the presence of Oxadiazon in water and bioassay responses as the main cause of the toxicity measured. This phenomenon affected both biological and conventional rice fields, due to the free circulation of water through irrigation canals. Therefore, the implementation of organic districts with water circulation isolated from conventional fields represents a simple and effective countermeasure to safeguard the agricultural practices of organic crops.

1. Introduction

For a long time, it has been known that the large amount of pesticides applied in paddy fields, in addition to the common practice of draining the paddy water in irrigation canals that flows into the freshwater system and eventually into the marine environment, is one of the major causes of pollution worldwide [1].
Plant-Protection Products (PPPs), because of their persistence, toxicity, and bioaccumulative properties, are of particular concern. They might have adverse ecological effects, causing both short-term (acute) and long-term (chronic), lethal or sub-lethal biological damage; in particular, changes in behavior, metabolism, development, alteration in the food chain or habitat of non-target organisms, such as amphibians and bats, and reduction in the populations of natural predators of insect pests [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]. In addition, pesticides are known to decrease the biodiversity and function of an ecosystem by promoting the dominance of undesired and invasive species [18].
The most common way to assess the biological effects of the use of PPPs derives from toxicological studies with the aim to define the existing relationship between doses of specific compounds and toxicity response in laboratory experiments. Relatively few field surveys were conducted with the purpose of correlating the in situ measurements of pesticide concentration with the bioassay responses, i.e., Daphnia magna [19]. Some studies were carried out in controlled microcosms on a range of representative soil organisms [20], while several studies used the Cornell University Environmental Impact Quotient (EIQ) calculator [21] and the ECOTOX Knowledgebase to determine the exposure risk associated with individual pesticides relative to their application rates and aquatic concentrations [22]. Moreover, Sánchez-Bayo and Goka [23] evaluated the ability of four community endpoints (species richness, abundance, diversity, and similarity indices) to assess the impacts of two insecticides (Imidacloprid and Etofenprox). The ecotoxicity of a mixture of PPPs and behavior of pesticide transformation products directly in the aquatic environment, such as paddies, are poorly understood.
The European Community established various actions to counter the impact of pesticides on biodiversity through the Directive 2009/128/CE [24], which was implemented in Italy with the Legislative Decree no. 150/2012 [25], compelling a minimization or prohibition of the use of pesticides in the areas designated by the Habitat (92/43/EEC) and Birds (2009/147/EC) Directives [26,27] and in the protected areas referred to in the Water Framework Directive (2000/60/EC) [28]. Europe and United States (USEPA) established Regulation (EC) No. 1107/2009 [29] and the Data Requirements for Conventional Chemicals for pesticide registration, requiring environment fate and ecotoxicity data to be provided [30].
These regulations are referred to as a tabular approach, which is limited to the compliance of the chemical thresholds with respect to the toxicity of pure individual pesticides. Moreover, this approach does not take into account the simultaneous action of the complex mixtures of contaminants commonly present in aquatic habitats that may result in antagonistic and more often synergistic effects on biota.
One of the most advanced pesticide actions is adopted in Japan by the Ministry of Environment [31]. In order to determine the eligibility of the product, acute toxicity tests must be conducted for fish (basically, Cyprinus carpio), daphnids (Daphnia magna), and algae (Raphidocelis subcapitata), and then the minimum value of the 50% effect concentration (EC50 or LC50) is divided by an uncertainty factor that considers the species sensitivity. The Japanese Agricultural Chemicals Regulation Law was revised in 2018, and the method of assessing pesticide registration criteria was also revised [31]. Toxicity tests using aquatic plants, such as Lemna sp., in addition to algae will be introduced in the setting of criteria for herbicides. The uncertainty factor applied to the algal EC50 was changed from 1 to 10 by default, which is then reduced depending on the number of algal species tested. However, registration criteria for eligibility of pesticides for the new method were not developed until 2021 [32].
A relatively common way to determine hazardous concentration for the protection of an ecosystem and to reveal ecological risk is the cumulative distribution function called SSD (Species Sensitivity Distribution) [32,33]. The 5th percentile of this distribution (called the 5% Hazardous Concentration, HC5) has been used by USEPA [34], the RIVM Institute (The Netherlands) [35], and the European Commission [36] for deriving threshold concentrations that protect almost all species in a community. Research based on SSD was conducted for Ecological Risk Assessment (ERA) of several paddy insecticides and herbicides applying the Potentially Affected Fraction (PAF) as an index for the magnitude of ecological risk capable of reducing diversity [32,37,38].
The aim of this work is to investigate the relationship between chemical and ecotoxicological Lines of Evidence (LOEs) in the assessment of environmental impact due to the use of PPPs in Italian paddies, supported by a statistical approach, in a more realistic way than the traditional tabular evaluation. The findings will provide scientifically useful indications on the environmental compatibility of the use of pesticides.
The chemical and ecotoxicological analyses were carried out in conjunction with the sampling of species and habitats linked to the agroecosystems, in a larger project aimed at testing the measures envisaged by the National Action Plan for the sustainable use of PPPs for Natura 2000 sites and in protected areas, through the comparison of the values in biological and conventional fields. These measures indicate the biological method as the one most compatible with the conservation of biodiversity.

2. Materials and Methods

2.1. Area of Study

The study area is located in the Po Valley, between Piedmont and Lombardy (Italy), a vast agricultural territory in which, for irrigation management needs, rice cultivation is mainly monoculture. The cultivation consists of an uninterrupted succession of rice-field chambers interwoven with a large and articulated irrigation network functional for the distribution of water.
Seven rice paddies, belonging to four organic and three conventional farms, located in the Vercelli plain, were investigated from 2018 to 2019 (Figure 1).
The zone includes five protected areas and Natura 2000 sites, as detailed in Table 1.

2.2. Sampling Strategy

The sampling campaign took place in 2018–2019 and was strongly influenced by the irrigation level and weather conditions. In 2018, soil sampling was performed on the paddy embankment (em) and inside the paddy chamber (ch), when the growing season was suitable. In 2019, because of long flooding events, samples were collected only in the paddy embankments and the paddy field RB4 was replaced with RB6 (Table 1).
The sampling strategy for water and soil samples was divided into two phases for each year: one sample was collected at the beginning of the growing season, before the phytosanitary treatment (t0), and a second one after the phytosanitary treatments (t1), for a total of four sampling events in two years. For each phase, in presence of water flow, water samples were taken at entry (in) and exit (out) of the paddy field. In the absence of water flow, two water samples were taken near the channels of entry (in) and exit (out) from the paddy field. This sampling strategy allowed for verifying the possible contributions of contaminants already present in the irrigation water before entering into the paddy field chamber.
Water samples were kept in decontaminated glass bottles and stored in the dark at −20 °C until further analyses.
Four aliquots of soil from each paddy embankment were collected, one on each side of the chamber; the aliquots were then pooled, homogenized, sieved (2 mm), and stored at −20 °C.

2.3. Chemical Characterization

2.3.1. Soil Samples

Soil samples were dried and homogenized by grinding with an IKA® mill equipped with a beater blade. Total Carbon (TC) and Total Nitrogen (TN) were determined with a CHNS analyzer Vario Micro Cube Elementar and referred to as dry weight. Inorganic carbon was removed by progressive additions of HCl, then evaporated at ~50 °C. The analyzer performed a controlled combustion (~900 °C), then a catalytic oxidation (Chromium oxide), and finally a reduction by metallic copper. The CO2 and N2 developed were determined by a thermal conductivity detector after their gas chromatographic separation. Quality control was performed by daily Acetanilide analysis and by repeated measurements of two standard soils (Boden Standard AIVA Analysentechnik and Soil standard 2.1 Elementar).
For Ca, K, and Mg analyses, the sample was dried at 35 °C for 48 h and homogenized in an agate mortar. It was then subjected to a microwave acid digestion, using a 1:3 mixture of nitric and hydrochloric acid [39]. Instrumental determination was performed by Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-OES) [40], at the following wavelengths: 766.492 nm for K; 422.673 nm for Ca; 279.800 nm for Mg, referring concentrations to dry weight. The limit of quantification (LOQ) for these elements was 0.010%. Quality control was performed analyzing certified reference materials (sediment PACS-2 and soil SRM 2709).
The following PPPs were determined in dried and homogenized soil samples: Chlorpyrifos, Penconazole, Metalaxyl, Metrafenone, Pendimenthalin, Metolachlor, Alpha Endosulfan, Beta Endosulfan, Oxadiazon, Boscalid, Deltamethrin, Lambda Cyhalothrin, Oxifluorfen, Tebuconazole, Folpet, Dimethomorph, and Kresoxim methyl.
PPPs were extracted by Pressurized Fluid Extraction, dehydrated, and concentrated under a nitrogen flow, filtered, and analyzed by Gas Chromatography coupled with a triple-quadrupole Mass Spectrometer (GC/MS/MS) in MRM mode. Identification was based on the presence of at least two characteristic transitions for each analyte. Quantification was performed using perdeuterated internal injection standards. The LOQ of the method was 0.1 mg Kg−1 for each analyte. The method was developed by recovery tests on spiked samples; the quality control involved the analysis, at each sequence, of blanks, replicates, and recovery tests.

2.3.2. Water Samples

The determination of Ca, K, Mg, and S was performed on whole samples, after acidification with 2% nitric acid, by Atomic Emission Spectroscopy technique with inductively coupled plasma [40] at the same wavelengths as for soil samples. Limits of quantification were 4 mg L−1 for Calcium and 1 mg L−1 for Mg, K, and S. Copper determination was carried out by AAS technique with a graphite furnace [41] at a wavelength of 327.4 nm, with a quantification limit (LOQ) of 1 μg L−1. The quality control was performed using different certified reference materials: NIST 1643f for Ca, K, Mg and Cu, and SRM2709 for S.
The same PPPs of soil were also detected in water samples. A solution of deuterated internal standards was added to 200 mL of aqueous sample, and an extraction with dichloromethane was performed three times by separation funnel shaking. The extract was dehydrated with anhydrous sodium sulfate to the final volume by nitrogen blowing and then analyzed by GC/MS in SIM mode. Identification was performed by comparison of retention indices with those of standards and comparison of the relative abundances of the ions. Quality control was ensured by the analysis, for each batch of samples, of standards, method blanks, and recovery tests on spiked samples. The limit of quantification of the method (LOQ) for Folpet and Deltamethrin was 0.20 μg L−1 and for all other substances was 0.10 μg L−1.

2.4. Ecotoxicological Characterization

Ecotoxicological assays were started simultaneously on all test species, using the same sample aliquot of chemical analyses.
Different batteries of bioassays were set up in order to optimize the ecological representativeness of test species with respect to the environmental characteristics of the crops under study. The main features of the toxicity tests are detailed in Table 2.
All bioassays were performed according to standardized ISO methods in the ISPRA certified laboratories, according to UNI EN ISO 9001:2015.

2.5. Ecotoxicological Hazard Assessment

ISPRA and Public University of Ancona (Italy) have already developed several synthetic indices for the ecological risk assessment in marine environments [42,43], transposed in 2016 in the Technical Annex of Ministerial Decree No. 173/2016, which regulates dredging activities in relation to dumping.
The model, by means of a modular structure corresponding to several LOEs, combines, in a weighted way, the chemical, biological, and ecotoxicological characteristics of sediments. Regarding the ecotoxicological LOE, the Hazard Quotient (HQeco) integrates the results of the bioassays based not only on the biological measured effects, but also on the severity of the endpoint (growth, bioluminescence, survival, embryonic development, etc.), the ecological relevance of the tested environmental matrix (pore water, elutriate, whole sediment, etc.) as well as the type of exposure (chronic or acute) [42]. In particular, weighted criteria to elaborate results from standardized ecotoxicological bioassays are based on specific thresholds and weights assigned to each bioassay depending on the biological endpoint, tested matrix, time of exposure, and the possibility of hormetic responses. The cumulative Hazard Quotient referred to as battery bioassays (HQBattery) is obtained by the summation (Σ) of the weighted effects (Ew), i.e., the variations measured for each test compared to specific thresholds, corrected for the statistical significance of the difference (w), the biological importance of the endpoint, and exposure conditions (w2) [42]:
H Q B a t t e r y = K = 1 N E f f e c t w ( k ) · w 2
The HQBattery is normalized to a scale ranging from 0 to 10, where 1 is the battery threshold when all the bioassays exhibit an effect equal to own toxicity threshold, while 10 indicates that all the assays exhibit 100% of the effect. The HQBattery is then assigned to one of five classes of hazard, from Absent to Severe [43,44].
The flexible structure of HQ index and weighted criteria can be easily adapted to other applications, simply by varying the test species and weights of the different variables considered environmentally relevant.
In the case of rice fields, the toxicity thresholds and weights assigned to variables considered in the integrated index for estimating ecotoxicological Hazard Quotient (HQeco) are detailed in Table S1 (Supplementary Materials). The weight of the different types of endpoint varies according to the severity of the biological effect, with a maximum value for mortality (2.2) and a minimum value for root elongation (1.1). According to the same principle, acute toxicity has greater weight (1.0) than chronic toxicity (0.8).

2.6. Statistical Analyses

In order to compare the ecotoxicological results obtained between organic and conventional paddies, a specific t-test was applied for inhomogeneous variance at a confidence level of 95% (p = 0.05). Finally, with the aim to assess possible relationships between the measured ecotoxicological responses and the chemical characteristics detected in the different matrices, a multivariate analysis (PCA) was carried out on the raw data. Statistical elaborations were performed using R (version 4.2.1) and RStudio (version 2022.07.1) with the installed packages lme4, vegan, factoextra, and ggplot2.

3. Results and Discussion

3.1. Chemical Characteristics of Soil

In Table S2 (Supplementary Materials) the raw data referred to as organic content (TOC, TC), macro-elements, and inorganic pesticides (Cu and S) sought for characterization of samples are detailed. The percentages of Total Nitrogen (TN), approximately in a range of 0.1–0.3, reflect a type of soil “well endowed”, according to the classification of Giardini [45], while the C/N ratios are generally > 11. The high percentage of organic component leads to nitrogen stabilization, making it less bioavailable. The ratios did not change substantially between organic and conventional paddies, with median values of 11.22 and 10.58, respectively. The amount of organic matter did not show significant differences between the organic and conventional fields (p = 0.335) throughout the study period, whereas conventional paddies resulted in significantly higher contents of Cu and S in 2018 (p = 0.005 and 0.004, respectively). For both metals, the highest values were always found in the conventionally treated rice fields, especially in RT5. The relatively high values measured could be associated with the large use of copper sulphate, as classic verdigris in the entire area.
Table S3 reports the detail of results for the pesticides for which at least one sample showed quantifiable concentrations (>LOQ). In particular, the PPP concentrations were <LOQ in all samples from organic crops, while a few cases of barely quantifiable concentrations of Oxadiazon, Oxyfluorfen, and Pendimethalin were found in samples from conventional paddies. The low number of quantifiable results does not allow for evaluating differences between the various types of samples from rice fields (chamber, embankment, soil inlet, or outlet) or even just between before (t0) and after treatment (t1).
Despite the many active ingredients sought, most measures were below the LOQ. The finding demonstrates that, despite frequent treatments with PPPs over the years, the irrigation cycles in rice fields favor soil leaching, which, therefore, retains little residue of contaminants in the chamber, promoting their transfer to the water.

3.2. Chemical Characteristics of Water

Table 3 shows the concentrations of macro elements (Ca, K, Mg) and inorganic and organic pesticides measured in the water samples for which at least one value was quantified, except for sample RT4_in_t0 because its amount was not sufficient to carry out the analyses.
The concentrations of the macro elements indicate that Ca is the dominant element (4.27–44.25 mg L−1), followed by Mg (up to 15.01 mg L−1) and K (up to 33.93 mg L−1) in similar quantities. No significant statistical differences were found between organic and conventional fields about the content of both macro elements and inorganic pesticides. However, a general trend is noticeable towards a lower Ca, Mg, and K content in organic rice fields, especially in 2019 (p = 0.076 for Ca; p = 0.097 for K; p = 0.099 for Mg).
No differences related to either the sampling point (in or out) or the two sampling campaigns following treatment were noticed. The concentrations of Cu and S measured in water are in very low ranges, especially for copper, as are the values found in soil samples from the same fields. This is likely due to the limitations on the use of copper sulphate in rice fields from 2012, according to Regulation 1107/2009/EC [29].
Oxadiazon was detected in almost all samples analyzed in concentrations up to 47.6 µg L−1 (RT5_out_t0), showing widespread contamination, although the organic paddies showed concentrations generally lower (0.45 ± 0.52 µg L−1) than conventional ones (2.36 ± 9.44 µg L−1), but still present, even before the first treatment of the culture cycle (t0). The overall variability among the rice fields is such that no statistically significant difference was identified (p = 0.343), considering both the individual years of investigation and the overall comparison among organic and conventional crops.
Other pesticides, such as Metolachlor and Oxyfluorfen, were detected in very low concentrations in sporadic cases (Table 3), whereas the other chemicals analyzed were all below the LOQ of the specific method, except for lambda-cyhalothrin in RB6_in_t0 (0.17 µg L−1).

3.3. Ecotoxicological Effects of Soil

The soil samples of conventional paddies showed significant toxicity for root elongation (Table 4) and for the bioluminescence inhibition assays on eluate (Table 5), in particular for RT5 and RT4 rice fields, making the difference between organic and conventional crops statistically significant for L. sativum (p = 0.047) and S. saccharatum (p = 0.001).
These results agree with the chemical data measured in conventional paddies regarding traces of Oxadiazon (RT1, RT4 and RT5) and Oxyfluorfen (RT1 and RT5). These PPPs are two herbicides, and the Oxadiazon could be used in Italy until 30 June 2020, whereas Oxyfluorfen is authorized until 2024 [46].

3.4. Ecotoxicological Effects of Water

Table 6 reports the results of bioassays expressed as bioluminescence inhibition for A. fischeri, average immobilization for D. magna, growth rate inhibition for R. subcapitata, and leaf growth inhibition for S. polyrhiza (only for 2018), respectively.
In 2018, with regard to A. fischeri, toxicity effects were found only in RB1, with particular reference to water leaving the organic rice chamber; regarding the 2019 campaign, most of the samples showed significant inhibition of bioluminescence, with no statistical differences between conventional and organic crops (p > 0.05).
Despite being the least sensitive among the organisms of the selected battery, D. magna showed some important toxicity effects on both conventional and organic paddies for both monitoring campaigns. It should be noted that the toxicity data on crustaceans are particularly relevant, since, from an ecotoxicological point of view, they are organisms characterized by relatively low sensitivity but, at the same time, high ecological value for the environment in question, as they live in the aquatic environments of rice paddies that adopt organic cultivation methods [47].
The green alga R. subcapitata was the most sensitive species in the battery of bioassays used: for all monitoring campaigns, water samples showed high toxicity, but without particular differences between inlet and outlet water samples, between conventional and organic rice fields. These ecotoxicological effects could be associated with the presence of at least three herbicides: Metolachlor, Oxadiazon, and Oxifluorfen. In particular, Oxadiazon is found in all organic and treated fields (Table 3), although these substances have not been officially used in organic rice fields. Oxadiazon is known to have important toxic effects on R. subcapitata [48]. The high toxicity could, therefore, be linked to the simultaneous presence of these herbicides in the water, considering other possible synergistic effects with other contaminants not researched in this study, but found in the past in paddy field waters in Piedmont (https://www.snpambiente.it/2017/07/11/monitoraggio-dei-fitofarmaci-delle-acque-piemontesi-nelle-risaie/, accessed on 30 November 2022).
Except for sample taken during t0 in paddy RB5, for S. polyrrhiza, no inhibition values were revealed, but a general tendency towards biostimulation, probably due to the presence of nutrients in the water was observed.

3.5. Ecotoxicological Hazard Index and Statistical Analysis

The results of bioassay applied to embankment, chamber soil, and water samples were elaborated by the ecotoxicological Hazard Index, producing the hazard levels shown in Table 7.
With regard to the soil samples, no ecotoxicological hazard was detected (HQ < 1), due to a general absence of toxicity, according with the low presence of pesticides (Table S3). Therefore, neither differences between organic fields and their corresponding conventional fields, nor differences between the two survey campaigns, can be highlighted. This is in perfect agreement with the general absence of pesticide residues in soil and sediment samples.
The situation is quite different for water samples. In general, important toxic effects were found on several organisms and for several endpoints, such as to determine an ecotoxicological hazard in all samples (HQ > 1) without exception, sometimes even “major” or “severe” (Table 7). The ecotoxicological hazard affects both conventional and organic paddy waters, with the highest hazard measured in the RB1 field, especially in 2019.
This ecotoxicological picture is also in agreement with the findings of the chemical analysis, with particular reference to the ubiquitous presence of Oxadiazon and Oxyfluorfen.
Figure 2 shows the plot of a Principal Component Analysis (PCA) applied to raw data that explains 67.7% of total variance. Almost all samples appear to be distributed along a gradient oriented with respect to the first component (PCA1), whose main variable contributing to this distribution is the parameter Ca (with a higher factor loading of −0.924) and secondly Mg and S in the same way (factor loading of −0.241 and −0.251, respectively), in accordance with the significantly higher content of these macro-elements in organic crops. With respect to the second component (PCA2), the sample RT5_t0_out_2019 is clearly distinguished in the factorial space due to the highest Oxadiazon concentration, associated to the highest factor loading (−0.940).
This herbicide is, therefore, the substance that most influences the characteristics of paddy water; its use is regularly reported by conventional paddies, while it has not been used by organic ones.

4. Conclusions

In conclusion, the study highlights the importance of assessing the impact of pesticides directly in paddy fields using an integrated approach, because the chemical or ecotoxicological approach alone is unable to consider numerous environmental factors that can occur as the photodecomposition responsible for rapid degradation [48], the volatilization favored by high temperatures [49,50], the hydrolysis influenced by the flooded conditions and pH of the paddy water [51,52,53], and microbial degradation [54].
Chemical analyses of paddy waters revealed the presence of several PPPs, such as Oxadiazon and Metolachlor, in conventional fields. Relatively high levels of Oxadiazon were also detected in the biological fields, demonstrating its ubiquitous presence in the paddy rice area examined in this study. Its homogeneous concentration in the fields might be derived by the cascade circulation pattern of irrigation water in the area. Indeed, it is quite likely that the Oxadiazon applied to conventional rice paddies moves to organic from ones through the circulation of water, determining a relevant ecotoxicological response in both conventional and in organic fields. Nevertheless, it is reasonable to assume that the pesticide level decreases rapidly over time in the runoff water, accounting for volatilization, degradation, leaching to groundwater, and sorption to soil [49].
The application of a multidisciplinary approach has allowed for the identification of the main cause of contamination in rice fields in relation to ecotoxicological effects, i.e., Oxadiazon freely circulating in waters and the definition of a suitable solution strategy: the realization of organic districts with water circulation isolated from conventional fields should limit the transfer of the PPPs through the fields, enhancing the benefits of organic farming practices.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14244136/s1, Table S1: Weights assigned to the variables considered in the integrated index for estimating ecotoxicological Hazard Quotient (HQeco); Table S2: Chemical characterization of soil samples from paddies; Table S3: Concentration of PPPs measured in soil samples of paddies.

Author Contributions

Conceptualization, F.O., A.T., A.P. and C.M.; Formal analysis, F.O.; Funding acquisition, S.D.; Investigation, A.T., A.P., G.S., M.T.B., E.R. and C.M.; Methodology, F.O., A.T., A.P. and C.M.; Project administration, F.O., C.M., S.D. and E.R.; Resources, A.T., A.P., G.S., M.T.B. and C.M.; Supervision, F.O., C.M. and S.D.; Validation, F.O., A.T., A.P. and C.M.; Visualization, F.O.; Writing—original draft, F.O.; Writing—review and editing, A.T. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Ecological Transition from 2015 to 2020 as part of the project for the experimentation of the Measures of the National Action Plan (PAN) for the sustainable use of pesticides for the protection of biodiversity, especially in “Natura 2000” sites, in accordance to 2009/128/CE Directive, transposed in Italy by legislative decree n. 150/2012.

Data Availability Statement

Data supporting reported results can be found in: La sperimentazione dell’efficiacia delle Misure del Piano d’Azione nazionale per l’uso sostenibile dei prodotti fitosanitari (PAN) per la tutela della biodiversità, ISPRA, Vol. 330/2020, 392 pp. (https://www.isprambiente.gov.it/it/pubblicazioni/rapporti/rapporto-330-2020-web-12-01-2021-rid-1.pdf; accessed on 5 December 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tejada, A.W. Pesticide residues in foods and the environment as a consequence of crop protection. Phillip. J. Agric. 1995, 78, 63–79. [Google Scholar]
  2. Lam, P.V.; Soon, B.H. Influence of broad-spectrum insecticides on main population of main predators in the rice field. Plant Prot. Bull. Vietnam 1994, 6, 7–12. [Google Scholar]
  3. Beketov, M.; Liess, M. Delayed effects of the neonicotinoid insecticide thiacloprid on seven freshwater arthropods. Environ. Toxicol. Chem. 2008, 27, 461–470. [Google Scholar] [CrossRef]
  4. Beketov, M.A.; Schäfer, R.B.; Marwitz, A.; Paschke, A.; Liess, M. Long-term stream invertebrate community alterations induced by the insecticide thiacloprid: Effect concentrations and recovery dynamics. Sci. Total Environ. 2008, 405, 96–108. [Google Scholar] [CrossRef]
  5. Botias, C.; David, A.; Hill, E.M.; Goulson, D. Quantifying exposure of wild bumblebees to mixtures of agrochemicals in agricultural and urban landscapes. Environ. Pollut. 2017, 222, 73–82. [Google Scholar] [CrossRef]
  6. Brühl, C.A.; Pieper, S.; Weber, B. Amphibians at risk? Susceptibility of terrestrial amphibian life stages to pesticides. Environ. Toxicol. Chem. 2011, 30, 2465–2472. [Google Scholar] [CrossRef]
  7. Brühl, C.A.; Schmidt, T.; Pieper, S.; Alscher, A. Terrestrial pesticide exposure of amphibians: An underestimated cause of global decline? Sci. Rep. 2013, 3, 1135. [Google Scholar] [CrossRef] [Green Version]
  8. Desnev, N.; Decourtye, A.; Delpuech, J.M. The sublethal effects of pesticides on beneficial arthropods. Annu. Rev. Entomol. 2007, 52, 81–106. [Google Scholar] [CrossRef]
  9. Diaz, S.; Settele, J.; Brondizio, E.S.; Ngo, H.T.; Agard, J.; Arneth, A.; Balvanera, P.; Brauman, K.A.; Butchart, S.H.M.; Chan, K.M.A.; et al. Pervasive human-driven decline of life on Earth points to the need for transformative change. Science 2019, 366, 6471. [Google Scholar] [CrossRef] [Green Version]
  10. EFSA (European Food safety Authority). Toxicological data analysis to support grouping of pesticide active substances for cumulative risk assessment of effects on liver, on the nervous system and on reproduction and development. EFSA Support. Publ. 2013, 10, 392E. [Google Scholar] [CrossRef] [Green Version]
  11. Howe, C.M.; Berrill, M.; Pauli, B.D.; Helbing, C.C.; Werry, K.; Veldhoen, N. Toxicity of Glifosate-based pesticides to four North American frog species. Environ. Toxicol. Chem. 2004, 23, 1928–1938. [Google Scholar] [CrossRef] [PubMed]
  12. Mohr, S.; Berghahn, R.; Schmiediche, R.; Hübner, V.; Loth, S. Macroinvertebrate community response to repeated short-term pulses of the insecticide imidacloprid. Aquat. Toxicol. 2012, 110–111, 25–36. [Google Scholar] [CrossRef] [PubMed]
  13. Slooff, W.; Canton, J.H.; Hermens, J.L.M. Comparison of the susceptibility of 22 freshwater species to 15 chemical compounds. I. (sub)acute toxicity tests. Aquat. Toxicol. 1983, 4, 113–128. [Google Scholar] [CrossRef]
  14. Stahlschmidt, P.; Brühl, C.A. Bats at risk? Bat activity and insecticide residue analysis of food items in an apple orchard. Environ. Toxicol. Chem. 2012, 31, 1556–1563. [Google Scholar] [CrossRef]
  15. Van Wijngaarden, R.P.A.; Brock, T.C.M.; Van Den Brink, P.J. Threshold levels for effects of insecticides in freshwater ecosystems: A review. Ecotoxicology 2005, 14, 355–380. [Google Scholar] [CrossRef]
  16. Van Swaay, C.A.M.; Van Strien, A.J.; Harpke, A.; Fontaine, B.; Stefanescu, C.; Roy, D.; Maes, D.; Kühn, E.; Õunap, E.; Regan, E.; et al. The European Grassland Butterfly Indicator: 1990–2011; European Environment Agency, Technical Report No. 11/2013; De Vlinderstichting: Wageningen, The Netherlands, 2013; p. 36. [Google Scholar]
  17. Santovito, A.; Audisio, M.; Bonelli, S. A micronucleus assay detects genotoxic effects of herbicide exposure in a protected butterfly species. Ecotoxicology 2020, 29, 1390–1398. [Google Scholar] [CrossRef]
  18. Zacharia, J.T. Ecological Effects of Pesticides. In Pesticides in the Modern World—Risks and Benefits; Stoytcheva, M., Ed.; IntechOpen: London, UK, 2011. [Google Scholar] [CrossRef] [Green Version]
  19. Liess, M.; Von Der Ohe, P.C. Analyzing effects of pesticides on invertebrate communities in streams. Environ. Toxicol. Chem. 2009, 24, 954–965. [Google Scholar] [CrossRef]
  20. Burrows, L.A.; Edwards, C.A. The use of integrated soil microcosm to predict effects of pesticides on soil ecosystems. Eur. J. Soil Biol. 2002, 38, 245–249. [Google Scholar] [CrossRef]
  21. Kovach, J.; Petzoldt, C.; Degni, J.; Tette, J. A method to measure the environmental impact of pesticides. New York’s Food Life Sci. Bull. 1992, 139, 1–8. [Google Scholar]
  22. Jayasiri, M.J.G.C.N.; Yadav, S.; Propper, C.R.; Kumar, V.; Dayawansa, N.D.K.; Singleton, G.R. Assessing potential environmental impacts of pesticides usage in paddy ecosystems: A case study in the Deduru Oya River basin, Sri Lanka. Environ. Toxicol Chem. 2022, 41, 343–355. [Google Scholar] [CrossRef]
  23. Sánchez-Bayo, F.; Goka, K. Evaluation of suitable endpoints for assessing the impacts of toxicants at the community level. Ecotoxicology 2012, 21, 667–680. [Google Scholar] [CrossRef] [PubMed]
  24. Directive 2009/128/EC of the European Parliament and of the Council of 21 October 2009 Establishing a Framework for Community Action to Achieve the Sustainable Use of Pesticides. Available online: http://data.europa.eu/eli/dir/2009/128/oj (accessed on 5 December 2022).
  25. Legislative Decree 14 August 2012, n. 150. Attuazione Della Direttiva 2009/128/CE che Istituisce un Quadro per L’azione Comunitaria ai fini Dell’utilizzo Sostenibile dei Pesticidi. Supplemento Ordinario n. 177/L, Gazzetta Ufficiale n. 202 del 30/08/2012. Available online: https://www.mite.gov.it/normative/dlgs-14-agosto-2012-n-150-attuazione-della-direttiva-2009128ce-che-istituisce-un-quadro (accessed on 5 December 2022).
  26. Council Directive 92/43/EEC of 21 May 1992 on the Conservation of Natural Habitats and of Wild Fauna and Flora. Available online: https://www.legislation.gov.uk/eudr/1992/43/contents (accessed on 5 December 2022).
  27. Directive 2009/147/EC of the European Parliament and of the Council of 30 November 2009 on the Conservation of Wild Birds. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2010:020:0007:0025:EN:PDF (accessed on 5 December 2022).
  28. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:32000L0060 (accessed on 5 December 2022).
  29. Regulation (EC) No 1107/2009 of the European Parliament and of the Council of 21 October 2009 Concerning the Placing of Plant Protection Products on the Market and Repealing Council Directives 79/117/EEC and 91/414/EEC. Available online: https://www.legislation.gov.uk/eur/2009/1107/contents (accessed on 28 November 2022).
  30. Environmental Protection Agency (EPA). Pesticides; Data Requirements for Conventional Chemicals. Federal Register, Vol. 72, No. 207/October 26, 2007/Rules and Regulations. Available online: https://www.federalregister.gov/documents/2007/10/26/E7-20826/pesticides-data-requirements-for-conventional-chemicals (accessed on 28 November 2022).
  31. Japan’s Ministry Of Environment. Setting of Pesticides Registration Criteria Regarding the Effect on Living Environmental Animals and Plants (Initial Report); Godochosha No. 5; Central Environment Council: Tokyo, Japan, 2019. (In Japanese)
  32. Nagai, T. Ecological effect assessment by species sensitivity distribution for 38 pesticides with various modes of action. J. Pestic. Sci. 2021, 46, 366–372. [Google Scholar] [CrossRef] [PubMed]
  33. Posthuma, L.; Suter, G.W.; Traas, T.P. (Eds.) Species Sensitivity Distributions in Ecotoxicology (Environmental and Ecological Risk Assessment); Lewis Publisher: Boca Raton, FL, USA, 2001; p. 587. [Google Scholar] [CrossRef]
  34. Environmental Protection Agency (EPA). Revised Deletion Process for the Site-Specific Recalculation Procedure for Aquatic Life Criteria; Environmental Protection Agency: Washington, DC, USA, 2013.
  35. National Institute for Public Health and the Environment (RIVM). Guidance for the Derivation of Environmental Risk Limits. Part 3. Derivation of ERLs for Freshwater and Marine Water; Version 1.0; National Institute for Public Health and the Environment: De Bilt, The Netherlands, 2015; p. 60.
  36. Scientific Committee on Health, Environmental and Emerging Risks (SCHEER). Scientific Advice on Guidance Document n. 27: Technical Guidance for Deriving Environmental Quality Standards; European Commission, Directorate General for Health and Food Safety: Brussels, Belgium, 2018; p. 31. [Google Scholar] [CrossRef]
  37. Nagai, T.; Yokoyama, A. Comparison of ecological risks of insecticides for nursery-box application using species sensitivity distribution. J. Pestic. Sci. 2012, 37, 233–239. [Google Scholar] [CrossRef] [Green Version]
  38. Nagai, T. Studies on ecological risk assessment of pesticide using species sensitivity distribution. J. Pestic. Sci. 2017, 42, 124–131. [Google Scholar] [CrossRef] [Green Version]
  39. Environmental Protection Agency (EPA). Method 3051A: Microwave assisted acid digestion of sediments, sludges, soils, and oils. Z. Für Anal. Chem. 2007, 111, 362–366. [Google Scholar]
  40. Environmental Protection Agency (EPA). Method 6010C (SW-846): Inductively Coupled Plasma-Atomic Emission Spectrometry. 2007. Available online: https://19january2017snapshot.epa.gov/sites/production/files/2015-07/documents/epa-6010c.pdf (accessed on 5 December 2022).
  41. Environmental Protection Agency (EPA). Method 7010: Graphite Furnace Atomic Absorption Spectrophotometry. 2007. Available online: https://www.epa.gov/sites/default/files/2015-12/documents/7010.pdf (accessed on 5 December 2022).
  42. Piva, F.; Ciaprini, F.; Onorati, F.; Benedetti, M.; Fattorini, D.; Ausili, A.; Regoli, F. Assessing sediment hazard through a weight of evidence approach with bioindicator organisms: A practical model to elaborate data from sediment chemistry, bioavailability, biomarkers and ecotoxicological bioassays. Chemosphere 2011, 83, 475–485. [Google Scholar] [CrossRef] [PubMed]
  43. Benedetti, M.; Ciaprini, F.; Piva, F.; Onorati, F.; Fattorini, D.; Notti, A.; Ausili, A.; Regoli, F. A multidisciplinary weight of evidence approach for classifying polluted sediments: Integrating sediment chemistry, bioavailability, biomarkers responses and bioassays. Environ. Int. 2012, 38, 17–28. [Google Scholar] [CrossRef]
  44. Regoli, F.; D’Errico, G.; Nardi, A.; Mezzelani, M.; Fattorini, D.; Benedetti, M.; Di Carlo, M.; Pellegrini, D.; Gorbi, S. Application of a weight of evidence approach for monitoring complex environmental scenarios: The case-study of off-shore platforms. Front. Mar. Sci. 2019, 6, 377. [Google Scholar] [CrossRef] [Green Version]
  45. Giardini, L. Agronomia Generale, 3rd ed.; Patron: Bologna, Italy, 1986; Available online: https://www.hoepli.it/libro/agronomia-generale/9788855526388.html (accessed on 5 December 2022).
  46. Ministry of Health. Withdrawal of Plant Protection Products Containing the Active Substance Oxadiazon Following the Expiry of Its Period of Community Approval under Regulation (EU) 1107/2009. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2009:309:0001:0050:en:PDF (accessed on 5 December 2022). (In Italian).
  47. Martínez-Jerónimo, F.; Villaseñor, R.; Rios, G.; Espinosa, F. Effect of food type and concentration on the survival, longevity, and Reproduction of Daphnia magna. Hydrobiologia 1994, 287, 207–214. [Google Scholar] [CrossRef]
  48. Deblois, C.P.; Qiu, B.; Juneau, P. Effect of Herbicides (Diuron and Oxadiazon) on Photosynthetic Energy Dissipation Processes of Different Species of Cyanobacteria and Two Green Algae. In Photosynthesis. Energy from the Sun; Springer: Dordrecht, The Netherlands, 2008; ISBN 978-1-4020-6707-5. [Google Scholar] [CrossRef]
  49. Comoretto, L.; Arfib, B.; Talva, R.; Chauvelon, P.; Pichaud, M.; Chiron, S.; Höhener, P. Runoff of pesticides from rice fields in the Ile de Camargue (Rhone River delta, France): Field study and modelling. Environ. Pollut. 2008, 151, 486–493. [Google Scholar] [CrossRef]
  50. Tanabe, S.; Ramesh, A.; Sakashita, D.; Iwata, H.; Mohan, D.; Subramaniam, A.N.; Tatsukawa, R. Fate of HCH in tropical paddy field: Application test in South India. Int. J. Environ. Anal. Chem. 1992, 45, 45–53. [Google Scholar] [CrossRef]
  51. Sethunathan, N.; Yoshida, T. Fate of diazinon in submerged soil: Accumulation of hydrolysis product. J. Agric. Food Chem. 1969, 17, 1192–1195. [Google Scholar] [CrossRef]
  52. Seiber, J.N.; Heinrichs, E.A.; Aquino, G.B.; Valencia, S.L.; Andrade, P.; Argente, A.M. Residues of Carbofuran Applied as a Systemic Insecticide in Irrigated Wetland Rice: Implications for Insect Control; IRRI Research Paper Series No. 17; FAO: Rome, Italy, 1978; pp. 1–28. [Google Scholar]
  53. Varca, L.M.; Magallona, E.D. Dissipation and degradation of DDT and DDE in Philippine soil under field conditions. J. Environ. Sci. Health 2008, 29, 25–35. [Google Scholar] [CrossRef]
  54. Venkateswarlu, K.; Sethunathan, N. Degradation of carbofuran by Azospirillum lipoferum and Streptomyces spp. isolated from flooded alluvial soil. Bull. Environ. Contam. Toxicol. 1984, 33, 556–560. [Google Scholar] [CrossRef]
Figure 1. Localization of crops involved in the research project.
Figure 1. Localization of crops involved in the research project.
Water 14 04136 g001
Figure 2. Plot of PCA applied to chemical and ecotoxicological characteristics of all water samples from rice paddies (t0 and t1 labels refer to the campaigns before and after PPPs treatment, respectively; in and out indicate the samples taken from streams entering and exiting the paddy water chamber, respectively, in 2018 and 2019. Only parameters (blue arrows) that contribute more than 15% in the variability of the samples are shown.
Figure 2. Plot of PCA applied to chemical and ecotoxicological characteristics of all water samples from rice paddies (t0 and t1 labels refer to the campaigns before and after PPPs treatment, respectively; in and out indicate the samples taken from streams entering and exiting the paddy water chamber, respectively, in 2018 and 2019. Only parameters (blue arrows) that contribute more than 15% in the variability of the samples are shown.
Water 14 04136 g002
Table 1. Identification code, the municipality, and type of agronomic management (RT = conventional; RB = organic) of the seven rice paddies considered in the research.
Table 1. Identification code, the municipality, and type of agronomic management (RT = conventional; RB = organic) of the seven rice paddies considered in the research.
CodeMunicipalityProtected Areas of the Natura 2000 NetworkAgronomic
Management
Year
RT1VillarboitZSC/ZPS IT1120014 Druma River marshConventional2018–2019
RT4RovasendaSIC IT1120026 Station of Isoetes malinvernianaConventional2018–2019
RT5Crescentino-Conventional2018–2019
RB1RovasendaZSC IT1120004 Baraggia of Rovasenda/EUAP0349 Natural Reserve of BaraggeOrganic2018–2019
RB4VillarboitZSC/ZPS IT1120014 Druma River marshOrganic2018
RB5Crescentino-Organic2018–2019
RB6CrovaZPS IT1120021 Paddies of VercelliOrganic2019
Table 2. Battery of bioassays used for ecotoxicological assessment of soil and water samples of rice paddies.
Table 2. Battery of bioassays used for ecotoxicological assessment of soil and water samples of rice paddies.
Sampling PointSpeciesCommon NameEnvironmental
Matrix
End-PointExpositionMethod
Embankment and paddy chamber soilLepidium sativumWatercressSoilGermination
and root
elongation
72 hISO 18763:2016
Sinapis albaMustard
Sorghum saccharatumSorghum
Aliivibrio fischeriBacteriumEluateBiolumin.30 minISO 11348-3:2019
WaterRaphidocelis subcapitataGreen algaeWaterGrowth rate72 hISO 8692:2012
Daphnia magnaWater fleaWaterImmobilization24 hISO 6341:2013
Spirodela polyrhizaDuckmeatWaterLeaves growth72 hISO 20227:2017
Aliivibrio fischeriBacteriumWaterBiolumin.30 minISO 11348-3:2019
Table 3. Chemical characterization of water samples from rice paddies (RB = organic cultures; RT = conventional cultures). In 2019 RB4 was substitute with RB6 (in = paddy field chamber entrance; out = paddy field chamber exit).
Table 3. Chemical characterization of water samples from rice paddies (RB = organic cultures; RT = conventional cultures). In 2019 RB4 was substitute with RB6 (in = paddy field chamber entrance; out = paddy field chamber exit).
SampleCa
(mg L−1)
K
(mg L−1)
Mg
(mg L−1)
Cu
(μg L−1)
S
(mg L−1)
Λ-Cyalothrin
(μg L−1)
Metolachlor
(μg L−1)
Oxadiazon
(μg L−1)
Oxyfluorfen
(μg L−1)
2018
RB1int04.811.071.123.601.79<0.1<0.10.15<0.1
t18.862.782.022.852.31<0.1<0.10.14<0.1
outt09.4419.793.088.072.44<0.1<0.10.13<0.1
t112.76<1.003.57<1.00<1.00<0.1<0.10.10<0.1
RB4int09.225.572.251.852.7<0.1<0.11.68<0.1
t111.004.003.11.941.17<0.1<0.10.67<0.1
outt013.853.813.316.863.13<0.1<0.10.21<0.1
t19.891.902.815.67<1.00<0.1<0.10.52<0.1
RB5int027.396.928.004.9210.08<0.1<0.11.72<0.1
t134.05<1.003.572.105.06<0.1<0.11.09<0.1
outt036.598.248.222.737.01<0.1<0.10.52<0.1
t136.591.539.993.0319.58<0.1<0.10.23<0.1
RT1int06.292.091.342.672.27<0.1<0.11.750.23
t110.92.792.565.582.44<0.1<0.10.61<0.1
outt05.632.401.095.563.12<0.1<0.11.240.12
t19.5233.933.056.691.50<0.1<0.10.35<0.1
RT4int0n.d.n.d.n.d.n.d.n.d.<0.1<0.10.13<0.1
t15.60<1.001.226.731.63<0.1<0.10.17<0.1
outt04.271.96<1.003.822.43<0.1<0.10.15<0.1
t16.141.311.857.22<1.00<0.1<0.10.42<0.1
RT5int042.663.588.722.1210.42<0.1<0.10.82<0.1
t130.733.8515.01<1.009.25<0.10.130.22<0.1
outt025.1316.96.684.569.31<0.1<0.10.26<0.1
t130.163.136.892.596.92<0.10.470.35<0.1
2019
RB1int08.026.41.741.031.88<0.1<0.11.56<0.1
t19.582.031.931.161.61<0.1<0.10.11<0.1
outt012.1913.353.192.752.81<0.1<0.1<0.1<0.1
t18.952.001.801.291.31<0.1<0.10.17<0.1
RB5int039.144.579.541.9216.36<0.1<0.10.26<0.1
t144.255.9911.32.3410.08<0.1<0.10.18<0.1
outt042.1716.748.112.237.41<0.1<0.1<0.1<0.1
t131.536.289.95<1.004.32<0.1<0.10.23<0.1
RB6int015.261.713.451.065.780.17<0.10.20<0.1
t117.701.084.73<1.004.85<0.1<0.1<0.1<0.1
outt013.881.983.11<1.006.61<0.1<0.10.43<0.1
t117.46<1.003.98<1.003.81<0.1<0.10.46<0.1
RT1int06.051.891.35<1.002.34<0.1<0.10.13<0.1
t16.582.651.371.05<1.00<0.1<0.10.37<0.1
outt05.433.321.116.863.2<0.1<0.1<0.1<0.1
t112.15<1.002.341.14<1.00<0.1<0.10.40<0.1
RT4int06.47<1.001.04<1.001.90<0.1<0.10.17<0.1
t16.371.401.411.232.06<0.1<0.1<0.1<0.1
outt05.111.571.062.633.01<0.1<0.1<0.1<0.1
t16.111.111.111.442.11<0.1<0.10.14<0.1
RT5int029.793.478.471.3015.28<0.1<0.10.29<0.1
t116.982.994.67<1.006.17<0.1<0.10.29<0.1
outt027.454.554.732.2412.97<0.1<0.147.6<0.1
t124.075.617.692.477.64<0.1<0.10.71<0.1
Table 4. Results of phytotoxicity tests with Lepidium sativum (Ls), Sinapis alba (Sa), Sorghum saccharatum (Ss) on samples of soil collected during campaigns t0 and t1 (in bold are shown the significative effects above the toxicity threshold). Negative sign indicates biostimulation of root growth in comparison to control (Ch = rise chamber; em = embankment).
Table 4. Results of phytotoxicity tests with Lepidium sativum (Ls), Sinapis alba (Sa), Sorghum saccharatum (Ss) on samples of soil collected during campaigns t0 and t1 (in bold are shown the significative effects above the toxicity threshold). Negative sign indicates biostimulation of root growth in comparison to control (Ch = rise chamber; em = embankment).
2018LsSaSs2018LsSaSs
Conventional(%)±sd(%)±sd%±sdOrganic(%)±sd(%)±sd(%)±sd
RT1cht017.817.611.842.75−48.627.88t0emRB430.4812.2421.809.8619.477.04
t138.204.8824.537.834.6612.09t126.9417.3715.5114.6221.0415.51
RT4emt025.8413.3636.2714.6340.813.88t0emRB118.1114.57−6.7910.25−12.4322.14
t1n.d.n.d.n.d.n.d.n.d.n.d.t123.672.34−1.5011.38−11.9713.09
cht030.9115.1829.9410.4235.8415.39t0ch−4.058.693.1315.74−20.4524.46
RT5emt194.340.5882.903.2532.6212.63t0emRB524.4110.997.698.2723.6014.84
cht025.156.458.286.3815.961.02
t159.719.3553.186.4824.627.46t157.6614.9723.596.3319.5723.58
2019LsSaSs2019LsSaSs
RT1emt024.185.9241.665.3317.5212.06t0emRB620.263.9620.7116.1414.022.83
t149.063.3759.856.0224.375.53t16.078.599.0411.63−10.855.62
RT4emt022.405.5221.142.961.6014.33t0emRB128.562.60−8.2510.645.543.60
t118.192.7934.9910.248.758.11t132.262.706.9720.025.905.22
RT5emt030.5510.4223.790.1319.442.39t0emRB521.1116.3010.6213.7614.0811.87
t172.784.8373.082.1858.795.15t115.859.8311.1910.0413.478.86
Table 5. Results of bioassays with Aliivibrio fischeri on samples of soil eluate collected during campaigns t0 and t1 (in bold are shown the effects above the toxicity threshold). Negative sign indicates biostimulation in comparison to control (Ch = rise chamber; em = embankment).
Table 5. Results of bioassays with Aliivibrio fischeri on samples of soil eluate collected during campaigns t0 and t1 (in bold are shown the effects above the toxicity threshold). Negative sign indicates biostimulation in comparison to control (Ch = rise chamber; em = embankment).
2018
Conventional crops%±sdOrganic crops(%)±sd
RT1cht0−15.607.92RB4emt021.509.50
t131.6720.65t1−11.3218.49
RT4emt0−3.5125.78RB1emt0−0.2415.68
t1−28.7036.55
cht0−11.0421.68cht025.3214.53
RT5emt0−18.6519.62RB5cht010.887.90
cht025.2614.42emt0−24.7739.92
t11.1621.12t12.694.60
2019
RT1emt010.0018.21RB6emt0−90.59107.65
t1−11.2832.56t1−11.8413.97
RT4emt030.5733.73RB1emt0−26.3358.08
t110.5614.81t1−0.3311.21
RT5emt0−24.1733.99RB5emt0−3.238.73
t1−55.4878.87t1−18.8911.82
Table 6. Bioassay results expressed as percent effect (%) with standard deviation (sd) on paddy water samples collected during campaigns t0 and t1 before and after PPPs treatment; in and out labels refer to samples taken from streams entering and exiting the paddy water chamber; boldface indicates inhibition effects above the toxicity thresholds; negative sign indicates biostimulation in comparison to the control.
Table 6. Bioassay results expressed as percent effect (%) with standard deviation (sd) on paddy water samples collected during campaigns t0 and t1 before and after PPPs treatment; in and out labels refer to samples taken from streams entering and exiting the paddy water chamber; boldface indicates inhibition effects above the toxicity thresholds; negative sign indicates biostimulation in comparison to the control.
Aliivibrio fischeri
Conventional20182019Organic20182019
(%)ds(%)ds(%)ds(%)ds
RT1int0−11.18.4−4.0010.00t0inRB4
(2018)
RB6
(2019)
−6.56.896.1310.24
out−8.97.495.9910.60out14.413.3−7.367.97
int1−5.412.8−4.2814.08t1in−3.18.5−19.6113.75
out−5.95.4−12.059.74out−7.45.5−26.9217.46
RT4int0−11.220.83.1117.97t0inRB1−8.211.994.0411.48
out−1.58.394.9511.51out21.226.990.2217.59
int1−7.110.1−5.356.57t1in−10.011.081.8818.99
out−0.810.197.167.52out60.124.746.9910.76
RT5int01.311.469.5313.28t0inRB54.79.132.3810.02
out7.74.377.9519.08out4.914.739.9921.93
int1−6.410.436.7710.96t1in−17.117.738.677.59
out−15.014.942.005.00out7.910.137.859.30
Daphnia magna
Conventional20182019Organic20182019
(%)ds(%)ds(%)ds(%)ds
RT1int037.051.80.00.0t0inRB4 (2018)
RB6 (2019)
25.027.80.00.0
out97.07.10.00.0out7.510.40.00.0
int10.00.00.00.0t1in0.00.00.00.0
out7.510.40.00.0out2.57.10.00.0
RT4int062.551.85.00.0t0inRB10.00.05.010.0
out15.023.30.00.0out97.57.155.019.1
int10.00.00.00.0t1in2.57.145.019.1
out50.053.50.00.0out15.09.375.019.1
RT5int025.020.715.010.0t0inRB50.00.015.019.1
out5.09.345.034.2out0.00.040.023.1
int10.00.025.019.1t1in0.00.020.016.3
out0.00.025.010.0out2.57.10.00.0
Raphidocelis subcapitata
Conventional20182019Organic20182019
(%)ds(%)ds(%)ds(%)ds
RT1int092.02.196.81.7t0inRB4 (2018)
RB6 (2019)
100.011.095.01.2
out91.12.997.21.9out100.04.395.33.1
int187.00.698.31.7t1in89.00.699.42.8
out100.03.897.64.4out86.02.199.42.0
RT4int0100.07.689.14.5t0inRB1100.010.388.53.9
out93.02.397.83.5out84.55.170.37.3
int187.03.394.51.7t1in100.07.2100.04.6
out87.32.480.53.9out90.011.977.18.3
RT5int097.52.955.18.3t0inRB584.12.2100.03.7
out96.33.836.64.5out93.02.378.43.7
int195.51.799.05.4t1in94.31.051.45.9
out94.81.475.35.7out98.72.276.02.6
Spirodela polyrhiza
Conventional2018Organic2018
(%)sd (%)(%)sd (%)
RT1int0−37.933.8t0inRB4/RB6−16.347.2
out−0.623.8out−5.135.6
int1−10.434.4t1in−19.526.8
out−8.034.5out−11.837.2
RT4int0−15.628.0t0inRB1−44.651.0
out−2.630.5out−35.932.9
int112.732.7t1in−19.545.4
out−34.427.2out−30.242.6
RT5int025.620.4t0inRB581.96.8
out10.027.7out51.510.2
int1−0.226.2t1in−1.727.8
out23.723.0out−3.331.9
Table 7. Ecotoxicological hazard (HQeco) applied to paddy field soil samples. The index is referred to as a bioassay battery with 7 test species. Colors refer to the level of hazard as follows: green = absent; orange = moderate; red = major; black = severe.
Table 7. Ecotoxicological hazard (HQeco) applied to paddy field soil samples. The index is referred to as a bioassay battery with 7 test species. Colors refer to the level of hazard as follows: green = absent; orange = moderate; red = major; black = severe.
20182019
SoilWater SoilWater
IDSampling Point
and Campaign
HQecoSampling Point
and Campaign
HQecoIDSampling Point
and Campaign
HQecoSampling Point
and Campaign
HQeco
RB1emt00.16int02.37RB1emt00.10int05.13
cht00.14out8.37out7.32
int12.41t10.17int16.29
emt10.09out4.64out7.42
RB4emt00.59int03.6RB5emt00.10int03.63
out3.22out5.96
t10.13int12.37t10.05int13.94
out2.41out3.69
RB5emt00.07int02.42RB6emt00.38int05.06
out3.01out2.37
int12.37t10.01int12.37
t10.30out2.59out2.37
RT1cht00.01int03.14RT1emt00.35int02.37
out7.16out5.06
t10.68int12.37t10.45int12.37
out2.55out2.37
RT4emt00.29int05.48RT4emt00.66int05.06
out2.66in2.46
cht10.25int12.37t10.21int12.37
out4.84out5.06
RT5emt00.71int03.61RT5emt00.26int04.62
cht00.57out2.69out3.49
t10.50int12.37t1 int14.15
out2.370.98out4.95
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Onorati, F.; Tornambé, A.; Paina, A.; Maggi, C.; Sesta, G.; Berducci, M.T.; Bellucci, M.; Rivella, E.; D’Antoni, S. Ecotoxicological and Chemical Approach to Assessing Environmental Effects from Pesticide Use in Organic and Conventional Rice Paddies. Water 2022, 14, 4136. https://doi.org/10.3390/w14244136

AMA Style

Onorati F, Tornambé A, Paina A, Maggi C, Sesta G, Berducci MT, Bellucci M, Rivella E, D’Antoni S. Ecotoxicological and Chemical Approach to Assessing Environmental Effects from Pesticide Use in Organic and Conventional Rice Paddies. Water. 2022; 14(24):4136. https://doi.org/10.3390/w14244136

Chicago/Turabian Style

Onorati, Fulvio, Andrea Tornambé, Andrea Paina, Chiara Maggi, Giulio Sesta, Maria Teresa Berducci, Micol Bellucci, Enrico Rivella, and Susanna D’Antoni. 2022. "Ecotoxicological and Chemical Approach to Assessing Environmental Effects from Pesticide Use in Organic and Conventional Rice Paddies" Water 14, no. 24: 4136. https://doi.org/10.3390/w14244136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop