What is a static head measurement? Evaluating structural error and implications for regional groundwater modeling

Most groundwater modelers avoid using static heads measured from active production wells because they can introduce a bias into model calibration. However, in the deep confined Cambrian-Ordovician Sandstone Aquifer System in the Central Midcontinent of North America, dedicated observation wells are sparse and remote from areas of most concentrated pumping. As a result, in areas where drawdown is the greatest and modeling is most needed, only static heads from production wells are available for calibration. This paper evaluates two leading sources of discrepancies in using production well data, spatial and temporal structural error (S.E.). A simple Theis solution is used to evaluate the potential magnitude of spatial S.E. when calibrating a regional MODFLOW model with coarse cell resolution. Despite theoretical analyses indicating that spatial S.E. could be significant, statistical analysis of the model results suggests that temporal S.E. is dominant. Long (ranging over decades) or frequent (monthly) head datasets are key in understanding temporal S.E., to better capture water-level variability. In this study, the range in static head observations impacted estimates of the remaining time a well could extract water from the aquifer by 0.1 to 16.0 years. This uncertainty in future water supply is highly relevant to stakeholders and must be assessed in hydrographs depicting risk.     

USE OF A GEOGRAPHIC INFORMATION SYSTEM TO ASSEMBLE INPUT-DATA SETS FOR A FINITE-DIFFERENCE MODEL OF GROUND-WATER FLOW1

ABSTRACT: A Geographic Information System (GIS) was used to expedite assembly of input-data sets for a model of ground-water flow in the Middle Patuxent River basin in Howard County, Maryland. The model grid was developed with GIS and used to select attributes from GIS data coverages. These attributes were then output from GIS into the input-data sets, which included model boundaries, regolith thickness, bedrock-surface altitude, stream locations, stream length, stream-bottom altitude, and trans-missivity. The ability to change large sets of spatial data quickly and accurately with GIS enhances the model-calibration process. (KEY TERMS: GIS; ground water; Piedmont; Maryland; MODFLOW.)     

The Transboundary Nature of the Allende–Piedras Negras Aquifer Using a Numerical Model Approach

The Allende–Piedras Negras (APN) aquifer is located between the states of Texas (United States [U.S.]) and Coahuila (Mexico). The Rio Grande crosses the aquifer, acting as a natural and political divide between the countries. However, it remains unclear whether the APN aquifer can be considered a truly transboundary aquifer flow system, which would potentially require joint management by two different administrative jurisdictions. The main purpose of this study was to evaluate the transboundary nature of this aquifer. This was achieved by developing a detailed hydrogeological model to analyze the direction of volumetric fluxes within the APN aquifer using Visual MODFLOW. The model simulated a spatially averaged cumulative drawdown of 0.76 m for the entire aquifer over an 18‐year modeling period (2000–2017). The flow convergence zone, previously located below the Rio Grande, has shifted to the U.S. side in most locations, driven by higher pumping rates of the wells located near the river. This shift of the convergence zone from one country to the other means that groundwater recharge from one side flows underneath the river to the other side. This qualifies the APN aquifer as a “transboundary groundwater flow system.” The procedure followed in this study may be applied to other aquifers that straddle the U.S.–Mexico border and may motivate future modeling studies on other poorly studied transboundary aquifers around the world and thereby enable bi‐national aquifer management.     

A Semianalytical Solution for Residual Drawdown at a Finite Diameter Well in a Confined Aquifer

After the end of pumping the water level in the observation well starts to recover and the reduced drawdown during the recovery period is named as the residual drawdown. Traditional approaches in analyzing the data of residual drawdown for estimating the aquifer hydraulic parameters are mostly based on the application of superposition principle and Theis equation. In addition, the effect of wellbore storage is commonly ignored in the evaluation even if the test well has a finite diameter. In this article, we develop a mathematical model for describing the residual drawdown with considering the wellbore storage effect and the existing drawdown distribution produced by the pumping part of the test. The Laplace‐domain solution of the model is derived using the Laplace transform technique and the time‐domain result is inverted based on the Stehfest algorithm. This new solution shows that the residual drawdown associated with the boundary and initial conditions are related to the well drawdown and the aquifer drawdown, respectively. The well residual drawdown will be overestimated by the Theis residual drawdown solution in the early recovery part if neglecting the wellbore storage. On the other hand, the Theis residual drawdown solution can be used to approximate the present residual drawdown solution in the late recovery part of the test.     

UTILIZING INDUCED RECHARGE FOR REGIONAL AQUIFER MANAGEMENT1

ABSTRACT: The deep aquifers of the Portland Basin are used as a regional water supply by at least six municipalities in Oregon and Washington. Maximum continuous use of the aquifers in 1998 was 13 mgd and peak emergency use was 55 mgd. Continuous use of the deep aquifers at a rate of 55 mgd has been proposed and inchoate water rights have been reserved for expansion of pumping to 121 mgd. A study was completed, using a calibrated ground water flow model, to evaluate the role of induced recharge from the Columbia River in mitigating aquifer drawdown from continuous‐use and expanded pumping scenarios in the center and eastern areas of the basin. The absolute average residual was less than 3.6 feet for steady‐state model calibrations, and less than 8.0 feet for transient calibration to a 42 mgd pumping event in 1987 with 170 feet of drawdown. Continuous use of the aquifers at a rate of 55 mgd is predicted to increase drawdown to 210 feet. Expansion of pumping to 121 mgd in the center basin is predicted to cause 400 feet of drawdown. However, expansion of pumping in the east basin is predicted to result in only 220 feet of drawdown because of induced recharge from the Columbia River.     

SWATMOD‐Prep: Graphical User Interface for Preparing Coupled SWAT‐MODFLOW Simulations

This article presents SWATMOD‐Prep, a graphical user interface that couples a SWAT watershed model with a MODFLOW groundwater flow model. The interface is based on a recently published SWAT‐MODFLOW code that couples the models via mapping schemes. The spatial layout of SWATMOD‐Prep guides the user through the process of importing shape files (sub‐basins, hydrologic response units [HRUs], river network) from an existing SWAT model, creating a grid, performing necessary geo‐processing operations to link the models, writing out SWAT‐MODFLOW files, and running the simulation. The option of creating a new single‐layer MODFLOW model for near‐surface alluvial aquifers is available, with the user prompted to provide groundwater surface elevation (through a digital elevation model), aquifer thickness, and necessary aquifer parameter values. The option of simulating nitrate transport in the aquifer also is available, using the reactive transport model RT3D. The interface is in the public domain. It is programmed in Python, with various software packages used for geo‐processing operations (e.g., selection, intersection of rasters) and inputting/outputting data, and is written for Windows. The use of SWATMOD‐Prep is demonstrated for the Little River Experimental Watershed, Georgia. SWATMOD‐Prep and SWAT‐MODFLOW executables are available with an accompanying user's manual at: http://swat.tamu.edu/software/swat-modflow/ . The user's manual also accompanies this article as Supporting Information.     

INTEGRATING GEOGRAPHIC INFORMATION SYSTEMS IN GROUND-WATER APPLICATIONS USING NUMERICAL MODELING TECHNIQUES1

ABSTRACT: The Sand and Gravel Aquifer is the sole source of potable water in Escambia County, Florida. In order to better understand the hydraulics of the aquifer, a numerical computer model of the aquifer was developed. The model applied a finite element technique which allowed for density-dependent transport and flow in three dimensions. The modeling technique was integrated with GIS to develop a system for optimal management of the resource. The GIS was the primary tool in the development of the model grid, as well as being the integral component in the modeling procedure. Multiple model grids were developed for simulating regional flow and local flow/transport phenomena. The model grids were generated by the GIS where nodal and element sequencing were recorded. The grid topology was stored in the GIS with the element numbers, node numbers, and the related hydrogeologic attributes. The Triangulated Irregular Network (TIN) module was used for transferring interpolated value between GIS coverages. TIN allowed a fit of the model grid to the physical dimensions of the aquifer and for interpolating boundary values for telescopically refined grids. Calculations between TIN surfaces provided the residuals of the dependent variable from observed TIN surfaces. Model calibration was conducted within the GIS environment through a combination of visual and relational querying. The GIS provided an integrated environment which facilitated model analyses and data storage and retrieval.     

SUMMARY OF MODIFICATIONS OF THE U.S. GEOLOGICAL SURVEY MODULAR,FINITE-DIFFERENCE,GROUND-WATER FLOW MODEL TO READ AND WRITE GEOGRAPHIC INFORMATION SYSTEM FILES1

ABSTRACT: The U.S. Geological Survey modular, three-dimensional, finite-difference, ground-water flow model, commonly called MODFLOW, has been modified so that it can read and write files used by a geographic information system (GIS). The modified model program is called MODFLOWARC. The design of MODFLOWARC parallels the design of the ground-water flow model program MODFLOW. The names of the variables, modules, and submodules used to explain the operations of MODFLOWARC were derived from the names used in MODFLOW. During the data input phase, MODFLOWARC reads array control records similar to the original control records of MODFLOW, except an additional variable is added. This additional variable is the name of the computer files containing array data in GIS format. Data output is achieved by setting record/input flags and by supplying a variable that is the name of the directory where the output data will be recorded. The modifications to MODFLOW were minimized so that MODFLOWARC will operate on an existing ground-water flow model without modifying array control records.     

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.     

RECHARGE TO ALLUVIAL VALLEY AQUIFERS FROM OVERBANK FLOW AND EXCESS INFILTRATION1

ABSTRACT: Recharge is an important parameter for models that simulate water and contaminant transport in unconfined aquifers. Unfortunately, measurements of actual recharge are not usually available causing recharge to be estimated or possibly added to the calibration procedure. In this study, differences between observed water-table elevations and water-table elevations simulated with a model based on the one-dimensional Boussinesq equation were used to identify both the timing and quantity of recharge to an alluvial valley aquifer. Observed water table elevations and river stage data were recorded during a five-year period from 1991 to 1995 at the Ohio Management Systems Evaluation Area located in south-central Ohio. Direct recharge attributed to overbank flow during and shortly after flood conditions accounted for 65 percent of the total recharge computed during the five-year study period. Recharge of excess infiltration to the aquifer was intermittent and occurred soon after large rainfall events and high river stage. Specification of constant recharge with time values in ground-water simulation models seems inappropriate for stream-aquifer systems given the strong influence of the river on water table elevations in these systems.     

Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas1

Mittelstet, Aaron R., Michael D. Smolen, Garey A. Fox, and Damian C. Adams, 2011. Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas. Journal of the American Water Resources Association (JAWRA) 1‐8. DOI: 10.1111/j.1752‐1688.2011.00524.x Abstract: We compared two approaches to administration of groundwater law on a hydrologic model of the North Canadian River, an alluvial aquifer in northwestern Oklahoma. Oklahoma limits pumping rates to retain 50% aquifer saturated thickness after 20 years of groundwater use. The Texas Panhandle Groundwater Conservation District’s (GCD) rules limit pumping to a rate that consumes no more than 50% of saturated thickness in 50 years, with reevaluation and readjustment of permits every 5 years. Using a hydrologic model (MODFLOW), we simulated river‐groundwater interaction and aquifer dynamics under increasing levels of “development” (i.e., increasing groundwater withdrawals). Oklahoma’s approach initially would limit groundwater extraction more than the GCD approach, but the GCD approach would be more protective in the long run. Under Oklahoma rules more than half of aquifer storage would be depleted when development reaches 65%. Reevaluation of permits under the Texas Panhandle GCD approach would severely limit pumping as the 50% level is approached. Both Oklahoma and Texas Panhandle GCD approaches would deplete alluvial base flow at approximately 10% development. Results suggest periodic review of permits could protect aquifer storage and river base flow. Modeling total aquifer storage is more sensitive to recharge rate and aquifer hydraulic conductivity than to specific yield, while river leakage is most sensitive to aquifer hydraulic conductivity followed by specific yield.     

A GIS-based DRASTIC vulnerability and net recharge reassessment in an aquifer of a semi-arid region (Metline-Ras Jebel-Raf Raf aquifer, Northern Tunisia)

This paper aims to elaborate new generic DRASTIC aquifer vulnerability maps of the coastal aquifer of Metline-Ras Jebel-Raf Raf (Northeast of Tunisia) using the GIS technique, making the data analyses easier to handle and providing better capabilities of dealing with large spatial data. A similar study was carried out in 1999 in the same aquifer using a method based on the soil water balance equation to determine the net recharge parameter. Unfortunately, the lack of data in the study area made the results unsatisfactory. By applying the Williams and Kissel equation and the Rao relationship, we intend to demonstrate that we could correctly evaluate the net recharge parameter. Moreover, new data related to the aquifer hydraulic conductivity, the soil cover and the vadose zone lithology have become available, and allowed us to elaborate suitable DRASTIC maps.     

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.     

DETERMINATION OF HYDRAULIC C.ONDUCTWITY TENSOR USING A NONLINEAR LEAST SQUARES ESTIMATOR1

ABSTRACT: Reliable and consistent estimation of the components of the hydraulic conductivity tensor provides information needed to make proper decisions regarding clean up and restoration of contaminated aquifers. In this study, the nonlinear least-squares estimation technique was applied to drawdown versus time data from three or more observation wells to determine a theoretical ellipse of equal drawdown. The angle of rotation of this ellipse with respect to the working coordinate axes was determined by a procedure based on contouring the drawdowns at a given time. This ellipse, in turn, was used to estimate the directions and magnitudes of the horizontal components of the hydraulic conductivity tensor. The technique is applicable to confined, as well as leaky, aquifers. Sources of error in this technique include nonhomogeneity of the aquifer and partial penetration of the pumping and observation wells into the aquifer. The procedure presented may be used as an additional tool to verify computations of hydraulic conductivity anisotropy based on other techniques.     

ESTIMATION OF RECHARGE FROM A RESERVOIR USING TWO WATER BUDGET MODELS1

ABSTRACT: Reservoir water levels, observation well data, and meteorological parameters were collected at a recharge dam site in Central Saudi Arabia. This data, along with other information on the reservoir and the underlying aquifer, were used to estimate the amounts of recharge through the reservoir bed by applying two water budget models. The first is a water budget model for the reservoir only, while the second is for an aquifer reach extending upstream and downstream from the reservoir. The results of the two approaches were discussed and compared.     

DETERMINATION OF UNCONFINED AQUIFER HYDRAULIC PROPERTIES FROM RECOVERY TEST DATA1

ABSTRACT: This paper analyzes the sensitivity of drawdown to four hydraulic parameters in unconfined aquifers: horizontal and vertical hydraulic conductivity K_r and K_z, storage coefficient S, and specific yield S_y. Sensitivity coefficients indicate that the sensitivity vanes with time for each aquifer parameter, and K_r, K_z, S, and S_y are identifiable from recovery test data. An inverse method was used to calculate the four parameters from residual drawdowns. Results of application examples demonstrate that residual data provide valid information in the determination of unconfined aquifer hydraulic parameters.     

DEVELOPMENT OF A WATER SUPPLY PROTECTION MODEL IN A GIS1

ABSTRACT: As part of a larger model to identify lands suitable for acquisition, a water supply protection model was developed using the Southwest Florida Water Management District's GIS. Several hydrologic and hydrogeologic data layers were overlaid to develop maps showing ground-water supply suitability, protection areas for surface-water supply, protection areas for major public supply wells, susceptibility to ground-water contamination, and recharge to the Floridan aquifer. These intermediate layers were combined into a final map to prioritize protection areas for water supply.     

CONTINENTAL SCALE SIMULATION OF THE HYDROLOGIC BALANCE1

This paper describes the application of a continuous daily water balance model called SWAT (Soil and Water Assessment Tool) for the conterminous U.S. The local water balance is represented by four control volumes; (1) snow, (2) soil profile, (3) shallow aquifer, and (4) deep aquifer. The components of the water balance are simulated using “storage” models and readily available input parameters. All the required databases (soils, landuse, and topography) were assembled for the conterminous U.S. at 1:250,000 scale. A GIS interface was utilized to automate the assembly of the model input files from map layers and relational databases. The hydrologic balance for each soil association polygon (78,863 nationwide) was simulated without calibration for 20 years using dominant soil and land use properties. The model was validated by comparing simulated average annual runoff with long term average annual runoff from USGS stream gage records. Results indicate over 45 percent of the modeled U.S. are within 50 mm of measured, and 18 percent are within 10 mm without calibration. The model tended to under predict runoff in mountain areas due to lack of climate stations at high elevations. Given the limitations of the study, (i.e., spatial resolution of the data bases and model simplicity), the results show that the large scale hydrologic balance can be realistically simulated using a continuous water balance model.     

EDWARDS AQUIFER EVALUATION: KINNEY COUNTY,TEXAS1

ABSTRACT: The Edwards Aquifer is one of the most studied and most prolific aquifers in the United States. The aquifer is a heavily fractured and faulted carbonate aquifer with transmissivities in excess of 100 ft²/s. The City of San Antonio relies upon the Edwards Aquifer as its sole source for water. Much work has been done on quantifying recharge to the aquifer and discharge from wells and acquiring aquifer characteristics from pumping tests, specific capacity tests, and geophysical logs. Although the aquifer has been well studied in Bexar County, much less is known about the Edwards Aquifer in Kinney County. This is partly due to the lower population within the county (approximately 3,500 people) relative to the eastern counties (Uvalde, Medina, Bexar, Comal, and Hays) and the great distance of Kinney County from high profile discharge areas such as the City of San Antonio and Comal and San Marcos Springs. Three key products resulted from this study: (1) exploratory well drilling and the largest aquifer test in the county that were conducted to evaluate the well yields within a 10,000 acre study area in which a drawdown of 2.5 ft approximately 1.2 miles away was observed while pumping at approximately 4,600 gpm; (2) a recharge estimate for the Edwards Aquifer within Kinney County of approximately 71,382 ac‐ft/yr; and (3) locating the Brackettville Groundwater Divide from an evaluation of ground water flow direction and hydrograph analysis. These results help evaluate the complex hydraulics occurring within Kinney County and aid in development of ground water modeling that will be used in managing the Edwards Aquifer.     

USING SENSITIVITY ANALYSIS TO ASSIST PARAMETER ZONATION IN GROUND WATER FLOW MODEL1

ABSTRACT: For numerical modeling of ground water movement in a real aquifer system, the aquifer is usually divided into hydrogeologically defined zones, each with its own parameter values. The responses of the system, such as head or drawdown, are often available only in some of the zones. The estimated parameters of all the zones are based on the measured response in these limited zones. However, the estimates for some of the zones may be very uncertain, and these zones are therefore not justified by the data. In this paper, an approach is presented to understand which zone may produce uncertain parameter values and should be lumped with its neighbor. This approach is demonstrated using a regional numerical model for pumping test analysis in the Nottinghamshire aquifer, UK. A step-by-step process is used in identifying the aquifer zones and estimating their parameters based on the principle of using the smallest possible numbers of zones and parameters for adequate representation of the drawdown response. After the parameters of each zone are estimated, the sensitivity features of these parameters are examined. The results show that the parameters in one zone can be estimated properly by the drawdown in another zone only when there is significant sensitivity. For transmissivity, sensitivity between zones occurs when there is significant flow between them. For storativity, sufficient sensitivity can occur without large flows between the zones, provided that one zone causes significant drawdown in the other. This idea can be extended to the flow model for a large aquifer system. If the aquifer is divided in such a way that aquifer responses are not sensitive to the parameters in some of the zones, the parameters in those zones cannot be estimated properly and should be lumped into their neighboring zones. In this way, a simple but more reasonable model can be built.