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.     

海啸波对近岸岛礁影响的数值模拟研究

基于Okada有限断层模型和非线性浅水波方程,结合高精度嵌套网格建立了越洋（中国近海）-局部-近岸岛礁的海啸生成与传播的数值模型。以三亚凤凰岛为例,首先针对2011日本地震海啸,模拟分析了海啸波沿中国沿海大陆架的传播特征及对凤凰岛的影响规律。在取得验证结果的基础上,进一步讨论了中国近海的马尼拉海沟和琉球海沟的潜在海啸源,以及环太平洋的21个潜在特大越洋海啸对凤凰岛的影响特征。依据海啸波在抵达凤凰岛的波浪特征,结合傅里叶频谱分析方法,探索了近岸岛礁对海啸波的放大效应。结果表明,中国近海一般震级的海啸和特大越洋海啸对凤凰岛存在一定影响,最大波幅接近1 m,传播时间从3 h到27 h不等。受三亚东南半岛的影响,琉球海沟激发的海啸和越洋海啸在凤凰岛的放大效应相对于马尼拉海沟较小,其频率集中在0.8×10-4~2×10-4 Hz。马尼拉海沟产生的海啸波在凤凰岛产生了较为显著的放大效应,对于凤凰岛是值得关注的高风险海啸源。     

“3.11”日本地震海啸及南海和琉球群岛假想地震海啸对福建近海影响的数值模拟

福建地处西北太平洋沿岸,在环太平洋地震带附近,是海啸灾害潜在风险区.＂3.11＂日本地震海啸,福建沿岸验潮站就监测到其海啸波.利用CTSU地震海啸数值模式,模拟了＂3.11＂日本地震海啸对福建近海的影响,模拟结果与实况较吻合.同时,利用该数值模式模拟分析了可能来自于琉球群岛和南海附近海域的地震海啸对福建近海的影响,分析表明,如果在琉球群岛海域（28.0°N,129.0°E）发生8.8级地震,引发的海啸波将在4.5 h左右抵达福建北部海岸,最大海啸波幅可达2 m;如果在马尼拉海沟附近海域（17.5°N,119.0°E）发生8.8级地震,引发的海啸波将在4 h左右抵达福建南部海岸,最大海啸波幅可达3 m,均会给福建沿海地区带来灾害性影响.为此,本文亦针对性提出了防范地震海啸的一些措施与建议,为福建省海洋防灾减灾提供参考.     

Scenarios of local tsunamis in the China Seas by Boussinesq model

The Okinawa Trench in the East China Sea and the Manila Trench in the South China Sea are considered to be the regions with high risk of potential tsunamis induced by submarine earthquakes. Tsunami waves will impact the southeast coast of China if tsunamis occur in these areas. In this paper, the horizontal two-dimensional Boussinesq model is used to simulate tsunami generation, propagation, and runnp in a domain with complex geometrical boundaries. The temporary varying bottom boundary condition is adopted to describe the initial tsunami waves motivated by the submarine faults. The Indian Ocean tsunami is simulated by the numerical model as a validation case. The time series of water elevation and runup on the beach are compared with the measured data from field survey. The agreements indicate that the Boussinesq model can be used to simulate tsunamis and predict the waveform and runup. Then, the hypothetical tsunamis in the Okinawa Trench and the Manila Trench are simulated by the numerical model. The arrival time and maximum wave height near coastal cities are predicted by the model. It turns out that the leading depression N-wave occurs when the tsunami propagates in the continental shelf from the Okinawa Trench. The scenarios of the tsunami in the Manila Trench demonstrate significant effects on the coastal area around the South China Sea.     

浙江沿海潜在区域地震海啸风险分析

采用COMCOT海啸模型建立三重网格模型模拟了2011年3月11日日本东北部9.0级地震引发的海啸发生、发展以及在我国东南沿海传播过程。震源附近浮标站以及浙江沿海的潮位站实测资料验证结果显示,大部分监测站首波到达时间和海啸波的计算值相差在15%以内,表明模型可较好的模拟海啸在计算域内的传播过程。研究表明日本南海海槽、冲绳海槽以及琉球海沟南部是影响浙江沿海主要的区域潜在震源,通过情景计算分别模拟3个潜在震源9.1级、8.0级和8.7级地震引发的海啸对浙江沿海的海啸风险,计算结果表明,海啸波产生后可在3~8h内传至浙江省沿岸,海啸波达1~3m,最大可达4m,此时浙江沿岸面临Ⅲ~Ⅳ级海啸风险,达到淹没至严重淹没等级。     

从2011 海啸演习看中国海啸危害

2011 年11 月, 联合国教科文组织政府间海洋学委员会(IOC/UNESCO)发起了代号为“Exercise Pacific Wave 11”的跨国海啸演习, 演习区域为整个太平洋海区。中国应邀参加了本次演习, 演习在我国海域地震带上假设了两个震源, 分别位于琉球海沟和马尼拉海沟。为了评估这两个潜在海啸对我国的影响, 本文采用数值模拟的方式, 对其进行了计算。计算结果表明: 这两处震源所引发的海啸均能对我国造成灾害性影响, 受影响严重的省市有江苏、上海、浙江、福建、广东和海南等; 从传播时间图上看, 若发生大海啸, 海啸波将在5 h 内, 陆续影响沿海各省市; 相比琉球海沟, 马尼拉海沟震源产生的海啸会更快到达我国沿岸。     

Worst scenarios of tsunami effects along the Mediterranean coast of Egypt

Scenarios of tsunami effects represent a very useful technique for the definition and evaluation of tsunami hazard and risk for the Egyptian coast. This paper is an attempt to develop different worst scenarios of tsunamigenesis toward the Egyptian Coast for five segment localities along three different sub-regions (Hellenic Arc, Cyprean Arc and Levantine Coast) in the eastern Mediterranean Sea. These segments are the southwest Hellenic Arc, southeast Hellenic Arc, northeast Hellenic Arc, west of Cyprean Arc and Levantine. For each of them, the scenario takes into account a seismic fault capable of generating an earthquake with magnitude equal to or larger than the highest magnitude registered in that region in historical times. Then the ensuing tsunamis are simulated numerically, highlighting the basic features of the wave propagation and roughly identifying the coastal sectors that are expected to suffer the heaviest tsunami effects. The output data indicated that the first wave of tsunamis from different segments attacked the nearest reference localities (city located nearest each segment) along the Egyptian shore between 28 and 50 min after an earthquake. Tsunamis from these earthquake scenarios produced maximum run-up heights ranging from 1.7 to 9.4 m at the shore. A Beirut Thrust scenario (Levantine segment) included the fact that only a small portion of the fault extended out into the sea, leading to a small effective tsunami source area. In contrast, the southwest Hellenic Arc segment (as in the A.D. 365 earthquake) has high displacement (15 m) and a long extensional fault, forming a highly effective tsunami source area.     

Tsunami Inundation Modelling for the Coast of Kerala,India

The disastrous tsunami of December 26, 2004, exposed the urgent need for implementing a tsunami warning system. One of the essential requirements of a tsunami warning system is the set up of tsunami inundation models which can predict inundation and run-up along a coastline for a given set of seismic parameters. The Tsunami Warning Centre and the State/District level Disaster Management Centres should have tsunami inundations maps for different scenarios of tsunami generation. In the event of a tsunamigenic earthquake, appropriate decisions on issue of warnings and/or evacuation of coastal population are made by referring to such maps. The nature of tsunami inundation and run-up along the Kerala coast for the 2004 Sumatra and 1945 Makran, and a hypothetical worst-case scenario are simulated using the TUNAMI N2 model and the results are presented in this paper. Further, scenarios of tsunami inundation arising out of possible rise in sea level as projected by the Intergovernmental Panel on Climate Change (IPCC 2001) are also simulated and analysed in the paper. For the study, three representative sectors of the Kerala coast including the Neendakara-Kayamkulam coast, which was the worst hit by the 2004 tsunami, are chosen. The results show that the southern locations and certain locations of central Kerala coast are more vulnerable for Sumatra when compared to Makran 1945 tsunami. From the results of numerical modelling for future scenarios it can be concluded that sea level rise can definitely make pronounced increase in inundation in some of the stretches where the backshore elevation is comparatively low.     

中国砂质海滩区域差异分布的构造成因及其堆积地貌研究

砂质海滩的形成与发育是联系地壳深部过程和地表过程的综合体现．本文根据中国沿岸砂质海滩分布状况调查结果,通过分析其分布规律与中国东部沿海中新生代构造地质分区不同海洋沉积环境的内在联系,探讨了砂质海滩区域性差异分布和堆积地貌的成因特点,从而为我国海岸资源利用的宏观配置和管理提供科学依据．研究得出：（1）中国沿岸砂质海滩的形成与发育受中新生代地质构造带控制,尤其受新构造期地壳断块升降活动差异性的影响,呈现出区域性富集的格局．（2）中国砂质海滩基本形成在构造隆起带岬湾沿岸,并且富集在其中新构造期处于上升区的海岸,这主要在于其区域性海岸环境具有形成砂质沉积的物质基础．（3）构造隆起带滨海泥沙的输沙途径及其再分配形式复杂多样,从而产生了各种不同的砂质海滩堆积地貌类型．（4）根据海岸地貌和沉积成因模式,将中国沿岸砂质海滩划分为岬湾海岸叠置沙体沉积形式、滨外坝沙体连续沉积形式和浪控三角洲平原岸改造沉积形式等3种基本模型．（5）沿海晚第四世纪古夷平地貌面与现代海平面相对高程分布的分析表明,山丘岬湾岸型砂质海滩主要形成在新构造运动下沉区段或构造相对稳定区段,而台地岬湾岸型和泻湖一沙坝岸型形成在上升区段,浪控三角洲平原岸型则分布在断块盆地．     

渤海海域地震海啸灾害概率性风险评估

渤海作为我国地震活动性最为活跃的近海,其地震海啸风险不可忽视。本文应用概率性海啸风险评估方法对渤海周边区域的海啸风险进行评估。根据历史地震目录建立了渤海区域的震级-频率关系,基于蒙特卡洛算法随机生成了一套10万年的地震目录,最终通过对地震事件的海啸数值模拟及最大波幅的统计分析给出了环渤海区域典型重现期的最大波幅分布以及重点城市的海啸波幅曲线。评估结果表明,渤海地区海啸风险主要集中在渤海湾和莱州湾周边,波幅可达到1~3 m,辽东湾地区海啸风险较低。     

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.     

越洋海啸的数值模拟及其对我国的影响分析

简要介绍了地震海啸产生的物理机制、海啸波在大洋中的传播特性以及海啸所具有的超强破坏力可能引发的巨大灾害;概述了全球地震海啸发生的频率和太平洋区域历史海啸的时空分布;整理分析了我国沿海发生海啸的频次和空间分布。针对越洋海啸传播的特点,采用基于波浪追逐原理和自适应网格加密技术的海啸数值模型对1960智利海啸进行了数值模拟,将模拟的结果与历史记录进行了对比,验证了模型的可靠性。通过对数值模拟结果的分析,初步讨论了我国沿海地区越洋海啸的危险性,并定量阐述了越洋海啸对我国各海区的影响。     

南海潜在海啸灾害的模拟

结合南海海域的地形条件、地质构造、地震学特征以及历史地震记录,在回顾总结国内外学者研究的基础上,分析了南海可能引发地震海啸的震源区域,并讨论了在我国南海沿岸发生海啸灾害的潜在可能性。采用目前在国际上广泛使用的COMCOT海啸模式,对马尼拉海沟的潜在地震引发的海啸进行了数值模拟计算,计算中包含了由地震参数到海面初始变形的转换、海啸的深水传播过程以及海啸的浅水传播过程。采用三重嵌套网格,外层网格对应于大范围的深水区域,使用球坐标系下的线性控制方程;第二层网格对应中等范围的较浅水区域,使用球坐标系下的非线性控制方程;第三层网格对应小范围的浅水区域,使用直角坐标系下的非线性控制方程。由模拟计算得到的海啸传时分布、近岸海面升降强度、四个特定点上海面高度随时间变化等的结果表明,我国南海沿岸遭受海啸袭击的可能性是存在的,应进一步对南海海啸进行监测、预警和研究。COMCOT模式性能良好,可用于对南海潜在地震海啸的进一步模拟研究。     

南海主要岛礁概率海啸风险评估

影响地震海啸的震源参数众多且具有很强的不确定性,充分评估海啸风险需要大量的情景模拟.本文基于建立的概率海啸风险模型,采用一种高效的海啸模拟方法,评估了南海主要岛礁的概率海啸风险.通过对历史地震数据的分析,综合考虑震级、震中位置、震源深度的随机性,形成了百万数量级的潜在地震情景集,并通过叠加近似方法实现了大量地震情景引发...     

天文潮与海啸耦合数学模型研究——以东中国海温州湾为例

在COMCOT海啸数学模型中加入潮汐边界条件,建立了东中国海天文潮与海啸耦合数学模型。在琉球海沟内侧设计震级为7.6级的海底地震,根据地震板块的错动方向不同,设计正波先行与负波先行两种海啸波,通过调整海啸波发生时间,使海啸波波峰遭遇温州湾天文高潮位。将天文潮与海啸耦合模型计算结果与线性叠加计算结果进行比较,结果表明:无论正波先行还是负波先行,天文潮与海啸耦合计算相比线性叠加的结果,海啸波的到达时间均有所提前;而从海啸波波高来看,线性叠加的计算结果则比耦合计算结果偏高。     

Numerical modeling of the propagation and strengthening of tsunami waves near the Crimean Peninsula and the Northeast Coast of the Black Sea

The linear model of long waves is used for the evaluation of the parameters of tsunami waves along the South Coast of Crimea, in the near-Kerch zone, and near the northeast coast of the Black Sea. Our numerical investigations are carried out for 24 probable locations of the elliptic zones of tsunami generation over the continental slope of the basin. The amplitude characteristics of tsunamis are computed for 27 sites of the Black-Sea coast. It is shown that significant strengthening of tsunami waves is possible in the course of their propagation toward the coast. The highest waves are formed at the sites of the coast closest to the seismic source. The dependence of the intensity of tsunami waves along the Black-Sea coast on the location of the seismic source and its magnitude is analyzed.     

Numerical analysis of the propagation and amplification of tsunami waves of seismic generation in the Sea of Azov

The numerical analysis of the evolution of tsunamis is performed for the Sea of Azov. Our calculations are carried out on a grid with steps of 500 m, as applied to seaquakes with magnitudes within the range 6–7 for 18 circular zones of generation covering almost the entire water area of the sea. It is shown that the oscillations of the sea level in the form of cellular waves are formed as a result of the wave reflections from the coasts. Small areas of the elevated activity of waves are formed in the zones of irregularity of the coastline on the north coast of the sea and in the zones of shoals in the southeast part of the basin. On the basis of the determined values of extreme elevations and lowerings of the sea level, we can make a conclusion that the tsunami hazard is quite low for the coast of the Sea of Azov.     

Strongest tsunamis in the World Ocean and the problem of marine coastal security

We consider the problems of assessing tsunami danger for sea coasts taking into account the risk of the strongest tsunamis of seismic origin. We identify a class of particularly dangerous transoceanic events characterized by extremely high runups (up to 40–50 m) at extended coastal areas (up to 500–1000 km). In most cases these transoceanic tsunamis are caused by underwater mega-earthquakes with a magnitude of 9.0 or more occurring with a period between 200–300 and 1000–1200 years in some areas of subduction zones. The possibility of these earthquakes in subduction zones directly threatening a given coast should be taken into account in creating maps of tsunami zoning of any scale.     

基于数值模拟的渤海海域地震海啸危险性定量化研究

根据地震海啸产生的条件,结合渤海海域的地形特征、地质构造、地震学特征和历史地震及海啸记录对渤海海域潜在的地震海啸进行了数值模拟研究。分析了渤海可能引发地震海啸的震源区域,讨论了渤海发生海啸灾害的可能性。文中通过数值模拟再现了渤海历史上几次规模较大的地震事件可能引发的海啸情景,研究分析了可能的地震海啸在渤海及周边海域的传播过程及波动特征.地震海啸传播模型采用基于四叉树原理的自适应网格加密技术,有效解决了局部分辨率与计算效率之间的矛盾。数值计算包括地震海啸产生及传播过程。利用该模型对渤海潜在的地震海啸进行了数值计算,基于数值计算结果定量阐述了渤海海域潜在地震海啸对渤海局部岸段及北黄海沿岸的影响,给出了渤海可能地震海啸危险性划分;研究结果将为我国海啸危险性分析和海啸预警技术研究工作提供技术支持。     

2015年9月16日智利Mw 8.3级地震:海啸波从近场到远场的特征分析

On September 16, 2015, an earthquake with magnitude of M_w 8.3 occurred 46 km offshore from Illapel, Chile,generating a 4.4-m local tsunami measured at Coquimbo. In this study, the characteristics of tsunami are presented by a combination of analysis of observations and numerical simulation based on sources of USGS and NOAA. The records of 16 DART buoys in deep water, ten tidal gauges along coasts of near-field, and ten coastal gauges in the far-field are studied by applying Fourier analyses. The numerical simulation based on nonlinear shallow water equations and nested grids is carried out to provide overall tsunami propagation scenarios, and the results match well with the observations in deep water and but not well in coasts closed to the epicenter. Due to the short distance to the epicenter and the shelf resonance of southern Peru and Chile, the maximum amplitude ranged from 0.1 m to 2 m, except for Coquimbo. In deep water, the maximum amplitude of buoys decayed from9.8 cm to 0.8 cm, suggesting a centimeter-scale Pacific-wide tsunami, while the governing period was 13–17 min and 32 min. Whereas in the far-field coastal region, the tsunami wave amplified to be around 0.2 m to 0.8 m,mostly as a result of run-up effect and resonance from coast reflection. Although the tsunami was relatively moderate in deep water, it still produced non-negligible tsunami hazards in local region and the coasts of farfield.