Hydrochemical and statistical studies of the groundwater salinization combined with MODPATH numerical model: case of the Sfax coastal aquifer,Southeast Tunisia

Shallow groundwater is one of the main water resources in the arid and semi-arid regions. However, it is threatened by not only the reduced rainfall and recharge capacity, but also the water table drawdown and seawater intrusion. Such factors could cause a deterioration of the water quality and consequently the loss of a valuable hydraulic resource. This study aimed to improve our knowledge on the groundwater chemical quality evolution of the Sfax shallow aquifer, one of the most vulnerable areas in Tunisia, by developing a geochemical study using statistical and numerical methods. Salinization was identified by factorial analysis, PCA, and hierarchical clustering analysis in addition to the numerical MODPATH model. These findings confirmed that the groundwater quality has deteriorated due to natural and anthropogenic processes with a different influence of mineralization factors. They also revealed the location of seawater intrusion by focusing on the most vulnerable areas which are Chaffar and Djbeniana. Methodologically, the use of MODPATH model for seawater intrusion determination might be considered as the backbone for future studies in Tunisian coastal aquifers. The numerical model supports the results obtained by the geochemical analysis. Both methods are valuable tools as they contribute to trend determinations, management, and recovery plans.