Clues from joint inversion of tsunami and geodetic data of the 2011 Tohoku-oki earthquake
http://www.nature.com/srep/2012/120427/srep00385/full/srep00385.html#/supplementary-information
F. Romano,
A. Piatanesi,
S. Lorito,
N. D'Agostino,
K. Hirata,
S. Atzori,
Y. Yamazaki
& M. Cocco
The 2011 Tohoku-oki (Mw 9.1) earthquake is so far the best-observed megathrust rupture, which allowed the collection of unprecedented offshore data. The joint inversion of tsunami waveforms (DART buoys, bottom pressure sensors, coastal wave gauges, and GPS-buoys) and static geodetic data (onshore GPS, seafloor displacements obtained by a GPS/acoustic combination technique), allows us to retrieve the slip distribution on a non-planar fault. We show that the inclusion of near-source data is necessary to image the details of slip pattern (maximum slip ~48 m, up to ~35 m close to the Japan trench), which generated the large and shallow seafloor coseismic deformations and the devastating inundation of the Japanese coast. We investigate the relation between the spatial distribution of previously inferred interseismic coupling and coseismic slip and we highlight the importance of seafloor geodetic measurements to constrain the interseismic coupling, which is one of the key-elements for long-term earthquake and tsunami hazard assessment.
A giant interplate earthquake of magnitude Mw 9.1 occurred on March 11, 2011 at 05:46 UTC (Japan Meteorological Agency, JMA) off the Tohoku coast (Japan) with epicentre at 142.68°E 38.19°N1, at the interface between the Pacific and Okhotsk plates (Fig. 1a). It generated a huge tsunami that devastated the eastern coasts of the Honshu Island, and locally inundated up to 5 km inland killing more than 15,000 people2. It also caused a severe nuclear accident at the Fukushima nuclear power plant, whose effects are still under investigation.
Figure 1: Location map of the 11 March 2011 Mw 9.1 Tohoku-oki earthquake.
a) Red star indicates the epicentre position1. Red and white “beach ball” represents the focal mechanism of this earthquake1. Red triangles indicate the DART stations used in the inversion; b) Yellow, pink and cyan shadow zones are approximated rupture areas of the 869 Jogan, 1896 and 1933 Sanriku-oki earthquakes respectively39. Cyan circles indicate GPS stations onshore, magenta circles the geodetic seafloor observation sites15, 16, red triangles the bottom pressure sensors and GPS-buoys (Table S3 in Supplementary Information). White arrow indicates the approximate convergence direction of the Pacific plate (estimated velocity of 9.2 cm/yr54).
Seismological, geodetic and marine observational networks provided unprecedented high-quality recordings of the 2011 Tohoku-oki earthquake. More important, this is the first megathrust event that allowed the collection of offshore near-field data recorded very close to the causative source and to the trench. Tsunami waves were recorded by coastal wave gauges3, GPS-buoys4 and bottom pressure sensors5, 6, 7, 8 close to the Japanese coast off Iwate, Miyagi and Fukushima Prefectures2, 9, 10, 11, 12, 13 and by Deep-ocean Assessment and Report Tsunami (DART) buoys14 in the open Pacific Ocean, whereas several seafloor instruments measured the coseismic deformation in the shallowest part of the Okhotsk plate15, 16; up to 31 m of horizontal coseismic displacement were measured in the epicentral zone by a GPS/acoustic combination technique15, 16.
Tsunami waveforms and onshore geodetic data (e.g. GPS, InSAR) are often used to retrieve the slip distribution of great subduction earthquakes17, 18, 19. Tsunami data improve the resolution on the offshore portion of the subduction interface, whereas onshore geodetic data well resolve the slip distribution on the landward portion of the fault plane. Their combined use thus helps to constrain the slip on the entire fault20, 21. In the case of the 2011 Tohoku-oki event the presence of near-source offshore instruments yields a more robust picture of the rupture2, 9, 10, 11, 12, 13, 22, 23, 24, with important implications both for long-term seismic and tsunami hazard assessment and tsunami warning11, 25.
The relation between interseismic coupling (plate convergence fraction that is not accommodated by aseismic sliding) and coseismic slip has been recently studied by many authors and for different great subduction earthquakes19, 26, 27, 28, 29, 30, 31, 32. Zones of possible future large earthquakes are those where stress accumulates during the interseismic period26, 27, 33 and they can be likely identified by mapping the degree of plate coupling. Highly coupled zones will eventually rupture in future earthquakes34 making the estimation of plate coupling a key ingredient in long-term seismic and tsunami hazard assessment. However, even assuming that the imaged coupling remains constant in time, forecasting which parts of the fault will rupture during the next seismic event is still an open issue, since this depends on the fault zone mechanical conditions determined by local heterogeneous pre-stress and frictional properties controlling the dynamic rupture evolution. For example, ruptures in large earthquakes are likely to nucleate in coupled zones but they may propagate also in relatively uncoupled portions of the fault plane extending all the way to the trench35. Large historical tsunamigenic events already occurred in the past off Tohoku (Fig.1b), such as the megathrust earthquakes in 869 (Jogan36, 37) and in 1896 (Sanriku38), and an outer-rise event in 1933 (Sanriku39, 40). Some authors argue that the 1896 Sanriku earthquake has been characterized by very shallow rupture with destructive effects locally similar to 2011 event. Most of the models of the 2011 Tohoku-oki earthquake slip distribution obtained up to now inverting tsunami, seismic, and geodetic data agree in imaging the rupture as characterized by extremely high slip values concentrated in a relatively small area, with the rupture extending up to the trench for most slip models2, 9, 10, 12, 13, 22, 23, 24, 31, 41, 42, 43, 44, 45, 46, 47, 48, 49.
Models of interseismic coupling obtained inverting onshore geodetic measurements of the interseismic deformation before the 2011 earthquake26, 27, 50 showed a large locked portion of the subduction zone offshore the Tohoku region. These models then could have suggested the occurrence of a giant megathrust event, even if they feature a relatively low coupling in the shallowest portion of the subduction zone near the trench. This last result can be due to the limited resolution of onshore geodetic data on the shallowest portion of the plate interface near the trench31. This bias can be removed by using data recorded in the very near-field by seafloor instruments that revealed to be very effective to measure the plate motion offshore13, 22, 23, 24, 51, 52. In other words, the amount of strain accumulation on the shallow part of the megathrust might be better imaged by means of seafloor near-trench observations.
In this study, we perform a joint inversion using tsunami waveforms (DART buoys, bottom pressure sensors, coastal wave gauges and GPS-buoys, Fig. 1) and geodetic data (GPS onshore and seafloor displacements obtained by a GPS/acoustic combination technique15, 16, Fig. 1b), to infer the slip distribution of the 2011 Tohoku-oki earthquake. We use realistic fault geometry with variable strike and dip and extend our fault model up to the trench (Table S1 in Supplementary Information). We then discuss the results in the context of long-term earthquake and tsunami hazard assessment along subduction zones, highlighting the importance of very near-source data25.
Results
We solve the inverse problem using a global search technique (the simulated annealing53) following previous studies that used this approach for inverting different geophysical data sets19, 21. The resolving power of different combinations of geodetic and tsunami data used here is evaluated with several synthetic checkerboard tests (Fig. S3 in Supplementary Information). The main result emerging from our study is that the use of seafloor observations significantly increases the control on the offshore slip pattern (Figs. S3b,c in Supplementary Information), and that the joint inversion of the whole data set provides a sufficiently homogeneous resolution on the entire fault plane (Fig. S3e in Supplementary Information). The solution of the inverse problem is inherently non-unique and therefore, instead of showing the best fitting model, we compute and show the average slip model calculated over the ensemble of models that satisfactorily fits the data19, 21. This will be done for each synthetic test as well as for each of the inversions for the 2011 Tohoku-oki earthquake presented below. The reported model errors (Fig. S6, Tables S2 in Supplementary Information for the slip distribution of the Tohoku-oki earthquake) are the standard deviations of the marginal distributions of the model parameter (slip value).
2011 Tohoku-oki slip distribution
The joint inversion of the whole geodetic and tsunami data sets shows that the rupture area with slip > 10 m extends over a region of about 250×180 km from about 36.5°N to 39.5°N and it is characterized by a main patch of slip (Fig. 2) located up-dip from the hypocenter (~143°E, 38°N), that features a peak slip of about 48 m at a depth of ~13 km (Fig. 2, Table S2 in Supplementary Information). The rupture extends eastward up to the trench with slip values of about 35 m and westward approximately near the Fukushima and Ibaraki regions at a depth of 25 km with slip values close to 10 m (Fig. 2, Table S2 in Supplementary Information).
Figure 2: 2011 Tohoku-oki slip model.
2011 Tohoku-oki slip distribution
The joint inversion of the whole geodetic and tsunami data sets shows that the rupture area with slip > 10 m extends over a region of about 250×180 km from about 36.5°N to 39.5°N and it is characterized by a main patch of slip (Fig. 2) located up-dip from the hypocenter (~143°E, 38°N), that features a peak slip of about 48 m at a depth of ~13 km (Fig. 2, Table S2 in Supplementary Information). The rupture extends eastward up to the trench with slip values of about 35 m and westward approximately near the Fukushima and Ibaraki regions at a depth of 25 km with slip values close to 10 m (Fig. 2, Table S2 in Supplementary Information).
Figure 2: 2011 Tohoku-oki slip model.
Slip distribution for the 2011 Tohoku-oki earthquake obtained from the joint inversion of tsunami and geodetic data. Orange arrows represent the slip direction (rake, Table S2 in Supplementary Information). Thin black contours above the fault plane indicate the interseismic coupling (from 10% to 100%, at 10% intervals) along the megathrust50. Black arrow indicates the approximate convergence direction of the Pacific plate (estimated velocity of 9.2 cm/yr54). Red star as of Figure 1.
The predicted tsunami waveforms show a generally good agreement with the observed records (Fig. 3c). Some discrepancies might suggest a more complex rupture history24, 41, 42, 43, 45, 47, 48, 49, 56, 58. In this study, we indeed assume a circular rupture front propagating with a velocity of 1 km/s relying on seismic estimates48, 59 because several performed tests (not shown) demonstrated that the data set used here does not allow a robust retrieval of rupture velocity. The overall slip pattern is consistent with the observations indicating a tsunami with two distinct phases. The two pressure gauges TM1-2 and most of the GPS-buoys recorded, after a flatter initiation likely due to rupture of a relatively deep patch12, a shorter period tsunami wave, with a maximum amplitude up to 5–6 meters (Fig. 3c), a dramatic anticipation of the devastating tsunami effects observed along Iwate, Miyagi and Fukushima Prefectures11. The very large slip (up to ~35 m) at shallow depths featured by the model appears sufficient to explain the impulsive character of the Tohoku tsunami12, though with a second order lack of peak amplitude (10% on the average), e.g. at TM1 and TM2 or at DART 21418. The horizontal coseismic displacements at GPS stations are well reproduced, whereas there is some local discrepancy between predicted and observed vertical displacements, which could be due to inaccuracies in the fault position or dip and to the homogeneous half space approximation, as rigidity contrasts may be important for modelling GPS60.
Figure 3: Comparison between observed and predicted data sets.
Comparison between observed (black) and predicted a) horizontal and b) vertical displacements at GPS (red) and geodetic seafloor observation sites (magenta); c) Comparison between the observed (black) and predicted (red) tsunami waveforms; peak value (predicted/observed) in meters is indicated for each station.
The large horizontal displacement (up to 31 m) seen in the offshore geodetic observations used in this work (Figs. 3a,b), including those relatively close to the trench, are predicted satisfactorily. However, even closer to the trench, there are two further sea bottom displacement measurements61 (Fig. S2 in Supplementary Information), indicating a greater horizontal displacement (up to ~70 m), which our model would largely underestimate. We do not attempt to include these data in the inversion for several reasons. We use relatively large subfaults in the elastic homogeneous half space approximation. Then, we are not able to model the inherently 3D fine scale structure of the trench region, with additional complexities such as rigidity contrasts, possible inelastic processes13, the eventual dislocation on a normal branching fault61, 62 with a likely consequent extensional regime of the sedimentary wedge, or the eventual contribution of the submarine mass failure evidenced by recent bathymetric surveys63, 64. We thus consider our model an accurate and robust long wavelength image of the slip on most of the fault, while it will need a refinement in a narrow band in the trench region, even with the aim to assess the relative importance of additional tsunamigenic mechanisms.
Interseismic coupling and coseismic slip during the Tohoku-oki earthquake
Our model (Fig. 2) suggests that the nucleation of the Tohoku-oki earthquake occurred within a strongly coupled zone (between 70 and 80 %). The overall rupture zone conversely has a limited overlap with the coupling itself, consistently with the results of other studies31. Slip distribution is more elongated in the SW-NE direction and its centroid is displaced toward the trench, featuring up to ~35 meters of slip from about 38°N to 39°N in a zone of relatively low to zero coupling (< 30%). Whether this rupture propagated in a relatively uncoupled zone or, alternatively, the coupling is poorly resolved and higher than expected close to the trench is a matter hard to address. It has been suggested that a classical elastic back-slip accumulation model may not be sufficient to explain the shallow rupture propagation for great subduction earthquakes65. However, the onshore GPS data cannot fully resolve the coupling distribution in the trench region31. This is also shown by our resolution tests (Fig. S3b in Supplementary Information), which further demonstrate that using the seafloor instruments increases the resolution of the offshore slip pattern (Figs. S3b,c in Supplementary Information). We therefore extend this result to the interseismic coupling and we suggest that, using offshore geodetic measurements collected during the interseismic period, would allow us to better constrain the spatial distribution of coupling25.This can be further investigated by inverting the different geodetic data sets of the Tohoku-oki earthquake. The slip model resulting from the inversion of only the onshore GPS data features a main slip patch concentrated around the hypocenter, and a maximum slip amplitude of ~25 meters at a depth of ~25 km (Fig. 4a), consistently to other GPS inversions66, 67. This slip model is also correlated to the coupling distribution. However, whereas the horizontal displacements at the GPS onshore stations are very well reproduced (Figs. S4a,b in Supplementary Information), the model is not suitable to explain the large coseismic deformations measured offshore15, 16 (Figs. S4a,b in Supplementary Information), and the observed tsunami is dramatically underestimated (Fig. S4c in Supplementary Information). On the other hand, inverting also seafloor coseismic displacements obtained by a GPS/acoustic combination technique15, 16 we obtain a shallower (about 15 km of depth) primary patch of slip (Fig. 4b) with larger slip amplitude (~46 meters) at about 38°N, and both the inland and offshore horizontal geodetic displacements are well matched (Figs. S5a,b in Supplementary Information). Therefore near-source seafloor instruments are essential to retrieve information about the slip offshore and GPS on land alone are not able to resolve. Even so, the impulsive tsunami component in the near-source (GPS-buoys and TMs) is underestimated (Fig. S5c in Supplementary Information). The present configuration of geodetic data has a sufficient control of the macroscopic (long-wavelength) offshore slip pattern, but can predict only the far-field tsunami (Fig. S5c in Supplementary Information), since the details of coseismic slip pattern are not retrieved. It thus appears also evident the primary role of offshore tsunami data (TMs and GPS buoys) to constrain some features of the slip distribution (Figs. 2,3) which cannot be well resolved by a few offshore geodetic data. This result is also consistent with those of the resolution tests in Supplementary Information (Figs. S3c,d,e). Near source tsunameters are also very important for tsunami warning as their records can be accessible in real-time, more easily than seafloor geodetic measurements, supporting their planned intensive deployment and their integration into the existing Japanese tsunami warning system11, 68.
Figure 4: Geodetic Slip models.
a) Slip distribution for the 2011 Tohoku-oki earthquake obtained from the inversion of GPS onshore only. Cyan circles within the upper left inset indicate GPS stations onshore used in the inversion; b) Slip distribution for the 2011 Tohoku-oki earthquake obtained by inverting the whole geodetic data set. Cyan circles within the upper left inset indicate GPS stations onshore, magenta circles the geodetic seafloor observation sites, both used in the inversion. Thin black contours and red star as of Figure 2.
The increase of near-trench resolution determined by offshore geodetic data (as shown in this study for the slip distribution), would extend to the estimation of interseismic coupling if offshore geodetic instruments were used for estimating displacement rates during the interseismic period, and of course deployment of additional instruments would further increase the accuracy of coupling estimation. Actually, some previous studies in Japan69, 70, though based only on single-point measurements, already showed the potentiality of offshore geodetic instruments. In particular, these coupling estimations are fairly consistent along strike direction with recent coupling distributions50, showing a decrease of coupling departing from Miyagi toward Fukushima. Conversely, coupling off Miyagi is found to increase in the updip direction at least up to a depth of about 20 km on the fault plane and in contrast with the deeper maximum coupling estimated by onshore measurements. This is analogous to our finding that the centroid of the slip distribution moves updip when using also offshore data (Fig. 4). However, getting a comprehensive image of the secular seafloor movement would require a significant improvement of observational capability and continuous measurements52 which would in turn improve our understanding of the physical processes generating megathrust earthquakes and associated tsunamis, as well as our ability to define tsunami scenarios for hazard assessment.
Discussion
The slip distribution of the 2011 Tohoku-oki earthquake is characterized by extremely high slip values (up to ~48 meters) concentrated, despite of the great magnitude (Mw 9.1), in a relatively small area, with the rupture extending with high slip values (~35 meters) up to the trench.
Our findings, resulting from a joint inversion of geodetic and tsunami data, account for the huge tsunami ensuing from the earthquake, and most of the very large and shallow coseismic deformation associated to this seismic event. A possible further contribution from the frontal wedge deformation or secondary faults will be further investigated using additional data and by means of more sophisticated models (e.g. finite elements models) in order to take into account the three-dimensional geometrical and rheological complexities.
We have analyzed the different resolving power of geodetic and tsunami data, stressing the importance to use offshore instruments in the very near-field of the source area to recover a robust picture of the rupture, which will yield important implications for long-term seismic and tsunami hazard assessment. Seafloor instruments, measuring the deformations along the megathrust, should be used to better estimate the interseismic plate coupling, which is a key element to infer the tectonic strain accumulation. Our results strongly support the implementation of offshore monitoring networks along subduction zones, along with near-source tsunameters which would significantly improve existing tsunami early warning systems.
© 2012 Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.partner of AGORA, HINARI, OARE, INASP, ORCID, CrossRef and COUNTER
