Keith J. Halford

A Computer Program for Flow-Log Analysis of Single Holes (FLASH)

A new computer program, FLASH (Flow-Log Analysis of Single Holes), is presented for the analysis of borehole vertical flow logs. The code is based on an analytical solution for steady-state multilayer radial flow to a borehole. The code includes options for (1) discrete fractures and (2) multilayer aquifers. Given vertical flow profiles collected under both ambient and stressed (pumping or injection) conditions, the user can estimate fracture (or layer) transmissivities and far-field hydraulic heads. FLASH is coded in Microsoft Excel(5) with Visual Basic for Applications routines. The code supports manual and automated model calibration.

Permeameter data verify new turbulence process for MODFLOW

Abstract A sample of Key Largo Limestone from southern Florida exhibited turbulent flow behavior along three orthogonal axes as reported in recently published permeameter experiments. The limestone sample was a cube measuring 0.2 m on edge. The published nonlinear relation between hydraulic gradient and discharge was simulated using the turbulent flow approximation applied in the Conduit Flow Process (CFP) for MODFLOW-2005 mode 2, CFPM2. The good agreement between the experimental data and the simulated results verifies the utility of the approach used to simulate the effects of turbulent flow on head distributions and flux in the CFPM2 module of MODFLOW-2005.

Arsenic Management Through Well Modification and Simulation

Arsenic concentrations can be managed with a relatively simple strategy of grouting instead of completely destroying a selected interval of well. The strategy of selective grouting was investigated in Antelope Valley, California, where groundwater supplies most of the water demand. Naturally occurring arsenic typically exceeds concentrations of 10 microg/L in the water produced from these long-screened wells. The vertical distributions of arsenic concentrations in intervals of the aquifer contributing water to selected supply wells were characterized with depth-dependent water-quality sampling and flow logs. Arsenic primarily entered the lower half of the wells where lacustrine clay deposits and a deeper aquifer occurred. Five wells were modified by grouting from below the top of the lacustrine clay deposits to the bottom of the well, which reduced produced arsenic concentrations to less than 2 microg/L in four of the five wells. Long-term viability of well modification and reduction of specific capacity was assessed for well 4-54 with AnalyzeHOLE, which creates and uses axisymmetric, radial MODFLOW models. Two radial models were calibrated to observed borehole flows, drawdowns, and transmissivity by estimating hydraulic-conductivity values in the aquifer system and gravel packs of the original and modified wells. Lithology also constrained hydraulic-conductivity estimates as regularization observations. Well encrustations caused as much as 2 microg/L increase in simulated arsenic concentration by reducing the contribution of flow from the aquifer system above the lacustrine clay deposits. Simulated arsenic concentrations in the modified well remained less than 3 microg/L over a 20-year period.

Advanced methods for modeling water-levels and estimating drawdowns with SeriesSEE, an Excel add-in

........................................................................................................................................................

Hydraulic Characterization of Overpressured Tuffs in Central Yucca Flat, Nevada Test Site, Nye County, Nevada

A sequence of buried, bedded, air-fall tuffs has been used extensively as a host medium for underground nuclear tests detonated in the central part of Yucca Flat at the Nevada Test Site. Water levels within these bedded tuffs have been elevated hundreds of meters in areas where underground nuclear tests were detonated below the water table. Changes in the ground-water levels within these tuffs and changes in the rate and distribution of land-surface subsidence above these tuffs indicate that pore-fluid pressures have been slowly depressurizing since the cessation of nuclear testing in 1992. Declines in ground-water levels concurrent with regional land subsidence are explained by poroelastic deformation accompanying ground-water flow as fluids pressurized by underground nuclear detonations drain from the host tuffs into the overlying water table and underlying regional carbonate aquifer. A hydraulic conductivity of about 3 x 10-6 m/d and a specific storage of 9 x 10-6 m-1 are estimated using ground-water flow models. Cross-sectional and three-dimensional ground-water flow models were calibrated to measured water levels and to land-subsidence rates measured using Interferometric Synthetic Aperture Radar. Model results are consistent and indicate that about 2 million m3 of ground water flowed from the tuffs to the carbonate rock as a result of pressurization caused by underground nuclear testing. The annual rate of inflow into the carbonate rock averaged about 0.008 m/yr between 1962 and 2005, and declined from 0.005 m/yr in 2005 to 0.0005 m/yr by 2300.

Detecting Drawdowns Masked by Environmental Stresses with Water-Level Models

Detecting and quantifying small drawdown at observation wells distant from the pumping well greatly expands the characterized aquifer volume. However, this detection is often obscured by water level fluctuations such as barometric and tidal effects. A reliable analytical approach for distinguishing drawdown from nonpumping water-level fluctuations is presented and tested here. Drawdown is distinguished by analytically simulating all pumping and nonpumping water-level stresses simultaneously during the period of record. Pumping signals are generated with Theis models, where the pumping schedule is translated into water-level change with the Theis solution. This approach closely matched drawdowns simulated with a complex three-dimensional, hypothetical model and reasonably estimated drawdowns from an aquifer test conducted in a complex hydrogeologic system. Pumping-induced changes generated with a numerical model and analytical Theis model agreed (RMS as low as 0.007 m) in cases where pumping signals traveled more than 1 km across confining units and fault structures. Maximum drawdowns of about 0.05 m were analytically estimated from field investigations where environmental fluctuations approached 0.2 m during the analysis period.

Predevelopment Water-Level Contours for Aquifers in the Rainier Mesa and Shoshone Mountain area of the Nevada Test Site, Nye County, Nevada

Contaminants introduced into the subsurface of the Nevada Test Site at Rainier Mesa and Shoshone Mountain by underground nuclear testing are of concern to the U.S. Department of Energy and regulators responsible for protecting human health and safety. Although contaminants were introduced into low-permeability rocks above the regional flow system, the potential for contaminant movement away from the underground test areas and into the accessible environment is greatest by ground-water transport. The primary hydrologic control on this transport is evaluated and examined through a series of contour maps developed to represent the water-level distribution within each of the major aquifers underlying the area. Aquifers were identified and their extents delineated by merging and analyzing multiple hydrostratigraphic framework models developed by other investigators from existing geologic information. The contoured water-level distribution in each major aquifer was developed from a detailed evaluation and assessment of available water-level measurements. Multiple spreadsheets that accompany this report provide pertinent water-level and geologic data by well or drill hole. Aquifers are mapped, presented, and discussed in general terms as being one of three aquifer types—volcanic aquifer, upper carbonate aquifer, or lower carbonate aquifer. Each of these aquifer types was subdivided and mapped as independent continuous and isolated aquifers, based on the continuity of its component rock. Ground-water flow directions, as related to the transport of test-generated contaminants, were developed from water-level contours and are presented and discussed for each of the continuous aquifers. Contoured water-level altitudes vary across the study area and range from more than 5,000 feet in the volcanic aquifer beneath a recharge area in the northern part of the study area to less than 2,450 feet in the lower carbonate aquifer in the southern part of the study area. Variations in water-level altitudes within any single continuous aquifer range from a few hundred feet in a lower carbonate aquifer to just more than 1,100 feet in a volcanic aquifer. Flow directions throughout the study area are dominantly southward with minor eastward or westward deviations. Primary exceptions are westward flow in the northern part of the volcanic aquifer and eastward flow in the eastern part of the lower carbonate aquifer. Northward flow in the upper and lower carbonate aquifers in the northern part of the study area is possible but cannot be substantiated because data are lacking. Interflow between continuous aquifers is evaluated and mapped to define major flow paths. These flow paths delineate tributary flow systems, which converge to form the regional ground-water flow system. The implications of these tributary flow paths in controlling transport away from the underground test areas at Rainier Mesa and Shoshone Mountain are discussed. The obvious data gaps contributing to uncertainties in the delineation of aquifers and development of water-level contours are identified and evaluated.

A unique approach to estimating lateral anisotropy in complex geohydrologic environments

Aquifers in fractured rock or karstic settings are likely to have anisotropic transmissivity distributions. Aquifer tests that are performed in these settings also are frequently affected by leakage from adjacent confining units. Finite-difference models such as MODFLOW are convenient tools for estimating the hydraulic characteristics of the stressed aquifer and adjacent confining units but are poor tools for the estimation of lateral anisotropy. This limitation of finite-difference methods can be overcome by application of the spin method, a technique whereby the positions of the observation wells are rotated about the production well to estimate anisotropy and orientation. Formal parameter estimation is necessary to analyze aquifer tests because of the number of parameters that are estimated. As a test, transmissivity, anisotropy, and orientation were successfully estimated for a simple hypothetical problem with known properties. The technique also was applied to estimate hydraulic properties of the Santee Limestone/Black Mingo (SL/BM) aquifer and a leaky confining unit beneath Charleston, South Carolina. A 9-day aquifer test with an average discharge of 6441/min was analyzed numerically. Drawdowns in the SL/BM aquifer and confining unit were simulated with a 12-layer MODFLOW model that was discretized into 81 rows of 81 columns. Simulated drawdowns at seven observation wells that ranged from 23 to 2700 m from the production well were matched to measured drawdowns. Transmissivity estimated along the minor axis ranged from 10 to 15 m2/day and along the major axis ranged from 80 to 10Om2/day. The major axis of transmissivity was oriented along compass heading 116° (degrees clockwise from north), which agrees with geologic interpretations. Vertical hydraulic conductivity and specific storage estimates for the overlying confining unit were 4 x 10-5 m/day and 2 x 10-4 1/m, respectively.

Scenario Evaluator for Electrical Resistivity Survey Pre-modeling Tool

Geophysical tools have much to offer users in environmental, water resource, and geotechnical fields; however, techniques such as electrical resistivity imaging (ERI) are often oversold and/or overinterpreted due to a lack of understanding of the limitations of the techniques, such as the appropriate depth intervals or resolution of the methods. The relationship between ERI data and resistivity is nonlinear; therefore, these limitations depend on site conditions and survey design and are best assessed through forward and inverse modeling exercises prior to field investigations. In this approach, proposed field surveys are first numerically simulated given the expected electrical properties of the site, and the resulting hypothetical data are then analyzed using inverse models. Performing ERI forward/inverse modeling, however, requires substantial expertise and can take many hours to implement. We present a new spreadsheet-based tool, the Scenario Evaluator for Electrical Resistivity (SEER), which features a graphical user interface that allows users to manipulate a resistivity model and instantly view how that model would likely be interpreted by an ERI survey. The SEER tool is intended for use by those who wish to determine the value of including ERI to achieve project goals, and is designed to have broad utility in industry, teaching, and research.

Evaluating Micrometeorological Estimates of Groundwater Discharge from Great Basin Desert Playas: T.R. Jackson et al. Groundwater xx, no. x: xx-xx

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg ) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2% of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.