Tsunami simulations on several scales

The tsunami event generated by the great Sumatra–Andaman earthquake on 26 December 2004 was simulated with the recently developed model TsunAWI. The model is based on the finite element method, which allows for a very flexible discretization of the model domain. This is demonstrated by a triangulation of the whole Indian Ocean with a resolution of about 14 km in the deep ocean but a considerably higher resolution of about 500 m in the coastal area. A special focus is put on the Banda Aceh region in the Northern tip of Sumatra. This area was heavily hit by the tsunami and the highest resolution in this area is about 40 m in order to include inundation processes in the model simulation. We compare model results to tide gauge data from all around the Indian Ocean, to satellite altimetry, and field measurements of flow depth in selected locations of the Aceh region. Furthermore, we compare the model results of TsunAWI to the results of a nested grid model (TUNAMI-N3) with the same initial conditions and identical bathymetry and topography in the Aceh region. It turns out that TsunAWI gives accurate estimates of arrival times in distant locations and in the same mesh gives good inundation results when compared to field measurements and nested grid results.

[1]  Eric Deleersnijder,et al.  An efficient Eulerian finite element method for the shallow water equations , 2005 .

[2]  Stephen Pond,et al.  A Numerical Model of the Circulation in Knight Inlet, British Columbia, Canada , 1995 .

[3]  Andrew Staniforth,et al.  Finite Elements for Shallow-Water Equation Ocean Models , 1998 .

[4]  Jens Schröter,et al.  A finite-element ocean model: principles and evaluation , 2004 .

[5]  R. Walters,et al.  ASSESSING TSUNAMI HAZARD ALONG THE NEW ZEALAND COAST , 2003 .

[6]  Lori Dengler,et al.  Northwest Sumatra and Offshore Islands Field Survey after the December 2004 Indian Ocean Tsunami , 2006 .

[7]  NUMERICAL MODELING OF THE GLOBAL TSUNAMI , 2005 .

[8]  Dinesh Manocha,et al.  Applied Computational Geometry Towards Geometric Engineering , 1996, Lecture Notes in Computer Science.

[9]  Richard E. Thomson,et al.  The 26 December 2004 Sumatra Tsunami: Analysis of Tide Gauge Data from the World Ocean Part 1. Indian Ocean and South Africa , 2007 .

[10]  R. Walters.  Design considerations for a finite element coastal ocean model , 2006 .

[11]  Chen Ji,et al.  Rupture Process of the 2004 Sumatra-Andaman Earthquake , 2005, Science.

[12]  Paul M. Whitmore,et al.  The Tsunami of 26 December, 2004: Numerical Modeling and Energy Considerations , 2007 .

[13]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[14]  William H. Frey,et al.  Mesh relaxation: A new technique for improving triangulations , 1991 .

[15]  Fred F. Pollitz,et al.  Implications of the 26 December 2004 Sumatra-Andaman earthquake on tsunami forecast and assessment models for great subduction-zone earthquakes , 2007 .

[16]  Rhett Butler,et al.  Global seismographic network records the Great Sumatra‐Andaman earthquake , 2005 .

[17]  S. Tinti,et al.  Finite-element simulations of the 28 december 1908 Messina Straits (Southern Italy) tsunami , 1999 .

[18]  C. Synolakis,et al.  Northern Sumatra Field Survey after the December 2004 Great Sumatra Earthquake and Indian Ocean Tsunami , 2006 .

[19]  V. R. Schneider,et al.  GUIDE FOR SELECTING MANNING'S ROUGHNESS COEFFICIENTS FOR NATURAL CHANNELS AND FLOOD PLAINS , 1989 .

[20]  B. G. McAdoo,et al.  Inundation distances and run‐up measurements from ASTER, QuickBird and SRTM data, Aceh coast, Indonesia , 2007 .

[21]  Jörn Behrens,et al.  TsunAWI - Unstructured Mesh Finite Element Model for the Computation of Tsunami Scenarios with Inundation , 2008 .

[22]  Stephan T. Grilli,et al.  Source Constraints and Model Simulation of the December 26 , 2007 .

[23]  Yutaka Hayashi,et al.  The 2004 Indian Ocean tsunami: Tsunami source model from satellite altimetry , 2006 .

[24]  Fred F. Pollitz,et al.  Coseismic Slip Distributions of the 26 December 2004 Sumatra–Andaman and 28 March 2005 Nias Earthquakes from gps Static Offsets , 2007 .

[25]  Vasily Titov,et al.  Implementation and testing of the Method of Splitting Tsunami (MOST) model , 1997 .

[26]  S. Tinti,et al.  A smoothing algorithm to enhance finite-element tsunami modelling: An application to the 5 February 1783 Calabrian case, Italy , 1995 .

[27]  Y. Okada.  Surface deformation due to shear and tensile faults in a half-space , 1985 .

[28]  Vasily Titov,et al.  The Global Reach of the 26 December 2004 Sumatra Tsunami , 2005, Science.

[29]  Franklin Liu,et al.  Modeling wave runup with depth-integrated equations , 2002 .

[30]  A. K. Lal,et al.  The Great Tsunami of 26 December 2004: A description based on tide-gauge data from the Indian subcontinent and surrounding areas , 2006 .

[31]  Extracting the 2004 Indian Ocean tsunami signals from sea surface height data observed by satellite altimetry , 2008 .

[32]  D. Kapoor.  General bathymetric chart of the oceans (GEBCO) , 1981 .

[33]  S. Kathiroli,et al.  Rupture process of the 2004 great Sumatra-Andaman earthquake estimated from tsunami waveforms , 2005 .