Contribution to the Hydrochemical Study of Groundwater from the Continental Terminal in Saloum ()
1. Introduction
Between the end of the 1960s and of the 1990s, West Africa was affected by transformations from a significant decrease in rainfall [1]. This decrease in rainfall is a result of climate change. In Senegal, new climate variations are perceived through precipitation and temperatures [2]. The decrease in rains and high temperatures in the southern area of Saloum lead to a decrease in surface runoff and freshwater inflows. The salinity of surface waters in this area is increasing. High salinities are observed (greater than 60 g/l) and may be (greater than 140 g/l) during the dry season, in the upstream part of the Saloum river estuary throughout the year [1]. Work carried out in the area [1] [3] [4], indicated significant values of salinity in the Saloum estuary, which may even exceed the salinity of seawater. As a result, surface waters are dominated by a high salt content except for temporary runoff in the rainy season. Groundwater is therefore the main source of water for the population. Thus, with the presence of brackish water in places, found in the Maastrichtian water table and high fluorine contents that may exceed WHO standards [1], the Continental Terminal (CT) water table is the major source of water for the inhabitants of the area. In this area of Saloum, the population’s water supply is mainly based on the exploitation of groundwater from the Continental Terminal aquifer [5]. Due to the decrease in rainfall, high temperatures and population growth in the area, groundwater recharge decreases, evaporation and pumping increases. Consequently, the Continental Terminal (CT) water table in its coastal area is exposed to saline intrusion. According to Fadili [6], high pumping in coastal areas decreases the potential of the existing aquifer, by penetrating marine waters, changing its quality, following an imbalance between fresh water and water from the sea. Today, the local population has given up operating wells because of their high chloride content (3195 mg/l) [7]. In many arid and semi-arid regions, the intrusion of salt water into coastal and continental aquifers is a widespread phenomenon that can sometimes irreversibly degrade the fate of these waters [3]. Several studies have been carried out in this southern zone of Saloum [8] [9] [10] [11] [12]. These studies focus on geology, hydrogeology, major and minor ion hydrogeochemistry, isotopes and hydrogeological modeling. They made it possible to identify the sources of salinity, the mechanisms and degree of salinization, the surface water-groundwater relationships in this southern area of Saloum [13] [14] [15] [16] [17]. These authors indicated that in addition to the appearance of salt water on the coastal part, the increase in the salinity of the waters of the Saloum River also contributes to the contamination of groundwater in this southern area of Saloum [8]. These studies aim to understand the Saloum system to allow its good management. The detection and monitoring of salinity are essential in the mechanisms of its intrusion into coastal aquifers, within the framework of sustainable management of water resources [3]. Our objective is to participate in the knowledge of the physicochemical quality of the water of the Continental Terminal (CT) groundwater and to contribute to understanding the origins and processes of the mineralization of these water sources.
2. Materials and Methods
2.1. Study Site
The study area (Figure 1) is located in central western Senegal, between latitudes 13˚35' and 14˚10'N, and between longitudes 15˚40' and 16˚30'W. This southern area of Saloum is limited to the north by the Saloum River, south by the border of Gambia, west by the Atlantic Ocean and east by the Kaffrine region. The climate in the area is Sudano-Sahelian, with the alternation of a wet season and a dry season. Altitudes in the area vary from 0 (western and northern part) to around 40 m (central, southern and eastern). The hydrographic network in the area is mainly made up of the Saloum river with its two branches (Diomboss and Bandiala), the Néma, the bolons, the Djikoye and the Baobolon.
Geologically, the study context is that of the Senegalese-Mauritanian basin. This basin is from Meso-Cenozoic to Quaternary with at its base a Precambrian to Paleozoic bedrock which rises to the surface on its eastern and south-eastern side [18].
For Lienou, two-thirds of the Senegalese-Mauritanian basin is made up of
Figure 1. Study area and location of measured and sampled points.
clayey, armored sands, generally variegated [19]. In the alternating phase of hot and humid climates, at the end and start of the Pliocene and Pleistocene respectively, a lateritic modification involves a deposit of armor that surpasses the Continental Terminal [3]. The Continental Terminal (CT) is made up of detrital terrain, characterized by significant ferralitic modifications accompanied by ferruginous deposits and armor-cladding, high variations of silica, neoformation of kaolinite [7]. Faye recognizes from stratigraphic correlations of borehole logs capturing the CT (Figure 2): a diversity of detrital formations, coarse sand, fine sand, clayey sand and sandy clay, irregular and alternated by variations of clay [3]. The Continental Terminal has a heterogeneous lithology, the deposits are detrital, clayey to gravels, and this is often encountered in continental dumping grounds, that is to say, continental sediments redistributed by marine sedimentation [8]. The sedimentation of the Continental Terminal is heterogeneous and consists mainly of fine and medium sand, clay, clayey sand, sandy clay, laterite, lateritic sand, silt, limestone, marl limestone and limestone-marly. The Continental aquifer Terminal occupies most of the territory of Senegal, in the southern zone (Sine-Saloum-Gambia), it is formed of intercalation of sands and clays, while in Ferlo it is sands, sandstones and clays [3]. This aquifer is contained in the sandy, sandy-clay and clay deposits of the Continental Terminal and the Quaternary. According to Diluca, its thickness increases from 20 m in Kaolack, to more than 100 m to the south [9]. To the west and north near the Saloum river, the thickness is 20 m, it becomes higher towards the center with 80 m in Keur Madiabel (thalweg zone) and 100 m in north of Keur Samba Gueye (thalweg zone) [3]. Based on a review of several studies carried out in the Sine-Saloum-Gambia zone [3] [7] [8] [9] [10] [15], and archives of the Direction of Water Resources Management and Planning (DGPRE), hydraulic conductivities vary between 9.0 × 10−7 and 1.0 × 10−1 m/s and transmissivity values (T) are between 1.8 × 10−4 and 7.0 m2/s. The storage coefficients (S) obtained from
Figure 2. Study area geological cross section. Source: Faye et al. [12].
pumping test and the coefficients obtained from drilling reports, in view of the drinking water distribution network in the zone are between 8.0 × 10−3 and 3.6 × 10−2 [9].
2.2. Sampling and Physicochemical Analysis
Temperature, pH and electrical conductivity (EC) were measured in situ at 46 water points (wells and boreholes) in February 2019 using a multiparameter probe. Ten (10) samples are taken in plastic bottles from wells and boreholes in April 2019. At these water points, the temperature, electrical conductivity (EC) and pH are measured in situ. The bottles are then stored in coolers for transport to the laboratory. The chemical analyzes of the elements were carried out in the hydrochemical laboratory of the Geology Department of Cheikh Anta Diop University of Dakar (UCAD) and in the sanitation laboratory of Polytechnic School of Thies. These analyses focused on magnesium (Mg2+), calcium (Ca2+), sodium (Na+), potassium (K+) ions, chlorides (Cl−), nitrates (
), bicarbonates (
), sulphates (
) and fluorides (F−). These ions analyses were carried out by ion chromatography (Dionex-AQUION). Analysis results of 14 water points (boreholes, wells and piezometers) were obtained from a database of Direction of Water Resources Management and Planning (DGPRE) in January 2018. These results relate to cations (Mg2+, Ca2+, Na+, K+) and anions (Cl−,
,
,
, F−).
The points were located with their geographic coordinates using GPS (Global Positioning System).
These coordinates were used to represent the spatial distribution of the water points sampled with ArcGIS software (Figure 1).
2.3. Data Processing
The methods applied to process the analysis results are based on Ion balance (BI) analysis, Piper diagram, Base Exchange indices (i.e.b.), Saturation Indices (SI), binary diagrams and multivariate statistical analysis.
2.3.1. Ion Balance Analysis
The study of ionic balance (BI) can help to verify the accuracy of test results. The results of analyses can be considered good if this ionic balance is less than 5%, for a value of the balance between 5% and 10%, the results can be maintained and when the balance is greater than 10%, the analyzes are to be rejected [20].
The ionic balance is given by the following Equation (1):
(1)
where
BI (Ionic Balance);
∑cations: sum of major cations in water in milliequivalents;
∑anions: sum of major anions in water in milliequivalents.
2.3.2. Piper’s Diagram
The Piper diagram allows a classification of water. This classification is based on a representation of the chemical facies of water on a triangle of cations (left), anions (right) and an associated diamond giving a synthesis of the overall facies. The approximate location of an analytical result on each of these two triangles gives the first clarification on the anionic and cationic dominance [21]. This Piper hydrochemical diagram is mainly used in the literature to display the chemical facies of a set of water samples [22] [23].
The Piper diagram is used with the help of Diagramme software for the processing of hydrochemical data.
2.3.3. Base Exchange Indice (i.e.b.)
Changes in the chemical quality of water during its movement in the aquifer can be highlighted by the base exchange index (i.e.b.) [24] [25]. The i.e.b. expresses the relationship between the ions exchanged and the ions of the same nature originally existing in water. This index gives the nature of the ionic exchanges between water and the surrounding land using Equation (2) [22] [25]:
(2)
where:
[Cl−]: concentration in meq/l of chloride ions;
([Na+] + [K+]): concentration in meq/l of the sum of sodium and potassium ions.
So:
- if i.e.b. is negative, the Ca2+ and Mg2+ ions of the water are replaced by the K+ and Na+ ions of the surrounding formations;
- if i.e.b. is positive, the Na+ and K+ ions of the water are substituted by the Mg2+ and Ca2+ ions of the surrounding lands;
- if i.e.b. = 0, we have equilibrium between the chemical constitution of the water and the constitution of the surrounding land.
2.3.4. Saturation Index (IS)
The influence of chemical elements in geochemical mineralization processes can be given by the state of saturation, determined by the saturation index. This saturation index provides information on the level of chemical equilibrium of water with the mineral in the aquifer matrix, in water-rock relationships, as a measure of the mechanism of dissolution and/or precipitation [26]. The saturation index is defined by the logarithm of the ionic activity ratio (PAI) to the equilibrium constant (K) (Equation (3)) [27]:
(3)
where:
PAI: ionic product of the solution;
K: equilibrium constant.
Three cases are possible:
IS < 0, the solution is undersaturated with respect to the mineral considered;
IS = 0, indicates saturation;
IS > 0, the solution is supersaturated.
Diagrams and PHREEQC Interactive software are used to calculate saturation indices with respect to minerals (calcite, aragonite, dolomite, anhydrite, gypsum and halite).
2.3.5. Binary Diagrams
The binary diagrams show on a 2D figure the variation of the concentrations of the different ions or the representation of ion concentrations as a function of the chlorides considered as conservative halogens and not impacted by redox phenomena and minerals with reduced solubility [28] [29]. These binary diagrams can be used to analyze mineralization from correlations between the concentrations of major elements. These linear correlations between some of the chemical variables also aid in the study of the origin of mineralization according to the type of relationship [1]. For Faye, binary diagrams make it possible to distinguish the geochemical mechanisms responsible for the chemical evolutions noted, in addition to specifying the types of water [10]. The chloride ion indicates the origin of the salt in the water, it is a mixing tracer absent in exchanges between water and rock and is not lost [29] [30] [31].
2.3.6. Multivariate Statistical Analysis
Multivariate statistical analysis applied in the processing of hydrochemical data is based on the techniques of Ascending Hierarchical Classification (AHC) and Principal Component Analysis (PCA). These analysis techniques reduce and classify the information resulting from hydrochemical data with the aim of explaining the process of water governance [22] [32] [33] [34]. They show similarities or dissimilarities by grouping chemical elements (variables) or individuals (water point) according to their chemical compositions.
2.3.7. Ascending Hierarchical Classification
The Ascending Hierarchical Classification (AHC) is an analysis assessing the similarity or difference between samples to give an indication of class association [35]. This analysis allows multiple data to be grouped and ordered into distinct classes on a graph to facilitate their interpretations. On this graph, the individuals are distributed in succession of nesting forming a classification tree, called a dendrogram. According to Ahmed, it is a method of hierarchical classification of parameters, widely used in earth sciences and in the processing of hydrochemical results [22]. The Ascending Hierarchical Classification (AHC) method is powerful for processing hydrochemical data and for establishing geochemical models [32] [33].
2.3.8. Principal Component Analysis
Principal Component Analysis (PCA) is a descriptive statistical technique that displays on a graph much more information forming a database [24] [36]. It thus makes it possible, from a graphical representation, to distribute the chemical elements on the basis of their contents in the sampled water by reducing the number of variables. PCA is a method which decreases variables by representing the point cloud in a two-dimensional subspace generated by factor axes [32] [37] [38] [39]; this analysis determines the essential elements making it possible to interpret the phenomena involved from correlation with the variables [40]. According to Ahmed, PCA is used to study the mechanisms governing the physico-chemical constitution of water and/or their origins, and it is a new descriptive method recommended by Farnham, Johannesson et al. [22] [41].
The Ascending Hierarchical Classification (CHA) and Principal Component Analysis (PCA) methods are processed using the XLSTAT software.
3. Results and Discussions
3.1. Physicochemical Characteristics of Ground Water
The temperature values vary between 23.20˚C and 32.00˚C with an average of 27.96˚C. The pH varies between 5.35 and 8.00, with an average value of 6.61 indicating acidic water overall. The electrical conductivity (EC) varies between 37.4 and 12,320 µS·cm−1 with an average of 729.3 µS·cm−1. Indeed the smallest values of EC are observed from the center to the South and South-East of the zone, while the higher electrical conductivities are noted in the West (near the Ocean) and North (near the Saloum river). Electrical conductivities increase as one approaches areas contaminated by salt water from the ocean and the Saloum River [3]. Based on WHO standards (2008), Dieng Ndao differentiates between 2 categories of water in the CT aquifer: fresh water (EC < 900 µs/cm) and those which are salty (EC > 900 µs/cm) [1].
The results of the analysis were subjected to an electric charge balance (ionic balance) between anions and cations to check their validity (Table 1). Samples which have an ionic balance greater than 10% in absolute value are discarded.
The concentrations of chlorides (Cl−) in water vary between 3.15 and 1152.60 mg/l with an average of 118.07 mg/l. High levels are measured near the Saloum river and the sea. These high values confirm possible contamination of the water table by the Saloum river and the sea [7]. The concentrations of potassium (K+) and Sodium (Na+) vary respectively between 0.33 and 18.79 mg/l with an average of 3.18 mg/l, and between 2.35 and 479.56 mg/l for an average of 51.78 mg/l. With regard to the calcium (Ca2+) and magnesium (Mg2+) contents, they vary respectively between 2.00 and 219.66 mg/l for an average of 36.66 mg/l, and between 1.20 and 93.36 mg/l with an average 10.81 mg/l. For bicarbonates (
), they are between 15.25 and 372.10 mg/l with an average value of 85.65 mg/l, while sulphates (
) vary between 0.51 and 55, 48 mg/l with an average of 6.15 mg/l. Nitrates (
) vary between 0.24 and 200.44 mg/l for an average value of 22.95 mg/l. Fluorides (F−) vary between 0.02 and 1.8 mg/l. Some of these
Table 1. Ionic balances of the samples.
aThe samples were taken from different types of water sources: Piezometer (Piézo), Borehole (F) and Well (Puits).
levels exceed the OMS standards of 50 mg/l (for nitrates) and 1.5 mg/l (for fluorides).
3.2. Hydrochemical Classification of Groundwater
An analysis by ascending hierarchical classification is carried out from the results of the analysis. This analysis showing the dendrograms (Figure 3(a) and Figure 3(b)) allowed us to classify the water points into three groups.
The representation of the concentrations of major elements on the Piper diagram (Figure 4) gives the main chemical facies of water.
Group 1 is represented by most of the structures sampled (Well6, Well8, Well2, Piezo.2, Piezo.6, Well7, Piezo.4, Well9, F1, Piezo.3, Piezo.5, F5, Well1, F8,
Figure 3. Ascending hierarchical classification dendograms.
F2, F6, Well5 and F7). This group is characterized by low mineralization with electrical conductivities varying between 48 and 653 µS/cm, for an average of 204.1 µS/cm. Group 1 does not have a particular geographical distribution. It is found in the center, North, South, East and West of the southern zone of Saloum. The distribution of the contents of the major elements is variable,
(between 15.25 and 91.5 mg/l), Cl− (between 3.15 and 171.30 mg/l), Ca2+ (between 2.00 and 43, 93 mg/l), Na+ (between 2.35 and 75.93 mg/l), Mg2+ (between 1.20 and 18.27 mg/l), K+ (between 0.33 and 4.55 mg/l),
(between 0.51 and 13.59 mg/l) and
(between 0.24 and 149.51 mg/l).
The main chemical facies in water are:
● the calcium bicarbonate facies (Ca-HCO3) (Piezo.2, Piezo.6, Well7, Well9, F1, F2 and F6);
● sodium bicarbonate facies (Na-HCO3) (Well6, Piezo. 4, Wells1, F8 and F7);
● he calcium chloride (Ca-Cl) and sodium (Na-Cl) facies (Well8, Piezo.5, Well2, Piezo.3, F5 and Well5).
Calcium bicarbonate and sodium bicarbonate facies dominate compared to the other facies. According to Dieng Ndao (2017), the low electrical conductivities characterizing freshwater of the Ca-Cl facies type show significant areas of groundwater recharge at the terminal continental [1]. Despite its diversity, this group consists mainly of sites in the area with fresh water. Generally, the water of the groundwater is soft, but near the sea and the river, it is salty [10]. However, the high levels of nitrates in some works suggest contamination due to anthropogenic origin. Nitrates show that groundwater is contaminated by recently infiltrated water [8]. The calcium chloride facies (Na-Cl) present in the center of the zone which is characterized by fresh water is thought to be due to brackish water. Pockets of brackish water are increasing in the aquifer, although they are currently not well demarcated [15].
Figure 4. Piper diagram of: (a) Water facies; (b) With groups of structures.
Group 2 is represented by the works (Piezo.1, Well3, F3 and Piezo.7). The electrical conductivity of this group varies between 296 and 656 µS/cm with an average of 469.5 µS/cm. Changes in the content of major ions give,
(131.15 - 305.00 mg/l), Ca2+ (23.00 - 75.00 mg/l), Cl− (19.31 - 50.40 mg/l), Na+ (5.74 - 46.29 mg/l), Mg2+ (4.89 - 22.48 mg/l), K+ (0.6 - 15.10 mg/l),
(1.02 - 9.12 mg/l) and
(0.32 - 8.82 mg/l). The facies of this group are mainly of the calcium bicarbonate (east, center and north of the area) and sodium bicarbonate (center-west) types.
Group 3 is characterized by an electrical conductivity varying between 2900 and 4100 µS/cm for an average of 3500 µS/cm. This group is represented by well 4 and F4 (Ndoffane forest and Passy). The waters of these localities correspond to the salty waters of the southern zone of Saloum. The concentrations of chloride (Cl−) and sodium (Na+) ions in these sites vary respectively between 960.80 and 1152.60 mg/l, and between 320.09 and 479.56 mg/l. The contents of anions and cations vary as follows: Cl− >
>
>
> F− and Na+ > Ca2+ > Mg2+ > K+. The facies of this group is a sodium chloride (Na-Cl) facies. These types of water show contamination by water from the sea and/or the Saloum River [1] [3] [7] [11] [15]. The high concentration of nitrate (
) in Ndoffane forest (200.44 mg/l) far exceeds the WHO standard (50 mg/l). This group 3 would highlight the contribution of human activities in the mineralization of the water of the terminal continental groundwater. The intensity of domestic and pastoral activities near wells without slabs and protective enclosures would promote percolation of rainwater to the water table [15].
3.3. Water Mineralization
3.3.1. Base Exchanges Indices
The base exchange indices (i.e.b.) of the sampled waters vary between −3.47 and 0.73 (Figure 5). Most of these water points have positive values with 58.3% of the indices. These positive clues indicate that Mg2+ and Ca2+ ions replace Na+
Figure 5. Change in base exchange index.
and K+ ions in water. In contrast, negative core trade indices represent 41.7%. These negative indices indicate that the surrounding terrain releases Mg2+ and Ca2+ to fix the Na+ and K+ ions. The dominance of kaolinite-type clay is remarkable, over this cation exchange mechanism through clay minerals [11].
3.3.2. Water Saturation Indices
Box whisker (Figure 6) are represented from the calculation of the saturation indices. These boxes show that the waters are undersaturated with respect to all the minerals (calcite, aragonite, dolomite, anhydrite, gypsum, halite). Carbonate minerals (calcite, aragonite, dolomite) tend to approach saturation, while evaporitic minerals (anhydrite, gypsum, halite) are far from the saturation line.
The dissolution of these minerals would therefore participate in the mineralization of the water.
3.4. Correlation Analysis of Physico-Chemical Parameters
3.4.1. Binary Diagrams
Figure 7 shows the cross diagrams of the concentrations of major elements (Na+, Ca2+, and Mg2+) with chlorides, (Ca2+,
) with Mg2+ and (
,
) with Ca2+. The distribution of points around the freshwater-saltwater mixing line can be adopted to determine other phenomena related to the mixing mechanism [31].
The representation of these major elements (Na+, Ca2+ and Mg2+) with (Cl−) indicates strong values of correlations. The correlation coefficients (r) are greater than 0.9 (Figures 7(a)-(c)). These correlation values would indicate a common source of these ions. This source could be marine and/or a dissolution of carbonate, evaporitic minerals. The high correlation between Na+ and Cl− would be due to the same origin of a dissolution of halite and/or the presence of salt water [8]. The halite saturation index, undersaturated with respect to the water samples (Figure 6) could corroborate this. The representation of Na+ as a function
Figure 6. Water saturation indices with respect to, calcite, aragonite, dolomite, gypsum, halite and anhydride.
Figure 7. Binary diagrams showing relationships between the contents: (a) Na+/Cl−, (b) Ca2+/Cl−, (c) Mg2+/Cl−, (d) Ca2+/Mg2+, (e)
/Ca2+, (f)
/Mg2+.
of Cl− (Figure 7(a)) (r = 0.98) indicates an important part of the points located on the mixing line and some placed above. The alignment of the points is made on the mixing line between two poles formed for a mixture of two waters, i.e. when the water is contaminated by salt water, on the sea water dilution line [42]. This representation shows a weak growth in sodium ion (Na+). Enrichment/depletion of Na+ in fresh water would result from exchange of bases ensuring the progression of facies from Ca-Cl or Na-Cl or Ca-HCO3 to Na-Cl or Na-HCO3 [10] [11]. The relationship between Ca2+ and Mg2+ (Figure 7(d)) shows a good correlation with r = 0.91. On this relation, the arrangement of the points with respect to the mixing line favors an enrichment in magnesium (Mg2+) compared to calcium (Ca2+). The movements of the freshwater and seawater compositions are ensured by ion exchange reactions largely resulting from the enrichments and depletion of magnesium and calcium, respectively [22]. Thus, the representation of Ca2+ with Cl− (Figure 7(b)) (r = 0.92) would show this depletion of Ca2+ by a large part of the points placed at the bottom of the mixing line. The diagrams
with Ca2+ and
with Mg2+ (Figure 7(e) and Figure 7(f)) indicate correlation coefficients respectively r = 0.55 and r = 0.65.
For these representations, calcium and magnesium are better correlated with chlorides (Figure 7(b) and Figure 7(c)) (r = 0.92 and r = 0.95) than with bicarbonates. These results seem to confirm a common, marine origin of Ca2+ and Mg2+, especially in the northern and western parts of the area. However, calcium and magnesium would come from a dissolution of dolomite and/or calcite, since the waters are undersaturated with respect to these minerals (Figure 6). The degradation of minerals would lead, from a dissolution of calcite and/or feldspar to the presence of Ca2+ [14]. The enrichment in calcium relative to bicarbonate (Figure 7(e)) would be due to a dissolution of calcite in the presence of CO2. In the group of fresh waters that are not contaminated by salty waters, the waters are often of calcium bicarbonate facies (Ca-HCO3) [8]. This type of Ca-HCO3 facies occurs largely in the center, south and east of the southern Saloum area and is a characterization of freshwater. The relationship between Ca2+ and
(Figure 7(g)) (r = 0.71) and the saturation index (Figure 6) of undersaturated water compared to gypsum indicate that calcium and sulfate ions would be bound to the dissolution of gypsum formations.
3.4.2. Principal Component Analysis (PCA)
Principal Component Analysis (PCA) is applied to major elements (
, Cl−,
,
, Ca2+, Mg2+, Na+, K+), conductivity (EC) and pH. This analysis gives the results of Table 2, Table 3 and Figure 8. The correlation matrix (Table 2)
Table 2. Correlation matrix (Pearson (n)) of elements, EC and pH.
Table 3. Eigenvalues, percent of variance, cumulative eigenvalue, cumulative percent of variance for the factor analysis of hydrochemical data.
Figure 8. Distribution of physicochemical variables in the correlation circle according to the factorial plane F1 - F2.
shows a strong correlation of electrical conductivity with major elements except
and
. This correlation varies between 0.716 and 0.988. These strong correlation values express the contribution of major elements in electrical conductivity. The correlations of Cl− with
, Ca2+, Mg2+, Na+ and K+ are strong.
The correlation coefficient varies between 0.714 and 0.988. These correlations corroborate the common origin of these major elements which could be marine and/or by a dissolution of carbonate, evaporitic minerals.
From the eigenvalues (Table 3), two principal components (F1 and F2) are retained. These components represent 85.49% of the total variance. The choice of the components (1 and 2) is based on the criterion of Kaiser (1960), according to which, factors with eigenvalues greater than or equal to 1 are considered as possible sources of variance. The factor 1 (F1) expresses 68.69% of the total variance, and is represented positively by EC, Na+, K+, Mg2+, Ca2+, Cl and
. This factor appears as a factor of mineralization. It would reflect the dissolution of carbonate, evaporite formations and/or a saline intrusion. On the other hand, the second factorial axis (F2) represents 16.80% of the total variance. This factor is characterized positively by
and pH, and negatively by
. The opposition of these two anions in the space of variables, reflects two different origins of these elements. This difference could be explained by an anthropogenic origin (domestic and agricultural activities) of nitrates (
), and bicarbonates (
) would come from the dissolution of carbonates.
4. Conclusion
The pH varies between 5.35 and 8.00 with an average of 6.61 indicating acidic water overall. The electrical conductivity values vary between 37.4 and 12,320 µS·cm−1 with an average of 729.3 µS·cm−1. The higher conductivities are recorded around the ocean and the Saloum River. These strong conductivities indicate more mineralized waters in this part of the zone. On the other hand, lower conductivities are recorded in the center, south, east and south-east of the southern zone of Saloum. Geochemistry and multivariate statistical analysis show three groups of the sampled structures. These groups are mainly composed of the calcium bicarbonate facies, the sodium bicarbonate facies, the calcium chloride facies and the sodium chloride facies. Group 1 largely occupies the study area with the presence of these four main facies. Calcium and sodium bicarbonate waters, less mineralized, are mostly located in the center and east of the southern zone of Saloum. On the other hand, the calcium chloride facies and the sodium chloride facies are for the most part located in the northern and western parts. This group 1 presents good quality water, with the exception of a few sites which have nitrate levels higher than the WHO standard. Group 2 is present in the center, east, west and north of the southern zone of Saloum. This group consists mainly of the calcium bicarbonate facies and the sodium bicarbonate facies, and its mineralization is moderate. Group 3 is formed from the sodium chloride facies and is located to the north and west of the zone, the southern zone of Saloum (near the river and the ocean). The conductivity and chloride concentration of groundwater increases the closer you get to the Saloum River and the ocean. These increases suggest saline intrusion into the system. This intrusion would come from the salty waters of the river and the ocean. Strong correlations are obtained from the results of the binary diagrams and the PCA. These strong correlations confirm the saline intrusion into the aquifer system south of Saloum. The calculated base exchange indices suggest a base exchange phenomenon between the waters and the formations of the terminal continental water table. The saturation indices suggest an undersaturation of the water in relation to carbonate and evaporitic minerals. Consequently, these minerals would contribute to the mineralization of the groundwater of the aquifer of the terminal continental south of Saloum.
Acknowledgements
The authors thank the sanitary engineering laboratory of Polytechnic School of Thies and the laboratory of hydrochemistry of Cheikh Anta Diop University of Dakar (UCAD) for the chemical analyses carried out. They also thank the Directorate General of Water Resources Planning (DGPRE) for making chemical data available.