USGS Science⌕ Search

Geology topics

Research about Death Valley

Source-linked reports with geographic coverage including Death Valley.

At least 19 recordsLinked to original sources
USGS2022

Hydroclimate response of spring ecosystems to a two-stage Younger Dryas event in western North America

The Younger Dryas (YD) climate event is the preeminent example of abrupt climate change in the recent geologic past. Climate conditions during the YD were spatially complex, and high-resolution sediment cores in the North Atlantic, western Europe, and East Asia have revealed it unfolded in two distinct stages, including an initial stable climatic period between ~ 12.9 and 12.2 ka associated with a weakened Atlantic Meridional Overturning Circulation (AMOC) and a second phase characterized by variable conditions until 11.7 ka as the AMOC recovered. Decades of investigations into the climate of western North America during the YD have failed to identify this stepped phenomenon. Here we present hydroclimate data from paleospring deposits in Death Valley National Park (California, USA) that demonstrate unequivocal evidence of two-stage partitioning within the YD event. High groundwater levels supported persistent and long-lived spring ecosystems between ~ 13.0 and 12.2 ka, which were immediately replaced by alternating wet and dry environments until ~ 11.8 ka. These results establish the mid-YD climate transition extended into western North America at approximately the same time it was recorded by hydrologic systems elsewhere in the Northern Hemisphere and show that even short-lived changes in the AMOC can have profound consequences for ecosystems worldwide.

California, Nevada↗
USGS2019

Geology of the Monte Blanco borate deposits, Furnace Creek Wash, Death Valley, California

The Monte Blanco borate deposits are located along the southern margin of Death Valley’s Furnace Creek Wash, south of Twenty Mule Team Canyon road in California. Topographic and geologic mapping by S. Muessig and F.M. Byers, Jr., in 1954 documented these deposits’ geologic settings, geometries, mineralogies, and chemical characteristics. They estimated borate resources at the time to be in excess of 550,000 tons B 2 O 3 . The borate bodies are composed of predominantly ulexite and colemanite. They lie beneath Monte Blanco itself and along a northwest-trending series of conspicuous, white hills and mounds formed by northeasterly dipping, fine-grained sedimentary beds and basaltic volcanic rocks of the Miocene and Pliocene Furnace Creek Formation. Geologic data suggest that in Miocene and Pliocene time, fine-grained sediments, volcanic debris and flows, and volcanically associated, boron-rich fluids gradually filled a fairly flat playa-like environment. At times, thick beds of felty crystals of ulexite developed and were interlayered as lenses in a thick series of mudstones as is seen today at the Eagle Borax works. After burial, the exterior of the ulexite deposit was altered to massive colemanite by ground water, which produced the “shell” of colemanite that typically surrounds the presently outcropping ulexite bodies.

California↗
USGS2018

Monogenetic origin of Ubehebe Crater maar volcano, Death Valley, California: Paleomagnetic and stratigraphic evidence

Paleomagnetic data for samples collected from outcrops of basaltic spatter at the Ubehebe Crater cluster, Death Valley National Park, California, record a single direction of remanent magnetization indicating that these materials were emplaced during a short duration, monogenetic eruption sequence ~ 2100 years ago. This conclusion is supported by geochemical data encompassing a narrow range of oxide variation, by detailed stratigraphic studies of conformable phreatomagmatic tephra deposits showing no evidence of erosion between layers, by draping of sharp rimmed craters by later tephra falls, and by oxidation of later tephra layers by the remaining heat of earlier spatter. This model is also supported through a reinterpretation and recalculation of the published age results (Sasnett et al., 2012) from an innovative and bold exposure-age study on very young materials. Their conclusion of multiple and protracted eruptions at Ubehebe Crater cluster is here modified through the understanding that some of their quartz-bearing clasts inherited from previous exposure on the fan surface (too old), and that other clasts were only exposed at the surface by wind and/or water erosion centuries after their eruption (too young). Ubehebe Crater cluster is a well preserved example of young monogenetic maar type volcanism protected within a National Park, and it represents neither a protracted eruption sequence as previously thought, nor a continuing volcanic hazard near its location.

California↗
USGS2016

Geology of the Greenwater Range, and the dawn of Death Valley, California—Field guide for the Death Valley Natural History Conference, 2013

Much has been written about the age and formation of Death Valley, but that is one—if not the last—chapter in the fascinating geologic history of this area. Igneous and sedimentary rocks in the Greenwater Range, one mountain range east of Death Valley, tell an earlier story that overlaps with the formation of Death Valley proper. This early story has been told by scientists who have studied these rocks for many years and continue to do so. This field guide was prepared for the first Death Valley Natural History Conference and provides an overview of the geology of the Greenwater Range and the early history (10–0 Ma) of Death Valley.

California, Nevada↗
USGS2013

Knowledge, transparency, and refutability in groundwater models, an example from the Death Valley regional groundwater flow system

This work demonstrates how available knowledge can be used to build more transparent and refutable computer models of groundwater systems. The Death Valley regional groundwater flow system, which surrounds a proposed site for a high level nuclear waste repository of the United States of America, and the Nevada National Security Site (NNSS), where nuclear weapons were tested, is used to explore model adequacy, identify parameters important to (and informed by) observations, and identify existing old and potential new observations important to predictions. Model development is pursued using a set of fundamental questions addressed with carefully designed metrics. Critical methods include using a hydrogeologic model, managing model nonlinearity by designing models that are robust while maintaining realism, using error-based weighting to combine disparate types of data, and identifying important and unimportant parameters and observations and optimizing parameter values with computationally frugal schemes. The frugal schemes employed in this study require relatively few (10–1000 s), parallelizable model runs. This is beneficial because models able to approximate the complex site geology defensibly tend to have high computational cost. The issue of model defensibility is particularly important given the contentious political issues involved.

California, Nevada↗
USGS2012

Modeling the Death Valley regional ground-water flow system

The development of a regional ground-water flow model of the Death Valley region in the southwestern United States is discussed in the context of the fourteen guidelines of Hill. This application of the guidelines demonstrates how they may be used for model calibration and evaluation, and to direct further model development and data collection.

California, Nevada↗
USGS2011

Comment on “An unconfined groundwater model of the Death Valley Regional Flow System and a comparison to its confined predecessor” by R.W.H. Carroll, G.M. Pohll and R.L. Hershey [Journal of Hydrology 373/3–4, pp. 316–328]

Carroll et al. (2009) state that the United States Geological Survey (USGS) Death Valley Regional Flow System (DVRFS) model, which is based on MODFLOW, is “conceptually inaccurate in that it models an unconfined aquifer as a confined system and does not simulate unconfined drawdown in transient pumping simulations.” Carroll et al. (2009) claim that “more realistic estimates of water availability” can be produced by a SURFACT-based model of the DVRFS that simulates unconfined groundwater flow and limits withdrawals from wells to avoid excessive drawdown. Differences in results from the original MODFLOW-based model and the SURFACT-based model stem primarily from application by Carroll et al. (2009) of head limits that can also be applied using the existing MODLOW model and not from any substantial difference in the accuracy with which the unconfined aquifer is represented in the two models. In a hypothetical 50-year predictive simulation presented by Carroll et al. (2009), large differences between the models are shown when simulating pumping from the lower clastic confining unit, where the transmissivity is nearly two orders of magnitude less than in an alluvial aquifer. Yet even for this extreme example, drawdowns and pumping rates from the MODFLOW and SURFACT models are similar when the head-limit capabilities of the MODFLOW MNW Package are applied. These similarities persist despite possible discrepancies between assigned hydraulic properties. The resulting comparison between the MODFLOW and SURFACT models of the DVRFS suggests that approximating the unconfined system in the DVRFS as a constant-saturated-thickness system (called a “confined system” by Carroll et al., 2009) performs very well.

Death Valley↗
USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗
USGS2003

Methods for using groundwater model predictions to guide hydrogeologic data collection, with application to the Death Valley regional groundwater flow system

Calibrated models of groundwater systems can provide substantial information for guiding data collection. This work considers using such models to guide hydrogeologic data collection for improving model predictions by identifying model parameters that are most important to the predictions. Identification of these important parameters can help guide collection of field data about parameter values and associated flow system features and can lead to improved predictions. Methods for identifying parameters important to predictions include prediction scaled sensitivities (PSS), which account for uncertainty on individual parameters as well as prediction sensitivity to parameters, and a new “value of improved information” (VOII) method presented here, which includes the effects of parameter correlation in addition to individual parameter uncertainty and prediction sensitivity. In this work, the PSS and VOII methods are demonstrated and evaluated using a model of the Death Valley regional groundwater flow system. The predictions of interest are advective transport paths originating at sites of past underground nuclear testing. Results show that for two paths evaluated the most important parameters include a subset of five or six of the 23 defined model parameters. Some of the parameters identified as most important are associated with flow system attributes that do not lie in the immediate vicinity of the paths. Results also indicate that the PSS and VOII methods can identify different important parameters. Because the methods emphasize somewhat different criteria for parameter importance, it is suggested that parameters identified by both methods be carefully considered in subsequent data collection efforts aimed at improving model predictions.

California, Nevada↗
USGS2003

Estimated ground-water withdrawals from the Death Valley regional flow system, Nevada and California, 1913-98

Ground-water withdrawals from 1913 through 1998 from the Death Valley regional flow system have been compiled to support a regional, three-dimensional, transient ground-water flow model. Withdrawal locations and depths of production intervals were estimated and associated errors were reported for 9,300 wells. Withdrawals were grouped into three categories: mining, public-supply, and commercial water use; domestic water use; and irrigation water use. In this report, groupings were based on the method used to estimate pumpage. Cumulative ground-water withdrawals from 1913 through 1998 totaled 3 million acre-feet, most of which was used to irrigate alfalfa. Annual withdrawal for irrigation ranged from 80 to almost 100 percent of the total pumpage. About 75,000 acre-feet was withdrawn for irrigation in 1998. Annual irrigation withdrawals generally were estimated as the product of irrigated acreage and application rate. About 320 fields totaling 11,000 acres were identified in six hydrographic areas. Annual application rates for high water-use crops ranged from 5 feet in Penoyer Valley to 9 feet in Pahrump Valley. The uncertainty in the estimates of ground-water withdrawals was attributed primarily to the uncertainty of application rate estimates. Annual ground-water withdrawal was estimated at about 90,000 acre-feet in 1998 with an assigned uncertainty bounded by 60,000 to 130,000 acre-feet.

California, Nevada↗
USGS2002

Hydrostructural maps of the Death Valley regional flow system, Nevada and California

The locations of principal faults and structural zones that may influence ground-water flow were compiled in support of a three-dimensional ground-water model for the Death Valley regional flow system (DVRFS), which covers 80,000 square km in southwestern Nevada and southeastern California. Faults include Neogene extensional and strike-slip faults and pre-Tertiary thrust faults. Emphasis was given to characteristics of faults and deformed zones that may have a high potential for influencing hydraulic conductivity. These include: (1) faulting that results in the juxtaposition of stratigraphic units with contrasting hydrologic properties, which may cause ground-water discharge and other perturbations in the flow system; (2) special physical characteristics of the fault zones, such as brecciation and fracturing, that may cause specific parts of the zone to act either as conduits or as barriers to fluid flow; (3) the presence of a variety of lithologies whose physical and deformational characteristics may serve to impede or enhance flow in fault zones; (4) orientation of a fault with respect to the present-day stress field, possibly influencing hydraulic conductivity along the fault zone; and (5) faults that have been active in late Pleistocene or Holocene time and areas of contemporary seismicity, which may be associated with enhanced permeabilities. The faults shown on maps A and B are largely from Workman and others (in press), and fit one or more of the following criteria: (1) faults that are more than 10 km in map length; (2) faults with more than 500 m of displacement; and (3) faults in sets that define a significant structural fabric that characterizes a particular domain of the DVRFS. The following fault types are shown: Neogene normal, Neogene strike-slip, Neogene low-angle normal, pre-Tertiary thrust, and structural boundaries of Miocene calderas. We have highlighted faults that have late Pleistocene to Holocene displacement (Piety, 1996). Areas of thick Neogene basin-fill deposits (thicknesses 1-2 km, 2-3 km, and >3 km) are shown on map A, based on gravity anomalies and depth-to-basement modeling by Blakely and others (1999). We have interpreted the positions of faults in the subsurface, generally following the interpretations of Blakely and others (1999). Where geophysical constraints are not present, the faults beneath late Tertiary and Quaternary cover have been extended based on geologic reasoning. Nearly all of these concealed faults are shown with continuous solid lines on maps A and B, in order to provide continuous structures for incorporation into the hydrogeologic framework model (HFM). Map A also shows the potentiometric surface, regional springs (25-35 degrees Celsius, D'Agnese and others, 1997), and cold springs (Turner and others, 1996).

Nevada;California↗
USGS2002

Three-dimensional hydrogeologic framework model for use with a steady-state numerical ground-water flow model of the Death Valley regional flow system, Nevada and California

The U.S. Geological Survey, in cooperation with the Department of Energy and other Federal, State, and local agencies, is evaluating the hydrogeologic characteristics of the Death Valley regional ground-water flow system. The ground-water flow system covers an area of about 100,000 square kilometers from latitude 35° to 38°15' North to longitude 115° to 118° West, with the flow system proper comprising about 45,000 square kilometers. The Death Valley regional ground-water flow system is one of the larger flow systems within the Southwestern United States and includes in its boundaries the Nevada Test Site, Yucca Mountain, and much of Death Valley. Part of this study includes the construction of a three-dimensional hydrogeologic framework model to serve as the foundation for the development of a steady-state regional ground-water flow model. The digital framework model provides a computer-based description of the geometry and composition of the hydrogeologic units that control regional flow. The framework model of the region was constructed by merging two previous framework models constructed for the Yucca Mountain Project and the Environmental Restoration Program Underground Test Area studies at the Nevada Test Site. The hydrologic characteristics of the region result from a currently arid climate and complex geology. Interbasinal regional ground-water flow occurs through a thick carbonate-rock sequence of Paleozoic age, a locally thick volcanic-rock sequence of Tertiary age, and basin-fill alluvium of Tertiary and Quaternary age. Throughout the system, deep and shallow ground-water flow may be controlled by extensive and pervasive regional and local faults and fractures. The framework model was constructed using data from several sources to define the geometry of the regional hydrogeologic units. These data sources include (1) a 1:250,000-scale hydrogeologic-map compilation of the region; (2) regional-scale geologic cross sections; (3) borehole information, and (4) gridded surfaces from a previous three-dimensional geologic model. In addition, digital elevation model data were used in conjunction with these data to define ground-surface altitudes. These data, properly oriented in three dimensions by using geographic information systems, were combined and gridded to produce the upper surfaces of the hydrogeologic units used in the flow model. The final geometry of the framework model is constructed as a volumetric model by incorporating the intersections of these gridded surfaces and by applying fault truncation rules to structural features from the geologic map and cross sections. The cells defining the geometry of the hydrogeologic framework model can be assigned several attributes such as lithology, hydrogeologic unit, thickness, and top and bottom altitudes.

California, Nevada↗
USGS2002

A three-dimensional numerical model of predevelopment conditions in the Death Valley regional ground-water flow system, Nevada and California

In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this "second-generation" regional model was to enhance the knowledge an understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-state representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration. The Death Valley regional ground-water flow system is situated within the southern Great Basin, a subprovince of the Basin and Range physiographic province, bounded by latitudes 35 degrees north and 38 degrees 15 minutes north and by longitudes 115 and 118 degrees west. Hydrology in the region is a result of both the arid climatic conditions and the complex geology. Ground-water flow generally can be described as dominated by interbasinal flow and may be conceptualized as having two main components: a series of relatively shallow and localized flow paths that are superimposed on deeper regional flow paths. A significant component of the regional ground-water flow is through a thick Paleozoic carbonate rock sequence. Throughout the flow system, ground water flows through zones of high transmissivity that have resulted from regional faulting and fracturing. The conceptual model of the Death Valley regional ground-water flow system used for this study is adapted from the two previous ground-water modeling studies. The three-dimensional digital hydrogeologic framework model developed for the region also contains elements of both of the hydrogeologic framework models used in the previous investigations. As dictated by project scope, very little reinterpretation and refinement were made where these two framework models disagree; therefore, limitations in the hydrogeologic representation of the flow system exist. Despite limitations, the framework model provides the best representation to date of the hydrogeologic units and structures that control regional ground-water flow and serves as an important information source used to construct and calibrate the predevelopment, steady-state flow model. In addition to the hydrogeologic framework, a complex array of mechanisms accounts for flow into, through, and out of the regional ground-water flow system. Natural discharges from the regional ground-water flow system occur by evapotranspiration, springs, and subsurface outflow. In this study, evapotranspiration rates were adapted from a related investigation that developed maps of evapotranspiration areas and computed rates from micrometeorological data collected within the local area over a multiyear period. In some cases, historical spring flow records were used to derive ground-water discharge rates for isolated regional springs. For this investigation, a process-based, numerical model was developed to estimate net infiltration. This result can be used to represent recharge by assuming that water that infiltrates past the “root zone” ultimately becomes ground-water recharge. Net infiltration based on this water-balance approach, however, suggests that recharge in the region exceeds the estimated ground-water discharge for the basin by almost a factor of four. The net infiltration model contains some assumptions and was limited in regard to soil moisture, bedrock permeability, and precipitation distributions that result in these higher than expected values. Despite this limitation, the net infiltration model does provide a good starting point for estimating regional ground-water recharge in the ground-water flow model. The Death Valley regional ground-water flow system was simulated using a three-dimensional steady-state model that incorporated a nonlinear least-squares regression technique to estimate selected model parameters. The numerical modeling program MODFLOW-2000 was used in creating a finite-difference model consisting of 194 rows, 160 columns, and 15 layers. The grid cells were oriented north-south and were of uniform size, with side dimensions of 1,500 meters. The layers span thicknesses of 50 to 300 meters. The model grid encompasses about 70,000 square kilometers. The required initial model parameter values were supplied by discretization of the three-dimensional hydrogeologic framework model and digital representations of the remaining conceptual model components. The three-dimensional simulation and corresponding sensitivity analysis supported the hypothesis of interactions between a relatively shallow local and subregional flow system and a deeper dominant regional system controlled by the carbonate aquifer. During calibration of the model, techniques available in MODFLOW-2000 allowed for estimation of a series of parameters that provided a best fit to observed hydraulic heads and flows. Numerous conceptual models were evaluated to test the validity of various interpretations about the flow system. Only those conceptual model changes contributing to a significant improvement in model fit, as indicated by a reduction in the sum of squared errors, were retained in the final optimized model. The final model was evaluated to assess the likely accuracy of simulated results by comparing measured and expected quantities with simulated values. Evaluation of the model indicates that although the model is clearly an improvement on previous representations of the flow system, there is an indication of important uncertainties and model error.

California, Nevada↗
USGS2002

Probability distributions of hydraulic conductivity for the hydrogeologic units of the Death Valley regional ground-water flow system, Nevada and California

The use of geologic information such as lithology and rock properties is important to constrain conceptual and numerical hydrogeologic models. This geologic information is difficult to apply explicitly to numerical modeling and analyses because it tends to be qualitative rather than quantitative. This study uses a compilation of hydraulic-conductivity measurements to derive estimates of the probability distributions for several hydrogeologic units within the Death Valley regional ground-water flow system, a geologically and hydrologically complex region underlain by basin-fill sediments, volcanic, intrusive, sedimentary, and metamorphic rocks. Probability distributions of hydraulic conductivity for general rock types have been studied previously; however, this study provides more detailed definition of hydrogeologic units based on lithostratigraphy, lithology, alteration, and fracturing and compares the probability distributions to the aquifer test data. Results suggest that these probability distributions can be used for studies involving, for example, numerical flow modeling, recharge, evapotranspiration, and rainfall runoff. These probability distributions can be used for such studies involving the hydrogeologic units in the region, as well as for similar rock types elsewhere. Within the study area, fracturing appears to have the greatest influence on the hydraulic conductivity of carbonate bedrock hydrogeologic units. Similar to earlier studies, we find that alteration and welding in the Tertiary volcanic rocks greatly influence hydraulic conductivity. As alteration increases, hydraulic conductivity tends to decrease. Increasing degrees of welding appears to increase hydraulic conductivity because welding increases the brittleness of the volcanic rocks, thus increasing the amount of fracturing.

California, Nevada↗
USGS2001

Preliminary evaluation of the importance of existing hydraulic-head observation locations to advective-transport predictions, Death Valley regional flow system, California and Nevada

When a model is calibrated by nonlinear regression, calculated diagnostic statistics and measures of uncertainty provide a wealth of information about many aspects of the system. This report presents a method of ranking the likely importance of existing observation locations using measures of prediction uncertainty. It is suggested that continued monitoring is warranted at more important locations, and unwarranted or less warranted at less important locations. The report develops the methodology and then demonstrates it using the hydraulic-head observation locations of a three-layer model of the Death Valley regional flow system. The predictions of interest are subsurface transport from beneath Yucca Mountain and 14 Underground Test Areas. The advective component of transport is considered because it is the component most affected by the system dynamics represented by the scale model being used. The problem is addressed using the capabilities of the U.S. Geological Survey computer program MODFLOW-2000, with its ADVective-Travel Observation (ADV) Package, and an additional computer program developed for this work. The methods presented in this report are used in three ways. (1) The ratings for individual observations are obtained by manipulating the measures of prediction uncertainty, and do not involve recalibrating the model. In this analysis, observation locations are each omitted individually and the resulting increase in uncertainty in the predictions is calculated. The uncertainty is quantified as standard deviations on the simulated advective transport. The increase in uncertainty is quantified as the percent increase in the standard deviations caused by omitting the one observation location from the calculation of standard deviations. In general, observation locations associated with larger increases are rated as more important. (2) Ratings for largely geographically based groups are obtained using a straightforward extension of the method used for individual observation locations. This analysis is needed where observations are clustered to determine whether the area is important to the predictions of interest. (3) Finally, the method is used to evaluate omitting a set of 100 observation locations. The locations were selected because they had low individual ratings and were not one of the few locations at which hydraulic heads from deep in the system were measured. The major results of the three analyses, when applied to the three-layer DVRFS ground-water flow system, are described in the following paragraphs. The discussion is labeled using the numbers 1 to 3 to clearly relate it to the three ways the method is used, as listed above. (1) The individual observation location analysis indicates that three observation locations are most important. They are located in Emigrant Valley, Oasis Valley, and Beatty. Of importance is that these and other observations shown to be important by this analysis are far from the travel paths considered. This displays the importance of the regional setting within which the transport occurs, the importance of including some sites throughout the area in the monitoring network, and the importance of including sites in these areas in particular. The method considered in this report indicates that the 19 observation locations that reflect hydraulic heads deeper in the system (in model layers 1, 2, and 3) are not very important. This appears to be because the locations of these observations are in the vicinity of shallow observation locations that also generally are rated as low importance, and because the model layers are hydraulically well connected vertically. The value of deep observations to testing conceptual models, however, is stressed. As a result, the deep observations are rated higher than is consistent with the results of the analysis presented, and none of these observations are omitted in the scenario discussed under (3) below. (2) The geographic grouping of the observations found one major area of importance not identified by the individual observation analysis. Five of the 49 groups are categorized as most highly important. The most important groups were those that, when omitted, produced mean increases greater than 10 percent at any UGTA site or Yucca Mountain. Four of the five groups were dominated by one individual observation. However, one group, located in Ash Meadows, had no individual observations ranked of high importance but collectively, when omitted, increased uncertainty substantially. Other groups also located in Ash Meadows, including intermediate depth observations, consistently ranked as more important than all other groups. (3) To demonstrate the importance of omitting a set of low-rated observations, one scenario is considered in which the 100 individually lowest-rated shallow and intermediate-depth observation locations are omitted. The measure of overall prediction uncertainty increased by just 0.59 percent, indicating that the wells associated with these observations probably could prudently be measured less frequently.

California, Nevada↗