The water resources are exhausted by the increasing demand related to the population growth. They are also affected by climate circumstances, especially in arid and semi-arid regions. These areas are already undergoing noticeable shortages and low annual precipitation rate. This paper presents a numerical model of the Sfax shallow aquifer system that was developed by coupling the geographical information system tool ArcGIS 9.3 and ground water modeling system GMS6.5’s interface, ground water flow modeling MODFLOW 2000. Being in coastal city and having an arid climate with high consumption rates, this aquifer is undergoing a hydraulic stress situation. Therefore, the groundwater piezometric variations were calibrated for the period 2003–2013 and simulated based on two scenarios; first the constant and growing consumption and second the rainfall forecast as a result of climate change scenario released by the Tunisian Ministry of Agriculture and Water Resources and the German International Cooperation Agency “GIZ” using HadCM3 as a general circulation model. The piezometric simulations globally forecast a decrease that is about 0.5 m in 2020 and 1 m in 2050 locally the decrease is more pronounced in “Chaffar” and “Djbeniana” regions and that is more evident for the increasing consumption scenario. The two scenarios announce a quantitative degradation of the groundwater by the year 2050 with an alarming marine intrusion in “Djbeniana” region.