Hydrochemical, Isotopic, and Geophysical Studies Applied to the Evaluation of Groundwater Salinization Processes in Quaternary Beach Ridges in a Semiarid Coastal Area of Northern Patagonia, Argentina

: Quaternary sea level fluctuations have led to the development of beach ridges on many South Atlantic coasts. The objective of this paper was to asses from lithological, hydrochemical, isotopic, and geophysical studies the salinization processes affecting groundwater stored in Pleistocene and Holocene beach ridges of the northern Patagonian coast. A hydrogeomorphological characterization of the area was performed using digital elevation models, the interpretation of satellite images, and field studies. Vertical electrical soundings were performed on transects running perpendicular to beach ridges in order to define variations in the freshwater ‐ saltwater interface position. The salinity, chemistry, and stable isotopes of the groundwater were analyzed. The results demonstrated that the groundwater salinization of Pleistocene ridges responds to processes associated with the geological ‐ geomorphological evolution of the area. The cementation of these surface sediments limits rainwater infiltration, which consequently prevents the development of freshwater lenses. This suggests that saline water is the result of ancient marine ingressions. Freshwater lenses develop in Holocene beach ridges; however, slight water salinization is detectable in the most populated areas as a result of intensive exploitation. The data provided are useful for freshwater resource prospection along the arid coast of Patagonia, where beach ridge deposits abound and populations experience serious drinking water supply problems.


Introduction
Quaternary sea level fluctuations that occur as a result of climate change have led to the formation of coastal plains, beach ridges, spits, and sandy bars in coastal areas in many parts of the world [1,2]. In coastal plains, groundwater is predominantly saline due to the marine ingressions that originated them and/or the current intrusion of seawater [3]. The low permeability of sediments in coastal plains prevents rainwater from infiltrating and washing salts [4,5]. By contrast, beach ridges and spits generally associated with dune fields constitute positive relief morphologies of high permeability where rainwater continuously infiltrates [6][7][8][9]. This preferential and localized recharge contributes to the development of a freshwater lens that displaces the surrounding saltwater downward by a density difference [10].
This determines that, in many coastal areas of the world, the only freshwater available sources are in the shape of lenses related to beach ridges and dune deposits [11][12][13][14][15][16]. These freshwater resources are usually relatively small, but significant for the development of coastal populations [17][18][19][20][21].
Freshwater lens salinization is a global issue that deteriorates water supply sources, therefore limiting the population development of coastal areas. This salinization may be due to seawater intrusion due to the intensive exploitation of freshwater lenses [22], by upward flowing of saline groundwater [23,24], or by past marine transgressions [25][26][27][28]. In arid and semi-arid areas, salinization problems are accentuated due to scarce and erratic rainfall and the high rate of evaporation that condition aquifers recharge [29][30][31]. To this end, understanding freshwater-saltwater distribution and the processes leading to water salinization is essential for water resource management in coastal areas [32,33]. Defining the origin is fundamental for groundwater management, since the problem is different if it is marine paleowaters or if it is a seawater intrusion due to anthropic actions. If the increase in salts is due to seawater intrusion, this can be mitigated with sustainable exploitation strategies and with actions that maximize freshwater lens recharge. In this sense, hydrogeochemical, isotopic, and geophysical tools are of importance in the study of the origin and mechanisms of groundwater salinization e.g., [34][35][36][37][38][39][40][41].
The Patagonian coast of Argentina constitutes an extensive coastal area, in which arid and semi-arid conditions predominate and coastal populations face serious freshwater supply problems [42,43]. Particularly on the coast of Patagonia, the uplift of the area resulted in the preservation of beach ridges generated during different Quaternary ingressions [44][45][46][47]. In the northern sector of the coast of Patagonia, in Bahía San Blas (Figure 1), beach ridge deposits can be found, separated by marsh and intertidal plain environments associated with the Jabalí channel ( Figure 1). Beach ridges located in a more continental position are associated with an ancient Pleistocene shoreline (MIS 5e), while those located between the Jabalí channel and the sea correspond to a Holocene shoreline (MIS 1) that is currently active [44,48,49]. In this arid zone, where annual rainfall and potential evapotranspiration average 300 and 730 mm [50], respectively, groundwater stored in the shape of freshwater lenses constitutes the only source of supply for local villagers. Previous studies identified the occurrence of predominantly saline groundwater in the Pleistocene ridges and the presence of freshwater lenses in the Holocene (Carol et al., 2021). However, the samples extracted from the Pleistocene ridges were located towards their edges, raising the question of whether the observed behavior occurs in the entire area of Pleistocene beach ridges or only at their edges, adjacent to the intertidal environments. From this arises the need to increase the number of sampling points in the area of Pleistocene ridges, as well as to carry out detailed studies of the sediments that compose both beach ridges, and to define the variation in water salinity in depth. In this context, the aim of this study was to evaluate from hydrochemical, isotopic, and geophysical studies the salinization processes affecting groundwater stored in Pleistocene and Holocene beach ridges of the Bahía San Blas on the northern Patagonian coast.

Materials and Methods
The geomorphological characteristics of the area were described by means of digital elevation models (DEM) performed by the National Geographic Institute of Argentina (IGN; MDE-Ar v2.0 30 m), the interpretation of satellite images (Google Earth software), and field surveys. The DEM MDE-Ar v2.0 offers a spatial resolution of 30 m and a vertical accuracy of approximately 2 m [51]. In addition, the data acquired from these DEMs were also cross-checked with the altimetric information collected in the field with a Trimble ® R2 differential GPS consisting of an integrated unit using RTX real-time positioning mode with a horizontal accuracy of 4 cm and vertical accuracy of 9 cm (40 points, coinciding with the sampling sites).
During the field surveys, soil trenches and drillings with manual auger (up to 2 m deep) were performed, and lithological profiles were described in quarries in order to determine the hydrolithological characteristics of the superficial sediments that condition rainwater infiltration.
Vertical electrical soundings (VES) were carried out on transects running perpendicular to the beach ridges to define variations in the freshwater-saltwater interface position. VES data were acquired using a Schlumberger array with a maximum electrode distance of 100 m. During field work, an apparent resistivity curve was obtained after a direct current circulated through the emission circuit and the potential difference generated between the receiver electrodes was measured. The inversion scheme followed here was proposed by [52]; the number of layers was then reduced using the Dar Zarrouk parameters [53] and the response curve was calculated with a linear filter [54]. Additional information was added to adjust the model, especially the thickness or depth of certain layers (manual auger), improving the resistivity determination.
Groundwater samples were obtained from shallow wells (4-8 m deep) and the electrical conductivity (EC), temperature, and pH of the water were measured in situ using a portable multiparameter (Lutron ® WA-2017SD, Lutron Electronics, Coopersburg, PA, USA). The water samples were obtained with polyethylene bottles and the analysis of the major ions was performed in the laboratory at the Centro de Investigaciones Geológicas (Geological Research Center) using standardized methods [55]. The carbonate (CO3 2− ), bicarbonate (HCO3 − ), chloride (Cl − ), calcium (Ca 2+ ), and magnesium (Mg 2+ ) were determined by titration; the sodium (Na + ) and potassium (K + ) by flame photometry; and the sulfate (SO4 2− ) by UV-Visible spectrophotometry. The determination of the saturation indexes (SI) was conducted using the PHREEQC software [56].
The analysis of stable isotopes in water (δ 2 H and δ 18 O) was performed in some samples at the Universidad Nacional de San Luis (National University of San Luis), following the Cavity Ring-Down Spectroscopy (CRDS) methods, using a Picarro L2120-i coupled to a high-precision vaporizer A0211. The isotopic results are expressed as δ‰, defined as δ = 1000 (Rs − Rr)/Rr‰, where δ is the isotopic deviation in ‰ relative to Vienna Standard Mean Ocean Water (V-SMOW) [57]; s: sample; r: International reference; R: isotopic ratio ( 2 H/ 1 H, 18 O/ 16 O). The analytical accuracy is ± 0.05‰ and ± 0.5‰, for δ 18 O and δ 2 H, respectively.
The freshening and salinization processes were analyzed using a hydrochemical facies evolution diagram, HFE-D [58], ionic and isotopic relationships graphs, and water mixing calculations. To evaluate the freshwater lens recharge from rainwater, isotopic values were compared with the local meteoric line δ 2 H = 7.73 ± 0.28 * δ 18 O + 8.89 ± 1.67 [59], these being the closest data for a coastal area. Theoretical water mixtures were obtained considering as end members the average chemical and isotopic composition of fresher water samples (Cfw) and the seawater sampled from the Jabalí channel (Csw).
where y is the percentage of seawater and (1 − y) the percentage of freshwater in the lenses.

Geomorphological and Hydrolithological Characterization
Field work supported by DEM and satellite image analysis enabled the geomorphological characterization of the Pleistocene and Holocene beach ridge deposits. Pleistocene beach ridges comprise a series of linear deposits with a prevailing N-S orientation separated by more depressed areas. These ridges reach heights of about 12 m a.s.l., while in depressed areas, the heights vary between 0.5 and 4 m a.s.l. Depressed areas could have constituted ancient coastal lagoons or ancient intertidal channels ( Figure 2a).
The more superficial sediments in these beach ridges comprise high-permeability gravels and sands (Figure 2b), observing from a depth of 50 cm levels of gravel cemented by carbonates (Figure 2c,d). In addition, in mills' large-diameter excavations, the presence of a cemented sandstone level was noticed at a 3 m depth, which constituted the base up to where the supply perforations were made.
Holocene beach ridges originated as a beach ridge plain developed in the N-S direction that diffracted in the N sector towards the NW due to the coastal drift that currently originates and models the current deposits of the barrier spit. The highest heights of Holocene beach ridges (8 m a.s.l.) are found in the central sector, decreasing both towards the north end of the spit, and the sea and Jabalí channel ( Figure 2a).
Within Holocene ridges, spatial variations in the hydrolithological characteristics of the superficial sediments are observed; all are associated with high-permeability lithologies. In the N sector at the tip of the spit where heights do not exceed 5 m a.s.l., the superficial sediments are composed of gravels with a sandy matrix (Figure 2e) passing to predominantly gravel deposits on the edges of the spit, limiting both the sea and the Jabalí channel. In the central part of the spit, sediments are gravelly, although the percentage of sand in the matrix increases; levels of interbedded sand in the gravel deposits are also observed ( Figure 2f). In the S sector, sandy deposits of coastal dunes develop on the surface (Figure 2g), mainly in areas close to the sea. On the other hand, towards the Jabalí channel, dune deposits are scarce and the superficial sediments comprise gravel with a sandy matrix.

Geophysical Data
The recorded curves demonstrate a great contrast in the apparent resistivity values, with differences of almost two orders of magnitude (1000 to 10 Ohm m). The general behavior of the curves starts with high values, which then follow a decreasing trend in the measured values. The proposed models feature an adjustment error (root mean square) of less than 5% between the observed and calculated curves. In the proposed resistivity models, it is recognized that the first layers are of high resistivity; later, a sequence of intermediate resistivity layers supported on a very low resistivity base can be observed.
The geophysical models proposed for Pleistocene ridges (VES P1, P2, and P3) demonstrate that the unsaturated zone presents resistivity values of 6 to 15 Ohm m ( Figure 3). In the VES P1, a sequence of layers of higher resistivity (200 and 250 Ohm m) is observed, probably due to a lithological change. The phreatic aquifer in this sector is characterized by layers of very low resistivity (less than 5 Ohm m) associated with the presence of brackish water. In the integrated geophysical models on the Holocene ridges, the unsaturated zone demonstrated higher resistivity values (more than 100 Ohm m) in the topographically higher sector (VES H3 and H4), while in the lowest one, it displayed low resistivity (8 and 9 Ohm m) ( Figure 4). The phreatic aquifer presented an intermediate resistivity (60 Ohm m) in the VES H2 and H3, attributable to freshwater and lower values (less than 5 Ohm m) in the VES H1 and H4, possibly caused by an increase in EC of the water (higher salinity) (Figure 4).

Hydrochemistry and Stable Isotopes
The analysis of EC values and major ions contents demonstrated sharp differences between the groundwater samples from the Pleistocene and Holocene beach ridges. The Pleistocene beach ridges were distinguished by high saline groundwater (with EC values between 5160 and 16,860 μS/cm), while in the Holocene ridges, groundwater salinity was lower in most samples (EC less than 2000 μS/cm). EC values between 2000 and 5150 μS/cm were recorded in the wells located in the most populated sector and the three wells in the S sector. On the other hand, only one groundwater sample from the Holocene ridges registered high salinity values (EC of 7960 μS/cm); it was located close to the tip of the spit ( Figure 5). The hydrochemical facies evolution diagram, HFE-D [58] also demonstrates a clear differentiation between groundwater samples from the Pleistocene and Holocene ridges ( Figure 6). In the Pleistocene beach ridges, groundwater is Na-Cl and is located along with seawater samples at the extreme associated with salinization processes of marine origin. Conversely, groundwater in the Holocene ridges moved from Na-Cl to Na-HCO3 facies; this variation was accompanied by a decrease in the EC (freshening trend). The salinization trends associated with seawater in Pleistocene ridges and those of freshening in the Holocene ridges were also reflected in the isotopic contents and major ions relations (Figure 7). Groundwater samples from Pleistocene ridges displayed isotopic values between −35.0 and −33.5‰ and between −5.2 and −4.6‰ of δ 2 H and δ 18 O, respectively. In the δ 2 H vs. δ 18 O plot (Figure 7a) they were located around a mixing line whose end member values corresponded to the average isotopic composition of the fresher water samples and the seawater sampled from the Jabalí channel. Within this theoretical mixing line, the isotopic values of the groundwater of the Pleistocene ridges were associated with mixing percentages, with seawater contributions of up to 20%. This trend towards salinization associated with seawater inputs was also observed in both δ 18 O vs. Cl − and Na + vs. Cl − ratios (Figure 7b,c). Note that in the Na + vs. Cl − ratio, when plotting other samples of Pleistocene ridges that did not featured isotopic values but instead featured major ion determinations, seawater contributions in the theoretical mixtures reached values of up to 32%. Conversely, the Na + -Cl − vs. (HCO3 − + SO4 2− ) − (Ca 2+ + Mg 2+ ) ratio demonstrated that samples were located around the 1:1 line, indicative of cation exchange processes, with the largest part of the samples located in the release Ca-fixation Na field (Figure 7d). Additionally, the Ca 2+ /Na + ratio in the samples from the Pleistocene ridges was low and similar to that of seawater; this trend was maintained with the increase in Cl − (Figure 7e).  (Figure 7a), they are located close to the local meteoric line [59], indicating recharge from rainwater infiltration. The scarce isotopic variation and the low Cl − content of these samples suggests that the samples with the least saline were located at the extreme representation of freshwater in the graph δ 18 O vs. Cl − , while those that presented EC values between 2000 and 5150 μS/cm moved towards the mixing line with seawater contributions lower than 5%, reaching 11% in the higher-saline sample (Figure 7b). This tendency and mixing of percentages were in turn evidenced in the Na + vs. Cl − ratio ( Figure  7c).

Discussion
Beach ridges deposited during the quaternary in the Bahía San Blas area constitute topographically positive geoforms, composed of coarse-grained sediments (gravel and sand), where freshwater lenses could be stored. Despite this, the presence of carbonate cementation recorded in Pleistocene ridges and the absence of cement in the Holocene ridges represents a hydrolithologically relevant difference conditioning the development of freshwater lenses.
In the Pleistocene beach ridges, layers of gravel and cemented sands were recorded both in shallow perforations and in road outcrops. Their presence was also displayed in the geophysical studies (Figure 3), which characterized the unsaturated zone of these beach ridges by low resistivity values (values of 6 to 15 Ohm m). This type of carbonatic cement filling spaces between gravels and sands, causing a strong decrease in the permeability of the sediments, has been described in coastal sectors of northern Patagonia as well [60]. The decrease in permeability due to cementation limits rainwater infiltration, consequently hindering freshwater lens development. Geophysical and hydrochemical studies evidence the absence of freshwater in the entire Pleistocene beach ridge area. Geophysical studies (Figure 3) demonstrated the presence of low resistivity levels (less than 5 Ohm m) below the electro-layer associated with the unsaturated zone, which was related to high-salinity groundwater. In addition, groundwater EC measurements in this area were in all cases higher than 5150 μS/cm ( Figure 5). Regarding the origin of this salinity, chemical and isotopic signals would indicate that it is primarily marine. In the hydrochemical facies evolution diagram, HFE-D [58], it is demonstrated that in Pleistocene beach ridges, Na-Cl facies dominate, with clear salinization trends associated with seawater ( Figure 6). Similar hydrochemical trends of salinization by seawater mixing are seen in isotopic and ionic ratios (Figure 7a-c). In the isotopic ratio δ 2 H vs. δ 18 O (Figure  7a), the analyzed samples register seawater contributions of up to 20%; however, the presence of more saline waters than those with isotopic data and that present similar Na + vs. Cl − trends (Figure 7c) displaying, in some sectors, higher mixing percentages (up to 32%).
In contrast to the Pleistocene ridges, the Holocene beach ridges were composed of loose sands and gravels without cementation ( Figure 2). The unsaturated zone in these ridges presented higher resistivity values (more than 100 Ohm m), while towards the Jabalí channel (VES H1) in marsh edge environments with saline soils, the resistivity was low (8 and 9 Ohm m). The high permeability of the sediments constituting these beach ridges enables rapid rainwater infiltration and the development of freshwater lenses. The presence of fresh groundwater was verified through geophysical and hydrochemical studies. In the central sectors of beach ridges, VES demonstrated that under the electrolayer corresponding to the unsaturated zone, another layer, with high resistivity corresponding to a groundwater level of low salinity, develops ( Figure 4). Additionally, it was observed that this lens is very thin and is salinized both towards the west and east, where it is bordered by the marsh and the sea, respectively. Towards the Jabalí channel in the topographically low sectors close to the marsh, VES indicates the existence of lowresistivity electro-layers associated with the proximity of shallow and saline groundwater. On the other hand, towards the sectors next to the sea, below the electro-layer corresponding to the unsaturated zone, the presence of an electro-layer with resistivities close to 25 Ohm, characteristic of brackish water, was observed. Thereby, freshwater lenses would develop in the central sectors, where beach ridge deposits present higher topographic heights. Similarly, it is important to note that the freshwater thicknesses estimated from the VES were less than 5 m and that VES H3 located in the central part of the ridges presented a depression in the water table. Considering that this VES is located in the most populated sector, where there are a large number of household wells, this decrease in the water level would suggest intensive exploitation of the aquifer. EC measurements carried out in different sectors of the Holocene ridges mainly presented low-salinity water (less than 2000 μS/cm). In the populated area and some topographically low sectors located to the south, EC tends to increase, registering saline water only in a well located at the north end of the spit ( Figure 5). The groundwater isotopic ratio δ 2 H vs. δ 18 O (Figure 7a), with similar values to the meteoric line, indicates that rainwater infiltrates rapidly in these ridges, recharging the freshwater lens. The processes of freshening and the associated cation exchange are evidenced in the hydrochemical facies evolution diagram, HFE-D [51], and in the ratios of Na + − Cl − vs. (HCO3 − +SO4 2− ) − (Ca 2+ + Mg 2+ ) and Ca 2+ /Na + (Figures 6 and 7d). In the HFE-D diagram and both the isotopic and Na + vs. Cl − ratios, slight salinization trends are recorded in the wells with the highest EC, primarily located in the populated area (Figures 6 and 7a-c).
Albeit with different magnitudes, groundwater salinization occurs in both Pleistocene and Holocene beach ridges. However, as explained above, changes in sediment permeability condition the development of freshwater lenses. The joint analysis of the results allowed us to identify that there are different salinization processes. These salinization processes would be geogenic in Pleistocene ridges and mainly anthropic in the Holocene ridges.
Groundwater salinization in Pleistocene ridges is geogenic and would respond to processes associated with the quaternary geological-geomorphological evolution itself. Quaternary is characterized by sea-level oscillations that gave rise to the deposition of beach ridges. In Bahía San Blas, beach ridges have been identified whose deposits represent at least MIS 5e and MIS 1 transgressive episodes resulting from the cortical uplift in Patagonia [44,46,49]. In addition, there were marked climatic changes during this period, with evidence that after MIS 5e, beach ridges were exposed on the surface affected by geomorphic and pedogenetic processes, the latter identified by the presence of calcretes, which mainly developed in the upper part of the gravel deposits, in the boundary with the overlying sandy Holocene deposits [61,62]. The presence of pedogenetic carbonates would indicate arid and semi-arid conditions [63]. Among the factors that favor cementation by carbonate precipitates, there are those associated with lithological changes. Sandy soils occurrence, which is in many cases associated with aeolian deposits above the gravel deposits of beach ridges, favors the precipitation of carbonates. This is due to the infiltrating water being retained before reaching the pressure necessary to fill the larger pores of the underlying horizon [64,65]. At this granulometric boundary, the solution becomes concentrated as water is removed by evapotranspiration; therefore, salts and pedogenic carbonate accumulations form in/or near upper gravel deposits [66]. In this way, different cemented levels may exist in some areas (Figure 2c,d), depending on the diverse processes that have led to their cementation. These cemented layers would later function as impermeable or very low-permeability levels conditioning rainwater infiltration, consequently preventing freshwater lenses development. On the other hand, it is interpreted that salinization by seawater intrusion evidenced in the chemical and isotopic signal of water could not be current. This is because the current sea level could not generate a marine intrusion that influences the entire area of Pleistocene ridges, which are topographically higher than Holocene ridges. Due to this, it is to be expected that the marine intrusion affecting Pleistocene ridges would have occurred during MIS 1 along with the marine ingression that led to the deposition of Holocene ridges, specifically during the climatic optimum or Hypsithermal (Hypsithermal ca. 6000 yr. BP). In this ingression, although the sea level does not cover Pleistocene ridges, it could cause seawater intrusion in them. During the climatic optimum, only Pleistocene ridges above approximately 7 m a.s.l. would have remained above sea level. Sea would have entered not only from the east, but also from the west due to the flooding of what is currently the Walker channel, and could also flood the lowlands that cross Pleistocene ridges ( Figure 2). Later, with sea level decrease, regardless of the existence of wet periods, the presence of cemented layers from MIS 5a would limit rainwater infiltration and consequently would prevent the formation of freshwater lenses in these ridges that can displace saline water. For this reason, groundwater currently reflects a chemical and isotopic signal of mixing with seawater, where scarce rainfall hinders freshwater lens development as a separate phase from that of seawater. A scheme of the quaternary evolution of the area is shown in Figure 8. Changes in coastline position as a result of sea-level change cause an imbalance between the groundwater and seawater interface position and the present-day sea level [67][68][69]. Aquifer salinization by Quaternary marine ingressions with chemical signals of mixing with seawater has been recorded in other coastal aquifers in the region [70][71][72]. This also occurs in other areas worldwide, where groundwater and seawater interface positions are a result of late Quaternary paleoenvironmental conditions [15,[73][74][75]. In other regions, these salinization processes due to paleo-intrusions are compounded by the dissolution of evaporites formed during arid periods, which causes the salinity of groundwater to be higher than that of seawater. (e.g., Laizhou Bay in China [76,77]). By contrast, in other regions of the world, the deglaciation that led to the sea level rise during the Pleistocene favored the recharge of aquifers with fresh water from melting ice, causing the freshening of seawater by dilution (e.g., in the Ledo-Paniselian aquifer in Flanders, Belgium [78]). Likewise, areas of extreme aridity (Sahara desert) currently feature freshwater aquifers that were recharged during rainy periods in the Pleistocene. All these examples demonstrate how the Quaternary geological and climatic evolution of an area conditions groundwater chemistry, which may not respond to current water balances [79].
Conversely, Holocene ridges in their most recent geological-geomorphological evolution have not undergone periods of aridity capable of forming carbonate cement, or have not formed for a sufficient amount of time. This makes its gravel and sandy sediments maintain the high permeability that allows the rapid infiltration of rainwater. Despite this, scarce water excesses in the area [50] determine the development of thin freshwater lenses (Figure 8d). In these lenses, incipient salinization processes were recorded in the central part of the ridges where the San Blas population is concentrated. As a consequence, water level depression was measured (VES H3 in Figure 4) that not only decreases the thickness of the freshwater lens but also causes an increase in the EC (Figure 3). Therefore, the groundwater salinization of Holocene ridges is associated with anthropically induced processes. Intensive exploitation is the main cause of coastal aquifer salinization on a global scale e.g., [80][81][82][83]. The increased rate of abstraction due to higher water demand has generated significant challenges to sustainable fresh groundwater management, especially to the effective prevention of coastal aquifer salinization [22,27]. This situation is even more serious in coastal aquifers in arid and semiarid areas, where the freshwater lenses stored are thin due to scarce water excess, which makes them very vulnerable [31]. In addition to this problem, there is the associated regional climatic trend of increasingly arid conditions, generating a decrease in groundwater recharge. Faced with this situation, the need to propose management guidelines to control freshwater lens exploitation arises in order to achieve sustainable pumping that preserves the only freshwater resource in the region. Furthermore, in order to assess the vulnerability of groundwater resources, it is essential to possess sufficient information on salinization processes. In this sense, the data provided in this work constitute a scientific basis that makes it possible to recognize geogenically and anthropogenically salinized areas. It is in the latter where guidelines can be generated to mitigate groundwater salinization. To this end, community participation is essential for the success of long-term management strategies. Communities, as users and beneficiaries of groundwater, must understand the fragility of the resource and become aware of its use [79].

Conclusions
Bahía San Blas constitutes a representative area of the Atlantic coast of Patagonia where beach ridges were deposited during the Pleistocene and Holocene as a result of sealevel fluctuations. Currently, these ridges emerge, forming positive geoforms composed of coarse-grained sediments that can potentially store freshwater lenses.
Increasing the sampling points in the area of Pleistocene ridges, as well as carrying out detailed studies of the sediments that compose both beach ridges, and defining the variation in water salinity in depth using geophysical methods, provides novel data that make it possible to improve the comprehension of the processes that cause groundwater salinization.
The different methodologies applied in this study allowed us to interpret the different factors and mechanisms responsible for groundwater salinization in the different studied environments. Sedimentological surveys provided data to understand the evolution of the area and the hydrolithological characteristics of the beach ridges. Groundwater EC measurements, along with VES interpretation, allowed us to characterize salinized areas and freshwater lenses morphologically. Similarly, the VES demonstrated changes in the resistivity of the unsaturated zone, responding to the presence, or non-presence, of cementation. In addition, the isotopic and hydrochemical data enabled us to interpret the origin of salinization and the existence of areas where freshening dominates.
Pleistocene beach ridges store brackish groundwater, while Holocene ridges mainly store low-salinity groundwater, with a slight increase in the most populated areas.
The groundwater salinization of Pleistocene ridges responds to processes associated with the geological-geomorphological evolution of the area. Sea level oscillations followed by periods of aridity during the Pleistocene caused the precipitation of carbonates in the gravels composing the Pleistocene ridges, causing their cementation. This cementation restricted rainwater infiltration during the wettest periods and consequently prevented the development of freshwater lenses, thus favoring saline water permanence as a result of ancient marine ingressions. Studies of trace elements or radioactive isotopes will allow the assessment of the presence of paleo-seawater in greater detail in the future; however, both the geological-geomorphological evolution of the area and the data collected support our hypothesis.
On the other hand, the sands and gravels that compose Holocene ridges are not cemented. This allows the high permeability to sediments favoring rainwater infiltration and the development of freshwater lenses. Nevertheless, slight water salinization was detected in the most populated areas associated with water levels decreasing as a result of intensive exploitation.
In this way, it was shown how, in the same area, groundwater salinization can be due both to factors depending on Quaternary geological-geomorphological evolution and to anthropic action. The data provided in this paper are useful for prospecting freshwater resources along the arid Patagonian coast, where Quaternary beach ridge deposits abound and populations experience serious freshwater supply problems.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.