东海海域潜在地震海啸的数值模拟初步研究

针对东海冲绳海槽地区的地震地质背景,对东海海域潜在的地震海啸进行了预研究. 假设了冲绳海槽在发生8.5级大地震,断层错距高达15 m的极端地震情况引发的海啸对中国东部沿海地区的影响. 初步数值模拟结果表明,该地震引发的海啸的最大初始波高为4.3 m,4小时左右传至浙江沿岸,近岸各处波高为1——2 m,其中局部地区波高为2.4 m;约7——8小时靠近上海海岸线（若震源在中冲绳海槽地区,海啸传到上海最快大约7小时）,近岸波高约为1 m. 近岸区域地形变化复杂,海岛密布,局部地形条件可能会很大地影响实际各地点海啸波高,加上海啸在岸边爬高及港湾效应,估计波高还会升高. 给出了冲绳海槽南、中、北部发生潜在地震海啸的传播等时图. 笔者在东海设置了3个地震及海啸监测站,基于海啸模拟结果绘制了监测站处的海啸随时间演化曲线,分析了预研究成果对海啸预警可能发挥的作用.      

设定琉球海沟发生罕遇地震评估我国东南沿海地区的海啸风险

根据构造相似条件分析,琉球海沟与日本海沟、智利海沟、印尼巽他海沟一样具备发生9级罕遇超巨大地震的可能。在对近几年来全球发生的超巨大地震参数及构造对比分析的基础上,设定琉球海沟9.0级地震参数,并将其引发的海啸进行数值模拟研究。结果表明,该地震可引发初始波高为8m的海啸,台湾东北部半小时后遭受10m以上海啸,3~4小时左右传至浙南、闽北沿岸,近岸各处波高在1~2m;5小时左右传至浙北、粤北沿岸,浙江近岸各处波高在2m左右,广东沿海、台湾海峡由于台湾岛的正面阻挡,海啸波高低于50cm;8小时后靠近上海海岸线,最大波高约1m。海啸的上岸高度与海岸附近的海深和海岸线的形态密切相关,我国东南海域地形变化复杂、海湾众多,对海啸波有放大作用,模拟结果可能比实际海啸偏小。我国沿海地区分布着不少已建和在建的核电厂,在核电设计时未考虑海啸,一旦发生这种罕遇地震海啸则影响不可忽视,尤其是若与风暴潮、天文大潮叠加则可能出现严重后果。由于核电安全要求万无一失,故须制订有效预警和应对措施。     

日本南海海槽发生罕遇地震情况下我国华东沿海的海啸危险性研究

基于日本南海海槽地震活动性和历史海啸事件记载的分析,本文对日本南海海槽发生M_W9.1罕遇地震情况下的海啸进行了数值模拟研究.结果表明:该地震可引发初始波幅约10 m的海啸,6个小时后传至浙江沿海,近岸各处波幅为1—2 m;8个小时后靠近上海海岸线,最大波幅约2 m,受地形影响局地爬高至近3 m;11个小时后抵达苏北黄海沿岸,预计波幅普遍在1 m左右.海啸的上岸高度与海岸附近的海深和海岸线的形态密切相关.我国近岸海域地形变化复杂,海湾众多,对海啸波有放大作用,该模拟结果可能比实际传播到近岸时偏小,因此综合评估日本海啸影响我国华东地区的规模m可达1—2级左右.一旦日本南海发生罕遇地震对我国的影响不容忽视,尤其遇上风暴潮与天文大潮叠加,则可能会造成一定程度的海啸灾害.      

地震海啸在东海大陆架上传播过程的数值模拟研究

在东海潜在震源区冲绳海槽假定了五个震源点,根据Steven地震海啸地震参数经验值作为初始条件,分别考虑6.5、7.0、7.5、8.0、8.5、9.0级地震条件下的30个震例,采用数值模拟的方法,对海啸在东海传播过程进行情境分析,特别是对上海沿岸地区可能会遭受的海啸灾害做了较为精细的研究.结果发现:小于8.0级的震例对上海地区几乎不会造成影响;8.0级震例只有最北端震源点震例会对上海地区有明显影响;8.5级以及9.0级震级基本上均会对上海沿岸地区造成较大的影响.特别是冲绳海槽北段9.0级震例可能会对上海沿岸局部地区造成危害,最大波高可达3.9m.     

冲绳海槽地震海啸的数值试验研究

在分析东海水深与地震情况基础上,我们认为冲绳海槽地震具有引发海啸的条件,并对该区海啸传播过程采用单侧破裂方式的有限移动源模型和高阶Boussinesq方程做了一维有限差分法数值模拟,给出了海啸传播过程波形和速度变化图.我们从模拟结果得到初步结论:冲绳海槽特大地震海啸传播到上海近海在5小时后,近海波高在几十厘米范围内.     

冲绳海槽地区地震潜在海啸对中国东南沿海的影响研究

采用球坐标系下非线性浅水波方程, 研究日本本州M9.0大地震引发的海啸对中国东南沿海的影响, 并计算了冲绳海槽构造带上3个不同段落可能发生潜在地震引发的海啸, 分析这些海啸与日本大海啸的浪高和走时关系. 结果表明, 日本地震海啸模拟结果与日本当地报道及中国东南沿海7个验潮站的报道结果相符. 冲绳海槽构造带中段可能发生的3次不同震级（M7.0, M7.5, M8.0）潜在地震引发的海啸到达中国东南沿海的时间比日本海啸提前约4个小时, 从震源区传播3个多小时即可到达华东沿海部分验潮站. 冲绳海槽M7.5潜在地震海啸在验潮站上计算的波高与日本海啸相当, 中冲绳海槽M8.0潜在地震海啸在大陈站的波高将超过0.9 m, 在坎门站波高将超过1.8 m. 北冲绳海槽的潜在地震海啸威胁主要集中在江苏盐城、 上海一带, 南冲绳海啸主要对台湾东北部和浙江沿海产生威胁. 本文对冲绳海槽构造带上潜在地震引发海啸的模拟结果, 可为中国东南沿海地区的防震减灾、 海啸预警提供有意义的参考.      

日本Mw9.0级地震海啸数值模拟与启示

2011年3月11日13时46分日本东北部海域发生Mw9.0级特大地震,地震诱发了海啸.本文依据USGS的震源机制解,进行了地震海啸的数值模拟,并与美国国家海洋与大气管理局布设的海底压力计记录的水深数据对比分析,结果表明数值模拟结果可信.同时,进一步分析了海啸造成巨大损失的原因,并对未来我国海啸防灾减灾工作给出了几点建...     

Modeling of the runup heights of the Hokkaido-Nansei-Oki tsunami of 12 July 1993

The Hokkaido-Nansei-Oki earthquake (M _w 7.7) of July 12, 1993, is one of the largest tsunamigenic events in the Sea of Japan. The tsunami magnitudeM _t is determined to be 8.1 from the maximum amplitudes of the tsunami recorded on tide gauges. This value is larger thanM _w by 0.4 units. It is suggested that the tsunami potential of the Nansei-Oki earthquake is large forM _w . A number of tsunami runup data are accumulated for a total range of about 1000 km along the coast, and the data are averaged to obtain the local mean heightsH _n for 23 segments in intervals of about 40 km each. The geographic variation ofH _n is approximately explained in terms of the empirical relationship proposed byAbe (1989, 1993). The height prediction from the available earthquake magnitudes ranges from 5.0–8.4 m, which brackets the observed maximum ofH _n , 7.7 m, at Okushiri Island.     

海啸传播模型与数值模拟研究进展

海啸在浅水大陆架的传播问题由于其非线性作用和浅水效应而变得十分复杂,然而目前成熟的海啸传播理论及数值模拟结果在这方面与实际并不一致.本文比较分析了可用来模拟大陆架海啸传播的浅水波模型和数值方法,并提出对我国东海陆架边缘可能发生的近海海啸需要开展数值试验研究.     

The July 15, 2009 Fiordland,New Zealand Tsunami: Real-Time Assessment

On 15 July 2009, a Mw 7.8 earthquake occurred off the New Zealand coast, which by serendipitous coincidence occurred while the International Tsunami Symposium was in session in Novosibirsk, Russia. The earthquake generated a tsunami that propagated across the Tasman Sea and was detected in New Zealand, Australia and as far away as the US West coast. Small boats close to the epicenter were placed in jeopardy, but no significant damage was observed despite a measured run-up height of 2.3 m in one of the Sounds in close proximity to the source (Wilson in GNS Science Report 46:62 2009). Peak-to-trough tsunami heights of 55 cm were measured at Southport, Tasmania and a height of 1 m was measured in Jackson Bay, New Zealand. The International Tsunami Symposium provided an ideal venue for illustration of the value of immediate real-time assessment and provided an opportunity to further validate the real time forecasting capabilities with the scientific community in attendance. A number of agencies with responsibility for tsunami forecast and/or warning, such as the NOAA Center for Tsunami Research, the Pacific Tsunami Warning Center, GNS Science in New Zealand, the Australian Bureau of Meteorology and the European Commission Joint Research Centre were all represented at the meeting and were able to demonstrate the use of state of the art numerical models to assess the tsunami potential and provide warning as appropriate.     

Field survey of the East Java earthquake and tsunami of June 3, 1994

A field survey of the June 3, 1994 East Java earthquake tsunami was conducted within three weeks, and the distributions of the seismic intensities, tsunami heights, and human and house damages were surveyed. The seismic intensities on the south coasts of Java and Bali Islands were small for an earthquake with magnitudeM 7.6. The earthquake caused no land damage. About 40 minutes after the main shock, a huge tsunami attacked the coasts, several villages in East Java Province were damaged severely, and 223 persons perished. At Pancer Village about 70 percent of the houses were swept away and 121 persons were killed by the tsunami. The relationship between tsunami heights and distances from the source shows that the Hatori's tsunami magnitude wasm=3, which seems to be larger for the earthquake magnitude. But we should not consider this an extraordinary event because it was pointed out byHatori (1994) that the magnitudes of tsunamis in the Indonesia-Philippine region generally exceed 1–2 grade larger than those of other regions.     

Field survey report on tsunami disasters caused by the 1993 Southwest Hokkaido earthquake

Detailed field work at Okushiri Island and along the southwest coast of Hokkaido has revealed quantitatively (1) the advancing direction of tsunami on land, (2) the true tsunami height (i.e., height of tsunami, excluding its splashes, as measured from the ground) and (3) the flow velocity of tsunami on land, in heavily damaged areas. When a Japanese wooden house is swept away by tsunami, bolts that tie the house to its concrete foundation resist until the last moment and become bent towards the direction of the house being carried away. The orientations of more than 850 of those bent bolts and iron pipes (all that can be measured, mostly at Okushiri Island) and fell-down direction of about 400 trees clearly display how tsunami behaved on land and caused serious damage at various places. The true tsunami height was estimated by using several indicators, such as broken tree twigs and a window pane. The flow velocity of tsunami on land was determined by estimating the hydrodynamic force exerted on a bent handrail and a bent-down guardrail by the tsunami throughin situ strength tests.Contrary to the wide-spread recognition after the tsunami hazard, our results clearly indicate that only a few residential areas (i.e., Monai, eastern Hamatsumae, and a small portion at northern Aonae, all on Okushiri Island) were hit by a huge tsunami, with true heights reaching 10 m. Southern Aonae was completely swept away by tsunami that came directly from the focal region immediately to the west. The true tsunami height over the western sea wall of southern Aonae was estimated as 3 to 4 m. Northern Aonae also suffered severe damage due to tsunami that invaded from the corner zone of the sand dune (8 m high) and tide embankment at the northern end of the Aonae Harbor. This corner apparently acted as a tsunami amplifier, and tide embankment or breakwater can be quite dangerous when tsunami advances towards the corner it makes with the coast. The nearly complete devastation of Inaho at the northern end of Okushiri Island underscored the danger of tsunami whose propagation direction is parallel to the coast, since such tsunami waves tend to be amplified and tide embankment or breakwater is constructed low towards the coast at many harbors or fishing ports. Tsunami waves mostly of 2 to 4 m in true height swept away Hamatsumae on the southeast site of Okushiri Island where there were no coastal structures. Coastal structures were effective in reducing tsunami hazard at many sites. The maximum flow velocity at northern Aonae was estimated as 10 to 18 m/s (Tsutsumi et al., 1994), and such a high on-land velocity of tsunami near shore is probably due to the rapid shallowing of the deep sea near the epicentral region towards Okushiri Island. If the advancing direction, true height, and flow velocity of tsunami can be predicted by future analyses of tsunami generation and progagation, the analyses will be a powerful tool for future assessment of tsunami disasters, including the identification of blind spots in the tsunami hazard reduction.     

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock     

日本宫崎县石卷平原海啸堆积物调查

2007年4月,为执行科技部科研院所社会公益课题 "中国沿海地区古海啸的定量研究",课题组部分成员访问了日本,并对日本东北地区宫崎县石卷海岸平原的海啸堆积物做了概略性的调查。2011年3月11日,该地区发生了9.0级大地震并引发海啸,造成了巨大的经济损失和人员伤亡。为了使人们对该地区过去的海啸历史有一个概略认识,并了解海啸堆积物的调查方法,对这次调查作一简要报道。这次调查揭示出3次海啸事件,最新一次发生在公元915年十和田火山喷发之前,对应于公元869年的贞观大海啸。     

华东沿海地震海啸模拟

利用数字地震资料及海洋地震地质基础数据, 基于Matlab平台在Windows系统下重新编译Philip liu等人编写的海啸模拟源程序, 将分层静态位移解结果与已建立的海啸传播初始场快速融合, 用交错网格蛙跳差分方法求解笛卡尔坐标系下非线性浅水波方程, 模拟华东沿海地区地震海啸传播过程, 对华东沿海地区可能遭受海啸袭击的危险做出评估, 结果显示华东沿海地区是发生地震海啸影响较弱地区。     

Faulting on the Anninghe fault zone, Southwest China in Late Quaternary and its movement model

The Anninghe fault is one of the significant earthquake-generating fault zones in the Southwest China. Local his-torical record shows that a M≥7 strong earthquake occurred in the year of 1536. On the basis of the detailed air-photographic interpretation and field investigation, we have acquired the following knowledge: 1 The average sinistral strike-slip rate since the Late Pleistocene is about 3~7 mm/a; 2 There is important reverse faulting along the fault zone besides the main left-lateral strike-slip motion, and the shortening rate across the Anninghe fault zone due to the reverse faulting is about 1.7~4.0 mm/a. If the Xianshuihe fault zone is simply partitioned into the Anninghe and Daliangshan faults, we can also get a slip rate of 3~7 mm/a along the Daliangshan fault zone, which is the same as that on the Anninghe fault zone. Moreover, on the basis of our field investigation and the latest knowledge concerning the active tectonics of Tibetan crust, we create a dynamic model for the Anninghe fault zone.     

Surface rupture zone of the 1303 Hongtong M=8 earthquake, Shanxi Province

Based on the latest displacement of Huoshan piedmont fault, Mianshan west-side fault and Taigu fault obtained from the beginning of 1990‘s up to the present, the characteristics of distribution and displacement of surface rupture zone of the 1303 Hongtong M = 8 earthquake, Shanxi Province are synthesized and discussed in the paper. If Taigu fault, Mianshan west-side fault and Huoshan piedmont fault were contemporarily active during the 1303 Hongtong M = 8 earthquake, the surface rupture zone would be 160 km long and could be divided into 3 segments, that is, the 50-km-long Huoshan piedmont fault segment, 35-km-long Mianshan west-side fault segment and 70-km-long Taigu fault segment, respectively. Among them, there exist 4 km and 8 km step regions. The surface rupture zone exhibits right-lateral features. The displacements of northern and central segments are respectively 6～7 m and the southern segment has the maximum displacement of 10 m. The single basin-boundary fault of Shanxi fault-depression system usually corresponds to M ≈ 7 earthquake, while this great earthquake (M = 8) broke through the obstacle between two basins. It shows that the surface rupture scale of great earthquake is changeable.     

The 26 December 2004 Sumatra Tsunami: Analysis of Tide Gauge Data from the World Ocean Part 1. Indian Ocean and South Africa

The M_w = 9.3 megathrust earthquake of December 26, 2004 off the northwest coast of Sumatra in the Indian Ocean generated a catastrophic tsunami that was recorded by a large number of tide gauges throughout the World Ocean. Part 1 of our study of this event examines tide gauge measurements from the Indian Ocean region, at sites located from a few hundred to several thousand kilometers from the source area. Statistical characteristics of the tsunami waves, including wave height, duration, and arrival time, are determined, along with spectral properties of the tsunami records.     

A Revised Evaluation of Tsunami Hazards along the Chinese Coast in View of the Tohoku-Oki Earthquake

Japan’s 2011 Tohoku-Oki earthquake and the accompanying tsunami have reminded us of the potential tsunami hazards from the Manila and Ryukyu trenches to the South China and East China Seas. Statistics of historical seismic records from nearly the last 4 decades have shown that major earthquakes do not necessarily agree with the local Gutenberg-Richter relationship. The probability of a mega-earthquake may be higher than we have previously estimated. Furthermore, we noted that the percentages of tsunami-associated earthquakes are much higher in major events, and the earthquakes with magnitudes equal to or greater than 8.8 have all triggered tsunamis in the past approximately 100 years. We will emphasize the importance of a thorough study of possible tsunami scenarios for hazard mitigation. We focus on several hypothetical earthquake-induced tsunamis caused by M _w 8.8 events along the Manila and Ryukyu trenches. We carried out numerical simulations based on shallow-water equations (SWE) to predict the tsunami dynamics in the South China and East China Seas. By analyzing the computed results we found that the height of the potential surge in China’s coastal area caused by earthquake-induced tsunamis may reach a couple of meters high. Our preliminary results show that tsunamis generated in the Manila and Ryukyu trenches could pose a significant threat to Chinese coastal cities such as Shanghai, Hong Kong and Macao. However, we did not find the highest tsunami wave at Taiwan, partially because it lies right on the extension of an assumed fault line. Furthermore, we put forward a multi-scale model with higher resolution, which enabled us to investigate the edge waves diffracted around Taiwan Island with a closer view.     

EFFECT OF TSUNAMIS GENERATED IN THE MANILA TRENCH ON CHINA MAINLAND

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 M_W7.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.