John Galetzka

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.

Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard

We document geodetic strain across the Nepal Himalaya using GPS times series from 30 stations in Nepal and southern Tibet, in addition to previously published campaign GPS points and leveling data and determine the pattern of interseismic coupling on the Main Himalayan Thrust fault (MHT). The noise on the daily GPS positions is modeled as a combination of white and colored noise, in order to infer secular velocities at the stations with consistent uncertainties. We then locate the pole of rotation of the Indian plate in the ITRF 2005 reference frame at longitude = − 1.34° ± 3.31°, latitude = 51.4° ± 0.3° with an angular velocity of Ω = 0.5029 ± 0.0072°/Myr. The pattern of coupling on the MHT is computed on a fault dipping 10° to the north and whose strike roughly follows the arcuate shape of the Himalaya. The model indicates that the MHT is locked from the surface to a distance of approximately 100 km down dip, corresponding to a depth of 15 to 20 km. In map view, the transition zone between the locked portion of the MHT and the portion which is creeping at the long term slip rate seems to be at the most a few tens of kilometers wide and coincides with the belt of midcrustal microseismicity underneath the Himalaya. According to a previous study based on thermokinematic modeling of thermochronological and thermobarometric data, this transition seems to happen in a zone where the temperature reaches 350°C. The convergence between India and South Tibet proceeds at a rate of 17.8 ± 0.5 mm/yr in central and eastern Nepal and 20.5 ± 1 mm/yr in western Nepal. The moment deficit due to locking of the MHT in the interseismic period accrues at a rate of 6.6 ± 0.4 × 10^(19) Nm/yr on the MHT underneath Nepal. For comparison, the moment released by the seismicity over the past 500 years, including 14 M_W ≥ 7 earthquakes with moment magnitudes up to 8.5, amounts to only 0.9 × 10^(19) Nm/yr, indicating a large deficit of seismic slip over that period or very infrequent large slow slip events. No large slow slip event has been observed however over the 20 years covered by geodetic measurements in the Nepal Himalaya. We discuss the magnitude and return period of M > 8 earthquakes required to balance the long term slip budget on the MHT.

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.

Deformation and Slip Along the Sunda Megathrust in the Great 2005 Nias-Simeulue Earthquake

Seismic rupture produced spectacular tectonic deformation above a 400-kilometer strip of the Sunda megathrust, offshore northern Sumatra, in March 2005. Measurements from coral microatolls and Global Positioning System stations reveal trench-parallel belts of uplift up to 3 meters high on the outer-arc islands above the rupture and a 1-meter-deep subsidence trough farther from the trench. Surface deformation reflects more than 11 meters of fault slip under the islands and a pronounced lessening of slip trenchward. A saddle in megathrust slip separates the northwestern edge of the 2005 rupture from the great 2004 Sumatra-Andaman rupture. The southeastern edge abuts a predominantly aseismic section of the megathrust near the equator.

Earthquake Supercycles Inferred from Sea-Level Changes Recorded in the Corals of West Sumatra

Records of relative sea-level change extracted from corals of the Mentawai islands, Sumatra, imply that this 700-kilometer-long section of the Sunda megathrust has generated broadly similar sequences of great earthquakes about every two centuries for at least the past 700 years. The moment magnitude 8.4 earthquake of September 2007 represents the first in a series of large partial failures of the Mentawai section that will probably be completed within the next several decades.

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.

Slip pulse and resonance of the Kathmandu basin during the 2015 Gorkha earthquake, Nepal

The bigger they are, the harder they fall The magnitude 7.8 Gorkha earthquake hit Nepal on 25 April 2015. The earthquake killed thousands and caused great damage. Galetzka et al. determined how the fault that caused this earthquake ruptured. The rupture showed a smooth slip pulse 20 km wide that moved eastward along the fault over about 6 s. The nature of the rupture limited damage to regular dwellings but generated shaking that collapsed taller structures. Science, this issue p. 1091 Continuous GPS and InSAR measurements record slip on the fault responsible for the 2015 Mw 7.8 Gorkha earthquake in Nepal. Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip onset, indicating a large (~5-meter) slip-weakening distance, caused moderate ground shaking at high frequencies (>1 hertz; peak ground acceleration, ~16% of Earth’s gravity) and minimized damage to vernacular dwellings. Whole-basin resonance at a period of 4 to 5 seconds caused the collapse of tall structures, including cultural artifacts.

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.

Paleogeodetic records of seismic and aseismic subduction from central Sumatran microatolls, Indonesia

We utilize coral microatolls in western Sumatra to document vertical deformation associated with subduction. Microatolls are very sensitive to fluctuations in sea level and thus act as natural tide gauges. They record not only the magnitude of vertical deformation associated with earthquakes (paleoseismic data), but also continuously track the long-term aseismic deformation that occurs during the intervals between earthquakes (paleogeodetic data). This paper focuses on the twentieth century paleogeodetic history of the equatorial region. Our coral paleogeodetic record of the 1935 event reveals a classical example of deformations produced by seismic rupture of a shallow subduction interface. The site closest to the trench rose 90 cm, whereas sites further east sank by as much as 35 cm. Our model reproduces these paleogeodetic data with a 2.3 m slip event on the interface 88 to 125 km from the trench axis. Our coral paleogeodetic data reveal slow submergence during the decades before and after the event in the areas of coseismic emergence. Likewise, interseismic emergence occurred before and after the 1935 event in areas of coseismic submergence. Among the interesting phenomenon we have discovered in the coral record is evidence of a large aseismic slip or “silent event” in 1962, 27 years after the 1935 event. Paleogeodetic deformation rates in the decades before, after, and between the 1935 and 1962 events have varied both temporally and spatially. During the 25 years following the 1935 event, submergence rates were dramatically greater than in prior decades. During the past four decades, however, rates have been lower than in the preceding decades, but are still higher than they were prior to 1935. These paleogeodetic records enable us to model the kinematics of the subduction interface throughout the twentieth century.

Coral evidence for earthquake recurrence and an A.D. 1390-1455 cluster at the south end of the 2004 Aceh–Andaman rupture

Coral records of relative sea level change provide a history of vertical interseismic and coseismic deformation along the coast of northern Simeulue Island, Sumatra, and reveal details about earthquakes in the 10th and 14th–15th centuries A.D. along the southern end of the December 2004 M_w 9.2 Sunda megathrust rupture. Over a 56 year period between A.D. 1390 and 1455, northern Simeulue experienced a cluster of megathrust ruptures, associated with total uplift that was considerably more than in 2004. Uplifted corals at two sites constrain the first event of the cluster to A.D. 1393 ± 3 and 1394 ± 2 (2σ). A smaller but well-substantiated uplift occurred in northern Simeulue in 1430 ± 3. An inferred third uplift, in A.D. 1450 ± 3, killed all corals on the reef flats of northern Simeulue. The amount of uplift during this third event, though confirmed only to have exceeded 28 and 41 cm at two sites, probably surpassed the 100 and 44 cm that occurred, respectively, at those sites in 2004, and it was likely more than in 2004 over all of northern Simeulue. The evidence for past earthquake clustering combined with the inference of considerably greater uplift in A.D. 1390–1455 than in 2004 suggests that strain may still be stored along the southernmost part of the 2004 rupture. Interseismic subsidence rates recorded by northern Simeulue coral microatolls have varied by up to a factor of 4 at some sites from one earthquake cycle to another.