TDEM Soundings as a Tool to Determine Seasonal Variations of Groundwater Salinity (Villaf á ﬁla Lakes, Spain)

: Interaction between groundwaters with different salinities and lakes show seasonal variations driven by changes in precipitation and evapotranspiration. In the vicinity of Villaf á ﬁla lakes, local fresh and brackish regional groundwaters feeds the lakes, forming a brine in the lake sediments aquitard. Two TDEM surveys (summer 2019 and winter 2020) were carried out. Five TDEM soundings were acquired at the same location for each survey, forming a proﬁle from the hills to the lake-shore. Simultaneously to the TDEM surveys, electric conductivity of lake water and groundwater was measured. The resistivity boundary between the local fresh (10–35 Ohm/m) and regional brackish groundwater (2–5 Ohm/m) is well marked at 600 m above sea level (masl) below the hills, and at 650 masl below the lowlands surrounding the lakes. During the summer, fresh-brackish groundwater interphase rises due to evaporative pumping occurring in the lowlands. This increases groundwater salinity close to the terrain surface favoring precipitation of halite efﬂorescences. Annual record of EC in a piezometer conﬁrms the summer ascendant of the brine contained in the lake aquitard. TDEM sounding is fast and simple technique to monitor seasonal variations in fresh-brackish groundwater interphase and to detect possible salinization of consumption wells and environmental changes.


Introduction
Many lakes show interactions between their waters and groundwaters which sometimes have different salinities, densities, and compositions [1][2][3][4]. Time domain electromagnetism (TDEM) is a geophysical technique that is able to determine the lateral and vertical (depth) variations in electric resistivity (ρ) of the subsoil. Since the beginning of its application, TDEM has proven useful in research fields such as hydrogeology [5,6], mining [7], underground environmental pollution [8], and archaeology [9]. TDEM technique is useful to identify the interface of groundwaters with different salinities due to its high capacity to discriminate between units with high electric conductivity (EC) [10][11][12][13]. The dynamics of groundwaters and the lake-groundwater interaction can show seasonal and multi-annual changes due to the variations in evapotranspiration and precipitation along the year. Such variability can result in lateral and vertical changes in ground-and lake-water salinity and circulation [14]. Villafáfila lakes in Zamora province (Spain) are brackish to saline lakes located in a saline grassland. These lakes are protected by several environmental protections such as Natural Park, area of special protection for birds, and a Site of Community Interest (SCI) of the European Union, and are included in the intergovernmental RAMSAR treaty for protection of wetlands. These lakes are wintering place for migrant birds such as ducks and geese, and host one of world's largest populations of Great Bustard (Otis tarda) [15,16]. The aim of this study is to characterize, by means of TDEM surveys, the margins of the Salina Grande lake in Villafáfila, during both the dry and wet seasons and to compare simultaneously the ρ data with the EC of lake and groundwaters. This characterization is of interest for correlating the variations in soil and water salinity with faunal and plant changes. The location of the saline-freshwater interface will serve to the local and regional stakeholders to identify the areas and maximum depths to pump fresh groundwaters for human consumption.

Study Area
Villafáfila lakes are located in the northwestern part of the Duero river hydrological basin in northern Spain ( Figure 1A), on Cenozoic sedimentary rocks that constitute a great aquifer system, with groundwaters flowing from the periphery of the basin toward the Duero river that flows westwards, to the Atlantic Ocean, along the central axis of the basin [17,18].   Figure 3 is marked with an orange line, and the different sounding sites with red squares. Note that Culebra Fm. crops out in the southwestern part of the map.

Climate
Villafáfila has a semi-arid, Continental Mediterranean climate, characterized by cold and wet winters (November-February) and hot and dry summers (June-September). During the winter, average monthly temperatures are about 4 • C, but minimum temperatures are commonly below zero. During the summer, average monthly temperatures are about 22 • C [19]. Mean annual evaporation is 720 mm and real evapotranspiration is 237 mm. Mean annual precipitation is 330 mm, largely concentrated from October to May. During this wetter period, when precipitation exceeds evapotranspiration, soil surplus occurs; the aquifer is recharged and the lakes fill up. On the other hand, during the summer, when evapotranspiration exceeds precipitation, the soils suffer water deficit and the lakes dry out [20].

Hydrogeology
Villafáfila lakes are brackish to saline lakes that occur in the Salado river valley, a tributary of the Valderaduey river, which discharges downstream into the Duero river. Villafáfila lakes occur in a flat area at 677 m above sea level (masl) surrounded by gently sloping hills that reach up to 730 masl [21]. The Miocene rocks outcropping in the area (Facies Tierra de Campos and Facies Aspariegos Units; mainly siliceous/litharenite sandstones and siltstones) constitute the main local aquifer, that belongs to the large Duero basin aquifer system. The clayey-sandy Quaternary lake deposits behaves as an aquitard at the bottom of the lakes ( Figure 1B).
The water of the lakes comes mainly from the rain, but the ions dissolved are provided by shallow fresh groundwaters coming from the hills and regional deep brackish groundwaters. In Villafáfila, the regional component of the groundwater flow comes from the northeast ( Figure 1A) and shows an increase in salinity southwestwardly [17]. Deep brackish groundwaters raise toward the lakes bottom and the lowland areas due to the presence of a near surface Paleozoic metamorphic basement elevation. This elevation is perpendicular to the general flow and locates in the southern zone of the Salina Grande Lake ( Figure 1B) [20]. The basement is constituted by the Culebra Fm. and the Salamanca sandstone Fm. which does not crop out in the Villafáfila basin.
Groundwaters below the hills are of the calcium-bicarbonate type, with low salinities (TDS = 250 mg/L). Deep groundwaters, however, are sodium-chlorine waters with high salinities all year round (TDS = 3200-5000 mg/L). Lake waters have sodium-chlorine composition and show salinity variations along the year, with values ranging from 800-6200 mg/L. The brine contained in the Quaternary lake aquitard has salinities of 12,000-27,000 mg/L TDS [20].

Methodology
Two TDEM surveys were carried out, one during the winter (February 2020) and the other during the summer (July 2019) at the same five sounding sites. Soundings were placed aligned in a NNW-SSE trend, from higher in the hills to the lowlands around Salina Grande lake, close to the Villafáfila village ( Figure 1B). The position of the different soundings was chosen to evaluate the distribution of salinities from the upper part of the hills (700 masl), where local groundwaters recharge, to the lowland areas (677 m), where the regional groundwaters discharge and mixing processes occur [20,22]. The soundings were placed 200-600 m apart from each other, at the same positions for the summer and winter surveys, to evaluate the ρ variations along the year. The surveys were carried out by the company Técnicas Geofísicas S.L. with a TerraTEM (MONEX GeoScope Co., Ltd., Melbourne, Australia) transient electromagnetic survey system. Time domain electromagnetic sounding (TDEM) consists on injecting a constant electric current into a cable loop (transmitter coil). The transmitter loop generates an electromagnetic wave that propagates into the subsoil through the subsurface when the electric current is suddenly interrupted [23][24][25]. As the electromagnetic energy finds different materials underground, it induces eddy currents that generate secondary electromagnetic fields. These secondary electromagnetic fields are measured on the surface by a receiver loop or magnetic antenna and recorded as the induced energy diffuses into the earth. The resistivity of the subsurface materials depends on the rate of diffusion. The array consisted on 50 × 50 m square loops in the mode of coincident loops. Several 10 ms measurements were performed in every site, varying the number of stacks (up to 3000 stacks) to minimize electromagnetic noise. Similar arrangements have been performed for hydrogeological research [20,26,27]. Management of the obtained data and its conversion from apparent to real resistivities were made with IX1D V3 (INTERPEX Co., Ltd., Golden, CO, USA). A smooth model was generated for every TDEM site measured. Occam's inversion [28] was used to obtain the 30-40 layers smooth model ( Figure 2). The obtained resistivities were used to construct the NNW-SSE profiles ( Figure 3). Alongside TDEM survey acquisition, hydrogeological field measurements (groundwater level, temperature, and electric conductivity) were made in the Salina Grande lake, springs, boreholes and dug wells located closest to the TDEM sites (see next sections for more details). EC and temperature were measured with a Hanna HI 9835 m and probe. SGR-3 piezometer was monitored at 8 m depth with a conductivity, temperature, depth data logger (CTD) (OTT Hydromet Inc., Loveland, CO, USA). Local piezometry was measured during the summer of 2018 in dug wells screened at depths lower than 40 m. The data were interpolated with a convergent interpolation method and later manually redrawn to avoid inconsistencies due to the lack of data in certain areas. The piezometer SGR-3 is screened below the aquitard, and overflows at 0.4 m above lake bottom.     Figure 1. EC values measured in lake waters, dug wells and springs close to the profile path have been projected (Table 1)  and U-4, in order to facilitate the description. U-1 appears in the upper part of the hills (700-680 masl), has ρ values ranging from 25 to 35 Ohm/m, and extends toward the toe of the hills. U-2 appears in the lowlands surrounding Salina Grande and in the lake shore and shows low resistivity values (ρ < 5 Ohm/m) in the upper part of the terrain. U-3 locates generally at an intermediate position between U-1 and U-4, below the upper parts of the hills (SEDT-6BIS and 7BIS), and between U-2 and U-4 in the lake shore. For the winter survey U-3 locates very close to the terrain surface (SEDT-3BIS and 4BIS). U-3 has ρ ranging from 10 to 25 Ohm/m and appears at 680-600 masl, below the hills and at 685-640 masl at the lowlands. U-4 appears at the lower part of the profile and is present at all sites; below the hills it is between 600-500 masl (SEDT-7BIS) and below the lake shore between 640-500 masl (SEDT-4BIS,) and has ρ lower than 5 Ohm/m ( Figure 3A). The boundary between the U-3 and U-4 is almost horizontal between SEDT-6BIS and 7BIS (600 masl) and has a gentle slope between SEDT-6BIS and 3BIS. Below SEDT-3BIS and SEDT-4BIS is horizontal again.

Summer Geoelectric Profile
Resistivities measured during the summer survey have similar values (2-35 Ohm/m) to those recorded during the winter. The lowest ρ is registered in the lowlands surrounding Salina Grande and at the lake shore (SEDT-5, 4, 3) and the highest is recorded in the upper parts of the hills (SEDT-6 and 7), in the upper zone ( Figure 3B). As in winter, in the upper parts of the hills there is a decrease in ρ from the upper parts of the terrain where geoelectric U-1 occurs toward the lower parts where U-4 appears. The boundary between U-1 and U-3 deepens slightly from the upper parts of the hills (SEDT-7) toward their toe (SEDT-6), but the boundary between U-3 and U-4 deepens 20 m from SEDT-7 (600 masl) to SEDT-6 (580 masl). The slope of this boundary is one of the main differences with the winter geoelectric profile, which is horizontal. In the lowland areas (SEDT-5, 4, 3) the upper part of the terrain (U-2) shows low ρ values (lower than 5 Ohm/m) and in comparison with the ρ data of the winter survey, U-2 extends from the lake shore (SEDT-4) to SEDT-5. U-3 has a similar thickness in SEDT-3, 4, and 5 and the ρ values are slightly lower (18 Ohm/m) towards the lake shore (SEDT-4). The interface between U-3 and U-4 in the lowland areas locates 20 m shallower in SEDT-5 (640 masl) relative to the winter survey, when it was located at 620 masl (SEDT-5BIS) ( Figure 3). The geometry of the interface between U-3 and U-4 in the 2D profile shows a steeper slope between SEDT-6 and SEDT-5, in the summer ( Figure 3B) than in the winter profile ( Figure 3A).

Piezometry
Piezometry in the Duero hydrological basin indicates that regional groundwater flow in the northern half of the basin has a southwestward component, from recharge areas located in the Cantabrian Mountains and the northern edge of the Duero basin to the valley of the Duero river ( Figure 1A) [17,18]. Local deep and brackish groundwater discharge occurs in Villafáfila due to the existence of an elevation of the Paleozoic impermeable basement [20].
Local piezometry also indicates a southwestern component of groundwater flow in the watershed of the Salado river, parallel to the river flow. The hills that surround the lakes have a local component flowing from the upper part of the hills toward the lowlands around the lakes (Figure 4). In the upper part of the hills the groundwater table locates about 700 masl, with depths between 5 and 10 m. In the lowlands around the lakes the groundwater table is at 680 masl, 10-30 cm below the surface. The hydraulic gradient in the upper part of the hills is about 0.1 and in the lowlands is about 0.001. The vertical component of the groundwater flow is downwards in the upper part of the hills, but it is upwards at its toe and in the lowlands (Figure 4). This is evidenced by the differences in hydraulic head between boreholes screened at more than 40 m depth and wells and dug wells screened at depths shallower than 15 m. In the upper part of the hills hydraulic head is lower in boreholes (screened > 40 m) than in dug wells (screened < 15 m). On the contrary, at the toe of the hills and in the lowlands hydraulic head in wells, some of which are flowing wells, is higher than in dug wells [20]. The occurrence of springs around Salina Grande lake shore is an additional evidence for the ascendant component of groundwater flow. Seasonal variations produce minor changes in groundwater level (about 10 cm) in the upper part of the hills. At the toe of the hills many springs dry out during the summer, and groundwater and Salina Grande lake level decrease by about 30-40 cm (Figure 4). Salina Grande lake is normally dry during the summer (from ≈ mid-July till October), although some years stays dry till December, when soils have surplus water and evapotranspiration is low. Lake level reaches maximum depth of 35 to 40 cm during the spring (April-May) ( Figure 5). In SGR-3 piezometer, located in the lake shore (Figure 4), hydraulic head overflows with a minimum discharge at 0.4 m above lake floor. Some springs and flowing boreholes, all of them with brackish Na-Cl regional waters, maintain a constant discharge and temperature all year round. Salina Grande lake is normally dry during the summer (from ≈ mid-July till October although some years stays dry till December, when soils have surplus water and evap transpiration is low. Lake level reaches maximum depth of 35 to 40 cm during the sprin (April-May) ( Figure 5). In SGR-3 piezometer, located in the lake shore (Figure 4), hydra lic head overflows with a minimum discharge at 0.4 m above lake floor. Some springs an flowing boreholes, all of them with brackish Na-Cl regional waters, maintain a consta discharge and temperature all year round.

Electric Conductivity
Electric conductivity (EC) of shallow groundwaters in the Villafáfila area show an increasing trend from the lowest values measured in the upper parts of the hills, toward the lowlands surrounding the lakes and to the lake itself (Figure 6), where a brine occurs within the lake sediment's aquitard [20]. EC data, measured in the summer of 2018, were used to construct with a convergent interpolation method the EC map ( Figure 6), that shows the lateral changes in this parameter. EC was also measured in lake waters, springs and dug wells in locations closest to the TDEM soundings at the time of the summer and winter surveys (in July 2019 and in February 2020) ( Table 1). EC in the shallow groundwaters occurring below the upper part of the hills are between 300-600 µS/cm, with no significant variations along the year (Table 1). Towards the toe of the hills and lowlands

Electric Conductivity
Electric conductivity (EC) of shallow groundwaters in the Villafáfila area show an increasing trend from the lowest values measured in the upper parts of the hills, toward the lowlands surrounding the lakes and to the lake itself (Figure 6), where a brine occurs within the lake sediment's aquitard [20]. EC data, measured in the summer of 2018, were used to construct with a convergent interpolation method the EC map ( Figure 6), that shows the lateral changes in this parameter. EC was also measured in lake waters, springs and dug wells in locations closest to the TDEM soundings at the time of the summer and winter surveys (in July 2019 and in February 2020) ( Table 1). EC in the shallow groundwaters occurring below the upper part of the hills are between 300-600 µS/cm, with no significant variations along the year (Table 1). Towards the toe of the hills and lowlands around the lakes, EC of shallow groundwaters increases to 1000-6000 µS/cm ( Figure 6). The brine within the lake sediment s aquitard has values reaching 42,000 µS/cm during the summer, when the lake is dry. Lake waters in Salina Grande vary along the year from about 1000 to 20,000 µS/cm ( Figure 6). Low EC values in lake waters are registered from February to May. On the other hand, high EC values in lake waters occur at the end of June or along July when the lake is concentrated by evaporation (Table 1). During December and January, there is an increase in lakes salinity, reaching values of 10,000 µS/cm. Groundwaters sampled from boreholes screened below 40 m has EC values ranging from 3000-6000 µS/cm, and show minor variations along the year [20].
Water 2022, 14, x FOR PEER REVIEW 1 around the lakes, EC of shallow groundwaters increases to 1000-6000 µS/cm (Figu The brine within the lake sediment´s aquitard has values reaching 42,000 µS/cm d the summer, when the lake is dry. Lake waters in Salina Grande vary along the yea about 1000 to 20,000 µS/cm ( Figure 6). Low EC values in lake waters are registered February to May. On the other hand, high EC values in lake waters occur at the e June or along July when the lake is concentrated by evaporation (Table 1). During D ber and January, there is an increase in lakes salinity, reaching values of 10,000 µ Groundwaters sampled from boreholes screened below 40 m has EC values ranging 3000-6000 µS/cm, and show minor variations along the year [20].    The electric conductivity record of SGR-3 piezometer along 2020-2021 has been measured at 8 m depth, below the aquitard ( Figure 5). It shows values that vary between 8760 µS/cm and 9030 µS/cm ( Figure 5). The lowest values occur during October while the highest ones are recorded in May and June ( Figure 5). During the summer there is a trend toward decreasing EC, till the end of September-October. From October to June there is an increasing trend in EC. Groundwater temperature in SGR-3 has little variations along the year, from 13.5 • C to 14.7 • C. The highest temperatures occur during December and January, when the lowest air temperatures occur. On the other hand, the lowest groundwater temperatures in SGR-3 occur during May and June ( Figure 5).

Aquifer System Characterization
TDEM resistivity data and electric conductivity (EC) values of lake waters and groundwaters show a consistent spatial distribution. Field geological mapping and well logs indicate that the upper part of the Cenozoic rock record is a sandstone aquifer [18,29,30].
Background resistivity values measured in the Miocene aquifer in the southern part of the Duero aquifer system are about 20 Ohm/m for similar lithologies [26]. Estimates of formation ρ using the measured groundwater EC (Table 1) and Archie s Law [31] for water saturated units with porosities of 0.35 for the Miocene unconsolidated sandstones and 0.2 for the Quaternary lake mudstones using Archie s standard parameters (a = 1.2; m = 1.3) are shown in Table 2. Estimated formation resistivity values are similar to those obtained in the TDEM surveys except for those measured in the proximities of dug well 30-P-006 (SEDT-7 and 7BIS), for which the calculated formation resistivity of 87.5 Ohm/m is notably higher than the ρ measured (25 Ohm/m) for U-1 in the SEDT-7 sounding. This is probably due to the existence of a 7 m thick non-saturated zone above the groundwater table, that if were considered in calculations (Sw-n = 0.3) would reduce the formation resistivity. U-1 is interpreted as the non-saturated part of the Miocene aquifer in the hills (Figure 7). Previous geophysical surveys have revealed the existence of an elevated threshold in the Paleozoic metamorphic basement that forces regional brackish groundwater flow to rise [20] (Figure 7). TDEM profiles perpendicular to the lake's axis (Figure 3) show the occurrence of brackish regional groundwaters (<8 Ohm/m) below 600-550 masl (U-4). The upper parts of the profiles with ρ ranging between 20 to 10 Ohm/m are interpreted as the local fresh groundwaters recharged on the hills (U-3). Piezometry indicates lateral groundwater flow from the hills to the lowlands and to the lake. TDEM profiles evidence the gradual decrease in ρ from the upper part of the hills towards the lowlands. This is consistent with the increase in EC of groundwaters and is interpreted as an increased contribution of the brackish regional groundwaters that mix with fresh groundwaters locally recharged in the hills. Table 2. Conversion of different EC values measured in waters into formation resistivities (ρf) using Archie's law [31]. Sw-n is the water saturation, ϕ is the porosity, and "a" and "m" are the parameters of Archie's formula that refers to tortuosity and cement respectively.   Figure 4. Geological reconstruction is based on the data of [20].

EC (µS/cm)
Underneath the lowlands and the Salina Grande lake bottom, the low ρ values in the upper part of the profile (U-2) are interpreted as a concentrated brine contained within the Quaternary lake sediment aquitard [20] and highly saline capillary waters and the halite efflorescences that form within the soils of the lake margins by evapotranspirative concentration (Figure 7).

Seasonal Changes in Groundwater
TDEM profiles show a similar ρ distribution pattern during summer and winter, but there are some differences that are related to the different hydrological processes occurring in Salina Grande lake and its surroundings. At the lake (SEDT-4BIS) ρ values in the upper part of the profile are low (r < 5 Ohm/m) during winter, due to the occurrence of brackish waters and the presence of the brine contained in the lake sediment aquitard (U-2 in Figure 3). During summer this low ρ zone (U-2) extends to the lowlands surrounding the lake. This is interpreted as an effect of groundwater concentration by evaporation in the upper part of the vadose-capillary zone in an area with mixing of local and regional  Figure 4. Geological reconstruction is based on the data of [20].
Underneath the lowlands and the Salina Grande lake bottom, the low ρ values in the upper part of the profile (U-2) are interpreted as a concentrated brine contained within the Quaternary lake sediment aquitard [20] and highly saline capillary waters and the halite efflorescences that form within the soils of the lake margins by evapotranspirative concentration (Figure 7).

Seasonal Changes in Groundwater
TDEM profiles show a similar ρ distribution pattern during summer and winter, but there are some differences that are related to the different hydrological processes occurring in Salina Grande lake and its surroundings. At the lake (SEDT-4BIS) ρ values in the upper part of the profile are low (r < 5 Ohm/m) during winter, due to the occurrence of brackish waters and the presence of the brine contained in the lake sediment aquitard (U-2 in Figure 3). During summer this low ρ zone (U-2) extends to the lowlands surrounding the lake. This is interpreted as an effect of groundwater concentration by evaporation in the upper part of the vadose-capillary zone in an area with mixing of local and regional groundwaters. This agrees with the occurrence of halite efflorescences on the surface of the lowland areas. The summer profile ( Figure 3B) also shows a rise of the interphase between the local groundwaters coming from the hills (U-3) and the regional brackish groundwaters (U-4), in the lowland areas. The annual record of SGR-3 piezometer (at 8 m depth) shows the highest EC values in May-June decreasing during the summer until reaching the minimum in October ( Figure 5). This decrease in EC ( Figure 5) is interpreted as the rise of the interphase between the brine (U-2) and the waters coming from the hills, that in the lake shore are mixed with the regional brackish groundwaters (U-3). The rise of the interphase is favored by evaporation pumping and capillary rise in the lake aquitard. This is consistent with the summer rise of the U-3/U-4 interphase observed in the TDEM profiles ( Figure 3). In contrast, at the beginning of winter the salinity increases not only in SGR-3 but in the lake waters as well due to the dissolution of halite efflorescences in the lowland areas and the aquifer recharge that favors the transport of the groundwaters and their dissolved salts to the lake area. This processes favor the lowering of the U-2/U-3 interphase in the lake area. Similar processes have been described in other saline lakes where the most soluble salts, mainly halite, can show precipitation-dissolution cycles [32,33]. Other less soluble mineralogical phases can remain, as it occur with the loose calcite cementation observed in the lacustrine mudstones.

Implications for Groundwater Use
TDEM identifies the fresh groundwater domain and is an important tool to locate its contact with brackish or saline groundwaters. Many research papers have highlighted the importance of TDEM soundings to identify groundwater resources [34][35][36]. In Villafáfila, where regional brackish groundwaters occur in many boreholes that then have to be abandoned, it is important to identify the mentioned interphase, especially for the construction of new boreholes for human consumption. TDEM results reported in this paper indicate that the best places to drill boreholes for human consumption of groundwater are the upper part of the hills, and that the lowlands surrounding the lakes should be avoided. The maximum depth of boreholes should not reach below 600 masl on the hills, which is the depth of the fresh (U-3) and brackish groundwaters (U-4) interphase ( Figure 3). Seasonal variations in groundwater dynamics do not affect severely the position of this boundary below the hills, but at the toe of the hills and in the lowlands this boundary can rise significantly during the summer. If intense groundwater pumping is planned in the area, monitoring the fresh-brackish groundwater boundary with TDEM is recommended. The collected data published in this paper and other previous hydrogeological, hydrochemical, and isotopic research in the area could be used to construct a hydrogeological model that will contribute to understand the role of groundwaters in the changes observed in this ecosystem.

Conclusions
TDEM profiles in the surroundings of Villafáfila lakes allowed identifying the lateral and vertical variations in terrain resistivity, which are mainly dominated by variations of groundwater salinity. The occurrence of high ρ zones below the hills (U-1 and U-3) show that these zones, dominated by fresh groundwaters, are the main recharge zones of the local aquifer. This is consistent with the piezometry which indicates that local groundwater flows go from the upper part of the hills toward the low lands and the lakes. The low ρ zones close to the surface toward the Salina Grande correspond to the brine contained in the Quaternary lake sediments aquitard (U-2). The deeper part of the terrain explored by the TDEM soundings shows low ρ that represent the regional brackish groundwater flows. The occurrence of a structural relief that elevates the metamorphic basement of the basin perpendicularly to regional groundwater flows forces the flow to increase, favoring the contribution of brackish groundwaters to the Villafáfila lakes [20,21].
Two TDEM surveys during summer (dry season) and winter (wet season), undertaken in the same location, from the hills toward the Salina Grande lake, have shown that there are some differences in ρ which are the response to the hydrological processes occurring in the lake vicinity. The summer expansion of the low ρ area (U-2) close to the surface in the lowland areas and in the lake is produced by the increase in salinity favored by the evapotranspiration of surficial groundwaters and from the capillarity zone. This evaporative groundwater pumping also produces the rise of the interphase between local fresh groundwaters and the regional brackish groundwaters. The continuous record of EC and temperature in the piezometer SGR-3 confirms the rise of the interphase between the saline brine in the aquitard and the local fresh groundwaters.
EC of groundwaters was measured alongside TDEM sounding acquisition. This has allowed correlating the ρ values of the terrain with those of water and groundwater salinity. The correlation between the measured EC of groundwaters and the ρf values obtained from Archie s law formula is consistent with the ρ values obtained from TDEM soundings. This means that TDEM sounding is an excellent tool to monitor the seasonal changes of groundwater salinity. This can help to prevent aquifer salinization problems. EC map reveals that surficial groundwaters in the low lands has values between 6000 and 9000 µS/cm and the brine below the Salina Grande Lake reaches 42,000 µS/cm.
The best places in the area to obtain fresh groundwaters are the hills that surround the Villafáfila lakes. Below 600 masl in the hills, brackish regional groundwaters are present. In the low lands, the brackish regional groundwaters can be found below 640 masl.