Geodetic imaging of coseismic slip and postseismic afterslip: Sparsity promoting methods applied to the great Tohoku earthquake

[1] Geodetic observations of surface displacements during and following earthquakes such as the March 11, 2011 great Tohoku earthquake can be used to constrain the spatial extent of coseismic slip and postseismic afterslip, and characterize the spectrum of earthquake cycle behaviors. Slip models are often regularized by assuming that slip on the fault varies smoothly in space, which may result in the artificial smearing of fault slip beyond physical boundaries. Alternatively, it may be desirable to estimate a slip distribution that is spatially compact and varies sharply. Here we show that sparsity promoting state vector regularization methods can be used to recover slip distributions with sharp boundaries, representing an alternative end-member result to very smooth slip distributions. Using onshore GPS observations at 298 stations during and in the ∼2 weeks following the Tohoku earthquake, we estimate a band of coseismic slip between 30 and 50 km depth extending 500 km along strike with a maximum slip of 64 m, corresponding to a minimum magnitude estimate of MW = 8.8. Our estimate of afterslip is located almost exclusively down-dip of the coseismic rupture, with a transition between 40 and 50 km depth and an equivalent moment magnitude MW = 8.2. This depth may be interpreted as coincident with the transition from velocity strengthening to velocity weakening frictional behavior, consistent with the upper limit of cold subduction estimates of the thermal structure of the Japan trench.

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

[2]  Sarah E. Minson,et al.  The 2011 Magnitude 9.0 Tohoku-Oki Earthquake: Mosaicking the Megathrust from Seconds to Centuries , 2011, Science.

[3]  R. Bilham,et al.  Postseismic deformation of the Andaman Islands following the 26 December, 2004 Great Sumatra–Andaman earthquake , 2007 .

[4]  Changrong He,et al.  Frictional sliding of gabbro gouge under hydrothermal conditions , 2007 .

[5]  Masataka Ando,et al.  Evidence of large scale repeating slip during the 2011 Tohoku‐Oki earthquake , 2011 .

[6]  James D. Byerlee,et al.  Frictional slip of granite at hydrothermal conditions , 1995 .

[7]  Motoyuki Kido,et al.  Trench‐normal variation in observed seafloor displacements associated with the 2011 Tohoku‐Oki earthquake , 2011 .

[8]  Kristine M. Larson,et al.  Coseismic and early postseismic slip for the 2003 Tokachi‐oki earthquake sequence inferred from GPS data , 2008 .

[9]  Tomokazu Kobayashi,et al.  Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake , 2011, Nature.

[10]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..

[11]  D. Donoho For most large underdetermined systems of equations, the minimal 𝓁1‐norm near‐solution approximates the sparsest near‐solution , 2006 .

[12]  P. Shearer,et al.  Compressive sensing of the Tohoku‐Oki Mw 9.0 earthquake: Frequency‐dependent rupture modes , 2011 .

[13]  Masaya Matsuura,et al.  Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip , 1992 .

[14]  P. Segall,et al.  Seismic and aseismic fault slip before and during the 2011 off the Pacific coast of Tohoku Earthquake , 2011 .

[15]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[16]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[17]  Chris Marone,et al.  On the mechanics of earthquake afterslip , 1991 .

[18]  Kenji Satake,et al.  Tsunami source of the 2011 off the Pacific coast of Tohoku Earthquake , 2011 .

[19]  David L. Donoho,et al.  Observed universality of phase transitions in high-dimensional geometry, with implications for modern data analysis and signal processing , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[20]  J. Claerbout,et al.  Robust Modeling With Erratic Data , 1973 .

[21]  John R. Rice,et al.  Crustal Earthquake Instability in Relation to the Depth Variation of Frictional Slip Properties , 1986 .

[22]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[23]  Nobuhito Mori,et al.  Survey of 2011 Tohoku earthquake tsunami inundation and run‐up , 2011 .

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

[25]  P. Segall,et al.  Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas Fault , 1987 .

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

[27]  Akira Asada,et al.  Displacement Above the Hypocenter of the 2011 Tohoku-Oki Earthquake , 2011, Science.

[28]  Kosuke Heki Space geodetic observation of deep basal subduction erosion in northeastern Japan , 2004 .

[29]  B. Meade,et al.  Geodetic imaging of plate motions, slip rates, and partitioning of deformation in Japan , 2010 .

[30]  K. Obara,et al.  Three‐dimensional crustal S wave velocity structure in Japan using microseismic data recorded by Hi‐net tiltmeters , 2008 .

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

[32]  R. Madariaga,et al.  The 2010 Mw 8.8 Maule Megathrust Earthquake of Central Chile, Monitored by GPS , 2011, Science.

[33]  Kurt L. Feigl,et al.  Crustal deformation and fault slip during the seismic cycle in the North Chile subduction zone, from GPS and InSAR observations , 2004 .

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

[35]  F. Santosa,et al.  Linear inversion of ban limit reflection seismograms , 1986 .

[36]  Gregory C. Beroza,et al.  Shallow Dynamic Overshoot and Energetic Deep Rupture in the 2011 Mw 9.0 Tohoku-Oki Earthquake , 2011, Science.

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

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

[39]  D. Pollard,et al.  Inverting for slip on three-dimensional fault surfaces using angular dislocations , 2005 .

[40]  Paul Segall,et al.  Rapid afterslip following the 1999 Chi‐Chi, Taiwan Earthquake , 2002 .

[41]  Kelin Wang,et al.  Seismic consequences of warm versus cool subduction metamorphism: examples from southwest and northeast japan , 1999, Science.

[42]  Brendan J. Meade,et al.  Spatial correlation of interseismic coupling and coseismic rupture extent of the 2011 MW = 9.0 Tohoku‐oki earthquake , 2011 .

[43]  Yuji Yagi,et al.  A unified source model for the 2011 Tohoku earthquake , 2011 .