日本地球惑星科学連合2015年大会

講演情報

インターナショナルセッション(口頭発表)

セッション記号 S (固体地球科学) » S-SS 地震学

[S-SS02] Frontier studies on subduction zone megathrust earthquakes and tsunamis

2015年5月25日(月) 16:15 〜 18:00 国際会議室 (2F)

コンビーナ:*金川 久一(千葉大学大学院理学研究科)、Demian Saffer(Dept. of Geosciences, The Pennsylvania State University, USA)、Michael Strasser(Geological Institute, Seiss Federal Insitute of Technology ETH Zurich)、山田 泰広(京都大学工学研究科都市社会工学専攻)、小平 秀一(海洋研究開発機構 地球内部ダイナミクス領域)、日野 亮太(東北大学災害科学国際研究所)、氏家 恒太郎(筑波大学生命環境系)、伊藤 喜宏(京都大学防災研究所)、座長:氏家 恒太郎(筑波大学生命環境系)、金川 久一(千葉大学大学院理学研究科)

16:50 〜 17:05

[SSS02-22] Rupture process of the 2014 Iquique, Chile earthquake estimated from tsunami waveforms, teleseismic, and GPS data

*Aditya GUSMAN1Satoko MUROTANI1Kenji SATAKE1Mohammad HEIDARZADEH1Shingo WATADA1Endra GUNAWAN2Bernd SCHURR3 (1.Earthquake Research Institute, the University of Tokyo、2.Graduate School of Environmental Studies, Nagoya University、3.GFZ Helmholtz Centre Potsdam, German Research Centre for Geosciences)

キーワード:Rupture process, Tsunami waveforms, GPS data, Teleseismic body waves, Tsunami dispersion, Joint inversion

We applied a new method to compute tsunami Green's functions for slip inversion of the 1 April 2014 Iquique earthquake using both near-field and far-field tsunami waveforms. Inclusion of the effects of the elastic loading of seafloor, compressibility of seawater, and the geopotential variation in the computed Green's functions reproduced the tsunami travel-time delay relative to long-wave simulation, and allowed us to use far-field records in tsunami waveform inversion. Multiple time window inversion (Satake et al., BSSA, 2013) was applied to tsunami waveforms iteratively until the result resembles the stable moment-rate function from teleseismic inversion. We also used GPS data to perform a joint inversion of tsunami waveforms and co-seismic crustal deformation (Gusman et al., EPSL, 2012). According to results, the major slip region with a size of 100 km × 40 km is located down-dip the epicenter at depth ~28 km, regardless of assumed rupture velocities. The total seismic moment from the slip distribution estimated by the joint inversion is 1.24 × 1021 Nm (Mw 8.0) (Gusman et al., GRL, 2015). This seismic moment is slightly smaller than 1.88 × 1021 Nm (Mw 8.1) from a teleseismic waveform inversion.
The tsunami arrival time and polarity reversal observed at far-field DART stations can be accurately reproduced by solving shallow water equations and applying the phase velocity correction to the simulated waveforms (Watada et al., JGR, 2014). The slip distribution of the 2014 Iquique earthquake from our joint inversion method can accurately explain the tsunami waveform in the near-field as well as in the far-field. We propose the tsunami phase velocity correction to be included as a standard procedure in inversion methods when using far-field tsunami waveforms.
The teleseismic inversion with different rupture velocities (1.5, 2.0, and 2.5 km/s) yielded similar moment rate functions which all peaked at ~35 s, but their spatial slip distributions are different. On the contrary, the joint inversion gives a stable spatial slip distribution for different rupture velocities. Among the slip distributions from the teleseismic inversions with the three different rupture velocities, the one for 1.5 km/s is most similar to the slip distribution from the joint inversion of tsunami waveforms and GPS data in terms of large slip area. Thus, the velocity of 1.5 km/s may better represent the rupture process of the 2014 Iquique earthquake (Gusman et al., GRL, 2015).