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.

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.

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

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.

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.

A Criterion for Tsunami Hazard Assessment at the Local Scale

A criterion for tsunami hazard assessment at the local scale is proposed. It is based on travel times and water level height, calculated by the tsunami numerical model, combined with the existence or not of an easy evacuation path from the shoreline to safely high ground and evaluated by field survey. Furthermore, the 1755 Lisbon Tsunami is considered as the worst case scenario, allowing evaluating the impact of a similar scenario at Figueira da Foz municipality, Portugal. The results show that all the beaches are inundated and should be evacuated within an hour after the earthquake. Since there is safely high ground nearby most areas leading to a local tsunami hazard of＂low＂. However, the presence of unstable sand dunes that has been showing signs of collapsing at the south of Cova and Leirosa allowing the tsunami to penetrate inland, inundating the residential areas. For that reason, the local tsunami hazard is ＂moderate＂. The other area which has ＂moderate＂ classification is Cabedelo, because it does not have any coastal protection from tsunami waves, and does not have easy access to the high ground. The marina and fishing port have ＂very low＂ classification, nevertheless it is recommended that vessels evacuate to art offshore area.

Evaluation of Height of Tsunami Induced by Submarine Landslide

Under the new regulatory requirements for nuclear power plants in Japan, which were enacted in response to the nuclear accident associated with the Great East Japan Earthquake Tsunami that occurred on 11 March 2011, it is a requirement to establish a site-specific ＂standard tsunami＂ based on numerical analysis considering non-seismic factors in addition to general seismic faults. It is necessary to establish a consistent evaluation scheme for estimation of tsunami height induced by submarine landslide, since a standard framework for evaluation has not yet been established even though several models for calculation have been proposed and applied in practice. In this study, we estimated the scale of submarine landslide from a literature survey and showed examples of tsunami height evaluation using multiple schemes. As a result of evaluation of tsunami height using three schemes, the Watts model, the KLS model, and the modified-KLS model, the result obtained by the KLS model was comparatively large for every case.

Long-Term Statistics of Extreme Tsunami Height at Crescent City

Historically, Crescent City is one of the most vulnerable communities impacted by tsunamis along the west coast of the United States, largely attributed to its offshore geography. Trans-ocean tsunamis usually produce large wave runup at Crescent Harbor resulting in catastrophic damages, property loss and human death. How to determine the return values of tsunami height using relatively short-term observation data is of great significance to assess the tsunami hazards and improve engineering design along the coast of Crescent City. In the present study, the extreme tsunami heights observed along the coast of Crescent City from 1938 to 2015 are fitted using six different probabilistic distributions, namely, the Gumbel distribution, the Weibull distribution, the maximum entropy distribution, the lognormal distribution, the generalized extreme value distribution and the generalized Pareto distribution. The maximum likelihood method is applied to estimate the parameters of all above distributions. Both Kolmogorov-Smirnov test and root mean square error method are utilized for goodness-of-fit test and the better fitting distribution is selected. Assuming that the occurrence frequency of tsunami in each year follows the Poisson distribution, the Poisson compound extreme value distribution can be used to fit the annual maximum tsunami amplitude, and then the point and interval estimations of return tsunami heights are calculated for structural design. The results show that the Poisson compound extreme value distribution fits tsunami heights very well and is suitable to determine the return tsunami heights for coastal disaster prevention.

An experimental study on tsunami inundation over complex coastal topography

In recent studies, the effects of complex costal topography on tsunami run-up has sparked heated discussion. This study mainly aims at investigating the effects of complex costal topography on the tsunami inundation distance and the effectiveness of sand dunes in dissipating tsunami wave energy. The experiments were carried out in a wave flume to investigate the potential reduction effects of wave run-up by non erodible sand dune like features. The results show that increasing dunes spacing could not significantly affect inundation distance. However, if the height of sand dunes is of the same order of magnitude as the incoming tsunami wave and the gaps between the dunes are large enough, successful tsunami mitigation could also be possible.

Inundation Studies Along the East Coast of India due to Tsunamigenic Earthquakes in North-Andaman and Car Nicobar Subduction zone

The Sumatra-Andaman arc is an active subduction zone and had generated several destructive Tsunamis in the past.In this paper we have analyzed two historical Tsunamigenic earthquakes from this region. One of the historical earthquake is the earthquake of 26th June 1941 in the North Andaman region,which was one of the strongest in the Andaman Sea and Bay of Bengal of magnitude M_w=7.7.This earthquake had triggered tsunami which affected the east coast of India.The other is the earthquake in Car Nicobar region on 31 st December 1881 of magnitude M_w=7.9. This submarine earthquake beneath the

A mathematical model of calculating local tsunami wave source constraints and propagation

This paper presents a local tsunami simulation, including the initial displacement field model of tsunami source and tsunami wave propagation model. We deduced the tsunami wave equation; applied the matching of interior and exterior solutions method and water mass method to determine the initial displacement field in different bottom topography. Tsunami wave propagation model was based on the Boussinesq equation. Difference format was based on the ADI method which discretized in alternating direction in the form of implicit scheme. The open boundary of ADI had been revised considering the influence of wave propagation in the equation of motion. The local tsunami mathematical model was used in the simulation of 2011 Japan tsunami, and the results and the observation data match well.

Effect of the Bottom Topography on Tsunami Propagation in the East (Japan) Sea

A study of tsunami events in the East （Japan） Sea using continuous Galerkin finite element model, aiming at reproducing tsunami waves generated by underwater earthquakes in 1983 and 1993 respectively has been performed focusing on the geographic extent of a topographic feature in the East （Japan） Sea. Numerical models can be the proper tools to study the combined effects of realistic topography. Subsequently, using the FEM based two-dimensional model we have simulated the smoothed and flattened topographic effects by removal of Yamato Rise and seamounts for the cases of tthe 1983 Central region earthquake tsunami and the 1993 southwestern Hokkaido earthquake tsunami. The results have shown that there will be higher tsunamis along the eastern coasts of Korea in general except some areas, like Sokcho with removal of topographic highs, thus providing complicated bottom topography of the East （Japan） Sea as effective tsunami energy scattering.

A Mega-Tsunami in the Baltic Sea 1171 BC: Geological Records with Special Reference to the Lake Mälaren Area in Sweden

At about 3000 C14-year BP or 1200 cal. yrs BC, the Baltic Sea experienced a mega-tsunami with a wave-height of 10 m or more, and a run-up height of up to 16.5 m. This event had significant geological and archaeological effects. We explore the records from the Lake M?laren area in Sweden. The tsunami event is linked to seismic ground shaking and methane venting tectonics at several sites. The triggering factor is proposed to be the Kaali meteor impact in Estonia of the same age. The documentation of a mega-tsunami in the middle of the Bronze Age has wide implications both in geology and in archaeology. The archaeological key sites at Annelund and Apalle are reinterpreted in terms of tsunami wave actions remodelling stratigraphy. By extensive coring, we are able to trace the tsunami effects in both off-shore and on-shore environment. At the time of the event, sea level was at +15 m (due to isostatic uplift). The tsunami wave erosion is traced 13.5 m below sea level. The tsunami run-up over land is traced to +29.5 m to +31.5 m (occasionally even higher), implying a run-up of 14.5 - 16.5 m. In ?ngermanland, the tsunami event was absolutely dated at 1171 varve years BC. Archaeologically, the tsunami event coincides well with the transition between Periods II and III of the South Scandinavian Bronze Age. Period III has traditionally been difficult to identify in the cultural materials of the Lake M?laren region.

Increase in Probability of the Correct Forecast of Tsunami Formation

The presented paper exhibits theory of the ＂gravitational＂ waves propagation near the Earth surface and in the ocean. There was determined an expression ＂gravitational＂ wave which was registered by the gravimeters being placed in several points of the Earth globe. Alteration of gravitational field was accompanied by alteration of the ＂gravitational＂ wave which has the velocity differing from the velocity of seismic waves. The theoretical model was proved by many experiments realized under registration of the underwater earthquake core by tens of gravimeters being placed in the Earth globe different points. The ＂gravitational＂ waves assist to increase the right forecast probability of the beginning tsunami to 50%.

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.

The Historical Review of the 1755 Lisbon Tsunami

The original accounts reporting the 1755 Lisbon Tsunami were compiled, focusing on the descriptions related to the tsunami parameters, damage and fatalities in the affected regions （Atlantic Northeast and Caribbean）. The accounts show the tsunami reached the southwest municipalities of Portugal in less than 30 minutes after the earthquake, and about one hour later the entire coasts of Portugal, Gulf of Cadiz （Spain） and Morocco were hit by the first tsunami waves. The tsunami took about five hours to reach Ireland and UK, and hit the Caribbean region about 9-10 hours after the earthquake. In addition, significant damage and fatalities were reported in Portugal （Peniche, the Lisbon Metropolitan Area and seven coastal areas in the south）, Spain （Cadiz） and in all Morocco coastline, including Marrakech which is located about 200 km inland. In Ireland, UK and the Caribbean minor damage was observed, and no fatalities were reported.

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.

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 absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.

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.

The September 16, 2015 Mw 8.3 Illapel, Chile Earthquake： characteristics of tsunami wave from near-field to far-field

On September 16, 2015, an earthquake with magnitude ofMw 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 from 9.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.