Submarine Landslides Research Articles

Submarine landslides and tsunami genesis in Sagami Bay, Japan, caused by the 1923 Great Kanto earthquake

The 1923 Great Kanto earthquake occurred on September 1, in Japan, and caused severe damage mainly in the Kanto region. Tsunamis were observed over wide regions from the east coast of the Izu Peninsula to the west coast of the Boso Peninsula, and particularly, the damage in Atami was devastating. Many earthquake fault models including those of Kanamori (1971) and Ando (1971) were proposed based on the records of land deformation. However, such fault models cannot sufficiently explain the tsunami elevation and its initial sea-level motion on the coasts. Hence, the detailed mechanisms remain elusive. This study examines the possibility that a leading mechanism of the tsunami in the 1923 Great Kanto earthquake was a large-scale submarine landslide that occurred at Sagami Bay and at the mouth of Tokyo Bay based on the records of depth data measured by the Imperial Japanese Navy (1924) before and after the earthquake. We first show that the tsunami calculated by each fault model was inconsistent with the waveform at Yokosuka, the coastal tsunami elevations, and initial sea-level motion. Then, based on statistical analysis of the depth changes at the 1923 Great Kanto earthquake, we found that the seafloor bathymetric changes represented large-scale submarine landslides that may correspond to long-runout submarine liquefied sediment flows. The seafloor gradient over a 40 km flow-out distance was equal to or less than 0.4∘\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^\\circ$$\\end{document}. Through the identification of the submarine landslide source by tsunami backpropagation analysis and utilizing an analytical solution of a high-density gravity flow and a sensitivity analysis, we conducted a range of numerical simulations of the 1923 Great Kanto earthquake tsunamis using a fault model and a submarine landslide tsunami source model due to a high-density liquefied gravity flow. The results quantitatively accounted for the discrepancy between the observed tsunami records with maximum tsunami elevations over 12 m and the fault-model–based simulations with maximum tsunami elevations of 2 to 5 m and explained consistently the maximum tsunami elevation distributions as well as the time-series tsunami waveforms. These results may thus facilitate and deepen our understanding of the earthquake-induced submarine landslide tsunami risk as cascading multi-geohazards.

A Mixture Model With Slip Velocity for Saturated Granular‐Liquid Free‐Surface Flows

AbstractIn this paper, a model is presented for modeling saturated granular‐liquid free‐surface flows, in which the volume‐averaged mixture bulk velocity is employed to derive the balance equations for the mass and momentum of mixture flow. Additionally, an evolution equation of the slip velocity between granular‐and liquid constituents is derived to describe the separation between these constituents. The frictional‐collisional constitutive relation for granular‐constituent is employed to determine the stress due to particles interaction. The governing equations for mixture flows are numerically solved by a finite difference two‐step projection method. The volume of fluid (VOF) method is employed to track the free surface of the mixture flow in the present numerical model. Good agreements between numerical results and experimental data are observed by modeling the dam‐break process of granular‐liquid mixture flow, dam‐break waves over the saturated erodible beds and surge waves induced by submarine landslides along an inclined plane. Furthermore, the difference between the volume‐averaged mixture bulk velocity and mass‐averaged mixture bulk velocity is found to vary as the instinct density ratio of granular‐constituent and liquid‐constituent and the volumetric concentration ns of the granular‐constituent, and the evolution in the slip velocity during the process of the settlement of sediments is numerically analyzed.

Can Glacial Sea‐Level Drop‐Induced Gas Hydrate Dissociation Cause Submarine Landslides?

AbstractWe conducted two‐dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea‐level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow‐water portion of the gas hydrate stability zone (GHSZ). This shallow‐water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low‐latitude regions. The deep BGHSZ (>1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large‐scale submarine landslides. Observed deep‐water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch‐out zone.

An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions, but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs. To overcome this limitation, a message passing interface (MPI) parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency. This framework utilises a static domain decomposition scheme, with the entire computation domain being decomposed into multiple subdomains according to predefined processors. A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation. In particular, a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces. Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework. Subsequently, the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides. The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.

Impact of initial slide on downslope seafloor with a weak layer

Submarine landslides are a major offshore geohazard posing increasing threat to offshore and coastal developments. A submarine landslide may usually involve multistage failure triggered by the first failure. During submarine landslide evolution, the emergent sliding mass may act as an impact load on the downslope seafloor potentially leading to surficial erosion of sediments or deeper ploughing failure with the presence of a weak layer. The study presents a comprehensive analysis of seafloor instability and post-failure behaviors subjected to the first slide mass impact, using a large deformation finite element (LDFE) method. Different failure patterns are investigated and the corresponding critical impact loads for triggering the seafloor instability are analyzed. The parameters controlling the strain softening behaviors of soils and the undrained shear strength in the weak layer have significant effects on the critical condition for the seafloor failure. The dynamic inertia of the slide mass, however, has little effect on the results. An empirical equation for assessing the seafloor instability subjected to the first slide impact, with respect to the impact load from the initial slide and the gravity stress relative to the shear strength in the weak layer, is proposed based on a parametric study.

An environmental-benign method to enhance hydrate production in clayey silty reservoir via splitting grouting

During the pivotal period of energy transition between carbon neutrality and carbon peaking, it is crucial to develop the natural gas hydrate (NGH) due to the advantages of high energy density and non-pollution. However, the NGH production also faces several significant environmental challenges that need to be addressed such as seawater intrusion, methane leakage, and submarine landslides. Therefore, an environmental-benign method to enhance hydrate production in clayey silty reservoir via splitting grouting was introduced. As a novel method for achieving the two goals of both reservoir skeleton reinforcement and permeability enhancement, is proposed to stimulate NGH reservoir. Herein, the formulation based on foamed ultrafine cement-based grout was established. The results showed that the grout has good fluidity with initial viscosity ∼500 mPa·s and consolidation time in the range of 28–50 min. Meanwhile, the consolidation performance evaluation was studied, and it showed that the consolidation is with porosity in the range of 9.41%–46.59%, permeability in the range of 51.23 mD-789.07 mD, and compressive strength in the range of 3.31 MPa–3.88 MPa. The internal pore channels are uniformly distributed with a large number of stomatal strings, which can be used as a multistage seepage channel for gas and water production from hydrate decomposition. To further examine the stimulating behavior, the effect of hydrate saturation, grouting well diameter, grouting flow rate, and in-situ stress on splitting pressure and grouting volume were analyzed, and it was found that the grout has good ability to split fracturing and grout diffusion under the conditions of axial pressure of 5 MPa, confining pressure of 4 MPa, and hydrate saturation not more than 0.30. In addition, an environmental-benign high-efficient reservoir stimulation method on a laboratory scale was proposed. The findings could provide fundamental insights and recommendations for the implementation of splitting grouting in NGH reservoirs.

A Comprehensive Assessment of Submarine Landslides and Mass Wasting Processes Offshore Southern California

AbstractIt is critical to characterize submarine landslide hazards near dense coastal populations, especially in areas with active faults, which can trigger slope failure, subsequent tsunamis, and damage seabed infrastructure during earthquake shaking. Offshore southern California, numerous marine geophysical surveys have been conducted over the past decade, and high‐resolution bathymetric and subsurface data now cover about 60 percent of the total region between Point Conception and the United States‐Mexico border from the California coast out to the base of Patton Escarpment ∼200 km offshore. In a comprehensive compilation and interpretive mapping effort, we find evidence of seafloor failure throughout offshore southern California with nearly 1,500 submarine landslide‐related features, including 63 discrete slide deposits with debris and >1,400 slide‐related scarps. In our analysis, we highlight new mapping of submarine landslides in Catalina Basin, the Del Mar slide, the San Gabriel slide complex, and the 232 km2 San Nicolas slide, the largest area of any known submarine landslide mass offshore southern California. Analysis of the spatial distribution of submarine landslide features suggests that most mapped slide features are located relatively near coastal sediment sources, particularly during sea‐level lowstand conditions, which underscores the importance of sediment supply and sediment accumulation on low‐gradient slopes as failure preconditioning processes. Tectonically driven uplift at shelf edges and along basin flanks is another key preconditioning factor, and our results also suggest that earthquakes along active faults trigger mass wasting, especially for repeated, small‐scale failures on tectonically steepened slopes.

How does rapid burial work? New insights from experiments with echinoderms

AbstractThis research explores the significance of rapid burial in preserving fossils, with a particular focus on free‐living echinoderms. Experiments were based on ophiuroids to simulate burial under different turbiditic flows. The results showed that a bed thickness of around 10 cm is a limit for the preservation of whole skeletons in most cases. The type of sediment can affect the integrity of the buried skeletons, with sand deposition resulting in higher rates of autotomy. However, mud deposition did not show any numbing effect, as previously believed for echinoderms. In contrast, freshwater‐rich sediments can play a critical role, paralysing specimens and preventing escape postures through rapid changes in salinity. From this, the study highlights the importance of extrabasinal turbidites, generated outside the marine basin, in the fossilization of marine invertebrates. Such sediments are rich in fresh water and can be more efficient burial traps compared to other intrabasinal deposits generated by storm waves or submarine landslides.

NUMERICAL MODELING OF THE 1998 PAPUA NEW GUINEA TSUNAMI USING THE COMCOT

The Papua New Guinea tsunami of 1998 is a unique phenomenon because the source of the tsunami propagation has been speculated. There was a 7.1-magnitude earthquake on July 17, 1998, at 18:49 WIT before the tsunami hit the Aitape area. However, previous studies have shown that the leading cause of the tsunami was not the earthquake but a submarine landslide. One of the steps to simulating the event is to do tsunami modeling. A tsunami propagation simulation will be conducted using Cornell Multi-grid Coupled Tsunami (COMCOT). This simulation was carried out with three scenarios to see which had the most significant effect on the tsunami event. The first scenario uses a tsunami source from a 7.1 magnitude earthquake, the following scenario is carried out using avalanche parameters, and the last scenario is a scenario with a combined source of earthquake and avalanche. The results of this study indicate that underwater landslides are the source of a tsunami similar to the original event.

Large-scale submarine landslides in the Barberton Greenstone Belt, southern Africa—Evidence for subduction and great earthquakes in the Paleoarchean

Abstract New mapping of the Barberton Greenstone Belt in South Africa shows that the central part is a pseudo-stratigraphy made of shallow-water and deep-water siliciclastic and volcanic slide blocks, with individual blocks ranging in size from tens of meters to >10 km in length. The outcrop pattern and scale are remarkably similar to those of large-scale Miocene to recent submarine landslides in New Zealand along the active Hikurangi subduction zone that are periodically triggered by earthquakes on the subduction megathrust, providing evidence for megathrust earthquakes in the Paleoarchean.

Offshore application of landslide susceptibility mapping using gradient-boosted decision trees: a Gulf of Mexico case study

Among natural hazards occurring offshore, submarine landslides pose a significant risk to offshore infrastructure installations attached to the seafloor. With the offshore being important for current and future energy production, there is a need to anticipate where future landslide events are likely to occur to support planning and development projects. Using the northern Gulf of Mexico (GoM) as a case study, this paper performs Landslide Susceptibility Mapping (LSM) using a gradient-boosted decision tree (GBDT) model to characterize the spatial patterns of submarine landslide probability over the United States Exclusive Economic Zone (EEZ) where water depths are greater than 120 m. With known spatial extents of historic submarine landslides and a Geographic Information System (GIS) database of known topographical, geomorphological, geological, and geochemical factors, the resulting model was capable of accurately forecasting potential locations of sediment instability. Results of a permutation modelling approach indicated that LSM accuracy is sensitive to the number of unique training locations with model accuracy becoming more stable as the number of training regions was increased. The influence that each input feature had on predicting landslide susceptibility was evaluated using the SHapely Additive exPlanations (SHAP) feature attribution method. Areas of high and very high susceptibility were associated with steep terrain including salt basins and escarpments. This case study serves as an initial assessment of the machine learning (ML) capabilities for producing accurate submarine landslide susceptibility maps given the current state of available natural hazard-related datasets and conveys both successes and limitations.

Origin, evolution and significance of giant buried sediment mounds near the Sahara Slide Complex, North‐west African margin

AbstractMixed turbidite–contourite depositional systems result from interactions between down‐slope turbidity currents and along‐slope bottom currents, comprising excellent records of past oceanographic currents. Modern and ancient systems have been widely documented along the continental margins of the Atlantic Ocean. Yet, few examples have so far been identified on the North‐west African continental margin, limiting understanding of the sedimentary and palaeoceanographic evolution in this area. This work uses two‐dimensional seismic reflection profiles to report, for the first time, the presence of three giant sediment mounds beneath the headwall region of the Sahara Slide Complex. The sediment mounds are elongated and separated by two broad canyons, showing a north‐west/south‐east orientation that is roughly perpendicular to the continental margin. These mounds are 24 to 37 km long and 12 to 17 km wide, reaching a maximum height of ca 1000 m. Numerous slide scarps are observed within and along the flanks of the mounds, hinting at the occurrence of submarine landslides during their development. Based on their geometries, external shapes, internal seismic architecture and stratigraphic stacking patterns, it is proposed that these sediment mounds comprise down‐slope elongated mounded drifts formed in a mixed turbidite–contourite system during four evolutionary stages: onset, growth, maintenance and burial. The significance of this work is that it demonstrates the gradual transition from a turbidite system to a full mixed turbidite–contourite system to be associated, in the study area, with the establishment of strong ocean currents along north‐west Africa.

Lithological features of the uryuk formation in the tolparovo section of the vendian (Southern Urals)

New results of lithological structural and textural observations and petrographic study of rocks of the Uryuk formation exposed in the basin of the Maly Tolpar River in the Bashkir meganticlinory are presented. For the first time, a wide development of consedimentary underwater landslide folds in the Uryuk deposits has been established, the sizes of which vary from the first tens of centimeters to the first meters. In sandstones, the presence of oblique, wavy, flasery and horizontal, with a dual character of thin clay layers, layering, as well as structureless layers with a massive texture, is recorded. Signs of shallow wave ripples have been established. Mechanoglyphs and textures similar to Arumberia banksi have been found, which presumably represent the lithified structures of bacterial mats that existed under certain conditions – shallow-sea environments and conditions of sandy-clay sedimentation. It is concluded that the Uryuk deposits, despite the absence of direct lithological signs of diamictites in them, form a single sedimentary sequence with the underlying Tolparov-Suirovsky marine glacial deposits, in which a shallow-sea sedimentation regime is established by the end of the Uryuk time. It is shown that the presence of ferruginous minerals in sandstone cement is secondary in nature and is associated with their epigenetic transformations. It is concluded that the red-colored color of the rocks of the Uryuk formation cannot be used for stratification of its sections and reconstruction of sedimentation conditions in the Uryuk time.

Submarine‐channel meandering reset by landslide filling, Taranaki Basin, New Zealand

AbstractLandslides are among the largest mass movements on Earth. As such, the deposits of landslides, also known as mass‐transport deposits, are significant architectural elements of continental margins, especially those receiving sediment from large deltas. Landslide dams have been shown to alter the courses of rivers and submarine channels. However, there are fewer examples of landslides completely filling submarine channels and examples of the subsequent stratigraphic evolution. A three‐dimensional seismic‐reflection dataset (<90 Hz) from the deep‐water (>1500 m) Taranaki Basin, offshore the North Island of New Zealand, was used to explore the response of a sequence of channel deposits to landslide filling. The basal channel system initially meandered like a river, with successive channel positions in close proximity, as it aggraded >250 ms two‐way travel time. This systematic, organised evolution is governed by the memory of early channel evolution, which sets the sea floor geomorphology that guides channel‐forming turbidity currents. Later, a channel approximately twice as wide as underlying channels cut off a number of channel bends, probably as a result of an increase in the discharge of channel‐forming turbidity currents. This last channel was filled with submarine landslides, which transported and deposited sediment as debris flows based on the presence of blocks within a matrix comprising chaotic, lower amplitude seismic facies. These debris‐flow deposits smoothed over the sea floor, effectively wiping the memory of channel evolution. As a result, the subsequent channel pattern bears no resemblance to the basal system. Submarine‐channel resetting by landslide filling is common in settings with frequent catastrophic basin‐margin collapses, like offshore New Zealand.

Applying cellular automata to dynamic simulation of a tsunamigenic submarine landslide in the South China Sea

This work presents the first known application of a new approach to tsunami modelling, by linking a discrete cellular automata (CA) dynamic model of a submarine landslide to an off-the-shelf tsunami modelling package. This has the advantage of simulating the temporal evolution of the submarine landslide rather than as a single rigid sliding mass. We first tested that the coupled model was able to satisfactorily reproduce observed tsunami wave heights resulting from the flank collapse of the Anak Krakatau volcano in the Sunda Strait (Indonesia) in December 2018, before applying it to a speculative area of interest on the North Borneo (Continental) Shelf (NBS). The NBS has one of the thickest near-shore clastic sediment wedges known globally (12 km thickness) and a very long pre-history of submarine landsliding. We modelled a small slide (5 km3 volume) because this is far more likely to occur than a recurrence of the gigantic Brunei slide (1200 km3 volume). Results indicate that a 5 km3 submarine landslide generates multiple tsunami waves, the largest of which reached 8 m at the coast of nearby Balabac Island (The Philippines). Small tsunami waves also arrived on the central Vietnam coast, more than 900 km distant. A larger potential failure on the NBS therefore appears capable of causing a tsunami with greater runup heights, which could pose a risk to much of central and southern Vietnam, northern Borneo (Sabah, Brunei) and Palawan. Finally, we note that the South China Sea is fringed by a necklace of other thick offshore clastic sediment accumulations, most notably on the southern China continental margin, offshore southern Taiwan, and offshore central Vietnam. Given the proof-of-concept of our coupled CA–tsunami modelling approach described here, we recommend this also be applied in these areas to assess the risk posed by a wide range of submarine landslides to coastal populations and infrastructural assets (including nuclear power plants).

Investigation of offshore landslides impact on bucket foundations using a coupled SPH–FEM method

BackgroundAs the exploitation of marine resources intensifies, the impact of submarine landslides on underwater structures has become a significant issue. Existing research primarily focuses on the impact on pipelines, often neglecting the actual deformation and mechanical response of underwater structures under impact loads in numerical simulations, thus complicating the evaluation of the reliability of these engineering structures in extreme conditions. Moreover, the dynamic response of bucket foundations, a common form of underwater base, under the effect of submarine landslide impacts remains unclear.MethodsTo address this knowledge gap, we have developed a fluid-structure coupling system that employs the coupled Smoothed Particle Hydrodynamics (SPH)-Finite Element Method (FEM) to investigate a single impact process and analyze the displacement response of bucket foundations within a water-offshore landslide-bucket foundationsubgrade context. The accuracy of this developed method has been systematically verified through comparisons with previous experimental and numerical results.ResultsDuring a submarine landslide impact event, the impact force demonstrates a distinct decrease followed by stabilization, and the displacement response of the bucket foundation exhibits a rebound effect after reaching its maximum value. Furthermore, we conducted an extensive analysis of different impact angles for underwater data centers equipped with multi-bucket foundations. Our study revealed that group-bucket foundations experience a combined translation-turnover failure when subjected to submarine landslide impacts, and the most unfavorable scenario for such impact is identified. The research introduces a novel numerical simulation approach for investigating the impact of submarine landslides on complex underwater structures.

National Archive of Marine Seismic Surveys Web Portal Provides Public Access to US Exclusive Economic Zone Marine Seismic Surveys

AbstractSeismic reflection data collected within the US Exclusive Economic Zone (EEZ) represent a scientific resource with inestimable value. The US Geological Survey (USGS), in partnership with the Bureau of Ocean Energy Management (BOEM), provide unrestricted public access to seismic data acquired since 1975 through the National Archive of Marine Seismic Surveys (NAMSS) portal. This map‐searchable web portal provides access to over 1,100 two‐ and three‐dimensional (2D, 3D) marine seismic data sets collected by the USGS and by the offshore resource industry. The industry data sets are released by BOEM under The Code of Federal Regulations 30 CFR Parts 551 and 580, allowing for the release of seismic reflection data after the completed 25‐year proprietary term from the issuance of the Geological and Geophysical (G&G) permit. The web portal provides tools to filter and identify surveys within the US EEZ and navigation and shapefiles can be downloaded with each survey. Seismic data sets downloaded from NAMSS have been instrumental in numerous scientific investigations along the US EEZ. These include studies focused on offshore marine geohazards, such as fault systems and submarine landslides, as well as gas hydrates and other geologic processes. Acquisition of such data sets, especially pneumatic‐sourced surveys, is becoming progressively more difficult due to environmental concerns and associated costs, thus preservation and public access to such high‐quality data sets is critical to maintaining the ability to conduct offshore seismic reflection scientific research along the US EEZ.

Study on the geomorphological changes of deep troughs under the influence of reclamation in the Caofeidian

Reclamation projects generally change the hydrodynamics and thus affect the topography of the seabed, which is the main threat to the safety of marine engineering. The Caofeidian Reclamation Project, the largest reclamation project in China to date, involved the reclamation of coastal regions of the Hebei Province. However, the project also created several issues and flaws due to the size of the reclaimed area. Through the comparative analysis in the Caofeidian sea area of multi-phase bathymetric topography, shallow stratigraphic changes, and numerical simulation of sediment dynamics, the characteristics of surface sediments and spatial distribution of erodibility were summarized. New marine geological hazards were discovered, such as submarine landslides and erosion depressions. It was concluded that after the reclamation of the Caofeidian, the south of the deep trough was still in an erosion state, and the erosion rate increased slightly. The shear stress formed by strong tidal currents in the deep trough was greater than the critical erosion shear stress of surface sediments, which caused the erosion transport of sediments. It was the primary dynamic mechanism leading to the long-term erosion state of the deep trough.

Smoothed particle hydrodynamics approach for modelling submerged granular flows and the induced water wave generation

In this study, we develop a Smoothed Particle Hydrodynamics (SPH) 2D-model for simulating fully submerged granular flows and their arising water waves. The granular particles are characterised by a non-Newtonian flow pattern, following a Casson constitutive law, generalised by applying the infinitesimal strain theory to avoid numerical singularities inherited from the original law. The implementation of this rheological model on the weakly compressible viscous Navier-Stokes equations enables the simultaneous modelling of the motion of granular flows and their resulting water waves, establishing a monolithic representation of fluid-structure coupling. The novelty of this model lies in the numerical continuity of the generalised rheological model based mainly on the yield stress criterion, which is computed purely from the mechanical properties of granular materials, including internal friction, cohesion, and viscosity coefficients. The proposed SPH model is validated through two benchmarks available in the literature, representing a submarine landslide along an inclined plane and an immersed granular column collapse. The outcomes of our study illustrate the effectiveness of the proposed model in accurately predicting the motions of submerged granular masses and their resulting water waves, which is crucial for accurately predicting the behaviour of underwater landslides and other natural hazards.

A STUDY ON THE APPLICATION CONDITIONS OF THE INITIAL WATER LEVEL PREDICTION EQUATIONS FOR TSUNAMI ASSOCIATED WITH SUBMARINE LANDSLIDE

A STUDY ON THE APPLICATION CONDITIONS OF THE INITIAL WATER LEVEL PREDICTION EQUATIONS FOR TSUNAMI ASSOCIATED WITH SUBMARINE LANDSLIDE