Modeling the near-field effects of the worst-case tsunami in the Makran subduction zone

As a first step towards the development of inundation maps for the northwestern Indian Ocean, we simulated the near-field inundation of two large tsunami in the Makran subduction zone (MSZ). The tsunami scenarios were based on large historical earthquakes in the region. The first scenario included the rupture of about 500 km of the plate boundary in the eastern MSZ, featuring a moment magnitude of M_w 8.6. The second scenario involved the full rupture of the plate boundary resulting from a M_w 9 earthquake. For each scenario, the distribution of tsunami wave height along the coastlines of the region is presented. Also, detailed runup modeling was performed at four main coastal cities in the region for the second scenario. To investigate the possible effect of splay fault branching on tsunami wave height, a hypothetical splay fault was modeled which showed that it can locally increase the maximum wave height by a factor of 2. Our results showed that the two tsunami scenarios produce a runup height of 12-18 m and 24-30 m, respectively. For the second scenario, the modeled inundation distance was between 1 and 5 km.     

A Note on Tsunami Bore Front

A note is presented on tsunami bore front. This tsunami bore front is an old dynamical problem but also a new problem to be understood. The tsunami event on 2004 December 26 has raised this is an urgent problem. The author introduces here a model in order to see a hydrodynamical specific property of the tsunami bore front. This modeling gives us a new understanding about what mechanics is for the interested tsunami bore front, especially, around a coastal zone. This work adds a new understanding about mechanics of water motions as the tsunamis generated by the earthquake undersea at a distant area from the coast. The model in this work points out a specific transitional pattern as a function of time and space of tsunami bore front. This modeling gives what is essential at considering tsunami bore front.     

A Note on Tsunami Bore Front

A note is presented on tsunami bore front. This tsunami bore front is an old dynamical problem but also a new problem to be understood. The tsunami event on 2004 December 26 has raised this is an urgent problem. The author introduces here a model in order to see a hydrodynamical specific property of the tsunami bore front. This modeling gives us a new understanding about what mechanics is for the interested tsunami bore front, especially, around a coastal zone. This work adds a new understanding about mechanics of water motions as the tsunamis generated by the earthquake undersea at a distant area from the coast. The model in this work points out a specific transitional pattern as a function of time and space of tsunami bore front. This modeling gives what is essential at considering tsunami bore front.     

Numerical assessment of tsunami attack on a rubble mound breakwater using OpenFOAM®

A numerical assessment study of tsunami attack on the rubble mound breakwater of Haydarpasa Port, located at the southern entrance of the Istanbul Bosphorus Strait in the Sea of Marmara, Turkey, is carried out in this study using a Volume-Averaged Reynolds-Averaged Navier-Stokes solver, IHFOAM, developed in OpenFOAM^® environment. The numerical model is calibrated with and validated against the data from solitary wave and tsunami overflow experiments representing tsunami attack. Furthermore, attack of a potential tsunami near Haydarpasa Port is simulated to investigate effects of a more realistic tsunami that cannot be generated in a wave flume with the present state of the art technology. Discussions on practical engineering applications of this type of numerical modeling studies are given focusing on pressure distributions around the crown-wall of the rubble mound breakwater, and the forces acting on the single stone located behind the crown-wall at the rear side of the breakwater. Numerical modeling of stability/failure mechanism of the overall cross-section is studied throughout the paper.The present study shows that hydrodynamics along the wave flume and over the breakwater can be simulated properly for both solitary wave and tsunami overflow experiments. Stability of the overall cross-section can only be simulated qualitatively for solitary wave cases; on the other hand, the effect of the time elapsed during tsunami overflow cannot be reflected in the simulations using the present numerical tool. However, the stability of the overall cross-section under tsunami overflow is assessed by evaluating forces acting on the rear side armor unit supporting the crown-wall of the rubble mound breakwater as a practical engineering application in the present paper. Furthermore, two non-dimensional parameters are derived to discuss the stability of this armor unit; and thus, the stability condition of the overall cross-section. Approximate threshold values for these non-dimensional parameters are presented comparing experimental and numerical results as a starting point for engineers in practice. Finally, investigations on the solitary wave and tsunami overflow experiments/simulations are extended to the potential tsunami simulation in the scope of both representation of a realistic tsunami in a wave flume and stability of the rubble mound breakwater.     

Analysis of the Tsunami of December 26, 2004, on the Kerala Coast of India—Part I: Amplitudes

The Indian Ocean tsunami of December 26, 2004, not only affected the Bay of Bengal coast of India but also part of the Arabian Sea coast of India. In particular, the tsunami caused loss of life and heavy damage on some parts of the Kerala coast in southwest India. The tsunami traveled west, south of Sri Lanka, and some of the tsunami energy was diffracted around Sri Lanka and the southern tip of India and moved northward into the Arabian Sea. However, tsunami, being a long gravity wave with a wave length of a few hundred kilometers, has to take a wide turn. In that process, it missed the very southern part of the Kerala coast and did not achieve large amplitudes there. However, further north, the tsunami achieved amplitudes of upto 5 m and caused loss of life and significant damage. Here we identify the physical oceanographic processes that were responsible for selective amplification of the tsunami in certain locations.     

Analysis of the Tsunami of December 26, 2004, on the Kerala Coast of India—Part I: Amplitudes

The Indian Ocean tsunami of December 26, 2004, not only affected the Bay of Bengal coast of India but also part of the Arabian Sea coast of India. In particular, the tsunami caused loss of life and heavy damage on some parts of the Kerala coast in southwest India. The tsunami traveled west, south of Sri Lanka, and some of the tsunami energy was diffracted around Sri Lanka and the southern tip of India and moved northward into the Arabian Sea. However, tsunami, being a long gravity wave with a wave length of a few hundred kilometers, has to take a wide turn. In that process, it missed the very southern part of the Kerala coast and did not achieve large amplitudes there. However, further north, the tsunami achieved amplitudes of upto 5 m and caused loss of life and significant damage. Here we identify the physical oceanographic processes that were responsible for selective amplification of the tsunami in certain locations.     

日本“3·11”海啸引起的乐清湾共振研究

本文基于TELEMAC-2D模型建立太平洋区域海啸传播模型,模拟2011年日本“3·11”海啸事件下海啸波的传播。使用实测数据对该模型进行验证,在模型验证良好的基础上分析日本“3·11”海啸事件对乐清湾的影响。通过频谱分析得到“3·11”海啸激发的乐清湾内240、180和103 min这3个主导模态的幅值及其相位。通过白噪声实验对乐清湾的固有共振特征进行估算,进一步支持了乐清湾在上述3个模态发生共振这一结论。白噪声实验还表明,海啸等海洋灾害发生时会在乐清湾湾顶及湾口处产生较大的增水,该结论对乐清湾内海洋灾害风险防范具有指导意义。     

Wave dispersion study for tsunami propagation in the Sea of Marmara

In this paper the aim is to investigate whether there are differences between the dispersion and non-dispersion solutions on tsunami propagation. For this purpose, two numerical models of tsunami propagation are compared. One of these numerical models is a nondispersive model that uses Saint Venant equations and the other is a dispersive model that uses Boussinesq equations. The tsunamis resulting from a submarine mass failure (SMF) which is settled at the bottom of the north eastern Sea of Marmara are examined. An analytical solution considering wave dispersion is developed for obtaining near-field tsunami amplitudes above the submarine mass failure. Numerical modeling is used at the sea surface from the common boundary called as liquid boundary with incident waves up to the coastal regions to get the tsunami amplitudes. The output of the analytical model is taken as the disturbances for the numerical method. In the numerical solutions TELEMAC-2D software system is used for both dispersive and nondispersive modeling. The results of the dispersive and nondispersive models are compared to each other. Both temporal and spatial differences in the amplitudes and wave shapes are examined. The obtained results demonstrate that there are no noticeable differences between the dispersion and non-dispersion solutions except some special cases and some special landslide velocities.     

Reconstructing tsunami run-up from the characteristics of tsunami deposits on the Thai Andaman Coast

Tsunamis can leave deposits on the land surface they inundate. The characteristics of tsunami deposits can be used to calculate tsunami run-up height and velocity. This paper presents a reconstruction of tsunami run-up from tsunami deposit characteristics in a simple mathematical model. The model is modified and applied to reconstruct tsunami run-ups at Ao Kheuy beach and Khuk Khak beach, Phangnga province, Thailand. The input parameters are grain-size and maximum run-up distance of the sediment. The reconstructed run-up heights are 4.16–4.91 m at Ao Kheuy beach and 5.43–9.46 m at Khuk Khak beach. The estimated run-up velocities (maximum velocity) at the still water level are 12.78–19.21 m/s. In the area located 70–140 m inland to the end of run-up inundation, estimated mean run-up velocities decrease from approximately 1.93 m/s to 0 m/s. Reasonably good agreements are found between reconstructed and observed run-up heights. The tsunami run-up height and velocity can be used for risk assessment and coastal development programs in the tsunami affected area. The results show that the area from 0 to 140 m inland was flooded by high velocity run-ups and those run-up energies were dissipated mainly in this area. The area should be designated as either an area where settlement is not permitted or an area where effective protection is provided, for example with flood barriers or forest.     

Tsunami Travel Time Atlas for the Atlantic Ocean

Compared to the Pacific Ocean, tsunamis are rare both in the Atlantic and Indian Oceans. However, the December 26, 2004, tsunami demonstrated that, no matter how rare they may be, when a major tsunami occurs, it could be very disastrous. The most basic information in tsunami warning center requires are charts showing tsunami travel times to various locations around the rim of the ocean. With this in mind, a tsunami travel time atlas for the Atlantic Ocean is in preparation. The Caribbean Sea is also included in this Atlas, as it is more or less a part of the Atlantic Basin.     

有限断层模型在2015年9月16日智利Mw8.3级地震海啸数值模拟中的应用与评估

2015年9月16日22时54分（当地时间）智利中部近岸发生M_w8.3级地震,震源深度25 km。同时,强震的破裂区长200 km,宽100 km,随之产生了中等强度的越洋海啸。海啸影响了智利沿岸近700 km的区域,局部地区监测到近5 m的海啸波幅和超过13 m的海啸爬坡高度。太平洋区域的40多个海啸浮标及200多个近岸潮位观测站详细记录了此次海啸的越洋传播过程,为详细研究此次海啸近场及远场传播及演化规律提供了珍贵的数据。本文选择有限断层模型和自适应网格海啸数值模型建立了既可以兼顾越洋海啸的计算效率又可以实现近场海啸精细化模拟的高分辨率海啸模型。模拟对比分析了海啸的越洋传播特征,结果表明采用所建立的模型可以较好地再现远场及近场海啸特征,特别是对近场海啸的模拟结果非常理想。表明有限断层可以较好地约束近场、特别是局部区域的破裂特征,可为海啸预警提供更加精确的震源信息,结合高分辨率的海啸数值预报模式实现海啸传播特征的精细化预报。本文结合观测数据与数值模拟结果初步分析了海啸波的频散特征及其对模型结果的影响。同时对观测中典型的海啸波特征进行的简要的总结。谱分析结果表明海啸波的能量主要分布在10~50 min周期域内。这些波特征提取是现行海啸预警信息中未涉及,但又十分重要的预警参数。进一步对这些波动特征的详细研究将为海啸预警信息及预警产品的完善提供技术支撑。     

Tsunami Travel Time Atlas for the Atlantic Ocean

Compared to the Pacific Ocean, tsunamis are rare both in the Atlantic and Indian Oceans. However, the December 26, 2004, tsunami demonstrated that, no matter how rare they may be, when a major tsunami occurs, it could be very disastrous. The most basic information in tsunami warning center requires are charts showing tsunami travel times to various locations around the rim of the ocean. With this in mind, a tsunami travel time atlas for the Atlantic Ocean is in preparation. The Caribbean Sea is also included in this Atlas, as it is more or less a part of the Atlantic Basin.     

Positioning requirements for the tsunami warning system

Abstract

Canada has increased the number of tsunami warning stations on the Pacific Coast from two to three. The last gauge was installed at the north end of Vancouver Island, thereby filling a large gap previously existing and providing full coverage along the coast. The record of gauges at two of the three locations is accessible either by telephone or by means of meteor burst communication, alleviating the difficulties experienced during the tsunami threat of May 6, 1986, when telephone communications were disrupted by heavy use. The gauge at Langara Island will be relocated in a more accessible and also a more tsunami‐responsive location at Rennell Sound in the Queen Charlotte Islands. All tsunami gauges also serve as tide gauges, recording the water level every 15 min. In the event of a tsunami, the recording interval can be altered to every 60 s. Suggestions have been made that Canada attempt deep‐sea recording of tsunamis off its Pacific Coast. Although this would be of great scientific value, no such program is contemplated at this time.     

Initial attack of large-scaled tsunami on ship motions and mooring loads

In this paper, we propose a numerical simulation procedure of moored ship motions due to initial attack of large-scaled tsunamis and investigate the effects on the motions and mooring loads. The effect of methodology on selection of tsunami wave components and of the drag forces are then considered by using the numerical simulation method, applying to several case studies for LNG-carrier. Large ship motions and excessive mooring loads beyond the safe working loads are induced by the resonant tsunami wave components in the sway and surge motions and drag forces.     

2016年全球地震海啸监测预警与数值模拟研究

回顾了国家海洋环境预报中心(国家海洋局海啸预警中心)2016年全球地震海啸监测预警的总体状况, 并基于震源生成模型和海啸传播数值模型的计算结果详细介绍了几次主要海啸事件及其影响特性。2016年全年国家海洋环境预报中心总共对全球6.5级(中国近海5.5级)以上海底地震响应了45次,发布海啸信息81期, 没有发生对我国有明显影响的海啸。结合精细化的数值模拟结果和浮标监测数据,重点介绍了苏门达腊7.8级地震海啸、厄瓜多尔7.8级地震海啸、新西兰7.1级和7.8级地震海啸, 以及所罗门7.8级地震海啸的波动特征和传播规律, 模拟结果与实测海啸波符合较好。针对厄瓜多尔7.8级地震海啸事件, 本文比较分析了均匀断层模型和多源有限断层模型对模拟结果的影响; 针对新西兰7.1级地震海啸, 探讨了色散效应对海啸波在大水深、远距离传播过程的影响规律。     

Tsunami amplitudes from local earthquakes in the pacific northwest region of North America part 1: The outer coast

Abstract

Maximum tsunami amplitudes that will result from major earthquakes in the Pacific Northwest region of North America are considered. The modeled region encompasses the coastlines of British Columbia in Canada, and Washington and Oregon in the United States. Three separate models were developed for the outer coast and one model for the system consisting of the Strait of Georgia, Juan de Fuca Strait, and Puget Sound (GFP model) (Part 2). Three different source areas were considered for the outer coast models and the resulting tsunami was propagated to the entrance of Juan de Fuca Strait. Using the output from the other models, the GFP model was run. The results showed that large tsunami amplitudes can occur on the outer coast, whereas inside the GFP system, unless the earthquake occurs in the system itself, no major tsunami will result (Part 2).     

Predicting tsunami arrivals: Estimates and policy implications

Tsunami is one of a few kinds of natural disasters that leave people some time for escape. This escape time, which is essentially the time for the giant wave to propagate from the epicentre to a coast, has to be estimated without delay upon the occurrence of the incident. With the advancement of water wave theories, much work has been done to model the propagation of tsunamis from deep oceans to shallow water. The authors argue that while much emphasis has been put on the expansion of the high-tech early warning system and the development of complicated tsunami models, a simple-to-use yet accurate predictive model is still wanting. This paper presents a handy linear wave model, which is capable of estimating the arrival time of a tsunami with very good accuracy, as has been verified by comparison with past incidents. With the availability of such a simple model, even local communities without access to a high-tech warning system can readily estimate the time left for emergency evacuation.     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

Assessment of tsunami hazard for coastal areas of Shandong Province,China

Shandong province is located on the east coast of China and has a coastline of about 3100 km. There are only a few tsunami events recorded in the history of Shandong Province, but the tsunami hazard assessment is still necessary as the rapid economic development and increasing population of this area. The objective of this study was to evaluate the potential danger posed by tsunamis for Shandong Province. The numerical simulation method was adopted to assess the tsunami hazard for coastal areas of Shandong Province. The Cornell multi-grid coupled tsunami numerical model (COMCOT) was used and its efficacy was verified by comparison with three historical tsunami events. The simulated maximum tsunami wave height agreed well with the observational data. Based on previous studies and statistical analyses, multiple earthquake scenarios in eight seismic zones were designed, the magnitudes of which were set as the potential maximum values. Then, the tsunamis they induced were simulated using the COMCOT model to investigate their impact on the coastal areas of Shandong Province. The numerical results showed that the maximum tsunami wave height, which was caused by the earthquake scenario located in the sea area of the Mariana Islands, could reach up to 1.39 m off the eastern coast of Weihai city. The tsunamis from the seismic zones of the Bohai Sea, Okinawa Trough, and Manila Trench could also reach heights of >1 m in some areas, meaning that earthquakes in these zones should not be ignored. The inundation hazard was distributed primarily in some northern coastal areas near Yantai and southeastern coastal areas of Shandong Peninsula. When considering both the magnitude and arrival time of tsunamis, it is suggested that greater attention be paid to earthquakes that occur in the Bohai Sea. In conclusion, the tsunami hazard facing the coastal area of Shandong Province is not very serious; however, disasters could occur if such events coincided with spring tides or other extreme oceanic conditions. The results of this study will be useful for the design of coastal engineering projects and the establishment of a tsunami warning system for Shandong Province.     

Optical measurements of tsunami inundation through an urban waterfront modeled in a large-scale laboratory basin

This paper presents optical measurements of tsunami inundation through an urban waterfront in a laboratory wave basin. The physical model was constructed at 1:50 scale and was an idealization of the town of Seaside, Oregon. The fixed-bed model was designed to study the initial inundation zone along an urban waterfront, such that the flow around several large buildings could be observed. This paper presents an analysis of the optical measurements made with two overhead video cameras, focusing on tracking the leading edge of the tsunami inundation through the urban waterfront and quantifies the accuracy of the algorithm used to track the edge. The results show that the methodology provides high-resolution information in both time and space of the leading edge position, and that these data can be used to quantify the influence of large macro-roughness features on the tsunami inundation processes in laboratory settings. The overall effect of the macro-roughness was to decrease the bore propagation speed relative to the control section with no macro-roughness. The bore speed could be reduced by as much as 40% due to the presence of the macro-roughness relative to the control section.