ALBERT

Notes: Israel and the Palestinian Authority share the southern Mediterranean coastal aquifer. Long-term overexploitation in the Gaza Strip has resulted in a decreasing water table, accompanied by the degradation of its water quality. Due to high levels of salinity and nitrate and boron pollution, most of the ground water is inadequate for both domestic and agricultural consumption. The rapid rate of population growth in the Gaza Strip and dependence upon ground water as a single water source present a serious challenge for future political stability and economic development. Here, we integrate the results of geochemical studies and numerical modeling to postulate different management scenarios for joint management between Israel and the Palestinian Authority. The chemical and isotopic data show that most of the salinity phenomena in the Gaza Strip are derived from the natural flow of saline ground water from Israel toward the Gaza Strip. As a result, the southern coastal aquifer does not resemble a classic “upstream-downstream” dispute because Israel's pumping of the saline ground water reduces the salinization rates of ground water in the Gaza Strip. Simulation of different pumping scenarios using a monolayer, hydrodynamic, two-dimensional model (MARTHE) confirms the hypothesis that increasing pumping along the Gaza Strip border combined with a moderate reduction of pumping within the Gaza Strip would improve ground water quality within the Gaza Strip. We find that pumping the saline ground water for a source of reverse-osmosis desalination and then supplying the desalinated water to the Gaza Strip should be an essential component of a future joint management strategy between Israel and the Palestinian Authority.

Notes: Recent legislation required regional grassroots water resources planning across the entire state of Texas. The Texas Water Development Board (TWDB), the state's primary water resource planning agency, divided the state into 16 planning regions. Each planning group developed plans to manage both ground water and surface water sources and to meet future demands of various combinations of domestic, agricultural, municipal, and industrial water consumers. This presentation describes the challenges in developing a ground water model for the Llano Estacado Regional Water Planning Group (LERWPG), whose region includes 21 counties in the Southern High Plains of Texas. While surface water is supplied to several cities in this region, the vast majority of the regional water use comes from the High Plains aquifer system, often locally referred to as the Ogallala Aquifer. Over 95% of the ground water demand is for irrigated agriculture. The LERWPG had to predict the impact of future TWDB-projected water demands, as provided by the TWDB, on the aquifer for the period 2000 to 2050. If detrimental impacts were noted, alternative management strategies must be proposed. While much effort was spent on evaluating the current status of the ground water reserves, an appropriate numerical model of the aquifer system was necessary to demonstrate future impacts of the predicted withdrawals as well as the effects of the alternative strategies. The modeling effort was completed in the summer of 2000. This presentation concentrates on the political, scientific, and nontechnical issues in this planning process that complicated the modeling effort. Uncertainties in data, most significantly in distribution and intensity of recharge and withdrawals, significantly impacted the calibration and predictive modeling efforts. Four predictive scenarios, including baseline projections, recurrence of the drought of record, precipitation enhancement, and reduced irrigation demand, were simulated to identify counties at risk of low final ground water storage volume or low levels of satisfied demand by 2050.

Notes: 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.

Notes: Intermontane basins in the Trans-Pecos region of westernmost Texas and northern Chihuahua, Mexico, are target areas for disposal of interstate municipal sludge and have been identified as possible disposal sites for low-level radioactive waste. Understanding ground water movement within and between these basins is needed to assess potential contaminant fate and movement. Four associated basin aquifers are evaluated and classified; the Red Light Draw Aquifer, the Northwest Eagle Flat Aquifer, the Southeast Eagle Flat Aquifer, and the El Cuervo Aquifer. Encompassed on all but one side by mountains and local divides, the Red Light Draw Aquifer has the Rio Grande as an outlet for both surface drainage and ground water discharge. The river juxtaposed against its southern edge, the basin is classified as a topographically open, through-flowing basin. The Northwest Eagle Flat Aquifer is classified as a topographically closed and drained basin because surface drainage is to the interior of the basin and ground water discharge occurs by interbasin ground water flow. Mountains and ground water divides encompass this basin aquifer on all sides; yet, depth to ground water in the interior of the basin is commonly 〉500 feet. Negligible ground water discharge within the basin indicates that ground water discharges from the basin by vertical flow and underflow to a surrounding basin or basins. The most likely mode of discharge is by vertical, cross-formational flow to underlying Permian rocks that are more porous and permeable and subsequent flow along regional flowpaths beneath local ground water divides. The Southeast Eagle Flat Aquifer is classified as a topographically open and drained basin because surface drainage and ground water discharge are to the adjacent Wildhorse Flat area. Opposite the Eagle Flat and Red Light Draw aquifers is the El Cuervo Aquifer of northern Chihuahua, Mexico. The El Cuervo Aquifer has interior drainage to Laguna El Cuervo, which is a phreatic playa that also serves as a focal point of ground water discharge. Our evidence suggests that El Cuervo Aquifer may lose a smaller portion of its discharge by interbasin ground water flow to Indian Hot Springs, near the Rio Grande. Thus, El Cuervo Aquifer is a topographically closed basin that is either partially drained if a component of its ground water discharge reaches Indian Hot Springs or undrained if all its natural ground water discharge is to Laguna El Cuervo.

Notes: Transboundary aquifers are as important a component of global water resource systems as are transboundary rivers; yet, their recognition in international water policy and legislation is very limited. Existing international conventions and agreements barely address aquifers and their resources. To rectify this deficiency, the International Association of Hydrogeologists and UNESCO's International Hydrological Programme have established the Internationally Shared (transboundary) Aquifer Resource Management (ISARM) Programme. This multiagency cooperative program has launched a number of global and regional initiatives. These are designed to delineate and analyze transboundary aquifer systems and to encourage riparian states to work cooperatively toward mutually beneficial and sustainable aquifer development. The agencies participating in ISARM include international and regional organizations (e.g., Organization of American States, United Nations Environment Programme, United Nations Economic Commission for Europe, Food and Agriculture Organization, and South African Development Community). Using outputs of case studies, the ISARM Programme is building scientific, legal, environmental, socioeconomic, and institutional guidelines and recommendations to aid sharing nations in the management of their transboundary aquifers. Since its start in 2000, the program has completed inventories of transboundary aquifers in the Americas and Africa, and several ISARM case studies have commenced.

Notes: The combination of flowmeter and depth-dependent water-quality data was used to evaluate the quantity and source of high-chloride water yielded from different depths to eight production wells in the Pleasant Valley area of southern California. The wells were screened from 117 to 437 m below land surface, and in most cases, flow from the aquifer into the wells was not uniformly distributed throughout the well screen. Wells having as little as 6 m of screen in the overlying upper aquifer system yielded as much as 50% of their water from the upper system during drought periods, while the deeper parts of the well screens yielded 15% or less of the total yield of the wells. Mixing of water within wells during pumping degraded higher-quality water with poorer-quality water from deeper depths, and in some cases with poorer-quality water from the overlying upper aquifer system. Changes in the mixture of water within a well, resulting from changes in the distribution of flow into the well, changed the quality of water from the surface discharge of wells over time. The combination of flowmeter and depth-dependent water quality data yielded information about sources of high-chloride water to wells that was not available on the basis of samples collected from nearby observation wells. Changing well design to eliminate small quantities of poor-quality water from deeper parts of the well may improve the quality of water from some wells without greatly reducing well yield.

Notes: This paper describes reactive transport simulations conducted to assess the impact of mineral fouling on the hydraulic behavior of continuous-wall permeable reactive barriers (PRBs) employing granular zero-valent iron (ZVI) in carbonate-rich alluvial aquifers. The reactive transport model included a geochemical algorithm for simulating corrosion and mineral precipitation reactions that have been observed in ZVI PRBs. Results of simulations show that porosity and hydraulic conductivity of the ZVI decrease over time and that flows are redistributed throughout the PRB in response to fouling of the pore space. Under typical conditions, only subtle changes occur within the first 10 years (i.e., duration of the current field experience record with PRBs), and the most significant changes do not occur until the PRB has operated for at least 30 years. However, changes can occur sooner (or later) if the rate at which mineral-forming ions are delivered to the PRB is higher (or lower) than that expected under typical conditions (i.e., due to higher/lower flow rate or inflowing ground water that has higher/lower ionic strength). When the PRB is more permeable than the aquifer, the median Darcy flux in the PRB does not change appreciably over time because the aquifer controls the rate of flow through the PRB. However, seepage velocities in the PRB increase, and residence times decrease, due to porosity reductions caused by accumulation of minerals in the pore space. When fouling becomes extensive, bypassing and reductions in flow rate in the PRB occur.

Notes: Heat was used as a natural tracer to characterize shallow ground water flow beneath a complex wetland system. Hydrogeologic data were combined with measured vertical temperature profiles to constrain a series of two-dimensional, transient simulations of ground water flow and heat transport using the model code SUTRA (Voss 1990). The measured seasonal temperature signal reached depths of 2.7 m beneath the pond. Hydraulic conductivity was varied in each of the layers in the model in a systematic manual calibration of the two-dimensional model to obtain the best fit to the measured temperature and hydraulic head. Results of a series of representative best-fit simulations represent a range in hydraulic conductivity values that had the best agreement between simulated and observed temperatures and that resulted in simulated pond seepage values within 1 order of magnitude of pond seepage estimated from the water budget. Resulting estimates of ground water discharge to an adjacent agricultural drainage ditch were used to estimate potential dissolved organic carbon (DOC) loads resulting from the restored wetland. Estimated DOC loads ranged from 45 to 1340 g C/(m2 year), which is higher than estimated DOC loads from surface water. In spite of the complexity in characterizing ground water flow in peat soils, using heat as a tracer provided a constrained estimate of subsurface flow from the pond to the agricultural drainage ditch.

Notes: A total of 176 wells in sand-and-gravel glacial aquifers in central Illinois were sampled for arsenic (As) and other chemical parameters. The results were combined with archived and published data from several hundred well samples to determine potential sources of As and the potential geochemical controls on its solubility and mobility. There was considerable spatial variability in the As concentrations. High concentrations were confined to areas smaller than 1 km in diameter. Arsenic and well depth were uncorrelated. Arsenic solubility appeared to be controlled by oxidation-reduction (redox) conditions, especially the presence of organic matter. Geochemical conditions in the aquifers are typically reducing, but only in the most reducing water does As accumulate in solution. In wells in which total organic carbon (TOC) was below 2 mg/L and sulfate (SO42−) was present, As concentrations were low or below the detection limit (0.5 μg/L). Arsenic concentrations 〉10 μg/L were almost always found in wells where TOC was 〉2 mg/L and SO42− was absent or at low concentrations, indicating post–SO42−reducing conditions. Iron (Fe) is common in the aquifer sediments, and Fe oxide reduction appears to be occurring throughout the aquifers. Arsenic is likely released from the solid phase as Fe oxide is reduced.

Notes: A nitrate-reactive porous media layer comprising wood particles with very high hydraulic conductivity (K∼ 1 cm/s) was used to successfully treat nitrate in a shallow sand-and-gravel aquifer in southern Ontario. Nitrate concentrations of 1.3 to 14 mg/L as N in the aquifer were attenuated to 〈0.5 mg/L as N in the reactive layer. Borehole dilution testing indicated that ground water velocities in the reactive layer, although variable, averaged five times higher than in the surrounding aquifer, suggesting that the layer was capturing ground water flow from deeper in the aquifer. The use of high-K reactive media opens up the possibility of installing permeable reactive barriers as horizontal layers in the shallow water table zone that do not necessarily have to penetrate the full depth of a contaminant plume to be effective. Model simulations show that the depth of capture of a high-K layer increases as the layer width in the direction of flow increases. Shallower emplacement could decrease barrier costs at some sites.

Notes: The Blue Ridge Province contains ubiquitous northeast-southwest–trending thrust faults or smaller thrust “slivers” that greatly impact the nature and character of ground water flow in this region. Detailed investigations at a field site in Floyd County, Virginia, indicate that high-permeability zones occur in the brittle crystalline rocks above these thrust faults. Surface and borehole geophysics, aquifer tests, and chlorofluorocarbon and geochemical data reveal that the shallow saprolite aquifer is separated from the deeper fault-zone aquifer by a low–fracture permeability bedrock confining unit, the hydraulic conductivity of which has been estimated to be six orders of magnitude less than the conductivity of the fault zones at the test site. Within the Blue Ridge Province, these fault zones can occur at depths of 300 m or more, can contain a significant amount of storage, and yield significant quantities of water to wells. Furthermore, it is expected that these faults may compartmentalize the deep aquifer system. Recharge to and discharge from the deep aquifer occurs by slow leakage through the confining unit or through localized breach zones that occur where quartz accumulated in high concentrations during metamorphism and later became extensively fractured during episodes of deformation. The results of this investigation stress the importance of thrust faults in this region and suggest that hydrogeologic models for the Blue Ridge Province include these ancient structural features. Faults in crystalline-rock environments may not only influence the hydrology, they may dominate the flow characteristics of a region.

Notes: Uncertainty of hydrogeologic conditions makes it important to consider alternative plausible models in an effort to evaluate the character of a ground water system, maintain parsimony, and make predictions with reasonable definition of their uncertainty. When multiple models are considered, data collection and analysis focus on evaluation of which model(s) is(are) most supported by the data. Generally, more than one model provides a similar acceptable fit to the observations; thus, inference should be made from multiple models. Kullback-Leibler (K-L) information provides a rigorous foundation for model inference that is simple to compute, is easy to interpret, selects parsimonious models, and provides a more realistic measure of precision than evaluation of any one model or evaluation based on other commonly referenced model selection criteria. These alternative criteria strive to identify the true (or quasi-true) model, assume it is represented by one of the models in the set, and given their preference for parsimony regardless of the available number of observations the selected model may be underfit. This is in sharp contrast to the K-L information approach, where models are considered to be approximations to reality, and it is expected that more details of the system will be revealed when more data are available. We provide a simple, computer-generated example to illustrate the procedure for multimodel inference based on K-L information and present arguments, based on statistical underpinnings that have been overlooked with time, that its theoretical basis renders it preferable to other approaches.

Notes: This paper presents a critical review of model-input parameters for transport of on-site wastewater treatment system (OWS) pollutants. Approximately 25% of the U.S. population relies on soil-based OWS for effective treatment and protection of public health and environmental quality. Mathematical models are useful tools for understanding and predicting the transport and fate of wastewater pollutants and for addressing water-budget issues related to wastewater reclamation from site to watershed scales. However, input parameters for models that simulate fate and transport of OWS pollutants are not readily obtained. The purpose of this analysis is to illustrate an objective, statistically supported method for choosing model-input parameters related to nitrogen (N) and phosphorus (P). Data were gathered from existing studies reported in the literature. Cumulative frequency distributions (CFDs) are provided for OWS effluent concentrations of N and P, nitrification and denitrification rates, and linear sorption isotherm constants for P. When CFDs are not presented, ranges and median values are provided. Median values for model-input parameters are as follows: total N concentration (44 mg/L), nitrate-N (0.2 mg/L), ammonium (60 mg/L), phosphate-P (9 mg/L), organic N (14 mg/L), zero-order nitrification rate (264 mg/L/d), first-order nitrification (2.9/d), first-order dentrification (0.025/d), maximum soil capacity for P uptake (237 mg/kg), linear sorption isotherm constant for P (15.1 L/kg), and OWS effluent flow rates (260 L/person/d).

Notes: Boreholes drilled through contaminated zones in fractured rock create the potential for vertical movement of contaminated ground water between fractures. The usual assumption is that purging eliminates cross contamination; however, the results of a field study conducted in a trichloroethylene (TCE) plume in fractured sandstone with a mean matrix porosity of 13% demonstrates that matrix-diffusion effects can be strong and persistent. A deep borehole was drilled to 110 m below ground surface (mbgs) near a shallow bedrock well containing high TCE concentrations. The borehole was cored continuously to collect closely spaced samples of rock for analysis of TCE concentrations. Geophysical logging and flowmetering were conducted in the open borehole, and a removable multilevel monitoring system was installed to provide hydraulic-head and ground water samples from discrete fracture zones. The borehole was later reamed to complete a well screened from 89 to 100 mbgs; persistent TCE concentrations at this depth ranged from 2100 to 33,000 μg/L. Rock-core analyses, combined with the other types of borehole information, show that nearly all of this deep contamination was due to the lingering effects of the downward flow of dissolved TCE from shallower depths during the few days of open-hole conditions that existed prior to installation of the multilevel system. This study demonstrates that transfer of contaminant mass to the matrix by diffusion can cause severe cross contamination effects in sedimentary rocks, but these effects generally are not identified from information normally obtained in fractured-rock investigations, resulting in potential misinterpretation of site conditions.

Notes: 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 μg/L were associated with dissolved oxygen.

Notes: The calculation of the discharge to a constant drawdown well or tunnel in the presence of an infinite linear constant head boundary in an ideal confined aquifer usually relies on the numerical inversion of a Laplace transform solution. Such a solution is used to interpret constant head tests in wells or to roughly estimate ground water inflow into tunnels. In this paper, a simple approximate solution is proposed. Its maximum relative error is on the order of 2% as compared to the exact analytical solution. The approximation is a weighted mean between the early-time and late-time asymptotes.

Notes: Submarine ground water discharge (SGD) rates were measured continuously by automated seepage meters to evaluate the process of ground water discharge to the ocean in the coastal zone of Suruga Bay, Japan. The ratio of terrestrial fresh SGD to total SGD was estimated to be at most 9% by continuous measurements of electrical conductivity of SGD. Semidiurnal changes of SGD due to tidal effects and an inverse relation between SGD and barometric pressure were observed. Power spectrum density analyses of SGD, sea level, and ground water level show that SGD near shore correlated to ground water level changes and SGD offshore correlated to sea level changes. SGD rates near the mouth of the Abe River are smaller than those elsewhere, possibly showing the effect of the river on SGD. The ratio of terrestrial ground water discharge to the total discharge to the ocean was estimated to be 14.7% using a water balance method.

Notes: Fractured rock aquifers are difficult to characterize because of their extremely heterogeneous nature. Developing an understanding of fracture network hydraulic properties in these aquifers is difficult and time consuming, and field testing techniques for determining the location and connectivity of fractures in these aquifers are limited. In the Clare Valley, South Australia, well interference is an important issue for a major viticultural area that uses a fractured aquifer. Five fracture sets exist in the aquifer, all dipping 〉25°. In this setting, we evaluate the ability of steady-state asymmetric dipole-flow tests (ADFTs) to determine the connections between a test well and a set of piezometers. The procedure involves dividing a test well into two chambers using a single packer and pumping fluid from the upper chamber to the lower chamber. By conducting a series of tests at different packer elevations, an “input” signal is generated in fracture zones connected to the test well. By monitoring the “output” response of the hydraulic dipole field at piezometers, the connectivity of the fractures between the test well and piezometers can be determined. Results indicate the test well used in this study is connected in a complex three-dimensional geometry, with drawdown occurring above and below areas of potentiometric buildup. The ADFT method demonstrates that the aquifer evaluated in this study cannot be modeled effectively on the well scale using continuum flow models.

Notes: An approximate analytical solution to the advection-dispersion equation was derived to describe solute transport during spherical-flow conditions in single-well push-pull tests. The spherical-flow case may be applicable to aquifer tests conducted in packed intervals or partially penetrating wells. Using results of two-dimensional numerical simulations, we briefly illustrate the applicability of the derived spherical-flow solution and provide a comparison with its cylindrical-flow counterpart. Good agreement between simulated extraction-phase breakthrough curves and the spherical-flow solution was found when the length of the injection/extraction region was small compared to both aquifer thickness and maximum solute frontal position at the end of the injection phase. On the other hand, discrepancies between simulated breakthrough curves and the spherical-flow solution increased with increasing anisotropy in hydraulic conductivities. Several inherent limitations embedded in its derivation such as assumptions of isotropy and homogeneity warrant the cautious use of the spherical-flow solution.

Notes: This paper describes a technique for applying the transition probability geostatistics method for stochastic simulation to a MODFLOW model. Transition probability geostatistics has some advantages over traditional indicator kriging methods including a simpler and more intuitive framework for interpreting geologic relationships and the ability to simulate juxtapositional tendencies such as fining upward sequences. The indicator arrays generated by the transition probability simulation are converted to layer elevation and thickness arrays for use with the new Hydrogeologic Unit Flow package in MODFLOW 2000. This makes it possible to preserve complex heterogeneity while using reasonably sized grids and/or grids with nonuniform cell thicknesses.

Notes: Improved surface-based geophysical technologies that are commercially available provide a new level of detail that can be used to guide ground water remediation. Surface-based multielectrode resistivity methods and tomographic seismic refraction techniques were used to image to a depth of ~30 m below the surface at the Natural and Accelerated Bioremediation Research Field Research Center. The U.S. Department of Energy (DOE) established the research center on the DOE Oak Ridge Reservation in Oak Ridge, Tennessee, to conduct in situ field-scale studies on bioremediation of metals and radionuclides. Bioremediation studies are being conducted on the saprolite, shale bedrock, and ground water at the site that have been contaminated with nitrate, uranium, technetium, tetrachloroethylene, and other contaminants (U.S. DOE 1997). Geophysical methods were effective in imaging the high-ionic strength plume and in defining the transition zone between saprolite and bedrock zones that appears to have a significant influence on contaminant transport. The geophysical data were used to help select the location and depth of investigation for field research plots. Drilling, borehole geophysics, and ground water sampling were used to verify the surface geophysical studies.

Notes: Calibration of ground water transport models is often performed using results of field tracer experiments. However, little attention is usually paid to the influence, on resulting breakthrough curves, of injection conditions and well-aquifer interactions, more particularly of the influence of the possible trapping of the tracer in the injection wellbore. Recently, a new mathematical and numerical approach has been developed to model injection conditions and well-aquifer interactions in a very accurate way. Using an analytical solution derived from this model, a detailed analysis is made of the evolution of the tracer input function in the aquifer. By varying injection conditions from one simulation to another, synthetic breakthrough curves are generated with the SUFT3D ground water flow and transport finite-element simulator. These tests show clearly that the shape of the breakthrough curves can be dramatically affected by injection conditions. Using generated breakthrough curves as “actual” field results, a calibration of hydrodispersive parameters is performed, neglecting the influence of injection conditions. This shows that neglecting the influence of actual injection conditions can lead to (1) errors on fitted parameters and (2) misleading identification of the active transport processes. Conclusions and guidelines are drawn in terms of proposed methodologies for better controlling the tracer injection in the field, in order to minimize risk of misinterpretation of results.

Notes: We describe the continuous time random walk (CTRW) MATLAB toolbox, a collection of MATLAB scripts and functions that compute breakthrough curves (BTCs) and one-dimensional/two-dimensional (1D/2D) resident concentration profiles for passive tracer dispersion. The transport model is based on the CTRW theory. CTRW includes as special cases the classical Fickian dispersion based advection-dispersion equation, multirate and mobile-immobile models, and the fractional-in-time derivative transport equation. Several models for treating the memory effects responsible for the anomalous character of dispersion have been implemented in the CTRW toolbox. In the current version of the toolbox, it is possible to solve explicitly for the forward problem (concentration prediction) in 1D and 2D and for the inverse problem (parameter identification from experimental BTC data) in 1D. Future extensions will include explicit treatment of sorbing tracers, simple subroutines for treating radial flow from wells, introduction of arbitrary initial conditions, treatment of heterogeneous domains by use of the Fokker-Planck with Memory equation, and treatment of transport in multidimensional systems.

Notes: We present a novel hybrid algorithm, integrating a genetic algorithm (GA) and constrained differential dynamic programming (CDDP), to achieve remediation planning for an unconfined aquifer. The objective function includes both fixed and dynamic operation costs. GA determines the primary structure of the proposed algorithm, and a chromosome therein implemented by a series of binary digits represents a potential network design. The time-varying optimal operation cost associated with the network design is computed by the CDDP, in which is embedded a numerical transport model. Several computational approaches, including a chromosome bookkeeping procedure, are implemented to alleviate computational loading. Additionally, case studies that involve fixed and time-varying operating costs for confined and unconfined aquifers, respectively, are discussed to elucidate the effectiveness of the proposed algorithm. Simulation results indicate that the fixed costs markedly affect the optimal design, including the number and locations of the wells. Furthermore, the solution obtained using the confined approximation for an unconfined aquifer may be infeasible, as determined by an unconfined simulation.

Notes: The water table in unconfined aquifers is often believed to be a subdued replica of the topography or land surface. However, this assumption has not been widely tested and in some cases has been found to be in error. An analysis of ground water rise in regional unconfined aquifers, using both a two-dimensional boundary element model and a one-dimensional Dupuit-Forchheimer model, reveals the conditions under which the water table does or does not closely follow the topography. A simple decision criterion is presented to estimate in advance under which conditions the water table is expected to be largely unrelated to the topography and under which conditions the topography controls the position of the water table.

Notes: Estimation of available ground water is a basic aspect of ground water management. Mathematical modeling is one of the methods that can be effectively used to obtain such estimates. A numerical model was used to calculate available ground water in the Zohor depression—an aquifer transcending national boundaries between the Slovak Republic and Austria. The aquifer, formed by Quaternary sediments overlying Neogene sequences, is composed of various clays interbedded with layers of sand, gravel, sandstones, and conglomerates. The AQUA computer model package was used to simulate flow in the aquifer. For model compilation, the following data were used: (1) effective precipitation; (2) surface water levels in surface water gauging profiles; and (3) withdrawal amounts. Hydraulic parameters of the aquifer were estimated based on information from 86 wells located in the area. The model was verified using data on ground water levels from a monitoring network. The simulation of the aquifer system permitted the estimation of the available ground water in the study area, showing that an additional 587 L/s can be abstracted. Ground water inflows to the Morava River, which flows through the region, range from 745 to 3100 L/s.

Notes: Ground water in deep confined aquifers is one of the major water resources for agricultural, industrial, and domestic uses in the North China Plain. Detailed information on ground water age and recharge is vital for the proper management of these water resources, and to this end, we used carbon 14 of dissolved inorganic carbon and tritium in water to measure the age and determine the recharge areas of ground water in the North China Plain. These isotopic data suggest that most ground water in the piedmont part of the North China Plain is 〈40 years old and is recharged locally. In contrast, ground water in the central and littoral portions of the North China Plain is 10,000 to 25,000 years old. The δ18O (δD) values of this ground water are 1.7‰ (11‰) less than that in the piedmont plain ground water and possibly reflect water recharged during a cooler climate during the last glaciation. The temperature of this recharge, based on δ18O values, ranges from 3.7°C to 8.4°C, compared to 12°C to 13°C of modern recharge water. The isotopic data set combined indicates that ground water in the central and littoral part of the North China Plain is being mined under non–steady state conditions.

Notes: Recharge events that deliver electron acceptors such as O2, NO3, SO4, and Fe3+ to anaerobic, contaminated aquifers are likely important for natural attenuation processes. However, the specific influence of recharge on (bio)geochemical processes in ground water systems is not well understood. The impact of a moderate-sized recharge event on ground water chemistry was evaluated at a shallow, sandy aquifer contaminated with waste fuels and chlorinated solvents. Multivariate statistical analyses coupled with three-dimensional visualization were used to analyze ground water chemistry data (including redox indicators, major ions, and physical parameters) to reveal associations between chemical parameters and to infer processes within the ground water plume. Factor analysis indicated that dominant chemical associations and their interpreted processes (anaerobic and aerobic microbial processes, mineral precipitation/dissolution, and temperature effects) did not change significantly after the spring recharge event of 2000. However, the relative importance of each of these processes within the plume changed. After the recharge event, the overall importance of aerobic processes increased from the fourth to the second most important factor, representing the variability within the data set. The anaerobic signatures became more complex, suggesting that zones with multiple terminal electron–accepting processes (TEAPs) likely occur in the same water mass. Three-dimensional visualization of well clusters showed that water samples with similar chemical associations occurred in distinct water masses within the aquifer. Water mass distinctions were not based on dominant TEAPs, suggesting that the recharge effects on TEAPs occurred primarily at the interface between infiltrating recharge water and the aquifer.

Notes: The validity of a previously proposed but untested modification to equations for flow to a horizontal well is assessed using a specially developed finite-difference model. This modification extends confined flow equations to allow the head in the well and the saturated depth at the well to be estimated in unconfined conditions. The study is limited to the case of two-dimensional flow with no flow in the direction parallel to the line of the well. The results show that the modified equations for both a finite unconfined aquifer and, by inference, an infinite unconfined aquifer are adequately accurate for practical application.

Notes: Artificial neural networks (ANNs) were developed for accurately predicting potentiometric surface elevations (monitoring well water level elevations) in a semiconfined glacial sand and gravel aquifer under variable state, pumping extraction, and climate conditions. ANNs “learn” the system behavior of interest by processing representative data patterns through a mathematical structure analogous to the human brain. In this study, the ANNs used the initial water level measurements, production well extractions, and climate conditions to predict the final water level elevations 30 d into the future at two monitoring wells. A sensitivity analysis was conducted with the ANNs that quantified the importance of the various input predictor variables on final water level elevations. Unlike traditional physical-based models, ANNs do not require explicit characterization of the physical system and related physical data. Accordingly, ANN predictions were made on the basis of more easily quantifiable, measured variables, rather than physical model input parameters and conditions. This study demonstrates that ANNs can provide both excellent prediction capability and valuable sensitivity analyses, which can result in more appropriate ground water management strategies.

Notes: The standard formulation of a block-centered finite-difference model, such as MODFLOW, uses the center of the cell as the location of a cell node. Simulations of a dynamic water table can be improved if the node of a cell containing the water table is located at the water table rather than at the center of the cell. The LPF package of MOD-FLOW-2000 was changed to position a cell's node at the water table in convertible cells with a water table. Improved accuracy in the upper regions of an unconfined aquifer is demonstrated for pumping from a partially penetrating well. The change introduces a nonlinearity into the solution of the flow equations that results in slightly slower convergence of the flow solution, 7% slower in the presented demonstration. Accuracy of simulations is improved where vertical flow is dominated by a moving water table, but not when a large water table gradient dominates over the water table movement.

Notes: Ground water quality networks for monitoring phreatic drinking water wellfields are generally established for two main purposes: (1) the short-term safeguarding of public water supply and (2) signaling and predicting future quality changes in the extracted ground water. Six monitoring configurations with different well locations and different screen depths and lengths were evaluated using a numerical model of the 3D ground water flow toward a partially penetrating pumping well in a phreatic aquifer. Travel times and breakthrough curves for observation and pumping wells were used to judge the effectiveness of different design configurations for three monitoring objectives: (1) early warning; (2) prediction of future quality changes; and (3) evaluation of protection measures inside a protection zone. Effectiveness was tested for scenarios with advective transport, first-order degradation, and linear sorption. It is shown that the location and especially the depth of the observation wells should be carefully chosen, taking into account the residence time from the surface to the observation well, the residual transit times to the extraction well, and the transformation and retardation rates. Shallow monitoring was most functional for a variety of objectives and conditions. The larger the degradation rates or retardation, the shallower should the monitoring be for effective early warning and prediction of future ground water quality. The general approach followed in the current study is applicable for many geohydrological situations, tuning specific monitoring objectives with residence times and residual transit times obtained from a site-specific ground water flow model.

Notes: In the upper Midwest, USA, elevated arsenic concentrations in public drinking water systems are associated with the lateral extent of northwest provenance late Wisconsin-aged drift. Twelve percent of public water systems located within the footprint of this drift (212 of 1764) exceed 10 μg/L arsenic, which is the U.S. EPA's drinking water standard. Outside of the footprint, only 2.4% of public water systems (52 of 2182) exceed 10 μg/L arsenic. Both glacial drift aquifers and shallow bedrock aquifers overlain by northwest provenance late Wisconsin-aged sediment are affected by arsenic contamination. Evidence suggests that the distinct physical characteristics of northwest provenance late Wisconsin-aged drift—its fine-grained matrix and entrained organic carbon that fosters biological activity—cause the geochemical conditions necessary to mobilize arsenic via reductive mechanisms such as reductive desorption and reductive dissolution of metal oxides. This study highlights an important and often unrecognized phenomenon: high-arsenic sediment is not necessary to cause arsenic-impacted ground water—when “impacted” is now defined as 〉10 μg/L. This analysis also demonstrates the scientific and economic value of using existing large but imperfect statewide data sets to observe and characterize regional-scale environmental problems.

Notes: 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.

Notes: Using heat as a tracer allows for estimation of ground water recharge rates based on subsurface temperature measurements. While possible in theory, it may be difficult in practice to discriminate the effects of climate from the effects of ground water advection. This study uses synthetic simulations to determine the influence of variability of ground surface temperature (GST) on the ability to estimate vertical specific discharge from temperature profiles. Results suggest that in cases where temperature measurements are sufficiently deep and specific discharge is sufficiently high, estimates of specific discharges will be reasonably accurate. Increasing the number of times temperatures are measured, or producing models that incorporate variations in GST, will increase the reliability of any studies using temperatures to estimate specific discharge. Furthermore, inversions of temperature measurements should be combined with other methods of estimating recharge rates to improve the reliability of recharge estimates.

Notes: One of the challenges in characterizing fractured-rock aquifers is determining whether the equivalent porous medium approximation is valid at the problem scale. Detailed hydrogeologic characterization completed at a small study site in a densely fractured dolomite has yielded an extensive data set that was used to evaluate the utility of the continuum and discrete-fracture approaches to aquifer characterization. There are two near-vertical sets of fractures at the site; near-horizontal bedding-plane partings constitute a third fracture set. Eighteen boreholes, including five coreholes, were drilled to a depth of ∼10.6 m. Borehole geophysical logs revealed several laterally extensive horizontal fractures and dissolution zones. Flowmeter and short-interval packer testing identified which of these features were hydraulically important. A monitoring system, consisting of short-interval piezometers and multilevel samplers, was designed to monitor four horizontal fractures and two dissolution zones. The resulting network consisted of 〉70 sampling points and allowed detailed monitoring of head distributions in three dimensions.Comparison of distributions of hydraulic head and hydraulic conductivity determined by these two approaches suggests that even in a densely fractured-carbonate aquifer, a characterization approach using traditional long-interval monitoring wells is inadequate to characterize ground water movement for the purposes of regulatory monitoring or site remediation. In addition, traditional multiwell pumping tests yield an average or bulk hydraulic conductivity that is not adequate for predicting rapid ground water travel times through the fracture network, and the pumping test response does not appear to be an adequate tool for assessing whether the porous medium approximation is valid.

Notes: Theis (1935) derived an exact solution for the residual drawdown in a well after the cessation of a pumping test by summing two drawdowns: one (s1), caused by imaginary continuation of the original pumping and the other (s2), due to an imaginary injection at the same constant rate. We approximated the Theis solution to obtain a simple linear relation for determining the transmissivity and storage coefficient from recovery data. Unlike other existing straight-line fitting methods, in our method, we applied different approximations to the well functions in the solutions of s1 and s2. We used the well-known Cooper-Jacob approximation for s1, truncating the expansion of the well function in s2 to its first three terms. For the same level of truncation errors, while the Cooper-Jacob approximation requires the argument u1≤ 0.01, the second approximation requires the argument u2≤ 0.2, thus dramatically improving the utility of short-term recovery data. Application of our method requires only recovery data from a single observation well and no knowledge of the drawdown at the moment of pumping cessation.

Notes: Consumption of aquifer Fe(III) during biodegradation of ground water contaminants may result in expansion of a contaminant plume, changing the outlook for monitored natural attenuation. Data from two research sites contaminated with petroleum hydrocarbons show that toluene and xylenes degrade under methanogenic conditions, but the benzene and ethylbenzene plumes grow as aquifer Fe(III) supplies are depleted. By considering a one-dimensional reaction front in a constant unidirectional flow field, it is possible to derive a simple expression for the growth rate of a benzene plume. The method balances the mass flux of benzene with the Fe(III) content of the aquifer, assuming that the biodegradation reaction is instantaneous. The resulting expression shows that the benzene front migration is retarded relative to the ground water velocity by a factor that depends on the concentrations of hydrocarbon and bioavailable Fe(III). The method provides good agreement with benzene plumes at a crude oil study site in Minnesota and a gasoline site in South Carolina. Compared to the South Carolina site, the Minnesota site has 25% higher benzene flux but eight times the Fe(III), leading to about one-sixth the expansion rate. Although it was developed for benzene, toluene, ethylbenzene, and xylenes, the growth-rate estimation method may have applications to contaminant plumes from other persistent contaminant sources.

Notes: This paper describes a methodology for resolving transboundary water disputes that arise when people/states/nations sharing a resource that crosses legal/political jurisdictions disagree about the use of the resource. Laws and treaties written in an attempt to settle disputes are frequently neither enforced nor effective, and disagreements continue. Crises, arising through resource overuse or shortages, worsen the conflict and typically result in further discord, lawsuits, depletion of the resource, and even open-armed hostility. Many water management experts call for either private/market-based or state/command-and-control resource management systems, but these eventually break down during crisis. The crises therefore necessitate the adoption of a more effective institutional arrangement to address and resolve present and future problems. A better alternative to management by private or state entities and the resolution of conflicts by the mere application of law is a cooperative approach. The Rowland-Ostrom Framework, introduced in this paper, incorporates Ostrom's eight design principles for sustainable common pool resource management within the context of crisis that involves an urgent threat to the quantity or quality of a resource such as water, as described by the author. This paper demonstrates that although established 15 years ago, Ostrom's design principles remain applicable today for effective, sustainable transboundary water management, and the Rowland-Ostrom Framework is a model for the equitable use of shared water resources throughout the world.

Notes: 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?

Notes: Access to fresh water is one of the major issues of northern and sub-Saharan Africa. The majority of the fresh water used for drinking and irrigation is obtained from large ground water basins where there is minor contemporary recharge and the aquifers cross national borders. These aquifers include the Nubian Aquifer System shared by Chad, Egypt, Libya, and Sudan; the Iullemeden Aquifer System, extending over Niger, Nigeria, Mali, Benin, and Algeria; and the Northwest Sahara Aquifer System shared by Algeria, Libya, and Tunisia. These resources are subject to increased exploitation and may be severely stressed if not managed properly as witnessed already by declining water levels. In order to make appropriate decisions for the sustainable management of these shared water resources, planners and managers in different countries need an improved knowledge base of hydrological information. Three technical cooperation projects related to aquifer systems will be implemented by the International Atomic Energy Agency, in collaboration with the United Nations Educational, Scientific and Cultural Organization and United Nations Development Programme/Global Environmental Facility. These projects focus on isotope hydrology studies to better quantify ground water recharge and dynamics. The multiple isotope approach combining commonly used isotopes 18O and 2H together with more recently developed techniques (chlorofluorocarbons, 36Cl, noble gases) will be applied to improve the conceptual model to study stratification and ground water flows. Moreover, the isotopes will be an important indicator of changes in the aquifer due to water abstraction, and therefore they will assist in the effort to establish a sustainable ground water management.

Notes: Variations in the concentrations and isotopic compositions (δ13CDIC) of dissolved inorganic carbon (DIC) reflect contamination and biogeochemical cycling of the carbon in ground water. In order to understand contamination and biogeochemical cycling of DIC, we carried out research on the geochemistry of ground water of Guiyang, the capital city of Guizhou Province, China. Results show that ground water is mainly characterized by SO4·HCO3-Ca·Mg and HCO3-Ca·Mg chemical compositions. The hydrochemical characteristics of these types of water are mainly controlled by lithology of the aquifers. HCO3− is the dominant species of DIC in ground water and has lower concentrations and more negative values of δ13CDIC in the high-flow (summer monsoon) season, as compared to the low-flow season. This indicates that DIC is relatively enriched in carbon of biological origin in the high-flow season as compared to the low-flow season and that biological activities are the predominant control on shifts of stable carbon isotope values. The evidence that the δ13CDIC values of ground water decrease with increasing concentrations of anthropogenic species shows that the carbon isotopic composition of DIC can be a useful tracer of contamination, in addition to biogeochemical cycling of inorganic carbon in ground water. Results from this study show that ground water is impacted by significant levels of contamination from human activities, especially in the urban areas, as well as the northeast and west suburbs, in Guiyang city, southwest China.

Notes: 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 δ18O, δD, 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°C to 160°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.

Notes: An exact, closed-form analytical solution is developed for calculating ground water transit times within Dupuit-type flow systems. The solution applies to steady-state, saturated flow through an unconfined, horizontal aquifer recharged by surface infiltration and discharging to a downgradient fixed-head boundary. The upgradient boundary can represent, using the same equation, a no-flow boundary or a fixed head. The approach is unique for calculating travel times because it makes no a priori assumptions regarding the limit of the water table rise with respect to the minimum saturated aquifer thickness. The computed travel times are verified against a numerical model, and examples are provided, which show that the predicted travel times can be on the order of nine times longer relative to existing analytical solutions.

Notes: Ground water plays an important role in water supply and the ecology of arid to semiarid areas such as Northwest China, where the landscape is fragile due to frequent drought in the past few decades. This paper discusses the role of ground water in these ecosystems, including the effect of condensation water and water table depth on the growth of plants and degree of soil salinity. The paper also discusses the controlling process for land desertification and soil salinization in Northwest China. Water table depth is a key factor controlling the water balance, ground water flow, and salt transport in the vadose zone. The suitable water table depth for vegetation growth, which can prevent land desertification and soil salinization, is within a range of 2 to 4 m; the optimal depth is ∼3 m. As examples, changes in ecosystems owing to water resources development in Tarim and Manas basins, Xinjiang, China, are discussed.

Notes: Hydrogeology is a field closely related to the needs of society. Many problems of current national and local interest require predictions of hydrogeological system behavior, and in a number of important cases, the period of prediction is tens to hundreds of thousands of years. It is argued that the demand for such long-term hydrogeological predictions casts a new light on the future needs of hydrogeological research. Key scientific issues are no longer concerned only with simple processes or narrowly focused modeling or testing methods but also with assessment of prediction uncertainties and confidence, couplings among multiple physicochemical processes occurring simultaneously at a site, and the interplay between site characterization and predictive modeling. These considerations also have significant implications for hydrogeological education. With this view, it is asserted that hydrogeological directions and education need to be reexamined and possibly refocused to address specific needs for long-term predictions.

Notes: Historically, specific capacity information has been used to calculate aquifer transmissivity when pumping test data are unavailable. This paper presents a simple computer program written in the MATLAB programming language that estimates transmissivity from specific capacity data while correcting for aquifer partial penetration and well efficiency. The program graphically plots transmissivity as a function of these factors so that the user can visually estimate their relative importance in a particular application. The program is compatible with any computer operating system running MATLAB, including Windows, Macintosh OS, Linux, and Unix. Two simple examples illustrate program usage.

Notes: Alternative fractional models of contaminant transport lead to a new travel time formula for arbitrary concentration levels. For an evolving contaminant plume in a highly heterogeneous aquifer, the new formula predicts much earlier arrival at low concentrations. Travel times of contaminant fronts and plumes are often obtained from Darcy's law calculations using estimates of average pore velocities. These estimates only provide information about the travel time of the average concentration (or peak, for contaminant pulses). Recently, it has been shown that finding the travel times of arbitrary concentration levels is a straightforward process, and equations were developed for other portions of the breakthrough curve for a nonreactive contaminant. In this paper, we generalize those equations to include alternative fractional models of contaminant transport.

Notes: Regulatory agencies are becoming increasingly interested in using young–ground water dating techniques, such as the 3H/3He method, in assessing the susceptibility of public supply wells (PSWs) to contamination. However, recent studies emphasize that ground water samples of mixed age may be the norm, particularly from long-screened PSWs, and tracer-based “apparent” ages can differ substantially from actual mean ages for mixed-age samples. We present age and contaminant data from PSWs in Salt Lake Valley, Utah, that demonstrate the utility of 3H and 3He measurements in evaluating well susceptibility, despite potential age mixing. Initial 3H concentrations (measured 3H + measured tritiogenic 3He) are compared to those expected based on the apparent 3H/3He age and the local precipitation 3H record. This comparison is used to determine the amount of modern water (recharged after ∼1950) vs. prebomb water (recharged before ∼1950) samples might contain. Concentrations of common contaminants were also measured using detection limits generally lower than those used for regulatory purposes. A clear correlation exists between the potential magnitude of the modern water fraction and both the occurrence and concentration of contaminants. For samples containing dominantly modern water based on their initial 3H concentrations, potential discrepancies between apparent 3H/3He ages and mean ages are explored using synthetic samples that are random mixtures of different modern waters. Apparent ages can exceed mean ages by up to 13 years for these samples, with an exponential age distribution resulting in the greatest discrepancies.

Notes: A common assumption of ground water models formulated using a block-centered finite-difference method is that a well is located at the center of a cell regardless of its actual location. Due to this assumption, errors are introduced in the spatial distribution of simulated heads. This paper presents an alternative approach for assigning the pumping rates of wells that are located off cell centers. This approach consists of assigning the pumping rate not only to the cell in which the well is located but also to adjacent cells, taking into account the length of the well screen, the hydraulic conductivity, and the distance from the well to the center of its cell. The advantage of this alternative approach over the conventional one is illustrated with a test problem of a synthetic aquifer. Statistical measures of error indicate a much better model fit when pumping rates of wells are distributed over several cells.

Notes: When fugitive methane migrates upward along boreholes of oil and gas wells, it may migrate into shallow ground water or pass through overlying soil to the atmosphere. Prior to this study, there was little information on the fate of fugitive methane that migrates into ground water. In a field study near Lloydminster, Alberta, Canada, we found hydrogeochemical evidence that fugitive methane from an oil well migrated into a shallow aquifer but has been attenuated by dissimilatory bacterial sulfate reduction at low temperature (∼5°C) under anaerobic conditions. Evidence includes spatial and temporal trends in concentrations of methane and sulfate in ground water and associated trends in concentrations of bicarbonate and sulfide. Within 10 m of the oil well, sulfate concentrations were low, and sulfate was enriched in both 34S and 18O. Sulfate concentrations had a strong positive correlation with δ13C values of bicarbonate, and sulfide was depleted in 34S compared to sulfate. These data indicate that bacterial sulfate reduction occurred near the production well. Near the oil well, elevated concentrations of bicarbonate were observed, and the bicarbonate was depleted in 13C. Modeling indicates that the main source of this excess 13C-depleted bicarbonate is oxidized methane. In concert with the sulfate concentration and isotope data, these results support an interpretation that in situ bacterial oxidation of methane has occurred, linked to bacterial sulfate reduction. Bacterial sulfate reduction may play a major role in bioattenuation of fugitive natural gas in ground water in western Canada.

Notes: Fracture trends (defined as kilometer-scale linear features interpolated between field observations of fractures along their strikes) often have a dominant orientation. Finding a correlation between this orientation and hydraulic data could shed light on their hydraulic influence. A significant correlation between head residuals from first-order regional drift and the orientation of 2- to 4-km-long fracture trends was found in a study site in the Negev, Israel, using the semivariogram cloud analysis. Correlation of head residuals rather than the head itself implies that the orientation of the fracture trends controls the anisotropy and heterogeneity at this scale, mainly because the fracture trends define the orientation of blocks, which differ in their hydraulic properties. Preferential transmissive pathways are probably shorter than the full extent of the fracture trends, causing a relatively high head difference along the trends on the 2- to 4-km scale. Fracture trend density and additional data from short-range hydraulic tests helped characterize two blocks separated by a fault zone. The identification of hydraulic features on a kilometer scale is necessary for better modeling of regional ground water flow and transport. Hydraulic tests at this scale are not feasible, thereby rendering combined analyses of head and structural data, such as the one presented here, essential.

Notes: Ground water nitrate concentrations on Nantucket Island, Massachusetts, were analyzed to assess the effects of land use on ground water quality. Exploratory data analysis was applied to historic ground water nitrate concentrations to determine spatial and temporal trends. Maximum likelihood Tobit and logistic regression analyses of explanatory variables that characterize land use within a 1000-foot radius of each well were used to develop predictive equations for nitrate concentration at 69 wells. The results demonstrate that historic nitrate concentrations downgradient from agricultural land are significantly higher than nitrate concentrations elsewhere. Tobit regression results demonstrate that the number of septic tanks and the percentages of forest, undeveloped, and high-density residential land within a 1000-foot radius of a well are reliable predictors of nitrate concentration in ground water. Similarly, logistic regression revealed that the percentages of forest, undeveloped, and low-density residential land are good indicators of ground water nitrate concentration 〉2 mg/L. The methodology and results outlined here provide a useful tool for land managers in communities with shallow water tables overlain with highly permeable materials to evaluate potential effects of development on ground water quality.

Notes: As part of its Wellhead Protection Program, the U.S. EPA mandates the delineation of “time-of-travel capture zones” as the basis for the definition of wellhead protection zones surrounding drinking water production wells. Depending on circumstances the capture zones may be determined using methods that range from simply drawing a circle around the well to sophisticated ground water flow and transport modeling. The simpler methods are attractive when faced with the delineation of hundreds or thousands of capture zones for small public drinking water supply wells. On the other hand, a circular capture zone may not be adequate in the presence of an ambient ground water flow regime. A dimensionless time-of-travel parameter T˜ is used to determine when calculated fixed-radius capture zones can be used for drinking water production wells. The parameter incorporates aquifer properties, the magnitude of the ambient ground water flow field, and the travel time criterion for the time-of-travel capture zone. In the absence of interfering flow features, three different simple capture zones can be used depending on the value of T˜. A modified calculated fixed-radius capture zone proves protective when T˜ 〈 0.1, while a more elongated capture zone must be used when T˜ 〉1. For values of T˜ between 0.1 and 1, a circular capture zone can be used that is eccentric with respect to the well. Finally, calculating T˜ allows for a quick assessment of the validity of circular capture zones without redoing the delineation with a computer model.

Notes: Ground water age is a fundamental, yet complex, concept in ground water hydrology. Discrepancies between results obtained through different modeling approaches for ground water age calculation have been reported, in particular, between ground water ages modeled by advection and direct simulation of ground water ages (e.g., age-mass approach), which includes effects of advection and dispersion. Here, through a series of two-dimensional (2D) simulations, the impact of water mixing through advection and dispersion on modeled 14C and directly simulated ground water ages is assessed. Impact of dispersion on modeled ages is systematically stronger in areas where water velocities are smaller and far more pronounced on 14C ages. This effect is also observed in one-dimensional models. 2D simulations show that longitudinal dispersion generally acts as a “source” of 14C, while vertical dispersion acts as a “sink,” leading to apparent younger or older modeled 14C ages as compared to advective and directly simulated ground water ages. The presence of permeable and impermeable faults provides an equally important source for discrepancies, leading to major differences in modeled ages among the three methods considered. Overall, our results show that a 14C modeling approach using a solute transport model for calculating ground water age appears to be more reliable in ground water systems without faults and where water velocities are relatively high than in systems that are relatively more heterogeneous and those where faults are present. Among the three modeling approaches considered here, direct simulation of ground water age seems to yield the most consistent results in complex, heterogeneous ground water flow systems, giving a vertical age structure consistent with ages expected from consideration of the flow system.

Notes: In this paper we present a hydrologic application of a new statistical learning methodology called support vector machines (SVMs). SVMs are based on minimization of a bound on the generalized error (risk) model, rather than just the mean square error over a training set. Due to Mercer's conditions on the kernels, the corresponding optimization problems are convex and hence have no local minima. In this paper, SVMs are illustratively used to reproduce the behavior of Monte Carlo–based flow and transport models that are in turn used in the design of a ground water contamination detection monitoring system. The traditional approach, which is based on solving transient transport equations for each new configuration of a conductivity field, is too time consuming in practical applications. Thus, there is a need to capture the behavior of the transport phenomenon in random media in a relatively simple manner. The objective of the exercise is to maximize the probability of detecting contaminants that exceed some regulatory standard before they reach a compliance boundary, while minimizing cost (i.e., number of monitoring wells). Application of the method at a generic site showed a rather promising performance, which leads us to believe that SVMs could be successfully employed in other areas of hydrology. The SVM was trained using 510 monitoring configuration samples generated from 200 Monte Carlo flow and transport realizations. The best configurations of well networks selected by the SVM were identical with the ones obtained from the physical model, but the reliabilities provided by the respective networks differ slightly.

Notes: Longitudinal dispersivity (α) data were compiled from 109 different authors for different types of geological media. The data were subdivided into different subsets. Dispersivity values for consolidated media were subdivided as basalts, granites, sandstones, and carbonate rocks, while unconsolidated sediments were subdivided into three reliability classes. The data sets provided here may provide ground water practitioners a preliminary guide to estimate dispersivity values at various scales and to guide and verify theories on scaling behavior. Based on the data set presented here, the relationship that empirically best described the dispersivity data in regard to scale of measurement was in the form of a power law. The scaling exponent for consolidated and unconsolidated geological media varied between 0.40 and 0.92, and 0.44 and 0.94, respectively. Higher reliability subsets of data for the unconsolidated sediments and more frequently tested rock formations indicate that the scaling exponent is at the lower end of the observed range, close to 0.5. No significant difference in scaling exponent was found among different media, and no clear evidence exists for the presence of an upper bound or asymptotic behavior on the relationship for any of the analyzed media.

Notes: A methodology using ordinal logistic regression is proposed to predict the probability of occurrence of heavy metals in ground water. The predicted probabilities are defined with reference to the background concentration and the maximum contaminant level. The model is able to predict the occurrence due to different influencing variables such as the land use, soil hydrologic group (SHG), and surface elevation. The methodology was applied to the Sumas-Blaine Aquifer located in Washington State to predict the occurrence of five heavy metals. The influencing variables considered were (1) SHG; (2) land use; (3) elevation; (4) clay content; (5) hydraulic conductivity; and (6) well depth. The predicted probabilities were in agreement with the observed probabilities under existing conditions. The results showed that aquifer vulnerability to each heavy metal was related to different sets of influencing variables. However, all heavy metals had a strong influence from land use and SHG. The model results also provided good insight into the influence of various hydrogeochemical factors and land uses on the presence of each heavy metal. A simple economic analysis was proposed and demonstrated to evaluate the cost effects of changing the land use on heavy metal occurrence.

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

Notes: Many people in sub-Saharan Africa have to rely on meager water resources within mudstones for their only water supply. Although mudstones have been extensively researched for their low permeability behavior, little research has been undertaken to examine their ability to provide sustainable water supplies. To investigate the factors controlling the occurrence of usable ground water in mudstone environments, an area of Cretaceous mudstones in southeastern Nigeria was studied over a 3 yr period. Transmissivity (T) variations in a range of mudstone environments were studied. The investigations demonstrate that within the top 40 m of mudstones, transmissivity can be sufficient to develop village water supplies (T 〉 1 m2/d). Transmissivity is controlled by two factors: low-grade metamorphism and the presence of other, subordinate, lithologies within the mudstones. Largely unaltered mudstones (early diagenetic zone), comprising mainly smectite clays, are mostly unfractured and have a low T of 〈 0.1 m2/d. Mudstones that have undergone limited metamorphism (late diagenetic zone) comprise mixed layered illite/smectite clays, and ground water is found in widely spaced fracture zones (T 〉 1 m2/d in large fracture zones; T 〈 0.1 m2/d away from fracture zones). Mudstones that have been further altered and approach the anchizone comprise illite clays, are pervasively fractured, and have the highest transmissivity values (T 〉 4 m2/d). Dolerite intrusions in unaltered, smectitic mudstones are highly fractured with transmissivity in the range of 1 〈 T 〈 60 m2/d. Thin limestone and sandstone layers can also enhance transmissivity sufficiently to provide community water supplies.