Contribution of Electric Tomography, Chemical and Isotopic Tools to the Understanding of the Resilience and Salinization of Wetlands along the Northern Coast of the Cape-Verde Peninsula: The Case of the Retba Lake

On the northern coast of the Cape Verde peninsula, Lake Retba is one of the few depressions that have preserved its water body due to the drought that plagued the Sahel in the 1960s. The aim of this article is to understand the resilience of the lake to the drought. For this purpose, a deep knowledge of interactions between the different components of the hydro-system, namely: lake, dune tablecloth, ponds and ocean is necessary. The use of electrical resistivity tomography (ETR) and chemical and isotopic tools yielded conclusive results. Electrical tomography of resistivity (ETR) shows throughout the coastal dune a superposition between fresh and salt water, with even a predominance of salt water at the level of the old canal. The chemistry of the major ions and isotopes confirms this contact between the saltwater bison of the sea and the lake by showing a superposition of the fresh/brackish waters on the salty waters of the bevel, which, in turn, are found in the very salty waters influenced by the lake. This translates a feeding of the lake by the sea through the bevel. At the canal level, this feeding is done, on the surface, by means of ponds whose water cha-racteristics therefore of the salt that settles on the bottom. This study shows that the sur-vival of the lake also depends on safeguarding the ponds located in the north of the lake.


Introduction
Lake Retba is located 35 km east of Senegal's capital, Dakar (Figure 1). Its watershed covers an area of 160 km 2 (Carn et al., 1976). The lake was established after post-Nouakchottian marine regression, by isolating the initial lagoon by a spire. The waters are over salted and deposit a salty crust at the bottom of the lake (Garnier, 1978), which has been exploited by the population since the 1970s.
This mining activity produces 110,000 tons of salt annually and generates profit margins ranging from CFAF 1.5 to 3 billion (Sarr, 1997). In addition to this mining activity, there are other activities such as tourism, market gardening and fruit arboriculture. As a result, it plays a significant socio-economic role for the country. It is the only lake, on the north coast, that maintains a permanent running water. However, in recent years, the level of the lake has continued to fall to alarming proportions, threatening all socio-economic activities.
Aware of the importance of this ecosystem, the Senegalese government has integrated the "scientific study of Lake Retba" component into the Mining Sector Support Program (PASMI) financed by the European Union. In fact, the protection of Lake Retba requires a thorough knowledge of the interactions between the various components of the hydro system, such as the lake, the groundwater of the Quaternary sands, the sea and the ponds.
These saline systems are very widespread in arid and semi-arid areas. Their existence requires mainly evaporation higher than precipitations (Hardie et al., 1978;Langbein, 1961). Their typology is based on geomorphology; it was established by Yechieli and Wood (2002). In fact, geomorphology distinguishes flat-surfaced systems without obvious depression or permanent running water (e.g., Sebkhas in the Arabic Gulf), from shallow water basins with seasonal expanse running water (e.g., Texas playas) and well-developed water basins containing permanent running water (e.g., Dead Sea) with the absence or presence of marine, underground and/or superficial influence. Retba Lake might belong to the latter system. For hydrogeologists, these salt lakes are the result of the same hydrogeological process influenced by fluctuations in shallow groundwater. As such, they may be active or representative of past hydro-geological conditions (Yechieli & Wood, 2002).
The first research specifically devoted to Retba Lake dated back to Elouard et al. (1975) and Elouard et al. (1977). This author highlighted and dated the  Sarr et al. (2009a) and Sarr et al. (2009b) highlight the paleo-geographic evolution of the lake by integrating biological markers such as molluscs, foraminifera, ostracods and diatoms. According to this author, the lake went through five stages before its current configuration: a slightly open lagoon stage, a small marine gulf stage, a closed lagoon stage in the process of drying out, a slightly open lagoon stage and a closed lagoon stage (lake). The work of Sow et al. (2008) shows an isolation of the lagoon of Retba Lake consecutive from a period of low sea level between 1200 and 250 B.P. They show the seasonal evolution of the ponds' salinity located in the old channel by the use of diatoms, in particular Paraliasulcata. Sow et al. (2006) mapped the former supplying routes of the lake by the ocean. This feeding was done in addition to the channel, by another route located North East of the lake. Those of Sarr (1997) point out the diachronic evolution of the water body characterized by a continuous reduction in its surface. The works of Martin (1970), Carn et al. (1976) and Garnier (1978) highlighted the hydrochemistry of lake water in relation to that of the groundwater. These authors consider the lake's hydro-system as fixed, fossilized and non-dynamic. They do not address the ocean-groundwater-lake relationship. They consider the lake's hydro-system as a relic of the salt bevel, with no hydrodynamic relationship with the ocean and evolving in an evaporative tank. This does not explain the lake's resilience to the drought of the 1970s as well as the salinization of its waters. The first studies suggesting that the lake is supplied by the ocean, and therefore having a dynamic system that is being renewed, are those of Gueye et al. (2016). This study corroborates with these first tioning of this ecosystem, which led to the salinization of the water and allowed the lake to maintain itself against the drought wave of the 1970s. This helps to identify strategies for managing and safeguarding the resource, that is the salt exploited at lake level by local populations. To achieve this objective, an experimental system composed of four piezometers has been set up on the coastal dune separating the lake and the ocean. A stratified sampling of the groundwater was carried out on these piezometers. The lake and ponds were also sampled.
Chemical analyses were performed on the various samples taken from the lake, groundwater and ponds. Three geophysical profiles were carried out to measure the variation of vertical resistivity as a function of depth. Two analyses, on the dune and another on the old channel of the lagoon currently occupied by sandy dunes. These resistivity measurements were performed using the Electric Resistivity Tomography (ETR) method.

Hydrogeological and Hydrological Context
The hydrogeology of the area is dominated by the quaternary sands. These formations, whether marine or continental, behave hydro geologically as a single reservoir directly in contact with the ocean to the west (Vallet, 1972), and its waterproof is variable and depends on the nature of the aquifer formations (Martin, 1970). According to Garnier (1978), the offshore groundwater flows on one hand to the ocean and on the other hand to the lakes.
Permeability is highly variable, as well as transmissivity which is between 5.2•10 −4 and 5.6•10 −3 m 2 /s. In Niakoul Rap and Bambilor located in the southern part of the catchment area ( Figure 1), it is respectively 2.3•10 −3 m 2 /s and 1.8•10 −3 m 2 /s. The storage coefficient is close to 1%, but can vary due to the presence of some of the more clayey horizons found in the sands. In some areas, the groundwater acquires a semi-captive behavior due to these clayey horizons (Vallet, 1972 (Ferry et al., 1998;Ouédraogo, 2001).
In consideration to the normal 1951-1980, the series can be divided into two periods ( Figure 2) according to the evolution of the rainfall index: a period from 1900 to 1969 and a period from 1970 to 2013:  The period from 1900 to 1969 The evolution of rainfall indices ( Figure 2) shows surplus rainfalls that vary between 1% (1900) and 83% (1918, 1951 Henry (1972), Vallet (1972) and Lenclud (2005) who estimate that evapotranspiration varies between 1500 mm and 2200 mm per year in the Cape Verde Peninsula.
From this analysis, we can conclude that our area has been marked by drought and high evaporation since the 1970s. These factors combined should lead to the disappearance of surface water bodies. This is not the case for Retba Lake, hence the importance of observing its functioning which allowed its resilience.

Theory and Basic Formulas
For four electrodes, for instance an AB-MN quadrupole of any geometry, the principle of potentials' additivity allows to write (Equation (1)): In the field, the resistivities' distribution is generally heterogeneous and anisotropic. Thus, an apparent resistivity ρ a (in ohm•m) is like the resistivity equivalent to a homogeneous and isotropic environment. This apparent resistivity ρ a (in ohm•m) is given by the following relationship (Equation (2)): K is therefore a factor that only depends on the configuration of the quadrupole geometry (ABMN): It is the geometric factor. Apparent resistivity is the variable measured in the field in direct current electrical prospecting and tomography. It does not directly reflect the actual resistivities or depths of the subsoil materials. As a result, the data measured in the field are processed using inversion software to calculate the true resistivities of the materials: it is inversion or deconvolution. It allows having a section of the basement in true resistivities that will be interpreted.

Application in the Coastal Dune
Using an LS ABEM terrameter, three ETR profiles were performed on the coastal dune between the ocean and the lake. The first two profiles are juxtaposed and have a line length of 630 m of cable from the beach to the approach of the lake ( Figure 3). The other profile was installed in the old channel of the lake that fed the lagoon (Figure 3). The profile is transverse to the channel and measures 210 m from the beach to the depression of the old channel that include pond 2. The device used is the Wenner. The topographic altitude of all electrodes in the three profiles was determined using a differential GPS (Global position system). These data were acquired after the wet season (September 2, 2015) in order to benefit from the moisture in the sands, which facilitates the passage of the current into the deep layers.

Data Processing and Inversion
The RES2DINV software version 3.54.44 was used for data processing and cutting of real resistivity. Once the data has been pre-processed, the modeling parameters are selected before loading the file into the RES2DINV program.  constraint (robust inversion), which implies a more restricted Laplacian distribution of data errors (Al-Chalabi, 1992). The use of a lower standard (L1) is required with data affected by non-random noise.
The quality of the acquisition is checked by visualizing the file where the different levels are represented with a reduced scale ( Figure 4). This representation highlights measurements with very low or very high apparent resistivity in relation to neighboring points. Since such a rapid change cannot be due to a geological phenomenon, such data must be eliminated (Marescot, 2008).
The data acquired are homogeneous since the contact problem did not arise due to the humidity of the dunes. This is illustrated by Figure 4, which shows the acquisition of profile 1 data. After saving the processed file, its modeling inversion is started. The result is subsurface sections representing the vertical and horizontal variation of resistivity. It is these sections that are presented in the result section in the form of profiles.

Isotope Analyses
The stable isotope contents of water 18 O and 2 H were measured as follows: 5 ml of sample were distilled at a temperature of 90˚C, and then at 105˚C for one hour in an Evapoclean, which is a closed evaporation device; this device avoids a high salt concentration for subsequent analysis of the stable isotopes. The oxygen-18 and deuterium contents, were determined using a Picarro L1102i WS-CRDS (Digital Wavelength Ring Cavity Spectroscopy System), as described by Halder and Decrouy (2013). Each sequence was calibrated using three differ- standards. The isotopic compositions measured are reported in delta units (deviation expressed in ‰) from the isotopic ratio of the sample to that of V-SMOW, standardized so that SLAP has a δ 18 O = −55.5‰ and a δ 2 H = −428.0‰ (Coplen, 1994;EOS, 1996). The uncertainties on repeated measurements, standards and samples are ±0.1‰ for oxygen −18‰ and ±1.4‰ for deuterium.

Electrical Resistivity Tomography (ERT)
Two transects were carried out: the Dune's transect, consisting of profiles 1 and 2, and the channel's transect, consisting of profile 3. The profiles are south (lake side) to north (ocean side) oriented. Profiles 1 and 2 are successive.

The Dune's Transect
Profile 1: The true resistivity model of this profile (Figure 6), obtained after inversion of the apparent resistivities, with an RMS error of 4.9%, shows five contrasting sets: 1) A superficial unit whose minimum resistivities are around 1000 Ohm•m and maximum resistivities exceeding 6000 Ohm•m. This set corresponds to the dry sands of the dune, which constitute the very resistant dune cover. Its thickness varies between 2 and 3 m.
2) A resistant unit, with resistivities of around 500 Ohm•m. This set is observed just after the very resistant superficial set. This sudden fall in resistivity suggests that this set corresponds to wet sands, i.e. the unsaturated area. This layer is very thin, about 0.5 m thick.
3) A moderately resistant unit, with resistivities between 40 and 200 Ohm•m. According to Archie's law, the resistivities of the imbibition water vary between 8 and 40 Ohm•m. These resistivities translated into conductivities give a range of conductivities that are from 250 to 1250 μS/cm. This set corresponds therefore to the sandy fresh water table of the dune's cordon. According to this profile, its maximum thickness is located at 560 m from the ocean and it is 18 m. 4) A moderately conductive unit, with resistivities between 10 and 40 Ohm•m. According to Archie's law, the resistivities of the imbibition water vary between 2 and 8 Ohm•m. These resistivities, translated into conductivities, give a range of conductivities ranging from 1250 to 5000 µS/cm. This set could therefore correspond to the transition zone between the underlying highly conductive set and set three.  Profile 1 illustrates well the contact between freshwater and salt water and shows that it is a slope that sinks into the dune and then rises towards the ocean.
The fresh water's highest thickness is 18 m at a point located at 560 m from the ocean. It also shows that fresh water forms a weak lens over salt water.  At the last level, lateral variation of the resistivity is noted. It may also be due to the variation in the proportion of clay in the aquifer sands or to changes in the resistivity of imbibition water or both.
On profile 2 (Figure 7), the freshwater-salt water contact occurs from the south to the north towards the ocean. This profile extends the first, so the recovery observed at the first level continues. This shows that the thickness of the freshwater sands aquifer is becoming increasingly weak by moving to the ocean, as shown in Figure 7.

The Channel Transect
The channel profile was reversed with a 20% RMS error. However, our knowledge of the field allowed us to interpret it as it was considered fairly consistent with our observations on the ground. Three contrasted sets have been identified ( Figure 8): 1) A very resistant set with resistivities above 1100 Ohm•m; this set corresponds to the sandy, very resistant cover of the dunes.
2) A moderately resistant set with resistivities between 50 and 250 Ohm•m. The resistivities of the imbibition water obtained from Archie's law therefore vary between 10 and 50 Ohm•m. These resistivities translated into conductivities give a range of conductivities ranging from 40 to 1000 µS/cm. This set corresponds to the freshwater/brackish r table. Its thickness decreases towards the ocean and is almost cancelled at the end of the profile (north side). Along the profile, small freshwater/brackish lenses are noted.
3) A very conductive set with very low resistivities, ranging from 0.024 to 11 Ohm•m. Water resistivities, obtained from Archie's Law, vary between 0.0048 and 5 Ohm•m. These resistivities translated into conductivity give a range of conductivities ranging between 200 and 2.08•10 5 mS/m. This set corresponds to the very salty water sands of the channel. These very high conductivities, suggesting a high salinity, could be explained by the fact that this area corresponds to the old channel and therefore there is an ancient salinity added to the current Journal of Geoscience and Environment Protection Figure 8. Section of true resistivities, obtained on the coastal dune between the ocean and the lake at the site of the old channel, reversal made the application of the stress L1-standardized. salinity due to the existence of the bevel. In the southern part it is noted an area of very low resistivity compared to the rest, these very low resistivities may be due to infiltration into the water table, the very salty pond or to clay enrichment, sands, next to the pond2.
On this profile (Figure 8), a slight north-south slope, i.e. from the ocean to the pond2, with a gradient of 0.0125 is observed. Pond2 is located at an altitude of −1.47 m, pond 1 to −3 m and the lake at −5.52 m. This shows a gradient from the ocean to the lake and promotes the flow of the bevel to the lake through both ponds. This result is consistent with the work of (Yechieli & Wood, 2002) on coastal saline hydrosystems and their relationship to the ocean. The freshwater/brackish water table does not appear immediately as a lens as seen at the profile level, which is probably due to its low thickness. However, small lenses are observed along the profile.
These results are almost consistent with those of (Comte, 2008), obtained by applying electric tomography to the Pikine-Dagoudane depression site, located west of the lake. Resistivity values also fall within the scale given by Debuisson (1965) and Moussu and Debuisson (1966). These geo-electric profiles show that between the lake and the ocean, there is a salty or even very salty water table that connects the lake to the ocean. This water table is overcome by the fresh/brackish water table of the dune cord. The maximum thickness of this fresh/brackish water table does not exceed 20 m. This maximum thickness is observed at a position 560 m from the ocean. From this position, the thickness of the table decreases by heading towards the lake and the ocean. The fresh water table floats on that of salt water table as a lens. However, it is not as powerful as it is to create a hydraulic barrier that prevents the lake from eating by the ocean, as reported (Garnier, 1978). Indeed, the lens is not capable of repelling the fresh water-salt water interface according to the proportions predicted by the law of Ghyben-Herzberg. Piezometric levels, measured on piezometers, show that the water table is below ocean level.
The ETR profile performed on the old channel shows that the fresh/brackish waters of the groundwater table at channel level is based on the salt waters of the

Vertical Evolution of Dune Groundwater Chemistry
During

Principal Component Analysis (PCA)
Choice of the main components We choose to study only the two main components. The choice of these main components is justified by the histogram of variances according to the different components ( Figure 9) on which the scree curve is drawn (red curve). This curve shows that the axes associated with the eigenvalues before the break represent the largest variance. These are axis 1 (Dim1) and axis 2 (Dim 2).
These two main components (axis1 and axis2) ( Figure 10, Figure 11) account for 79.33% of the total variability, of which axis1 explains nearly 58.14% of the variations and axis2 25.24%. This reflects the quality of the data since only 20.67% of the information is lost. The factorial design will therefore consist of these two axes. The analysis will be carried out according to these two axes.

Study of the cloud of variables
The representation of the variables on the correlation circle ( Figure 10    The group of variables whose vector ends representing them is close to the Journal of Geoscience and Environment Protection origin (i.e. 0). These variables are not significant and are therefore poorly represented. These are nitrite ( 2 NO − ), phosphate oxide ( 2 4 PO − ), pH and fluorine (F − ). These variables will not be included in the analysis since they do not contribute to the variance of the chemical composition of the water.
The variables: conductivity (cond), bicarbonate ( 3 HCO − ), ammonium ( 4 NH + ), bromine (Br − ), deuterium ( 2 H), oxygen-18 ( 18 O), magnesium (Mg 2+ ), sodium (Na + ), calcium (Ca 2+ ), chlorine (Cl − ) and sulfate ( 2 4 SO − ) are strongly correlated with axis 1 (Dim 1). These variables therefore contribute to the formation of this axis. They are strongly correlated with conductivity. On the other hand, the variables altitude and distance (from the sea) are not correlated with this axis. This axis (Dim 1) therefore contrasts the variables that contribute to salinity with those that do not. It thus defines a salinity gradient.

Study of the cloud of individuals
The projection of the individuals (samples) on the factorial plane consisting of axis 1 and axis 2 gives Figure 11. It shows that the individuals (sample waters) can be separated into five groups: G1 (sample waters from PB and PC piezometers), G2 (surface waters from the first three samples from the PD piezometer), G3 (intermediate waters from the two samples from the PD piezometer), G4 (sample waters from the PA piezometer) and G5 (deep waters from the PD piezometer). Individuals in groups G1 and G2 have negative coordinates on axis 1 (Dim 1) and those in groups G3, G4 and G5 have positive coordinates on the same axis. Samples of groups G3, G4 and G5 therefore have higher values for the following variables: conductivity (cond), bicarbonate ( 3 HCO − ), ammonium ( 4 NH + ), bromine (Br − ), deuterium ( 2 H), oxygen-18 ( 18 O), magnesium (Mg 2+ ), sodium (Na + ), calcium (Ca 2+ ), chlorine (Cl − ) and sulfate ( 2 4 SO − ) than those of groups G3 and G2. Axis 1 therefore contrasts individuals in groups G1 and G2 who have low to zero salinity with individuals in groups G3, G4 and G5 who have high salinity. This axis therefore represents salinity. It thus defines a salinity gradient, ranging from fresh/brackish water, consisting of the waters of the PC and PB piezometers and the waters of the first three surface samples of the PD piezometer, to salt water, consisting of the waters of the deep and intermediate samples of the PD piezometer and the PA piezometer. The factorial design shows that the surface water of the PD piezometer is closer to the water of the PB and PC piezometers. These waters are representative of the freshwater/brackish sand table of the barrier boom. The waters of the intermediate samples of the PD piezometer (G3) are closer to the waters of the PA piezometer samples (these samples have coordinates between 2 and 4 on axis 1) which are under marine influence due to their proximity to the ocean. These intermediate waters are therefore influenced by the bevel of the sea. The deep sample waters of the PD piezometer stand out completely from the other groups with coordinates between 5 and 6 on axis 1 (Dim 1). These waters therefore show a high salinity higher than that of the waters of the PA piezometer, which is closer to the sea. This salinity is therefore justified by a lacustrine influence. The water from the lake seeps into the dune, increasing the salinity of these waters.
From this analysis, it appears that at the level of the coastal dune, there are three types of water: the fresh/brackish water table represented by groups G1 and G2, the salt water table of the sea bevel represented by groups G3 and G4 and the very salt water table represented by group G5. Groups G2, G3, G5 are waters of the piezometer PD. Therefore, at the level of this piezometer we note a stratification of these three types of water. This piezometer contains water from the water table of the coastal dune (fresh/brackish water) which is found on the waters of the sea bevel, which rests on the very salty water probably coming from the lake bevel. This shows a continuity of the sea bevel at the level of this piezometer and a superposition of ocean waters on the lake waters.  Figure 12). The slope of the regression right is equal to 6.5; it is 8 for the world meteoric waters and 7 for the local one. This shows that these waters have slightly evaporated. The sample distribution includes the five groups determined from the CPA: G1 and G2 groups: the waters of the samples have oxygen-18 levels, which vary between −5.6 and −4.3. Deuterium content varies between −33.9 and −30.6. The waters of the samples of these groups plow on the local meteorological right, traced through the isotopes of local rainwater by Diouf (2012). They are therefore an isotopic signal close to that of rainwater. They are representative of the freshwater table of the coastal dune.  to that of seawater, which is 0.73‰. The fraction of sea water in water under marine influence may be calculated by the following formula (Equation (3)) (Appelo & Postma, 2005):

Isotopic Study of the Waters of the Groundwater of the
With f seawater = fraction of seawater in percentage, M is the mmol/l concentration of the sample chloride, m is the concentration of chloride ions in mmol/l freshwater, S is the concentration of chloride ions in mmol/l of seawater.
For the groups G3 and G4 groups the fraction of seawater varies between 31 and 51%. This shows that the groundwater at the level of the PA piezometer is under marine influence as well as the water at the intermediate level of the PD piezometer.
The G5 group consists of the waters of the deep samples of the PD piezometer. Their oxygen-18 content varies between −1‰ and 0‰ and those of deuterium between −10‰ and −1.5‰ while the oxygen-18 levels content of PB and PC piezometers vary between −5.6‰ and −4.9‰ and deuterium of these piezometers varies between −36.9‰ and −33.2‰. The lake has 18 oxygen levels ranging from 0.7‰ and 1‰ and deuterium levels ranging between 0.6‰ and 1‰. The deuterium content of seawater is 4.2‰ and that of oxygen-18 is 0.7‰. This shows, on the one hand, that the waters of the lake have an isotopic signature near those of seawater and that the deep waters of the PD piezometer are close to those of the lake. These waters are influenced by lake waters that infiltrate the dune.
This isotopic study confirms the distribution of the groups displayed by the CPA and the relationship between G1 and G2 groups and between G3 and G4 groups on the other. It also confirms that the G5 group's waters of do not belong to the same type as the other representative groups of the groundwater (fresh and salty); these waters are closer to the lake. This confirms the stratification of the waters at the level of the PD piezometer where the water of the fresh/brackish dune groundwater (G1 and G2) which is based on water from the sea bevel (G2 and G3) which in turn is based on very salty water (mixing water of the sea and water of the lake) (G5).

Deuterium/Chloride Relationship
The "heavy isotopes/chlorides" relationship and the discretization of the samples based on their depth of sampling also the identification of sample groups identified by the PCA: It is a question of: G1 and G2 groups: G1 group has low deuterium levels (ranging from 36.9‰ to 33.2‰), with chloride concentrations varying between 1 and 3 meq/l. The G2 group exhibited low deuterium levels (ranging from −33 to −30‰) and chloride concentrations varying from 40 to 55 meq/l. Their deuterium levels are therefore close to that of the G1 group but they have a high salinity.
These waters represent the fresh water table of the coastal dune. The waters of G2 group (surface waters of the PD piezometer) are part of this aquifer but they have been influenced by deep salted waters resulting in relatively large chloride concentrations.
G3 and G4 groups: They have deuterium levels ranging between −25.6 and −15‰. Their chloride concentrations vary between 200 and 300 meq/l. They are more enriched in deuterium than the waters of the PB and PC piezometers. Figure 15 confirms the proximity, already seen at the CPA level, between these two groups (G3 and G4) and between these two groups and the seawater sample on the other hand. This shows that these waters are under marine influences and their chloride concentrations corroborate this.
In addition, G5 consisting of the waters of the deep samples of the PD piezometer. They show deuterium that varies between −5‰ and −1.5‰ and chloride levels ranging from 300 to 500 meq/l. Their deuterium isotopic signature approximates those of lake water and their chloride concentrations are closer to that of seawater ( Figure 13). This confirms the results already reported on the origin of the deep waters of the PD piezometer. These waters have marine origin, but have a lake influenced from the infiltration of lake water into the dune.
According to Figure 13, the ponds have an isotopic signature and chloride concentrations close to the seawater sample. This confirms the geophysical results already mentioned: These ponds are an outcrop of the sea's biswater. In addition, the lake's feeding by the sea, at the channel level, passes through these surface ponds, and by the deep biswater.

Relationship between the Br − /Cl − Ratio and Chlorides (Cl − )
The evolution of bromides (Br − ) based on chlorides (Cl − ) ( Figure 14) showed a strong correlation between Br − and ions (Cl − ), with a correlation coefficient equal to 0.959. This means that both ions have a common origin or come from the same process.
The Br − /Cl − ratio ( Figure 15) highlights the influence of seawater on the coastal groundwater. This ratio, bromide-chloride, is constant and equal to 3.47 × 10 −3 for seawater (Marjoua et al., 1997). When the values of the Br − /Cl − ratio of fresh water are close to or greater than 3.47 × 10 −3 this translates into a seawater flow into the groundwater. This report is therefore used to identify the saline origin of water (Faye et al., 2005;Marjoua et al., 1997;Richter & Kreitler, 1993;Viviana et al., 2010). The surface waters of the PD piezometer (G2) and the wa-    varies between 1 and 3 meq/l. By comparing these values with the seawater ratio, it can be said that these piezometers are not or slightly influenced by seawater. The lower chloride values observed at the surface water level of the PD piezometer were explained by its proximity to the lake. The waters of the PA piezometer (G3) and the waters of the intermediate sample of the PD piezometer have a Br − /Cl − ratio ranging from 0.2723 to 0.4552 and chloride concentrations varying between 230 and 300 meq/l. These ratios Br − /Cl − higher than that of seawater indicates that the waters of this piezometer are influenced by sea water through its salty biswater.
The deep waters of the PD piezometer have a Br − /Cl − ratio that varies between 0.00153 and 0.00157. Their chloride concentration varies between 300 and 500 meq/L. Their Br − /Cl − ratios are lower than that of seawater but their chloride concentrations are close to the latter (550 meq/l). This shows that the deep waters of the PD piezometer also have a marine influence by the arrival of the bis-water at the level of this piezometer. Their isotopic signature also advocates for lake influence.
The water of the lake has a Br − /Cl − ratio equal to 0.00344, thus substantially equal to that of seawater. Its chloride content is 2000 meq/l while that of seawater is 550 meq/l. This clearly shows that lake waters are marine waters that have undergone evaporative concentration. These isotopic tools show that this diet passes through the dune via the salt bevel of the sea and from the old channel through the ponds.

Conclusion
From these developments, it can be concluded that isotopes, water chemistry, and geophysics have elucidated the operation of the lake's hydro-system. The ocean and the water table of the coastal dune influence the lake. During the dry season, it receives continental inputs. The superposition of "fresh water/salt water table" is also attested by geo-electric models. The same is applied to water chemistry, such as the Br − /Cl − ratio, which confirms that the lake's waters are indeed of marine origin. The old channel is also a witness to this lake feeding by the ocean. There is a strong resurgence of salted bevel of the sea in the ponds that mark it. The lake is thus an evaporative pond, the permanent source of which is the ocean and the freshwater lens of the coastal dune. There is such a concentration of water there, by evaporation, that the salinity of its waters reaches 9 times that of seawater. During the rainy season, there is also runoff, through the hydrographic network, which provides a slight softening of the water. In the dry season, the concentration is again triggered. This operation explains the salinization of these waters as well as its resilience to the drought of the 1970s because it has a permanent source of food that is the sea. Protection of the ponds must therefore be taken to prevent the lake from drying out.