In situ measurement and geochemical data of groundwater samples.
Multivariate statistical techniques were applied to improve the understanding of the aquifers hydrodynamic and to identify the natural and anthropogenic processes that control the BouHafna and Haffouz groundwater quality. Some other parameters, such as sodium adsorption ratio (SAR), percent sodium (%Na), residual sodium carbonate (RSC), and permeability index (PI), were used to examine the suitability of groundwater for irrigation applications. Groundwater samples are classified into Ca-Mg-HCO3 and Ca-Mg-SO4 water-type. The statistical investigation permits to identify three different groups. The first group reflects the influence of water-rock interaction in relation with the dissolution of evaporatic minerals, the cation exchange process with phyllosilicates and the dedolomitization. The second and third groups, including the weakly mineralized groundwater samples, suggest, firstly, that the return flow of irrigation waters has a small, but not negligible contribution to the groundwater contamination, and secondly, the reduction of nitrate (NO3) to nitrogen gas (N2). Furthermore, it has been demonstrated that the majority of the groundwater samples are suitable for irrigation uses.
- Central Tunisia
- statistical techniques
- water-rock interaction
- cation exchange
- irrigation uses
In most parts of Central Tunisia, and particularly in Haffouz and BouHafna regions, groundwater has played a fundamental role in shaping the social and economic development. Due to rapid demographic growth and agricultural progress, exploitation of groundwater has increased dramatically to provide drinking water to rural community, support irrigation, and maintain ecosystems. During the past 30 years, the exponential increase in water abstraction from Haffouz and BouHafna aquifers has caused ground water depletion, a term often defined as long-term water-level declines, and an increase in salinity of groundwater pumped from wells situated mainly downstream. On the other hand, recent changes in agricultural land use and irrigation may result in groundwater contamination throughout agricultural fertilizers and pesticides applied to fields. Therefore, information about irrigation groundwater quality is critical to the understanding of necessary management changes for long-term productivity . Proper assessment of groundwater sustainability requires understanding and quantification of human effects on water resources using analysis, application of management practices, and revision.
Within this framework, a combined hydrogeological and hydro-chemical data were examined, using statistical methods, to determine (i) natural and anthropogenic processes that control the groundwater mineralization; (ii) the origin of different water bodies and their sources of recharge; and (iii) to assess the suitability of groundwater for agricultural purposes.
2. General features
The Haffouz study areas, which make up part of the Kairouan plain in Central Tunisia, cover about 1192 Km2 and lie between longitudes 39G 55′ and 39G 70′ north and latitudes 7G 88′ and 8G 33′ east (Figure 1). This region is bounded by Ouesselat and Jebil Mountains in the north, the Trozza Mountain in the south, and the plain of El Ala in the west. The altitude of the study area ranges from 200 m a.m.s.l. at Haffouz plain to 997 m at Trozza Mountain. The study area has a semi-arid climate with mild, wet winters and warm, dry summers . The average monthly temperature varies between a minimum of 10.15°C measured in January and a maximum of 33.07°C measured in August. It receives an annual rainfall ranging between 300 and 500 mm/year. The annual total evaporated exceeds 1720 mm Piche .
The drainage network consists of several nonperennial Wadis such as Mourra, 55 Zabbes, and MsilahWadis that drain toward the Merguellil Wadi, the most important Wadi of the Kairouan plain  (Figure 5).
2.1. Geological setting
Several geologic studies [4, 5, 6, 7, 8, 9, 10, 11] have focused on the stratigraphy, structure, and surficial and subsurface deposits characteristics in the study district. As shown in the geologic map and in the lithostratigraphic column, the geological formations, which outcrop in the study area, are represented by sedimentary series ranging from the Trias to the Quaternary (Figures 1 and 2). The Trias deposits, which outcrop in Cherichira Mountain, consist mainly of an alternation of gypsum, clay, and carbonate levels . The Cretaceous units that outcrop in the foot of Trozza Mountain are constituted of the Yprisian limestone of AlagAbiod and El Metlaoui formations. These formations are covered by the Lutetian limestone and clay of Cherahil and Souar formations, respectively. On the Lutetian deposits, the sandy and sandstone units of the Oligocene repose with an average thickness of about 350 m. These deposits belong to the Fortuna formation out crop in the central part of the study area. The Miocene deposits are represented by thick calcareous and clayey layers, which correspond to A in Grab and Mahmoud formations, respectively. These formations are overlaid by the sandy clay of the Beglia Formation and the clay Saouaf Formation. The Quaternary sediments that are represented by sand, clay, and conglomerate deposits occur in the plains of El Ala and Haffouz and in the Merguellil Wadi depression.
Structurally, the BouHafna basin is bordered to the east by the major fault of Ouesselat Mountain and corresponds to a syncline structure, which consists of the sandy deposits of Fortuna Formation. However, the Haffouz region, which is limited by the faults F1 in the east and F2 in the west, corresponds to a Graben structure filled by the sand, sandstone, sandy clay, and clay of Mio-plio-quaternary (Figure 3).
2.2. Hydrogeological setting
The BouHafna aquifer is lodged in the sand and sandstone deposits of the Oligocene. This aquifer, which is unconfined overall the basin, has a thickness varying largely from 50 to 300 m. The bedrock of this aquifer is constituted of the Lutetian marly deposits of Souar Formation (Figure 3). Nevertheless, the Haffouz unconfined aquifer is constituted of the Mio-Plio-82 Quaternary detrital deposits with a thickness of about 1000 m.
The piezometric map (Figure 4) shows that the general directions of groundwater flow are N-S, NW-SE, and W-E in the case of BouHafna aquifer, and they are oriented NW and NE in the case of Haffouz aquifer. This may indicate that the recharge of BouHafna aquifer occurs in the Western part of the basin and in the pediment of Ouesselat and Jebil Mountains, where the Oligocene sediments outcrop largely. For the Haffouz region, rainwater infiltrates and recharges the aquifer, locally, in the pediment of Ouesselat Mountain. The hydraulic gradient of BouHafna and Haffouz aquifers varies from 0.009 to 0.003, highlighting the importance of localized recharge in the pediment of mountains, in the northern part, and the linear recharge by the Merguellil Wadi, in the southern part of the region. This variation of hydraulic gradient can also be explained by the lateral variation of lithology caused by abundance of clayey layer, especially in the south of this aquifer.
3. Sampling and analytical methods
A total of 22 boreholes were sampled for chemical analysis, 15 samples collected from the BouHafna aquifer and 7 samples from the Haffouz aquifer. The physicochemical parameters such as Temperature (T°C), pH, and electrical conductivity (EC) were measured directly in situ using portable meters. Groundwater samples were analyzed for major ions in the laboratory of Radio-Analyses and Environment of the school of engineers of Sfax, Tunisia. Major cation concentrations were determined by Waters Ion chromatograph using IC-PakTM CM = D columns. Major anion concentrations were measured using a Metrohm ion chromatograph equipped with CI SUPER-SEP columns. The analysis of bicarbonate was undertaken by titration to the methyl orange endpoint.
4. Results and discussion
4.1. In situ measurements interpretation
The physicochemical parameters and overall hydrochemical data of Haffouz and BouHafna groundwater samples are presented in Table 1. Groundwater temperature values are relatively homogenous and vary from 19.9 to 25.3°C. The pH measurements range between 7.39 and 7.73 for Haffouz aquifer and between 7.01 and 7.74 for BouHafna aquifer. This may indicate a neutral to slightly alkaline pH but suitable for drinking and agricultural purposes. Electrical conductivity (EC) values, which vary in a wide range between 422 and 1623 μS/cm, are probably related to the signification variation of the TDS values that range from 203 to 1188 mg/l, respectively.
4.2. Water type
The data plotted in the Piper diagram  show that Haffouz and BouHafna groundwater samples have the same Ca-Mg-HCO3water-type, except for some samples collected from BouHafna aquifer, which are distinguished by Ca-Mg-SO4 water-type (Figure 5).
4.3. Correlation matrix
The correlation matrix of 12 parameters (Cl, NO3, SO4, HCO3, Na, K, Mg, Ca, TDS, EC, pH, and T°C) was computed in order to calculate the contribution degree of each hydrochemical parameters to the groundwater mineralization [15, 16] (Table 2).
Overall groundwater samples display positive and strong correlations between Cl, SO4, Na, Mg, and Ca versus TDS (>0.9), providing insight into the large contribution of these ions to BouHafna and Haffouz groundwater salinization. There is well positive correlation between Na and Cl (r = 0.84), indicating the same origin of these elements likely related to the halite dissolution. The strong and positive relationship between Ca and SO4 (r = 0.98) suggests that these ions derive from the same origin probably in relation with the dissolution of gypsum and/or anhydrite. The moderate correlation between NO3 and the majority of ions can be explained by an anthropogenic effect related to the pollution through return flow of irrigation water.
4.4. Principal component analysis
The principal component analysis (PCA), which is widely used in environmental studies, exhibits complex associations among several variables and individuals [17, 18]. Factors analysis was applied to the hydrochemical data set (Na, K, Mg, Ca, Cl, SO4, HCO3, NO3, EC, and TDS) of BouHafna and Haffouz aquifers in order to precisely specify the main processes controlling the groundwater mineralization (Figure 6a). The PCA approach has preserved only the first two factors, which represent 74.63% of total samples variance (62.77% for F1 and 11.86% for F2). In the variables space, the F1 factor displays strong positive loadings for Na, K, Mg, Ca, Cl, SO4, NO3, and TDS. The strong correlations for the referred major ions suggest that the groundwater mineralization is acquired through water-rock interaction processes. The positive loading for NO3 may reflect the influence of the return flow of irrigation water as a potential source of contamination related to the application of fertilizers. The F2 factor takes positive loadings for pH and negative loading for NO3 and HCO3. The inverse relationship between pH and NO3 with respect to F2 factor lends support to the implication of the denitrification process in the groundwater salinization.
The projection of the 22 samples on the same first factorial plan (F1, F2) permits to classify groundwater samples into three groups (Figure 6b). The first group, which lies in the positive side of the first factor (F1), is mainly represented by the BouHafna groundwater samples. This group is mainly influenced by natural processes of mineralization in relation with the dissolution of evaporatic minerals, dedolomitization, and cation exchange. These processes are summarized in the schematic model, which indicates that all groundwater samples display an undersaturation state with respect to gypsum and anhydrite, suggesting that dissolution of these minerals takes place there (Figure 7). Supplementary concentration of Ca2+deriving from gypsum dissolution can cause calcite precipitation . On the other hand, the decrease of bicarbonate concentration related to the precipitation of calcite causes the groundwater to be undersaturated with respect to dolomite and promotes the incongruent dissolution of this mineral known as dedolomitization [20, 21, 22].
The second group is mainly composed of BouHafna groundwater samples, which are strictly associated with intensive agricultural activities. This group, which plots in the negative side of the second factor (F2), is strongly associated with nitrate and highlights the considerable contribution of the return flow of irrigation waters to the contamination of groundwater. The third group comprises mainly the groundwater samples collected from the Haffouz aquifer. This group, which falls in the negative side of F1 and the positive side of F2, is extremely related to pH and negatively correlated with nitrate element. This position lends support to the predominance of denitrification processes, which corresponds to the biological reduction of nitrate (NO3) to nitrogen gas (N2).
4.5. Gibbs plot
Gibbs plot is used in the present investigation to confirm the significant role played by the natural hydrochemical processes already cited and their effect on groundwater quality . Plot of (Na + K)/(Na + K + Ca)
4.6. Suitability of groundwater for irrigation
Several parameters can be used to determine the suitability of groundwater for irrigation, that is, the electrical conductivity (EC), sodium adsorption ratio (SAR), and percent sodium.
The plot of analytical data on the Wilcox (1955) diagram shows that 83% of groundwater samples collected from the Haffouz aquifer belong to excellent category. All groundwater samples collected from the BouHafna aquifer, and only one sample of the Haffouz aquifer falls in the field of good to permissible (Figure 9). These very low to low SAR and low to medium salinity suggest that the studied groundwaters are suitable to moderately suitable for irrigation purposes without any threat of imposition of any hazard. Therefore, the application of these groundwaters in irrigation will be very advantageous as it will increase the agricultural yield.
The correlation of sodium adsorption ratio (SAR)
4.7. Permeability index (PI)
The permeability index (PI) parameter was used to assess the suitability of groundwater for irrigation. Indeed, the long-term irrigation with relatively enriched Na+, Ca2+, Mg2+, and HCO3− groundwater can affect soil permeability . Groundwater in the study area displays PI value ranging from 43 to 70% with an average value of 55% (PI greater than 25%) indicating that are good and suitable for irrigation purposes. Moreover, Doneen plot shows that groundwater samples in the study area fall in the fields of Class I and II, highlighting excellent to good permeability (Figure 11).
The hydrochemistry of BouHafna and Haffouz aquifer was investigated employing multivariate statistical approach in order to identify different processes that control groundwater mineralization. This investigation reveals the predominance of Ca-Mg-HCO3 and Ca-Mg-SO4 water-types. These water facies are derived mainly from water-rock interaction processes, i.e., the dissolution of halite, gypsum, the dedolomitization, and the cation exchange. On the other hand, return flow of irrigation water has resulted in elevated nitrate concentrations in groundwater especially in the agricultural zones, which are characterized by an excessive use of fertilizer. Thus, it is important to protect the aquifers against overexploitation and groundwater quality deterioration related to the evaporate dissolution and agricultural contamination. For these reasons, in the BouHafna and Haffouz regions where groundwater resources are under the great development stress and environmental pressure, some preventive measures should be taken. These are (1) control the exploitation groundwater; (2) the definition of special groundwater resources protection zones; (3) control the potential processes and sources of salinization; and (4) improvement of diffuse source groundwater pollution.