skip to main content

Title: Diverse Slip Behavior of the Banyak Islands Subsegment of the Sunda Megathrust in Sumatra, Indonesia

Nearly every type of slip behavior of which megathrusts are known to be capable of has been observed within a narrow ~75‐km‐wide strip of the Sunda megathrust—the Banyak Islands subsegment. The diverse list of recorded slip events starts with the great 1861MW ~ 8.5 earthquake and includes the 1907MW ~ 8.2 shallow tsunami earthquake, the 1966–1981 slow slip event, the great 2005MW8.6 Nias‐Simeulue earthquake and associated afterslip, many small (M < 6) to moderate (6 M < 7) earthquakes, and the large 6 April 2010MW7.8 Banyak Islands earthquake. In this paper we map out the spatial and temporal relationships between these slip events to search for a pattern and possible controls on slip behavior. We use GPS and coral geodetic data to derive a coseismic slip distribution for the relatively understudied 2010 event and find that this event fits like a puzzle piece into a low slip patch left by the 2005 event. We compare the recent slip sequence with the historical sequence associated with the 1861 event and find that the modern slip pattern could be a repeat of the historical slip pattern, possibly suggesting the importance of stationary fault properties in controlling slip. The diverse set of slip events mostly fit our expectations within depth‐dependent rupture zones; however, the overlaps between events of different slip types may complicate this conceptual model. Along strike, we suggest that a subducting fracture zone could act to diversify the slip behavior, and we explore possible mechanisms for this slip control.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    An important manifestation of the background or coseismic deformation surrounding megathrust earthquakes is the offshore microseismicity, which are difficult to be directly detected by land‐based seismic instruments. Here, we improve the capability of detecting offshore events by combining two popular techniques: backprojection (BP) imaging and matched‐filter (MF) detection. The BP method is effective in retrieving offshore seismicity (Mw> 4.5) buried in the coda wave of large earthquakes. The capability of the MF method depends on the availability of the template pool; therefore, the BP‐inferred events can be used as additional templates to expand the MF detections (BP‐MF). We performed the BP‐MF approach in the period within 600 days after the 2011Mw9.0 Tohoku earthquake. We find overall 44.2% more offshore events than those listed in the Japan Meteorological Agency (JMA) catalog. In the near‐trench area, we detect 213% more events. Among the newly detectedMw> 4events, the BP template contributions are twice more than those matched by the JMA templates. Based on the spatial consistency between aftershock‐depleted zones and large coseismic slip, we identify a possible large coseismic slip zone in the near‐trench region offshore Fukushima. Largebvalues (1.2) are found close to large aftershocks, possibly indicating localized pockets of small differential stresses. At several locations close to the trench,pvalues (0.93–1.11) are higher than those in the inland area (0.64–0.85). This may be due to the larger coseismic slip and hence larger stress drop of the outer‐rise normal‐faulting events compared to the deeper thrust‐faulting events.

    more » « less
  2. Abstract

    We conduct a detailed study of the foreshock sequence preceding the 2010Mw6.7 Yushu, Qinghai earthquake in the Tibetan plateau by examining continuous waveforms recorded at a seismic station near the mainshock rupture zone. By using a deep learning phase picker—EQTransformer and a matched‐filter technique, we identify 120 foreshocks with magnitude ranging from −0.7 to 1.6, starting with aMw4.6 foreshock approximately 2 hr before theMw6.7 Yushu mainshock. Our analyses show that the foreshock sequence follows a typical Omori's law decay with ap‐value of 0.73 and the Gutenberg‐Richer frequency‐magnitudeb‐value of 0.66. We do not find any evidence of accelerating events leading up to the Yushu mainshock. Hence, they could be considered as aftershocks of theMw4.6 earthquake. We further invert for the focal mechanisms and rupture directions for both the largest foreshock and the mainshock. TheMw4.6 foreshock likely occurred on a NE‐SW trending fault conjugating to the NW‐SE trending fault of the mainshock. Coulomb stress analysis suggests theMw4.6 foreshock induces negative stress on the mainshock source area. These observations do not support either the pre‐slip or the cascade triggering model for foreshock generation. The occurrence of the foreshock, mainshock and large aftershocks appear to be modulated by the Earth's tidal forces, likely reflecting the role of high pore‐fluid pressures. Our observations, together with other recent studies, suggest that extensional step‐overs and conjugate faults along major strike‐slip faults play an important role in generating short‐term foreshock sequences.

    more » « less
  3. Abstract

    We quantify sliding stability and rupture styles for a horizontal interface between an elastic layer and stiffer elastic half‐space with a free surface on top and rate‐and‐state friction on the interface. This geometry includes shallowly dipping subduction zones, landslides, and ice streams. Specific motivation comes from quasiperiodic slow slip events on the Whillans Ice Plain in West Antarctica. We quantify the influence of layer thickness on sliding stability, specifically whether a steadily loaded system produces steady sliding or stick‐slip sequences. We do this using both linear stability analysis and nonlinear earthquake sequence simulations. We restrict our attention to the 2‐D antiplane shear problem but anticipate that our findings generalize to more complex 2‐D in‐plane and 3‐D problems. Steady sliding with velocity‐weakening rate‐and‐state friction is linearly unstable to Fourier mode perturbations having wavelengths greater than a critical wavelength (λc). We quantify the dependence ofλcon the rate‐and‐state friction parameters, elastic properties, loading, and the layer thickness (H). Confirming previous studies, we find thatλc ∝ H1/2for smallHand is independent ofHfor largeH. The linear stability analysis provides insight into nonlinear earthquake sequence dynamics of a nominally velocity‐strengthening interface containing a velocity‐weakening region of widthW. Sequence simulations reveal a transition from steady sliding at smallWto stick‐slip events whenWexceeds a critical width (Wcr), withWcr ∝ H1/2for smallH. Overall, this study demonstrates that the reduced stiffness of thin layers promotes instability, with implications for sliding dynamics in thin layer geometries.

    more » « less
  4. Abstract

    The 2010MW7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events (M≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus onPnand its depth phasespPnandsPn, which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stackPnwaveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of thePn, directP, and directSphases recorded for event pairs. Focal depths of 93 out of 122M≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.

    more » « less
  5. Abstract

    Foreshocks can provide valuable information about possible nucleation process of a mainshock. However, their physical mechanisms are still under debate. In this study, we present a comprehensive analysis of the earthquake sequence preceding the 2010 Mw7.2 El Mayor‐Cucapah mainshock, including waveform detection of missing smaller events, relative relocation, and source parameter analysis. Based on a template matching method, we find a tenfold increase in the number of earthquakes than reported in the Southern California Seismic Network catalog. The entire sequence exhibits nearly continuous episodes of foreshocks that can be loosely separated into two active clusters. Relocated foreshocks show several seismicity streaks at depth, with a consistently active cluster at depths between 14 and 16 km where the mainshock was nucleated. Stress drop measurements from a spectral ratio approach based on empirical Green's functions show a range between 3.8 and 41.7 MPa with a median of 13.0 MPa and no clear temporal variations. The relocation results, together with the source patches estimated from earthquake corner frequencies, revealed a migration front toward the mainshock hypocenter within last 8 hr and a chain of active burst immediately 6 min prior to the mainshock. Our results support combined effects of aseismic slip and cascading failure on the evolution of foreshocks.

    more » « less