M. Chlieh

Frictional afterslip following the 2005 Nias-Simeulue earthquake, Sumatra

Continuously recording Global Positioning System stations near the 28 March 2005 rupture of the Sunda megathrust [moment magnitude (Mw) 8.7] show that the earthquake triggered aseismic frictional afterslip on the subduction megathrust, with a major fraction of this slip in the up-dip direction from the main rupture. Eleven months after the main shock, afterslip continues at rates several times the average interseismic rate, resulting in deformation equivalent to at least a Mw 8.2 earthquake. In general, along-strike variations in frictional behavior appear to persist over multiple earthquake cycles. Aftershocks cluster along the boundary between the region of coseismic slip and the up-dip creeping zone. We observe that the cumulative number of aftershocks increases linearly with postseismic displacements; this finding suggests that the temporal evolution of aftershocks is governed by afterslip.

Plate-boundary deformation associated with the great Sumatra-Andaman earthquake

The Sumatra–Andaman earthquake of 26 December 2004 is the first giant earthquake (moment magnitude Mw > 9.0) to have occurred since the advent of modern space-based geodesy and broadband seismology. It therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events. Here we report estimates of the ground displacement associated with this event, using near-field Global Positioning System (GPS) surveys in northwestern Sumatra combined with in situ and remote observations of the vertical motion of coral reefs. These data show that the earthquake was generated by rupture of the Sunda subduction megathrust over a distance of >1,500 kilometres and a width of <150 kilometres. Megathrust slip exceeded 20 metres offshore northern Sumatra, mostly at depths shallower than 30 kilometres. Comparison of the geodetically and seismically inferred slip distribution indicates that ∼30 per cent additional fault slip accrued in the 1.5 months following the 500-second-long seismic rupture. Both seismic and aseismic slip before our re-occupation of GPS sites occurred on the shallow portion of the megathrust, where the large Aceh tsunami originated. Slip tapers off abruptly along strike beneath Simeulue Island at the southeastern edge of the rupture, where the earthquake nucleated and where an Mw = 7.2 earthquake occurred in late 2002. This edge also abuts the northern limit of slip in the 28 March 2005 Mw = 8.7 Nias–Simeulue earthquake.

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. In March 2005, the Sunda megathrust ruptured again, producing an event 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). 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.

Heterogeneous coupling of the Sumatran megathrust constrained by geodetic and paleogeodetic measurements

Geodetic and paleogeodetic measurements of interseismic strain above the Sumatran portion of the Sunda subduction zone reveal a heterogeneous pattern of coupling. Annual banding in corals provides vertical rates of deformation spanning the last half of the 20th century, and repeated GPS surveys between 1991 and 2001 and continuous measurements at GPS stations operated since 2002 provide horizontal velocities. Near the equator, the megathrust is locked over a narrow width of only a few tens of kilometers. In contrast, the locked fault zone is up to about 175 km wide in areas where great interplate earthquakes have occurred in the past. Formal inversion of the data reveals that these strongly coupled patches are roughly coincident with asperities that ruptured during these events. The correlation is most spectacular for rupture of the M_w 8.7 Nias-Simeulue earthquake of 2005, which released half of the moment deficit that had accumulated since its previous rupture in 1861, suggesting that this earthquake was overdue. Beneath the Mentawai islands, strong coupling is observed within the overlapping rupture areas of the great earthquakes of 1797 and 1833. The accumulated slip deficit since these events is slowly reaching the amount of slip that occurred during the 1833 earthquake but already exceeds the slip that occurred during the 1797 earthquake. Thus, rerupture of part of the Mentawai patch in September 2007 was not a surprise. In contrast, coupling is low below the Batu islands near the equator and around Enggano island at about 5°S, where only moderate earthquakes (M_w < 8.0) have occurred in the past two centuries. The correlation of large seismic asperities with patches that are locked during the interseismic period suggests that they are persistent features. This interpretation is reinforced by the fact that the large locked patches and great ruptures occur beneath persistent geomorphologic features, the largest outer arc islands. Depth- and convergence-rate-dependent temperature might influence the pattern of coupling, through its effect on the rheology of the plate interface, but other influences are required to account for the observed along-strike heterogeneity of coupling. In particular, subduction of the Investigator Fracture Zone could be the cause for the low coupling near the equator.

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 Mw = 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 Mw = 8.0 in the Pisco area is estimated to be 250 years.

Source parameters of the great Sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls

Large uplifts and tilts occurred on the Sumatran outer arc islands between 0.5° and 3.3°S during great historical earthquakes in 1797 and 1833, as judged from relative sea level changes recorded by annually banded coral heads. Coral data for these two earthquakes are most complete along a 160-km length of the Mentawai islands between 3.2° and 2°S. Uplift there was as great as 0.8 m in 1797 and 2.8 m in 1833. Uplift in 1797 extended 370 km, between 3.2° and 0.5°S. The pattern and magnitude of uplift imply megathrust ruptures corresponding to moment magnitudes (M_w) in the range 8.5 to 8.7. The region of uplift in 1833 ranges from 2° to at least 3.2°S and, judging from historical reports of shaking and tsunamis, perhaps as far as 5°S. The patterns and magnitude of uplift and tilt in 1833 are similar to those experienced farther north, between 0.5° and 3°N, during the giant Nias-Simeulue megathrust earthquake of 2005; the outer arc islands rose as much as 3 m and tilted toward the mainland. Elastic dislocation forward modeling of the coral data yields megathrust ruptures with moment magnitudes ranging from 8.6 to 8.9. Sparse accounts at Padang, along the mainland west coast at latitude 1°S, imply tsunami runups of at least 5 m in 1797 and 3–4 m in 1833. Tsunamis simulated from the pattern of coral uplift are roughly consistent with these reports. The tsunami modeling further indicates that the Indian Ocean tsunamis of both 1797 and 1833, unlike that of 2004, were directed mainly south of the Indian subcontinent. Between about 0.7° and 2.1°S, the lack of vintage 1797 and 1833 coral heads in the intertidal zone demonstrates that interseismic submergence has now nearly equals coseismic emergence that accompanied those earthquakes. The interseismic strains accumulated along this reach of the megathrust have thus approached or exceeded the levels relieved in 1797 and 1833.

Interseismic coupling and seismic potential along the Central Andes subduction zone

We use about two decades of geodetic measurements to characterize interseismic strain build up along the Central Andes subduction zone from Lima, Peru, to Antofagasta, Chile. These measurements are modeled assuming a 3-plate model (Nazca, Andean sliver and South America Craton) and spatially varying interseismic coupling (ISC) on the Nazca megathrust interface. We also determine slip models of the 1996 M_w = 7.7 Nazca, the 2001 M_w = 8.4 Arequipa, the 2007 M_w = 8.0 Pisco and the M_w = 7.7 Tocopilla earthquakes. We find that the data require a highly heterogeneous ISC pattern and that, overall, areas with large seismic slip coincide with areas which remain locked in the interseismic period (with high ISC). Offshore Lima where the ISC is high, a M_w∼8.6–8.8 earthquake occurred in 1746. This area ruptured again in a sequence of four M_w∼8.0 earthquakes in 1940, 1966, 1974 and 2007 but these events released only a small fraction of the elastic strain which has built up since 1746 so that enough elastic strain might be available there to generate a M_w > 8.5 earthquake. The region where the Nazca ridge subducts appears to be mostly creeping aseismically in the interseismic period (low ISC) and seems to act as a permanent barrier as no large earthquake ruptured through it in the last 500 years. In southern Peru, ISC is relatively high and the deficit of moment accumulated since the M_w∼8.8 earthquake of 1868 is equivalent to a magnitude M_w∼8.4 earthquake. Two asperities separated by a subtle aseismic creeping patch are revealed there. This aseismic patch may arrest some rupture as happened during the 2001 Arequipa earthquake, but the larger earthquakes of 1604 and 1868 were able to rupture through it. In northern Chile, ISC is very high and the rupture of the 2007 Tocopilla earthquake has released only 4% of the elastic strain that has accumulated since 1877. The deficit of moment which has accumulated there is equivalent to a magnitude M_w∼8.7 earthquake. This study thus provides elements to assess the location, size and magnitude of future large megathurst earthquakes in the Central Andes subduction zone. Caveats of this study are that interseismic strain of the forearc is assumed time invariant and entirely elastic. Also a major source of uncertainty is due to fact that the available data place very little constraints on interseismic coupling at shallow depth near the trench, except offshore Lima where sea bottom geodetic measurements have been collected suggesting strong coupling.

Tsunami inundation modeling for western Sumatra

A long section of the Sunda megathrust south of the great tsunamigenic earthquakes of 2004 and 2005 is well advanced in its seismic cycle and a plausible candidate for rupture in the next few decades. Our computations of tsunami propagation and inundation yield model flow depths and inundations consistent with sparse historical accounts for the last great earthquakes there, in 1797 and 1833. Numerical model results from plausible future ruptures produce flow depths of several meters and inundation up to several kilometers inland near the most populous coastal cities. Our models of historical and future tsunamis confirm a substantial exposure of coastal Sumatran communities to tsunami surges. Potential losses could be as great as those that occurred in Aceh in 2004.

Interseismic deformation above the Sunda Megathrust recorded in coral microatolls of the Mentawai islands, West Sumatra

past half century vary from 2 to 14 mm yr 1 and increase southwestward, toward the subduction trench. The pattern is consistent with rates of subsidence measured by a sparse network of continuously recording Global Positioning System (cGPS) stations and with locking of a 400-km-long section of the underlying subduction megathrust, between about 1S and 4S. This record of subsidence and tilting, extending nearly a century into the past, implies that the region is advancing toward the occurrence of another giant earthquake. However, evidence of episodic rather than steady subsidence reflects a behavior that is more complex than simple elastic strain accumulation and relief. Most prominent of these episodes is an extensive emergence/subsidence couplet in about 1962, which may be the result of rapid, aseismic slip on the megathrust, between the islands and the trench. Lower subsidence rates recorded by the corals since about 1985 may reflect failure on many small patches within the locked section of the megathrust.

Postseismic relocking of the subduction megathrust following the 2007 Pisco, Peru, earthquake

Characterizing the time evolution of slip over different phases of the seismic cycle is crucial to a better understanding of the factors controlling the occurrence of large earthquakes. In this study, we take advantage of interferometric synthetic aperture radar data and 3.5 years of continuous Global Positioning System (GPS) measurements to determine interseismic, coseismic, and postseismic slip distributions in the region of the 2007, Mw 8.0 Pisco, earthquake, Peru, using the same fault geometry and inversion method. Our interseismic model, based on pre-2007 campaign GPS data, suggests that the 2007 Pisco seismic slip occurred in a region strongly coupled before the earthquake while afterslip occurred in low coupled regions. Large afterslip occurred in the peripheral area of coseismic rupture in agreement with the notion that afterslip is mainly induced by coseismic stress changes. The temporal evolution of the region of maximum afterslip, characterized by a relaxation time of about 2.3 years, is located in the region where the Nazca ridge is subducting, consistent with rate-strengthening friction promoting aseismic slip. We estimate a return period for the Pisco earthquake of about 230 years with an estimated aseismic slip that might account for about 50% of the slip budget in this region over the 0–50 km seismogenic depth range. A major result of this study is that the main asperity that ruptured during the 2007 Pisco earthquake relocked soon after this event.