Despite the high concentrations of heavy metals (HMs) in the water bodies of northern Mexico, current monitoring efforts and regulations are insufficient. This is especially troubling, given the indigenous and vulnerable communities and the overall lack of public health information. To characterize exposure, we conducted an exhaustive review of HM studies in drinking water in northern Mexico. This information was used to estimate hazard quotient (HQ) and carcinogenic risk (CR) values by age for all of northern Mexico. In total, 316 concentrations of HMs were reported in 55 studies from 1975 to 2021, with concentration ranges of 0–11,100, 0–5,250, 0–350, 0–25, and 0–9,470 μg/L for arsenic, cadmium, chromium, mercury, and lead, respectively. The probable maximum concentrations (95th percentile) of metals in drinking water were 244.55, 10, 56, 1.68, and 90.35 μg/L for arsenic, cadmium, chromium, mercury, and lead, respectively. The HQ for all HMs ranged from 0.0011 to 404.62. Children younger than 2 years had an extremely high risk (HQ > 40) of presenting adverse health effects from arsenic consumption. Children younger than 6 years had the highest risk of developing cancer, with CR values from 2.24 × 10−5 to 7.77 × 10−1, emphasizing the need for continuous HM monitoring in northern Mexico.

  • Heavy metal concentrations in water in northern Mexico varied over time; almost half of the reported concentrations exceeded WHO standards.

  • Arsenic hazard quotients for children under 6 years were more than 100 times higher than the safe level.

  • Arsenic represents an extremely high non-carcinogenic hazard for northern Mexicans.

  • Carcinogenic risks related to As and Cr consumption varied from medium to extremely high.

Aquifers and rivers are essential drinking water sources for many communities in Mexico, especially in the northern and central areas of the country (Pacheco-Treviño & Manzano-Camarillo 2024; González-Horta et al. 2015). Unfortunately, population growth has induced hydrological changes and introduced contamination from minerals, heavy metals, and microorganisms, deteriorating the quality of these water sources and posing risks to human health (WHO 2017; Zhang et al. 2023).

Heavy metals are recognized as persistent pollutants that come from the Earth's crust and are distributed in the environment (Jomova et al. 2024; Singh & Kostova 2024). Unlike other metallic elements, heavy metals have a relatively high density (>4 g/cm3) (Jomova et al. 2024) and are toxic even in small concentrations (Zhang et al. 2023; Singh & Kostova 2024). Heavy metals can replace the functional group OH in enzymes, disrupting metabolic functions (Jomova et al. 2024) and often leading to enzyme inactivation, suppression of antioxidant defenses, and the generation of reactive oxygen species (ROS) in cells, which produce oxidative stress and lead to health conditions including congenital disorders, immune system problems, and tumor formation and progression (Balali-Mood et al. 2021; Zhang et al. 2023; Singh & Kostova 2024). Globally, the list of the most toxic heavy metals includes arsenic (As), cadmium (Cd), chromium (Cr), mercury (Hg), and lead (Pb) (Balali-Mood et al. 2021; Jomova et al. 2024). These metals contribute notably to ecological contamination (Singh & Kostova 2024) and are among the 13 metals included in the priority pollutant list of the United States (USEPA 2014); this classification is assigned to pollutants that have at least a 2.5% frequency of occurrence in water (USEPA 2024b). With the exception of Cr, these metals are also found in the global list of the ten chemicals of public health concern (WHO 2020).

Heavy metal concentrations in the environment have increased due to inputs from natural (e.g., soil mineral composition, hydrothermal activity, and volcanic activity) and anthropogenic (e.g., sewage and industrial waste) sources, although anthropogenic activities constitute the main source of these metals (Zhang et al. 2023). In northern Mexico, agriculture, aquaculture, mining, and other industrial activities are crucial sources of heavy metal environmental inputs and also influence their dispersion in the region (Jara-Marini et al. 2013). In addition to these sources, arsenic-rich minerals like arsenopyrite are responsible for natural As occurrence in groundwater in some parts of northern Mexico (Armienta & Segovia 2008); these minerals are widely distributed throughout Mexico and are susceptible to erosion by hot springs (Armienta & Segovia 2008; González-Horta et al. 2015).

In Mexico, chronic exposure to heavy metals in groundwater was first described in 1958 (McClintock et al. 2012). This exposure type has been associated with several adverse health effects, including nephropathy (Chávez-Gómez et al. 2017), a reduction of the intelligence quotient in infants (Thatcher et al. 1982), neurotoxicity, alterations in the cardiovascular and central nervous systems, cancer, and death (Londoño-Franco et al. 2016). Despite the health risks, people still use and consume water contaminated with heavy metals from household connections, public fountains, and wells (Sobsey et al. 2008). The quality of all water sources in northern Mexico cannot be evaluated due to the paucity of research. This has resulted in a lack of action or the implementation of ineffective measures, such as activated sludge systems (Guzmán-Colis et al. 2011) in certain sites where natural enrichment or heavy metal contamination has been detected. Therefore, the objective of this study was to compile reports of heavy metals in drinking water in northern Mexico from 1975 to 2021 and undertake a review of the associated carcinogenic and non-carcinogenic risks faced by the inhabitants of northern Mexico.

The continental territory of Mexico covers 1,959,248 km2 and borders the United States of America, Guatemala, and Belize (CONAGUA 2018). Mexico also has 11,122 km of coastline, with 15,000 km2 of coastal lagoons and 29,000 km2 of inland water bodies (SEMARNAT 2020). Mexico is ranked 94th in renewable water (renewed via the hydrological cycle) per capita among 152 countries, with 3,656 m3 of available water/inhabitant/year (CONAGUA 2018). However, according to the National Water Commission (Comisión Nacional del Agua [CONAGUA]), five hydrological-administrative regions located in northern Mexico exhibit high hydrological stress (40–100%), which is the amount of offstream water that is used out of the total amount available (CONAGUA 2018). In addition, the northwestern and central areas of Mexico cover two-thirds of the country and receive less than 500 mm of annual rainfall due to their arid and semi-arid climates; the northern, central, and northwestern regions host four-fifths of the population but only have access to one-third of the renewable water in the country (CONAGUA 2018).

At present, 115 of the 653 aquifers in the country are overexploited (SEMARNAT 2020). In addition, 168 of the 757 hydrological basins exhibit water deficits (negative water availability), with most of these basins located in northern Mexico (Salinas-Rodríguez & Martínez Pacheco 2024). This has notable implications for water management. To date, Mexico has established 147 decrees that ban groundwater use in certain areas, 3 aquifer regulations, 3 zone regulations, and 3 declarations of reserve zones for urban public use, which collectively protect approximately 55% of the national territory (CONAGUA 2018). Despite these measures, water management has not been effective, partly due to a lack of monitoring.

The National Water Quality Monitoring Network (Red Nacional de Medición de Calidad del Agua [RENAMECA]) has 5,028 monitoring stations; however, none of the RENAMECA stations monitor concentrations of hazardous elements such as heavy metals (CONAGUA 2018). The surface waters in the country are contaminated due to untreated municipal and industrial wastewater discharges, as well as agrochemical inputs (SEMARNAT 2020). Water contamination has led to a loss of ecosystem services, which have particularly affected rural communities and indigenous people (SEMARNAT 2020), given that rural populations, especially in arid northern areas, depend on groundwater (CONAGUA 2018).

Literature selection strategy

We performed a literature search in Google Scholar, Scopus, ScienceDirect, JSTOR, PubMed, ResearchGate, and Mexican governmental databases, using the following keywords: heavy metals, arsenic, cadmium, chromium, mercury, lead, water, and Mexico. In total, 950 studies were identified. In addition, 243 dissertations (Master's and doctoral theses) were identified in the digital repositories of the Autonomous University of Sinaloa, National Polytechnic Institute, University of Sonora, University of Nuevo León, and University of Arizona using the same keywords. Exclusion criteria were used to eliminate articles that contained information that was unrelated to the nine states of northern Mexico (Baja California, Baja California Sur, Chihuahua, Coahuila, Durango, Nuevo León, Sinaloa, Sonora, and Tamaulipas). However, studies from the central states of Aguascalientes, Zacatecas, and San Luis Potosí were not excluded, given that they share aquifers with Coahuila, Nuevo León, and Tamaulipas. The states selected for this review are shown in Figure 1.
Figure 1

States in northern Mexico (Baja California [1], Baja California Sur [2], Sonora [3], Chihuahua [4], Coahuila [5], Nuevo León [6], Tamaulipas [7], Sinaloa [8], Durango [9]) and central Mexico (Zacatecas [10], Aguascalientes [11], San Luis Potosí [12]) included in the bibliographic review.

Figure 1

States in northern Mexico (Baja California [1], Baja California Sur [2], Sonora [3], Chihuahua [4], Coahuila [5], Nuevo León [6], Tamaulipas [7], Sinaloa [8], Durango [9]) and central Mexico (Zacatecas [10], Aguascalientes [11], San Luis Potosí [12]) included in the bibliographic review.

Close modal
Studies were included if the authors stated that the water sources in the study were used for human consumption, if the studies compared heavy metal concentrations with the concentration limits listed in the Official Mexican Standard NOM-127-SSA1-1994 established by the General Directorate of Standards (Dirección General de Normas [DGN] 1994), or if the studies compared heavy metal concentrations with the concentration limits established by the World Health Organization (WHO 2017) (Table 1). On 2 May 2022, the Official Mexican Standard NOM-127-SSA1-2021 replaced NOM-127-SSA1-1994. However, this review included articles from 1975 to 2021; thus, NOM-127-SSA1-2021 was not considered in the inclusion criteria. The screening process and the retrieved studies are summarized in the flowchart (Figure 2), which was based on PRISMA guidelines (Page et al. 2021).
Table 1

Official concentration limits established by the Official Mexican Standard NOM-127-SSA1-1994 and by the World Health Organization (WHO) for heavy metals in drinking water

Heavy metalNOM-127-SSA1-1994 concentration limit (μg/L)NOM-127-SSA1-2021 concentration limit (μg/L)WHO concentration limit (μg/L)
Arsenic (As) 25 10* 10 
Cadmium (Cd) 3* 
Chromium (Cr) 50 50 50 
Mercury (Hg) 
Lead (Pb) 10 10 10 
Heavy metalNOM-127-SSA1-1994 concentration limit (μg/L)NOM-127-SSA1-2021 concentration limit (μg/L)WHO concentration limit (μg/L)
Arsenic (As) 25 10* 10 
Cadmium (Cd) 3* 
Chromium (Cr) 50 50 50 
Mercury (Hg) 
Lead (Pb) 10 10 10 

*The concentration limit does not yet apply to localities with less than 500,000 inhabitants (Diario Oficial de la Federación [DOF] 2022).

Figure 2

Flow diagram describing literature selection (diagram based on Page et al. (2021)).

Figure 2

Flow diagram describing literature selection (diagram based on Page et al. (2021)).

Close modal

After applying the inclusion and exclusion criteria, 55 studies that contained 316 concentrations of heavy metals were selected from 1975 to 2021. If a report presented multiple concentrations at the same metal for the same site, then the maximum and minimum values across these concentrations were considered. Similarly, if a report presented concentrations for multiple sites, the maximum and minimum values reported for each site were considered. If reports did not present multiple concentrations, then the value reported (usually a mean) was used. Some reports only presented a range of values. When results were represented solely by maps or graphs, concentration values ​​were extracted manually. Only a few reports included Supplementary data.

The following information related to the sampling area was also extracted from the studies: region (city or town), sample origin (e.g., well or river), sample collection year, and additional comments related to the selection of the report (Supplementary Table S1).

Statistical analysis of heavy metal concentrations

We estimated the quartiles of metal concentrations and represented their distributions with interquartile ranges (box plots). Different methods have been used to establish maximum concentrations of metals in water, soil, and air while minimizing the impact of extreme values (Frank et al. 2019; da Silva et al. 2020; USEPA 2021); among them, the 90th and 95th percentiles have been widely used (Spektor 1998; Frank et al. 2019; da Silva et al. 2020; USEPA 2021). Given that the 95th percentile has been used to establish the maximum concentration of metals in water (UKTAG 2020; Ravalli et al. 2022; Proshad et al. 2025), it was considered as the probable maximum concentration present in the water of northern Mexico (excluding extreme values). The 95th percentiles were calculated according to the National Institutes of Health guidelines (National Institutes of Health 2017). For this calculation, the difference between the 3rd quantile and the 1st quantile of heavy metal concentrations was determined (Supplementary Table S2). Then, the outliers to be eliminated were determined (Supplementary Table S3), and the 95th percentiles were calculated from the remaining data using the ‘PERCENTILE.INC’ function in Excel.

In risk characterizations, contaminant exposure levels are compared against chemical-specific toxicity criteria (established through research reports) to determine if the contamination levels constitute a health concern (USEPA 2001). The health risk associated with heavy metal exposure varies according to the metabolic state of individuals, the contamination source, environmental concentration, exposure time, chemical characteristics, and metal toxicity (Echeverry et al. 2015). In general, health risks can be classified as carcinogenic or non-carcinogenic and are calculated according to institutional guidelines such as those established by the United States Environmental Protection Agency (USEPA) (USEPA 2001, 2005). This agency recommends that reasonable exposure frequency estimates should be considered when there are no statistical data (USEPA 2001). Given that drinking water is a vital liquid associated with long-term exposure, we assumed a daily exposure frequency to calculate health risks. Some authors prefer to state that hazard assessments are ‘theoretical’ because mathematical formulas are used to calculate approximate measures (Hirano 1990) that do not represent the exact level of risks. However, for reasons of clarity in the text, the word ‘theoretical’ will not be included when mentioning risks.

We calculated the non-carcinogenic risk using the hazard quotient (HQ). To calculate the HQ, it is necessary to calculate the exposure dose (D) (mg/kg/day), which is recommended when exposure to contaminants occurs by ingesting drinking water (ATSDR 2022b); D was calculated with Equation (1):
(1)
where C is the environmental concentration of the metal in drinking water (mg/kg), IR is the daily ingestion rate (water consumption) from drinking water sources in northern Mexico (kg/day) by age (years) (Supplementary Table S4) (1 L of water was considered 1 kg), EF (unitless) is the exposure factor (1 [water is consumed daily, thus exposure is daily]) (ATSDR 2022b), and BW is the age-specific body weight (kg). To obtain weight information, we consulted the anthropometric database of the National Institute of Statistics and Geography (Instituto Nacional de Estadística y Geografía [INEGI]) (INEGI 2018); weights were classified by state and age (years) (Supplementary Table S5).
After calculating D, we calculated HQ (unitless) using Equation (2) (ATSDR 2022a):
(2)
where RfD is the reference dose established by USEPA (1995, 2000, 2012) or the tolerable daily intake (TDI) established by Health Canada (2010) (Table 2). If HQ > 1, then the heavy metal concentration exceeds the proposed safety limits, and consumers are at risk of developing adverse health effects (USEPA 2001, 2005; ATSDR 2022a).
Table 2

Proposed maximum reference doses for heavy metals in drinking water

MetalRfD, TDI (mg/kg/day)CSF, OSF (mg/kg/day)
Lead (Pb) 0.0036a 0.0085b 
Chromium (Cr) III 1.5c – 
Chromium (Cr) VI 0.0009d  0.27d 
Arsenic (As) 0.00006e 32e 
Cadmium (Cd) 0.0005f – 
Inorganic mercury (Hg) 0.0003f – 
MetalRfD, TDI (mg/kg/day)CSF, OSF (mg/kg/day)
Lead (Pb) 0.0036a 0.0085b 
Chromium (Cr) III 1.5c – 
Chromium (Cr) VI 0.0009d  0.27d 
Arsenic (As) 0.00006e 32e 
Cadmium (Cd) 0.0005f – 
Inorganic mercury (Hg) 0.0003f – 

Reference doses: RfD, reference dose; TDI, tolerable daily intake; CSF, cancer slope factor; OSF, oral slope factor.

Table 3

Classification of non-carcinogenic risks

Risk classificationValue
Slightly high >1–3 
Medium high >3–10 
High >10–25 
Very high >25–40 
Extremely high >40 
Risk classificationValue
Slightly high >1–3 
Medium high >3–10 
High >10–25 
Very high >25–40 
Extremely high >40 

Due to the lack of references that established a classification of HQ values, the HQ classification (Table 3) was arbitrarily chosen according to the HQ values calculated in this study.

The prevailing metric used to assess carcinogenic risk (CR) has been the lifetime average daily dose (LADD) (USEPA 2005); however, considering that heavy metal concentrations in northern Mexico have varied over time, a short-term exposure estimate was considered to be more appropriate than the LADD. We based this decision on USEPA recommendations (USEPA 2005), which state that when the dose rates are important in carcinogenic processes, a short-term exposure estimate may be more appropriate than the LADD when conducting risk assessments, given that using LADD in less-than-lifetime human exposure scenarios may underestimate or overestimate cancer risk. Furthermore, northern Mexicans were exposed to fluctuating heavy metal concentrations. Thus, it could not be assumed that older adults in northern Mexico had been exposed to the same concentrations for 78 years, which is the exposure period usually considered in the LADD (Ramirez-Andreotta & Superfund 2016). For this reason, we used Equation (3), which was derived from the cancer risk formula established by USEPA (2005) and ATSDR (ATSDR 2022a). Equation (3), which calculates the chronic daily intake (CDI), has been used to calculate carcinogenic risks from drinking water with heavy metals in states in central Mexico (Puebla and Colima) (Mendoza-Cano et al. 2017; Pérez-Castresana et al. 2019):
(3)
where C is the environmental concentration of the metal in drinking water (mg/kg), IR is the daily water consumption from drinking water sources in northern Mexico (kg/day) by age (years) (Supplementary Table S4), and BW is the age-specific body weight (kg) (Supplementary Table S5).
Then, to calculate CR, CDI was multiplied by the cancer slope factor (CSF) or the oral slope factor (OSF) established for each metal (Table 2) (USEPA 1995; OEHHA 2011a, b); CR < 1 × 10−6 is considered low (ATSDR 2022a). Cadmium and Hg do not have CSF values, as they are not considered to be conclusively carcinogenic (OEHHA 2006). The CR value was calculated with Equation (4) and classified by structure as described by Archundia et al. (2021) and Ramal et al. (2021) (Supplementary Table S6):
(4)

Heavy metals in the drinking water of northern Mexico

The review identified 55 studies with 316 heavy metal concentrations. Twelve states were considered for the review. Heavy metal contaminations in the water resources of northern Mexico varied by state. Pollution in Aguascalientes has been mainly due to aquifer overexploitation (González et al. 2012; CONAGUA 2019). In addition, the main water resource in the state, the San Pedro River, has been seriously polluted by wastewater discharges and the municipal slaughterhouse (Guzmán-Colis et al. 2011). The presence of heavy metals in Baja California Sur has been strongly influenced by mining activity, regional geothermalism, and the presence of chloro-fluoro green rocks and arsenopyrite (Cassassuce et al. 2005). In an exhaustive water quality study of 500 wells in Baja California Sur (2004–2005), it was reported that the wells used for human consumption contained 10–500 μg As/L, which was likely due to mineral residues and mine tailings; some residents, such as those in the community of San Antonio, drank water purified from municipal treatment plants, but in other ranches and communities, such as El Rosario, Valle Perdido, and Texcalama, people drank water directly from contaminated sources (Cassassuce et al. 2005).

Chihuahua has been recognized for having high As concentrations in its water bodies (Westerhoff et al. 2004), which are mainly due to volcanic processes, arsenopyrite deposits in hot spring areas, and anthropogenic contributions related to mining and metal smelting (Camacho et al. 2011). This type of contamination has been reported since 1995, when the central water well of La Casita contained As concentrations between 7,000 and 11,000 μg/L, which caused chronic arsenicosis in 8% of the inhabitants of the town (Ochoa-Reyes et al. 2009). Years later, in 2009, the water from the well in La Casita contained 2 μg As/L; however, prolonged exposure to high As concentrations has caused the inhabitants of the town to develop diseases (Ochoa-Reyes et al. 2009).

Numerous wells in the Guadiana Valley (Durango) were also studied, with As concentrations ranging from 5 to 167 μg/L and suspected contamination of geological origin (Alarcón-Herrera et al. 2001). Recently, Alcázar-Medina et al. (2020) reported that the community wells of San José de Avino (Panuco de Coronado municipality) contained 5,250 μg Cd/L, 2139.2 μg Pb/L, and 60.5 μg Cr/L; in general, contamination was considered to be of geological origin (Gamboa-Loira et al. 2020). A similar situation was observed in Nuevo León, where the most important water resource of this state, the Río San Juan basin, has exhibited variable heavy metal concentrations over the last two decades (Flores-Laureano 1997). Indeed, it is possible that ∼37% of the inhabitants of Nuevo León are drinking water with As concentrations greater than 10 μg/L (maximum of 34.54 μg As/L) (Gamboa-Loira et al. 2020). Similarly, children from communities in San Luis Potosí may have suffered adverse effects due to high As levels in drinking water (Díaz-Barriga et al. 1993). López-Álvarez et al. (2013) suggested the presence of cemeteries, mine tailings, wastewater discharges, solid waste disposal, and gas stations as possible pollution sources.

Contamination has also affected the indigenous communities of Sonora. Several towns in the Yaqui Valley (southern Sonora) were exposed to 4.71–49.26 μg As/L from 2001 to 2002 (Meza et al. 2004, 2005). The most recent study in the Yaqui and Mayo Valleys reported 57 wells with As concentrations ranging from 2.7 to 98.7 μg As/L (García-Rico et al. 2019). Similar contamination levels have been observed in Coahuila; wells in the southwest had Pb and CD concentrations of 8–51 and 1–18 μg/L, respectively (Saldarriaga et al. 2014).

Zacatecas is abundant in metallic minerals. Mining activity in the state is high, which has resulted in the pollution of many aquatic resources in the San Ramón system (González-Dávila 2013). The municipal water of Ojocaliente City contained 40 μg As/L (Sandoval-Alvarado et al. 2020), exceeding the Official Mexican Standard (DGN 1994). Fortunately, the Cr concentrations in the water sources of the Mazapil municipality did not exceed the Official Mexican Standard (<10 μg/L) (Escamilla-Rodríguez et al. 2021).

There are few studies on the quality of freshwater in Sinaloa despite the vast amount of land designated for agriculture (Frías-Sarmiento 2007) and mining (SGM 2019). Ayala-Rodríguez (2010) established that the water from the Sinaloa River near the towns of Las Quemazones and Nío had 4–130 μg Pb/L and 10–40 μg Cd/L, respectively. However, a few years later, all rivers in Sinaloa exhibited low Cd and Pb concentrations (Frías-Espericueta et al. 2014). Finally, the water wells in the Culiacán Valley exhibited heavy metal concentrations of 0.01–14.3 μg As/L, 0–0.41 μg Cd/L, 0–5.5 μg Cr/L, and 0.01–6.84 μg Pb/L (Rivera-Hernández et al. 2021). Lastly, only one report was available on drinking water contamination in Tamaulipas. López et al. (2020) reported that surface water from the state, sampled near Victoria City, presented 30–1,000 μg Pb/L and 0–140 μg Cd/L, respectively, with averages of 400 μg Pb/L and 140 μg Cd/L. This surface water provides drinking water to the population of Tamaulipas. A visual summary of the cumulative maximum concentrations reported in drinking water (surface, river, well, dam, underground, and bottled water) from 1990 to 2019 is shown in Figure 3.
Figure 3

Cumulative maximum concentrations of heavy metals reported in drinking water (surface, river, well, dam, underground, and bottled water) in northern Mexico from 1975 to 2021. The maps are arranged according to the year (horizontal) and metal (vertical). The concentration scale appears at the top for each metal. Individual maps were created with Microsoft Excel.

Figure 3

Cumulative maximum concentrations of heavy metals reported in drinking water (surface, river, well, dam, underground, and bottled water) in northern Mexico from 1975 to 2021. The maps are arranged according to the year (horizontal) and metal (vertical). The concentration scale appears at the top for each metal. Individual maps were created with Microsoft Excel.

Close modal

Descriptive statistics of heavy metal studies

Overall, 316 concentrations of heavy metals were reported in 55 studies from 1975 to 2021. Of these, 174 concentrations were of As, with values ranging from 0 to 11,100 μg/L; 55 concentrations were of Cd, with values ranging from 0 to 5,250 μg/L; 24 concentrations were of Cr, with values ranging from 0 to 350 μg/L; 8 concentrations were of Hg, with values ranging from 0 to 25 μg/L; and 55 concentrations were of Pb, with values ranging from 0 to 9,470 μg/L. The averages and standard deviations of the metal concentrations are shown in Table 4.

Table 4

Summary statistics of the 316 heavy metal concentrations (μg/L) in 55 studies from 1975 to 2021

MetalMinimumMaximumAverageStandard deviationN
Arsenic (As) 11,100 264.84 1287.21 174 
Cadmium (Cd) 5,250 159.64 755.53 55 
Chromium (Cr) 350 42.00 84.77 24 
Mercury (Hg) 25 3.52 8.71 
Lead (Pb) 9,470 282.19 1305.45 55 
MetalMinimumMaximumAverageStandard deviationN
Arsenic (As) 11,100 264.84 1287.21 174 
Cadmium (Cd) 5,250 159.64 755.53 55 
Chromium (Cr) 350 42.00 84.77 24 
Mercury (Hg) 25 3.52 8.71 
Lead (Pb) 9,470 282.19 1305.45 55 

All reported concentrations were highly variable over time and study area. On average, 43.04% of heavy metal concentrations exceeded the concentration limits established by the Official Mexican Standard NOM-127-SSA1-1994, and 50.63% exceeded the concentration limits established by the WHO (Table 5). Many of the As concentrations reported in the studies exceeded the Official Mexican Standard, with some exceeding the concentration limit by 200 times.

Table 5

Percentage of heavy metal values that exceeded the limits established for drinking water

Arsenic (As)Cadmium (Cd)Chromium (Cr)Lead (Pb)Mercury (Hg)
NOM-127-SSA1-1994 52.3% 25.5% 25% 41.82% 25% 
WHO standard 65.52% 27.3% 25% 41.82% 25% 
Arsenic (As)Cadmium (Cd)Chromium (Cr)Lead (Pb)Mercury (Hg)
NOM-127-SSA1-1994 52.3% 25.5% 25% 41.82% 25% 
WHO standard 65.52% 27.3% 25% 41.82% 25% 

The distribution of heavy metal concentrations in the drinking water of northern Mexico (Supplementary Figure S1) indicates that Coahuila and Nuevo León exhibited the highest As contamination over the years; this metalloid comes from natural compounds, such as arsenite and arsenate, and its concentration in water has increased mainly due to hydrothermalism, arsenical pesticides, and arsenical chemical waste (Jaishankar et al. 2014). Durango and Tamaulipas contained the highest Pb concentrations; this contamination may have come from fertilizers and pesticides, ore smelting, factory chimneys, soil wastes, and wastes from the battery industry (Jaishankar et al. 2014).

The probable maximum concentrations (95th percentile) for As, Cd, Cr, Hg, and Pb data were 244.55, 10.00, 56.00, 1.68, and 90.35 μg/L, respectively. All probable maximum concentrations exceeded the concentration limits established by the Official Mexican Standard NOM-127-SSA1-1994 and the concentration limits established by the WHO. The box plot with the original concentrations is shown in Supplementary Figure S2. The box plot with the transformed concentrations (Log10 scale) is shown in Figure 4.
Figure 4

Box plot of the heavy metal concentrations in northern Mexico from 1975 to 2021. The box represents 50% of the data. The first line of the box represents the first quartile. The median is indicated by the line in center line of the box. The third line represents the third quartile. The whiskers indicate the maximum and minimum values. Data points that fall outside the whiskers are outliers (exceeded the interquartile range from the first and third quartile by more than 1.5 times). The concentrations were transformed using the Log10 function to ensure the whiskers were visible. Values of 0 were omitted. Created with XLSTAT.

Figure 4

Box plot of the heavy metal concentrations in northern Mexico from 1975 to 2021. The box represents 50% of the data. The first line of the box represents the first quartile. The median is indicated by the line in center line of the box. The third line represents the third quartile. The whiskers indicate the maximum and minimum values. Data points that fall outside the whiskers are outliers (exceeded the interquartile range from the first and third quartile by more than 1.5 times). The concentrations were transformed using the Log10 function to ensure the whiskers were visible. Values of 0 were omitted. Created with XLSTAT.

Close modal

Heavy metal exposure risk

The maximum heavy metal concentrations were used to calculate the highest non-carcinogenic HQ that was most likely to occur in northern Mexico (Table 6). From the 54 HQ values, 28 ​​exceeded the value of 1 (Supplementary Figure S3). The highest HQ values were found in children younger than 6 years. The population that may be at greatest risk was comprised of children between 1 and 2 years in age; these values exceeded those of children under 1 year, given that children aged 1–2 years tended to drink more water.

Table 6

Non-carcinogenic hazard quotient (HQ) most likely to occur in northern Mexico by age

HQ value
AgeAsCdCr IIICr VIHgPb
<1 374.03 1.84 0.0034 5.71 0.51 2.30 
1–2 404.62 1.99 0.0037 6.18 0.56 2.49 
3–5 345.78 1.70 0.0032 5.28 0.48 2.13 
6–7 265.94 1.30 0.0024 4.06 0.37 1.64 
8–11 206.67 1.01 0.0019 3.16 0.28 1.27 
12–17 146.53 0.72 0.0013 2.24 0.20 0.90 
18–35 127.83 0.63 0.0012 1.95 0.18 0.79 
36–59 119.06 0.58 0.0011 1.82 0.16 0.73 
≥60 131.61 0.65 0.0012 2.01 0.18 0.81 
HQ value
AgeAsCdCr IIICr VIHgPb
<1 374.03 1.84 0.0034 5.71 0.51 2.30 
1–2 404.62 1.99 0.0037 6.18 0.56 2.49 
3–5 345.78 1.70 0.0032 5.28 0.48 2.13 
6–7 265.94 1.30 0.0024 4.06 0.37 1.64 
8–11 206.67 1.01 0.0019 3.16 0.28 1.27 
12–17 146.53 0.72 0.0013 2.24 0.20 0.90 
18–35 127.83 0.63 0.0012 1.95 0.18 0.79 
36–59 119.06 0.58 0.0011 1.82 0.16 0.73 
≥60 131.61 0.65 0.0012 2.01 0.18 0.81 

As, arsenic; Cd, cadmium; Cr, chromium; Hg, mercury; Pb, lead.

If the Cr ion that is consumed is Cr+3, then the consumption of this heavy metal would not represent a health risk; however, studies specifying Cr valence are scarce. Based on the HQ categories, the HQ values for Cd and Hg in drinking water indicate that there is either no risk or a slightly high health risk. The HQ values for Pb indicated no risk to a slightly high health risk. The HQ values for Cr+6 indicated a slightly high to medium-high health risk. Arsenic was responsible for the highest health risks, with extremely high hazard values (>40).

The maximum heavy metal concentrations were also used to calculate the most likely CR values in people in northern Mexico by age. According to the CR classification (Supplementary Table S6), the carcinogenic risks related to As for all ages were extremely high (>1 × 10−3) (Table 7). The CR values for Cr+6 varied from medium-high to extremely high risk. It should be noted that most international standards do not establish a limit for the concentration of Cr+6 in water, although they do establish limits for total Cr; Cr+6 is the only form of Cr reported to have carcinogenic effects (USEPA 2000). The CR values for Pb reflected low-medium risk and medium risk.

Table 7

Carcinogenic risk (CR) most likely to occur in people from northern Mexico by age

CR value
AgeArsenic (As)Chromium (Cr6)Lead (Pb)
<1 7.18 × 10−1 1.39 × 10−3 7.05 × 10−5 
1–2 7.77 × 10−1 1.50 × 10−3 7.62 × 10−5 
3–5 6.64 × 10−1 1.28 × 10−3 6.52 × 10−5 
6–7 5.11 × 10−1 9.87 × 10−4 5.01 × 10−5 
8–11 3.97 × 10−1 7.67 × 10−4 3.89 × 10−5 
12–17 2.81 × 10−1 5.44 × 10−4 2.76 × 10−5 
18–35 2.45 × 10−1 4.74 × 10−4 2.41 × 10−5 
36–59 2.29 × 10−1 4.42 × 10−4 2.24 × 10−5 
≥60 2.53 × 10−1 4.88 × 10−4 2.48 × 10−5 
CR value
AgeArsenic (As)Chromium (Cr6)Lead (Pb)
<1 7.18 × 10−1 1.39 × 10−3 7.05 × 10−5 
1–2 7.77 × 10−1 1.50 × 10−3 7.62 × 10−5 
3–5 6.64 × 10−1 1.28 × 10−3 6.52 × 10−5 
6–7 5.11 × 10−1 9.87 × 10−4 5.01 × 10−5 
8–11 3.97 × 10−1 7.67 × 10−4 3.89 × 10−5 
12–17 2.81 × 10−1 5.44 × 10−4 2.76 × 10−5 
18–35 2.45 × 10−1 4.74 × 10−4 2.41 × 10−5 
36–59 2.29 × 10−1 4.42 × 10−4 2.24 × 10−5 
≥60 2.53 × 10−1 4.88 × 10−4 2.48 × 10−5 

Descriptive statistics of heavy metal studies

The variation in heavy metal concentrations in drinking water in northern Mexico was mainly due to geological origin, with metals originating naturally from soil minerals (e.g., silicate, carbonate, oxide, hydroxide, and sulfide structures) and being present in dissolved, particulate, or complex (bonded to other metals) forms in aquatic ecosystems due to natural processes (e.g., weathering, volcanic eruptions, erosion, and runoff) or anthropogenic waste (Zhang et al. 2023). The distributions may also be influenced by the presence of arsenopyrite in northern Mexico and the extensive generation sites of arsenic trioxide (arsenolite), a by-product of gold and silver extraction, especially in the mining regions of Baja California Sur (Carrillo-Chávez et al. 2000). The anthropogenic sectors that may have contributed the most to heavy metal concentrations in drinking water are likely those that have received the greatest amounts of water from the Mexican government: agriculture (76.04%), electric (non-hydroelectric) power (4.72%), industry (4.86%), and public supply (14.38%) (CONAGUA 2018). The variability of heavy metal concentrations was also due to the mobility and bioavailability of metals being affected by chemical processes, such as sorption, desorption, and precipitation–dissolution reactions (Caporale & Violante 2016).

In northern Mexico, 35.6% of the concentrations of As in drinking water exceeded 50 μg/L. This is similar to what was reported in Bangladesh, where 50% of the tube wells contained As concentrations greater than 50 μg/L (Khan et al. 2003). While the outlook is better for northern Mexico than for Bangladesh, the situation in northern Mexico is still worrying. If members of a community drink 1 L of water with 50 μg/L of a heavy metal each day, then 13 out of every 1,000 people could develop a malignant neoplasm (Aryan et al. 2024), such as skin, liver, or other types of cancer, due to oxidative stress, chromosomal anomalies, and altered growth factor expression (Balali-Mood et al. 2021).

It is important to highlight that the range of As concentrations found in drinking water sources was greater than the concentration range found in the groundwater of South Asia: India (0–3,700 μg As/L, 2008–2022), Pakistan (0–4,730 μg As/L, 2009–2022), Bangladesh (0.5–4,730 μg As/L, 2010–2022), Nepal (0–2,620 μg As/L, 2006–2021), and Sri Lanka (0–400 μg As/L, 2013–2019) (Aryan et al. 2024). This situation is extremely worrying, as more than 100,000 deaths have been attributed to groundwater contamination in India, and in Bangladesh, ∼43,000 people die every year for the same reason (Aryan et al. 2024). Thus, a relatively high number of northern Mexicans may be dying due to high concentrations of As in drinking water. Furthermore, in Antofagasta, Chile, from 1958 to 1965, people were exposed to drinking water with 860 μg As/L, which was responsible for 18–24% of all infant deaths during that period (Ferreccio & Sancha 2006). Considering that children in northern Mexico have been exposed to similar and even higher concentrations in drinking water, As-contaminated drinking water may be responsible for a comparable infant mortality rate in Mexico.

Lead and Cd concentrations in water sources in northern Mexico were higher than those in the streams, wells, and stormwater in Anambra, Nigeria (Pb: 0–370 μg/L; Cd: 100–350 μg/L) (Emmanuel et al. 2022). However, Hg concentrations were lower than those reported in Anambra (50–380 μg/L) (Emmanuel et al. 2022). The mean concentrations of As, Cd, and Pb in drinking water from northern Mexico (264.84, 159.64, and 282.19 μg/L) were also higher than the average concentrations found in 71 dug-well water samples from houses from Polonnaruwa, Sri Lanka (As: 0.37 μg/L; Cd: 0.03 μg/L; Pb: 0.09 μg/L), in which people with chronic kidney disease (CKD) of unknown etiology (CKDu) were living. A disease similar to CKDu, CKD, is currently a public health problem in Mexico; it is one of the top 10 causes of mortality among patients of the Mexican Social Security Institute (∼377 cases per million inhabitants) (Chávez-Gómez et al. 2017). Although the presence of Pb, Cd, As, and Hg in water could not be associated with CKDu incidence (Kulathunga et al. 2022), CKD has been shown to be related to exposure to heavy metals (Chávez-Gómez et al. 2017).

Heavy metal exposure risk

The HQ results are highly concerning, as they were many orders of magnitude greater than 1. When heavy metals enter the human body, they bind to blood proteins and are distributed to various organs, especially the small intestine, liver, kidneys, brain, and bones; a portion of the ingested concentration is stored in organs, and the rest is excreted through urine and feces (Chávez-Gómez et al. 2017). In addition, pathological symptoms due to heavy metal exposure can be up to 5-fold more severe in children than in adults and may be exacerbated by iron deficiencies (Chávez-Gómez et al. 2017). Thus, many people in northern Mexico are at risk of organ malfunction and other diseases related to heavy metal exposure.

The maximum HQ values calculated for adults and children from the probable maximum concentrations of As, due to the As concentrations in this study, were found to carry a higher risk than that posed by the groundwater samples from India, Pakistan, Bangladesh, Nepal, and Sri Lanka (Badeenezhad et al. 2023); groundwater from Shiraz, Iran (0.426–0.183) (Badeenezhad et al. 2023); and groundwater and surface water from Maharashtra, India (5.5108) (Mawari et al. 2022) (Supplementary Table S7). Most of the estimated non-carcinogenic risks for children and adults exposed to Cd, Pb, and Hg in drinking water in this study were lower than those reported for the same elements in streams, wells, and stormwater in Anambra, Nigeria (1.57–107.56) (Emmanuel et al. 2022). In addition, the Cd, Pb, and Hg risks were higher than those established for these elements in the groundwater of Shiraz, Iran (0.004–0.066) (Badeenezhad et al. 2023) and for groundwater and surface water in Maharashtra, India (0.0033–0.4724) (Mawari et al. 2022). In addition, the HQ values for As, Cd, and Pb in this study were also higher than those for the same elements in the water collected from wells that had been dug outside the houses of CKDu patients from Sri Lanka (0.00169–0.04) (Kulathunga et al. 2022).

The CR values related to As ​​were similar to those reported in groundwater from India, Pakistan, Bangladesh, Nepal, and Sri Lanka (4.95 × 10−4 to 2.10 × 10−1) (Aryan et al. 2024) (Supplementary Table S8). This is especially worrying considering that As has been reported in various parts of northern Mexico. The CR values for Pb were similar to the incremental lifetime cancer risks established for streams, wells, and stormwater in Anambra, Nigeria (5.34 × 10−6 to 5.88 × 10−5) (Emmanuel et al. 2022).

It is crucial to highlight that in northern Mexico, all CR values ​​for As fell within category VII of extremely high risk. In 2022, 89,574 deaths from malignant tumors were reported in Mexico, accounting for 10.6% of all deaths that year (INEGI 2024). In addition, Sonora, Chihuahua, and Nuevo León are the states in northern Mexico with the highest death rates due to malignant tumors (74.71–95.96 deaths per 100,000 inhabitants) (INEGI 2024). In addition, indigenous children in northern Mexico who were exposed to high As through drinking water exhibited significant increases in DNA damage (Maldonado-Escalante et al. 2018); thus, it is possible that the high concentrations of As in drinking water in northern Mexico are related to the high incidence of cancer in the region. It is imperative that adequate monitoring and management measures are implemented for aquatic resources; in addition, water treatment processes must be established, and the quality of public water pipes must be improved.

The present study relied on summary data of varying quality, including minimum, mean, and maximum values, which precluded any distributional development or trend analyses of metal concentrations over the analyzed time period. Comparisons between HQ and CR values were conducted in the context of the major differences in data sources; we recognize that they might not truly reflect the real differences. Possible relationships were conjectured between the results and disease incidence; these relationships were based on results and background, without considering other environmental variables. This article does not provide a probabilistic model that specifies the exact relationships between metal concentrations and disease incidence.

In all, 316 heavy metal concentrations were reported in 55 studies from 1975 to 2021: 174 As concentrations (0 to 11,100 μg/L), 55 Cd concentrations (0 to 5,250 μg/L), 24 Cr concentrations (0 to 350 μg/L), 8 Hg concentrations (0–25 μg/L), and 55 Pb concentrations (0 to 9,470 μg/L). The states of Coahuila and Nuevo León had the highest levels of As contamination over the years included in the study; in addition, Durango and Tamaulipas had the highest recorded Pb concentrations. All reported concentrations exhibited high variability with respect to time and the study area.

Almost half of the reported concentrations (43.04%) exceeded the limits established by the Official Mexican Standard NOM-127-SSA1-1994, and 50.63% exceeded the concentration limits established by the WHO. Moreover, 65.5% of the As concentrations exceeded the concentration limit established by the WHO. The analyses indicated that the probable maximum concentrations were 244.55 μg As/L, 10.00 μg Cd/L, 56.00 μg Cr/L, 1.68 μg Hg/L, and 90.35 μg Pb/L.

The non-carcinogenic risk calculated for As, Cd, Cr, Hg, and Pb resulted in HQ values that ranged from 0.0011 to 404.62. Arsenic exhibited the highest health risks, with extremely high HQ values (>40). The CR values for Cr+6 varied from medium-high (>1 × 10−4 to 5 × 10−4) to extremely high (>1 × 10−3) risk, and CR values for Pb varied from low-medium (>1 × 10−5 to 5 × 10−5) to medium-high (>5 × 10−5 to 1 × 10−4) risk. Furthermore, all CR values for As represented an extremely high risk, with CR values from 2.29 × 10−1 to 7.77 × 10−1. Based on the HQ and CR values, the highest risks were found in children under 6 years of age.

Although not all concentrations exceeded those established by the Official Mexican Standard, most of the associated non-carcinogenic HQ and CR were high, indicating that some standards do not consider the established RfD or CSF values. Therefore, it is necessary to revise the current guidelines and lower the maximum permissible concentrations. In addition, financial support must be increased for research on the distribution, concentrations, and dispersal of heavy metals and their consequences in Mexico. Management actions should also be improved. It is necessary to conduct massive sampling campaigns of various aquatic resources in northern Mexico to determine heavy metal concentrations and evaluate the real risks faced by the population. In particular, heavy metal monitoring programs for water bodies must be expanded, especially those led by RENAMECA. Awareness campaigns must also be created and implemented that educate the public on the risks associated with heavy metal exposure. Lastly, adsorption instruments should be installed in factory effluent outlets, and public purification and distribution systems for drinking water must be improved to remove heavy metals. Collectively, these actions will serve to safeguard the drinking water sources in northern Mexico, lowering the potential health risks associated with the consumption of heavy metals and improving the quality of life of people in the region.

The authors are thankful to Dr Andrea Lievana MacTavish for professionally editing this article.

This study was partially supported by the Universidad Autónoma de Sinaloa through PROFAPI [PRO_A2_025].

M.R.-C.: conceptualization, design, data curation, formal analysis, visualization, and writing – original draft. O.E.-S.: conceptualization, funding acquisition, supervision, validation, and writing – review & editing. M.B.-L.: conceptualization, supervision, validation, and writing – review & editing. M.G.F.-E.: validation and writing – review & editing. N.Y.Z.-A.: validation and writing – review & editing. C.C.O.-M.: validation and writing – review & editing. All authors approved the final version.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Agency for Toxic Substances and Disease Registry (ATSDR)
(
2022a
)
Calculating Hazard Quotients and Cancer Risk Estimates
.
Atlanta, GA, USA
:
U.S. Department of Health and Human Services, Public Health Service
. .
Agency for Toxic Substances and Disease Registry (ATSDR)
(
2022b
)
Estimating Site-Specific Ingestion and Dermal Exposure Doses
.
Atlanta, GA, USA
:
U.S. Department of Health and Human Services, Public Health Service
. .
Alarcón-Herrera
M. T.
,
Montenegro
I.
,
Romero-Navar
P.
,
Martín-Domínguez
I. T.
&
Vázquez
R.
(
2001
)
Contenido de arsénico en el agua potable del Valle del Guadiana en México
,
Ingeniería Hidráulica en México
,
16
(
4
),
63
70
.
Alcázar-Medina
F.
,
Núñez-Núñez
C.
,
Villanueva-Fierro
I.
,
Antileo
C.
&
Proal-Nájera
J.
(
2020
)
Removal of heavy metals present in groundwater from a northern Mexico mining community using Agave tequilana Weber extracts
,
Revista Mexicana De Ingeniería Química
,
19
(
3
),
1187
1199
.
https://doi.org/10.24275/rmiq/Bio1049
.
Archundia
D.
,
Prado-Pano
B.
,
González-Méndez
B.
,
Loredo-Portales
R.
&
Molina-Freaner
F.
(
2021
)
Water resources affected by potentially toxic elements in an area under current and historical mining in northwestern Mexico
,
Environmental Monitoring and Assessment
,
193
(
4
),
236
.
https://doi:10.1007/s10661-021-08998-z
.
Armienta
M. A.
&
Segovia
N.
(
2008
)
Arsenic and fluoride in the groundwater of Mexico
,
Environmental Geochemistry and Health
,
30
(
4
),
345
353
.
https://doi.org/10.1007/s10653-008-9167-8
.
Aryan
Y.
,
Pon
T.
,
Panneerselvam
B.
&
Dikshit
A. K.
(
2024
)
A comprehensive review of human health risks of arsenic and fluoride contamination of groundwater in the south Asia region
,
Journal of Water and Health
,
22
(
2
),
235
267
.
https://doi:10.2166/wh.2023.082
.
Ayala-Rodríguez
L.
(
2010
)
Metales pesados en agua y sedimentos del río Sinaloa
.
MSc thesis
.
National Polytechnic Institute
. .
Badeenezhad
A.
,
Soleimani
H.
,
Shahsavani
S.
,
Parseh
I.
,
Mohammadpour
A.
,
Azadbakht
O.
,
Javanmardi
P.
,
Faraji
H.
&
Nalosi
K. B.
(
2023
)
Comprehensive health risk analysis of heavy metal pollution using water quality indices and Monte Carlo simulation in R software
,
Scientific Reports
,
13
,
15817
.
https://doi.org/10.1038/s41598-023-43161-3
.
Balali-Mood
M.
,
Naseri
K.
,
Tahergorabi
Z.
,
Khazdair
M. R.
&
Sadeghi
M.
(
2021
)
Toxic mechanisms of five heavy metals: mercury, lead, chromium, cadmium, and arsenic
,
Frontiers in Pharmacology
,
12
,
643972
.
https://doi.org/10.3389/fphar.2021.643972
.
Camacho
L.
,
Gutierrez
M.
,
Alarcón-Herrera
M.
,
Villalba
L.
&
Deng
S.
(
2011
)
Occurrence and treatment of arsenic in groundwater and soil in northern Mexico and southwestern USA
,
Chemosphere
,
83
(
3
),
211
225
.
https://doi.org/10.1016/j.chemosphere.2010.12.067
.
Caporale
A. G.
&
Violante
A.
(
2016
)
Chemical processes affecting the mobility of heavy metals and metalloids in soil environments
,
Current Pollution Reports
,
2
,
15
27
.
https://doi.org/10.1007/s40726-015-0024-y
.
Carrillo-Chávez
A.
,
Drever
J. I.
&
Martínez
M.
(
2000
)
Arsenic content and groundwater geochemistry of the San Antonio-El Triunfo, Carrizal and Los Planes aquifers in southernmost Baja California, Mexico
,
Environmental Geology
,
39
(
2000
),
1295
1303
.
https://doi.org/10.1007/S002540000153
.
Cassassuce
F.
,
Olvera
J.
&
Loera-Pizarro
E.
(
2005
)
Estudio del Arsénico en Baja California Sur
.
La Paz, México
:
Sociedad de Historia Natural Niparajá, A.C. en Colaboración con la CONAGUA
.
Sociedad de Historia Natural Niparajá AC. Available at: https://defiendelasierra.org/wp-content/uploads/2015/09/ESTUDIO-ARSENICO-NIPARAJA.pdf [Accessed 21st November 2021]
.
Chávez-Gómez
N. L.
,
Cabello-López
A.
,
Gopar-Nieto
R.
,
Aguilar-Madrid
G.
,
Marín-López
K. S.
,
Aceves-Valdez
M.
,
Jiménez-Ramírez
C.
,
Cruz-Angulo
M. C.
&
Juárez-Pérez
C. A.
(
2017
)
Enfermedad renal crónica en México y su relación con los metales pesados
,
Revista Médica del Instituto Mexicano del Seguro Social
,
55
(
6
),
725
734
.
CONAGUA [Comisión Nacional del Agua]
(
2018
)
Atlas del Agua en México
.
Mexico City: Comisión Nacional del Agua
. .
CONAGUA [Comisión Nacional del Agua]
(
2019
)
Estadísticas del Agua en México 2019
.
Mexico City: Comisión Nacional del Agua
.
Available at: http://sina.conagua.gob.mx/publicaciones/EAM_2019.pdf [Accessed 23rd March 2021]
.
da Silva
E. B.
,
Gao
P.
,
Xu
M.
,
Guan
D.
,
Tang
X.
&
Ma
L. Q.
(
2020
)
Background concentrations of trace metals As, Ba, Cd, Co, Cu, Ni, Pb, Se, and Zn in 214 Florida urban soils: different cities and land uses
,
Environmental Pollution
,
264
,
114737
.
https://doi.org/10.1016/j.envpol.2020.114737
.
Diario Oficial de la Federación (DOF)
(
2022
)
NOM-127-SSA1-2021. Agua para uso y consumo humano. Límites permisibles de la calidad del agua
.
Mexico City: Norma Oficial Mexicana
. .
Díaz-Barriga
F.
,
Santos
M. A.
,
Mejía
J. J.
,
Batres
L.
,
Yáñez
L.
,
Carrizales
L.
,
Vera
E.
,
del Razo
L. M.
&
Cebrián
M. E.
(
1993
)
Arsenic and cadmium exposure in children living near a smelter complex in San Luis Potosí, Mexico
,
Environmental Research
,
62
(
2
),
242
250
.
https://doi.org/10.1006/enrs.1993.1109
.
Dirección General de Normas (DGN)
(
1994
)
NOM-127-SSA-1994. Norma oficial mexicana. Salud ambiental, agua para uso y consumo humano, límites permisibles de calidad y tratamientos a que debe someterse el agua para su potabilización
.
Mexico City: Diario Oficial de la Federación
.
Available at: http://www.ordenjuridico.gob.mx/Documentos/Federal/wo69543.pdf [Accessed 19th November 2021]
.
Echeverry
G.
,
Zapata
A.
,
Páez
M.
,
Méndez
F.
&
Peña
M.
(
2015
)
Valoración del riesgo en salud en un grupo de población de Cali, Colombia, por exposición a plomo, cadmio, mercurio, ácido 2,4-diclorofenoxiacético y diuron, asociada al consumo de agua potable y alimentos
,
Biomédica
,
35
(
2
),
110
119
.
https://dx.doi.org/10.7705/biomedica.v35i0.2464
.
Emmanuel
U. C.
,
Chukwudi
M. I.
,
Monday
S. S.
&
Anthony
A. I.
(
2022
)
Human health risk assessment of heavy metals in drinking water sources in three senatorial districts of Anambra State, Nigeria
,
Toxicology Reports
,
9
,
869
875
.
https://doi.org/10.1016/j.toxrep.2022.04.011
.
Escamilla-Rodríguez
A.
,
Carlos-Hernández
S.
&
Díaz-Jiménez
L.
(
2021
)
Evidence of resistance of heavy metals from bacteria isolated from natural waters of a mining area in Mexico
,
Water
,
13
(
19
),
2766
.
http://dx.doi.org/10.3390/w13192766
.
Ferreccio
C.
&
Sancha
A. M.
(
2006
)
Arsenic exposure and its impact on health in Chile
,
Journal of Health, Population, and Nutrition
,
24
(
2
),
164
175
.
Flores-Laureano
J. S.
(
1997
)
Evaluación de la Calidad del Agua del Rio San Juan, en el Estado de Nuevo León
.
MSc thesis
.
University of Nuevo León
.
Available at: http://eprints.uanl.mx/7699/1/1020118308.PDF (Accessed: 11 September 2021)
.
Frank
J. J.
,
Poulakos
A. G.
,
Tornero-Velez
R.
&
Xue
J.
(
2019
)
Systematic review and meta-analyses of lead (Pb) concentrations in environmental media (soil, dust, water, food, and air) reported in the United States from 1996 to 2016
,
Science of the Total Environment
,
694
,
133489
.
https://doi.org/10.1016/j.scitotenv.2019.07.295
.
Frías-Espericueta
M. G.
,
Mejía-Cruz
R.
,
Osuna-López
I.
,
Muy-Rangel
M. D.
,
Rubio-Carrasco
W.
,
Aguilar-Juárez
M.
&
Voltolina
D.
(
2014
)
Metal discharges by Sinaloa rivers to the coastal zone of NW Mexico
,
Bulletin of Environment Contamination and Toxicology
,
92
(
2
),
132
136
.
https://doi: 10.1007/s00128-013-1158-5
.
Frías-Sarmiento
E.
(
2007
)
Financiamientos para la agricultura comercial de Sinaloa: 1932–1949. El creciente papel de los actores privados regionales y estadounidenses
,
Región y Sociedad
,
19
(
39
),
135
158
.
Gamboa-Loira
B.
,
Cebrián
M. E.
&
López-Carrillo
L.
(
2020
)
Arsenic exposure in northern Mexican women. Exposición a arsénico en mujeres del norte de México
,
Salud Pública de México
,
62
(
3
),
262
269
.
https://doi.org/10.21149/11085
.
García-Rico
L.
,
Meza-Figueroa
D.
,
Jay-Gandolfi
A.
,
Del Rivero
C. I.
,
Martínez-Cinco
M. A.
&
Meza-Montenegro
M. M.
(
2019
)
Health risk assessment and urinary excretion of children exposed to arsenic through drinking water and soils in Sonora, Mexico
,
Biological Trace Element Research
,
187
(
1
),
9
21
.
https://doi.org/10.1007/s12011-018-1347-5
.
González
F. J. A.
,
López
E. M. R.
,
Saldaña
M. C. M.
,
Barrera
A. L. G.
,
Juárez
F. J.
&
Sánchez
J. L. R.
, (
2012
)
Water quality in the state of Aguascalientes and its effects on the population's health
. In:
Oswald Spring
Ú
(ed.)
Water Resources in Mexico. Hexagon Series on Human and Environmental Security and Peace
, Vol.
7
.
Heidelberg, Berlin: Springer
, pp.
217
218
.
González-Dávila
O.
(
2013
)
Groundwater Contamination and Contingent Valuation of Safe Drinking Water in Guadalupe, Zacatecas, Mexico
.
Working Papers 180
.
UK
:
Department of Economics, SOAS, University of London
. .
González-Horta
C.
,
Ballinas-Casarrubias
L.
,
Sánchez-Ramírez
B.
,
Ishida
M. C.
,
Barrera-Hernández
A.
,
Gutiérrez-Torres
D.
,
Zacarias
O. L.
,
Saunders
R. J.
,
Drobná
Z.
,
Mendez
M. A.
,
García-Vargas
G.
,
Loomis
D.
,
Stýblo
M.
&
Del Razo
L. M.
(
2015
)
A concurrent exposure to arsenic and fluoride from drinking water in Chihuahua, Mexico
,
International Journal of Environmental Research and Public Health
,
12
(
5
),
4587
4601
.
https://doi.org/10.3390/ijerph120504587
.
Guzmán-Colis
G.
,
Thalasso
F.
,
Ramírez-López
E. M.
,
Rodríguez-Narciso
S.
,
Guerrero-Barrera
A. L.
&
Avelar-González
F. J.
(
2011
)
Evaluación espacio-temporal de la calidad del agua del río San Pedro en el Estado de Aguascalientes, México
,
Revista Internacional de Contaminación Ambiental
,
27
(
2
),
89
102
.
Health Canada
(
2010
)
Toxicological Review and Recommended Toxicological Reference Values for Environmental Lead Exposure in Canada
. Ontario: Health Canada. .
Hirano
T.
(
1990
)
A theoretical approach to hazard assessments
,
Journal of Loss Prevention in the Process Industries
,
3
(
3
),
321
326
.
https://doi.org/10.1016/0950-4230(90)80026-7
.
Instituto Nacional de Estadística y Geografía
(
2018
)
Censo Nacional 2018
.
Aguascalientes City: Gobierno Federal
.
Available at: https://www.inegi.org.mx/programas/cngf/2018/ [Accessed 19th November 2021]
.
Instituto Nacional de Estadística y Geografía
(
2024
)
Estadísticas a propósito del día mundial contra el cáncer
.
Comunicado 78/24
.
Aguascalientes City: Gobierno Federal
. .
Jaishankar
M.
,
Tseten
T.
,
Anbalagan
N.
,
Mathew
B. B.
&
Beeregowda
K. N.
(
2014
)
Toxicity, mechanism and health effects of some heavy metals
,
Interdisciplinary Toxicology
,
7
(
2
),
60
72
.
https://doi.org/10.2478/intox-2014-0009
.
Jara-Marini
M. E.
,
Tapia-Alcaraz
J. N.
,
Dumer-Gutiérrez
J. A.
,
García-Rico
L.
,
García-Hernández
J.
&
Páez-Osuna
F.
(
2013
)
Distribution and accumulation of Cd, Cu, Hg, Pb and Zn in the surface sediments of El Tobari Lagoon, central-East Gulf of California: an ecosystem associated with agriculture and aquaculture activities
,
Journal of Environmental Science and Health. Part A, Toxic/Hazardous Substances & Environmental Engineering
,
48
(
14
),
1842
1851
.
https://doi.org/10.1080/10934529.2013.823341
.
Jomova
K.
,
Alomar
S. Y.
,
Nepovimova
E.
,
Kuca
K.
&
Valko
M.
(
2024
)
Heavy metals: toxicity and human health effects
,
Archives of Toxicology
,
98
(
12
), 153–209.
https://doi.org/10.1007/s00204-024-03903-2
.
Khan
M. M.
,
Sakauchi
F.
,
Sonoda
T.
,
Washio
M.
&
Mori
M.
(
2003
)
Magnitude of arsenic toxicity in tube-well drinking water in Bangladesh and its adverse effects on human health including cancer: evidence from a review of the literature
,
Asian Pacific Journal of Cancer Prevention
,
4
(
1
),
7
14
.
Kulathunga
M. R. D. L.
,
Wijayawardena
M. A. A.
,
Naidu
R.
,
Wimalawansa
S. J.
&
Rahman
M. M.
(
2022
)
Health risk assessment from heavy metals derived from drinking water and rice, and correlation with CKDu
,
Frontiers in Water
,
3
,
786487
.
https://doi.org/10.3389/frwa.2021.786487
.
Londoño-Franco
L. F.
,
Londoño-Muñoz
P. T.
&
Muñoz-García
F. G.
(
2016
)
Los riesgos de los metales pesados en la salud humana y animal
,
Rev Bio Agro
,
14
(
2
),
145
153
.
https://doi.org/10.18684/BSAA(14)145-153
.
López
E.
,
Reynaldo-Patiño
M. L.
,
Vázquez-Sauceda
R.
,
Pérez-Castañeda
L. U.
,
Arellano-Méndez
L.
,
Ventura
R.
&
Heyer
L.
(
2020
)
Water quality and ecological risk assessment of intermittent streamflow through mining and urban areas of San Marcos River sub-basin, Mexico
,
Environmental Nanotechnology, Monitoring & Management
,
14
(
2020
),
100369
.
https://doi.org/10.1016/j.enmm.2020.100369
.
López-Álvarez
B.
,
Ramos-Leal
J. A.
,
Moran-Ramírez
J.
,
Cardona-Benavides
A.
&
Hernández-García
G.
(
2013
)
Origen de la calidad del agua del acuífero colgado y su relación con los cambios de uso de suelo en el Valle de San Luis Potosí
,
Boletín de la Sociedad Geológica Mexicana
,
65
(
1
),
9
26
.
Maldonado-Escalante
J. F.
,
Meza-Figueroa
D.
,
Figueroa
A. G.
,
García-Rico
L.
,
Burgess
J. L.
,
Lantz
R. C.
,
Yáñez-Estrada
L.
,
Martínez-Cinco
M. A.
,
Balderas-Cortés
J. D. J.
,
Mondaca-Fernández
I.
&
Meza-Montenegro
M. M.
(
2018
)
An integrated health risk assessment of indigenous children exposed to arsenic in Sonora, Mexico
,
Human and Ecological Risk Assessment
,
25
(
3
),
706
772
.
https://doi.org/10.1080/10807039.2018.1449098
.
Mawari
G.
,
Kumar
N.
,
Sarkar
S.
,
Frank
A. L.
,
Daga
M. K.
,
Singh
M. M.
,
Joshi
T. K.
&
Singh
I.
(
2022
)
Human health risk assessment due to heavy metals in ground and surface water and association of diseases with drinking water sources: a study from Maharashtra, India
,
Environmental Health Insights
,
16
,
11786302221146020
.
https://doi.org/10.1177/11786302221146020
.
McClintock
T. R.
,
Chen
Y.
,
Bundschuh
J.
,
Oliver
J. T.
,
Navoni
J.
,
Olmos
V.
,
Lepori
E. V.
,
Ahsan
H.
&
Parvez
F.
(
2012
)
Arsenic exposure in Latin America: biomarkers, risk assessments and related health effects
,
Science of the Total Environment
,
429
(
2012
),
76
91
.
https://doi.org/10.1016/j.scitotenv.2011.08.051
.
Mendoza-Cano
O.
,
Sánchez-Piña
R. A.
,
Barrón-Quintana
J.
,
Cuevas-Arellano
H. B.
,
Escalante-Minakata
P.
&
Solano-Barajas
R.
(
2017
)
Riesgos potenciales de salud por consumo de agua con arsénico en Colima, México
,
Salud Pública de México
,
59
(
1
),
34
40
.
https://doi.org/10.21149/8413
.
Meza
M. M.
,
Kopplin
M. J.
,
Burgess
J. L.
&
Gandolfi
A. J.
(
2004
)
Arsenic drinking water exposure and urinary excretion among adults in the Yaqui Valley, Sonora, Mexico
,
Environmental Research
,
96
(
2
),
119
126
.
https://doi.org/10.1016/j.envres.2003.08.010
.
Meza
M. M.
,
Yu
L.
,
Rodríguez
Y. Y.
,
Guiad
M.
,
Thompson
D.
,
Gandolfi
J.
&
Klimecki
W. T.
(
2005
)
Developmentally restricted genetic determinants of human arsenic metabolism: association between urinary methylated arsenic and CYT19 polymorphisms in children
,
Environmental Health Perspectives
,
113
(
6
),
775
781
.
https://doi.org/10.1289/ehp.7780
.
National Institutes of Health
(
2017
)
OPA Excel Tips: Creating a Box and Whisker Plot in Excel
.
Bethesda, MD: Office of Portfolio Analysis
. .
Ochoa-Reyes
J.
,
Vidal-Lozano
A.
,
Lerma-Santana
A.
,
Gómez-Reza
L.
,
Reta-Sánchez
A.
&
Fernández-Escajeda
D.
(
2009
)
Arsenicismo subagudo y crónico en una población rural
,
Dermatología Revista Mexicana
,
53
(
2
),
63
69
.
Office of Environmental Health Hazard Assessment (OEHHA)
(
2006
)
Update of the Public Health Goal for Cadmium
.
Sacramento, CA: OEHHA
. .
Office of Environmental Health Hazard Assessment (OEHHA)
(
2011a
)
Lead and Lead Compounds
.
Sacramento, CA: OEHHA
.
Available at: https://oehha.ca.gov/chemicals/lead-and-lead-compounds [Accessed 23rd March 2021]
.
Office of Environmental Health Hazard Assessment (OEHHA)
(
2011b
)
Public Health Goals for Chemicals in Drinking Water: Hexavalent Chromium (Cr VI)
.
Sacramento, CA: OEHHA
. .
Pacheco-Treviño, S. & Manzano-Camarillo, M. G. F. (2024) Review of water scarcity assessments: Highlights of Mexico's water situation, WIREs Water, 11 (4), e1721. https://doi.org/10.1002/wat2.1721
Page
M. J.
,
McKenzie
J. E.
,
Bossuyt
P. M.
,
Boutron
I.
,
Hoffmann
T. C.
,
Mulrow
C. D.
&
Moher
D.
(
2021
)
The PRISMA 2020 statement: an updated guideline for reporting systematic reviews
,
BMJ
,
10
(
1
),
1
11
.
Pérez-Castresana
G.
,
Castañeda-Roldán
E.
,
García-Suastegui
W. A.
,
Morán-Perales
J. L.
,
Cruz-Montalvo
A.
&
Handal-Silva
A.
(
2019
)
Evaluation of health risks due to heavy metals in a rural population exposed to Atoyac river pollution in Puebla, Mexico
,
Water
,
11
(
2
),
277
.
https://doi.org/10.3390/w11020277
.
Proshad
R.
,
Chandra
K.
,
Akash
P. B.
,
Ritu
S. A.
,
Khan
M. S. U.
,
Dey
H. C.
,
Sarker
J.
,
Baroi
A.
&
Idris
A. M.
(
2025
)
An overview of metal(oid) pollution, sources, and probabilistic health risk evaluations based on a Monte Carlo simulation of surface river water in a developing country
,
Water
,
17
(
5
),
630
.
https://doi.org/10.3390/w17050630
.
Ramal
M. M.
,
Abdulhameed
U. H.
&
Jalal
A. D.
(
2021
)
Trace elements risk assessment in taps drinking water of Ramadi City, Anbar Province, Iraq
,
International Journal of Safety and Security Engineering
,
11
(
6
),
623
634
.
https://doi.org/10.18280/ijsse.110602
.
Ramirez-Andreotta
M.
&
Superfund Research Program (Superfund)
(
2016
)
Understanding Environmental Health Risk Assessment
.
Tucson, AZ: Arizona Department of Health Services – Web
. .
Ravalli
F.
,
Yu
Y.
,
Bostick
B. C.
,
Chillrud
S. N.
,
Schilling
K.
,
Basu
A.
,
Navas-Acien
A.
&
Nigra
A. E.
(
2022
)
Sociodemographic inequalities in uranium and other metals in community water systems across the USA, 2006-11: a cross-sectional study
,
The Lancet Planetary Health
,
6
(
4
),
e320
e330
.
https://doi.org/10.1016/S2542-5196(22)00043-2
.
Rivera-Hernández
J. R.
,
Green-Ruiz
C. R.
,
Pelling-Salazar
L. E.
&
Russell-Flegal
A.
(
2021
)
Monitoring of As, Cd, Cr, and Pb in groundwater of Mexico's agriculture Mocorito river aquifer: implications for risks to human health
,
Water, Air, & Soil Pollution
,
232
,
291
.
https://doi.org/10.1007/s11270-021-05238-5
.
Saldarriaga
H.
,
De la Garza
I.
,
Waliszewski
S.
,
Colunga
E.
,
Amador
O.
,
Moreno-Dávila
M.
&
Morales-Cueto
R.
(
2014
)
Chemical evaluation of groundwater from supply wells in the state of Coahuila, México
,
Journal of Water Resource and Protection
,
6
,
49
54
.
https://doi.org/10.4236/jwarp.2014.61008
.
Salinas-Rodríguez
S. A.
&
Martínez Pacheco
A. I.
(
2024
)
Water reserves for the environment: a strategic and temporal analysis (2012–2022) for the implementation of environmental flows in Mexico
,
Diversity
,
16
(
3
),
190
.
https://doi.org/10.3390/d16030190
.
Sandoval-Alvarado
M.
,
Conejo-Flores
R.
,
Guzmán-Pantoja
J.
&
García
J.
(
2020
)
Remoción de arsénico presente en el agua potable del municipio de Ojocaliente, Zac. empleando perlas de gel de quitosano
,
Revista de Desarrollo Urbano y Sustentable
,
6
(
18
),
18
27
.
https://doi.org/10.35429/JUSD.2020.18.6.18.27
.
Secretaría de Medio Ambiente y Recursos Naturales (SEMARNAT)
(
2020
)
Decreto por el que se Aprueba el Programa Nacional Hídrico 2020–2024
.
Mexico City: Diario Oficial de la Federación – Web
. .
Servicio Geológico Mexicano (SGM)
(
2019
)
Panorama Minero del Estado de Sinaloa
.
Hidalgo: SGMWeb
. .
Singh
A.
&
Kostova
I.
(
2024
)
Health effects of heavy metal contaminants vis-à-vis microbial response in their bioremediation
,
Inorganica Chimica Acta
,
568
,
122068
.
https://doi.org/10.1016/j.ica.2024.122068
.
Sobsey
M. D.
,
Stauber
C. E.
,
Casanova
L. M.
,
Brown
J. M.
&
Elliott
M. A.
(
2008
)
Point of use household drinking water filtration: a practical, affective solution for proving sustained access to safe drinking water in the developing world
,
Environmental Science & Technology
,
42
,
4261
4267
.
https://doi: 10.1021/es702746n
.
Spektor
D. M.
(
1998
)
A Review of the Scientific Literature as It Pertains to Gulf War Illnesses
.
Oil Well Fires. Chapter Two: Environmental Measurements
.
Santa Monica, CA: National Defense Research Institute
.
Available at: https://www.gulflink.osd.mil/library/rowl/mr1018ch2.html [Accessed 20th February 2025]
.
Thatcher
R. W.
,
Lester
M. L.
,
McAlaster
R.
&
Horst
R.
(
1982
)
Effects of low levels of cadmium and lead on cognitive functioning in children
,
Archives of Environmental & Occupational Health
,
37
(
3
),
159
166
.
https://doi.org/10.1080/00039896.1982.10667557
.
U.K. Technical Advisory Group on the Water Framework Directive (UKTAG)
(
2020
)
Paper 11b(i). Groundwater Chemical Classification for the Purposes of the Water Framework Directive and the Groundwater Daughter Directive
.
Strasbourg: Water Framework Directive Web
. .
U.S. Environmental Protection Agency
(
1995
)
Arsenic, Inorganic; CASRN 7440-38-2
.
Washington, DC: USEPAWeb
.
Available at: https://iris.epa.gov/static/pdfs/0278_summary.pdf [Accessed 15th March 2021]
.
U.S. Environmental Protection Agency
(
2000
)
Chromium Compounds. Hazard Summary
. Washington, DC:
USEPAWeb
. .
U.S. Environmental Protection Agency
(
2001
)
Risk Assessment Guidance for Superfund. Volume I. Part A, Process for Conducting Probabilistic Risk Assessment
.
Washington, DC: USEPAWeb
. .
U.S. Environmental Protection Agency
(
2005
)
Guidelines for Carcinogen Risk Assessment
.
Washington, DC: USEPAWeb
. .
U.S. Environmental Protection Agency
(
2012
)
Basic Information About the Integrated Risk Information System. Appendix Q – Human Health Benchmarks
.
Washington, DC: USEPAWeb
.
Available at: https://archive.epa.gov/epawaste/hazard/web/pdf/appndxq-w.pdf [Accessed 15th September 2021]
.
U.S. Environmental Protection Agency
(
2014
)
Priority Pollutant List
. Washington, DC: USEPAWeb .
U.S. Environmental Protection Agency
(
2021
)
Targeted National Sewage Sludge Survey (TNSSS): Summary Statistics and Estimates of 95th Percentiles for 84 Additional Analytes
.
Washington, DC: USEPAWeb Available at: https://www.epa.gov/sites/default/files/2021-04/documents/tnsss-appendix-95th-percentiles-report.pdf [Accessed 20th February 2025]
.
U.S. Environmental Protection Agency
(
2024a
)
Chromium (VI)
.
Integrated Risk Information System
. Washington, DC: USEPAWeb
Available at: https://iris.epa.gov/ChemicalLanding/&substance_nmbr=144 [Accessed 20th February 2025]
.
U.S. Environmental Protection Agency
(
2024b
)
Toxic and Priority Pollutants Under the Clean Water Act
. Washington, DC: USEPAWeb .
U.S. Environmental Protection Agency
(
2025
)
Arsenic, Inorganic
.
Integrated Risk Information System
.
Washington, DC: USEPAWeb Available at: https://iris.epa.gov/ChemicalLanding/&substance_nmbr=278 [Accessed 20th February 2025]
.
Westerhoff
P.
,
Esparza-Soto
M.
,
Mata
P.
,
Parry
W.
&
Johnson
W.
(
2004
)
Drinking Water Quality in the US-Mexico Border Region
. Tempe, AZ: Arizona State University. .
World Health Organization
(
2017
)
Guidelines for Drinking-Water Quality, Incorporating the 1st Addendum
, 4th ed
.
Geneva: WHOWeb
.
Available at: https://www.who.int/publications/i/item/9789241549950 [Accessed 15th June 2021]
.
World Health Organization
(
2020
)
10 Chemicals of Public Health Concern
.
Geneva: WHOWeb
. .
Zhang
P.
,
Yang
M.
,
Lan
J.
,
Huang
Y.
,
Zhang
J.
,
Huang
S.
,
Yang
Y.
&
Ru
J.
(
2023
)
Water quality degradation due to heavy metal contamination: health impacts and eco-friendly approaches for heavy metal remediation
,
Toxics
,
11
(
10
),
828
.
https://doi.org/10.3390/toxics11100828
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data