Adaptive Mesh Refinement Applied to Tsunami Modeling: TsunaFLASH

After the devastating Andaman-Sumatra Tsunami of December 2004, the international community has strived to introduce measures to prevent hazards from future tsunamis as much as possible. In this endeavor numerical modeling plays a key role, since forecasts as well as inundation mapping efforts rely on numerical simulation. In this presentation, we introduce a new triangle-based adaptive mesh finite element numerical model for tsunami propagation (and inundation) simulations. TsunaFLASH combines numerical methods developed in the framework of the unstructured triangular, yet non-adaptive, tsunami model TsunAWI with adaptive mesh refinement capabilities provided by the mesh refinement library amatos. Adaptive methods are well suited for accurate resolution of localized features, maintaining computational efficiency in terms of the number of computations and the required memory. This presentation will introduce the first developments of TsunaFLASH, which use bathymetry and topography data derived from ETOPO5 and focus on the Indian Ocean region. The finite element based discretization scheme with conforming and non-conforming linear elements will be introduced, as well as the coupling with the software RuptGen (Babeyko, 2007), which generates the initial uplift function.

[1]  G. Mellor.  USERS GUIDE for A THREE-DIMENSIONAL, PRIMITIVE EQUATION, NUMERICAL OCEAN MODEL , 1998 .

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

[3]  António M. Baptista,et al.  An Efficient and Robust Tsunami Model on Unstructured Grids. Part I: Inundation Benchmarks , 2008 .

[4]  Teh-Ru Alex Song.  Broadband Modeling of Earthquake Sources and Mantle Structures , 2007 .

[5]  Hai Cheng,et al.  Source parameters of the great Sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls , 2006 .

[6]  B. C. Papazachos,et al.  GLOBAL RELATIONS BETWEEN SEISMIC FAULT PARAMETERS AND MOMENT MAGNITUDE OF EARTHQUAKES , 2004 .

[7]  Stephen Roberts,et al.  Hydrodymamic modelling of coastal inundation , 2005 .

[8]  Harry Yeh,et al.  Advanced numerical models for simulating tsunami waves and runup , 2008 .

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

[10]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

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

[12]  Natalja Rakowsky,et al.  amatos: Parallel adaptive mesh generator for atmospheric and oceanic simulation , 2005 .

[13]  Barbara Romanowicz,et al.  Slip of the 2004 Sumatra-Andaman Earthquake from Joint Inversion of Long-Period Global Seismic Waveforms and GPS Static Offsets , 2007 .

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

[15]  Jens Zimmermann,et al.  Parallelizing an Unstructured Grid Generator with a Space-Filling Curve Approach , 2000, Euro-Par.

[16]  Costas E. Synolakis,et al.  NOAA Technical Memorandum OAR PMEL-135 STANDARDS, CRITERIA, AND PROCEDURES FOR NOAA EVALUATION OF TSUNAMI NUMERICAL MODELS , 2007 .

[17]  Nobuo Shuto,et al.  Numerical simulation of tsunamis — Its present and near future , 1991 .

[18]  Richard Tiley,et al.  The elastic thickness of the British Isles , 2003, Journal of the Geological Society.

[19]  Paul D. Williams.  A Proposed Modification to the Robert–Asselin Time Filter* , 2009 .

[20]  Michael L. Gittings,et al.  Two- and Three-Dimensional Simulations of Asteroid Ocean Impacts , 2002 .

[21]  H. Kanamori.  The energy release in great earthquakes , 1977 .

[22]  Jens Schröter,et al.  Simulation of the minor tsunami generated by the September 30/2009 earthquake near Padang, Sumatra , 2009 .

[23]  Hermann M. Fritz,et al.  Extreme runup from the 17 July 2006 Java tsunami , 2007 .

[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]  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 .

[26]  Chen Ji,et al.  Coseismic Slip and Afterslip of the Great Mw 9.15 Sumatra–Andaman Earthquake of 2004 , 2007 .

[27]  Tokutaro Hatori.  Tsunami magnitudes in Taiwan, the Philippines and Indonesia , 1996 .

[28]  Wolfgang Hiller,et al.  Evolution of Small-Scale Filaments in an Adaptive Advection Model for Idealized Tracer Transport , 2000 .

[29]  Jörn Behrens.  An adaptive semi-lagrangian advection scheme and its parallelization , 1996 .

[30]  Jörn Behrens,et al.  Comparison of approaches with unstructured meshes and nested grids , 2022 .

[31]  Jeffrey Park,et al.  Sumatra-Andaman earthquake , 2006 .

[32]  R. Asselin,et al.  Frequency Filter for Time Integrations , 2003 .

[33]  Chen Ji,et al.  Source Description of the 1999 Hector Mine, California, Earthquake, Part II: Complexity of Slip History , 2002 .

[34]  Pierre Henry,et al.  The Sumatra subduction zone: A case for a locked fault zone extending into the mantle , 2004 .

[35]  Quentin F. Stout,et al.  Adaptive Grids for Weather and Climate Models , 2004 .

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

[37]  Kenji Satake,et al.  Tsunami Source of the 2004 Sumatra–Andaman Earthquake Inferred from Tide Gauge and Satellite Data , 2007 .

[38]  Widodo Setiyo Pranowo,et al.  Modelling of the Bengkulu Minor Tsunami Event, September 12, 2007, West of Sumatra, Indonesia , 2008 .

[39]  Guan-Yu Chen,et al.  Numerical Study on Tsunamis Excited by 2006 Pingtung Earthquake Doublet , 2008 .

[40]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[41]  Natasha Flyer,et al.  Transport schemes on a sphere using radial basis functions , 2007, J. Comput. Phys..

[42]  E. Asphaug,et al.  Impact tsunami Eltanin , 2002 .

[43]  Natalja Rakowsky,et al.  A parallel adaptive barotropic model of the atmosphere , 2007, J. Comput. Phys..

[44]  Fumihiko Imamura,et al.  Source characteristics of the Nicaraguan Tsunami Earthquake of September 2, 1992 , 1993 .

[45]  Martin Vallée,et al.  Rupture Properties of the Giant Sumatra Earthquake Imaged by Empirical Green's Function Analysis , 2007 .

[46]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Transport Scheme on the Cubed Sphere , 2005 .

[47]  Roberto Tonini,et al.  Preliminary Numerical Simulations of the September 12, 2007 Southern Sumatra Tsunami: Forward Modeling and Comparison With Publicly Available Tsunami Data , 2007 .

[48]  Naoyuki Kato,et al.  Hypocenter depths of large interplate earthquakes and their relation to seismic coupling , 2003 .

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

[50]  Kenji Satake,et al.  Source of the July 2006 West Java tsunami estimated from tide gauge records , 2006 .

[51]  Jrn Behrens.  Adaptive Mesh Generator for Atmospheric and Oceanic Simulations - amatos , 2004 .

[52]  C. Chapman.  The hazard of near-Earth asteroid impacts on earth , 2004 .

[53]  Michael Griebel,et al.  Hash based adaptive parallel multilevel methods with space-filling curves , 2002 .

[54]  Jörn Behrens,et al.  Credible Worst Case Tsunami Scenario Simulation for Padang , 2009 .

[55]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[56]  Roy A. Walters,et al.  A semi‐implicit finite element model for non‐hydrostatic (dispersive) surface waves , 2005 .

[57]  D. H. Peregrine,et al.  Surf and run-up on a beach: a uniform bore , 1979, Journal of Fluid Mechanics.

[58]  B. Keating,et al.  INTRODUCTION TO A TSUNAMI-DEPOSITS DATABASE , 2008 .

[59]  H. Kanamori.  Mechanism of tsunami earthquakes , 1972 .

[60]  Hiroo Kanamori,et al.  The physics of earthquakes , 2001 .

[61]  Naoshi Hirata,et al.  Destructive Earthquake and Disastrous Tsunami in the Indian Ocean, What Next? , 2005 .

[62]  R. Dmowska,et al.  Local Tsunamis and Distributed Slip at the Source , 1999 .

[63]  Chen Ji,et al.  Source Description of the 1999 Hector Mine, California, Earthquake, Part I: Wavelet Domain Inversion Theory and Resolution Analysis , 2002 .

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

[65]  Anthony Sladen,et al.  Numerical Modeling of the Great 2004 Indian Ocean Tsunami: Focus on the Mascarene Islands , 2007 .

[66]  André Robert,et al.  The Integration of a Low Order Spectral Form of the Primitive Meteorological Equations , 1966 .

[67]  Wm. G. Van Dorn.  A low‐frequency microbarograph , 1960 .

[68]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[69]  Hiroo Kanamori,et al.  The diversity of the physics of earthquakes , 2004, Proceedings of the Japan Academy. Series B, Physical and Biological Sciences.

[70]  Kenji Satake,et al.  Effects of bathymetry on tsunami propagation: Application of ray tracing to tsunamis , 1988 .

[71]  John M. Dennis,et al.  A Comparison of Two Shallow-Water Models with Nonconforming Adaptive Grids , 2008 .

[72]  Efim Pelinovsky,et al.  Study of tsunami propagation in the Ligurian Sea , 2001 .

[73]  S. Yabushita,et al.  On the possible hazard on the major cities caused by asteroid impact in the Pacific ocean , 1994 .

[74]  Mader,et al.  TSUNAMI FROM ASTEROID AND COMET IMPACTS : THE VULNERABILITY OF EUROPE , 1999 .

[75]  A. Gusman,et al.  Source Model of the 2007 Bengkulu Earthquake Determined from Tsunami Waveform Analysis , 2008 .

[76]  B. Sindhu,et al.  Improved bathymetric datasets for the shallow water regions in the Indian Ocean , 2007 .

[77]  Robert J. Geller,et al.  Scaling relations for earthquake source parameters and magnitudes , 1976 .

[78]  Ralf Ludwig,et al.  Hazard analysis and estimation of people exposure as contribution to tsunami risk assessment in the West Coast of Sumatra, the South Coast of Java and Bali , 2010 .

[79]  A. Snider,et al.  Atmospheric and ocean modeling with an adaptive finite element solver for the shallow-water equations , 1986 .

[80]  A. Bennett,et al.  TOPEX/POSEIDON tides estimated using a global inverse model , 1994 .

[81]  C. Synolakis,et al.  The Runup of Long Waves , 1986 .

[82]  Stefano Lorito,et al.  Rupture Process of the 2004 Sumatra-Andaman Earthquake from Tsunami Waveform Inversion , 2007 .

[83]  C. E. Synolakis,et al.  Validation and Verification of Tsunami Numerical Models , 2008 .

[84]  E. Asphaug,et al.  Asteroid impact tsunami of 2880 March 16 , 2003 .

[85]  R. LeVeque,et al.  FINITE VOLUME METHODS AND ADAPTIVE REFINEMENT FOR GLOBAL TSUNAMI PROPAGATION AND LOCAL INUNDATION. , 2006 .

[86]  Andreas Hoechner,et al.  Enhanced GPS inversion technique applied to the 2004 Sumatra earthquake and tsunami , 2008 .

[87]  Louis W. Uccellini,et al.  The Effects of the Asselin Time Filter on Numerical Solutions to the Linearized Shallow-Water Wave Equations , 1983 .

[88]  H. Kanamori.  Mechanics of Earthquakes , 1994 .

[89]  Costas E. Synolakis,et al.  Tsunami inundation modeling for western Sumatra , 2006, Proceedings of the National Academy of Sciences.

[90]  D. Wells,et al.  New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement , 1994, Bulletin of the Seismological Society of America.

[91]  Yannis N. Krestenitis,et al.  Modelling the cohesive sediment transport in the marine environment: the case of Thermaikos Gulf , 2006 .

[92]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

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

[94]  David Ronald Henderson.  Local Source Tsunami Inundation Modelling for Poverty Bay, Gisborne , 2008 .

[95]  Fumihiko Imamura,et al.  Field survey of the East Java earthquake and tsunami of June 3, 1994 , 1995 .

[96]  Remko Scharroo,et al.  Satellite Altimeters Measure Tsunami , 2005 .

[97]  P. Lipman.  Natural hazards , 1993, Nature.

[98]  Fumihiko Imamura,et al.  Tsunami Catalog and Zones in Indonesia , 2000 .

[99]  Monte Kietpawpan,et al.  Method of calculating tsunami travel times in the Andaman Sea region , 2008, Natural hazards.

[100]  E. Asphaug,et al.  Asteroid Impact Tsunami: A Probabilistic Hazard Assessment , 2000 .

[101]  Tom De Groeve,et al.  Modelling Disaster Impact for the Global Disaster Alert and Coordination System , 2006 .

[102]  Michael Bader,et al.  Efficiency considerations in triangular adaptive mesh refinement , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[103]  Robert I. Tilling,et al.  This Dynamic Earth: The Story of Plate Tectonics , 1996 .

[104]  Frederick G. Shuman,et al.  NUMERICAL METHODS IN WEATHER PREDICTION: II. SMOOTHING AND FILTERING* , 1957 .

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

[106]  Fumihiko Imamura,et al.  Source models for the 1993 Hokkaido Nansei-Oki earthquake tsunami , 1995 .

[107]  Stefan Dech,et al.  "Last-Mile" preparation for a potential disaster - Interdisciplinary approach towards tsunami early warning and an evacuation information system for the coastal city of Padang, Indonesia , 2009 .

[108]  Isamu Aida,et al.  29. Numerical Experiments for the Tsunami Propagation : the 1964 Niigata Tsunami and the 1968 Tokachi-oki Tsunami , 1969 .

[109]  Roy A. Walters,et al.  Removing Tidal-Period Variations from Time-Series Data Using Low-Pass Digital Filters , 1982 .

[110]  Jiwei Tian,et al.  Estimates of global M2 internal tide energy fluxes using TOPEX/POSEIDON altimeter data , 2009 .

[111]  Y. Okada.  Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.