| Submarine landslides can generate tsunamis, and the generated waves can be catastrophic when a large volume of landslide material is involved. Moreover, large earthquakes are often accompanied by submarine landslides that can enhance the magnitude of the resulting tsunamis. In this thesis, numerical schemes are developed to solve the wave propagation problems generated by submarine landslides. Assuming the landslides in a flow regime, depth-averaged models are studied, and finite volume methods are extended to the fully coupled multi-layer shallow water equations. From the fully coupled model, an efficient simplified approach is derived that is often appropriate for tsunamis generated by submarine landslides. These waves can have relatively short wavelength, and another class of equations may be necessary that can handle the dispersion of waves. Several types of the Boussinesq equations have been reviewed and implemented with a hybrid of high-resolution finite volume and finite difference methods. Stability analysis and convergence tests have been performed for the hybrid scheme. The develpoment has been done in the context of the textsc{Geoclaw} framework, a code designed to handle the single-layer shallow water equations, that uses adaptive mesh refinement to model tsunami propagation on a global scale with inundation of specific regions on a fine grid. The newly developed methods, tested on the exact solutions, are validated by comparing to laboratory experiments and by applying to historic events such as the Papua New Guinea 1998 and Storegga slides. Possible scenarios of submarine landslides and resulting tsunamis on the Washington coast were investigated. |
