Hydraulic recharge and element dynamics during salinization in an overexploited coastal aquifer of the world's driest zone: Atacama Desert.

The management of water resources in hyper-arid coastal regions is a challenging task because proper information regarding groundwater recharge and water budget is needed for maintaining the hydraulic balance in optimal conditions, avoiding salinization and seawater intrusion. Thus, this article deals with the estimation of the hydraulic recharge and the study of the effects of salinization on the dynamics of major and trace elements in an alluvial aquifer located in the world's driest zone, the northern Atacama Desert. The result of stable water isotopes (δD and δ18O) and tritium (3H) indicated that groundwater in the area is not recent, whereas 14C results estimated a groundwater residence time ranging between 11,628 and 16,067 yBP. The estimation of the artificial recharge coming from the urban water-supply-system leaks and wastewater/river-water/groundwater infiltration during irrigation was about 19.84 hm3/year, which represents an annual negative water balance of 177 hm3/year for the aquifer. The groundwater salinization triggered by seawater intrusion (up to 32.6 %) has caused the enrichment of Li, Rb, Ca, Ba, and Sr in groundwater by cationic exchange, where the excess of aqueous Na is exchanged by these elements in the aquifer sediments. Other elements such as B, Se, Si, and Sb are enriched in groundwater by ionic strength and/or anionic exchange during salinization. The heightened B concentrations derived from the B-rich alluvial sediments were higher than the limit suggested by international guidelines, representing a risk to consumers. Vanadium seems to be unaffected by salinization, whereas Pb, Mo, As, U, and Zr did not show a clear behavior during saline intrusion. Finally, this article highlights the consequences of conducting improper water management in coastal hyper-arid regions with exacerbated agriculture.

FlowSOM clustering - A novel pattern recognition approach for water research: Application to a hyper-arid coastal aquifer system.

Monitoring and understanding of water resources have become essential in designing effective and sustainable management strategies to overcome the growing water quality challenges. In this context, the utilization of unsupervised learning techniques for evaluating environmental tracers has facilitated the exploration of sources and dynamics of groundwater systems through pattern recognition. However, conventional techniques may overlook spatial and temporal non-linearities present in water research data. This paper introduces the adaptation of FlowSOM, a pioneering approach that combines self-organizing maps (SOM) and minimal spanning trees (MST), with the fast-greedy network clustering algorithm to unravel intricate relationships within multivariate water quality datasets. By capturing connections within the data, this ensemble tool enhances clustering and pattern recognition. Applied to the complex water quality context of the hyper-arid transboundary Caplina/Concordia coastal aquifer system (Peru/Chile), the FlowSOM network and clustering yielded compelling results in pattern recognition of the aquifer salinization. Analyzing 143 groundwater samples across eight variables, including major ions, the approach supports the identification of distinct clusters and connections between them. Three primary sources of salinization were identified: river percolation, slow lateral aquitard recharge, and seawater intrusion. The analysis demonstrated the superiority of FlowSOM clustering over traditional techniques in the case study, producing clusters that align more closely with the actual hydrogeochemical pattern. The outcomes broaden the utilization of multivariate analysis in water research, presenting a comprehensive approach to support the understanding of groundwater systems.

Hydrogeochemical Characterization and Identification of Factors Influencing Groundwater Quality in Coastal Aquifers, Case: La Yarada, Tacna, Peru.

The coastal aquifer La Yarada has anthropogenic and geogenic contamination that adversely affect the quality of groundwater for population and agricultural use. In this scenario, multivariate statistical methods were applied in 20 physicochemical and isotopic parameters of 53 groundwater pumping wells in October 2020, with the aim of characterizing the hydrogeochemical processes that dominate the groundwater of the coastal aquifer and the factors that cause them to optimize the effective management of water resources, delimiting areas affected by more than one salinization process. The samples were grouped into three clusters (C1, C2, and C3) with cluster analysis, the spatial distribution of C2 and C3 (reclassified in stiff diagrams), evidenced hydrogeochemical facies associated with the flow and recharge directions governed by the structural lineaments (NE-SO), favoring some areas more than others, arising different facies and hydrogeochemical processes. Factor analysis was applied from three different approaches: (1) main elements, (2) trace elements, and (3) physicochemical and isotopic parameters; exposing 6 distinguishable hydrogeochemical processes in the aquifer and factors that cause them: (i) salinization-marine intrusion, (ii) fertilizer leaching and dissolution of (Ca2+, Mg2+), (iii) wastewater mixture (NO3-), (iv) reducing conditions (Fe, Mn, Al), (v) contributions of (B, Sr), (vi) conservative mixtures and dissolution (As, F). It was validated with water quality indices (WQI) according to the national limits, delimiting 67 km2 parallel to the coast with "bad" to "very bad" quality for human consumption and unsuitable for irrigation according to the Wilcox diagram thus pre-treatment in this area is indispensable.

Predicting adverse scenarios for a transboundary coastal aquifer system in the Atacama Desert (Peru/Chile).

The Caplina/Concordia transboundary coastal aquifer system, located in the Atacama Desert, is the primary source of water supply for domestic use and irrigation for La Yarada-Los Palos (Peru) and Concordia (Chile) agriculture districts, and to a lesser extent, for Tacna province public supply use (Peru). Despite the scarce amount of rainfall (<20 mm/year) in the area and the limited recharge coming from the Andean highlands, this transboundary aquifer system has been overexploited mainly for agriculture since before the 2000s on the Peruvian side. Consequently, this has caused groundwater depletion and seawater intrusion. In this study, comprehensive hydrogeological information was integrated to understand the aquifer system's behavior and the effects to which it has been subjected to groundwater overexploitation. To that end, a 3D hydrogeological framework was developed using the LEAPFROG software and a constant-density groundwater flow model with equivalent heads was generated in FEFLOW software, which was adjusted with Monte Carlo analysis and conventional automated calibration. Finally, eight scenarios, considering various water resource management options proposed by the authority and potential climatic trends (CMIP6), were simulated from 2020 to 2040. The results showed that between 2002 and 2020, the increase in the seawater wedge and the average groundwater level decline were 216 hm3/year and 7 m, respectively. It is expected that the depletion will continue with a groundwater level decline between 5 and 8 m and an increase in the seawater wedge between 1120 hm3/year and 1175 hm3/year for the forecast period. The study concludes that the aquifer system will remain unsustainable for the next 20 years, regardless of the selected scenarios, and suggests that any mitigation measure requires the participation of stakeholders from Peru, Chile, and Bolivia.