MTBE and gasoline hydrocarbons in ground water of the United States

The occurrence of methyl tert-butyl ether (MTBE) and gasoline hydrocarbons was examined in three types of studies of ground water conducted by the U.S. Geological Survey: major aquifer surveys, urban land-use studies, and agricultural land-use studies. The detection frequency of MTBE was dependent on the study type, with the highest detection frequency in urban land-use studies. Only 13 ground water samples from all study types, or 0.3%, had concentrations of MTBE that exceeded the lower limit of the U.S. EPA's Drinking-Water Advisory. The detection frequency of MTBE was highest in monitoring wells located in urban areas and in public supply wells. The detection frequency of any gasoline hydrocarbon also was dependent on study type and generally was less than the detection frequency of MTBE. The probability of detecting MTBE in ground water was strongly associated with population density, use of MTBE in gasoline, and recharge. Ground water in areas with high population density, in areas where MTBE is used as a gasoline oxygenate, and in areas with high recharge rates had a greater probability of MTBE occurrence. Also, ground water from public supply wells and shallow ground water underlying urban land-use areas had a greater probability of MTBE occurrence compared to ground water from domestic wells and ground water underlying rural land-use areas. The probability of detecting MTBE in ground water was weakly associated with the density of leaking underground storage tanks, soil permeability, and aquifer consolidation, and only concentrations of MTBE >0.5 microg/L were associated with dissolved oxygen.     

Variations of bromide in potable ground water in the United States

Concentrations of bromide in potable ground water that has <10 mg/L chloride range from 0.0032 to 0.058 mg/L with a median value of 0.016 mg/L. The chloride/bromide mass ratio for the same water ranges from 43 to 285 with a median value of 101. The ratios, which resulted from screening approximately 165 analyses of water from 32 locations in 24 states in the United States, show a distinct geographic variation with highest values near the coast and trending toward a value of approximately 50 in the continental interior.     

Michigan basin regional ground water flow discharge to three Great Lakes

Ground water discharge to the Great Lakes around the Lower Peninsula of Michigan is primarily from recharge in riparian basins and proximal upland areas that are especially important to the northern half of the Lake Michigan shoreline. A steady-state finite-difference model was developed to simulate ground water flow in four regional aquifers in Michigan's Lower Peninsula: the Glaciofluvial, Saginaw, Parma-Bayport, and Marshall aquifers interlayered with the Till/"red beds," Saginaw, and Michigan confining units, respectively. The model domain was laterally bound by a continuous specified-head boundary, formed from lakes Michigan, Huron, St. Clair, and Erie, with the St. Clair and Detroit River connecting channels. The model was developed to quantify regional ground water flow in the aquifer systems using independently determined recharge estimates. According to the flow model, local stream stages and discharges account for 95% of the overall model water budget; only 50% enters the lakes directly from the ground water system. Direct ground water discharge to the Great Lakes' shorelines was calculated at 36 m3/sec, accounting for 5% of the overall model water budget. Lowland areas contribute far less ground water discharge to the Great Lakes than upland areas. The model indicates that Saginaw Bay receives only approximately 1.13 m3/sec ground water; the southern half of the Lake Michigan shoreline receives only approximately 2.83 m3/sec. In contrast, the northern half of the Lake Michigan shoreline receives more than 17 m3/sec from upland areas.     

Heat flow and ground water movement in the Bohemian cretaceous basin (Czechoslovakia)

Sub-surface temperature fields may be considerably affected by active ground water systems, thereby seriously hampering the interpretation of heat flow data. Quantitative evaluation of the convective component of heat transfer is thus very important in cases such as large sedimentary basins with vast underground water circulation. We propose in this study a simple model of horizontal aquifer. This model was used to examine the effect of the lateral convection on the surface heat flow near the recharge zone of basinal margins. The perturbation of the heat flow field above the aquifer was calculated for various aquifer geometry and various flow velocities and the regional scale dependence of the perturbation on the hydraulic properties of the aquifer was demonstrated. The model was applied to the Bohemian Cretaceous Basin and it was shown that within a few kilometres from the recharge zones the observed surface heat flow may be underestimated by up to several tens of percent. The procedure was tested in two locations in this area, in an attempt to make hydrogeological corrections to the measured heat flow values in several boreholes.     

Intercomparison of snow water equivalent observations in the Northern Great Plains

In the Northern Great Plains, melting snow is a primary driver of spring flooding, but limited knowledge of the magnitude and spatial distribution of snow water equivalent (SWE) hampers flood forecasting. Passive microwave remote sensing has the potential to enhance operational river flow forecasting but is not routinely incorporated in operational flood forecasting. We compare satellite passive microwave estimates from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E) to the National Oceanic and Atmospheric Administration Office of Water Prediction (OWP) airborne gamma radiation snow survey and U.S. Army Corps of Engineers (USACE) ground snow survey SWE estimates in the Northern Great Plains from 2002 to 2011. AMSR‐E SWE estimates compare favourably with USACE SWE measurements in the low relief, low vegetation study area (mean difference = ?3.8 mm, root mean squared difference [RMSD] = 34.7 mm), but less so with OWP airborne gamma SWE estimates (mean difference = ?9.5 mm, RMSD = 42.7 mm). An error simulation suggests that up to half of the error in the former comparison is potentially due to subpixel scale SWE variability, limiting the maximum achievable RMSD between ground and satellite SWE to approximately 26–33 mm in the Northern Great Plains. The OWP gamma versus AMSR‐E SWE comparison yields larger error than the point‐scale USACE versus AMSR‐E comparison, despite a larger measurement footprint (5–7 km² vs. a few square centimetres, respectively), suggesting that there are unshared errors between the USACE and OWP gamma SWE data.     

Heat flow maps and deep ground water circulation: Examples from Switzerland

The effect of regional and local ground water circulation systems on the Heat Flow Density (HFD) field is demonstrated by two examples from Switzerland, one near St. Gall in an area at the northern border of the Alps, and the other northwest of Zurich along the eastern end of the Jura mountains. Detailed HFD maps of both areas slow pronounced high heat flow zones which are attributed to discharge of subsurface water which has migrated laterally over several 10 km. Seepage velocities on the order of several mm/yr have been calculated. Geothermal information is not available about the infiltration zones where low HFD values are expected. Geochemical and isotopic analysis of water samples from springs and drillholes indicates the recharge zones and demonstrates the effect of extensive regional systems. These results indicate that in regions with significant topographic relief HFD mapping can be seriously biased if drillholes are positioned in valleys which correspond to discharge areas with relatively high HFD, whereas the low heat flow zones remain undetected.     

Heat as a ground water tracer

Heat carried by ground water serves as a tracer to identify surface water infiltration, flow through fractures, and flow patterns in ground water basins. Temperature measurements can be analyzed for recharge and discharge rates, the effects of surface warming, interchange with surface water, hydraulic conductivity of streambed sediments, and basin-scale permeability. Temperature data are also used in formal solutions of the inverse problem to estimate ground water flow and hydraulic conductivity. The fundamentals of using heat as a ground water tracer were published in the 1960s, but recent work has significantly expanded the application to a variety of hydrogeological settings. In recent work, temperature is used to delineate flows in the hyporheic zone, estimate submarine ground water discharge and depth to the salt-water interface, and in parameter estimation with coupled ground water and heat-flow models. While short reviews of selected work on heat as a ground water tracer can be found in a number of research papers, there is no critical synthesis of the larger body of work found in the hydrogeological literature. The purpose of this review paper is to fill that void and to show that ground water temperature data and associated analytical tools are currently underused and have not yet realized their full potential.     

Securing the future of ground water resources in the Great Lakes basin

Ground water is a vital, but underappreciated, natural resource in the Great Lakes basin. It meets many human needs and contributes significantly to the hydrology of the Great Lakes and the health of ecosystems. This paper provides an overview of ground water in the Great Lakes and the institutional and legal setting that governs the use, protection, diversion, and removal of water from the basin and proposes a citizen-centered vision for management of ground water in the 21st century.     

Fracture control of ground water flow and water chemistry in a rock aquitard

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.     

Role of ground water in geomorphology,geology, and paleoclimate of the Southern High Plains,USA

Study of ground water in the Southern High Plains is central to an understanding of the geomorphology, deposition of economic minerals, and climate change record in the area. Ground water has controlled the course of the Canadian and Pecos rivers that isolated the Southern High Plains from the Great Plains and has contributed significantly to the continuing retreat of the westward escarpment. Evaporative and dissolution processes are responsible for current plateau topography and the development of the signature 20,000 small playa basins and 40 to 50 large saline lake basins in the area. In conjunction with eolian processes, ground water transport controls the mineralogy of commercially valuable mineral deposits and sets up the distribution of fine efflorescent salts that adversely affect water quality. As the water table rises and retreats, lunette and tufa formation provides valuable paleoclimate data for the Southern High Plains. In all these cases, an understanding of ground water processes contributes valuable information to a broad range of geological topics, well beyond traditional interest in water supply and environmental issues.     

Grid resolution study of ground water flow and transport

Three-dimensional grids representing a heterogeneous, ground water system are generated at 10 different resolutions in support of a site-scale flow and transport modeling effort. These grids represent hydrostratigraphy near Yucca Mountain, Nevada, consisting of 18 stratigraphic units with contrasting fluid flow and transport properties. The grid generation method allows the stratigraphy to be modeled by numerical grids of different resolution so that comparison studies can be performed to test for grid quality and determine the resolution required to resolve geologic structure and physical processes such as fluid flow and solute transport. The process of generating numerical grids with appropriate property distributions from geologic conceptual models is automated, thus making the entire process easy to implement with fewer user-induced errors. The series of grids of various resolutions are used to assess the level at which increasing resolution no longer influences the flow and solute transport results. Grid resolution is found to be a critical issue for ground water flow and solute transport. The resolution required in a particular instance is a function of the feature size of the model, the intrinsic properties of materials, the specific physics of the problem, and boundary conditions. The asymptotic nature of results related to flow and transport indicate that for a hydrologic model of the heterogeneous hydrostratigraphy under Yucca Mountain, a horizontal grid spacing of 600 m and vertical grid spacing of 40 m resolve the hydrostratigraphic model with sufficient precision to accurately model the hypothetical flow and solute transport to within 5% of the value that would be obtained with much higher resolution.     

Attenuation relations of ground motion for rock and soil sites in eastern United States

This paper presents the attenuation relations of peak ground acceleration and spectral accelerations for rock and soil sites in the central and eastern United States (CEUS). For the bedrock site condition, 56 pairs of moment magnitude M and epicentral distance R are used to simulate ground motion, and for each pair of M and R, 550 samples of ground motion parameters are generated using a seismological model together with random vibration theory and distribution of extreme values. From the regression analyses of these data, the attenuation relations of ground motion parameters for the bedrock site are established. With the aid of appropriate site coefficients, these attenuation relations are modified for the site categories specified in the 1994 NEHRP Provisions. These attenuation relations are appropriate for the assessment of seismic hazards at far-field rock and soil sites in the CEUS.     

Relating nitrogen sources and aquifer susceptibility to nitrate in shallow ground waters of the United States

Characteristics of nitrogen loading and aquifer susceptibility to contamination were evaluated to determine their influence on contamination of shallow ground water by nitrate. A set of 13 explanatory variables was derived from these characteristics, and variables that have a significant influence were identified using logistic regression (LR). Multivariate LR models based on more than 900 sampled wells predicted the probability of exceeding 4 mg/L of nitrate in ground water. The final LR model consists of the following variables: (1) nitrogen fertilizer loading (p-value = 0.012); (2) percent cropland-pasture (p < 0.001); (3) natural log of population density (p < 0.001); (4) percent well-drained soils (p = 0.002); (5) depth to the seasonally high water table (p = 0.001); and (6) presence or absence of a fracture zone within an aquifer (p = 0.002). Variables 1-3 were compiled within circular, 500 m radius areas surrounding sampled wells, and variables 4-6 were compiled within larger areas representing targeted land use and aquifers of interest. Fitting criteria indicate that the full logistic-regression model is highly significant (p < 0.001), compared with an intercept-only model that contains none of the explanatory variables. A goodness-of-fit test indicates that the model fits the data well, and observed and predicted probabilities of exceeding 4 mg/L nitrate in ground water are strongly correlated (r2 = 0.971). Based on the multivariate LR model, vulnerability of ground water to contamination by nitrate depends not on any single factor but on the combined, simultaneous influence of factors representing nitrogen loading sources and aquifer susceptibility characteristics.     

International borders, ground water flow, and hydroschizophrenia

A substantial body of research has been conducted on transboundary water, transboundary water law, and the mitigation of transboundary water conflict. However, most of this work has focused primarily on surface water supplies. While it is well understood that aquifers cross international boundaries and that the base flow of international river systems is often derived in part from ground water, transboundary ground water and surface water systems are usually managed under different regimes, resulting in what has been described as "hydroschizophrenia." Adding to the problem, the hydrologic relationships between surface and ground water supplies are only known at a reconnaissance level in even the most studied international basins, and thus even basic questions regarding the territorial sovereignty of ground water resources often remain unaddressed or even unasked. Despite the tensions inherent in the international setting, riparian nations have shown tremendous creativity in approaching regional development, often through preventive diplomacy, and the creation of "baskets of benefits," which allow for positive-sum, integrative allocations of joint gains. In contrast to the notion of imminent water wars, the history of hydropolitical relations worldwide has been overwhelmingly cooperative. Limited ground water management in the international arena, coupled with the fact that few states or countries regulate the use of ground water, begs the question: will international borders serve as boundaries for increased "flows" of hydrologic information and communication to maintain strategic aquifers, or will increased competition for shared ground water resources lead to the potential loss of strategic aquifers and "no flows" for both ground water users?     

The role of hand calculations in ground water flow modeling

Most ground water modeling courses focus on the use of computer models and pay little or no attention to traditional analytic solutions to ground water flow problems. This shift in education seems logical. Why waste time to learn about the method of images, or why study analytic solutions to one-dimensional or radial flow problems? Computer models solve much more realistic problems and offer sophisticated graphical output, such as contour plots of potentiometric levels and ground water path lines. However, analytic solutions to elementary ground water flow problems do have something to offer over computer models: insight. For instance, an analytic one-dimensional or radial flow solution, in terms of a mathematical expression, may reveal which parameters affect the success of calibrating a computer model and what to expect when changing parameter values. Similarly, solutions for periodic forcing of one-dimensional or radial flow systems have resulted in a simple decision criterion to assess whether or not transient flow modeling is needed. Basic water balance calculations may offer a useful check on computer-generated capture zones for wellhead protection or aquifer remediation. An easily calculated "characteristic leakage length" provides critical insight into surface water and ground water interactions and flow in multi-aquifer systems. The list goes on. Familiarity with elementary analytic solutions and the capability of performing some simple hand calculations can promote appropriate (computer) modeling techniques, avoids unnecessary complexity, improves reliability, and is likely to save time and money. Training in basic hand calculations should be an important part of the curriculum of ground water modeling courses.     

Market transfers of water for environmental protection in the western United States

Buying and leasing water rights for environmental protection is becoming an important method for protecting river and stream flows in the western United States. The region has experienced an increasing number of market transfers of water to protect water quality and fish and wildlife habitat. From 1990 to 1997, more than US$37 million was spent to lease 2 million acre-feet of water and US$23.8 was spent to purchase 132,000 acre-feet of water for environmental protection. State and federal agencies are responsible for most of these market transfers, but there is increasing activity on the part of private organizations to acquire water for instream needs. This paper examines recent developments of instream flow marketing in the western United States. The paper reviews recent transfers and acquisitions programs and discusses the growing entrepreneurial efforts by private organizations to acquire water rights for instream use. Market information including price and quantity of water traded was collected from market participants. The average purchase and lease prices for the region are US$397 and US$30 per acre-foot, respectively.     

Influence of calibration methodology on ground water flow predictions

We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).     

Using hydraulic head measurements in variable-density ground water flow analyses

The use of hydraulic head measurements in ground water of variable density is considerably more complicated than for the case of constant-density ground water. A theoretical framework for dealing with these complications does exist in the current literature but suffers from a lack of awareness among many hydrogeologists. When corrections for density variations are ignored or not properly taken into account, misinterpretation of both ground water flow direction and magnitude may result. This paper summarizes the existing theoretical framework and provides practical guidelines for the interpretation of head measurements in variable-density ground water systems. It will be argued that, provided that the proper corrections are taken into account, fresh water heads can be used to analyze both horizontal and vertical flow components. To avoid potential confusion, it is recommended that the use of the so-called environmental water head, which was initially introduced to facilitate the analysis of vertical ground water flow, be abandoned in favor of properly computed fresh water head analyses. The presented methodology provides a framework for determining quantitatively when variable-density effects on ground water flow need to be taken into account or can be justifiably neglected. Therefore, we recommend that it should become part of all hydrogeologic analyses in which density effects are suspected to play a role.     

Heat flow anomalies caused by water circulation

The practically important part of geothermal systems belongs to the convective type where the thermal energy is transported by movement of water or steam. Both geothermics and hydrology should be in very close cooperation at the interpretation of convective geothermal anomalies.In the first part of the study the parameters required for the calculation of water- and thermal-balance will be enumerated and their obtainable accuracy will be discussed based mainly on the praxis used in Hungary.In the second part, heat convection problems connected to subterranean water movement will be discussed, divided into three cases which have importance in praxis:

• — regional water-flow systems with great inflow and outflow areas;

• — mountainous — mainly karstic — areas of infiltration with springs at the foot of the mountain;

• — closed convective systems of circulation.

For illustrating the conceptual examples given above, Hungarian case histories with characteristic data will be presented: The Transdanubian Middle Range, Spa of Budapest, Spa of Héviz, the Great Hungarian Plain and the Thermal Anomaly at Tiszakécske.     

Implications of ground water chemistry and flow patterns for earthquake studies

Ground water can facilitate earthquake development and respond physically and chemically to tectonism. Thus, an understanding of ground water circulation in seismically active regions is important for earthquake prediction. To investigate the roles of ground water in the development and prediction of earthquakes, geological and hydrogeological monitoring was conducted in a seismogenic area in the Yanhuai Basin, China. This study used isotopic and hydrogeochemical methods to characterize ground water samples from six hot springs and two cold springs. The hydrochemical data and associated geological and geophysical data were used to identify possible relations between ground water circulation and seismically active structural features. The data for delta18O, deltaD, tritium, and 14C indicate ground water from hot springs is of meteoric origin with subsurface residence times of 50 to 30,320 years. The reservoir temperature and circulation depths of the hot ground water are 57 degrees C to 160 degrees C and 1600 to 5000 m, respectively, as estimated by quartz and chalcedony geothermometers and the geothermal gradient. Various possible origins of noble gases dissolved in the ground water also were evaluated, indicating mantle and deep crust sources consistent with tectonically active segments. A hard intercalated stratum, where small to moderate earthquakes frequently originate, is present between a deep (10 to 20 km), high-electrical conductivity layer and the zone of active ground water circulation. The ground water anomalies are closely related to the structural peculiarity of each monitoring point. These results could have implications for ground water and seismic studies in other seismogenic areas.