Satellite Data for a Rapid Assessment of Tsunami Inundation Areas after the 2011 Tohoku Tsunami

The M _w = 9.0 earthquake that occurred off the coast of Japan’s Tohoku region produced a great tsunami causing catastrophic damage and loss of life. Within hours of the tsunami event, satellite data were readily available and massive media coverage immediately circulated thousands of photographs and videos of the tsunami. Satellite data allow a rapid assessment of inundated areas where access can be difficult either as a result of damaged infrastructure (e.g., roads, bridges, ports, airports) or because of safety issues (e.g., the hazard at Nuclear Power Plant at Fukushima). In this study, we assessed in a day tsunami inundation distances and runup heights using satellite data (very high-resolution satellite images from the GeoEye1 satellite and from the DigitalGlobe worldview, SRTM and ASTER GDEM) of the Tohoku region, Northeast Japan. Field survey data by Japanese and other international scientists validated our results. This study focused on three different locations. Site selection was based on coastal morphologies and the distance to the tsunami source (epicenter). Study sites are Rikuzentakata, Oyagawahama, and Yagawahama in the Oshika Peninsula, and the Sendai coastal plain (Sendai City to Yamamoto City). Maximum inundation distance (6 km along the river) and maximum runup (39 m) at Rikuzentakata estimated from satellite data agree closely with the 39.7 m inundation reported in the field. Here the ria coastal morphology and horn shaped bay enhanced the tsunami runup and effects. The Sendai coastal plain shows large inundation distances (6 km) and lower runup heights. Natori City and Wakabayashi Ward, on the Sendai plain, have similar runup values (12 and 16 m, respectively) obtained from SRTM data; these are comparable to those obtained from field surveys (12 and 9.5 m). However, at Yagawahama and Oyagawahama, Miyagi Prefecture, both SRTM and ASTER data provided maximum runup heights (41 to 45 m and 33 to 34 m, respectively), which are higher than those measured in the field (about 27 m). This difference in DEM and field data is associated with ASTER and SRTM DEM’s pixel size and vertical accuracy, the latter being dependent on ground coverage, slope, aspect and elevation. Countries with less access to technology and infrastructure can benefit from the use of satellite imagery and freely available DEMs for an initial, pre-field surveys, rapid estimate of inundated areas, distances and runup, and for assisting in hazard management and mitigation after a natural disaster.     

Modeling of the 2011 Japan Tsunami: Lessons for Near-Field Forecast

During the devastating 11 March 2011 Japanese tsunami, data from two tsunami detectors were used to determine the tsunami source within 1.5 h of earthquake origin time. For the first time, multiple near-field tsunami measurements of the 2011 Japanese tsunami were used to demonstrate the accuracy of the National Oceanic and Atmospheric Administration (NOAA) real-time flooding forecast system in the far field. To test the accuracy of the same forecast system in the near field, a total of 11 numerical models with grids telescoped to 2 arcsec (~60 m) were developed to hindcast the propagation and coastal inundation of the 2011 Japanese tsunami along the entire east coastline of Japan. Using the NOAA tsunami source computed in near real-time, the model results of tsunami propagation are validated with tsunami time series measured at different water depths offshore and near shore along Japan’s coastline. The computed tsunami runup height and spatial distribution are highly consistent with post-tsunami survey data collected along the Japanese coastline. The computed inundation penetration also agrees well with survey data, giving a modeling accuracy of 85.5 % for the inundation areas along 800 km of coastline between Ibaraki Prefecture (north of Kashima) and Aomori Prefecture (south of Rokkasho). The inundation model results highlighted the variability of tsunami impact in response to different offshore bathymetry and flooded terrain. Comparison of tsunami sources inferred from different indirect methods shows the crucial importance of deep-ocean tsunami measurements for real-time tsunami forecasts. The agreement between model results and observations along Japan’s coastline demonstrate the ability and potential of NOAA’s methodology for real-time near-field tsunami flooding forecasts. An accurate tsunami flooding forecast within 30 min may now be possible using the NOAA forecast methodology with carefully placed tsunameters and large-scale high-resolution inundation models with powerful computing capabilities.     

Japan's 2011 Tsunami: Characteristics of Wave Propagation from Observations and Numerical Modelling

We use a numerical tsunami model to describe wave energy decay and transformation in the Pacific Ocean during the 2011 Tohoku tsunami. The numerical model was initialised with the results from a seismological finite fault model and validated using deep-ocean bottom pressure records from DARTs, from the NEPTUNE-Canada cabled observatory, as well as data from four satellite altimetry passes. We used statistical analysis of the available observations collected during the Japan 2011 tsunami and of the corresponding numerical model to demonstrate that the temporal evolution of tsunami wave energy in the Pacific Ocean leads to the wave energy equipartition law. Similar equipartition laws are well known for wave multi-scattering processes in seismology, electromagnetism and acoustics. We also show that the long-term near-equilibrium state is governed by this law: after the passage of the tsunami front, the tsunami wave energy density tends to be inversely proportional to the water depth. This fact leads to a definition of tsunami wave intensity that is simply energy density times the depth. This wave intensity fills the Pacific Ocean basin uniformly, except for the areas of energy sinks in the Southern Ocean and Bering Sea.     

Manifestation of the 2011 Great Tohoku Tsunami on the Coast of the Kuril Islands: A Tsunami with Ice

The tsunami generated by the 2011 Tohoku Earthquake (M _w = 9.0) reached maximum heights of about 5 m along the coast of the Kuril Islands. The most essential feature of this event was sea ice about 0.5 m thick moved by the ocean water. The tsunami did not cause any essential damage on the Kuril Islands, but significantly affected coastal zones and produced interesting effects. The problem of a tsunami accompanied by marine ice is discussed and illustrated with photos.     

利用海啸数值模拟结果进行海底地震有限断层模型验证

基于地震有限断层模型进行海啸的数值模拟通常被用来估计海啸的到时、 波高等, 另一方面, 海啸数值模拟的结果也可以作为限定条件用来考察同一地震的不同断层模型之间的相对合理性。 采用国际上各地震研究机构在震后各自得出的不同的有限断层模型作为海啸源, 使用基于二维浅水波方程的海啸传播模型对2011年日本东北地震海啸的传播过程进行模拟, 以海啸模拟所得到的沿岸浪高分布、 平均波高、 最大波高等与实际观测值相比较, 进而判断由各有限断层模型所计算得到的海啸中哪个结果与实际的海啸情况更为符合, 由此推断断层模型的相对合理性。     

W Phase Inversion and Tsunami Inundation Modeling for Tsunami Early Warning: Case Study for the 2011 Tohoku Event

Centroid moment tensor solutions for the 2011 Tohoku earthquake are determined by W phase inversions using 5 and 10 min data recorded by the Full Range Seismograph Network of Japan (F-net). By a scaling relation of moment magnitude to rupture area and an assumption of rigidity of 4 × 10¹⁰ N m^?2, simple rectangular earthquake fault models are estimated from the solutions. Tsunami inundations in the Sendai Plain, Minamisanriku, Rikuzentakata, and Taro are simulated using the estimated fault models. Then the simulated tsunami inundation area and heights are compared with the observations. Even the simulated tsunami heights and inundations from the W phase solution that used only 5 min data are considerably similar to the observations. The results are improved when using 10 min of W phase data. These show that the W phase solutions are reliable to be used for tsunami inundation modeling. Furthermore, the technique that combines W phase inversion and tsunami inundation modeling can produce results that have sufficient accuracy for tsunami early warning purposes.     

Observations, Effects and Real Time Assessment of the March 11, 2011 Tohoku-oki Tsunami in New Zealand

The great Tohoku-oki earthquake of March 11, 2011 generated a devastating tsunami in the near field as well as substantial far-field effects throughout the Pacific Ocean. In New Zealand, the tsunami was widely observed and instrumentally recorded on an extensive array of coastal tidal gauges and supplemented by current velocity data from two sites. While the tsunami's first arrival was on the morning of March 12 in New Zealand, the strongest effects occurred throughout that afternoon and into the following day. Tsunami effects consisted primarily of rapid changes in water level and associated strong currents that affected numerous bays, harbors, tidal inlets and marine facilities, particularly on the northern and eastern shores of the North Island. The tsunami caused moderate damage and significant overland flooding at one location. The tsunami signal was clearly evident on tide gauge recordings for well over 2 days, clearly illustrating the extended duration of far field tsunami hazards. Real time analysis and modelling of the tsunami through the night of March 11, as the tsunami crossed the Pacific, was used as a basis for escalating the predicted threat level for the northern region of New Zealand. A comparison to recorded data following the tsunami shows that these real time prediction models were accurate despite the coarse near-shore bathymetry used in the assessment, suggesting the efficacy of such techniques for future events from far-field sources.     

利用PETSc和FEPG编制海啸数值模拟程序的研究

海啸的数值模拟是海啸研究的一个重要领域, 它对于帮助理解海啸的基本物理特性和预防减灾具有重要意义。 海啸数值模拟程序的编制是一项繁杂的工作, 该文介绍了利用两种功能强大的通用软件(PETSc和FEPG)来进行海啸数值模拟程序编写的研究。 PETSc和FEPG采用有限差分、 有限元、 有限体积等多种离散方法, 可以对数值问题给出稳定的求解。 该文采用了直角坐标系下的非线性浅水波方程作为海啸波传播的控制方程进行离散求解, 并将其结果与TUNAMI N1模型进行了比较, 表明这两种方法方便而且有效。     

The 2011 Tohoku-oki Earthquake Tsunami: Similarities and Differences to the 869 Jogan Tsunami on the Sendai Plain

A post-tsunami field survey following the 2011 Tohoku-oki Earthquake Tsunami was carried out to asses inundated area in Sendai Plain, Northeast Japan. The type of inundation was classified into two categories (major and minor) according to the amount of accumulated debris, garbage and sediment. Major and minor inundations were identified up to 4 and 5 km from the coastline, respectively. Many artificial geomorphological features, such as roadway embankments and canals, were believed to have affected the run-up process of the tsunami. The inundation area of the 2011 tsunami on the Sendai Plain is compared with that of the 869 Jogan tsunami, which was reconstructed using numerical modeling based on available historical and geological records. The inundation area of the 2011 Tohoku-oki tsunami is comparable to that of the 869 Jogan tsunami, although a direct comparison is difficult due to differences in geomorphological contexts between the paleo period and the present.     

Source Fault Model of the 1771 Yaeyama Tsunami, Southern Ryukyu Islands, Japan, Inferred from Numerical Simulation

The 1771 Yaeyama tsunami is successfully reproduced using a simple faulting model without submarine landslide. The Yaeyama tsunami (M 7.4), which struck the southern Ryukyu Islands of Japan, produced unusually high tsunami amplitudes on the southeastern coast of Ishigaki Island and caused significant damage, including 12,000 casualties. Previous tsunami source models for this event have included both seismological faults and submarine landslides. However, no evidence of landslides in the source has been obtained, despite marine surveying of the area. The seismological fault model proposed in this study, describing a fault to the east of Ishigaki Island, successfully reproduces the distribution of tsunami runup on the southern coast of the Ryukyu Islands. The unusual runup heights are found through the numerical simulation attributable to a concentration of tsunami energy toward the southeastern coast of Ishigaki Island by the effect of the shelf to the east. Thus, the unusual runup heights observed on the southeastern coast of Ishigaki Island can be adequately explained by a seismological fault model with wave-ray bending on the adjacent shelf.     

From 3-Hz P Waves to 0 S 2: No Evidence of A Slow Component to the Source of the 2011 Tohoku Earthquake

In the hours following the 2011 Honshu event, and as part of tsunami warning procedures at the Laboratoire de Géophysique in Papeete, Tahiti, the seismic source of the event was analyzed using a number of real-time procedures. The ultra-long period mantle magnitude algorithm suggests a static moment of 4.1 × 10²⁹ dyn cm, not significantly different from the National Earthquake Information Center (NEIC) value obtained by W-phase inversion. The slowness parameter, $\Uptheta = -5.65, $ is slightly deficient, but characteristic of other large subduction events such as Nias (2005) or Peru (2001); it remains significantly larger than for slow earthquakes such as Sumatra (2004) or Mentawai (2010). Similarly, the duration of high-frequency (2–4 Hz) P waves in relation to seismic moment or estimated energy, fails to document any slowness in the seismic source. These results were confirmed in the ensuing weeks by the analysis of the lowest-frequency spheroidal modes of the Earth. A dataset of 117 fits for eight modes (including the gravest one, ₀ S ₂, and the breathing mode, ₀ S ₀) yields a remarkably flat spectrum, with an average moment of 3.5 × 10²⁹ dyn cm (*/1.07). This behavior of the Tohoku earthquake explains the generally successful real-time modeling of its teleseismic tsunami, based on available seismic source scaling laws. On the other hand, it confirms the dichotomy, among mega-quakes (M ₀ > 10²⁹ dyn cm) between regular events (Nias, 2005; Chile, 2010; Sendai, 2011) and slow ones (Chile, 1960; Alaska, 1964; Sumatra, 2004; and probably Rat Island, 1965), whose origin remains unexplained.     

The 15 August 2007 Peru Earthquake and Tsunami: Influence of the Source Characteristics on the Tsunami Heights

The tsunami caused by the 2007 Peru earthquake (M_w 8.0) provoked less damage than by the seismic shaking itself (numerous casualties due to the earthquake in the vicinity of Pisco). However, it propagated across the Pacific Ocean and small waves were observed on one tide gauge in Taiohae Bay (Nuku Hiva, Marquesas, French Polynesia). We invert seismological data to recover the rupture pattern in two steps. The first step uses surface waves to find a solution for the moment tensor, and the second step uses body waves to compute the slip distribution in the source area. We find the slip distribution to consist of two main slip patches in the source area. The inversion of surface waves yields a scalar moment of 8.9 10²⁰ Nm, and body-wave inversion gives 1.4 10²¹ Nm. The inversion of tsunami data recorded on a single deep ocean sensor also can be used to compute a fault slip pattern (yielding a scalar moment of 1.1 10²¹ Nm). We then use these different sources to model the tsunami propagation across the Pacific Ocean, especially towards Nuku Hiva. While the source model taken from the body-wave inversion yields computed tsunami waves systematically too low with respect to observations (on the central Pacific Ocean DART buoy as on the Polynesian tide gauge), the source model established from the surface-wave inversion is more efficient to fit the observations, confirming that the tsunami is sensitive to the low frequency component of the source. Finally we also discuss the modeling of the late tsunami arrivals in Taiohae Bay using several friction coefficients for the sea bottom.     

Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami

— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented.     

基于应变能变化的芦山强震同震效应的数值模拟

芦山地震发生后,川滇地区主要活动断裂的地震活动趋势变化值得关注.本文基于川滇地区的三维非线性有限元模型,利用地震矩张量反演的断层节面性质及震级与破裂长度、同震位错量的经验公式,从应变能积累/释放角度分析地震的发生对川滇地区主要断裂的同震加卸载效应.初步结果表明,在33条断裂段中,有17条表现为同震加载效应,15条表现为卸载效应,1条影响不明显.地震的发生造成马尔康断裂、岷江断裂加载效应较强.但需要注意的是,滇西、滇西南地区主要活动断裂表现出群体性的加载效应.