Kinematic fault slip evolution source models of the 2008 M7.9 Wenchuan earthquake in China from SAR interferometry, GPS and teleseismic analysis and implications for Longmen Shan tectonics

The M_w 7.9 2008 Wenchuan earthquake ruptured about 280 km of faults in the Longmen Shan of Sichuan province, China, at the eastern edge of the Tibetan Plateau. We use teleseismic waveforms with geodetic data from Global Positioning System, synthetic aperture radar interferometry and image amplitude correlation to produce a source model of this earthquake. The model describes evolution of fault slip during the earthquake. The geodetic data constrains the spatial distribution of fault slip and the seismic waveforms constrain mostly the time evolution of slip. We find that the earthquake started with largely thrust motion on an imbricate system of faults beneath the central Longmen Shan, including the Beichuan Fault and Pengguan Fault, with fault slip at depth extending up to 50 km northwest of the mountain front. The fault ruptures continued northeast along the Beichuan Fault with more oblique slip (right-lateral and thrust) and the proportion of lateral motion increasing in the northern Longmen Shan. The northernmost fault segment has a much steeper dip, consistent with nearly pure strike-slip motion. The kinematic source model shows that the rupture propagated to the northeast at about 2.5–3.0 km s^(−1), producing a cascade of subevents with a total duration of about 110 s. The complex fault ruptures caused shortening and uplift of the extremely steep central Longmen Shan, which supports models where the steep edge of the plateau is formed by thrusting over the strong crust of the Sichuan Basin

A detailed source model for the M_w9.0 Tohoku-Oki earthquake reconciling geodesy, seismology, and tsunami records

The 11 March 2011 M_w9.0 Tohoku-Oki earthquake was recorded by an exceptionally large amount of diverse data offering a unique opportunity to investigate the details of this major megathrust rupture. Many studies have taken advantage of the very dense Japanese onland strong motion, broadband, and continuous GPS networks in this sense. But resolution tests and the variability in the proposed solutions have highlighted the difficulty to uniquely resolve the slip distribution from these networks, relatively distant from the source region, and with limited azimuthal coverage. In this context, we present a finite fault slip joint inversion including an extended amount of complementary data (teleseismic, strong motion, high-rate GPS, static GPS, seafloor geodesy, and tsunami records) in an attempt to reconcile them into a single better resolved model. The inversion reveals a patchy slip distribution with large slip (up to 64 m) mostly located updip of the hypocenter and near the trench. We observe that most slip is imaged in a region where almost no earthquake was recorded before the main shock and around which intense interplate seismicity is observed afterward. At a smaller scale, the largest slip pattern is imaged just updip of an important normal fault coseismically activated. This normal fault has been shown to be the mark of very low dynamic friction allowing extremely large slip to propagate up to the free surface. The spatial relationship between this normal fault and our slip distribution strengthens its key role in the rupture process of the Tohoku-Oki earthquake

Joint Inversion of Coseismic and Early Postseismic Slip to Optimize the Information Content in Geodetic Data: Application to the 2009 M_w 6.3 L'Aquila Earthquake, Central Italy

When analyzing the rupture of a large earthquake, geodetic data are often critical. These data are generally characterized by either a good temporal or a good spatial resolution, but rarely both. As a consequence, many studies analyze the coseismic rupture with data that also include one or more days of early postseismic deformation. Here, we invert simultaneously for the coseismic and postseismic slip with the condition that the sum of the two models remains compatible with data covering the two slip episodes. We validate the benefits of this approach with a toy model and an application to the 2009 M_w 6.3 L'Aquila earthquake, using a Bayesian approach and accounting for epistemic uncertainties. For the L'Aquila earthquake, we find that if early postseismic deformation is not an explicitly acknowledged coseismic signal, coseismic slip models may overestimate the peak amplitude while long‐term postseismic models may largely underestimate the total postseismic slip amplitude. This example illustrates how the proposed approach could improve our comprehension of the seismic cycle, fault frictional properties, and the spatial and temporal relationship between seismic rupture, afterslip, and aftershocks

Superficial simplicity of the 2010 El Mayor–Cucapah earthquake of Baja California in Mexico

The geometry of faults is usually thought to be more complicated at the surface than at depth and to control the initiation, propagation and arrest of seismic ruptures. The fault system that runs from southern California into Mexico is a simple strike-slip boundary: the west side of California and Mexico moves northwards with respect to the east. However, the M_w 7.2 2010 El Mayor–Cucapah earthquake on this fault system produced a pattern of seismic waves that indicates a far more complex source than slip on a planar strike-slip fault. Here we use geodetic, remote-sensing and seismological data to reconstruct the fault geometry and history of slip during this earthquake. We find that the earthquake produced a straight 120-km-long fault trace that cut through the Cucapah mountain range and across the Colorado River delta. However, at depth, the fault is made up of two different segments connected by a small extensional fault. Both segments strike N130° E, but dip in opposite directions. The earthquake was initiated on the connecting extensional fault and 15 s later ruptured the two main segments with dominantly strike-slip motion. We show that complexities in the fault geometry at depth explain well the complex pattern of radiated seismic waves. We conclude that the location and detailed characteristics of the earthquake could not have been anticipated on the basis of observations of surface geology alone

Coseismic Slip and Afterslip of the Great M_w 9.15 Sumatra–Andaman Earthquake of 2004

We determine coseismic and the first-month postseismic deformation associated with the Sumatra–Andaman earthquake of 26 December 2004 from near- field Global Positioning System (GPS) surveys in northwestern Sumatra and along the Nicobar-Andaman islands, continuous and campaign GPS measurements from Thailand and Malaysia, and in situ and remotely sensed observations of the vertical motion of coral reefs. The coseismic model shows that the Sunda subduction megathrust ruptured over a distance of about 1500 km and a width of less than 150 km, releasing a total moment of 6.7–7.0 x 10^(22) N m, equivalent to a magnitude M_w 9.15. The latitudinal distribution of released moment in our model has three distinct peaks at about 4° N, 7° N, and 9° N, which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated T waves. Our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very-long-period surface waves. The tsunami predicted from this model fits relatively well the altimetric measurements made by the JASON and TOPEX satellites. Neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave. The near-field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500-sec-long seismic rupture. The postseismic geodetic moment of about 2.4 x 10^(22) N m (M_w 8.8) is equal to about 30 ± 5% of the coseismic moment release. Evolution of postseismic deformation is consistent with rate-strengthening frictional afterslip

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

Geophysical observations from the 2011 moment magnitude (M_w) 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 M_w 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

Partial rupture of a locked patch of the Sumatra megathrust during the 2007 earthquake sequence

The great Sumatra–Andaman earthquake and tsunami of 2004 was a dramatic reminder of the importance of understanding the seismic and tsunami hazards of subduction zones [1,2,3,4]. In March 2005, the Sunda megathrust ruptured again, producing an event [5] of moment magnitude (Mw) 8.6 south of the 2004 rupture area, which was the site of a similar event in 1861 (ref. 6). Concern was then focused on the Mentawai area, where large earthquakes had occurred in 1797 (Mw = 8.8) and 1833 (Mw = 9.0) [6,7]. Two earthquakes, one of Mw = 8.4 and, twelve hours later, one of Mw = 7.9, indeed occurred there on 12 September 2007. Here we show that these earthquakes ruptured only a fraction of the area ruptured in 1833 and consist of distinct asperities within a patch of the megathrust that had remained locked in the interseismic period. This indicates that the same portion of a megathrust can rupture in different patterns depending on whether asperities break as isolated seismic events or cooperate to produce a larger rupture. This variability probably arises from the influence of non-permanent barriers, zones with locally lower pre-stress due to the past earthquakes. The stress state of the portion of the Sunda megathrust that had ruptured in 1833 and 1797 was probably not adequate for the development of a single large rupture in 2007. The moment released in 2007 amounts to only a fraction both of that released in 1833 and of the deficit of moment that had accumulated as a result of interseismic strain since 1833. The potential for a large megathrust event in the Mentawai area thus remains large

Seismic and aseismic slip on the Central Peru megathrust

Slip on a subduction megathrust can be seismic or aseismic, with the two modes of slip complementing each other in time and space to accommodate the long-term plate motions. Although slip is almost purely aseismic at depths greater than about 40 km, heterogeneous surface strain suggests that both modes of slip occur at shallower depths, with aseismic slip resulting from steady or transient creep in the interseismic and postseismic periods. Thus, active faults seem to comprise areas that slip mostly during earthquakes, and areas that mostly slip aseismically. The size, location and frequency of earthquakes that a megathrust can generate thus depend on where and when aseismic creep is taking place, and what fraction of the long-term slip rate it accounts for. Here we address this issue by focusing on the central Peru megathrust. We show that the Pisco earthquake, with moment magnitude M_w = 8.0, ruptured two asperities within a patch that had remained locked in the interseismic period, and triggered aseismic frictional afterslip on two adjacent patches. The most prominent patch of afterslip coincides with the subducting Nazca ridge, an area also characterized by low interseismic coupling, which seems to have repeatedly acted as a barrier to seismic rupture propagation in the past. The seismogenic portion of the megathrust thus appears to be composed of interfingering rate-weakening and rate-strengthening patches. The rate-strengthening patches contribute to a high proportion of aseismic slip, and determine the extent and frequency of large interplate earthquakes. Aseismic slip accounts for as much as 50–70% of the slip budget on the seismogenic portion of the megathrust in central Peru, and the return period of earthquakes with M_w = 8.0 in the Pisco area is estimated to be 250  years