Search USGSSearch

Geology topics

Keith J. Halford

Publications and source records attributed to Keith J. Halford.

At least 19 recordsLinked to original sources

Hydraulic characterization of carbonate-rock and basin-fill aquifers near Long Canyon, Goshute Valley, northeastern Nevada

Understanding groundwater flow and pumping effects near pending mining operations requires accurate subsurface hydraulic characterization. To improve conceptual models of groundwater flow and development in the complex hydrogeologic system near Long Canyon Mine, in northwestern Goshute Valley, northeastern Nevada, the U.S. Geological Survey characterized the hydraulic properties of carbonate rocks and basin-fill aquifers using an integrated analysis of steady-state and stressed aquifer conditions informed by water chemistry and aquifer-test data. Hydraulic gradients and groundwater-age data in northern Goshute Valley indicate carbonate rocks in the Pequop Mountains just west and south of the Long Canyon Mine project area constitute a more permeable and active flow system than saturated rocks in the northern Pequop Mountains, western Toano Range, and basin fill. Permeable carbonate rocks in the northern Pequop Mountains, in part, discharge to the Johnson Springs wetland complex (JSWC), where mean groundwater ages range from 500 to 2,400 years and samples all contain a small fraction of modern waters, relative to mean ages of 8,600 to more than 22,000 years for most groundwater sampled to the north and east. Recharge to the JSWC occurs from a roughly 27-square-mile area in the upgradient Pequop Mountains to the west, composed mostly of permeable carbonate rock and fractured quartzite, and bounded by low-permeability shales and marbleized and siliclastic rocks. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from 7,000 to 400,000 feet squared per day (ft 2 /d) in carbonate rocks and from 2,000 to 80,000 ft 2 /d in basin fill near the Long Canyon Mine. Water-level drawdown from multiple-well aquifer testing and rise from unintentional leakage into the overlying basin-fill aquifer were estimated and distinguished from natural fluctuations in 93 pumping and monitoring sites using analytical water-level models. Leakage of disposed aquifer-test pumpage occurred south of the aquifer test area through an unlined irrigation ditch. Drawdown was detected at distances of as much as 3 miles (mi) from pumping wells at all but one carbonate-rock site, at basin-fill sites on the alluvial fan immediately downgradient from pumping wells, and in Big Spring and spring NS-05. Similar drawdowns in carbonate rocks within the drawdown detection area suggest all wells penetrate a highly transmissive zone (HTZ) that is bounded by low-permeability rocks. Drawdown was not detected in carbonate rocks to the west of Canyon fault, in any basin-fill sites on the valley floor east of the Hardy fault, or at volcanic sites to the north, indicating that these major fault structures and (or) permeability contrasts between hydrogeologic units impeded groundwater flow or obscured pumping signals. Alternatively, unintentional leakage might have obscured drawdown at basin-fill sites on the valley floor, where water-level rise was detected at nine sites over 3 mi. Consistent hydraulic properties were estimated by simultaneously interpreting steady-state flow during predevelopment conditions and changes in groundwater levels and springflows from the 2016 carbonate-rock aquifer test with an integrated groundwater-flow model. Hydraulic properties were distributed across carbonate rocks, basin fill, volcanic rocks, and siliciclastic rocks with a hydrogeologic framework developed from geologic mapping and hydraulic testing. Estimated transmissivity distributions spanned at least three orders of magnitude in each rock unit. In the HTZ, simulated transmissivities ranged from 10,000 to 23,000,000 ft 2 /d, with the most transmissive areas occurring around Big Spring. Comparatively low carbonate-rock transmissivities of less than 10,000 ft 2 /d were estimated in the northern Pequop Mountains and poorly defined values of less than 1,000 ft 2 /d were estimated in the western Toano Range. Transmissivities in basin fill ranged from less than 10 to 80,000 ft 2 /d and were minimally constrained by the 2016 carbonate-rock aquifer test because poorly quantified leakage affected water levels more so than pumping. The most transmissive areas were informed by single-well aquifer tests along the eastern edge of the Pequop Mountains near Long Canyon Mine and could be indicative of a hydraulic connection between basin fill and more transmissive underlying carbonate rocks. Simulated transmissivities of volcanic and low-permeability rocks mostly are less than 1,000 ft 2 /d. The estimated hydraulic-property distributions and informed interpretation of hydraulic connections among hydrogeologic units improved the characterization and representation of groundwater flow near the Long Canyon Mine.

Nevada

Groundwater characterization and effects of pumping in the Death Valley regional groundwater flow system, Nevada and California, with special reference to Devils Hole

Groundwater flow and development were characterized in four groundwater basins of the Death Valley regional flow system in Nevada and California with calibrated, groundwater-flow models. Natural groundwater discharges in the Furnace Creek, Lower Amargosa, and Saratoga Spring areas were defined and distributed consistently with a revised hydrogeologic framework. This simplified hydrogeologic framework was limited to four hydraulically unique, hydrogeologic units: (1) basin fill; (2) carbonate rocks; (3) volcanic rocks; and (4) low-permeability granitic and siliciclastic rocks. Hydrogeologic units and division of carbonate and volcanic rocks between shallow and deep were supported by results from 271 aquifer tests and specific-capacity estimates. Greater than 90 percent of field-estimated transmissivity occurred within 1,600 feet (ft) of the water table. Pumping in the study area from 1960 to 2010 averaged 46,000 acre-feet per year (acre-ft/yr), which is 80 percent of the predevelopment discharge. The central Amargosa Desert and Pahrump Valley were the two primary pumping centers and measurably affected water levels across 900 square miles in 2018. Water levels in Devils Hole were a special focus because endangered Devils Hole pupfish ( Cyprinodon diabolis ) are affected by water-level declines. Pumping 42,100 acre-ft by Cappaert Enterprises, formerly Spring Meadows, Inc., caused a 2.3-ft water-level decline in Devils Hole , which temporarily reduced habitat of Devils Hole pupfish by 85 percent in 1972. If no pumping occurred, water levels in Devils Hole would have risen naturally about 1 ft between 1973 and 2018 from temporal variations in recharge. The 2.6-ft range of measured water-level changes in Devils Hole was simulated with a root-mean-square error of 0.2 ft during the 70-year period of record. Simulated water-level declines from pumping totaled 1.4 ft in 2018, with 25 and 34 percent attributed to pumping by Cappaert Enterprises and the central Amargosa Desert, respectively. Water levels in Devils Hole will decline at rates of 0.1–0.2 ft per decade if pumping from Ash Meadows groundwater basin and the central Amargosa Desert continue at current rates. Effects of future natural water-level fluctuations remain unknown. Ash Meadows and Alkali Flat–Furnace Creek Ranch groundwater basins are hydraulically connected near well AD-4 , about 5 miles south of the town of Amargosa Valley, Nevada. About 40 percent of the discharge from the Furnace Creek area is recharged in the Ash Meadows groundwater basin. Basin fill in the central Amargosa Desert hydraulically connects carbonate rocks east of well AD-4 with saturated carbonate rocks in the Funeral Range. About 7 percent of the 960,000 acre-ft pumped from Ash Meadows and Alkali Flat–Furnace Creek Ranch groundwater basins prior to 2019 was captured discharge from springs and phreatophytes. Greater than 40 percent of the 2,080,000 acre-ft pumped from Pahrump Valley between 1910 and 2019 was capture that primarily discharged from Bennetts and Manse Springs. Simulated advective-flow distances and velocities from underground nuclear tests are within the range of advective transport calculations from tritium data and previous radionuclide transport investigations. Boundary conditions and flow rates from the regional model in this study are plausible for local-scale flow and radionuclide transport models. Simulated 165-year groundwater-flow paths do not extend into pumping areas and effects of regional pumping on advective transport are negligible.

California, Nevada

Estimation of groundwater flow through Yucca Flat based on a multiple-well aquifer test at well ER-6–1–2 main , Nevada National Security Site, southern Nevada

The rate of groundwater flow past underground nuclear testing areas in Yucca Flat at the Nevada National Security Site, southern Nevada, was estimated using results from the ER-6-1-2 main multiple-well aquifer test (MWAT), done during February 5–July 23, 2004. Drawdowns in 13 observation wells were evaluated from pumping in well ER-6-1-2 main and used as observations in a regional groundwater-flow model to estimate the groundwater-flow rate. Flow from southern Yucca Flat regionally was constrained by detecting drawdown in Tracer Well 3 , which is 33 miles south-southwest of pumping well ER-6-1-2 main . Predevelopment flow, the ER-6-1-2 main MWAT, and regional groundwater pumping were simulated in separate models that jointly informed estimates of transmissivity and specific-yield distributions. Predevelopment flow constrained groundwater-flow estimates from southern Yucca Flat because the steady-state model correctly simulated regional groundwater discharges in the Ash Meadows groundwater basin and the surrounding basins. Less than 600 acre-feet per year of groundwater moves from underground nuclear testing areas and flows south of well ER-6-1-2 main in southern Yucca Flat.

Nevada

Evaluation of groundwater-flow models for estimating drawdown from proposed groundwater development in Tule Desert, Nevada

At the request of the Bureau of Land Management (BLM), the U.S. Geological Survey (USGS) is releasing with this open-file report (OFR) a previously unpublished review and comparison of two numerical models for Tule Desert, Nevada. The original review was performed in spring 2013, and only minor editorial revisions were made in the current (2019) OFR for clarity and to reformat the original interagency correspondence to the USGS OFR template. No revisions have been made to the technical content of the original review for this OFR release. Report content presented in the purpose and scope statement, and all subsequent sections of the OFR, are original content submitted to BLM in May 2013. Model review and comparisons described in the following paragraphs are based on, in part, results of a long-term (more than 2 years) aquifer test mandated by Nevada State Engineer Order 1169. Additional information on Order 1169 and associated aquifer test results can be found at the State of Nevada Division of Water Resources website (State of Nevada, 2019).

Nevada

Documentation of single-well aquifer tests and integrated borehole analyses, Pahute Mesa and Vicinity, Nevada

Single-well aquifer testing has been carried out at Pahute Mesa in southern Nevada since 1962. These tests include single-well pumping and slug tests to estimate geologic formation hydraulic properties. Initially, aquifer tests focused on identifying low-permeability rocks suitable for testing large-yield nuclear devices, whereas later hydrologic investigations focused on potential subsurface transport of radionuclides away from Pahute Mesa. The number of single-well aquifer tests analyzed for Pahute Mesa and vicinity was much greater than the number of actual tests because each response to a water-level displacement was interpreted multiple times using different analytical models by different investigators, which were reported as distinct, individual aquifer tests. Some aquifer tests also were reinterpreted multiple times because the wells had not been defined uniquely. Previous databases and reports identified wells by different borehole names, reported various depths to top and bottom of open intervals, and used inconsistent wetted-aquifer thicknesses, which were coupled with inappropriate analytical methods to interpret aquifer-test data. This work standardized the reporting of single-well aquifer-test results at Pahute Mesa and vicinity. A hydraulic-parameter database was created that uniquely identifies tested wells, reports all replicated aquifer tests in each well, and identifies the best transmissivity estimate for each well. Datasets, interpreted models, and results from 1,459 analyses are presented and evaluated for 360 unique wells. Integrated borehole analyses were used to estimate hydraulic conductivity by depth in Pahute Mesa because vertical variations in hydraulic conductivity greatly affect groundwater velocities, which directly affect subsurface transport rates and directions. The integrated analyses reconciled differences between transmissivity estimated from pumping tests and the summed transmissivity from straddle-packer slug tests done at different depth intervals in a borehole. An integration borehole database was created that presents integrated aquifer-test results from 17 boreholes.

Nevada

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

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 percent 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.

Ground Water

A process to estimate net infiltration using a site-scale water-budget approach, Rainier Mesa, Nevada National Security Site, Nevada, 2002–05

This report documents a process used to estimate net infiltration from precipitation, evapotranspiration (ET), and soil data acquired at two sites on Rainier Mesa. Rainier Mesa is a groundwater recharge area within the Nevada National Security Site where recharged water flows through bedrock fractures to a deep (450 meters) water table. The U.S. Geological Survey operated two ET stations on Rainier Mesa from 2002 to 2005 at sites characterized by pinyon-juniper and scrub-brush vegetative cover. Precipitation and ET data were corrected to remove measurement biases and gap-filled to develop continuous datasets. Net infiltration (percolation below the root zone) and changes in root-zone water storage were estimated using a monthly water-balance model. Site-scale water-budget results indicate that the heavily-fractured welded-tuff bedrock underlying thin (<40 centimeters) topsoil is a critical water source for vegetation during dry periods. Annual precipitation during the study period ranged from fourth lowest (182 millimeters [mm]) to second highest (708 mm) on record (record = 55 years). Annual ET exceeded precipitation during dry years, indicating that the fractured-bedrock reservoir capacity is sufficient to meet atmospheric-evaporative demands and to sustain vegetation through extended dry periods. Net infiltration (82 mm) was simulated during the wet year after the reservoir was rapidly filled to capacity. These results support previous conclusions that preferential fracture flow was induced, resulting in an episodic recharge pulse that was detected in nearby monitoring wells. The occurrence of net infiltration only during the wet year is consistent with detections of water-level rises in nearby monitoring wells that occur only following wet years.

Nevada

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.

Groundwater

Hydraulic characterization of volcanic rocks in Pahute Mesa using an integrated analysis of 16 multiple-well aquifer tests, Nevada National Security Site, 2009–14

An improved understanding of groundwater flow and radionuclide migration downgradient from underground nuclear-testing areas at Pahute Mesa, Nevada National Security Site, requires accurate subsurface hydraulic characterization. To improve conceptual models of flow and transport in the complex hydrogeologic system beneath Pahute Mesa, the U.S. Geological Survey characterized bulk hydraulic properties of volcanic rocks using an integrated analysis of 16 multiple-well aquifer tests. Single-well aquifer-test analyses provided transmissivity estimates at pumped wells. Transmissivity estimates ranged from less than 1 to about 100,000 square feet per day in Pahute Mesa and the vicinity. Drawdown from multiple-well aquifer testing was estimated and distinguished from natural fluctuations in more than 200 pumping and observation wells using analytical water-level models. Drawdown was detected at distances greater than 3 miles from pumping wells and propagated across hydrostratigraphic units and major structures, indicating that neither faults nor structural blocks noticeably impede or divert groundwater flow in the study area. Consistent hydraulic properties were estimated by simultaneously interpreting drawdown from the 16 multiple-well aquifer tests with an integrated groundwater-flow model composed of 11 well-site models—1 for each aquifer test site. Hydraulic properties were distributed across volcanic rocks with the Phase II Pahute Mesa-Oasis Valley Hydrostratigraphic Framework Model. Estimated hydraulic-conductivity distributions spanned more than two orders of magnitude in hydrostratigraphic units. Overlapping hydraulic conductivity ranges among units indicated that most Phase II Hydrostratigraphic Framework Model units were not hydraulically distinct. Simulated total transmissivity ranged from 1,600 to 68,000 square feet per day for all pumping wells analyzed. High-transmissivity zones exceeding 10,000 square feet per day exist near caldera margins and extend along the northern and eastern Pahute Mesa study area and near the southwestern edge of the study area. The estimated hydraulic-property distributions and observed hydraulic connections among geologic structures improved the characterization and representation of groundwater flow at Pahute Mesa.

Nevada

Testing the suitability of geologic frameworks for extrapolating hydraulic properties across regional scales

The suitability of geologic frameworks for extrapolating hydraulic conductivity ( K ) to length scales commensurate with hydraulic data is difficult to assess. A novel method is presented for evaluating assumed relations between K and geologic interpretations for regional-scale groundwater modeling. The approach relies on simultaneous interpretation of multiple aquifer tests using alternative geologic frameworks of variable complexity, where each framework is incorporated as prior information that assumes homogeneous K within each model unit. This approach is tested at Pahute Mesa within the Nevada National Security Site (USA), where observed drawdowns from eight aquifer tests in complex, highly faulted volcanic rocks provide the necessary hydraulic constraints. The investigated volume encompasses 40 mi 3 (167 km 3 ) where drawdowns traversed major fault structures and were detected more than 2 mi (3.2 km) from pumping wells. Complexity of the five frameworks assessed ranges from an undifferentiated mass of rock with a single unit to 14 distinct geologic units. Results show that only four geologic units can be justified as hydraulically unique for this location. The approach qualitatively evaluates the consistency of hydraulic property estimates within extents of investigation and effects of geologic frameworks on extrapolation. Distributions of transmissivity are similar within the investigated extents irrespective of the geologic framework. In contrast, the extrapolation of hydraulic properties beyond the volume investigated with interfering aquifer tests is strongly affected by the complexity of a given framework. Testing at Pahute Mesa illustrates how this method can be employed to determine the appropriate level of geologic complexity for large-scale groundwater modeling.

Nevada

T-COMP — A suite of programs for extracting transmissivity from MODFLOW models

Simulated transmissivities are constrained poorly by assigning permissible ranges of hydraulic conductivities from aquifer-test results to hydrogeologic units in groundwater-flow models. These wide ranges are derived from interpretations of many aquifer tests that are categorized by hydrogeologic unit. Uncertainty is added where contributing thicknesses differ between field estimates and numerical models. Wide ranges of hydraulic conductivities and discordant thicknesses result in simulated transmissivities that frequently are much greater than aquifer-test results. Multiple orders of magnitude differences frequently occur between simulated and observed transmissivities where observed transmissivities are less than 1,000 feet squared per day. Transmissivity observations from individual aquifer tests can constrain model calibration as head and flow observations do. This approach is superior to diluting aquifer-test results into generalized ranges of hydraulic conductivities. Observed and simulated transmissivities can be compared directly with T-COMP, a suite of three FORTRAN programs. Transmissivity observations require that simulated hydraulic conductivities and thicknesses in the volume investigated by an aquifer test be extracted and integrated into a simulated transmissivity. Transmissivities of MODFLOW model cells are sampled within the volume affected by an aquifer test as defined by a well-specific, radial-flow model of each aquifer test. Sampled transmissivities of model cells are averaged within a layer and summed across layers. Accuracy of the approach was tested with hypothetical, multiple-aquifer models where specified transmissivities ranged between 250 and 20,000 feet squared per day. More than 90 percent of simulated transmissivities were within a factor of 2 of specified transmissivities.

Nevada

Delineation of the Pahute Mesa–Oasis Valley groundwater basin, Nevada

This report delineates the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, where recharge occurs, moves downgradient, and discharges to Oasis Valley, Nevada. About 5,900 acre-feet of water discharges annually from Oasis Valley, an area of springs and seeps near the town of Beatty in southern Nevada. Radionuclides in groundwater beneath Pahute Mesa, an area of historical underground nuclear testing at the Nevada National Security Site, are believed to be migrating toward Oasis Valley. Delineating the boundary of the PMOV groundwater basin is necessary to adequately assess the potential for transport of radionuclides from Pahute Mesa to Oasis Valley. The PMOV contributing area is defined based on regional water-level contours, geologic controls, and knowledge of adjacent flow systems. The viability of this area as the contributing area to Oasis Valley and the absence of significant interbasin flow between the PMOV groundwater basin and adjacent basins are shown regionally and locally. Regional constraints on the location of the contributing area boundary and on the absence of interbasin groundwater flow are shown by balancing groundwater discharges in the PMOV groundwater basin and adjacent basins against available water from precipitation. Internal consistency for the delineated contributing area is shown by matching measured water levels, groundwater discharges, and transmissivities with simulated results from a single-layer, steady-state, groundwater-flow model. An alternative basin boundary extending farther north than the final boundary was rejected based on a poor chloride mass balance and a large imbalance in the northern area between preferred and simulated recharge.

California, Nevada

Evaluating connection of aquifers to springs and streams, Great Basin National Park and vicinity, Nevada

Federal agencies that oversee land management for much of the Snake Range in eastern Nevada, including the management of Great Basin National Park by the National Park Service, need to understand the potential extent of adverse effects to federally managed lands from nearby groundwater development. As a result, this study was developed (1) to attain a better understanding of aquifers controlling groundwater flow on the eastern side of the southern part of the Snake Range and their connection with aquifers in the valleys, (2) to evaluate the relation between surface water and groundwater along the piedmont slopes, (3) to evaluate sources for Big Springs and Rowland Spring, and (4) to assess groundwater flow from southern Spring Valley into northern Hamlin Valley. The study focused on two areas&mdash;the first, a northern area along the east side of Great Basin National Park that included Baker, Lehman, and Snake Creeks, and a second southern area that is the potential source area for Big Springs. Data collected specifically for this study included the following: (1) geologic field mapping; (2) drilling, testing, and water quality sampling from 7 test wells; (3) measuring discharge and water chemistry of selected creeks and springs; (4) measuring streambed hydraulic gradients and seepage rates from 18 shallow piezometers installed into the creeks; and (5) monitoring stream temperature along selected reaches to identify places of groundwater inflow. The Snake Range was formed by a generally normal-faulted uplift, where late Proterozoic and Cambrian siliciclastic rocks and metamorphic rocks are present at the highest altitudes and younger Paleozoic carbonate rocks are exposed along the flanks. The consolidated rocks are intruded by Jurassic to Tertiary age plutons, which are most common between the Lehman and Snake Creek drainage basins. Older Cenozoic rocks, including Oligocene volcanic rocks and Miocene sedimentary rocks, crop out locally and fill the basins that underlie Snake, Spring, and Hamlin Valleys. Younger Tertiary and Quaternary sedimentary (basin-fill) deposits overlie the older Cenozoic rocks. The rocks and deposits can be divided into three distinct aquifers. These aquifers include (1) basin-fill aquifers that consist of the permeable parts of the Cenozoic basin fill and some fractured or jointed Cenozoic volcanic rocks, (2) an upper carbonate-rock aquifer that consists of upper Paleozoic carbonate rocks overlying a regionally extensive middle Paleozoic siliciclastic confining unit, and (3) a lower carbonate-rock aquifer that consists of lower Paleozoic carbonate rocks. Secondary openings created by faults, shear zones, fractures, and, in the carbonate rocks, karst solution features, largely determine the water-transmitting properties of the volcanic- and carbonate-rock aquifers. The basin-fill aquifers are composed of a wide variety of rock types and have highly variable hydraulic properties. The three aquifers are stratigraphically and structurally heterogeneous, causing large variations in the ability to store and transmit water. The aquifers are separated by confining units in some areas and are in contact with each other in other areas, yet function as a single, composite aquifer system. Basin-fill aquifers most often overlie or adjoin the lower and upper carbonate-rock aquifers. Baker, Lehman and Snake Creek drainage basins were divided into five hydrologic zones on the basis of climate, geology, and topography. The five zones, from highest to lowest altitudes, are the mountain-upland, karst-limestone, upper-piedmont, lower-piedmont, and valley-lowland zones. The primary hydrologic connection between the mountain-upland and the valley-lowland zones is streamflow. Much of the streamflow from the mountain-upland zone is generated above tree line. Groundwater flow increases in the karst-limestone zone because of increased permeability caused by dissolution, which results in increased streamflow losses. Most of the increased groundwater flow is to springs near faults that form the boundary with the upper-piedmont zone. Thus, groundwater flow from the karst-limestone zone to the upper-piedmont zone was only 10 percent of the combined flow of streams and springs that exit the karst-limestone zone. About 60 percent of the water flowing from Rowland Spring in the Lehman Creek drainage basin was from streamflow losses along Baker Creek. The remaining flow from Rowland Spring comes from local recharge in the karst-limestone zone. In the upper-piedmont zone, the water table by Baker, Lehman and Snake Creeks was near the water level in the creeks for several hundred feet downstream from the karst-limestone zone. Water levels in piezometers along Snake Creek downstream from its confluence with Spring Creek were far below the streambed, indicating gravity drainage beneath this section of the creek. Estimated vertical hydraulic conductivity along a 3-mile reach of Snake Creek downstream of this confluence was 0.5 foot per day, which was an order of magnitude less than that estimated for Baker and Lehman Creeks. The low vertical hydraulic conductivity in the streambed along the lower reaches of Snake Creek results from chemical precipitation of calcite caused by off-gassing of carbon dioxide derived from springs at the end of the karst-limestone zone. The younger alluvial deposits thicken rapidly across faults that form the upper boundary of the lower-piedmont zone. The absence of springs or groundwater flow to the creeks upstream of these faults indicates they are not a complete barrier to groundwater flow. The water table was shallow in the valley-lowland zone in the Baker and Lehman Creek drainage basins, whereas the water table was more than 50 feet below land surface in the Snake Creek drainage basin. In contrast to thick basin fill in the valley-lowland zone in the Baker and Lehman Creek drainage basins, fractured and karst limestone underlie basin fill at relatively shallow depths in Snake Creek drainage basin. The underlying limestone acts as a drain for groundwater in the basin fill beneath Snake Creek. A groundwater divide in southern Spring Valley south of Baking Powder Flat separates groundwater flow to the flat from southeastward flow into northern Hamlin Valley. Groundwater flow from southern Spring Valley south of the groundwater divide into northern Hamlin Valley was estimated to range from 6,000 to 11,000 acre-feet per year. This groundwater does not flow to Big Springs in southern Snake Valley; rather, the source of water to Big Springs is groundwater recharge in the Big Spring Wash drainage basin and in nearby smaller drainage basins at the south end of the Snake Range. Groundwater flow from southern Spring Valley continues through the western side of Hamlin Valley before being directed northeast toward the south end of Snake Valley. This flow is constrained by southward-flowing groundwater from Big Spring Wash and northward-flowing groundwater beneath central Hamlin Valley. The redirection to the northeast corresponds to a narrowing of the width of flow in southern Snake Valley caused by a constriction formed by a steeply dipping middle Paleozoic siliciclastic confining unit exposed in the flanks of the mountains and hills on the east side of southern Snake Valley and shallowly buried beneath basin fill in the valley. The narrowing of groundwater flow could be responsible for the large area where groundwater flows to springs or is lost to evapotranspiration between Big Springs in Nevada and Pruess Lake in Utah.

Nevada

Potential effects of groundwater pumping on water levels, phreatophytes, and spring discharges in Spring and Snake Valleys, White Pine County, Nevada, and adjacent areas in Nevada and Utah

Assessing hydrologic effects of developing groundwater supplies in Snake Valley required numerical, groundwater-flow models to estimate the timing and magnitude of capture from streams, springs, wetlands, and phreatophytes. Estimating general water-table decline also required groundwater simulation. The hydraulic conductivity of basin fill and transmissivity of basement-rock distributions in Spring and Snake Valleys were refined by calibrating a steady state, three-dimensional, MODFLOW model of the carbonate-rock province to predevelopment conditions. Hydraulic properties and boundary conditions were defined primarily from the Regional Aquifer-System Analysis (RASA) model except in Spring and Snake Valleys. This locally refined model was referred to as the Great Basin National Park calibration (GBNP-C) model. Groundwater discharges from phreatophyte areas and springs in Spring and Snake Valleys were simulated as specified discharges in the GBNP-C model. These discharges equaled mapped rates and measured discharges, respectively. Recharge, hydraulic conductivity, and transmissivity were distributed throughout Spring and Snake Valleys with pilot points and interpolated to model cells with kriging in geologically similar areas. Transmissivity of the basement rocks was estimated because thickness is correlated poorly with transmissivity. Transmissivity estimates were constrained by aquifer-test results in basin-fill and carbonate-rock aquifers. Recharge, hydraulic conductivity, and transmissivity distributions of the GBNP-C model were estimated by minimizing a weighted composite, sum-of-squares objective function that included measurement and Tikhonov regularization observations. Tikhonov regularization observations were equations that defined preferred relations between the pilot points. Measured water levels, water levels that were simulated with RASA, depth-to-water beneath distributed groundwater and spring discharges, land-surface altitudes, spring discharge at Fish Springs, and changes in discharge on selected creek reaches were measurement observations. The effects of uncertain distributed groundwater-discharge estimates in Spring and Snake Valleys on transmissivity estimates were bounded with alternative models. Annual distributed groundwater discharges from Spring and Snake Valleys in the alternative models totaled 151,000 and 227,000 acre-feet, respectively and represented 20 percent differences from the 187,000 acre-feet per year that discharges from the GBNP-C model. Transmissivity estimates in the basin fill between Baker and Big Springs changed less than 50 percent between the two alternative models. Potential effects of pumping from Snake Valley were estimated with the Great Basin National Park predictive (GBNP-P) model, which is a transient groundwater-flow model. The hydraulic conductivity of basin fill and transmissivity of basement rock were the GBNP-C model distributions. Specific yields were defined from aquifer tests. Captures of distributed groundwater and spring discharges were simulated in the GBNP-P model using a combination of well and drain packages in MODFLOW. Simulated groundwater captures could not exceed measured groundwater-discharge rates. Four groundwater-development scenarios were investigated where total annual withdrawals ranged from 10,000 to 50,000 acre-feet during a 200-year pumping period. Four additional scenarios also were simulated that added the effects of existing pumping in Snake Valley. Potential groundwater pumping locations were limited to nine proposed points of diversion. Results are presented as maps of groundwater capture and drawdown, time series of drawdowns and discharges from selected wells, and time series of discharge reductions from selected springs and control volumes. Simulated drawdown propagation was attenuated where groundwater discharge could be captured. General patterns of groundwater capture and water-table declines were similar for all scenarios. Simulated drawdowns greater than 1 ft propagated outside of Spring and Snake Valleys after 200 years of pumping in all scenarios.

Scientific Investigations Report

The source, discharge, and chemical characteristics of water from Agua Caliente Spring, Palm Springs, California

Agua Caliente Spring, in downtown Palm Springs, California, has been used for recreation and medicinal therapy for hundreds of years and currently (2008) is the source of hot water for the Spa Resort owned by the Agua Caliente Band of the Cahuilla Indians. The Agua Caliente Spring is located about 1,500 feet east of the eastern front of the San Jacinto Mountains on the southeast-sloping alluvial plain of the Coachella Valley. The objectives of this study were to (1) define the geologic structure associated with the Agua Caliente Spring; (2) define the source(s), and possibly the age(s), of water discharged by the spring; (3) ascertain the seasonal and longer-term variability of the natural discharge, water temperature, and chemical characteristics of the spring water; (4) evaluate whether water-level declines in the regional aquifer will influence the temperature of the spring discharge; and, (5) estimate the quantity of spring water that leaks out of the water-collector tank at the spring orifice.

California

Assessing hydraulic connections across a complex sequence of volcanic rocks - Analysis of U-20 WW multiple-well aquifer test, Pahute Mesa, Nevada National Security Site, Nevada

Groundwater beneath Pahute Mesa flows through a complexly layered sequence of volcanic rock aquifers and confining units that have been faulted into distinct structural blocks. Hydraulic property estimates of rocks and structures in this flow system are necessary to assess radionuclide migration near underground nuclear testing areas. The U.S. Geological Survey (USGS) used a 12 month (October 1, 2008&mdash; October 1, 2009) intermittent pumping schedule of well U-20 WW and continuously monitored water levels in observation wells ER-20-6 #3, UE-20bh 1, and U-20bg as a multi-well aquifer test to evaluate hydraulic connections across structural blocks, bulk hydraulic properties of volcanic rocks, and the hydraulic significance of a major fault. Measured water levels were approximated using synthetic water levels generated from an analytical model. Synthetic water levels are a summation of environmental water-level fluctuations and a Theis (1935) transform of the pumping signal from flow rate to water-level change. Drawdown was estimated by summing residual differences between measured and synthetic water levels and the Theis-transformed pumping signal from April to September 2009. Drawdown estimates were used in a three‑dimensional numerical model to estimate hydraulic properties of distinct aquifers, confining units, and a major fault. A maximum water-level drawdown of nearly 0.4 foot in well UE-20bh 1, which is more than 1 mile from the pumping well, was detected across a major fault. Drawdown estimates in the observation well nearest to (ER-20-6 #3, less than 1 mile) and within the same structural block as the pumping well were less than detection (<0.1 foot). Minimal drawdown within the same structural block indicates that lava units are likely separated by bedded tuff confining units. Hydraulic property estimates indicate that wells U-20 WW, UE-20bh 1, and ER-20-6 #3 produce water from moderately permeable fractured lava, as hydraulic conductivity and specific storage estimates average 4.8 feet per day and 2.1&times;10 &ndash;6 per foot, respectively, and transmissivity estimates range from 1,200 to 3,600 feet squared per day. Sensitivity analyses indicate that the major fault is hydraulically similar to the permeable host rock and connects flow between structural blocks.

Nevada

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 (mu or u)g/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 (mu or u)g/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 hydraulicconductivity 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 (mu or u)g/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 (mu or u)g/L over a 20-year period.

Ground Water

Interpretation of Flow Logs from Nevada Test Site Boreholes to Estimate Hydraulic Conductivity Using Numerical Simulations Constrained by Single-Well Aquifer Tests

Hydraulic conductivities of volcanic and carbonate lithologic units at the Nevada Test Site were estimated from flow logs and aquifer-test data. Borehole flow and drawdown were integrated and interpreted using a radial, axisymmetric flow model, AnalyzeHOLE. This integrated approach is used because complex well completions and heterogeneous aquifers and confining units produce vertical flow in the annular space and aquifers adjacent to the wellbore. AnalyzeHOLE simulates vertical flow, in addition to horizontal flow, which accounts for converging flow toward screen ends and diverging flow toward transmissive intervals. Simulated aquifers and confining units uniformly are subdivided by depth into intervals in which the hydraulic conductivity is estimated with the Parameter ESTimation (PEST) software. Between 50 and 150 hydraulic-conductivity parameters were estimated by minimizing weighted differences between simulated and measured flow and drawdown. Transmissivity estimates from single-well or multiple-well aquifer tests were used to constrain estimates of hydraulic conductivity. The distribution of hydraulic conductivity within each lithology had a minimum variance because estimates were constrained with Tikhonov regularization. AnalyzeHOLE simulated hydraulic-conductivity estimates for lithologic units across screened and cased intervals are as much as 100 times less than those estimated using proportional flow-log analyses applied across screened intervals only. Smaller estimates of hydraulic conductivity for individual lithologic units are simulated because sections of the unit behind cased intervals of the wellbore are not assumed to be impermeable, and therefore, can contribute flow to the wellbore. Simulated hydraulic-conductivity estimates vary by more than three orders of magnitude across a lithologic unit, indicating a high degree of heterogeneity in volcanic and carbonate-rock units. The higher water transmitting potential of carbonate-rock units relative to volcanic-rock units is exemplified by the large difference in their estimated maximum hydraulic conductivity; 4,000 and 400 feet per day, respectively. Simulated minimum estimates of hydraulic conductivity are inexact and represent the lower detection limit of the method. Minimum thicknesses of lithologic intervals also were defined for comparing AnalyzeHOLE results to hydraulic properties in regional ground-water flow models.

Scientific Investigations Report