Rigid block characteristics on subaerial landslide-tsunamis using a 3D coupled Eulerian-Lagrangian model

To quantitatively reveal how rock blocks falling into water affect the impulse waves,the influence of a rigid block on induced first wave and second wave is systematically investigated.The block characteristics include the initial velocity,density,volume,and incident angle,and the investigated wave behavior characteristics include the maximum kinetic energy of the water,the transformation ratio of the kinetic energy from the block to the waves,the duration of the waves,the maximum movement speed,and the maximum height and width of the waves.The coupled Eulerian-Lagrangian method(CEL)is introduced to establish the numerical models of the fluid-solid coupling,and a laboratory test of a rigid wedge sliding into water demonstrates that it can reasonably describe the dynamic behavior of a landslide-induced wave.A typical process of a block entering water and its energy variation are described and analyzed in detail.Further,the relationship between each characteristic parameter of the block and the waves is quantitatively investigated and fitted.The simulation results show that energy exchange between the block and the water is very rapid after the block collides with the water.The maximum kinetic energy,maximum velocity,duration,and side dimension of the waves mainly increase non-linearly with the above characteristic parameters of the block.The transformation ratio of the kinetic energy from the block to the water,the first wave,and the second wave are usually below 60%,45%,and 30%,respectively.The velocity of the block first decreases and then maintains a constant speed after entering the water.The displacement of the block increases linearly with the initial velocity and density of the block and exponentially increases with the block volume at different times.With the increase in the incident angle of the block,the kinetic energy and scale of the second wave increase correspondingly.

Simulation on 2D underwater landslide-induced tsunamis

The Green-Naghdi （GN） theory is used here to simulate two-dimensional （2D） underwater landslide-induce tsunamis. Finite difference method is used to solve the GN equations. GN theory has different levels. There are GN-1, GN-2, .--, GN-K theory in GN theory. When K goes up, the GN equations will be more complicated, and the results will be more accurate. For the case simulated here, results of GN-5 theory are better than results of GN-3 theory. GN-7 results are almost the same as GN-5 results. That means GN-5 results are the converged results from GN theory. GN-5 results agree well with other＇s experimental results.

Scenarios of Local Tsunamis in the China Seas by Boussinesq Model

The Okinawa Trench in the East China Sea and the Manila Trench in the South China Sea are considered to be the regions with high risk of potential tsunamis induced by submarine earthquakes. Tsunami waves will impact the southeast coast of China if tsunamis occur in these areas. In this paper, the horizontal two-dimensional Boussinesq model is used to simulate tsunami generation, propagation, and runnp in a domain with complex geometrical boundaries. The temporary varying bottom boundary condition is adopted to describe the initial tsunami waves motivated by the submarine faults. The Indian Ocean tsunami is simulated by the numerical model as a validation case. The time series of water elevation and runup on the beach are compared with the measured data from field survey. The agreements indicate that the Boussinesq model can be used to simulate tsunamis and predict the waveform and runup. Then, the hypothetical tsunamis in the Okinawa Trench and the Manila Trench are simulated by the numerical model. The arrival time and maximum wave height near coastal cities are predicted by the model. It turns out that the leading depression N-wave occurs when the tsunami propagates in the continental shelf from the Okinawa Trench. The scenarios of the tsunami in the Manila Trench demonstrate significant effects on the coastal area around the South China Sea.

Comparison of two great Chile tsunamis in 1960 and 2010 using numerical simulation

Great differences in hazard and losses were shown from two tsunamis, both generated in Chile, one in 1960 and the other in 2010. Numerical simulation was applied to the tsunami analysis. The fault dislocation of the seafioor was assumed to equal to the initial tsunami wave field, which can be calculated by the formula of fault dislocation in the elastic isotropic half-space. The linear long wave theory was used as the tsunami hydrodynamic model, and the finite difference method and leap-frog scheme were selected for solving the equations. The accuracy of the simulated results was verified by the observed data in five tide gauges. By means of two scenario tsunamis, the analytical results show that the earthquake magnitude, bathymetry in rupture zone and rapid release of warning information in 2010 tsunami are the main explanations of the aforementioned great difference.

Will oscillating wave surge converters survive tsunamis?

With an increasing emphasis on renewable energy resources, wave power technology is becoming one of the realistic solutions. However, the 2011 tsunami in Japan was a harsh reminder of the ferocity of the ocean. It is known that tsunamis are nearly undetectable in the open ocean but as the wave approaches the shore its energy is compressed, creating large destructive waves. The question posed here is whether an oscillating wave surge converter （OWSC） could withstand the force of an incoming tsunami. Several tools are used to provide an answer： an analytical 3D model developed within the framework of linear theory, a numerical model based on the non-linear shallow water equations and empirical formulas. Numerical results show that run-up and draw-down can be amplified under some circumstances, leading to an OWSC lying on dry ground t

Extraction of two tsunamis signals generated by earthquakes around the Pacific rim

As one of the ocean sudden natural disasters,the tsunami is not easily to differentiate from the ocean variation in the open ocean due to the tsunami wave amplitude is less than one meter with hundreds of kilometers wavelength. But the wave height will increases up to tens of meters with enormous energy when the tsunami arrives at the coast. It would not only devastate entire cities near coast,but also kill millions of people. It is necessary to forecast and make warning before the tsunami arriving for many countries and regions around the Pacific rim. Two kinds of data were used in this study to extract the signals of 2011 Tohoku tsunami and 2014Iquique tsunami. Wave undulations from DART（ Deep-ocean Assessment and Reporting of Tsunamis） buoys and SLA from altimetry could extract the tsunami signals generated by this two earthquake. The signals of Tohoku tsunami were stronger than that of Iquique tsunami probably due to the 2011 Tohoku tsunami was generated by a magnitude 9. 0 earthquake and the 2014 Iquique tsunami was triggered by a magnitude 8. 2 earthquake.

Numerical Study on the Effect of the Manila Seismic Tsunami on the Guangdong Coast——The Nonlinear Effects of Tides and Tsunamis

The coupling numerical model between astronomical tide and tsunamis was built based on the data of the COMCOT model, which was verified by the Japan "3·11" tsunami. The tsunami source of the Manila seismic tsunami was designed and computed to analyze the risk of tsunami which will happen in the Guangdong coast. According to the results, the maximum quantity of water increasing and the time of the tsunami arriving at the Guangdong seacoast were calculated. The coupling simulation between astronomical tide and tsunamis has significant meaning for evaluating the risk and early-warning of tsunamis in the Guangdong coastal area.

Numerical simulation of landslide-generated tsunamis in lakes:A case study of the Xiluodu Reservoir

China has planned and built several world-class cascade high dams in Jinsha River,Dadu River,Lancang River,and Yarlung Zangbo River.The complex geological conditions in the reservoir area and numerous large-scale landslide bodies make the potential disaster risk of overflowing and cascade dam failure caused by landslide-generated tsunami under increasing severe situations.However,the study on describing and predicting the complex dynamic processes of generation,propagation,overflowing,wave setup,and the interaction between tsunami and lakeshore has not been systematically carried out.Based on the high-order Boussinesq-type equations,the development of the dynamic system of tsunamis in lakes coupled with the landslide process is realized using the finite volume method in this paper.To verify the accuracy and reliability of the study,the Xiluodu Reservoir is selected as the object to simulate the potential landslide-generated tsunamis.The factors such as the generation and propagation of tsunamis,dam overflowing,and wave setup in the downstream river are quantitatively evaluated and analyzed.The constructed landslide with a total volume of 24×10^(6)m^(3)generates a near-field wave amplitude of about 28 m.The maximum wave run-up height is about 95 m,the volume of the dam overflowing water up to 2.13×10^(6)m^(3),and the maximum wave height above the dam crest presents an M-shaped distribution.This LGWs event raises the downstream water level by nearly 40 m.The results show that the risk of landslide-generated tsunamis in the reservoir area in China cannot be ignored.The developed Boussinesq-type equations coupled with the landslide dynamics can simulate the whole process of generation,propagation,runup,and estimating the overflowing water volume of the tsunamis in the lake,laying a foundation for the quantitative risk assessment of tsunamis in lakes of high cascade dams in China.

Tsunamis and tsunami warning:Recent progress and future prospects

Tsunamis are one of the most destructive disasters in the ocean.Large tsunamis are mostly generated by earthquakes,and they can propagate across the ocean without significantly losing energy.During the shoaling process in coastal areas,the wave amplitude increases dramatically,causing severe life loss and property damage.There have been frequent tsunamis since the 21 st century,drawing the attention of many countries on the study of tsunami mechanism and warning.Tsunami records also play an essential role in deriving earthquake rupture models in subduction zones.This paper reviews the recent progress and limitations of tsunami research,from the aspects of tsunami generation,propagation,inversion and warning.Potential tsunami warning strategies are discussed and future prospects on tsunami research are provided.

Numerical simulation of submarine landslide tsunamis using particle based methods

This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics（SPH）, as a particle based method, is for solving problems of fast moving boundaries in the field of continuum mechanics. Other particle based methods, like the discrete element method（DEM）, are suitable for modeling the displacement and the collision related to the rigid landslides. In the present work, we use the SPH and the DEM to simulate tsunamis generated by rigid and deformable landslides with consideration of submerged conditions. The viscous free-surface flows are solved by a weakly compressible SPH and the displacement and the rotation of the rigid body slides are calculated using a multi-sphere DEM allowing for modeling solids of arbitrarily complex shapes. The fluid-solid interactions are simulated by coupling the SPH and the DEM. A rheology model combining the Papanastasiou and the Herschel-Bulkley models is applied to represent the viscoplastic behavior of the non-Newtonian flow in the submarine deformable landslide cases. Submarine landslide tsunamis due to rigid and deformable landslides are both simulated as typical landslide cases in this investigation. Our simulated results and the previous experimental results in the literatures are in good agreement, which shows that the proposed particle based methods are capable of modeling the submarine landslide tsunamis.

Risk Assessment of Tsunamis Along the Chinese Coast Due to Earthquakes

China’s coastal areas are densely populated,economically developed, and located in close proximity to several potential tsunami sources;therefore, tsunami risk cannot be ignored. This study assessed tsunami risk in coastal areas of China by developing a framework for tsunami risk assessment from the perspectives of hazards,vulnerability, and exposure. First, a probabilistic tsunami hazard assessment(PTHA) model was applied to estimate the potential tsunami sources in both local crustal faults and circum-Pacific subduction zones based on numerical simulations. The output of the PTHA includes tsunami wave height distributions along the coast. Then, an indicator system reflecting exposure and vulnerability to tsunamis in the coastal areas of China was established by using the entropy method and analytic hierarchy process.The PTHA findings show that the tsunami wave height is close to 3 m on the southern coast of the Bohai Sea, the Pearl River Estuary, and the Yangtze River Delta and exceeds 2 m near the Taiwan Strait for the 2000-year return period. The results of the tsunami risk assessment show that the cities at the highest risk level(level I) include Tangshan, Yantai, and Hong Kong, while cities at the high risk level(level II) include Fuzhou, Xiamen, and Quanzhou near the Taiwan Strait and many cities on the Yangtze River Delta, the Pearl River Estuary, and the southern coast of the Bohai Sea. Our findings can provide an understanding of differences in tsunami risk between Chinese coastal cities that may be affected by tsunamis in the future.

散粒体滑坡涌浪运动特征与能量转化规律研究

滑坡涌浪灾害造成的伤亡和损失一般远超滑坡灾害本身,已引起广泛关注。针对该类灾害的预测评价一直是研究的难点,精细刻画滑坡涌浪运动特征和能量转化规律是解决问题的关键和前提。因此,以天然碎石模拟散粒体滑坡,建立缩尺的矩形水槽滑坡涌浪三维物理试验模型,研究滑坡体积、速度和水深对堆积体形态和涌浪特征的影响以及滑体与水体能量转化规律。研究结论包括:(1)滑体通过冲击、挤压和抬升水体形成涌浪的特征受滑体规模和形态影响显著;厚度大且速度快的滑体趋向产生非线性过渡波,厚度小且速度低的滑体一般产生非线性震荡波;(2)基于试验提出的体积效应和速度效应揭示了一定水深的涌浪规模与滑坡体积和速度呈正相关规律,水深效应则解释了一定波能的涌浪在波形和速度上的分布差异,结果表明浅水比深水工况涌浪高度平均高出67%、波速平均高出51.17%;(3)基于试验的数值反演总结了滑体势能与波能转化率为1.00%~3.07%,由于三维试验中滑体与水体在水平和环向的扩散造成更多能量耗散,相较二维试验转化率较低。通过探讨散粒体滑坡运动过程和首浪产生、传播及爬高特征,揭示滑体-水体能量转化基本规律,研究成果对滑坡涌浪防灾减灾工作具有一定的理论价值和意义。

全球海啸相关研究进展综述

With the advancement of the global economy,the coastal region has become heavily developed and densely populated and suffers significant damage potential considering various natural disasters,including tsunamis,as indicated by several catastrophic tsunami disasters in the 21st century.This study reviews the up-to-date tsunami research from two different viewpoints:tsunamis caused by different generation mechanisms and tsunami research applying different research approaches.For the first issue,earthquake-induced,landslide-induced,volcano eruption-induced,and meteorological tsunamis are individually reviewed,and the characteristics of each tsunami research are specified.Regarding the second issue,tsunami research using post-tsunami field surveys,numerical simulations,and laboratory experiments are discussed individually.Research outcomes from each approach are then summarized.With the extending and deepening of the understanding of tsunamis and their inherent physical insights,highly effective and precise tsunami early warning systems and countermeasures are expected for the relevant disaster protection and mitigation efforts in the coastal region.

Experimental study on the mitigation effect of mangroves during tsunami wave propagation

Mangroves are crucial for protecting coastal areas against extreme disasters such as tsunamis and storm surges.An experimental study was conducted to determine how mangroves can mitigate the tsunami wave propagation.The test was performed in a flume, where mangrove models were installed on a slope, and dam-burst waves were used to simulate tsunami waves. To study how mangrove forests reduce the impact of tsunamis, this paper measured the heights of the incoming waves under different initial conditions(tsunami wave intensity and initial water depth) and plant factors(arrangement and distribution density) and described the reduction process. The results show that, after passing through the mangrove, the tsunami bore height will decrease within a certain range as the initial water depth increases. However, there is no correlation between the increase of inundation level and the drop of water level. The bore height attenuation is more significant at higher density of mangroves,but after tsunami passing through the mangroves, the relative bore height will decrease. When the distribution density of mangroves is constant, the wave attenuation at different locations(before, on and after the slope)shows different relationships with the initial water depth and wave height for different models. The transmission coefficient(K_(i)) shows a parabolic correlation with its density. The proportion of the energy loss caused by the mangrove resistance to the total energy(E_(b)) is defined as C_(m2). The variation trend of C_(m2) corresponds to the tsunami wave energy attenuation rate(C_(a)) and K_(i).

Radio Direction Finding Method to Mitigate Tsunami Risk in Sierra Leone

In this study,the Radio Direction Finding method is proposed for the detection of electromagnetic signals,in the VLF band,to try to anticipate the occurrence of potentially destructive geophysical events.The experimentation concerns the interception of electromagnetic anomalies in Sierra Leone,in the five-day time window,associated with seismic events that could potentially generate tsunamis.The area of investigation is Sierra Leone,whose coastline is subjected to tidal wave hazards triggered by earthquakes generated in the Mid-Atlantic Ridge.Although Sierra Leone is not affected by recurrent earthquakes,there is nevertheless a low probability,estimated at 2 percent,of the occurrence of destructive earthquakes in the next 50 years.Also in estimates,the risk of rogue and potentially damaging waves is estimated to strike the Sierra Leone coast at least once in the next 10 years.The Radio Direction Finding experiment carried out continuously 24/7,has shown a close relationship between increased radio-anomalies,in the frequencies of 6,000 Hz,a time window between electromagnetic anomaly detection and the imminence of an earthquake,and higher frequency times for the risk of earthquake occurrence in the Mid-Atlantic Ridge.

Tsunami Hazard Assessment on Qatar Coastline from Makran Earthquakes Considering Tidal Effect and Coastal Landslides Scenarios

To assist the analysis of tsunami hazards for Qatar coastal areas were using numerical model. By Tsunamis waves created from submarine earthquakes of magnitude of (Mw) 8.6 and 9.0 in Richard scale along the Makran Subduction Zone (MSZ) as well as coastal landslides with soil volume of 1.25 to 2.0 km3 along Iranian coast inside the Arabian Gulf is considered. TUNAMI-N2KISR model (Al-Salem) was applied in this study to predict the tsunami propagation and magnitude of Tsunami induced wave heights. The model adopts to solve shallow water equations describing nonlinear long-wave theory. The model also incorporate tidal effect inside the Arabian Gulf as a tsunami travel time from Makran Subduction to Qatar coastline takes more than 9 hours with the tidal range of about 1.6 m during Spring Tide event. For coastal landslides, tsunami generation was simulated using a two-layer numerical model, developed by solving nonlinear long-wave equations. Two-layer model was used to determine initial wave deformation generated by a landslide case. Then TUNAMI-N2KISR was use to simulate tsunami wave propagation. Tsunami waves from landslide scenario arrived after 2.5 - 3 hr with maximum tsunami amplitudes along coasts of Ras laffan-Qatar were 0.8 to 1.0 m. Incorporation of ocean tide is found to impose some small effect on tsunami amplitude at Qatar coastline and nearby areas for the Mw 9.0 earthquake due to small tidal range in this area. In addition, it is found that the tsunami arrival time has become shorter.

Converting Tsunami Wave Heights to Earthquake Magnitudes

There is a fairly strict relation between maximum tsunami wave heights and causation earthquake magnitudes. This provides a new tool for estimating the magnitude of past earthquakes from the observed wave heights of related paleo-tsunami events. The method is subjected to a test versus two paleoseismic events with multiple independent estimates of corresponding earthquake magnitude. The agreement to the tsunami wave height conversion is good, confirming very high magnitudes of M 8.5 - 9.0 and M 8.4 - 8.5. Applying the same method to two Late Holocene events of methane venting tectonics indicates a ground shaking of forces equivalent to a M 8.0 earthquake, seriously changing previous long-term crustal hazard assessments.

Trapping Mechanism of Submerged Ridge on Trans-oceanic Tsunami Propagation

Based on the linear shallow water equations,an analytic solution of trapped waves over a symmetric parabolicprofile submerged ridge is derived.The trapped waves act as propagating waves along the ridge and as standing waves across the ridge.The amplitude gets the maximum at the ridge top and decays gradually towards both sides.The decaying rate gets more gently with higher modes.Besides,an explicit first-order approximate dispersion relation is derived to simplify transcendental functions in the exact solution,which is useful to describe trapped waves over shallowly submerged ridges in reality.Furthermore,the trapping mechanism of the submerged ridge waveguides on the trans-oceanic tsunami propagation can be explained by the ray theory.A critical incident angle exists as a criterion to determine whether the wave is trapped.Besides,a trapped parameter γ is proposed to estimate the ratio of the energy trapped by the oceanic ridge if a tsunami is generated at its top.

Laboratory Generation of Solitary Waves:An Inversion Technique to Improve Available Methods

Solitary waves are often used in laboratory experiments to study tsunamis propagation and interaction with coasts. However, the experimental shape of the waves may differ from the theoretical one. In this paper, a correction technique aiming at minimizing the discrepancies between the two profiles is presented. Laboratory experiments reveal their effectiveness in correcting the experimental shape of solitary waves, mainly for low nonlinearities.

Feasibility Analysis of Performance Validation for Satellite Altimeters Using Tide Gauge and Deep-ocean Bottom Pressure Recorder

Independent of traditional approach of satellite altimeter calibration, the feasibility of altimeter validation using tide gauge located on solitary island at open sea (TGSI) and deep-ocean bottom pressure recorder (DBPR) separately is initially studied. Bias of Jason-3 sea surface height (SSH) and relative SSH bias (Δbias) between Jason-2 and Jason-3 is calculated using the data of tide gauge on Harvest oil platform, tide gauge No. 1890000 and DBPR No. 21419. The standard deviations of calculated SSH bias sequence are 3.98 cm, 2.87 cm and 8.61 cm respectively, and Δbias (Jason-3—Jason-2) is -3.62± 2.17 cm , -2.58±1.97 cm and -2.60±1.30 cm respectively. Comparing to the results reported by international calibration sites, the results show that Jason-3 SSH is 3.0 cm lower than that of Jason-2, the selected DBPR is appropriate to the calculation of relative SSH bias between Jason-2 and Jason-3, but it is not suitable for calibration or validation of single satellite, TGSI is appropriate to both.