2017年9月8日墨西哥沿岸Mw8.2级地震海啸观测数据分析与模拟

2017年9月8日4时49分（UTC）,墨西哥瓦哈卡州沿岸海域（15.21°N,93.64°W）发生M_w8.2级地震,震源深度30 km。强震在该海域引发海啸,海啸对震源附近数百千米范围内造成了严重影响。位于太平洋上的多个海啸监测网络捕捉到了海啸信号并详细记录了此次海啸的传播过程。本文选用了近场2个DART浮标和6个验潮站的水位数据,通过潮汐调和分析和滤波分离出海啸信号,对近场海啸特征值进行了统计分析,并采用小波变换分析方法进一步分析了海啸的波频特征。基于Okada弹性位错理论断层模型计算得到了强震引发的海底形变分布,并采用MOST海啸模式对本次海啸事件近场传播特征进行了模拟,模拟结果与观测吻合较好。最后,基于实测和模拟结果,详细分析了此次地震海啸的近场分布特征,发现除受海啸源的强度和几何分布特征影响外,近岸海啸波还主要受地形特征控制,在与特定地形相互作用后波幅产生放大效应,会进一步加剧海啸造成的灾害。     

印尼附近海域地震海啸发生的构造背景综述

2004年12月26日在印度尼西亚苏门答腊岛以西海域发生的里氏9级地震,引发了历史上第五大地震海啸,引起了地学界的高度关注。印度尼西亚西部位于环太平洋地震带和地中海-喜马拉雅地震带结合部位,处于洋壳和陆壳的汇聚边界以及弧状压缩构造格局中。苏门答腊陆缘NW—SE向延伸约1600km,苏门答腊俯冲带源于印度-澳大利亚和欧亚板块3°N的汇聚。苏门答腊弧前区沿着平行于俯冲海沟的两大走滑断裂——苏门答腊大断裂和明打威断裂向北运动。本文对印度尼西亚苏门答腊岛附近海域的地质背景、大地构造特征和地震活动评估等进行了综述,对该海域产生大地震的深部动力学机制进行了分析,并从地质、地球物理资料的获取与分析角度出发,对区域灾害模型进行了探讨。     

局地海啸波传播数学模型的研究

推导了地震海啸波传播的理论方程,并时局地地震海啸情况下的理论方程进行了求解.基于Boussinesq方程出发,建立了二维局地海啸渡传播数学模型;对局地海啸进行了数值计算,计算方法采用有限差分方法.差分格式采用交替方向隐格式,即ADI方法.利用该模型对发生在台湾海峡的一次地震海啸进行了模拟计算.将计算结果与理论方程的计算结果以及实际情况进行了对比.计算结果较好.为局地海啸波传播提供了一种模拟方法.     

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

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

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

The role of horizontal impulses of the faulting continental slope in generating the 26 December 2004 tsunami

For a long time, people have believed that the vertical displacement of seafloor due to undersea earthquakes is the primary cause of tsunami genesis. However, seismically-inverted seafloor deformation of the 2004 Sumatra–Andaman earthquake shows that the total vertical displacement is not enough to have generated the powerful Indian Ocean tsunami. Based on the seismically-inverted data and a three-dimensional ocean general circulation model (OGCM), we show that the momentum force, transferred by the horizontal impulses of the faulting continental slope in that earthquake, has accounted for two thirds of the satellite-observed tsunami height and generated kinetic energy 5 times larger than the potential energy due to the vertical displacement. The asymmetric tsunami pattern, recorded by tide-gauges showing leading-elevation waves toward Sri Lanka and leading-depression waves toward Thailand, is best explained by the horizontally-forced mechanism. The same mechanism has also explained the March 2005 Nias earthquake and tsunami data, suggesting that the horizontal motions of faulting have played more important roles in tsunami genesis than previously thought.     

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

2011年日本Mw 9.0级地震的GPS响应研究

2011年3月11日日本宫城县以东太平洋海域发生Mw9.0级特大地震,造成了地表的严重错位并引发海啸。文中利用位于日本及周边国家的IGS站和国家海洋局GPS业务站观测数据,采用作者研制的精密单点定位(PPP)软件UniP,对此次地震的GPS数据响应进行了研究。结果表明:(1)GPS观测数据能清晰、连续地记录震时地表形变的过程,我国CHAN,NCST等站点水平方向的震时最大位移在10 cm以内,高程方向的震时最大位移在15 cm以内,且形变以可恢复性的弹性形变为主。(2)我国距震中较远,受此次日本地震的影响较小,且大部分站点是在东坐标方向出现不同程度的震后永久性位移。其中CHAN站点的震后位移最为明显,东向形变量为(1.8±0.11)cm;NCST、NLHT站点次之,东向形变量分别为(1.1±0.26)cm和(1.0±0.18)cm。(3)地震波传输到国家海洋局GPS业务站NCST、NLHT等的时间约为10 min,比海啸在深海的传播速度快约14倍,可为海啸预警提供所需的时间差。这些结果显示出GPS能够为地震监测和动力学特征研究提供有价值的基础资料,也表明中国沿海GPS业务观测系统在海底地震监测、海啸预警服务中的应用潜力。     

Energetics of the Tsunami of 26 December 2004 in the Indian Ocean: A Brief Review

The energetics of the most destructive tsunami in historical time, and that of the under ocean earthquake that triggered this tsunami of 26 December 2004 in the Indian Ocean have been briefly reviewed. This latest tsunami has several other unique characteristics besides being one of the worst natural disasters in human history. It is the first truly global tsunami after modern seismographic and sea level monitoring networks have been put in place. It was the first tsunami on record detected by a satellite, even though at present, global satellite coverage of the oceans for real time tsunami detection is not adequate. Finally, the energy associated with the tsunami and the earthquake that triggered it is so large that speculation has been made about the normal modes of oscillation of the earth, that were triggered by the earthquake as well as some suggestions, that some of the earth's rotational characteristics may have temporarily changed to a discernible degree. Here, we briefly review the energetics of the tsunami and the earthquake that triggered it.     

越洋海啸波传播数学模型的研究

应用内外解匹配的方法和驻相法推导了柱坐标系下地震引起的水面波动方程的解,即推导了地震海啸波生成与传播的理论方程,并对越洋地震海啸理论方程进行了求解。基于Boussinesq方程出发,建立了二维越洋海啸波传播数学模型,并对越洋海啸进行了数值计算,计算方法采用有限差分方法,差分格式采用交替方向隐格式(Alternating Direction Implicit即ADI方法)。利用越洋海啸计算模式对发生在大洋中的地震海啸进行了模拟,将数值模拟结果与地震海啸波理论方程的计算结果进行了比较,两种计算结果吻合较好。     

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.     

The possibility of tsunami in the Sea of Okhotsk caused by deep-focus earthquakes

The earthquake that occurred on May 24, 2013, in the basin of the Sea of Okhotsk with a magnitude of 8.3 was the strongest in this region. We have modeled a possible tsunami caused by such an earthquake. The simulations confirm that the wave heights were sufficiently small because the earthquake epicenter depth was 640 km. We analyze the oscillations of the DART buoys in the vicinity of the earthquake source and show that they were not associated with the tsunami waves. Analysis of the available pressure gauge records at different points of the Sea of Okhotsk show that only in one case (Iturup Island) can the observed oscillations of the sea level with a height of approximately 4 cm be classified as tsunami waves.     

The Kuril Earthquakes and tsunamis of November 15, 2006, and January 13, 2007: Observations,analysis, and numerical modeling

Major earthquakes occurred in the region of the Central Kuril Islands on November 15, 2006 (M _w = 8.3) and January 13, 2007 (M _w = 8.1). These earthquakes generated strong tsunamis recorded throughout the entire Pacific Ocean. The first was the strongest trans-Pacific tsunami of the past 42 years (since the Alaska tsunami in 1964). The high probability of a strong earthquake (M _w ≥ 8.5) and associated destructive tsunami occurring in this region was predicted earlier. The most probable earthquake source region was investigated and possible scenarios for the tsunami generation were modeled. Investigations of the events that occurred on November 15, 2006, and January 13, 2007, enabled us to estimate the validity of the forecast and compare the parameters of the forecasted and observed earthquakes and tsunamis. In this paper, we discuss the concept of “seismic gaps,” which formed the basis for the forecast of these events, and put forward further assumptions about the expected seismic activity in the region. We investigate the efficiency of the tsunami warning services and estimate the statistical parameters for the observed tsunami waves that struck the Far Eastern coast of Russia and Northern Japan. The propagation and transformation of the 2006 and 2007 tsunamis are studied using numerical hydrodynamic modeling. The spatial characteristics of the two events are compared.     

Correlation between the intensity of tsunami in the Black Sea and the magnitude of underwater earthquakes

We perform the analysis of tsunami waves in the shelf zone of the Crimean peninsula generated by underwater earthquakes whose epicentres are located near the lower boundary of the continental slope. For this purpose, we use a one-dimensional nonlinear dissipative numerical model of long waves. The investigated distributions of the depth of the basin correspond to four points of the south coast of the Crimean peninsula. We use the empirical dependences of parameters of the sources of tsunami waves on the magnitude of the earthquake obtained earlier for the Pacific Ocean. It is shown that the height, vertical climb, and duration of tsunami waves increase with the magnitude of the earthquake. For equal magnitudes of the earthquake, the highest tsunamis on the south coast of the Crimea are observed between Alushta and Yalta. We also deduced a generalized regression dependence of the height of tsunami waves near the coast on the magnitude of the earthquake. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Observations of Tsunami Impact on the Coast of Kerala, India

The tsunami generated by the December 2004 Sumatra-Andaman earthquake had a devastating effect on some parts of Kerala coast, which is a coast located in southwest India. Results of post-tsunami field surveys carried out to understand the changes in coastal morphology and sediment characteristics in the worst affected Kayamkulam region of Kerala coast are documented in this study. Analysis of offshore bathymetric data indicates the shifting of depth contours towards shore, indicating erosion of sediments and deepening of innershelf due to the tsunami. Depth measurement along the backwater (T-S canal) in the hinterland region indicates siltation due to the inundation of the canal.     

On the near-bottom pressure variations in the region of the 2003 Tokachi-Oki tsunami source

The Tokachi-Oki earthquake was the strongest seismic event in 2003. The tsunami caused by the earthquake reached a height of four meters at the northeastern coast of Hokkaido. The JAMSTEC successfully recorded the variations of the near-bottom pressure in the region of the tsunami source. An analysis of the data reveals low-frequency (～ 0.15 Hz) elastic vibrations of the water layer. Estimates of the amplitude, velocity, and duration of the bottom deformation at the tsunami source were obtained.     

Observations of Tsunami Impact on the Coast of Kerala,India

The tsunami generated by the December 2004 Sumatra-Andaman earthquake had a devastating effect on some parts of Kerala coast, which is a coast located in southwest India. Results of post-tsunami field surveys carried out to understand the changes in coastal morphology and sediment characteristics in the worst affected Kayamkulam region of Kerala coast are documented in this study. Analysis of offshore bathymetric data indicates the shifting of depth contours towards shore, indicating erosion of sediments and deepening of innershelf due to the tsunami. Depth measurement along the backwater (T-S canal) in the hinterland region indicates siltation due to the inundation of the canal.     

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.     

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

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

Propagation Mechanisms of Incident Tsunami Wave in Jiangsu Coastal Area,Caused by Eastern Japan Earthquake on March 11,2011

At 13:46 on March 11, 2011(Beijing time), an earthquake of Mw=9.0 occurred in Japan. By comparing the tsunami data from Guanhekou marine station with other tsunami wave observation gathered from southeast coastal area of China, it was evident that, only in Guanhekou, the position of the maximum wave height appeared in the middle part rather than in the front of the tsunami wave train. A numerical model of tsunami propagation based on 2-D nonlinear shallow water equations was built to study the impact range and main causes of the special tsunami waveform discovered in Jiangsu coastal area. The results showed that nearly three-quarters of the Jiangsu coastal area, mainly comprised the part north of the radial sand ridges, reached its maximum tsunami wave height in the middle part of the wave train. The main cause of the special waveform was the special underwater topography condition of the Yellow Sea and the East China Sea area, which influenced the tsunami propagation and waveform significantly. Although land boundary reflection brought an effect on the position of the maximum wave height to a certain extent, as the limits of the incident waveform and distances between the observation points and shore, it was not the dominant influence factor of the special waveform. Coriolis force's impact on the tsunami waves was so weak that it was not the main cause for the special phenomenon in Jiangsu coastal area. The study reminds us that the most destructive wave might not appear in the first one in tsunami wave train.