Lingsen Meng

The 2011 Magnitude 9.0 Tohoku-Oki Earthquake: Mosaicking the Megathrust from Seconds to Centuries

Detailed geophysical measurements reveal features of the 2011 Tohoku-Oki megathrust earthquake. Geophysical observations from the 2011 moment magnitude (Mw) 9.0 Tohoku-Oki, Japan earthquake allow exploration of a rare large event along a subduction megathrust. Models for this event indicate that the distribution of coseismic fault slip exceeded 50 meters in places. Sources of high-frequency seismic waves delineate the edges of the deepest portions of coseismic slip and do not simply correlate with the locations of peak slip. Relative to the Mw 8.8 2010 Maule, Chile earthquake, the Tohoku-Oki earthquake was deficient in high-frequency seismic radiation—a difference that we attribute to its relatively shallow depth. Estimates of total fault slip and surface secular strain accumulation on millennial time scales suggest the need to consider the potential for a future large earthquake just south of this event.

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.

Earthquake in a Maze: Compressional Rupture Branching During the 2012 Mw 8.6 Sumatra Earthquake

Earthquake in a Maze The 11 April 2012 magnitude 8.6 earthquake offshore of Sumatra was the largest measured earthquake along a strike-slip boundary that modern seismological instruments have ever recorded. Despite its size and proximity to a large population, there was no subsequent tsunami and there were no reported fatalities. Meng et al. (p. 724, published online 19 July) used teleseismic data from seismological networks in Japan and Europe to image the source of high-frequency radiation generated by the earthquake to understand the mechanics of this unique event. The resultant back projections showed that the earthquake slowly ruptured along a complex series of faults. The deeper-than-usual rupture path and large stress drop are both features that may not be unique to this earthquake, suggesting that regions in a similar tectonic environment may have the potential for more complex—or larger—intraplate earthquakes than might have been expected. The mechanics of the largest strike-slip earthquake ever recorded give clues about how intraplate earthquakes rupture. Seismological observations of the 2012 moment magnitude 8.6 Sumatra earthquake reveal unprecedented complexity of dynamic rupture. The surprisingly large magnitude results from the combination of deep extent, high stress drop, and rupture of multiple faults. Back-projection source imaging indicates that the rupture occurred on distinct planes in an orthogonal conjugate fault system, with relatively slow rupture speed. The east-southeast–west-northwest ruptures add a new dimension to the seismotectonics of the Wharton Basin, which was previously thought to be controlled by north-south strike-slip faulting. The rupture turned twice into the compressive quadrant, against the preferred branching direction predicted by dynamic Coulomb stress calculations. Orthogonal faulting and compressional branching indicate that rupture was controlled by a pressure-insensitive strength of the deep oceanic lithosphere.

High-resolution backprojection at regional distance: Application to the Haiti M7.0 earthquake and comparisons with finite source studies,

A catastrophic M_w7 earthquake ruptured on 12 January 2010 on a complex fault system near Port-au-Prince, Haiti. Offshore rupture is suggested by aftershock locations and marine geophysics studies, but its extent remains difficult to define using geodetic and teleseismic observations. Here we perform the multitaper multiple signal classification (MUSIC) analysis, a high-resolution array technique, at regional distance with recordings from the Venezuela National Seismic Network to resolve high-frequency (about 0.4 Hz) aspects of the earthquake process. Our results indicate westward rupture with two subevents, roughly 35 km apart. In comparison, a lower-frequency finite source inversion with fault geometry based on new geologic and aftershock data shows two slip patches with centroids 21 km apart. Apparent source time functions from USArray further constrain the intersubevent time delay, implying a rupture speed of 3.3 km/s. The tips of the slip zones coincide with subevents imaged by backprojections. The different subevent locations found by backprojection and source inversion suggest spatial complementarity between high- and low-frequency source radiation consistent with high-frequency radiation originating from rupture arrest phases at the edges of main slip areas. The centroid moment tensor (CMT) solution and a geodetic-only inversion have similar moment, indicating most of the moment released is captured by geodetic observations and no additional rupture is required beyond where it is imaged in our preferred model. Our results demonstrate the contribution of backprojections of regional seismic array data for earthquakes down to M ≈ 7, especially when incomplete coverage of seismic and geodetic data implies large uncertainties in source inversions.

A dynamic model of the frequency-dependent rupture process of the 2011 Tohoku-Oki earthquake

We present a 2D dynamic rupture model that provides a physical interpretation of the key features of the 2011 Tohoku-Oki earthquake rupture. This minimalistic model assumes linear slip-weakening friction, the presence of deep asperities and depth-dependent initial stresses. It reproduces the first-order observations of the along-dip rupture process during its initial 100 s, such as large static slip and low-frequency radiation up-dip from the hypocenter, and slow rupture punctuated by high-frequency radiation in deeper regions. We also derive quantitative constraints on the ratio of shallow versus deep radiation from teleseismic back-projection source imaging. This ratio is explained in our model by the rupture of deep asperities surrounded by low stress drop regions, and by the decrease of initial stresses towards the trench.

Mitigating artifacts in back-projection source imaging with implications for frequency-dependent properties of the Tohoku-Oki earthquake

Comparing teleseismic array back-projection source images of the 2011 Tohoku-Oki earthquake with results from static and kinematic finite source inversions has revealed little overlap between the regions of high- and low-frequency slip. Motivated by this interesting observation, back-projection studies extended to intermediate frequencies, down to about 0.1 Hz, have suggested that a progressive transition of rupture properties as a function of frequency is observable. Here, by adapting the concept of array response function to non-stationary signals, we demonstrate that the “swimming artifact”, a systematic drift resulting from signal non-stationarity, induces significant bias on beamforming back-projection at low frequencies. We introduce a “reference window strategy” into the multitaper-MUSIC back-projection technique and significantly mitigate the “swimming artifact” at high frequencies (1 s to 4 s). At lower frequencies, this modification yields notable, but significantly smaller, artifacts than time-domain stacking. We perform extensive synthetic tests that include a 3D regional velocity model for Japan. We analyze the recordings of the Tohoku-Oki earthquake at the USArray and at the European array at periods from 1 s to 16 s. The migration of the source location as a function of period, regardless of the back-projection methods, has characteristics that are consistent with the expected effect of the “swimming artifact”. In particular, the apparent up-dip migration as a function of frequency obtained with the USArray can be explained by the “swimming artifact”. This indicates that the most substantial frequency-dependence of the Tohoku-Oki earthquake source occurs at periods longer than 16 s. Thus, low-frequency back-projection needs to be further tested and validated in order to contribute to the characterization of frequency-dependent rupture properties.

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.

The 2015 Mw 8.3 Illapel, Chile, earthquake: direction-reversed along-dip rupture with localized water reverberation

National Science Foundation of China 41374040 41090294 Hellman Fellowship University of California, Los Angeles (UCLA) Faculty Research Grant National Aeronautics and Space Administration (NASA) NAS7-03001 JPL Award 1468977

Application of Seismic Array Processing to Earthquake Early Warning

Earthquake early warning (EEW) systems that issue warnings prior to the arrival of strong shaking are essential in mitigating earthquake hazard. Currently operating EEW systems work on point‐source assumptions and are of limited effectiveness for large events, for which ignoring finite‐source effects result in magnitude underestimation. Here, we explore the concept of characterizing rupture dimensions in real time for EEW using small‐aperture seismic arrays located near active faults. Back tracing array waveforms allow estimation of the extent of the rupture front (as a proxy of the rupture size) and directivity in real time, providing complementary EEW capabilities for M>7 earthquakes to existing EEW systems. We implement it in a simulated real‐time environment and analyze the 2004 M 6 Parkfield, California, earthquake recordings by the U.S. Geological Survey Parkfield dense Seismograph ARray (UPSAR) array and the 2010 M 7.2 El Mayor–Cucapah earthquake recordings by strong‐motion sensors in San Diego, California. We find it important to correct for the bias in back azimuth induced by dipping structures beneath the UPSAR array, based on data from smaller events. Our estimated rupture length is 30% shorter than those inferred from other studies but still reasonable for EEW purposes. We attribute this difference to rupture directivity effects and the limited field of view of a single array. The accuracy of the approach may be improved with a network of arrays with overlapping fields of view. We demonstrate this by tracking the 2011 Tohoku earthquake rupture with two clusters of Hi‐net stations in Kyushu and northern Hokkaido. The obtained results are consistent with teleseismic back‐projection results and yield reasonable estimates of rupture length and directivity. Compared with other proposed finite‐fault EEW approaches, the array method is less affected by the coarseness of a Global Positioning System or seismic network and provides a high‐frequency characterization of the rupture that yields more suitable predictors of ground shaking for certain structures.

Transpressional Rupture Cascade of the 2016 Mw 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.