This study presents a planning model to assist decision-makers in implementing proposed sustainable development policies based on the government's development schemes. It focuses on the exploitation of groundwater from the Eocene aquifer in the West-West Minya region in the Western desert of Egypt as one of the new development areas. The visual MODFLOW model is applied to serve as a base model for the study area's local modeling and assess the impact of operating scenarios on the groundwater aquifer. An optimal scenario for groundwater sustainability is achieved by considering the water meter consumption in the newly reclaimed area with a value of 6 m3/d for each acre to irrigate 3953.68 acres. The results reveal a lower decline in groundwater levels than the other proposed scenarios, and dynamic groundwater levels reached a steady-state after an estimated 9700 days. The GIS Multi-Criteria Analysis Model is used to assess the impact of groundwater deficiency as a result of reclamation. Results of the groundwater model are merged as input layers in Multicriteria Decision Analysis, which include groundwater salinity, groundwater levels, aquifer transmissivity, and aquifer storativity. The prospective groundwater zones are categorized into Most suitable (1584 sq km, 9.9 %), Suitable (4592 sq km, 28.7%), Good (4784 sq Km, 29.9%), Moderate (1488 sq Km, 9.3%) and Unsuitable (3552 sq Km, 22.2%). It is recommended to optimize levels of groundwater withdrawal from the current or future drilling wells to achieve a balance between the use and the protection of water potential.
Research Article
Planning for groundwater management using Visual MODFLOW Model and Multi-Criteria Decision Analysis, West- West Minia, Egypt.
https://doi.org/10.21203/rs.3.rs-1319259/v1
This work is licensed under a CC BY 4.0 License
Version 1
posted
You are reading this latest preprint version
Groundwater management
Groundwater Modeling
Visual MODFLOW Model
Multi-Criteria Decision Analysis
On a global scale, groundwater is the primary source of freshwater. It accounts for about one-third of freshwater withdrawals (Siebert et al., 2010; Famiglietti, 2014). An estimated 2.5 billion people worldwide depend on groundwater for their use (Dhar Chakrabarti 2017). This condition is a fact in arid and semi-arid regions because of the scarcity of rainfall and the limitation of surface water masses which are sometimes even absent. Due to high population density, global climatic change, and seawater intrusion in coastal aquifers, the demand for groundwater resources is expected to increase about 110 billion cubic meters per year (BCM/Y) (El Alfy, M., and Abdalla, F., 2021). In Egypt, the surface water from the Nile River is a primary source of the daily drinking water production about 22.1 million m3 with 82.1% (HCWW 2017); and the groundwater is an integral part of its national policy where it is used for drinking, domestic and industrial purposes. With the continuous deficiency of fresh water and the importance of groundwater, the Egyptian government started a new development project to reclaim 1.5 million acres to increase agricultural areas. Groundwater is the only source of irrigation water in the newly reclaimed areas; it is also a supplementary source of irrigation in the old, cultivated areas. Two main aquifer systems are relevant in Egypt. The first is the granular aquifer system represented by Nile Valley and Delta, Nubian sandstone, coastal and Moghra systems. The second is the fissured and the karstified aquifer system represented by limestone and hard igneous and metamorphic rocks (Hefny and Shata 2004, El Alfy and Abdalla, 2021). These aquifer systems comprise five hydrological zones: the Nile Valley and Delta, Western Desert, Eastern Desert, Coastal region, and the Sinai Peninsula; they are referred by many authors (RIGW, 1989, 1992; MPWWR 1988; RIGW-IWACO 1988,1993,1999; Abdel Moneim 2009). The purpose of the present work is to present a planning model for decision-makers for groundwater management and to predict the environmental impact of long-term extraction schemes on groundwater. And to achieve the optimum groundwater exploitation plans for sustainable development projects (domestic water supply, agriculture, and industry) with further rationing of the existing wells in line with groundwater potential and detect the suitability of groundwater after reclamation proceeded. The Eocene fractured limestone aquifer, occupying the West-West Minya area in the Western desert of Egypt, has been investigated as a case study of one of the new agricultural development areas. The West-West Minya agricultural development area extends further west from the West Minya agricultural zone.
The study area lies at the eastern border of the Western Desert in the central part of Egypt between latitudes 27ᵒ00ʹ-29ᵒ00ʹ and longitudes 29ᵒ30ʹ-31ᵒ00ʹ with an area of about 16000 Square Kilometers (Fig. 1). Its climate is typically characterized by hot, dry summers and mild winters. According to the NCEI (2021) and EMA (2016), the average climatic records for the period (2000-2019) are of temperature ranging from 31.9ᵒC to 14.1ᵒC, and average moisture percentage ranges from 30–40%. Annual precipitation intensity is around 0 to 5 mm/year. For the implementation of this work, two effective techniques have been employed. In the first technique, a three-dimensional groundwater flow model based on GIS-Database integration is applied to simulate the hydrological conditions of the Eocene aquifer in the study area and interpolate the aquifer response and groundwater level variation under different pumping stresses. Model development scenarios involving government operating plans have been suggested. In the second technique, the obtained output results of the optimal proposed scenario have been merged under the GIS framework and used as input criteria for GIS-Multi-Criteria Decision Analysis (MCDA) model to delineate the groundwater suitability zones after groundwater withdrawal during reclamation processes. El Alfy, M., (2014), Abdel Moneim, A. et al. (2016), and Abdelhalim, A. et al. (2019) applied the numerical groundwater model to assess the responses of the aquifer under different pumping scenarios. It is proved the efficiency of this technique in managing groundwater resources. Sayed et al. (2019) investigated the feasibility of using solar energy to exploit groundwater supplies from the Moghra Oasis in the Western Desert of Egypt using a multicriteria analysis model. El-Hadidy and Morsy (2022) integrated the Visual MODFLOW model and GIS model and applied the weighted overlay approach to producing the groundwater potentiality map in the Nile Valley region.
Studying the surface and subsurface geological characteristics of the area concerned is vital for detecting the hydrogeological configuration in terms of aquifer properties, water potentiality, and water flow. The geomorphology of the El Minya province resembles the prevailing conditions in the Nile Valley. The limestone plateau borders the valley from both sides along its length; causes a steep slope of the earth's surface on both sides of the Nile Valley to the east and west. Three geomorphological units are studied (Fig. 2) (Said 1981):
a. Young alluvial plains (Holocene silty clay) form the cultivated lands bordering the course of the Nile River on both east and west sides and dissected irrigation canals and drains.
b. Old alluvial plains (sand and gravel from the Pliocene and the Quaternary) are exposed as terraces at different heights above the young alluvial plains.
c. Limestone Calcareous Structural Plateau (Middle Eocene limestone) is structurally formed and bounds the Nile Valley east and west. The surface lithologic units exposed in the study area range from Middle Eocene limestone to Oligocene-Pleistocene gravel and sand and Quaternary deposits (Fig. 3) (Abd El-Aziz, 1994; Abou Heleika and Niesner, 2009). The subsurface stratigraphic sequence in the study area ranges from Pre-Cambrian to the Quaternary ages; it is built-up up from the base to top as follows:
a. Pre-Cambrian Basement Complex's surface elevations in the study area range from -1500m to -4000 (CEDARE, 2002).
b. The Nubia Sandstone Sequence rests unconformably on the rugged surface of the basement complex. It comprises all clastic sediments ranging from the Upper Cretaceous to the Pre-Cambrian Basement consisting of alternating sandstone beds, shale, and clay with thickness ranges from 500n to 2500m in the study area (LaMoreaux et al. 2008).
c. Upper Cretaceous shale mainly consists of shale, clay, and mudstone with thin limestone interbeds with phosphates bands forming the confining formation to the Nubia Sandstone northward El-Kharga and El-Dakhla Oases. The thickness of these variegated shales varies between 126 m in El-Dakhla and 265 m in El-Kharga (MPWWR, 1998).
d. Tertiary carbonate rocks of Eocene age consist of limestone dolomitic limestone and cover the plateau area surrounding the depressions of south-Kharga, El-Kharga, El-Dakhla El-Farafra, El-Bahariya, Siwa, and Qattara. Their thickness ranges from 200m to 1200m (RIGW;1992, 2015).
e. The Oligocene sequence consists of alternating classics, siltstone, and claystone. It overlays the Eocene limestone and occupies the western part of the Western Plateau. (Abdel Moneim et al.,2016)
f. Quaternary deposits are widely distributed in the eastern side, including the Nile silt, dunes, and fanglomerates (Korany et al., 2008).
The structural framework is composed of a system of wrench faults parallel either to the NW-SE trend of the Gulf of Suez or to the NE-SW trend of the Gulf of Aqaba. These faults play a vital role in the groundwater aquifers in the study area (Youssef,1968; Shabana, 2010; El Kashouty et al., 2010; Embabi, 2018).
Four aquifers are present in the study area: Quaternary alluvial deposits, Oligocene sandstone, Eocene limestone, and Nubian sandstone aquifers. The groundwater of the Quaternary and Oligocene sandstone is recorded under unconfined conditions, and Pleistocene sediments overlie the fractured Eocene limestone. In contrast, the Nubian sandstone is registered as a confined aquifer (Yousif et al.,2018). Regionally, groundwater exhibits a northward trending gradient in this area. The recharge of the Upper Eocene water-bearing units may occur by seepage from the Nile and by upward percolation of groundwater in hydraulic charge from the deep NSA through the upper confining Cenomanian shale. No recharge of the upper Pleistocene aquifer is expected to occur from the average 4 mm of the annual rainfall occurring as very limited precipitations and believed to be lost in evaporation (Tantawi, 2005). The Middle Eocene aquifer has an isotopic signature of the modern Nile. (Yousif et al.,2018). The investigated groundwater flow pattern in the Eocene aquifer shows the flow direction from the southwest. The total hydraulic head value is 50m and reaches 30m northeast (Fig. 4). The present estimation of groundwater quality of the Eocene aquifer in the new development study area is based on the results of the chemical analysis of groundwater samples collected from the boreholes. The Total Dissolved Solids (TDS) content varies between 2000 to 3000 mg/l (Fig. 5). The spatial extension of the Eocene aquifer in the study area and its approximated thickness is calculated using spatial interpolation Inverse Distance Weighted (IDW) in GIS Tools (Fig. 6).
The procedures of this work are based on applying and integrating two techniques, the Visual MODFLOW groundwater model and Multi-Criteria Decision Analysis (MCDA) model. The operational steps and modeling processes are illustrated in a flow chart (Fig. 7).
3.1. Numerical modeling for groundwater system
3.1.1 Designing the conceptual model
Groundwater flow models are used to calculate the rate and direction of groundwater movement through aquifers and confining units in the subsurface; these calculations are referred to as simulations (Mandle, 2002). The model is also used to simulate possible future changes in hydraulic head or groundwater flows due to future changes in aquifer system constraints. MODFLOW 2005 is used as a numerical model. It is a standard USGS software developed by McDonald & Harbaugh (1988) and Harbaugh (2005). The integral finite difference method is employed to resolve the transient flux equation. There are virtually no hydrological studies specifically focused on the study area. Currently, the drilling of test wells is being carried out under the Ministry of Water Resources and Irrigation (MRWI) supervision. Available and updated field data from wells in the Eocene aquifer are collected and organized by integrating the attributes and spatial databases with GIS tools. GIS is used to manage the model's spatially distributed input parameters and outputs. The database constructed in the GIS technique is imported into the conceptual model, including the SRTM (Dem) 1-ARC resolution with 30 m, potentiometric water levels, hydraulic properties, and recharge and discharge components.
3.1.2. Grid system and boundary conditions
The model domain’s dimensions in the study area are discretized with a uniform, regular 500x500 m grid, the number of rows and columns is 200 and 150, respectively. To obtain results that reflect the hydrological conditions of the aquifer, it is essential to describe its hydraulic boundaries with known hydraulic heads. Boundary conditions are based on the hydraulic head map and hydrological conditions. Due to the absence of major branches or banks in or near the area at the northern, southern, and western boundaries, the initial hydraulic conditions were used to identify the sources of recharge of the study area. It is assumed that the aquifer is extended and represented by constant hydraulic pressure limits. The flow through these boundaries is determined as constant head boundaries (Dirichlet Conditions), where the derivatives of the groundwater head are constant and do not change with time. The eastern boundary of the model is a specified head boundary to represent the Nile River. For other areas, the no-flow boundary was specified (Neumann conditions) in which the derivatives of the head (flux) are set to zero. All boundaries were set distance from the area of main interest and simulated pumping wells to avoid the impact of uncertainty in subsurface boundary inflows and outflows and drawdown effect during transient simulation (Fig. 8).
3.1.3. Initial conditions
When studying the dynamic behavior of the aquifer to determine the amount of change in groundwater levels, it is necessary to decide on the initial conditions of the groundwater levels. Due to lack of available data for the study area, current static groundwater levels of the investigated wells have been used as the initial distribution head for the model (Fig. 9)
3.1.4. Aquifer hydraulic parameters
The present estimation of aquifer characteristics, including transmissivity, hydraulic conductivity, and storativity, is based on step-drawdown, constant discharge, and recovery tests for the investigated wells in the modeled area. Six wells were selected to represent field measurement data; their locations are shown in Fig. (4). The data of step drawdown tests are plotted on a linear scale (Fig. 10) to estimate the general well equation using Jacob (1947) and Rorabaugh (1953) techniques.
$${S}_{w }=BQ+C{Q}^{2}$$………………… 1
Where:
\({S}_{w:}\)is the total drawdown in the pumping well (m)
\(BQ\): is the formation loss (m) caused by laminar flow through the aquifer.
\(C{Q}^{2}\) is the well loss (m) caused by turbulent flow inside the well and its filter and damage zones. From this equation, the specific capacity curve is drawn. The initial Specific Capacity value is used to estimate the transmissivity (T) of the aquifer by using Driscoll (1986) equation:
$$T=2000\times \frac{Q}{S}$$…………………2
Where: \(T\) = Transmissivity (gpd/ft) and \(\frac{Q}{S}\) = Specific Capacity (gpm/ft). By using Formula/Conversion Table for Water Treatment and Water Distribution (FDEP, 2020), equivalent units for transmissivity (m2/day) and specific capacity (m2/hr.) values are obtained.
The records of pumping and recovery tests were represented graphically to estimate the transmissivity (T), storativity (S), and hydraulic conductivity(k) values by using Jacob’s straight-line method (Cooper and Jacob, 1946) and Theis’s recovery methods (Theis, 1935). (Fig. 11)
Table (1): Hydraulic parameters of the Eocene aquifer, West-West Minia, Egypt.
|
Well No. |
Step drawdown test |
Jacob straight-line method |
Theis recovery method |
||||
|---|---|---|---|---|---|---|---|
|
Specific capacity (m2/hr.) |
T (m2/d) |
T (m2/d) |
S |
K (m/d) |
T (m2/d) |
K (m/d) |
|
|
44 |
71.43 |
5952 |
7000 |
5.8×104− |
3.88 |
||
|
247 |
17.14 |
1428.6 |
2900 |
6.1×104− |
2.41 |
||
|
182 |
11.26 |
938.3 |
1100 |
0.916 |
|||
|
43 |
42.2 |
3516.6 |
5000 |
4.166 |
|||
|
42 |
13.2 |
1100 |
1500 |
1.25 |
|||
|
53 |
26.5 |
2208.33 |
3200 |
2.13 |
|||
3.1.5. Steady-State calibration
The objective of model calibration is to assure that it can produce field-measured heads, which are calibration values. The calibration process has been carried out for the boundary conditions and the hydraulic conductivity spatial distribution modification to match between the observed and simulated groundwater heads (Figs. 12&13).
3.1.6. Transient simulation
Transient simulation is needed to solve time-dependent groundwater problems. The calibrated steady-state model is run under transient calibration for one year to adjust the storativity values (Fig. 14). The model is subjected to sensitivity analysis to quantify the uncertainty in the calibrated parameters. It revealed the high sensitivity of the model toward the boundary conditions, hydraulic conductivity, and storativity, respectively. The new agricultural development project in the study area covers 4000 acres. Proposed pumping wells were distributed with a spacing of 1000 meters to decrease the pressure on the aquifer and reduce the interference between the wells. Virtual observation wells were distributed among the extraction wells to monitor the groundwater levels under the impact of each development plan.
3.2. Model results
Based on the future governmental development plan in the West-West Minia area, model scenarios have been proposed to support the decision-makers with the results and recommendations. The first scenario has been assumed to interpolate the sustainability of the groundwater aquifer in the newly reclaimed area under the current operation scheme by applying the proposed pumping wells with pumping rates of 150 m3/d for irrigation of 3953.68 acres. In the second proposed scenario, the sustainability of the groundwater aquifer in the newly reclaimed area is interpolated by considering the water meter consumption in the newly reclaimed area with a value of 6 m3/d for each acre (MALR 2007). The transient simulations have been carried out for 100 years divided into 10-time steps. Boundary conditions and wells pumping rates were constant during transient simulation.
3.2.1. Impact of the first proposed scenario
By applying the first proposed scenario, a significant depletion in hydraulic heads results from the end of each time step. The cone of depression is formed gradually at the center of the modeled area, reaching a maximum drawdown of about 100 m. The dynamic groundwater level got steady-state conditions after an estimated period of 5000 days (Fig. 15). Negative values of groundwater levels are recorded; it reached -80 m.b.s.l.
3.2.2. Impact of the second proposed scenario
This scenario results show a smaller decrease in groundwater level than the first proposed scenario. The lowest groundwater level is recorded at -20 m.b.m.s.l. The cone of depression with a maximum drawdown of about 60 m is formed. The dynamic groundwater level reached steady-state conditions after an estimated period of 9700 days (Fig. 16)
3.3. Multi-criteria Decision modeling
Based on the obtained results of the optimal scenario for groundwater sustainability of the applied Visual MODFLOW model, Multi-Criteria Decision Analysis (MCDA) model technique (Malczewski, 1999&2000) is used to assess the impact of groundwater consumption by detecting the groundwater suitability zones after long run reclamation in the concerned development area. The (MCDA) model is not a substitute for cost-benefit analysis; it is a complement analysis it gives ranges of decisions. In this work, the MCDA model integrates geospatial data of the most effective parameters in the groundwater potentials both for aquifer hydraulic characteristics and groundwater quality by merging all multicriteria layers into one prospective layer for groundwater exploitation. The contributing criteria are groundwater levels (GWL) which are obtained from the outputs of the optimal simulated scenario (second scenario) of groundwater flow model, groundwater salinity (TDS), aquifer transmissivity (T), and aquifer storativity (S). Arc-GIS is used to generate a thematic map for each criterion. All-controlling criteria thematic maps are reclassified in ARC-GIS to use in the model as a layer (Fig. 17). Each criterion is ranked according to its priority in the groundwater assessment and integrated into a GIS database model using the weighted overlay method. Weighted-Sum Overlay analysis provides the ability to obtain weight for each ranked criterion. It combines all inputs by overlaying all rasters, multiplying each by its given weight, and summing them together to calculate the overall suitably of each cell and create one output raster layer that defines the appropriate groundwater zones. The ranks, rates, and attributes of input criteria are indicated in the table (2). The equation used in this method is defined as the following:
………………… 3Where,
-
WI: is the Relative weight of the I criterion
-
n: is the number of criteria
-
Wi: is the rank of the criterion
Table (2): Categorizing the contributing criteria and their rank, weight, and rate affecting the groundwater assessment.
|
Thematic layer |
Rank |
Weight |
Relative weight |
Attribute |
Rate |
|---|---|---|---|---|---|
|
Groundwater salinity |
1 |
4 |
0.4 |
< 500 2000-3200 |
1 2 |
|
Groundwater levels |
2 |
3 |
0.3 |
< 0 7-0 7-15 |
3 2 1 |
|
Transmissivity |
3 |
2 |
0.2 |
<1000 1,000-3000 3,000-6000 6,000-10,000 10,000-23,300 |
5 4 3 2 1 |
|
Storativity |
4 |
1 |
0.1 |
< 0.00552 0.00552 - 0.02443 0.02443 - 0.05911 0.05911 - 0.11980 0.11980 - 0.200991 |
5 4 3 2 1 |
3.4. Assessing groundwater prospects in the new development area
The raster layer generated from Multi-Criteria Decision Analysis (MCDA) is shown in Fig. (18). It defines the prospective groundwater zones which categorized into Most suitable (1584 sq km, 9.9%), Suitable (4592 sq km, 28.7%), Good (4784 sq Km, 29.9%), Moderate (1488 sq Km, 9.3%), and Unsuitable (3552 sq Km, 22.2%). It is observed that nearly 80% of the study area has moderate to good potential for groundwater prospecting. There is a matching between the Unsuitable zone of groundwater and the higher salinity zone of groundwater, which indicates the first rank and the highest weight of the salinity factor in groundwater development. In the rest of the area, where the salinity is almost constant, the Most suitable and Suitable groundwater zones match with high groundwater level zones, indicating the second rank and the high weight of the groundwater level factor in groundwater development.
-
The Eocene aquifer in the West-West Minya area in the Western desert of Egypt is one of the new development areas with a total area of 80,000 acres. It has been investigated to assess the aquifer for sustainable development.
-
The Visual MODFLOW model is a powerful tool to estimate the impact of operating scenarios on the groundwater aquifer and assist decision-makers in implementing proposed policies for sustainable development. The model's output results indicated the high efficiency of the Eocene aquifer as a source of groundwater for irrigation of 4000 acres, and the daily withdrawal rate from wells should not exceed 23722.2 m3/day according.
-
Applying the Multi-Criteria Decision Analysis (MCDA) model revealed the groundwater suitability zones after long-run reclamation. About 80% of the study area has moderate to good potential for groundwater prospecting.
-
In sustainable development projects, it is necessary to assess the continuity of groundwater resources under current and future operating conditions considering the expected population increase and increased water needs of industry and agriculture.
-
Optimizing levels of groundwater withdrawal from the current or future drilling wells to achieve a balance between the use and the protection of water potential and address their impacts on aquifer storage and depletion trajectory management is the main target of groundwater assessment policies.
Financial interests: The author declares that no funds, grants, or other support were received during the preparation of this manuscript.
Declaration of Competing Interest: The author declares that she has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Author Contributions: The author has read and agreed to the published version of the manuscript.
- Abdel Aziz RS (1994) : Geological and sedimentological studies in west El Minia – Beni Mazar area, Egypt. M.Sc. Thesis
- Abdel Moneim A, A (2009) Overview of water resources and requirements in Egypt; the factors controlling its management and development. Journal of Environmental Studies 2:85–97
- Abdel Moneim A, Álvarez A, Abu El Ella JPF, Masoud A, M (2016) Groundwater Management at West El-Minia Desert Area, Egypt Using Numerical Modeling. Journal of Geoscience and Environment Protection 4:66–76. http://dx.doi.org/10.4236/gep.2016.47008
- Abdelhalim A, Sefelnasr A, Ismail E (2019) Numerical modeling technique for groundwater management in Samalut city, Minia Governorate, Egypt. Arab J Geosci 12:124. https://doi.org/10.1007/s12517-019-4230-6
- Abu Heleika M, Niesner E (2009) Configuration of the limestone aquifers in the central part of Egypt using electrical measurements. Hydrogeol J 17(2):433–446. https://doi.org/10.1007/s10040-008-0360-8
- Centre for Environment & Development for the Arab Region and Europe (CEDARE) and International Fund for Agricultural Development (IFAD) (2002) : Regional Strategy for Utilisation of The Nubian Sandstone Aquifer System. Online: http://water.cedare.int/cedare.int/files15/File2919.pdf
- CONOCO (1987) Geologic map of Egypt. Egyptian General Authority for Petroleum (UNESCO joint map project), 20 sheets, scale 1:500 000. Cairo, Egypt
- Cooper HH, Jacob CE (1946) : A generalized graphical method for evaluating formation constants and summarizing well history.Am. Geophys. Union.Trans.;27526–27534
- Dhar Chakrabarti (2017) : People, Planet, and progress beyond 2015. The Energy and Resources Institute (TERI), New Delhi, p 462. "MENA Report, Albawaba (London) Ltd., Aug. 2020
- Driscoll FG (1986) : Groundwater and wells. 2nd edition, Johnson Division, St. Paul, Minnesota, 1089 pp
- El Alfy M (2014) Numerical groundwater modeling as an effective tool for the management of water resources in arid areas. Hydrol Sci J 59(6):1259–1274. https://doi.org/10.1080/02626667.2013.836278
- El Alfy M, Abdalla F (2021) : Groundwater Resources in Egypt, pp. 699-711. In: Hamimi, Z., Arai, S., Fowler, A., El-Bialy, M.Z. (Eds.); The Geology of the Egyptian Nubian Shield, Part of the Regional Geology Reviews book series (RGR), Springer Nature Switzerland, 711p https://doi.org/10.1007/978-3-030-49771-2
- El Kashouty M, El Sayed E, Kamel A (2010) The hydrochemical characteristics and evolution of groundwater and surface water in the western part of the River Nile, El Minia district, Upper Egypt. Arab J Geosci. http://dx.doi.org/10.1007/s12517-010-0219-x
- El-Hadidy S, Morsy S (2022) M., : Expected Spatio-temporal variation of groundwater deficit by integrating groundwater modeling, remote sensing, and GIS techniques. The Egyptian Journal of Remote Sensing and Space Sciences 25 (2022) 97–111 https://doi.org/10.1016/j.ejrs.2022.01.001
- Egyptian M, Authority EMA (2016) Annual Meteorological reports. Cairo, Egypt
- Embabi NS (2018) The Geologic Setting. Landscapes and Landforms of Egypt. World Geomorphological Landscapes. Springer, Cham. https://doi.org/10.1007/978-3-319-65661-8_2
- Famiglietti JS (2014) The global groundwater crisis. Nat Clim Change 4:945–948
- Florida Department of Environmental Protection (FDEP) (2020) : Water Treatment and Water Distribution Formula/Conversion Table, https://floridadep.gov/sites/default/files/DW-DS%20Formula%20Sheet%20June%202020.pdf
- Harbaugh AW (2005) : MODFLOW-2005, the U.S. Geological Survey Modular Ground-Water Model—The Ground-Water Flow Process. U.S. Geological Survey Techniques and Methods 6-A16.USGS TM 6-A16. https://pubs.usgs.gov/tm/2005/tm6A16/
- HCWW (2017) : Holding Company for Water and Wastewater, Cairo. Internal report
- Hefny K, Shata A (2004) Underground water in Egypt. Ministry of water supplies and irrigation. Cairo, Egypt, p 295. (Unpublished report)
- Jacob CE (1947) : Drawdown test to determine the effective radius of artesian well: Trans. Amer. Soc. of Civil Engrs., Vol.112, paper 2321, pp.1047-1064
- Korany E, Sakr S, Darwish M, Morsy S, World C (2008) Univ., p. 703-711
- LaMoreaux PE, James WL, Mostafa MS, Bashir AM, Fakhry AA (2008) : Environmental Hydrogeology, Second Edition, CRC Press, 375 pages
- McDonald MG, Harbaugh AW (1988) : A modular three-dimensional finite-difference groundwater flow model. U.S. Geological Survey Techniques of Water-Resources Investigations Report 06-A1
- Malczewski J (1999) GIS and Multicriteria Decision Analysis. Wiley, New York
- Malczewski J (2000) : On the use of weighted linear combination method in GIS
- common and best practice approaches.Trans. GIS4 (1),5–22
- Mandle RJ (2002) Groundwater Modeling Guidance. Groundwater Modeling Program, Michigan Department of Environmental Quality
- MALR (2007) Ministry of Agriculture and Land Reclamation. The possibility of achieving self-sufficiency in Egyptian crops. Internal report
- MPWWR (1988) : Minister of Public Works and Water Resources. Rehabilitation and improvement of water delivery systems in old lands. Proj. No. EGY/85/012. Final report
- MPWWR (1998) : Minister of Public Works and Water Resources. Hydrogeology of deep aquifers in the Western Desert and Sinai. Report No. 10
- NCEI (2021) : National Centers for Environmental Information https://www.ncei.noaa.gov/
- RIGW (1989) Research institute for groundwater. Hydrogeological map of Egypt, Cairo, scale 1:100,000, 1st edn. El-Qanatir
- RIGW (1992) Research institute for groundwater. Hydrogeological map of the Nile Delta (1st edn) project development and management of groundwater resources in the Nile Valley and Delta. RIGW, Cairo
- RIGW/IWACO (1988) 1993): Research institute for groundwater, Hydrogeological map of Egypt, scale 1:2,000,000. Cairo
- RIGW/IWACO (1999) Research institute for groundwater, Contributions to environmental management of Egypt’s groundwater resources. Final report. EMGR Project, Cairo
- RIGW (2015) : Research institute for groundwater, technical report, internal report
- Rorabaugh MJ (1953) : Graphical and theoretical analysis of step-drawdown test of artesian well: Proc.Amer.Soc. Civil. Engrs., Vol.79, separate no. 362,23p
- Said R (1981) The Geological Evolution of the River Nile. Springer-Verlag
- Sayed E, Riad P, Elbeih S, Hagras M, Hassan A (2019) A., : Multi-criteria analysis for groundwater management using solar energy in Moghra Oasis, Egypt. 22 (2019) 227–235 https://doi.org/10.1016/j.ejrs.2019.04.001
- Siebert S, Burke J, Faures JM, Frenken K, Hoogeveen J, Doll P, Portmann FT (2010) Groundwater use for irrigation: A global inventory. Hydrol Earth Syst Sci 14:1863–1880
- Shabana AR (2010) Hydrogeological studies on the area West Deir Mouas-Mallawi, El Minia governorate-Egypt. Egypt J Geol 54:61–78
- Tantawi M, Elsayed E, Setto I, Abuelsebaa E (2005) Hydrochemical evaluation of the groundwater resources in the area east of Minia District, Egypt. Geophysical Research Abstracts 7:00360
- Thesis CV (1935) The relation between the lowering of the piezometric surface and the rate and duration of discharge of well-using groundwater storage. Trans Amer Geophys Union 16:519–524
- Youssef MI (1968) Structural pattern of Egypt and its interpretation. Bull Amer Petrol Geol 52(4):601–614
- Yousif M, Sabetb HS, Ghoubachia SY, Aziz A (2018) Utilizing the geological data and remote sensing applications for investigation of groundwater occurrences. Western Desert of Egypt, West El Minia
- NRIAG Journal of Astronomy and Geophysics Volume 7 (2018) - Issue2. https://doi.org/10.1016/j.nrjag.2018.07.002