USGS ScienceSearch

SEARCH · USGS Science

Results for “Geochemistry: Exploration, Environment, Analysis”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

633 records · Page 34Linked to original sources

Evaluating machine learning approaches to identify and predict oil and gas produced water lithium concentrations

Recently, the demand for battery-grade lithium has substantially increased, largely due to electrification of the transportation sector. The search for new lithium sources has turned to produced waters (frequently brines), a large-volume wastewater by-product of oil and gas extraction. Geochemical analysis indicates the presence of varying concentrations of lithium from produced water samples collected across the United States and represented in the U.S. Geological Survey’s National Produced Water Geochemical Database, as well as mixtures of Marcellus Shale produced water included in the Pennsylvania Department of Environmental Protection’s Oil and Gas Well Waste Reports. We first examined whether the geochemical signature of the lithium-bearing produced waters is sufficiently distinct so that machine learning (ML) can be used to correctly classify samples to the formation of origin. The produced water sample data used to assess classification accuracy were from the Marcellus Shale, Utica Shale and Point Pleasant Formation (Utica), and Smackover Formation oil and gas wells. Further, we evaluated the potential for ML to accurately classify Marcellus Shale produced water spatially (i.e., northeast versus southwest Pennsylvania). We then investigated whether ML algorithms applied to a suite of geochemical concentration data (i.e. Ba, Br, Cl, K, Mg, Sr) may be used to predict the lithium concentration of an unknown sample. Finally, we applied an estimated economic lithium grade cutoff of 150 milligrams per liter (mg/l) and assessed the utility of ML to predict whether a produced water sample would fall above or below the grade cutoff based on the suite of geochemical parameters. Four machine learning algorithms—Random Forest (RF), Gradient Boosting Trees (GBT), Extreme Boosting (XGBoost), and Deep Neural Networks (DNN) were assessed. This study successfully demonstrates that all four machine learning methods can precisely and accurately estimate lithium concentrations and geologic formation classification. The products of this study contribute to the growing body of knowledge aimed at expanding the lithium resource base within the United States.

Alabama, Arkansas, Florida, Georgia, Louisiana, Mi

Land cover change within wetland complexes at Dixie Meadows, Churchill County, Nevada: 2015 – 2023

Dixie Meadows, Nevada, is a system of geothermal springs and seeps that feed a complex of marshes and wetland meadows that are located within lands managed by the Bureau of Land Management (BLM) and the Department of Defense (DOD). A previous U.S. Geological Survey report documented variability in satellite imagery-based land cover classifications for seven wetland complexes at near monthly time intervals between October 2015 and January 2022. This report presents additional data, extending analysis to November 2023. Land cover classifications between October 2015 and November 2023 demonstrated an association between vegetation cover characteristics and surface moisture, with Class 1 having dry, bare soil or sparse upland vegetation, Class 2 having moist, bare soil or sparse to small vegetation, Class 3 having dense green vegetation with potentially saturated soil conditions, Class 4 having a mix of shallow surface water, saturated soil, and dense green vegetation, and Class 5 having open surface water. Most of the wetland complexes occur close to spring outflows primarily within land managed by the DOD, though portions are also within BLM lands. The intervening and surrounding landscape outside of the wetland complexes assessed in this study are managed by the BLM. As a result, Class 1 land covers had the largest areal coverage for BLM managed lands. Classes 2 and 3 land covers were primarily mapped inside the wetland complexes and thus had the largest area coverage within DOD managed lands. Class 4 was almost exclusively mapped within the wetland complexes and thus was largely contained within DOD managed lands. Class 5 (open water) was exclusively mapped in and adjacent to a single wetland complex with catchment ponds on land managed by the BLM. The distribution of these land cover classes over the study period was seasonally and annually variable. Land cover areas of Classes 1 and 2 were larger during the spring months. Conversely, land cover areas of Classes 3 and 4 tended to be greatest during the summer or fall. These patterns might be influenced by differences in seasonal water sources and phenology.

Nevada

Comparison of hydrologic data and water budgets between 2003–08 and 2018–23 for the eastern part of the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer is divided spatially into three parts (eastern, central, and western). The largest groundwater withdrawals are from the eastern part of the Arbuckle-Simpson aquifer, which provides water to approximately 39,000 people in Ada and Sulphur, Oklahoma, and surrounding areas. The Arbuckle-Simpson aquifer, including the eastern part, is designated a sole source aquifer for its service area. Based primarily on data collected between 2003 and 2008, a series of comprehensive hydrologic studies of the Arbuckle-Simpson aquifer was published to provide the information necessary to perform groundwater-flow model simulations so that the Oklahoma Water Resources Board could determine how much water could be withdrawn from the aquifer while maintaining flow to springs and streams. As part of the Phase 1 studies, an aquifer water budget was developed from a numerical model for the period 2003–08. For this report, Phase 1 refers to the 2003–08 data collection period, although for some of the analyses, data collected prior to 2003 were used to inform model development work. Allocation of water from this aquifer was then established by the Oklahoma Water Resources Board in 2013. Additional well-spacing rules were also established by the Oklahoma Water Resources Board for sensitive sole source groundwater basins. To determine how the water budget for the eastern part of the Arbuckle-Simpson aquifer has changed over time, recently collected hydrologic data (2018–23) were compared to data collected during 2003–08. The analysis of changes in the aquifer water budget from 2003–08 to 2018–23 could help resource managers better understand changes in the overall balance of water in storage and the potential effects on streamflow, changes in groundwater levels, and the effects of different water uses in the aquifer area on available water in the eastern part of the Arbuckle-Simpson aquifer and streams overlying the eastern part of the Arbuckle-Simpson aquifer.

Oklahoma

Oligocene–Miocene northward growth of the Tibetan Plateau: Insights from intermontane basins in the West Qinling Belt, NW China

Growth of the Tibetan Plateau, Earth’s broadest and highest elevation collisional system, shapes orographic barriers, reorganizes drainage networks, and influences surface erosion and sediment delivery, whose changes in space and provenance feed back to intracontinental tectonic processes. Studies of interior basins within the northern Tibetan Plateau provide new sediment accumulation, provenance, paleodrainage, and deformation timing data that enable a reconstruction of the far-field tectono-geomorphic evolution of the rising Tibetan Plateau. Along the northern plateau margin, topographic growth in the West Qinling Belt is inferred to have initiated in the Eocene, nearly coeval with the India-Asia collision, as well as in the late Miocene. However, geological knowledge about the intervening period remains at present enigmatic, and the kinematics and dynamics are uncertain. This study presents a multidisciplinary data set from the intermontane Anhua-Huicheng Basin (AHB; Gansu Province, China) to fill this gap. Magnetostratigraphic dating, regional mapping, and sedimentological analysis imply that contractional deformation and thrust-top basin systems formed within the West Qinling Belt in the Oligocene (not later than ca. 24 Ma). A combination of observations including paleocurrent changes, detrital zircon U-Pb age variations, and appearance of growth strata along the Anhua-Huicheng Basin reveal the rapid uplift of the West Qinling Belt at ca. 15 Ma. Sedimentation in the intermontane basins ended after the late Miocene (ca. 8 Ma), when the region experienced intrabasinal deformation, uplift, and erosion with the establishment of an external drainage system. Since the late Miocene, the growth of the West Qinling Belt reached a climax with the lack of substantial contractional deformation in Cenozoic sequences heralding the onset of the modern kinematic regime and attainment of high elevation. Observed transitions in the tectonostratigraphy and paleodrainage define different phases of deformation and plateau-wide shifts in stress reorganization, which led to the northward growth and later lateral expansion of the Tibetan Plateau.

GSA Bulletin

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay

Monitoring marine eruptions

Introduction Submarine volcanoes produce much of the same seismicity and eruptive activity as subaerial volcanoes and can pose hazards to society. Although they can be monitored with similar techniques and methods as described in other chapters of this volume, their submerged location brings unique challenges. This chapter addresses these challenges and provides recommendations for monitoring volcanoes fully or partly in marine environments to meet the capabilities described in other chapters of this volume. The United States and its territories host dozens of submarine volcanoes with most (around 60) in the Commonwealth of the Northern Mariana Islands. Approximately 20 of the Northern Mariana Islands submarine volcanoes are known to be hydrothermally active, and 10 have confirmed eruptions since the 1950s (for example, Baker and others, 2008; Tepp and others, 2019a). Nine of those volcanoes were considered by the National Volcanic Threat Assessment (Ewert and others, 2018) to have a combination of eruptive type and summit depth that poses a higher risk of hazardous eruptions, although only one was listed as a moderate (level 3) threat. Other notable submarine volcanoes of interest to the United States that have historically erupted are Axial Seamount off the Washington State coast, Kamaʻehuakanaloa in Hawaiʻi, and Vailuluʻu seamount in American Samoa. All of these, however, have a low risk of hazards because of their depth (greater than 600 meters below sea level) and eruptive type and so are not included in the National Volcanic Threat Assessment. In addition to submarine volcanoes, the submerged flanks of island volcanoes can also be a source of hazardous submarine eruptions—for example, the 1877 eruption of Mauna Loa, Hawai‘i, in Kealakekua Bay (Wanless and others, 2006). The most notable submarine eruption in recent times was the 2022 eruption of Hunga Tonga–Hunga Haʻapai in Tonga, which was one of the largest eruptions on Earth in the past 100 years. It created a massive volcanic plume, unprecedented shock waves, and far-reaching tsunami (Lynett and others, 2022). Other recent submarine eruptions in the Pacific Ocean Basin have produced subaerial plumes that reached aircraft heights (Carey and others, 2014) and large pumice rafts that can affect marine traffic and harbors (for example, Jutzeler and others, 2014; Kornei, 2019). These examples illustrate the potential hazards of major submarine eruptions. Yet, submarine volcanoes are largely unmonitored, and many eruptions occur that are unnoticed or only identified hours or days afterward. Within U.S. territory, submarine volcanoes in the Northern Mariana Islands have been known to produce eruptive activity that can affect society. Reports from fishermen and other marine vessels in the Northern Mariana Islands have noted underwater explosions, sea-surface discoloration, and bubbling water, all of which are known to be signs of submarine volcanic activity. South Sarigan seamount, located about 160 kilometers (km) north of Saipan, erupted in 2010 from greater than 150 meters below the sea surface, resulting in a gas and ash plume that reached more than 11.9 km into the atmosphere (for example, Searcy, 2013; Embley and others, 2014), high enough to affect international air traffic. Precursory and co-eruptive seismicity was detected on the regional Northern Mariana Islands seismic network (Searcy, 2013) and on global monitoring instruments (Green and others, 2013). Monitoring of submarine volcanoes is best accomplished with marine-based instrumentation, which is also useful for monitoring small island volcanoes that may not have the land area necessary for comprehensive subaerial monitoring. The primary marine-based instrumentation used for submarine volcanoes includes ocean-bottom pressure sensors to assess sea-floor deformation, ocean-bottom seismometers (OBSs) to detect seismicity, and both moored and ocean-bottom hydrophones to detect submarine explosions. Other sensors offer important monitoring data, such as turbidity, temperature, and chemistry of hydrothermal emissions. Marine-based instruments are typically deployed in campaign-style networks with no real-time telemetry owing to cost considerations and technical limitations. However, when necessary, marine instruments can be operated in real time using cables to transmit data to land-based facilities; other technologies for this purpose are in use or in development, such as acoustic transmission from the instrument to a moored buoy (Matsumoto and others, 2016) and a winch-based system with a satellite antenna that is part of the instrument mooring (Matsumoto and others, 2019). Emerging technologies for marine-based monitoring may be considered as part of a future monitoring plan. These technologies include ocean gliders and floats with on-board hydrophones that have been used to record earthquakes and submarine eruptions (for example, Matsumoto and others, 2013; Sukhovich and others, 2015) and fiber-optic cables that have been used as strainmeters to detect earthquakes (for example, Marra and others, 2018; Lindsey and others, 2019). Land-based instruments and satellites can also provide some capability for monitoring submarine volcanoes, but they provide more limited observations than marine-based instrumentation.

Scientific Investigations Report

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

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.

EarthArXiv

Technique for estimating magnitude and frequency of peak flows in Maryland

A convenient and reliable technique for estimating flood magnitudes is required for effective flood-plain management and for the efficient design of bridges, culverts, embankments, and flood-protection structures. Methods are presented for estimating peak-flow magnitudes of selected frequencies, ranging from 2 to 500 years, for all nontidal drainage basins in Maryland. The methods were developed by generalized least-squares regression techniques using data from 219 gaged basins in and near Maryland. The State is divided into five hydrologic regions: the Appalachian Plateaus and Allegheny Ridges region, the Blue Ridge and Great Valley region, the Piedmont region, the Western Coastal Plain region, and the Eastern Coastal Plain region. These regions correspond to the physiographic provinces of the State, with the exceptions that (1) the Coastal Plain Province is divided into two hydrologic regions, and (2) there is no distinct hydrologic region corresponding to the Valley and Ridge Province as it is divided into its constituent Allegheny Ridges and Great Valley subdivisions. Sets of equations for calculating peak discharges based on physical basin characteristics are provided for each of the regions. Based on the peak-flow equations, methods for estimating peak flows are presented for ungaged and gaged streams in Maryland. The methods and equations are supported by generalized least-squares analysis of basin and flood-frequency characteristics data from 219 drainage basins in and near Maryland. estimates for each of the five regions are calculated using combinations of the fol-lowing basin characteristics: drainage area, forest cover, basin relief, carbonate rock coverage, storage, and runoff-curve number. Drainage area contributes to the estimate in all five study regions. Carbonate rock coverage is used only in the Blue Ridge and Great Valley region. Storage and runoff-curve number are used solely in the Eastern Coastal Plain region. All other basin characteristics are used in two or more regions. Standard errors of estimate for the regression equations range from 19 to 31 percent in the Appalachian Plateaus and Allegheny Ridges region, 34 to 47 percent in the Blue Ridge and Great Valley region, 33 to 48 percent in the Piedmont region, 45 to 64 percent in the Western Coastal Plain region, and 36 to 42 percent in the Eastern Coastal Plain region.

Maryland

Linking distribution and return-on-investment models to optimize woody management for prairie grouse in Nebraska

Grasslands in Nebraska, USA, face threats from agricultural conversion, urban development, and woody encroachment, all of which negatively affect prairie grouse ( Tympanuchus spp.) populations. To optimize conservation planning, Nebraska wildlife agencies developed probabilistic area-based surveys for greater prairie-chicken ( T. cupido ) and sharp-tailed grouse ( T. phasianellus ) to sample landscapes across a range of environmental conditions. This design improves historical surveys and enables the development of distribution models that quantitatively define habitat associations and support scenario-based conservation planning. Using survey data collected during 2020–2022, we modeled prairie grouse occurrence and abundance as functions of land cover, topography, and climate using Bayesian logistic and zero-inflated negative binomial models with regularized horseshoe priors. We then conducted a maximum potential return-on-investment analysis of woody cover treatments, assuming sustained treatment success, relative to projected impacts of woody encroachment on prairie grouse populations by 2050. Among modeled associations were a positive association with grasslands having low woody cover and a negative association with grasslands having high woody cover. Across the 3-year period, median estimated annual populations were 142,380 for greater prairie-chicken (range of 95% CIs across years = 68,821–277,615) and 64,154 for sharp-tailed grouse (range of 95% CIs across years = 27,550–146,559). Under projected woody encroachment, mean predicted population declines were 10% for greater prairie-chicken (range of 95% CIs = 7–14%) and 6% for sharp-tailed grouse (range of 95% CIs = 5–7%). Areas with high prairie grouse density and low treatment costs in 2021, and high projected woody encroachment and population loss by 2050, offered the greatest return on investment for woody management. Return on investment was greatest in the northwestern Shortgrass Prairie ecoregion (northwestern Nebraska) for sharp-tailed grouse and the eastern Sandhills ecoregion (central Nebraska) for both species. These models underscore the value of evidence-based, quantitative approaches for prioritizing conservation actions on working lands. Scenario-based modeling could be extended to guide other treatments, such as optimizing restoration (e.g., Conservation Reserve Program) or incentivizing grassland persistence in areas with predicted climate resilience.

Nebraska

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–23, California Groundwater Ambient Monitoring and Assessment Program, Priority Basin Project

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–23. 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 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 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 total dissolved solids 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 methyl tert -butyl ether are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

California

Trace elements in streambed sediment and fish liver at selected sites in the Upper Colorado River Basin, Colorado 1995-96

Trace elements were analyzed in streambed-sediment samples collected from 16 sites and in fish-liver samples collected from 14 sites in the Upper Colorado River Basin in Colorado as part of the National Water-Quality Assessment program. Sites sampled represented agricultural, mining, mixed, and urban/recreation land uses and background conditions. The results for 15 trace elements in streambed-sediment and in fish-liver samples are presented in this report. Fourteen of the selected trace elements were detected in streambed-sediment samples collected at all sites. Twelve of the selected trace elements were detected in fish liver at more than 50 percent of the sites. Cadmium, copper, selenium, and zinc were selected for a more detailed analysis. Cadmium, copper, and zinc concentrations in streambed sediment were highest at mining land-use sites in the Southern Rocky Mountains physiographic province. Selenium concentrations in streambed sediment were highest at an agricultural land-use site in the Colorado Plateau physiographic province. The concentration of trace elements in streambed sediment generally increased as particle size decreased. Concentrations of trace elements in fish liver generally did not follow the same relation to land use as concentrations in streambed sediment; however, cadmium concentrations in fish liver were highest at a mining land-use site in the Southern Rocky Mountains physiographic province, and selenium concentrations in fish liver were highest at an agricultural land-use site in the Colorado Plateau physiographic province. Copper and zinc concentrations in fish liver were highest at mixed land-use sites. Comparison of streambed-sediment and fish-liver concentrations to two other similar NAWQA studies in the Rocky Mountain region generally indicated similar patterns in relation to land use for streambed sediment, but not for fish liver. Cadmium, copper, and zinc concentrations in streambed sediment were highest at sites affected by mining in all three study units. Selenium concentrations in streambed sediment did not indicate relations among the three study units when compared to land use. Cadmium in fish liver was highest at sites affected by mining in all three study units. Copper, selenium, and zinc in fish liver did not indicate relations among the three study units when compared to land use.

Colorado

Metal fingerprints of Eocene rhyolite magmas coincident with Carlin-type gold deposition in Nevada USA

Eocene magmatic systems contemporaneous with world-class Carlin-type Au deposits in Nevada (USA) have been proposed by some researchers as a key ingredient for Au mineralization, though evidence conclusively demonstrating their genetic relationship remains tenuous. This study provides the first direct evidence of the pre-eruptive metal budget of volatile- and metal-charged silicic magmas coincident in time (~41 to 34 Ma) and space (within 5 km) with Carlin-type Au deposits. We characterize the pre-eruptive metal fingerprints of these diverse magmatic systems to assess their potential as sources of metals for Carlin-type Au mineralization. Metal abundances from quartz-hosted melt inclusions (Au, Te, Ag, Sb, Tl, Mo, W, Sn, As, Pb, Co, Cu, Ni, and Zn) characterized in situ by SHRIMP-RG and LA-ICP-MS represent our best (and only) estimates for the pre-eruptive metal budget in these systems. Median metal concentrations are generally within one order of magnitude of average upper crust and average continental rhyolite values. But there are two notable exceptions, with median Au contents extending >1 order of magnitude higher than average upper crust and median Cu contents ranging >1 order of magnitude lower than upper crust. Despite this, melts contain lower Au/Cu (<0.1), Au/Ag (<5), and Au/Tl (<0.3) than most ore-grade Carlin-type rock samples and quartz-hosted fluid inclusions, regardless of their age and timing relative to nearby Carlin-type Au mineralization. The metal fingerprints of these magmatic systems, de-fined both by traditional and multivariate compositional data analysis techniques, are distinct from one another. Yet none are particularly specialized, e.g., high Au/Cu, in terms of being ideal ingredients as postulated by magmatic models for Carlin-type Au mineralization. Magmatic Au contents do not appear to be correlated with rhyolite “flavors” in the way that Cu, Sn, and Nb contents are. Fluid/melt partitioning modeling and magma volume estimates support the idea that a diverse array of non-specialized silicic magmas could feasibly contribute some or potentially all of the Au, Ag, and Cu in Carlin-type systems. The compositional diversity among contemporaneous magmatic systems could possibly contribute to some of the diversity observed across Carlin-type Au districts in Nevada.

Nevada

Potential for leakage among principal aquifers in the Memphis area, Tennessee

The principal aquifers in the Memphis area consist primarily of sand or sand and gravel, and the confining beds consist of clay, silt, sand, and lignite. The Jackson Formation and upper part of the Claiborne Group serve as the confining bed separating the water table aquifers from the Memphis Sand; the Flour Island Formation separates the Memphis Sand from the Fort Pillow Sand. Differences in total hydraulic head among the principal aquifers in the Memphis urban area result in vertical hydraulic gradients which create a potential for inter-aquifer exchange of water. Throughout this area, the gradient is downward from the water table aquifers to the Memphis Sand. In the central part of the Memphis urban area, the vertical hydraulic gradient is upward from the Fort Pillow Sand to the Memphis Sand, and in the eastern and western parts, it is downward from the Memphis Sand to the Fort Pillow Sand. The vertical distribution of carbon-14 data for water from the fluvial deposits, Memphis Sand, and Fort Pillow Sand shows an increase in the relative age of the water with depth. The areal distribution of carbon-14 data for water from the upper part of the Memphis Sand indicates that relatively recent water has been brought into the major cone of depression in the potentiometric surface of the Memphis Sand, either by horizontal movement or from downward vertical leakage. The normal, near-surface geothermal gradient in the Memphis area was determined to be 0.6 C/100 ft. Deviations from the normal geothermal gradient, in areas affected by intense pumping from the Memphis Sand, indicate that downward vertical leakage occurs from the water table aquifers through the Jackson-upper Claiborne confining bed to the Memphis Sand. The velocity of downward vertical leakage of water from the Memphis Sand through the Flour Island confining bed to the Fort Pillow Sand was determined to be 0.0066 ft/day by analysis of borehole temperature data from an observation well in the northeastern part of the Memphis area. From this velocity and the head difference between the Memphis Sand and the Fort Pillow Sand at this locality, the hydraulic conductivity of the Flour Island confining bed was determined to be 0.00114 ft/day. (Lantz-PTT)

Tennessee

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho

Water use across the conterminous United States, water years 2010–20

Withdrawals of water for human use are fundamental to the evaluation of the Nation’s water availability. This chapter provides an analysis of public supply, crop irrigation, and thermoelectric power water use for the conterminous United States (CONUS) during water years 2010–20. These three categories account for about 90 percent of water withdrawals in the Nation. The values presented here are based on modeling approaches that estimate water use at temporal (monthly) and spatial scales (12-digit hydrologic unit code—small watersheds sized 50–100 square kilometers) compatible for integration into a broader national assessment of water availability. Models also provide an understanding of factors that influence water use. An estimated 244,817 million gallons per day (Mgal/d; 28,677 million cubic meters per month [Mm 3 /mo]) were withdrawn on average within the CONUS during water years 2010–20 from fresh water and saline water for crop irrigation, public supply, and thermoelectric power, with shares of 43, 14.5, and 42.5 percent for each of these categories, respectively. In the same period, estimated withdrawals and consumptive use (1) for public supply were 35,400 and 4,219 Mgal/d (4,081 and 486 Mm 3 /mo), respectively; (2) for crop irrigation were 105,497 and 75,698 Mgal/d (12,147 and 8,716 Mm 3 /mo), respectively; and (3) for thermoelectric power from fresh water were 82,656 and 2,904 Mgal/d (9,952 and 345 Mm 3 /mo), respectively. Withdrawals for these categories of water use are highly spatially variable, with western States dominated by crop irrigation and eastern States dominated by thermoelectric-power water use. Public supply accounts for the largest percentage of water use in several heavily populated northeastern States. Reliance on groundwater compared to surface water depends on the availability of water sources and the type of water use. For public supply, withdrawals from groundwater are greater than withdrawals from surface water in the Western aggregated hydrologic regions, whereas the balance shifts to more surface water for the rest of the CONUS. In all aggregated hydrologic regions, the predominant source of water for crop irrigation is groundwater. Most thermoelectric power facilities in the eastern half of the CONUS use surface water from freshwater and saline sources; most thermoelectric power facilities in the western half of the CONUS use groundwater.

conterminous United States