skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Raton Basin Induced Seismicity Is Hosted by Networks of Short Basement Faults and Mimics Tectonic Earthquake Statistics
Abstract The Raton Basin has been an area of injection induced seismicity for the past two decades. Previously, the reactivated fault zone structures and spatiotemporal response of seismicity to evolving injection have been poorly constrained due to sparse publicly available seismic monitoring. The application of a machine‐learning phase picker to 4 years of continuous seismic data from a local array enables the detection and location of ∼38,000 earthquakes. The events from 2016 to 2020 are ∼2.5–6 km below sea level and range from ML < −1 to 4.2. Most earthquakes occur within previously identified ∼N‐S zones of seismicity, however our new catalog illuminates that these zones are composed of many short faults with variable orientations. The two most active zones, the Vermejo Park and Tercio zones, are potentially linked by small intermediate faults. In total, we find ∼60 short (<3 km long) basement faults with strikes from WNW to NNE. Faulting mechanisms are predominantly normal but some variability, including reverse dip‐slip and oblique‐slip, is observed. The Trinidad fault zone, which previously hosted a Mw5.3 earthquake in 2011, is quiescent during 2016–2020, likely in response to both slow accumulation of tectonic strain after the 2011 sequence, and the significant decrease (80% reduction) in nearby wastewater injection from 2012 to 2016. Unlike some other regions, where induced seismicity was triggered in response to higher injection rates, the Raton Basin's frequency‐magnitude and spatiotemporal statistics are not distinguishable from tectonic seismicity. The similarity suggests that seismicity in the Raton Basin is predominantly releasing tectonic stress.  more » « less
Award ID(s):
1554908
PAR ID:
10372288
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
11
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract We report on the tectonic framework, seismicity, and aftershock monitoring efforts related to the 31 March 2020 Mw 6.5 Stanley, Idaho, earthquake. The earthquake sequence has produced both strike-slip and dip-slip motion, with minimal surface displacement or damage. The earthquake occurred at the northern limits of the Sawtooth normal fault. This fault separates the Centennial tectonic belt, a zone of active seismicity within the Basin and Range Province, from the Idaho batholith to the west and Challis volcanic belt to the north and east. We show evidence for a potential kinematic link between the northeast-dipping Sawtooth fault and the southwest-dipping Lost River fault. These opposing faults have recorded four of the five M≥6 Idaho earthquakes from the past 76 yr, including 1983 Mw 6.9 Borah Peak and the 1944 M 6.1 and 1945 M 6.0 Seafoam earthquakes. Geological and geophysical data point to possible fault boundary segments driven by pre-existing geologic structures. We suggest that the limits of both the Sawtooth and Lost River faults extend north beyond their mapped extent, are influenced by the relic trans-Challis fault system, and that seismicity within this region will likely continue for the coming years. Ongoing seismic monitoring efforts will lead to an improved understanding of ground shaking potential and active fault characteristics. 
    more » « less
  2. Abstract The Calaveras Fault (CF) branches from the San Andreas Fault (SAF) near San Benito, extending sub‐parallel to the SAF for about 50 km with only 2–6 km separation and diverging northeastward. Both the SAF and CF are partially coupled, exhibit spatially variable aseismic creep and have hosted moderate to large earthquakes in recent decades. Understanding how slip partitions among the main fault strands of the SAF system and establishing their degree of coupling is crucial for seismic hazard evaluation. We perform a timeseries analysis using more than 5 years of Sentinel‐1 data covering the Bay Area (May 2015–October 2020), specifically targeting the spatiotemporal variations of creep rates around the SAF‐CF junction. We derive the surface creep rates from cross‐fault InSAR timeseries differences along the SAF and CF including adjacent Sargent and Quien Sabe Faults. We show that the variable creep rates (0–20 mm/yr) at the SAF‐CF junction are to first order controlled by the angle between the fault strike and the background stress orientation. We further examine the spatiotemporal variation of creep rates along the SAF and CF and find a multi‐annual coupling increase during 2016–2018 the subparallel sections of both faults, with the CF coupling change lagging behind the SAF by 3–6 months. Similar temporal variations are also observed in bothb‐values inferred from declustered seismicity and aseismic slip rates inferred from characteristic repeating earthquakes. 
    more » « less
  3. Abstract Oceanic transform faults connect spreading centers and are imprinted with previous tectonic events. However, their tectonic interactions are not well understood due to limited observations. The Discovery transform fault system at 4°S, East Pacific Rise (EPR), represents a young transform system, offering a unique opportunity to study the interplay between faulting and other tectonic events at an early phases of an oceanic transform system. Discovery regularly hostsM5–6 characteristic earthquakes, and the seafloor north of Discovery includes a 35 km‐long rift zone that records a complex history of rifting, faulting and volcanism, suggesting that the transform faults likely interact with regional tectonic activity. We apply a machine‐learning enabled workflow to locate 21,391 earthquakes recorded during a 1‐year ocean bottom seismometer experiment in 2008. Our results indicate that seismicity on the western Discovery fault is separated into seven patches with distinct aseismic and seismic slip modes. Additionally, we observe a patch of off‐fault seismicity near where seafloor abyssal hills intersect the rift zone. This seismicity may have been caused by varying opening rates as spreading rate decreases from north to south in the rift zone. Our findings suggest that the Discovery system is still evolving, and that system equilibrium has not been reached between rifting and faulting. These results reflect the complex yet rarely observed interactions between fault slip, plate rotation, and rifting which are likely ubiquitous at oceanic transform systems. 
    more » « less
  4. Abstract Mature strike‐slip faults are usually surrounded by a narrow zone of damaged rocks characterized by low seismic wave velocities. Observations of earthquakes along such faults indicate that seismicity is highly concentrated within this fault damage zone. However, the long‐term influence of the fault damage zone on complete earthquake cycles, that is, years to centuries, is not well understood. We simulate aseismic slip and dynamic earthquake rupture on a vertical strike‐slip fault surrounded by a fault damage zone for a thousand‐year timescale using fault zone material properties and geometries motivated by observations along major strike‐slip faults. The fault damage zone is approximated asan elastic layer with lower shear wave velocity than the surrounding rock. We find that dynamic wave reflections, whose characteristics are strongly dependent on the width and the rigidity contrast of the fault damage zone, have a prominent effect on the stressing history of the fault. The presence of elastic damage can partially explain the variability in the earthquake sizes and hypocenter locations along a single fault, which vary with fault damage zone depth, width and rigidity contrast from the host rock. The depth extent of the fault damage zone has a pronounced effect on the earthquake hypocenter locations, and shallower fault damage zones favor shallower hypocenters with a bimodal distribution of seismicity along depth. Our findings also suggest significant effects on the hypocenter distribution when the fault damage zone penetrates to the nucleation sites of earthquakes, likely being influenced by both lithological (material) and rheological (frictional) boundaries. 
    more » « less
  5. Abstract Strain partitioning in oblique convergent margins results in margin‐parallel shear in the overriding plate. Margin‐parallel shear is often accommodated by margin‐parallel strike‐slip faults proximal to active volcanic arcs. Along the Nicaraguan segment of the Central American Forearc (CAFA) in the Cocos‐Caribbean plate convergent margin, there are no well‐expressed right‐lateral faults that accommodate CA‐CAFA relative motion. Instead, historical earthquakes and mapped fault orientations indicate that the ∼12 mm/yr of dextral motion is accommodated on arc‐normal, left‐lateral faults (i.e., bookshelf faults). We investigate three upper‐plate earthquakes; the 10 April 2014 (Mw6.1), 15 September 2016 (Mw5.7), and 28 September 2016 (Mw5.5), using Global Position System co‐seismic displacements and relocated earthquake aftershocks. Our analyses of the three earthquakes indicate that the 10 April 2014 earthquake ruptured an unmapped margin‐parallel right‐lateral fault in Lago Xolotlán (Managua) and the September 2016 earthquakes ruptured arc‐normal, left‐lateral and oblique‐slip faults. These earthquakes represent a triggered sequence whereby the 10 April 2014 earthquake promoted failure of the faults that ruptured in September 2016 by imparting a static Coulomb stress change (ΔCFS) of 0.02–0.07 MPa. Likewise, the 15 September 2016, earthquake additionally promoted failure (ΔCFS of 0.08–0.1 MPa) on sub‐parallel faults that ruptured in two subsequent earthquakes. We also present an instance of magma‐tectonic interaction whereby the 10 April 2014 earthquake dilated (10s of μStrain) the shallow magmatic system of Momotombo Volcano, which led to magma injection, ascent, and eruption on 1 December 2015, after ∼100 years of quiescence. 
    more » « less