Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake

Most large earthquakes occur along an oceanic trench, where an oceanic plate subducts beneath a continental plate. Massive earthquakes with a moment magnitude, M(w), of nine have been known to occur in only a few areas, including Chile, Alaska, Kamchatka and Sumatra. No historical records exist of a M(w) = 9 earthquake along the Japan trench, where the Pacific plate subducts beneath the Okhotsk plate, with the possible exception of the ad 869 Jogan earthquake, the magnitude of which has not been well constrained. However, the strain accumulation rate estimated there from recent geodetic observations is much higher than the average strain rate released in previous interplate earthquakes. This finding raises the question of how such areas release the accumulated strain. A megathrust earthquake with M(w) = 9.0 (hereafter referred to as the Tohoku-Oki earthquake) occurred on 11 March 2011, rupturing the plate boundary off the Pacific coast of northeastern Japan. Here we report the distributions of the coseismic slip and postseismic slip as determined from ground displacement detected using a network based on the Global Positioning System. The coseismic slip area extends approximately 400?km along the Japan trench, matching the area of the pre-seismic locked zone. The afterslip has begun to overlap the coseismic slip area and extends into the surrounding region. In particular, the afterslip area reached a depth of approximately 100?km, with M(w) = 8.3, on 25 March 2011. Because the Tohoku-Oki earthquake released the strain accumulated for several hundred years, the paradox of the strain budget imbalance may be partly resolved. This earthquake reminds us of the potential for M(w)?≈?9 earthquakes to occur along other trench systems, even if no past evidence of such events exists. Therefore, it is imperative that strain accumulation be monitored using a space geodetic technique to assess earthquake potential.     

Seismology: earthquake risk on the Sunda trench

On 28 March 2005 the Sunda megathrust in Indonesia ruptured again, producing another great earthquake three months after the previous one. The rupture was contiguous with that of the December 2004 Sumatra-Andaman earthquake, and is likely to have been sparked by local stress, although the triggering stresses at its hypocentre were very small - of the order of just 0.1 bar. Calculations show that stresses imposed by the second rupture have brought closer to failure the megathrust immediately to the south, under the Batu and Mentawai islands, and have expanded the area of increased stress on the Sumatra fault. Palaeoseismologic studies show that the Mentawai segment of the Sunda megathrust is well advanced in its seismic cycle and is therefore a good candidate for triggered failure.     

Tracking the rupture of the Mw = 9.3 Sumatra earthquake over 1,150 km at teleseismic distance

On 26 December 2004, a moment magnitude Mw = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region. The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave-height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN), located approximately 9,000 km from the event. Here we present a modification of a standard array-seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.     

Simultaneous teleseismic and geodetic observations of the stick-slip motion of an Antarctic ice stream

Long-period seismic sources associated with glacier motion have been recently discovered, and an increase in ice flow over the past decade has been suggested on the basis of secular changes in such measurements. Their significance, however, remains uncertain, as a relationship to ice flow has not been confirmed by direct observation. Here we combine long-period surface-wave observations with simultaneous Global Positioning System measurements of ice displacement to study the tidally modulated stick-slip motion of the Whillans Ice Stream in West Antarctica. The seismic origin time corresponds to slip nucleation at a region of the bed of the Whillans Ice Stream that is likely stronger than in surrounding regions and, thus, acts like an 'asperity' in traditional fault models. In addition to the initial pulse, two seismic arrivals occurring 10-23 minutes later represent stopping phases as the slip terminates at the ice stream edge and the grounding line. Seismic amplitude and average rupture velocity are correlated with tidal amplitude for the different slip events during the spring-to-neap tidal cycle. Although the total seismic moment calculated from ice rigidity, slip displacement, and rupture area is equivalent to an earthquake of moment magnitude seven (M(w) 7), seismic amplitudes are modest (M(s) 3.6-4.2), owing to the source duration of 20-30 minutes. Seismic radiation from ice movement is proportional to the derivative of the moment rate function at periods of 25-100 seconds and very long-period radiation is not detected, owing to the source geometry. Long-period seismic waves are thus useful for detecting and studying sudden ice movements but are insensitive to the total amount of slip.     

Plateau 'pop-up' in the great 1897 Assam earthquake

The great Assam earthquake of 12 June 1897 reduced to rubble all masonry buildings within a region of northeastern India roughly the size of England, and was felt over an area exceeding that of the great 1755 Lisbon earthquake. Hitherto it was believed that rupture occurred on a north-dipping Himalayan thrust fault propagating south of Bhutan. But here we show that the northern edge of the Shillong plateau rose violently by at least 11 m during the Assam earthquake, and that this was due to the rupture of a buried reverse fault approximately 110 km in length and dipping steeply away from the Himalaya. The stress drop implied by the rupture geometry and the prodigious fault slip of 18 +/- 7 m explains epicentral accelerations observed to exceed 1g vertically and surface velocities exceeding 3 m s-1 (ref. 1). This quantitative observation of active deformation of a 'pop-up' structure confirms that faults bounding such structures can penetrate the whole crust. Plateau uplift in the past 2-5 million years has caused the Indian plate to contract locally by 4 +/- 2 mm yr-1, reducing seismic risk in Bhutan but increasing the risk in northern Bangladesh.     

Earthquake nucleation by transient deformations caused by the M = 7.9 Denali, Alaska, earthquake

The permanent and dynamic (transient) stress changes inferred to trigger earthquakes are usually orders of magnitude smaller than the stresses relaxed by the earthquakes themselves, implying that triggering occurs on critically stressed faults. Triggered seismicity rate increases may therefore be most likely to occur in areas where loading rates are highest and elevated pore pressures, perhaps facilitated by high-temperature fluids, reduce frictional stresses and promote failure. Here we show that the 2002 magnitude M = 7.9 Denali, Alaska, earthquake triggered widespread seismicity rate increases throughout British Columbia and into the western United States. Dynamic triggering by seismic waves should be enhanced in directions where rupture directivity focuses radiated energy, and we verify this using seismic and new high-sample GPS recordings of the Denali mainshock. These observations are comparable in scale only to the triggering caused by the 1992 M = 7.4 Landers, California, earthquake, and demonstrate that Landers triggering did not reflect some peculiarity of the region or the earthquake. However, the rate increases triggered by the Denali earthquake occurred in areas not obviously tectonically active, implying that even in areas of low ambient stressing rates, faults may still be critically stressed and that dynamic triggering may be ubiquitous and unpredictable.     

Insight into the 2004 Sumatra-Andaman earthquake from GPS measurements in southeast Asia

Data collected at approximately 60 Global Positioning System (GPS) sites in southeast Asia show the crustal deformation caused by the 26 December 2004 Sumatra-Andaman earthquake at an unprecedented large scale. Small but significant co-seismic jumps are clearly detected more than 3,000 km from the earthquake epicentre. The nearest sites, still more than 400 km away, show displacements of 10 cm or more. Here we show that the rupture plane for this earthquake must have been at least 1,000 km long and that non-homogeneous slip is required to fit the large displacement gradients revealed by the GPS measurements. Our kinematic analysis of the GPS recordings indicates that the centroid of released deformation is located at least 200 km north of the seismological epicentre. It also provides evidence that the rupture propagated northward sufficiently fast for stations in northern Thailand to have reached their final positions less than 10 min after the earthquake, hence ruling out the hypothesis of a silent slow aseismic rupture.     

Seismic reflection imaging of two megathrust shear zones in the northern Cascadia subduction zone

At convergent continental margins, the relative motion between the subducting oceanic plate and the overriding continent is usually accommodated by movement along a single, thin interface known as a megathrust. Great thrust earthquakes occur on the shallow part of this interface where the two plates are locked together. Earthquakes of lower magnitude occur within the underlying oceanic plate, and have been linked to geochemical dehydration reactions caused by the plate's descent. Here I present deep seismic reflection data from the northern Cascadia subduction zone that show that the inter-plate boundary is up to 16 km thick and comprises two megathrust shear zones that bound a >5-km-thick, approximately 110-km-wide region of imbricated crustal rocks. Earthquakes within the subducting plate occur predominantly in two geographic bands where the dip of the plate is inferred to increase as it is forced around the edges of the imbricated inter-plate boundary zone. This implies that seismicity in the subducting slab is controlled primarily by deformation in the upper part of the plate. Slip on the shallower megathrust shear zone, which may occur by aseismic slow slip, will transport crustal rocks into the upper mantle above the subducting oceanic plate and may, in part, provide an explanation for the unusually low seismic wave speeds that are observed there.     

Earthquakes as beacons of stress change

Aftershocks occurring on faults in the far-field of a large earthquake rupture can generally be accounted for by changes in static stress on these faults caused by the rupture. This implies that faults interact, and that the timing of an earthquake can be affected by previous nearby ruptures. Here we explore the potential of small earthquakes to act as 'beacons' for the mechanical state of the crust. We investigate the static-stress changes resulting from the 1992 Landers earthquake in southern California which occurred in an area of high seismic activity stemming from many faults. We first gauge the response of the regional seismicity to the Landers event with a new technique, and then apply the same method to the inverse problem of determining the slip distribution on the main rupture from the seismicity. Assuming justifiable parameters, we derive credible matches to slip profiles obtained directly from the Landers mainshock. Our results provide a way to monitor mechanical conditions in the upper crust, and to investigate processes leading to fault failure.     

Rupture of the 2004 Sumatra-Andaman earthquake inferred from direct P-wave imaging

The Sumatra-Andaman earthquake on December 26, 2004 is the first well recorded gigantic earthquake (moment magnitude MW 9.3) by modern broadband seismic and Global Positioning System networks. The rich seismic and geodetic recordings have documented unprecedented details about the earthquake rupture, coseismic and postseismic deformations. This is a report of detailed images of the rupture process using the first-arriving compressional waves recorded by the China National Digital Seismic Network (CNDSN). An improved imaging condition was employed to account for the sparse distribution of the CNDSN stations. The resulting images are consistent with the major rupture features reported by previous seismic and geodetic studies. It is found that the earthquake rupture initiated at offshore of northwestern Sumatra and propagated in the north northwest direction at a speed of 2.7 ± 0.2 km/s. The rupture continued for at least 420 s and extended about 1200-1300 km along the Andaman trough with two bursts of seismic energy.     

The potential for giant tsunamigenic earthquakes in the northern Bay of Bengal

The great Sumatra-Andaman earthquake and Indian Ocean tsunami of 2004 came as a surprise to most of the earth science community. Although it is now widely recognized that the risk of another giant earthquake is high off central Sumatra, just east of the 2004 earthquake, there seems to be relatively little concern about the subduction zone to the north, in the northern Bay of Bengal along the coast of Myanmar. Here I show that similar indicators suggest a high potential for giant earthquakes along the coast of Myanmar. These indicators include the tectonic environment, which is similar to other subduction zones that experience giant megathrust earthquakes, stress and crustal strain observations, which indicate that the seismogenic zone is locked, and historical earthquake activity, which indicates that giant tsunamigenic earthquakes have occurred there in the past. These are all consistent with active subduction in the Myanmar subduction zone and I suggest that the seismogenic zone extends beneath the Bengal Fan. I conclude therefore that giant earthquakes probably occur off the coast of Myanmar, and that a large and vulnerable population is thereby exposed to a significant earthquake and tsunami hazard.     

Predicting the endpoints of earthquake ruptures

The active fault traces on which earthquakes occur are generally not continuous, and are commonly composed of segments that are separated by discontinuities that appear as steps in map-view. Stress concentrations resulting from slip at such discontinuities may slow or stop rupture propagation and hence play a controlling role in limiting the length of earthquake rupture. Here I examine the mapped surface rupture traces of 22 historical strike-slip earthquakes with rupture lengths ranging between 10 and 420 km. I show that about two-thirds of the endpoints of strike-slip earthquake ruptures are associated with fault steps or the termini of active fault traces, and that there exists a limiting dimension of fault step (3-4 km) above which earthquake ruptures do not propagate and below which rupture propagation ceases only about 40 per cent of the time. The results are of practical importance to seismic hazard analysis where effort is spent attempting to place limits on the probable length of future earthquakes on mapped active faults. Physical insight to the dynamics of the earthquake rupture process is further gained with the observation that the limiting dimension appears to be largely independent of the earthquake rupture length. It follows that the magnitude of stress changes and the volume affected by those stress changes at the driving edge of laterally propagating ruptures are largely similar and invariable during the rupture process regardless of the distance an event has propagated or will propagate.     

Static slip model of the M w 9.0 Tohoku (Japan) earthquake: Results from joint inversion of terrestrial GPS data and seafloor GPS/acoustic data

Based on co-seismic displacements recorded by terrestrial GPS stations and seafloor GPS/acoustic stations, the static slip model of the 2011 Mw 9.0 Tohoku earthquake was determined by inverting the data using a layered earth model. According to a priori information, the rupture surface was modeled with a geometry that is close to the actual rupture, in which the fault dip angle increases with depth and the fault strike varies with the trend of the trench. As shown by the results inferred from the joint inversion, the "geodetic" moment is 3.68 × 10 22 Nm, corresponding to Mw 9.01, and the maximum slip is positioned at a depth of 13.5 km with a slip magnitude of 45.8 m. Rupture asperities with slip exceeding 10 m are mainly distributed from 39.6 to 36.97°N, over a length of almost 240 km along the trench. The slip was mostly concentrated at depths shallower than 40 km, up-dip of the hypocenter. "Checkerboard" tests reveal that a joint inversion of multiple datasets can resolve the slip distribution better than an inversion with terrestrial GPS data only-especially when aiming to resolve slip at shallow depths. Thus, the joint inversion results obtained by this work may provide a more reliable slip model than the results of other studies that are only derived from terrestrial GPS data or seismic waveform data.     

Slip zone and energetics of a large earthquake from the Taiwan Chelungpu-fault Drilling Project

Determining the seismic fracture energy during an earthquake and understanding the associated creation and development of a fault zone requires a combination of both seismological and geological field data. The actual thickness of the zone that slips during the rupture of a large earthquake is not known and is a key seismological parameter in understanding energy dissipation, rupture processes and seismic efficiency. The 1999 magnitude-7.7 earthquake in Chi-Chi, Taiwan, produced large slip (8 to 10 metres) at or near the surface, which is accessible to borehole drilling and provides a rare opportunity to sample a fault that had large slip in a recent earthquake. Here we present the retrieved cores from the Taiwan Chelungpu-fault Drilling Project and identify the main slip zone associated with the Chi-Chi earthquake. The surface fracture energy estimated from grain sizes in the gouge zone of the fault sample was directly compared to the seismic fracture energy determined from near-field seismic data. From the comparison, the contribution of gouge surface energy to the earthquake breakdown work is quantified to be 6 per cent.     

Resonant slow fault slip in subduction zones forced by climatic load stress

Global Positioning System (GPS) measurements at subduction plate boundaries often record fault movements similar to earthquakes but much slower, occurring over timescales of approximately 1 week to approximately 1 year. These 'slow slip events' have been observed in Japan, Cascadia, Mexico, Alaska and New Zealand. The phenomenon is poorly understood, but several observations hint at the processes underlying slow slip. Although slip itself is silent, seismic instruments often record coincident low-amplitude tremor in a narrow (1-5 cycles per second) frequency range. Also, modelling of GPS data and estimates of tremor location indicate that slip focuses near the transition from unstable ('stick-slip') to stable friction at the deep limit of the earthquake-producing seismogenic zone. Perhaps most intriguingly, slow slip is periodic at several locations, with recurrence varying from 6 to 18 months depending on which subduction zone (or even segment) is examined. Here I show that such periodic slow fault slip may be a resonant response to climate-driven stress perturbations. Fault slip resonance helps to explain why slip events are periodic, why periods differ from place to place, and why slip focuses near the base of the seismogenic zone. Resonant slip should initiate within the rupture zone of future great earthquakes, suggesting that slow slip may illuminate fault properties that control earthquake slip.     

Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array

The disastrous Sumatra-Andaman earthquake of 26 December 2004 was one of the largest ever recorded. The damage potential of such earthquakes depends on the extent and magnitude of fault slip. The first reliable moment magnitude estimate of 9.0 was obtained several hours after the Sumatra-Andaman earthquake, but more recent, longer-period, normal-mode analyses have indicated that it had a moment magnitude of 9.3, about 2.5 times larger. Here we introduce a method for directly imaging earthquake rupture that uses the first-arriving compressional wave and is potentially able to produce detailed images within 30 min of rupture initiation. We used the Hi-Net seismic array in Japan as an antenna to map the progression of slip by monitoring the direction of high-frequency radiation. We find that the rupture spread over the entire 1,300-km-long aftershock zone by propagating northward at roughly 2.8 km s(-1) for approximately 8 minutes. Comparisons with the aftershock areas of other great earthquakes indicate that the Sumatra-Andaman earthquake did indeed have a moment magnitude of approximately 9.3. Its rupture, in both duration and extent, is the longest ever recorded.     

Reflection signature of seismic and aseismic slip on the northern Cascadia subduction interface

At the northern Cascadia margin, the Juan de Fuca plate is underthrusting North America at about 45 mm x yr(-1) (ref. 1), resulting in the potential for destructive great earthquakes. The downdip extent of coupling between the two plates is difficult to determine because the most recent such earthquake (thought to have been in 1700) occurred before instrumental recording. Thermal and deformation studies indicate that, off southern Vancouver Island, the interplate interface is presently fully locked for a distance of approximately 60 km downdip from the deformation front. Great thrust earthquakes on this section of the interface (with magnitudes of up to 9) have been estimated to occur at an average interval of about 590 yr (ref. 3). Further downdip there is a transition from fully locked behaviour to aseismic sliding (where high temperatures allow ductile deformation), with the deep aseismic zone exhibiting slow-slip thrust events. Here we show that there is a change in the reflection character on seismic images from a thin sharp reflection where the subduction thrust is inferred to be locked, to a broad reflection band at greater depth where aseismic slip is thought to be occurring. This change in reflection character may provide a new technique to map the landward extent of rupture in great earthquakes and improve the characterization of seismic hazards in subduction zones.     

Near-field surface movement during the Wenchuan M s8.0 earthquake measured by high-rate GPS

Based on high-rate (1 Hz) GPS data from the Sichuan GPS Continuous Observation Network on the footwall of the Longmenshan Fault, we have characterized the near-field surface movement process during the 5.12 Wenchuan Earthquake. Results show that the maximum deformation near the fault is larger than the deformation set. Stations on the northern segment of the fault moved towards the epicenter first, and then turned toward the vertical orientation of the fault. Deformation of stations on the southern segment is smaller and recovered. The initial motion at all of the stations was downward followed by periodic up and down movements. Comparing the displacement from high-rate GPS and accelerograph data, we can see that they are consistent prior to the arrival of the principal shock, but afterwards a 10 cm difference is found, even though they are synchronized and in phase. More work is yet to be done to explain this. This is the first time that actual the real near-field surface movements of an earthquake of M >7.0 have been determined in china. These measurements are therefore of high value for studies of surface rupture processes and the analysis of seismic wave travels paths in this region.     

Rupture process of the Wenchuan <Emphasis Type="Italic">M</Emphasis>8.0 earthquake of Sichuan China: the segmentation feature

Moment tensor solution, rupture process and rupture characteristics of the great Wenchuan M8.0 earthquake are studied by using 39 long-period P and SH waveforms with evenly azimuth coverage of stations. Our results reveal that the Wenchuan M8.0 event consisted of 5 sub-events of Mw≥7.3 occurring succesively in time and space. Rupture started with a Mw7.3 introductory strike-slip faulting in the first 12 s, then within 12?40 s, two sub-events with Mw7.6 and Mw7.4 occurred within 80 km northeast from the init...     

Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit

Our understanding of the earthquake process requires detailed insights into how the tectonic stresses are accumulated and released on seismogenic faults. We derive the full vector displacement field due to the Bam, Iran, earthquake of moment magnitude 6.5 using radar data from the Envisat satellite of the European Space Agency. Analysis of surface deformation indicates that most of the seismic moment release along the 20-km-long strike-slip rupture occurred at a shallow depth of 4-5 km, yet the rupture did not break the surface. The Bam event may therefore represent an end-member case of the 'shallow slip deficit' model, which postulates that coseismic slip in the uppermost crust is systematically less than that at seismogenic depths (4-10 km). The InSAR-derived surface displacement data from the Bam and other large shallow earthquakes suggest that the uppermost section of the seismogenic crust around young and developing faults may undergo a distributed failure in the interseismic period, thereby accumulating little elastic strain.