Hydrochemical and geoelectrical investigation of the coastal shallow aquifers in El-Omayed area, Egypt

Monitoring and assessment of the coastal aquifers are becoming a worldwide concern for the need of additional and sustainable water resources to satisfy demographic growth and economic development. A hydrochemical and geoelectrical investigation was conducted in the El-Omayed area in the northwestern coast of Egypt. The aim of the study was to delineate different water-bearing formations, provide a general evaluation of groundwater quality, and identify the recharge sources in aquifers. Thirty-seven water samples were collected and chemically analyzed from the sand dune accumulations and oolitic limestone aquifers. Fifteen profiles of vertical electrical soundings (VESs) were obtained in the oolitic limestone aquifer to examine the variations of subsurface geology and associated groundwater chemistry. The groundwater reserves in the El-Omayed area are mainly contained in sand dune accumulations and oolitic limestone aquifers. The aquifer of sand dune accumulations contains freshwater of low salinity (average total dissolved solids (TDS)?=?974 mg/l). Groundwater of oolitic limestone aquifer is slightly brackish (average TDS?=?1,486 mg/l). Groundwater of these aquifers can be used for irrigation under special management for salinity control, and regular leaching as indicated by electrical conductivity and sodium adsorption ratio. Results of VES interpretation classified the subsurface sequence of oolitic limestone aquifer into four geoelectric zones, with increasing depth, calcareous loam, gypsum, oolitic limestone, and sandy limestone. Oolitic limestone constitutes the main aquifer and has a thickness of 12–32 m.     

Arsenic-safe alternate aquifers and their hydraulic characteristics in contaminated areas of Middle Ganga Plain, Eastern India

Arsenic groundwater contamination exceeding 0.05 mg/l affecting the Newer Alluvial tracts of Patna and Bhojpur, the two worst affected districts located in the Middle Ganga Plain in the Bihar state, has been studied The area is underlain by two-tier Quaternary aquifer systems within a depth of 300 m below ground level, separated by a 15?C32-m-thick clay/sandy clay aquitard. The upper part (<50 m depth) of the shallow aquifer system is arsenic-contaminated. The deeper aquifer system (lying below 120?C130 m depth) exhibits low arsenic load (max 0.0035 mg/l), having hydraulic conductivity between 64.88 and 82.04 m/day. Groundwater in the deeper aquifer occurs under semi-confined to confined condition due to poor hydraulic conductivity of the middle clay (4.7 × 10^???2???7.2 × 10^???3 m/day). Hydraulic head of the deeper aquifer remains close to the surface than the shallow aquifer. The two aquifer systems in the Newer Alluvium are replaced by a thick single aquifer system in the adjoining Older Alluvium, within a depth of 330 m below ground. In the arsenic-contaminated area, deeper aquifer is protected by a middle clay, which may be developed for community drinking water supply by deep tube wells having a yield capacity of 150 m³/h.     

Groundwater potential zoning of a peri-urban wetland of south Bengal Basin, India

Demand for groundwater for drinking, agricultural, and industrial purposes has increased due to rapid increase in population. Therefore, it is imperative to assess the groundwater potential of different areas, especially in a fragile wetland ecosystem to select appropriate sites for developing well fields to minimize adverse environmental impacts of groundwater development. This study considers East Calcutta Wetlands (ECW)??a freshwater peri-urban inland wetland ecosystem located at the lower part of the deltaic alluvial plain of South Bengal Basin and east of Kolkata city. This wetland is well known over the world for its resource recovery systems developed by local people through ages, using wastewater of the city. The subsurface geology is completely blanketed by the Quaternary sediments comprising a succession of silty clay, sand of various grades, and sand mixed with occasional gravels and thin intercalations of silty clay. Groundwater occurs mostly under confined condition except in those places where the top aquitard has been obliterated due to scouring action of past channels. The groundwater in the study area is being over-extracted at the rate of 65 × 10³ m³/day. Overlay analysis in Geographic Information System platform using multiple criteria such as water quality index, hydraulic conductivity, groundwater velocity, and depth to piezometric surface reveals that in and around ECW, there are five groundwater potential zones. About 74% of the aquifer of this area shows very poor to medium groundwater potential. Management options such as minimization of groundwater abstraction by introducing the treated surface water supply system and the implementation of rainwater harvesting and artificial recharge in high-rise buildings and industries are suggested for different potential zones.     

Groundwater vulnerability assessment using DRASTIC and Pesticide DRASTIC models in intense agriculture area of the Gangetic plains,India

Delineating areas susceptible to contamination from anthropogenic sources form an important component of sustainable management of groundwater resources. The present research aims at estimating vulnerability of groundwater by application of DRASTIC and Pesticide DRASTIC models in the southern part of the Gangetic plains in the state of Bihar. The DRASTIC and Pesticide DRASTIC models have considered seven parameters viz. depth to water level, net recharge, aquifer material, soil material, topography, impact of vadose zone and hydraulic conductivity. A third model, Pesticide DRASTIC LU has been adopted by adding land use as an additional parameter, to assess its impact on vulnerability zonation. The DRASTIC model indicated two vulnerable categories, moderate and high, while the Pesticide DRASTIC model revealed moderate, high and very high vulnerable categories. Out of the parameters used, depth to water level affected the vulnerability most. The parameter caused least impact was topography in DRASTIC, while in case of Pesticide DRASTIC and Pesticide DRASTIC LU models, the parameter was hydraulic conductivity. A linear regression between groundwater NO₃ concentrations and the vulnerability zonation revealed better correlation for Pesticide DRASTIC model, emphasising the effectiveness of the model in assessing groundwater vulnerability in the study region. Considering all three models, the most vulnerable areas were found to be concentrated mainly in two zones, (i) in the south-western part along Ekangarsarai-Islampur patch and (ii) around Biharsharif-Nagarnausa area in the central part. Both zones were characterised by intensive vegetable cultivation with urban areas in between.     

Neural Network Based Automated Interpretation Algorithm for Combined Geophysical Soundings in Coastal Zones

The fresh water availability in coastal aquifers is an important problem faced by a major part of world’s population dwelling there. For in situ and dynamic characterization of seawater encroachment into coastal aquifers, electrical geophysical methods are better suited. However, the resolving power of such data in distinguishing saline sands from moist clays in the subsurface is very poor. To meet this aspect and also the problem of analyzing voluminous data sets, we propose a feed forward back-propagation neural network (BPNN) based approach for the analysis of combined vertical electrical and induced polarization soundings. Our method is tested on synthetic data computed from available geo-electric sections and prevailing subsurface geological information of coastal aquifers of East Coast of India.The synthetic data are comprised of 18 combined soundings spread over five profiles. 15 out of 18 are used for training the BPNN, while 3 are used for testing. The trained BPNN (one node each in each of the input and output layers and 18 hidden nodes) showed 84.85% accuracy in testing phase for distinguishing clays from saline sands.Our method is also tested on real data concerning a shaly aquifer in Bahia, Brazil yielding an overall accuracy of 84.9%, comparable to that of synthetic case; thereby validating our approach.     

An appraisal of subsurface geology and groundwater resources of Owerri and environs based on electrical resistivity survey and borehole data evaluation

The research was aimed at determining the depth to the watertable, aquifer thickness and subsurface geology of the studyarea thus revealing its groundwater distribution as well as its potential as a substitute to the surface water resources.Vertical electrical soundings were carried out in the study area with maximum electrode spread. The Schlumberger electrode configuration technique was adopted. VES data were processedusing Schlumberger analysis package. Lithologic logs of alreadyexisting boreholes in the study area were collected, evaluatedand comparison were carried out. The results reveal alternatinglayers of sands, sandstones, gravel and clay. The lithologiclogs revealed that the study area is underlain by coastal sands (Benin formation). The water table varies from 10–64 m and thickness of the aquifer ranges from 20–80 m. Results show that the study area is underlain by a thick extensive aquifer that has a transmissivity of 2.8 × 10^-2m² s^-1 to 3.3 × 10^-1m² s^-1 and storativity 1.44 × 10^-4 to 1.68 × 10^-3m s^-1 values. The specific yield is about 0.31.The sandy component of the study area forms more than 90% ofthe sequence, therefore the permeability, the transmissivity andthe storage coefficient are high with an excellent source of groundwater resources.     

Aquifer Transmissivity and Basement Structure Determination Using Resistivity Sounding at Jos Plateau State Nigeria

A geoelectric investigation involving twenty-six vertical electrical soundings was carried out at Jos, Plateau State, North Central Nigeria. The survey was aimed at determining the structure of the underlying bedrock, as well as computing the transmissivity for the aquifer in the area. The basement geometrymapproduced from the results of the survey indicates that the bedrock is undulating lying at depths between 30 m to 6.5 m. There is evidence of faulting and fracturing within the area. Computation of aquifer transmissivity values based on the results obtained made it possible to demarcate regions with good ground water potential in the area.     

Groundwater Modeling in Agricultural Watershed under Different Recharge and Discharge Scenarios for Quaternary Aquifer Eastern Nile Delta,Egypt

Different scenarios of recharge and discharge were assessed for sustainable management of groundwater in Quaternary aquifer east of Nile Delta. MODFLOW was utilized to investigate the effect of land use change and damming construction in the upstream of the Nile River on the current and short-term groundwater management strategies. The interpretive transient simulation was performed between 2004 and 2016 after steady-state calibration in 2004, and transient state from 2004 to 2013 with different irrigation recharges associated with land use change in this period. Sensitivity analysis was performed for hydraulic conductivities, recharge, and conductance parameters. The predictive transient simulation was run till 2023 under three scenarios of increasing pumping rates by 15, 30, and 50% for agriculture expansion and specified head reduction of Port Said Canal by 0.2, 0.4, and 0.6 m associated with the reduction of Nile water levels after Grand Ethiopian Residence Dam, GERD operation in 2017. Results from the in- and out-flow budgets showed that groundwater aquifer is stable at the current rate of pumping till 2023. Groundwater heads decreased by 0.2 and 0.42 m in the southern section, and a slight increase in the northern part was noticed for the first and second scenarios, respectively. When additional pumping stress is applied (50% increase), groundwater head dropped by 0.66 m, and the storage is no longer able to maintain the aquifer capacity after 2020 (worst-case scenario).     

Determining the depositional pattern by resistivity–seismic inversion for the aquifer system of Maira area,Pakistan

Velocity and density measured in a well are crucial for synthetic seismic generation which is, in turn, a key to interpreting real seismic amplitude in terms of lithology, porosity and fluid content. Investigations made in the water wells usually consist of spontaneous potential, resistivity long and short normal, point resistivity and gamma ray logs. The sonic logs are not available because these are usually run in the wells drilled for hydrocarbons. To generate the synthetic seismograms, sonic and density logs are required, which are useful to precisely mark the lithology contacts and formation tops. An attempt has been made to interpret the subsurface soil of the aquifer system by means of resistivity to seismic inversion. For this purpose, resistivity logs and surface resistivity sounding were used and the resistivity logs were converted to sonic logs whereas surface resistivity sounding data transformed into seismic curves. The converted sonic logs and the surface seismic curves were then used to generate synthetic seismograms. With the utilization of these synthetic seismograms, pseudo-seismic sections have been developed. Subsurface lithologies encountered in wells exhibit different velocities and densities. The reflection patterns were marked by using amplitude standout, character and coherence. These pseudo-seismic sections were later tied to well synthetics and lithologs. In this way, a lithology section was created for the alluvial fill. The cross-section suggested that the eastern portion of the studied area mainly consisted of sandy fill and the western portion constituted clayey part. This can be attributed to the depositional environment by the Indus and the Kabul Rivers.     

Delineating shallow saline groundwater zones from Southern India using geophysical indicators

A geophysical survey was conducted over an industrial belt encompassing 80 functional leather factories in Southern India. These factories discharge untreated effluents which pollute shallow groundwater where electrical conductivity (EC) value had a wide range between 545 and 26,600 μS/cm (mean, 3, 901 μS/cm). The ranges of Na⁺ and Cl^? ions were from 46 to 4,850 mg/L (mean, 348 mg/L) and 25 to 10,390 mg/L (mean, 1,079 mg/L), respectively. Geoelectrical layer parameters of 37 vertical electrical soundings were analyzed to demarcate fresh and saline water zones. However, the analysis not did lead to a unique resolution of saline and fresh waters. It was difficult to assign a definitive value to the aquifer resistivity of any area. Thus, geophysical indicators, namely longitudinal unit conductance (S), transverse unit resistance (T), and average longitudinal resistivity (R _s), were calculated for identifying fresh and saline waters. Spatial distributions of S, T, and R _s reflected widely varying ranges for the saline and fresh water zones. Further, the empirical relation of formation factor (F) was established from pore-water resistivity and aquifer resistivity for fresh and saline aquifers, which may be used to estimate local EC values from the aquifer resistivity, where well water is not available.     

Delineation of groundwater development potential zones in parts of marginal Ganga Alluvial Plain in South Bihar, Eastern India

A part of the Gangetic Alluvial Plain covering 2,228 km², in the state of Bihar, is studied for demarcating groundwater development potential zones. The area is mainly agrarian and experiencing intensive groundwater draft to the tune of 0.12 million cubic metre per square kilometres per year from the Quaternary marginal alluvial deposits, unconformably overlain northerly sloping Precambrian bedrock. Multiparametric data on groundwater comprising water level, hydraulic gradient (pre- and post-monsoon), aquifer thickness, permeability, suitability of groundwater for drinking and irrigation and groundwater resources vs. draft are spatially analysed and integrated on a Geographical Information System platform to generate thematic layers. By integrating these layers, three zones have been delineated based on groundwater development potential. It is inferred that about 48% of the area covering northern part has high development potential, while medium and low development potential category covers 41% of the area. Further increase in groundwater extraction is not recommended for an area of 173 km², affected by over-exploitation. The replenishable groundwater resource available for further extraction has been estimated. The development potential enhances towards north with increase in thickness of sediments. Local deviations are due to variation of—(1) cummulative thickness of aquifers, (2) deeper water level resulting from localised heavy groundwater extraction and (3) aquifer permeability.     

A groundwater flow model for water resources management in the Ismarida plain,North Greece

This paper presents the development of a regional flow simulation model of the stream–aquifer system of Ismarida plain, northeastern Greece. It quantifies the water budget for this aquifer system and describes the components of groundwater and the characteristics of this system on the basis of results of a 3-year field study. The semiconfined aquifer system of Ismarida Lake plain consists of unconsolidated deltaic clastic sediments, is hydraulically connected with Vosvozis River, and covers an area of 46.75 km². The annual precipitation ranges in the study area from 270 to 876 mm. Eighty-seven irrigation wells are densely located and have been widely used for agricultural development. Groundwater flow in this aquifer was simulated with MODFLOW. Model calibration was done with observed water levels, and match was excellent. To evaluate the impacts of the current pumping schedule and propose solutions, four management scenarios were formulated and tested with the model. Based on model results, the simulated groundwater budget indicates that there must be approximately 33% decrease of withdrawals to stop the dramatic decline of groundwater levels. The application of these scenarios shows that aquifer discharge to the nearby river would be very low after a 20-year period.     

Assessment of the statistical significance of seasonal groundwater quality change in a karstic aquifer system near Izmir-Turkey

The main objective of this study was to statistically evaluate the significance of seasonal groundwater quality change and to provide an assessment on the spatial distribution of specific groundwater quality parameters. The studied area was the Mount Nif karstic aquifer system located in the southeast of the city of Izmir. Groundwater samples were collected at 57 sampling points in the rainy winter and dry summer seasons. Groundwater quality indicators of interest were electrical conductivity (EC), nitrate, chloride, sulfate, sodium, some heavy metals, and arsenic. Maps showing the spatial distributions and temporal changes of these parameters were created to further interpret spatial patterns and seasonal changes in groundwater quality. Furthermore, statistical tests were conducted to confirm whether the seasonal changes for each quality parameter were statistically significant. It was evident from the statistical tests that the seasonal changes in most groundwater quality parameters were statistically not significant. However, the increase in EC values and aluminum concentrations from winter to summer was found to be significant. Furthermore, a negative correlation between sampling elevation and groundwater quality was found. It was shown that with simple statistical testing, important conclusions can be drawn from limited monitoring data. It was concluded that less groundwater recharge in the dry period of the year does not always imply higher concentrations for all groundwater quality parameters because water circulation times, lithology, quality and extent of recharge, and land use patterns also play an important role on the alteration of groundwater quality.     

The depositional setting of the Late Quaternary sedimentary fill in southern Bannu basin,Northwest Himalayan fold and thrust belt,Pakistan

Geostatistical variogram and inversion techniques combined with modern visualization tools have made it possible to re-model one-dimensional electrical resistivity data into two-dimensional (2D) models of the near subsurface. The resultant models are capable of extending the original interpretation of the data to depict alluvium layers as individual lithological units within the 2D space. By tuning the variogram parameters used in this approach, it is then possible to visualize individual lithofacies and geomorphological features for these lithologic units. The study re-examines an electrical resistivity dataset collected as part of a groundwater study in an area of the Bannu basin in Pakistan. Additional lithological logs from boreholes throughout the area have been combined with the existing resistivity data for calibration. Tectonic activity during the Himalayan orogeny uplifted and generated significant faulting in the rocks resulting in the formation of a depression which subsequently has been filled with clay-silt and dirty sand facies typical of lacustrine and flood plain environments. Streams arising from adjacent mountains have reworked these facies which have been eroded and replaced by gravel-sand facies along channels. It is concluded that the sediments have been deposited as prograding fan shaped bodies, flood plain, and lacustrine deposits. Clay-silt facies mark the locations of paleo depressions or lake environments, which have changed position over time due to local tectonic activity and sedimentation. The Lakki plain alluvial system has thus formed as a result of local tectonic activity with fluvial erosion and deposition characterized by coarse sediments with high electrical resistivities near the mountain ranges and fine sediments with medium to low electrical resistivities towards the basin center.     

Time Series Analysis to Monitor and Assess Water Resources: A Moving Average Approach

An understanding of the behavior of the groundwater body and its long-term trends are essential for making any management decision in a given watershed. Geostatistical methods can effectively be used to derive the long-term trends of the groundwater body. Here an attempt has been made to find out the long-term trends of the water table fluctuations of a river basin through a time series approach. The method was found to be useful for demarcating the zones of discharge and of recharge of an aquifer. The recharge of the aquifer is attributed to the return flow from applied irrigation. In the study area, farmers mainly depend on borewells for water and water is pumped from the deep aquifer indiscriminately. The recharge of the shallow aquifer implies excessive pumping of the deep aquifer. Necessary steps have to be taken immediately at appropriate levels to control the irrational pumping of deep aquifer groundwater, which is needed as a future water source. The study emphasizes the use of geostatistics for the better management of water resources and sustainable development of the area.     

A Compartmental Model for Stable Time-Dependent Modeling of Surface Water and Groundwater

Water resource system planners make decisions that guide water management policy. The fundamental tools for assessing management and infrastructure strategies are combined hydro-economic models of river basins (RBHE models). These models have improved the economic efficiency of water use in situations of competition for scarce water resources. In RBHE models, a groundwater model is coupled with surface water models of the various water resources. Today, the groundwater models used in an RBHE model can be of two types: cell models or numerical models. Cell models are easy to use, but they are too simple to realistically describe the geology and hydrology of the area under investigation. Numerical models, in contrast, are closer to the physical behavior of the aquifer. However, the vast quantity of data to be analyzed makes them impractical for many management scenarios. Moreover, the calibrations of these high-resolution models are generally difficult and sensitive to the variation of parameters, especially when boundary conditions are dynamic. This is the case when dynamic river data or dynamic surface lake data are present. In this work, a compartmental cell model is built on the hydrogeology of the aquifer. In this approach, the hydrogeology of the aquifer and the dynamic boundary conditions are treated with separate models. A general mathematical formulation is presented where the calibration stage, the validation stage, and the prediction stage are formulated as a series of sub-model calibrations and solved using a general least squares routine. With this approach, it becomes possible to treat both the water level and the pumping rate in each cell as variables to be predicted. In most of the models, the pumping rates are known and the goals of the computation are to estimate the groundwater level. However, when for political or technical reasons access to some of the wells is difficult, the pumping rates are only partially known. Then, both groundwater levels and pumping rates are variables to be predicted by the groundwater model. A computer program was developed using MATLAB, with a Visual Basic graphical user interface using COM technology to access the advanced mathematical libraries. The approach was implemented with a real case study of the Yarkon–Taninim aquifer in Israel. The results indicate that the method is more stable than the classical approach.     

Hydrochemical Changes Over Time in the Zahedan Aquifer, Iran

Groundwater in the Zahedan Aquifer located in the arid southeast of Iran was chemically characterized to understand both the nature of the alluvium aquifer and the effect of human activities, specifically sewage disposal on groundwater quality. Concentrations of major cations and anions in the Zahedan Aquifer are much higher than concentrations observed in groundwater of similar settings. Although the nature of the aquifer and climatic parameters have affected the chemistry of groundwater, human impacts on the groundwater quality are more significant. The electrical conductivity in some areas of the Zahedan aquifer increased up to 7,500μS/cm in 25 years. Human influences in some areas are so prevalent, that the groundwater type changed from a Na⁺-HCO₃ ^- in 1976 to a Na⁺-Cl^- type in 2000. The impact of human influences on the groundwater quality is also indicated through observed nitrate concentration up to 4.81 meq/L in the urban area.     

Hydrochemical and water quality assessment of groundwater in Doon Valley of Outer Himalaya, Uttarakhand, India

The present study discusses ion sources and assesses the chemical quality of groundwater of Doon Valley in Outer Himalayan region for drinking and irrigational purposes. Valley is almost filled with Doon gravels that are main aquifers supplying water to its habitants. Recharged only by meteoric water, groundwater quality in these aquifers is controlled essentially by chemical processes occurring between water and lithology and locally altered by human activities. Seventy-six water samples were collected from dug wells, hand pumps and tube wells and were analysed for their major ion concentrations. The pH is varying from 5.6 to 7.4 and electrical conductivity from 71 to 951 μmho/cm. Groundwater of Doon valley is dominated by bicarbonate contributing 83% in anionic abundance while calcium and magnesium dominate in cationic concentrations with 88%. The seasonal and spatial variation in ionic concentration, in general, is related to discharge and lithology. The high ratio of (Ca + Mg)/(Na + K), i.e. 10, low ratio of (Na + K)/TZ+, i.e.0.2 and also the presence of carbonate lithology in the northern part of valley, is indicative of carbonate dissolution as the main controlling solute acquisition process in the valley. The low abundance of silica content and high HCO?/H?SiO? ratio also supports carbonate dissolution and less significant role of silicate weathering as the major source for dissolved ions in Doon Valley. The analytical results computed for various indices show that water is of fairly good quality, although, hard but have moderate dissolved solid content. It is free from sodium hazard lying in C?-S? and C?-S? class of USSL diagram and in general suitable for drinking and irrigation except few locations having slightly high salinity hazard.     

Integrated approach for demarcating subsurface pollution and saline water intrusion zones in SIPCOT area: a case study from Cuddalore in Southern India

This paper deals with a systematic hydrogeological, geophysical, and hydrochemical investigations carried out in SIPCOT area in Southern India to demarcate groundwater pollution and saline intrusion through Uppanar River, which flows parallel to sea coast with high salinity (average TDS 28, 870 mg/l) due to back waters as well as discharge of industrial and domestic effluents. Hydrogeological and geophysical investigations comprising topographic survey, self-potential, multi-electrode resistivity imaging, and water quality monitoring were found the extent of saline water intrusion in the south and pockets of subsurface pollution in the north of the study area. Since the area is beset with highly permeable unconfined quaternary alluvium forming potential aquifer at shallow depth, long-term excessive pumping and influence of the River have led to lowering of the water table and degradation of water quality through increased salinity there by generating reversal of hydraulic gradient in the south. The improper management of industrial wastes and left over chemicals by closed industries has led surface and subsurface pollution in the north of the study area.     

Groundwater modeling of Saq Aquifer Buraydah Al Qassim for better water management strategies

Saudi Arabia is an arid country. It has limited water supplies. About 80?C90% of water supplies come from groundwater, which is depleting day by day. It needs appropriate management. This paper has investigated groundwater modeling of Saq Aquifer in Buraydah Al Qassim to estimate the impact of its excessive use on depletion of Saq Aquifer. MODFLOW model has been used in this study. Data regarding the aquifer parameters was measured by pumping tests. Groundwater levels and discharge of wells in the area for the year 2008 and previous record of year 1999 have been collected from Municipal Authority of Buraydah. Location of wells was determined by Garmin. The model has been run for different sets of pumping rates to recommend an optimal use of groundwater resources and get prolonged life of aquifer. Simulations have been made for a long future period of 27?years (2008?C2035). Model results concluded that pumping from the Saq Aquifer in Buraydah area will result into significant cones of depression if the existing excessive pumping rates prevail. A drawdown up to 28?m was encountered for model run for 27?years for existing rates of pumping. Aquifer withdrawals and drawdowns will be optimal with the conservation alternative. The management scheme has been recommended to be adopted for the future protection of groundwater resources in Kingdom of Saudi Arabia.