VARIABLE DISCHARGE TYPE CURVE SOLUTIONS FOR CONFINED AQUIFERS1

ABSTRACT: Almost all of the existing solutions of aquifer test analysis assume constant discharge rate boundary conditions. In many actual pumping tests, constant discharge cannot be maintained. This paper presents a general solution for an aquifer test with variable discharge. The exact solution of drawdown distribution around an infinitesimally small diameter well in a uniform, horizontal, extensive, homogeneous, and isotropic confined aquifer is presented when the discharge changes with time during the aquifer test period. A general equation for the type curves resulting from any discharge variability is given, and its application for the exponential changes is presented in detail. A simple straight line procedure is proposed for field applications by considering late time drawdown data. The consideration of constant discharge leads to overestimation of transmissivity but storativity is underestimated. Practical application of the discharge variability is illustrated by a field example.     

SELECTED METHODS OF AQUIFER TEST ANALYSIS1

Changes in groundwater levels due to a well discharging at a constant rate are used with various formulas to determine hydraulic properties of aquifers and their confining bed and to detect the presence of aquifer boundaries. These formulas are generally solved by graphical methods.     

ESTIMATING AQUIFER PARAMETERS FROM A HORIZONTAL WELL PUMPING TEST IN AN UNCONFINED AQUIFER1

ABSTRACT: A pumping test on a phreatic glacial till aquifer was performed near Ames, Iowa, in November 1990. The head in a horizontal well was pumped down rapidly and then held constant for the duration of the 30-hour test. Throughout the test, the flow rate at the pumped well and the head at an adjacent vertical observation well were recorded. The pumping test data were used to determine the hydraulic conductivity and specific yield of the aquifer. The results indicate a hydraulic conductivity of 2.23 × 10^?5 cm/s and a representative specific yield of 0.03. Hydraulic conductivity was calculated by a simple integration of Darcy's Law after extrapolating the data to steady state. Specific yield was determined by use of several methods from the literature and a new method proposed by the author. The results show that specific yield increases with time, and that each method is within an order of magnitude of the others.     

NON-STEADY SPHERICAL FLOW TO A CAVITY WELL IN AN INFINITE NON-LEAKY AQUIFER1

ABSTRACT: The non-steady drawdown distribution near a cavity well discharging from an infinite non-leaky artesian aquifer is presented. The variation of drawdown with time and distance caused by a cavity well of constant discharge in a confined aquifer of uniform thickness and uniform permeability is obtained. The solution is expressed in a series form which converges rapidly so that only two terms of the series are needed to obtain an accuracy of more than 95 percent. A simplified approach has been suggested to find the aquifer characteristics.     

AQUIFER PARAMETERS FROM CONSTANT DRAWDOWN NONSTEADY-LEAKY TYPE CURVES1

ABSTRACT: A drain function and set of type curves were defined for the mathematical solution that represents one-dimensional flow under nonsteady conditions in a leaky aquifer for the constant drawdown boundary condition. A match point procedure was developed for determining the aquifer parameters transmissivity, storage coefficient, and leakance based on the drain function and type curves. Use of the procedure is illustrated by an example that utilizes simulated aquifer drawdowns and flowrate data. The drain function and type curves developed in this investigation include the effects of leakage for the constant drawdown boundary condition, which is not included in the existing drain function and type curve found in the literature. Thus, a new set of type curves was developed that can be used to analyze drawdowns for one-dimensional flow in a leaky aquifer with constant drawdown at a line sink. Applications would include flow to a canal or river, drainage of agricultural lands, and dewatering associated with strip mining operations.     

CONTROL OF NATURAL BRINE SPRINGS IN BRAZOS RIVER BASIN PART II: BRINE DISPOSAL1

ABSTRACT: Water quality in the Brazos River of Texas is seriously degraded by natural salt pollution. Two thousand tons/day of total dissolved solids emanate from brine springs and seeps in the Upper Brazos River drainage. Approximately 45 percent of the total salt load comes from a relatively small flow in the Dove Creek area. The companion paper demonstrates that a system of wells pumping brine at a constant rate of about 2 cfs from the near surface aquifer should eliminate the brine springs in this area. In this paper, injection into deep brine aquifers is shown to be a feasible brine disposal alternative. Four brine aquifers were determined from the literature to be possible injection zones. Accurate net aquifer thickness maps were generated in a 23 by 14 mile area centered on the Dove Creek area for three of the aquifers from an interpretation of 41 well logs. Constant injection for a project life of 100 years was simulated using the SWIFT/486 software. Modeling suggests that one well would be sufficient to inject the entire disposal volume into either the Strawn or Ellenburger Formation.     

NONSTEADY FLOW TO A WELL WITH TIME DEPENDENT DRAWDOWN1

Confined flow toward a single well of finite radius in an extensive aquifer of uniform transmissibility is studied under the assumption of time-dependent drawdown. Three particular cases are considered: (a) linear drawdown (including constant drawdown); (b) exponential drawdown; (c) periodic (sinusoidal) drawdown. The differential equation governing unsteady axial symmetric flow toward a single well in a confined aquifer is solved for the three different situations by the use of the Laplace transform method. The resulting expressions are integrated by adapting a modified Gemant scheme. General computer programs have been developed and operated for several combinations of characteristics. The results are plotted to show the effect of time dependent drawdown on the variation of the well discharge and the piezometric head distribution.     

MODELING CONCENTRATION VARIATIONS IN HIGH-CAPACITY WELLS: IMPLICATIONS FOR GROUNDWATER SAMPLING1

ABSTRACT: High-capacity wells are used as a convenient and economical means of sampling groundwater quality. Although the inherent limitations of using these wells are generally recognized, little has been done to investigate how these wells actually sample groundwater. A semi-analytical particle tracking model is used to illustrate the influence of variable vertical contaminant distributions and aquifer heterogeneity on the composition of water samples from these wells during short pumping periods. The hypothetical pumping well used in the simulations is located in an unconfined, alluvial aquifer with a shallow water table and concentration gradients of nitrate-nitrogen contamination. This is a typical setting for many irrigated areas in the United States. The main conclusions are: (1) high-capacity wells underestimate the average amount of contamination within an aquifer; (2) shapes of concentration-time curves for high-capacity wells appear to be governed by the distribution of the contaminant and travel times to the well; (3) variables such as well construction, pumping rate, and hydrogeologic properties contribute to the magnitude of the concentration-time curves at individual high-capacity wells; and (4) a sampling strategy using concentration-time curves based on the behavioral characteristics of the well rather than individual samples will provide a much better framework for interpreting spatial contaminant distributions.     

Assessing the economic feasibility of regional deep saline aquifer CO2 injection and storage: A geomechanics-based workflow applied to the Rose Run sandstone in Eastern Ohio, USA

Typical top-down regional assessments of CO₂ storage feasibility are sufficient for determining the maximum volumetric capacity of deep saline aquifers. However, they do not reflect the regional economic feasibility of storage. This is controlled, in part, by the number and type of injection wells that are necessary to achieve regional CO₂ storage goals. In contrast, the geomechanics-based assessment workflow that we present in this paper follows a bottom-up approach for evaluating regional deep saline aquifer CO₂ storage feasibility. The CO₂ storage capacity of an aquifer is a function of its porous volume as well as its CO₂ injectivity. For a saline aquifer to be considered feasible in this assessment it must be able to store a specified amount of CO₂ at a reasonable cost per ton of CO₂. The proposed assessment workflow has seven steps that include (1) defining the storage project and goals, (2) characterizing the geology and developing a geomechanical model of the aquifer, (3) constructing 3D aquifer models, (4) simulating CO₂ injection, (5,6) evaluating CO₂ injection and storage feasibility (with and without injection well stimulation), and (7) determining whether it is economically feasible to proceed with the storage project. The workflow was applied to a case study of the Rose Run sandstone aquifer in the Eastern Ohio River Valley, USA. We found that it is feasible in this region to inject 113 Mt CO₂/year for 30 years at an associated well cost of less than US $1.31/t CO₂, but only if injectivity enhancement techniques such as hydraulic fracturing and injection induced micro-seismicity are implemented.     

GROUND-WATER TEMPERATURE FLUCTUATIONS AT LYONS FERRY FISH HATCHERY,WASHINGTON1

ABSTRACT: The well field serving the Lyons Ferry Fish Hatchery has experienced reduced water temperatures following continued, periodic withdrawal of large volumes of water. In January 1985, the well field temperature was 49°F, which is less than the optimal 52°F for raising salmon and steelhead trout. The aquifer supplying the hatchery is in hydraulic and thermal connection with the Snake River and a flooded embayment of the Palouse River. Ground-water temperatures in the well field cycle on an annual basis in response to changes in surface water temperature and pumping rate. Numerical simulation of the well field, using a simplified mixing cell model, demonstrates the coupling of well field hydraulics and aquifer thermal response. Alternative pumping schedules indicate that it is feasible to adjust ground-water pumping to effectively store heat in the aquifer during the summer months when surface water temperatures are elevated. Sensitivity analysis of this model indicated that the primary controls of the system's thermal response are the volume of the aquifer assumed to contribute to the well field and temperature of the overlying surface water body.     

A Stochastic Modeling Approach to Address Hydrogeologic Uncertainties in Modeling Wellhead Protection Boundaries in Karst Aquifers1

Abstract: Development of any numerical ground‐water model is dependent on hydrogeologic data describing the subsurface. These data are obtained from geologic core analyses, stratigraphic analyses, aquifer performance tests, and geophysical studies. But typically in remote areas, these types of data are very sparse and site‐specific in terms of the aerial extent of the resource to be modeled. Uncertainties exist as to how well the available data from a few locations defines a heterogeneous surficial aquifer such as the Biscayne Aquifer in Miami‐Dade County, Florida. This is particularly the case when an exceptionally conductive horizontal flow zone is detected at one site due to specialized testing that was not historically conducted at the other at sites that provided data for the model. Not adequately accounting for the potential effect of the high flow zone in the aquifer within a ground‐water numerical model, even though the zone may be of very limited thickness, might underpredict the well field protection capture boundaries. Applied Stochastic ground‐water modeling in determining well field protection zones is steadily becoming important in addressing the uncertainty of the hydrogeologic subsurface parameters, specifically in karstic heterogeneous aquifers. This is particularly important in addressing the uncertainty of a 60‐day travel time capture zone in the Northwest Well Field, Miami‐Dade County, where a predominantly high flow zone controls much of the flow in the production wells. A stochastic ground‐water modeling application along with combination of pilot points and regularization technique is presented to further consolidate the uncertainty of the subsurface.     

Modeling the Potential Impact of Seasonal and Inactive Multi‐Aquifer Wells on Contaminant Movement to Public Water‐Supply Wells1

Johnson, R.L., B.R. Clark, M.K. Landon, L.J. Kauffman, and S.M. Eberts, 2011. Modeling the Potential Impact of Seasonal and Inactive Multi‐Aquifer Wells on Contaminant Movement to Public Water‐Supply Wells. Journal of the American Water Resources Association (JAWRA) 47(3):588‐596. DOI: 10.1111/j.1752‐1688.2011.00526.x Abstract: Wells screened across multiple aquifers can provide pathways for the movement of surprisingly large volumes of groundwater to confined aquifers used for public water supply (PWS). Using a simple numerical model, we examine the impact of several pumping scenarios on leakage from an unconfined aquifer to a confined aquifer and conclude that a single inactive multi‐aquifer well can contribute nearly 10% of total PWS well flow over a wide range of pumping rates. This leakage can occur even when the multi‐aquifer well is more than a kilometer from the PWS well. The contribution from multi‐aquifer wells may be greater under conditions where seasonal pumping (e.g., irrigation) creates large, widespread downward hydraulic gradients between aquifers. Under those conditions, water can continue to leak down a multi‐aquifer well from an unconfined aquifer to a confined aquifer even when those multi‐aquifer wells are actively pumped. An important implication is that, if an unconfined aquifer is contaminated, multi‐aquifer wells can increase the vulnerability of a confined‐aquifer PWS well.     

UTILIZATION OF A GEOGRAPHIC INFORMATION SYSTEM TO IDENTIFY THE PRIMARY AQUIFER PROVIDING GROUND WATER TO INDWIDUAL WELLS IN EASTERN ARKANSAS1

ABSTRACT: A Geographic Information System (GIS) was used to develop an automated procedure for identifying the primary aquifers supplying ground water to individual wells in eastern Arkansas. As mandated by state law, water-use data are reported by ground-water withdrawers annually to the Arkansas Soil and Water Conservation Commission, and stored in the Arkansas Site-Specific Water-Use Data System provided and supported by the U.S. Geological Survey. Although most withdrawers are able to provide the amount of water withdrawn and the depth of their wells, very few are able to provide the name of the aquifer from which they withdraw water. GIS software was used to develop an automated procedure for identifying the primary aquifers supplying ground water to individual wells in eastern Arkansas. The software was used to generate a spatial representation of the bottom boundary for the Mississippi River Valley alluvial aquifer (the shallowest aquifer) in eastern Arkansas from well log-data collected by the U.S. Geological Survey. The software was then used to determine the depth of the aquifer bottom at reported well locations to ascertain whether the Mississippi River Valley alluvial aquifer or a deeper aquifer was the primary aquifer providing water to each well. The alluvial aquifer was identified as the primary aquifer for about 23,500 wells.     

ANALYSIS OF PUMPING TESTS: SIGNIFICANCE OF WELL DIAMETER,PARTIAL PENETRATION,AND NOISE1

ABSTRACT: The nonlinear least squares (NLS) method was applied to pumping and recovery aquifer test data in confined and unconfined aquifers with finite diameter and partially penetrating pumping wells, and with partially penetrating piezometers or observation wells. It was demonstrated that noiseless and moderately noisy drawdown data from observation points located less than two saturated thicknesses of the aquifer from the pumping well produced an exact or acceptable set of parameters when the diameter of the pumping well was included in the analysis. The accuracy of the estimated parameters, particularly that of specific storage, decreased with increases in the noise level in the observed drawdown data. With consideration of the well radii, the noiseless drawdown data from the pumping well in an unconfined aquifer produced good estimates of horizontal and vertical hydraulic conductivities and specific yield, but the estimated specific storage was unacceptable. When noisy data from the pumping well were used, an acceptable set of parameters was not obtained. Further experiments with noisy drawdown data in an unconfined aquifer revealed that when the well diameter was included in the analysis, hydraulic conductivity, specific yield and vertical hydraulic conductivity may be estimated rather effectively from piezometers located over a range of distances from the pumping well. Estimation of specific storage became less reliable for piezometers located at distances greater than the initial saturated thickness of the aquifer. Application of the NLS to field pumping and recovery data from a confined aquifer showed that the estimated parameters from the two tests were in good agreement only when the well diameter was included in the analysis. Without consideration of well radii, the estimated values of hydraulic conductivity from the pumping and recovery tests were off by a factor of four.     

Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge1

Boggs, Kevin G., Robert W. Van Kirk, Gary S. Johnson, Jerry P. Fairley, and P. Steve Porter, 2010. Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge. Journal of the American Water Resources Association (JAWRA) 46(6):1116–1132. DOI: 10.1111/j.1752-1688.2010.00479.x Abstract: There is a need to develop a general understanding of how variations in aquifer recharge are reflected in discharge. Analytical solutions to the linearized Boussinesq equation governing flow in an unconfined aquifer provide a unified mathematical framework to quantify relationships among lag time, attenuation and distance between aquifer recharge and discharge and the effect of an up-gradient no-flow boundary. We applied this framework to three types of recharge: (1) instantaneous, (2) periodic, and (3) constant rate for a finite duration. When the temporal scale of recharge exceeds the diffusive aquifer time scale, recharge will be reflected in discharge quickly and with little attenuation. When aquifer time scale is large, most recharge events are shorter in scale than that of the aquifer, resulting in large attenuation. Attenuation is more sensitive to boundary effects than lag time, and boundary effects increase as recharge time scale increases. Boundary effects can often be ignored when the recharge source is farther than 1/3 of the domain length away from the no-flow boundary. We illustrate analytical results with application to the economically critical Eastern Snake River Plain Aquifer in Idaho. In this aquifer, detectable annual and decadal cycles in discharge can result from recharge no farther than 20 and 60 km away from the discharge point, respectively. The effects of more distant, long-term recharge can be detected only after a time lag of several decades.     

STREAMFLOW DEPLETION: MODELING OF REDUCED BASEFLOW ANI INDUCED STREAM INFILTRATION FROM SEASONALLY PUMPED WELLS1

ABSTRACT: Numerical modeling techniques are used to analyze streamflow depletion for stream‐aquifer systems with baseflow. The analyses calculated two flow components generated by a pumping well located at a given distance from a river that is hydraulically connected to an unconfined aquifer. The two components are induced stream infiltration and reduced baseflow; both contribute to total streamflow depletion. Simulation results suggest that the induced infiltration, the volume of water discharged from the stream to the aquifer, has a shorter term impact on streamflow, while the reduced baseflow curves show a longer term effect. The peak impacts of the two hydrologic processes on streamflow occur separately. The separate analysis helps in understanding the hydrologic interactions between stream and aquifer. Practically, it provides useful information about contaminant transport from stream to aquifer when water quality is a concern, and for areas where water quantity is an issue, the separate analysis offers additional information to the development of water resource management plan.     

SENSITIVITY OF GROUNDWATER FLOW MODELS TO VERTICAL VARIABILITY OF AQUIFER CONSTANTS1

ABSTRACT: The Ogallala aquifer in the Oklahoma Panhandle is in need of better management because of increased groundwater demand which has caused declines in static water levels at an alarming rate. A groundwater management computer model was developed for the Ogallala aquifer in the Texas Panhandle and treats the aquifer as a homogeneous system. In this study, the computer model has been modified in order to evaluate the effects of vertical layering on semi-static water level changes which occur during the dewatering of a single unconfined aquifer. The modified model was applied to a study area near Guymon, Oklahoma, using both the homogeneous and the multilayered cases. The aquifer is characterized by a saturated thickness of 400 feet. The accumulated drawdown values of the homogeneous and the multilayered cases demonstrate that an average difference of approximately 22% of the original saturated thickness occurs between the two cases before the base of the aquifer is encountered. Approximately 25% more time is required to dewater the layered aquifer. Thus, vertical variations of lithology in an aquifer such as the Ogallala should be considered when prediction is made relative to groundwater management.     

NUMERICAL ANALYSIS OF PUMPING FROM CONFINED-UNCONFINED AQUIFERS1

ABSTRACT: A numerical method is presented for the analysis of a pumped well in a homogeneous aquifer with allowance made for the decrease in saturated depth, vertical components of flow, the possibility of regions of the aquifer changing between the confined and unconfined states and the effect of different outer boundaries. The method is based on a discrete space, backward difference time, approximation. A particular example considered in detail concerns heavy pumping from one of a regular array of wells in an unconfined aquifer until the drawdown in the well reaches a critical value. Non-dimensional curves are presented relating the time and volume dewatered to the quantity discharged from the well. A further example investigates the effect of an initial confining pressure on the aquifer behaviour.     

THE ECONOMIC VALUE OF GROUND WATER RECHARGE FOR IRRIGATION USE1

ABSTRACT: Artificial recharge as a means of augmenting water sup plies for irrigation is a management alternative which policy makers in ground water decline areas are beginning to consider seriously. A conceptual model is developed to evaluate the economic benefits from ground water recharge under conditions where the major water use is irrigation. The methodology presented separates recharge benefits into two components: pumping cost savings and aquifer extension benefits. This model is then applied to a Nebraska case to approximate the value of recharge benefits as a function of aquifer response. discount rate, and commodity prices. It was found that recharge benefits vary from less than $2 to over $6 an acre foot recharged.     

HYDROGEOLOGY AND PATHS OF FLOW IN THE CARBONATE BEDROCK AQUIFER,NORTHWESTERN INDIANA1

ABSTRACT: The U.S. Geological Survey (USGS) is assessing the ground-water resources of the carbonate bedrock aquifers in Indiana and Ohio as part of their Regional Aquifer Systems Analysis program. Part of this assessment includes the determination of unknown aspects of the hydraulic characteristics, boundaries, and flow paths of the carbonate aquifer. To accomplish this, the USGS drilled three wells through the carbonate aquifer near the Kankakee River in northwestern Indiana. Geophysical logs were used to help determine depths and thicknesses for testing and to help describe geology at the three wells. Packer tests were used to determine direction of ground-water flow and to provide data for an analysis of the distribution of transmissivity in the carbonate aquifer. Transmissivity of the carbonates is associated with two physical characteristics of the rocks: fractures and interconnected porosity. Almost all of the transmissivity is derived from horizontal fracturing; however, only a few of the fractures present in the carbonate are transmissive. Some transmissivity is associated with a zone of fossiliferous, vuggy dolomite, which yields water from the rock matrix. Most of the transmissivity is associated with large fractures and solution crevices in the upper 30 feet of the bedrock; less transmissivity is associated with the deeper vuggy reef material, even where extensively fractured. Transmissivity of individual fractures and fossiliferous zones ranges from 300 to 27,000 feet squared per day. The aquifer bottom is defined by a lack of transmissive fractures and an increased shale content near the contact of the Silurian and Ordovician sections. Water-level data from the three wells indicate that flow is horizontal at well site 1 north of the Kankakee River, upward at well site 2 near the river, and downward at well site 3 south of the river. Most of the flow occurs in the upper part of the carbonate bedrock where fracturing and solution-enlarged crevices are most developed. Water levels indicate the the Kankakee River is a hydrologic boundary for the regional carbonate bedrock aquifer.