Parallel Implementation of Dispersive Tsunami Wave Modeling with a Nesting Algorithm for the 2011 Tohoku Tsunami

Because of improvements in offshore tsunami observation technology, dispersion phenomena during tsunami propagation have often been observed in recent tsunamis, for example the 2004 Indian Ocean and 2011 Tohoku tsunamis. The dispersive propagation of tsunamis can be simulated by use of the Boussinesq model, but the model demands many computational resources. However, rapid progress has been made in parallel computing technology. In this study, we investigated a parallelized approach for dispersive tsunami wave modeling. Our new parallel software solves the nonlinear Boussinesq dispersive equations in spherical coordinates. A variable nested algorithm was used to increase spatial resolution in the target region. The software can also be used to predict tsunami inundation on land. We used the dispersive tsunami model to simulate the 2011 Tohoku earthquake on the Supercomputer K. Good agreement was apparent between the dispersive wave model results and the tsunami waveforms observed offshore. The finest bathymetric grid interval was 2/9 arcsec (approx. 5 m) along longitude and latitude lines. Use of this grid simulated tsunami soliton fission near the Sendai coast. Incorporating the three-dimensional shape of buildings and structures led to improved modeling of tsunami inundation.

[1]  T. Sagiya,et al.  Slip distribution of the 2011 off the Pacific coast of Tohoku Earthquake inferred from geodetic data , 2011 .

[2]  Nobuo Shuto,et al.  The Nihonkai-Chubu Earthquake Tsunami on the North Akita Coast , 1985 .

[3]  T. Furumura,et al.  Tsunami waveform inversion including dispersive waves: the 2004 earthquake off Kii Peninsula, Japan , 2010 .

[4]  Dae-Hong Kim,et al.  Nested and multi-physics modeling of tsunami evolution from generation to inundation , 2011 .

[5]  Nobuhito Mori,et al.  Nationwide Post Event Survey and Analysis of the 2011 Tohoku Earthquake Tsunami , 2012 .

[6]  Wave Dispersion Effect in the Indian Ocean Tsunami , 2006 .

[7]  Marcelo H. Kobayashi,et al.  Shock-capturing Boussinesq-type model for nearshore wave processes , 2010 .

[8]  G. Pedersen,et al.  Modeling propagation and inundation of the 11 March 2011 Tohoku tsunami , 2012 .

[9]  Yoshinori Shigihara,et al.  Wave Dispersion Study in the Indian Ocean-Tsunami of December 26, 2004 , 2006 .

[10]  Narumi Takahashi,et al.  Tsunami Inundation Modeling of the 2011 Tohoku Earthquake Using Three-Dimensional Building Data for Sendai, Miyagi Prefecture, Japan , 2014 .

[11]  Fumihiko Imamura,et al.  Tsunami waveform inversion incorporating permanent seafloor deformation and its application to tsunami forecasting , 2012 .

[12]  William H. Press,et al.  Numerical recipes , 1990 .

[13]  Kazuhiro Kimura,et al.  Tsunami source area of the 2011 off the Pacific coast of Tohoku Earthquake determined from tsunami arrival times at offshore observation stations , 2011 .

[14]  T. Haraguchi,et al.  Field measurements and numerical modeling for the run-up heights and inundation distances of the 2011 Tohoku-oki tsunami at Sendai Plain, Japan , 2012, Earth, Planets and Space.

[15]  John D. Jakeman,et al.  Towards spatially distributed quantitative assessment of tsunami inundation models , 2010 .

[16]  Y. Ito,et al.  Tsunami source of the 2011 Tohoku‐Oki earthquake, Japan: Inversion analysis based on dispersive tsunami simulations , 2011 .

[17]  Isamu Aida,et al.  RELIABILITY OF A TSUNAMI SOURCE MODEL DERIVED FROM FAULT PARAMETERS , 1978 .

[18]  Kenji Satake,et al.  Time and Space Distribution of Coseismic Slip of the 2011 Tohoku Earthquake as Inferred from Tsunami Waveform Data , 2013 .

[19]  Stephan T. Grilli,et al.  Dispersive tsunami waves in the ocean: Model equations and sensitivity to dispersion and Coriolis effects , 2013 .

[20]  Yuji Yagi,et al.  Rupture process of the 2011 Tohoku‐oki earthquake and absolute elastic strain release , 2011 .

[21]  H. Suito,et al.  Interplate fault slip along the Japan Trench before the occurrence of the 2011 off the Pacific coast of Tohoku Earthquake as inferred from GPS data , 2011 .

[22]  Hiroshi Munekane,et al.  The 2011 off the Pacific coast of Tohoku Earthquake and its aftershocks observed by GEONET , 2011 .

[23]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

[24]  Stephan T. Grilli,et al.  Numerical Simulation of the 2011 Tohoku Tsunami Based on a New Transient FEM Co-seismic Source: Comparison to Far- and Near-Field Observations , 2013, Pure and Applied Geophysics.

[25]  T. Sakakiyama,et al.  A Study of Tsunami Wave Fission in an Undistorted Experiment , 2007 .

[26]  Chen Ji,et al.  Focal mechanism and slip history of the 2011 Mw 9.1 off the Pacific coast of Tohoku Earthquake, constrained with teleseismic body and surface waves , 2011 .

[27]  Hiroo Kanamori,et al.  A rupture model of the 2011 off the Pacific coast of Tohoku Earthquake , 2011 .

[28]  Per A. Madsen,et al.  On the solitary wave paradigm for tsunamis , 2008 .

[29]  D. Peregrine Equations for Water Waves and the Approximation behind Them , 1972 .

[30]  V. Titov,et al.  A nested-grid Boussinesq-type approach to modelling dispersive propagation and runup of landslide-generated tsunamis , 2011 .

[31]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[32]  T. Furumura,et al.  Significant tsunami observed at ocean-bottom pressure gauges during the 2011 off the Pacific coast of Tohoku Earthquake , 2011 .

[34]  P. Lynett,et al.  Parallel computation of a highly nonlinear Boussinesq equation model through domain decomposition , 2005 .