Modeling of the Groundwater Flow and Saltwater Intrusion in the Coastal Aquifer of Fum Al Wad , Province of Laayoun , Morocco

Saltwater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. The Fum Al Wad aquifer is located between Atlantic Ocean in the West and Laayoun in the East. This aquifer covers an area of 250 Km, and represents an essential water resource for Laayoun city and the periphery regions. It is heavily exploited for water supply, agriculture and industry. The freshwater-saltwater interface is affected by groundwater extraction by public supplies, irrigation wells, and domestic wells in the coastal of this aquifer. The position of the interface is controlled by several factors: these include precipitation, recharge rate, dryness, evapotranspiration, hydraulic conductivity and hydraulic head. Landward migration of the interface freshwater-saltwater often results in a significant decrease in the water resources available for coastal communities. The volume pumped by public for irrigation and the domestic usage in 2010-2011 is estimated 2.5 Mm/year, and in 2015 about 2.91 Mm/year only for domestic usages. The objectives of this work are to model the groundwater flow and saltwater intrusion in the coastal aquifer of Fum Al Wad, by SEAWAT-2000 program which coupled both the version of MODFLOW-2000 and MT3DMS. They are designed to simulate variable-density groundwater flow and solute transport in three dimensions. The model is calibrated for hydraulic conductivity, specific yield, porosity, and recharge rate as well as dispersivity coefficient. The result of simulation of the hydrodynamic model during the period of 1986-2015 has revealed a piezometric drawdown with 2.3 m approximately at the level of the pumping zone. Furthermore, this piezometric depression is caused by excessive pumping of the various uses of water and its corresponding that has reHow to cite this paper: El Mokhtar, M., Chibout, M., El Mansouri, B., Chao, J., Kili, M. and El Kanti, S.M. (2018) Modeling of the Groundwater Flow and Saltwater Intrusion in the Coastal Aquifer of Fum Al Wad, Province of Laayoun, Morocco. International Journal of Geosciences, 9, 71-92. https://doi.org/10.4236/ijg.2018.91005 Received: December 28, 2017 Accepted: January 27, 2018 Published: January 30, 2018 Copyright © 2018 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/


Introduction
Globally, the saltwater intrusion is manifested particularly in coastal and insular aquifers.They are sensitive to the salinization of the freshwater in natural and anthropogenic-influenced conditions [1].
Salinization of coastal groundwater is driven by factors of the over-exploitation, the climatic situation and the geometry of aquifers.Thus, the interface of transition between freshwater-saltwater develops spatially and temporally based on the development of the urban areas and the increase of the pumping activities and climate change.
The problem of the saltwater intrusion has become the main focus of many researchers.Hence, many works have been published by hydrogeologists, geophysicists and hydrochimistes authors [2] [3] [4].In this respect, they have conducted many studies to provide suggestive solutions.These specialists have set and worked on many hypotheses and perspectives [5] related to the issue being studied.Other recent activities similar to our saltwater intrusion theme use the same seawat-2000 modeling technique that is integrated into the GMS platform [6] [7] [8] [9].Lately, the climate change and human intervention have aggravated the situation of saltwater intrusion in the: unconfined coastal aquifers [10], confined aquifers [11] and also in estuaries.
The results of many experiments carried out by some laboratories made it possible to direct the science towards the optimization of the intensive exploitation of the coastal water [12].The desalination technique encourages the economic aspect and reinforces the optimization and perseverance of the water-table [13].
The piezometric, geophysical, hydrochemical and numerical approaches are certainly the most used to determine the state of coastal groundwater and the degree of their contamination by seawater or oceans water [14].
The numerical modeling of saltwater intrusion determines the space-time evolution of the variable-density-flow by using different simulation codes of the saltwater into the freshwater in the coastal groundwater [15].It serves as a good management method of groundwater resources.The hydrogeological model represents a decision support system that allows managers to conduct works in the future International Journal of Geosciences to improve the exploitation of this resource.
The objective of this work then is to realize the hydrodynamic and the solute transport modeling of the groundwater of Fum Al Wad aquifer.It also aims at presenting their methodology and the procedure for constructing complete conceptual and mathematical models.The construction of such models is the ultimate goal of this study.In fact, these models are used as essential tools for the planning and management for sustainable use of groundwater resources.This use of models is based on our belief that the physical reality can be represented by mathematical models, albeit with acceptable approximations.

The Governing Mechanism of the Interface between Freshwater and Saltwater
The presence of the saltwater into coastal aquifers is governed by factors such as: the geometry of the aquifer and the topography of its impermeable bottom.
Other factors can also contribute to the increase of salinity in freshwater as the geographical location of the water resource, the geology of the aquifer, climate and over-exploitation of groundwater.
The phenomenon of saltwater intrusion in coastal aquifers porous has been studied by Ghyben and Herzberg at the end of the 19th century [2].Saltwater is denser than the freshwater; it is naturally in coastal aquifers to a certain depth (z) depending on the hydraulic load of freshwater (h) [16] (Figure 1).
The freshwater and the saltwater are miscible in all proportions.The flow is a justifiable monophasic aspect of modeling by the dispersion theory [17].Indeed, these two fluids are in hydrostatic equilibrium in homogeneous and porous environment in a unique aquifer.In addition, these fluids' contact is supposed to Figure 1.Ground-water flow patterns and the freshwater-saltwater transition zone in an idealized coastal aquifer.A circulation of saltwater from the sea to the transition zone and then back to the sea is induced by mixing of freshwater and saltwater in the transition zone [19].International Journal of Geosciences be represented by a net interface [3].Equal pressures from the freshwater and saltwater are expressed by the following expression of the Ghyben-Herzberg: where (z) is the depth of the interface below the level of the sea; (h f ) represents the piezometric head of freshwater; ρ f and ρ s are densities of freshwater and saltwater [18].These last are equal respectively 1 and 1.025 kg/m 3 [3].The above equation can be simplified and thus becomes z = 40 h f .

Study Area
The Fum Al Wad aquifer covers an area which is about 250 km 2 ; it is situated between the East of Laayoun city (the up-rising of impermeable formations of the marls of the Upper Cretaceous), and the West of the Atlantic Ocean.However, the North and the South limits are not exact, but these boundaries were traced from wells' database [20] [21] as in (Figure 2 and Figure 3).
Figure 2. Extract of the geological map of Laayoun in the 1/100,000 scale [34] and the geological section AB oriented NNE -SSW.D: Sandhills; A: Silts and evaporates; q1: Alluvial terrace (named in the Sahara by "recent reg"); q2: Spreading loamy and stony average terrace, "old reg"; Q5: Slab made of carbonated sandstone and conglomerate; PVlh3: Carbonated whitish marl; PVlh2: Clear lumachelle within shells of gastropods, Brachiopods and conglomerate in the base; PVlh1: Redhead lumachelle mainly within oyster shells and gastropods; Mla2: whitish marl sandstone, sometimes with flint; Mla1: Lumachelle, sandy yellowish marls, sandstone, conglomerate, micro alternation of the whitish marls and the yellowish with a bench lumachellique; Com: Grey marl and gypsum.International Journal of Geosciences It is characterized by a coastal Saharian climate type that is warm in the winter and more or less sunny (cloudiness) with strong winds from the NNE.According to the Hassan I airport weather station of Laayoun city, this climate is softened by the ocean along the coast, while it becomes more and more hostile as one penetrates to the inside of the Saharian lands.
The temperature is somehow high during the year with maximum values in August about 33.7˚C except a few wet months.The area of the study receives minimum rainfall rate which do not exceed 10.78 mm as recorded in the previous month of January.In contrast, the maximum raining years are recorded during the periods 1989-1991 and 1995-1997 with a height of 155 mm in 1989 as annual rainfall rate.The Saquia Hamra Wadi (river) represents the main hydrographic element of the region of Laâyoune city.It is characterized by a discontinuous flow depending on the intensity of rainfall and exceptional floods.This river has created a kind of canyon which has distorted the surface floor; it has formed a slight slope towards the Atlantic Ocean.The morphology of this river became less obvious near the coast; this is due to sand-hills that have invaded partially its traces.In the mouth of this river, a form of delta has been created because of the swell of the Atlantic Ocean.The elevation of the natural terrain varies gradually from the International Journal of Geosciences Atlantic Ocean to Laayoun city, from 0 to 60 m.
The geology of the study area is part of the Meso-Cenozoic basin of Tarfaya-Laayoun-Dakhla which has developed during the Mesozoic and the Cenozoic [22] [23].Several authors have highlighted a general description of the structural and regional geology of this basin [24]- [33] from the Trias to the recent Quater-

The Mathematical Model of the Variable-Density of Groundwater Flow
To simulate the variable-density of groundwater flow in water-table, a SEAWAT-2000 numerical model was used which coupled the version of MODFLOW 2000 [35] and MT3DMS [36] [37].This model is designed to simulate three-dimensional of the variable-density of groundwater flow and multispecies transport.The variable density flow (VDF) process uses the familiar and well-established MODFLOW methodology to solve the variable-density of ground water flow equation.The MT3DMS which is part of SEAWAT that is referred to as the Integrated MT3DMS Transport (IMT) process can be used to solve the solute transport equation of the process of saltwater intrusion in coastal aquifers.
Guo and Langevin [38] derive the governing equation for variable-density ground-water flow, in terms of equivalent freshwater head [39], as: where: x, y, z are orthogonal coordinate axes, aligned with the principal directions of permeability; h f : equivalent freshwater head (m); Z: is the elevation at the measurement point (m); f ρ : is the density of freshwater (kg/m 3 ); s ρ : is the density of saltwater (kg/m 3 ); K (x,y,z) : is equivalent freshwater hydraulic conductivity International Journal of Geosciences (m/s); S s : is equivalent freshwater specific storage (1/m); t: time (s); e ω : is effec- tive porosity (dimensionless); C: is solute concentration (Kg/m 3 ); w ρ is fluid density source or sink water.W: is the volumetric flow rate of sources and sinks per unit volume of aquifer (1/s).
The conversions between h f (freshwater head) and h s (saltwater head) can be made using the following equations [38]: For a coupled variable-density flow and solute-transport simulation, a fluid density is assumed to be a function only of solute concentration; the effects of pressure and temperature on fluid density are not considered.A linear equation of state is used to represent fluid density as a function of solute concentration [38]: where: f ρ : is the density of freshwater (kg/m 3 ); s ρ : is the density of saltwater (kg/m 3 ); C: is the solute concentration (Kg/m 3 );

Geographical Situation of the Control Piezometric and Pumping Wells
The geographical situation of the pumping wells and the control piezometers are illustrated by the following Figure 3.The used coordination system of the all maps of this work is: WGS84, UTM zone 28 Northing.

Geometry of the Coastal of Fum Al Wad Aquifer
The geometry of the aquifer is built using the stratigraphic model which works by the borhole datas of the GMS (Groundwater Modeling System) software [15].
The distribution of these borhole datas is based on the choice of the grids of wells covering the study area (Figure 4).This geometry is retained in the modeling of the saline intrusion phase.
Lithological analysis of the three-dimensional model of the Fum Al Wad aquifer has allowed viewing seven hydrogeological formations which overlap with each other above the impermeable bottom formed by the gray marls of the International Journal of Geosciences

The Discretization and Boundary Conditions of Fum Al Wad Aquifer
The study area of the Fum Al Wad aquifer is subdivided into eight layers each with a thickness of 20 meters, except the first layer which has a varying size from 22 to 83 meters in Eastern zone (Figure 5).These three points are clearly illustrated below as in Figure 7.

Piezometry Calculated by the Steady-State Model in 1986
The mathematical model was established in the piezometric reference-state in International Journal of Geosciences The analysis of the piezometric map (Figure 6) below can give the following information:  The hydraulic gradient: it is highly manifested in the East and marked a maximum value which is about 0.4%; this is due to the groundwater recharge by the lateral exchange with other water-tables of Laayoun city in the East.
However, at the level of the ground-water pumping in the center of Fum Al Wad aquifer, it became lower and marked a minimum value which is about 0.004%.This indicates that the hydraulic gradient is influenced by the intensive pumping wells by the different cultivators and water supplies. The direction of the flow: it is coming from the East towards the West of the Fum Al Wad water-table. Piezometry: it varies from 0 meter from the Atlantic Ocean side (West) to 15 meters from Laayoun City side (East).
The permeability map is obtained by creating six (06) areas with hydraulic conductivity values ranging from 2 × 10 -4 m/s to 9.8 × 10 −4 m/s.The analysis of the sensitivity of the hydraulic conductivity parameter allowed selecting the most sensitive areas in the present model.Thus, the model is sensitive in zones 4, 5 and 6 while it is less sensitive in the zones 1, 2 and 3 as it is illustrated in Figure 7 which also shows the distribution of permeability per area.

Diagram of the steady-state calibration model of 1986
After several tests of the calibration process in the steady-state regime, the piezometric values simulated by the model and those observed in the land at the level of 40 controlled piezometric points (Figure 3) are almost diagonally aligned in Figure 8.International Journal of Geosciences

Piezometry of Fum Al Wad Obtained by the Transient
Model in 2015  9).This depression is displayed at the zone of the pumping wells in the center of Fum Al Wad water-table and it reflected the deficit flow budget of it.This is due to the overexploitation of this aquifer compared to the 1986 piezometric data.
The calibration of the transient model is performed with six controlling piezometers of Fum Al Wad piezometry since they contain a sufficient history of the piezometric head for the period 1986-2015.Thus, we have used these piezometers to calibrate the second hydrodynamic parameter (coefficient of storage); its value is selected carefully for better calibration which is about 10%.
The results of this calibration are illustrated by six wells as in Figures 10-15.
They show the evolution of the observed and simulated piezometric head by the transient model.

2) Flow budget of the Fum Al Wad water-table
The flow budget of the Fum Al Wad water-table has been obtained by simula-       3) The predicted results of the piezometric head for 2020 and 2030 In order to perform the predicted piezometric simulations results of 2020 and 2030, we have used extracted volumes that refer to 2015, namely the ones utilized for water-supply.Furthermore, the data of the other volumes, related to agriculture and industry, are not used since they are absent from the available data.
The predicted piezometric head vision at the horizon of 2020 shows a piezometric depression which is about 2.4 m compared to 1986 piezometric data.In 2030, the Fum Al Wad water-table will experience a piezometric head decrease simulated by the prediction transient model which is about 2.5 m in comparison to 1986 piezometric data.
The following Figure 16 represents the predicted piezometric heads.

Simulation of the Fum Al Wad Saltwater Intrusion in 2015
The saltwater intrusion modeling of the Fum Al Wad water-table has elaborated  retained is 30%; the density ratio is 1.5 (dimensionless).The Figure 17  Wad aquifer with 4.3 km approximately, as it is illustrated in section B in Figure 16.
In order to perform the predicted salt-wedge simulations results of 2020 and 2030, we have used extracted volumes that refer to 2015, namely the ones utilized for water-supply.Furthermore, the data of the other volumes, related to agriculture and industry, are not used since they are absent from the available data.
The two predicted simulations of the years 2020 and 2030 are obtained by simulating saltwater intrusion in freshwater by the transient model.The Figure 16 shows the advancement of saltwater towards the pumping wells as it is illustrated in sections A, B and C in the West-East direction.
The results of saltwater-intrusion modeling found in 2015 in the Fum Al Wad aquifer has highlighted that the pumping wells located in less than 4.3 km from the Atlantic Ocean were contaminated by sea-water (see Figure 17 below).
Furthermore, the future modeling of saltwater intrusion has highlighted that the salt wedge has immigrated into the freshwater in pumping zone with 4.4 km approximately in 2020 and 4.7 km in 2030.This indicates that the pumping wells will be partially contaminated.The following Figure 18 shows the state of salt

Discussion
The conceptual model consists of the construction of the geometry of the Fum Al Wad aquifer, which is based on the use of a sufficient history of the lithological data of the boreholes, wells, piezometers and other geophysical (geoelectric, seismic) data.The initialization of the hydrogeological model of this aquifer was done by defining the boundary conditions: the null flux in North and South limits, the imposed flow around 0.18 m 3 /s in the East and the imposed potential H = 0 m in West (Atlantic Ocean).Furthermore, the initial steady-state model was calibrated in 1986 using 40 observation piezometers.Thus, the main calibrated parameters are: the hydraulic conductivity, the effective porosity and the source and sink recharge rate (the Eastern limit).
The calibration of the transient model was done during a period of 29 years (1986-2015).The calibrated main parameter of the storage coefficient is done using 6 controlling piezometers with a complete piezometric history covering the entire simulation period of the transient model.
In addition, saltwater intrusion model of Fum Al Wad Aquifer was initiated by the same boundary conditions of the hydrodynamic model, except the western boundary which was defined as 35g/l of salinity concentration.
The model is calibrated with salinity data of the 5 salinometers (named: 737/120, 815/120, 826/120, 1151/120 and 1156/120) which have a sufficient salinity history .They are located in the pumping wells zone of Fum Al Wad.However, the main hydrodispersive parameters (longitudinal and vertical dispersivity, molecular diffusion) are calibrated in a period of 29 years (1986-2015) taking into account the previous results obtained from hydrodynamic modeling in the steady-state and transient model.
In this work, the form of the freshwater-saltwater interface modeled by the SEAWAT-2000 program confirms the result of the geophysical study which explains the presence of the salt-wedge phenomenon in the pumping wells zone of Fum Al Wad Aquifer [41].However, the salt-wedge moves diagonally from the West to the East with a gradual decrease in salinity concentration (35 g/l) from the coast to the continent (pumping wells side).
The present study shows that this interface is not in natural equilibrium since its temporal-spatio evolution reflects advancement towards the continent; this occurs in parallel with the increase of the extraction volumes from the water-table and the decrease of the piezometric head.This situation is aggravated at the level of the extraction zone of the water-supply (National Office of Drinking Water) where the salinization goes up gradually at the local level.
The predicted simulations have shown that piezometry is always in a dynamic state in the extraction zone.This reflects that the piezometric head will decrease in 2020 and 2030 accompanied by the contamination of the pumping wells with International Journal of Geosciences seawater.This requires a quarterly artificial recharge rate with a necessary volume which is around 2 m 3 /s at the pumping zone.
By giving this critical situation, it is highly crucial to start acting so as to stop this serious trend and avoid the irreversible degradation of land.This requires raising farmers' awareness to use the new economical mode of drip irrigation with respect to groundwater resources.
The seawater desalination alternative which is already put in place by the Agency of Hydraulic Basin of the Sakia Al Hamra and Oued Eddahab (ABHSHOD), and it can be developed to reduce the impact of pumping wells on the salt-wedge.
The elaborated models in this work gave a preliminary vision of the qualitative and quantitative aspects of the Fum Al Wad water-table.These models require other recent hydrogeological data to cover the entire area of the aquifer in order to further improve the quality and credibility of the results in this work.

Conclusions
The  Wad aquifer (Sections B in Figure 18).

Figure 3 .
Figure 3.The geographical situation of the used points in calibration models.
nary.The main geological formations constituting the aquifer of Fum Al Wad field are grouped into four families namely: lower sands, dolomitic limestone and sandstone lower, the upper sands and alluvium Plio-Quaternary of the Saquia Al Hamra Wad (Figure2).The groundwater of Fum Al Wad aquifer flows into the permeable formations of Plio-Quaternary which rest on an impermeable bottom that is formed by the gray marls of Miocene.It is alluvial phreatic water that is limited to the East by the aquifers of the right bank and the left bank of Laayoun city which are geologically continuous and unrelated hydraulic since their bottom (clays and impermeable marls of the Upper Cretaceous) appear in the middle of the Wadi (Saguia Al Hamra) giving birth to some sources of low flow at the level of Laayoun city.
for the range of expected densities and concentrations.For example, if meters and kilograms are used for the simulation, ∂ρ/∂C is set to a value of 0.7143, which approximately equals the change in fluid density divided by the change in solute concentration for freshwater and seawater.

Figure 4 .
Figure 4.The Geometry and the four lithological cross sections (West-East) of the Fum Al Wad aquifer represented in three dimensions.
The model grid consists of 106,476 active cells with an approximate cell dimension of 100 m × 100 m per each cell in the horizontal plane.The top surface of the model grid is interpolated to the terrain elevation which ranges from 0 to 60 m and the bottom elevation is plunged downward towards the Atlantic Ocean to reach −200 m/sea.The boundary conditions of this aquifer system are highlighted as the following:  The eastern limit of the Fum Al Wad aquifer is considered as a recharge limit, to take account of the waters coming from the Laayoun aquifer by lateral entry; it is represented by an imposed flow condition;  The western limit (limit with the Atlantic Ocean) is represented by a condition of imposed potential (H = 0 m);  The South and North limits are fictitious and are represented by a null flow.

Figure 5 .
Figure 5. Grid model and boundary conditions of the Fum Al Wad aquifer.

Figure 6 .
Figure 6.Fum Al Wad piezometry obtained by the calibration of the steady-state regime in 1986.

Figure 7 .
Figure 7. Hydraulic conductivity map of Fum Al Wad aquifer obtained by the calibration the steady-state in 1986.

Figure 8 .
Figure 8. Diagram of the observed and simulated values of the piezometric head obtained by the steady-state model simulation in 1986.

Figure 9 .
Figure 9. Fum Al Wad piezometry obtained by the calibration transient model in 2015.

Figure 10 .
Figure 10.Calibration results of the observed and calculated head values by the transient model (Control piezometer 277/120).

Figure 11 .
Figure 11.Calibration results of the observed and calculated head values by the transient model (Control piezometer 735/120).

Figure 12 .
Figure 12.Calibration results of the observed and calculated head values by the transient model (Control piezometer 815/120).

Figure 13 .
Figure 13.Calibration results of the observed and calculated head values by the transient model (Control piezometer 737/120).

Figure 14 .
Figure 14.Calibration results of the observed and calculated head values by the transient model (Control piezometer 707/120).

Figure 15 .
Figure 15.Variations of the flow budget of Fum Al Wad water-table calculated by transient regime during 1986-2015.

annually throughout 29
years, which started from 1986 to 2015.The results obtained from the hydrodynamic transient model (see Figure 9) are retained for the new simulation of the saltwater intrusion in the water-table.The new databases used for this simulation are mainly: the freshwater density (1 g/l) and the seawater density (1.025 g/l).Thus, the hydrodispersifs parameters obtained by the calibration transient model of the saltwater intrusion are mainly: longitudinal and vertical dispersivity and molecular diffusion.Their values are retained successively: 5 m, 5 m and 10 −9 m 2 /s.Hence, the model has been initialized with a salinity condition on the coastal ocean which is 35 g/l.The effective porosity International Journal of Geosciences
below illustrates the results of saltwater intrusion that is obtained by SEAWAT-2000 simulation program in 2015.The immigration of saltwater to the zone of pumping wells is usually caused by excessive water-extraction and the high permeability of sandy subsoil; the salt-wedge has advanced towards the pumping wells of the freshwater in Fum Al
borate some piezometric maps which testify the recent evolution of the aquifer state.The piezometry obtained by the hydrodynamic transient model in 2015 marked a piezometric decrease at the pumping wells zone with about 2.3 m, in comparison to the 1986 piezometric data.This is due to the increase of water-extraction rate for irrigation during the past 15 years.In addition, the calibration of the transient regime is carried out by using six piezometric controlling piezometers.The simulated piezometry of 2015 should be verified with the measured piezometry in the study area to elaborate a new piezometric reference of that year.This operation could constitute a fundamental element for validating the mathematical model of the Fum Al Wad aquifer.The flow budget calculated by the transient model during the 29 years (1986-2015) has showed an increase of the pumped volumes which reached 7.02 Mm 3 in 2015.International Journal of Geosciences This increase in the volume of the pumped water is accompanied by a decrease of the outflow discharge in the Atlantic Ocean.The numerical modeling of the saltwater intrusion has highlighted the lateral extension of the salt-wedge in the Fum Al Wad aquifer.The simulated results showed the location of the interface between saltwater and freshwater, as in section B of the transient model in 2015 (Figure 17), with 4.3 km approximately at the level of the pumping zone.The results of the simulated scenarios for the years 2020 and 2030 helped in predicting the future states of the piezometric head and of the saltwater intrusion in Fum Al Wad water-table.In this respect, the predicted results showed that the piezometric values decrease successively 2.4 and 2.5 m in 2020 and 2030.This decrease of the piezometric head, which is due to the over-exploitation of groundwater, will accompany an advancement of the salt-wedge, with 4.4 km in 2020 and 4.7 km in 2030, towards the freshwater in the pumping zone of Fum Al

Table 1 .
The five salinometers to controlling salinity of the Fum Al Wad aquifer (the coordinate system is: WGS84 UTM zone 28 northing).
Fum Al Wad aquifer represents an important reservoir for groundwater exploitation which is used for water supply, agriculture and industry.It supplies the Laayoun city, which has 217,732 citizens (2014 statistics), and it ensures the sustainability of the various economic activities in the region.This aquifer is located in a coastal area, west of the Laâyoune city; it extends on a surface of approximately 250 km2.Its limits are drawn based on facts and studies in the East and not identified in North and South.This resource belongs to a region with a Saharan climate where the average annual rainfall between 1993 and 2007 did not exceed a few tens of millimeters.Thus, the infiltration of rainwater is very low and insufficient for recharging the Fum Al Wad water-table.The annual water volume subtracted from the aquifer during the periodaccording to the data from the National Office of Drinking Water is constantly growing and reached approximately 2.91 Mm 3 in 2015.This exploitation was accompanied by a gradual decline of the water-table and a considerable advancement of the salt-wedge.The numerical modeling of the Fum Al Wad water-table has permitted to ela-