Next Article in Journal
Hydrological Response of the Kunhar River Basin in Pakistan to Climate Change and Anthropogenic Impacts on Runoff Characteristics
Previous Article in Journal
The Influence of Fracture on the Permeability of Carbonate Reservoir Formation Using Lattice Boltzmann Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamics, Hydrochemistry, and Stable Isotope Geochemistry to Assess Temporal Behavior of Seawater Intrusion in the La Yarada Aquifer in the Vicinity of Atacama Desert, Tacna, Peru

1
Departamento de Ingenieria Civil, Universidad Nacional Jorge Basadre Grohmann, Tacna 23000, Peru
2
Geociencia, Universidad Politécnica de Nochixtlán “Abraham Castellanos”, Oaxaca 69600, Mexico
3
Departamento de Geologia-Geotecnia, Universidad Nacional Jorge Basadre Grohmann, Tacna 23000, Peru
4
Departamento de Recursos Hídricos, Universidad Nacional Agraria La Molina, Lima 15012, Peru
5
Escuela de Ingeniería y Ciencias, Tecnológico de Monterrey, Campus Monterey, Monterrey 64849, Mexico
6
Escuela de Ingeniería y Ciencias, Tecnológico de Monterrey, Campus Puebla, Puebla 72453, Mexico
*
Author to whom correspondence should be addressed.
Water 2021, 13(22), 3161; https://doi.org/10.3390/w13223161
Submission received: 31 July 2021 / Revised: 29 October 2021 / Accepted: 2 November 2021 / Published: 9 November 2021
(This article belongs to the Section Oceans and Coastal Zones)

Abstract

:
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.

1. Introduction

In coastal aquifers, seawater intrusion is a challenging problem affecting groundwater management [1]. However, the situation is more complex in the arid and semi-arid areas, resulting in aquifer heads declination, groundwater quality deterioration, lower crop yields, and ecosystem degradation [2]. With rapid industrial growth and agriculture intensification, groundwater sustainability has become a significant issue in Tacna, Peru [3,4]. The seawater infiltration into the aquifer due to overexploitation with high pumping wells [5,6] had forced the abandonment of wells in the region with increased salinity [7,8]. The encroachment of saline water into fresh groundwater has been a subject of extensive study for well over a century [9]. Since 1920, various graphic and statistical techniques have facilitated the homogeneous grouping of groundwater [10,11]. Multivariate statistical methods, cluster analysis (CA), and factor analysis (FA) effectively manipulate, interpret, and represent groundwater pollutants and geochemistry data [12,13].
Villegas et al. [14] used the principal component analysis to evaluate the dissolution and cation-exchange processes controlling the groundwater chemistry in the Uraba aquifer. Black et al. [15] assessed the thermal water origin in the Dublin carboniferous basin (Ireland) by implementing CA and FA. Werner and Simmons [16] illustrated the first-order assessment of seawater intrusion through a simple conceptual framework, classifying aquifers as flux-controlled and head-controlled systems. Binda et al. [17] argued for applying an integrated interdisciplinary approach and illustrated it by evaluating potentially toxic element sources in a mountainous watershed. The same procedure defined the natural and anthropogenic contributions in the winter snow of the Dolomites, Italy [18].
The present study implemented the integrated interdisciplinary approach [17] to characterize the spatiotemporal variations of water level and salt content in the La Yarada coastal aquifer in the vicinity of the Atacama Desert, Tacna, Peru. The study compiled the hydrodynamic, hydrochemical, and isotopic data (δ18O and δ2H) from the literature. Water samples of ice, rivers, springs, and groundwater wells in the basin were collected in two campaigns. A multivariate statistical analysis was performed to characterize the dominant factors to understand the mechanism of seawater intrusion in the aquifer. Likewise, an isotope–hydrochemical mixing method was employed to identify the salinization processes due to seawater intrusion in the aquifer.

2. Study Area

2.1. Location and Climate

The Atacama Desert, a cold arid region in northern Chile, is approximately 1100 km long from north to south, and the north part of the desert is in Peru. Figure 1 shows the study area of the Caplina hydrographic basin, which encompasses a delta that caused the formation of the La Yarada coastal aquifer. Desertification started 14 Ma during the global climatic desiccation [19]. Similarly, the Barroso mountain range separates the Caplina and Maure basins [20]. The coast and the entire western slope of the mountain range are in a rainy region of around 400 mm. The annual rainfall in the Caplina basin is only up to 10 mm, except for the years 2019 and 2020 [21].
The main economic activity in the area is agriculture, covering more than 40,000 ha cultivated on the surface of the La Yarada aquifer in 2020 [21,22]. Consequently, the groundwater has a severe imbalance problem in the coastal zone to meet the population irrigation and domestic use demand [20]. The principal causes of the aquifer’s progressive salinization process are quality deterioration due to seawater intrusion and depletion of the aquifer [22,23]. Consequently, there is a crisis of governability and governance in groundwater use [23]. The excessive and uncontrolled groundwater pumping fulfilled the local water resources demand in the last century [19,21]. In 2018, illegal pumping wells exploited more than 160 Hm3 per year [22], far exceeding the aquifer capacity.

2.2. Hydrogeological Settings

The La Yarada aquifer is mainly of quaternary alluvial origin. Its shape is a rectangular polygon with a flat bottom and steep and abrupt flanks. Downstream, the ejected cone of the Caplina river (Peru) constitutes a physiographic unit that begins in the Magollo gorge and progressively widens downward in a delta unit, reaching the beach line [24]. The aquifer is delimited from Calientes to the beach line from the northeast to the southwest by rocky outcrops with incipient wind cover and volcanic ash deposits [25].
During the last two decades, the water level has declined due to groundwater overuse, and consequently, there was seawater intrusion in the aquifer. In the coming years, the water level will continue to drop from 0.23 to 0.38 m/year, and saline intrusion will be five times higher (56–59 hm3/year), increasing groundwater quality deterioration [23,26].

3. Materials and Methods

3.1. Monitoring Network and Sampling

The hydrochemical and isotopic monitoring network, distributed in the Caplina basin, was first established in an IAEA project in 1988 [27]. In this study, the first campaign was performed in August 2020 (dry season), while the second was in November 2020 (wet season) [28,29]. One hundred sampling points were selected, distributed in three basins: 86 points from the Caplina basin, 9 points from the NE Maure basin, and 5 points in the NW Sama basin for the chemical and isotopic analyses (Figure 1). Of the total samples, 76 are for agriculture use, 9 for domestic use, 5 for cattle-ranch use, 3 for urban use, 1 for industrial use, and 6 for no use due to high contamination.
Each sample was conducted in 4 days from the highest part of the Maure basin at 4500 m above sea level (masl) to the beach at 0 masl by two brigades. The sample points of the Maure and Sama basins were selected to observe the influence of interaction or contribution from the aquifers to the La Yarada aquifer of the Caplina basin.
For the present study, we also compiled the isotopic data of 44 groundwater and three from the Caplina and Uchusuma rivers from the report of the IAEA project [24].

3.2. Analysis and Data Quality

Cation (Ca2+, Mg2+, Na+, K+) analysis of water samples was conducted by inductively coupled plasma mass spectrometry (ICP-MS) following the method ISO 17294-2. The anions (SO42−, Cl, NO3, F) were analyzed on the ion chromatograph, using EPA methods 300.0. The carbonic species (HCO3) and alkalinity were determined by the acid-base titration method. The δ18O and δ2H of water samples were measured by the analyzer ABB-LGR model 912-0008 with analytical precision ± 0.1‰ and ± 1.0‰, respectively.

3.3. Statistical Interpretation Procedure

The univariate, bivariate, and multivariate analyses of the chemical and isotopic data were performed with the computer code written during this study in Python. The following 14 variables were considered: electrical conductivity (EC), pH, total dissolved solids (TDS), Ca2+, Mg2+, Na+, K+, HCO3, SO42−, Cl, NO3, F, δ18O, and δ2H to identify chemical groups in groundwater.

3.4. Groundwater and Seawater Mixing Model

A mixing model was proposed to quantify the seawater component in the Caplina basin. It consists of establishing the chemical and isotopic compositions of end members: (1) local fresh groundwater (Xgnd) and (2) seawater (Xsea) [30]. The fraction of seawater component ( γ ) in each groundwater well is calculated as
γ = X m i x X g n d X s e a X g n d    
where X is the concentration of chemical ( C l , N a + ) and isotopic ( δ 18 O ) species in the saline groundwater samples. The subscript (mix) indicates the mixture component (i.e., the La Yarada aquifer) groundwater.
Figure 2 presents the conceptual model of the negative effect on the hydraulic barriers [31] and geometry of the freshwater–seawater contact after 13 years of operation of extraction wells [32]. On defining the end members (i.e., local fresh groundwater and seawater), Equation (1) is used to calculate the seawater fraction in each groundwater well in the aquifer. Thus, the mixing model permits in defining the fresh–seawater contact geometry.

4. Results and Discussion

4.1. Spatial and Temporal Univariate Trends in Groundwater Chemistry

Table 1 presents the data statistics of chemical and isotopic analyses of samples collected in 2020. The physical chemistry of saline groundwater has more significant variability in EC, which is related to TDS. Figure 3 shows the spatial and temporal distribution of Na+ and Cl in the aquifer for 1998 and 2020. These parameters are generally conservative during the mixing of different types of water [30]. The iso-value curves indicate that the high Na+ and Cl concentration zones had evolved from W to SW along the aquifer edge in contact with seawater.
Similarly, the saline seawater front is moving slowly inside the coastal area. The phreatic level inside the Caplina basin was decreasing due to extensive pumping up to 2010. However, high precipitation in 2019 and 2020 was favorable for improving the phreatic level profile [33]; however, there was no significant change in the chemical and isotopic compositions of the La Yarada groundwater.

4.1.1. Stable Isotope of Caplina Basin

Figure 4 shows the isotopic values of different types of waters in the Craig diagram [34]. Ninety percent of the water samples lie around the GMWL, and the data were fitted to the LMWL, δ 2 H = 7.06   δ 16 O + 0.97 , with R2 = 0.97 in 2020. There is a high scattering in the values of 1988; however, the values are around the GMWL.
The snow samples (δ18O = −13.9‰ and δ2H = −99.6‰) were collected at the elevation (Lat. = 17°37′ and Long. = 69°48′, Elev. = 5164 masl). The highest peak of the mountain is at the elevation (Elev. = 5800 masl). So, the isotopic composition of rainwater at the mountain range may still be lighter and close to that of the Maure aquifer. The Caplina basin is probably recharged by the river waters and evaporated during the river flow trajectory. The isotopic values of Caplina and Maure Basin are on the trend line. The isotopic values of snow suggest the Barroso Mountain range as a separation boundary for the basins. The groundwaters of the La Yarada aquifer are the recharge of the evaporated mountain rain waters.

4.1.2. Isotopic Composition Interpolation

The isotopic composition of groundwater in arid regions may be different from the composition of the recharge zone [28]. Figure 5a presents isotopic composition iso-lines that separate the isotopic system of the Sama basin in 1988 in the littoral part; however, the iso-line of −11‰ stands out at 2 km from the coastline and −12‰ at 5 km, generating an isotopic enrichment anomaly due to seawater influence in 2020 (Figure 5b).

4.2. Statistical Correlation

Table 2 shows the correlation matrix among the 14 variables. The correlation coefficients between the physicochemical parameters and isotopes of δ2H and δ18O have low values except for fluorine. Likewise, the strong positive correlations in pairs TDS and Cl (0.98), TDS and Ca2+ (0.92), TDS and Mg2+ (0.93) reflect the dependence of salinity increase with the primary ions of seawater Cl (0.98), Ca2+ (0.92), and Mg2+. However, the increase in TDS is not necessarily due to seawater intrusion. It is crucial to consider the solid correlation for salts corresponding to seawater [35,36], highlighting the exclusive non-correspondence of magnesium and calcium coming from calcareous rocks. Similarly, the weak correlations of F, NO3, and HCO3 with electrical conductivity (EC) and TDS reveal a low contribution in groundwater salinity compared to the major ions Cl, Ca2+, and Mg2+.

4.3. Multivariate Analysis

4.3.1. Cluster Analysis

The cluster analysis is a statistical tool to classify the principal data groups according to a Euclidean distance dissimilarity measure [36]. According to their parameters, the sample classifications are known as Q mode classifications. The Ward method uses Euclidean distances to measure similarity and has a small distortion effect in space [12,13]. The 14 chemical parameters (EC, pH, TDS, Ca2+, Mg2+, Na+, K+, HCO3, SO42−, Cl, NO3, F, δ18O, and δ2H) of 53 groundwater samples of the La Yarada aquifer were grouped into four clusters of hydrogeochemical facies according to their similarity in chemical signature (Figure 6).
Table 3 explains the relative abundance of ions. Cl and SO42− are the dominant ions in all clusters except for cluster 2, which shows the domination of Cl and Na+. A relationship (<0.05) of the Na+/Cl ratio reveals the influence of very salty water from seawater intrusion [36]. Figure 7a presents the spatial distribution of the chemical signature of the four clusters. Cluster 1 is parallel to the coastline and has high TDS, confirming the high salinity of the coastline.
The chloride abundance in clusters 3 and 4 indicates a possible zone of seawater intrusion. Still, it may also result from clay dissolution, anhydrite from sulfide leaching, sewage infiltration, agricultural inputs such as sulfate, and ammonium return of irrigation [18,23]. The spatial distribution of these groups covers most of the aquifer, which indicates a double source of sulfates as irrigation return and seawater intrusion.
The primary domain of cations is Na+ > Ca2+ >> Mg2+, and in anions: SO42− >> Cl > HCO3 (i.e., calcium sulfate waters). However, the occurrence of more chlorinated sodium samples is observed in clusters 3 and 4. The distribution of group 4 (Cl >> Na+) towards the coastline confirms that a large volume of saltwater intrusion is estimated at the SW of the aquifer [26].
For clusters 1 and 2, the water samples close to the coastline experience Cl–SO42− ion exchange (Figure 7b). These waters contain the highest TDS of the entire sampling network and are influenced by very salty water as a product of seawater intrusion [23,26,32]. However, cluster 2 (SO42− ≥ Cl >> Ca2+) waters reappear in the middle part of the aquifer, with hydrochemical values similar to cluster 3 (SO42− >> Cl). Sulfates predominantly salinize the water of cluster 2 due to the SW underground flow that recharges the aquifer [20]. The spatial distribution of water of cluster 3 covers most of the aquifer. There were two sources of sulfate: sewage infiltration and agricultural fertilizer inputs (i.e., irrigation water infiltration) [14,15].

4.3.2. Factorial Analysis

Factor analysis permits the identification of the correlations between the physicochemical components of the groundwater samples. Table 4 presents the three independent factors, classified as strong, moderate, and weak [37], corresponding to values (strong, moderate, and weak correlation) as >0.75, 0.75–0.50, and 0.50–0.30, respectively.

4.4. Conceptual Model for Seawater Intrusion in Groundwater Salinization

The water mixing model water delimited the variability of the chemical concentration in the mixture of seawater and freshwater [38]. The La Yarada aquifer groundwaters were clustered in four groups. The factorial analysis converted the 14 parametric chemical and isotopic data into three critical factors: FA1, FA2, and FA3. Figure 8 shows the spatial variation of the first two factors, explaining 44.38% and 21.37% of the total variance (Table 4). The spatial distribution of the chemical signature of the four clusters is also plotted in Figure 8. The factor FA1 is characterized by a positive correlation of Cl, TDS, Mg2+, Ca2+, K+, and Na+ and moderately by SO42− and EC. FA1 was named the “seawater chemical front”, because it represents the primary ions of seawater. Figure 8a shows the spatial distribution of the factor scores. The high-value regions represent high salinization regions (Cl, TDS, Mg2+, Ca2+, K+, and Na+), and the low-value regions denote the discharge regions for these chemical species.
The FA2 factor characterized the positive correlation of isotopes (δ2H, δ18O) and ion F (Table 4); therefore, it was named the “seawater isotopic front”. Figure 8b represents the spatial distribution of FA2 scores.
The spatial distributions of 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 (Figure 8b). 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 FA3 factor explains 12.38% of the total variance. Moderate loads of HCO3 and NO3 characterized with a negative correlation of pH and NO3. It means an entry of municipal wastewater infiltrated from the forest in the center of the aquifer since the forest was irrigated with secondary treated wastewater. The HCO3 may also be originated from the dissolution of carbonates. The negative correlation between HCO3 and pH would indicate the little significant impact of pH in the dissolution of carbonates.
There is a good relationship between clusters and factors. The data points corresponding to cluster 4 are located in the region where the FA2 seawater isotopic front process dominates. Similarly, cluster 4 represents the process corresponding to the FA1 seawater chemical front.

4.4.1. Sama Basin Underground Connection

The Sama basin has the δ18O values (−10 to −12‰). The multitemporal iso-value maps in Figure 5 show the variation from −11‰ to −10‰ in 32 years (i.e., isolines in the SE region of the aquifer). This theory is reinforced when observing the distribution of cluster 2 (Figure 7) in superposition to the higher scores in FA2 (Figure 8b). It may be inferred that there is an underground flow to recharge the wells to the east of the aquifer. Under the above postulation, the isotopic composition of the water, once in the aquifer, does not change quickly, except for mixing with another source of different isotopic values [35].

4.4.2. Rock Dissolution

The δ18O value tends towards more positive values that are influenced by the isotopic exchange in the water-rock reactions, favored by the low slope (3–5%) in the aquifer [39,40]. The soil’s primary source of fluoride, emphasized by the factor FA2 and cluster 4, is clay minerals [41]. Similarly, paleochannels [39] were identified in piezometric records of the aquifer, formed of silt and clay banks. These paleochannels were formed due to weathering, and the leaching process occurred during the groundwater flow. When fluoride-rich minerals are in contact with highly alkaline water, they release fluoride into groundwater through hydrolysis, replacing the hydroxyl ion (OH) [42]. It is worth highlighting the considerable positive correlation between fluoride and pH in factor FA2 (Table 4). The rock dissolution process is not only present in the aquifer strata but throughout the Caplina basin of the sulfated calcium family from the oxidation of sulfides due to the hydrothermal alteration in the mineral occurrences of the intrusive Yarabamba and Challaviento [40].

4.5. Seawater Salinization Processes

Figure 9 presents the conceptual model of the negative effect on the hydraulic barriers [32] and geometry of the freshwater–seawater contact after 13 years of operation of extraction wells [33].
The γ-factor for parameters (Clsub, Na+sub, and δ18Osub) was calculated by Equation (1), considering the distance as their position perpendicular to the coastline in the NW direction. The factor represents the fraction of seawater origin in the groundwater samples. The data were fitted in the logarithmic functions (Figure 9a). Despite having a positive correlation, the samples are very dispersed to the trend consequence of seawater salinization anomalies. In contrast, the rest of the samples with γ = 0–0.05 depends on the aquifer’s hydrogeological and paleo-hydrogeological regimes. In this sense, in Figure 9a, emphasis was placed on relating the anomalies of 1988 and 2020 with factors γ = 0.20–0.30, configuring boundary lines of the facies with the highest seawater proportion (Figure 9b).
A typical mixing interface was obtained for the data of the year 1988. For the year 2020, at 11 km, samples that differ from the behavior of a salinized interface are distributed, generating a behavior similar to that obtained after applying hydraulic barriers [32]. Based on these already established models of the effects of negative hydraulic barriers [31,32], a 3D model was generated to visualize the “pseudo hydraulic barriers” (Figure 9b).
The extensive exploitation would drive this behavior at 11 km from the coastline that overlaps in the range of exploited volumes (Figure 9c). This anomaly maintains the isotopic signature in the aquifer (Figure 3) since exploitation is continuous every month of the year. Therefore, they do not allow the majority salinization of Cl and Ca2+ in more than 30% of the aquifer. The graph of the γ fraction versus littoral distance (Figure 9a) appears to be a process called saltwater up-coning [43] as a result of the intense pumping of water in the area (Figure 9c).

5. Conclusions

The multivariate analysis and chemical–isotopic trends studying different types of water in the semi-arid Tacna, Peru, to investigate seawater intrusion and hydrogeological processes affecting water quality are summarized in the following outcomes:
  • No change in the chemical and isotopic analysis of groundwater of two samplings, performed in the dry (August 2020) and wet (November 2020), indicate practically no direct influence of precipitation on the La Yarada aquifer characteristics.
  • The hydrochemical and isotopic mixing model supports 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 factor analysis reduces the chemical and isotopic parameters into the FA1 seawater chemical front and the FA2 seawater isotopic front.
  • In the La Yarada aquifer, the seawater and freshwater mixing is dominant at the SW of the aquifer; however, the chemical species migration is along the sea-land boundary. 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.
  • The intrusion of the mixing zone into groundwater behaves differently in the unconfined zone due to the variability of exploitation volumes. In the area proximal to the coastline, the extracted volumes do not exceed 0.5 Hm3. In contrast, more wells than those registered exist to pump 2.5 Hm3 water from the Caplina aquifers nearby the wedge formation.

Author Contributions

Conceptualization, E.P.-V., M.P.V. and A.V.; software, S.C. and A.V.; validation, E.P.-V.; formal analysis, E.C. and M.C.; writing—original draft preparation, A.V.; writing—review and editing, J.A.T.-M., A.M. and J.M.; supervision; project administration, E.P.-V. All authors have read and agreed to the published version of the manuscript.

Funding

Universidad Nacional Jorge Basadre Grohmann, Tacna, Perú.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable, part of the data corresponds to reports from Peruvian public institutions that under agreement provided such information with academic purposes, but the majority of data is in Table 2.

Acknowledgments

This work was developed within the framework of the research project, “Integration of hydrodynamic, hydrochemical and isotopic methods to specify the operation and sustainable management of the La Yarada aquifer, Tacna, Peru” founded by “Fondos de canon, sobrecanon y regalías mineras de la Universidad Nacional Jorge Basadre Grohmann, Tacna, Perú”. The constructive comments of two anonymous reviewers helped to improve the content and presentation of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Datta, P.S. Water Harvesting for Groundwater Management: Issues, Perspectives, Scope, and Challenges; John Wiley & Sons Ltd: Hoboken, NJ, USA, 2019. [Google Scholar]
  2. Chen, J.; Wu, H.; Qian, H.; Li, X. Challenges and Prospects of Sustainable Groundwater Management in an Agricultural Plain along the Silk Road Economic Belt, North-West China. Int. J. Water Resour. Dev. 2018, 34, 354–368. [Google Scholar] [CrossRef]
  3. Grosjean, M.; Veit, H. Water Resources in the Arid Mountains of the Atacama Desert (Northern Chile): Past Climate Changes and Modern Conflicts. In Global Change and Mountain Regions; Springer: Dordrecht, The Netherlands, 2005; pp. 93–104. [Google Scholar] [CrossRef]
  4. Jordan, T.E.; Herrera, L.C.; Godfrey, L.V.; Colucci, S.J.; Gamboa, P.C.; Urrutia, M.J.; González, L.G.; Paul, J.F.; Jordan, T.E.; Herrera, L.C.; et al. Características Isotópicas e Implicaciones Paleoclimáticas Del Evento de Precipitación Extrema de Marzo de 2015 En El Norte de Chile. Andean Geol. 2019, 46, 1–31. [Google Scholar] [CrossRef]
  5. Naciones Unidas La Agenda 2030 y Los Objetivos de Desarrollo Sostenible: Una Oportunidad Para América Latina y El Caribe; CEPAL: Santigo, Chile, 2018.
  6. Pino, V.E.; Chávarri, V.E.; Ramos, F.L. Crisis de Gobernanza y Gobernabilidad y Sus Implicancias En El Uso Inadecuado Del Agua Subterránea, Caso Acuífero Costero de La Yarada, Tacna, Perú. Idesia 2018, 36, 77–85. [Google Scholar] [CrossRef] [Green Version]
  7. Post, V.E.A. Fresh and Saline Groundwater Interaction in Coastal Aquifers: Is Our Technology Ready for the Problems Ahead? Hydrogeol. J. 2005, 13, 120–123. [Google Scholar] [CrossRef]
  8. Mahlknecht, J.; Merchán, D.; Rosner, M.; Meixner, A.; Ledesma-Ruiz, R. Assessing Seawater Intrusion in an Arid Coastal Aquifer under High Anthropogenic Influence Using Major Constituents, Sr and B Isotopes in Groundwater. Sci. Total Environ. 2017, 587–588, 282–295. [Google Scholar] [CrossRef] [PubMed]
  9. Cheng, A.H.-D.; Ouazar, D.; Bear, J.; Sorek, S.; Herrera, I. (Eds.) Seawater Intrusion in Coastal Aquifers-Concepts, Methods and Practices; Springer Science & Business Media: Dordrecht, The Netherlands, 1999; Volume XXXIII. [Google Scholar]
  10. Chai, Y.; Xiao, C.; Li, M.; Liang, X. Hydrogeochemical Characteristics and Groundwater Quality Evaluation Based on Multivariate Statistical Analysis. Water 2020, 12, 2792. [Google Scholar] [CrossRef]
  11. Güler, C.; Thyne, G.D.; McCray, J.E.; Turner, K.A. Evaluation of Graphical and Multivariate Statistical Methods for Classification of Water Chemistry Data. Hydrogeol. J. 2002, 10, 455–474. [Google Scholar] [CrossRef]
  12. Menció, A.; Mas-Pla, J. Assessment by Multivariate Analysis of Groundwater-Surface Water Interactions in Urbanized Mediterranean Streams. J. Hydrol. 2008, 352, 355–366. [Google Scholar] [CrossRef]
  13. Belkhiri, L.; Boudoukha, A.; Mouni, L. A Multivariate Statistical Analysis of Groundwater Chemistry Data. Int. J. Environ. Res. 2011, 5, 537–544. [Google Scholar]
  14. Villegas, P.; Paredes, V.; Betancur, T.; Ribeiro, L. Assessing the Hydrochemistry of the Urabá Aquifer, Colombia by Principal Component Analysis. J. Geochem. Explor. 2013, 123, 120–129. [Google Scholar] [CrossRef]
  15. Blake, S.; Henry, T.; Murray, J.; Flood, R.; Muller, M.R.; Jones, A.G.; Rath, V. Compositional Multivariate Statistical Analysis of Thermal Groundwater Provenance: A Hydrogeochemical Case Study from Ireland. Appl. Geochem. 2016, 75, 171–188. [Google Scholar] [CrossRef] [Green Version]
  16. Werner, A.D.; Simmons, C.T. Impact of Sea-Level Rise on Sea Water Intrusion in Coastal Aquifers. Ground Water 2009, 47, 197–204. [Google Scholar] [CrossRef]
  17. Binda, G.; Pozzi, A.; Livio, F. An Integrated Interdisciplinary Approach to Evaluate Potentially Toxic Element Sources in a Mountainous Watershed. Environ. Geochem. Health 2020, 42, 1255–1272. [Google Scholar] [CrossRef] [PubMed]
  18. Gabrielli, P.; Cozzi, G.; Torcini, S.; Cescon, P.; Barbante, C. Trace Elements in Winter Snow of the Dolomites (Italy): A Statistical Study of Natural and Anthropogenic Contributions. Chemosphere 2008, 72, 1504–1509. [Google Scholar] [CrossRef] [PubMed]
  19. Ritter, B.; Wennrich, V.; Medialdea, A.; Brill, D.; King, G.; Schneiderwind, S.; Niemann, K.; Fernández-Galego, E.; Diederich, J.; Rolf, C.; et al. Climatic Fluctuations in the Hyperarid Core of the Atacama Desert during the Past 215 Ka. Sci. Rep. 2019, 9, 1–13. [Google Scholar] [CrossRef] [Green Version]
  20. Pino, V.E.; Tacora, P.; Steenken, A.; Alfaro, L.; Valle, A.; Chavarri, E.; Ascencios Templo, D.R.; Mejía, J. Efecto de Las Características Ambientales y Geológicas Sobre La Calidad Del Agua En La Cuenca Del Río Caplina, Tacna, Perú. Tecnología y Ciencias del Agua 2017, 8, 77–99. [Google Scholar] [CrossRef]
  21. Pino, V.E.; Ramos, F.L.; Mejía, M.J.; Chávarri, V.E.; Ascensios, T.D. Medidas de Mitigación Para El Acuífero Costero La Yarada, Un Sistema Sobreexplotado En Zonas Áridas. Idesia 2020, 38, 21–31. [Google Scholar] [CrossRef]
  22. Ahmed, I.; Umar, R. Groundwater Flow Modelling of Yamuna-Krishni Interstream, a Part of Central Ganga Plain Uttar Pradesh. J. Earth Syst. Sci. 2009, 118, 507–523. [Google Scholar] [CrossRef]
  23. Pino, E.; Ramos, L.; Avalos, O.; Tacora, P.; Chávarri, E.; Angulo, O.; Ascensios, D.; Mejía, J. Factors Affecting Depletion and Pollution by Marine Intrusion in the La Yarada’s Coastal Aquifer, Tacna, Peru. Tecnol. Y Cienc. Del Agua 2019, 10, 177–213. [Google Scholar] [CrossRef]
  24. ANA. Autoridad Nacional del Agua Caracterización Hidrogeológica Del Acuífero Del Valle Del Caplina-La Yarada; ANA: Lima, Perú, 2009. [Google Scholar]
  25. ANA. Autoridad Nacional del Agua Caracterización Hidrogeoquímica Del Acuífero Caplina; ANA: Lima, Perú, 2010. [Google Scholar]
  26. Narvaez-Montoya, C.; Torres-Martínez, J.A.; Pino-Vargas, E.; Cabrera-Olivera, F.; Loge, F.J.; Mahlknecht, J. Predicting Adverse Scenarios for a Transboundary Coastal Aquifer System in the Atacama Desert (Peru/Chile). Sci. Total. Environ. 2022, 806, 150386. [Google Scholar] [CrossRef]
  27. Aranyossy, J. Estudio Isótoico Del Aquifero La Yarada (Perú). Informe de Interpretación; Proyecto PER/8/005; PET-INADE: Dakara, Senegal, 1988. [Google Scholar]
  28. Valenzuela, L.; Ramírez-Hernández, J.; Palomares, R.B. Composición Isotópica Del Agua Subterránea y Su Relación Con La Salinidad En El Valle de San Luis Río Colorado En Sonora, México. Inf. Tecnológica 2013, 24, 57–66. [Google Scholar] [CrossRef]
  29. Cruz Falcón, A.; Troyo Diéguez, E.; Murillo Jiménez, J.M.; García Hernández, J.L.; Murillo Amador, B. Groundwater Types and Total Dissolved Solids Distribution in the Aquifer of La Paz Baja California Sur, México. Terra Latinoam. 2018, 36, 31–38. [Google Scholar] [CrossRef] [Green Version]
  30. Henley, R.W.; Truesdell, A.H.; Barton, P.B. Fluid-Mineral Equilibria in Hydrothermal Systems; Society of Economic Geologists, Inc.: Littleton, CO, USA, 1984. [Google Scholar]
  31. Pool, M.; Carrera, J. Dynamics of Negative Hydraulic Barriers to Prevent Seawater Intrusion. Hydrogeol. J. 2010, 18, 95–105. [Google Scholar] [CrossRef]
  32. Jorreto, S.; Sola, F.; Vallejos, Á.; Sánchez-Martos, F.; Gisbert, J.; Molina, L.; Rigol, J.P.; Pulido-Bosch, A. Evolución de La Geometría de La Interfase Agua Dulce-Agua de Mar En Un Acuífero Costero Sometido a Un Intenso Bombeo de Agua de Mar. Geogaceta 2017, 62, 87–90. [Google Scholar]
  33. Pino, E.; Chávarri, E. Evidence of Climate Change in the Hyper-Arid Region of the Southern Coast of Peru, Head of the Atacama Desert. Tecnol. Cienc. Agua 2022, 12. in press. [Google Scholar]
  34. Craig, H. Isotopic Variations in Meteoric Waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef]
  35. Arslan, H. Application of Multivariate Statistical Techniques in the Assessment of Groundwater Quality in Seawater Intrusion Area in Bafra Plain, Turkey. Environ. Monit. Assess. 2013, 185, 2439–2452. [Google Scholar] [CrossRef]
  36. Moujabber, M.E.; Samra, B.B.; Darwish, T.; Atallah, T. Comparison of Different Indicators for Groundwater Contamination by Seawater Intrusion on the Lebanese Coast. Water Resour. Manag. 2006, 20, 161–180. [Google Scholar] [CrossRef]
  37. Liu, C.W.; Lin, K.H.; Kuo, Y.M. Application of Factor Analysis in the Assessment of Groundwater Quality in a Blackfoot Disease Area in Taiwan. Sci. Total Environ. 2003, 313, 77–89. [Google Scholar] [CrossRef]
  38. Appelo, C.A.J.; Postma, D. Groundwater Geochemistry. In A Practical Guide to Modelling of Natural and Contaminated Aquatic Systems; Cambridge University Press: Cambridge, UK, 1993; Volume 536, ISBN 978-84-9732-487-8. [Google Scholar]
  39. Acosta, H.; Alván, A.; Cutipa, M.; Mamani, M.; Rodríguez, J. Geología de los Cuadrángulos de La Yarada, Tacna y Huaylillas, Hojas 37-u, 37-v y 37-x, Escala 1: 50,000–[Boletín A 145]; INGEMMET: Lima, Peru, 2012. [Google Scholar]
  40. Acosta, H.; Alván, A.; Cutipa, M.; Mamani, M.; Rodríguez, J. Geología, de los Cuadrángulos de Pachía y Palca, hojas 36-v y 36-x, Escala 1: 50,000; Serie A: C.; INGEMMET: Lima, Peru, 2011. [Google Scholar]
  41. Robinson, W.O.; Edginton, G. Fluorine in soils. Soil Sci. 1946, 61, 341–354. [Google Scholar] [CrossRef]
  42. Arif, M.; Hussain, I.; Hussain, J.; Sharma, M.K.; Kumar, S.; Bhati, G. GIS-Based Inverse Distance Weighting Spatial Interpolation Technique for Fluoride Distribution in South West Part of Nagaur District, Rajasthan. Cogent Environ. Sci. 2015, 1, 1038944. Available online: http://www.editorialmanager.com/cogentenv (accessed on 4 October 2021). [CrossRef]
  43. Motevalli, A.; Moradi, H.R.; Javadi, S. A Comprehensive Evaluation of Groundwater Vulnerability to Saltwater Up-Coning and Sea Water Intrusion in a Coastal Aquifer (Case Study: Ghaemshahr-Juybar Aquifer). J. Hydrol. 2018, 557, 753–773. [Google Scholar] [CrossRef]
Figure 1. Location of the study area. (a) Atacama Desert in Chile and Peru, (b) Barroso Mountain range, a separation boundary between Caplina and Maure basins, (c) La Yarada coastal aquifer in the Caplina basin.
Figure 1. Location of the study area. (a) Atacama Desert in Chile and Peru, (b) Barroso Mountain range, a separation boundary between Caplina and Maure basins, (c) La Yarada coastal aquifer in the Caplina basin.
Water 13 03161 g001
Figure 2. (a) Conceptual model of the negative effect on the hydraulic barriers, modified after [31] and (b) geometry of the freshwater–seawater contact, predicted from the mixing model [32].
Figure 2. (a) Conceptual model of the negative effect on the hydraulic barriers, modified after [31] and (b) geometry of the freshwater–seawater contact, predicted from the mixing model [32].
Water 13 03161 g002
Figure 3. Iso-value plot in the groundwater of the La Yarada aquifer (this study) for 1988: (a) Na+, (b) Cl−, and for 2020: (c) Na+, (d) Cl. The high concentration zones of Na+ and Cl (i.e., seawater intrusion) have moved from W to SW along the sea-land boundary during the period.
Figure 3. Iso-value plot in the groundwater of the La Yarada aquifer (this study) for 1988: (a) Na+, (b) Cl−, and for 2020: (c) Na+, (d) Cl. The high concentration zones of Na+ and Cl (i.e., seawater intrusion) have moved from W to SW along the sea-land boundary during the period.
Water 13 03161 g003
Figure 4. Distribution of δ18O and δ2H values of samples collected at the study area with the global meteoric water line (GMWL). The local meteoric water lines (LMWL) for 1988 and 2020 were plotted with respective equations. The isotopic values of Caplina and Maure Basin are on the trend line.
Figure 4. Distribution of δ18O and δ2H values of samples collected at the study area with the global meteoric water line (GMWL). The local meteoric water lines (LMWL) for 1988 and 2020 were plotted with respective equations. The isotopic values of Caplina and Maure Basin are on the trend line.
Water 13 03161 g004
Figure 5. (a) Spatial distribution of δ18O in the Caplina, Sama, and Maure hydrological basins in 1988. (b) Update of the isotopic composition in the aquifer for the year 2020.
Figure 5. (a) Spatial distribution of δ18O in the Caplina, Sama, and Maure hydrological basins in 1988. (b) Update of the isotopic composition in the aquifer for the year 2020.
Water 13 03161 g005
Figure 6. Dendrogram derived using 14 parameters of 53 groundwater samples the log-transformed and standardized (z-score) water chemistry and isotopic ratios dataset for 2020 sampling campaign. A dashed horizontal line is used to define four groups of water samples.
Figure 6. Dendrogram derived using 14 parameters of 53 groundwater samples the log-transformed and standardized (z-score) water chemistry and isotopic ratios dataset for 2020 sampling campaign. A dashed horizontal line is used to define four groups of water samples.
Water 13 03161 g006
Figure 7. (a) Spatial distribution of four clusters, obtained from the grouping of 53 groundwaters sampled in 2020. (b) Piper diagram for samples from the years 2020 and 1988.
Figure 7. (a) Spatial distribution of four clusters, obtained from the grouping of 53 groundwaters sampled in 2020. (b) Piper diagram for samples from the years 2020 and 1988.
Water 13 03161 g007
Figure 8. Spatial distribution of factors: (a) FA1 (seawater chemical front) and (b) FA2 (seawater isotopic front). The clusters were drawn from Figure 7 to compare the different multivariate approaches.
Figure 8. Spatial distribution of factors: (a) FA1 (seawater chemical front) and (b) FA2 (seawater isotopic front). The clusters were drawn from Figure 7 to compare the different multivariate approaches.
Water 13 03161 g008
Figure 9. (a) Logarithmic correlation of the proportion factor (γ) versus littoral distance along with anomalous distribution configuration. (b) 3D model of the interface (γ = 0.20–0.30) represented in percentages. (c) Map of exploitation volume in the La Yarada aquifer and piezometric levels for 2018.
Figure 9. (a) Logarithmic correlation of the proportion factor (γ) versus littoral distance along with anomalous distribution configuration. (b) 3D model of the interface (γ = 0.20–0.30) represented in percentages. (c) Map of exploitation volume in the La Yarada aquifer and piezometric levels for 2018.
Water 13 03161 g009
Table 1. Data statistical of 14 physicochemical and isotopic parameters in 2020, sampled in the dry (August 2020) and wet (November 2020) seasons.
Table 1. Data statistical of 14 physicochemical and isotopic parameters in 2020, sampled in the dry (August 2020) and wet (November 2020) seasons.
StatisticsCEpHTDSCa2+Mg2+Na+K+HCO3SO42−ClNO3F δ 18 O δ 2 H
μS/cm mg/L
Dry Season (August 2020)
Mean1772.87.61248.6164.934.3161.817.487.8410.3362.13.60.2−12.1−88.0
Minimum25.55.863.510.70.847.91.553.445.465.50.020.1−13.5−96.2
Maximum8190.08.75594.0874.0223.0872.051.7196.01030.02417.020.40.6−6.7−49.9
Wet Season (November 2020)
Mean2075.17.61381.1185.635.2166.515.186.6394.7369.94.10.2−12.2−86.4
Minimum653.07.0425.012.10.951.82.846.242.870.20.60.1−13.8−95.8
Maximum8290.08.35918.0899.6209.1871.044.3170.0969.02266.024.10.6−4.7−31.2
Table 2. Correlation coefficient table for 14 parameters of water sample in 2020. EC stands for electrical conductivity (μS/cm), TDS is total dissolved solids (mg/L). The significant correlation values are marked as bold.
Table 2. Correlation coefficient table for 14 parameters of water sample in 2020. EC stands for electrical conductivity (μS/cm), TDS is total dissolved solids (mg/L). The significant correlation values are marked as bold.
ECpHTDSCa2+Mg2+Na+K+HCO3SO42−ClNO3F δ 18 O δ 2 H
EC1.00−0.400.700.780.730.450.650.250.500.660.350.310.100.12
pH−0.401.00−0.52−0.62−0.55−0.29−0.30−0.47−0.52−0.47−0.510.220.460.48
TDS0.70−0.521.000.920.930.860.840.190.710.980.430.270.030.08
Ca2+0.78−0.620.921.000.970.600.760.240.720.880.480.14−0.11−0.09
Mg2+0.73−0.550.930.971.000.660.820.140.770.890.330.11−0.13−0.07
Na+0.45−0.290.860.600.661.000.780.160.520.870.300.430.240.32
K+0.65−0.300.840.760.820.781.000.060.700.800.230.280.100.21
HCO30.25−0.470.190.240.140.160.061.000.280.130.530.03−0.06−0.13
SO42−0.50−0.520.710.720.770.520.700.281.000.580.27−0.16−0.40−0.31
Cl0.66−0.470.980.880.890.870.800.130.581.000.400.350.120.17
NO30.35−0.510.430.480.330.300.230.530.270.401.000.15−0.02−0.11
F0.310.220.270.140.110.430.280.03−0.160.350.151.000.740.75
δ 18 O 0.100.460.03−0.11−0.130.240.10−0.06−0.400.12−0.020.741.000.96
δ 2 H 0.120.480.08−0.09−0.070.320.21−0.13−0.310.17−0.110.750.961.00
Table 3. Relative ion abundance of the four clusters.
Table 3. Relative ion abundance of the four clusters.
ClusterRelative Abundance
1Cl>> SO42−>> Ca2+> Na+>> Mg2+> HCO3> K+> NO3
2Cl>> Na+> SO42−>> Ca2+> HCO3>> K+> Mg2+> NO3
3SO42−>> Cl> Ca2+> Na+> HCO3> Mg2+> K+> NO3
4SO42−Cl>> Ca2+> Na+> HCO3>> Mg2+> K+> NO3
Table 4. Values of total variance for the parameters in each factor. Var stands for variance and Acum for accumulative percentage variance.
Table 4. Values of total variance for the parameters in each factor. Var stands for variance and Acum for accumulative percentage variance.
FA1FA2FA3Communality
EC0.6760.1240.2820.551
pH−0.4310.423−0.5910.715
TDS0.9690.0790.2190.993
Ca2+0.895−0.0850.3190.910
Mg2+0.963−0.1130.1340.957
Na+0.7660.3010.1160.690
K+0.8960.135−0.0130.821
HCO30.064−0.0390.7020.499
SO42−0.736−0.3580.1650.696
Cl0.9250.1760.1870.922
NO30.2600.0280.7250.594
F0.2040.7940.0840.678
δ18O−0.0350.971−0.0440.946
δ2H0.0620.964−0.1780.964
Var6.2132.9911.733
%Var0.4440.2140.124
Accum0.4440.6570.781
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vera, A.; Pino-Vargas, E.; Verma, M.P.; Chucuya, S.; Chávarri, E.; Canales, M.; Torres-Martínez, J.A.; Mora, A.; Mahlknecht, J. Hydrodynamics, Hydrochemistry, and Stable Isotope Geochemistry to Assess Temporal Behavior of Seawater Intrusion in the La Yarada Aquifer in the Vicinity of Atacama Desert, Tacna, Peru. Water 2021, 13, 3161. https://doi.org/10.3390/w13223161

AMA Style

Vera A, Pino-Vargas E, Verma MP, Chucuya S, Chávarri E, Canales M, Torres-Martínez JA, Mora A, Mahlknecht J. Hydrodynamics, Hydrochemistry, and Stable Isotope Geochemistry to Assess Temporal Behavior of Seawater Intrusion in the La Yarada Aquifer in the Vicinity of Atacama Desert, Tacna, Peru. Water. 2021; 13(22):3161. https://doi.org/10.3390/w13223161

Chicago/Turabian Style

Vera, Alissa, Edwin Pino-Vargas, Mahendra P. Verma, Samuel Chucuya, Eduardo Chávarri, Miguel Canales, Juan Antonio Torres-Martínez, Abrahan Mora, and Jürgen Mahlknecht. 2021. "Hydrodynamics, Hydrochemistry, and Stable Isotope Geochemistry to Assess Temporal Behavior of Seawater Intrusion in the La Yarada Aquifer in the Vicinity of Atacama Desert, Tacna, Peru" Water 13, no. 22: 3161. https://doi.org/10.3390/w13223161

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop