The Wadi-Fira region in eastern Chad, adjacent to Sudan (Figure 1), is located in an arid environment, exposed to the climate changes, where crystalline basement is mainly exposed. This region is regularly faced with recurring problems of water supply. This situation is linked to the climatic context of this region with a rather low rainfall, but also to the lack of knowledge about the basement aquifer.
However, the basement aquifer of Wadi-Fira constitutes the main source of water supply for the populations of this region. The aquifer system has two main types of reservoirs, the alterites reservoir and the fractured reservoir. Most of the underground water resources exploited by deep wells and hand dug wells are located in alterites and basement fractured rocks. Sedimentary formations play a negligible role.
In 2003, the Wadi-Fira region hosted »201,000 refugees from Darfur out of a total of »924,000 residents that is 22% of the population. This migratory flow led to an expansion of agglomerations (temporary new towns), an increase in agricultural and pastoral practices. The groundwater of the crystalline system is thus now heavily exploited to meet the dramatically increasing water requirements in recent years.
The rational exploitation of this vital resource requires a prior assessment of
Figure 1. Location and geological setting of the study site. Wadi-Fira region (Eastern Chad).
the processes of mineralization of the groundwater in order not to degrade its quality by a disorderly overexploitation.
Within this framework, the present investigation aims at better understanding and identifying hydrogeochemical processes and their relations to groundwater quality within this complex environment, and groundwater recharge mechanisms.
The geological formations encountered in the study area are mainly granitoids, schists, migmatites and volcano-sedimentary formations (basalts and alluvium) covering the crystalline basement. Two regions were sampled to conduct this work, the Am Zoer region and the Guereda-Iriba region.
2. Geological and Hydrogeological Framework
The study area lies between latitudes 14˚ and 16˚North, and longitudes 21˚ and 23˚East. It covers an area of approximately 50,000 km². It is characterized by a Sahelo-Saharan climate, with a great spatial and temporal variability in precipitation, i.e. 150 to 400 mm of rain per year, between the North and the South. Most of the rainfall is concentrated over two months, between July and August. The Am Zoer site is located in the south-west, and that of Guereda-Iriba in the northeast (Figure 1). The particularities of these two contexts will be analyzed from the hydrogeochemical point of view related to the aspect of the water resource.
The geological context of the study area (Figure 1) was elaborated from the geological map of Niéré at 1/500,000  .
The formations of the Precambrian basement identified in Chad, of which the Ouaddaï massif forms part, covers between 15% and 20% of the territory of the country  . The formations encountered are granitoids, metamorphic rocks (schists, quartzites, micaschists, gneiss, green rocks and migmatites) and volcano-sedimentary formations. The granitoids are composed of granodiorites and granites that are exposed on most parts of the basement of the Ouaddaï. The volcanic formations occur in basaltic flows along the wadi Absounout near Iriba, or in isolated pitons with circular base about 200 m in diameter and with steep slope. The sedimentary formations consist of the Palaeozoic (Cambrian) sandstones, the Continental Terminal sandstones and alluvial formations (sand, clay, silt).
From recent results, it is observed that in basement hard rocks, aquifers are located within the first tens of meters in the rocks    . The typical vertical profile comprises from top to bottom: i) The laterite layer with a few meters thickness, ii) The alterite (or regolith) layer, a few tens of meters thick, derived from the weathering of the basement rocks; iii) The fractured layer which consists of the massive rocks. The fracture network mainly results from the tectonic activity in the region, but also from the weathering process itself   .
In the Wadi-Fira study area, groundwater is exploited by shallow large-di- ameter wells (average depth »10 m) dug in the alterites and deeper wells penetrating the fractured aquifer (mean depth »45 m).
The hydrodynamic properties are very variable in space, showing the strong heterogeneity of this system. Transmissivities estimated by pumping tests on deep wells range between 1.9 × 10−6 to 8.8 × 10−4 m2/h. Discharge rates are generally low. More than 80% of wells have discharge rates between 0.5 and 5 m3/h. Flow rates exceeding 10 m3∙h−1 were obtained at Diker (Am Zoer site) and on two wells at Guereda. Such high discharges are linked respectively to a thick sedimentary layer (20 to 40 m) and a well-developed (10 and 26 m) weathered and/or fractured horizon.
The piezometric map is reported in Figure 2. This map does not differentiate the laterite reservoir from the fractured rocks reservoir. Given the scale of the study area, a continuous system can be assumed. Interpretation of the map indicates that groundwater is mainly drained towards rivers beds.
3. Sampling and Analytical Procedures
Two sampling campaigns were conducted in 2011 and 2015 for a total of 31
Figure 2. Piezometric map of the study area.
groundwater samples. The first campaign involved 15 deep wells. The analyses were carried out at the National Laboratory of Water Analysis (LNAE) in Ndjamena (Chad). The second field mission in July 2015 resulted in the sampling of 16 water samples from boreholes (9) and dug wells (7). The chemical and isotopic analyzes were carried out at the Center of Hydrogeology and Geothermics (CHYN) of the University of Neuchâtel in Switzerland.
During both field missions, the physical parameters (electrical conductivity, pH and temperature) of the groundwater were measured in situ. The data of the conductimeter during these two campaigns proved unfortunately wrong.
Other complementary parameters were also measured, the location and altitude of the water points using a Garmin GPS, the piezometric level and the depth of the wells respectively by a piezometric probe and a decametric tape.
Analyzes of the major ions, cations and anions were carried out respectively by ion chromatography Dionex DX-120 and Dionex-ICS-1600. The determination of the bicarbonate ion ( ) was carried out by the titrimetry method using an apparatus of the Metrohm 805 Dosimat type. The quality of the analyses was verified by the calculation of the ion balance.
The isotopic analyzes were carried out only on the samples collected during the 2nd campaign in 2015. The stable isotopes contents of the water molecule are measured by the Picarro L2140-i analyzer, which performs simultaneous measurements of δ18O and δ2H with an accuracy of ±0.0015‰.
All stable isotopic composition (18O and ²H) are reported in standard δ‰ notation (part per thousand versus VSMOW (Vienna Standard Mean Ocean Water)) as follows:
where Rsample and Rstandard represent the ratio of heavy to light isotopes of the samples and standard, respectively  .
4. Interpretative Methods
The results of the physico-chemical and isotopic analyzes carried out on the 31 groundwater samples are interpreted using statistical, graphical and hydrogeochemical methods. Statistical interpretations were performed using box-plots, correlation matrice and Principal Component Analysis (PCA).
The hydrochemical classification of water is developed using Piper's trilinear diagram  . The study of hydrochemical processes is based on the calculation of the saturation indices and scatter diagrams of chemical elements.
The Phreeqc geochemical program allowed to evaluate the saturation indices, to characterize the geochemical processes at the origin of the mineralization of groundwater.
Isotopic data on stable isotopes (oxygen-18 and deuterium) in groundwater have been interpreted using conventional IAEA methods  .
The application of all these different methods to the hydrochemical and isotope data allowed to characterize the mechanisms of mineralization of groundwater and the recharge processes in the basement aquifer of the Wadi-Fira area.
5. Groundwater Chemistry
5.1. Descriptive Statistics
The physico-chemical data of the groundwaters samples are reported in Table 1.
Table 1. Hydrochemical data from samples collected in the field at the 2 sites, Am-Zoer (AZ) and Guereda-Iriba (GI). P: hand dug wells; F: deep wells.
The descriptive statistics (average, variance, coefficient of variation) are reported in Table 2. The groundwater samples are of low overall mineralization. The average mineralization is 336 mg/l and is in the range of 106 to 628 mg/l, showing that the groundwater in the study area is fresh.
The pH of the groundwater, for both hand dug wells and deep wells ranges between 4.6 and 8.9 with an average of 7.0. More than half of the pH values (51%) are in the range 5.5 - 7.5. Such intermediate natural-water pH values implies a balance between the production of H+ ions from the dissociation of weak acids, and H+ (and CO2) consumption by weathering reactions, as the weathering of albite (Na-feldspar) to kaolinite  .
The weak acids involved in these processes are carbonic acid and organic acids such as fulvic acid  . These weak acids are the cause which provokes the weathering of silicate and alumino-silicate minerals within such basement environment.
Higher pH values are also observed in the study area with values in the range 7.5 - 8.5 common in groundwaters in hard rocks (basalt and granite). The maximum observed pH value for groundwaters is 8.9. High pH values occur in waters whose composition is dominated by minerals such as silicates and aluminosilicates, which tend to raise the pH  . Groundwaters with pH values greater than 9.0 are classified as’rock-dominated’  , meaning that they are likely not to be recharged very fast. All pH values higher that 7.5 characterize deep wells in the study area, showing that groundwater recharge in such sites is slower.
The groundwater temperature varies from 16˚C to 33˚C with an average of 28˚C. These temperatures are close to the air temperature in the region of 29˚C, showing that recharge takes place by infiltration of rainfall, in the absence of
Table 2. Descriptive statistics of the hydrochmical parameters. CV: coefficient of variation; SD: standard deviation.
very deep circulations.
An analysis of the hydrochemical data, carried out using Box-plots, is shown in Figure 3. The purpose of this plot is to demonstrate a possible differentiation between subsurface groundwater sampled in wells dug in alterites and the deep groundwater sampled in wells penetrating the fractured crystalline reservoir. The mean values of the variables are also reported in Table 3.
The box-plot (Figure 3) compares the parameters of waters from the boreholes and dug wells. This results in significant differences in the values of T˚C, pH, Ca2+, Na+, K+, Cl−, , , and F−.
Figure 3. Box plots representation of the hydrochemical parameters sampled in hand dug wells and in deep wells.
Table 3. Hydrochemical average values in deep wells (D.W.) and in hand dug wells (H.D.W.) Depth in m. All ions concentration in mg/l.
All box-plots show that concentrations of the elements are higher in deep wells than in wells dug in the alterites except for K+ and . This highlights the increasing mineralization of infiltrating waters and that mineralization processes are more pronounced in fractured rock. The average mineralization of dug wells water is lower (average TDS = 235 mg/l) than that of deep wells (average TDS = 357 mg/l).
5.2. Correlation Matrice
Pearson correlation matrice  was calculated between hydrochemical variables to display their respective statistical relationships (Table 4). Significant correlation coefficients, at a level of 0.05 (p < 0.05), are highlighted in the table.
Strong correlations are found between Ca2+ and (R = 0.77), Mg2+ and (R = 0.68), Na+ and (R = 0.78), Ca2+ and Cl− (R = 0.73), Na+ and F− (R = 0.73). Strong relationships between Ca2+, Mg2+, indicate contribution from dissolution of or . Relations between Ca2+ and Cl−, Na+ and would result from complex processes discussed below.
Fluoride content is rather an important issue regarding potability of groundwater. According to the World Health Organization  , water with a content higher than 1.5 mg/lis unsuitable for drinking. In the Wadi-Fira area, F− content ranges between 0.1 mg/l and 2.9 mg/l, with an average of 1.0 mg/l, indicating that at some places, the groundwater is unsuitable for drinking. It is thus important to understand the origin of F− in groundwater.
The high correlation between F− and Na+ should be noticed, providing an indication on the possible source of F−. Dissolution of fluoride-bearing minerals
Table 4. Correlation matrix between variables.
such as cryolite (Na3AlF6) was depicted as a source of F- in hard rocks groundwater  . However further investigations are still required to deepen this issue in the Wadi-Fira area.
Cations that can exchange through ion exchange processes exhibit moderate correlation (Na+-Ca2+: R = 0.45; Na+-Mg2+: R = 0.51). This indicates that cation exchange processes can affect the groundwater chemical composition, but are not dominant processes.
5.3. Principal Component Analysis
Principal Component analysis is a well-established and widely used multivariate method in hydrochemical studies. As groundwater mineralization is a multivariate concept, analyzing the hydrochemical parameters pairwise, may lead to ignore relevant relationships  . Thus, a Principal Component analysis (R mode factor analysis,   was undertaken in this work in order to extract significant information about analyzed data set and better display the inter-relations between the variables. A PCA may also contribute to understand the underlying geochemical reaction and hydrochemical evolution processes.
The data processed by the ACP consists of 12 hydrochemical elements measured on 31 samples from the Am Zoer and Guéréda-Iriba sites. The eigenvalues and the percentages of variance explained by each Principal Component (PC) are shown in Table 5.
The PC1 component with an expressed variance of 37.2% is the most important of all the factorial axes. The PC2 component expresses the greatest residual variance, i.e. 18.6% of the total variance. The factor plane PC1-PC2 expresses 58.8% of the cumulative variance. As can be seen on the screeplot (Figure 4), the first 2 components are sufficient to explain the variability of the data set. Thus, this PC1-PC2 plan incorporates enough significant information to highlight and explain the general trends of the studied system.
Table 5. Physicochemical variables loadings on principal components.
Figure 4. Principal components analysis (PCA) screeplot.
The component PC1, which has the highest variance (37.2%), is strongly associated with major chemical elements (Ca2+, Mg2+, Na+, , Cl−). Correlations with the other elements (F− and ) also remain significant. It is noted that the physical parameters (temperature, pH, depth) are also well associated with this axis. and K+ are poorly represented on this axis. The component PC1 thus expresses the mineralization of the water in the study region.
The second component PC2 is positively associated with F− (0.71), and negatively with potassium (−0.86).
The analysis of the projection of the variables on the factor plane (PC1-PC2) reveals that the factorial axis PC1 gathers almost all the elements in the positive pole with the exception of potassium (K+) and . These two elements have no significant correlation with the set of variables represented in the correlation matrix. The axis CP2 is negatively correlated with the variables; Ca2+, Cl−, , pH, T˚C, and K+, the positive pole is represented by the variables; , Na+, Mg2+, F−.
The spatial analysis of the variables of the PC1-PC2 plane makes it possible to classify, according to the affinities of the ions, two groups of variables. The first group is determined by , Na+ and F−. The second group contains the variables pH, T˚C, and Cl−. The Ca2+, Mg2+ and Depth variables are intermediate to these 2 groups.
Figure 5. Principal Components Analysis (PCA). Hydrochemical variables loading plot on PC1 vs. PC2.
Figure 6. Plot of groundwater samples on PC1 vs. PC2.
These links between the chemical elements are in agreement with the chemical facies derived from the Piper diagram, which are , , and . This shows the acid hydrolysis of the silicate minerals present in the aquifer system and at the origin of the different ions. This cluster reflects a mineralization of groundwater due mostly to the water-rock interaction.
Nitrate and potassium do not correlate with the other elements, which presuppose that they are dissolved by different mechanisms. Nitrates have a superficial origin and are evidence of anthropogenic pollution.
The projection of the samples on the plane PC1-PC2 makes it possible to differentiate 2 very distinct clusters according to their mineralization. Waters flowing through the fractured reservoir are more mineralized than waters in the alterites. The latter also have a lower temperature. It is interesting to note that waters in the alterites remain undifferentiated on both sites (Am Zoer and Guereda-Iriba).
On the other hand, certain discrimination can be observed on the deep waters between the two sites. Accordingly, it can be noted that most of the deep wells at Am Zoer site have warmer, more basic waters, differenciated by and Cl− ions. Guereda-Iriba waters are more associated with Na+ and elements, and especially with F−.
The alkaline earths (Ca2+, Mg2+) are present in both sites and do not discriminate between these two sites.
5.4. Groundwater Types
The analysis of the data from the 31 water samples shows the order of abundance of major ions as follows: calcium (51%) and bicarbonate (80%) are respectively the main cation and anion in the groundwater. The relative abundance in descending order for cations in groundwater is Ca2+ > Na+ > Mg2+ > K+ while that of the anions is . Since nitrates are not significantly correlated with any other element (indicating anthropogenic origin), they ae not considered in the Piper diagram (Figure 7).
The Piper diagram shows that the dominant anion is at both sites. For the waters of most deep wells at Am Zoer site, the dominant cation is Ca2+. Waters from dug wells and some boreholes at Guereda-Iriba, are marked by the poles Ca2+ and Mg2+. A single well (F23) at Guereda-Iriba, is characterized by the Na+ pole.
The Piper diagram thus shows the following hydrochemical types:
- bicarbonate-calcium type secondary bicarbonate-sodium and bicarbonate- magnesium
- bicarbonate calcium type secondarily magnesium
- bicarbonate sodium type secondary bicarbonate-calcium and magnesium
These hydrochemical types are also represented on a semi-logarithmic Schoeller diagram  (Figure 8) on 10 selected samples. This diagram also shows that the Ca2+ and ions are preponderant in groundwater chemical composition and reveals more clearly that Cl− and ions play a minor role in the mineralization of the basement aquifer groundwater.
6. Hydrogeochemical Processes
6.1. Scatter Plots
Interpretation of ionic relations in the form of scatter plots can greatly help to understand the origin of solutes and the various hydrogeochemical processes involved in the groundwater chemistry.
Figure 7. Trilinear Piper plot of the groundwater samples.
Figure 8. Plot on the Schoeller semi-logarithmic diagram of some selected groundwater samples.
shows that in all samples these elements are very weakly correlated. All samples display an enrichment with Na+ compared to Cl−, indicating that these elements are not derived from the same source. The Na+/Cl− molar ratio is higher than 1 (1.2 < Na+/Cl− < 13.8), reflecting that silicate weathering can be the main contributor for Na+  . For example albite weathering contributes Na+ (see equation 01).
The relation of Na+ associated to base exchanges process is discussed below. The presence of Cl− can be attributed to another source for this ion than the dissolution of halite and could be associated, at least partly, to anthropogenic sources (agricultural fertilizers).
The source for alkaline earth elements (Ca2+ and Mg2+) may be dissolution of silicate minerals (e.g., plagioclase feldspar, chlorite, or biotite), carbonates (dolomite or calcite), gypsum, and/or cation exchange of Na+ for Ca2+ and Mg2+ on clay minerals. vs. Ca2+ (Figure 9(b)) displays a significant linear trend (R = 0.75). Samples from Am Zoer site display an excess of Ca2+ while those from Guereda-Iriba plot close to 1:1 line. An essential point is that mineralization processes in the study area, are not homogeneous but display a noteworthy spatial variability.
The scatter plot vs. Mg2+ (Figure 9(c)) is characterized by a significant correlation (R = 0.68), but also shows a deficit in Mg2+. These relations indicate that dissolution of carbonates (dolomite or calcite) is not the main source for these elements. Other sources for these ions are discussed below.
In crystalline rocks, Ca2+ is found in the feldspar network, associated with sodium and potassium, while Mg2+ is associated with iron in ferromagnesian materials (micas, amphiboles, pyroxenes, peridots). Thus excess of Ca2+ in granitic rocks (Figure 9(d)) may be caused by means of hydrolysis of plagioclase feldspars, as anorthite (Equation (2)):
Figure 9. Scatter plots between chemical elements.
Weathering of ferromagnesian minerals such as biotite may be responsible for the release of Mg2+ in groundwater in basement rocks (reaction 3):
Additional geochemical processes such as the Ca2+ ↔ Mg2+ exchanges during precipitation reactions or calcite recrystallization  may contribute to excess of Ca2+ compared to Mg2+ (Figure 9(d)). Figure 9(e) shows the scatter plot vs. Ca2+. This plot displays a large deficit in . The same observation is reached when comparing vs. Mg2+. The absence of any significant correlation indicates that the main source of these elements is not the dissolution of gypsum or any Mg-sulfate mineral. Concentrations of K+ are low in Wadi-Fira groundwater (average 3.5 mg/l), which can be expected as there are no important mineral sources of K+ in these aquifers.
6.2. Mineral Saturation Indices
Saturation indices of minerals indicate to which extent water mineralization is controlled by equilibrium with solid phases  . When this indice is negative, the water is undersaturated with respect to the concerned mineral. A positive SI means that water is oversaturated regarding the mineral. When the SI value is close to 0, water is in equilibrium with the mineral.
The software PHREEQC  was used to determine mineral saturation indices using the specific ionic concentrations in the water.
The saturation index of a mineral is obtained from Equation (5)  .
where IAP is the ion activity product of the dissociated chemical species in solution, Kt is the equilibrium solubility product for the chemical involved at the sample temperature.
The saturation indices of the following minerals in all water samples were evaluated: calcite, dolomite, gypsum and halite. The average values of the SI for each mineral are given in Table 6.
The saturation indices regarding calcite and dolomite have high and negative values, expressing that the groundwater is undersaturated with these minerals.
The evaporitic minerals (gypsum, halite) also have negative SI values. These results show that the minerals in question (calcite, dolomite, gypsum and halite) are not the main source of the groundwater mineralization. This is in agreement with the conclusions derived above.
6.3. Role of Evaporation
As the study is located in a semi-arid climatological context, it may expected that evaporation could also affect groundwater mineralization. The role of evapora
Table 6. Average Saturation Indices (SI) of principal minerals.
tion has been assessed using the Gibbs diagrams  .
Gibbs diagrams represent the ratios of (Na+ + K+)/(Na+ + K+ + Ca2+) and Cl−/(Cl− + ) as a function of TDS. They assume that evaporation increases salinity by increasing Na+ and Cl- contents in relation with increasing TDS. The diagram differentiates 3 areas highlighting various mechanisms: precipitation dominance, rocks dominance, evaporation dominance.
The plots of the hydrochemical data in Gibbs diagram are reported in Figure 10.
This figure shows that chemical weathering of rock-forming minerals is the most important functional source of dissolved elements in almost all groundwater samples. A few samples fall in the evaporation dominance area. These samples were collected from hand dug wells.
We may thus conclude that evaporation does not represent a major hydrochemical process in the study area.
6.4. Ions Exchanges Processes
Ions exchanges processes between the groundwater and its host rocks during its flow, can be analyzed using the plot of (Ca2+ + Mg2+) versus ( + ). In this scatter diagram, the points close to the 1:1line (Ca2+ + Mg2+ = + ) suggest that these ions have resulted from minerals dissolutions  .
When groundwater has been in contact with clayey materials able to release the interchangeable alkaline-earth ions easily, then Na+ contained in groundwater can exchange with these alkaline-earth ions (Ca2+ and Mg2+) sorbed on the exchangeable sites. In such cases, groundwater is enriched with Ca2+ and Mg2+ and depleted in Na+. This process is called reverse ion exchange. Direct ion exchange may also occur (exchange between adsorbed Na+ or K+ with Mg2+ or Ca2+ in the groundwater).
Figure 11 shows the plot of (Ca2+ and Mg2+) vs. ( + ) for all the samples. Most of the samples are close to the 1:1 line, demonstrating that the ion exchanges process is minor. A few samples display an excess of ( + ) suggesting the occurrence of direct ion exchange, but with limited extent.
The fact that ions exchanges processes between groundwater and the crystalline aquifer are relatively limited, is well in agreement with the previous findings.
7. Isotopes Analysis
The environment isotopes of oxygen d18O and hydrogen d2H are excellent tracers for determining the origin of groundwater and widely used in studying the natural water circulation and groundwater movement.
Figure 10. Gibbs diagrams.
Figure 11. Diagram (Ca + Mg) vs. (SO4 + HCO3).
The groundwater has δ18 O values in the range of −5.06‰ to −0.39‰ and δ²H from −30.18‰ to +6.0‰ (Table 7). The average values of δ18O and δ²H in the groundwater are −3.6‰ and −18.9‰ respectively.
The local meteoric water line (LMWL) is defined at the Geneina station (capital of West Darfur in Sudan). Its equation is:
The stable isotopes values are plotted in Figure 12. The global meteoric water line (GMWL)  is also given as reference. The GMWL is written as:
This diagram shows that the groundwater data are close to the GMWL meaning that the groundwater is in relation to modern rainfall.
This finding indicates that the groundwater recharge resource is mainly from precipitation and is weakly affected by evaporation.
Suc result is quite important since we can deduce that the renewal of groundwater is sufficiently rapid. Dating using C isotopes (13C, 14C) would help to determine the residence time of water in the aquifer system.
Table 7. Isotope composition of groundwater in the study area.
Figure 12. Stable isotopes diagram.
8. Summary and Conclusions
Groundwater is a vital resource for the Wadi-Fira region in eastern Chad. The water needs of this region of Chad have increased dramatically in recent years, given the increase in population and the development of agricultural and urban activities caused by an influx of refugees from neighboring Darfur.
The aquifer system exploited in this region is represented by the rocks of the crystalline basement. Groundwater is exploited using wells dug in the alterites and deep wells penetrating the fractured rocks. However, the available knowledge about the hydrogeology and the water quality of this system is still very low.
The objective of this work was to study the processes of groundwater mineralization and recharge mechanisms. Hydrochemical and isotopic analyses were carried out on water samples collected during two field campaigns organized in 2011 and 2015. The results of the analyses and their interpretations have shown that groundwater is on the whole fairly gentle. However, the mineralization of the deep waters of the fractured crystalline reservoir is more important than that of the waters of the alterites. The predominant process of mineralization is the interaction between water and rock and the solution of minerals. Other processes intervene (basic exchanges, precipitation, evaporation), but their impact is moderate. The isotopic data have shown that the groundwater of this system is related to modern precipitation. Recharging the system is relatively fast. This study has allowed to significantly improve knowledge on this complex system and laid the bases to consider a more rational exploitation of the resource.