SEARCH · USGS Science

Results for “Water-Resources Investigations Report”

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

SourceResearch connectionFrom yearThrough year

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

1,206 records · Page 32Linked to original sources

USGS2026

Standardized method for logging drill core at the Idaho National Laboratory, Idaho

The U.S. Geological Survey’s (USGS) Lithologic Core Storage Library (CSL) at the Idaho National Laboratory stores more than 120,000 feet of drill core that is accessible to the public for research and sampling. To effectively convey the physical and descriptive properties of the drill core, USGS staff at the Idaho National Laboratory Project Office log the drill core and publish the lithologic logs as data releases. The logs provide essential data on the lithology, texture, mineralogy, alteration, and other physical properties of the core, which serve as valuable information for researchers to guide their research and sampling efforts. To ensure consistent, quality, and dependable lithologic logs, this document outlines the procedures and expectations for logging drill core at the CSL. This document describes the processes for storing, photographing, and logging core, and includes a variety of resources, reference materials, and appendixes designed to standardize and aid the logging process. Following the procedures outlined in this document will produce consistent, detailed logs that facilitate dependable observations and serve as an easy reference for researchers and other interested parties.

Idaho↗

USGS2026

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

Mathematical Geosciences↗

USGS1976

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report↗

USGS2026

USGS Geochron—A database of geochronological and thermochronological dates and data—Technical documentation

Geochronological and thermochronological data are essential for constraining the timing and rates of geological processes, supporting geologic mapping, natural hazard assessment, and resource exploration. The U.S. Geological Survey (USGS) Geochron Database is a centralized, relational database that integrates USGS and State geological survey data in accordance with the National Geologic Mapping Act of 1992 and its reauthorizations. This report documents the structure and development of the database, including its standardized schema, controlled vocabularies, and protocols for compiling legacy and newly published data. The database is designed to adhere to FAIR (Findable, Accessible, Interoperable, and Reusable) data principles and currently (2026) includes data from a wide range of geochronological and thermochronological methods. Public access is provided through the USGS Geochron Database Explorer, a browser-based geographic information system (GIS) interface, and through versioned data releases available on ScienceBase, which provide sample-level summary data and detailed analytical information. The USGS Geochron Database is structured to be extensible, supporting the inclusion of additional methods and data types as they are compiled. This format ensures long-term utility and aligns with the Mapping Act’s directive to create a national archive of geochronological information that adds interpretive value to geologic map data. Its design reflects a commitment to data transparency, scientific reproducibility, and national geoscience priorities.

Data Report↗

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗

USGS1977

Acid-rain weathering of a metasedimentary rock basin, Herkimer County, New York

A geochemical balance of major chemical constituents of precipitation inflow with those of surface-water outflow was used to evaluate the effect of acid rainfall (average pH = 4.3) on weathering rates of metasedimentary rocks underlying a surface-water basin in New York State. Precipitation samples and surface-water samples show that approximately 45 percent of the net cation outflow is due to weathering by precipitation acidity. The remaining 55 percent of the cation load in the outflow is attributed to interactions between atmosphere, rock, and soil water in the unsaturated zone. Comparison of the acidity of acid precipitation and air-equilibrated distilled water suggests that since preindustrial time, precipitation acidity has probably increased weathering rates for this basin approximately 75 percent. The increased acidity of the precipitation may also increase the rate of dissolution of heavy metals from water-piping systems.

New York↗

USGS2024

Environmental monitoring of groundwater, surface water, and soil at the Ammonium Perchlorate Rocket Motor Destruction Facility at the Letterkenny Army Depot, Chambersburg, Pennsylvania, 2021

Letterkenny Army Depot in Chambersburg, Pennsylvania, built an Ammonium Perchlorate Rocket Motor Destruction (ARMD) Facility in 2016 to centralize rocket motor destruction and contain all waste during the destruction process. The U.S. Geological Survey has collected environmental samples from groundwater, surface water, and soils at ARMD since 2016. During 2021, samples were collected from four groundwater wells in September, one surface-water site in October, and five soil sites in November near the facility. Samples were analyzed for nutrients, trace metals, major ions, total volatile organic compounds, and perchlorate. Perchlorate was not detected in any 2021 samples. Groundwater results showed no constituents exceeded any U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL). Dissolved arsenic (As) was detected in one well above the reporting detection level (RDL) of 3 micrograms per liter (μg/L) at 5.4 μg/L but below its MCL of 10 μg/L. Dissolved iron (Fe) was the only inorganic constituent measured above an EPA secondary maximum contaminant level (SMCL). All groundwater samples collected in 2021 exceeded the Fe SMCL of 300 μg/L, with concentrations ranging from 390 μg/L to 3,500 μg/L. Surface-water data collected during 2021 showed no measured constituents in the surface-water sample that exceeded any EPA MCL or SMCL. Soil samples collected from 2016 through 2021 showed all concentrations of As exceeded the EPA soil screening levels of 3 milligrams per kilogram (mg/kg) but did not exceed the Pennsylvania medium-specific concentrations for As of 61 mg/kg. Arsenic concentrations in 2021 ranged from 9.1 mg/kg to 12.9 mg/kg. The 2021 results for the ARMD Facility indicate no increases in concentrations of reported compounds compared to data from 2016 to 2020. The contained burn treatment facility for demilitarization of rocket motors during 2021 appears to have operated without elevating concentrations of target compounds compared to previous years.

Pennsylvania↗

USGS2025

The U.S. Geological Survey National Water Quality Network—Groundwater, 2024

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, 8 of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,089 wells in the network. All wells are sampled for physical properties, nutrients, major ions, trace elements, per- and polyfluoroalkyl substances (PFAS), and a subset of wells are sampled for pesticides, volatile organic compounds, radionuclides, and microbiological contaminants.

Contiguous United States↗

USGS2025

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗

USGS2024

Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida

The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.

Florida↗

USGS2025

The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi‐year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage‐limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

Colorado↗

USGS2025

Scoping decision-maker needs and science availability to support regional natural capital accounting in the U.S. Colorado River Basin

Natural capital accounting has the potential to yield important policy insights at multiple scales, but there remains a disconnect between regional-scale natural capital accounts and their use for informing policy. In this paper, we propose a roadmap that could lead to the creation of policy-relevant regional accounts, with steps split across an initial scoping phase and a subsequent development phase. We demonstrate the scoping steps in action with an application to the Colorado River Basin (“Basin”), a large watershed in the southwestern United States (U.S.) that has faced aridification and substantial high-profile tradeoffs around the use of its water and other natural resources. Drawing on prior U.S. Geological Survey science co-production efforts, we conducted a series of eight discussion sessions with 41 scientists and science representatives whose work is relevant to Basin water, riparian and riverine ecosystems, upland ecosystems and energy and minerals. We summarise participants' thoughts on key topics and economic linkages, their insights and questions of interest and their recommendations on existing scientific data sources and gaps. We evaluate the suitability of the available data for construction of System of Environmental-Economic Accounting (SEEA) Central Framework and SEEA Ecosystem Accounting accounts, including those for land, water, forests, energy and minerals and ecosystems (covering extent, condition and ecosystem services). We present a series of lessons learned during the scoping phase, as well as lessons that could be relevant for future practitioners engaging in the development phase. The information can help guide the development of timely and relevant regional-scale environmental-economic accounts in the U.S. and beyond.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2026

Satellite time series analysis to quantify changing climax ciénegas using a state and transition model approach

Ciénegas are rare wetlands in arid landscapes of the North American Southwest, historically providing critical ecological and hydrological functions but increasingly threatened by changing climate and land use pressures. This study quantifies changes in ciénega condition and floodplain dynamics using a state-and-transition model (STM) informed by expert knowledge and remote sensing. Key factors include woody plant encroachment, water availability, and soil aggradation. We mapped 31 ciénegas with high-resolution imagery and analyzed Landsat data (1985–2023) to assess vegetation health and moisture using the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Infrared Index (NDII). Results show substantial interannual variability in phenology, water stress, and soil moisture, with regional drying and elevation strongly influencing ciénega resilience. We classified ciénegas into three functional states—healthy, desiccated, and dormant—and mapped their 2023 condition. Trend analyses indicate most ciénegas exhibit greening despite drought, though localized variability underscores the need for site-specific management. None are in a stable climax (reference) state; rather, they transition among states in response to external drivers. Increasing woody plant cover and surface drying, likely linked to declining regional water tables, favor deep-rooted species over wetland grasses—a pattern mirrored in adjacent control plots. Spatially explicit analysis revealed intra-ciénega variability often masked by aggregated data, highlighting the importance of high-resolution monitoring. Seasonal and long-term trends provide context for understanding ciénega dynamics, including degradation and restoration pathways. This study emphasizes the importance of groundwater conservation and demonstrates how remote sensing supports long-term monitoring. The STM framework offers a practical tool for adaptive management to sustain freshwater resources in arid environments.

Arizona, New Mexico↗

USGS2024

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper↗

USGS2026

Understanding the resource potential of natural hydrogen on Earth: Scientific gaps, uncertainties and recommendations

A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.

Earth-Science Reviews↗

USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗

USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗