Hg Pollution in Groundwater of Andean Region of Ecuador and Human Health Risk Assessment

: In Ecuador, illegal gold mining has led to significant environmental impacts, with the release of harmful elements such as mercury (Hg) into the environment. Mercury, due to its physical–chemical characteristics and the transport elements involved between different environmental matrices, can easily percolate through the soil and reach groundwater. The purpose of this study was to evaluate the mercury concentration levels in the Andean region in order to perform a human health risk assessment. For this purpose, 175 water samples were analyzed, of which 9.71% exceeded the maximum permissible limit (MPL) established for drinking water in accordance with Ecuadorian regulations. The risk analysis was conducted by applying two approaches: deterministic and probabilistic. The deterministic approach involves a specific analysis based on the calculation of the risk quotient (HQ) and risk index (HI) for both receptors (adults and children) and scenarios (residential and recreational) considered; the probabilistic approach is based on the use of stochastic simulation techniques. The results obtained from the two approaches show a discrepancy, with the deterministic analysis providing more conservative results; however, they coincide in showing higher risk for the child population; decision-makers could use these results to identify areas to be monitored and plan more detailed investigation plans.


Introduction
Mining is the main source of precious metals and minerals worldwide [1].Currently, the world market is experiencing exponential growth, so the quantities extracted have increased exponentially [2].In particular, gold (Au) is one of the metals most in demand in global trade and is, therefore, the most exploited metal in mining production on a global scale [3].Gold mining can be conducted at multiple scales, but ASGM is widespread in developing countries.In fact, many South American countries are significant producers and exporters of gold in the world.Among them are Colombia, Peru, Chile, Brazil, and Argentina.In recent years, Ecuador has also become competitive in this sector as a result of a strategy to diversify its production output [4,5].In 2018, gold (Au) production in Ecuador reached 11.5 tons, placing it in tenth place for production among South American countries [6].
The central problem with ASGM is that it is often carried out outside the existing regulations and applies unsafe techniques that involve the use of Hg in Au recovery processes [7], despite such techniques being banned in 2010 [8].The extensive use of Hg for the extraction of precious metals can be traced to the ease with which it amalgamates with them; in particular, in Au mining processes, Hg allows its separation from other heavy minerals: in Au mining, Hg is spread on its surface, forming a substrate that is removed at a later stage, allowing Au of high purity to be obtained via distillation.In terms of its general properties, Hg is classified in the category of transition metals, is a water-insoluble heavy metal, is the only one known to date that occurs in a liquid state under standard conditions, and tends to bioaccumulate, and is therefore significantly persistent.The World Health Organization (WHO) and US EPA [9], such as other states and federal agencies, have defined Hg as one of the most harmful metals for human health and the environment; although it is an element that occurs naturally in the Earth's crust, it is estimated that between 50 and 75% of Hg found in the environment comes from anthropogenic sources [10][11][12].Hg exists in three different forms: elemental Hg, organic Hg compounds, and inorganic Hg compounds; each form is toxic and can result in different types of human health effects.Although it often enters the environment in inorganic form from anthropogenic sources, natural processes can occur that can convert it into an organic form, among which methylmercury and dimethylmercury are considered the most dangerous to human health [13].
According to the methodology of risk analysis, the possible migration routes by which a contaminant may reach receptors are influenced by the characteristics of the analyzed environment, and also by the dominant species present there [14].In addition, the reference scenario also has an influence.The scenarios that are considered in a risk analysis make it possible to identify the types of human receptors present according to land use and the main exposure routes according to the activities/uses and habits that are carried out.The residential scenario refers to living contexts and home environments, where both adult and child human receptors are present.In the residential scenario, receptors are in frequent contact with the contaminants, so exposure is daily and long-term.
The main activities that can bring receptors into contact with the contaminant, in a residential-type scenario, are therefore the accidental ingestion of contaminated water [15] and also the consumption of some particular fish species.The recreational scenario identifies all those spaces intended for playful activities such as hunting, fishing, sports, etc.These are activities that can, in general, be considered occasional.Therefore, the human receptors present in these environments, which may be either adults or children, are exposed to contamination for limited intervals of time.All these activities involve contact with the environment, so one of the most significant modes of exposure is dermal contact.However, the recreational scenario encompasses a wide range of different activities, so a more specific analysis could detect, with reference to some of these activities, a higher and more frequent exposure of the child population.Finally, it should be noted that with reference to work contexts in which Hg amalgam processes take place, the most frequent route of exposure is inhalation of Hg vapors [13,14].
The Andean region of Ecuador has historically been one of the most exploited for mining.Taking into account current mining concessions and ongoing projects, it is estimated that more than 28% of the territory could be occupied for mining activities in the near future [8].Currently, the Andean region has around 797 gold mining concessions distributed in nine provinces, as registered in the Mining Cadastre Web Geoportal of the Mining Regulation and Control Agency (ARCOM), a company that operates in the mining and metals sector, and in Ecuador manages operations related to regulation, mining concessions, production, and illegal mining.These concessions correspond to different levels of activity: artisanal, small, medium, and large-scale mining, and there are also mining activities belonging to the general regime.However, in the Andean region of Ecuador, there is a marked presence of illegal mining operations.Following several complaints, between January and October 2019, the ARCOM, in collaboration with the Mining Crimes Unit of the National Police and the Armed Forces, carried out a total of 418 operations to combat illegal mining, achieving a 60 percent reduction compared to 2018.However, the health emergency from COVID-19 diverted the government's attention from this issue [2].From a social point of view, the Andean population shows a divide between people in favor and people against; in fact, the region's economy is based on activities such as agriculture, livestock, trade, and fishing, so gold mining is still a mainstay of the economy [15].The purpose of this study was to investigate the concentration of (Hg) in groundwater in the Andean region of Ecuador to assess the potential risks to the environment and human health by applying an assessment approach that combined deterministic and probabilistic methods.

Study Area
The research focused on the provinces of Carchi, Imbabura, Pichincha, Cotopaxi, Chimborazo, Bolívar, Cañar, Azuay, and Loja (Figure 1), the populations of which amount to approximately 6,081,342 inhabitants.The economy in this region is mainly based on primary production (agriculture, livestock, etc.) and mining.According to the Mining Cadaster prepared by ARCOM, mining is the main source of contamination in the region because of the uncontrolled use of Hg in gold mining [2].The Andean region takes its name from the Andes, an imposing mountain barrier characterized by extremely steep slopes and altitudes of up to 3500-4000 m, decreasing from the north of the region to the south.In the north, the landscape consists of two parallel mountain ranges with average altitudes between 4000 and 4500 m, separated by a depression with altitudes between 1600 and 3000 m.The Western Cordillera consists of the Cretaceous volcanic complex, while the Eastern Cordillera develops on metamorphic facies formed by the Andean orogeny.These chains are crowned by two parallel rows of large volcanoes that form the "Avenue of Volcanoes", culminating in the active Cotopaxi at 5897 m and the extinct Chimborazo at 6310 m [16].Due to its great geomorphological variety, along with the presence of volcanoes, the region has great advantages in terms of mineral deposits [8].

Sampling and Laboratory Analysis
The analysis started with field data collection work performed by researchers along the territory during the period between March and July 2022, and the choice of sampling points was based on the distribution, along the territory, of mining activities.For this reason, the sampling network is denser in some areas than in others, in particular, the provinces with the most sampling were Azuay with just under 26 percent and Loja with just under 22 percent.For the purpose of the analysis, water samples were taken from wells intended for human consumption, and then the study focused on the analysis of Hg concentrations in groundwater in the region.To collect and transport samples, 250 mL amber flasks were used and acidified with 0.10 mL of nitric acid.The samples were then delivered to the scientific laboratory at the Escuela Superior Politécnica de Chimborazo in Orellana, Ecuador, adhering to a strict custody protocol.Mercury concentrations were measured using Atomic Absorption Spectrophotometry, with a precision of ±5 and a measurement range of 0.0005 to 10 mg/L.The analysis followed the Standard Methods, 23rd Edition, 2017, 3112B-Acid Digestion: EPA Method 3015, 2007 [17,18].Samples were prepared using the nitric acid digestion procedure outlined in EPA Method 7473 [19].The entire workflow, from sample collection to laboratory processing, adhered to international standards for quality, confidentiality, and ethical practices.
Application of the method for Hg concentration determination was preceded by a quality assurance and quality control procedure in accordance with EPA Method 7473 [20].Analyses were carried out with scrupulous use of the appropriate safety devices to avoid any possible contact between operators and Hg compounds, and quality control of the water samples was also carried out to detect organic content.Initially, the laboratory demonstrated its competence by following the sample preparation and analysis procedures according to the established method.For each group of samples processed, a sample

Sampling and Laboratory Analysis
The analysis started with field data collection work performed by researchers along the territory during the period between March and July 2022, and the choice of sampling points was based on the distribution, along the territory, of mining activities.For this reason, the sampling network is denser in some areas than in others, in particular, the provinces with the most sampling were Azuay with just under 26 percent and Loja with just under 22 percent.For the purpose of the analysis, water samples were taken from wells intended for human consumption, and then the study focused on the analysis of Hg concentrations in groundwater in the region.To collect and transport samples, 250 mL amber flasks were used and acidified with 0.10 mL of nitric acid.The samples were then delivered to the scientific laboratory at the Escuela Superior Politécnica de Chimborazo in Orellana, Ecuador, adhering to a strict custody protocol.Mercury concentrations were measured using Atomic Absorption Spectrophotometry, with a precision of ±5 and a measurement range of 0.0005 to 10 mg/L.The analysis followed the Standard Methods, 23rd Edition, 2017, 3112B-Acid Digestion: EPA Method 3015, 2007 [17,18].Samples were prepared using the nitric acid digestion procedure outlined in EPA Method 7473 [19].The entire workflow, from sample collection to laboratory processing, adhered to international standards for quality, confidentiality, and ethical practices.
Application of the method for Hg concentration determination was preceded by a quality assurance and quality control procedure in accordance with EPA Method 7473 [20].Analyses were carried out with scrupulous use of the appropriate safety devices to avoid any possible contact between operators and Hg compounds, and quality control of the water samples was also carried out to detect organic content.Initially, the laboratory demonstrated its competence by following the sample preparation and analysis procedures according to the established method.For each group of samples processed, a sample blank and a laboratory control sample were carried out.The method blank is prepared with reagent water in the specified volume and processed according to the appropriate steps.If the method blank does not contain the target analyte in an amount that affects the quality of the design, it is considered acceptable.In the absence of concrete evidence, this implies that a maximum level of 10% of the lowest concentration in the sample is considered acceptable.Laboratory control sampling should be performed at the project level for each analyte of interest.
The application of the method for the determination of the concentration of (Hg) was carried out through the quality assurance process and a quality control protocol according to EPA Method 7473 [16].The operators performed the measurements while wearing the necessary safety equipment to prevent any form of potential contact with mercury compounds (ingestion, inhalation, and skin contact).In addition, a quality check was carried out on the samples to be analyzed to detect any samples with high organic content, as these require special attention during analysis and need to be reduced in size to avoid ignition in the decomposition tube.The laboratory initially demonstrated its competence by adhering to the sample preparation and analysis procedures as per the established method.For each batch of processed samples, a method blank and a laboratory control sample were included.The method blank was prepared using a reagent in water at the specified volume and was considered acceptable if it did not contain the target analyte at a level that would compromise analysis quality.In the absence of specific guidelines, a maximum level of 10% of the lowest sample concentration was deemed acceptable.Generally, it is necessary to collect laboratory control samples for each analyte of interest at the project level.Acceptance criteria can be based on historical data or, if such data are unavailable, the limit can be set at ±20% of the spiked value.Finally, for quality assurance/quality control (QA/QC) blanks, standards, and duplicated samples were analyzed.The recovery rates were in the range of 75-110%.

Risk Assessment and Characterization
For the health risk assessment, specifically for human health, two possible scenarios were considered: residential and recreational.In both scenarios, the main ways in which people can come into contact with contamination that should be considered are water consumption and direct skin contact.The potential risk associated with human health was assessed through the Hazard Quotient (HQ), which is a risk parameter used for noncarcinogenic substances and expresses how many times the maximum tolerable daily intake of contaminant, per unit of body weight, is exceeded.For each form of exposure, it is estimated and is defined as the ratio of the average daily dose ADD to the reference dose RfD.The average daily dose was determined using Equations ( 1) and (2) below, USEPA 2001 and 2004, respectively.Organic Hg is taken as the reference for RfD because it is the most water-soluble Hg compound and its value was obtained as reported on the Risk Assessment Information System website [16].Since two scenarios were considered, the risk was assessed in terms of cumulative risk through the so-called HI Hazard Index, given by the sum of the HQs for the two scenarios.If one of the HQs or the HI tends to take values close to or even above unity, this means that the safe exposure threshold has been exceeded and, therefore, adverse effects related to Hg exposure may occur.To mitigate the uncertainties inherent in the risk assessment, the analysis was carried out using the traditional deterministic method [21] and a probabilistic approach [22].
Equation (1)-Average daily dose by the route of ingestion, USEPA 2001 where C GW (mg/L) is the concentration of Hg detected in the analyzed samples; EF (days/year) is the frequency of exposure; IR (L/days) is the ingestion rate; and ED (years) is the duration of exposure.The parameter "AT" represents the average exposure time.For non-carcinogenic toxic substances, the AT value is considered to coincide with the exposure duration (ED).The latter is expressed in years, while AT is expressed in days.Therefore, AT is obtained by multiplying ED by the number of days present in the year (365).In this way, the ADD ingestion that is obtained is expressed in (mg/kg-day).Finally, BW (kg) is the individual's body weight.
Equation (2)-Average daily dose by the route of dermal contact, USEPA 2004 where CGW (mg/L) is the concentration of Hg detected in the analyzed samples; EF (days/year) is the frequency of exposure; ET (hours/event) is the exposure time; ED (years) is the exposure time; SA (cm 2 ) is the exposed skin area of the individual; and kp (cm/hours) is the skin permeability constant.AT (days) is the average exposure time, which for noncarcinogenic toxic substances can be assumed to be equivalent to the exposure duration ED (years).Therefore, AT is obtained by converting ED from years to days, multiplying it by the number of days present in the year, i.e., 365.BW (kg) represents the body weight of the individual.The meaning of each parameter is summarized below:  • Body Weight (BW) signifies the average body weight of the receptor and thus varies accordingly; • Chronic Reference Dose (RfD): This represents the maximum dose of toxic contaminant considered acceptable.Essentially, it is the concentration value of the pollutant for which no adverse health effects have been documented in the literature.Specifically, for the purposes of this analysis, reference was made to RfD values for elemental mercury (metallic).In reality, the reference dose (RfD) is the result of a rough estimate and may have an uncertainty of up to an order of magnitude.The value of the reference dose (RfD) is decisive for the outcome of the analysis results.The reference dose (RfD) is usually determined, for different compounds and exposure modes according to the following equation: where: • NOAEL stands for "No Observed Adverse Effect Level", and among several possible values, the most conservative one is chosen.In the absence of the NOAEL, the LOAEL is used; • UF is the uncertainty factor; • MF is the modifying factor.This section may be divided by subheadings.It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn;

•
For the purpose of this analysis, the EPA-recommended values were used for all the parameters indicated [23].

Calculation of Hazard Quotient (HQ) and Hazard Index (HI)
The Hazard Quotient HQ is calculated as the ratio between the Average Daily Dose and the Reference Dose, as follows: Equation (3)-Hazard Quotient (HQ) Spatial data processing and map creation were conducted using the open-source software QGIS (version 3.16).Probabilistic calculations were executed using the R programming language, involving 10,000 iterations to enhance the reliability of the analysis.The parameters utilized in the analysis, along with their probabilistic distributions, are presented in Table 1.

Mercury (Hg) Accumulation in Groundwater
Table 2 shows a statistical summary of samples with mercury accumulation in groundwater present in selected provinces.Of the total number of samples studied, 54.8% were below the permissible limit according to the equipment used.(LoD of 0.0005 mg/L).For the purpose of the analysis, a concentration of half of the same instrument detectability threshold, or 0.00025 mg/L, was assumed at these points.On the other hand, 9.71 percent of the samples exhibited concentrations of 0.006 mg/L imposed by Ecuadorian regulations with reference to standards for safe drinking water.Of the data exceeding the Maximum Permissible Limit (MPL) value, 32.60% were found in Azuay province alone, and 7.69% in Loja province.In Chimborazo province, 100% of the tests studied exhibited Hg concentrations under the instrument's identification threshold (LoD = 0.0005 mg/L).This is followed by the provinces of Cotopaxi at 88.23%, Imbabura at 68.42%, Bolivar at 62.50%, Loja at 61.53%, Pichincha at 45.45%, Azuay at 36.95%, and Carchi, and Cañar at the same value of 28.57%.Figure 2

Spontaneous Perspective
Taking into account the philosophical doctrine of determinism, the impacts on human welfare are cataloged in Figures 3-6 by means of a punctual damage map.The numerical data of HI by sector (Table 3) are presented below: Loja > Azuay > Bolivar > Imbabura > Cotopaxi > Pichincha > Cañar > Chimborazo > Carchi.The risk analysis performed using the traditional deterministic approach showed a risk of 7.43% for the adult population in the residential scenario and 10.86% in the recreational scenario.Meanwhile, for the child population, a risk of 18.29% and 17.71% respectively emerged according to the above scenarios.The risk analysis makes it possible to estimate the risk associated with different receptors, according to their vulnerability.The latter is influenced, among other factors, by age.In risk analysis, unless specific data are provided, all individuals aged 0-6 years are considered 'children'.In a more detailed risk analysis, surveys are carried out to allow more accurate quantification of the type of targets present in the survey area.Therefore, the results of the analysis suggest a higher vulnerability in both scenarios for the child population, which is more exposed.Finally, thanks to the results obtained from the residential context, the main route of exposure was found to be ingestion.On the other hand, the following result demonstrates that the highest risk mode of exposure is dermal contact.3) are presented below: Loja > Azuay > Bolivar > Imbabura > Cotopaxi > Pichincha > Cañar > Chimborazo > Carchi.The risk analysis performed using the traditional deterministic approach showed a risk of 7.43% for the adult population in the residential scenario and 10.86% in the recreational scenario.Meanwhile, for the child population, a risk of 18.29% and 17.71% respectively emerged according to the above scenarios.The risk analysis makes it possible to estimate the risk associated with different receptors, according to their vulnerability.The latter is influenced, among other factors, by age.In risk analysis, unless specific data are provided, all individuals aged 0-6 years are considered 'children'.In a more detailed risk analysis, surveys are carried out to allow more accurate quantification of the type of targets present in the survey area.Therefore, the results of the analysis suggest a higher vulnerability in both scenarios for the child population, which is more exposed.Finally, thanks to the results obtained from the residential context, the main route of exposure was found to be ingestion.On the other hand, the following result demonstrates that the highest risk mode of exposure is dermal contact.To evaluate the determination of mercury (Hg) concentration in groundwater in different regions, a detailed analysis was carried out considering different receptors and scenarios.This study was carried out in order to identify mercury levels that could represent a risk to public health and the environment.Table 4 below shows the determinism values (p95) of Hg in groundwater, broken down by the different receptors and scenarios evaluated.This information is crucial for understanding variations in Hg concentration and for developing appropriate mitigation strategies.Figure 3 shows the danger index for adults in the residential scenario, where the highest index was found in the provinces of Azuay and Loja.
The figure provides a detailed overview of the Hazard quotient (HQ) for children in residential areas.It highlights that the provinces of Azuay and Loja show the most significant levels of risk.This crucial analysis provides a clear understanding of the areas where children may face the greatest dangers in their residential environment.The evidence underscores the need for targeted interventions in these regions to ensure the safety and well-being of the very young.
The analysis of the Hazard Index (HQ) for adults in recreational settings reveals significant concerns, particularly in the provinces of Azuay and Loja, where the maximum permissible limit is exceeded.
Analysis of the Danger Index (Hq) in recreational settings for children reveals areas of concern, particularly in the provinces of Carchi, Pichincha, Cotopaxi, Bolivar, Azuay, and Loja.These results highlight the urgency of addressing the specific risks faced by children in these recreational settings.It is crucial to implement preventive measures and effective safety policies to ensure their protection during their recreational activities.

Azuay
HQ_dermal contact 8.62 × 10     The analysis of the Hazard Index (HQ) for adults in recreational settings reveals significant concerns, particularly in the provinces of Azuay and Loja, where the maximum permissible limit is exceeded.The analysis of the Hazard Index (HQ) for adults in recreational settings reveals significant concerns, particularly in the provinces of Azuay and Loja, where the maximum permissible limit is exceeded.Analysis of the Danger Index (Hq) in recreational settings for children reveals areas of concern, particularly in the provinces of Carchi, Pichincha, Cotopaxi, Bolivar, Azuay, and Loja.These results highlight the urgency of addressing the specific risks faced by children in these recreational settings.It is crucial to implement preventive measures and effective safety policies to ensure their protection during their recreational activities.

Probability Perspective
Process hazard studies use the perspective of traditional determinism under a regime of unwanted events; in analyses of this type, the larger the sample of data available, the greater the reliability of the results.However, it may also be useful to conduct a risk analysis using a different approach, in particular a probability of comparing corresponding

Probability Perspective
Process hazard studies use the perspective of traditional determinism under a regime of unwanted events; in analyses of this type, the larger the sample of data available, the greater the reliability of the results.However, it may also be useful to conduct a risk analysis using a different approach, in particular a probability of comparing corresponding responses.In particular, with the help of the R programming language [28].This technique makes it possible to assess the risk associated with a given process or phenomenon through the construction of a series of simulated scenarios.The construction of the simulated scenarios consists of generating sequences of random numbers using a special algorithm, according to a given probabilistic distribution law, to be assigned to the variables of the phenomenon to be analyzed.The application of the method first requires the preliminary definition of a mathematical model suitable for representing the phenomenon under analysis.Then, all the variables, their probability distribution, and the relationships between them must be identified and defined.Based on the constructed model, using a suitable software tool, sequences of random numbers can be generated and assigned to the variables involved: That is, the more simulations, the more precise the final risk assessment will be.
By applying the Monte Carlo simulation, a number of manifestations of the phenomenon (process) analyzed can be determined.In the present study, the random determinations of the phenomenon under investigation were obtained by synthetically generating random numbers by writing a special algorithm in R. In particular, 10,000 iterations were performed for each simulated quantity.Accurate modeling of a process with the Monte Carlo method requires the generation of high-quality pseudo-random numbers.The main limitation of generating pseudo-random numbers with a computer is that the values begin to repeat themselves after a certain generation sequence.However, due to the great interest in the development of the Monte Carlo technique, many scientists have developed algorithms to solve this problem.In particular, R implements the Mersenne-Twister algorithm by default.The Mersenne-Twister algorithm, developed by Matsumoto and Nishimura, guarantees a sequence of 2 19937 −1 elements before repetition and this property has been proven by its developers and more, this algorithm is faster than most other algorithms and has passed many statistical tests.
Monte Carlo simulation offers a full range of possible outcomes by implementing these stochastic simulation techniques [24].This type of method is considered valid by many scholars for the analysis of problems requiring risk evaluation or in decision-making problems involving many variables [29].Once the random sample has been determined, it is possible to determine the quantities of interest and obtain very explanatory graphical representations.The probability distributions and intervals used for each parameter are shown in Table 5.The R software (version 4.4.1)has special libraries that allow sequences to be generated in accordance with each probability distribution.
Figure 7 presents four box-and-whisker plots, corresponding to each exposure route and receptor category.Each of these plots illustrates the distribution of the Hazard Quotient (HQ) for the exposure pathway and receptor type analyzed.The lower and upper boundaries of the 'box' symbolize the 25th and 75th percentiles, respectively.The 'lines' or 'whiskers' that extend from the box display the variability in the expected data.The points located outside these whiskers constitute atypical or extreme observations.

Discussion
The study's findings paint a troubling picture regarding Hg pollution in groundwater in the Andean region of Ecuador.A discrepancy was found between the results obtained by applying the two approaches; based on the deterministic approach, certain locations present risks associated with exposure.This inconsistency in the findings derived from the two methodologies has been observed in other research as well, such as a study by Mestanza-Ramón et al. in 2023 that analyzed mercury levels in superficial water bodies within the Andean area and encountered similar discrepancies in the results obtained through the different techniques applied [22].This divergence lies in the inherent methodology of each of the two approaches considered, since the deterministic approach involves a specific and concrete assessment, whereas the probabilistic perspective consists of a stochastic-type analysis that considers the variability of all the parameters involved in the phenomenon.However, in qualitative terms, both approaches suggest greater exposure to the child population.
There was a non-uniform distribution of Hg concentrations across the samplings, and with a wide range of variability: in some places concentrations were at levels below the minimum detectable limit of the measuring device, in others they exceeded the levels surpassed the maximum permissible limit by Ecuadorian standards for water intended for human consumption.This spatial pattern may be influenced by depends on multiple variables, only some of which are described below.First, the sampling points were chosen according to "reasoned location", that is, based on the location of the mining areas and the information collected in the field.Nevertheless, the mining activities that have the most significant impact on Hg contamination are illegal activities, so they are not easy to track.In addition, one must consider the particular inclination of mercury to disperse very rapidly: when in its liquid form, it fractionates into tiny spheres that exhibit a significant ability to move, and additionally, contingent upon environmental conditions, it has the potential to evaporate and subsequently condense again with swiftness.These particular characteristics mean that mercury can disperse very rapidly by diffusing from its point of origin.
In any case, it appears that the most affected provinces are those located in the south of the region: Loja and Azuay.The network of sampling points is not uniformly distributed over the territory but was determined according to a "reasoned location".In other words, based on the information collected, sampling sites were selected in the areas closest to mining activities.At mining locations, especially ASGM, other potentially toxic substances are used in addition to Hg, such as arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), nickel (Ni), lead (Pb), tin (Sn) and aluminum (Al).It should be considered that multiple toxic substances in combination can dramatically exacerbate potentially hazardous conditions, and additionally, as mining activities continue to increase, Hg (and also of any other toxic substances) concentrations could continue to rise over time, resulting in a higher risk status for populations in these areas.
Figure 9 shows the choice of sampling points and the locations where pollution from mining activities has been officially recognized in the Andean region.It is possible to see that most of the sampling points fall in the vicinity of the points where pollution was detected.The prospect of contamination affecting groundwater makes the overall picture even more worrisome.Indeed, from groundwater, the chances of the contaminant reaching human receptors or entering the food chain are very high.Actually, underground aquifers constitute the primordial water resource for supplying the population.Previous research maintains that a primary means of mercury exposure is the unintentional consumption of contaminated water.In addition, from groundwater, contamination can reach surface water bodies (especially during lean periods), contaminating the fish population.In particular, as attested by the "Exposure to Mercury" document prepared by WHO (2021), the foremost source of contact with this organometallic form of mercury is through the dietary intake of aquatic animals such as finfish and shellfish, especially for populations where the diet is based on the consumption of predatory fish [27].A more specific analysis of the dietary habits of the populations inhabiting the analyzed area could allow a more accurate assessment of the impact of this source of contamination.Contamination can also accumulate on the ground, both shallow and deep, posing a danger, especially to children in recreational scenarios.
The essential purpose of a risk analysis consists of offering an estimation of the danger posed to the health of the population and ecosystems by the presence of one or more contaminants in different environmental components.Accurate risk analysis requires very thorough site-specific investigations, which require a significant deployment of financial funds, specialized knowledge, and support infrastructure.Consequently, it is invariably advisable for such analyses to be preceded by general analyses, which, considering vast territories (such as the entire Andean region), allow one to identify which areas are most affected by the presence of the contaminant(s) in question, and therefore potentially contaminated.For these reasons, first-level risk analyses are conducted on the basis of very conservative assumptions.Therefore, total mercury was used as the analysis criterion because potentially all mercury in the environment could potentially be transformed into its most toxic form.In fact, mercury introduced into the environment by anthropogenic activities often occurs in an inorganic form, however, if favorable conditions are established, natural processes may occur that lead to its conversion into methylmercury [CH 3 Hg] + , which has a high persistence capacity that allows it to bioaccumulate in the environment even for centuries [10].By applying this precautionary approach, results can be obtained to the benefit of safety and the most affected areas can be easily identified, enabling decision-makers to plan and conduct more detailed analyses.specific analysis of the dietary habits of the populations inhabiting the analyzed area could allow a more accurate assessment of the impact of this source of contamination.Contamination can also accumulate on the ground, both shallow and deep, posing a danger, especially to children in recreational scenarios.The essential purpose of a risk analysis consists of offering an estimation of the danger posed to the health of the population and ecosystems by the presence of one or more contaminants in different environmental components.Accurate risk analysis requires very thorough site-specific investigations, which require a significant deployment of financial funds, specialized knowledge, and support infrastructure.Consequently, it is invariably advisable for such analyses to be preceded by general analyses, which, considering vast territories (such as the entire Andean region), allow one to identify which areas are most affected by the presence of the contaminant(s) in question, and therefore potentially contaminated.For these reasons, first-level risk analyses are conducted on the basis of very conservative assumptions.Therefore, total mercury was used as the analysis criterion because potentially all mercury in the environment could potentially be transformed into its most toxic form.In fact, mercury introduced into the environment by anthropogenic activities often occurs in an inorganic form, however, if favorable conditions are established, natural processes may occur that lead to its conversion into methylmercury [CH3Hg] + , which has a high persistence capacity that allows it to bioaccumulate in the environment even for centuries [10].By applying this precautionary approach, results can be obtained to the benefit of safety and the most affected areas can be easily identified, enabling decision-makers to plan and conduct more detailed analyses.A research based on the thorough analysis of 1573 scientific publications found consistent results regarding the relationship between exposure to heavy metallic elements (HM) and the repercussions on neurocognitive abilities in adult individuals [28].The conclusions indicated a correlation between contact with HM and the deterioration of neurocognitive functions in the adult population.Specifically, certain observations have detected a decrease in the aptitude to retain information in the short term as exposure to mercury in its elemental form increases [29].Likewise, another two-cycle cross-sectional study, carried out in Korea with a total sample of 14,682 participants, discovered that Hg levels (together with two other metals) showed a strong association with alcohol-related liver disease [30].

Conclusions
This research focused on analyzing Hg levels within the Andean area of Ecuador.The provinces located in the southern region were found to be the most impacted.However, in certain zones, the identified levels of concentration were lower than the maximum permissible limit set by Ecuador's standards for water meant for human use.The results from the two approaches applied, deterministic and probabilistic, show different results: the deterministic approach provided higher levels of risk.Nevertheless, both showed an increased vulnerability of the child population compared to adults.Therefore, it might be appropriate to carry out a more detailed analysis, especially in the areas most at risk.In fact, in order to conduct a more detailed risk analysis, it is necessary to dispose of a set of data to contextualize the risk to the territories of reference.
The contaminant in the pathway from the source to the point of exposure may be partially attenuated, depending on the transport processes involved.To take these aspects into account, it is necessary to have a complete mapping of the structure and physical characteristics of the individual site analyzed.However, the aim of this study is to realize a general screening of concentrations; generally, in studies of this type, the analysis is carried out with very conservative assumptions, so transport mechanisms are neglected, and the detected concentrations are assumed to be valid.Where decision-makers deem it appropriate, they may use the results summarized in this study as a guide and proceed to gather more information to conduct a more detailed analysis.

Figure 1 .
Figure 1.Study area and location of the sampled sites.

Figure 1 .
Figure 1.Study area and location of the sampled sites.
Hg concentration in groundwater mg L EF = exposure f requency days year IR = ingestion rate o f water L day ET = exposure time hours event ED = li f etime exposure duration (years) SA = skin sur f ace area exposed cm 2 k p = skin permeability constant cm hour AT = averaging time (days) BW = body weight (kg) 2.3.1.Insight into the Significance of Parameters below shows the Hg concentrations at the locations where the samples were taken.

Figure 2 .
Figure 2. Hg concentrations at the sampled sites.

Figure 2 .
Figure 2. Hg concentrations at the sampled sites.

3. 2 .
Harm to People's Wellbeing 3.2.1.Spontaneous Perspective Taking into account the philosophical doctrine of determinism, the impacts on human welfare are cataloged in Figures 3-6 by means of a punctual damage map.The numerical data of HI by sector (Table

Figure 3
Figure3shows the danger index for adults in the residential scenario, where the highest index was found in the provinces of Azuay and Loja.

Figure 6 .
Figure 6.Hazard index (HI) for children in recreational scenarios.

Figure 6 .
Figure 6.Hazard index (HI) for children in recreational scenarios.

Figure 7 .
Figure 7. Boxplots of the hazard index (HI) for both receptors.The prefixes C and A refer to children and adult receptors respectively, while the suffixes -res and -ricr refer to the residential and recreational scenarios.

Figure 7 .
Figure 7. Boxplots of the hazard index (HI) for both receptors.The prefixes C and A refer to children and adult receptors respectively, while the suffixes -res and -ricr refer to the residential and recreational scenarios.

Figure 8
Figure 8 similarly shows the development of the probability density function for the different scenarios and receptors.Specifically, the x-axis shows the values of HI on a logarithmic scale, and the corresponding probability densities are shown on the y-axis.

Figure 7 .
Figure 7. Boxplots of the hazard index (HI) for both receptors.The prefixes C and A refer to children and adult receptors respectively, while the suffixes -res and -ricr refer to the residential and recreational scenarios.

Figure 8 .
Figure 8. Histograms of probability densities for different receptors and scenarios.Figure 8. Histograms of probability densities for different receptors and scenarios.

Figure 8 .
Figure 8. Histograms of probability densities for different receptors and scenarios.Figure 8. Histograms of probability densities for different receptors and scenarios.

Figure 9 .
Figure 9. Location of sampling points with respect to areas contaminated by mining activities (general pollution).

Figure 9 .
Figure 9. Location of sampling points with respect to areas contaminated by mining activities (general pollution).

Table 1 .
Parameters used in the risk assessment.
*: Data with minimum results according to the measuring device (LoD = 0.0005 mg/L) are presented below; for analysis purposes, a so-called present value (LoD/2 = 0.00025 mg/L) was assumed.Percentiles are referred to as minimum data over a certain percentage present.That is to say, p50 translates as a maximum value below which 50% falls.

Table 5 .
Parameters used in probabilistic risk assessment.