A nonlinear parametric model based on a power law relationship for predicting the coastal tsunami height

When a subduction-zone earthquake occurs, the tsunami height must be predicted to cope with the damage generated by the tsunami. Therefore, tsunami height prediction methods have been studied using simulation data acquired by large-scale calculations. In this research, we consider the existence of a nonlinear power law relationship between the water pressure gauge data observed by the Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) and the coastal tsunami height. Using this relationship, we propose a nonlinear parametric model and conduct a prediction experiment to compare the accuracy of the proposed method with those of previous methods and implement particular improvements to the extrapolation accuracy.

     

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.     

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.     

Tsunami flooding extension for coastal zones of Mexico

Abstract

Joint probabilities of runup height due to astronomical tides and tsunamis are calculated for five bays on the Pacific Ocean coast of Mexico. Runup heights equaled or exceeded on the average of once per 100 years, and once per 200 years, for those bays, are evaluated. Calculations are based on wave height probabilities obtained from sets of harmonic constituents and tsunami heights recorded during the last 28 years at tidal stations on each location.

Extension of potential tsunami terrestrial flooding for populated coastal zones around those bays can be determined from the previous results. One case is illustrated showing endangered features in the hazardous zones, and recommendations to prevent casualties and property damages are given.     

How the tide influences dangerous level rises on the coast of the Sea of Okhotsk and adjacent areas in cases of tsunami and storm surges

A long-term sea level series were analyzed, recorded at 12 coastal tide gauges located on the coast of the Sea of Okhotsk and the Pacific coast of the Kuril Islands and Kamchatka. Estimates for the maximum heights of the tidal level, storm surges, and tsunamis were obtained separately, as well as for the rare recurrence of the total sea level height with the probability of these individual components superposition. The maximum total height of the sea level without a tsunami were obtained for the Magadan station, where the main factor is anomalously large tides, as well as for Iturup and Matua islands, where the highest storm surges were recorded. The minimum values were obtained for Ust’-Kamchatsk and Malokuril’sk (Shikotan Island) on different flanks of the study area. When a tsunami is included, the maximum values of possible total sea level rises were observed on the Pacific coast of the Kuril Ridge and the influence of tides and meteorologically induced oscillations are small. On the east coast of Kamchatka adjacent to the considered closed area, the role of tsunamis is much smaller. At the Kuril’sk station, where the height of the largest tsunami (Chilean, May 1960) was about a half the strongest surge height, the contribution of the tsunami scarcely affected the resulting estimates. As a rule, the contribution of a tsunami becomes significant at other stations on the coast of the Sea of Okhotsk with a recurrence period of 100 years.     

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海啸模式对本次海啸事件近场传播特征进行了模拟,模拟结果与观测吻合较好。最后,基于实测和模拟结果,详细分析了此次地震海啸的近场分布特征,发现除受海啸源的强度和几何分布特征影响外,近岸海啸波还主要受地形特征控制,在与特定地形相互作用后波幅产生放大效应,会进一步加剧海啸造成的灾害。     

The detectability of shallow slow earthquakes by the Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET) in Tonankai district,Japan

In order to understand the characteristics of shallow very low-frequency (VLF) events as revealed by recent ocean-floor observation studies, we perform a trial simulation of earthquake cycles in the Tonankai district by taking the characteristics of the 1944 Tonankai earthquake and assuming that slow earthquakes occur on numerous small asperities. Our simulation results show that the increase of moment release rate of shallower VLF events in the pre-seismic stage of a megathrust earthquake is higher than that of deeper VLF events. This increase may make leveling change due to VLF swarms detectable at Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET). We also introduce the time series of hydraulic pressure data at DONET, comparing with the leveling change expected from our numerical simulation. Since leveling change due to shallower VLF swarms is so local as to be incoherent, removal of the moving-averaged data from the data stacked by four nearby observation points in the same node may be useful to detect the short-term local leveling change.     

南海潜在海啸灾害的模拟

结合南海海域的地形条件、地质构造、地震学特征以及历史地震记录,在回顾总结国内外学者研究的基础上,分析了南海可能引发地震海啸的震源区域,并讨论了在我国南海沿岸发生海啸灾害的潜在可能性。采用目前在国际上广泛使用的COMCOT海啸模式,对马尼拉海沟的潜在地震引发的海啸进行了数值模拟计算,计算中包含了由地震参数到海面初始变形的转换、海啸的深水传播过程以及海啸的浅水传播过程。采用三重嵌套网格,外层网格对应于大范围的深水区域,使用球坐标系下的线性控制方程;第二层网格对应中等范围的较浅水区域,使用球坐标系下的非线性控制方程;第三层网格对应小范围的浅水区域,使用直角坐标系下的非线性控制方程。由模拟计算得到的海啸传时分布、近岸海面升降强度、四个特定点上海面高度随时间变化等的结果表明,我国南海沿岸遭受海啸袭击的可能性是存在的,应进一步对南海海啸进行监测、预警和研究。COMCOT模式性能良好,可用于对南海潜在地震海啸的进一步模拟研究。     

Run-Up of a Solitary Tsunami Wave on a Sloping Coast

Within the framework of the nonlinear theory of long waves, we perform the numerical analysis of the one-dimensional run-up of solitary tsunami waves upon a plane sloping coast. We study the dependences of the run-up heights on the parameters of waves at the entrance of the shelf zone and on the slope of the coast. The run-up heights of tsunami waves are estimated for the bottom topography typical of the south coast of the Crimean Peninsula. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 4, pp. 11–18, July–August, 2005.     

Statistical Simulation of Boxing Day Tsunami of The Indian Ocean and a Predictive Equation for Beach Run up Heights Based on Work-Energy Theorem

The tsunami similar to the one that has occurred in December 26, 2004 (Boxing Day Tsunami) in the Indian Ocean is simulated using the expression derived from Modified Weibull Distribution (for maximum wave height simulation) for extreme wave height predictions. The tuning coefficient plays a significant role in estimating the tsunami heights at various stages. It follows well defined mathematical laws at different stages. It is time dependent in the first three stages and depth dependent in the last two stages. The beach run-up heights estimated by the expression derived from the work-energy relation are comparable with observed values with reasonable accuracy.     

Statistical Simulation of Boxing Day Tsunami of The Indian Ocean and a Predictive Equation for Beach Run up Heights Based on Work-Energy Theorem

The tsunami similar to the one that has occurred in December 26, 2004 (Boxing Day Tsunami) in the Indian Ocean is simulated using the expression derived from Modified Weibull Distribution (for maximum wave height simulation) for extreme wave height predictions. The tuning coefficient plays a significant role in estimating the tsunami heights at various stages. It follows well defined mathematical laws at different stages. It is time dependent in the first three stages and depth dependent in the last two stages. The beach run-up heights estimated by the expression derived from the work-energy relation are comparable with observed values with reasonable accuracy.     

Temporal variability of tsunami arrival detection distance revealed by virtual tsunami observation experiments using numerical simulation and 1-month HF radar observation

The combination of a high-frequency ocean surface radar and a tsunami detection method should be assessed as the onshore-offshore distribution of tsunami detection probability, because the probability will vary in accordance with the signal-to-noise ratio (SNR) and the tsunami magnitude in addition to the radar system specifications. Here, we statistically examine the tsunami detection distance based on virtual tsunami observation experiments by using signals received by a high-frequency radar in February 2014 installed on the southern coast of Japan and numerically simulated velocities induced by a Nankai Trough earthquake. In the experiments, the Doppler frequencies associated with the simulated velocities were superimposed on the receiving signals of the radar, and the radial velocities were calculated from the synthesized signals by the fast Fourier transform. Tsunami arrival was then detected based on the temporal change in the cross-correlation of the velocities, before and after tsunami arrival, between two points 3 km apart along a radar beam. We found that the possibility of tsunami detection primarily depends on the kinetic energy ratio between tsunami current and background current velocities. The monthly average detection probability is over 90% when the energy ratio exceeds 5 (offshore distance: 9 km ≤ L ≤ 36 km) and reduces to 50% when the energy ratio is approximately 1 (L = 42 km) over the shelf slope. The ratio varied with the background current physics and SNR, which was mainly affected by ocean surface wave heights and ionospheric electron density.     

2015年9月16日智利Mw 8.3级地震:海啸波从近场到远场的特征分析

On September 16, 2015, an earthquake with magnitude of M_w 8.3 occurred 46 km offshore from Illapel, Chile,generating a 4.4-m local tsunami measured at Coquimbo. In this study, the characteristics of tsunami are presented by a combination of analysis of observations and numerical simulation based on sources of USGS and NOAA. The records of 16 DART buoys in deep water, ten tidal gauges along coasts of near-field, and ten coastal gauges in the far-field are studied by applying Fourier analyses. The numerical simulation based on nonlinear shallow water equations and nested grids is carried out to provide overall tsunami propagation scenarios, and the results match well with the observations in deep water and but not well in coasts closed to the epicenter. Due to the short distance to the epicenter and the shelf resonance of southern Peru and Chile, the maximum amplitude ranged from 0.1 m to 2 m, except for Coquimbo. In deep water, the maximum amplitude of buoys decayed from9.8 cm to 0.8 cm, suggesting a centimeter-scale Pacific-wide tsunami, while the governing period was 13–17 min and 32 min. Whereas in the far-field coastal region, the tsunami wave amplified to be around 0.2 m to 0.8 m,mostly as a result of run-up effect and resonance from coast reflection. Although the tsunami was relatively moderate in deep water, it still produced non-negligible tsunami hazards in local region and the coasts of farfield.     

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.     

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.     

Time-Dependent Tsunami Source Following the 2018 Anak Krakatau Volcano Eruption Inferred from Nearby Tsunami Recordings

The eruption of the Anak Krakatau volcano,Indonesia,on 22 December 2018 induced a destructive tsunami(the Sunda Strait tsunami),which was recorded by four nearby tidal gauges.In this study we invert the tsunami records and recover the tsunami generation process.Two tsunami sources are obtained,a static one of instant initial water elevation and a time-dependent one accounting for the continuous evolution of water height.The time-dependent results are found to reproduce the tsunami recordings more satisfactorily.The complete tsunami generation process lasts approximately 9 min and features a two-stage evolution with similar intensity.Each stage lasts about 3.5 min and elevates a water volume of about 0.13 km3.The time,duration and volume of the volcano eruption in general agree with seismic records and geomorphological interpretations.We also test different sizes of the potential source region,which lead to different maximum wave height in the source area,but all the results of time-dependent tsunami sources show the robust feature of two stages of wave generation.Our results imply a time-dependent and complex process of tsunami generation during the volcano eruption.     

Registration of tsunamis in the open ocean

Abstract

The small tsunami of February 23, 1980, which originated near the southern part of the Kuril Islands was recorded by the bottom vibrotron sensor installed on the shelf near Shikotan Island and connected to the island observatory by underwater cable. The marigram of the tsunami is given and compared with marigrams obtained at shore tide gauges. Some spectral features of the records are discussed.     

Seismic monitoring of western Pacific typhoons

Typhoons inflict large damage to societies, but are usually difficult to monitor in close proximity in real-time without expensive instruments. Here we study the possibility of using seismic waveforms on the seafloor and on land to monitor the turning of a far away or approaching typhoon. Up to 67% of the typhoons making landfall in Taiwan come from the eastern shore, so that we deployed broadband ocean-bottom seismometers (OBSs) offshore eastern Taiwan in 2006 to study ground motion in close proximity to a typhoon. Typhoons generate ocean waves, which generate pressure signals in the water column before being transmitted to the seafloor as seismic waves and recorded by the OBSs. The ground motions on the seafloor correlate with locally increased (ocean) wave heights and wave periods, suggesting that the ground motions are mostly induced by in situ or nearby pressure fields, as shown by coherence function analyses. When a typhoon turns and changes wave-wave interaction near the source region, a new set of en echelon patterns develops which can be observed by OBSs and land stations. Similar features occur when a typhoon crosses a landmass and re-enters the ocean. The energy level ratio between the single-frequency and double-frequency microseisms also changes abruptly when the typhoon turns. These features can potentially help near real-time early warning with little cost to complement other conventional typhoon early warning methods.     

The buckling of a corrugated carbon fibre cylinder under external hydrostatic pressure

A theoretical and an experimental investigation was carried out, where a carbon fibre corrugated circular cylinder was tested to destuction under external hydrostatic pressure. The theoretical investigation was via the finite element method, where the structure was modelled with several orthotropic axisymmetric thin-walled shell elements. The experimental observations were aided with strategically placed strain gauges. Comparison between theory and experiment showed that the experimentally observed buckling pressure was a little lower than the theoretical prediction. This may have been due to the fact that the model had slight initial geometrical imperfections in the circumferenential direction.