Study and Evaluation of Surface Water Resources Affected by Ancient and Illegal Mining in the Upper Part of the Milluni Micro-Basin, Bolivia

One of the factors with the most significant impact due to pollution caused by mining activities is water. Surface water resources are limited, and in many cases, they are the primary water source for consumption; thus, its quality is a fundamental subject of study. The traditional way of studying water quality is to measure parameters and compare the results with tables of values established in the regulations. The quality of drinking water is characterized and controlled. However, this traditional comparison method may be subjective and insufficient to classify water quality as safe. Therefore, it is necessary to develop and complement the conventional monitoring, linking physicochemical analysis with remote perception. The different variables analyzed can be integrated and result in much more complete water quality indices, achieving the evolution of the monitoring programs of the water quality of the surface bodies. This research combines remote sensing with traditional monitoring to obtain a more robust and affordable indicator, resulting in a study of the central water bodies of Milluni. The main finding—the alarming contamination in the Milluni Chico and Milluni Grande lagoons—provides more accurate information to decision-makers to improve water management in Bolivia.


Introduction
The sustainable and affordable supply of clean, safe, and adequate quality water is one of the world's most challenging problems [1]. As is known, the relationship between mining and sustainability is complex and controversial [2]. That is why many developing countries whose main economic activity is mining are highly vulnerable to the negative impacts on different environmental factors, especially water. According to the World Health Organization, diseases associated with the lack of drinking water, adequate sanitation, and hygiene continue to be one of the leading causes of mortality for millions of developing countries [3]. The current water problem is reflected as a priority in Goal 6, "Clean Water and Sanitation," within the Sustainable Development Goals (SDG), which frame the 2030 agenda [4].
Water resources in developing regions of the world have high social and ecological value, and many vulnerable human communities depend on these resources for their • Pata Khota: A natural lagoon of irregular shape that receives water from the glacier of the Huayna Potosí. It is located at 4665 m.a.s.l. • Jankho Khota: Located at 4560 m.a.s.l., an irregular-shaped natural lagoon that receives water from the Pata Khota lagoon. • Milluni Chico: An artificial lagoon at an altitude of 4550 m.a.s.l., having an irregular shape and receiving water from natural springs, Jankho Khota lagoon, and mines. This artificial lagoon aims to capture the acid drainage of the mines to prevent it from entering the Milluni Grande lagoon. • Milluni Grande: Located at 4530 m.a.s.l., the Milluni Grande lagoon receives effluents from natural springs, mines' drainage, and the Milluni Chico lagoon, specifically located 2.7 km from this lagoon; it also receives water from the Jankho Khota lagoon through a bypass system. There is a water storage dam for public supply in this lagoon, which has a high probability of contamination by heavy metals due to its tributaries.
The upper part of the Milluni micro-basin is relevant because its lagoons are a source of water for the public supply. This area has a high probability of contamination by heavy metals, which was exposed in the study in 2008 [19] and was ratified in 2019 [26]. This last study carried out a review of historical data on water quality in the area, where the presence of heavy metals such as arsenic, cadmium, iron, zinc, and manganese was identified, which were found with concentrations higher than those values established in Bolivian Standard 512 for drinking water [27]. For this reason, the Milluni area must be closely monitored to avoid risk to public health. For this purpose, a monitoring program for surface water bodies was designed, which considers the Milluni area's specific characteristics [22].
The Water Quality Monitoring Program (WQMP) for Milluni was based on two international standards: the UNEP/WHO monitoring program design guide [28] and Bolivian Norm NB, which adopts ISO 5667-1 [29], the latter specific to bodies' surface water. The objective of the WQMP is to determine the quality of the water of Milluni to become a control program that generates representative data and contributes to the water management of the place. There were three designated monitoring points; their locations can be seen in Figure 1, and the image of each of them can be seen in Figure 2.
Resources 2022, 11, x FOR PEER REVIEW 3 of 16 • Milluni Chico: An artificial lagoon at an altitude of 4550 m.a.s.l., having an irregular shape and receiving water from natural springs, Jankho Khota lagoon, and mines. This artificial lagoon aims to capture the acid drainage of the mines to prevent it from entering the Milluni Grande lagoon.

•
Milluni Grande: Located at 4530 m.a.s.l., the Milluni Grande lagoon receives effluents from natural springs, mines' drainage, and the Milluni Chico lagoon, specifically located 2.7 km from this lagoon; it also receives water from the Jankho Khota lagoon through a bypass system. There is a water storage dam for public supply in this lagoon, which has a high probability of contamination by heavy metals due to its tributaries.
The upper part of the Milluni micro-basin is relevant because its lagoons are a source of water for the public supply. This area has a high probability of contamination by heavy metals, which was exposed in the study in 2008 [19] and was ratified in 2019 [26]. This last study carried out a review of historical data on water quality in the area, where the presence of heavy metals such as arsenic, cadmium, iron, zinc, and manganese was identified, which were found with concentrations higher than those values established in Bolivian Standard 512 for drinking water [27]. For this reason, the Milluni area must be closely monitored to avoid risk to public health. For this purpose, a monitoring program for surface water bodies was designed, which considers the Milluni area's specific characteristics [22].
The Water Quality Monitoring Program (WQMP) for Milluni was based on two international standards: the UNEP/WHO monitoring program design guide [28] and Bolivian Norm NB, which adopts ISO 5667-1 [29], the latter specific to bodies' surface water. The objective of the WQMP is to determine the quality of the water of Milluni to become a control program that generates representative data and contributes to the water management of the place. There were three designated monitoring points; their locations can be seen in Figure 1, and the image of each of them can be seen in Figure 2.   (a) (b) (c) Point 1 (E: 592549; N: 8197635) has the purpose of determining the reference conditions in the watercourse system; for this reason, it is located in the headwaters of the lagoon system in the Pata Khota lagoon. Point 2 (E: 589887; N: 8193022) seeks to identify any sign of deterioration in the quality of the water in the water system due to anthropogenic activities and is located in the Milluni Chico lagoon. Finally, Point 3 (E: 589147; N: 8190424) aims to expose whether the body of water in the course system meets the desired quality standards when located in a storage dam.

Methodology
The study was carried out in 3 stages, as shown in Figure 3, and developed in the following sections.   Point 1 (E: 592549; N: 8197635) has the purpose of determining the reference conditions in the watercourse system; for this reason, it is located in the headwaters of the lagoon system in the Pata Khota lagoon. Point 2 (E: 589887; N: 8193022) seeks to identify any sign of deterioration in the quality of the water in the water system due to anthropogenic activities and is located in the Milluni Chico lagoon. Finally, Point 3 (E: 589147; N: 8190424) aims to expose whether the body of water in the course system meets the desired quality standards when located in a storage dam.

Methodology
The study was carried out in 3 stages, as shown in Figure 3, and developed in the following sections. Point 1 (E: 592549; N: 8197635) has the purpose of determining the reference conditions in the watercourse system; for this reason, it is located in the headwaters of the lagoon system in the Pata Khota lagoon. Point 2 (E: 589887; N: 8193022) seeks to identify any sign of deterioration in the quality of the water in the water system due to anthropogenic activities and is located in the Milluni Chico lagoon. Finally, Point 3 (E: 589147; N: 8190424) aims to expose whether the body of water in the course system meets the desired quality standards when located in a storage dam.

Methodology
The study was carried out in 3 stages, as shown in Figure 3, and developed in the following sections.

Climatological Analysis
The hydrological analysis is essential for monitoring the quality of water resources [30,31]. Through the climatological analysis of the area of interest, it will be possible to know the distribution of dry and wet seasons by month in a hydrological year (July-June). The above is used to interpret the incidence of hydrological behavior on contamination levels in the central bodies of water in the Milluni micro-basin. For the development of this analysis, the following procedure was carried out: 1. Compilation of meteorological information from the Alto Achachicala station of SENHAMI for the climatological analysis [32]: precipitation (mm) and evapotranspiration (mm) with a return rate of 18 years (1999-2018), in the framework of a hydrological year;

Climatological Analysis
The hydrological analysis is essential for monitoring the quality of water resources [30,31]. Through the climatological analysis of the area of interest, it will be possible to know the distribution of dry and wet seasons by month in a hydrological year (July-June). The above is used to interpret the incidence of hydrological behavior on contamination levels in the central bodies of water in the Milluni micro-basin. For the development of this analysis, the following procedure was carried out: 1.
Compilation of meteorological information from the Alto Achachicala station of SEN-HAMI for the climatological analysis [32]: precipitation (mm) and evapotranspiration (mm) with a return rate of 18 years (1999-2018), in the framework of a hydrological year;

2.
Statistical processing of the information collected to obtain total data and averages per month of the precipitation and evapotranspiration parameters; 3.
Generation of a graph with monthly data of total precipitation, mean, and total evapotranspiration parameters; 4.
Analysis of the results obtained to identify dry and wet seasons.

Study of Water Quality of the Milluni Lagoons
This section presents the two parts developed for the Milluni water quality study. The first part is traditional monitoring, and the second is remote sensing of contaminants.

Traditional Field Monitoring
Fieldwork to monitor Milluni was completed in 2019 following the designed WQMP. The monitoring frequency was determined according to the temporality of the study area. The samplings were carried out three times in the dry season and three times in the rainy season with intervals of 1 month between samples. Fundamental in situ control parameters (pH, dissolved oxygen, conductivity, turbidity, and temperature) were measured with an HQ40D-Hach Handheld Multimeter and Hach 2100Q Handheld Turbidimeter. Ex situ parameters (heavy metals) were analyzed with an Agilent Model 7700x Inductively Coupled Plasma Mass Spectrometer (ICP-MS). The sampling protocol adopted was ISO 5667-4:2016 Water quality-Sampling-Part 4: Guidance on sampling of lakes, natural and artificial [33].

Remote Sensing
This part applied sequential processes developed based on remote sensing through the Google Earth Engine platform. Thanks to its adequate spatial resolution, the satellite images selected for this purpose were captured by the Sentinel-2 satellite. The temporality that was handled was monthly during 2019; it was fitted to the availability of images captured by the satellite and the monitoring dates in the field.
To study the water quality in the central water bodies of the Milluni micro-basin, the manipulation of satellite bands was carried out using scripts in Java language within the Google Earth Engine platform. Two leading indicators were obtained: the Normalized Water Index (NDWI) and the Normalized Chlorophyll-a Index (NDCl). It should be noted that the NDCI is an indicator of the presence of eutrophication. The equations to determine the values for the indices are presented below.

•
NDWI. To obtain the pixels that represent water, the work was carried out with two spectral bands, green and shortwave infrared (SWIR), following Equation (1).
To obtain the pixels representing water, the work was completed with two spectral bands, Red edge one and Red, following Equation (2).
It was possible to identify the representative pixels of areas with water coverage and their quality level. Figure 4 details the sequence of processes that have been followed in this phase of the investigation. The satellite images obtained from the previous process are presented in Section 3.2.2. The lagoons will show three colors related to the concentration levels determined for the NDCI. It is essential to analyze the three levels of the indicator to assess the risk of eutrophication in each lagoon. Low concentration, represented by yellow, means little chance of eutrophication. The medium concentration implies a moderate risk of eutrophication with an orange coloration. Finally, a high concentration, presented in red, suggests a high risk of eutrophication.

Analysis and Interpretation of Results
The last phase of this research interprets the correlation of the results obtained from remote sensing and the values obtained through implementing the first year of the Water Quality Monitoring Program (WQMP), where in situ and ex situ parameters were measured.

Results
In this section, the results will be presented for the different stages of the study.

Climatological Analysis
Information was collected from SENHAMI meteorological stations on precipitation and evapotranspiration from 1999 to 2018 [32]. Table 1 below shows the average data for all months of the year in terms of total rainfall and total and mean evapotranspiration.  The satellite images obtained from the previous process are presented in Section 3.2.2. The lagoons will show three colors related to the concentration levels determined for the NDCI. It is essential to analyze the three levels of the indicator to assess the risk of eutrophication in each lagoon. Low concentration, represented by yellow, means little chance of eutrophication. The medium concentration implies a moderate risk of eutrophication with an orange coloration. Finally, a high concentration, presented in red, suggests a high risk of eutrophication.

Analysis and Interpretation of Results
The last phase of this research interprets the correlation of the results obtained from remote sensing and the values obtained through implementing the first year of the Water Quality Monitoring Program (WQMP), where in situ and ex situ parameters were measured.

Results
In this section, the results will be presented for the different stages of the study.

Climatological Analysis
Information was collected from SENHAMI meteorological stations on precipitation and evapotranspiration from 1999 to 2018 [32]. Table 1 below shows the average data for all months of the year in terms of total rainfall and total and mean evapotranspiration. Note: Total precipitation and total evapotranspiration are referred to as the full water considered through these climatic processes in a given period, such as a month or year. Mean evapotranspiration refers to the average data calculated from all the data collected in a certain period; for example, by collecting daily data, we can obtain an average in a particular month. Source: Data from SENAMHI, 2021.
With the data in Table 1, Figure 5 was made, and this area's dry and humid periods can be appreciated. On the one hand, the dry period starts in May and last until October. On the other hand, the wet period begins in November and ends at the end of April. Note: Total precipitation and total evapotranspiration are referred to as the full water considered through these climatic processes in a given period, such as a month or year. Mean evapotranspiration refers to the average data calculated from all the data collected in a certain period; for example, by collecting daily data, we can obtain an average in a particular month. Source: Data from SENAMHI, 2021.
With the data in Table 1, Figure 5 was made, and this area's dry and humid periods can be appreciated. On the one hand, the dry period starts in May and last until October. On the other hand, the wet period begins in November and ends at the end of April. From the climatological analysis of the Milluni area, it is observed that the maximum precipitation peak occurs during the month of January. This can cause the lagoons to increase volume and the pollutant dilution phenomenon to occur naturally. In addition, it should be considered that the driest months are June and July. In the dry period, a concentration of pollutants could occur due to reducing the volume of water in the lagoons. The months of January, June, and July should be carefully examined later when interpreting the results.

Study of the Water Quality of Milluni Lagoons
This section presents the results of traditional monitoring where physicochemical control parameters were measured and the results were obtained by remote sensing of contaminants. From the climatological analysis of the Milluni area, it is observed that the maximum precipitation peak occurs during the month of January. This can cause the lagoons to increase volume and the pollutant dilution phenomenon to occur naturally. In addition, it should be considered that the driest months are June and July. In the dry period, a concentration of pollutants could occur due to reducing the volume of water in the lagoons. The months of January, June, and July should be carefully examined later when interpreting the results.

Study of the Water Quality of Milluni Lagoons
This section presents the results of traditional monitoring where physicochemical control parameters were measured and the results were obtained by remote sensing of contaminants. In the three points defined by the Water Quality Monitoring Program for the upper part of Milluni, 4 in situ control parameters and 33 ex situ control parameters were monitored. The results of the in situ and ex situ parameters are presented below.

Results of In-Situ Parameters
The pH and conductivity are related to the presence of metal ions and dissolved oxygen, and turbidity is directly associated with NDCI, which is an indicator of eutrophication. As for the four parameters measured in situ, point 1 is the only one with values below what is stipulated in the Bolivian standard 512 [27], indicating that it is free of contamination. However, points 2 and 3 are negatively impacted by anthropogenic activities in the area, presenting acidic pH values, conductivities above 700 uS/cm, and high turbidity in January for point 2 and October for point 3. Figure 6 shows the measured values for the in situ parameters at the three monitored points. In the three points defined by the Water Quality Monitoring Program for the upper part of Milluni, 4 in situ control parameters and 33 ex situ control parameters were monitored. The results of the in situ and ex situ parameters are presented below.

Results of In-Situ Parameters
The pH and conductivity are related to the presence of metal ions and dissolved oxygen, and turbidity is directly associated with NDCI, which is an indicator of eutrophication. As for the four parameters measured in situ, point 1 is the only one with values below what is stipulated in the Bolivian standard 512 [27], indicating that it is free of contamination. However, points 2 and 3 are negatively impacted by anthropogenic activities in the area, presenting acidic pH values, conductivities above 700 uS/cm, and high turbidity in January for point 2 and October for point 3. Figure 6 shows the measured values for the in situ parameters at the three monitored points. These four parameters will be correlated with the indicator of the presence of eutrophication in Section 3.3.

Results of Ex Situ Parameters
At the same time, 33 heavy metals were monitored. Seven metal ions were found with concentrations outside the current environmental regulations: arsenic, cadmium, iron, zinc, manganese, nickel, and aluminum. Some of these metals were also identified These four parameters will be correlated with the indicator of the presence of eutrophication in Section 3.3.

Results of Ex Situ Parameters
At the same time, 33 heavy metals were monitored. Seven metal ions were found with concentrations outside the current environmental regulations: arsenic, cadmium, iron, zinc, manganese, nickel, and aluminum. Some of these metals were also identified in previous studies [19,22]. In this study, only three of these will be analyzed due to their toxicity (arsenic), their high concentration (zinc), and their persistence in the area (manganese). Figure 7 shows the attention of the three metals analyzed. in previous studies [19,22]. In this study, only three of these will be analyzed due to their toxicity (arsenic), their high concentration (zinc), and their persistence in the area (manganese). Figure 7 shows the attention of the three metals analyzed. From the graphs presented above, we can see that point 1 continues to be the point without heavy metal contamination. However, points 2 and 3 are contaminated by the three metals studied. It should be noticed that point 2 is the most affected. The three metals were previously exposed. These parameters will also be correlated with the indicator of the presence of eutrophication in Section 3.3.

Results of Remote Sensing
The satellite images selected for the study include the months of February, April, June, July, August, September 2019, and January 2020. The months of January, March, May, October, November, and December 2019 could not be considered due to the low visibility caused by the dense clouds of those months. Figure 8 presents the satellite views of Milluni. From the graphs presented above, we can see that point 1 continues to be the point without heavy metal contamination. However, points 2 and 3 are contaminated by the three metals studied. It should be noticed that point 2 is the most affected. The three metals were previously exposed. These parameters will also be correlated with the indicator of the presence of eutrophication in Section 3.3.

Results of Remote Sensing
The satellite images selected for the study include the months of February, April, June, July, August, September 2019, and January 2020. The months of January, March, May, October, November, and December 2019 could not be considered due to the low visibility caused by the dense clouds of those months. Figure 8 presents the satellite views of Milluni. From Figure 8, it can be seen that there are no high concentrations of NDCI in any of the months studied, which shows that the risk of extreme eutrophication in the surface water bodies of Milluni is low. However, the medium concentrations of the NDCI are prevalent in some months, such as February, June, and August 2019, and January 2020. From Figure 8, it can be seen that there are no high concentrations of NDCI in any of the months studied, which shows that the risk of extreme eutrophication in the surface water bodies of Milluni is low. However, the medium concentrations of the NDCI are prevalent in some months, such as February, June, and August 2019, and January 2020. February 2019 and January 2020 are months in the wet season according to the area's climatological analysis, which would indicate that an increase or decrease in precipitation would not affect the medium concentration prevalent in the Milluni lagoons. This would imply that time does not fully influence the presence of organic pollutants in the Milluni lagoon system.

Analysis and Interpretation of the Results
The correlation of the results of the parameters monitored in situ, ex situ, and the NDCI for each lagoon belonging to the upper part of the Milluni micro-basin is shown below.

Pata Khota Lagoon
The Pata Khota Lagoon is the first lagoon of the water system, with waters coming from the natural glacier called Huayna Potosí. Monitoring point 1 was located in this lagoon. Point 1 is the least polluted according to the contrast of the in situ and ex situ parameters with the Bolivian Regulation 512. The dynamics of the Pata Khota lagoon are based on the climatological analysis where wet and dry times are denoted. The organic load based on what is registered by the NDCI is in the function of the season. Figure 9 shows concentrations of NDCI in the lagoon.
Resources 2022, 11, x FOR PEER REVIEW climatological analysis, which would indicate that an increase or decrease in precip would not affect the medium concentration prevalent in the Milluni lagoons. This imply that time does not fully influence the presence of organic pollutants in the M lagoon system.

Analysis and Interpretation of the Results
The correlation of the results of the parameters monitored in situ, ex situ, a NDCI for each lagoon belonging to the upper part of the Milluni micro-basin is below.

Pata Khota Lagoon
The Pata Khota Lagoon is the first lagoon of  When comparing the predominant chlorophyll-a concentration in the Pata Kh goon, which is the low one, with the turbidity and dissolved oxygen values, they s that there would not be a pollutant increasing the organic load of the lagoon. On t trary, if the dissolved oxygen values decrease and the turbidity values rise, a eutro tion process would occur in the Pata Khota lagoon, and the predominant concentra chlorophyll-a would be medium or high range. In addition, the monitored point neutral pH and a low conductivity throughout the year without being affected by or wet season. This would indicate the absence of metal ions, which is confirmed heavy metals' low concentrations.
In summary, the results of the monitoring and what can be seen in the satellite indicate that the Pata Khota lagoon would not be negatively affected by contam from human activities. The water in this body of surface water is of sufficient qual can be used for human consumption. When comparing the predominant chlorophyll-a concentration in the Pata Khota lagoon, which is the low one, with the turbidity and dissolved oxygen values, they showed that there would not be a pollutant increasing the organic load of the lagoon. On the contrary, if the dissolved oxygen values decrease and the turbidity values rise, a eutrophication process would occur in the Pata Khota lagoon, and the predominant concentration of chlorophyll-a would be medium or high range. In addition, the monitored point 1 has a neutral pH and a low conductivity throughout the year without being affected by the dry or wet season. This would indicate the absence of metal ions, which is confirmed by the heavy metals' low concentrations.
In summary, the results of the monitoring and what can be seen in the satellite images indicate that the Pata Khota lagoon would not be negatively affected by contamination from human activities. The water in this body of surface water is of sufficient quality and can be used for human consumption.

Jhanko Khota Lagoon
This lagoon was not monitored because, visually, it is very similar to the first Pata Khota lagoon. However, satellite images show that its water quality is different from the first lagoon. In the months of January and February, belonging to the wet season according to the climatological analysis carried out, the most notable differences are presented. This exposes the existence of a polluting agent that would be altering the water quality of this lagoon, the source of which must be studied. Figure 10 reveals the concentrations of NDCI in the second lagoon.
Resources 2022, 11, x FOR PEER REVIEW 1 This exposes the existence of a polluting agent that would be altering the water qua this lagoon, the source of which must be studied. Figure 10 reveals the concentrat NDCI in the second lagoon. The second monitored point was located in the Milluni Chico lagoon. Accord the monitored parameters, point 2 is the most visually and analytically contam point. This is attributed to the mining activity around it, and the indicators of the pr of heavy metals corroborate it. The Milluni Chico lagoon has an acid character throu the year without being affected by the time. It also has conductivities above 700 which suggests the presence of metal ions in the water. The lowest conductivities o the season with the most significant precipitation, which would present the dilution of pollutants. The concentration of metal ions (As, Zn, and Mn) is outside the perm limits. Likewise, there is a decrease in turbidity in May, which is congruent with t cline in the average concentration of the Normalized Chlorophyll-A Index in M Chico. Regarding the dissolved oxygen values, it can be said that they are acceptabl a tendency to decrease from May when the dry season begins. Figure 11 shows th centrations of NDCI in the lagoon.

Milluni Chico Lagoon
The second monitored point was located in the Milluni Chico lagoon. According to the monitored parameters, point 2 is the most visually and analytically contaminated point. This is attributed to the mining activity around it, and the indicators of the presence of heavy metals corroborate it. The Milluni Chico lagoon has an acid character throughout the year without being affected by the time. It also has conductivities above 700 uS/cm, which suggests the presence of metal ions in the water. The lowest conductivities occur in the season with the most significant precipitation, which would present the dilution effect of pollutants. The concentration of metal ions (As, Zn, and Mn) is outside the permissible limits. Likewise, there is a decrease in turbidity in May, which is congruent with the decline in the average concentration of the Normalized Chlorophyll-A Index in Milluni Chico. Regarding the dissolved oxygen values, it can be said that they are acceptable, with a tendency to decrease from May when the dry season begins. Figure 11 shows the concentrations of NDCI in the lagoon. Figure 11 reveals that the Milluni Chico lagoon has a more significant amount of chlorophyll-a in medium concentration in February. This decreases throughout the year, with a minimum peak in September. Contrary to the low chlorophyll-a concentration, which has its most down season in February, it tends to rise throughout the year with a maximum growing peak in September. It denotes that the dynamics of organic pollutants in this lagoon are not related to the dry or wet season detected in Milluni.
limits. Likewise, there is a decrease in turbidity in May, which is congruent with cline in the average concentration of the Normalized Chlorophyll-A Index in M Chico. Regarding the dissolved oxygen values, it can be said that they are acceptabl a tendency to decrease from May when the dry season begins. Figure 11 shows t centrations of NDCI in the lagoon.

Milluni Grande Lagoon
The third monitored point is located in the water storage dam located in the Milluni Grande lagoon. The water at this point should be of high quality, considering that it comes from a natural glacier and is used for public supply. Milluni Grande has an acid character throughout the year without being affected by the wet or dry season. The conductivities are above 1000 uS/cm. Once again, the lowest conductivity occurs at the time of most significant precipitation, denoting the dilution effect of the pollutants present. This is corroborated with high conductivities during the dry season. Both the pH and the conductivity suggested the presence of metal ions in the water, which is confirmed by the concentrations of arsenic, zinc, and manganese outside the permissible limits. In this sense, it is verified that this water body would be negatively impacted by the mining in the area.
Regarding turbidity, it is observed that this decreases in the rainy season and increases from May, which exposes the effect of dilution of pollutants generated by the rise in water in the dam. The dissolved oxygen values at this point are acceptable, with a tendency to decrease from May when the dry season begins. Figure 12 presents the concentrations of NDCI in the last lagoon.
Resources 2022, 11, x FOR PEER REVIEW maximum growing peak in September. It denotes that the dynamics of organic pol in this lagoon are not related to the dry or wet season detected in Milluni.

Milluni Grande Lagoon
The third monitored point is located in the water storage dam located in the M Grande lagoon. The water at this point should be of high quality, considering that it from a natural glacier and is used for public supply. Milluni Grande has an acid ch throughout the year without being affected by the wet or dry season. The conduc are above 1000 uS/cm. Once again, the lowest conductivity occurs at the time of m nificant precipitation, denoting the dilution effect of the pollutants present. This is c orated with high conductivities during the dry season. Both the pH and the condu suggested the presence of metal ions in the water, which is confirmed by the con tions of arsenic, zinc, and manganese outside the permissible limits. In this sen verified that this water body would be negatively impacted by the mining in the a Regarding turbidity, it is observed that this decreases in the rainy season a creases from May, which exposes the effect of dilution of pollutants generated by t in water in the dam. The dissolved oxygen values at this point are acceptable, with dency to decrease from May when the dry season begins. Figure 12 presents the c trations of NDCI in the last lagoon.   season. This could suggest that the factor influencing the low and medium concentrations of chlorophyll-a is independent of the Milluni season.

Discussion
Many previous studies have shown that mining is one of the activities with the most significant negative impact on the environment, especially on water. The ecological problem generated by unsustainable mining has been seen worldwide, but developing countries currently have the most significant vulnerabilities. This is due to many factors but in particular to the lack of resources to carry out complete controls and to be able to identify the sources and the degree of contamination.
Access to safe water is the current challenge for humanity, since it is directly linked to public health. This initial study exposes the water quality of the lagoons located in the upper part of the Milluni micro-basin. On the one hand, the first lagoon of the system (P1) does not present contamination according to international water quality standards, which implies that its water could be suitable for consumption after simple treatment. On the other hand, the Milluni Chico (P2) and Milluni Grande (P3) lagoons show alarming contamination. For P2 and P3, the measured pH values are around three and have high conductivities. In addition, toxic heavy metals were found for both points with values outside the permissible limits by the EPA and the WHO, exposing significant health risks if the water is used for public supply. The problem of water with compromised quality must be addressed in the first instance by reinforcing the monitoring of water bodies and subsequently seeking treatment alternatives to guarantee access to safe water for the population.
Although the traditional monitoring program designed for Milluni is functional and the results of its first year of implementation revealed the contamination of its lagoons, there are parameters within the permitted limits that are not a problem. However, over time, these parameters, heavy metals, could bioaccumulate and biomagnify within living organisms, leading to an ecological or public health disaster. Thus, this study exposes a future problem to be considered by decision-makers within water resources management.

Conclusions
The study presented the first approximation to the quality of surface water bodies in the upper part of the Milluni micro basin, which is the area affected by ancient and illegal mining. This is by correlating the results of the first year of the water quality monitoring program with satellite images obtained through remote sensing methodologies and programming on the "Google Earth Engine" platform. The four lagoons in the studied water system presented different qualities. Based on this, it can be seen that anthropogenic activities are not negatively impacting the first part, but they are affecting the second.
The part without contamination includes the first two lagoons, Pata Khota and Jhanco Khota. Still, as an important finding, the quality of the Pata Khota and Jhanko Khota lagoons do not have similar organic compositions (Figure 8). This would mean that a control point should be added within the Jhanko Khota lagoon for future monitoring to obtain adequate control of its quality. As for the most polluted part of the system, it was observed that it is the Milluni Chico lagoon, and the Milluni Grande lagoon is receiving its effluents. This compromises the water quality inside the water storage dam located at the end of the Milluni Grande lagoon, presenting a risk to public health.
Finally, combining the remote sensing tool with traditional monitoring was a more robust and affordable indicator of water quality. It is possible to constantly monitor surface water bodies without requiring additional financial resources because several remote sensing methodologies use resources such as free software or platforms. In addition, it was seen that remote sensing served to observe the bodies of water in a given period in a global manner, reducing the risk of error when evaluating the quality of surface water bodies that depend on a single control point.
Future research related to the analysis of surface water quality should consider the application of contaminant remote sensing methodologies as a mechanism to validate their results. In addition, monitoring programs looking for optimization could incorporate remote sensing as was mentioned above; these methodologies would not represent extra costs for the program but would help to carry out a constant and global control of surface water bodies.