USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

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.

1,666 records · Page 89Linked to original sources
USGS2026

Relating glassy rind thicknesses to ambient air temperatures at the Lost Jim flow field in the Imuruk Lake volcanic field, Alaska

The Lost Jim flow field, in the Imuruk Lake volcanic field, Alaska, extends west ~34 km from a single vent, crossing subarctic tundra and currently touches several lakes and streams. The weighted mean of five 36Cl cosmogenic exposure ages from the Lost Jim pāhoehoe flow is 7.73 ± 0.37 ka, indicating this eruption occurred substantially after the eruption of the underlying Camille flow, which was emplaced at 39.7 ± 1.3 ka. Paleoclimate records indicate the period when the Lost Jim flow field was emplaced was after deglaciation, and the climate was similar to today. We propose that the emplacement of lava in these cold subarctic conditions can lead to faster cooling of the lava surface compared to lava emplaced in warmer locations such as mid- latitude cold deserts. Glass abundance in the outermost rinds at the Lost Jim flow field was on average 74 % with 6.4 mm thick rims, compared to 60 % with 2.9 mm rims for cold mid-latitude desert samples. We interpret increased glass content as a proxy for rapid cooling likely occurring partly during winter. Glassiness values varied less across vent, margin, and mid-flow locations when compared to the mid-latitude flows suggesting the Lost Jim flow field was broadly impacted by the subarctic climate as opposed to responding to local microclimates. Our results indicate that lava glassiness may be a useful environmental indicator of cooler (in this case subarctic) conditions.

Alaska↗
USGS2023

WellSTIC: A cost-effective sensor for performing point dilution tests to measure groundwater velocity in shallow aquifers

Many individual measurement points are required to characterize groundwater velocity within an aquifer. Groundwater velocity is most commonly measured using a network of >5 cm diameter monitoring wells, which, if not already present at a site, are expensive and labor-intensive to install. Drive-point piezometers—simple, cost-effective wells that can be installed by hand—are a common tool for sampling groundwater in shallow, alluvial aquifers, but most groundwater velocity measurement techniques require equipment that is too large for these narrow (usually <2 cm inside diameter) piezometers. In this technical note, we introduce a low-cost sensor and well packer system (<$90 USD) for performing point dilution tests in narrow piezometers. Field data show that the magnitude of groundwater velocity measured with this technique agrees with velocities computed from natural gradient tracer tests. Additionally, with proper calibration, these sensors can be used to continuously monitor in-well specific conductance, either during inter-well tracer tests with saline tracers or for water quality monitoring. This system is a viable tool for rapid assessment of the magnitude of groundwater velocity in shallow aquifers.

Water Resources Research↗
USGS2026

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota↗
USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗
USGS2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi↗
USGS2026

Emergent headwater types for watershed monitoring, research, and planning

Headwaters—that is headwater streams and the watersheds they drain—strongly influence water quality, ecosystem services, and hydrologic connectivity, yet remain poorly mapped, monitored, and understood at large spatial scales. Using available spatial data to classify headwater “types” could therefore enhance the utility of sparse monitoring data to advance research and watershed planning—and help update literature-based perceptual models of headwaters as steep, forested, and minimally disturbed watersheds. Here, we analyze headwaters across the contiguous United States to quantify how their attributes (climate, topography, land cover, soils) vary regionally and differ from their downstream watersheds. Our analysis (a) identified seven dominant headwater types with distinct combinations of biophysical and climatic attributes and (b) revealed systematic contrasts with downstream watersheds, with strongest differences in attributes representing potential flowpath depth. Our findings offer a transferable framework for classifying, modeling, and managing these critical yet understudied components of the hydrologic network.

conterminous United States↗
USGS1996

Pesticides in surface water from three agricultural basins in south-central Georgia, 1993-97

Twenty-two of 43 pesticides analyzed were detected in 128 water samples collected from the Tucsawhatchee Creek, the Little River, and the Withlacoochee River. These streams drain agricultural basins in south-central Georgia and were sampled from March 1993 through June 1995. Herbicides were detected more frequently than insecticides. The most frequently detected herbicides were atrazine and metolachlor and the most frequently detected insecticide was carbaryl. Pesticide concentrations in the three streams were low and did not exceed U.S. Environmental Protection Agency drinking water standards. The maximum pesticide concentration was 2.6 mg/L (micrograms per liter) for propargite, a miticide detected in only one sample. The maximum concentrations of the remaining 21 pesticides were less than 0.25 mg/L. The median concentrations were equal to the method detection limit for all pesticides except atrazine (0.008 mg/L) and metolachlor (0.012 mg/L). The ratio of herbicide detections to nondetections was largest in the planting season, smaller in the harvest season and smallest in the fallow season for the three basins. The same pattern existed for the insecticide ratios in the Little River and the Withlacoochee River. Pairwise correlations between concentrations of atrazine and metolachlor, and four parameters (discharge, and concentrations of dissolved organic carbon, suspended organic carbon, and suspended sediment) were evaluated for each stream. The strongest correlations existed between metolachlor and mean daily discharge, and metolachlor and sediment in the Withlacoochee River. The only significant correlation for the Little River was between atrazine and suspended sediment.

Georgia↗
USGS2026

PFAS mixture composition and internal exposure profiles shape biological responses under field-realistic exposure

Per- and polyfluoroalkyl substances (PFAS) occur as complex mixtures, yet mixture-dependent biological effects under environmentally realistic exposure conditions remain poorly understood. We conducted multiyear (2018, 2019, 2021) continuous-flow, field-based exposures of male fathead minnows ( Pimephales promelas ) using a low-PFAS reference well (REF; sum of measured PFAS (∑PFAS) 0.1–0.2 μg L –1 ) and PFAS-contaminated groundwater from a fire-training area (FTA) at Joint Base Cape Cod, Massachusetts. Dilution treatments enabled separation of concentration and mixture effects. Groundwater from well FTA1 was perfluorooctanesulfonate (PFOS) dominated (∑PFAS 10–31 μg L –1 ), whereas groundwater from well FTA2 had higher concentrations and was enriched in perfluorooctanoate (PFOA) and diverse precursors (∑PFAS ∼ 80 μg L –1 ). Despite comparable plasma ∑PFAS and PFOS in FTA1–100% and FTA2–50% on day-7, cumulative mortality reached approximately 35% in FTA2–50% and 17% in FTA2–100%, and secondary sex trait expression was reduced by approximately 65–80% relative to REF, whereas sperm motility effects were mixture- and time-dependent. Plasma PFAS profiles were dominated by perfluorohexanesulfonate (PFHxS), PFOS, and sulfonamide precursors. Liver transcriptomics from REF-acclimated fish revealed robust disruption of metabolic, mitochondrial, and endocrine pathways that link PFAS mixture chemistry and internal exposure profiles to organismal outcomes, and testis transcriptomics provided complementary insight into reproductive impairment. These results indicate that PFAS mixture composition and internal exposure profiles, including precursor-associated differences, are important determinants of ecotoxicological outcomes under field-realistic conditions.

Massachusetts↗
USGS2025

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona↗
USGS2026

Measuring stress In high pressure deformation experiments with high speed fiber-optics

High‐pressure, high‐temperature rock deformation experiments are essential for understanding deep Earth processes, but accurately measuring stress is challenging due to the inaccuracy introduced by seal friction within the apparatus and large inertia, which affects earthquake process measurements made far from the experimental fault. To overcome these limitations, we developed and implemented a simple, compact fiber‐optic sensor based on an External Cavity Fabry–Perot Interferometer for in situ load measurement on a piston inside the pressure vessel of a Griggs‐type high‐pressure apparatus. The sensor can be used at sample temperatures up to 800°C for both slow creep and fast rupture testing with bandwidth capability from DC to 6 MHz. Two important experimental results are described: (1) the first direct measurements of seal friction in this apparatus type; these measurements record seal friction approximately one‐tenth of the confining pressure under both low (150 MPa) and high (1 GPa) pressure conditions; and (2) high‐temporal resolution capture of dynamic stress drops during stick‐slip faulting events at high‐confining pressure, revealing near‐fault wave propagation details, high‐frequency oscillations, and implying high‐transient slip rates (6–12 m/s) previously inaccessible with external sensors.

The Seismic Record↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗
USGS2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

Comptes Rendus. Géoscience↗
USGS2026

Extreme precipitation variability and soil texture controls on water-table response

Extreme precipitation events (EPEs), a key class of hydrometeorological extremes, are intensifying globally under climate change; however, their effects on water-table dynamics across varying soil textures remain poorly understood. To better understand the impacts of EPEs, we conducted one-dimensional modeling to evaluate water-table response time, displacement, recession time, and total recharge under EPEs of 0.20 m, 0.40 m, and 0.60 m amounts, applied over 1-, 7-, and 20-day durations across twelve soil textures. The results show that coarse soils (i.e., sand) respond within days, while fine soils (i.e., clay) may take over 200 days. Water-table displacement ranged from 0.30 to 1.64 m and increased with EPE magnitude. The time it took for water tables to recede ranged from 1.2 to 3.0 years. A first-order estimate of total possible recharge, calculated from porosity and displacement, ranged from 17% (clay) to 97% (sand), averaging ~63% across soil textures. These findings highlight that recharge is primarily governed by EPE magnitude and soil properties, not event duration. This modeling effort provides new insight into how soil texture modulates groundwater response to extreme precipitation, informing future water budget and resilience assessments.

Water↗
USGS1998

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report↗
USGS2024

Telecommunications fiber for sensing earthquake aftershocks: Progress and hurdles

Aftershocks offer valuable clues to earthquake behavior. The challenge: quickly deploying sensors to capture the early details of earthquake ruptures within the zone of aftershocks. Telecommunication fibers might be an answer, providing denser networks in otherwise difficult areas, potentially faster than traditional methods.

Conference Paper↗
USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗
USGS2025

Preventing overfitting when using tree-based methods for mapping hydrothermal favorability

Ensemble tree-based algorithms are robust tools for estimating sparsely distributed resources with non-linear dependencies (e.g., hydrothermal systems). These algorithms naturally accommodate the threshold conditions necessary to enable and support hydrothermal systems (e.g., having sufficient heat and permeability) and are simpler than many other non-linear machine learning strategies (e.g., artificial neural networks), which is an advantage when working with few labeled examples from which to learn. In previous work, we used eXtreme Gradient Boosting (XGBoost) to produce regional prediction and uncertainty maps of hydrothermal favorability; however, recent studies suggest that, even when properly applied, XGBoost has some risk of overfitting when there are few labeled examples from which to learn. To evaluate overfitting when constructing hydrothermal favorability maps with tree-based methods, we compare XGBoost with Extremely Randomized Trees (ExtraTrees), another ensemble tree-based algorithm that has the potential to underfit when using few labeled examples. We hold all other modeling parameters constant, resulting in two contrasting favorability maps of conventional geothermal resources for the Great Basin. Our results indicate that ExtraTrees demonstrably reduces overfitting compared with XGBoost. After considering overall performance, we conclude that ExtraTrees provides a more suitable modeling approach than XGBoost for the purposes of conventional hydrothermal resource assessments.

Conference Paper↗
USGS2026

A practical framework for identifying genetic subpopulations and ESUs: Insights for IUCN assessments and broader management

Species conservation assessments evaluate extinction risk, and recovery potential, advancing species persistence through guiding resource prioritization and planning. Assessment frameworks, including the International Union for Conservation of Nature Red List and Green Status of Species, typically focus on species as a whole. Importantly, they do not routinely account for genetically distinct units or do not have standardized methods of unit delineation. This limits the representation of genetically distinct components, including adaptive genetic diversity that underpins long-term resilience and recovery. Incorporating standardized within-species units like subpopulations and Evolutionarily Significant Units (ESUs) into species assessments could help address this oversight. However, identifying and delineating such units remain challenging, particularly when molecular data are limited. Here, we propose a flexible framework that integrates molecular and non-molecular evidence to identify both subpopulations and ESUs across taxa, providing a practical tool to incorporate within-species diversity into conservation assessments.

BioScience↗