Groundwater Resources Potential in Hard Rock Terrain: A Multivariate Approach

Groundwater resources are limited and difficult to predict in crystalline bedrock due to heterogeneity and anisotropy in rock fracture systems. Municipal‐level governments often lack the resources for traditional hydrogeological tests when planning for sustainable use of water resources. A new methodology for assessing groundwater resources potential (GRP) based on geological and topographical factors using principal component analysis (PCA) and analysis of variance (ANOVA) was developed and tested. ANOVA results demonstrated statistically significant differences in classed variable groups as well as in classed GRP scores with regard to hydrogeological indicators, such as specific capacity (SC) and transmissivity. Results of PCA were used to govern the weight of the variables used in the prediction maps. GRP scores were able to identify 79% of wells in a verification dataset, which had SC values less than the total dataset median. GRP values showed statistically significant correlations using both parametric (using transformed datasets) and non‐parametric methods. The method shows promise for municipal or regional level planning in crystalline terrains with high levels of heterogeneity and anisotropy as a hydrogeologically and statistically based tool to assist in assessing groundwater resources. The methodology is executed in a geographic information systems environment, and uses often readily available data, such as geological maps, feature maps and topography, and thus does not require expensive and time‐consuming aquifer tests.     

识别地下水超采区井水位异常性质的理论与方法

长期过量开采地下水,使地下水位持续下降、水质发生变化,动水位观测井断流;地面沉降造成井管上窜,观测管路系统被损坏等,这些现象对地震地下流体观测地震前兆异常的正确判断带来很大困难。应用水文地质理论与方法,分析含水层的水均衡状态、应力-应变状态及其与水位动态的关系,初步探讨了超采区井水位异常性质的理论与方法。结果表明,根据井孔所在区水位下降漏斗的扩散特征,结合以上所提到的理论和方法,依据资料多年变化特征,可以较准确地判断异常的性质。研究结果有助于区分单一集中抽水与长期地下水超采对水位观测的影响,有助于正确识别超采区水位前兆异常,有助于地震分析预报水平的提高     

Estimating Aquifer Properties Using Derivative Analysis of Water Level Time Series from Active Well Fields

Water level time series from groundwater production wells offer a transient dataset that can be used to estimate aquifer properties in areas with active groundwater development. This article describes a new parameter estimation method to infer aquifer properties from such datasets. Specifically, the method analyzes long‐term water level measurements from multiple, interacting groundwater production wells and relies on temporal water level derivatives to estimate the aquifer transmissivity and storativity. Analytically modeled derivatives are compared to derivatives calculated directly from the observed water level data; an optimization technique is used to identify best‐fitting transmissivity and storativity values that minimize the difference between modeled and observed derivatives. We demonstrate how the consideration of derivative (slope) behavior eliminates uncertainty associated with static water levels and well‐loss coefficients, enabling effective use of water level data from groundwater production wells. The method is applied to time‐series data collected over a period of 6 years from a municipal well field operating in the Denver Basin, Colorado (USA). The estimated aquifer properties are shown to be consistent with previously published values. The parameter estimation method is further tested using synthetic water level time series generated with a numerical model that incorporates the style of heterogeneity that occurs in the Denver Basin sandstone aquifers.     

Reverse Water‐Level Fluctuations Associated with Fracture Connectivity

Reverse water‐level fluctuations (RWFs), a phenomenon in which water levels rise briefly in response to pumping, were detected in monitoring wells in a fractured siliciclastic aquifer system near a deep public supply well. The magnitude and timing of RWFs provide important information that can help interpret aquifer hydraulics near pumping wells. A RWF in a well is normally attributed to poroelastic coupling between the solid and fluid components in an aquifer system. In addition to revealing classical pumping‐induced poroelastic RWFs, data from pressure transducers located at varying depths and distances from the public supply well suggest that the RWFs propagate rapidly through fractures to influence wells hundreds of meters from the pumping well. The rate and cycling frequency of pumping is an important factor in the magnitude of RWFs. The pattern of RWF propagation can be used to better define fracture connectivity in an aquifer system. Rapid, cyclic head changes due to RWFs may also serve as a mechanism for contaminant transport.     

Analytical and Semi‐Analytical Tools for the Design of Oscillatory Pumping Tests

Oscillatory pumping tests—in which flow is varied in a periodic fashion—provide a method for understanding aquifer heterogeneity that is complementary to strategies such as slug testing and constant‐rate pumping tests. During oscillatory testing, pressure data collected at non‐pumping wells can be processed to extract metrics, such as signal amplitude and phase lag, from a time series. These metrics are robust against common sensor problems (including drift and noise) and have been shown to provide information about aquifer heterogeneity. Field implementations of oscillatory pumping tests for characterization, however, are not common and thus there are few guidelines for their design and implementation. Here, we use available analytical solutions from the literature to develop design guidelines for oscillatory pumping tests, while considering practical field constraints. We present two key analytical results for design and analysis of oscillatory pumping tests. First, we provide methods for choosing testing frequencies and flow rates which maximize the signal amplitude that can be expected at a distance from an oscillating pumping well, given design constraints such as maximum/minimum oscillator frequency and maximum volume cycled. Preliminary data from field testing helps to validate the methodology. Second, we develop a semi‐analytical method for computing the sensitivity of oscillatory signals to spatially distributed aquifer flow parameters. This method can be quickly applied to understand the “sensed” extent of an aquifer at a given testing frequency. Both results can be applied given only bulk aquifer parameter estimates, and can help to optimize design of oscillatory pumping test campaigns.     

A Model to Simulate Pumping Tests in Wells with High Loss and Backflow

Traditional methods of analyzing pumping tests in single wells fail when the well loss is very high due to a low transmissivity skin. Because of the restricted rate at which water can enter a high loss well from the aquifer, well casing storage becomes a significant factor. Additionally, if a slug of water enters the well from the pump column immediately after the pump is switched off, it has a long‐lasting significant effect on the recovering water level in the well because it cannot be absorbed rapidly by the aquifer. A theoretical model is derived here that simulates the water level in a well in these circumstances. In the model, the continuously changing rate of water inflow from the aquifer to the well is approximated by a step function with a finite difference time step. It is demonstrated by a real example that the model can be applied easily to analyze pumping tests, including tests with a varying pumping rate. The analysis confirms suspected high well loss, calculates the unknown rate of backflow, and determines the aquifer's transmissivity.     

Lithological control of land subsidence induced by groundwater withdrawal in new urban AREAS (Granada Basin,SE Spain). Multiband DInSAR monitoring

Ground subsidence in the southeastern border of the Granada Basin (SE Spain) has been studied using remote sensing techniques. Over the last decades, the region has experienced a huge urban expansion, which has caused a substantial increase in water supply requirements. Water needs are exclusively met by groundwater by means of numerous pumping wells, which exploit a confined detrital aquifer of alluvial fan deposits with a heterogeneous facies distribution. A general piezometric level decline (up to 50 m) has been recorded in the aquifer during the past 30 years that has induced the generation of a subsiding area with oval shape oriented WNW‐ESE just where the new urban areas and pumping wells are located. Subsidence has been monitored by exploiting synthetic aperture radar (SAR) images from ENVISAT (2003–2009) and Cosmo‐SkyMed (2011–2014). A new approach, which combines A‐DInSAR and small‐area persistent scatterer interferometry (PSI) analysis, has been applied obtaining a good accuracy regarding temporal and spatial dimension of the subsidence. ENVISAT data (2003–2009) reveal subsidence rates up to 10–15 mm/year, and Cosmo‐SkyMed (2011–2014) values slightly lower; up to 10 mm/year. Temporal variations in the subsidence velocity are in accordance with the rainfall pattern and piezometric fluctuations in the aquifer. The sector with highest rates of subsidence does not correspond to the area with more intense groundwater exploitation but to the area with greater presence of clays in the confining layer of the aquifer. There is a clear lithological control in the spatial distribution of the ground subsidence. This work integrates detailed geological and hydrogeological data with differential SAR interferometry monitoring with the aim to better understand subsidence processes in detrital aquifers with small‐scale heterogeneity. Copyright © 2016 John Wiley & Sons, Ltd.     

Interpretation of Water Level Changes in the High Plains Aquifer in Western Kansas

Water level changes in wells provide a direct measure of the impact of groundwater development at a scale of relevance for management activities. Important information about aquifer dynamics and an aquifer's future is thus often embedded in hydrographs from continuously monitored wells. Interpretation of those hydrographs using methods developed for pumping‐test analyses can provide insights that are difficult to obtain via other means. These insights are demonstrated at two sites in the High Plains aquifer in western Kansas. One site has thin unconfined and confined intervals separated by a thick aquitard. Pumping‐induced responses in the unconfined interval indicate a closed (surrounded by units of relatively low permeability) system that is vulnerable to rapid depletion with continued development. Responses in the confined interval indicate that withdrawals are largely supported by leakage. Given the potential for rapid depletion of the unconfined interval, the probable source of that leakage, it is likely that large‐scale irrigation withdrawals will not be sustainable in the confined interval beyond a decade. A second site has a relatively thick unconfined aquifer with responses that again indicate a closed system. However, unlike the first site, previously unrecognized vertical inflow can be discerned in data from the recovery periods. In years of relatively low withdrawals, this inflow can produce year‐on‐year increases in water levels, an unexpected occurrence in western Kansas. The prevalence of bounded‐aquifer responses at both sites has important ramifications for modeling studies; transmissivity values from pumping tests, for example, must be used cautiously in regional models of such systems.     

Hydrogeological Studies to Identify the Trend of Concealed Section of the North Tabriz Fault (Iran)

The North Tabriz Fault (NTF) is the predominant regional‐scale tectonic structure in the northwest of Iran. In the east side of the city of Tabriz, a portion of the fault trend has been completely concealed by recent sediments and urbanization. In this paper, some hydrogeological methods are used to locate the concealed sector. As is clear from the pumping tests results, despite the fact that the northern observation wells were closer to the pumping wells than the southern ones, they have not been affected by pumping. Conversely, all southern wells were affected by pumping and displayed decline of the water table. In addition, obvious differences in groundwater levels combined with clear differences in groundwater quality within a short distance across the probable fault trend are sufficient reasons for the presence of the fault that behaves as a barrier to groundwater lateral flows. Significant change in the elevation of the bedrock base of the aquifer over less than 200 m suggests that the fault has near vertical dip. These results indicate that the inferred trend of the NTF closely conforms to its actual trend. Therefore, the hydrogeological studies can be complementary tools to determine the position and trend of concealed faults.     

Analysis of steady ground water flow toward wells in a confined-unconfined aquifer

A confined aquifer may become unconfined near the pumping wells when the water level falls below the confining unit in the case where the pumping rate is great and the excess hydraulic head over the top of the aquifer is small. Girinskii's potential function is applied to analyze the steady ground water flow induced by pumping wells with a constant-head boundary in a mixed confined-unconfined aquifer. The solution of the single-well problem is derived, and the critical radial distance at which the flow changes from confined to unconfined condition is obtained. Using image wells and the superposition method, an analytic solution is presented to study steady ground water flow induced by a group of pumping wells in an aquifer bounded by a river with constant head. A dimensionless function is introduced to determine whether a water table condition exists or not near the pumping wells. An example with three pumping wells is used to demonstrate the patterns of potentiometric surface and development of water table around the wells.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Identifying Aquifer Type in Fractured Rock Aquifers using Harmonic Analysis

Determining aquifer type, unconfined, semi‐confined, or confined, by drilling or performing pumping tests has inherent problems (i.e., cost and complex field issues) while sometimes yielding inconclusive results. An improved method to cost‐effectively determine aquifer type would be beneficial for hydraulic mapping of complex aquifer systems like fractured rock aquifers. Earth tides are known to influence water levels in wells penetrating confined aquifers or unconfined thick, low‐porosity aquifers. Water‐level fluctuations in wells tapping confined and unconfined aquifers are also influenced by changes in barometric pressure. Harmonic analyses of water‐level fluctuations of a thick (～1000 m) carbonate aquifer located in south‐central Oklahoma (Arbuckle‐Simpson aquifer) were utilized in nine wells to identify aquifer type by evaluating the influence of earth tides and barometric‐pressure variations using signal identification. On the basis of the results, portions of the aquifer responded hydraulically as each type of aquifer even though there was no significant variation in lithostratigraphy. The aquifer type was depth dependent with confined conditions becoming more prevalent with depth. The results demonstrate that harmonic analysis is an accurate and low‐cost method to determine aquifer type.     

Detecting Drawdowns Masked by Environmental Stresses with Water‐Level Models

Detecting and quantifying small drawdown at observation wells distant from the pumping well greatly expands the characterized aquifer volume. However, this detection is often obscured by water level fluctuations such as barometric and tidal effects. A reliable analytical approach for distinguishing drawdown from nonpumping water‐level fluctuations is presented and tested here. Drawdown is distinguished by analytically simulating all pumping and nonpumping water‐level stresses simultaneously during the period of record. Pumping signals are generated with Theis models, where the pumping schedule is translated into water‐level change with the Theis solution. This approach closely matched drawdowns simulated with a complex three‐dimensional, hypothetical model and reasonably estimated drawdowns from an aquifer test conducted in a complex hydrogeologic system. Pumping‐induced changes generated with a numerical model and analytical Theis model agreed (RMS as low as 0.007 m) in cases where pumping signals traveled more than 1 km across confining units and fault structures. Maximum drawdowns of about 0.05 m were analytically estimated from field investigations where environmental fluctuations approached 0.2 m during the analysis period.     

Simulation of Flow in a Complex Aquifer System Subjected to Long-Term Well Network Growth

In west-central Lower Peninsula of Michigan, population growth and expanded agricultural activities over recent decades have resulted in significant increases in distributed groundwater withdrawals. The growth of the extensive well network and anecdotes of water shortages (dry wells) have raised concerns over the region's groundwater sustainability. We developed an unsteady, three-dimensional (3D) groundwater flow model to describe system dynamics over the last 50 years and evaluate long-term impacts of groundwater use. Simulating this large aquifer system was challenging; the site is characterized by strong, spatially distributed, and statistically nonstationary heterogeneity, making it difficult to avoid over-parameterization using traditional approaches for conceptualizing and calibrating a flow model. Moreover, traditional pumping and water level data were lacking and prohibitively expensive to collect given the large-scale and long-term nature of this study. An integrated, stochastic-deterministic approach was developed to characterize the system and calibrate the flow model through innovative use of high-density water well datasets. This approached allowed (1) implementation of a “zone-based,” nonstationary stochastic approach to conceptualize complex spatial variability using a small set of geologic material types; (2) modeling the spatiotemporal evolution of many water well withdrawals across several decades using sector-based parameterization; and (3) critical analysis of long-term water level changes at different locations in the aquifer system for characterizing the system dynamics and calibrating the model. Results show the approach is reasonably successful in calibrating a complex model for a highly complex site in a way that honors complex distributed heterogeneity and stress configurations.     

Poroelastic Effects on Fracture Characterization

Reverse water‐level fluctuations have been widely observed in aquitards or aquifers separated from a pumped confined aquifer (Noordbergum effect) immediately after the initiation of pumping. This same reverse fluctuation has been observed in a fractured crystalline‐rock aquifer at the Coles Hill uranium site in Virginia in which the reverse water‐level response occurs within a pumped fracture and results from an instantaneous strain response to pumping that precedes the pore‐pressure response in observation wells of sufficient distance from the pumped well. This response is referred to as the Mandel‐Cryer effect. The unique aspect of this water level rise during a controlled 24 h pumping test was that the reverse water levels lasted for approximately 100 min and reached a magnitude of nearly 1 cm prior to a typical drawdown response. The duration and magnitude of the response reflects the poromechanical properties of the fractured host rock and hydraulic properties of the pumped fracture. An axisymmetric flow and deformation model were developed using Biot2 in an effort to simulate the observed water‐level response along an assumed 0.5 to 1.0 cm aperture horizontal fracture 176 m from the pumping well and to identify the importance of the poroelastic effect. Results indicate that traditional aquifer‐testing methods that ignore the poromechanical response are not significantly different than results that include the response. However, the poroelastic effect allows for more accurate and efficient parameter calibration.     

Empirical relationships for estimating stream depletion by a well pumping near a gaining stream

Siting wells near streams requires an accurate estimate of the quantity of water derived from the river due to pumping. A number of hydrogeological and hydraulic parameters influence this value. This study estimates stream depletion under steady-state conditions for a variety of hydrogeological systems. A finite differences model was used to analyze several hydrogeological situations, and for each of these the stream depletion was estimated using an advective transport method. An empirical equation for stream depletion was obtained for the case of a stream that partially penetrates the aquifer and a pumping well that is screened over a portion of the aquifer. The derived equation, which is valid for both isotropic and anisotropic conditions, expresses stream depletion as a function of the unit inflow to the river, the discharge of the pumping well, the well screen length, the distance between the river and pumping well, the wetted perimeter, and a new parameter called "overlap," which is defined to be the distance between the riverbed and the top of well screen. The overlap parameter makes it possible to consider indirectly the vertical component of flow, which is accentuated when the well is screened below the streambed. The formula proposed here should be useful in deciding where to locate a pumping well and to decide the appropriate length of its screen.     

Fecal Indicator Bacteria Variability in Samples Pumped from Monitoring Wells

The detection of microbiological contamination in drinking water from groundwater wells is often made with a limited number of samples that are collected using traditional geochemical sampling protocols. The objective of this study is to examine the variability of fecal indicator bacteria, as observed using discrete samples, due to pumping. Two wells were instrumented as multilevel piezometers in a bedrock aquifer, and bacterial enumeration was conducted on a total of 166 samples (for total coliform, fecal coliform, Escherichia coli, and fecal streptococci) using standard membrane filtration methods. Five tests were conducted using pumping rates ranging from 0.3 to 17 L/min in a variety of purging scenarios, which included constant and variable (incremental increase and decrease) flow. The results clearly show a rapid and reproducible, 1 to 2 log‐unit decrease in fecal indicator bacteria at the onset of pumping to stabilized, low‐level concentrations prior to the removal of three to five well volumes. The pumping rate was not found to be correlated with the magnitude of observed bacterial counts. Based on the results, we suggest sampling protocols for fecal indicator bacteria that include multiple collections during the course of pumping, including early‐time samples, and consider other techniques such as microscopic enumeration when assessing the source of bacteria from the well‐aquifer system.     

Transient analysis of the source of water to wells: Cape Cod, Massachusetts

A transient flow modeling analysis for potential public-supply wells on western Cape Cod, Massachusetts, demonstrates the difference between transient and steady-state recharge areas can have important implications for wellhead protection. An example of a single pumping well illustrates that commonly used steady-state time-related capture areas do not represent the recharge area and travel times of water being pumped from the well until sufficient time has elapsed for steady-state flow conditions to be established. Until that time, transient recharge areas are needed to account for the portion of water discharging from the well that entered the aquifer before pumping started. An example of two pumping wells demonstrates the same area at the water table cannot supply water to more than one well under steady-state flow conditions. Transient recharge areas to multiple wells can overlap, however, until steady-state flow conditions are established. The same area can, therefore, be a source of water to more than one well during early pumping times, and the water pumped from a given well may derive from source areas, including contaminated areas, that do not lie within the well's steady-state recharge area.     

Aquifer modelling of the Ganga–Mahawa sub‐basin,a part of the Central Ganga Plain,Uttar Pradesh,India

The Ganga–Mahawa sub‐basin, which has an area of 1280 km² forms the western part of the Central Ganga Plain in the Moradabad and Badaun districts of western Uttar Pradesh, India. The Bundelkhand granite forms the basement complex, overlain unconformably by the upper Vindhyan sequence, which is further overlain by the Neogene (Middle and Upper) Siwaliks and finally by Quaternary alluvium. Four geomorphological units, the Varanasi older alluvial plain, Aligarh older alluvial plain, terrace zones and the Ganga recent floodplain, abandoned channels, channel scars and meander scars represent various landforms. The hydrogeological cross‐sections indicate the occurrence of a single aquifer down to 120 m. Some influent seepage from the River Ganga could be seen around Gangeswari, but the rest of the River Ganga is effluent. Groundwater‐flow modelling was carried out to assess the degree of Ganga river and aquifer interaction. The River Ganga marks the western boundary; boundaries to the northeast and southeast are set as fixed heads to simulate lateral inflow into and outflow from the sub‐basin respectively. The eastern boundary is simulated as a no‐flow condition. The Mahawa and Badmar rivers are considered to be effluent. The area modelled is covered by a grid of 34 rows×46 columns with three layers, viz., an unconfined aquifer, an aquitard which is underlain by a semi‐confined to confined aquifer. The permeability distribution was inferred from morphometric analysis and pumping tests. Natural recharge due to monsoon rainfall forms the main input. The River Ganga stage data at Ahar, Naora and Ramghat has been used for assigning surface water levels and river bed elevations in the model. Abstraction from all existing deep and shallow tube wells has been assigned as output at various cells. A steady state flow simulation was carried out and calibrated against the June 1986 water level; subsequent transient conditions were calibrated up to May 1995. The computed groundwater balance was comparable to that estimated from field investigations. The aquifer modelling study has attempted to integrate all available information and provided a tool that could be used for predictive simulation. Copyright © 2000 John Wiley & Sons, Ltd.