Application of Numerical Modeling for Groundwater Flow System Analysis in the Akaki Catchment,Central Ethiopia

A three dimensional steady-state finite difference groundwater flow model is used to quantify the groundwater fluxes and analyze the subsurface hydrodynamics in the Akaki catchment by giving particular emphasis to the well field that supplies water to the city of Addis Ababa. The area is characterized by Tertiary volcanics covered with thick residual and alluvial soils. The model is calibrated using head observations from 131 wells. The simulation is made in a two layer unconfined aquifer with spatially variable recharge and hydraulic conductivities under well-defined boundary conditions. The calibrated model is used to forecast groundwater flow pattern, the interaction of groundwater and surface water, and the effect of pumping on the well field under different scenarios. The result indicates that the groundwater flows regionally to the south converging to the major well field. Reservoirs and rivers play an important role in recharging the aquifer. Simulations made under different pumping rate indicate that an increase in pumping rate results in substantial regional groundwater level decline, which will lead to the drying of springs and shallow hand dug wells. Also, it has implications of reversal of flow from contaminated rivers into productive aquifers close to main river courses. The scenario analysis shows that the groundwater potential is not enough to sustain the ever-growing water demand of the city of Addis Ababa. The sensitivity and scenario analysis provided important information on the data gaps and the specific sites to be selected for monitoring, and may be of great help for transient model development. This study has laid the foundation for developing detailed predictive groundwater model, which can be readily used for groundwater management practices.     

Postaudit evaluation of conceptual model uncertainty for a glacial aquifer groundwater flow and contaminant transport model

Numerical groundwater flow and contaminant transport modeling incorporating three alternative conceptual models was conducted in 2005 to assess remedial actions and predict contaminant concentrations in an unconfined glacial aquifer located in Milford, Michigan, USA. Three alternative conceptual models were constructed and independently calibrated to evaluate uncertainty in the geometry of an aquitard underlying the aquifer and the extent to which infiltration from two manmade surface water bodies influenced the groundwater flow field. Contaminant transport for benzene, cis-DCE, and MTBE was modeled for a 5-year period that included a 2-year history match from July 2003 to May 2005 and predictions for a 3-year period ending in July 2008. A postaudit of model performance indicates that predictions for pumping wells, which integrated the transport signal across multiple model layers, were reliable but unable to differentiate between alternative conceptual model responses. In contrast, predictions for individual monitoring wells with limited screened intervals were less consistent, but held promise for evaluating alternative hydrogeologic models. Results of this study suggest that model conceptualization can have important practical implications for the delineation of contaminant transport pathways using monitoring wells, but may exert less influence on integrated predictions for pumping wells screened over multiple numerical model layers.     

A conceptual and numerical model for groundwater management: a case study on a coastal aquifer in southern Tuscany, Italy

Ongoing hydrogeological research aims to develop a correct management model for the Plio-Pleistocene multi-aquifer system of the Albegna River coastal plain (southern Tuscany, Italy); overexploitation of this aquifer for irrigation and tourism has caused seawater intrusion. The conceptual model is based on field and laboratory data collected during the 1995–2003 period. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer. Natural outflow occurs through a diffuse flow into the sea and river; artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage, closing of aquitard) and artificially (multiscreened wells). The aquifer was represented by a three-dimensional finite element model using the FEFLOW numerical code. The model was calibrated for steady-state and transient conditions by matching computed and measured piezometric levels (February 1995–February 1996). The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it occurs above all in the Osa-Albegna sector, as well as along the river that at times feeds the aquifer. The effects of hypothetical aquifer exploitation were assessed in terms of water budget and hydraulic head evolution.     

The dependence of stream depletion by seasonal pumping on various hydraulic characteristics and engineering factors

Compensation pumping is used to alleviate deficiencies in streamflow discharge during dry seasons. Short-term groundwater pumping can use aquifer storage instead of catchment-zone water until the drawdown reaches the edge of the stream. The capacitance is a complex, dimensionless parameter of an aquifer system that defines the delayed effect on streamflow when there is groundwater pumping. This parameter is a function of aquifer hydraulic characteristics, pumping time, and distance between the well and stream edge; the latter can involve stream leakance and vertical leakance of an associated aquitard. Three typical hydraulic cases of combined water systems (major catchment-zone wells close to the stream and compensation pumping wells) were classified depending on their capacitance structure (i.e. the relationship between surface water and groundwater): (1) perfect hydraulic connection between the stream and aquifer; (2) imperfect hydraulic connection between the stream and aquifer; and (3) essentially imperfect hydraulic connection between the stream and the underlying confined aquifer. The impact of various hydraulic characteristics and engineering factors on stream depletion was examined by conceptual and numerical modeling. To predict the suitability and efficiency of a combined water system application, regression tests were undertaken on unit stream depletion and capacitance, and power dependencies were defined.     

Evaluation and numerical modeling of seawater intrusion in the Gaza aquifer (Palestine)

A numerical assessment of seawater intrusion in Gaza, Palestine, has been achieved applying a 3-D variable density groundwater flow model. A two-stage finite difference simulation algorithm was used in steady state and transient models. SEAWAT computer code was used for simulating the spatial and temporal evolution of hydraulic heads and solute concentrations of groundwater. A regular finite difference grid with a 400 m² cell in the horizontal plane, in addition to a 12-layer model were chosen. The model has been calibrated under steady state and transient conditions. Simulation results indicate that the proposed schemes successfully simulate the intrusion mechanism. Two pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and solute concentrations in the groundwater for a planning period of 17 years. The results show that seawater intrusion would worsen in the aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area, to moderate pumpage rates from water supply wells far from the seashore and to increase the aquifer replenishment by encouraging the implementation of suitable solutions like artificial recharge, may limit significantly seawater intrusion and reduce the current rate of decline of the water levels.     

Regional hydrostratigraphy and groundwater flow modeling in the arsenic-affected areas of the western Bengal basin, West Bengal, India

The first documented interpretation of the regional-scale hydrostratigraphy and groundwater flow is presented for a ～21,000-km² area of the arsenic-affected districts of West Bengal [Murshidabad, Nadia, North 24 Parganas and South 24 Parganas (including Calcutta)], India. A hydrostratigraphic model demonstrates the presence of a continuous, semi-confined sand aquifer underlain by a thick clay aquitard. The aquifer thickens toward the east and south. In the south, discontinuous clay layers locally divide the near-surface aquifer into several deeper, laterally connected, confined aquifers. Eight 22-layer model scenarios of regional groundwater flow were developed based on the observed topography, seasonal conditions, and inferred hydrostratigraphy. The models suggest the existence of seasonally variable, regional, north–south flow across the basin prior to the onset of extensive pumping in the 1970s. Pumping has severely distorted the flow pattern, inducing high vertical hydraulic gradients across wide cones of depression. Pumping has also increased total recharge (including irrigational return flow), inflow from rivers, and sea water intrusion. Consequently, downward flow of arsenic contaminated shallow groundwater appears to have resulted in contamination of previously safe aquifers by a combination of mechanical mixing and changes in chemical equilibrium.     

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

Groundwater flow and capture zone analysis of the Central Passaic River Basin, New Jersey

Delineating capture zones of pumping wells is an important part of safe drinking water and well protection programs. Capture zones or contributing areas of a groundwater extraction well are the parts of the aquifer recharge areas from which the wells draw their water. Their extent and location depend on the hydrogeologic conditions such as groundwater recharge, pumping scenario and the aquifer properties such as hydraulic conductivity, porosity, heterogeneity of the medium and hydraulic gradient. Different methods of delineation can be used depending on the complexity of the hydrogeologic conditions. In this study, a 3-dimensional transient numerical MODFLOW model was developed for the Central Passaic River Basin (CPRB), and used with a MODPATH particle tracking code to determine 3-dimensional transient capture zones. Analytically calculated capture zones from previous studies at the site were compared with the new numerically simulated capture zones. The study results revealed that the analytical solution was more conservative, estimating larger capture zones than the numerical models. Of all the parameters that can impact the size, shape and location of a capture zone, the hydraulic conductivity is one of the most critical. Capture zones tend to be smaller in lower hydraulic conductivity areas.     

Groundwater modeling in semiarid Central Sudan: adequacy and long-term abstraction

The present study assesses groundwater resources in the semiarid central Sudan, where 20 deep productive wells were installed to supply a major city in the region, El Obeid. The wells, which has an average 20 L/s discharge each, are taping a deep semiconfined to confined aquifer of fluvial silisiclastics deposited in the Tertiary–Pleistocene. Groundwater modeling was used as a technique to interpret the hydrologic system in arid to semiarid central Sudan and to simulate the future influence of the project on the hydrogeologic system. The simulation confirmed that steady-state flow conditions have been currently reached as indicated by consistency of computed heads. It also calibrated the values of the conductivity and recharge and ensured the sustainability of the El Obeid water supply project. A total of 3.5 × 10⁷ m³/year can be continually extracted from the deep aquifer to supply El Obeid city without endangering the groundwater resources in the region. The decline in water level will not exceed 25 m during the first 10 years, while indefinite continuous pumping will affect only the vicinity of the wells in a circle of 30 km diameter. Therefore, aquifer storage capacity and hydraulic properties encourage further groundwater exploitation. The present use of groundwater is extremely lower than the present demand, and it can potentially cover future demands without introducing significant changes to the system. The increase of pumping cost due to the decline in head subsequent to project operation was found to be minimal and of local effect.     

水资源和环境工程中水平井研究简介

简单介绍了水平井在水资源和环境工程中的近期研究进展, 重点介绍了在不同含水层(如承压含水层、潜水含水层、越流含水层及河下含水层) 中, 水平井抽水条件下降深的半解析解, 同时分析了降深随抽水时间变化的标准曲线和微分标准曲线.这些半解析解可用于分析小流量水平井在中长时段的降深特性.分析了用于排水和供水的大流量水平井的水力学特征, 并介绍了求解渗流-管流耦合井流系统的一种新方法.同时介绍了在非饱和含水层中水平井抽取气体的动力学特征, 并分析了地面覆盖和未覆盖2种情况下的气体清除效率, 报告了在不同饱和含水层条件下水平井的捕获区和捕获时间的计算.      

河北省滦县李家屿灰场2#副坝区灰水运动的水文地质评价

韩家哨村位于河北滦县李家屿灰场北部谷地,地势低于灰场,距离约800m。灰场运行期间,该区地下水出现水位升高及水质变差等问题。本文从水文地质条件角度,详细研究灰水的运移途径及其对该区地下水的影响。韩家哨村区域地下水系统分为上层滞水和孔隙一裂隙（溶隙）潜水含水层系统。潜水含水层系统水力传导系数为172．6～203．4m^2／d。地下水补给主要来源于韩家哨村以南坝后地形较高区域及灰水。灰水渗漏优势区韩家哨村南部宽约200m的带状区域。灰场堆灰形成新的地下水分水岭致使灰水在20副坝坝前垂直下渗进入风化带,沿白云岩裂隙（溶隙）补给潜水含水层,并向韩家哨村区域流动补给孔隙含水层,改变后的潜水部分通过民井向外排泄。     

基于解析法和数值法的非稳定流抽水试验参数反演

含水层的水文地质参数是进行地下水资源计算、地下水污染防控等所必需的基础数据,结合数值模拟技术进行含水层参数反演很有必要.按照1:5万水文地质调查规范在江汉平原仙桃市杨林尾镇复兴水厂不同含水层位开展抽水试验,包括深层含水层单孔抽水试验以及浅层含水层中群孔(2孔)抽水试验.对于单孔抽水试验,应用第1类越流系统井流理论进行参数反演;对于群孔抽水试验,推导了特定综合井函数,并利用特定标准曲线匹配法和直线图解法求解了含水层参数.随后利用FEFLOW软件建立了相应数值模型,拟合了含水层参数.结果表明:浅层含水层的渗透系数变化范围为21.66~54.00 m/d,贮水率变化范围为1.28×10-5~8.00×10-4 m-1;深层含水层渗透系数变化范围为1.27~7.00 m/d,贮水率变化范围为3.90×10-6~5.00×10-6 m-1.对于深层承压含水层而言,越流补给量较大.采用数值模拟方法结合抽水试验数据求参,综合考虑了含水层结构,拟合效果好,所得结果更加可靠.      

Effects of model layer simplification using composite hydraulic properties

The effects of simplifying hydraulic property layering within an unconfined aquifer and the underlying confining unit were assessed. The hydraulic properties of lithologic units within the unconfined aquifer and confining unit were computed by analyzing the aquifer-test data using radial, axisymmetric two-dimensional (2D) flow. Time-varying recharge to the unconfined aquifer and pumping from the confined Upper Floridan aquifer (USA) were simulated using 3D flow. Conceptual flow models were developed by gradually reducing the number of lithologic units in the unconfined aquifer and confining unit by calculating composite hydraulic properties for the simplified lithologic units. Composite hydraulic properties were calculated using either thickness-weighted averages or inverse modeling using regression-based parameter estimation. No significant residuals were simulated when all lithologic units comprising the unconfined aquifer were simulated as one layer. The largest residuals occurred when the unconfined aquifer and confining unit were aggregated into a single layer (quasi-3D), with residuals over 100% for the leakage rates to the confined aquifer and the heads in the confining unit. Residuals increased with contrasts in vertical hydraulic conductivity between the unconfined aquifer and confining unit. Residuals increased when the constant-head boundary at the bottom of the Upper Floridan aquifer was replaced with a no-flow boundary.     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

基于多源数据的弱透水层水文地质参数反演研究

弱透水层是含水层系统的重要组成部分,其水文地质参数的准确获取一直以来都是研究难点,传统室内试验难以克服应力状态及环境变化对土层参数的影响,因此如何在原位状态下评价其水文地质参数有着重要意义。选取呼和浩特盆地某淤泥层为研究对象,基于抽水试验和溶质运移试验获得的多源数据,采用基于控制体积有限差分法的非结构化网格（Unstructured Grid）方法,建立三维地下水流-溶质运移耦合模型,对弱透水层水文地质参数进行反演。结果显示:模型计算的抽水井中水位和溶质浓度变化规律与实际观测数据能够较好拟合,模型识别出的抽水含水层渗透参数（4.8 m/d）与通过解析法计算得出的结果（4.17 m/d）相近,建立的数值模型符合实际水文地质条件;反演得到淤泥质弱透水层垂向渗透系数为1.2×10?4 m/d、储水率为1.0×10?5 m?1。基于多源数据的参数反演方法可为弱透水层参数研究提供借鉴。     

Integrated approach of geospatial visualization and modeling for groundwater management in hard rock terrains in Nagpur Urban Area,India

Groundwater management requires complete visualization of aquifer characteristics to understand scientific aspects and hence remains a challenge, especially in hard rock terrains. In the present research paper, a comprehensive approach using 3D stratigraphic model, fence diagrams and numerical modeling has been proposed to understand the groundwater status for effective recharge. The stratigraphy, groundwater flow, and groundwater fluctuations for the period 1999–2010 were analyzed. The total volume of formations, volume of voids, storage capacity, and quantities of recharge of unconfined aquifer system in the Nagpur urban area were estimated. The steady state groundwater flow model of Basalt formation was calibrated to evaluate the subsurface system using Processing Modflow (PMWIN 5.3.2). The calibrated hydraulic head is compared with field observed head. The comparative spatial analysis presents a simple integrated approach in identifying zones with falling groundwater trends suitable for groundwater recharge in hard rock terrain in Nagpur urban area.     

Seawater intrusion in Syrian coastal aquifers, past, present and future, case study

Overexploitation of shallow aquifers on the Syrian coast, north of Latakia (Damsarkho) for irrigation and tourism has caused an intrusion of seawater. The seawater intrusion into this aquifer has been presented by a three-dimensional finite element model using the FEFLOW numerical code. This conceptual model is based on field and laboratory data collected during the period 1966–2003. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer; natural outflow occurs through a diffuse flow into the sea; and artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage) and artificially (multi-screened wells). The model was calibrated for transient conditions. The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it mainly occurs in the Damsarkho plain. The effects of hypothetical aquifer exploitation were assessed in terms of salt budget.     

Groundwater recharge in suburban areas of Hanoi,Vietnam: effect of decreasing surface-water bodies and land-use change

Over-exploited groundwater is expected to remain the predominant source of domestic water in suburban areas of Hanoi, Vietnam. In order to evaluate the effect on groundwater recharge, of decreasing surface-water bodies and land-use change caused by urbanization, the relevant groundwater systems and recharge pathways must be characterized in detail. To this end, water levels and water quality were monitored for 3 years regarding groundwater and adjacent surface-water bodies, at two typical suburban sites in Hanoi. Stable isotope (δ¹⁸O, δD of water) analysis and hydrochemical analysis showed that the water from both aquifers and aquitards, including the groundwater obtained from both the monitoring wells and the neighboring household tubewells, was largely derived from evaporation-affected surface-water bodies (e.g., ponds, irrigated farmlands) rather than from rivers. The water-level monitoring results suggested distinct local-scale flow systems for both a Holocene unconfined aquifer (HUA) and Pleistocene confined aquifer (PCA). That is, in the case of the HUA, lateral recharge through the aquifer from neighboring ponds and/or irrigated farmlands appeared to be dominant, rather than recharge by vertical rainwater infiltration. In the case of the PCA, recharge by the above-lying HUA, through areas where the aquitard separating the two aquifers was relatively thin or nonexistent, was suggested. As the decrease in the local surface-water bodies will likely reduce the groundwater recharge, maintaining and enhancing this recharge (through preservation of the surface-water bodies) is considered as essential for the sustainable use of groundwater in the area.

     

The relationship between the water table and the surface flow of a losing stream, lower Medano Creek, Great Sand Dunes National Monument, Colorado

 Proposed groundwater withdrawals in the San Luis Valley of Colorado may lower the water table in Great Sand Dunes National Monument. In response, the National Park Service initiated a study that has produced a generalized conceptual model of the hydrologic system in order to assess whether a lowering of the water table might decrease the surface flow of lower Medano Creek. Based upon information obtained during the drilling of several boreholes, there appear to be five important hydrostratigraphic units underlying lower Medano Creek within the upper 30 m of the ground surface: 1. a perched aquifer overlying an aquitard located between about 5 and 6 m below the ground surface; 2. the aquitard itself; 3. an unconfined aquifer located between the upper and lower aquitards; 4. an aquitard located between about 27 and 29 m below the ground surface; and 5. a confined underlying the lower aquitard. Because the areal extent of the aquitards cannot be determined from the borehole data, a detailed conceptual model of the hydrogeologic system underlying lower Medano Creek cannot be developed. However, a generalized conceptual model can be envisioned that consists of a complex system of interlayered aquifers and leaky aquitards, with each aquifer having a unique hydraulic head. Water levels in the perched aquifer rise rapidly to their annual maximum levels in response to the arrival of the flow terminus of Medano Creek during the spring runoff event, and the location of the flow terminus is directly dependent upon the discharge of the creek. Water levels in the deeper, non-perched aquifers do not appear to fluctuate significantly in response to the arrival of the flow terminus, demonstrating that it is unlikely that the proposed groundwater withdrawals will decrease the surface flow of lower Medano Creek. Received: 27 December 1995 · Accepted: 20 February 1996     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.