USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Reports”

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.

751 records · Page 40Linked to original sources
USGS2025

Introduction to the special section on improving measurements of earthquake source parameters

Earthquake source parameters such as magnitude, seismic moment, source dimension, stress drop, and radiated energy are fundamental to understanding earthquake physics, and are also key ingredients in earthquake ground‐motion modeling, rupture simulation, and statistical seismology. However, the uncertainties in these parameters estimated from the radiated seismic wavefield are large due to variability in approaches, including site and attenuation characterizations, and so estimates for an individual earthquake made by different studies can vary greatly. Estimating spectral source parameters remains a popular topic, due to a combination of their intrinsic importance and their apparent ease of measurement, but also a controversial one, due to many sources of variability and large uncertainty. The available methods coupled with necessary parameter choices and assumptions in the analysis make it challenging for researchers to apply methods or understand the reliability in results or reported source parameters. This Special Section on Improving Measurements of Earthquake Source Parameters showcasing comparisons between methods and studies seeks to alleviate some of these difficulties to help the community identify the important components and trade‐offs of decomposing recorded seismograms into their source, path, and site components.

Bulletin of the Seismological Society of America↗
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

Triple-oxygen isotopic evidence of prolonged direct bioleaching of pyrite with O2

Sulfate is often touted as containing atmospheric oxygen whose isotopic signature can constrain redox, environmental conditions, and biological activity. Yet, the amount and isotopic fractionation associated with air-O 2 incorporation during sulfate formation is still debated, making its verification difficult. In this study, we identify a distinct, microbially dominated environment with the potential to preserve maximum signals of air-O 2 in sulfate. We report triple-oxygen isotope data for sulfate produced from pyrite oxidation in microbial and abiotic experiments, and from natural dissolved sulfate from the Rio Tinto, Spain, an acid mine drainage site. The oxygen isotope systematics of sulfate in these environments define a unique kinetic isotope effect associated with initial stage pyrite oxidation by Acidithiobacillus ferrooxidans that preserves >80 % oxygen from air-O 2 in sulfate. Unlike experiments, which evolve toward water-oxygen dominated sulfate on short time scales, Rio Tinto, Spain hosts a microbe rich environment with distinct geochemistry that maintains high O 2 -oxygen in sulfate. Therefore, in addition to containing isotopic records from water and air, sulfates can also contain a biosignature that is promising for understanding conditions on Mars and early Earth.

Earth and Planetary Science Letters↗
USGS2024

Uncertainty in ground-motion-to-intensity conversions significantly affects earthquake early warning alert regions

We examine how the choice of ground‐motion‐to‐intensity conversion equations (GMICEs) in earthquake early warning (EEW) systems affects resulting alert regions. We find that existing GMICEs can underestimate observed shaking at short rupture distances or overestimate the extent of low‐intensity shaking. Updated GMICEs that remove these biases would improve the accuracy of alert regions for the ShakeAlert EEW system for the West Coast of the United States. ShakeAlert uses ground‐motion prediction equations (GMPEs), which calculate spatial distributions of peak ground acceleration (PGA) and peak ground velocity (PGV) from earthquake source estimates, combined with GMICEs to translate GMPE output into modified Mercalli intensity (MMI). We find significant epistemic uncertainty in alert distances; near‐source MMI estimates from different GMICEs can differ by over 1 MMI unit, and MMI extents used for public EEW alerts can differ by hundreds of kilometers for larger magnitude earthquakes ( M ∼6.5+). We use a catalog of “Did You Feel It?” shaking reports to evaluate how well GMICEs predict observed shaking. Our preferred GMICE is the one that computes MMI using PGV for high intensities and transitions to using PGA for nondamaging intensities. These results motivate updating GMICE relationships more generally, including in ShakeMap applications.

The Seismic Record↗
USGS2024

Perspectives on transportable array Alaska background noise levels

Background seismic noise fundamentally sets a lower bound on our ability to record signals arising from earthquakes. The background noise spectrum at a station is a combination of cultural noise, ocean-generated microseism noise, intrinsic instrument self-noise, and the sensitivity of the instrument to nonseismic noise sources. The USArray-Transportable Array Alaska deployed 195 stations across Alaska and parts of Canada (Yukon, British Columbia, and Northwest Territories). These stations were all installed using similar techniques and made use of instruments with similar self-noise levels. As such, this network provides an opportunity to look at how geographic location influences seismic background. Using these broadband stations, we report background noise levels from 0.2 to 75 s period in six discrete bands. By constructing “noise maps,” we depict both spatial and temporal changes in the background noise field. Using these maps, combined with targeted analysis, we infer sources and contributing factors to noise levels in these different period bands. These include cultural noise, the formation of sea ice, seasonal changes in permafrost and wave activity in the Gulf of Alaska, and magnetic field variability. We use this study as an opportunity to review several previous studies examining seismic noise in Arctic regions.

Book chapter↗
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↗
USGS2026

Remote sensing enables basin-scale inventories of coal mine methane

Underground coal mines are important global sources of methane, but emission estimates are uncertain. We show that emission estimates for individual mines from aircraft remote-sensing surveys in the United States agree within 40% with direct measurements used for national emission reporting (IPCC Tier 3 estimate). Such direct measurements are unavailable in most countries, which rely on estimated emission factors (EFs) applied to coal-production rates. We find that EFs from IPCC Tier 1 and the Model for Calculating Coal Mine Methane (MC2M) methods overestimate U.S. emissions 3-fold due to incorrect dependence on mine depth. An IPCC Tier 2 method using measured basin-specific mine gas content agrees with direct emission measurements but does not account for gob well emissions and requires gas content data that are generally unavailable. We show that aircraft remote sensing for a small sample of mines can successfully estimate basin-specific EFs for ventilation shafts and gob wells, enabling estimates of basin- and national-scale emissions. We discuss how the method can be applied with satellite remote sensing to quantify coal emissions worldwide.

Alabama, Colorado, Kentucky, New Mexico, Ohio, Pen↗
USGS2025

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv↗
USGS2025

Over, under, and through: Hydrologic connectivity and the future of coastal landscape salinization

Seawater intrusion (SWI) affects coastal landscapes worldwide. Here we describe the hydrologic pathways through which SWI occurs - over land via storm surge or tidal flooding, under land via groundwater transport, and through watersheds via natural and artificial surface water channels—and how human modifications to those pathways alter patterns of SWI. We present an approach to advance understanding of spatiotemporal patterns of salinization that integrates these hydrologic pathways, their interactions, and how humans modify them. We use examples across the East Coast of the United States that exemplify mechanisms of salinization that have been reported around the planet to illustrate how hydrologic connectivity and human modifications alter patterns of SWI. Finally, we suggest a path for advancing SWI science that includes (a) deploying standardized and well-distributed sensor networks at local to global scales that intentionally track SWI fronts, (b) employing remote sensing and geospatial imaging techniques targeted at integrating above and belowground patterns of SWI, and (c) continuing to develop data analysis and model-data fusion techniques to measure the extent, understand the effects, and predict the future of coastal salinization.

Water Resources Research↗
USGS2025

Borehole geophysical time-series logging to monitor passive ISCO treatment of residual chlorinated-ethenes in a confining bed, NAS Pensacola, Florida

In-situ chemical oxidation (ISCO) is a common method to remediate chlorinated ethene contaminants in groundwater. Monitoring the effectiveness of ISCO can be hindered because of insufficient observations to assess oxidant delivery. Advantageously, potassium permanganate, one type of oxidant, provides the opportunity to use its strong electrical signal as a surrogate to track oxidant delivery using time-series borehole geophysical methods, like electromagnetic (EM) induction logging. Here we report a passive ISCO (P-ISCO) experiment, using potassium permanganate cylinders emplaced in boreholes, at a chlorinated ethene contamination site, Naval Air Station Pensacola, Florida. The contaminants are found primarily at the base of a shallow sandy aquifer in contact with an underlying silty-clay confining bed. We used results of the time-series borehole logging collected between 2017 and 2022 in 4 monitoring wells to track oxidant delivery. The EM-induction logs from the monitoring wells showed an increase in EM response primarily along the contact, likely from pooling of the oxidant, during P-ISCO treatment in 2021. Interestingly, concurrent natural gamma-ray (NGR) logging showed a decrease in NGR response at 3 of the 4 wells possibly from the formation of manganese precipitates coating sediments. The coupling of time-series logging and well-chemistry data allowed for an improved assessment of passive ISCO treatment effectiveness.

Florida↗
USGS2025

Kiloyear cycles of carbonate and Mg-silicate replacement at Von Damm hydrothermal vent field

The Von Damm vent field (VDVF) on the Mid-Cayman Rise in the Caribbean Sea is unique among modern hydrothermal systems in that the chimneys and mounds are almost entirely composed of talc. We analyzed samples collected in 2020 and report that in addition to disordered talc of variable crystallinity, carbonates are a major class of mineral at VDVF. The carbonate minerals include aragonite, calcite, magnesium-rich calcite, and dolomite. Talc and carbonate mineral textures indicate that, rather than replacing volcanic host rock, they precipitate from the mixing of hydrothermal fluids and seawater at the seafloor, occurring in chimneys and surrounding rubble. Alternating precipitation of this mineral assemblage is pervasive, with carbonate minerals typically being succeeded by talc, and with indications that in some cases talc and carbonate minerals replace one another. Stable carbon isotopic data indicate the carbonate minerals originate from the mixing of seawater and hydrothermal fluid, which is supported by U-Th data. Radiocarbon calcite ages and talc 234 U- 230 Th isochron ages indicate mineral ages spanning over thousands to tens of thousands of years. Analyses of these samples illustrate a dynamic system that transitions from carbonate-dominated to Mg-silicate−dominated precipitation over time scales of thousands of years. Our observations raise questions regarding the eventual fate of seafloor precipitates and whether carbonate and silicate minerals in such settings are sequestered and represented in the rock record.

Geology↗
USGS2024

A conceptual framework to assess post-wildfire water quality: State of the science and knowledge gaps

Wildfire substantially alters aquatic ecosystems by inducing moderate to catastrophic physical and chemical changes. However, the relations of environmental and watershed variables that drive those effects are complex. We present a Driver-Factor-Stressor-Effect (DFSE) conceptual framework to assess the current state of the science related to post-wildfire water-quality. We reviewed 64 peer-reviewed papers using the DFSE framework to identify drivers, factors, stressors, and effects associated with each study. A total of five drivers were identified and ranked according to their frequency of occurrence in the literature: atmospheric processes > fire characteristics > ecologic processes and characteristics > land surface characteristics > soil characteristics. Commonly reported stressors include increased nutrients, runoff, and sediment transport. Furthermore, although several different factors have been used at least once to explain water-quality effects, relatively few factors outside of precipitation and fire characteristics are frequently studied. We identified several gaps indicating the need for long-term monitoring, multi-factor studies, consideration of organic contaminants, consideration of groundwater, and inclusion of soil characteristics. This assessment expands on other reviews and meta-analyses by exploring causal linkages between influential variables and overall effects in post-wildfire watersheds. Information gathered from our assessment and the framework itself can be used to inform future monitoring plans and as a guide for modeling efforts focused on better understanding specific processes or to mitigate potential risks of post-wildfire water quality.

Water Resources Research↗
USGS2025

Forecast, monitor, adapt: A multi-agency strategy to protect people from postfire debris flows

In 2020, a wildfire burned across Glenwood Canyon in Colorado, USA. A history of postfire debris flows in the region and a hazard assessment for the burn area indicated that potentially life-threatening debris flows could be triggered by rainfall within months of a wildfire. As a result, four government agencies evaluated strategies to help mitigate hazards, including the loss of human life, that may be associated with debris-flow events. After the fire, 26 large debris flows occurred in the summer of 2021 and three sediment-laden flows occurred in the summer of 2023, but there were no major injuries or fatalities reported. We found that integrating hazard assessment/ forecasting, monitoring, and adaptation scenarios was a successful strategy for reducing postfire debris-flow risks to human life (including injuries and fatalities). Weather forecasts and estimates of debris-flow triggering rainfall thresholds, likelihood, and volume were used to anticipate the timing, location, and magnitude of debris-flow events. Rainfall monitoring and detailed recordkeeping of storms that triggered debris flows were used to validate and update debris-flow warning thresholds that varied with time following the wildfire. Although the governmental agencies working in this burn area had distinct and differing agency mandates, they were able to integrate information to reduce the risk of debris-flow events to human life.

Colorado↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2025

Statistical approaches for modeling correlated grade and tonnage distributions and applications for mineral resource assessments

Correlations between grade and tonnage exist in mineral resource data compiled from published reports, but they are not always addressed during quantitative assessment of undiscovered mineral resources. Failure to account for correlated grade and tonnage distributions can result in geologically unrealistic assessment results. Current software tools simulate univariate ore tonnage and multivariate resource grades of undiscovered deposits independently. As a result, analysts are forced to rely on ad-hoc solutions to minimize the correlation issues by: 1) creating subsets of data with restricted criteria; 2) truncating grade and tonnage distributions; and 3) testing model robustness using exploratory data analysis. While these methods represent pragmatic solutions, the statistical solutions presented here provide additional options to address real correlations in grade and tonnage data used for mineral resource assessments. We present a modified version of the MapMark4 package in R that introduces two alternatives for modeling grade and tonnage distributions, consisting of a multivariate solution that accounts for correlations between ore tonnage and metal grades and an empirical solution that utilizes simple random sampling with replacement to reproduce coupled grades and tonnages from the input data. We present simulations for contained ore and metal for three case studies representing tungsten skarn, komatiite-hosted nickel, and sediment-hosted carbonate amagmatic zinc-lead (Mississippi Valley-type) deposits. Employing the methods presented here yields quantitative mineral resource assessment results that more closely reflect the empirical distributions of grades and tonnages observed in nature and expands the applicability of these tools for ongoing critical mineral resource assessments.

Applied Computing and Geosciences↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2026

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗