USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Investigations Report”

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.

923 records · Page 29Linked to original sources
USGS2025

Examining the role of elevated and sustained strain in dynamically triggering earthquakes on the Anza section of the San Jacinto fault

Microearthquakes can be dynamically triggered in southern California by remote earthquakes. However, directly connecting dynamic triggering mechanisms with observational data remains challenging. One proposed failure mechanism suggests that both the amplitude and duration of cyclic fatigue caused by the passing seismic wave contribute to triggering occurrence. Here, we measure dynamic strains recorded by borehole strainmeters in the Anza section of the San Jacinto fault zone from 710 earthquakes that occurred over 300 km away between 2008 and 2017 to systematically investigate the role of elevated and sustained strain in controlling dynamic triggering. We design a suite of tests to evaluate whether specific amplitude thresholds and durations of strain can predict dynamic triggering cases. We further test whether the peak dynamic strain (PDS) can predict triggering occurrence in combination with the strain amplitude and duration. Based on these tests, there is no strain amplitude–duration threshold that can distinguish triggering occurrence in Anza. Dynamic triggering is more likely to occur if a remote earthquake causes a PDS above 100 nanostrain, though many cases were triggered at smaller PDSs. The lack of clear correlation between triggering and characteristics of the dynamic strain field suggests that the tested features of the incoming waves do not determine triggering occurrence and local fault conditions and slip processes are more important in controlling dynamic triggering in Anza.

California↗
USGS2025

Review and synthesis of the applications of machine learning to coalbed methane recovery

Over the last 30 years, a substantial literature has evolved on the use of machine learning (ML) to assess, predict, and improve the efficiency of coalbed methane (CBM) recovery. In the United States, the production of CBM declined as shale gas production matured, but CBM continues to be an important energy resource in other parts of the world. ML applications that have the potential to improve CBM reservoir management and production forecasts, and to increase exploration and operational efficiency, are still of significant interest. The integration of geostatistical techniques into the CBM ML applications has been largely absent but represents an opportunity for improvement. The literature demonstrates the widespread interest in, and applicability of, ML algorithms applied to CBM problems, and that they continue to result in improvements in predictive performance. However, (1) much of the research is more academic than operational, (2) many results are based on simulations, or small or proprietary datasets, (3) ML performance information can be inconsistent and sometimes entirely omitted, (4) most methodologies are unique to the specific CBM situation and likely not generalizable, (5) no standard data repositories are available to directly compare the performance of competing algorithms, and (6) the spatial component is often omitted. Finally, relatively new ML protocols involving causality analysis and reinforced learning, as well as hybrid workflows combining both supervised and unsupervised learning, are anticipated to dominate the future investigations. Integration of geostatistical and geospatial analysis with ML should enhance performance.

Book chapter↗
USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2025

Insights into widespread landsliding in southern Appalachia from Hurricane Helene

Between 23 and 27 September 2024, antecedent rain followed by Hurricane Helene produced one of the most damaging weather events in southern Appalachia history. The back-to-back storm events resulted in a maximum cumulative rainfall of 848 mm and hurricane-force wind gusts over 170 km/h in western North Carolina, eastern Tennessee, and southwestern Virginia. The resulting regional flooding, landslides, and tree blowdown caused over 100 fatalities, damaged or destroyed critical infrastructure and thousands of structures, and severed connectivity across the region. Over the next several weeks, a multi-agency landslide response produced a rapid hazard assessment and mapped 2217 landslides, 55% of which damaged infrastructure or property. Orographic uplift enhanced rainfall, resulting in concentrated landsliding along the ~250 km swath of the Blue Ridge escarpment in western North Carolina. Landslides initiated predominantly on windward-facing (southeast-facing) slopes, and localized clustering of initiation points indicated a strong influence of hillslope-scale meteorological and geomorphic factors. Many shallow landslides mobilized into larger, highly mobile, and damaging debris flows that graded into floods. Here, we put our preliminary observations in the context of historical storm-driven landslide events and open new avenues for investigating the nature and extent of landslides and their effects in southern Appalachia and similar environments.

Georgia, North Carolina, South Carolina, Tennessee↗
USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗
USGS2026

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗
USGS2025

Wetter winters, drier summers: Quantifying the change in hydrological response around the Puget Sound area using the wflow_sbm hydrological model and CMIP6 projections

Climate change is expected to impact hydrological regimes worldwide, including the Pacific Northwest of the United States. This study investigates how climate change will affect river discharge in the Puget Sound region of the State of Washington, with a focus on King and Pierce Counties. We simulated river discharge under historical and future conditions using the physically based, spatially distributed wflow_sbm hydrological model, which was calibrated and validated against U.S. Geological Survey discharge records. Future forcing was based on an ensemble of six high-resolution CMIP6 climate models, which were bias corrected using empirical quantile mapping. The results indicate a decrease in summer discharges (5–10%) and an increase in winter discharges (5–10%) across the study region. The high discharges (90th percentile) are projected to increase in winter, and the low discharges are projected to decrease in summer, due to shifts in precipitation regimes, snowpack hydrology, and evapotranspiration. However, variability between individual CMIP6 models often exceeds the magnitude of ensemble mean changes, underscoring substantial uncertainty in climate projections and the importance of including multiple climate models in climate change analysis. Furthermore, model consensus increased with elevation, which could be the result of the higher elevation areas being driven by less diverse hydrological processes. These findings highlight potential challenges for regional water management, ecosystem health, and flood risk mitigation in the Puget Sound region under future climate conditions.

EarthArXiv↗
USGS2025

Factors influencing landslide occurrence in low-relief formerly glaciated landscapes: Landslide inventory and susceptibility analysis in Minnesota, USA

In landscapes recently impacted by continental glaciation, landslides may occur where topographic relief has been generated by the drainage of glacial lakes and ensuing post-glacial fluvial network development into unconsolidated glacially derived sediments and exhumed bedrock. To investigate linkages among environmental variables, post-glacial landscape development, and landslides, we created a landslide inventory of nearly 10,000 landslides in five regions of the formerly glaciated low-relief state of Minnesota, USA. Multivariate logistic regression indicates the importance of slope angle, lithology, and the development of stream valleys to landslide distribution. Areas underlain by fine-grained glaciolacustrine and nearshore deposits that are incised by streams are particularly prone to shallow (<1-2 m depth) landslides. Landslides also occur in a wide range of glacial and fluvial deposits, and as rockfall in layered Paleozoic sedimentary rocks in central and southern Minnesota and Precambrian igneous and sedimentary rocks in northeastern Minnesota. Although no more than 1-2% of the studied regions are susceptible to landslides, they can pose risk to life and safety, damage infrastructure, and impact water quality. The combination of recently generated low-relief steep slopes, extensive unconsolidated sediments, and layered sedimentary bedrock make this formerly glaciated landscape more susceptible to landslides than current national-scale models indicate.

Minnesota↗
USGS2026

Preface to the focus section on intraplate earthquakes

More than a half century after plate tectonics provided an overarching framework to explain earthquakes along active plate boundaries, numerous theories have been proposed to explain where, why, and how often earthquakes occur well away from active plate boundaries, but a paradigm remains elusive. Even the classification of earthquakes away from active plate boundaries as "intraplate" raises issues, with potentially important distinctions between Stable Continental Regions and more actively deforming regions including passive margins and failed rifts. Some of the largest known intraplate earthquakes themselves remain enigmatic, having occurred before the modern instrumental era. Hazard assessments are often data-limited: low fault-slip rates relative to landscape modification rates result in poor discoverability of fault sources, challenging the characterization of source zones and earthquake recurrence; the completeness and homogenization of instrumental earthquake catalogs using uncertain magnitude conversions can lead to uncertainties in earthquake recurrence; and, limited strong-motion observations for large-magnitude events at near-source distances leads to uncertainties in the selection and development of ground-motion models for seismic hazard studies. Data from recent intraplate earthquakes around the world—from the moment magnitude M 7.7 2001 Bhuj, India, earthquake 25 years ago to the 2024 M 4.8 Tewksbury, New Jersey earthquake—have yielded both new insights and new questions. The papers in this special focus discuss many of the long-standing challenges involved with intraplate earthquake investigations and provide a snapshot of the state of the art with current research to advance our understanding.

Seismological Research Letters↗
USGS2024

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface↗
USGS2024

Groundwater-Surface water interactions research: Past trends and future directions

Interactions between groundwater and surface water sustain groundwater-dependent ecosystems and regulate river temperature and biogeochemical cycles, amongst many other processes. These interactions occur in freshwater environments including rivers, springs, lakes, and wetlands, and in coastal environments via tidal pumping, submarine groundwater discharge, and seawater intrusion. Here, we explore groundwater-surface water interactions research using bibliometric analyses of titles, abstracts, and keywords from 20,275 journal papers published between 1970 and 2023 extracted from Scopus. Analyses show that research into groundwater-surface water interactions is highly multi-disciplinary, with growing contributions from the social and biological sciences. The number of groundwater-surface water interactions papers is rapidly increasing with over 1200 papers published per year since 2020. Drawing on our data-driven approach and expert knowledge, we synthesise current research trends and identify critical future research directions. Despite the thousands of papers on groundwater-surface water interactions, important processes are still difficult to quantify or predict at meaningful spatial scales to inform water-resources management. We see benefits in future groundwater-surface water interactions research focusing on: (1) using new technologies including internet-of-things-based sensors, uncrewed vehicles, and remote-sensing approaches for data collection to inform groundwater-surface water interactions at large scales, (2) seeking approaches to upscale site-specific findings to better inform management, and (3) continuing the movement towards multi-disciplinary investigations to better inform the understanding of groundwater-surface water interactions and processes that will enable better management outcomes.

Journal of Hydrology↗
USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗
USGS2026

Stepovers and beyond: Structural control of The Geysers geothermal system and the broader Clear Lake region

Fault geometry exerts a first-order control on geothermal systems by governing stress localization, fracture development, and permeability, yet in complex fault networks or broader shear zones, the relative influence of individual geometric features is often difficult to resolve. In the northern California Coast Ranges, The Geysers geothermal field is commonly interpreted to occur within a releasing stepover, although no single, clearly defined stepover is identified in published studies. To investigate the structural controls on The Geysers and the broader Clear Lake region, a two-dimensional elastic boundary element model is developed to evaluate spatial patterns of dilational strain associated with progressively more complete fault geometries. Model results show that dilation in the region is not controlled by a single structure but instead reflects the combined effects of multiple interacting fault elements. Three primary controls are identified: (1) opposing bends in the regional strike-slip fault system, including a releasing bend along the Maacama fault; (2) the southern fault tip of the Collayomi fault, which generates a prominent dilational lobe beneath the southern Geysers; and (3) a releasing stepover between the Collayomi fault and the Geyser Peak–Mercuryville–Big Sulphur Creek fault system, inferred to collectively behave as a right-lateral shear zone bounding the western margin of The Geysers. Predicted dilational strain magnitudes are sufficient to localize permeability between faults. These results highlight that incorporating complete fault networks and bedrock geological mapping can enhance geothermal assessments and provide a transferable framework for evaluating structurally controlled permeability in tectonically active regions.

California↗
USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗
USGS2024

Managing basin-scale carbon sequestration: A tragedy of the commons approach

The Tragedy of the Commons is a well studied problem in the literature of ecology, economics, and environmental policy which illustrates the deleterious consequences of managing common pool resources when individual and social incentives are misaligned. In this work, we apply a simple model of carbon sequestration in a deep saline aquifer by two neighboring geologic CO 2 storage (GCS) operators to begin investigating if a Tragedy of the Commons framework applies to GCS. Specifically, we consider the pressure space as a “commons” because the injection by each firm at its own well increases the downhole injection pressure at both wells. We assume that a firm will decrease its injection rate if the downhole pressure at its well exceeds a predefined maximum (i.e., exceeds the “pressure limit”). With this assumption in place, we find that the same injection flowrates are optimal for both wells, regardless of whether they are owned by the same firm or competing firms. This suggests that GCS may not be best represented by a pure Tragedy of the Commons framework under our initial assumptions. However, there could be economic incentives or contractual obligations that may result in either or both GCS operators being unwilling to reduce their injection rates. Thus, we conclude the conference paper with a discussion of future extensions of our approach that may demonstrate closer alignment with the Tragedy of the Commons, including explicit definitions of pore-space rights, firm uncertainty regarding the parameters of the Theis equation, and the potential role of unitization.

Conference Paper↗
USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗
USGS2024

Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in-situ distribution of fault damage remain limited and along-strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub-bottom profiles. The combination of high-resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi-trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

California↗