The El Salvador and Philippines Tsunamis of August 2012: Insights from Sea Level Data Analysis and Numerical Modeling

We studied two tsunamis from 2012, one generated by the El Salvador earthquake of 27 August (Mw 7.3) and the other generated by the Philippines earthquake of 31 August (Mw 7.6), using sea level data analysis and numerical modeling. For the El Salvador tsunami, the largest wave height was observed in Baltra, Galapagos Islands (71.1 cm) located about 1,400 km away from the source. The tsunami governing periods were around 9 and 19 min. Numerical modeling indicated that most of the tsunami energy was directed towards the Galapagos Islands, explaining the relatively large wave height there. For the Philippines tsunami, the maximum wave height of 30.5 cm was observed at Kushimoto in Japan located about 2,700 km away from the source. The tsunami governing periods were around 8, 12 and 29 min. Numerical modeling showed that a significant part of the far-field tsunami energy was directed towards the southern coast of Japan. Fourier and wavelet analyses as well as numerical modeling suggested that the dominant period of the first wave at stations normal to the fault strike is related to the fault width, while the period of the first wave at stations in the direction of fault strike is representative of the fault length.     

A fault model of the 1946 Nankaido earthquake derived from tsunami data

A fault model of the 1946 Nankaido earthquake (M = 8.2) is determined by the use of tsunami records of Uwajima, Shimotsu and Hososhima which were located within or near the area of major coseismic crustal deformation. Synthetic tsunamis computed for various fault models are matched with the observed tsunamis to determine the fault parameters. A low-angle thrust model slightly revised from a previous model by Ando is consistent with the observed tsunamis. The duration of faulting is constrained as less than 10 min based upon the tsunami. The fault is divided into an eastern and a western segment corresponding to areas associated with and without aftershocks, respectively. The fault area and dislocation for the western segment are 150 × 70 km² and 6 m, and those for the eastern segment are 150 × 70 km² and 3 m, respectively. The total seismic moment is 4.7 × 10²⁸ dyn·cm, significantly smaller than that obtained from a geodetic model by Fitch and Scholz, but still larger than that of the seismic model by Kanamori. The discrepancy in seismic moment between the seismic and the present models (RAN2) could be interpreted in terms of a slow dislocation on the fault, but this interpretation does not match the seismic intensity distribution and damage pattern, and the slow-slip model for the Nankaido earthquake is rejected. The discrepancy between the two seismic moments is considered insignificant within error involved in data and modeling assumptions. If the revised geodetic model (RAN2) is modified, the seismic moment required to explain the observed tsunamis would be reduced further by ～30%. If we consider the uncertainties involved in the fault model of Kanamori and the fault-finiteness effect affecting the amplitude of seismic waves, the seismic moment required to interpret the seismic-wave data could be increased, possibly being more than twice that of Kanamori. Thus, the two seismic moments from the different data sets could be close to each other within allowable tolerance. This implies that the rise time of the Nankaido earthquake was short enough to generate short-period seismic waves from both the western and the eastern fault segments.     

The source of the 1722 Algarve earthquake: evidence from MCS and Tsunami data

The 27 December 1722 Algarve earthquake destroyed a large area in southern Portugal generating a local tsunami that inundated the shallow areas of Tavira. It is unclear whether its source was located onshore or offshore and, in any case, what was the tectonic source responsible for the event. We analyze available historical information concerning macroseismicity and the tsunami to discuss the most probable location of the source. We also review available seismotectonic knowledge of the offshore region close to the probable epicenter, selecting a set of four candidate sources. We simulate tsunamis produced by these candidate sources assuming that the sea bottom displacement is caused by a compressive dislocation over a rectangular fault, as given by the half-space homogeneous elastic approach, and we use numerical modeling to study wave propagation and run-up. We conclude that the 27 December 1722 Tavira earthquake and tsunami was probably generated offshore, close to 37°01′N, 7°49′W.     

Introduction to “Historical and Recent Catastrophic Tsunamis in the World: Volume II. Tsunamis from 1755 to 2010”

Eighteen papers on past and recent destructive tsunamis are included in Volume II of the PAGEOPH topical issue “Historical and Recent Catastrophic Tsunamis in the World.” Three papers discuss deep-sea (DART) and coastal tsunami observations, warning systems and risk management in the Pacific Ocean. Four papers examine the 1755 Lisbon, 1964 Alaska, 2003 Algeria, and 2011 Haiti tsunamis. Four more papers, as well as some papers in Volume I, report on various aspects of the 2010 Chile tsunami. Two papers present some results of field survey and modelling investigation of the 2010 Mentawai, Indonesia, tsunami. Three papers report on modelling efforts of tsunami generation by earthquake and landslide, and of tsunami propagation. Finally, two papers discuss hazard assessment using a probabilistic approach.     

Numerical simulation of the 1992 Flores tsunami: Interpretation of tsunami phenomena in northeastern Flores Island and damage at Babi Island

Numerical analysis of the 1992 Flores Island, Indonesia earthquake tsunami is carried out with the composite fault model consisting of two different slip values. Computed results show good agreement with the measured runup heights in the northeastern part of Flores Island, except for those in the southern shore of Hading Bay and at Riangkroko. The landslides in the southern part of Hading Bay could generate local tsunamis of more than 10 m. The circular-arc slip model proposed in this study for wave generation due to landslides shows better results than the subsidence model, It is, however, difficult to reproduce the tsunami runup height of 26.2 m at Riangkroko, which was extraordinarily high compared to other places. The wave propagation process on a sea bottom with a steep slope, as well as landslides, may be the cause of the amplification of tsunami at Riangkroko. The simulation model demonstrates that the reflected wave along the northeastern shore of Flores Island, accompanying a high hydraulic pressure, could be the main cause of severe damage in the southern coast of Babi Island.     

海啸波近岸共振响应的数值模拟及分析

应用有限差分方法求解非线性浅水长波方程,建立了海啸波产生和传播的二维数值模型;对太平洋地震引起的夏威夷群岛海啸波进行模拟,并将模拟结果和测站实测值进行比较,验证了模型的正确性。利用快速傅里叶变换对数值模拟结果进行谱分析,得到整个计算区域的能量谱分布,并给出了发生能量聚集的位置及相应的谱峰周期。讨论了海啸波传播到近岸时可能产生的共振响应现象,发现海啸波和近岸的共振响应不仅与近岸复杂地形有关,还与海啸波传播到近岸时的波浪入射方向有关。     

Study on the relation of river morphology and tsunami propagation in rivers

The relations of river morphology and tsunami propagation in rivers were studied at several rivers in the Tohoku region during The Great Chilean Tsunami of 2010 and The Great East Japan Tsunami of 2011. It was found that river mouth morphological features play an important role in the intrusion of low magnitude tsunamis in which the geological and geographical conditions are an important factor. Nevertheless, the effects of these features were not found in the case of an extreme tsunami wave. As the wave enters the river, the propagation depends on other factors. It was found that the intrusion distance correlates well to the riverbed slope. The measurements of water level and riverbed slope were analyzed to propose an empirical method for estimating the damping coefficient for the tsunami propagation in rivers based on the tsunami of 2011. The proposed empirical method was used to approximate the length of the tsunami intrusion into a river by assuming that the furthest distance is given for the ratio of local tsunami wave height to the tsunami wave height at the river entrance of 0.05 (5 %). The estimated intrusion length from the proposed method in this study shows a good comparison with measurement data.     

Some Aspects of Energy Balance and Tsunami Generation by Earthquakes and Landslides

— Tsunamis are generated by displacement or motion of large volumes of water. While there are several documented cases of tsunami generation by volcanic eruptions and landslides, most observed tsunamis are attributed to earthquakes. Kinematic models of tsunami generation by earthquakes — where specified fault size and slip determine seafloor and sea-surface vertical motion — quantitatively explain far-field tsunami wave records. On the other hand, submarine landslides in subduction zones and other tectonic settings can generate large tsunamis that are hazardous along near-source coasts. Furthermore, the ongoing exploration of the oceans has found evidence for large paleo-landslides in many places, not just subduction zones. Thus, we want to know the relative contribution of faulting and landslides to tsunami generation. For earthquakes, only a small fraction of the minimum earthquake energy (less than 1% for typical parameter choices for shallow underthrusting earthquakes) can be converted into tsunami wave energy; yet, this is enough energy to generate terrible tsunamis. For submarine landslides, tsunami wave generation and landslide motion interact in a dynamic coupling. The dynamic problem of a 2-D translational slider block on a constant-angle slope can be solved using a Green's function approach for the wave transients. The key result is that the largest waves are generated when the ratio of initial water depth above the block to downslope vertical drop of the block H ₀ /W sin δ is less than 1. The conversion factor of gravitational energy into tsunami wave energy varies from 0% for a slow-velocity slide in deep water, to about 50% for a fast-velocity slide in shallow water and a motion abruptly truncated. To compare maximum tsunami wave amplitudes in the source region, great earthquakes produce amplitudes of a few meters at a wavelength fixed by the fault width of 100 km or so. For submarine landslides, tsunami wave heights — as measured by b, block height — are small for most of the parameter regime. However, for low initial dynamic friction and values of H ₀ /W sin δ less than 1, tsunami wave heights in the downslope and upslope directions reach b and b/4, respectively.Wavelengths of these large waves scale with block width. For significant submarine slides, the value of b can range from meters up to the kilometer scale. Thus, the extreme case of efficient tsunami generation by landslides produces dramatic hazards scenarios.     

Slip Distribution and Seismic Moment of the 2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms and Coastal Geodetic Data

The slip distribution and seismic moment of the 2010 and 1960 Chilean earthquakes were estimated from tsunami and coastal geodetic data. These two earthquakes generated transoceanic tsunamis, and the waveforms were recorded around the Pacific Ocean. In addition, coseismic coastal uplift and subsidence were measured around the source areas. For the 27 February 2010 Maule earthquake, inversion of the tsunami waveforms recorded at nearby coastal tide gauge and Deep Ocean Assessment and Reporting of Tsunamis (DART) stations combined with coastal geodetic data suggest two asperities: a northern one beneath the coast of Constitucion and a southern one around the Arauco Peninsula. The total fault length is approximately 400 km with seismic moment of 1.7 × 10²² Nm (Mw 8.8). The offshore DART tsunami waveforms require fault slips beneath the coasts, but the exact locations are better estimated by coastal geodetic data. The 22 May 1960 earthquake produced very large, ~30 m, slip off Valdivia. Joint inversion of tsunami waveforms, at tide gauge stations in South America, with coastal geodetic and leveling data shows total fault length of ~800 km and seismic moment of 7.2 × 10²² Nm (Mw 9.2). The seismic moment estimated from tsunami or joint inversion is similar to previous estimates from geodetic data, but much smaller than the results from seismic data analysis.     

The July 1998 Papua New Guinea Earthquake: Mechanism and Quantification of Unusual Tsunami Generation

— The unusual tsunami generated by the July 17, 1998 Papua New Guinea earthquake was investigated on the basis of various geophysical observations, including seismological data, tsunami waveform records, and on-land and submarine surveys. The tsunami source models were constructed for seismological high-angle and low-angle faults, splay fault, and submarine slumps. Far-field and near-field tsunamis computed from these models were compared with the recorded waveforms in and around Japan and the measured heights along the coast around Sissano Lagoon, respectively. In order to reproduce the far-field tsunami waveforms, small sources such as splay fault or submarine slump alone were not enough, and a seismological fault model was required. Relocated aftershock distribution and observed coastal subsidence were preferable for the low-angle fault, but the low-angle fault alone could not reproduce the large near-field tsunamis. The low-angle fault with additional source, possibly a submarine slump, is the most likely source of the 1998 tsunami, although other possibilities cannot be excluded. Computations from different source models showed that the far-field tsunami amplitudes are proportional to the displaced water volume at the source, and the comparison with the observed tsunami amplitudes indicated that the displaced water volume at the 1998 tsunami source was ～0.6 km³. The near-filed tsunami heights, on the other hand, are determined by the potential energy of displaced water, and the comparison with the observed heights showed that the potential energy was ～2 × 10¹² J.     

Tide Gauge Observations of 2004–2007 Indian Ocean Tsunamis from Sri Lanka and Western Australia

Tide gauge data collected from Sri Lanka (three stations) and Western Australia (eleven stations) during the Indian Ocean tsunamis, which occurred in December 2004, March 2005, July 2006, and September 2007, and incorporated five tsunamis, were examined to determine tsunami behaviour during these events. During the December 2004 tsunami, maximum wave heights of 3.87 m and 1.75 m were recorded at Colombo (Sri Lanka) and Bunbury (Western Australia), respectively. The results indicated that although the relative magnitudes of the tsunamis varied, the tsunami behaviour at each station was similar. This was due to the effect of the local and regional topography. At all tide gauges, the spectral energy corresponding to periods between 20 and 85 minutes increased during the tsunami. The sea-level data obtained from the west and south coasts of Sri Lanka (Colombo and Kirinda) indicated the importance of wave reflections from the Maldives Island chain, which produced the maximum wave two to three hours after the arrival of the first wave. In contrast, Trincomalee on the east coast did not show evidence of a reflected wave. Similarly, along the west coast of Australia, the highest waves occurred 15 hours after the arrival of the first wave. Here, based on travel times, we postulated that the waves were reflected from the Mascarene Ridge and/or the Island of Madagascar. Reflected waves were not present in the 2006 tsunami, where the primary waves propagated away from topographic features. One of the main influences of the tsunami was to set up oscillations at the local resonance frequency. Because Sri Lanka and Western Australia have relatively straight coastlines, these oscillations were related to the fundamental period of the shelf oscillation. For Colombo, this corresponded to 75-minute period, whereas in Geraldton and Busselton (Australia), the four-hour period was most prominent; at Jurien Bay and Fremantle, the resonance period was 2.7 hours.     

A Revised Evaluation of Tsunami Hazards along the Chinese Coast in View of the Tohoku-Oki Earthquake

Japan’s 2011 Tohoku-Oki earthquake and the accompanying tsunami have reminded us of the potential tsunami hazards from the Manila and Ryukyu trenches to the South China and East China Seas. Statistics of historical seismic records from nearly the last 4 decades have shown that major earthquakes do not necessarily agree with the local Gutenberg-Richter relationship. The probability of a mega-earthquake may be higher than we have previously estimated. Furthermore, we noted that the percentages of tsunami-associated earthquakes are much higher in major events, and the earthquakes with magnitudes equal to or greater than 8.8 have all triggered tsunamis in the past approximately 100 years. We will emphasize the importance of a thorough study of possible tsunami scenarios for hazard mitigation. We focus on several hypothetical earthquake-induced tsunamis caused by M _w 8.8 events along the Manila and Ryukyu trenches. We carried out numerical simulations based on shallow-water equations (SWE) to predict the tsunami dynamics in the South China and East China Seas. By analyzing the computed results we found that the height of the potential surge in China’s coastal area caused by earthquake-induced tsunamis may reach a couple of meters high. Our preliminary results show that tsunamis generated in the Manila and Ryukyu trenches could pose a significant threat to Chinese coastal cities such as Shanghai, Hong Kong and Macao. However, we did not find the highest tsunami wave at Taiwan, partially because it lies right on the extension of an assumed fault line. Furthermore, we put forward a multi-scale model with higher resolution, which enabled us to investigate the edge waves diffracted around Taiwan Island with a closer view.     

地震海啸在东海大陆架上传播过程的数值模拟研究

在东海潜在震源区冲绳海槽假定了五个震源点,根据Steven地震海啸地震参数经验值作为初始条件,分别考虑6.5、7.0、7.5、8.0、8.5、9.0级地震条件下的30个震例,采用数值模拟的方法,对海啸在东海传播过程进行情境分析,特别是对上海沿岸地区可能会遭受的海啸灾害做了较为精细的研究.结果发现:小于8.0级的震例对上海地区几乎不会造成影响;8.0级震例只有最北端震源点震例会对上海地区有明显影响;8.5级以及9.0级震级基本上均会对上海沿岸地区造成较大的影响.特别是冲绳海槽北段9.0级震例可能会对上海沿岸局部地区造成危害,最大波高可达3.9m.     

冲绳海槽地区地震潜在海啸对中国东南沿海的影响研究

采用球坐标系下非线性浅水波方程, 研究日本本州M9.0大地震引发的海啸对中国东南沿海的影响, 并计算了冲绳海槽构造带上3个不同段落可能发生潜在地震引发的海啸, 分析这些海啸与日本大海啸的浪高和走时关系. 结果表明, 日本地震海啸模拟结果与日本当地报道及中国东南沿海7个验潮站的报道结果相符. 冲绳海槽构造带中段可能发生的3次不同震级（M7.0, M7.5, M8.0）潜在地震引发的海啸到达中国东南沿海的时间比日本海啸提前约4个小时, 从震源区传播3个多小时即可到达华东沿海部分验潮站. 冲绳海槽M7.5潜在地震海啸在验潮站上计算的波高与日本海啸相当, 中冲绳海槽M8.0潜在地震海啸在大陈站的波高将超过0.9 m, 在坎门站波高将超过1.8 m. 北冲绳海槽的潜在地震海啸威胁主要集中在江苏盐城、 上海一带, 南冲绳海啸主要对台湾东北部和浙江沿海产生威胁. 本文对冲绳海槽构造带上潜在地震引发海啸的模拟结果, 可为中国东南沿海地区的防震减灾、 海啸预警提供有意义的参考.      

Tsunami Assessment for Risk Management at Nuclear Power Facilities in Japan

The present study focuses on evaluation of the maximum and minimum water levels caused by tsunamis as risk factors for operation and management at nuclear power facilities along the coastal area of Japan. Tsunamis generated by submarine earthquakes are examined, basing literature reviews and databases of information on historical tsunami events and run-up heights. For simulation of water level along the coast, a numerical calculation system should be designed with computational regions covering a particular site. Also the calculation system should be verified by comparison of historical and calculated tsunami heights. At the beginning of the tsunami assessment, the standard faults, their locations, mechanisms and maximum magnitudes should be carefully estimated by considering historical earthquake-induced tsunamis and seismo-tectonics at each area. Secondly, the range of errors in the model parameters should be considered since earthquakes and tsunamis are natural phenomena that involve natural variability as well as errors in estimating parameters. For these reasons, uncertainty-induced errors should be taken into account in the process of tsunami assessment with parametric study of the tsunami source model. The element tsunamis calculated by the standard fault models with the errors would be given for the design. Then, the design tsunami can be selected among the element tsunamis with the most significant impact, maximum and minimum water levels, on the site, bearing in mind the possible errors in the numerical calculation system. Finally, the design tsunami is verified by comparison with the run-up heights of historical tsunamis, ensuring that the design tsunami is selected as the highest of all historical and possible future tsunamis at the site.     

Analysis of seismological and tsunami data from the 1993 Guam earthquake

The fault parameters of the Guam earthquake of August 8, 1993 are estimated from seismological analyses, and the possibility of identifying the actual fault plane from tsunami waveforms is tested. The Centroid Moment Tensor solution of long-period surface waves shows one nodal plane shallowly dipping to the north and the other nodal plane steeply dipping to the south. The seismic moment is 3.5×10²⁰ Nm and the corresponding moment magnitude is 7.7. The Moment Tensor Rate Function inversion ofP waves also yields a similar focal mechanism and seismic moment. The point source depth is estimated as 40–50 km.This earthquake generated tsunamis that propagated toward the Japanese coast along the Izu-Bonin-Mariana ridge system. The tsunamis are recorded on ocean bottom pressure gauges and tide gauges. Numerical computation of tsunamis shows that the computed waveforms from the two possible fault planes match well with the observed tsunami waveforms. The numerical computation also shows that the tsunami waveforms at Guam Island, just above the fault, should contain useful information regarding the identification of the actual fault plane. However, the current sampling rate of the tide gauges is so small that the records cannot help the identification.     

Forecasting tsunamis in Poverty Bay, New Zealand, with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.     

Forecasting tsunamis in Poverty Bay,New Zealand,with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.

     

Could a 1755-Like Tsunami Reach the French Atlantic Coastline? Constraints from Twentieth Century Observations and Numerical Modeling

The tsunami generated by the 1 November, 1755 earthquake off the coast of Portugal affected mainly the coastlines of the Iberian Peninsula and Northwest Morocco, but was also observed in some places along the North Atlantic coasts. To determine whether the event could have effected the French coastline, we conducted a study to search for signs of the tsunami in historical records from all tide gauge stations off the French Atlantic coast during the twentieth century, specifically for the 28 February, 1969 and the 26 May, 1975 tsunamis that were recorded by the Portuguese tide gauge network. Because many recordings are available in La Rochelle (located on the southwest coast of France), we focused our study on this harbor. The analysis of the tide gauge data shows no evidence for tsunamis in La Rochelle, neither in 1969 nor in 1975. To confirm this lack of tsunami signals, we used nonlinear, shallow water equations to compute the tsunami propagation to the French Atlantic coastline for both 1969 and 1975 events. Results obtained from these simulations confirm otherwise unnoticeable wave amplitudes at La Rochelle harbor. In a second step, tsunamis from three different scenarios for the 1755 earthquake were modeled to estimate the impact of such a tsunami on the French Atlantic coast, with a focus on La Rochelle harbor. A comparison of the functions of tide configuration was made in order to analyse the difference in impact. The results show that, while the harbor is poorly impacted, several areas (western part of the island of Ré and northern coast of the island of Oléron) may have experienced a moderate impact from 0.5 to 1 m, especially since the tide was high at the time of arrival, possibly causing local inundations in lowland areas.     

A Theoretical Comparison of Tsunamis from Dislocations and Landslides

— We use simple physical models to evaluate and compare the orders of magnitude of the energy generated into a tsunami wave by seismic dislocations and underwater slumps. We conclude that the two sources can generate tsunamis of comparable total energy. However, the slumping source is shown to be fundamentally dipolar in nature, which results in a low-frequency deficiency in the far-field. These simple conclusions corroborate the interpretation of the 1998 Papua New Guinea tsunami as being generated by an underwater slump.