Modelling Flow in Aquifer
Modelling Flow in Aquifer
Modelling Flow in Aquifer
2, 2007
Associate Professor, Department of Civil Engineering, Jordan University of Science and Technology,
Irbid, Jordan, E-mail: wail@just.edu.jo
2)
PhD Student, Department of Civil Engineering, The Hashemite University, Zarqa, Jordan,
E-mail: rhatamleh@gmail.com
ABSTRACT
Numerical modeling has emerged as an effective tool for managing groundwater resources and predicting future
responses, especially when dealing with complex aquifers systems and heterogeneous formations. Among these
models, MODFLOW and MT3D are the most commonly used simulators for groundwater flow and solute transport
in subsurface systems, respectively. These models have been used herein as a management tool for the Azraq basin,
one of the most important groundwater resources for domestic and agricultural sectors in Jordan. Groundwater
extraction from this basin already exceeded the safe yield of the aquifer, and a sharp drop in the water table, a dry
out of the springs at the center of the basin and the problem of increased salinity in many parts of the aquifers have
been reported. Currently, more than 600 wells including governmental, private and unauthorized wells are operating
within the basin boundary. In its attempts to restore and sustain the aquifer, the Ministry of Water and Irrigation
(MWI), the official entity in charge of water resources in Jordan, is considering several scenarios of controlled
pumping. In this study, five suggested scenarios of pumping with different abstraction rates for years 2005 through
2020 have been explored by using the three dimensional finite difference flow model [MODFLOW (PM5)] to
simulate the flow system, and the solute transport model (MT3D) to predict the transport of total dissolved solids
given in terms of Electric Conductivity (EC). These scenarios include: first, maintaining the current pumping rate of
57 MCM for the study period; second, reducing the current pumping rates by half; third, increasing the pumping
rate by half; fourth, reducing the pumping rates in public wells by half and maintaining the current rates for other
wells; and finally reducing the pumping rates by half for the farm wells (private) and maintaining the rates at the
other wells. Results indicate that the first and fourth scenarios have similar effect on the drawdown. Also, the
second and fifth scenarios have similar effects and provide the lowest drawdown values. The third scenario gives
the worst drawdown. The transport of Total Dissolved Solids (TDS) given in terms of Electric Conductivity (EC),
has also been explored. Different parameters including EC, recharge, model boundary and advection parameters
were adjusted to run the model. Simulation results indicated that the effect of the different scenarios on the values of
EC is less profound than the effects on the drawdown values. The third scenario caused a slight increase of EC
values over the values simulated by other scenarios.
Keywords: Numerical Models, Modflow, MT3D, Groundwater Flow, Transport, Management
Scenarios, Azraq Basin.
-153-
INTRODUCTION
With nearly half the population of Jordan drawing its
drinking water from groundwater resources, issues of
groundwater flow and contaminant transport in
subsurface systems have elicited considerable interest
from both public and private sectors. Among the twelve
groundwater basins in Jordan, Azraq basin is one of the
most important since it provides the three main cities of
Jordan, namely: Amman, Irbid and Zarqa with significant
portions of its domestic water supplies. It consists of a
three aquifer system, the upper (unconfined), middle and
lower aquifers, separated by two confining formations.
The water is mainly pumped from the upper aquifer and
only few wells penetrated the middle aquifer and the deep
aquifer, respectively.
In the beginning, water was pumped from the basin
to supply the city of Irbid only, and the increased water
demand due mainly to the explosive population growth in
the cities of Amman and Zarqa, led to drill new well
fields and draw groundwater at high rates. Furthermore,
private wells in the basin increased dramatically in the
past two decades; some of these wells are not even
authorized. The total pumping from the basin was
estimated to reach 57.7 MCM/yr in the year 2001 (AlHadidi and Subah, 2001), which is far beyond the safe
yield of the basin (30 MCM/yr.). Absence of full control
on private wells and its pumping rates in addition to
pumping from the water authority well fields resulted in
overexploitation of water storage and deterioration of
groundwater quality in some parts of the aquifer.
The Azraq basin received great attention, and
numerous studies were conducted for estimating the
potential safe yield and for assessing current and potential
future scenarios (Agrar and Hydrotechnick, 1977;
Humphry and Sons, 1978 and 1982; Rimawi, 1985; AlMomani, 1993; Al-Kharabsheh, 1995; Ayed, 1996;
UNDP- Azraq Oasis Conservation Project, 1996;
Jordanian Consulting Engineering, 1997; Al-Khatib,
1999; Hatamleh, 2001; Al-Hadidi and Subah, 2001;
Rihani, 2003; Moqbel, 2004). All computer modeling
studies of Azraq basin (except for Rihani, 2003 and
- 154 -
- 155 -
Figure 1: Map of Jordan showing the location of Azraq basin (Al-Hadidi and Subah, 2002).
The Upper Aquifer System is an unconfined aquifer
except at the middle of the basin (Qaa Azraq). The Qaa
Azraq deposits form a lens up to 15 meters thick
overlying the aquifer and give rise to leaky artesian
- 156 -
- 157 -
(1)
- 158 -
h
h h
Kx + Ky + Kz w = 0.0
x
x y
y z
z
(2)
t xi xi xi
(3)
k =1
Figure 2: Steady state flow calibration matching map between calculated and
measured water levels for the upper aquifer.
Steady State Calibration
Steady state calibration for the flow model was
achieved by comparing the hydraulic heads obtained from
available groundwater level contour maps of the first and
second model layer and the calculated hydraulic heads of
the MODFLOW simulation in order to simulate the flow
- 159 -
(3), respectively.
As shown in Figure (4), the horizontal hydraulic
conductivity values for the upper aquifer ranged from
0.10-90.00 m/ day for the first model layer (B4/B5),
Figure 3: Steady state flow calibration matching map between calculated and
measured water levels for the middle aquifer.
Transient Calibration
For transient simulation, the time was divided into
fifteen-stress period. The first one stared from the
beginning of the 1970s to the end of 1984 with fifteen
time steps. This period represented the steady state
period where there was no significant drawdown in the
water level. The second starts from the beginning of
1985 till the end of that year with one time step. And the
- 160 -
- 161 -
- 162 -
- 163 -
Table 1: Prediction results for the flow and transport model with the same current
discharge 57 MCM/year (The first scenario).
Year
Observation well
Max.
F1022
F1043
F1060
F1280
Draw down (m)
25.7
5.3
13.6
11.1
11.7
2005
EC (s/cm)
1882
1069
853
1776
820
Draw down (m)
28.0
6.2
14.7
12.9
12.8
2010
EC (s/cm)
1886
1099
867
1781
924
Draw down (m)
32.4
7.8
16.8
15.9
15
2020
EC (s/cm)
1905
1122
873
1804
929
Year
2005
2010
28.0
20.5
35.6
28.0
20.5
2020
32.4
20.8
43.7
32.4
20.8
The EC values for the year 2020 are shown for the first
and third scenarios in Figures (13) and (14), respectively.
There is no significant change in the values of EC through
the different scenarios for the same year.
SENSITIVITY ANALYSIS
CONCLUSION
- 164 -
Figure 7: Comparison of observed and simulated drawdown for well F1043 (Governmental well area).
Figure 8: Comparison of observed and simulated drawdown for well F1060 (Farm well area).
- 165 -
- 166 -
Figure 10: Measured iso-EC (S/cm) for the upper aquifer system
(JCE, 1997).
- 167 -
- 168 -
- 169 -
- 170 -
REFERENCES
Agrar, Hydrotechnik and GTZ. 1977. National Water Master
Plan of Jordan. 8 Volume, Essen-Hanover.
Al-Hadidi K. and Subah, A. 2001. Jordan Badia Research
and Development Program: Integrated Studies of Azraq
Basin for Optimum Utilization of the Natural
- 171 -
- 172 -