Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Article”

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.

979 records · Page 13Linked to original sources
USGS2025

High-precision earthquake catalog for Minto Flats fault zone, central Alaska, reveals complex and conjugate faulting

The Minto Flats fault zone (MFFZ) in central Alaska is a left‐lateral strike‐slip fault system situated between the continental‐scale right‐lateral Denali and Kaltag‐Tintina faults. The MFFZ has the potential to generate magnitude 7 earthquakes, and it hosted a magnitude 6 earthquake in 1995. It has also produced exotic events, such as very‐low‐frequency earthquakes and nucleation signals. We use network‐matched filtering and relative earthquake relocation techniques to derive a detailed catalog of earthquake locations for the MFFZ. The catalog spans from August 2014 to December 2019, a time period including 13 temporary seismic stations in the region. Our results provide the most complete catalog for the MFFZ and include deeper events, clusters of shallow seismicity, and a complex and segmented fault structure not observed in the original regional catalog. We document right‐lateral strike‐slip faulting, conjugate to the main northeast‐striking left‐lateral faults of the MFFZ. Below Nenana basin, the relocated seismicity reveals northwest‐dipping left‐lateral faults, supporting the inference that deep crustal active faulting is associated with recent basin deformation.

Alaska↗
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↗
USGS2025

Divide and conquer: Separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process – the probability of the existence of a seismic phase (detection probability) and the probability of correctly identifying the phase arrival time (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labeled phase arrival time, and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks) and the height of the peaks are taken as “probability” of the picks. Here, we show that this “probability” represents neither the detection nor the timing probabilty because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multi-class classifier to obtain a probability mass function to address the timing probability. This new approach makes the deep learning-based phase picking process more interpretable and gives us options to logically control seismic monitoring workflows.

Geophysical Journal International↗
USGS2024

Mitigating climate change by abating coal mine methane: A critical review of status and opportunities

Methane has a short atmospheric lifetime compared to carbon dioxide (CO 2 ), ∼decade versus ∼centuries, but it has a much higher global warming potential (GWP), highlighting how reducing methane emissions can slow the rate of climate change. When considering the contribution of greenhouse gas (GHG) emissions to current global warming (2010–2019) relative to the industrial revolution (1850–1900) levels, methane contributes 0.5 °C or ∼ a third of the total. The most recent post-2023 global estimates of methane emissions by bottom-up (BU) and top-down (TD) approaches for the coal mining sector are in the range of ∼41 ± 3 Tg yr −1 and 33 ± 5 Tg yr −1 , respectively. This divergence, notwithstanding overlapping confidence intervals, is a result of differences between applied TD global inversion models and BU emission inventories. Further research can help to better refine emissions from the various contributing coal mine methane (CMM) emissions sources. The coal mining sector accounts for over 10 % of global anthropogenic methane emissions. The contribution of CMM emissions to the global budget have increased since 2000, although upward and downward regional trends have been observed.

International Journal of Coal Geology↗
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↗
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↗
USGS2025

Disturbance is the primary determinant of food chain length when the top predator is constant

Food chain length (FCL) is a primary determinant of food web structure and is hypothesized to be influenced by habitat size, productivity, and disturbance. Understanding the environmental characteristics that determine food chain length can assist in understanding how food webs may be impacted due to changes in habitats and environmental characteristics. This study examines the impact of hydrologic disturbance on stream food webs when the top predator is constant. We analyzed FCL in less disturbed groundwater flashy streams and more disturbed runoff flashy streams using stable isotopes. Despite no difference in species richness or fish density, food chains in more disturbed streams had a lower FCL compared to food chains in more stable streams. Assemblage analysis showed that flow regime and drainage area significantly impacted individual species abundances. The more disturbed runoff flashy streams had higher proportions of primary consumer fish, such as the algivorous Campostoma sp. (Stonerollers), which likely drives the reduced FCL. Drainage area and land cover had non-significant relationships with FCL. Shifting community structure due to hydrologic variability likely leads to differences in diet of Micropterus dolomieu (Smallmouth Bass), and thus a difference in FCL.

Arkansas, Missouri, Oklahoma↗
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↗
USGS2026

A geographic dataset of rocky reefs habitat areas of particular concern for the United States West Coast

The United States National Marine Fisheries Service determines “essential fish habitat (EFH)” for federally managed species in coordination with regional fishery management councils, considers adverse effects to those habitats, and provides information to further habitat conservation and enhancement. Identifying discrete subsets of EFH as “habitat areas of particular concern (HAPC)” can help focus conservation, management, and research efforts. In 2006, the Pacific Fishery Management Council designated rocky reefs along the United States (U.S.) West Coast as HAPCs for groundfishes because of their ecological significance, sensitivity to human impacts, and relative rarity. To better understand where rocky reefs occur, we (1) located rocky reef areas that were not included in the 2006 rocky reef dataset, and (2) incorporated best available data into a refined geographic dataset that enables visualization. Our update shows that rocky reefs are distributed throughout the U.S. West Coast continental margin, are patchier than previously known, and comprise 8% of the extent of all data inputs. This updated dataset will inform resource management decisions in coastal and marine environments.

California, Oregon, Washington↗
USGS2026

New constraints on northeast Seattle basin structure from converted seismic waves

The Seattle basin is a deep sedimentary basin in the Seattle–Bellevue, Washington metropolitan area within the Puget Lowland of Washington State. We determine the structure of a portion of the basin and the underlying basement using analysis of P waves converted from direct S incident from below. A deep local crustal event beneath Monroe, about 35 km northeast of Seattle, was recorded by a 100‐station nodal array deployed in 2019. The event produced a variety of coherent seismic phases, including converted waves from the sediment—basement boundary, internal structure within the basin, and additional crustal discontinuities. Using observed Sp converted waves, we apply an adjoint‐based full waveform inversion (FWI) method to determine the amplitude and extent of seismic discontinuities at depth. We find the strongest source of converted waves for this event lies ∼6 to 7 km depth below northern Lake Washington, interpreted to be the local depth to basement rock. The newly imaged shallow basement structure may be part of a deformation zone associated with the Siletzia eastern boundary. Our results highlight the utility of converted seismic waves recorded by a dense array, combined with an FWI method, to illuminate crustal structure.

Washington↗
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↗
USGS2008

Geology and geochemistry of Jurassic plutonic rocks, Baboquivari Mountains, south-central Arizona

Among the plutons of the Jurassic magmatic arc segment in southern Arizona, northern Sonora, and southern California, two distinctive rock types dominate, occurring together in numerous mountain ranges: porphyritic, titanite-bearing hornblende-biotite granodiorite; and biotite leucogranite. These characteristic rocks are particularly well preserved and well exposed around Kitt Peak, in the Baboquivari Mountains, southern Arizona. The Middle to early Late Jurassic Kitt Peak Plutonic Suite (KPPS) comprises three units: Aguirre Peak Quartz Diorite (APQD;<170 Ma, >=165 Ma, U-Pb), chiefly hornblende mesodiorite; Kitt Peak Granodiorite (KPGD; 165 Ma), hornblende-biotite granodiorite and monzogranite; and Pavo Kug Granite (PKG; 159 Ma), equigranular biotite leucogranite. This range also hosts a fourth, regionally unusual, Jurassic plutonic unit: the Baboquivari Peak Perthite Granite, related to post-arc, Late Jurassic crustal extension. Binary variation diagrams, incompatible element patterns, and REE spectra indicate APQD and KPGD are consanguineous. APQD basaltic mesodiorite evidently approximates the parental magma for the APQD–KPGD series. More evolved rocks of this series evidently formed by combinations of fractional crystallization and mixing with or assimilation of crustal material. PKG is less directly related to APQD and KPGD than these two are to one another.

Arizona, Sonora↗
USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗
USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗
USGS2026

A 481 m-high landslide-tsunami in a cruise ship-frequented Alaska fjord

Early in the morning of 10 August 2025, a >64 × 10 6 –cubic meter landslide struck Tracy Arm fjord in Alaska. The landslide was preconditioned by glacial retreat caused by climate change. The resulting 481-meter runup megatsunami followed an initial 100-meter-high breaking wave traveling at >70 meters per second. The landslide was preceded by several days of microseismicity, which increased in rate and magnitude until ~1 hour before failure. The landslide produced globally observed long-period seismic waves equivalent in size to a moment magnitude 5.4 earthquake. A long-period (~66 second) global seismic signal, produced by a landslide-induced seiche trapped within the fjord, persisted for up to 36 hours, the second time a days-long seiche had thus been observed. With fjord regions increasingly visited by cruise ships, and climate change making similar events more likely, this unanticipated, near-miss event highlights the growing risk from landslides and tsunamis in coastal environments.

Alaska↗
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗
USGS2024

Reduced injection rates and shallower depths mitigated induced seismicity in Oklahoma

The proximity of wastewater disposal to the Precambrian basement is a critical factor influencing induced earthquake rates in the Central United States, but the impact of reducing injection depths has not been widely demonstrated. Beginning in 2015, state regulatory efforts in Oklahoma and Kansas mandated that wells injecting into the lower Arbuckle Group, a basal sedimentary unit, be backfilled with cement (i.e. “plugged back”) so that they inject into shallower formations. This plug back activity gives us a unique opportunity to investigate the relationship between injection depth and induced seismicity rate. To evaluate the impact that decreased injection rates and plug backs had on the seismicity rates, we create a suite of rate-state earthquake models. Observed seismicity rates are best fit when only lower Arbuckle volumes are considered, suggesting the lower Arbuckle injectors were primarily responsible for the seismicity and that plug backs were effective at isolating the injected volumes to shallower formations. Our models demonstrate that if these wells had not been plugged back, seismicity rates would be multiple times larger than they are today. We find that the combination of well plug backs and injection volume decreases can be effective strategies for reducing induced seismicity rates.

Oklahoma↗
USGS2025

Analyzing multi-year nitrate concentration evolution in Alabama aquatic systems using a machine learning model

Rising nitrate contamination in water systems poses significant risks to public health and ecosystem stability, necessitating advanced modeling to understand nitrate dynamics more accurately. This study applies the long short-term memory (LSTM) modeling to investigate the hydrologic and environmental factors influencing nitrate concentration dynamics in rivers and aquifers across the state of Alabama in the southeast of the United States. By integrating dynamic data such as streamflow and groundwater levels with static catchment attributes, the machine learning model identifies primary drivers of nitrate fluctuations, offering detailed insights into the complex interactions affecting multi-year nitrate concentrations in natural aquatic systems. In addition, a novel LSTM-based approach utilizes synthetic surface water nitrate data to predict groundwater nitrate levels, helping to address monitoring gaps in aquifers connected to these rivers. This method reveals potential correlations between surface water and groundwater nitrate dynamics, which is particularly meaningful given the lack of water quality observations in many aquifers. Field applications further show that, while the LSTM model effectively captures seasonal trends, limitations in representing extreme nitrate events suggest areas for further refinement. These findings contribute to data-driven water quality management, enhancing understanding of nitrate behavior in interconnected water systems.

Alabama↗