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

See allHide authors and affiliations

Science  17 Jun 2011:
Vol. 332, Issue 6036, pp. 1421-1425
DOI: 10.1126/science.1206731


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.

The 2011 Tohoku-Oki earthquake occurred on the megathrust where the Pacific Plate subducts below Japan at an average rate of about 8 to 8.5 cm/year (Fig. 1) (1). Historically, many moment magnitude (Mw) 7 to Mw 8 earthquakes have occurred on the Japan Trench megathrust (2). Geodetic observations of crustal strain during the interseismic period have been used to infer spatial variations in the degree of plate coupling (i.e., regions of the megathrust expected to produce large earthquakes) for this section of the Japan Trench (3). Generally, these models infer high coupling in regions where earthquakes were known to have already occurred (Fig. 1 and fig. S1), with only partial or even no coupling from the trench to a point approximately midway between the trench and the coastline—precisely the region where the 2011 Tohoku-Oki earthquake occurred. It is fundamentally difficult to use land-based data to assess the state of coupling on distant portions of a megathrust. The Jogan earthquake of 13 July 869 may be the only documented event to have occurred with a possible magnitude and location similar to that of the 2011 earthquake (4).

Fig. 1

Map of central and northern Honshu, Japan. Vectors indicate the horizontal component of the GPS displacements for the mainshock (yellow) and the Mw 7.9 aftershock (orange). Approximate locations of historical megathrust earthquakes are indicated by closed curves colored by region: pink, Tokachi (1968 M8.2); green, Sanriku (1896 M8.5, 1901 M7.4, 1931 M7.6, 1933 M7.6); purple, Miyagi (1897 M7.4, 1936 M7.4, 1978 M7.4, 2005 M7.2); brown, Fukushima (1938 Mw 7.4, 1938 Mw 7.7, 1938 Mw 7.8) [modified from (2, 3, 33)]. Yellow and orange moment tensors indicate the W-phase centroid for the mainshock (34) and the GCMT (global centroid moment tensor) location for the M7.9 aftershock. The closed yellow curve indicates the outline of the Mw 9.0 mainshock (8 m slip contour). The region of inferred slip deficit or high plate coupling is indicated by dark blue nested contour lines for 35%, 70%, and 100% coupling (3). Barbed lines indicate subduction plate boundaries. The white arrow indicates the direction of convergence between the Pacific Plate and northeast Japan (1). T, S, and K indicate the cities of Tokyo, Sendai, and Kamaishi. The yellow box in the inset reference map shows the region of this figure and the locations of deep-sea bottom pressure gauges used in this study, all superimposed on the peak tsunami wave heights predicted by our preferred earthquake source model.

Observations of the 2011 Tohoku-Oki earthquake from a dense regional geodetic network and globally distributed broadband seismographic networks, as well as open-ocean tsunami data, allow the construction of a family of models that describe the distribution and evolution of subsurface fault slip. Surface displacements due to the Tohoku-Oki earthquake were observed by more than 1200 continuously recording Global Positioning System (GPS) sites installed and operated by the Geodetic Survey of Japan (GSI). Here, we used data sampled at 5-min intervals to produce individual three-component positional time series, from which we isolated coseismic displacements (Fig. 1) (5). Significant quasi-permanent displacements due to the mainshock occurred over the entire northern half of Honshu, with peak GPS-measured offsets exceeding 4.3 m horizontally and 66 cm of subsidence (Figs. 1 and 2). We also isolated surface displacements associated with an Mw 7.9 aftershock that occurred about 30 min after the mainshock (Fig. 1). The spatial extent and the azimuth of the horizontal displacement vectors indicate that the aftershock was located to the south of the mainshock in the Ibaraki segment. Peak horizontal GPS-measured displacements for this aftershock were approximately 44 cm (Fig. 1). We used observations of open-ocean tsunami wave heights measured by deep sea-bottom pressure gauges to constrain the distribution of coseismic slip on the shallowest portions of the megathrust (Fig. 1). On the basis of their spatial and azimuthal distribution, we selected 12 sensors in the Pacific ocean east of the Japan trench. The closest of these pressure gauges detected a maximum tsunami wave height of more than 1.9 m (Fig. 2A).

Fig. 2

(A) Observed (green) and predicted (white) deep ocean tsunami records of the Tohoku-Oki earthquake. The predicted records correspond to a model constructed using the mean of each fault slip parameter in the Bayesian inversion. These waveforms are superimposed on the map of maximum model-predicted tsunami height. (B) GPS vertical coseismic surface displacements (circles colored and scaled with amplitude) as well as model-predicted vertical sea-floor displacements (filled contours). Other overlay features are as in Fig. 1.

We first describe static coseismic slip models based on the GPS observations of coseismic offsets and the sea-floor pressure gauge data (Fig. 3 and fig. S2). Static models constrain the final distribution of slip for the event but not its temporal evolution. We adopt a novel fully Bayesian probabilistic formalism requiring no a priori spatial regularization (5, 6). We conservatively define the section of the megathrust directly involved with the earthquake by considering only the areas where inferred slip exceeds 8 m (approximately 15% of the maximum slip value, depending on the model). The model predicts maximum seafloor subsidence of about 2 m located 50 km offshore Sendai and Kamaishi, and maximum sea-floor uplift just under 9 m about 50 km from the trench (Fig. 2). This model fits the GPS and tsunami data with variance reductions of 99.7% and 90.1%, respectively. Residual GPS displacements after removal of the model predictions are shown in fig. S3.

Fig. 3

Inferred distribution of subsurface fault slip (color and black contours with a contour interval of 8 m). Fault slip associated with the Mw 7.9 aftershock is indicated by nested 1-m orange contours. Historic earthquake ellipses are as in Fig. 1. Location of points of HF radiation estimated using back-projection methods with data from the European Union seismic array and the USArray are indicated by squares and circles, respectively, with color intensity indicating time of the activity relative to the beginning of the event and with symbol size proportional to amplitude of the HF radiation normalized to the peak value. The star indicates the location of the Japan Meteorological Agency epicenter. See fig. S2 for a plot of the slip model without overlays.

The spatial distribution of slip (Fig. 3) can be divided into several sections. The central section contains the highest estimated slip values, with peak displacement of around 60 m. The up-dip limit of the forearc is an active accretionary prism that extends about 50 km landward of the trench. Generally, the majority of the fault slip does not extend below this zone, with the exception being just up-dip of the region of maximum fault slip where estimated slip values near the trench range from about 5 to 15 m. A tendril of slip extends more than 100 km north from the central slip zone and just down-dip from the inferred source of the 1896 Mw 8.0 Sanriku earthquake. Average slips in this region are approximately 5 to 10 m—similar to those inferred for the 1896 earthquake (7). A lobe of about 10 m of fault slip extends down-dip toward the Oshika Peninsula east-northeast of Sendai. This lobe overlaps several of the inferred historical Miyagi-Oki rupture areas. Slip in the up-dip portion of these rupture areas exceeds 20 m. Another tendril of significant slip (5 to 10 m) extends southward of the main high-slip asperity. This tendril clearly overlaps the inferred locations of the 1938 Fukushima earthquake sequence.

We estimated probability distributions for derived scalar rupture quantities including rupture area, potency, scalar seismic moment, and static stress drop (fig. S4). Estimates of moment magnitude range from 8.8 to 9.2. We note that static slip models are relatively insensitive to absolute scaling of the elastic moduli; thus, estimates of moment are less certain than estimates of potency. Estimates of static stress drop vary between 2 and 10 MPa, depending on the area of fault considered. These values are high relative to previous estimates for megathrust events, which typically lie in the 1 to 5 MPa range (8) and reflect the relatively small area over which there is high slip. As a point of comparison, relative to our model for the Tohoku-Oki event, models of the 2010 Mw 8.8 Maule, Chile earthquake typically find twice the along-strike extent of slip and half the peak slip (9, 10).

We also developed two kinematic finite fault models based on one and two fault planes, respectively, and constrained by broadband seismic data and GPS observations (but no tsunami data) (5, 11). Examples of the displacement and velocity waveform fits are shown in fig. S5. The inferred moment rate function suggests that most of the rupture occurred in a little over 3 min (fig. S6). We find a low average rupture velocity of about 1.2 km/s (fig. S6). This model explains 99% of the variance of the GPS data. The peak estimated fault slip in this model is about 45 m, a little less than found in the static model. This difference is due to the imposed smoothing in the kinematic model, which is absent in the static model. The extension of slip to the south (offshore Fukushima) is evident in the kinematic model; however, it is located up-dip relative to the static model, presumably because of the absence of the tsunami constraints in the kinematic model. The Mw 7.9 aftershock occurred just beyond the southernmost extent of the mainshock slip area with an estimated maximum slip of about 4 m for this event (Fig. 3) (5).

Observations at the high-frequency (HF) end of the seismic spectrum (2 to 4 Hz) can constrain rupture direction and duration (12); 90% of the energy release in this frequency band occurs within about 3 min (fig. S7). Unilateral rupture propagation would result in an azimuthally dependent duration. This earthquake displays uniform durations at most azimuths with slightly shortened durations in the down-dip direction, suggesting bilateral along-strike rupture with some down-dip propagation (fig. S7).

We developed an image of the rupture process at high frequencies (between 0.5 and 1 Hz) using back-projection of teleseismic array waveforms (13) based on high-resolution array processing techniques (5, 14, 15). The most energetic of the HF sources as a function of time during the rupture are systematically down-dip of the regions of largest fault slip (Fig. 3). The robustness of the locations of the HF radiators relative to the assumed hypocenter is supported by the consistency between the results obtained with USArray data and those from the European array (Fig. 3 and fig. S8).

HF radiation is usually assumed to be spatially correlated to seismic slip (16) or uniformly distributed over the fault (17). This assumption contrasts with the spatial complementarity between HF and low-frequency (LF) source properties observed for this earthquake. Such a relationship has been inferred for other earthquakes, although not systematically (18), and might reflect the general lack of correlation between HF and LF in ground motions (19, 20). Dynamic rupture models generate HF radiation mostly during sudden changes of rupture speed along sharp contrasts of fault rheology or geometry, or along remnant stress concentrations from previous earthquakes (21), which can define the boundaries of the slip area.

Instead of relying on characteristics of the rupture dynamics to explain the predominance of HF radiators down-dip of the region of significant coseismic slip, an alternative explanation is that they occur at the transition between brittle and ductile regions. Recent observations of medium-sized earthquakes down-dip from aseismic sections as well as of slow-slip and nonvolcanic tremor phenomena suggest the presence of frictional heterogeneities in the brittle-ductile transition regions of subduction megathrusts (22). Coseismic triggering of compact brittle asperities embedded in the ductile fault matrix could explain the relative locations of HF and lower-frequency slip and the apparent lack of HF radiators at shallow depths.

Similarly, we find that large aftershocks with thrust mechanisms are located outside the region of coseismic slip (fig. S2), as has been documented for previous events (23, 24). Interestingly, for this event, aftershocks are dominantly down-dip of the regions of major coseismic slip—consistent with the idea that fault slip at shallow depths is relatively slow because of higher fracture energy and thus radiates energy less efficiently than at greater depths (8). This latter interpretation is supported by a comparison of seismic excitation between 0.5 and 4 Hz from the 2010 Maule and 2011 Tohoku-Oki events (fig. S9). Although the value of peak slip inferred for the Tohoku-Oki earthquake was larger than for the Maule earthquake by a factor of about 2 to 3 and the area of appreciable slip for the Tohoku-Oki is approximately half that of the Maule earthquake, we find that the former produced much less HF radiation. We suggest that this difference is because the Maule earthquake rupture is on average much deeper than the Tohoku-Oki rupture (fig. S10). Thus, it may be that shallow ruptures generally produce large displacements with relatively weak HF excitation. Such an interpretation is also broadly consistent with the general observations that tsunami earthquakes rupture very slowly, with slip concentrated in the shallowest part of the megathrust (25).

Our inferred fault slip model suggests high static stress drop with large amounts of slip in a small region. An explanation for this behavior theorizes the existence of barriers that require much more stress accumulation than other regions before they rupture. Such barriers may pin the fault locally, limiting the amount of seismic slip occurring on neighboring areas that have lower thresholds for failure. When the strongest barrier finally ruptures, then surrounding areas can catch up. Subducted seamounts are the most obvious candidates for such barriers (26). Indeed, several seamounts are known to have subducted in this segment of the Japan Trench (27). The distribution of slip in the Tohoku-Oki earthquake suggests that small areas can have high effective yield stresses that serve to limit rupture propagation during some earthquakes but then eventually rupture with large slip during others.

The extent to which the 2011 earthquake was unexpected suggests that we should consider the potential for similar large events elsewhere on the Japan Trench megathrust. The secular interseismic velocity field for the Miyagi segment shows more than 3 cm/year of relative convergence across Honshu (Fig. 4). On average, faults in the interior and off the western coast of Honshu, are believed to account for between 1 and 2 cm/year (28)—leaving 1 to 2 cm/year associated with interseismic strain accumulation on the subduction interface. We adopt 1100 years as a representative time period because this corresponds to the last large event that is inferred to have occurred in this region (4). Thus, over this time period, we must still account for 11 to 22 m of relative motion across Honshu. Similarly, at about 8.5 cm/year of convergence, we must account for more than 90 m of fault slip on the megathrust.

Fig. 4

Secular interseismic surface deformation (blue vectors) as observed by the GEONET continuous GPS network using the GSI F2 solution. The upper left inset shows the decrease of this deformation field at the latitude of the Miyagi (purple profile) and Ibaraki (orange profile) segments, with distance measured along the profile. The coseismic slip distribution is indicated by the yellow contours at 8-m intervals. The inset uses vectors within 100 km of the profile location. The question mark indicates a region of possible high seismic hazard. Other features are as in Fig. 1.

Earthquake activity offshore of Miyagi (Fig. 1) has been suggested to be dominated by M7+ earthquakes recurring every 30 to 40 years (2). However, it has already been established that the historical events are not exact repeats of one another (29). Further, such M7+ events only produce 3 to 4 m of fault slip and 5 to 20 cm of surface displacement per event. In the 2011 mainshock, fault slip in the region of the historical M7+ events ranged from 5 to 25 m, which suggests that the concept of a characteristic subduction earthquake with approximately the same slip per event at a given location may be of limited use (30). That the 2011 event produced approximately 50 m of slip up-dip of the historical Miyagi M7+ events is roughly consistent with a 500- to 1000-year potential slip accumulation period. However, there is no basis on which to assume that the aforementioned interval of 1100 years is representative of the recurrence interval of great earthquakes in this segment—it could be shorter by a factor of 2 and still be consistent with the surface displacement budget and the peak slip inferred in this recent earthquake.

The only previously recorded large events offshore Fukushima and Ibaraki occurred as a sequence in 1938, which taken together correspond to about an Mw 8.1 event (31). At about 8.0 cm/year convergence, the 73 years since those earthquakes imply about 6 m of accumulated potential fault slip—surprisingly similar to the estimates of fault slip for this region during the 2011 event. However, this agreement is probably coincidental because there is no record of an equivalent set of earthquakes preceding the 1938 sequence. Similarly, the 2011 Mw 7.9 aftershock offshore of Ibaraki (Fig. 1) produced about 4 m of fault slip, implying a 50-year recurrence if these events are characteristic of this segment of the megathrust. There is no documentation of large (M8+) events prior to 1938 (31).

The 2 cm/year GPS-observed onshore convergence in this region (Fig. 4) implies that, in terms of the 1100-year budget, there is between 0 and 10 m of surface convergence that cannot be plausibly associated with faults farther to the west (28). The combination of the 2011 mainshock and the M7.9 aftershock only produced about 2 m of surface displacement.

There is no record of a large event up-dip of the 1938 Fukushima and Ibaraki sequence. Thus, the slip budget on the megathrust and the surface velocity data suggest that an earthquake similar to the 2011 event is possible offshore Ibaraki and Fukushima just south of the most recent event (Fig. 4). During such an event, the 1938 asperities and the M7.9 aftershock rupture area could experience much greater slip than has been documented for previous events, similar to what just occurred offshore Miyagi. However, if this region is in fact not strongly coupled, simple mechanical models (32) would predict high rates of postseismic afterslip—as is inferred to have occurred after several large recent earthquakes (24). Thus, it is essential to monitor this region to quantify the extent of any postseismic slip in order to further understand the long-term fault slip budget and associated seismic hazard.

Supporting Online Material

SOM Text

Figs. S1 to S12

Tables S1 and S2

References (3547)

References and Notes

  1. “Seismic Activity in Japan—Regional Perspectives on the Characteristics of Destructive Earthquakes (Excerpt)” (Earthquake Research Committee, Headquarters for Earthquake Research Promotion, Prime Minister’s Office, Tokyo, 1998).
  2. See supporting material on Science Online.
  3. S. E. Minson, M. Simons, J. L. Beck, paper presented at the AGU 2010 fall meeting, San Francisco, 13 to 17 December 2010, abstr. G12A-02.
  4. Z. Duputel et al., (2011).
  5. Acknowledgments: Supported in part by the Gordon and Betty Moore Foundation. M.S. and S.E.M. are supported by NSF grant CDI-0941374; J.P.A. and L.S.M. are supported by NSF grant EAR-1015704 and the Southern California Earthquake Center, which is funded by NSF Cooperative Agreement EAR-0106924 and USGS Cooperative Agreement 02HQAG0008. A portion of the research was carried out at the Jet Propulsion Laboratory (JPL), California Institute of Technology, under a contract with the National Aeronautics and Space Administration and funded through the internal Research and Technology Development program. We acknowledge the Geospatial Information Authority (GSI) of Japan for kindly providing all the GEONET RINEX data. Raw RINEX data are available directly from GSI. Processed GPS time series are provided through the ARIA project. We thank T. Ito for providing the interseismic velocity estimates used in Fig. 4. This paper is Caltech Tectonics Observatory contribution 165 and Caltech Seismological Laboratory contribution 10059.
View Abstract

Stay Connected to Science

Navigate This Article