Geology topics

Andrew J. Barbour

Publications and source records attributed to Andrew J. Barbour.

At least 19 recordsLinked to original sources

USGS2026

Rapid earthquake magnitude classification via P-wave strains from borehole strainmeters and Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) offers a promising approach for earthquake early warning (EEW) in settings where seismic networks are costly to maintain. By repurposing fiber-optic cables as dense strainmeter arrays, DAS enables real-time earthquake detection wherever those fibers are accessible. However, poor azimuthal coverage and challenges in estimating magnitude from strain measurements remain key hurdles in applying for earthquake monitoring. Here, we develop a machine learning method to distinguish large (M≥5.4) earthquakes from smaller ones within the first 4 seconds of a strain waveform after a P-wave arrival without determining location. Using ensemble decision tree models trained on borehole strainmeter data (3.5≤M≤7.1) and tested on onshore DAS waveforms (including the 2024 M7 Offshore Cape Mendocino earthquake), we find that low-frequency (0.2–0.5 Hz) continuous wavelet transform coefficients are the strongest predictors of magnitude, in addition to strain amplitude. Both DAS and borehole strainmeters effectively capture long-period strain signals, making these findings valuable for EEW systems. Our method shows high precision compared to the real-time EEW system, ShakeAlert®, supporting the position that DAS is a viable technology for earthquake monitoring and magnitude classification.

California↗

USGS2026

Compact seismicity bursts have different characteristics from regional seismicity

Earthquakes tend to cluster, developing into sequences driven by stress perturbations and transient fault-zone processes. Depending on the driving process, earthquake sequences show differing behaviors. This variability challenges our ability to observe or distinguish these driving processes in high resolution. Here we systematically identify seismicity bursts throughout southern California using new statistical methods and examine their causes with other independent geophysical observations. Seismicity bursts are defined as statistically significant seismicity rate anomalies. Our approach has the advantage of being data driven, depending on few earthquake occurrence assumptions. We find abundant seismicity bursts across southern California, most frequently occurring along the San Jacinto Fault and in the Salton Sea and Coso geothermal fields. These seismicity bursts are highly compact in space and time, often encompassed by a 5 km radius and have durations less than 10 hr. Many of the seismicity bursts have their largest earthquake near the beginning of the sequence, but the largest earthquake is usually not an obvious mainshock. We utilize a variety of independent geophysical data sets to analyze the characteristics of the seismicity bursts, finding that the seismicity bursts have low b-values, low spectral stress drops, and varied stress ratios compared to regional seismicity. These differences suggest that seismicity bursts are driven by transient processes acting frequently across fault networks.

California↗

USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗

USGS2025

S/P amplitude ratios with Distributed Acoustic Sensing and application to earthquake focal mechanisms

Distributed acoustic sensing (DAS), which transforms a fiber optic cable into an array of high frequency strainmeters, has the potential to help us characterize earthquakes with a dense sampling of measurements. While earthquake focal mechanisms are frequently determined using P -wave polarities and S / P amplitude ratios with inertial seismometers, the dense sampling of DAS over potentially large portions of the focal sphere can aid our solutions. Here, we consider ∼200 regional earthquakes and thousands of S / P measurements on a DAS cable colocated with a network of inertial seismometers near Arcata, California. We demonstrate the S / P ratio measurements made on DAS are similar to those made on the inertial seismometers and can be used to constrain focal mechanism solutions.

California↗

USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗

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↗

USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗

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

Fiber-imaged supershear dynamics in the 2024 Mw 7 Mendocino Fault earthquake

Fault structure and rupture physics are deeply intertwined, and observations of this coupling are critical for understanding earthquake behavior. Rupture propagation is observable at fine scales using dense seismic networks. Fiber-optic sensing allows for long-term deployments of ultradense arrays that enable high-resolution measurements of infrequent, large earthquakes. We recorded the 2024 moment magnitude ( M w ) 7 Mendocino Fault earthquake with a nearby fiber-optic array and imaged its behavior with seismic beamforming. The rupture propagated to the east at subshear velocity; stagnated near the Mendocino Triple Junction, a zone of structural complexity; and subsequently transitioned to supershear velocity. The correlation between source physics and structure shows how lithospheric heterogeneity affects first-order characteristics of earthquake ruptures. Our results also demonstrate the potential for fiber-optic sensing to improve real-time estimation of key parameters for early warning.

California↗

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↗

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

Relatively stable pressure effects and time-increasing thermal contraction control Heber geothermal field deformation

Due to geological complexities and observational gaps, it is challenging to identify the governing physical processes of geothermal field deformation including ground subsidence and earthquakes. In the west and east regions of the Heber Geothermal Field (HGF), decade-long subsidence was occurring despite injection of heat-depleted brines, along with transient reversals between uplift and subsidence. These observed phenomena contradict current knowledge that injection leads to surface uplift. Here we show that high-yield production wells at the HGF center siphon fluid from surrounding regions, which can cause subsidence at low-rate injection locations. Moreover, the thermal contraction effect by cooling increases with time and eventually overwhelms the pressure effects of pressure fluctuation and poroelastic responses, which keep relatively stable during geothermal operations. The observed subsidence anomalies result from the siphoning effect and thermal contraction. We further demonstrate that thermal contraction dominates long-term trends of surface displacement and seismicity growth, while pressure effects drive near-instantaneous changes.

California↗

USGS2024

Induced seismicity strategic vision

Executive Summary The U.S. Geological Survey has a long history of contributions to the understanding and resolution of various scientific questions related to earthquakes associated with human activities, referred to as induced seismicity. Work started with the Rocky Mountain Arsenal studies in the 1960’s (Healy and others, 1968) when it was first discovered that fluid waste-disposal operations can cause earthquakes. U.S. Geological Survey work on induced seismicity continued in the intervening years but expanded dramatically in the early 2010s. Disposal of large volumes of wastewater, a byproduct of oil and gas production, into the subsurface caused an exponential increase in earthquakes in the central United States, including earthquakes that caused damage to buildings and infrastructure in nearby communities. Established in 2012 within the Earthquake Hazards Program (EHP), the Induced Seismicity Project (ISP) examines earthquakes that are induced by human activities to assess and mitigate the hazards associated with induced earthquakes as well as understand the conditions and processes controlling induced and natural earthquake generation and recurrence. The ISP examines instances of suspected induced earthquakes, in real-time and retrospectively, to assess the probabilistic hazard and investigate possibilities for reducing that hazard, something currently not possible with natural earthquake activity. The ISP has many synergies with work across multiple areas of EHP, including earthquake monitoring, hazard mitigation, and fundamental research into earthquake processes. The overarching questions that guided the development of our research strategies and approaches are as follows: How can seismicity that is induced by industrial activities be better distinguished from natural processes, particularly in regions with higher rates of natural seismicity? Can faults that may be at increased hazard for induced seismicity be identified based on existing fault structures, past earthquake activity, local stress conditions, or other information? What geologic and operational conditions control the occurrence, spatial and temporal evolution, maximum magnitude, ground shaking, and other characteristics of induced seismic sequences? What roles do surface and subsurface deformation, aseismic slip, and other processes play in the occurrence of injection-induced earthquakes? How can induced seismicity be better forecasted to support the Nation’s energy production and climate mitigation goals while limiting related earthquake hazards? In this report, background on induced seismicity and an overview of the state of knowledge is provided in three topical task areas: (1) oil and gas, (2) geothermal, and (3) geologic CO2 sequestration. This report then presents the EHP’s goals, strategies, and approaches—primarily focused on work in the ISP—for understanding the conditions that lead to induced seismicity, mitigating the associated hazards, and assisting the Nation in meeting its future needs in energy production and mitigation of climate change. Finally, the key limitations and opportunities for innovation are outlined at the end of the report.

Circular↗

USGS2024

Telecommunications fiber for sensing earthquake aftershocks: Progress and hurdles

Aftershocks offer valuable clues to earthquake behavior. The challenge: quickly deploying sensors to capture the early details of earthquake ruptures within the zone of aftershocks. Telecommunication fibers might be an answer, providing denser networks in otherwise difficult areas, potentially faster than traditional methods.

Conference Paper↗

USGS2022

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗

USGS2022

Very low frequency earthquakes in between the seismogenic and tremor zones in Cascadia?

Megathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.

Oregon, Washington↗

USGS2022

Fast rupture of the 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California dynamically triggers seismicity in California

In the Gulf of California, Mexico, the relative motion across the North America-Pacific boundary is accommodated by a series of marine transform faults and spreading centers. About 40 M>6 earthquakes have occurred in the region since 1960. On 3 August 2009, an Mw 6.9 earthquake occurred near Canal de Ballenas in the region. The earthquake was a strike-slip event with a shallow hypocenter that is likely close to the seafloor. In contrast to an adjacent M7 earthquake, this earthquake triggered a ground-motion-based earthquake early warning algorithm being tested in southern California (∼600 km away). This observation suggests that the abnormally large ground motions and dynamic strains observed for this earthquake relate to its rupture properties. To investigate this possibility, we image the rupture process and resolve the slip distribution of the event using a P-wave back-projection approach and a teleseismic, finite-fault inversion method. Results from these two independent analyses indicate a relatively simple, unilateral rupture propagation directed along-strike in the northward direction. However, the average rupture speed is estimated around 4 km/s, suggesting a possible supershear rupture. The supershear speed is also supported by a Rayleigh wave Mach cone analysis, although uncertainties in local velocity structure preclude a definitive conclusion. The Canal de Ballenas earthquake dynamically triggered seismicity at multiple sites in California, with triggering response characteristics varying from location-to-location. For instance, some of the triggered earthquakes in California occurred up to 24 hours later, suggesting that nonlinear triggering mechanisms likely have modulated their occurrence.

Baja California↗

USGS2022

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗