South Lahontan

Death Valley regional ground-water flow system, Nevada and California—Hydrogeologic framework and transient ground-water flow model

U.S. Geological Survey (USGS) | July 1st, 2010

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 Ener

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.

Death Valley regional groundwater flow system, Nevada and California—Hydrogeologic framework and transient ground-water flow model

U.S. Geological Survey (USGS) | August 23rd, 2010

Decadal Shifts in Groundwater Age Detected by Environmental Tracers Across California

American Geophysical Union (AGU) | March 23rd, 2026

Groundwater age tells us how long water has been underground since it last infiltrated the land surface. Monitoring changes in groundwater age helps identify aquifer stre

Groundwater age tells us how long water has been underground since it last infiltrated the land surface. Monitoring changes in groundwater age helps identify aquifer stress from pumping or shifts toward new hydrologic conditions—key factors for sustainable water management. In this study, changes in groundwater age were tested by comparing chemical tracer data collected about 10 years apart from hundreds of wells across California. Chemical tracers are dissolved substances with known historical records and decay rates that can be used to infer groundwater age. Using models that estimate expected tracer concentrations under steady-age conditions, this study found that nearly one-third of the wells showed significant changes in groundwater age over time. These changes were most common in areas with heavy agricultural use, where pumping and irrigation have altered how water moves underground. These shifts in age can also affect water quality, especially when younger water carries contaminants like nitrate. By using a statistical modeling approach, this study shows how changes in groundwater age can help identify aquifers under stress and support efforts to manage groundwater sustainably.

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–2023, California groundwater ambient monitoring and assessment program, priority basin project

U.S. Geological Survey (USGS) | February 25th, 2026

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring site

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–2023. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban- and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert-butyl ether (MTBE) showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids (TDS) statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in TDS by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and MTBE are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

Decline of the world's saline lakes

Nature Portfolio (Springer Nature) | October 23rd, 2017

Many of the world's saline lakes are shrinking at alarming rates, reducing waterbird habitat and economic benefits while threatening human health. Saline lakes are long-t

Many of the world's saline lakes are shrinking at alarming rates, reducing waterbird habitat and economic benefits while threatening human health. Saline lakes are long-term basin-wide integrators of climatic conditions that shrink and grow with natural climatic variation. In contrast, water withdrawals for human use exert a sustained reduction in lake inflows and levels. Quantifying the relative contributions of natural variability and human impacts to lake inflows is needed to preserve these lakes. With a credible water balance, causes of lake decline from water diversions or climate variability can be identified and the inflow needed to maintain lake health can be defined. Without a water balance, natural variability can be an excuse for inaction. Here we describe the decline of several of the world's large saline lakes and use a water balance for Great Salt Lake (USA) to demonstrate that consumptive water use rather than long-term climate change has greatly reduced its size. The inflow needed to maintain bird habitat, support lake-related industries and prevent dust storms that threaten human health and agriculture can be identified and provides the information to evaluate the difficult tradeoffs between direct benefits of consumptive water use and ecosystem services provided by saline lakes.

Delta Chromium-53/52 isotopic composition of native and contaminated groundwater, Mojave Desert, USA

Applied Geochemistry (Elsevier) | April 20th, 2012

Chromium(VI) concentrations in groundwater sampled from three contaminant plumes in aquifers in the Mojave Desert near Hinkley, Topock and El Mirage, California, USA, wer

Chromium(VI) concentrations in groundwater sampled from three contaminant plumes in aquifers in the Mojave Desert near Hinkley, Topock and El Mirage, California, USA, were as high as 2600, 5800 and 330 μg/L, respectively. δ53/52Cr compositions from more than 50 samples collected within these plumes ranged from near 0‰ to almost 4‰ near the plume margins. Assuming only reductive fractionation of Cr(VI) to Cr(III) within the plume, apparent fractionation factors for δ53/52Cr isotopes ranged from εapp = 0.3 to 0.4 within the Hinkley and Topock plumes, respectively, and only the El Mirage plume had a fractionation factor similar to the laboratory derived value of ε = 3.5. One possible explanation for the difference between field and laboratory fractionation factors at the Hinkley and Topock sites is localized reductive fractionation of Cr(VI) to Cr(III), with subsequent advective mixing of native and contaminated water near the plume margin. Chromium(VI) concentrations and δ53/52Cr isotopic compositions did not uniquely define the source of Cr near the plume margin, or the extent of reductive fractionation within the plume. However, Cr(VI) and δ53/52Cr data contribute to understanding of the interaction between reductive and mixing processes that occur within and near the margins of Cr contamination plumes. Reductive fractionation of Cr(VI) predominates in plumes having higher εapp, these plumes may be suitable for monitored natural attenuation. In contrast, advective mixing predominates in plumes having lower εapp, the highly dispersed margins of these plumes may be difficult to define and manage.

Design and methods of the California stream quality assessment (CSQA), 2017

U.S. Geological Survey (USGS) | April 21st, 2020

During 2017, as part of the National Water-Quality Assessment Project, the U.S. Geological Survey conducted the California Stream Quality Assessment to investigate the qu

During 2017, as part of the National Water-Quality Assessment Project, the U.S. Geological Survey conducted the California Stream Quality Assessment to investigate the quality of streams in the Central California Foothills and Coastal Mountains ecoregion, United States. The goal of the California Stream Quality Assessment study was to assess the health of wadeable streams in the region by characterizing multiple water-quality factors that are stressors to aquatic biota and by evaluating the relation between these stressors and biological indicators of stream health. Urbanization, agriculture, and modifications to streamflow are anthropogenic changes that affect water quality in the region; consequently, the study design primarily targeted sites and specific stressors associated with these activities. For the study, 85 stream sites were selected to represent the types and intensity of land use in the watershed; categories of site types were undeveloped, urban (low, medium, high), agriculture (low, high), and mixed (urban and agriculture). Most sites (about 70 percent) represent a gradient of urbanization from undeveloped to 99-percent urbanized. At most of the sites, streamgages or pressure transducers were used to monitor stream discharge and stage, as well as temperature. Water-quality samples were collected routinely at all sites and were analyzed for major ions, organic contaminants, nutrients, and suspended sediment. Sampling frequency varied on the basis of site type and location. Discrete water samples were collected weekly and generally 6 times per site, except for 11 undeveloped sites that were sampled only 4 times (during the last 4 weeks). Water sampling began at sites in the southern part of the study on March 13, 2017, and at sites in the northern part of the study on April 3, 2017. Passive samplers were deployed at most sites for measurement of polar organic contaminants (pesticides and pharmaceuticals). In May 2017, coincident with completion of water-quality sampling, an ecological survey was conducted at each site to assess benthic algal and macroinvertebrate communities and instream habitat. During the ecological surveys, a single composite streambed-sediment sample was collected for chemical analysis and toxicity testing. In addition, a few focused studies were done at subsets of sites, namely, measuring pesticides using small-volume automated samplers, measuring pesticides in biofilms, and sampling suspended sediments using passive samplers. This report describes the various study components and methods of the California Stream Quality Assessment, including measurements of water quality, sediment chemistry, habitat assessments, and ecological surveys, as well as procedures for sample analysis, quality assurance and quality control, and data management.

Detection and Measurement of Land Subsidence Using Interferometric Synthetic Aperture Radar and Global Positioning System, San Bernardino County, Mojave Desert, California

U.S. Geological Survey (USGS) | July 1st, 2003

Detection of aquifer system compaction and land subsidence using interferometric synthetic aperture radar, Antelope Valley, Mojave Desert, California

U.S. Geological Survey (USGS) | April 15th, 1998

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with h

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (10 mm in change of land surface elevation) were developed for a groundwater basin (103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

Development of a Mono Lake Water Level Model and Projections across the 21st Century

University of California, Los Angeles (UCLA) | March 17th, 2026

The UCLA Mono Lake Model (UCLA-MLM) was developed to assess how climate-driven changes and existing and alternative export criteria influence Mono Lake’s water level th

The UCLA Mono Lake Model (UCLA-MLM) was developed to assess how climate-driven changes and existing and alternative export criteria influence Mono Lake’s water level throughout the 21st century.