马尼拉海沟地震海啸对中国大陆的影响
石峰1, 何宏林1, 周本刚1, 魏占玉1, 毕丽思2
1中国地震局地质研究所, 活动构造与火山重点实验室, 北京 100029
2广东省地震局, 广州 510070

〔作者简介〕 石峰, 男, 1984年生, 2007年本科毕业于北京大学, 2014年于中国地质局地质研究所获博士学位, 助理研究员, 主要从事活动构造和海啸方面的研究, E-mail: Shifeng@ies.ac.cn

摘要

分析可能影响中国的地震海啸发生的构造环境, 查明中国大陆沿海历史上有没有遭受过海啸的袭击, 对于预测和预防将来可能发生的海啸灾害具有重要意义。中国大陆以东受宽阔的大陆架和一系列岛弧保护, 只有南海马尼拉海沟俯冲带具备发生可能引发海啸的逆冲型地震的条件。文中在分析史料和马尼拉海沟俯冲带构造环境的基础上, 利用数值模拟技术, 认为对中国沿海可以产生最大达4.0m浪高的海啸威胁。历史文献记录也支持这一结论。这些结果都表明, 马尼拉海沟潜在地震海啸对中国大陆影响很大, 值得我们重视和预防。

关键词: 马尼拉海沟; 数值模拟; 海啸
中图分类号:P317.9 文献标志码:A 文章编号:0253-4967(2018)03-579-11
EFFECT OF TSUNAMIS GENERATED IN THE MANILA TRENCH ON CHINA MAINLAND
SHI Feng1, HE Hong-lin1, ZHOU Ben-gang1, WEI Zhan-yu1, BI Li-si2
1)Key Laboratory of Active Tectonics and Volcano, Institute of Geology, China Earthquake Administration, Beijing 100029, China
2)Earthquake Administration of Guangdong Province, Guangzhou 510070, China
Abstract

Tsunami is one of the most devastating natural coastal disasters. Most of large tsunamis are generated by submarine earthquakes occurring in subduction zones. Tsunamis can also be triggered by volcano eruptions and large landslides. There are many records about “sea-overflow” in Chinese ancient books, which are not proved to be tsunamis. Tectonics and historical records analysis are import to forecast and prevention of tsunami. Consider the tectonic environment of the China sea, the possibility of huge damage caused by the offshore tsunami is very small. And the impact of the ocean tsunami on the Bohai sea, the Yellow sea, and the East China sea is also small. But in the South China Sea, the Manila subduction zone has been identified as a high hazardous tsunamigenic earthquake source region. No earthquake larger than MW7.6 has been recorded in the past 100a in this region, suggesting a high probability for larger earthquakes in the future. If a tsunamigenic earthquake were to occur in this region in the near future, a tragedy with the magnitude similar to the 2004 Indian Ocean tsunami could repeat itself. In this paper, based on tectonics and historical records analysis, we have demonstrated that potential for a strong future earthquake along the Manila subduction zone is real. Using a numerical model, we have also shown that most countries in the South China Sea will be affected by the tsunamis generated by the future earthquake. For China, it implies that the maximum wave height over 4.0 meter on China mainland, especially the Pearl River Estuary. But the island, local relief maybe influence the maximum wave. But it takes nearly 3 hours to attack China mainland, if there is the operational tsunami warning system in place in this region, should be greatly reduced losses. And the simulated results are conformable to historical records. It indicates that the tsunami hazards from Manila trench to China mainland worthy of our attention and prevention.

Keyword: Manila Trench; Numerical simulation; Tsunami
0 引言

海啸(tsunami)是一种巨大的海浪, 是以重力为恢复力的重力波。海底大规模的、 突然的上下变动, 包括海底火山喷发, 海底或海岸滑坡、 崩塌、 滑塌, 陨星或彗星的撞击以及海底地震都会激发海啸(Lay et al., 1995; Satake, 2000)。但是, 在激发海啸的诸多原因中, 最主要的原因还是海底地震, 特别是以沿着断层面上下错动为其特征的 “ 倾滑型” 地震(陈运泰等, 2005)。在20世纪, 全球发生了498次海啸, 有66次引起了灾难。最终识别出海啸分别起源于地震事件(86%)、 火山活动(5%)、 滑坡(4%), 以及这些事件的综合过程(5%)(Rhodes et al., 2006)。地震海啸(earthquake-generated tsunami)系指地震激发(产生)的海啸。通过对海啸特点的分析, 可知影响地震激发海啸的主要因素有: 1)地震的大小(以地震矩M0或矩震级MW量度); 2)地震机制; 3)震源深度; 4)震源破裂过程(陈运泰等, 2005); 此外还与区域的海水深度有很大关系。2004年印度尼西亚苏门塔腊大海啸和2011年的日本东北大海啸, 造成了巨大的伤亡和经济损失, 引起了各国对防范海啸灾害的重视。中国拥有漫长的海岸线, 并且占全国GDP绝大比例的发达地区都集中在沿海, 因此, 分析可能影响中国的地震海啸发生的构造环境, 查明中国大陆沿海历史上或史前有没有遭受过海啸的袭击, 对于预测和预防将来可能发生的海啸灾害具有重要意义。

中国大陆以东受宽阔的大陆架和一系列岛弧保护, 包括日本列岛、 朝鲜半岛、 琉球群岛、 台湾岛、 菲律宾群岛、 印度尼西亚群岛、 中南半岛等, 太平洋的远洋海啸波传播到中国近海, 经过岛弧和浅海大陆架, 能量极大衰减; 这种天然地理屏障, 可以阻挡太平洋传播来的远洋海啸。例如1960年智利大海啸, 对菲律宾和日本都造成了灾害, 但传到中国沿岸, 香港记录到的浪高才38cm, 长江口记录的浪高也仅20cm, 不至于产生破坏。

此外, 中国4大海域, 渤海和黄海的平均水深都< 50m, 东海也只有370m, 都不具备发生破坏性地震海啸的条件。但是南海的平均水深≥ 1i000m, 而且在南海有马尼拉海沟俯冲带, 这里具备发生可能引发海啸的逆冲型地震的条件和传播海啸的条件(图1)。因此, 本文在马尼拉海沟构造环境分析的基础上, 结合潜在地震的海啸模拟来讨论未来可能对中国大陆产生影响的地震海啸。

图 1 中国海域地形图
黑点为海啸记录点, 数据来自Lau et al., 2010; 红心圆为马尼拉海沟附近7级以上仪器地震, 数据来自CSNDMC; 黄心圆为马尼拉海沟附近7级以上历史地震, 数据来自PHIVOLCS; 绿心圆为马尼拉海沟附近7级以上历史地震, 数据来自Soloviev et al., 1974
Fig. 1 Map showing topography of the floor of Chinese seas.

1 南海历史海啸分析

中国是世界上最早记载地震的国家, 也是世界上最早记录了海啸的国家。几千年来中国记载与海啸或 “ 海溢” 相关的事件超过100次, 几乎在中国沿海从南到北都有记录, 但大部分集中于南海沿岸(王晓青等, 2006)。其中大部分并非现代意义上的地震海啸, 而是台风等气候原因引起的风暴潮(王锋等, 2005; 石峰等, 2012)。例如, 中国南海最早的海啸记录为1076年10月31日, 于广东一带观测到异常高的海浪。但是对于产生的原因, 大部分研究者认为可能是台风(Keimatsu, 1963; Cox, 1970; 李善邦, 1981)。1604年12月29日福建泉州以东海域发生8级大地震, 这是有史料记载以来台湾海峡最强烈的1次地震, 也是中国东南沿海地区最大的地震。这次地震使福建沿海尤其是古城泉州遭受了严重破坏, 震感范围可达中国东部的江苏、 上海、 安徽、 湖北、 广西等十余个省市(黄昭等, 2006)。截至目前, 尚未发现一处文献提到这次地震近海或海岸有明显的海啸异常现象(吴维汕等, 1988)。但是也有比较一致认为是海啸事件的记录。例如1640年侵扰澄海、 揭阳等地的广东海啸(Keimatsu, 1963; Berning-Hausen, 1969; Cox, 1970; 李善邦, 1981)和1781年的台湾海啸(Soloviev et al., 1974; Li et al., 2015)。1867年12月8日, 台湾地区狂浪倾覆船只, 又越过堤岸, 冲入街市, 冲毁屋宇, 浪高5m左右, 也是中国比较典型的1次地震海啸(雷土成等, 1991)。

此外调查发现在广东汕尾市海丰县梅陇镇流传着 “ 3个涌盖梅陇” 的传说, 发生在 1718 年(康熙五十七年)。《清乾隆志· 卷十· 邑事》记载: “ 康熙五十七年夏五月廿八日, 飓风滛雨。半夜海水泛溢, 浪高数丈, 杨安都村落民畜淹没殆尽” 。《惠州府志· 卷十八》中也发现了与之相同时段的灾难记录: “ 康熙五十七年(1718 年)夏五月飓风霪雨, 覆船数百艘, 夜半海水泛溢, 浪高数丈, 濒海居民漂溺无数。知府余毓浩与归善令欧加意赈恤灾民被泽, 龙川雨雹” 。初步判定这是1次造成较大灾难的海啸。

对于距离马尼拉海沟更近的菲律宾, 虽然文献记录时间不及中国悠久。但是菲律宾沿岸也记载了大量的海啸破坏情况, 多达几十次; 所有的记载皆来自于各国的历史文献资料, 且每一笔资料都有2个以上的相关研究文献, 大部分内容包含了发生的时间、 地点、 原因、 海啸的浪高记录、 造成的财产损失和人民伤亡(Soloviev et al., 1974)。

这些关于南海或马尼拉海沟地震与海啸的历史文献记载, 说明了南海地区可能存在潜在的海啸威胁。因此, 更不能轻视马尼拉海沟地震而引发海啸的可能性。

2 马尼拉海沟构造环境分析

南海处于欧亚板块、 太平洋板块(菲律宾板块)和印度洋板块的交会地带。3大板块的相互作用, 使南海处在1个复杂的构造应力环境中, 并且形成了性质各异的四周边界。在这些边界中尤其以东部边界— — 马尼拉海沟俯冲带最为复杂。马尼拉海沟位于吕宋岛西侧, 北至北纬21.5° , 南至南纬13° , 长990km。在这里欧亚板块向E俯冲至菲律宾海板块下方而形成马尼拉海沟。

马尼拉海沟的形成与南海的形成和俯冲密切相关。在6MaiBP时南海古扩张洋脊已经接近马尼拉海沟, 但并未发生碰撞; 在台湾岛和其南部地区欧亚板块和北吕宋岛弧的碰撞也即将开始, 而此时欧亚板块沿马尼拉海沟的俯冲则形成了西火山岛链。在4~5MaiBP时, 南海古扩张脊已经与海沟接触并开始俯冲, 而且形成台湾岛的欧亚板块和北吕宋岛弧的碰撞也在此阶段发生, 但是由于古扩张洋脊对俯冲的阻碍和形成台湾岛的碰撞作用, 使得西火山链北部的火山停止喷发, 而已经俯冲的板块所具有的重力和动能, 引起板块在深部的撕裂。2MaiBP左右, 由于古洋脊阻碍作用而停止的俯冲重新开始, 形成东火山岛链。俯冲的古扩张洋脊会造成俯冲板块倾角的改变, 并因此引发南海洋陆过渡带的撕裂, 这个破裂的形成可以解释东火山岛链喷发的岩浆中幔源成分的富集(Yang et al., 1996)。GPS测量数据显示马尼拉俯冲带聚敛速率达80mm/a(Yu et al., 1999), 其中在吕宋岛北部速率最大, 向南略微减少(图2)。马尼拉海沟附近地震分布在17.5° N以南相对密集, 而且150~300km深度的地震在南部更靠近海沟, 俯冲震源深度剖面资料显示, 马尼拉海沟附近俯冲板块倾角自北向南加深, 南端倾角可以达到30° 以上(图2)。

图 2 马尼拉海沟及邻区地震构造图Fig. 2 Seismotectonics of the Manila trench and adjacent areas.

马尼拉海沟俯冲带1个最大的特点是, 自北向南俯冲板块分别为陆壳、 洋壳、 古扩张洋脊、 洋壳, 其差异性使得俯冲带具有明显的分段和火山活动差异。自北向南大体可以分为: ①巴士海峡段, 该段俯冲板块为陆壳; ②北吕宋海槽段, 该段俯冲板块为洋壳, 发育北吕宋海槽; ③南海古扩张洋脊段, 该段为古洋脊向下俯冲, 没有发育海槽, 该段北部NW向菲律宾断裂与之交会, 并可能使海沟左旋错列; ④南吕宋海槽段, 该段俯冲板块为洋壳, 发育吕宋海槽; 在该段以南又转变为陆壳。俯冲板块最大的不连续发生在古扩张洋脊, 由于俯冲板块的横向撕裂, 使得向S板块俯冲倾角突然增加, 形成 “ 板片窗” (刘再峰等, 2007)。此外, 海沟附近最大主应力方向分布资料也反映了海沟存在自北向南的差异(刘再峰等, 2007)。基于前述地震构造资料分析和前人对马尼拉海沟海啸破裂源的划分(Kirby, 2006; 黄惠绢, 2008), 我们将马尼拉海沟破裂源分为5段, 其中破裂段1位于巴士海峡段, 破裂段3和破裂段4的分界点就是南海古扩张洋脊。由于古扩张洋脊的构造作用, 导致马尼拉海沟南北差异很大, 马尼拉海沟不太可能发生全段联合破裂。其次, 破裂段1和破裂段2由于分属于陆壳和洋壳, 也不太可能发生联合破裂; 破裂段2和3虽同属洋壳, 但是倾角相差较大, 也不太可能发生联合破裂; 只有破裂段4和5同属洋壳, 倾角也接近, 有一定的可能发生联合破裂。在此基础上, 给出了海沟地震破裂面的划分(图2)。

3 潜在地震的海啸模拟与结果分析
3.1 数值模拟方法简介

海啸事件发生次数少, 且不易观测, 因此研究海啸的1个重要方法是数值模拟。数值模拟可以重现海啸在无法量测及观测之地的波高、 了解海啸传播的动力过程、 探讨海啸波对海岸工程的破坏, 以及研究地震震源的几何特性等。目前常见的海啸数值模拟方法有: 美国南加利福尼亚州大学开发的MOST(Method of SplittingTsunami)模型(Titov et al., 1998), 美国康乃尔大学开发的COMCOT(Cornell Multigrid Coupled Tsunami)模型(Liu et al., 1994)和日本东北大学开发的TIME模型(Imamura, 1996)。

本研究采用康乃尔大学开发的COMCOT模型, 它是由康乃尔大学土木工程学系研发的基于Fortran语言的数值仿真模型, 专门模拟受扰动的海水波浪传递的过程, 适用于长波运动, 并可用于海啸波对海洋沿岸的损害及工程评估, 可以模拟出海啸的波高、 溯上高度以及淹溢范围等。该模型较为成熟, 其准确性已经通过与卫星资料、 潮位资料和实地调查资料比对获得了验证, 如1960年智利海啸(Titov et al., 1998)、 1986 及2002年花莲地震海啸(Wang et al., 2005)、 2003年Algerian 海啸以及2004和2005年印度洋海啸事件(Wang et al., 2006)等。

COMCOT模型的主要特色有以下几个方面: 1)可选用Okada(1985)的断层模型或Mansinha等(1971)的断层模型计算破裂面错动引起的海床垂直位移量; 2)可输入初始波高资料进行波的传播模拟; 3)可使用滑坡与人造波进行海啸数值模拟; 4)可计算线性或非线性的浅水波方程; 5)可选择使用球坐标或卡式坐标系统; 6)可进行多层网格套迭以及计算淹溢范围(Liu et al., 1998)。

3.2 地形数据选择及潜在地震的参数选择原则

地形数据资料来源于美国地球物理资料中心(National Geophysical Data Center, NOAA)下载的ETOPO2资料, 本研究选用的资料格式为“ .xyz” 。

地震海啸的模拟需要精确的震源参数, 图3的断层模型表示了应用COMCOT模型所需要的震源参数。由于是对马尼拉海沟潜在地震的模拟, 所以我们应用以下几个原则进行处理:

图 3 断层参数示意图Fig. 3 Sketch of fault parameters.

(1)震级的判定主要依据基于全球俯冲带倾滑数据给出的震级与破裂长度统计关系式计算获得, 该关系式在海啸震源参数评价中得到了广泛的应用(Papazachos et al., 2004); 同时根据破裂段上的最大历史地震进行修正。

(2)震源参数通过以下关系式得到(Papazachos et al., 2004)。

震级与破裂长度关系式:

LogL=0.55M-2.19, σ=0.18, 6.7M9.3(1)

震级与破裂面宽度关系式:

LogW=0.31M-0.63, 6.7M9.2(2)

震级与平均滑动量关系式:

Logu=0.64M-2.78, 6.7M9.2(3)

(3)破裂面倾角和滑动角主要依据俯冲带50km以浅的震源深度剖面估计, 在没有深度剖面时, 根据相邻段倾角推断获得。

(4)根据张国民等(2002)统计的1970年1月至2000年5月中国东部ML≥ 2.0地震的震源深度绝大部分位于5~24km, 而马尼拉海沟附近地震也大都集中在5~25km(图2), 所以本文采用其平均值15km计算。

根据上述原则, 给出了不同破裂源潜在的最大地震的相关参数(表1)。

表1 马尼拉海沟强震震源断层的相关参数及其模拟浪高 Table1 Parameters of source faults and simulated height of waves in the Manila trench
3.3 模拟结果

当破裂段1(M1)完全破裂时, 地震的矩震级约为8.1, 造成的海啸最大浪高达到了2.2m。约3h到达中国大陆沿海, 到岸时的最大浪高约为0.8m, 主要集中在广东沿岸的深圳— 汕头一带(图4a)。这样的浪高会对最沿海有一定的破坏。

图 4 马尼拉海沟潜在地震的海啸数值模拟结果
左为海啸走时图; 中为最大浪高; 右为局部最大浪高
a M1; b M2; c M3; d M4; e M5; f M4+5
Fig. 4 Tsunami simulation of possible Manila trench earthquakes.

当破裂段2(M2)完全破裂时, 地震的矩震级约为8.6, 造成的海啸最大浪高达到了6.3m。约3h到达中国大陆沿海, 到岸时的最大浪高约4m, 几乎整个广东沿岸和福建沿岸的北部浪高都达到了1m左右(图4b)。这样的浪高会对最沿海产生一定程度的破坏。由于需要3h左右才能到达中国大陆沿岸, 而在20min以内就可到达菲律宾西北部沿岸, 所以可以根据地震的发生和菲律宾沿岸的破坏来进行有效的海啸预警。根据模拟结果这次地震时震级最大也是破坏最大的1次, 而且到达中国大陆沿海也只需要不到3h, 所以对破裂段2的地震要尤其注意。

当破裂段3(M3)完全破裂时, 地震的矩震级约为8.4, 造成的海啸最大浪高达到了7.7m。约3h到达中国广东和福建沿海, 到岸时的最大浪高约0.9m; 约3~4h到达浙江一带, 到岸浪高约0.5m(图4c)。这样的浪高会对最沿海产生一定的破坏。

当破裂段4(M4)完全破裂时, 地震的矩震级约8.0, 造成的海啸最大浪高达到了1.7m。约3h到达中国大陆沿海, 到岸时的最大浪高约0.5m(图4d)。这样的浪高会对沿海产生一定的影响。

当破裂段5(M5)完全破裂时, 地震的矩震级约7.9, 造成的海啸最大浪高达到了2.1m。约3~4h到达中国大陆沿海, 到岸时的最大浪高< 0.3m(图4e)。这样的浪高几乎不会对沿海造成破坏。

当破裂段4和5(M4+5)完全破裂时, 地震的矩震级约为8.5, 造成的海啸最大浪高达到了4.5m。约3~4h到达中国大陆沿海, 到岸时的最大浪高约1.0m(图4f)。这样的浪高会对最沿海产生一定的破坏。

根据模拟结果, 除了破裂段5(M5)以外其他破裂段完全破裂时产生的海啸都会对中国沿海产生一定的破坏。尤其是破裂段2(M2)完全破裂时, 可能造成的破坏最大, 这可能与破裂段长度、 倾角与滑动量等参数相关。最大到岸浪高到达了4.0m, 主要集中在珠江口附近。考虑到近岸区域地形变化复杂, 海岛密布, 局部地形条件可能会很大地影响海啸浪高。历史文献中记录的1718年的海啸浪高数丈也在这个地方附近, 说明模拟结果显示的受灾最严重的地方与历史文献记录较一致。至于历史文献记录的浪高达到数丈也就是十几m, 考虑到古文记录中存在一定的夸大和局部地形条件可能会很大地影响模拟海啸浪高, 所以具有一定的合理性。而且造成的海啸到达中国海岸的时间只有3~4h, 给人们的反应时间也因此而短得多。但是如果产生海啸, 会在约20min内先到达菲律宾, 而且在菲律宾的浪高普遍是中国海岸的2倍以上, 所以如果能及时得到菲律宾西北海岸的破坏程度和浪高, 可以很好地指示我们对中国沿海海啸的预警和防护。

当然, 如果马尼拉海沟全段完全破裂, 会造成的更大的影响, 可能在广东沿海产生10m高的海啸抵达波(黄慧绢, 2008)。但是, 由于马尼拉海沟的构造特征, 这种可能性几乎不存在。

4 结论

通过对中国附近海域的构造环境分析, 认为近海海啸造成破坏的可能性很小。而远洋海啸对渤海、 黄海和东海的影响也较小。最有可能造成破坏性海啸的是马尼拉海沟俯冲带对南海的影响。在对马尼拉海沟俯冲带构造环境分析的基础上, 利用数值模拟技术, 认为对中国沿海可以产生最大达4.0m浪高的海啸威胁, 主要集中在珠江口附近。考虑到近岸区域地形变化复杂, 海岛密布, 局部地形条件可能会很大地影响海啸浪高。文献记录也支持这一结论。这些结果都表明, 马尼拉海沟潜在地震对中国大陆的影响很大, 值得我们重视和预防。

The authors have declared that no competing interests exist.

参考文献
[1] 陈运泰, 杨智娴, 许力生. 2005. 海啸、 地震海啸与海啸地震[J]. 物理, 34(12): 864872.
CHEN Yun-tai, YANG Zhi-xian, XU Li-sheng. 2005. Tsunamis, earthquake-generated tsunamis and tsunamigenic earthquakes[J]. Physics, 34(12): 864872(in Chinese). [本文引用:2]
[2] 黃惠絹. 2008. 馬尼拉海溝地震引發海嘯的潛勢分析 [D]. 中国台湾: 國立中央大學.
HANG Hui-juan. 2008. Analysis of the potential tsunami generated by the earthquakes along the Manila sub-duction zone[D]. National Central University, Taiwan, China(in Chinese). [本文引用:2]
[3] 黄昭, 王善雄, 王喜年. 2006. 1604年泉州海外大地震及其海啸影响分析[J]. 地震, 26(4): 94102.
HUANG Zhao, WANG Shan-xiong, WANG Xi-nian. 2006. Analysis of the 1604 Quanzhou offshore earthquake-tsunami effects[J]. Earthquake, 26(4): 94102(in Chinese). [本文引用:1]
[4] 雷土成, 欧秉松. 1991. 台湾海峡及其邻近地区的地震海啸与海溢[J]. 台湾海峡, 10(3): 264270.
LEI Tu-cheng, OU Bing-song. 1991. Seaquake and sea overflow in Taiwan Strait and its adjacent regions[J]. Journal of Oceanography in Taiwan Strait, 10(3): 264270(in Chinese). [本文引用:1]
[5] 李善邦. 1981. 中国地震 [M]. 北京: 地震出版社.
LI Shan-bang. 1981. Earthquake in China [M]. Seismological Press, Beijing(in Chinese). [本文引用:2]
[6] 刘再峰, 詹文欢, 张志强. 2007. 台湾-吕宋岛双火山弧的构造意义[J]. 大地构造与成矿学, 31(2): 145150.
LIU Zai-feng, ZHAN Wen-huan, ZHANG Zhi-qiang. 2007. Tectonic implications of the double island arc between Taiwan and Luzon[J]. Geotectonica et Metallogenia, 31(2): 145150(in Chinese). [本文引用:2]
[7] 石峰, 毕丽思, 谭锡斌, . 2012. 渤海海域历史上发生过地震诱发海啸吗?[J]. 地球物理学报, 55(9): 30973104.
SHI Feng, BI Li-si, TAN Xi-bin, et al. 2012. Did earthquake tsunami occur in Bohai Sea in history?[J]. Chinese Journal of Geophysics, 55(9): 30973104(in Chinese). [本文引用:1]
[8] 王锋, 刘昌森, 章振铨. 2005. 中国古籍中的地震海啸记录[J]. 中国地震, 21(3): 437443.
WANG Feng, LIU Chang-sen, ZHANG Zhen-quan. 2005. The earthquake-generated tsunami records in Chinese ancient books[J]. Earthquake Research in China, 21(3): 437443(in Chinese). [本文引用:1]
[9] 王晓青, 吕金霞, 丁香. 2006. 中国地震海啸危险性初步探讨[J]. 华南地震, 26(1): 7680.
WANG Xiao-qing, Jin-xia, DING Xiang. 2006. A preliminary study on the risk of tsunami in China[J]. South China Journal of Seismology, 26(1): 7680(in Chinese). [本文引用:1]
[10] 吴维灿, 许东亚, 丁祥焕, . 1988. 1604年福建泉州8级地震 [M]∥国家地震局. 中国特大地震研究. 北京: 地震出版社: 5269.
WU Wei-chan, XU Dong-ya, DING Xiang-huan, et al. 1988. 1604 Fujian Quanzhou M8 earthquake [M]∥State Seismological Bureau. China Great Earthquake Research. Seismological Press, Beijing: 5269(in Chinese). [本文引用:1]
[11] 张国民, 汪素云, 李丽, . 2002. 中国大陆地震震源深度及其构造含义[J]. 科学通报, 47(9): 663668.
ZHANG Guo-min, WANG Su-yun, LI Li, et al. 2002. Focal depth research of earthquakes in Mainland China: Implication for tectonics[J]. Chinese Science Bulletin, 47(12): 969974. [本文引用:1]
[12] Berning-Hausen W H. 1969. Tsunamis and seismic seiches of Southeast Asia[J]. Bulletin of the Seismological Society of America, 59(1): 289297. [本文引用:1]
[13] Cox D C. 1970. Discussion of “Tsunamis and seismic seiches of Southeast Asia” by William H. Berninghausen[J]. Bulletin of the Seismological Society of America, 60(1): 281287. [本文引用:2]
[14] Imamura F. 1996. Review of tsunami simulation with a finite-difference method [M]∥Yeh H, Liu P, Synolakis C. Long-Wave Runup Models. World Scientific Publishing Co, Hackensack. [本文引用:1]
[15] Keimatsu M. 1963. On the historical tidal waves in China[J]. Zishin, 16: 149160(in Japanese). [本文引用:2]
[16] Kirby S. 2006. Tsunami source characterization for western Pacific subduction zone [R]. USGS Tsunami Source Working Group. [本文引用:1]
[17] Lau A Y A, Switzer A D, Dominey-Howes D, et al. 2010. Written records of historical tsunamis in the northeastern South China Sea-Challenges associated with developing a new integrated database[J]. Natural Hazards and Earth System Science, 10(9): 17931806. [本文引用:1]
[18] Lay T, Wallace T C. 1995. Modern Global Seismology [M]. San Diego: Academic Press. [本文引用:1]
[19] Li L L, Switzer A D, Wang Y, et al. 2015. What caused the mysterious eighteenth century tsunami that struck the southwest Taiwan coast?[J]. Geophysical Research Letters, 42(20): 84988506. [本文引用:1]
[20] Liu P, Cho Y, Yoon S et al. Yoon S 1994. Numerical simulations of the 1960 Chilean tsunami propagation and inundation at Hilo, Hawaii[M]∥El-Sabh M I. Recent Development in Tsunami Research. Kluwer Academic, Dordrecht. [本文引用:1]
[21] Liu P L F, Woo S B, Cho Y S. 1998. Computer programs for tsunami propagation and inundation [R]. Cornell University. [本文引用:1]
[22] Mansinha L, Smylie D E. 1971. The displacement fields of inclined faults[J]. Bulletin of the Seismological Society of America, 61(5): 14331440. [本文引用:1]
[23] Okada Y. 1985. Surface deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 75(4): 11351154. [本文引用:1]
[24] Papazachos B C, Scordilis E M, Panagiotopoulos D G, et al. 2004. Global relations between seismic fault parameters and moment magnitude of earthquakes[J]. Bulletin of the Geological Society of Greece, 36: 14821489. [本文引用:2]
[25] Rhodes B, Tuttle M, Horton B, et al. 2006. Paleotsunami research[J]. Eos, Transactions American Geophysical Union, 87(21): 205209. [本文引用:1]
[26] Satake K. 2000. Tsunami [M]∥Lee W H K, Jennings P, Kisslinger C, et al. International Hand book of Earthquake and Engineering Seismology, Part A. Academic Press, San Diego. [本文引用:1]
[27] Soloviev S L, Go C N. 1974. A Catalogue of Tsunami on the Western Shore of the Pacific Ocean[M]. Translated by Canada Institute for Scientific and Technical Information, National Research Council Ottawa, Nauka House, Moscow. [本文引用:2]
[28] Titov V V, Synolakis C E. 1998. Numerical modeling of tidal wave runup[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4): 157171. [本文引用:2]
[29] Wang X, Liu P L F. 2005. Preliminary simulation of 1986 & 2002 Taiwan Hualien Tsunami [R]. Cornell University. [本文引用:1]
[30] Wang X M, Liu P L F. 2006. An analysis of 2004 Sumatra earthquake fault plane mechanisms and Indian Ocean tsunami[J]. Journal of Hydraulic Research, 44(2): 147154. [本文引用:1]
[31] Yang T F, Lee T, Chen C H, et al. 1996. A double island arc between Taiwan and Luzon: consequence of ridge subduction[J]. Tectonophysics, 258(1-4): 85101. [本文引用:1]
[32] Yu S B, Kuo L C, Punongbayan R S, et al. 1999. GPS observation of crustal deformation in the Taiwan-Luzon region[J]. Geophysical Research Letter, 26(7): 923926. [本文引用:1]