Validation of linearity assumptions for using tsunami waveforms in joint inversion of kinematic rupture models: Application to the 2010 Mentawai Mw 7.8 tsunami earthquake

Tsunami observations have particular importance for resolving shallow offshore slip in finite‐fault rupture model inversions for large subduction zone earthquakes. However, validations of amplitude linearity and choice of subfault discretization of tsunami Green's functions are essential when inverting tsunami waveforms. We explore such validations using four tsunami recordings of the 25 October 2010 Mentawai Mw 7.8 tsunami earthquake, jointly inverted with teleseismic body waves and 1 Hz GPS (high‐rate GPS) observations. The tsunami observations include near‐field and far‐field deep water recordings, as well as coastal and island tide gauge recordings. A nonlinear, dispersive modeling code, NEOWAVE, is used to construct tsunami Green's functions from seafloor excitation for the linear inversions, along with performing full‐scale calculations of the tsunami for the inverted models. We explore linearity and finiteness effects with respect to slip magnitude, variable rake determination, and subfault dimensions. The linearity assumption is generally robust for the deep water recordings, and wave dispersion from seafloor excitation is important for accurate description of near‐field Green's functions. Breakdown of linearity produces substantial misfits for short‐wavelength signals in tide gauge recordings with large wave heights. Including the tsunami observations in joint inversions provides improved resolution of near‐trench slip compared with inversions of only seismic and geodetic data. Two rupture models, with fine‐grid (15 km) and coarse‐grid (30 km) spacing, are inverted for the Mentawai event. Stronger regularization is required for the fine model representation. Both models indicate a shallow concentration of large slip near the trench with peak slip of ~15 m. Fully nonlinear forward modeling of tsunami waveforms confirms the validity of these two models for matching the tsunami recordings along with the other data.

[1]  T. Lay.  The surge of great earthquakes from 2004 to 2014 , 2015 .

[2]  L. Ye,et al.  Tsunami surges around the Hawaiian Islands from the 1 April 2014 North Chile Mw 8.1 earthquake , 2014 .

[3]  T. Lay.  A GLOBAL SURGE OF GREAT EARTHQUAKES FROM 2004-2014 AND IMPLICATIONS FOR CASCADIA , 2014 .

[4]  T. Lay,et al.  Localized fault slip to the trench in the 2010 Maule, Chile Mw = 8.8 earthquake from joint inversion of high‐rate GPS, teleseismic body waves, InSAR, campaign GPS, and tsunami observations , 2014 .

[5]  Anthony Sladen,et al.  A detailed source model for the Mw9.0 Tohoku‐Oki earthquake reconciling geodesy, seismology, and tsunami records , 2014 .

[6]  F. Romano,et al.  Structural control on the Tohoku earthquake rupture process investigated by 3D FEM, tsunami and geodetic data , 2014, Scientific Reports.

[7]  K. Sieh,et al.  Rupture process of the 2010 Mw 7.8 Mentawai tsunami earthquake from joint inversion of near‐field hr‐GPS and teleseismic body wave recordings constrained by tsunami observations , 2014 .

[8]  C. An,et al.  Tsunami source and its validation of the 2014 Iquique, Chile, earthquake , 2014 .

[9]  Chao An,et al.  The 1 April 2014 Iquique, Chile, Mw 8.1 earthquake rupture sequence , 2014 .

[10]  H. Kanamori.  The Diversity of Large Earthquakes and Its Implications for Hazard Mitigation , 2014 .

[11]  K. Satake,et al.  Traveltime delay and initial phase reversal of distant tsunamis coupled with the self‐gravitating elastic Earth , 2014 .

[12]  S. Allgeyer,et al.  Numerical tsunami simulation including elastic loading and seawater density stratification , 2014 .

[13]  J. Avouac,et al.  The 2013, Mw 7.7 Balochistan earthquake, energetic strike-slip reactivation of a thrust fault , 2014 .

[14]  H. Kanamori,et al.  The February 6, 2013 Mw 8.0 Santa Cruz Islands earthquake and tsunami , 2013 .

[15]  Yehuda Bock,et al.  Near‐field tsunami models with rapid earthquake source inversions from land‐ and ocean‐based observations: The potential for forecast and warning , 2013 .

[16]  S. Owen,et al.  The 5 September 2012 Nicoya, Costa Rica Mw 7.6 earthquake rupture process from joint inversion of high‐rate GPS, strong‐motion, and teleseismic P wave data and its relationship to adjacent plate boundary interface properties , 2013 .

[17]  H. Kanamori,et al.  The October 28, 2012 Mw 7.8 Haida Gwaii underthrusting earthquake and tsunami: Slip partitioning along the Queen Charlotte Fault transpressional plate boundary , 2013 .

[18]  K. Cheung,et al.  Dispersion and nonlinearity of multi-layer non-hydrostatic free-surface flow , 2013, Journal of Fluid Mechanics.

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

[20]  T. Lay,et al.  Source Rupture Models for the Mw 9.0 2011 Tohoku Earthquake from Joint Inversions of High‐Rate Geodetic and Seismic Data , 2013 .

[21]  H. Kanamori,et al.  Estimating the effect of Earth elasticity and variable water density on tsunami speeds , 2013 .

[22]  Yefei Bai,et al.  Depth‐integrated free‐surface flow with parameterized non‐hydrostatic pressure , 2013 .

[23]  H. Latief,et al.  Tsunami Source of the 2010 Mentawai, Indonesia Earthquake Inferred from Tsunami Field Survey and Waveform Modeling , 2013, Pure and Applied Geophysics.

[24]  Pedro Elosegui,et al.  The 2010 Mw 7.8 Mentawai earthquake: Very shallow source of a rare tsunami earthquake determined from tsunami field survey and near‐field GPS data , 2012 .

[25]  N. D’Agostino,et al.  Clues from joint inversion of tsunami and geodetic data of the 2011 Tohoku-oki earthquake , 2012, Scientific Reports.

[26]  H. Kanamori,et al.  The 2011 Mw 9.0 off the Pacific coast of Tohoku Earthquake: Comparison of deep-water tsunami signals with finite-fault rupture model predictions , 2011 .

[27]  S. Singh,et al.  Seismic images of the megathrust rupture during the 25th October 2010 Pagai earthquake, SW Sumatra: Frontal rupture and large tsunami , 2011 .

[28]  E. Engdahl,et al.  The 25 October 2010 Sumatra tsunami earthquake: Slip in a slow patch , 2011 .

[29]  T. Kanazawa,et al.  Joint inversion of strong motion, teleseismic, geodetic, and tsunami datasets for the rupture process of the 2011 Tohoku earthquake , 2011 .

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

[31]  Fred F. Pollitz,et al.  Geodetic slip model of the 2011 M9.0 Tohoku earthquake , 2011 .

[32]  K. Cheung,et al.  Modeling near‐field tsunami observations to improve finite‐fault slip models for the 11 March 2011 Tohoku earthquake , 2011 .

[33]  Thorne Lay,et al.  Inversion of high‐rate (1 sps) GPS data for rupture process of the 11 March 2011 Tohoku earthquake (Mw 9.1) , 2011 .

[34]  Gavin P. Hayes,et al.  The 25 October 2010 Mentawai tsunami earthquake, from real‐time discriminants, finite‐fault rupture, and tsunami excitation , 2011 .

[35]  John McCloskey,et al.  Limited overlap between the seismic gap and coseismic slip of the great 2010 Chile earthquake , 2011 .

[36]  H. Kanamori,et al.  The 25 October 2010 Mentawai tsunami earthquake (Mw 7.8) and the tsunami hazard presented by shallow megathrust ruptures , 2011 .

[37]  Slip distribution of the 2003 Tokachi‐oki Mw 8.1 earthquake from joint inversion of tsunami waveforms and geodetic data , 2010 .

[38]  David T. Sandwell,et al.  Coseismic slip model of the 2008 Wenchuan earthquake derived from joint inversion of interferometric synthetic aperture radar, GPS, and field data , 2010 .

[39]  V. Popov.  Earthquakes and Friction , 2010 .

[40]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model with grid nesting for tsunami generation, propagation, and run‐up , 2010 .

[41]  H. Kanamori,et al.  Effects of Kinematic Constraints on Teleseismic Finite-Source Rupture Inversions: Great Peruvian Earthquakes of 23 June 2001 and 15 August 2007 , 2009 .

[42]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model for wave breaking and run‐up , 2009 .

[43]  Peizhen Zhang,et al.  Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake , 2009 .

[44]  P. Rosen,et al.  Interferometric Synthetic Aperture Radar Geodesy , 2007 .

[45]  Rongjiang Wang,et al.  Erratum to: "Computation of deformation induced by earthquakes in a multi-layered elastic crust - FORTRAN programs EDGRN/EDCMP": [Computers & Geosciences, 29(2) (2003) 195-207] , 2006, Comput. Geosci..

[46]  Kristine M. Larson,et al.  Modeling the rupture process of the 2003 September 25 Tokachi‐Oki (Hokkaido) earthquake using 1‐Hz GPS data , 2004 .

[47]  WangRongjiang,et al.  Computation of deformation induced by earthquakes in a multi-layered elastic crust , 2003 .

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

[49]  P. Rosen,et al.  SYNTHETIC APERTURE RADAR INTERFEROMETRY TO MEASURE EARTH'S SURFACE TOPOGRAPHY AND ITS DEFORMATION , 2000 .

[50]  Rongjiang Wang,et al.  A simple orthonormalization method for stable and efficient computation of Green's functions , 1999, Bulletin of the Seismological Society of America.

[51]  C. Scholz.  Earthquakes and friction laws , 1998, Nature.

[52]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[53]  Kenji Satake,et al.  Tsunami generation by horizontal displacement of ocean bottom , 1996 .

[54]  Kenji Satake,et al.  The 1964 Prince William Sound earthquake: Joint inversion of tsunami and geodetic data , 1996 .

[55]  P. Segall,et al.  The co-seismic slip distribution of the Landers earthquake , 1994, Bulletin of the Seismological Society of America.

[56]  R. Goldstein,et al.  Satellite Radar Interferometry for Monitoring Ice Sheet Motion: Application to an Antarctic Ice Stream , 1993, Science.

[57]  Masayuki Kikuchi,et al.  Source complexity of the 1988 Armenian Earthquake: Evidence for a slow after‐slip event , 1993 .

[58]  K. Feigl,et al.  The displacement field of the Landers earthquake mapped by radar interferometry , 1993, Nature.

[59]  K. Satake.  Depth distribution of coseismic slip along the Nankai Trough, Japan, from joint inversion of geodetic and tsunami data , 1993 .

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

[61]  Masayuki Kikuchi,et al.  Inversion of complex body waves—III , 1991, Bulletin of the Seismological Society of America.

[62]  K. Satake.  Inversion of tsunami waveforms for the estimation of heterogeneous fault motion of large submarine earthquakes: The 1968 Tokachi‐oki and 1983 Japan Sea earthquakes , 1989 .

[63]  Kenji Satake,et al.  Inversion of tsunami waveforms for the estimation of a fault heterogeneity: Method and numerical experiments. , 1987 .

[64]  Thomas H. Heaton,et al.  Inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 Imperial Valley, California, earthquake , 1983 .

[65]  Masayuki Kikuchi,et al.  Inversion of complex body waves , 1982 .

[66]  Kinjiro Kajiura,et al.  The Leading Wave of a Tsunami , 1963 .