The 2019–2020 Southwest Puerto Rico Earthquake Sequence: Seismicity and Faulting

Author(s):  
Blaž Vičič ◽  
Seyyedmaalek Momeni ◽  
Alessandra Borghi ◽  
Anthony Lomax ◽  
Abdelkrim Aoudia

Abstract The 2019–2020 Southwest Puerto Rico earthquake sequence ruptured multiple faults with several moderate magnitude earthquakes. Here, we investigate the seismotectonics of this fault system using high-precision hypocenter relocation and inversion of the near-field strong motions of the five largest events in the sequence (5.6≤Mw≤6.4) for kinematic rupture models. The Mw 6.4 mainshock occurred on a northeast-striking, southeast-dipping normal fault. The rupture nucleated offshore ∼15 km southeast of Indios at the depth of 8.6 km and extended southwest–northeast and up-dip with an average speed of 1.55 km/s, reaching the seafloor and shoreline after about 8 s. The 6 January 2020 (10:32:23) Mw 5.7 and the 7 January 2020 (11:18:46) Mw 5.8 events occurred on two east–southeast-striking, near-vertical, left-lateral strike-slip faults. However, the 7 January 2020 (08:34:05) Mw 5.6 normal-faulting aftershock, which occurred only 10 min after the Mw 6.4 normal-faulting mainshock, ruptured on a fault with almost the same strike as the mainshock but situated ∼8 km farther east, forming a set of parallel faults in the fault system. On 11 January 2020, an Mw 6.0 earthquake occurred on a north–northeast-striking, westing-dipping fault, orthogonal to the faults hosting the strike-slip earthquakes. We apply template matching for the detection of missed, small-magnitude earthquakes to study the spatial evolution of the main part of the sequence. Using the template-matching results along with Global Positioning System analysis, we image the temporal evolution of a foreshock sequence (Caja swarm). We propose that the swarm and the main sequence were a response to a tectonic transient that most affected the whole Puerto Rico Island.

Author(s):  
Copeland W. Cromwell ◽  
Kevin P. Furlong ◽  
Eric A. Bergman ◽  
Harley M. Benz ◽  
Will L. Yeck ◽  
...  

Abstract We present a new catalog of calibrated earthquake relocations from the 2019–2020 Puerto Rico earthquake sequence related to the 7 January 2020 Mw 6.4 earthquake that occurred offshore of southwest Puerto Rico at a depth of 15.9 km. Utilizing these relocated earthquakes and associated moment tensor solutions, we can delineate several distinct fault systems that were activated during the sequence and show that the Mw 6.4 mainshock may have resulted from positive changes in Coulomb stress from earlier events. Seismicity and mechanisms define (1) a west–southwest (∼260°) zone of seismicity comprised of largely sinistral strike-slip and oblique-slip earthquakes that mostly occurs later in the sequence and to the west of the mainshock, (2) an area of extensional faulting that includes the mainshock and occurs largely within the mainshock’s rupture area, and (3) an north–northeast (∼30°)-striking zone of seismicity, consisting primarily of dextral strike-slip events that occurs before and following the mainshock and generally above (shallower than) the normal-faulting events. These linear features intersect within the Mw 6.4 mainshock’s fault plane in southwest Puerto Rico. In addition, we show that earthquake relocations for M 4+ normal-faulting events, when traced along their fault planes, daylight along east–west-trending bathymetric features offshore of southwest Puerto Rico. Correlation of these normal-faulting events with bathymetric features suggests an active fault system that may be a contributor to previously uncharacterized seismic hazards in southwest Puerto Rico.


1995 ◽  
Vol 85 (3) ◽  
pp. 705-715
Author(s):  
Mark Andrew Tinker ◽  
Susan L. Beck

Abstract Regional distance surface waves are used to study the source parameters for moderate-size aftershocks of the 25 April 1992 Petrolia earthquake sequence. The Cascadia subduction zone had been relatively seismically inactive until the onset of the mainshock (Ms = 7.1). This underthrusting event establishes that the southern end of the North America-Gorda plate boundary is seismogenic. It was followed by two separate and distinct large aftershocks (Ms = 6.6 for both) occurring at 07:41 and 11:41 on 26 April, as well as thousands of other small aftershocks. Many of the aftershocks following the second large aftershock had magnitudes in the range of 4.0 to 5.5. Using intermediate-period surface-wave spectra, we estimate focal mechanisms and depths for one foreshock and six of the larger aftershocks (Md = 4.0 to 5.5). These seven events can be separated into two groups based on temporal, spatial, and principal stress orientation characteristics. Within two days of the mainshock, four aftershocks (Md = 4 to 5) occurred within 4 hr of each other that were located offshore and along the Mendocino fault. These four aftershocks comprise one group. They are shallow, thrust events with northeast-trending P axes. We interpret these aftershocks to represent internal compression within the North American accretionary prism as a result of Gorda plate subduction. The other three events compose the second group. The shallow, strike-slip mechanism determined for the 8 March foreshock (Md = 5.3) may reflect the right-lateral strike-slip motion associated with the interaction between the northern terminus of the San Andreas fault system and the eastern terminus of the Mendocino fault. The 10 May aftershock (Md = 4.1), located on the coast and north of the Mendocino triple junction, has a thrust fault focal mechanism. This event is shallow and probably occurred within the accretionary wedge on an imbricate thrust. A normal fault focal mechanism is obtained for the 5 June aftershock (Md = 4.8), located offshore and just north of the Mendocino fault. This event exhibits a large component of normal motion, representing internal failure within a rebounding accretionary wedge. These two aftershocks and the foreshock have dissimilar locations in space and time, but they do share a north-northwest oriented P axis.


2019 ◽  
Vol 220 (1) ◽  
pp. 585-597 ◽  
Author(s):  
Maria C Araya ◽  
Juliet Biggs

SUMMARY Tectonic slivers form in the overriding plate in regions of oblique subduction. The inner boundaries of the sliver are often poorly defined and can consist of well-defined faults, rotating blocks or diffuse fault systems, which pass through or near the volcanic arc. The Guanacaste Volcanic Arc Sliver (GVAS) as defined by Montero et al., is a segment of the Central American Forearc Sliver, whose inner boundary is the ∼87-km-long Haciendas-Chiripa fault system (HCFS), which is located ∼10 km behind the volcanic arc and consists of strike slip faults and pull apart steps. We characterize the current ground motion on this boundary by combining earthquake locations and focal mechanisms of the 2016 Bijagua earthquake sequence, with the surface ground deformation obtained from Interferometric Synthetic Aperture Radar (InSAR) images from the ALOS-2 satellite. The coseismic stack of interferograms show ∼6 cm of displacement towards the line of sight of the satellite between the Caño Negro fault and the Upala fault, indicating uplift or SE horizontal surface displacement. The largest recorded earthquake of the sequence was Mw 5.4, and the observed deformation is one of the smallest earthquakes yet detected by InSAR in the Central American region. Forward and inverse models show the surface deformation can be partially explained by slip on a single fault, but it can be better explained by slip along two faults linked at depth. The best-fitting model consists of 0.33 m of right lateral slip on the Caño Negro fault and 0.35 m of reverse slip on the Upala fault, forming a positive flower structure. As no reverse seismicity was recorded, we infer the slip on the Upala fault occurred aseismically. Observations of the Bijagua earthquake sequence suggests the forearc sliver boundary is a complex and diffuse fault system. There are localized zones of transpression and transtension and areas where there is no surface expression suggesting the fault system is not yet mature. Although aseismic slip is common on subduction interfaces and mature strike-slip faults, this is the first study to document aseismic slip on a continental tectonic sliver boundary fault.


Active and recent faulting along the main north—south road in Tibet is dominated by normal faulting occurring on northerly-trending planes and by strike-slip faulting, both of which reflect an east-west extension of the plateau. Normal faulting is prevalent in the southern half of the plateau, but we saw no evidence for any major graben in the northern half. Strike-slip faulting on roughly easterly-trending structures is m ore prevalent in the northern half, but conjugate faulting, with right-lateral slip on northwesterly-trending planes and left-lateral slip on northeasterly-trending planes, is common in the southern half. In two areas, we also observed components of thrust faulting, apparently in association with young strikeslip faulting. Our most important results are bounds on the rates of slip on the two main strands of the Kunlun strike-slip fault system, which trends east-w est through the Kunlun range. Ground moraine containing boulders of pyroxenite is separated by 30 km from the nearest outcrop of such rock, implying that amount of displacement in the last 1.5 to 3 M a. Therefore the average rate of slip during the Quaternary period has been between 10 and 20 mm/a , with a likely value of 13 mm/a . Abundant fresh tension cracks and mole tracks imply continued slip on the main strand, the Xidatan -Tuosuohu-Maqu fault, and the likely occurrence of a major earthquake in the last few hundred years. Consistent offsets of gullies and dry stream channels of about 10 m may reflect slip of that amount during such an earthquake, and possible multiple offsets at one site suggest that slip may occur by large displacements of 10 m during infrequent great earthquakes. Along the other strand, the Kunlun Pass fault, offsets of roughly 50 to 150 m of, apparently, post-glacial valleys and of one glacier and its terminal moraine suggest a Holocene rate of slip between 5 and 20 mm/a , and most likely about 10 mm/a , on this fault. These rapid rates of displacement imply that Tibet is being extruded rapidly eastward, at a rate com parable to the rate at which India is penetrating into Eurasia, and therefore that, at present, a substantial fraction of this penetration is being absorbed by the eastward extrusion of Tibet.


2018 ◽  
Vol 40 (3) ◽  
pp. 1125 ◽  
Author(s):  
A. Kiratzi ◽  
C. Benetatos ◽  
Z. Roumelioti

Nearly 2,000 earthquake focal mechanisms in the Aegean Sea and the surroundings for the period 1912- 2006, for 1.5 <M<7.5, and depths from 0 to 170 km, indicate a uniform distribution and smooth variation in orientation over wide regions, even for the very small magnitude earthquakes. ~ 60% of the focal mechanisms show normal faulting, that mainly strikes ~E-W. However, a zone ofN-S normal faulting runs the backbone of Albanides-Hellenides. Low-angle thrust and reverse faulting is confined in western Greece (Adria-Eurasia convergence) and along the Hellenic trench (Africa-Eurasia). In the central Aegean Sea the effect of the propagating tip of the North Anatolian Fault into the Aegean Sea is pronounced and strike-slip motions are widely distributed. Shearing does not cross central Greece. Strike-slip motions reappear in the Cephalonia-Lefkada Transform Fault zone and in western Péloponnèse, which shows very complex tectonics, with different types of faulting being oriented favourably and operating under the present stress-field. Moreover, in western Péloponnèse the sense of the observed shearing is not yet clear, whether it is dextral or sinistral, and this lack of data has significant implications for the orientation of the earthquake slip vectors compared to the GPS obtained velocity vectors.


Author(s):  
Fred F. Pollitz ◽  
William C. Hammond ◽  
Charles W. Wicks

Abstract The 2020 M 6.5 Stanley, Idaho, earthquake produced rupture in the north of the active Sawtooth fault in the northern basin and range at depth, without any observable surface rupture. Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data yield several millimeters of static offsets out to ∼100  km from the rupture and up to ∼0.1  m of near-field crustal deformation. We combine the GPS and InSAR data with long-period regional seismic waveforms to derive models of kinematic slip and afterslip. We find that the coseismic rupture is complex, likely involving up to 2 m combined left-lateral strike slip and normal slip on a previously unidentified ∼south-southeast-striking fault. This slip is predominantly left-lateral strike slip, different from the dominant east-northeast–west-northwest normal faulting of the region. At least one ∼northeast-trending fault, likely associated with the Trans-Challis fault system, is inferred to have accommodated a few decimeters of right-lateral afterslip, consistent with vigorous aftershock activity at depth along northeast-trending lineations.


2021 ◽  
Author(s):  
◽  
Vasiliki Mouslopoulou

<p>The North Island of New Zealand sits astride the Hikurangi margin along which the oceanic Pacific Plate is being obliquely subducted beneath the continental Australian Plate. The North Island Fault System1 (NIFS), in the North Island of New Zealand, is the principal active strike-slip fault system in the overriding Australian Plate accommodating up to 30% of the margin parallel plate motion. This study focuses on the northern termination of the NIFS, near its intersection with the active Taupo Rift, and comprises three complementary components of research: 1) the investigation of the late Quaternary (c. 30 kyr) geometries and kinematics of the northern NIFS as derived from displaced geomorphic landforms and outcrop geology, 2) examination of the spatial and temporal distribution of  paleoearthquakes in the NIFS over the last 18 kyr, as derived by fault-trenching and displaced landforms, and consideration of how these distributions may have produced the documented late Quaternary (c. 30 kyr) kinematics of the northern NIFS, and 3) Investigation of the temporal stability of the late Quaternary (c. 30 kyr) geometries and kinematics throughout the Quaternary (1-2 Ma), derived from gravity, seismic-reflection, drillhole, topographic and outcrop data. The late Quaternary (c. 30 kyr) kinematics of the northern NIFS transition northward along strike, from strike-slip to oblique-normal faulting, adjacent to the rift. With increasing proximity to the Taupo Rift the slip vector pitch on each of the faults in the NIFS steepens gradually by up to 60 degrees, while the mean fault-dip decreases from 90 degrees to 60 degrees W. Adjustments in the kinematics of the NIFS reflect the gradual accommodation of the NW-SE extension that is distributed outside the main physiographic boundary of the Taupo Rift. Sub-parallelism of slip vectors in the NIFS with the line of intersection between the two synchronous fault systems reduces potential space problems and facilitates the development of a kinematically coherent fault intersection, which allows the strike-slip component of slip to be transferred into the rift. Transfer of displacement from the NIFS into the rift accounts for a significant amount of the northeastward increase of extension along the rift. Steepening of the pitch of slip vectors towards the northern termination of the NIFS allows the kinematics and geometry of faulting to change efficiently, from strike-lip to normal faulting, providing an alternative mechanism to vertical axis rotations for terminating large strike-lip faults. Analyses of kinematic constraints from worldwide examples of synchronous strike-lip and normal faults that intersect to form two or three plate configurations, within either oceanic or continental crust, suggest that displacement is often transferred between the two fault systems in a similar manner to that documented at the NIFS - Taupo Rift fault intersection. The late Quaternary (c. 30 kyr) change in the kinematics of the NIFS along strike, from dominantly strike-slip to oblique-normal faulting, arises due to a combination of rupture arrest during individual earthquakes and variations in the orientation of the coseismic slip vectors. At least 80 % of all surface rupturing earthquakes appear to have terminated within the kinematic transition zone from strike-slip to oblique-normal slip. Fault segmentation reduces the magnitudes of large surface rupturing earthquakes in the northern NIFS from 7.4-7.6 to c. 7.0. Interdependence of throw rates between the NIFS and Taupo Rift suggests that the intersection of the two fault systems has functioned coherently for much of the last 0.6-1.5 Myr. Oblique-normal slip faults in the NIFS and the Edgecumbe Fault in the rift accommodated higher throw rates since 300 kyr than during the last 0.6-1.5 Myr. Acceleration of these throw rates may have occurred in response to eastward migration of rifting, increasing both the rates of faulting and the pitch of slip vectors. The late Quaternary (e.g. 30 kyr) kinematics, and perhaps also the stability, of the intersection zone has been geologically short lived and applied for the last c. 300 kyr.</p>


2021 ◽  
Author(s):  
Michael Rudolf ◽  
Joscha Podlesny ◽  
Esther Heckenbach ◽  
Matthias Rosenau ◽  
Anne Glerum ◽  
...  

&lt;p&gt;The release of elastic energy along an active fault is accommodated by a wide range of slip modes. It ranges from long-term slow slip events (SSEs) and creep to short-term tremors and earthquakes. They vary not only in their characteristic duration but also in their magnitude, spatial exten&lt;span&gt;&lt;span&gt;t&lt;/span&gt;&lt;/span&gt; and slip velocities. The exact relationship is unclear, as in some regions many slip modes occur simultaneously (e.g. Tohoku-Oki) and in others certain slip modes are completely absent (e.g. Cascadia).&lt;/p&gt;&lt;p&gt;One of the driving factors in the generation of this large variety of slip modes is the interplay of fault heterogeneity and geometrical complexity of the fault system. We test various settings in terms of fault heterogeneity and geometrical complexity with a scaled physical model. The experimental results are then validated and benchmarked through multi-scale numerical simulations. We describe &lt;span&gt;&lt;span&gt;the&lt;/span&gt;&lt;/span&gt; system using &lt;span&gt;&lt;span&gt;a&lt;/span&gt;&lt;/span&gt; rate-and-state frictional framework and introduce on-fault heterogeneity with variable frictional properties. All properties are the same for analogue and numerical simulation as far as they can be determined or realized experimentally (a-b, v&lt;sub&gt;load&lt;/sub&gt;, S&lt;sub&gt;hmax&lt;/sub&gt;, S&lt;sub&gt;hmin&lt;/sub&gt;, etc...). As analogue material we use segmented, decimetre sized neoprene foam blocks in multiple configurations (e.g. biaxial shear at forces &lt;1 kN) to simulate the elastic upper crust. The contact surfaces are spray-painted with acrylic paint to generate velocity weakening characteristics in between the blocks which is similar to the frictional behaviour of natural faults. We add heterogeneity to the fault surface by varying the fault area that is velocity weakening using grease. Geometrical complexity is implemented using conjugated or parallel sets of additional faults with the same characteristics.&lt;/p&gt;&lt;p&gt;We are able to reliably generate frequent stick-slip events of variable size and recurrence intervals. The slip characteristics, such as slip distribution, are in good agreement with analytical solutions of fault slip in elastic media. In a geometrically simple strike-slip model the recurrence behaviour and magnitude follows straightforward scaling relations in accordance with existing studies. If geometrical complexity is added to the model we observe clustering and variable recurrence that differ from the simpler geometry. Additionally, we are going to give an outlook on the interaction behaviour of multiple faults in dependence of their geometric configuration and the generation of power-law type magnitude scaling relations.&lt;/p&gt;


Author(s):  
Marjolein Blasweiler ◽  
Matthew W. Herman ◽  
Fenna Houtsma ◽  
Rob Govers

Abstract An historically unprecedented seismic moment was released by crustal events of the 2019–2020 earthquake sequence near southwest Puerto Rico. The sequence involved at least two, and perhaps three interacting fault systems. The largest Mw 6.4 event was likely triggered by left lateral strike-slip events along the eastern extension of the North Boquerón Bay-Punta Montalva fault zone. The mainshock occurred in a normal fault zone that extends into a region where previous studies documented extensional deformation, beyond the Ponce fault and the Bajo Tasmanian fault. Coulomb stress changes by the mainshock may have triggered further normal-faulting aftershocks, left lateral strike-slip events in the region where these two fault zones interacted, and possibly right lateral strike-slip aftershocks along a third structure extending southward, the Guayanilla fault zone. Extension directions of the seismic sequence are consistently north-northwest–south-southeast-oriented, in agreement with the Global Navigation Satellite Systems-inferred motion direction of eastern Hispaniola relative to western Puerto Rico, and with crustal stress estimates for the overriding plate boundary zone.


Sign in / Sign up

Export Citation Format

Share Document