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.     

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.     

CHARACTERIZATION OF GROUND WATER FLOW FROM SPRING DISCHARGE IN A CRYSTALLINE ROCK ENVIRONMENT1

ABSTRACT: Recent investigations describing the hydrogeology of the Blue Ridge Province of Virginia suggest the occurrence of multiple aquifers and flow paths that may be responsible for the variable flow behavior of springs and seeps appearing throughout the region. Deep, confined aquifers associated with ubiquitous faults and shallow, variably confined saprolite aquifers may contribute water to spring outlets resulting in significantly different quantities of discharge and water quality. Multiple analyses are required to adequately identify the flow paths to springs. In this investigation, hydrograph analyses, surface electrical resistivity surveys, aquifer tests, and nitrate concentrations are used in conjunction with previously reported analyses from borehole logs and age dating of ground water to identify two distinct flow paths. Results indicate that base flow occurs from a deep fault zone aquifer and such discharge can be maintained even during prolonged periods of drought, while increased discharge identified on hydrograph peaks suggests the occurrence of rapid flow through the saprolite aquifer within a radius of about 25 meters of the spring orifice. Springflow hydrograph analysis is suitable for rapid characterization of flow paths leading to spring outlets. Rapid characterization is important for evaluation of potential water quality problems arising from contamination of shallow and deep aquifers and for evaluation of water resource susceptibility to drought. The techniques evaluated here are suitable for use in other locations in fractured crystalline rock environments.     

HYDROLOGIC ANALYSIS OF DISCHARGE SUSTAINABILITY FROM AN ABANDONED UNDERGROUND COAL MINE1

ABSTRACT: Discharge from flooded abandoned subsurface coal mines is considered a potential source for water supplies where other acceptable water sources are not available. The objective of this study was to develop procedures for determining sustainability of mine‐water discharge using rainfall and discharge data for a case study site. The study site is located in southwest Virginia where Late Paleozoic sequences of sandstone, coal, and shale predominate. A rain gauge and a flow rate monitoring system were installed at the site and data were collected for a period of 100 days. The recording period corresponded with one of the driest periods in recent years and, therefore, provided valuable information regarding the flow sustainability during baseflow conditions. From available data on underground mining patterns, geology, and ground water flow regimes, it was determined that a coal mine aquifer exhibits hydraulic characteristics very similar to the extremely heterogeneous systems observed in karst aquifers, and the mine discharge is analogous to springflow. Thus, techniques commonly used in karst‐water systems and springflow analysis were used to develop rainfall/mine‐discharge relationships. Springflow recession analysis was performed on five rainfall recessions and the coefficient for each recession was compared and interpreted in light of known geologic information. It was found that the recession coefficients described the mine discharge adequately and the mine aquifer response to a rainfall pulse was very similar to the response from certain types of karst aquifers. A cross‐correlation analysis was performed to verify the results of the recession analysis and to develop a “black box” statistical model for discharge data. The correlation analysis proved the validity of springflow recession analysis for mine discharge. The recorded data length was not adequate to create a statistical model, however, but a procedure was proposed for a statistical model that could be used with large flow records. For the study site, the mine discharge was found to be sustainable for a prolonged period of time.     

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.     

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.     

PARAMETER ESTIMATION FOR AQUIFER EVALUATION1

ABSTRACT A discussion of the various tools and techniques available for evaluation of aquifers is presented. Data for evaluation studies can be obtained from laboratory experiments, geological maps, well logs and field studies. Such data can be used in analytic or electric analog models to determine the sustained yield of the aquifer under various operating conditions.     

A SYSTEMS STUDY OF FUTURE GROUNDWATER DEVELOPMENT COSTS IN THE CHICAGO REGION1

The cost of developing groundwater resources in northeastern Illinois from 198cL2020 is estimated for the purpose of providing a basis for comparing alternative sources. Demands for each township in the study area are estimated at 10-year increments and are satisfied, where the supply is sufficient, in such a way as to minimize the cost subject to constraints on supply. Sources of water are two shallow aquifers with known potential yields and a series of deep aquifers treated as a single unit and modeled on a digital computer. For each township the costs of wells, pumps, power and rehabilitation is estimated for each aquifer on a per million gallons of water per day basis. In addition the cost of groundwater treatment necessary to raise the quality to that of treated Lake Michigan water is considered. Raw water costs are found to vary from 2 to 14 cents per 1000 gallons depending upon the depth to the deep aquifer water. Treated water costs vary from 22 to 53 cents per 1000 gallons, the lower costs applying to the largest users because of the economy of scale. It is found that with proper distribution of pumpage there is sufficient water in storage in the deep aquifers to meet groundwater demands through 2020.     

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.     

EVALUATING THE INFLUENCE OF SOURCE BASINS ON DOWNSTREAM WATER QUALITY IN THE MISSISSIPPI RIVER1

ABSTRACT: Chemical variability in the Mississippi River during water years 1989 to 1998 was evaluated using stream discharge and water‐quality data in conjunction with the DAFLOW/BLTM hydraulic model. Model simulations were used to identify subbasin contributions of water and chemical constituents to the Mississippi River upstream from its confluence with the Ohio and the Mississippi River and at the Atchafalaya Diversion in Louisiana. Concentrations of dissolved solids, sodium, and sulfate at the Thebes site showed a general decreasing trend, and concentrations of silica and nitrate showed a general increasing trend as the percentage of discharge from the Mississippi River upstream from Grafton increased. Concentrations of most chemical constituents in the Mississippi River at the Atchafalaya Diversion exhibited a decreasing trend as the percentage of water from the Ohio River increased. Regression models were used to evaluate the importance of the source of water to the water chemistry in the Mississippi River at Thebes and the Atchafalaya Diversion. The addition of terms in regression equations to account for the percent of water from sub‐basins improved coefficients of determination for predicting chemical concentrations by as much as nine percent at the Thebes site and by as much as 48 percent at the Atchafalaya Diversion site. The addition of source‐water terms to regression equations increased the estimated annual loads of nitrate and silica delivered from the Mississippi River Basin to the Gulf of Mexico by as much as 14 and 13 percent, respectively.     

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.     

OPTIMAL GROUND WATER MANAGEMENT IN TWO-AQUIFER SYSTEMS1

ABSTRACT; This paper presents a numerical model for the prediction of optimal ground water withdrawal from a two-aquifer system by observing a set of constraints determined by the ecological conditions of the ground water basin. The aquifer system consists of an upper unconfined and a lower confined aquifer with a leaky stratum between them. It is assumed that water is withdrawn from the confined aquifer only, but the unconfined aquifer will also be affected due to the leakiness of the layer separating the upper and lower aquifers. Simulation and linear programming are employed for developing a computer model for the optimal management of such systems, with the objectives of determining withdrawal rates for predetermined ground water levels.     

CHARACTERIZING GROUND WATER FLOW IN THE MUNICIPAL WELL FIELDS OF CEDAR RAPIDS,IOWA, WITH SELECTED ENVIRONMENTAL TRACERS1

ABSTRACT: Cedar Rapids obtains its municipal water supply from a shallow alluvial aquifer along the Cedar River in east-central Iowa. Water samples were collected and analyzed for selected isotopes and chlorofluorocarbons to characterize the ground-water flow system near the municipal well fields. Analyses of deuterium and oxygen-18 indicate that water in the alluvial aquifer and in the underlying carbonate bedrock aquifer was recharged from precipitation during modern climatic conditions. Analyses of tritium indicate modern, post-1952, water in the alluvial aquifer and older, pre-1952, water in the bedrock aquifer. Mixing of the modern and older waters occurs in areas where (1) the confining layer between the two aquifers is discontinuous, (2) the bedrock aquifer is fractured, or (3) pumping of supply wells induces the flow of water between aquifers. Analyses of chlorofluorocarbons were used to determine the date of recharge of water samples. Water in the bedrock aquifer likely was recharged prior to the 1950s. Water in the alluvial aquifer likely was recharged from the 1960s to 1990s. Biodegradation or sorption probably affected some of the ground water analyzed for chlorofluorocarbons. These processes reduce the concentrations of CFCs, which results in older than actual calculated dates of recharge.     

HISTORICAL CHANGES IN UPPER MISSISSIPPI RIVER WATER AREAS AND ISLANDS1

ABSTRACT: Periodic surveys of the upper Mississippi River since 1866 and a discharge record of nearly equal length provided an opportunity to learn more about the magnitudes and rates of geomorphic processes at work in large stream systems. Furthermore, geomorphic and hydrologic adjustments could be evaluated in relation to watershed land use changes, small‐scale climate fluctuations, and considerable modifications to the channel and floodplain during the period of record. The present study uses GIS mapping to quantitatively compare historical changes in mapped land and water phenomena in the upper Mississippi River Pool 10, located along southwest Wisconsin's border. Modest channel widening and decreases in island area throughout the study reach during the last century are detectable. Flood magnitudes and frequencies also have varied during this time, and stages and low flow discharges have increased since the 1940s. The latter hydrologic change appears to be closely associated with the reach's geomorphic adjustments. Results are representative of a valley reach where a major tributary contributes a large sand bedload, forming an alluvial fan of considerable size in the floodplain.     

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.     

HYDROLOGIC STAGE PERIODICITY OF THE MISSISSIPPI RIVER BEFORE AND AFTER SYSTEMATIC CHANNEL MODIFICATIONS1

ABSTRACT: Periodic flood disturbance is a well known controlling factor of in channel and floodplain ecosystem function. However, channel manipulations during the last century have potentially altered hydrologic fluctuations, and thus ecosystem function. We examined temporal river stage hydrology, through autocorrelation analysis, at seven gauges along the Mississippi River to quantify flow periodicity and effects of systematic channel modifications on flow periodicity. Intraannual variation follows a strong one‐year cycle of six months higher flow and six months lower flow for the entire Mississippi River drainage, with precipitation as a driving force. Interannual hydrologic variation differs between the upper and lower river segments. A clear quasi‐biennial oscillation pattern was evident throughout the lower river section. The effect of channel alterations was a decreased magnitude of differences between lower and higher flows. The upper section, however, suggests a 12‐to 14‐year periodicity prior to alterations and a decreased duration of lower flow years following systematic modifications. Interannual variograms clearly depict very different temporal hydrology between the upper Mississippi River and the lower Mississippi River, suggesting the simple transfer of knowledge from one segment to the other oversimplifies the complexity of a large river system.     

AN EXAMINATION OF THE MECHANISMS CONTROLLING GROUNDWATER GRADIENTS IN HYPER-ARID REGIONAL SEDIMENTARY BASINS1

ABSTRACT: Although evidence of modern recharge in the North African and Arabian sedimentary basin aquifers exists, it is difficult to determine the volume of recharge. Also, from the evidence of regional groundwater gradients, the flow within the aquifers seems to be appreciably greater than one would intuitively expect. A hypotehtical model embodying the characteristics of the aquifers has been used to investigate the likely significance of various possible flow mechanisms. It is shown that while dewatering in the unconfined area can possibly contribute to flows for a considerable period of time, the maintenance of water levels in the unconfined zone must be the result of modern recharge. It is also shown that recharge depths of less than 10 mm per annum are sufficient given suitable aquifer parameters. Results for various combinations of aquifer parameters and configurations are given, including layered aquifers and the effects of restricted oufflows. Comparisons are made using a “bench mark” example. The work indicates that there is little point in carrying out conventional hydrological balance studies in hyper-arid areas and that, instead, more emphasis should be placed upon good groundwater hydrographic data and modeling.     

A VERIFICATION OF THE QUAL-1 WATER QUALITY MODEL FOR THE LOWER MISSISSIPPI RIVER1

ABSTRACT: This paper describes the verification of the QUAL-1 mass transport model for the lower Mississippi River between St. Francisville and Point a la Hache using dye studies conducted by the U. S. Geological Survey. QUAL-1 is a one-dimensional steady-state model for rivers and is capable of predicting longitudinal profiles of soluble materials entering rivers from point sources. Both conservative and nonconservative parameters of water quality can be considered. The major problems surmounted were the determination of a diffusion coefficient and the use of transient data to verify a steady state model.     

NUTRIENT LOAD CHARACTERIZATION FROM INTEGRATED SOURCE DATA FOR THE LOWER MISSISSIPPI RWER1

ABSTRACT: Nutrient data from all available sources for the lower Mississippi River were examined for potential differences among sampling agencies and geographic locations for the period between 1960 and 1998. Monthly means grouped by parameter, sampling location and agency, were calculated and compared as paired sets, excluding those months where data were not available for both sets. Some significant differences were found between various agencies collecting nutrient data on the river, as well as between various stretches of river, especially in the case of phosphorus nutrient data. Results were used to synthesize data sets from which a history of nutrient loading in the Mississippi River was determined. General trends in nitrate+nitrite, total Kjeldahl nitrogen, orthophosphate, total phosphorus and silica loads, as well as changes in nutrient proportions and the specific limiting nutrient (by month) are reported. This study provides a useful summary of contemporary and historical nutrient data that may assist in the evaluation of Mississippi River water quality and its potential effect on the Gulf of Mexico.     

Mapping the vulnerability of groundwater to the contamination of four carbonate aquifers in Europe

The vulnerability of four European aquifers with different hydrogeological and climatic characteristics was evaluated using the COP method. The results obtained were statistically analyzed by determination coefficients to measure which factor has greater importance in the vulnerability index. Furthermore, a new parameter has been designed to measure the vulnerability for the whole of the aquifer. The results demonstrate that COP is a useful method to assess the vulnerability of the test sites under consideration. The results obtained are coherent with the conceptual model of each pilot aquifer and the available hydrogeological information (hydrographs, isotopic data, tracer tests). Fissured carbonate aquifers (diffuse flow systems) are less vulnerable than karst aquifers (conduit flow systems) and the vulnerability index is more positively correlated with the O factor (unsaturated zone protection capacity) in the first case. The karst aquifers are more vulnerable than fissured aquifers and they show a higher correlation between the C factor (karst features) and the vulnerability index. Climatic variation (precipitation for example) influences the final vulnerability index of the aquifers according to the weight in the index and the spatial distribution.