A window into the complexity of the dynamic rupture of the 2011 Mw 9 Tohoku-Oki earthquake

The 2011 Mw 9 Tohoku-Oki earthquake, recorded by over 1000 near-field stations and multiple large-aperture arrays, is by far the best recorded earthquake in the history of seismology and provides unique opportunities to address fundamental issues in earthquake source dynamics. Here we conduct a high resolution array analysis based on recordings from the USarray and the European network. The mutually consistent results from both arrays reveal rupture complexity with unprecedented resolution, involving phases of diverse rupture speed and intermittent high frequency bursts within slow speed phases, which suggests spatially heterogeneous material properties. The earthquake initially propagates down-dip, with a slow initiation phase followed by sustained propagation at speeds of 3 km/s. The rupture then slows down to 1.5 km/s for 60 seconds. A rich sequence of bursts is generated along the down-dip rim of this slow and roughly circular rupture front. Before the end of the slow phase an extremely fast rupture front detaches at about 5 km/s towards the North. Finally a rupture front propagates towards the south running at about 2.5 km/s for over 100 km. Key features of the rupture process are confirmed by the strong motion data recorded by K-net and KIK-net. The energetic high frequency radiation episodes within a slow rupture phase suggests a patchy image of the brittle-ductile transition zone, composed of discrete brittle asperities within a ductile matrix. The high frequency is generated mainly at the down-dip edge of the principal slip regions constrained by geodesy, suggesting a variation along dip of the mechanical properties of the mega thrust fault or their spatial heterogeneity that affects rise time

The 2013 Okhotsk deep-focus earthquake: Rupture beyond the metastable olivine wedge and thermally controlled rise time near the edge of a slab

The 2013 M8.3 Okhotsk earthquake involves two primary mechanisms of deep-focus earthquake rupture, mineral phase transformation of olivine to spinel and thermal shear instability. Backprojection imaging of broadband seismograms recorded by the North American and European networks indicates bilateral rupture toward NE and SSE. The rupture paths of the NE segment and other regional M7 earthquakes are confined in narrow regions along the slab contours, consistent with the phase transformation mechanism. However, the SSE rupture propagates a long distance across the slab and aftershocks are distributed across a ~60 km wide zone, beyond the plausible thickness of the metastable olivine wedge, favoring thermal shear weakening. While the NE rupture is only visible at high frequencies, the SSE rupture is consistently observed across a broad-frequency range. This frequency-dependent rupture mode can be explained by lateral variations of rise time controlled by thermal thinning of the slab near its northern end

Localized water reverberation phases and its impact on back-projection images

Coherent radiators imaged by backprojections (BP) are commonly interpreted as part of the rupture process. Nevertheless, artifacts introduced by structure related phases are rarely discriminated from the rupture process. In this study, we use a calibration event to discriminate between rupture and structure effects. We reexamine the waveforms and BP images of the 2012 M_w 7.2 Indian Ocean earthquake and a calibration event (M_w 6.2). The P wave codas of both events present similar shape with characteristic period of approximately 10 s, which are backprojected as coherent radiators near the trench. S wave BP does not image energy radiation near the trench. We interpret those coda waves as localized water reverberation phases excited near the trench. We perform a 2-D waveform modeling using realistic bathymetry model and find that the steep near-trench bathymetry traps the acoustic water waves forming localized reverberation phases. These waves can be imaged as coherent near-trench radiators with similar features as that in the observations. We present a set of methodologies to discriminate between the rupture and propagation effects in BP images, which can serve as a criterion of subevent identification

Transpressional Rupture Cascade of the 2016 M_w 7.8 Kaikoura Earthquake, New Zealand

Large earthquakes often do not occur on a simple planar fault but involve rupture of multiple geometrically complex faults. The 2016 M_w 7.8 Kaikoura earthquake, New Zealand, involved the rupture of at least 21 faults, propagating from southwest to northeast for about 180 km. Here we combine space geodesy and seismology techniques to study subsurface fault geometry, slip distribution, and the kinematics of the rupture. Our finite‐fault slip model indicates that the fault motion changes from predominantly right‐lateral slip near the epicenter to transpressional slip in the northeast with a maximum coseismic surface displacement of about 10 m near the intersection between the Kekerengu and Papatea faults. Teleseismic back projection imaging shows that rupture speed was overall slow (1.4 km/s) but faster on individual fault segments (approximately 2 km/s) and that the conjugate, oblique‐reverse, north striking faults released the largest high‐frequency energy. We show that the linking Conway‐Charwell faults aided in propagation of rupture across the step over from the Humps fault zone to the Hope fault. Fault slip cascaded along the Jordan Thrust, Kekerengu, and Needles faults, causing stress perturbations that activated two major conjugate faults, the Hundalee and Papatea faults. Our results shed important light on the study of earthquakes and seismic hazard evaluation in geometrically complex fault systems

Dual megathrust slip behaviors of the 2014 Iquique earthquake sequence

The transition between seismic rupture and aseismic creep is of central interest to better understand the mechanics of subduction processes. A Mw 8.2 earthquake occurred on April 1st, 2014 in the Iquique seismic gap of northern Chile. This event was preceded by a long foreshock sequence including a 2-week-long migration of seismicity initiated by a Mw 6.7 earthquake. Repeating earthquakes were found among the foreshock sequence that migrated towards the mainshock hypocenter, suggesting a large-scale slow-slip event on the megathrust preceding the mainshock. The variations of the recurrence times of the repeating earthquakes highlight the diverse seismic and aseismic slip behaviors on different megathrust segments. The repeaters that were active only before the mainshock recurred more often and were distributed in areas of substantial coseismic slip, while repeaters that occurred both before and after the mainshock were in the area complementary to the mainshock rupture. The spatiotemporal distribution of the repeating earthquakes illustrates the essential role of propagating aseismic slip leading up to the mainshock and illuminates the distribution of postseismic afterslip. Various finite fault models indicate that the largest coseismic slip generally occurred down-dip from the foreshock activity and the mainshock hypocenter. Source imaging by teleseismic back-projection indicates an initial down-dip propagation stage followed by a rupture-expansion stage. In the first stage, the finite fault models show an emergent onset of moment rate at low frequency (0.5 Hz>0.5 Hz). This indicates frequency-dependent manifestations of seismic radiation in the low-stress foreshock region. In the second stage, the rupture expands in rich bursts along the rim of a semi-elliptical region with episodes of re-ruptures, suggesting delayed failure of asperities. The high-frequency rupture remains within an area of local high trench-parallel gravity anomaly (TPGA), suggesting the presence of subducting seamounts that promote high-frequency generation. Our results highlight the complexity of the interactions between large-scale aseismic slow-slip and dynamic ruptures of megathrust earthquakes

Rupture Heterogeneity and Directivity Effects in Back-Projection Analysis

The back projection method is a tremendously powerful technique for investigating the time dependent earthquake source, but its physical interpretation is elusive. We investigate how earthquake rupture heterogeneity and directivity can affect back-projection results (imaged location and beam power) using synthetic earthquake models. Rather than attempting to model the dynamics of any specific real earthquake, we use idealized kinematic rupture models, with constant or varying rupture velocity, peak slip rate, and fault-local strike orientation along unilateral or bilateral rupturing faults, and perform back-projection with the resultant synthetic seismograms. Our experiments show back-projection can track only heterogeneous rupture processes;homogeneous rupture is not resolved in our synthetic experiments. The amplitude of beam power does not necessarily correlate with the amplitude of any specific rupture parameter (e.g., slip rate or rupture velocity) at the back-projected location. Rather, it depends on the spatial heterogeneity around the back-projected rupture front, and is affected by the rupture directivity. A shorter characteristic wavelength of the source heterogeneity or rupture directivity toward the array results in strong beam power in higher frequency. We derive an equation based on Doppler theory to relate the wavelength of heterogeneity with synthetic seismogram frequency. This theoretical relation can explain the frequency- and array-dependent back-projection results not only in our synthetic experiments but also to analyze the 2019 M7.6 bilaterally rupturing New Ireland earthquake. Our study provides a novel perspective to physically interpret back-projection results and to retrieve information about earthquake rupture characteristics

Transpressional Rupture Cascade of the 2016 M_w 7.8 Kaikoura Earthquake, New Zealand

Large earthquakes often do not occur on a simple planar fault but involve rupture of multiple geometrically complex faults. The 2016 M_w 7.8 Kaikoura earthquake, New Zealand, involved the rupture of at least 21 faults, propagating from southwest to northeast for about 180 km. Here we combine space geodesy and seismology techniques to study subsurface fault geometry, slip distribution, and the kinematics of the rupture. Our finite‐fault slip model indicates that the fault motion changes from predominantly right‐lateral slip near the epicenter to transpressional slip in the northeast with a maximum coseismic surface displacement of about 10 m near the intersection between the Kekerengu and Papatea faults. Teleseismic back projection imaging shows that rupture speed was overall slow (1.4 km/s) but faster on individual fault segments (approximately 2 km/s) and that the conjugate, oblique‐reverse, north striking faults released the largest high‐frequency energy. We show that the linking Conway‐Charwell faults aided in propagation of rupture across the step over from the Humps fault zone to the Hope fault. Fault slip cascaded along the Jordan Thrust, Kekerengu, and Needles faults, causing stress perturbations that activated two major conjugate faults, the Hundalee and Papatea faults. Our results shed important light on the study of earthquakes and seismic hazard evaluation in geometrically complex fault systems

The 2012 Sumatra great earthquake sequence

The equatorial Indian Ocean is a well known place of active intraplate deformation defying the conventional view of rigid plates separated by narrow boundaries where deformation is confined. On 11 April 2012, this region was hit in a couple of hours by two of the largest strike-slip earthquakes ever recorded (moment magnitudes Mw=8.6 and 8.2). Broadband seismological observations of the Mw=8.6 mainshock indicate a large centroid depth (∼30 km) and remarkable rupture complexity. Detailed study of the surface-wave directivity and moment rate functions clearly indicates the partition of the rupture into at least two distinct subevents. To account for these observations, we developed a procedure to invert for multiple-point-source parameters. The optimum source model at long period consists of two point sources separated by about 209 km with magnitudes Mw=8.5 and 8.3. To explain the remaining discrepancies between predicted and observed surface waves, we can refine this model by adding directivity along the WNW–ESE axis. However, we do not exclude more complicated models. To analyze the Mw=8.2 aftershock, we removed the perturbation due to large surface-wave arrivals of the Mw=8.6 mainshock by subtracting the corresponding synthetics computed for the two-subevent model. Analysis of the surface-wave amplitudes suggests that the Mw=8.2 aftershock had a large centroid depth between 30 km and 40 km. This major earthquake sequence brings a new perspective to the seismotectonics of the equatorial Indian Ocean and reveals active deep lithospheric deformation

Early and persistent supershear rupture of the 2018 Mw 7.5 Palu earthquake

The speed at which an earthquake rupture propagates affects its energy balance and ground shaking impact. Dynamic models of supershear earthquakes, which are faster than the speed of shear waves, often start at subshear speed and later run faster than Eshelby’s speed. Here we present robust evidence of an early and persistent supershear rupture at the sub-Eshelby speed of the 2018 magnitude 7.5 Palu, Indonesia, earthquake. Slowness-enhanced back-projection of teleseismic data provides a sharp image of the rupture process, along a path consistent with the surface rupture trace inferred by subpixel correlation of synthetic-aperture radar and satellite optical images. The rupture propagated at a sustained velocity of 4.1 km s^(–1) from its initiation to its end, despite large fault bends. The persistent supershear speed is further validated by seismological evidence of far-field Rayleigh Mach waves. The unusual features of this earthquake probe the connections between the rupture dynamics and fault structure. An early supershear transition could be promoted by fault roughness near the hypocentre. Steady rupture propagation at a speed unexpected in homogeneous media could result from the presence of a low-velocity damaged fault zone