Groundwater balance in the Khor Arbaat basin, Red Sea State, eastern Sudan

The Khor Arbaat basin is the main source of potable water supply for the more than 750,000 inhabitants of Port Sudan, eastern Sudan. The variation in hydraulic conductivity and storage capacity is due to the heterogeneity of the sediments, which range from clay and silt to gravely sand and boulders. The water table rises during the summer and winter rainy seasons; it reaches its lowest level in the dry season. The storage capacity of the Khor Arbaat aquifer is estimated to be 21.75?×?10⁶ m³. The annual recharge through the infiltration of flood water is about 1.93?×?10⁶ m³. The groundwater recharge, calculated as underground inflow at the ‘upper gate’, is 1.33?×?10⁵ m³/year. The total annual groundwater recharge is 2.06?×?10⁶ m³. The annual discharge through underground outflow at the ‘lower gate’ (through which groundwater flows onto the coastal plain) is 3.29?×?10⁵ m³/year. Groundwater discharge due to pumping from Khor Arbaat basin is 4.38?×?10⁶ m³/year on average. The total annual groundwater discharge is about 4.7?×?10⁶ m³. A deficit of 2.6?×?10⁶ m³/year is calculated. Although the total annual discharge is twice the estimated annual recharge, additional groundwater flow from the fractured basement probably balances the annual groundwater budget since no decline is observed in the piezometric levels.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

An assessment of the potential and impacts of winter water banking in the Sokh aquifer, Central Asia

The dynamics of artificial recharge of winter surface flows coupled with increased summer groundwater use for irrigation in the Sokh aquifer (Central Asia) have been investigated. Water release patterns from the giant Toktogul reservoir have changed, as priority is now given to hydropower generation in winter in Kyrgyzstan. Winter flows have increased and summer releases have declined, but the Syr Darya River cannot pass these larger winter flows and the excess is diverted to a natural depression, creating a 40?×?10⁹m³ lake. A water balance study of all 18 aquifers feeding the Fergana Valley indicated the feasibility of winter groundwater recharge in storage created by summer abstraction. This modeling study examines the dynamics of the process in one aquifer over a 5-year period, with four scenarios: the current situation; increased groundwater abstraction of around 625 million (M) m³/year; groundwater abstraction with an artificial recharge of 144 Mm³/year, equivalent to the volume available in low flow years in the Sokh River; and with a larger artificial recharge of 268 Mm³/year, corresponding to high flow availability. Summer surface irrigation diversions can be reduced by up to 350 Mm³ and water table levels can be lowered.     

雄安新区地下水资源概况、特征及可开采潜力

地下水资源在中国社会经济发展中发挥重要作用,特别是在地表水资源相对匮乏的北方地区。掌握一个地区地下水资源状况、动态变化特征及可开采潜力,对该地区的供水安全保障至关重要。本文选择雄安新区,在近年来开展的区域水文地质调查、监测及综合研究等成果基础上,结合前人研究,对雄安新区区域水文地质条件、地下水动态变化特征等进行分析总结;以恢复地下水降落漏斗为地下水可持续开采利用方案的目标,从白洋淀流域平原区尺度,设置现状开采条件、河流补水、工农业节水及地下水禁(限)采等不同情景方案,采用地下水数值模拟技术,综合分析不同情景30年后的预测结果,提出白洋淀流域平原区地下水可持续开采利用方案;在流域地下水可持续开采利用方案基础上,分析雄安新区地下水可开采的最大资源量,进而评价雄安新区地下水可开采潜力。结果显示,雄安新区区域水文地质条件相对简单,浅层富水性中等,深层富水性较强;地下水位为多年下降状态,近年来,浅、深层地下水整体呈企稳或回升状态,局部地区仍有所下降;地下水质量总体良好,且较为稳定。根据评价结果,雄安新区地下水可开采潜力约为1.80×10~8m~3/a,其中,浅层地下水可开采潜力约为1.50×10~8m~3/a,深层地下水可开采潜力约为0.30×10~8m~3/a。     

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin,northwest China

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?10⁴ m³/year in 2017 to 11,875?×?10⁴ m³/year in 2021, and to 17,039?×?10⁴ m³/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?10⁴ m³/year in 2017 to 1857?×?10⁴ m³/year in 2021, and to 510?×?10⁴ m³/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.     

Groundwater budget for the upper and middle parts of the River Gash Basin, eastern Sudan

The River Gash Basin is filled by the Quaternary alluvial deposits, unconformably overlying the basement rocks. The alluvial deposits are composed mainly of unconsolidated layers of gravel, sand, silt, and clays. The aquifer is unconfined and is laterally bounded by the impermeable Neogene clays. The methods used in this study include the carry out of pumping tests and the analysis of well inventory data in addition to the river discharge rates and other meteorological data. The average annual discharge of the River Gash is estimated to be 1,056?×?10⁶ m³ at El Gera gage station (upstream) and 587?×?10⁶ m³ at Salam-Alikum gage station (downstream). The annual loss mounts up to 40% of the total discharge. The water loss is attributed to infiltration and evapotranspiration. The present study proofs that the hydraulic conductivity ranges from 36 to 105 m/day, whereas the transmissivity ranges from 328 to 1,677 m²/day. The monitoring of groundwater level measurements indicates that the water table rises during the rainy season by 9 m in the upstream and 6 m in the midstream areas. The storage capacity of the upper and middle parts of the River Gash Basin is calculated as 502?×?10⁶ m³. The groundwater input reach 386.11?×?10⁶ m³/year, while the groundwater output is calculated as 365.98?×?10⁶ m³/year. The estimated difference between the input and output water quantities in the upper and middle parts of the River Gash Basin demonstrates a positive groundwater budget by about 20?×?10⁶ m³/year     

Application of the water balance model J2000 to estimate groundwater recharge in a semi-arid environment: a case study in the Zarqa River catchment, NW-Jordan

Pollution and overexploitation of scarce groundwater resources is a serious problem in the Zarqa River catchment, Jordan. To estimate this resource’s potential, the amount and spatial distribution of groundwater recharge was calculated by applying the hydrological model J2000. The simulation period is composed of daily values gathered over a 30-year period (July 1977 to June 2007). The figure finally obtained for estimated groundwater recharge of the Zarqa River catchment is 105 × 10⁶ m³ per year (21 mm a^?1). This is 19 % higher than the value previously assumed to be correct by most Jordanian authorities. The average ratio of precipitation to groundwater recharge is 9.5 %. To directly validate modelled groundwater recharge, two independent methods were applied in spring catchments: (1) alteration of stable isotope signatures (δ¹⁸O, δ²H) between precipitation and groundwater and (2) the chloride mass balance method. Recharge rates determined by isotopic investigations are 25 % higher, and recharge rates determined by chloride mass balance are 9 % higher than the modelled results for the corresponding headwater catchments. This suggests a reasonably modelled safe yield estimation of groundwater resources.     

Comparison of groundwater recharge estimation methods for the semi-arid Nyamandhlovu area, Zimbabwe

The Nyamandhlovu aquifer is the main water resource in the semi-arid Umguza district in Matebeleland North Province in Zimbabwe. The rapid increase in water demand in the city of Bulawayo has prompted the need to quantify the available groundwater resources for sustainable utilization. Groundwater recharge estimation methods and results were compared: chloride mass balance method (19–62 mm/year); water-table fluctuation method (2–50 mm/year); Darcian flownet computations (16–28 mm/year); ¹⁴C age dating (22–25 mm/year); and groundwater modeling (11–26 mm/year). The flownet computational and modeling methods provided better estimates for aerial recharge than the other methods. Based on groundwater modeling, a final estimate for recharge (from precipitation) on the order of 15–20 mm/year is believed to be realistic, assuming that part of the recharge water transpires from the water table by deep-rooted vegetation. This recharge estimate (2.7–3.6% of the annual precipitation of 555 mm/year) compares well with the results of other researchers. The advantages/disadvantages of each recharge method in terms of ease of application, accuracy, and costs are discussed. The groundwater model was also used to quantify the total recharge of the Nyamandhlovu aquifer system (20?×?10⁶–25?×?10⁶ m³/year). Groundwater abstractions exceeding 17?×?10⁶ m³/year could cause ecological damage, affecting, for instance, the deep-rooted vegetation in the area.     

Characterisation of groundwater recharge conditions and flow mechanisms in bedrock aquifers of the Johannesburg area,South Africa

Groundwater is the major source of water and a critical resource for socioeconomic development in semi-arid environments like the Johannesburg area. Environmental isotopes are employed in this study to characterise groundwater recharge and flow mechanisms in the bedrock aquifers of Johannesburg, which is known for polluted surface water. With the exception of boreholes near the Hartbeespoort Dam, groundwater in the study area was derived from meteoric water that has undergone some degree of evaporation before recharge, possibly via diffuse mechanisms. Boreholes that tap groundwater from the Transvaal Supergroup Formation show depletion in δ¹⁸O and δ²H values. This is attributed to diffuse recharge through weathering fractures at high elevation that are undergoing deep circulation or recharge from depleted rainfall from the high-latitude moisture sources. The influence of focused recharge from the Hartbeespoort Dam was observed in the boreholes north of the dam, possibly as a result of the north–south trending fault lines and the north-dipping fractures in the bedding planes of quartzites. This is also supported by a reservoir water budget method which indicated a mean annual net flux of 2,084,131 m³ from Hartbeespoort Dam recharging groundwater per annum. Using tritium in the dam and boreholes located at 750 m and 5400 m downstream, average groundwater flow velocity was estimated as 202 m/year. An open system was observed in shale, andesite and granitic-gneiss aquifers indicating soil CO₂ as a dominant source of carbon (δ¹³C) in groundwater. A closed system was also observed in dolomitic aquifers indicating carbonate dissolution as the predominant source of carbon.     

Assessment of village-wise groundwater draft for irrigation: a field-based study in hard-rock aquifers of central India

Extracted groundwater, 90% of which is used for irrigated agriculture, is central to the socio-economic development of India. A lack of regulation or implementation of regulations, alongside unrecorded extraction, often leads to over exploitation of large-scale common-pool resources like groundwater. Inevitably, management of groundwater extraction (draft) for irrigation is critical for sustainability of aquifers and the society at large. However, existing assessments of groundwater draft, which are mostly available at large spatial scales, are inadequate for managing groundwater resources that are primarily exploited by stakeholders at much finer scales. This study presents an estimate, projection and analysis of fine-scale groundwater draft in the Seonath-Kharun interfluve of central India. Using field surveys of instantaneous discharge from irrigation wells and boreholes, annual groundwater draft for irrigation in this area is estimated to be 212 × 10⁶ m³, most of which (89%) is withdrawn during non-monsoon season. However, the density of wells/boreholes, and consequent extraction of groundwater, is controlled by the existing hydrogeological conditions. Based on trends in the number of abstraction structures (1982–2011), groundwater draft for the year 2020 is projected to be approximately 307 × 10⁶ m³; hence, groundwater draft for irrigation in the study area is predicted to increase by ～44% within a span of 8 years. Central to the work presented here is the approach for estimation and prediction of groundwater draft at finer scales, which can be extended to critical groundwater zones of the country.     

Estimation of groundwater recharge in arid, data scarce regions; an approach as applied in the El Hawashyia basin and Ghazala sub-basin (Gulf of Suez, Egypt)

In this study, an approach for runoff and recharge estimations that can be applied in arid regions which suffer from lack of data is presented. Estimating groundwater recharge in arid regions is an extremely important but difficult task, the main reason is the scarcity of data in arid regions. This is true for the Eastern Egyptian Desert where groundwater is used for irrigation purposes in agricultural reclamation along the Red Sea coast line. As a result of the scarcity of hydrologic information, the relation between rainfall and runoff was calculated depending on the paleo-flood hydrology information. Two models were used to calculate the rainfall–runoff relationships for El Hawashyia basin and Ghazala sub-basin. Two computer programs known as Gerinne (meaning channel in German) and SMADA6 (Stormwater Management and Design Aid, version 6) were conjunctively used for this purpose. As a result of the model applied to El Hawashyia basin, a rainfall event of a total of 18.3 mm with duration 3 h at the station of Hurghada, which has an exceedance probability of 5–10 %, produces a discharge volume of 10.2 × 10⁶ m³ at the delta, outlet of the basin, as 4.7 mm of the rainfall infiltrates (recharge). For the Ghazala sub-basin, the model yields a runoff volume of 3.16 × 10⁶ m³ transferred from a total rainfall of 25 mm over a period of 3 h, as 3.2 mm of it was lost as infiltration.     

Groundwater resources assessment for irrigable agricultural lands in the Wadi Araba area, southern Jordan

The main target of this research paper was to the hydrogeological assessment of the groundwater resources to irrigate 600 ha of irrigable agricultural lands, distributed along the Dead Sea–Aqaba Highway in Umm, Methla, Wadi Musa, Qa’ Saideen and Rahma, southern Jordan. Therefore, a comprehensive groundwater study was commenced by drilling eight new wells which can be used to supply irrigable areas with the existing groundwater that would be enriched by the yield of three proposed recharge dams on Wadi Musa, Wadi Abu-Burqa, and Wadi Rahma. The evaluation of the pumping test data of the drilled was carried out using the standard methods of pumping test interpretation. This was based on the available water table measurements at well locations and knowledge of water flow in the general. The sustainable yield of each well was calculated based on the pumping test parameters. The obtained results indicate that pumping out of Beer Mathkor wells should not exceed 1,100 m³/day in the case of continuous pumping and 8,700 m³/day in the case of intermittent pumping. Since the water table did not significantly change with small changes in pumping (it took eightfolds of magnitude increase in pumping from approximately 1,100 to 8,700 m³/day to show a significant drop in the water table equivalent to about 5.5 MCM per year from the aquifer.     

Radiocarbon age distribution of groundwater in the Konya Closed Basin, central Anatolia, Turkey

Annual abstraction of 2.6?×?10⁹ m³ of groundwater in the 53,000 km² Konya Closed Basin of central Turkey has caused a head decline of 1 m/year over the last few decades. Therefore, understanding the hydrogeology of this large endorheic basin, in a semi-arid climate, is important to sustainable resource management. For this purpose, the groundwater’s radiocarbon age distribution has been investigated along a 150-km transect parallel to regional flow. Results show that the groundwater ranges in age from Recent at the main recharge area of the Taurus Mountains in the south, to about 40,000 years around the terminal Salt Lake located in the north. In this predominantly confined flow system, radiocarbon ages increase linearly by distance from the main recharge area and are in agreement with the hydraulic ages. The mean velocity of regional groundwater flow (3 m/year) is determined by the rate of regional groundwater discharge into the Salt Lake. Calcite dissolution, dedolomitization and geogenic carbon dioxide influx appear to be the dominant geochemical processes that determine the carbon isotope composition along the regional flow path. The groundwater’s oxygen-18 content indicates more humid and cooler paleorecharge. A maximum drop of 5°C is inferred for the past recharge temperature.     

Aquifer system response to intensive pumping in urban areas of the Gangetic plains,India: the case study of Patna

A significant component of domestic demand for water of urban areas located in the Gangetic plains is met by heavy pumping of groundwater. The present study is focused on the Patna municipal area, inhabited by 17 million people and spanning over 134 km², where entire urban water demand is catered from pumping by wells of various capacities and designs. The present study examines the nature of the aquifer system within the urban area, the temporal changes in the water/piezometric level and the recharge mechanism of the deeper aquifers. The aquifer system is made up of medium-to-coarse unconsolidated sand, lying under a ~40-m-thick predominantly argillaceous unit holding 8- to 13-m-thick localised sand layers and continues up to 220 m below ground. Groundwater occurs under semi-confined condition, with transmissivity of aquifers in 5,500–9,200 m² day^?1 range. Hydraulic head of the deeper aquifer remains in 9–19 m range below ground, in contrast to 1–9 m range of that of the upper aquitard zone. The estimated annual groundwater extraction from the deeper aquifer is ~212.0 million m³, which has created a decline of 3.9 m in the piezometric level of the deeper aquifer during the past 30 years. Unregulated construction of deep tube wells with mushrooming of apartment culture may further exacerbate the problem. The sand layers within the aquitard zone are experiencing an annual extraction of 14.5 million m³ and have exhibited stable water level trend for past one and half decades. This unit is recharged from monsoon rainfall, besides contribution from water supply pipe line leakage and seepage from unlined storm water drains.     

Methodology to evaluate the renewal period of carbonate aquifers: a key tool for their management in arid and semiarid regions,with the example of Becerrero aquifer,Spain

A methodological procedure is proposed for determining the renewal period (RP), which expresses the ratio of total storage to recharge of carbonate aquifers, and it was applied to the overexploited moderate-size Becerrero carbonate aquifer (southern Spain). To this end, geological and subsurface data—time domain electromagnetic (TEM) soundings and borehole logs—were integrated to construct a three-dimensional (3D) geological model of the aquifer. The interconnected porosity was estimated by analyzing 73 rock samples. The resulting 3D geometrical model makes it possible to quantify the fractions of the aquifer having a confined or unconfined behaviour. Based on the total storage capacity (179?·?10⁶–514?·?10⁶ m³) and available aquifer recharge estimation (4.8?·?10⁶–6.4?·?10⁶ m³/year), an RP between 37 and 106 years is obtained. In view of the RP, an exploitation rate slightly lower than the average recharge of the system is recommended, so that the piezometric level will be stable but below the discharge head that is produced through the springs in natural conditions. The proposed methodology to obtain an aquifer RP and the management strategies designed accordingly are of broad interest, especially for carbonate aquifers, which are abundant in arid and semiarid regions.     

In situ infiltration test using a reclaimed abandoned river bed: managed aquifer recharge in Shijiazhuang City,China

Managed aquifer recharge (MAR) is necessary for water resources management in arid and semiarid regions. Infiltration rate is often a decisive limiting factor in site selection for MAR. In order to avoid scale effects in the application of infiltration rate parameters, the largest in situ infiltration test in China was undertaken between August 19 and August 30, 2009 to measure the infiltration rate of the field selected for MAR in Shijiazhuang City, China. The in situ test lasted for 10 days, and about 1.82 × 10⁷ m³ of water was introduced into the infiltration field. Groundwater level variations were monitored during the test. Monitoring showed that the infiltration rate of surface water was 1.5 m/day, which means that about 10–15 × 10⁸ m³/a of water could be injected into the target aquifer. Also, groundwater level variations showed that the northern part of the infiltration field had a higher infiltration rate, as predicted, and the test result supplied a sound foundation for validation of the groundwater numerical simulation, which will be of benefit for future predictions of the response of the groundwater level to artificial recharge engineering. Finally, an artificial recharge plan was proposed based on the infiltration test results and the water source conditions, which would be useful for the development of MAR programs and management of local water resources.     

Hydrogeology of the Kabul Basin (Afghanistan), part I: aquifers and hydrology

Shallow groundwater represents the main source for water supply in Kabul, Afghanistan. Detailed information on the hydrogeology of the Kabul Basin is therefore needed to improve the current supply situation and to develop a sustainable framework for future groundwater use. The basin is situated at the intersection of three major fault systems of partially translational and extensional character. It comprises three interconnected aquifers, 20–70 m thick, consisting of coarse sandy to gravely detritus originating from the surrounding mountains. The aquifers were deposited by three rivers flowing through the basin. The coarse aquifer material implies a high permeability. Deeper parts are affected by cementation of pore spaces, resulting in formation of semi-diagenetic conglomerates, causing decreased well yields. Usually the aquifers are covered by low-permeability loess which acts as an important protection layer. The main groundwater recharge occurs after the snowmelt from direct infiltration from the rivers. The steadily rising population is estimated to consume 30–40 million m³ groundwater per year which is contrasted by an estimated recharge of 20–45 million m³/a in wet years. The 2000–2005 drought has prevented significant recharge resulting in intense overexploitation indicated by falling groundwater levels.     

Large-scale replenishment of groundwater depletion by long-lived extreme rainstorms in arid regions: case study of the 2019 floods in Iran

The Dena rainstorm in Iran in March and April 2019 caused about US$ 8.3?×?10⁹ damage in the country; however, it resulted in the replenishment of half of the dam reservoirs and 35% of ponds and lakes. Also, it increased the volume of groundwater stored in aquifers by 3.6?×?10⁹ m³. In arid and semiarid regions such as most parts of Iran, which usually face water scarcity, getting water from rainstorms is essential for replenishing water resources. This research aims to quantify the direct and indirect effects of the Dena rainstorm on the replenishment of Iran’s groundwater storage using the groundwater balance method and water-table fluctuation method. Studies showed that the main mechanisms for replenishment of groundwater storage due to the rainstorm included increases in precipitation recharge, surface runoff recharge, and artificial recharge, and reductions in irrigation withdrawal and evapotranspiration, while the contribution of each factor is estimated to be about 23, 28, 2, 15, and 32%, respectively.

     

Residence times of shallow groundwater in West Africa: implications for hydrogeology and resilience to future changes in climate

Although shallow groundwater (<50 mbgl) sustains the vast majority of improved drinking-water supplies in rural Africa, there is little information on how resilient this resource may be to future changes in climate. This study presents results of a groundwater survey using stable isotopes, CFCs, SF₆, and ³H across different climatic zones (annual rainfall 400–2,000 mm/year) in West Africa. The purpose was to quantify the residence times of shallow groundwaters in sedimentary and basement aquifers, and investigate the relationship between groundwater resources and climate. Stable-isotope results indicate that most shallow groundwaters are recharged rapidly following rainfall, showing little evidence of evaporation prior to recharge. Chloride mass-balance results indicate that within the arid areas (<400 mm annual rainfall) there is recharge of up to 20 mm/year. Age tracers show that most groundwaters have mean residence times (MRTs) of 32–65 years, with comparable MRTs in the different climate zones. Similar MRTs measured in both the sedimentary and basement aquifers suggest similar hydraulic diffusivity and significant groundwater storage within the shallow basement. This suggests there is considerable resilience to short-term inter-annual variation in rainfall and recharge, and rural groundwater resources are likely to sustain diffuse, low volume abstraction.     

Groundwater Flow Simulation and its Application in Groundwater Resource Evaluation in the North China Plain,China

The purpose of this study is to establish a 3D groundwater flow modelling for evaluating groundwater resources of the North China Plain.First,the North China Plain was divided into three aquifers vertically through a characterization of hydrogeological conditions.Groundwater model software GMS was used for modeling to divide the area of simulation into a regular network of 164 rows and 148 lines.This model was verified through fitting of the observed and the simulated groundwater flow fields at deep and shallow layers and comparison between the observed and simulated hydrographs at 64 typical observation wells.Furthermore,water budget analysis was also performed during the simulation period(2002-2003).Results of the established groundwater flow model showed that the average annual groundwater recharge of the North China Plain during 1991 to 2003 was 256.68×10~8 m~3/yr with safe yield of groundwater resources up to 213.49×10~8 m~3/yr,in which safe yield of shallow groundwater and that of deep groundwater was up to 191.65×10~8 m~3/yr and 22.64×10~8 m~3/yr respectively.Finally,this model was integrated with proposal for groundwater withdrawal in the study area after commencement of water supply by South-North Water Transfer Project,aiming to predict the changing trend of groundwater regime.As indicated by prediction results,South-North Water Transfer Project,which is favorable for effective control of expansion and intensification of existing depression cone,would play a positive role in alleviation of short supply of groundwater in the North China Plain as well as maintenance and protection of groundwater.