ALBERT

Notes: A convolution integral is developed to evaluate transient, drilling speed–dependent discharge rates into a tunnel gradually excavated in a homogeneous infinite aquifer. Comparison with the classical instantaneous-drilling evaluation commonly used in practice reveals drastically reduced early-time discharge rates, higher maximum rates, and similar long-term rates. Dimensionless-type curves are provided to help assess total discharge sensitivity to drilling time and predict safer maximum flow rates.

Notes: Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with 〉2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.

Notes: Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models—the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)—to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain.

Notes: System dynamics is a computer-aided approach to evaluating the interrelationships of different components and activities within complex systems. Recently, system dynamics models have been developed in areas such as policy design, biological and medical modeling, energy and the environmental analysis, and in various other areas in the natural and social sciences. The Idaho National Engineering and Environmental Laboratory, a multipurpose national laboratory managed by the Department of Energy, has developed a system dynamics model in order to evaluate its utility for modeling large complex hydrological systems. We modeled the Bear River basin, a transboundary basin that includes portions of Idaho, Utah, and Wyoming. We found that system dynamics modeling is very useful for integrating surface water and ground water data and for simulating the interactions between these sources within a given basin. In addition, we also found that system dynamics modeling is useful for integrating complex hydrologic data with other information (e.g., policy, regulatory, and management criteria) to produce a decision support system. Such decision support systems can allow managers and stakeholders to better visualize the key hydrologic elements and management constraints in the basin, which enables them to better understand the system via the simulation of multiple “what-if” scenarios. Although system dynamics models can be developed to conduct traditional hydraulic/hydrologic surface water or ground water modeling, we believe that their strength lies in their ability to quickly evaluate trends and cause-effect relationships in large-scale hydrological systems, for integrating disparate data, for incorporating output from traditional hydraulic/hydrologic models, and for integration of interdisciplinary data, information, and criteria to support better management decisions.

Notes: This paper presents a general model of optimal water management for a transboundary aquifer under three different management approaches: cooperative, noncooperative, and myopic. Comparing the results from the approaches, we find the cooperative solution, where a single management plan is executed for all parts of the aquifer, results in the highest level of net social welfare, followed by the noncooperative and then the myopic. The trade-offs we find for the higher levels of welfare are lower use levels in the earlier periods. We present a short discussion of factors that can increase the inefficiencies of the models and suggest directions for future research.

Notes: The goal for any property rights system is to achieve equity, efficiency, and certainty. Trying to achieve these goals for ground water is difficult because a ground water right is not exclusive. To make matters more complicated, ground water is often under the jurisdiction of more than one political unit. The result is transboundary conflicts. Two critical elements must be included in any system of ground water rights. The system must define how the ground water can be used and define the relationships that each user and each use has with the other users and uses in the system. Unfortunately, these relationships are seldom completely defined and are made more complex by the transboundary scales at which they operate. As ground water moves horizontally across boundaries, different users or different jurisdictions have sequential control, creating conflicts between the first users and subsequent ones. Other problems occur because of vertical relationships, with more than one person or entity having control over ground water at the same time. This simultaneous exercise of authority can create conflicts between an individual who possesses a right to use ground water and a state or federal agency that regulates the same water. Transboundary conflicts occur at different scales and include conflicts between neighboring property owners as well as conflicts between countries. Scale, the property rights structure, and the nature of the relationship between users influence the way transboundary ground water conflicts are resolved.

Notes: In Europe, a long history of cooperation over transboundary rivers—most notably the Rhine and Danube rivers—exists. To help foster cooperation and communication vis-à-vis transboundary ground water, the United Nations Economic Commission for Europe (UNECE), as part of its ground water program, conducted a survey on transboundary aquifers in Europe. The survey produced 25 responses from 37 countries and identified 89 transboundary aquifers. Respondents reported on the degree of ground water use within their own boundaries, transboundary aspects (agreements, joint commissions, etc.) of ground water, and transboundary aquifers themselves. The inventory proved useful, but a number of problems were identified: different map scales and symbols, difficulty in identifying transboundary aquifers, inconsistent labeling of aquifers, and data discrepancies. The UNECE ground water program also drafted guidelines for monitoring and assessment of transboundary ground water. These guidelines are not legally binding but have been adopted by 25 countries, deal mainly with monitoring and assessment, and are being implemented through a number of pilot projects. Other organizations—the United Nations Scientific, Educational and Cultural Organization, the Food and Agriculture Organization, the International Association of Hydrogeologists, and the European Union—are all supporting the investigation of transboundary aquifers in an effort to facilitate data sharing and coordinated management of these valuable resources.

Notes: More than one-half of the world's population is dependent on ground water for everyday uses such as drinking, cooking, and hygiene. In fact, it is the most extracted natural resource in the world. As a result of growing populations and expanding economies, many aquifers today are being depleted while others are being contaminated. Notwithstanding the world's considerable reliance on this resource, ground water resources have long received only secondary attention as compared to surface water, especially among legislatures and policymakers. Today, while there are hundreds of treaties governing transboundary rivers and lakes, there is only one international agreement that directly addresses a transboundary aquifer. Given that many of the aquifers on which humanity so heavily relies cross international borders, there is a considerable gap in the sound management, allocation, and protection of such resources. In order to prevent future disputes over transboundary aquifers and to maximize the beneficial use of this resource, international law must be clarified as it applies to transboundary ground water resources. Moreover, it must be defined with a firm basis in sound scientific understanding. In this paper we offer six conceptual models is which ground water resources can have transboudary consequences. The models are intended to help in assessing the applicability and scientific soundness of existing and proposed rules governing transboundary ground water resources. In addition, we consider the development of international law as it applies to ground water resources and make recommendations based on the models and principles of hydrogeology. The objective is the development of clear, logical, and science-based norms of state conducts as they relate to aquifers that traverse political boundaries.

Notes: A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation.

Notes: Previous studies have used flowmeters in environments that are within the expectations of their published ranges. Electromagnetic flowmeters have a published range from 0.1 to 79.0 m/min, and impeller flowmeters have a published range from 1.2 to 61.0 m/min. Velocity-log data collected in five long-screened production wells in the Pleasant Valley area of southern California showed that (1) electromagnetic flowmeter results were comparable within ±2% to results obtained using an impeller flowmeter for comparable depths; (2) the measured velocities from the electromagnetic flowmeter were up to 36% greater than the published maximum range; and (3) both data sets, collected without the use of centralizers or flow diverters, produced comparable and interpretable results. Although either method is acceptable for measuring wellbore velocities and the distribution of flow, the electromagnetic flowmeter enables collection of data over a now greater range of flows. In addition, changes in fluid temperature and fluid resistivity, collected as part of the electromagnetic flowmeter log, are useful in the identification of flow and hydrogeologic interpretation.

Notes: While the discussion of model uncertainty has centered on spatial heterogeneity, it is possible that ground water models have not enjoyed much success as predictive tools often because the sources that were eventually imposed in the field differed from those represented in the simulations. This is because deterministic prediction of future conditions is often inaccurate due to the random nature of contaminant sources, in terms of their timing, location, and magnitude. This paper presents a stochastic framework for accommodating random contaminant sources in conventional, deterministic advection-dispersion transport models. The contaminant sources are first classified into two types: those occurring continuously with a deterministic component and random variations and those occurring randomly at instantaneous discrete-time intervals. For the first type, the governing partial differential equation (PDE) is replaced by a stochastic PDE. The random variations are modeled by Gaussian noise or Brownian motion, and the solution is obtained by using Ito's integration technique. For the second type, Markovian analysis is used for discrete-time contamination events. Both approaches use a deterministic transport model to generate response functions at any observation location and time. The response functions are then integrated to yield probabilistic description of contaminant transport, from which key statistical properties such as mean, standard deviation, and confidence interval can be drawn.

Notes: A model of porosity development has been developed to investigate general relationships between simple fracture aperture growth laws and fracture porosity in evolved fracture arrays in aquifers. The growth of fracture apertures in two-dimensional orthogonal arrays with initially spatially uncorrelated lognormal aperture distributions has been studied, where aperture growth rate is proportional to an exponent of the flow rate through each fracture. The evolved arrays show geometrical phase changes as a function of the aperture growth rate exponent, e, and the standard deviation of the initial aperture distribution, σz. Low values of e and σz lead to bimodal aperture distributions, where apertures parallel to flow are preferentially enlarged. At moderate values of e and σz, there is a transition to a regime of more complex geometries consisting of networks of channel-like structures of preferentially enlarged apertures. At larger values of e, array-spanning channel-like paths of preferentially enlarged apertures develop, where the tortuosity of the channel-like paths is a linear function of σz. Following an initial growth phase, during which dynamically stable aperture configurations develop, arrays undergo simple amplification. The geometry of the evolved aperture fields is diverse and they can be highly complex; consequently, parameterization and prediction of their evolution in terms of the initial aperture distributions and growth rate laws is not trivial.

Notes: A field-scale demonstration project was conducted to evaluate the capability of eastern cottonwood trees (Populus deltoides) to attenuate trichloroethene (TCE) contamination of ground water. By the middle of the sixth growing season, trees planted where depth to water was 〈3 m delivered enough dissolved organic carbon to the underlying aquifer to lower dissolved oxygen concentrations, to create iron-reducing conditions along the plume centerline and sulfate-reducing or methanogenic conditions in localized areas, and to initiate in situ reductive dechlorination of TCE. Apparent biodegradation rate constants for TCE along the centerline of the plume beneath the phytoremediation system increased from 0.0002/d to 0.02/d during the first six growing seasons. The corresponding increase in natural attenuation capacity of the aquifer along the plume centerline, from 0.0004/m to 0.024/m, is associated with a potential decrease in plume-stabilization distance from 9680 to 160 m. Demonstration results provide insight into the amount of vegetation and time that may be needed to achieve cleanup objectives at the field scale.

Notes: The hydraulic conductivity (K) variation has important ramifications for ground water flow and the transport of contaminants in ground water. The delineation of the nature of that variation can be critical to complete characterization of a site and the planning of effective and efficient remedial measures. Site-specific features (such as high-conductivity zones) need to be quantified. Our alluvial field site in the Kansas River valley exhibits spatial variability, very high conductivities, and nonlinear behavior for slug tests in the sand and gravel aquifer. High-resolution, multilevel slug tests have been performed in a number of wells that are fully screened. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been used to analyze the data, employing an automated processing system that runs within the Excel spreadsheet program. It is concluded that slug tests can provide the necessary data to identify the nature of both horizontal and vertical K variation in an aquifer and that improved delineation or higher resolution of K structure is possible with shorter test intervals. The gradation into zones of higher conductivity is sharper than seen previously, and the maximum conductivity observed is greater than previously measured. However, data from this project indicate that well development, the presence of fines, and the antecedent history of the well are important interrelated factors in regard to slug-test response and can prevent obtaining consistent results in some cases.

Notes: To improve understanding of property measurements in heterogeneous media, an energy-based weighting function concept is developed. In (assumed) homogeneous media, the instrument spatial weighting function (ISWF) depends only on the energy dissipation distribution set up by the measurement procedure and it reduces to simply inverse sample volume (uniform weighting) for 1-D parallel flow case (ideal permeameter). For 1-D transient flow in homogeneous media, such as with slug tests, the ISWF varies with position and time, with 95% of the total weighting contained within 115 well radii, even late in the test. In the heterogeneous case, the determination of the ISWF is connected to the problem of determining an equivalent hydraulic conductivity (K), where the criterion for equivalence is based on equal energy dissipation rate rather than equal volume discharge. The discharge-based equivalent K (KE) and the energy-based equivalent K in heterogeneous media (Keh) are not equal in general, with Keh typically above the nodal arithmetic mean K. The possibly more fundamental problem is that as one makes K measurements in heterogeneous media at different locations or on different cores of heterogeneous materials, the ISWF will be heterogeneity dependent, implying that the averaging process resulting in the equivalent K value also varies with position. If the testing procedure is transient, then the averaging process varies with time. This suggests a fundamental ambiguity in the interpretation of hydraulic conductivity measurements in heterogeneous media that may impact how we approach modeling and prediction in a practical sense (Molz 2003). Further research is suggested.

Notes: At a number of sites, a plume of methyl tert-butyl ether (MTBE) in ground water has dived below the screen of conventional monitoring wells and escaped detection. Techniques are needed to predict the vertical extent of MTBE in an aquifer. Two techniques that are emerging in the site characterization market are electrical conductivity logging and pneumatic slug testing performed in temporary push wells. These techniques were evaluated at a diving plume of MTBE in the aquifer that supplies water to the village of East Alton, Illinois. The plume stayed near the water table for the first 100 m from the potential sources and then dived below conventional monitoring over the next 100 m. At the location where the MTBE plume dived, the depth to water was 9.1 m below land surface. The first 10 m of material below the water table had an electrical conductivity near 100 mS/m, indicating silts and clays. An electrical conductivity near 25 mS/m, indicating sands or gravels, was encountered at a depth of 10.6 m below the water table, and the sands and gravel extended to a depth of at least 15.2 m below the water table. Pneumatic slug tests measured low hydraulic conductivity in the interval of silt and clay (0.34 and 0.012 m/d) and higher hydraulic conductivity in the interval with sands and gravels (12.5, 11.6, and 11.3 m/d). Ground water with the highest concentration of MTBE was produced just below the contact between the silt and clay and the sands and gravel.

Notes: Steam injected below the water table tends to move upward because of buoyancy. This limits the horizontal steam zone development, which determines the optimal spacing between injection wells. In this study, observations indicating steam override at a full-scale remediation of an unconfined aquifer are analyzed by numerical modeling using the code T2VOC. A simplified three-dimensional (3D) numerical model is set up, which qualitatively shows the same mechanisms as observed at the site. By means of the model, it is found that it will be possible to achieve a larger horizontal extent of the steam zone in a layered geology compared to the homogeneous case. In the homogeneous case, the steam injection rate increases dramatically when the injection pressure is increased, which is necessary to achieve a larger horizontal development. The development of the steam zone under unconfined conditions is found to be a complex function of the geologic layering, the water table at steady-state extraction, and the injection/extraction system. Because of this complexity, it will be difficult to predict steam behavior without site-specific 3D numerical modeling.

Notes: This work describes the integration of data acquisition (DAQ) hardware and software for the purpose of acquiring not only data but real-time transport model parameter estimates in the context of subsurface flow and transport problems. Integrated DAQ parameter estimation systems can be used to reduce data storage requirements, trigger event recognition and more detailed sampling actions, and otherwise enhance remote monitoring capabilities. The contaminant transport problem is posed here as the analogous heat transfer problem in a three-dimensional, intermediate-scale physical aquifer model. A constant source of warm water is fed into a sandy aquifer undergoing steady, unidirectional flow. The spatial distribution of temperature in the medium is monitored over time using 17 thermocouples embedded in the medium. These sensors log temperatures via conventional analog-to-digital conversion hardware driven by commercially available DAQ software (LabVIEW™). Parameter estimation routines programmed in MATLAB™-based M-files are embedded in the LabVIEW DAQ routine and access parameter estimation libraries, such as the descent method employed here, via the Internet. The integrated DAQ parameter estimation system is demonstrated for the estimation of (1) the thermal dispersion coefficients (analogous to mass dispersion coefficients), given a known heat source; and (2) the location of a heat source, given known thermal dispersion coefficients. In both cases, the parameter estimation procedure is executed repeatedly as the data are acquired. For the case of source location, the effect of the number of sensors on the parameter estimation procedure is also demonstrated. Reasonable parameter estimates are provided rapidly during both the transient and steady-state phases of the experiments, with accuracy increasing with time and with the number of observations employed.

Notes: Herbicide and herbicide degradate compounds were measured in tile drains, monitoring wells, domestic wells, and in the receiving stream at a farm in south-central Michigan, in which corn, soybeans, and livestock are the major products. Immunoassay screens of two herbicide groups (triazines and alachlor) in 2000 showed a decrease in the number of detections from June to September, after the well-known spring flush, of 48% to 26% for alachlor and 62% to 22% for the triazines. Gas chromatography/mass spectrometry and liquid chromatography/mass spectrometry analyses of 2001 samples documented a fall flush as well, mobilizing relatively high levels of both parent and degradate compounds to the tile drain system. The compounds included both metolachlor and atrazine, which were applied on the farm earlier that year, as well as the degradates of alachlor, which had not been applied on the farm for several years. In four tile drain samples, metolachlor was detected at concentrations ranging from 18 to 25 μg/L, whereas atrazine ranged from nondetect to 2.25 μg/L. Contrary to previous studies, metolachlor concentrations were higher than metolachlor ethane sulfonic acid, by a factor of 3 to 7. The high soil permeability, high water table, and shallow screened intervals of monitoring wells directly below the point of application may explain the lesser degree of breakdown at the point of sampling. Although most detections were below maximum contaminant level or health advisory levels, the unknown effects of degradates and the aggregate levels of compounds measured in many samples may pose higher health risks than for individual compounds.

Notes: Monitoring data collected over a 6-year period show that a plume of chlorinated ethene–contaminated ground water has contracted significantly following treatment of the contaminant source area using in situ oxidation. Prior to treatment (1998), concentrations of perchloroethene (PCE) exceeded 4500 μg/L in a contaminant source area associated with a municipal landfill in Kings Bay, Georgia. The plume emanating from this source area was characterized by vinyl chloride (VC) concentrations exceeding 800 μg/L. In situ oxidation using Fenton's reagent lowered PCE concentrations in the source area below 100 μg/L, and PCE concentrations have not rebounded above this level since treatment. In the 6 years following treatment, VC concentrations in the plume have decreased significantly. These concentration declines can be attributed to the movement of Fenton's reagent–treated water downgradient through the system, the cessation of a previously installed pump-and-treat system, and the significant natural attenuation capacity of this anoxic aquifer. While in situ oxidation briefly decreased the abundance and activity of microorganisms in the source area, this activity rebounded in 〈6 months. Nevertheless, the shift from sulfate-reducing to Fe(III)-reducing conditions induced by Fenton's treatment may have decreased the efficiency of reductive dechlorination in the injection zone. The results of this study indicate that source-area removal actions, particularly when applied to ground water systems that have significant natural attenuation capacity, can be effective in decreasing the areal extent and contaminant concentrations of chlorinated ethene plumes.

Notes: Two small-diameter (≤0.5-inch/13-mm O.D.) bladder pumps were evaluated under the U.S. EPA Environmental Technology Verification Program to assess their ability to provide representative water quality samples for several volatile organic compounds (VOCs) and major element cations. Both a pneumatic bladder pump (PnBP) and a newly designed mechanical bladder pump (MBP) were evaluated. Since the MBP does not require compressed air to operate, the field logistics are simplified and cost of conducting low-flow sampling can be reduced. Replicate samples were collected from an above-ground standpipe spiked with selected VOCs and cations. Low-flow sampling protocol was followed during collection of replicate samples from six direct push (DP)–installed monitoring wells at Tyndall Air Force Base, Florida. Simple statistical methods are used to evaluate the sample quality provided by the pumps as compared to reference samples. F-tests are used to evaluate pump sample precision, and t-tests are applied to evaluate comparability of the pump sample results to the colocated reference or port sample results. Statistical analysis reveals that the samples collected with the pumps provide precision comparable to that observed in the reference samples. Significantly greater variability was observed between the pump and reference samples collected from the DP wells than was observed during the standpipe trials. The results and statistics from the trials at the standpipe indicate that both the PnBP and MBP provide water quality samples that are comparable to the reference samples, for the analytes investigated.

Notes: Previous site-specific studies designed to assess the impacts of unsewered subdivisions on ground water quality have relied on upgradient monitoring wells or very limited background data to characterize conditions prior to development. In this study, an extensive monitoring program was designed to document ground water conditions prior to construction of a rural subdivision in south-central Wisconsin. Previous agricultural land use has impacted ground water quality; concentrations of chloride, nitrate-nitrogen, and atrazine ranged from below the level of detection to 296 mg/L, 36 mg/L, and 0.8 μg/L, respectively, and were highly variable from well to well and through time. Seasonal variations in recharge, surface topography, aquifer heterogeneities, surficial loading patterns, and well casing depth explain observed variations in ground water chemistry. This variability would not have been detected if background conditions were determined from only a few monitoring wells or inferred from wells located upgradient of the subdivision site. This project demonstrates the importance of characterizing both ground water quality and chemical variability prior to land-use change to detect any changes once homes are constructed.

Notes: From the mid-1940s through the 1980s, large volumes of waste water were discharged at the Hanford Site in southeastern Washington State, causing a large-scale rise (〉20 m) in the water table. When waste water discharges ceased in 1988, ground water mounds began to dissipate. This caused a large number of wells to go dry and has made it difficult to monitor contaminant plume migration. To identify monitoring wells that will need replacement, a methodology has been developed using a first-order uncertainty analysis with UCODE, a nonlinear parameter estimation code. Using a three-dimensional, finite-element ground water flow code, key parameters were identified by calibrating to historical hydraulic head data. Results from the calibration period were then used to check model predictions by comparing monitoring wells’ wet/dry status with field data. This status was analyzed using a methodology that incorporated the 0.3 cumulative probability derived from the confidence and prediction intervals. For comparison, a nonphysically based trend model was also used as a predictor of wells’ wet/dry status. Although the numerical model outperformed the trend model, for both models, the central value of the intervals was a better predictor of a wet well status. The prediction interval, however, was more successful at identifying dry wells. Predictions made through the year 2048 indicated that 46% of the wells in the monitoring well network are likely to go dry in areas near the river and where the ground water mound is dissipating.

Notes: A method is proposed to estimate the uncertainty of the location of pathlines in two-dimensional, steady-state confined or unconfined flow in aquifers due to the uncertainty of the spatially variable unconditional hydraulic conductivity or transmissivity field. The method is based on concepts of the semianalytical first-order theory given in Stauffer et al. (2002, 2004), which allows estimates of the lateral second moment (variance) of the location of a moving particle. However, this method is reformulated in order to account for nonuniform recharge and nonuniform aquifer thickness. One prominent application is the uncertainty estimation of the catchment of a pumping well by considering the boundary pathlines starting at a stagnation point. In this method, the advective transport of particles is considered, based on the velocity field. In the case of a well catchment, backtracking is applied by using the reversed velocity field. Spatial variability of hydraulic conductivity or transmissivity is considered by taking into account an isotropic exponential covariance function of log-transformed values with parameters describing the variance and correlation length. The method allows postprocessing of results from ground water models with respect to uncertainty estimation. The code PPPath, which was developed for this purpose, provides a postprocessing of pathline computations under PMWIN, which is based on MODFLOW. In order to test the methodology, it was applied to results from Monte Carlo simulations for catchments of pumping wells. The results correspond well. Practical applications illustrate the use of the method in aquifers.

Notes: Extremely alkaline ground water has been found underneath many shuttered steel mills and slag dumps and has been an impediment to the cleanup and economic redevelopment of these sites because little is known about the geochemistry. A large number of these sites occur in the Lake Calumet region of Chicago, Illinois, where large-scale infilling of the wetlands with steel slag has created an aquifer with pH values as high as 12.8. To understand the geochemistry of the alkaline ground water system, we analyzed samples of ground water and the associated slag and weathering products from four sites. We also considered several potential remediation schemes to lower the pH and toxicity of the water. The principal cause of the alkaline conditions is the weathering of calcium silicates within the slag. The resulting ground water at most of the sites is dominated by Ca2+ and OH− in equilibrium with Ca(OH)2. Where the alkaline ground water discharges in springs, atmospheric CO2 dissolves into the water and thick layers of calcite form. Iron, manganese, and other metals in the metallic portion of the slag have corroded to form more stable low-temperature oxides and sulfides and have not accumulated in large concentrations in the ground water. Calcite precipitated at the springs is rich in a number of heavy metals, suggesting that metals can move through the system as particulate matter. Air sparging appears to be an effective remediation strategy for reducing the toxicity of discharging alkaline water.

Notes: Artificial neural networks (ANNs) were developed to accurately predict highly time-variable specific conductance values in an unconfined coastal aquifer. Conductance values in the fresh water lens aquifer change in response to vertical displacements of the brackish zone and fresh water–salt water interface, which are caused by variable pumping and climate conditions. Unlike physical-based models, which require hydrologic parameter inputs, such as horizontal and vertical hydraulic conductivities, porosity, and fluid densities, ANNs can “learn” system behavior from easily measurable variables. In this study, the ANN input predictor variables were initial conductance, total precipitation, mean daily temperature, and total pumping extraction. The ANNs were used to predict salinity (specific conductance) at a single monitoring well located near a high-capacity municipal-supply well over time periods ranging from 30 d to several years. Model accuracy was compared against both measured/interpolated values and predictions were made with linear regression, and in general, excellent prediction accuracy was achieved. For example, although the average percent change of conductance over 90-d periods was 39%, the absolute mean prediction error achieved with the ANN was only 1.1%. The ANNs were also used to conduct a sensitivity analysis that quantified the importance of each of the four predictor variables on final conductance values, providing valuable insights into the dynamics of the system. The results demonstrate that the ANN technology can serve as a powerful and accurate prediction and management tool, minimizing degradation of ground water quality to the extent possible by identifying appropriate pumping policies under variable and/or changing climate conditions.

Notes: A new distance-drawdown method for aquifers with anisotropy on the horizontal plane is presented. The method uses scalar transformation to convert to an equivalent, isotropic medium, thus permitting application of the Cooper-Jacob Method. The method is applicable to cases where at least one ellipse of equal drawdown can be delineated but can also be applied where no ellipse can be discerned from the data. In the latter case, a least-squares regression approach can be employed to estimate the orientation and magnitude of the anisotropy. The regression R2 value provides a quantitative assessment of the degree to which the drawdown data are indicative of a systematic areal anisotropy in the aquifer or whether the data simply reflect natural aquifer heterogeneity. In addition to confined aquifers, this methodology, like the Cooper-Jacob Method, is also applicable to unconfined aquifers either before the onset of delayed drainage or following the completion of delayed drainage provided that the u value meets the recommended criterion.

Notes: A new multilayer approach is presented for the modeling of ground water flow to radial collector wells. The approach allows for the inclusion of all aspects of the unique boundary condition along the lateral arms of a collector well, including skin effect and internal friction losses due to flow in the arms. The hydraulic conductivity may differ between horizontal layers within the aquifer, and vertical anisotropy can be taken into account. The approach is based on the multilayer analytic element method, such that regional flow and local three-dimensional detail may be simulated simultaneously and accurately within one regional model. Horizontal flow inside a layer is computed analytically, while vertical flow is approximated with a standard finite-difference scheme. Results obtained with the proposed approach compare well to results obtained with three-dimensional analytic element solutions for flow in unconfined aquifers. The presented approach may be applied to predict the yield of a collector well in a regional setting and to compute the origin and residence time, and thus the quality, of water pumped by the collector well. As an example, the addition of three lateral arms to a collector well that already has three laterals is investigated. The new arms are added at an elevation of 2 m above the existing laterals. The yield increase of the collector well is computed as a function of the lengths of the three new arms.

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.