Landscape-Scale Mining and Water Management in a Hyper-Arid Catchment: The Cuajone Mine, Moquegua, Southern Peru

1. Introduction

Developmental processes that impact the quality and availability of water in hyper-arid regions is an overriding concern when assessing environmental impact of any proposal to develop resources [1]. Determining the character of water chemistry and identifying the relative importance of both natural and anthropogenic factors in the contamination of water resources in such regions is a global priority, especially with predicted climate change impacts [2]. This is particularly important on the arid and mineral-rich western flank of the Andes, where two thirds of Peru’s population share less than 2% of the country’s water and where the world’s largest copper deposits reside [3,4,5,6]. Increasing pressure from mining on water resources poses a serious challenge to the ecology of these regions and regularly produces social conflict [7]. The projected nine-fold expansion of copper production in southern Peru, already the second largest producer after Chile, can only accentuate these challenges. Since trace metals are a major factor in the toxicity of contaminated water [8,9], mapping their presence in potentially affected river systems and understanding their baseline geochemistry is an important preparatory challenge. The direct impacts of increased mining can be partly mitigated by re-routing rivers using canals and diverting waste water to tailings dams, but the increase in mining inevitably increases the human population in remote mineral-rich regions, with a concomitant increase in urban and agricultural development, both of which are often associated with contamination from trace-metal and persistent residues [7,10,11,12].

Trace metal contamination in rivers systems is, in part, the product of natural acid rock drainage (ARD), in which the oxidation of sulfide-bearing minerals increases the acidity of run-off waters and enhances the further chemical breakdown of the underlying lithology. Hydrological systems dominated by acidic waters show increased trace metal mobilisation and metal-hydroxide precipitation [13]. The downstream eco-toxicity of these contaminants depends on their bio-availability and consequent accumulation in the food chain [14,15]. Polymetallic mining is a well-documented cause of increased ARD where it is more commonly known as acid mine drainage (AMD). Sulfide minerals and ores, in which the dominant metal ion can be iron, zinc, copper or nickel, are commonly associated with AMD. Chalcopyrite (a copper-iron-sulfide) is the most commonly mined copper bearing ore, and it is frequently found alongside other sulfides, making copper mines particularly prone to AMD [7,11,16,17]. In geologically active regions such as the Andes, ARD as a consequence of acidic epithermal waters entering the river catchment systems, lead to elevated chemical weathering rates and increased addition of ions to the rivers, the exact chemistry dependent upon the geology of the source region [18]. In combination, the effect of geology and anthropogenic elevation of weathering at disturbed mine sites in regions such as southern Peru may make water resources unsuitable for human consumption and agricultural irrigation, and therefore, water quality needs careful monitoring. Pesticide and fertilzer effluents and industrial and urban waste waters are also well-documented sources of trace-metal contamination. Arsenic is used in herbicides and pesticides; cadmium in batteries and plastics; chromium in dyes and tanning processes; lead in batteries, wiring and cabling; mercury and manganese in pesticides and batteries; and zinc is found in pharmaceuticals, dyes and batteries. All of these elements can leach into the wider environment from agricultural and urban run-off. The uptake of trace metals depends upon soil characteristics, plant species and the metal concerned [11,19]. The trophic transfer of trace metals, even at low concentrations, is known to have mutagenic effects [20].

In addition to the more global anthropogenic and natural geogenic factors that contribute to the determination of water quality, there are factors specific to areas in which mining and catchment transfer projects have been developed [16,21]. The Torata river sub-catchment and lower reaches of the Moquegua river drainage basin represent a region in which both a landscape-scale mine and a regional-scale catchment transfer project could reasonably be thought to impact catchment dynamics, including water resource quality and availability. (Figure 1) The exploitation of mineral resources increasingly necessitates environmental impact statements (EIS) with concomitant impact mitigation strategies which may become advisory or statutory obligations [22]. Development of landscape-scale mines increasingly requires such mitigation strategies to avoid the otherwise inevitable environmental consequences and political issues, although these may not be properly identified in many impact assessment processes [7,23,24,25,26]. The development of the Cuajone, Quellaveco and Toquepala mines in southern Peru have all involved the implementation of significant impact mitigation, with both river channel and waste material diversion with a complex network of sealed and open channels, tunnels and existing river channels [22]. Waste materials from the Cuajone and Quellaveco mines do not enter the greater Moquegua river drainage basin but are transferred to the tailings dams at Quebrada Honda (17 28 19.52S, 70 41 51.85W) at 1135m asl and Cortaderas 2 (17 12 02.49S, 70 41 51.85W) at 3086m asl, all of which can be readily seen on satellite imagery. Additionally, a tunnel diverting Torata river channel flow to the north of the principal mine pit at Cuajone was constructed to ensure that water quality from the headwaters is not affected by mine activity. Water scarcity in the region, combined with population growth and agricultural development has increased the need for water resources and driven the development of the catchment transfer Pasto Grande project. This a project transfers water from the headwaters of the Tambo river into the greater Moquegua river catchment in 262km of open channels, tunnels and existing river channels, that include the Sajena, Torata and Moquegua rivers. The long term impacts of this transfer on catchment ecology, hydrological dynamics and water chemistry have not been studied. Very few studies have been carried out on the impacts of catchment transfer projects but they are widely considered to have drastic impacts on hydrological regimes, water chemistry and aquatic biota [21]. The effect of the Pasto Grande project is the subject of further study.

As such, the Torata river was selected for several significant reasons: (i) the presence of the landscape-scale Cuajone copper mine, owned by Southern Copper, located in the headwaters, and the proposed expansion of mining in the region, gives the opportunity to look at the efficiency of water management measures associated with large-scale mining infrastructure; (ii) there is a global need for detailed studies of the controls on water chemistry in hyper-arid regions to guide resource management; (iii) the nature of the topography provides an ideal context in which to identify the sequential impacts of trace metal dispersion and concentration; (iv) the presence of two well spaced urban settlements at 2200m asl (Torata) and 1400m asl (Moquegua), allowing the identification of distinct chemical signals; and (v) the presence of two distinct and well separated agricultural areas associated with the above settlements Figure 1.

The general aim of this paper is to establish geogenic and anthropogenic determinants of water chemistry in the Torata river sub-catchment and lower Moquegua river, lying in a hyper-arid region in southern Peru. The catchment is periodically affected by El Niño–Southern Oscillation (ENSO) events with Pacific-derived moisture and regional atmospheric instability producing heavy precipitation and catastrophic flooding in the Moquegua region, as in 1998 [27] and more recently in 2016 [28].

Topographic conditions in the region provide an ideal transect from the high Andes to the Pacific coast with the Torata river and the lower Moquegua river crossing four distinct geological units, passing the landscape-scale Cuajone mine, the urban settlements of Torata and Moquegua and two agricultural areas (Figure 1). A sample-sites network was established throughout the Torata river sub-catchment, above and below the Cuajone copper mine, below the village of Torata and associated agricultural areas. Additionally, sites were adopted above and below the city of Moquegua and along the lower Moquegua river to the Pacific ocean. We describe the spatial and temporal variation in the chemistry of the river and its main tributaries, and use this to investigate the various contributors to the chemistry. Spatial patterns can be linked to processes specific to the upper, middle and lower catchments, particularly the underlying geology, any influence from mining and related acid rock interactions, agriculture, and hydro-fluvial processes such as ground flow, surface run off, distance from the ocean and elevation. Secondly, we discuss temporal patterns, linked to dry/wet seasons and El Niño and how these affect the baseline. Where the data allow we have indicated qualitative source apportionment based on concentration data. In addition to the overall aim we have, where the data allow, indicated qualitative source apportionment based on concentration data. In this we have attempted to achieve four specific objectives: 1) To identify the impact of the Cuajone mine on the geochemistry of the Torata river. 2) To identify and describe the hydro-chemical impact of the rural settlement of Torata and associated agricultural lands on the Torata river. 3) To identify and describe the hydro-chemical impact of the urban settlement of Moquegua on the Moquegua river. 4) To identify the the hydro-chemical impacts of the agricultural area lying along the lower Moquegua river valley.

2. Materials and Methods

2.1. Study Area

Hyper-arid regions are characterised by long periods of very low precipitation [18]. In the Peruvian Coastal Desert, the development of hyper-aridity co-evolved with the emergence of the episodic El Niño phenomena, creating a stable xeric environment that supports little vegetation, but which is affected by periodic catastrophic flooding and erosion [18,27,30]. The World’s largest supergene copper deposits are also located in this hyper-arid zone, stretching from southern Peru to northern Chile [31,32]. Currently there are three landscape-scale copper mines in the region; the Cuajone mine, the newly developed Quellaveco mine and the Toquepala mine. The study area focuses on the Torata river in the Moquegua drainage. This river rises in the headwater bofedales at 4590m asl at Pampa Titijones, and flows to the Cuajone copper mine at 3500m asl, before descending through the settlement of Torata, approximately 11km below the mine at 2200m asl. The channel then passes through a narrow canyon to join the Sajena river and, subsequently, the Moquegua river at 1290m asl below the city of Moquegua. Thereafter, the river flows through a narrow, intensively farmed, floodplain before reaching the Pacific Ocean at Ilo.

2.2. Study Sites

The sampling network was designed to cover the principal tributaries in the Torata subcatchment and lower Moquegua river catchment. Data discussed here are from 16 sites within the catchment, 11 along the Torata river, four lying at or below the confluence of the Torata river with the Moquegua and Sajena rivers, and one site close to the ocean (Figure 1). Precise GPS locations and elevations of the sites are given in Table 1. Site 0A is in the arid tropical zone, near Ilo on the coast. Site 1 is located in the upper tropics, 20 km below both the city of Moquegua and an extensive area of irrigated agriculture. The sites 1B, 1C and 1D are located at the confluence of the Sajena, Torata and Moquegua rivers. Site 1E, lies above all urban development and records the quality of the upper Asana River. Site 2, lies above all agricultural influence from the Moquegua city district, emerging from a canyon that starts in the Torata district. Sites 3 and 4 lie below the small rural town of Torata, and sites 5 and 5B above Torata, with 5B located immediately beneath the Cuajone mine. Directly above the Cuajone mine are sites 16, 17 and 18, and finally, sites 19 and 20 represent the source bofedale (peat wetland) and spring respectively.

2.3. Sample Collection

Three field investigations were undertaken: January 2017 (denoted 17-R), July 2017 (denoted 17-D) and January 2018 (denoted 18-R). The 17-R sampling campaign followed directly after an El Niño episode in which there was considerable flooding in the lower catchment. Despite this, it did not rain during the 17-R visit. During 17-D, the weather was dry and had been dry since the 17-R visit. Visit 18-R was nominally during the wet season but data collection preceded the onset of rains in the headwaters and as such many of the rivers retained dry season characteristics. Field visits were arranged to collect water quality data in-situ and water samples for later analysis. Prior to sample and data collection all field equipment was calibrated at the Universidad Nacional Autonoma de Moquegua (UNAM) laboratories. Owing to weather conditions and logistical constraints, it was not possible to collect samples at all sites on all field visits.

2.3.1. Water Quality Parameters and River Data Collection

At each site the following data were collected: GPS coordinates, atmospheric and water temperature, altitude, dissolved oxygen (DO mg/l and %), pH, conductivity, total dissolved solids (TDS) and oxygen reduction potential (ORP). A simple assessment of nearby land usage was noted.

2.3.2. Water Sample Collection

At each site, 1 litre composite samples were collected in pre-washed polyethylene bottles. Sub samples for Mercury (Hg) analysis were preserved with sulphuric acid and potassium dichromate to pH < 2, and samples for cyanide ($CN$) analysis were preserved with sodium hydroxide to pH > 12. Samples for total metal analysis were decanted into 50 ml tubes and preserved to pH < 2 with HNO₃ (99.9% trace metal basis) in the UNAM laboratory. All samples were packed in cool boxes for transport and stored at 4°C until analysis.

2.4. Analysis of Water Samples

Chemical analyses for 44 components including major cations, major anions, trace metals, nutrients and biological components (see Table 2) were conducted by Northumbrian Water Scientific Services (Newcastle, UK), Greenwich University (Kent, UK) and the University of Cambridge (Cambridge, UK) laboratories. Samples for metal analysis were treated to acid digestion at 105°C and analysed by ICP-MS (7500ce, Agilent, UK). The anions were measured using ion chromatography (Dionex IC with Ionpac analytical column¸ Thermo Scientific, UK) and nitrates were measured using discrete automated colorimetry (Aquakem 600, Thermo Scientific, UK). Hg was measured using atomic florescence (Millenium Merlin, PS Analytical, LOD: 0.003 $\mu$g/l) and total cyanide (CN) was measured using a segmented flow analyser (Skalar, Nederland, LOD: 0.02 mg/l). Total phosphorous (TP) was measured by ICP-OES (iCap 6500, Thermo Scientific, UK) and total nitrogen (TN) was measured using chemiluminescence (Shimadzu TOC-V, UK). Chemical oxygen demand (COD) was measured by digestion followed by colorimetry (Merck COD cell test kit, Germany, COD 5 – 80 mg/l). E.coli and coliform were quantified using multiple-tube fermentation techniques. All analytical methods were carried out according to EPA and ISO protocols in quality-controlled labs. All samples were analysed in duplicate for precision, and certified reference standards and internal calibration were used for accuracy measurements [33].

2.5. Statistical Analysis

Cluster analysis (CA) and principal component analysis (PCA) were used to quantify the relative dominance of various physico-chemical components of the river water with the aim of discriminating the main controls over the river quality. These methods are widely used in the analysis of correlations in multivariate data [9,16]. Both CA and PCA were applied to the z-score normalised value of the measured parameters [34]. z-score normalisation facilitates discussion of correlations between parameters with vastly different magnitudes or amplitudes of fluctuations. CA provides a framework for creating subsets (clusters) of samples with high similarities. The internal structure of the data is ordered according to a dissimilarity measure, grouping similar samples. The analysis presented here (Figure 7 and Figure 8), uses Ward’s hierarchical agglomerative CA method, together with a Euclidean dissimilarity measure [9,35,36]). This method finds the two closest samples (the data points with the smallest sum of squares of Euclidean distance), producing a secondary pseudo-sample. It then iterates this procedure until all the samples have been grouped, recording the dissimilarity distances, $\left\{d_i\right\}$, where i indicates the $i^{\mathrm{th}}$ link. In the CA dendrograms shown in this paper, the dissimilarity distance was reported as a fraction of the largest distance: $d_i/\max \left(\left\{d_i\right\}\right)$. The cophenetic correlation coefficient (CCC) was used to measure how well the cluster dendrograms preserve the pairwise dissimilarity distances between original, un-clustered, data [37]. A CCC of 100% corresponds to a dendrogram that perfectly preserves the original distances.

CA was employed in two ways in this study. Q-mode CA was used to find how the sample sites cluster with respect to the chemical parameters. This reveals spatial correlations of parameters over the set of sample sites. R-mode CA was used to find how the parameters cluster with respect to the sample sites. Clusters of similar parameters indicate that underlying factors affect them alike [38].

PCA is a methodology for the dimensional reduction of the number of parameters needed to describe a set of data [34,39]. The analysis begins with the calculation of the eigenvectors of the covariance matrix of the data, and the arrangement of them in descending order with respect to decreasing eigenvalues. If some eigenvalues are significantly larger (above some threshold) than the others, their corresponding eigenvectors are called “principal components”. By picking out the principal components with the largest eigenvalues, it is possible to describe the covariance of the data to a level given by the ratio between the sum of their eigenvalues and the sum of all the eigenvalues. In this paper, PCA was employed on clusters with similarly behaving parameters to understand the intra-cluster correlation.

2.6. Contextual Materials and Analyses

In addition to the above methodologies applied to our sample site data, we also refer to data collected under our longer-term campaign, and data collated from published papers and technical reports on water management in the region. This provides a useful context in which to fully understand both the geogenic and anthropogenic determinants of water quality and availability, while underlining the need for further studies.

3. Results

Water chemistry data from the three field seasons (17-R, 17-D and 18-R) are presented in Figure 2. Measured chemical parameters with their minimum and maximum values are in Table 2. The following sections discuss the temporal and spatial variations in river water chemistry during the three sampling campaigns.

3.1. Water Quality Parameters and Major Ion Composition

3.1.1. Headwater Sites Above the Cuajone Mine (Sites 16, 17, 18, 19, 20)

These sites were neutral or slightly alkaline (pH 7.84 to 9.45) with most sites recording a pH suitable for aquatic life. Site-19 showed the most alkaline conditions with pH values up to 9.45 and water temperature readings $\sim 5^{\circ }$C higher than ambient temperature. This site is a bofedal (Andean peat wetland) just below the source spring, site-20, and water was collected at midday in 17-R but early morning in 17-D. Both dissolved oxygen (DO) levels and pH at this site were higher at midday (12.4 mg/l and pH 9.45 in 17-R) than in the early morning (6.62 mg/l and pH 8.52 in 17-D) and, when taken together, are consistent with photosynthetic activity consuming CO₂ [40,41] and adding free oxygen to the water during daylight hours. However, vegetative growth was more evident during the warmer wetter season (17-R) than the drier colder season (17-D), as to be expected, and therefore increased photosynthesis in season 17-R may also have contributed to the higher pH and DO levels. Similarly, the high quantity of decomposing organic matter related to the rich floral composition of the bofedal, explains why the highest COD measurements in the Torata river (17 mg/l) system were observed at this site.

The headwater sites have the lowest concentrations of major cations (Na⁺≳ Ca²⁺≫ Mg²⁺, K⁺) and anions (${\mathrm{SO}}_4^{2-}$≫ Cl⁻> F⁻, ${\mathrm{NO}}_3^{-}$), with total concentrations in all ions increasing downstream (Figure 2) as expected. Na/Cl equivalent molar ratios at the headwater sites, Figure 3a, were significantly higher (10.6 - 3.4, 17-R, 14.4 - 5.3 17-D) when compared to sites downstream of the Cuajone mine (2.2 - 1.4 17-R, 2.7 - 1.9 17-D), with the highest Na/Cl ratios recorded at site-19 (the bofedal, 10.6 during 17-R) and site-18 (14.4, during 17-D). ${\mathrm{SO}}_4^{2-}$/Cl⁻ equivalent molar ratios are also significantly higher in the headwater sites, with site-18 recording the highest values (13.45 – 17-R, 18.02 – 17-D). Na/Cl equivalent molar ratio can be used as an indicator of ground water contributions to water bodies [42]. High Na/Cl ratio is an indication of a significant contribution from mineral-rich ground water sources rather than precipitation, which usually has Na/Cl ratios between 0.8 – 1.0 [42,43]. The high Na/Cl and ${\mathrm{SO}}_4^{2-}$/Cl⁻ ratios, when combined with the higher water temperature, are consistent with a geothermal contribution to source waters.

Following the El Niño, the bofedal, site-19, had the highest Fe and Al concentrations of the headwater sites (1.5 mg/l and 0.09 mg/l in 17-R respectively, Figure 4). This Fe concentration exceeds the safe limit of 1 mg/l for aquatic life [44]. Both the high pH and the organic processes driving reducing conditions in the waters of the bofedal encourage dissolution of iron, and aluminium solubility also increases with higher pH and in the presence of sulphate in the water.

3.1.2. Foothill Sites Below the Cuajone mine (Sites 5B, 5, 4, 3, 2, 1B, 1C, 1D, 1E)

pH measurements and field observations at sites below the Cuajone mine did not show evidence of AMD. pH values at sites 5, 5B, below the Cuajone mine, ranged from 7.98 – 8.56, little different from site 16 above the mine where pH ranged from 7.84 – 7.97. A healthy riparian eco-system dominates the valley at this altitude. This highly-sensitive environmental indicator suggests that the water quality in this area is good. Foothill sites-2, 5, 5B recorded elevated Al, Fe and Mn concentrations compared to other sites, and levels were particularly high at these sites and site-4 during 17-R, after the El Niño episode. Cretaceous (91-65 Ma) volcanic and sedimentary rocks, with Eocene intrusions (55-45 Ma), underlie the area (Figure 1) and represent the oldest geological units in the sample area. Chemical weathering leads to formation of clay minerals and insoluble oxides which are relatively high in Al and Fe. Changes in run off rate will have elevated the amount of very fine particulate matter in the rivers and the observed levels of iron and aluminium probably reflect this. The Asana river (sites-1B and 1E) also shows relatively higher concentrations of these elements than equivalent sites on the Torata and Sajena rivers. Figure 4 shows the spatial distribution of Al, Fe, Mn along the Torata and Moquegua rivers for the three seasons. Interestingly, although Al and Mn levels were not detectable and Fe only found in moderate concentrations (0.18 – 0.24 mg/l) above the Cuajone mine (sites 16-18), all three metals were found in higher concentrations below the mine (site-5, Al, Fe and Mn, 2.0, 1.5 and 0.33mg/l, respectively). USEPA guidelines stipulate that recoverable Al exceeding 0.087 mg/l (pH 6.5 – 9.0) and Fe exceeding 1 mg/l is unsuitable for aquatic life and the secondary drinking water standards specify that Al, Fe and Mn exceeding, 0.2, 0.3 and 0.05 mg/l, respectively, is undesirable in drinking water. Mn concentrations exceed the safe limits for irrigation use (0.2 mg/l). Although these metals are associated with both AMD and neutral mine drainage (NMD) [14,15,45,46,47], the concentrations detected were lower than those normally found under AMD or NMD conditions [48,49,50]. Under alternating flooding and draining conditions such as produced during an El Niño event and following dry seasons, changes in the pH of the water will greatly affect Al, Mn and Fe solubility via influencing either reductive dissolution or carbonate formation [51]. Therefore, it seems that seasonal variability controls the concentrations recorded in the data, and suggests that these metal concentrations probably derive from surface run-off, soil erosion and throughflow during the recent El Niño event [52,53,54]. Observed concentrations were thus considered to be the result of differential weathering of the underlying bedrock and trace metal mobilisation during the episodic rainfall during the recent El Niño event rather than mine leachate.

The metal concentrations decrease downstream from site-5, increase again at site-2 irrespective of the season (Figure 4). Such an increase in metal concentration could reflect evaporation as the river passes through the desert landscape, however the concentrations of conservative ions such as Cl⁻ did not show a similar spike, pointing to other origins. Cation enrichment from atmospheric aerosols is a possibility. Regional high pressure and sinking air in the foothills concentrates atmospheric aerosols which then settle out in the lower sections of the Moquegua, Sajena and Torata rivers, possibly affecting water chemistry. However, data from other studies [55] shows the river chemistry here is more significantly determined by contaminated groundwater. The aquifer is recharged by filtration of surface waters through the porous alluvial deposits of the Moquegua valley floor. The chemistry of these surface waters south of the city provide numerous pollutants from agriculture and the city itself. Therefore, the elevated cations (and bicarbonate) in the river at this altitude, are most likely due to recycled contaminated surface waters through the porous alluvial deposits.

Nitrate (${\mathrm{NO}}_3^{-}$) concentration in water can be used to assess chemical addition to the river from local agriculture. We observed raised ${\mathrm{NO}}_3^{-}$ concentrations, in the wet seasons 17-R and 18-R (0 - 16.64 mg/l and 0 – 7.9 mg/l), compared to the dry season 17-D (0 – 3.4 mg/l). The highest ${\mathrm{NO}}_3^{-}$ concentrations were recorded at site-3, located in the agricultural area downstream of Torata. The concentration drops at site-2 after the river has passed through 8km of desert canyon. Interestingly, the ${\mathrm{NO}}_3^{-}$ concentration at site-1C in the farmland around Moquegua, is similar to that at site 2 in both 17-D and 17R, possibly deriving from the greater emphasis on arable versus pastoral farming in this region. ${\mathrm{NO}}_3^{-}$ concentrations are controlled by increased surface run-off and leaching during the wet season, and relative concentration is linked to land use, being elevated in areas used for arable farming.

3.1.3. Sites Downstream of Moquegua City (Sites 1 and 0A)

The lower stretches of the Moquegua/Osmore river below Moquegua city record slightly alkaline (site-1, pH 7.81 – 8.68) conditions during all seasons. Water composition at site-1 (20km outside Moquegua and below an extensive area of more intensively farmed land) was dominated by Ca²⁺ (180 – 250 mg/l), Na⁺ (160 - 230 mg/l), ${\mathrm{SO}}_4^{2-}$ (370 - 470 mg/l) and Cl⁻ (260 - 360 mg/l). The concentrations of Na⁺, ${\mathrm{SO}}_4^{2-}$ and Cl⁻ were above or close to unsafe levels for drinking water (EPA and WHO standards, 200 mg/l, 250 mg/l respectively) for all seasons. The tributaries (sites 1B, 1C, 1D, 1E) feeding into the river above the confluence have lower concentrations and slightly more alkaline pH values (e.g. site 1B – 9.09, site 1E – 8.10). At sites 1C (Torata river) and 1D (Sajena river), above the confluence with the Sajena and Torata rivers, metal concentrations were below EPA and WHO safe limits set for drinking water, although relative concentrations are seasonally affected, and Moquegua sites 1B and 1E (Asana river) recorded lower cation and anion concentrations, in all seasons. Mn concentrations at site-1, during both wet and dry seasons (0.65 mg/l and 0.11 mg/l, respectively) are also higher than those recorded from the tributaries (see Figure 4): site-1B (0.064 mg/l), site-1D (< 0.05 mg/l) and site-1C (< 0.05 mg/l) with the exception of data from 18-R, when Mn concentration was 0.12 mg/l at site-1D, showing a significantly elevated contribution of this metal from the Sajena river. These results suggest that water from the Torata, Sajena and Asana/Moquegua tributaries does not significantly contribute to the elevated cation and anion concentrations at sites 1 and 0A, and the data is consistent with an extraneous source for the elevated chemical levels measured at site 1.

Cl⁻ concentration was above that considered safe for aquatic life over long periods of time during all seasons (chronic, >230 mg/l) and for irrigation of sensitive crops such as avocados [56], a common crop in the Moquegua region. The highest K⁺ concentration (8.7 – 11.0 mg/l) for all study sites was also found at site-1. As a result of these high concentrations of cations and anions, specific electrical conductivity (EC, 2040 – 2711 $\mu$S/cm) and total dissolved solids (TDS, 1020 – 1833 ppm) at site-1 were close to or above the maximum limits set by Peruvian water authority for agricultural use (category 3: 2500 $\mu$S/cm and 1800 ppm, respectively) [57]. EC and TDS values exceeded the safe limits for agriculture use in the dry season 17-D.

In general, cation and anion concentrations will gradually increase downstream along river systems due to accumulation from natural mineral weathering, evaporation and ground water discharge. However, the concentrations at site-1 were 2-3 times higher than those recorded at sites-1C and 2 in all seasons (Figure 2: line graphs) which suggests an additional process causing increased concentration. It is unlikely that evaporation alone (see Figure 5) could account for this dramatic increase so it is likely to be a result of changes in natural mineral weathering or the composition of groundwater discharge. Softer, more friable, sedimentary rocks of the Moquegua Group, which are more susceptible to physical and therefore chemical erosion, underlie the land between sites 1 and 2, in contrast to the more competent volcanic extrusive units found above Moquegua. Site 1 also lies downstream of the agricultural land below Moquegua. It is not possible to separate the relative influence of these two factors with the current data, but given the low precipitation rates (and therefore reduced weathering potential), such a dramatic increase in concentration probably reflects the influence of urban and agricultural activities in the lower valley.

${\mathrm{SO}}_4^{2-}$, Cl⁻ and K₊ were among the main contaminants at site-1 and they are commonly found in agricultural fertilizers, in the form of potassium sulphate (K₂SO₄, 50% K and 17% S), potassium chloride (KCl, 60% K) and magnesium potassium sulphate (MgK(SO₄), 22% S) [60]. These anions may also derive from urban waste [61].

Site-1 also showed the highest ${\mathrm{NO}}_3^{-}$ concentrations in the study site network (19.65 mg/l 17-R, 17.7 mg/l 18-R, 6.1 mg/l 17-D). ${\mathrm{NO}}_3^{-}$ concentrations from the tributaries were much lower throughout the sampling campaign (site-1D, 2.2 - 2.4 mg/l; 1B and 1E 0-0.96 - 0-1.5 mg/l). ${\mathrm{NO}}_3^{-}$ is highly soluble and commonly associated with the use of agricultural fertilizers and domestic waste [9]. To analyse the effect of the latter, samples from site-1 were further analysed for E.coli and thermal coliform, which are good indicators of sewage contamination. Both E.coli and thermal coliform were found to be below detection levels (<1.8 MPN) indicating little or no domestic sewage enters the river above site-1. Commonly used nitrogen-based fertilizers in the Moquegua district [60] are: urea (CO(NH₂)₂, which has the highest nitrogen content (46% N), ammonium nitrate (NH₄NO₃, 33%N) and ammonium sulphate ((NH₄)₂SO₄, 21% N). Another nitrogen-phosphoric fertilizer, commonly used is monoammonium phosphate (NH₄H₂PO₄, 10% N and 48% P), also a source of highly soluble phosphate [60]. Similar increases in total nitrogen (TN) and total phosphate (TP) were observed at site-1 ($r^2=0.99$, $p<0.05$), where aquatic plant growth was prolific. Correspondingly, an increase in Mn concentration at this site could also arise from agriculture since Mn is less readily complexed by organic ligands, resulting in high mobility in agriculture arras [62].

These data and the discussion above are consistent with urban effluent and agricultural run-off and leaching, especially during the wet season, as a major source of raised cationic and anionic concentrations at sites downstream of both the city of Moquegua and the associated downstream industrial agriculture area. However, Figure 2 shows that there are no significant correlations between the trace metals and cation-anion concentrations, suggesting that the actual system is complex with multivariate determination.

3.1.4. Endmember Controls on Water Chemistry in the Moquegua river Catchment

Water chemistry in rivers is controlled by precipitation, absorption, evaporation, crystallization and weathering. The Gibbs plot of TDS against the weight ratio of Na⁺/(Na⁺ + Ca²⁺) is a useful tool to identify these mechanisms in a river system [59,63,64]. The Gibbs plot for the Torata river sites is shown in Figure 5 (a). This shows that, irrespective of season, site-16 was the most Na-abundant (0.59 – 0.67 Na/Ca ratio) and site-5 was the most Ca-abundant (0.35 – 0.39 Na/Ca ratio) while the TDS remains between 98 – 115 mg/l at site-16 and 179 – 288 mg/l at site-5. This indicates a major change in water chemistry determinants between site-16 and site-5. The Gibbs plot confirms the dominance of precipitation and silicate weathering at headwater sites on water chemistry, and the change from Na abundance to Ca abundance reflects the change in underlying geology between sites 16 and 5. TDS naturally increases downstream and site-1 records concentration levels a factor of 1000 higher than the rest of the river system in all seasons. Evaporation becomes increasingly important in the lower parts of the river system but the data indicate that both agricultural leachates and evaporation play significant roles in the chemistry of the lower catchment.

3.2. Trace Metals in the Moquegua River System

Total concentration of trace metals (Cd, Cu, Co, Cr, Pb, Li, Mo, Ni, Se, U, V, Zn) along the Torata and Moquegua rivers for each season are shown in Figure 2 (bar charts). Samples were analysed for trace metals and other compounds associated with AMD and the natural geochemistry of the region (Table 2). Trace metals can have significant health impacts and are therefore of particular interest. Water analysis showed that Hg, Bi, Th and Te concentrations in all water samples were below detection limits (<0.06, < 2, <4 and <1 $\mu$g/l, respectively). The study also confirmed that cyanide (CN) remains below detection limits for all study sites (< 20 $\mu$g/l). In all water samples, Be, Cd, Co, Cr, Cu, Li, Pb, Mo, Ni, Se, U, V and Zn concentrations were well below safe limits for drinking and irrigation purposes, established by USEPA, WHO, FAO and the Peruvian water authority [44,57,65,66]. Generally, the concentrations of trace metals increase downstream. In our study catchment Be, Co, Pb, Se, U and Zn were only detectable below the Cuajone mine.

For both wet and dry seasons, a notable increase in total trace metal (Cd, Cu, Co, Cr, Pb, Li, Mo, Ni, Se, U, V, Zn) concentration was recorded below the Cuajone mine, with a 5-fold increase in the wet season and a 3.5-fold increase in the dry season measured between site-16 and site-5. The highest total metal concentration at site-5 was recorded after the El Niño episode (96 $\mu$g/l, 17-R). At site-5, significant levels of Li were recorded in 17-D and 18-R (33, 20 $\mu$g/l) while Cu was dominant in 17-R (53 $\mu$g/l: maximum recorded level in this study at any site, but still within safe limits). Pb was detected only at site-5 and site-4 in 17-R (2.9 – 3.9 $\mu$g/l). These trace metals occur naturally in river waters arising from typical Andean geological units but their concentration is generally increased by AMD. However, apart from a slight peak at sites 4 and 5 after the El Niño episode in 17-R, trace metal concentrations are not dramatically elevated above site 2 but follow a natural concentration profile concurrent with the cations and anions.

Spatially, total trace metal concentration peaked below Moquegua at site-1 during 17-R and at site-2 during 17-D. Trace metals are used in a wide variety of agricultural chemicals and domestic goods and, as such these peaks may derive from increased run-off and domestic waste in the lower valley, or increased weathering processes and the dominance of evaporation in the hyper-arid climes of the tropics. It is also true that trace metals can be easily transported by Fe, Mn and humic substances [53,67] which were found abundantly throughout the sample-site network. The trace metals data back up the conclusions from the major cations and anions above – that the dominant endmember in the lower catchment is the agricultural and urban runoff.

It should be noted that the hardness of water, determined by Ca²⁺ and Mg²⁺ ion concentrations, has a major influence on the toxicity of trace metals [68,69]. Water samples from the Torata headwaters (sites 16-20) showed limited hardness (17 – 44 mg/l CaCO3) while at lower elevations (sites-5 to 1C) water was permanently hard (75 – 378 mg/l CaCO3). Irrespective of the season, below Moquegua, water became extremely hard (548 – 764 mg/l CaCO₃). Soft water (<17 mg/l CaCO₃) is generally associated with increased trace metal ecotoxicity [68,69,70]. High levels of water hardness in the region almost certainly reduce the bioavailability, and therefore ecotoxicity, of trace metals in the Moquegua region. The elevated hardness could be an artefact of hyper aridity.

3.3. Arsenic in the Torata and Lower Moquegua River System

Arsenic (As) concentrations are shown in Figure 6. As levels along the Torata river were frequently above recommended levels for human consumption and irrigated agriculture (10 $\mu$g As/l, WHO and USEPA standards). The highest As concentrations were recorded at upper Torata valley sites 17 and 18 (14 - 21 $\mu$g As/l) where aresenic containing sulfidic minerals deposited by epithermal and hydrothermal processes associated with the active arc are more prevalent. Site-1D along the Sajena river also recorded elevated As in 17-D (11 $\mu$g As/l). By contrast, the Asana river sites 1B and 1E showed consistently low concentrations of As (2.1 – 4.3 $\mu$g As/l) for both 17-D and 17-R seasons. The highest As concentration observed in the Moquegua river, below the confluence of the Sajena, Asana and Torata rivers, was at site-1 during 17-R (10 $\mu$g As/l) which followed extensive flooding owing to the 2017 El Niño episode.

Generally, Peruvian rivers are known to contain As arising from geological sources [71]. The high Andean sites in the study area (16-20), are underlain by relatively young volcanoclastic deposits that are undergoing low temperature hydrothermal alteration and localised sulfidisation. Previous studies have reported that As is predominantly mobilised from such rocks and their weathering products [71]. Therefore, the elevated As levels found in the upper Torata sites are most likely to be geogenic in origin.

The dominant minerals found at the Cuajone Copper Mine are sulfides, chiefly chalcopyrite (CuFeS₂), chalcosine (Cu₂S), molybdenite (MoS₂) and pyrite (FeS₂), as well as arsenic-bearing minerals, including enargite (Cu₃AsS₄) and tennantite (Cu₁₂As₄S₁₃) [72]. These arsenic-bearing minerals, frequently associated with copper-porphyry deposits and generally marginal to the main orebody, develop during sulfide mineralization [72]. Mineralisation is very localised, concentrated within a 3km alternation halo, and does not extend significantly beyond this. It has been shown that weathering of porphyry Cu deposits and associated As-bearing minerals results in the dispersion of As into the groundwater and surface water bodies, depending on the redox state and pH of the system [73]. Site 5B sits just below the mine and does not show elevated levels of arsenic. However, sites further downstream show variably elevated arsenic levels depending on the season, consistent with remobilisation via groundwater and periods of increased run off and recharge through the shallow aquifer.

It is therefore likely, that the pattern of arsenic concentrations observed, Figure 6, reflects both local geochemistry and groundwater dispersion.

3.4. Multivariate Analysis of Water Parameters

Multivariate analysis was used to give an unbiased assessment of statistical correlations in our data. Correlations between water quality parameters and correlations between sites have been analysed using both CA and PCA. Our analysis took into account the 33 water quality parameters listed in Table 2. Parameters NO₂, CN, Be, Bi, Te, Th, Hg, Pb were not included because their concentrations were below measurement accuracy. Parameters CO₃, and DO were omitted because the datasets were incomplete.

Figure 7a presents a dendrogram of an R-mode CA for all sites visited in 17-R, 17-D and 18-R. It shows how the parameters cluster according to sites. It was produced with a CCC of 82%. Parameters that have a similar distribution across the sites cluster. At the highest level, we see three clusters. The right cluster contains major ions, and trace metals. These parameters arise from the natural water chemistry and dynamics in combination with the local geology. The left cluster contains nitrates, phosphates, trace metals, pH and COD. These parameters are indicative of agricultural or biological impact. The central cluster contains Al, As, Fe, Si, Mn and trace metals associated with of weathering of geological units with economic potential and / or urban effluent [74,75]. Peaks in these elements are also associated with periods of increased turbulence and greater suspended load (see above).

Figure 7b,c present a dendrogram from a Q-mode CA and corresponding Beck map (similar to a subway map with sites as stops) for the same 33 water quality parameters. They show how the sites cluster according to parameters. It was produced with a CCC of 84%. Sites with similar patterns of parameters cluster. At the highest level we see three clusters: the left cluster contains the two lowland sites downstream of Moquegua city, 0A and 1; the central cluster contains foothill sites; and the right cluster mostly contains the highland sites reinforcing the conclusion that chemistry is controlled by geographical location and regional inputs.

PCA analysis carried out on the same data shows that 74% of the variance can be accounted for by the first four components. The first component PC1 accounts for 46% of the variance, PC2 accounts for 12%, PC3 accounts for 9% and PC4 for 7%. Results are shown in Figure 8a - biplot of PC1 against PC2, and Figure 8b - PC1 against PC3. The parameter vectors in the biplots are coloured to indicate which parameter cluster they belong to in Figure 7a. As can be seen in Figure 8, the right cluster (red, orange; those controlled by natural hydrological processes and underlying geology) aligns well with PC1, the central cluster (light blue, grey; those controlled by the suspended load and increased runoff) aligns with PC2 and the left cluster (green, blue; those controlled principally by anthropogenic activity) aligns with PC3. As and Si align with PC4. Coloured circles are used in the biplot to indicate the site clusters from Figure 7b,c. The colour of the outer rings indicates the season.

3.4.1. Headwater Sites

In general, the headwater sites 16-20 cluster with site 5 (and site 4 in 17-R) below the Cuajone mine. This precludes chemical addition to the river through mining activity. PC1 and PC2 are low for the headwater sites (in contrast to sites 4 and 5 see below) suggesting the main controls on river chemistry up to this point are hydrological. Contribution from PC3 (chemicals associated with organics) is also low except for site-19, the bofedal, and this is interpreted as being naturally high levels of biogeochemical components in an Andean peatland. Data from the Asana and Sajena rivers was collected during 17-D and site 1E, on the Asana river, also clusters with headwater sites suggesting this river has similar chemistry and dynamics to the Torata river to this point. In 18-R headwater sites were inaccessible.

3.4.2. Foothill Sites

As mentioned above, sites 4 and 5 (5B) (between the mine and Torata) show strong alliance with the headwater sites but share some affinity with sites 2 and 3, below Torata. All four sites showed an average contribution from PC1 as expected from their location in the middle of the river system. During 17-R, Sites (4, 5) showed some contribution from PC2 whereas the headwater sites and sites (2, 3) showed average values. Elevated contribution from PC2 is consistent with increased runoff creating turbulent flow, following El Niño. All foothill sites showed an average value for PC3. Samples collected from the confluence of the Torata with the Osana and Sajena rivers (sites 1B, 1C, 1D) during 17-D, cluster with the other foothill sites (2, 3, 4). These sites showed average values for PC1, PC2 and PC3, consistent with the season being dry and the hydrographic location of these sites. The specific chemistry of the tributary sites in this area is seasonally regulated, sometimes sites share characteristics with the foothill sites and other times share characteristics with headwater sites. This is consistent with variable runoff rates in different seasons changing the exact chemical elements added to the river system at various sites. In general though, hydrographic processes determine the chemistry of the foothill sites (see Figure 5), and the multivariate analysis for these sites shows that the Cuajone mine had little or no effect on the river system over the observation period. The fact that this period included an El Niño episode shows that the measures to isolate the mine from the Torata river are effective.

3.4.3. Lowland Sites

The lowland sites (sites 1, 0) formed an isolated cluster irrespective of season. PC1 was high in all seasons consistent with these sites being towards the lower reaches of a river running through a hyper-arid zone with high levels of evaporation. In addition, elemental concentrations are elevated by accumulation of dissolved ions through chemical weathering of all three geological zones, from several agricultural regions and the two main urban centres; Torata and Moquegua. In 17-R PC2 was high indicating elevated wash through caused by the El Niño episode. In 17-D and 18-R PC3 was high, consistent with the agricultural activities along the river at that time.

3.5. Analysis of ANA and INGEMMET Data

In Figure 9 and Figure 10 we present a collation of data for the Torata river system including our surface water data, available survey data collected by ANA, the National Water Authority in Peru, from 2013 and 2014 [76,77,78] and groundwater analysed by INGEMMET around the same time [55]. Figure 9 shows the downstream concentration of a selection of major ions (calcium, sodium, sulfate) plotted against altitude. The data show an increase in concentration from source to sink, consistent with conservative behaviour of the ions controlled predominantly by cumulative addition of elements along the water course through weathering, and further concentration through evaporation. Points plotting to the right of the main cluster at just below 1000m asl are surface waters affected by contaminated groundwater (c.f. conclusions drawn by INGEMET [55]). The Cuajone mine is marked by a grey box with diagonal lines. If the mine had a significant impact on the surface waters we would expect to see a similar deviation of the surface waters to the right in the vicinity of the mine. In fact the two grey data points at 1100m asl, are samples of runoff analysed by ANA from the base of the tailing fan in 2013 and 2014. These data clearly show elevated concentrations of major ions. However, these elevated concentrations are not spatially distributed and data gathered only a few hundred meters from the tailing fan show normal concentrations. The ions are quickly scavenged from the surface water through sorption and precipitation and do not travel far from source.

Figure 10 shows the concentration of arsenic, cadmium and copper plotted against altitude. The concentration of these three trace elements in the water is a major concern for the people who live in the area. Data are not controlled by conservative processes as in the major ions, but are more dispersed. Of note are the digressions to higher concentrations around 4250m asl in all three elements. The higher altitude headwaters are underlain by young volcanic arc units, with very high permeability. There is extensive epithermal and hydrothermal activity associated with the active arc. This allows mobilisation of these cations via natural acidic chemical weathering of the mantle derived volcanic material, and rapid recharge into the surface waters as a result of the highly porous nature of the lithology (c.f. findings by INGEMMET [55]). The Cuajone mine is once again represented by a grey box with dashed lines. There is no significant increase in concentration of any of the three trace metals at the altitude of the mine. In contrast, the green lines just below 2500 masl and 1000m asl represent the start and end of the irrigated land between Torata and Moquegua. The data in all three plots is very dispersed in this region. Surface waters contaminated by the presence of the urban settlements and agricultural practices, are rapidly mixed into the ground water through the porous alluvial valley fill at this altitude, and recharge leads to increased concentrations in surface waters down stream. The concentrations of cadmium and copper are below guidelines for safe water throughout the river network, both for drinking and use in agriculture. The levels of arsenic are generally below the limits quoted by USEPA guidelines [65] except a few sites in the headwaters, where aggressive chemical weathering takes place, and also at times of elevated flow following heavy rains and/or flash floods in the agricultural plains.

4. Discussion

Hyper-arid regions rely heavily on river systems as the principal source of water. Changes in both river water availability and chemistry therefore have a direct and immediate impact on dependant systems in such regions. This is particularly true in the Andean transect in this study, where the Moquegua river system passes through the hyper-arid Atacama Desert and is used for human consumption and agricultural and industrial use throughout its length.

The location of the Cuajone mine above Torata required that both channel and waste material diversion mitigation strategies were essential, to avoid impacting water quality in heavily settled areas beneath the mine. To this end a diversionary channel was constructed around the projected mine pit at the Cuajone mine. Additionally, the construction of tunnels and concrete lined channels enables waste materials to be diverted to the tailings dams at Quebrada Honda and Cortaderas 2. The expansion of urban areas and the extension and modernisation of agriculture in the region necessitated the development of a regional irrigation project. The Pasto Grande catchment transfer project diverts water from the Tambo catchment, to the north, to sub-catchment tributaries in the greater Moquegua drainage basin. The catchment transfer project has significantly increased channel flow in the Sajena, Otora, Torata and Moquegua rivers and will undoubtedly add physico-chemical characteristics of the upper Tambo river to the Moquegua-Torata river system. Our data indicate a clear urban and agricultural signal in the lower parts of the Torata and Moquegua river sample sites. While data from Pasto Grande samples sites have not been analysed to identify a clear signal from that source in the lower Moquegua drainage basin.

Water availability in the Moquegua river system arises from seasonal rains and geothermal sources. In the month prior to the 17-R investigation, an El Niño episode had led to catastrophic flooding and a temporary superabundance of turbid river flow. River courses were altered, farmland destroyed and bridges and roads damaged. However, climatic conditions during the study period were generally dry and the rivers had already returned to dry season flow rates during the sampling campaign. The same was observed for 17-D and 18-R. The cycles of flooding and drought in this region, which disrupt water availability and complicate water resource allocation, are irregular and therefore difficult to predict.

Irregular seasonality affecting water chemistry, through the erosion and concentration of trace metal containing silts, further complicates water resource allocation and treatment. Elevated concentrations of Al in drinking water is often considered to increase the risk of Alzheimer’s disease, with prolonged exposure leading to systematic toxicity and renal failure; increased Fe and Mn is associated with aesthetic problems such as staining of laundry, unpleasant smell and taste [79,80]; elevated As concentrations in drinking water is associated with skin, liver and kidney cancers [81,82]. This study has shown significant seasonal variability in the levels of Al, Fe, Mn and As at any one site at any one time, making it hard to determine the social and economic cost associated with potentially harmful hydrochemistry.

The study indicates multi-factoral determination of river chemistry at many sites, including more significantly seasonality. This has significant implications for the management of water resources in hyper-arid regions. These areas are particularly prone to episodic and localised precipitation patterns, as well as distinct seasonality which may suffer increased climatic volatility in the future [83,84]. The data shows increased mobilisation of trace metals into the river system following episodes of more intense rainfall. These have a relatively short residence time and are scavenged from the water by sorption and precipitation, within a matter of months (17-R c.f. 17-D). However, management of the purity of the water in the lower catchment requires consideration of agricultural practises and supervision of the disposal of anthropogenic waste as these contaminated superficial waters are rapidly recycled into the river through the porous aquifer of the alluvial plain [55].

The desert micro-habitats in the region, although depauperate, have a high level of biological endemicity [85]. Hyper-arid environments are particularly sensitive to perturbation and recuperation is slow. Thus, the impacts on biodiversity and ecological integrity from irregular water contamination are disproportionately high; making the monitoring of water quality and control of contamination all the more consequential, for what is a scarce and very limited resource. The impact of elevated metal, metalloid and anion concentrations in this river system on ecologies is not yet known. It is of note that the Peruvian National Water authority (ANA) monitors only 98 of the 159 river catchments identified in Peru. Only 60% of these monitored rivers meet nationally agreed environmental quality standards [57].

This study finds three distinct geochemical and anthropogenic determinants to the water chemistry of the Torata River system; natural hydrological factors, addition of elements through increased run off and turbulence, and biogeochemical parameters associated with elevated organics, both natural and anthropogenic. Sites (16, 17, 18, 19, 20) fall within the arid temperate zone from 2800 – 4600m asl and are all above the Cuajone mine. Water quality here is determined by geological conditions and weathering processes. Cation, anion and trace metal levels were well below the safe limits for human consumption, with the exception of Arsenic concentrations, which exceeded safe limits for human consumption at sites 17 and 18. This can be attributed to geochemical weathering processes, but water chemistry also suggests the river is fed by Na⁺ and ${\mathrm{SO}}_4^{2-}$ rich ground waters, which can also be a source of As contamination [71]. Sites 3, 4, 5 and 5B, fall within the arid sub-tropics from 1400-2800m asl below the Cuajone Mine, town of Torata, and its associated agricultural areas. These sites showed moderate water quality with significant seasonal variation. While agricultural leachates are still detectable, naturally weathered contaminant loading is more prevalent. Increased surface runoff following El Niño, led to elevated Al, Fe and Mn in the water, where concentrations were recorded above recommended limits for domestic consumption and agricultural use. An increase in total trace metals, Li, Cu and Pb concentrations was also found below the Cuajone mine, but no indication of AMD was seen. Cation-anion composition at these sites showed water chemistry undergoing a transition between the Cuajone mine and Torata village, in which water composition changes from Na-dominant to Ca-dominant. This reflects a change in dominance from precipitative to evaporative processes in determining the hydrochemistry of the river. Sites 0A, 1, 1B, 1C, 1D, 1E and 2, fall within a hot, hyper-arid desert extending from sea level to 1400 m asl. Sites 1B, 1C, 1D, 1E and 2 show more chemical affinity with the sites in the arid subtropical zone, with water chemistry controlled though natural hydrological processes and minor concentration of some elements during periods of high runoff. Sites-1 and 0A, below Moquegua, exhibited the worst water quality irrespective of season. The data presented here is consistent with elevated levels of these chemicals from urban waste and agricultural leachates rather than natural weathering or hydrological processes. Significant levels of cations (Ca²⁺, Na⁺, K⁺) and anions (${\mathrm{SO}}_4^{2-}$, Cl⁻) as well as high EC and TDS levels were recorded. Intermittently, trace metals and metalloids (Mn, As) exceeded guideline values. Water quality in this section of the catchment did not meet strict criteria for safe domestic consumption and, on occasion, agricultural use.

The preceding discussion of chemical analyses and cluster analysis (CA) of site data show that the impact of the Cuajone mine is less significant than agricultural and urban activities in the contamination of this river system. Owing to hyper-aridity in the region, the low weathering regime reduces the mobilization of soil minerals and it is difficult to conclude that ARD occurs to any significant extent in this catchment, or that surface water quality has been greatly affected. Direct contamination as a result of AMD from the mine itself is reduced by the bypass of the Torata river around the Cuajone mine through a system of canals and diversions, and an upstream dam that largely prevents flooding from affecting the lower catchment. However, although AMD from open pit mining may be controlled during operation periods, it can cause long-term environmental problems after closure [14] and it is therefore important to develop a strategy to maintain good control over the mine drainage and tailing disposal to protect this drainage basin long after mine closure.

The data shows that the impact of agriculture and urban activities are far more influential than mining in determining cation-anion concentrations in the hyper-arid river system in this study. These rivers, in combination with the local climate and alluvial soils, provide ideal conditions for the development of highly-productive small and medium scale agriculture, which has been practised for centuries [66]. The modernization of agriculture has led to the rapid increase in use of pesticides and fertilizers, which in sensitive environments will have consequences for the quality of water resources and their value for domestic consumption and agricultural irrigation. Additionally, the use of river water and river channels to dispose of urban waste and industrial effluents is clearly deleterious to water quality and environmental integrity. Effective regulation of agrochemical use and tight control of urban and industrial effluent and sewage disposal in Peru will be key for the long-term environmental stability and management of a very scarce resource and therefore the long-term sustainability of development.

5. Conclusions

Our research shows that the geochemistry of the Torata river subcatchment is substantially determined by the underlying geological units and, in combination with weathering processes, characteristic of conditions found along the transect from the high Puna at 4500m asl to the lower coastal deserts. There does appear to be an indication of seasonal variation, perhaps linked to the episodic El Niño phenomena. Our data show that the Cuajone mine does not affect water quality parameters in the Torata river nor significantly alter the geochemistry of the river. We do not have data prior to the onset of mining at the site, and so cannot constrain any quantitative changes to the geological contribution to river chemistry as a result of the mine, but the data are consistent with chemical weathering being the dominant control on the fluvial geochemistry, both above and below the mine site during the period of investigation. It seems, at least in the time period concerned, that mine activity has had no adverse effect on water qualities necessary for aquatic life. Data does indicate a clear signal from the rural settlement of Torata and associated agricultural lands beneath the Cuajone mine. A clear geochemical signal is also discernible below the city of Moquegua in the Moquegua river, and also below the confluence of the Sajena and Torata with the Moquegua river. The data are consistent with contamination by urban and agro-chemical usage in the built up and agricultural areas lying below the city of Moquegua. This is the result of rapid recharge of contaminated surface waters from the urban and agricultural areas through the highly porous alluvial fill in the Moquegua valley. Additionally, significant increases in organic matter and sulphates were identified at El Conde, 18km below the city.

This study represents the time-limited findings of research undertaken from 2015 to 2020 in the Torata and greater Moquegua river systems in southern Peru and offers an insight into the determinants of water quality and availability in hyper-arid catchments in the region and, in furtherance of our understanding of geochemical dynamics in hyper-arid river systems, a unique case study. Our findings do indicate the need for further investigation into the significance of both mine operations and water management strategies, including catchment transfer projects, for water resources and environmental integrity. In southern Peru, this would involve investigating the wider environmental and resource implications associated with landscape-scale mine operations, including impact mitigation strategies employed prior to, during and long after mining operations cease, as well as the impact of the mine activity on groundwater flow and storage, headwater ecology and waste material diversion.