SciELO - Scientific Electronic Library Online

 
vol.38Manejo de plaguicidas y percepción de los agricultores sobre los problemas ambientales y de salud por el uso de plaguicidas en el estado de Yucatán, México: estudio de casoFluoruro en aguas subterráneas de la región centro-norte de México y su posible origen índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista internacional de contaminación ambiental

versión impresa ISSN 0188-4999

Rev. Int. Contam. Ambient vol.38  Ciudad de México  2022  Epub 25-Nov-2022

https://doi.org/10.20937/rica.54180 

ArticleS

Pesticide distribution and ecotoxicological risk in the Ayuquila-Armería river

Distribución de pesticidas y riesgo ecotoxicologico en el río Ayuquila-Armería

Brian Arturo Rodríguez-Aguilar1 

Luis Manuel Martínez-Rivera1  * 

Roberto Muñiz-Valencia2  3 

Norman Mercado-Silva4 

Luis Ignacio Íñiguez-Dávalos1 

Alejandro Aarón Peregrina-Lucano5 

1Departamento de Ecología y Recursos Naturales, Centro Universitario de la Costa Sur, Universidad de Guadalajara. Independencia Nacional 151, 48900 Autlán de Navarro, Jalisco, México.

2Facultad de Ciencias Químicas, Universidad de Colima, km 9 carretera Colima-Coquimatlán, 28400, Coquimatlán, Colima, México.

3Centro de Investigación en Recursos Naturales y Sustentabilidad, Universidad Bernardo O’Higgins, Fábrica 1990, 2º piso, Santiago, Chile.

4Centro de Investigación en Biodiversidad y Conservación, Universidad Autónoma del Estado de Morelos, Av. Universidad 1001, Col. Chamilpa, 62209 Cuernavaca, Morelos, México.

5Departamento de Farmacobiología, Centro Universitario de Ciencias Exactas e Ingeniería, Universidad de Guadalajara, Blvd. Marcelino García Barragán 1421, Calzada Olímpica, ,44430 Guadalajara, Jalisco, México.


ABSTRACT

Pesticides make an important contribution to increased global agricultural production; however, their intensive use represents a risk to biota in both the short and long terms. The physical and chemical characteristics of these compounds, as well as their versatility in agricultural, domestic, and public health use, have led them to become widely present in the environment. The objectives of this study were to evaluate the spatiotemporal distribution of pesticides in the surface water of the Ayuquila-Armería river and to conduct an ecotoxicological risk assessment in algae, macroinvertebrates, and fish. The presence of 11 analytes was determined, being λ-cyhalothrin, ametrine, and malathion the pesticides with highest concentrations, and carbendazim, malathion, and glyphosate those with the highest frequency of detection (> 70 %). The number of pesticides detected per sample presented no significant temporal differences, and there was no significant effect of the crops or vegetation adjacent to the study sites on the number of pesticides detected per sample. The sampling conducted in November 2019 (wet season) presented the greatest difference in concentrations of most pesticides, compared to the other samplings. The ecotoxicological risk assessment showed that the macroinvertebrates (Daphnia magna) are the most vulnerable to the concentrations of pesticides in surface water, followed by fish and algae. Regulation of the importation, trade, and management of malathion and λ-cyhalothrin is crucial in order to reduce the presence of their residues and the risks associated with them.

Key words: water pollution; agrochemical; intensive farming; pesticide monitoring; solid-phase extraction

RESUMEN

Los plaguicidas contribuyen de manera importante al aumento de la producción agrícola mundial; sin embargo, su uso intensivo representa un riesgo para la biota tanto a corto como a largo plazo. Las características físicas y químicas de estos compuestos, así como su versatilidad en el uso agrícola, doméstico y de salud pública, los han llevado a tener una amplia presencia en el ambiente. Los objetivos de este estudio fueron evaluar la distribución espaciotemporal de plaguicidas en agua superficial del río Ayuquila-Armería, y realizar una evaluación de riesgo ecotoxicológico en algas, macroinvertebrados y peces. Se determinó la presencia de 11 analitos, siendo la λ-cihalotrina, la ametrina y el malatión los plaguicidas con mayores concentraciones, en tanto que el carbendazim, el malatión y el glifosato fueron los de mayor frecuencia de detección (> 70 %). El número de plaguicidas detectados por muestra no presentó diferencias temporales significativas, y no hubo efecto significativo de los cultivos o la vegetación adyacente a los sitios de estudio sobre el número de plaguicidas detectados por muestra. El muestreo realizado en noviembre de 2019 (temporada de lluvias) presentó la mayor diferencia en las concentraciones de la mayoría de los plaguicidas, en comparación con los demás muestreos. La evaluación de riesgos ecotoxicológicos mostró que los macroinvertebrados (Daphnia magna) son los más vulnerables a las concentraciones de plaguicidas en el agua superficial, seguidos de los peces y las algas. La regulación de la importación, el comercio y la gestión del malatión y la λ-cihalotrina es fundamental para reducir la presencia de sus residuos y los riesgos asociados a ellos.

Palabras clave: contaminación del agua; agroquímicos; agricultura intensiva; monitoreo de plaguicidas; extracción en fase sólida

INTRODUCTION

The global production of synthetic chemical pesticides is approximately two million tons per year, for use in intensive agriculture for food production (Sharma et al. 2019). However, their intensive use and inefficient application impact the environment and represent a potential risk for living beings (Sumon et al. 2018). Mexico is a producer and exporter of foodstuffs at a global level; however, its production system is based on monocultures and excessive use of pesticides (SAGARPA 2017). The use of synthetic pesticides in Mexico began in the 1940s with the adoption of the technological paradigm for the modernization of agriculture known as the “green revolution” (Bejarano-González 2017).

Aquatic ecosystems are very vulnerable to contamination by pesticides, which can occur via specific or diffuse routes, such as surface runoff, erosion, leaching, and poor waste management (Sumon et al. 2018). Pesticides present high variation in surface waters due to environmental (precipitation, temperature, solar radiation) and anthropic (periods of application, productive activities) variables, which determine the presence and concentration of these chemicals (Belenguer et al. 2014, Vryzas 2018).

Ecotoxicological risk assessment is a useful tool for evaluating potential adverse effects on an ecosystem. Toxicity reference values depend on the response capacity of the species to one or a combination of different stressors, which reflects its sensitivity at a trophic level (Kuzmanović et al. 2016). The potential effect on the biota of stressors can be estimated using two models: Independent Action (IA) or Concentration Addition (CA) (Ccanccapa et al. 2016). Considering that surface water contamination due to the use of pesticides is characterized by the presence of a mixture of compounds, the CA model is considered the first approach to risk assessment (Ginebreda et al. 2014, Castro-Catalá et al. 2016). In this sense, the ecotoxicological risks associated with concentrations of pesticide can be evaluated using Risk Quotients (RQs) and Toxic Units (TUs) (Tsaboula et al. 2016, Carazo-Rojas et al. 2018). RQs are used to evaluate the ecotoxicological risk posed by the concentration of a single pesticide, in terms of chronic exposure, while TUs determine the toxic effects of a combination of pesticides, expressed as a single value per sampling site. These tools are useful for the determination of potential risks to different biotic components of an ecosystem. However, one disadvantage is the current lack of knowledge concerning the antagonistic and synergistic relationships that occur in a mixture of pesticides (Palma et al. 2014, Kuzmanović et al. 2016).

In Mexico, the Ayuquila-Armería basin is found between the states of Jalisco and Colima and has an area of 9864 km2. It has about 450 000 inhabitants and is considered one of the most important basins of western Mexico (Santana et al. 1993). The main economic activities practiced within the basin are intensive and subsistence agriculture, forestry, and fishing, with sugarcane monoculture one of the most economically important productive activities, and one that is dependent on the use of pesticides. The river is recognized as a provider of fishery resources and water at a regional level; however, these have been affected over time by activities that have modified the integrity of the aquatic ecosystem, such as agricultural, domestic, industrial, and livestock production discharges, as well as diversion and retention of water, overexploitation of fishery resources and deforestation of the riparian zones (Meza-Rodríguez 2006, Mercado-Silva et al. 2011).

The Ayuquila-Armería river is one of the 15 most important rivers among the 100 found on the Pacific slope of the country and presents important biodiversity in the two states it occupies, with endemic and threatened species according to Mexican Official Standard NOM-059-SEMARNAT-2010 (CONABIO 2020). In the basin, at least 143 active ingredients are available for sale, along with an even higher number of commercial products. Commonly used pesticides in agricultural activities within the basin include glyphosate, carbofuran, malathion, ametrine, paraquat, cypermethrin, and atrazine (Rodríguez-Aguilar et al. 2019). The aims of this study were (1) to evaluate the spatiotemporal distribution of the pesticides in the surface water of the Ayuquila-Armería river, and (2) to conduct an ecotoxicological risk assessment in algae, macroinvertebrates, and fish, based on pesticide residues in the surface water samples and adopting the RQ and TU methods.

METHODS

Sampling sites

Twenty-four sites were chosen based on the main sources of contamination by domestic, industrial, and agricultural discharges within the Ayuquila-Armería basin. These included the Melchor Ocampo sugar mill; the agave production area of Tuxcacuesco, Tonaya, Tolimán and San Gabriel Jalisco; the greenhouse zone of San Gabriel, Jalisco, and residual discharges from the cities of Autlán, El Grullo and Colima (Fig. 1). Three reference sites were established: one in the state of Colima (Zacualpan creek) and two in the state of Jalisco (La Taza spring and Manantlán creek). These reference sites were selected because they are close to or within the polygon of the Sierra de Manantlán Biosphere Reserve (SMBR) and present null or very low agricultural activity.

Fig. 1 Sampling sites and main sources of contamination in the Ayuquila-Armería basin in western Mexico. 

The sites are found in a surrounding manner, with different land uses associated with various productive activities, including irrigated and seasonal agriculture and the sugarcane and agave distillate production industries. The exceptions to this were the reference sites (8, 10, and 18), and those found mainly in zones with canyons (1, 9, 17, and 19), which generally present tropical low deciduous forest vegetation (Table I). The choice of pesticides was based on their sale, the extent, and frequency of application to crops. These pesticides comprise a mixture of active ingredients that control different types of organisms, such as fungicides, insecticides, and herbicides (Rodríguez-Aguilar et al. 2019).

TABLE I PHYSICAL MORPHOLOGY AND LAND USE AT THE SURFACE WATER SAMPLING SITES IN THE AYUQUILA-ARMERÍA RIVER BASIN. 

Study sites Land use Geomorphic type
Corcovado; Manantlán creek and Zenzontla in Jalisco; Alseseca Zacualpan creek and Zacualpan in Colima Deciduous tropical dry forest Canyon
El Chacalito; El Grullo bridge; Autlán and El Grullo sewage discharges; Palo Blanco and El Aguacate in Jalisco; San Miguel; San Buenaventura; El Chical; Arroyo Seco; Colima river; Asmoles sewage discharge and Armería in Colima Sugarcane Valley
La Taza spring in Jalisco Deciduous tropical dry forest Valley
Ayuquila river; Tuxcacuesco town; Tuxcacuesco river Agave and greenhouses Valley
La Croix in Jalisco Agave and greenhouses Canyon

Sample collection

Four water samples were taken at each site over the course of a year. Two were taken during the dry season (May 2018 and February 2019) and two during the rainy season (September and November 2018). Water samples were taken with 1 L jars, previously washed with methanol and ultrapure water prior to being used, in zones with a flowing current, from the middle of the water column and in the center of the river flow, based on the protocol published in NMX-AA-003-1980 (SECOFI 1980). During transport and fieldwork, the samples were placed at 4 ºC in order to increase the time of preservation and reduce microbial activity. Once in the laboratory, the samples were stored at -20 °C until analysis.

Materials and equipment

The pesticides under study were 2,4-D, acetochlor, ametrine, atrazine, carbendazim, carbofuran, diazinon, dimethoate, emamectin, glyphosate, imazalil, λ-cyhalothrin, malathion, methomyl, metoxuron, molinate, parathion, picloram, pyraclostrobin, and thiabendazole. Formic acid, acetonitrile, sodium chloride, ammonium formate, and methanol (MeOH), all analytical grade and sourced from Sigma-Aldrich, were used in the extraction of pesticides and preparation of the mobile phase of the chromatograph. Solid-phase extraction (SPE) cartridges Supel-select HLB 500 mg/12 mL (Supelco) were used.

A range of equipment was used during the process of extraction of pesticides from the water samples, including Milli-Q water purification system of Millipore, Vacuum pump EV-40 of EVAR, 24-port SPE Manifold system of Phenomenex, Orion 3 Star pH meter, and Cimarec stirring hotplate, both of Thermo Scientific. For quantification of the pesticides, a liquid chromatograph model 1200 was used, coupled to a mass spectrometry detector model 6430B, of Agilent Technologies. The software packages used for data acquisition were MassHunter Workstation Acquisition Software v. B.02.01 and MassHunter Workstation Quantitative Analysis Software v. B.03.02. A chromatographic column Zorbax Eclipse XDB-C18 Rapid Resolution 2.1 mm in diameter × 50 mm in length and particle size 3.5 µm of Agilent Technologies was used.

Extraction and analysis of pesticides

Water samples were filtered in a vacuum system with filter paper of pore size 0.45 µm. A volume of 400 mL of each sample was used for the extraction of the pesticides. Prior to extraction, 5 g of sodium chloride was added per 100 mL of sample. The SPE cartridges were then conditioned with 5 mL of MeOH and 10 mL of ultrapure water. Samples were passed through the cartridge at a flow rate of 2 mL/min. The cartridges were washed with 10 mL of ultrapure water followed by 5 min of vacuum drying. Elution of the analytes was conducted with 6 mL of MeOH. The samples were evaporated at 1 mL at ambient temperature, according to the procedure published by El-Osmani et al. (2014). The final extract was placed in 2 mL vials for injection into the liquid chromatography. The samples were analyzed based on Sierra-Díaz et al. (2019).

Risk quotient and toxic units for water

The toxicity values of effective concentration (EC50), lethal concentration (LC50), and no observable effect concentration (NOEC) are used in ecotoxicological risk assessments as points of comparison for the concentrations of pesticides in water and to establish their possible effect at a trophic level in an ecosystem (Carazo-Rojas et al. 2018). Ecotoxicological risk assessment was conducted using the approaches of RQs and TUs (Palma et al. 2014, Castro-Català et al. 2016) at three trophic levels of the aquatic ecosystem: (a) Raphidocelis subcapitata, Pseudokirchneriella subcapitata, and Scenedesmus subspicatus (algae), (b) Daphnia magna (aquatic invertebrates), and (c) Oncorhynchus mykiss, Lepomis macrochirus, and Pimephales promelas (fishes).

The toxicity reference values EC50 and CL50 were used to evaluate acute effects (Mu et al. 2016, Yuniari et al. 2016), while NOEC was used for the analysis of chronic effects (Warne and van Dam 2008). To evaluate the acute effects, acute toxicological data EC50 of 72 h were used for algae; EC50 of 48 h for the aquatic invertebrates, and LC50 of 96 h for the fishes. For the chronic effects, of each pesticide, data of 96 h NOEC were used for algae, 21 days NOEC for aquatic invertebrates, and 21 days NOEC for the fishes. This information was obtained from the Pesticide Properties Database developed by the Agriculture and Environment Research Unit (AERA) of the University of Hertfordshire (Lewis et al. 2016).

TUi = Ci/EC50 or LC50 (1)

where TUi is the toxic unit of the pesticide i; Ci is the quantified concentration of a pesticide (µg/L) in the water samples, and EC50 or LC50 (µg/L) are the effective and lethal concentrations that affect 50 % of the individuals when exposed to a given pesticide (Ccanccapa et al. 2016). The TUi was obtained for each of the pesticides in a sample and later summed to obtain the ecotoxicological risk level per site. Sites that present a value of 0.1-1 suggest moderate acute effects, and ≥ 1 suggest a high acute risk (Carazo-Rojas et al. 2018).

The toxic units per site (TU site ) were subsequently determined using a cumulative approach, conducted by summing the toxic units (TUi) of each pesticide quantified in a sample (Sánchez-Bayo et al. 2002). The TU site was obtained for each of the samples collected during the different samplings.

Sum TUsite = i=1nTUi (2)

The RQ values were calculated according to the following equation:

RQ = EC/PNEC (3)

where EC is the mean and maximum quantified concentration of each pesticide in the water samples, with the aim of evaluating the general and worst-case scenarios, respectively (Palma et al. 2014). The value of the predicted no-effect concentration (PNEC) was estimated through the NOEC. Where the NOEC value was unavailable in the database, EC50 was used. The PNEC values were then divided by 1000 as an assessment factor (AF) when at least a short-term assay was available (LC50 or EC50) at a trophic level. These were divided by 100 when data were available from long-term assays (EC10 or NOEC) with fish or aquatic invertebrates, and by 50 and 10 when two or three long-term assays (EC10 or NOEC) were available, respectively (EC 2011, Tsaboula et al. 2016). AF is an arbitrary value that considers the inherent uncertainty of the acquisition of toxicological laboratory data, according to the Technical Guidance for Deriving Environmental Quality Standards of the European Commission (EC 2011). Finally, if the RQ value was >1, it was considered that chronic adverse effects could be expected, if the value was 0.1-1, it was considered an intermediate risk and if it was < 0.1, the risk was considered low.

Statistical analysis

To evaluate the temporal distribution of the number of pesticides detected per sample, a statistical analysis was conducted using a generalized linear model (GLM) of Poisson type, in order to determine the significant differences between the number of pesticides detected per sample and the samplings conducted; likewise, to determine differences between the number of pesticides detected per sample and the land use adjacent to each study site.

To evaluate the spatiotemporal distribution of each pesticide, statistical analyses were performed in order to determine significant differences per season and per sampling. The statistical tests for each pesticide were conducted through left-censored data analysis, with the non-parametric Kaplan-Meier test and a parametric test with maximum likelihood estimation, using R project v. 3.6.2 and the packages Survival v. 3.1-8 and NADA v. 1.6-1. This type of analysis allows working with data found below the limit of quantification and conferring greater robustness to the results obtained, in place of truncating or substituting the results below the limit, whether for a zero or a fraction of the limit, causing bias in the results and increasing the probability of committing errors in the conclusions and subsequent decision-making (Helsel 2012). According to the theory of censored data analysis, it is not possible to conduct analysis with this statistical test for analytes that present censorship, or more than 80 % of results below the quantification limit, in the total number of samples (Hewett and Ganser 2007, Fox 2015). For this reason, the analytes acetochlor, dimethoate, emamectin, metoxuron, parathion, pyraclostrobin, thiabendazole, carbofuran, imazalil, and methomyl were excluded from the statistical analyses.

Finally, a temporal analysis was conducted of the TUs results produced by the ecotoxicological risk assessment for each level (algae, aquatic invertebrates, and fish). An analysis was conducted to determine differences among the results of the TUs of each sampling site per season, using the Wilcoxon Rank Sum test. A further analysis was conducted to determine differences among the results of the TUs of each sampling site per sample, using the Kruskal-Wallis test. The use of non-parametric tests was necessary since none of the results behaved normally according to the Shapiro-Wilk test.

RESULTS AND DISCUSSION

Occurrence and concentration of pesticides in surface water

The concentrations of 20 pesticides were analyzed in the samples (Fig. 2). Of these, six pesticides were invariably found below the limit of quantification (0.025 µg/L) (Table II). The maximum individual concentration detected was for λ-cyhalothrin at 50.92 µg/L in the site El Grullo bridge in the sample from September (wet season), followed by ametrine at 18.12 µg/L in the site El Grullo sewage discharge in the sample from May (dry season) and by malathion at 16.82 µg/L in the site La Croix in the sample from September.

Fig. 2 Pesticides detected by each site during the four samplings in the surface water of the Ayuquila-Armería river. The gray color in the box represents the pesticide positive detection. AME: ametrine, ATZ: atrazine, CBZ: carbendazim, CBF: carbofuran, DZN: diazinon, IMZ: imazalil, λCY: λ-cyhalothrin, MAL: malathion, MTM: methomyl, MLT: molinate, PIC: picloram, GPS: glyphosate, and 24D: 2,4-D. Thiabendazole was only detected at site 17. Acetochlor, dimethoate, emamectin, metoxuron, parathion, and pyraclostrobin were not detected. 

TABLE II MEAN AND MAXIMUM CONCENTRATION, AND STANDARD DEVIATION OF ALL DETECTED CONCENTRATIONS OF PESTICIDES IN SURFACE WATER SAMPLES FROM THE AYUQUILA-ARMERÍA RIVER. 

Pesticide* Mean concentration (µg/L) Maximum concentration (µg/L) SD Pesticide Mean (µg/L) Maximum concentration (µg/L) SD
2,4-D 0.14 0.81 0.19 Imazalil 0.19 0.93 0.25
Ametrine 0.58 18.12 2.36 l-cyhalothrin 1.75 50.92 9.13
Atrazine 0.47 2.59 0.46 Malathion 3.73 16.82 3.28
Carbendazim 0.26 1.86 0.37 Methomyl 0.09 0.12 0.01
Carbofuran 0.12 0.25 0.08 Molinate 0.18 0.40 0.09
Diazinon 0.73 10.08 2.06 Picloram 0.71 2.02 0.39
Glyphosate 0.73 3.34 0.64 Thiabendazole 0.15 0.24 **

SD: standard deviation.

*Acetochlor, dimethoate, emamectin, metoxuron, parathion and pyraclostrobin were below the limit of quantification in all samples; **pesticide quantification on only two occasions.

The pesticides with the greatest frequency of detection in the total number of samples were carbendazim, malathion, and glyphosate, being found in more than 70 % of the total number of samples analyzed. Their high frequencies of detection are associated with their wide sale and use in different crops in the study area, causing their residues to be distributed across most of the basin (Rodríguez-Aguilar et al. 2019).

The broad spectrum of use of malathion, probably with different periods of application in agriculture, has led to its permanence in the environment and thus its continuous detection (Deknock et al. 2019). However, the frequent detection of this pesticide isn’t only related to agricultural practices but also to its use in urban areas (Masiá et al. 2013); it isn’t only used to control insects in different crops, but also to control vectors of diseases and for domestic use. This leads to its permanent presence through continuous use and incorporation in the environment, which coincides with the definition of “pseudo-persistent” pollutants provided by Daughton (2003).

All the samples presented at least one pesticide. For this reason, the reference sites also presented concentrations of mixed pesticides in all their samples, apart from the sample taken from the Manantlán creek in September 2018 (wet season), in which only one analyte was detected. However, the Zacualpan creek presented a mean of five pesticides during the four samplings, the Manantlán creek a mean of 3.2 and La Taza spring presented a mean of 2.5, with the latter presenting the lowest mean relative to all the study sites.

According to various authors (Bailly-Comte et al. 2008, Guo et al. 2010), karst zones are vulnerable systems that are prone to contamination by pesticides, since their hydrogeological condition favors the movement of pollutants found on the surface towards the subterranean water. These zones present a minimum natural autodepuration of pollutants that, added to the rapid infiltration of precipitation, short residence times of the water, and diminution of photolysis and thermal decomposition of the pesticides, contribute to maintaining these pollutants in a stable condition and distributing them widely within these systems (Alam et al. 2014). The SMBR comprises Sierra Manantlán, Cerro de Enmedio, and Cerro Grande, the latter two of which are karstic mountainous massifs. The reference sites La Taza spring and Zacualpan creek are located on the slopes of Cerro de Enmedio and Cerro Grande, respectively, for which reason their hydrogeological situation favors the transfer of pesticides even though Cerro de Enmedio presents no agricultural activities, while in Cerro Grande the main activities are forestry and agriculture.

This indicates that some pesticides could already be contaminating subterranean waters used to supply the human population (Polanco-Rodríguez et al. 2018). The Zacualpan creek, one of the reference sites, supplies water to the conurbation of the municipalities of Colima and Villa de Álvarez and is considered the main source of potable water for both municipalities.

Some triazines and organophosphates have the capacity to move throughout the system of caverns and represent a risk of contamination for the aquifers (Lorenzo-Flores et al. 2017, Modrá et al. 2017). The results obtained not only revealed the presence of pesticides belonging to the groups of the triazines and organophosphates, but also pyrethroids, carbamates, thiocarbamates, benzimidazole, pyridine, and phosponoglycine, which not only increases the risk of the effects that these concentrations can have on the biota but also has the adverse effects on the other uses of the water from the springs, such as is human consumption.

Temporal and spatial distribution of pesticides in surface water

The analytes with the greatest frequency of detection during the dry season were malathion, glyphosate, and atrazine while, in the wet season, these were malathion, carbendazim, and glyphosate. Mean concentrations of ametrine, atrazine, carbofuran, malathion, picloram, and glyphosate were higher during the samplings conducted in the dry season compared to those of the wet season; however, carbendazim, diazinon, imazalil, λ-cyhalothrin, methomyl, molinate, 2,4-D, and thiabendazole presented the highest mean concentrations during the wet season. The highest concentration detected of ametrine during the dry season was 18.115 µg/L, while in the wet season this was 50.917 µg/L (Table III).

TABLE III RANGE, MEAN AND DETECTION FREQUENCY OF PESTICIDES QUANTIFIED DURING THE DRY AND WET SEASONS IN THE SURFACE WATER OF THE AYUQUILA-ARMERÍA RIVER. 

Analyte Dry season Wet season
Range (µg/L) Mean (µg/L) D.F. (%) Range (µg/L) Mean (µg/L) D.F. (%)
2,4-D LQ-0.771 0.12 31 LQ-0.806 0.16 40
Ametrine LQ-18.115 0.96 60 LQ-0.708 0.20 62
Atrazine LQ-2.587 0.60 67 LQ-0.904 0.31 55
Carbendazim LQ-0.464 0.11 60 LQ-1.857 0.37 81
Carbofuran LQ-0.248 0.15 6 LQ-0.120 0.08 6
Diazinon LQ-4.560 0.47 29 LQ-10.076 1.00 28
Glyphosate LQ-3.340 0.85 87 LQ-1.730 0.59 77
Imazalil LQ-0.236 0.16 10 LQ-0.928 0.21 15
λ-cyhalothrin LQ-0.336 0.10 31 LQ-50.917 3.30 34
Malathion LQ-12.140 4.38 94 LQ-16.819 2.97 81
Methomyl LQ-0.093 0.09 21 LQ-0.120 0.09 13
Molinate LQ-0.397 0.16 42 LQ-0.369 0.21 28
Picloram LQ-1.585 0.76 44 LQ-2.017 0.68 68
Thiabendazole LQ-0.06 * 2 LQ-0.241 * 2

DF: detection frequency; LQ: 0.025 µg/L.

*Only one quantification.

There are no temporal differences between the number of pesticides found in the study sites and the samplings (P = 0.371). These results differ from those presented by Masiá et al. (2013) since, in the case of the Ayuquila-Armería basin, the number of pesticides detected per sample over the entire period presented no temporal variation. The May 2018 sampling coincided with the period of least rainfall, which facilitated the process of pesticide concentration due to evaporation of the water and diminution of river flows. However, this same condition of low river flow facilitates greater contact between sunlight and the pesticides, causing even greater photolysis in the intertropical zones where the solar incidence angle is almost perpendicular, which could explain why the mean value for May was the lowest (Narváez-Valderrama et al. 2012).

The results of this study don’t coincide with that indicated by Belenguer et al. (2014) and Deknock et al. (2019), since there were no differences between the number of pesticides detected in the samples with respect to the crops and vegetation types adjacent to each sampling site (P = 0.271). Nevertheless, the results coincide with those of Carazo-Rojas et al. (2018) in that the frequency and concentration of the pesticides are not a function of the areas of intensive agriculture, since the reference sites of the Manantlán creek and the La Taza and Zacualpan creek, where agricultural activity is low or null, presented concentrations and presence of pesticides that were similar to other sites located in the agricultural valleys. The Zacualpan creek presented a mean of five different pesticides during the sampling period, which is equal or greater than the number detected in the sites San Miguel, Tuxcacuesco, El Chical, Arroyo Seco, and Armería, which are all located in agricultural zones where pesticides are intensively applied.

The wide spatial distribution of pesticides in the study area is due to a lack of training in their application related to the dosage and the pesticide mixing, the use of equipment that is obsolete or in poor condition, and the absence of regulation of the sale and adequate management of the containers of these agrochemicals (Rodríguez-Aguilar et al. 2019). In addition, the highly diffuse characteristics of certain pesticides also favor their dispersion across most of the surface water of the Ayuquila-Armería River (Rasmussen et al. 2016).

Of the 10 pesticides included in the statistical analysis with the censored data in order to evaluate their temporal distribution, only atrazine, carbendazim, malathion and glyphosate presented significant differences per season (Table IV). These pesticides present different intrinsic characteristics that range from high to low water solubility and from a moderate to low octanol-water partition coefficient (K ow ). For these reasons, their residuality in the surface water must respond to other intrinsic aspects of the pesticides, to the quantity or frequency of use, or to environmental factors within the study area.

TABLE IV RESULTS OF THE STATISTICAL TEST BY LEFT-CENSORED DATA ANALYSIS OF PESTICIDE CONCENTRATIONS PER SEASON AND SAMPLING IN THE SURFACE WATER OF THE AYUQUILA-ARMERÍA RIVER. 

Analyte Statistical test Censorship percentage (%) Season (p value) Sampling (p value)
2,4-D MLE 64.2 0.170 0.190
Ametrine K-M 38.9 0.600 0.100
Atrazine K-M 38.9 0.030* 0.060
Carbendazim K-M 29.5 5 × 10-5* 1 × 10-4*
Diazinon MLE 71.6 0.930 0.009*
Glyphosate K-M 17.9 0.040* 4e-10*
l-cyhalothrin MLE 67.4 0.42 3.8 × 10-8*
Malathion K-M 12.6 0.001* 2 × 10-14*
Molinate MLE 65.3 0.240 0.041*
Picloram K-M 44.2 0.090 9 × 10-5*

K-M: Kaplan-Meier, MLE: maximum likelihood estimation.

*Statistically significant.

Except for carbendazim, the pesticides malathion, atrazine, and glyphosate presented greater concentrations during the dry season, coinciding with the findings of Palma et al. (2014) and Papadakis et al. (2015) who indicated that periods of drought promote the concentration of pesticides due to reduced river flow. However, the dry season also coincides with the period of pesticide application in the sugarcane crop.

The results showed that the concentrations of carbendazim, diazinon, λ-cyhalothrin, malathion, molinate, picloram, and glyphosate in the surface water presented significant differences per sample (Table IV). The analytes that presented no differences per sample suggest that their wide temporal distribution could be the result of their frequent use in concentrations greater than those recommended, leading to their constant presence in the surface water, which coincides with that published by Carazo-Rojas et al. (2018). However, the results do not coincide with the findings of Belenguer et al. (2014), who indicate that the compounds that present the greatest persistence were those with relatively more stable values.

The analysis per sample provided more detail regarding the temporal behavior of the pesticides in the surface water of the Ayuquila-Armería river, compared to the analysis per season. The reduced concentrations of glyphosate, molinate, malathion, and diazinon during November could be a function of the increased frequency and intensity of hydrometeorological events due to the rainy season, facilitating the dilution of the analytes in the surface water, thus reducing the capacity for quantification (Carazo-Rojas et al. 2018).

Ecotoxicological risk assessment of pesticides in surface water

The results of the TU site values evidenced that pesticide concentrations can be associated with acute toxic effects for D. magna in all the study sites. However, this does not mean that the sites present no live invertebrates, since the effect generated by a pesticide in live organisms is a function of the concentration to which they are exposed, the toxicity of the pesticide, and the species involved, which results in different responses to the presence of pesticides. This situation is even more complex when pesticides are found in combination due to the antagonistic, synergistic, and cumulative relationships involved in their toxicity (Hoffman 2002, Rand 2008).

According to the results, the sites El Grullo bridge, San Miguel, La Croix, and Tuxcacuesco presented the highest TU site values for D. magna, standing out for its high values in the three trophic levels of El Grullo bridge. These results don’t coincide with those obtained by Ginebreda et al. (2014), Ccanccapa et al. (2016), and Carazo-Rojas et al. (2018), since they show that the risk associated with the combination of pesticides is distributed in neither a diffuse nor specific manner, nor does it respond to the loss of ecological quality. This study shows that the processes of pesticide dispersion have led to the presence of these chemicals across most of the basin, regardless of the degree of conservation of the sites or their surrounding productive activities (Table V).

TABLE V TOXIC UNITS (TU) FOR SITES IN ALGAE, AQUATIC INVERTEBRATES AND FISHES BASED ON THE CONCENTRATIONS OF ALL PESTICIDES DETECTED (EXCLUDING THOSE OF MALATHION AND λ-CYHALOTHRIN) IN EACH OF THE SAMPLING SITES OF THE AYUQUILA-ARMERÍA RIVER. 

Sites All pesticides Excluding malathion and λ-cyhalothrin
Algae Aquatic invertebrates Fishes Algae Aquatic invertebrates Fishes
Corcovado 0.05 7.72 0.64 0.03 0.13 < 0.01
El Chacalito 0.1 4.55 0.79 0.07 0.65 < 0.01
El Grullo bridge 10.35 232.05 242.88 0.17 0.1 < 0.01
Autlán sewage discharge 0.27 15.16 0.59 0.27 0.14 < 0.01
Palo Blanco 0.24 8.71 0.67 0.22 0.15 < 0.01
El Grullo sewage discharge 5.11 7.29 0.92 5.08 0.63 0.01
El Aguacate 0.35 9.73 1.14 0.31 0.26 < 0.01
Manantlán creek 0.11 16.97 1.11 0.08 4.56 < 0.01
Zenzontla 0.31 6.46 0.25 0.31 0.12 < 0.01
La Taza spring 0.02 8.99 0.77 < 0.01 0.01 < 0.01
Ayuquila river 0.35 6.61 0.85 0.32 0.01 < 0.01
San Miguel 0.02 17.58 0.3 0.02 10.08 0.01
San Buenaventura 0.08 10.14 0.39 0.08 0.1 < 0.01
La Croix 0.1 24.63 1.24 0.09 0.33 0.01
Tuxcacuesco town 0.14 5.11 0.89 0.11 0.15 < 0.01
Tuxcacuesco river 0.08 17.5 2 0.02 < 0.01 < 0.01
Alseseca 0.22 11.68 0.87 0.2 0.26 < 0.01
Zacualpan creek 0.06 15.95 1.26 0.03 0.01 < 0.01
Zacualpan 0.07 10.73 0.82 0.05 0.2 0.01
El Chical 0.1 6.81 0.6 0.09 0.1 < 0.01
Arroyo Seco 0.3 10.85 0.82 0.28 0.77 0.01
Colima river 0.07 14.22 1.2 0.05 0.05 < 0.01
Asmoles sewage discharge 0.28 15.55 1.27 0.25 0.11 < 0.01
Armería 0.03 10.27 0.9 0.01 0.25 < 0.01

Numbers in bold denote high risk of acute effects.

It is worth noting that most of the risk is caused by concentrations of malathion and λ-cyhalothrin. Excepting the concentrations of malathion and λ-cyhalothrin, only three sites could present a high potential risk (Table V). Of the three sites, Manantlán creek is the only one found in the SMBR and presents little agricultural activity; the high TU site value is due to the concentration of diazinon, from the group of the insecticides. This pesticide is highly toxic for aquatic invertebrates and could be used in different settings, such as agricultural, livestock production, and urban activities (EPA 2004). For this reason, its presence could be related to activities ongoing in the communities upstream of the sampling site, or to atmospheric transport and deposition with subsequent entry into the aquatic medium.

Malathion and λ-cyhalothrin have the characteristics of insecticides, which are authorized and recommended by the Mexican National Center of Preventive Programs and Disease Control (CENAPRECE, by its Spanish acronym) for the control of insect vectors of disease (CENAPRECE 2019). The dengue virus is a public health problem for the populations found within the Ayuquila-Armería basin, for which reason campaigns of control of the mosquito Aedes aegypti are common practice and are conducted through the spraying of insecticides in the urban zones (Anguiano-Moreno et al. 2011). The use of these chemicals in the control of vectors is a situation that contributes to the number of residues dispersed in the environment.

In most of the samples that presented high toxicity, this was due to the concentration of one or more compounds that contribute in greatest measure to the risk, and not to the set of pesticides in combination, which coincides with Stenström (2013). Regarding the TU site results for algae, El Grullo bridge and El Grullo sewage discharge were the sites that presented a potentially high risk. In El Grullo bridge, the concentration of λ-cyhalothrin presented the greatest effect on the growth of algae P. subcapitata, while the herbicide ametrine had the highest impact on the growth of the algae R. subcapitata in El Grullo sewage discharge. El Grullo bridge, El Aguacate, Manantlán creek, La Croix, Tuxcacuesco, Zacualpan creek, Colima river, and the Asmoles sewage discharge were the study sites that presented the greatest potential acute risk in fish (P. promelas and O. mykiss), along with the concentrations of malathion and λ-cyhalothrin, making the greatest contribution to this risk.

The result of the statistical analysis to evaluate the temporal variation of the TU site results per season through the Wilcoxon test showed that aquatic invertebrates are the only biotic component that presented differences (P = 0.005). The temporal analysis of the TU site values per sampling, through the Kruskal-Wallis test, showed that all the biotic components (algae, P = 0.015; aquatic invertebrates, P = 3.72-8; fishes, P = 0.002) presented differences. These differences in the three components are due to the diminution of the risk values obtained during the samplings conducted in the wet season (September and November 2018).

The analytes ametrine, carbendazim, and λ-cyhalothrin are associated with a chronic risk (RQ > 1) in algae (R. subcapitata, S. subspicatus, and P. subcapitata, respectively), at both their maximum and mean concentrations. Carbofuran, diazinon, λ-cyhalothrin, malathion, and methomyl also represent a danger for D. magna in both their maximum and mean concentrations. Finally, the concentrations of carbendazim, carbofuran, imazalil, l-cyhalothrin, malathion, and thiabendazole presented RQ values of moderate to high risk in both their maximum and mean concentrations (Table VI). The similar concentrations of a wide variety of pesticides represent a high risk, which leads to a high probability of potential affectation in the biotic communities of the aquatic ecosystem, mainly due to concentrations of ametrine, carbendazim, carbofuran, l-cyhalothrin, and malathion (Ccanccapa et al. 2016, Carazo-Rojas et al. 2018).

TABLE VI RESULTS OF THE RISK ASSESSMENT BY RISK QUOTIENT (RQ) IN ALGAE, AQUATIC INVERTEBRATES AND FISHES.* 

Pesticides Chronic 96 h NOEC and 72 h EC50 in algae Chronic 21 days NOEC and 48 h EC50 in aquatic invertebrates Chronic 21 days NOEC in fishes
PNEC (µg/L) RQ mean RQ max PNEC (µg/L) RQ mean RQ max PNEC (µg/L) RQ mean RQ max
Ametrine 0.072 8.0 252 6.4 <0.1 2.8 14 <0.1 1.3
Atrazine 10 <0.1 0.3 25 <0.1 0.1 200 <0.1 <0.1
Carbendazim 0.003 85.3 619 154 <0.1 <0.1 0.064 4.0 29.0
Carbofuran 320 <0.1 <0.1 0.8 0.1 0.3 0.22 0.5 1.1
Diazinon 1000 <0.1 1.0 0.056 13.0 179.9 70 <0.1 0.1
Imazalil 8.7 <0.1 0.1 35 <0.1 <0.1 0.43 0.4 2.2
λ-cyhalothrin 0.5 3.5 102 0.0002 7964 231441 0.003 565.2 16424.8
Malathion 260 <0.1 <0.1 0.0012 3112 14016 1.82 2.1 9.2
Methomyl 2000 <0.1 <0.1 0.032 2.8 3.8 1.52 <0.1 <0.1
Picloram 1204 <0.1 <0.1 135.8 <0.1 <0.1 11 <0.1 0.2
Thiabendazole 320 <0.1 <0.1 4.2 <0.1 <0.1 1.2 0.1 0.2

*The mean and maximum RQ for each pesticide represent the general and the worst-case scenarios, respectively.

NOEC: non-observable effect concentration, EC: effective concentration, PNEC: predicted no-effect concentration.

Numbers in bold denote high potential risk of chronic effects; 2,4-D, glyphosate and molinate presented RQ values lower than 0.1 in the three trophic levels.

Pesticide combination is common in the environment worldwide; however, the relationships that occur among them when in combination are generally unknown (Xie et al. 2019, De Souza et al. 2020). Organophosphate and carbamate concentrations can be managed using an additive approach since both are recognized as inhibitors of acetylcholinesterase; however, the effect of the mixture with other pesticide families remains largely unknown (Laetz et al. 2009).

Risk assessments of combinations of pesticides through an additive model of concentration would allow to determe the most problematic sites and evaluate the influence of concentrations of certain pesticides on the risk level, with the aim of directing actions for the regulation of the sale and use of these pesticides (Zhang et al. 2008). The results of the ecotoxicological risk assessment associated with concentrations of pesticides in the surface water of the Ayuquila-Armería basin should not be ignored, since the concentrations present are associated with both acute and chronic effects, mainly in the sites located within the valley of Autlán-El Grullo, in which there is a high potential of acute affectation at different trophic levels. Insecticides were found in greater measure with RQ values > 1, compared to the values obtained for herbicides and fungicides, which coincides with the findings of Palma et al. (2014).

CONCLUSIONS

Pesticides in the water samples from the Ayuquila-Armería river are current-use, their intensive application, in both quantity and frequency, leads to the presence of these compounds in combination, regardless of the season. However, concentrations of certain pesticides decreased during the November sampling in response to the general increase of the flow in the stretches of the river under study, with the consequent dilution effect.

Attention to the use and management of malathion, carbendazim, and glyphosate is important, given that they presented a frequency of detection of more than 70 % in the total samples. This suggests that their use is common and excessive in the basin, leading to the presence of their residues in most of the sites sampled.

This study used the model of concentration addition for the risk assessment per site. This procedure does not consider the antagonistic or synergistic relationships that can occur with a mixture of pesticides, but for practical terms, would allow determining the sites under greatest ecological risk. This is crucial in order to focus efforts on the regulation of certain pesticides and management measures for the protection of species associated with the aquatic ecosystem.

The results with the TU site for the determination of potential acute risk showed a high risk for aquatic invertebrates (D. magna) in all the study sites, being the most sensitive group. Followed by the group of fish, the sites that presented a high risk were El Grullo bridge (P. promelas), El Aguacate, Manantlán creek, La Croix, Tuxcacuesco, Zacualpan creek, Colima river and Asmoles sewage discharge (P. promelas and O. mykiss, both species in the previous sites). And finally, the least sensitive group was algae, being the sites El Grullo bridge (P. Subcapitata) and El Grullo sewage discharge (R. subcapitata) the ones that presented a high risk. The concentrations of malathion and l-cyhalothrin influenced to the greatest extent the TU site values for most of the study sites.

The use of RQ as a tool to determine chronic risk showed that all three trophic levels are potentially influenced by chronic effects, mainly due to the concentrations of ametrine, carbendazim, carbofuran, λ-cyhalothrin, and malathion. The use of TU site and RQ to determine acute and chronic risk, respectively, are viable tools to evaluate the risk in different trophic levels, with the aim of protecting the integrity of an ecosystem and generating suitable information for natural resources managers and for the development of public policies.

Regulation of the importation, sale, and management of synthetic pesticides has become a substantial aspect in reducing the concentration of their residues in surface waters of the country’s rivers. Special attention should be paid to malathion and λ-cyhalothrin, which generate the greatest influence on the associated risks in this study; however, it is possible that the reduction of these two pesticides alone will not resolve the situation of pesticides in the basin, because more active ingredients may be present, which generate greater impacts. For this reason, it is recommended to generate a long-term monitoring plan for the surface waters, considering a greater number of pesticides that are currently used within the basin and were not considered in this study.

REFERENCES

Alam M.J., Daoxian Y., Jiang Y.J., Yuchuan S., Yong L. and Xin X. (2014). Sources and transports of organochlorine pesticides in the Nanshan underground river, China. Environmental Earth Science 71, 1977-1987. https://doi.org/10.1007/s12665-013-2919-5 [ Links ]

Anguiano-Moreno M., Aguayo-López M.A., Álvarez-Lucas C.H., Torres-Ornelas P. and Lara-Esqueda A. (2011). Estrategia estatal de combate al dengue en Colima. Medicina Interna de México 27, 131-140 [ Links ]

Bailly-Comte V., Jourde H., Roesch A., Pistre S. and Batiot-Guilhe C. (2008). Time series analyses for karst/river interactions assessment: case of the Coulazou river (southern France). Journal of Hydrology 349, 98-114. https://doi.org/10.1016/j.jhydrol.2007.10.028 [ Links ]

Bejarano-González F. (2017). Los plaguicidas altamente peligrosos nuevo tema normativo internacional y su perfil nacional en México. Red de Acción sobre Plaguicidas y Alternativas en México, Estado de México, Mexico, 351 pp. [ Links ]

Belenguer V., Martinez-Capel F., Masiá A. and Picó Y. (2014). Patterns of presence and concentration of pesticides in fish and waters of the Júcar river (eastern Spain). Journal of Hazardous Materials 265, 271-279. https://doi.org/10.1016/j.jhazmat.2013.11.016 [ Links ]

Carazo-Rojas E., Pérez-Rojas G., Pérez-Villanueva M., Chinchilla-Soto C., Chin-Pampillo J., Aguilar-Mora P., Alpízar-Marín M., Masís-Mora M., Rodríguez-Rodríguez C.E. and Vryzas Z. (2018). Pesticide monitoring and ecotoxicological risk assessment in surface water bodies and sediments of a tropical agro-ecosystem. Environmental Pollution 241, 800-809. https://doi.org/10.1016/j.envpol.2018.06.020 [ Links ]

Castro-Català N., Kuzmanovic M., Roig N., Sierra J., Ginebreda A., Barceló D., Pérez S., Petrovic M., Picó Y., Chuhmacher M. and Muñoz I. (2016). Ecotoxicity of sediments in rivers: invertebrate community, toxicity bioassays and the toxic unit approach as complementary assessment tools. Science of the Total Environment 540, 297-306. https://doi.org/10.1016/j.scitotenv.2015.06.071 [ Links ]

Ccanccapa A., Masiá A., Navarro-Ortega A., Picó Y. and Barceló D. (2016). Pesticides in the Ebro river basin: occurrence and risk assessment. Environmental Pollution 211, 414-424. https://doi.org/10.1016/j.envpol.2015.12.059 [ Links ]

CENAPRECE (2019). Productos recomendados por el CENAPRECE para el combate de insectos vectores de enfermedades a partir de 2019. Centro Nacional de Programas Preventivos y Control de Enfermedades [online]. https://www.gob.mx/cms/uploads/attachment/file/469289/Lista_de_Insumos_Recomendados_por_el_CENAPRECE.pdf 13/11/20 [ Links ]

CONABIO (2020). 25. Ríos Purificación y Armería. Comisión Nacional para el Conocimiento y Uso de la Biodiversidad [online]. https://www.conabio.gob.mx/conocimiento/regionalizacion/doctos/rhp_025.html 18/06/20 [ Links ]

Daughton C.G. (2003). Cradle-to-cradle stewardship of drugs for minimizing their environmental disposition while promoting human health. I. Rationale for and avenues toward a green pharmacy. Environmental Health Perspectives 111, 757-774. https://doi.org/10.1289/ehp.5947 [ Links ]

De Souza R.M., Seibert D., Quesada H.B., Bassetti F.J., Fagundes-Klen R. and Bergamasco R. (2020). Occurrence, impacts and general aspects of pesticides in surface water: a review. Process Safety and Environmental Protection 135, 22-37. https://doi.org/10.1016/j.psep.2019.12.035 [ Links ]

Deknock A., Troyer N.D., Houbraken M., Dominguez-Granda L., Nolivos I., Echelpoel W.V., Forio M.A.E., Spanoghe P. and Goethals P. (2019). Distribution of agricultural pesticides in the freshwater environment of the Guayas river basin (Ecuador). Science of the Total Environment 646, 996-1008. https://doi.org/10.1016/j.scitotenv.2018.07.185 [ Links ]

El-Osmani R., Net S., Dumoulin D., Bigan M., Ouddane B. and Baroudi M. (2014). An experimental design approach to the optimization of pesticide extraction from water. Analytical Methods 6, 6514-6521. https://doi.org/10.1039/c4ay00610k [ Links ]

EPA (2004). Interim reregistration eligibility decision for diazinon EPA case No. (0238). Environmental Protection Agency [online]. https://archive.epa.gov/pesticides/reregistration/web/pdf/diazinon_red.pdf 18/06/20 [ Links ]

EC (2011). Guidance document No. 27. Technical guidance for deriving environmental quality standards. European Commission [online]. https://circabc.europa.eu/sd/a/0cc3581b-5f65-4b6f-91c6-433a1e947838/TGD-EQS%20CIS-WFD%2027%20EC%202011.pdf 18/06/20 [ Links ]

Fox G.A. (2015). What you don’t know can hurt you: Censored and truncated data in ecological research. In: Ecological statistics: Contemporary theory and application (G.A. Fox, S. Negrete-Yankelevich and V.J. Sosa, Eds.). Oxford University Press, New York, USA, pp. 106-130 [ Links ]

Ginebreda A., Kuzmanovic M., Guasch H., López de Alda M., López-Doval J., Muñoz I., Ricart M., Romaní A.M., Sabater S. and Barcerló D. (2014). Assessment of multi-chemical pollution in aquatic ecosystems using toxic unit: compound prioritization, mixture, characterization and relationships with biological descriptors. Science of the Total Environment 468, 715-723. https://doi.org/10.1016/j.scitotenv.2013.08.086 [ Links ]

Guo F., Yuan D. and Qin Z. (2010). Groundwater contamination in karst areas of southwestern China and recommended countermeasures. Acta Carsologica 39, 389-399. https://doi.org/10.3986/ac.v39i2.107 [ Links ]

Helsel D.R. (2012). Statistics for censored environmental data using Minitab and R. John Wiley and Sons, New Jersey, USA, 324 pp. [ Links ]

Hewett P. and Ganser G.H. (2007). A comparison of several methods for analyzing censored data. Annals of Occupational Hygiene 51, 611-632. https://doi.org/10.1093/annhyg/mem045 [ Links ]

Hoffman D.J. (2002). Wildlife toxicity testing. In: Handbook of ecotoxicology (D.J. Hoffman, B.A. Rattner, G.A. Burton and J. Cairns, Eds.). Lewis Publishers, Florida, USA, pp. 75-110. [ Links ]

Kuzmanović M., López-Doval J.C., De Castro-Català N., Guasch H., Petrović M., Muñoz I., Ginebreda A. and Barceló D. (2016). Ecotoxicological risk assessment of chemical pollution in four Iberian river basins and its relationship with the aquatic macroinvertebrate community status. Science of the Total Environment 540, 324-333. https://doi.org/10.1016/j.scitotenv.2015.06.112 [ Links ]

Laetz C.A., Baldwin D.H., Collier T.K., Hebert V., Stark J.D. and Scholz N.L. (2009). The synergistic toxicity of pesticide mixtures: implications for risk assessment and the conservation of endangered Pacific Salmon. Environmental Health Perspectives 117: 348-353. https://doi.org/10.1289/ehp.0800096 [ Links ]

Lorenzo-Flores A., Giácoman-Vallejos G., Ponce-Caballero M.C. and Ghoveisi H. (2017). Adsorption of organophosphorus pesticides in tropical soils: the case of karst landscape of northwestern Yucatan. Chemosphere 166, 292-299. https://doi.org/10.1016/j.chemosphere.2016.09.109 [ Links ]

Lewis K.A., Tzilivakis J., Warner D.J. and Green A. (2016). An international database for pesticide risk assessment and management. Human and Ecological Risk Assessment 22 (4), 1-15. https://doi.org/10.1080/10807039.2015.1133242 [ Links ]

Masiá A., Campo J., Vázquez-Roig P., Blasco C. and Picó Y. (2013). Screening of currently used pesticides in water, sediments and biota of the Guadalquivir river basin (Spain). Journal of Hazardous Materials 263, 96-104. https://doi.org/10.1016/j.jhazmat.2013.09.035 [ Links ]

Mercado-Silva N., Santana-Castellón E., Martínez-Rivera L.M., Lyons J. and Moermond T. (2011). Subsistence fisheries in the Sierra Manantlán Biosphere Reserve (Jalisco/Colima, Mexico). E-Gnosis 9, 1-19. [ Links ]

Meza-Rodríguez D. (2006). Caracterización hidrográfica de la cuenca Ayuquila-Armería mediante la aplicación de un sistema de información geográfica. B.Sc. Thesis. University of Guadalajara, Autlán de Navarro, Jalisco, Mexico, 134 pp. [ Links ]

Modrá H., Gruberová E., Konečný O., Ulmann V., Kaucká P., Vlková M., Tůma A., Halešová T., Kudělka J., Geršl M. and Pavlík I. (2017). Influx and concentration of triazine pesticides in the Amaterska cave system, Moravian karst, Czech Republic. Journal of Soils and Sediments 18, 640-647 https://doi.org/10.1007/s11368-017-1831-0 [ Links ]

Mu C., Wu X., Ma H., Tao W., Zhang G., Xia X., Shen J., Mai J., Sun T., Sun X., Arlinghaus R.B. and Shen H. (2016). Effective concentration of a multikinase inhibitor within bone marrow correlates with in vitro cell killing in therapy-resistant chronic myeloid leukemia. Molecular Cancer Therapeutics 15, 899-910. https://doi.org/10.1158/1535-7163.MCT-15-0577-T [ Links ]

Narváez-Valderrama J.F., Palacio-Baena J.A. and Molina-Pérez F.J. (2012). Persistencia de plaguicidas en el ambiente y su ecotoxicidad: una revisión de los procesos de degradación natural. Gestión y Ambiente 15, 27-38. [ Links ]

Palma P., Köck-Schulmeyer M., Alvarenga P., Ledo L., Barbosa I.R., López de Alda M. and Barceló D. (2014). Risk assessment of pesticides detected in surface water of the Alqueva reservoir (Guadiana basin, southern of Portugal). Science of the Total Environment 488, 208-219. https://doi.org/10.1016/j.scitotenv.2014.04.088 [ Links ]

Papadakis E.N., Vryzas Z., Kotopoulou A., Kintzikoglou K., Makris K.C. and Papadopoulou-Mourkidou E. (2015). A pesticide monitoring survey in rivers and lakes of northern Greece and its human and ecotoxicological risk assessment. Ecotoxicology and Environmental Safety 116, 1-9. https://doi.org/10.1016/j.ecoenv.2015.02.033 [ Links ]

Polanco-Rodríguez A.G., Riba-López M.I., Delvalss-Casillas A., Araujo-León J.A. and Banik S.D. (2018). Impact of pesticides in karst groundwater. Review of recent trends in Yucatan, Mexico. Groundwater for Sustainable Development 7, 20-29. https://doi.org/10.1016/j.gsd.2018.02.003 [ Links ]

Rand G.M. (2008). Fish toxicity studies. In: The toxicology of fishes (R.T. Di Giulio and D.E. Hinton, Eds.). CRC Press, Florida, USA, pp. 659-682. [ Links ]

Rasmussen J.J., Reiler E.M., Carazo E., Matarrita J., Muñoz A. and Cedergreen N. (2016). Influence of rice field agrochemicals on the ecological status of a tropical stream. Science of the Total Environment 542, 12-21. https://doi.org/10.1016/j.scitotenv.2015.10.062 [ Links ]

Rodríguez-Aguilar B.A., Martínez-Rivera L.M., Peregrina-Lucano A.A., Ortiz-Arrona C.I. and Cárdenas-Hernández O.G. (2019). Analysis of pesticide residues in the surface water of the Ayuquila-Armería river watershed, Mexico. Terra Latinoamericana 37, 151-161. https://doi.org/10.28940/terra.v37i2.462 [ Links ]

Sánchez-Bayo F., Baskaran S. and Kennedy I.R. (2002). Ecological relative risk (EcoRR): another approach for risk assessment of pesticides in agriculture. Agriculture, Ecosystems and Environment 91, 37-57. https://doi.org/10.1016/S0167-8809(01)00258-4 [ Links ]

Santana E., Navarro S., Martínez L.M., Aguirre A., Figueroa P. and Aguilar C. (1993). Contaminación, aprovechamiento y conservación de los recursos acuáticos del río Ayuquila, Reserva de la Biosfera Sierra de Manantlán, Jalisco-Colima. Tiempos de Ciencia 30, 29-38. [ Links ]

Sharma A., Kumar V., Shahzad B., Tanveer M., Sidhu G.P.S., Handa N., Kohli S.K., Yadav P., Bali A.S., Parihar R.D., Dar O.I., Singh K., Jasrotia S., Bakshi P., Ramakrishnan M., Kumar S., Bhardwaj R. and Thurkral A.K. (2019). Worldwide pesticide usage and its impacts on ecosystem. SN Applied Sciences 1, 1446. https://doi.org/10.1007/s42452-019-1485-1 [ Links ]

SAGARPA (2017). Planeación agrícola nacional 2017-2030. Secretaría de Agricultura, Ganadería Desarrollo Rural, Pesca y Alimentación [online]. https://www.gob.mx/cms/uploads/attachment/file/255627/Planeaci_n_Agr_cola_Nacional_2017-2030-_parte_uno.pdf 08/03/20 [ Links ]

SECOFI (1980). Norma Mexicana NMX-AA-003-1980. Aguas residuales. Muestreo. Secretaría de Comercio y Fomento Industrial, Mexico [online]. https://www.gob.mx/cms/uploads/attachment/file/166762/NMX-AA-003-1980.pdfLinks ]

Sierra-Díaz E., Celis-de la Rosa A.J., Lozano-Kasten F., Trasande L., Peregrina-Lucano A.A., Sandoval-Pinto E. and González-Chávez H. (2019). Urinary pesticide levels in children and adolescents residing in two agricultural communities in Mexico. International Journal of Environmental Research and Public Health 16, 562. https://doi.org/10.3390/ijerph16040562 [ Links ]

Stenström J.R. (2013). Mixture toxicity of pesticides and biological effects in agricultural streams. B.Sc. Thesis. Swedish University of Agricultural Sciences, Upsala, Swedish, 43 pp. [ Links ]

Sumon K.A., Rashid H., Peeters E.T., Bosma R.H. and Van den Brink P.J. (2018). Environmental monitoring and risk assessment of organophosphate pesticides in aquatic ecosystems of north-west Bangladesh. Chemosphere 206, 92-100. https://doi.org/10.1016/j.chemosphere.2018.04.167 [ Links ]

Tsaboula A., Papadakis E.N., Vryzas Z., Kotopoulou A., Kintzikoglou K. and Papadopoulou-Mourkidou E. (2016). Environmental and human risk hierarchy of pesticides: A prioritization method, based on monitoring hazard assessment and environmental fat. Environment International 91, 78-93. https://doi.org/10.1016/j.envint.2016.02.008 [ Links ]

Vryzas Z. (2018). Pesticide fate in soil-sediment-water environment in relation to contamination preventing actions. Current Opinion in Environmental Science and Health 4, 5-9. https://doi.org/10.1016/j.coesh.2018.03.001 [ Links ]

Warne M.J. and van Dam R. (2008). NOEC and LOEC data should no longer be generated or used. Australasian Journal of Ecotoxicology 14, 1-5. [ Links ]

Xie H., Wang X., Chen J., Li X., Jia G., Zou Y., Zhang Y. and Cui Y. (2019). Occurrence, distribution and ecological risk of antibiotics and pesticides in coastal waters around Liaodong peninsula, China. Science of the Total Environment 656, 946-951. https://doi.org/10.1016/j.scitotenv.2018.11.449 [ Links ]

Yuniari S.H., Hertika A.M.S. and Leksono A.S. (2016). Lethal concentration 50 (LC50 - 96 hours) Nile tilapia (Oreochromis niloticus) exposed cypermethrin-based pesticide. Journal of Experimental Life Science 6, 58-62. https://doi.org/10.21776/ub.jels.2016.006.02.01 [ Links ]

Zhang Y.H., Liu S.S. and Ge H.L. (2008). Prediction for the mixture toxicity of six organophosphorus pesticides to the luminescent bacterium Q67. Ecotoxicology and Environmental Safety 71, 880-888. https://doi.org/10.1016/j.ecoenv.2008.01.014 [ Links ]

Received: November 21, 2020; Accepted: November 25, 2021

* Author for correspondence: lmartinez@cucsur.udg.mx

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License