Evaluation of methods for determination of hydraulic properties in an aquifer–aquitard system hydrologically connected to a river

Field and laboratory methods have been used to determine the hydraulic properties in a multiple-layer aquifer–aquitard system that is hydrologically connected to a river. First, hypothetical pumping tests in aquifer–aquitard systems were performed to evaluate the feasibility of MODFLOW-PEST in determining these parameters. Sensitivity analyses showed that: the horizontal hydraulic conductivity in the aquifer has the highest composite sensitivity; the vertical hydraulic conductivity has higher composite sensitivity than the horizontal hydraulic conductivity in the aquitard; and a partial penetration pumping well in an aquifer layer can improve the quality of the estimated parameters. This inverse approach was then used to analyze a pumping-recovery test conducted near the Platte River in southeastern Nebraska, USA. The hydraulic conductivities and specific yield were calculated for the aquitard and aquifer. The direct-push technique was used to generate sediment columns; permeameter tests on these columns produced the vertical hydraulic conductivities that are compatible with those obtained from the pumping-recovery test. Thus, the combination of the direct-push technique with permeameter tests provides a new method for estimation of vertical hydraulic conductivity. The hydraulic conductivity, determined from grain-size analysis, is smaller than the horizontal one but larger than the vertical one determined by the pumping-recovery test.     

Mapping an aquitard breach using shear-wave seismic reflection

In multi-layered hydrostratigraphic systems, aquitard breaches caused by faulting or paleo-erosion can allow substantial quantities of water of differing quality to be exchanged between aquifers. Seismic reflection technology was used to map the extent and orientation of an aquitard breach connecting a shallow alluvial aquifer to the deeper semi-confined Memphis aquifer in southwestern Tennessee, USA. Geophysical well logs indicate the presence of the aquitard at borehole locations that define the beginning and end points on two seismic survey lines, which intersect at a borehole where the aquitard is absent. A SE–NW-oriented paleochannel, 350 m wide and approximately 35–40 m deep, is interpreted from the seismic reflection surveys. The paleochannel cuts through the aquitard and into the upper part of the Memphis aquifer, thus creating a hydraulic connection between the shallow unconfined and deeper, semi-confined aquifers. The results indicate the potential of the shear-wave seismic reflection methods to resolve shallow breaches through fine-grained aquitards given availability of sufficient well control.     

Vertical movement of water in a High Plains Aquifer induced by a pumping well

Field observation and numerical simulations were carried out to evaluate the hydraulic relationship between the shallow and deep aquifer of a High Plains Aquifer system, in which shallow and deep aquifers are separated by an aquitard. Pumping from the lower aquifer resulted in a small drawdown in the upper aquifer and a larger drawdown in the aquitard; pumping from the shallow aquifer caused a small drawdown in the aquitard and the deep aquifer. Analysis of pumping test data gives the values of the hydraulic conductivity of the aquitard and the deep aquifer. Long-term observation of groundwater levels in the shallow and deep aquifers showed that a strong downward hydraulic gradient was maintained during an irrigation season. Numerical simulations were used to calculate the induced leakage of water from the shallow to the deep aquifer. Water budget analyses suggested that after pumping continues for a couple of days, the leakage from the overlying layers begins to supply the majority of the withdrawal from the deep aquifer. However, the induced leakage from the upper shallow aquifer can travel only a few meters into the aquitard, and it can not reach the lower aquifer during a 90 day pumping period. The major portion of the induced leakage occurred during the pumping period, but a small leakage can continue as a residual effect after the pumping period. The vertical hydraulic conductivity of the aquitard plays a major role in partitioning the ratio of the induced leakage for the pumping and after-pumping periods.     

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.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Groundwater response to dual tidal fluctuations in a peninsula or an elongated island

Groundwater flow in a peninsula or an elongated island is influenced by tidal fluctuations on both sides of the peninsula or island, which is named as dual tidal fluctuations. In this study, semianalytical solutions of transient groundwater flow in response to dual tidal fluctuations in an aquifer–aquitard system were presented for cases with and without the aquitard storage. These solutions were first derived using the Laplace transform and subsequently computed by the Fourier series numerical inverse Laplace transform. The derived solutions were found to agree very well with the results of numerical simulations by MODFLOW. The solution ignoring the aquitard storage approached the quasi‐steady state solution quickly when the mean sea level initial condition was used. The solutions with and without the aquitard storage were nearly the same at the early time and were separated from each other during the intermediate time, and the difference of solutions became constant at late time for small aquifer/aquitard storativity ratio and large tidal frequency. The propagation bias, which is the departure from the theoretical ratio of tidal attenuation to tidal lag, was enhanced not only with increase of the dimensionless specific leakage (aquitard/aquifer hydraulic conductance ratio) but also with decrease of the aquifer/aquitard storativity ratio and with the increase of the dimensionless tidal frequency. The solution with the aquitard storage was more sensitive to these three parameters. The newly developed solutions were capable of handling realistic initial conditions that might be approximated by piecewise linear functions. Copyright © 2012 John Wiley & Sons, Ltd.     

气压式振荡试验确定地层渗透系数适宜性

为评价滨海地区海（咸）水入侵区不同类型地层渗透性,提出振荡试验适宜模型选择方法。利用Geoprobe系统对不同平面位置地层实施了气压式振荡试验,部分点位进行了分层测试;利用配线法计算地层渗透系数。结果表明:当试验目标层为承压含水层时,利用Butler模型分析欠阻尼振荡试验,KGS模型分析过阻尼振荡试验,均可获得良好的计算结果;当试验目标层为弱透水层时,KGS模型可用于过阻尼振荡试验的数据分析,但计算结果会受到相邻高渗透性含水层以及钻孔编录不准确的影响。地层渗透性等级未知条件下,需根据地层类型、介质岩性、测试井结构和水位恢复过程曲线特征选择适宜的振荡试验分析模型,以准确获取地层渗透系数值。     

Influence of Hydraulic Conductivity and Wellbore Design in the Fate and Transport of Nitrate in Multi-aquifer Systems

Nitrate concentrations in multi-aquifer systems are heavily affected by the presence of wellbores (active or abandoned) that are screened in several aquifers. The spatial variability of hydraulic conductivity in the confining layers has also an important impact on the concentrations. A synthetic three-dimensional flow and transport exercise was carried in a multi-aquifer system consisting of two aquifers separated by an aquitard in which 100 vertical wellbores had been drilled. To model the wellbores and the flow and transport connection between aquifers that they may induce, we assign a high vertical hydraulic conductivity and a low effective porosity to the cell blocks including the wells. With these parameters, a solute will travel quickly from one aquifer to the other without being stored in the well itself. The wellbores will act as preferential pathways, and the solute will move quickly between aquifers according to the hydrodynamic conditions. Not considering these preferential pathways could induce erroneous interpretations of the solute distribution in an aquifer. We also noted that when there are vertical wellbores that connect aquifers in a multi-aquifer system, low conductivity in the aquitard enhances the flow of solute through the wellbores. Time-varying pumping rates induce important fluctuations in nitrate concentrations; therefore, any estimate of the water quality of the aquifer will depend on the moment when the data has been recorded. Consequently, concentration maps obtained by interpolation of point samples are seldom a good indicator of the chemical status of groundwater bodies; alternatively, we recommend complementing the usual interpolated maps with numerical models to gain a true understanding of the spatial distribution of the solute concentration.     

Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard

Recent work in southern Ontario, Canada, demonstrates anomalously high vertical groundwater flow velocities (>1 m/year) through a thick (as much as 60 m), sandy silt till aquitard (Northern till), previously assumed to be of very low permeability (hydraulic conductivity <10^–10 m/s). Rapid recharge is attributed to the presence of fractures and sedimentary heterogeneities within the till, but the field-scale flow regime is poorly understood. This study identifies the nature of physical groundwater pathways through the till and provides estimates of the associated groundwater fluxes. The aquitard groundwater flow system is characterized by integrating details of the outcrop and subsurface sedimentary characteristics of the till with field-based hydrogeologic investigation and numerical modeling. Outcrop and subsurface data identify a composite internal aquitard stratigraphy consisting of tabular till beds (till elements) separated by laterally continuous sheet-like sands and gravels (interbeds) and boulder pavements. Individual till elements contain sedimentary heterogeneities, including discontinuous sand and gravel lenses, vertical sand dikes, and zones of horizontal and vertical fractures. Hydrogeologic field investigations indicate a three-layer aquitard flow system, consisting of upper and lower zones of more hydraulically active and heterogeneous till separated by a middle unit of relatively lower hydraulic conductivity. Groundwater pathways and fluxes in the till were evaluated using a two-dimensional aquitard/aquifer flow model which indicates a step-wise flow mechanism whereby groundwater moves alternately downward along vertical pathways (fractures, sedimentary dikes) and laterally along horizontal sand interbeds within the till. This model is consistent with observed hydraulic-head and isotope profiles, and the presence of tritiated pore waters at various depths throughout the till. Simulations suggest that a bulk aquitard vertical hydraulic conductivity on the order of 1×10^–9 m/s is required to reproduce observed hydraulic-head and tritium profiles. Electronic Publication     

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology, aquitard heterogeneity is often less considered compared to aquifers, despite its significant impact on groundwater hydraulics and groundwater resources evaluation. A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity (K) and specific storage (Ss). The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or Ss. Two scenarios of analytical investigations are explored: one is the presence of aquitard interlayers with distinct K or Ss values, a common field-scale occurrence; another is an exponentially depth-decaying aquitard Ss, a regional-scale phenomenon supported by statistical analysis. Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion; a high-Ss interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion. Locations of low-K or high-Ss interlayers also significantly impact well hydraulics and groundwater budget. In the context of an exponentially depth-decaying aquitard Ss, a larger decay exponent can enhance aquifer drawdown. When using current models with a vertically homogeneous aquitard, half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard Ss should be used to calculate aquitard depletion and unconfined aquifer leakage.     

Tritium�Chelium groundwater age used to constrain a groundwater flow model of a valley-fill aquifer contaminated with trichloroethylene (Quebec, Canada)

Tritium?Chelium groundwater dating was carried out in a trichloroethylene (TCE)-contaminated valley-fill aquifer system in Quebec, Canada, where a numerical groundwater flow model was developed. Forty seven discrete groundwater and dissolved gas samples were obtained along two flow paths originating from known TCE source zones whose related plumes converge down gradient to form a single plume. Sampling points in monitoring wells were projected onto vertical sections showing particle tracks along the two flow paths. At these points, simulated advective ages obtained from particle tracking were matched to tritium?Chelium ages using different porosity values; the best match was 0.35. Ages were also obtained above and below a prodeltaic silty aquitard in a portion of the aquifer where some source zones are located, which provide groundwater and TCE transit times through the aquitard as well as a mean vertical hydraulic conductivity that agrees with previous estimates used in the model. In certain locations, anomalously old ages associated with high terrigenic ⁴He indicate areas where groundwater from the underlying proglacial unit flows upward into the deltaic sand aquifer through aquitard windows. Upflow locations correspond with increased TCE concentrations, suggesting significant TCE provenance through the proglacial unit originating from a previously unrecognized TCE source zone.     

Redox evolution in glacial drift aquifers: role of diamicton units in reduction of Fe(III)

High iron concentrations create water quality problems for municipal use in glacial drift aquifer units. The chemical evolution of oxic groundwater in shallow aquifer units to anoxic groundwater in deeper aquifer units, in which soluble Fe(II) is stable, is attributed to coupled reduction of Fe(III) on aquifer solids with oxidation of organic carbon. The objective of this study was to characterize the distribution of organic carbon in aquifer and aquitard sediments to determine the availability of potential electron donors to drive these reactions. To do this, four complete rotasonic cores in a glacial aquifer/aquitard system were sampled at close intervals for analyses of grain-size distribution and organic carbon content. The results indicate significantly higher organic carbon concentrations in diamicton (till) units that function as aquitards, relative to coarse-grained aquifer units. In addition, readily reducible iron content in the diamicton units and lower aquifer unit materials is sufficient to produce far more dissolved iron than is present in the aquifer. Groundwater evolves to the level of iron reduction as a terminal electron-accepting process as it moves downward through aquitard units along flow paths from upland recharge areas to downgradient discharge areas. Deeper aquifer units are therefore unlikely to contain groundwater with low iron concentration.     

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     

Geostatistical mapping of leakance in a regional aquitard, Oak Ridges Moraine area, Ontario, Canada

The Newmarket Till forms a regionally extensive aquitard separating two major aquifer systems in the Greater Toronto area, Canada. The till is incised, and sometimes eroded entirely, by a network of sand- and gravel-filled channels forming productive aquifers and, locally, high-conductivity windows between aquifer systems. Leakage through the till may also be substantial in places. This study investigates the spatial variability of aquitard leakance in order to assess the relative importance of recharge processes to the lower aquifers. With a large database derived from water-well records and containing both hard and soft information, the Sequential Indicator Simulation method is used to generate maps of aquitard thickness and window probability. These can be used for targeting channel aquifers and for identifying potential areas of recharge to the lower aquifers. Conductivities are modeled from sparse data assuming that their correlation range is much smaller than the grid spacing. Block-scale leakances are obtained by upscaling nodal values based on simulated conductivity and thickness fields. Under the "aquifer-flow' assumption, upscaling is performed by arithmetic spatial averaging. Histograms and maps of upscaled leakances show that heterogeneities associated with aquitard windows have the largest effect on regional groundwater flow patterns. Electronic Publication     

Logging in deep water wells in central Jutland, Denmark

During a search for unpolluted resources, new wells have been drilled to a deep confined aquifer in central Jutland. Since little is known about the extent and vulnerability of the deep reservoir, geophysical logging was used to supplement hydraulic well-testing. Gamma-ray, SP and resistivity logs were recorded in the boreholes. Gamma-logs were also recorded in two cased wells. Analysis of the logs, combined with information from samples, shows an unconfined aquifer of quartz and mica sand to a depth of about 55 m at Grindsted. Below this level there is a 40 m thick aquitard with an apparent southerly dip. The aquitard is composed of interbedded sand, silt and clay. The deep confined aquifer extends from 95 to 110 m below the surface at Grindsted. Both the aquitard and the deep aquifer are provisionally dated as Early/Middle Miocene. The electrical logs discriminate well between sand, silt and clay beds. If a suitable mud is used in the well it is possible to calculate the formation water resistivities in good agreement with water samples obtained after completion. From the resistivity logs formation factors were calculated in the range 3–5 using no-invasion departure curves. The high content of heavy mineralsin the Tertiary deposits makes it impossible to obtain precise lithological information from the gamma-ray log. On the other hand this log provides clear markers for identifying the hydrogeological units in the well-bores.     

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.     

山东省巨野煤田梁宝寺井田水文地质条件及充水因素分析

通过对巨野煤田梁宝寺井田含水层、隔水层及断层导水规律的分析,指出了对煤矿生产所造成的影响及生产中应注意的问题,认为梁宝寺井田开采上组煤的水文地质类型为裂隙、岩溶类简单至中等类型,开采下组煤的水文地质类型为岩溶裂隙类中等至复杂型。     

地面沉降对含水层参数及给水能力的影响研究

地面沉降主要由深层地下水开采造成,含水层的压缩释水是深层地下水开采量主要构成来源,同时随其压缩变形孔隙比减小造成储水系数、渗透系数的减小,对弱透水层非弹性释水量、越流补给量造成影响。地面沉降过程中的水文地质参数非线性变化及对承压含水层系统的反馈作用也成为水文地质领域的前沿问题。为对地下水开采量、沉降量、地质参数变化以及给水能力变化之间的关系做一个较为定量定性的探究,以含水层压缩过程中的物理机制为依据,并基于沧州地区深层承压含水层的地下水位变化过程和水文地质参数,采用数学手段构建储水系数和渗透系数变化的一维非线性沉降模型。模拟结果显示随承压水头下降,储水系数最高可减小77%,含水层系统的给水能力和储水能力会随压缩变形减小50%甚至更多。研究成果为深入认识水文地质参数与应力变化相关关系、科学评价承压含水层地下水储水调节能力有重要借鉴意义。     

基于现场双管试验确定弱透水层水力参数的方法

弱透水层水力参数对基础工程设计、建设和地下水资源评价、地下水污染预测具有重要意义。现提出一种现场确定含水层系统中弱透水层水力参数的钻孔双管试验方法。利用麻花钻钻孔,到达地下水位线以下揭露到目的层位后布设外管,以起到护壁防止土层塌孔的作用。在外管内通过锤击的方式将内管嵌入到目的层中,使其下部管中有一定厚度的土样作为现场固定试验土柱。待内外管中地下水位稳定后迅速加高内管水头,并通过马里奥特瓶装置保持定水头。建立钻孔双管试验土柱的数学模型,推导了土柱顶底面单位面积流量随时间变化的解析解,将其无量纲化后作出标准曲线。作为该方法应用的实例,在长江下游三角洲苏锡常冲积平原的常州天宁区德安医院附近的空旷场地进行了专门的钻孔试验。确定了该地区的细砂质粉土和粉质黏土的渗透系数和贮水率分别为：K=1.66×10?6 m/s,?s =8.21×10?3 m?1;K=6.06×10?7 m/s,?s =2.26×10?3 m?1。结果表明,该方法对于现场确定弱透水层的水力参数有较好的适用性。     

A generalized solution for groundwater head fluctuation in a tidal leaky aquifer system

A new analytical solution is developed for describing groundwater level fluctuations in a coupled leaky confined aquifer system which consists of an unconfined aquifer, confined aquifer, and an aquitard in between. The aquifer system has a tidal boundary at the seashore, a no flow boundary at remote inland side, and a confined aquifer extending under the sea and terminated with an outlet-capping. This new solution has shown to be a generalisation of most existing analytical solutions for a tidal aquifer system which includes single confined and leaky confined aquifers. In addition, the solution is used to explore the influences of the dimensionless leakance of the outlet-capping, the dimensionless hydraulic diffusivities, and the leakages of the inland and offshore aquitards on the head responses in the leaky confined aquifer.