In Mauritania, 4/5 of the territory is subject to an arid climate, where groundwater resources are the main sources of drinking water supply for populations and livestock.
The study area is located in southeast Mauritania, and covers the southern part of the Taoudenni Basin in Mauritania. It is characterized by a high concentration of livestock (49%)  , and a high spatio-temporal variability in rainfall and increasing population pressure on water resources, which cause water shortages and decrease in boreholes yield.
The implementation of sustainable water resources management plan is essential in this context which requires a better knowledge of the system functioning. The present work aims to investigate the hydrodynamic and hydrochemical characteristics of Aioune sandstone aquifer, the most important generalized aquifer in the Taoudenni transboundary basin in Mauritania, and to identify and quantify the recharge at this Aioune sandstone aquifer. Hydrodynamic, hydrochemical and isotopic data used in this study are powerful tools to identify the origin, the age of the water, mineralization processes as well as location of recharge areas in order to establish measures for sustainability of these valuable water resources  .
2. Study Area Description
The study area called “Afolle” is located in the southern part of the Taoudenni Transboundary Basin in Mauritania. It has a triangular shape and covers an area of 21,000 km2. In the center of the basin, culminates a cliff of 350 m IGN which is strongly affected by weathering and its borders are covered in the North, East and West by recent sedimentary layers (Figure 1).
The climate in the study area is generally Sahelian. Rainfall is irregular, and occurs mainly between July and October, with an average of 200 mm at Aioune station (Figure 2(a)). 70% of annual rainfall events have an intensity between 10 - 30 mm  . The series of rainfall data (1939-2006) is characterized by a deficit cycle (1969-1990), followed by slightly increase relative to annual mean (Figure 2(b)). The average temperature is 30˚C, with peaks of heat in May and June. Mean annual PICHE evaporation at Aioune zone is 5.32 m  .
Apart the Senegal River, there is no perennial flow in Mauritania. In the study area, the hydrographic network is more active in the center, with a coefficient of flow between 1.5% - 10%  .
The tributaries of the Senegal River in the study area are mainly the “Kolimbine” Wadi to the east of Afolle, and to the West the Wadi “Karakoro”. These Wadis are currently cut by sandy cover (endoreism), which is often present in the study area  (Figure 3).
2.4. Geology and Hydrogeological Context
Several French geologists, studied this zone, from the 19th century to the 1960s (H. HUBERT, R. HUDEAU, R. FURON and C. BENSE) and had identified an anticlinal structure inferred by tectonics, with the North, East and West edges plunging under recent sedimentation layers (Figure 1). The Afolle relief is characterized by a succession of tabular formations affected by numerous diaclases. The oldest sedimentary formations which outcrop are the sandstones of the
Figure 1. Geology, section cross AB and localization map of the Aioune aquifer area.
Figure 2. (a) Pluviometric and evaporation average per month at aioune station (1970-2009); (b) Pluviometric variation average per year to the Aioune station (1970-2009).
Figure 3. Piezometric head map of aioune groundwater (Auguste 2013).
infracambrian series. This sedimentary series unconformably overlay the granitic basement of “Birrimien”. The infracambrian series begin at the base with basic compact clayey quartzite sandstones of 400 - 1000 m thickness  , covered by ferruginous sandstones which extents a large area and have oblique stratification, called “Aioune sandstones”. This unit is between 80 and 120 m in thickness, and its oblique stratification characteristic could be due to strong marine currents during the past  . At the top, the infracambrian series occur with the sandstone compact quartzite “brain sandstone”. The thickness of the infracambrian series is in the order of 2000 m  .
In the north, the previous series are covered in angular unconformity  , by the lower Cambrian outcrop series, which begin at the base with coarse sandstones, interspersed by a granite cobble clay conglomerate of so-called glacial origin, namely Tillite. The thickness of this series in the west is 210 m  . To the West and East, the upper Cambrian outcrop is characterized by the presence of the Pelite, jasper and dolomitic limestone minerals.
The study area is characterized by a large extension of the unconsolidated and interstratified sandstone of Aioune  . This formation becomes aquifer in area where are favorable conditions occur  , such as the primary and secondary permeabilities of the Aioune sandstone formation occur. The superficial part of this formation mostly exhibits these characteristics. The thickness of the one layer aquifer of “Grès d’Aioune” is variable, with an average of 15 m  . The static level varies between 2 and 35 m below ground surface. The average transmissivity is 6 × 10−4 m2/s. The porosity in the order of 5% and flow rates can reach 15 m3/h  , and the renewable reserve is estimated at 2.1 million m3/year  .
A total of 23 groundwater samples, two springs water samples, and two surface water samples, were collected from the Aioune sandstone aquifer, in December 2011 and September 2013. In addition rainwater samples covering only the rainy events of August 2013 were collected during the couese of this stady. Physico-chemical parameters (pH, conductivity, temperature) were measured on site. The determination of the concentration of the major elements (28 samples) was carried out by high performance liquid chromatography (HPLC). The stable isotopes were analyzed by the laser spectrometer (28 samples). The tritium content was performed by liquid scintillator after electrolytic enrichment (28 samples). Only ten samples were analyzed for 14C and 13C. The radiocarbon analyses were performed by accelerator mass spectrometry (AMS). All analyses were carried out at the Laboratory of Radio-Analyzes and Environment (LRAE) of ENIS of Sfax (Tunisia) except for 14C and 13C, carried out at the isotopic laboratory of the University of Groningen (Holland). The analytical errors on the measurements are ±0.2‰ for oxygen-18, ±0.6‰ for deuterium, for Tritium ±0.3‰ and carbon-13 (Table 1).
Only the chemical balance data of ±5% were considered.
The rainfall samples collected at the Aioune station are representative (((2H-(18O *8))-10): 0.7), located on the global meteorological water line, and are in the range of stable isotope values recorded at the nearby GNIP station (Bamako (Mali) and Niamey (Niger) (Table 2).
The correct Bamako GNIP data (low deviation from mean ((2H-(18O *8)-10)), and extrapolated from the Vienna GNIP series, were used in tritium tracing.
4. Results and Discussions
Head data measured from 16 piezometers were used to establish piezometric map of Aioune groundwater. The piezometric surface morphology of the aquifer follows the soil topography, characterizing a low depth water table (water table). Head values range between 238 and 143 m of altitude and the flow direction is generally from the cliff towards the plains (Figure 4). The hydraulic gradient calculated for this water table varies between 1 and 5‰, the effective velocity is 94 m/year (Figure 3).
The data collected in 2013, at the National Hydraulic Database (SIPPE2)  , were additionally used and extrapolated to evaluate the volume taken in 2006, which is estimated at 18,000 m3/d (represents 6.5 million m3/year).
4.2. Hydrochemical Characteristics
4.2.1. Surface Water and Rainfall
In December 2011 and August 2013, two samples of surface water were collected from the study area, and all the August 2011 rain events greater than 5 mm, at the Aioune station. The average temperatures are 23˚C and 14˚C, respectively for surface water and rainwater. The average pH is 8.2 and 7.1, respectively for surface water and rain. Salinity ranges from 318 to 755 mg/l for surface water and 18 mg/l for rainwater. The surface water and rain samples exhibit HCO3-Ca water type.
1) Salinity and chemical facies
The chemical and physicochemical data from the 2011 campaign samples were assessed to evaluate the chemical quality of the water in the study area and identify the geochemical processes controlling mineralization of the groundwater.
Indeed, depending on the sampling points, the waters were classified into three groups: wells, boreholes and springs (Figure 1).
The temperatures range from 33.1˚C (Borehole F7) to 27.2˚C (P8 Well), with average 30˚C, and show some increase in depth, except for El Ghouds borehole (F6).
The pH varies between 6.5 and 8.6, with an average of 7.5. The electric
Table 2. Stable isotopes of precipitation in the study area and near regional stations sample name.
Figure 4. (a) Major elements of the Aioune sandstone water; (b) Piper diagram of the analysed Aioune groundwater samples.
conductivity varies from 140 μs/cm to 3200 μs/cm, and the average is 451 μs/cm. The spatial distribution of the electric conductivity shows higher values in the shallow depth of the alluvial spreading zones and in particular near the major agglomerations (eg Gounguel, Aioune, etc.).
The general order of abundance of the cations in meq is Ca2+ > Na+ > Mg2+ > K+, with the exception of few water points where Na > Ca. For anions, the general order is > > Cl− > , with high concentrations of Nitrates, reaching 1.1 g/l, mainly in the shallow level, reflects the vulnerability of the water table to the pollution (Figure 4(a)).
The representation of the chemical data in the Piper diagram shows waters whose chemical facies evolve mainly between two poles, a HCO3-Ca (low mineralized) pole, and a Cl-Ca/Na (generally more mineralized) pole. Some other facies types occur such as, HCO3-Ca/Na, HCO3-Na/Ca, HCO3-Na and SO4-Na, reflecting concommittent other chemical processes (Figure 4(b)).
2) Geochemical processes
In order to determine the origin of the mineralization of the Aioune sandstone groundwaters samples, some correlations between the main major elements have been studied, notably with Cl which is considered as conservative.
The positive correlation between Ca/Cl vs. Cl (Figure 5(a)) generally indicates the contribution of this element to mineralization, and similarly for Na/Cl vs. Cl (Figure 5(c)), indicates the significance of the contribution of these elements. These elements generally result from the dissolution of carbonates and evaporites and the basic exchange phenomenon. The K and Mg elements do not show a clear trend (Figure 5(b) and Figure 5(d)).
The positive correlation between Ca, Mg, HCO3 and SO4 (Figures 5(e) and Figures 5(f)), probably reflects carbonates and gypsum dissolution phenomenon. Figure 5(g) shows a good correlation between Na, K and HCO3, up to high values (3 meq-l−1) of HCO3, where there is a stability of the HCO3 concentration, the increase in Na + K, can be explained by the achievement of the equilibrium of water with respect to the carbonates, therefore the presence of the phenomenon of precipitation-dissolution.
Figure 5(h) indicates the presence of both types of cationic exchange with type I clays through the release in water of Ca2+ and or Mg2+ ions against fixation Na+ ion (P 3), and type II, where the Na+ ions of the complex are released and are replaced by the Ca2+ (F8) ions  . This cationic exchange mainly affects the mineralized waters, which are close to the alluvial spreading (Gounguel Zone, Douerare etc.).
Saturation index (SI) vs. major elements (Figure 6) show waters mostly saturated to supersaturated in carbonates, which explain by the presence of the phenomena of precipitation/dissolution of carbonates (Figures 6(a)-(c)). The waters are generally under-saturated in evaporite, which explains by the presence of the phenomenon of dissolution of evaporites (Figures 6(d)-(f))  .
The high exposure of the aquifer to contamination is probably due to variation in the unsaturated zone thickness (<30 m), characterizing a water shallow.
Figure 5. (a) Cl vs. Ca/Cl; (b) Cl vs. K/Cl; (c) Cl vs. Na/Cl; (d) Cl vs. Mg/Cl; (e) (HCO3 + SO4) vs. (Ca + Mg); (f) SO4 vs. (Ca + Mg); (g) HCO3 vs. (Na + K); (h) ((Na + K) − Cl) vs. ((Ca + Mg) − HCO3 + SO4).
Figure 6. Correlation between the saturation indices and major elements of the Aioune sandstone water.
4.3. Isotopic Characterization
4.3.1. Surface Water and Rainfall
For rainwater, the mean value content are −5.7‰ and −36.3‰ V-SMOW, for δ 18O and δ 2H respectively, and are aligned on the Meteoric Water Line (MWL). For surface waters, the δ 18O contents vary between 23.3 and 29.7‰ V-SMOW, with an average of 26.5‰ and for the δ 2H, values range between 99.6 to 130.2‰ V-SMOW, with an average of 114.9‰. Plotted in the conventional δ 18O vs δ 2H, data scatter and deviate the meteoric line evidencing evaporation phenomenon and probably mixing.
1) Stable Isotopes and Tritium
The δ 18O contents vary between −5.7 and −2.0‰ V-SMOW, with an average of −4.6‰ and for δ 2H, values range between −42.6 to −20.4‰ V-SMOW, with an average of −35.3‰. This large range probably reflects the multiple phenomena present in the system.
The δ 18O/δ 2H correlation characterizes two groups, Group A, with indirect evaporated recharge, occuring mainly at the wells and springs. Secondly, Group B, reflects mixing water between current and older recharge, located mainly at the boreholes (Figure 7). A significant variation of the stable isotopes of water, is recorded between 2011 and 2013, that characterize a water table with a fast response time, and also reflect the state of a superficial aquifer.
The δ 18O-Depth Correlation shows a water table affected by evaporation mainly at shallow depths between 0 - 30 m.
Correlation Cl-δ 18O identifies the presence of the phenomenon of dissolution, and to a less a extent evaporation phenomenon. This latter is located mainly at the shallow level.
Tritium values vary between 5.2 TU and 0.1 TU, with an average of 1.2 TU. The spatial variation of the tritium content shows the presence of local recharge, localized at the spreading zones, probably due to the topographic and favorable lithological conditions to the convergence and the water infiltration (Figure 8).
The average tritium contents of 1.5, 1 and 0.7 TU, respectively for wells, springs and boreholes, characterize the stratification of the aquifer. For some
Figure 7. 2H versus 18O of Aioune Sandstone waters.
Figure 8. Tritium map of aioune aquifer (December 2011).
wells the recharge occurred after 1961  , while all the boreholes and springs were recharged before 1961, except the F1 borehole, which is located at Tintane City, which suffered very strong inundations in 2007, causing the evacuation of the city (Figure 9).
2) Carbon Isotopes
Correlation 13C-14C exhibits large variation in 13C, from −11.9 to −7.4, indicating the presence of several carbon sources (the atmospheric reservoir and the carbonate matrix). The Aioune sandstone aquifer is characterized by a deep water enriched in δ 13C, identifying an exchange with the carbonate matrix  ; and surface water depleted in δ 13C (Figure 10), indicating the presence of isotopic exchange with soil CO2 or organic matter  .
The 13C equilibrium, representative of soil CO2, defines an open to intermediate system aquifer.
For the Aioune sandstone aquifer, activities in carbon 14 vary from 60 to 104 pcm, and an average of 84 pcm, confirming the contribution of recent waters to recharge.
3) Groundwater Dating
The variation in radiocarbon contents makes it difficult to identify the initial activity A0. Several models were used to correct this activity taking into account the different reactions. These models evolve from the simplest (model dilution), to the more complex (exchange isotopic).
In the presence of isotopic exchange between the water and the two carbon tanks, identified previously, the F and G, IAEA and Mook models are the most appropriate since they give close values; but in the presence of an open to intermediate system, the exchange phenomena with the gaseous phase and matrix, the presence of a calcareous medium  , the IAEA model is most indicated,
Figure 9. Datation Tritium of Aioune sandstone groundwater (2011) (pluviometric station Bamako(Mali).
Figure 10. Relationship between activity carbone 14―carbone 13.
and gives an age closer to the piezometric age.
The age ranges from the present to 1060 years, which is in accordance with the Tritium data, which characterize recent recharge (wells) to mixing water (boreholes) (Table 3).
Quantification of the recharge was discussed using the chloride masse balance and the carbon 14 method.
The chloride masse balance consists of using the Allison and Hughes formula, 1978, for rainfall events greater than 5 mm sampled at the study area and groundwater samples from the recharge zones ( at TU > 1).
where P: the average height of rain in mm・year−1, R: the recharge in mm・year−1, Cl rain: the chloride content in precipitation in mg・L−1, Claq: the content of chloride in the sample in mg・L−1.
The 14C method uses the formula adapted to free groundwater  , which takes into account the depth of structures and the thickness collected.
with: R: Recharge rate, H: Total aquifer thickness, Z: Depth of the strained area, θ: Porosity, λ: Radioactive decay constant, A: Activity measured in carbon-14, A0: Activity initial carbon-14.
The computed average recharges are 0.4 and 5 mm/year, respectively for the 14C method and the chloride masse balance, which represents respective volumes of 0.4 to 5 million m3/year.
The volume of water extracted from 494 water points (175 Wells and 319 boreholes), from the Aioune sandstone aquifer, is in the order of 6.5 million m3/year by considering a rate of 1 m3/hr for each well and 5 m3/hr for borehole  .
Table 3. 14C Age of Aioune groundwater.
The comparison between the annual recharge and the annual pompung rate, evidences an imbalance in the water budget, which results in the recurrence of water shortages during summer, and makes it urgent to take appropriate measures to make sustainable and rational this water resourses.
This pioneering work in the Aioune sandstone one layer aquifer, located in the southern part of the Taoudenni Basin in Mauritania, is based on the interpretation of hydrodynamic (piezometry, porosity, etc.), hydrochemical (elements major) and isotopic (δ 18O, δ 2H, δ 3H, δ 13C and δ 14C) data to better understand the behavior of this aquifer.
Piezometrie configuration of the water table reveals effective velocity of 94 m/year and direction flow generally from the cliff to the plains, spreading areas. The chemical data have identified a layer whose chemical facies evolving between two poles, a calcium-magnesium bicarbonated pole (low mineralized), and a second pole chloride calcium-sodium. However, the shallow water table in part of aquifer reveals high vulnerability to pollution, especially at the shallow horizons (Douerare, Aioune). The dissolution of carbonate and evaporitic minerals is likely the main process controlling mineralization, and to a lesser extent cation exchange,
The isotopic tool allowed to identify waters classified in two groups. Group A, with indirect evaporated moderne recharge is located at the shallow level of the spraying areas; and Group B consists of mixing water, between old and modern recharge water, located deep within the water table. The sandstone Aioune aquifer is characterized by an open to intermediate system, whose water age ranges from the present day to 1060 years.
The calculated average recharge is 0.4 and 5 mm/year, respectively for the 14C method and the chloride masse balance.
The 14C method takes into account the depth of structures and the thickness collected, and the chloride masse balance consists of using the rainfall events greater than 5 mm and groundwater samples from the recharge area (at TU > 1).
Recharge compared to water pumping characterizes a water table at the limit of imbalance. In order to be able to provide a sustainable solution, it is necessary to exploit these results to arrive at a model of the water table.
The authors gratefully acknowledge the contributions of the staff members of National Center of Water Resources/Mauritania Hydraulic Ministry, for their help during field-work. We also thank the member team of Department of Geology, University Cheikh Anta Diop of Dakar, and technical staff of MAU 7003 IAEA, and also staff of the Laboratory of Radio-Analyses and Environment of the National Engineering School of Sfax (ENIS) for their constant help and assistance during laboratory analyses. Finally, we wish to thank the editors, and reviewers for their constructive criticism, which helped to improve the paper.