Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Water Research”

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.

638 records · Page 6Linked to original sources
USGS2026

Synthesizing a twelve-year sediment trap time series of planktic foraminiferal flux in the Gulf of America (Mexico)

Sediment trap time series provide powerful frameworks for testing hypotheses about planktic foraminiferal assemblage composition, seasonality, and geochemical responses to environmental variability, all of which are central to improving paleoceanographic reconstructions. We present high resolution foraminiferal assemblage data from a long-running sediment trap (2008–2020) in the northern Gulf of America (Mexico). This study summarizes the species composition, seasonality, and size distribution of the fifteen most abundant species of planktic foraminifera, which account for 98% of total flux. Foraminiferal flux peaks in winter and reaches a minimum in summer, following the seasonal pattern of primary production in the northern Gulf. Winter assemblages are dominated by non-spinose taxa, whereas spinose taxa prevail during summer. Across nearly all species, average monthly test size covaries with temperature, independent of seasonal flux trends. Notably, Trilobatus sacculifer and Neogloboquadrina dutertrei show a significant decline in relative abundance and flux during 2017–2020, nearly disappearing from the assemblage.

Journal of Foraminiferal Research↗
USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗
USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗
USGS2025

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗
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

A case study of maximum depositional ages from terrestrial sandstones near the Cretaceous—Paleogene transition, western Williston Basin, USA

We present a new Bayesian method for deriving maximum depositional ages from detrital K-feldspar using total-fusion 40 Ar/ 39 Ar data. Individual analyses weighted by K/Ca ratio, age uncertainty, and percent radiogenic argon prioritize a result to come from accurately measured volcanic sources. Results from this method applied to sandstones from the Hell Creek region of northeastern Montana show that detrital K-feldspar maximum depositional ages align closely with detrital U-Pb zircon ages as well as tephra based chronostratigraphic constraints from both 40 Ar/ 39 Ar and U-Pb data. An age model informed by maximum depositional ages and available tephra data provides an estimate of 1.74 ± 1.04 Ma for the duration of deposition for the Hell Creek Formation. Combined age model and maximum depositional age data suggest ≤ 1.84 Ma of missing time is represented by the Hell Creek and Fox Hills formational disconformity, coinciding with the final regression of the Western Interior seaway.

Montana↗
USGS2024

(Re)discovering the seismicity of Antarctica: A new seismic catalog for the southernmost continent

We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on‐continent temporary and long‐term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (⁠Mw ⁠) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice‐covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Seismological Research Letters↗
USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗
USGS2025

Climate-driven sulfate export in alpine watersheds may stimulate methylmercury production

Climate change is increasing sulfate export and changing wetland extent in mountain regions. These changes may increase microbially mediated production of the neurotoxic substance methylmercury due to enhanced sulfate metabolism in mountain environments. Here, we assess methylmercury concentrations and formation rates across high-elevation wetlands in the Colorado Rocky Mountains. We also investigate sulfate controls on methylmercury production within subalpine peatlands by amending soils with sulfate to mimic increased stream export of sulfate from the alpine zone and measuring methylmercury formation rates for different sulfate treatments. We found that subalpine peatlands have statistically significant higher methylmercury concentrations and formation rates compared to alpine, mineral-soil wetlands. Methylmercury production in subalpine peatlands also increased significantly ( p < 0.05) following sulfate additions; the highest rates occurred in sediments with intermediate extractable sulfate concentrations (∼0.60–1.4 mg sulfate g −1 dry soil). Our study is the first to identify soil sulfate-related thresholds for methylmercury production and sulfate-limitation of methylmercury production in subalpine peatlands. These findings highlight important linkages between climate-driven mineral weathering and mercury cycling in mountain regions globally.

Colorado↗
USGS2025

Geochemical evidence for the origin of late Quaternary loess, Seward Peninsula, Alaska

Loess is the most widespread surficial deposit in the state of Alaska. Although loess of last glacial age is common in mid-continental North America, records of last glacial loess in Alaska have been elusive. Here we report a record of last glacial loess on the Seward Peninsula, along with studies of a transect of loess localities also thought to date to the last glacial period. Highly diverse hypotheses on the origin of loess on the Seward Peninsula include local volcanic rocks and mountain ranges to the north and south, as well as offshore deposits. We present particle size and mineralogical and geochemical data that indicate that the most likely sources for loess in this part of eastern Beringia are glaciogenic silt particles of the Kobuk River and Noatak River, both of which drain the Brooks Range to the northeast. Northeasterly paleowinds were therefore responsible for loess on the Seward Peninsula. Such winds probably occurred during late summer or early autumn, when river discharge was at a minimum but temperatures were still above freezing. Northeasterly winds during the last glacial period are in good agreement with a recent paleoclimate model.

Alaska↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗
USGS2024

Joint pilot fish habitat framework

This story map will take you through the process of exploring and testing methods necessary for a higher resolution, seamless fish habitat assessment across both inland and estuarine waters through the lens of our joint pilot assessment Fish habitat assessments attempt to relate past, current, or future landscape conditions to the state of fish species occurrence, distribution, abundance, or community and habitat condition in streams, rivers, or estuaries. Previous fish habitat assessments, such as the National Fish Habitat Assessment, conducted separate and disconnected assessments for inland waters and estuaries. In this project, National Oceanic and Atmospheric Administration ( NOAA ) and U.S. Geological Survey ( USGS ) researchers created a seamless spatial framework to allow assessments that integrate influences on fish habitat from headwaters to the estuary. This effort began when the Chesapeake Bay Program Fish Habitat Action Team expressed interest in a Baywide fish habitat assessment spanning tidal salt, tidal fresh, warm non-tidal, and cold non-tidal waters. However, the complexity of the myriad of implementation details to consider when developing such an assessment necessitated the need for a tributary-specific pilot assessment. To conduct this pilot assessment, a NOAA/USGS joint partnership was formed with cooperation and support from the Chesapeake Bay Agreement and Chesapeake Bay Fish Habitat Action Team (FHAT).

Maryland↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS2026

Digitizer Suite: The Albuquerque Seismological Laboratory Digitizer Testing Suite

Laboratory testing of digitizers and seismometers helps ensure that prior to deployment the instrumentation can produce high quality data and is operating within specifications. In this work we detail the software package called: the Albuquerque Seismological Laboratory (ASL) Digitizer Test Suite. This Java software package provides several algorithms to verify various performance parameters of digitizers commonly used for recording analog seismic instruments. The goal of these tests is not to be exhaustive, but to identify common failures that could compromise the integrity of seismic data being recorded on the digitizer. For example, Sandia National Laboratories (e.g., Slad and Merchant, 2018) routinely do comprehensive testing of digitizers for various monitoring missions. While these tests reports are valuable for comprehensively characterizing a recording system, it would be resource intensive to conduct such tests on every seismic recorder used in a network. We focus on tests that include ways to estimate the sensitivity, timing, self-noise, and clip-level of the digitizer, as well as the fidelity of the signal being recorded. The software is publicly available and provides a way for the community to verify the integrity of a digitizer using a minimum amount of outside equipment.

Seismological Research Letters↗
USGS2024

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 2 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km 2 ) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 2, the subject of this scientific investigations map, has a mapped area of 209 km 2 and has water depths that range from about 19 meters (m) on the Stellwagen Bank crest to about 68 m in the Stellwagen Basin. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates across the southwestern part of Stellwagen Bank, in Stellwagen Basin to the west and southwest of the bank, and in Little Stellwagen Basin and the western part of Race Point Channel to the south of the bank. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 733 stations were analyzed, including 656 sediment samples. The geologic substrate maps of quadrangle 2 show the distribution of 19 geologic substrates that represent a wide range of textures, such as rippled and immobile sand, immobile muddy sand and sandy mud, sand that partially veneers gravel, and a boulder ridge. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 2. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that utilize the various substrates, and support seabed management in the region.

Massachusetts↗
USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗
USGS2025

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗
USGS2026

Restoring dryland water cycles for precipitation feedback and climate stability; a review

Drylands across the globe are experiencing intensifying water scarcity, land degradation, and hydroclimatic extremes. This review integrates evidence from multidecadal field studies, hydrologic monitoring, geomorphic and ecological assessments, remote sensing, and land–atmosphere science to evaluate how restoration influences key components of the terrestrial water cycle. Low-tech natural infrastructure in dryland streams (NIDS)—including check dams, leaky weirs, one-rock dams, and gabions—has emerged as a promising but under-synthesized nature-based solution for restoring hydrologic function in these environments. We describe the mechanisms through which these interventions modify runoff detention, infiltration, sediment and alluvial storage, shallow-groundwater recharge, vegetation recovery, and surface-energy partitioning, and we summarize outcomes across diverse dryland settings. Findings consistently show increased water residence time, enhanced soil-moisture storage, expanded riparian vegetation, extended flow duration, and shifts toward greater latent-heat flux—producing localized cooling and strengthened ecohydrological feedbacks. Building on these localized effects, we articulate a hypothesis that links the spatial extent of restoration, the density of NIDS per unit drainage area, and the magnitude of the latent-to-sensible-heat contrast generated by wetter post-rainfall conditions. Specifically, we hypothesize that when NIDS are implemented at densities permitted by topography and across areas large enough to maintain elevated soil moisture after storm events, the resulting increases in latent heat flux, surface cooling, and boundary-layer moistening may enhance moisture convergence and boundary-layer development, potentially increasing the likelihood or stability of convective precipitation, analogous to how reductions in these processes have contributed to regional drought intensification. These land–atmosphere feedbacks remain untested at scale but represent an important research Frontier. By integrating hydrologic, geomorphic, ecological, and atmospheric perspectives, this review provides a comprehensive framework for considering how low-tech, landscape-scale interventions can strengthen watershed resilience and contribute to climate-relevant nature-based solutions.

Arizona↗