Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geotechnical and Geoenvironmental Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,723 records · Page 8Linked to original sources
USGS2026

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States↗
USGS2024

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗
USGS2024

A new water temperature modeling approach to predict thermal habitat suitability for nonnative cichlids in Florida rivers

As global temperatures increase, the spatiotemporal arrangement of thermal habitats in Florida rivers may shift, creating the potential for greater dispersal and establishment of nonnative tropical freshwater fishes. To understand how water temperature changes may affect the spatial distribution of these nonnative species, more effective water temperature prediction models are necessary. Currently, most models employ either a generalized air–water temperature relationship or require expensive and complicated tools to measure hydrometeorological factors (e.g. groundwater input). Thus, we developed a novel modeling approach that is accurate, accessible, and cost-effective in allowing fisheries managers to project water temperatures in rivers across Central and North Florida. To characterize the potential for nonnative fishes to spread northward, we evaluated two hardy and abundant species currently found primarily in South Florida: Mayan Cichlid ( Mayaheros urophthalmus ) and Oscar ( Astronotus ocellatus ). Our results show an increase in thermally suitable winter days for both species in 10 of 11 rivers studied, consistent with predicted water temperature warming under 16 climate-change scenarios spanning different levels of air temperature warming (+1 °C, +2 °C, +3 °C, +4 °C) and precipitation/groundwater thermal sensitivity (0, 0.33, 0.66, 1). Considering resource limitations, fisheries managers can use our water temperature modeling approach to predict effects of climate change on Mayan Cichlid and Oscar survival, growth, and dispersal and take actions to manage potential northward movement of these species.

Florida↗
USGS2025

Systematic approach to prioritize wells for effective groundwater monitoring and management in the Arkansas Headwaters Basin, Colorado, USA

Study region The Arkansas Headwaters Basin, an intermountain basin in the Southern Rocky Mountains of North America. Study focus Our specific focus is choosing a set of wells to support a possible future regional groundwater-surface water model that would support water management. We present a three-step process using multiple criteria to score, predict, and choose prioritized wells that capture the full distribution of data including extremes. The three-step process provides accessible visualizations, fiscally efficient well prioritization, and screening useful for subsequent groundwater modeling. The novelty of the proposed methodology is the systematic approach integrating a scoring and a predictive approach to support a selection path. The systematic approach may be broadly adapted for other basins. New hydrological insights for the region Understanding regional hydrology hinges on efficient collection of hydrologic data that captures the relevant dynamics including extremes. The present study, a case study for a particular basin in the Southern Rocky Mountains, is the first use of a scripted (R software) strategy to select an economical and representative set of monitoring wells. Our findings suggest caution when using proximity as a proxy for correlation, because proximal wells in the same geologic formation and similar depths are not always correlated. In the Arkansas Headwaters Basin, subsurface geology may be less influential on groundwater elevations than broader hydrologic influences, such as regional drought.

Colorado↗
USGS2024

Cross-fade sampling: Extremely efficient Bayesian inversion for a variety of geophysical problems

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density function (PDF). Secondly, it is much faster to evaluate the posterior probability of each sample than conventional sampling methods for simulating Bayesian posterior PDFs. Conventional methods require evaluating the prior probability (which describes your a priori constraints) and data likelihood (which describes the fit between the observations and the predictions of the model) for each sample model. However, cross-fading does not require evaluating the data likelihood, meaning that ‘big data’ can be fit with zero additional computational cost. Further, the cross-fading approach can be used to calculate the marginal likelihood associated with a model design, facilitating model comparison and Bayesian model averaging. Topics covered in this paper include derivation of the cross-fade approach and how it can be used to simulate Bayesian posterior PDFs and compute the marginal likelihood, discussion of the class of problems to which cross-fading can be applied (with examples from earthquake statistics, earthquake ground motion modelling, volcanic eruption forecasting, and finite fault slip modelling), demonstration of efficiency relative to existing sampling methods and discussion of how cross-fading can be used to account for prediction errors (i.e. epistemic errors) as part of the geophysical inverse problem.

Geophysical Journal International↗
USGS2026

Decadal-scale characteristics of natural and anthropogenic dune morphology along North Carolina Barrier Islands (SE Atlantic Coast)

Dunes offer protection to coastal communities from storms, yet quantitative research characterizing human influences on dune systems across decades and 10–100 km spatial scales constitutes a known data gap. This study analyzed a multidecadal archive of dune crest elevations and positions to map mesoscale spatiotemporal trends in foredune morphology across the barrier island coastline of North Carolina, SE Atlantic coast. Results are summarized within four categories of coastline representing a range of coastal management intensities, integrating beach width and dune lateral accommodation space to assess the comparative influence of differing intensities of anthropogenic intervention on dune morphology. These categories account for the presence of structures, roads, and alongshore connectivity to anthropogenic sediment sources (beach nourishment). These categories show distinct dune and beach system characteristics. Dunes are largest in the areas with the longest legacy of coastal management. As categories decrease in development intensity, dune and beach system characteristics generally depart from the large, static dunes and narrow beaches consistent with hold-the-line management. Where small or transgressive dunes exist on heavily developed coasts, they coincide with relatively strong hydrodynamic forcing or sediment-limiting geology. There is evidence of indirect coastal management influence on undeveloped barrier islands with alongshore connectivity to developed barriers. While this analysis was conducted at the mesoscale, the results provide researchers with the framework to make insights at smaller scales than those demonstrated in this study.

North Carolina↗
USGS2026

Assessing future hydrologic extremes using an integrated hydrology and river operations model in the Russian River watershed

Study region The Russian River watershed, situated in coastal, northern California, experiences hydrologic extremes, including periodic droughts and flooding. Water managers are working to maintain sustainable water supplies and environmental flows, while mitigating flood risks. Study focus This paper introduces an integrated hydrology and river operations model for the Russian River watershed. This model is distinct from models in previous studies because it represents surface-groundwater interactions and uses climate forcings to estimate dynamic water use demands that are superimposed onto both reservoir operations and water supply constraints. The model was used to examine three historical (1990–2015) and eight future (2016–2099) water use and climate change scenarios. New hydrological insights for the region The direct connection between streams and aquifers facilitated both annual aquifer replenishment by high winter streamflows and streamflow depletion by groundwater wells (19 % of pumped groundwater in alluvial aquifers from stream leakage) during critical low flow periods. Simulated streamflow changes included 59 % longer and 54 % more severe streamflow droughts, 26 % lower seasonal low streamflows, and up to 125 % higher peak streamflows, averaged over future climate and water use scenarios, suggesting increased future flood and water availability risks. Results showed the importance of reservoir operations for mitigating the impacts of increased hydroclimatic volatility, despite a decrease in reservoir reliability at Lake Mendocino, suggesting that reservoir management may be used to decrease future risks.

California↗
USGS2026

Evaluating the U.S. Geological Survey’s earthquake shaking hazard forecasts and their implications on seismic risk

We analyze the last six update cycles of the long-term probabilistic earthquake shaking hazard forecast from the U.S. Geological Survey (USGS) and discuss the changes in hazard estimates from the 1996 to the latest 2023 update of the National Seismic Hazard Model (NSHM) for the conterminous U.S. We summarize how our understanding of earthquake hazards has evolved over the last two decades and quantify the implications of changing hazards on people, buildings, and lifeline infrastructure. The net effect of changes in hazard estimates, with reference to the mean hazard from the 2023 NSHM update that has 2% probability of exceedance in 50 years, suggests an overall increase in total geographic areas facing very strong shaking levels (modified Mercalli Intensity (MMI) of VII or more) when compared with any of the previous cycles. This increase puts ~159.2 million people (an increase from 111.8 million when using the 2018 NSHM), ~47.3 million residential buildings (an increase from 32.5 million with the 2018 NSHM), and ~$25.6 trillion of economic exposure of all buildings (an increase from $17.4 trillion when using the 2018 NSHM) at risk from earthquakes. We demonstrate that small changes in hazard estimates do not necessarily imply small changes in risk estimates. Therefore, changes in risk estimates can be used to highlight key USGS NSHM updates for improving risk mitigation efforts.

Earthquake Spectra Journal↗
USGS2025

Hydrologic variability and groundwater age of springs in eastern Oregon and northern Nevada, USA

The ecological importance of springs in semiarid regions is far greater than their small size and sparse distribution, yet little is known about the hydrologic functioning of these systems. During 2016–22, 261 springs were visited in the volcanic terrane of eastern Oregon and northern Nevada. When conditions were suitable, measurements of discharge, water temperature, and specific conductance were made, and samples for the analysis of carbon-14, tritium, and water stable isotopes (WSI) were collected. A subset of 60 springs was revisited during different seasons in the same year and during the dry season in multiple years to evaluate variability in discharge, chemistry, and groundwater age. Specific conductance and WSI varied considerably among springs across the study area but were unexpectedly stable across seasons and years at individual springs. Seasonal and interannual variability in spring discharge was related to the residence time of the discharging groundwater. Springs discharging older groundwater (10 3 –10 4 years) had significantly less variability in their discharge compared to springs discharging younger groundwater (10 0 –10 1 years). Variability among springs discharging younger groundwater included cessation of late-summer discharge at 18 % of the repeat-visit springs. A logistic regression model predicted the age of discharging spring water with 89 % accuracy using only the spring latitude, longitude, elevation, and δ 2 H value. This study framework provides a simple, inexpensive, and robust method to provisionally assess the hydrologic behavior of springs having little or no prior information in understudied, semiarid regions across the globe.

Nevada, Oregon↗
USGS2026

Ringtail (Bassariscus astutus) survival in southwestern Oregon

The effective conservation and management of small carnivore populations requires understanding species’ life-history traits and identifying important vital rates that drive population trajectories. However, many of these species are rare or elusive and of state or federal conservation concern, and demographic information is often lacking and difficult to obtain. At the northern limit of their range in Oregon, USA, ringtail ( Bassariscus astutus ) occupy mid-elevation forests and are a species of conservation concern because of their limited distribution and suspected low density. We initiated a radio-telemetry study in 2020 to estimate monthly and annual survival of ringtail in southwest Oregon. We monitored 26 ringtail from November 2020 to October 2022 and estimated survival rates using a known-fate framework and Program MARK. Model-averaged monthly survival estimates ranged from lows of 0.963 (SE = 0.022, 95% CI = 0.887–0.988) in spring (Feb–May) of 2021 to highs of 0.980 (SE = 0.018, 95% CI = 0.889–0.997) in summer (Jun–Oct) of 2022. Model-averaged estimates of annual survival were 0.695 (SE = 0.176, 95% CI = 0.310–0.920) during 2020–2021 and 0.728 (SE = 0.168, 95% CI = 0.336–0.934) during 2021–2022. Predation was the leading cause of mortality, but notably, no mortalities could be attributed to avian predation. Survival rates were much higher within forested landscapes in Oregon than for populations in the southwestern United States, which may have implications for regional conservation and management strategies.

Oregon↗
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

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2025

Per- and polyfluoroalkyl substances (PFAS) mass flux and mass balance at an aqueous film-forming foam release site in semiarid eastern New Mexico, USA

Passive flux meters (PFMs) directly measure groundwater chemistry mass flux and Darcy flux, providing insight into contaminant source-zone architecture and transport properties. This study uses PFMs to characterize PFAS flux in groundwater at a semiarid site with a thick (greater than 90-m) unsaturated zone where groundwater has been contaminated with per - and polyfluoroalkyl substances (PFAS) related to the use of aqueous film-forming foam (AFFF) for fire training and fire suppression. PFAS mass discharge (PFAS mass flux integrated over a control plane) in groundwater downgradient from several PFAS release areas is calculated using PFM results. In groundwater downgradient from fire-training areas, total PFAS mass discharge (summed across 14 compounds) was estimated to be between 6.0 and 31 g per day in 2020 and between 5.9 and 23 g per day in 2021. Site-specific documentation, generic information on AFFF properties, and literature values of PFAS concentration in AFFF are used to estimate site-specific PFAS-application rates at fire-training areas. These PFAS-application rates are compared to groundwater PFAS-discharge rates. Results suggest that transformation processes (exact pathways unknown) have led to increased discharge of measured PFAS in groundwater relative to initial AFFF formulations. The mass balance approach has broad applicability as a high-level approach that can provide insight into PFAS transport at AFFF sites.

New Mexico↗
USGS2026

Evaluating the effects of a multi-modal deterrent on upstream passage by invasive carp and native species

Invasive carp, including Silver Carp ( Hypophthalmichthys molitrix ) and Grass Carp ( Ctenopharyngodon idella ), are increasing their range within the Mississippi River Basin, threatening native fishes and ecosystem function. We conducted a 3-year field test of a multi-modal deterrent that uses sound, bubbles, and light to test its ability to impede upstream passage by invasive carp. The test was performed at a lock and dam on the Cumberland River, Kentucky, USA, using acoustic telemetry to track the movements of Silver Carp, Grass Carp, and several native fish species. When the deterrent was operating, the estimated probability of upstream lock passage was 53% lower (95% CI: 41%–63%) for Silver Carp but 74% higher (95% CI: 3%–190%) for Grass Carp compared to times when the deterrent was not operating. Paddlefish ( Polyodon spathula ) passage was reduced by 43% (95% CI: 71% reduction to 9% increase), although this was not statistically significant. We detected no significant effects on Freshwater Drum ( Aplodinotus grunniens ) or Smallmouth Buffalo ( Ictiobus bubalus ). However, sample sizes were small for these native species, limiting our ability to identify effects. Water temperature, tailwater elevation, and vessel lockages affected passage rates, although this varied among species.

Kentucky↗
USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗
USGS2025

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal↗
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↗