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

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

海啸地震震源机制解统计（英文）

近年来,全球海啸灾害频发,给人们带来了巨大损失。为了探索海啸地震发震时的规律,本文统计分析了1976-2010年的美国全球CMT地震矩心矩张量数据,并结合美国海洋与大气局的海啸灾害数据库,从海啸事件和较大海啸事件两个方面对震源机制解展开了对比统计分析,给出了海啸地震在震源方面的统计特征。     

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

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

新西兰海啸沉积研究进展

新西兰长期受到海啸影响，历史和地质学记录都表明新西兰受到多次海啸洪水的影响，新西兰对海啸沉积的研究历史只有十多年，但发展很快，成为国际海啸沉积研究的热点区域。我们简要地小结了新西兰海啸沉积的分布、诊断特征以及海啸沉积研究方法的新进展，并对海啸沉积与风暴沉积的区别进行了总结。     

海啸地震震源机制解统计

为探索海啸地震发震时的规律,统计分析了1976-2010年的美国全球CMT地震矩心矩张量数据,并结合美国海洋与大气局的海啸灾害数据库,从海啸事件和较大海啸事件两个方面对震源机制解展开了对比统计分析,给出了海啸地震在震源方面的统计特征。     

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

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

基于海啸单位源数据库的南中国海区域海啸数值预报方法

基于全球俯冲带板片模型Slab 2.0将马尼拉海沟划分为80个海啸单位源,构建了海啸单位源数据库,建立了基于单位源数据库方法预报海啸的业务化流程。经过与现有海啸数值模型计算结果的比对,单位源数据库预报的最大海啸波幅平均预报一致性可以达到88%,能够满足业务化需求。     

基于完整目录数据的全球海啸时空分异研究

海啸是自然灾害中对人类生命财产安全有严重威胁的灾难之一。随着全球气候变化和全球化贸易日益增强,越来越多的人口和经济暴露于海啸灾害。历史海啸灾害的时空分异分析可以帮助我们认识海啸灾害的演变规律,为灾害预警、灾害防控等提供有益参考。本文通过提取具有完整性和同质性的数据（爬高高度（RH））进行全球海啸的时空分异规律研究,结果表明:（1）对于0.1 m≤RH<0.5 m、0.5 m≤RH<1 m、1 m≤RH<5 m、5 m≤RH<10 m、10 m≤RH<20 m和20 m≤RH的间隔,海啸目录分别自1963年、1940年、1950年、1946年、1922年和1885年以来可以被认为是完整的;（2）全球海啸发生有一定的增加趋势,大约每年会多观测到7次波浪爬高事件。在0.1 m≤RH<5 m区间内,海啸发生呈现一定的周期性。当RH大于5 m时,表现出明显的增加趋势;（3）西北太平洋区域、南太平洋区域、东南太平洋区域、印度洋区域海啸发生有一定的增加趋势,而在北美区域则呈减少趋势,东北大西洋区域无显著变化;（4）除北美区域外,其他区域的海啸发生遵循一定的自组织临界行为,相比来说,东北大西洋区域更容易发生小的海啸事件,而西北太平洋区域和印度洋区域更容易发生各种强度的海啸事件。     

基于全球海洋模型GOCTM的海啸传播过程数值模拟——以2004印尼大海啸为例

文中利用在有限体积近海模型FVCOM基础上拓展开发的全球海洋模型GOCTM(Global Ocean Circulation and Tide Model)进行了海啸波传播过程模拟,GOCTM采用全球无结构三角形网格,避免了开边界条件引入带来的误差,利用德国AWI研究所提供的海啸源作为初始水位场,模拟了2004年12月26日苏门答腊-安达曼Mw 9.2地震引发的海啸传播过程。通过模拟结果与印度沿岸潮位站数据以及海啸发生过程观测到的卫星高度计数据进行了对比,发现模拟结果与观测值相近,相关系数最高达0.82,相关性较好。模拟的海啸波到达苏门答腊岛北部的时间与日本的TUNAMI 模型和德国AWI研究所的TsunAWI模型的模拟结果相符,时间相差不到30 min,证明GOCTM全球模型可以较好地对海啸传播过程进行模拟,结果令人满意,希望本工作可以为我国海啸预报和预警提供参考。     

20世纪夏威夷最大海啸的地质证据

夏威夷历史海啸记录包含远距离和当地产生的海啸,由这些事件引起的人员伤亡和对建筑、环境造成的破坏与损失是惊人的,但是在近来才有一些对历史海啸沉积事件的综合性的研究.夏威夷多岩石的海岸线的演变主要是海啸和风暴的贡献,然而区分风暴和海啸是很难的,特别是相关的小的历史海啸事件.因此夏威夷海啸的研究往往关注于大的古海啸沉积事件.研究表明远距离产生的事件,如1946年阿留申海啸对瓦胡岛岩石岸线的演变具有重要影响.活火山也有助于岩石岸线的形成,且容易引发当地海啸.现代活动的Kilauea火山南侧海域的逐步沉降并没有引起太平洋范围内的破坏性的大海啸,但引发了局部一系列重要的海啸事件,这些事件的周期被估计认为是200年.     

中国的海啸灾害危险性及海啸预警系统——解读2011年太平洋海啸演习

新世纪以来频发的海啸灾害引起了国际社会的广泛关注,各滨海国家不仅加大了在海啸预警建设和海啸灾害危险性、海啸基础理论研究方面的投入,更进一步加强了国际社会在应对海啸灾害方面交流与合作.为测试太平洋各国海啸预警系统的有效性以及政府职能部门的应急管理能力,切实提高各国区域和局地海啸的应对水平,促进国家和地区间海啸预警的交流与合作,联合国教科文组织政府间海洋学委员会(IOC/UNESCO)决定2011年11月9-10日,在整个太平洋地区举行一次代号为“Exercise Pacific Wave 11”的海啸演习.此次演习恰逢日本“3.11”大地震海啸发生后的8个月,通过本次演习对督促太平洋各国进一步检验本国的海啸预警系统、评估本国的海啸危险性均有着重要的意义.中国作为IOC和太平洋海啸预警系统的成员国,积极组织实施了我国历史上第一次涉及当地人员疏散的海啸演习.本文将在本次演习所涉及的海啸源评估、海啸数值计算的基础上,应用新的海啸灾害分级标准对我国沿海的海啸危险性和海啸预警系统进行重新评估分析,期望本文的研究将为今后的海啸预警及海啸灾害评估工作提供科学的决策依据标准.     

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.     

南海定量海啸预警系统

我国面临着来自于马尼拉海沟、台湾岛以及琉球海沟等海域可能形成的局地和区域海啸风险,及时、准确地预警海啸是一个难题。2004年印度洋海啸之后,基于海量数据库和GIS技术,国家海洋环境预报中心开发了南海定量海啸预警系统,可以对潜在海啸进行快速定量化的预警,并利用GIS软件进行预警结果可视化。利用该系统对2006年台湾南部地震海啸进行模拟预报,预报结果与实测结果基本吻合。     

Deep water signature of a tsunami

The dynamics of tsunamis can be divided, for convenience, into three parts: tsunami generation, tsunami propagation, and the coastal problems. Out of these three, the problem of tsunami propagation is probably better understood than the other two. One of the main hindrances to the quantitative prediction of tsunami amplitudes at various coastal locations is a lack of detailed knowledge about the deep water signature of a tsunami. Here, the present understanding of this problem is discussed.     

On local probability of invasive tsunami

Abstract

Tsunami occurrence and invasive tsunami at a local area in the circum‐Pacific seismic zone were studied as a Poisson process. The tsunami height at Osaka, Japan, was related to tsunami magnitude. The exceedence frequency of invaded tsunami at Osaka showed a good fitness to the Poisson process. However, an adapted process should be introduced for exceedence frequency of tsunami occurrence in the western Pacific. The exceedence probability of invasive tsunami at Osaka was shown as a function of tsunami magnitude on a diagram with a parameter of the time period. The obtained result might be useful for the planning of coastal area, warning of invasive tsunami, and designing coastal structures as protection within the scope of tsunami economics.     

The tsunami of 3 February 1605

Whether the origins of the tsunami of 3 February 1605 were separated ones or a joined one, is one of the most important problems for Japanese society. When the marine knowledge has not been popularized, it has been considered that the damage pattern on Hachijo Island (including Kojima) was an important key to solving this question.It resembles the situation that the reports concerning the tsunami of 1 April 1946 were introduced to Japanese society under an internal disturbance just after the war, and we overlooked the most important consideration on the combination of tsunami and storm waves caused by the trade wind, and received this tsunami as an extraordinarily huge one.We studied the statistics of recent wind directions around Hachijo Island in January and February 1973, 1974, 1977, 1983, 1984, 1985, 1986 and 1987. And we believed that the probability was high that we could explain well that the damage pattern on Hachijo Island in 1605 was due to the combination of tsunami and storm waves caused by the monsoon.We need not, therefore, look back the damage pattern on Hachijo Island when we consider the fact that the necrologies of many temples in Shizuoka Prefecture record no death on 3 February 1605, proves to be the key to infer that the origins were separated ones.     

A new objective tsunami magnitude scale

Abstract

Tsunami research has suffered considerably due to the lack of an objective magnitude scale. The Imamura‐Iida Grade Scale, although it has been used until now as a magnitude scale, is not objective, and it is more appropriate to call it a tsunami intensity scale. The tsunami magnitude scale proposed by Abe (1979) is also based on coastal tide gauge records, and this is not truly an objective magnitude scale. A new magnitude scale is proposed here, based on total tsunami energy, and it is demonstrated that this new scale will provide an adequate representation of the whole spectrum of tsunamis, starting from the negligible ones and including those that have devastated whole coastlines.     

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.     

美国夏威夷州海啸灾害研究概述

夏威夷是世界上最易遭受海啸洪水灾害的地区之一。本文介绍了美国夏威夷群岛在海啸实时监测和数值模拟方面的研究进展,并重点介绍了夏威夷海啸预警系统的组成与作用。     

Coastal forest effects on tsunami run-up heights

An experimental study was carried out to determine the effects of a coastal forest on tsunami run-up heights. The beach was built as a natural sandy beach at laboratory scale. The coastal forest model was constructed using artificial trees (FM–I) and cylindrical timber sticks (FM–II). Artificial trees were placed on a 1:5 slope in three different layouts: rectilinear, staggered, and dense rectilinear. It was shown that in the case when the trees were placed in the dense rectilinear pattern and close to the still water level (SWL), the run-up height was reduced by approximately 45% compared with the case without trees. After evaluation of the experimental results, the parameters that affect the run-up height were determined. These parameters were written as a dimensionless group using Buckingham's Pi theorem. An extensive regression analysis was carried out and equations proposed. Furthermore, all experiments were repeated with a slope of 1:3.5 to verify the proposed equations. The experimental results were compared with the results of the proposed equations, and it was shown a good agreement between the results.