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.

The La Yarada aquifer is the primary water resource for municipal, irrigation, and industrial uses in the semi-arid Tacna, Peru. Presently, over-pumping has caused severe groundwater management problems, including the abandonment of saline water wells. This study presents multivariate analysis and chemical–isotopic trends in water to investigate seawater intrusion and hydrogeological processes affecting water quality. The chemical and isotopic analysis of water samples, collected in two campaigns in the dry (August 2020) and wet (November 2020) seasons, together with the 1988 data, were evaluated with a mixing model, cluster, and factor analysis. The hydrochemical and isotopic mixing model suggested the formation of a wedge with 20% seawater intrusion. The heterogeneity of piezometric map isolines corroborates the wedge formation associated with the groundwater movement. The spatial distributions of factors, FA1 and FA2, suggest two processes of seawater front movement: dispersion (diffusion) of chemical elements and different types of water mixing, respectively. At the edge of the La Yarada aquifer, the water head was relatively low, permitting seawater and freshwater mixing. On the other hand, along the sea-land boundary, the water head of the La Yarada aquifer was relatively high, avoiding seawater and freshwater mixing; however, the chemical species were migrating from the seawater to the groundwater due to the diffusion processes. The cluster 4 samples are in the region corresponding to the isotopic mixing process represented by the FA2, while cluster 4 describes the chemical diffusion process represented by the FA2. Thus, the integrated approach is helpful to assess the seawater intrusion mechanisms in coastal aquifers in a semi-arid region.

One of the most notable problems in the Atacama desert is the low occurrence of rainfall, which leads to a shortage of surface and groundwater for different users in the region. Therefore, the task of carrying out new exploration studies of potential groundwater zones (GWPZs) is of vital importance for decision-makers in water resources. The main objective of this research is to determine potential sources of groundwater using a Multi-Criteria Decision-Making technique with remote sensors. A method of exploration using the Analytical Hierarchy Process (AHP) techniques applied to remote sensing data is provided. The AHP method allows calculating the influence of multiple factors, and along with the GIS environment, a map of groundwater exploitation potential can be produced. The results indicating GWPZs showed four classifications of groundwater potential. The distribution shows 15.02%, 23.93%, 59.80%, and 1.25% of the total area with high, moderate, low, and very low potential, respectively. The results were validated with existing wells in the study area, offering an acceptance of 86.9%. This reveals the effectiveness and accuracy of the AHP-based GIS approach as a strategy for analyzing groundwater potential in arid zones. Similarly, the tested high GWP areas are helpful for the development and management of water resources in the Caplina basin.