Global tsunami modelling on a spherical multiple-cell grid

Abstract

A model of shallow water equations (SWEs) on a spherical multiple-cell (SMC) grid of 4-level (2.5–5–10–20 km) spatial resolutions is used to simulate tsunami propagation on global ocean surfaces. The unstructured SMC grid retains rectangular cells of the longitude-latitude grid so that efficient finite-difference schemes could be used. It also supports multi-resolutions like mesh refinement to resolve small islands and coastline details while keeping models compact enough to fit into available computers. Two earthquake-induced tsunami cases are simulated and compared with available observations. Results indicate that the modelled tsunami arrival times agree well with observations while tsunami wave heights are underestimated, particularly the observed runups on remote coastal lands. Possible reasons for this underestimation include the smoothing scheme used to suppress numerical oscillations and the missing of initial kinetic energy input from the earthquakes. Another reason is the limitation of the SWEs to describe coastal bores and breaking waves in coastal waters. A possible tsunami scenario induced by a landslide in the Canary Islands is also simulated to assess its potential impact on Atlantic coastal regions. Model results indicate that this kind of tsunami may cause severe damage to local areas but its effects on far fields, like the UK and American coastal regions are small. As the initial landslide disturbance is overly simplified, this study does not give a true representation of a real landslide tsunami but rather a qualitative assessment of its impact on the Atlantic Ocean. More realistic initial condition and improved model representation of coastal processes are needed for further studies of this possible landslide hazard.