日本南海海槽地震海啸对我国东海沿岸影响

本文利用数值模拟技术重现了1707年宝永地震海啸的传播过程,定量分析了我国东海沿岸海啸时空分布特征。计算结果表明,地震发生2.5小时后海啸波传至东海陆架,震后6小时浙江沿海地区遭到海啸的袭击,沿岸最大海啸波高为0.8米。通过海啸波在东海大陆架传播时海底地形与波幅的关系,研究分析了东海陆架缓变地形下海啸放大效应,为及时判断沿海可能的海啸强度和受灾程度提供了便捷的估算方法。此外,本文还评估了南海海槽发生极端地震时,中国东海沿岸的海啸危险性,为东海区域针对日本南海海槽进行海啸预警和减灾评估提供定量科学的参考。     

日本南海海沟潜在地震海啸对中国东部沿海影响的数值模拟研究

2011-03-11日本东北部地震海啸发生后,围绕日本南海海沟发生潜在9级地震的可能性也进一步增大。采用高精度高分辨率有限体积浅水波数值模型GeoClaw对日本东北大地震海啸的产生、传播过程进行数值模拟研究,模拟结果对比海上浮标观测数据及中国东部沿海验潮站记录数据可知首波波峰偏差小于10%,表明模型GeoClaw可以很好地模拟海啸在大洋中及东海大陆架的传播过程。利用模型GeoClaw对日本南海海沟9级地震断层模型进行海啸数值模拟研究。研究表明,地震引发的海啸波能量巨大,同时向各个方向传播。地震发生5h后海啸波到达中国东部福建沿海地区,进而影响浙江、上海、江苏等沿海地区,海啸波高可达到1 m以上。鉴于日本南海海沟大地震发生的历史性规律及对中国东部沿海地区造成的潜在灾害,需要对其进行相应的数值模拟研究。     

精细化地震海啸波流实时预警系统研究与应用

基于Okada模型和非线性浅水波模型,结合高精度多层嵌套网格针对我国浙江沿海的温州和台州地区建立了越洋–近海–局部的精细化地震海啸波流实时预警系统,近岸的分辨率为900 m。该预警系统包括了并行化的数值计算模块,基于Python 2D绘图库的计算结果可视化处理模块,以及通过Python语言将所有经过数值计算的图形与动画产品集成在一个网页上的产品集成模块。一旦地震发生,该系统可根据地震的震源参数信息在10 min内完成数值计算、可视化处理,以及产品集成。选取2011年日本东北9.0级地震海啸结合实测数值对该系统进行模拟验证,进一步应用该系统模拟计算了日本南海海槽和琉球海沟潜在极端海啸的影响规律。结果表明,该预警系统可有效地提高地震海啸实时预警的时效性和准确度,为海啸的预警、减灾,以及辅助决策提供科学依据。     

马尼拉海沟潜在地震海啸对我国华南沿海危险性研究

采用COMCOT模型建立南海多重嵌套网格对马尼拉海沟潜在震源进行数值模拟,评估马尼拉海沟地震海啸对我国华南附近海域的影响。敏感性分析结果表明:马尼拉海沟地震震源深度为33~40 km时,地震引发的海啸危险性最大;在此基础上,设置5种马尼拉海沟潜在地震海啸情景,其中情景二、四和五,其地震震级分别为Mw8.6、Mw8.8和Mw9.0,引发的海啸会严重威胁到我国华南沿海大部分海域,危险等级为Ⅲ-Ⅳ级,有淹没至严重淹没危险。     

2014年4月2日智利海啸数值模拟与分析

北京时间2014年4月2日智利北部近海发生8.2级地震。地震引发了海啸,南美智利、秘鲁等国沿岸的多个海洋站监测到了明显的海啸波动。文章利用国家海洋环境预报中心开发的CTSU海啸数值模式对这次智利海啸事件进行了数值模拟。模拟结果显示距离震源最近的智利北部受灾严重,秘鲁以及智利南部等海域的海啸波相对较小。沿海站点的第一波海啸波的数值模拟曲线与实测曲线基本吻合。由于数值模型的理想化和近岸水深地形数据分辨率不够,后续海啸波部分模拟结果与实测值存在一定误差。     

“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,均会给福建沿海地区带来灾害性影响.为此,本文亦针对性提出了防范地震海啸的一些措施与建议,为福建省海洋防灾减灾提供参考.     

马尼拉俯冲带潜在地震海啸对华南沿海的影响

马尼拉俯冲带潜在地震海啸对我国南部沿海城市构成巨大威胁,利用情景式数值模拟技术重构灾害过程并评估危险等级有助于理解南海海啸传播规律并指导预警预报和防灾减灾工作。根据美国太平洋海洋环境研究中心(Pacific Marine Environmental Laboratory, PMEL)发布的马尼拉俯冲带断层参数设计M_w 7.5、M_w 8.1和M_w 8.5三个震级下共19个震源,应用非静压海啸数值模型(Non-hydrostatic Evolution of Ocean WAVE, NEOWAVE)模拟各震源激发海啸在南海海盆的传播过程,通过最大波辐和测点时间序列发现海啸波能量传输分布并评估代表区域危险等级。研究表明, M_w 7.5级地震海啸对我国南部沿海的影响较低,波幅一般不超过30 cm; M_w 8.1级地震海啸对华南沿海主要造成太平洋海啸预警中心定义的Ⅱ或Ⅲ级海啸危险等级,海啸影响范围和能量分布特征由震源位置决定; M_w 8.5级地震海啸主要对中国沿海构...     

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

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

2010 年智利和2011 年日本地震海啸在我国东南沿岸诱导的波流特征及危险性分析

本文基于有限断层模型和OKADA 位错模型计算海表变形场作为初始条件,利用MOST 海啸数值模型模拟分析了2010年智利和2011 年日本地震海啸在我国东南沿海地区的海啸传播特征,海啸波模拟结果与观测数据吻合较好。重点研究分析了沈家门港口海域的海啸波流特征及其诱导的涡旋结构。研究结果表明：尽管两次事件的海啸源位置及破裂特征完全不同,但海啸波流在我国东南沿岸的分布特征大致相似;另外相对于海啸波幅而言,港湾中海啸流具有更强的空间差异性,港口入口、岬角地形处和岛屿间水道中往往会有强流存在。尽管这两次越洋海啸均未能在我国东南沿海引发淹没情形,但通过数值计算发现局部均存在超过3 m/s 以上的强流,因此进行海啸预警及风险管理时应综合考虑海啸波流的影响。     

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

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

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.     

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.     

Tsunamis in the central part of the Caspian Sea

This paper describes the geotectonics of the Caspian Sea basin and the seismicity of its central part. The seismicity analysis enables us to identify the most probable zones of tsunami generation. We also present a brief review of the historical records of tsunamis in the Caspian Sea. In order to estimate the tsunami risk, we used the method of numerical hydrodynamic simulation while taking into account the real topography of the Caspian Sea. The computation of the wave field for the possible tsunamis occurring in the central part of the Caspian Sea allowed us to estimate the maximum expected heights of the waves along the coast of the CIS countries (Russia, Azerbaijan, Kazakhstan, and Turkmenistan). On the basis of the earthquake statistics in the region and the results of numerical experiments, we show that the extreme wave heights can reach 10 m at certain parts of the coast. Such extreme events correspond to extended (up to 200 km) seismic sources with M _S ～ 8 and a recurrence period of T ≈ 1600 years. The tsunami wave heights are expected to be as high as 3 m for sources of lesser extent (<50 km) with earthquake magnitudes of M _S ～ 7 and a recurrence period of 200 years.     

Numerical Simulation of Tsunamis on the Tamil Nadu Coast of India

The State of Tamil Nadu was the most affected region in India during the tsunami of December 26, 2004, in the Indian Ocean, in terms of loss of life and damage. Numerical simulation was made for three tsunamis, the December 26, 2004, event, the Sumatra tsunami of 1833, and a hypothetical tsunami originating in the Andaman-Nicobar region. Since inundation is not included in these simulations, the tsunami amplitudes were deduced at the 10m depth contour in the ocean, off several locations on the coast of Tamil Nadu. The computed amplitudes appear reasonable as compared to known tsunami amplitudes from past events.     

有限断层模型在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周期域内。这些波特征提取是现行海啸预警信息中未涉及,但又十分重要的预警参数。进一步对这些波动特征的详细研究将为海啸预警信息及预警产品的完善提供技术支撑。     

The reduction effects of mangrove forest on a tsunami based on field surveys at Pakarang Cape,Thailand and numerical analysis

Using an integrated approach including satellite imagery analysis, field measurements, and numerical modeling, we investigated the damage to mangroves caused by the 2004 Indian Ocean tsunami at Pakarang Cape in Pang Nga Province, Thailand. Comparing pre- and post-tsunami satellite imagery of the study area, we found that approximately 70% of the mangrove forest was destroyed by the tsunami. Based on field observations, we found that the survival rate of mangroves increased with increasing stem diameter. Specifically, we found that 72% of Rhizophora trees with a 25–30 cm stem diameter survived the tsunami impact, whereas only 19% with a 15–20 cm stem diameter survived. We simulated the 2004 Indian Ocean tsunami using the nonlinear shallow-water wave theory to reproduce the tsunami inundation flow and investigated the bending moment acting on the mangrove trees. Results of the numerical model showed that the tsunami inundated areas along the mangrove creeks, and its current velocity reached 5.0 m s⁻¹. Based on the field measurements and numerical results, we proposed a fragility function for mangroves, which is the relationship between the probability of damage and the bending stress caused by the maximum bending moment. We refined the numerical model to include the damage probability of mangrove forests using the obtained fragility function to investigate the tsunami reduction effect of mangrove forest. Under simple numerical conditions related to the mangrove forest, ground level, and incident wave, the model showed that a mangrove forest of Rhizophora sp. with a density of 0.2 trees m⁻² and a stem diameter of 15 cm in a 400 m wide area can reduce the tsunami inundation depth by 30% when the incident wave is assumed to have a 3.0 m inundation depth and a wave period of 30 min at the shoreline. However, 50% of the mangrove forest is destroyed by a 4.5 m tsunami inundation depth, and most of the mangrove forest is destroyed by a tsunami inundation depth greater than 6 m. The reduction effect of tsunami inundation depth decreased when the tsunami inundation depth exceeded 3 m, and was mostly lost when the tsunami inundation depth exceeded 6 m.     

Forecasting of tsunami wave heights at the Russian coast of the Black Sea

Forecasting of tsunami wave heights at the Russian coast of the Black Sea is discussed. Prognostic numerical calculations of tsunamis were carried out for the tsunami sources uniformly distributed in the Black Sea basin (a total of 55 events). Their results are compared with the results of numerical modeling of the historical events (in 1939 and 1966) and the data of not numerous measurements. A preliminary forecast is made on this basis for the tsunami wave heights along the Russian coast of the Black Sea.