The 2013 Mw 7.7 Balochistan Earthquake: Seismic Potential of an Accretionary Wedge

Great earthquakes rarely occur within active accretionary prisms, despite the intense long‐term deformation associated with the formation of these geologic structures. This paucity of earthquakes is often attributed to partitioning of deformation across multiple structures as well as aseismic deformation within and at the base of the prism (Davis et al., 1983). We use teleseismic data and satellite optical and radar imaging of the 2013 M_w 7.7 earthquake that occurred on the southeastern edge of the Makran plate boundary zone to study this unexpected earthquake. We first compute a multiple point‐source solution from W‐phase waveforms to estimate fault geometry and rupture duration and timing. We then derive the distribution of subsurface fault slip from geodetic coseismic offsets. We sample for the slip posterior probability density function using a Bayesian approach, including a full description of the data covariance and accounting for errors in the elastic structure of the crust. The rupture nucleated on a subvertical segment, branching out of the Chaman fault system, and grew into a major earthquake along a 50° north‐dipping thrust fault with significant along‐strike curvature. Fault slip propagated at an average speed of 3.0  km/s for about 180 km and is concentrated in the top 10 km with no displacement on the underlying decollement. This earthquake does not exhibit significant slip deficit near the surface, nor is there significant segmentation of the rupture. We propose that complex interaction between the subduction accommodating the Arabia–Eurasia convergence to the south and the Ornach Nal fault plate boundary between India and Eurasia resulted in the significant strain gradient observed prior to this earthquake. Convergence in this region is accommodated both along the subduction megathrust and as internal deformation of the accretionary wedge.

[1]  Steven G. Wesnousky,et al.  Predicting the endpoints of earthquake ruptures , 2006, Nature.

[2]  K. Priestley,et al.  Surface waveform tomography of the Turkish–Iranian plateau , 2005 .

[3]  Marian Werner,et al.  Shuttle Radar Topography Mission (SRTM) Mission Overview , 2001 .

[4]  Luis Rivera,et al.  A note on the dynamic and static displacements from a point source in multilayered media , 2002 .

[5]  Y. Klinger continental strike-slip earthquake segmentation and thickness of the crust , 2010 .

[6]  Lin Ding,et al.  Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard , 2011 .

[7]  Rémi Michel,et al.  Horizontal coseismic deformation of the 1999 Chi-Chi earthquake measured from SPOT satellite images: Implications for the seismic cycle along the western foothills of central Taiwan , 2003 .

[8]  R. Yeats,et al.  Thrust and strike slip fault interaction along the Chaman transform zone, Pakistan , 1981, Geological Society, London, Special Publications.

[9]  N. Lapusta,et al.  Three‐dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip , 2009 .

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

[11]  Bertrand Meyer,et al.  Coseismic and early post-seismic slip associated with the 1999 Izmit earthquake (Turkey), from SAR interferometry and tectonic field observations , 2003 .

[12]  Yann Klinger,et al.  Coseismic deformation of the 2001 Mw = 7.8 Kokoxili earthquake in Tibet, measured by synthetic aperture radar interferometry , 2005 .

[13]  Mustafa Aktar,et al.  Extended Nucleation of the 1999 Mw 7.6 Izmit Earthquake , 2011, Science.

[14]  Jeffrey T. Freymueller,et al.  Coseismic slip distribution of the 2002 MW7.9 Denali fault earthquake, Alaska, determined from GPS measurements , 2003 .

[15]  Hiroo Kanamori,et al.  Uncertainty estimations for seismic source inversions , 2011 .

[16]  L. Rivera,et al.  Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations , 2002 .

[17]  Urs Wegmüller,et al.  Gamma SAR processor and interferometry software , 1997 .

[18]  R. Bilham,et al.  Earthquakes and Associated Deformation in Northern Baluchistan 1892-2001 , 2003 .

[19]  Fuk K. Li,et al.  Synthetic aperture radar interferometry , 2000, Proceedings of the IEEE.

[20]  M. Sambridge Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[21]  J. Tromp,et al.  Theoretical Global Seismology , 1998 .

[22]  James L. Beck,et al.  Bayesian inversion for finite fault earthquake source models I—theory and algorithm , 2013 .

[23]  Daniel M. Davis,et al.  GREAT THRUST EARTHQUAKES AND ASEISMIC SLIP ALONG THE PLATE BOUNDARY OF THE MAKRAN SUBDUCTION ZONE , 1992 .

[24]  H. Kanamori,et al.  The 2012 Sumatra great earthquake sequence , 2012 .

[25]  T. Farr,et al.  Shuttle radar topography mission produces a wealth of data , 2000 .

[26]  Hiroo Kanamori,et al.  Source inversion of W phase: speeding up seismic tsunami warning , 2008 .

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

[28]  M. Mokhtari,et al.  The diffuse transition between the Zagros continental collision and the Makran oceanic subduction (Iran): microearthquake seismicity and crustal structure , 2007 .

[29]  Jonathan M. Bull,et al.  The structure and fault activity of the Makran accretionary prism , 2012 .

[30]  D. Sandwell,et al.  Three-dimensional deformation caused by the Bam, Iran, earthquake and the origin of shallow slip deficit , 2005, Nature.

[31]  Demitris Paradissis,et al.  GPS constraints on continental deformation in the Africa‐Arabia‐Eurasia continental collision zone and implications for the dynamics of plate interactions , 2005 .

[32]  J. Beck,et al.  Accounting for prediction uncertainty when inferring subsurface fault slip , 2014 .

[33]  S. Lallemant,et al.  Impact of Sedimentation on Convergent Margin Tectonics: Example of the Makran Accretionary Prism (Pakistan) , 2007 .

[34]  J. Suppe,et al.  Mechanics of fold-and-thrust belts and accretionary wedges , 1983 .

[35]  Marie-Pierre Doin,et al.  Shallow creep on the Haiyuan Fault (Gansu, China) revealed by SAR Interferometry , 2012 .

[36]  Zhaohui Yang,et al.  Stick–slip advance of the Kohat Plateau in Pakistan , 2012 .

[37]  Sebastien Leprince,et al.  The 2005, Mw 7.6 Kashmir earthquake: Sub-pixel correlation of ASTER images and seismic waveforms analysis , 2006 .

[38]  M. Sambridge Geophysical inversion with a neighbourhood algorithm—II. Appraising the ensemble , 1999 .

[39]  Shuhab D. Khan,et al.  Slip-rates along the Chaman fault: Implication for transient strain accumulation and strain partitioning along the western Indian plate margin , 2013 .

[40]  Hiroo Kanamori,et al.  W phase source inversion for moderate to large earthquakes (1990–2010) , 2012 .

[41]  Rowena B. Lohman,et al.  Some thoughts on the use of InSAR data to constrain models of surface deformation: Noise structure and data downsampling , 2005 .

[42]  W. Szeliga,et al.  Interseismic strain accumulation along the western boundary of the Indian subcontinent , 2012 .