Geodetically constrained models of viscoelastic stress transfer and earthquake triggering along the North Anatolian fault

Over the past 80 years, 8 MW > 6.7 strike-slip earthquakes west of 40° longitude have ruptured the North Anatolian fault (NAF) from east to west. The series began with the 1939 Erzincan earthquake in eastern Turkey, and the most recent 1999 MW = 7.4 Izmit earthquake extended the pattern of ruptures into the Sea of Marmara in western Turkey. The mean time between seismic events in this westward progression is 8.5 ± 11 years (67% confidence interval), much greater than the timescale of seismic wave propagation (seconds to minutes). The delayed triggering of these earthquakes may be explained by the propagation of earthquake-generated diffusive viscoelastic fronts within the upper mantle that slowly increase the Coulomb failure stress change ( ΔCFS) at adjacent hypocenters. Here we develop three-dimensional stress transfer models with an elastic upper crust coupled to a viscoelastic Burgers rheology mantle. Both the Maxwell (ηM = 4 × 1018−1 × 1019 Pa s) and Kelvin (ηK = 1 × 1018−1 × 1019 Pa s) viscosities are constrained by studies of geodetic observations before and after the 1999 Izmit earthquake. We combine this geodetically constrained rheological model with the observed sequence of large earthquakes since 1939 to calculate the time evolution of ΔCFS changes along the North Anatolian fault due to viscoelastic stress transfer. Apparent threshold values of mean ΔCFS at which the earthquakes in the eight decade sequence occur are between ∼0.02 to ∼3.15 MPa and may exceed the magnitude of static ΔCFS values by as much as 177%. By 2023, we infer that the mean time-dependent stress change along the northern NAF strand in the Marmara Sea near Istanbul, which may have previously ruptured in 1766, may reach the mean apparent time-dependent stress thresholds of the previous NAF earthquakes.

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

[2]  Fred F. Pollitz,et al.  Stress Triggering of the 1999 Hector Mine Earthquake by Transient Deformation Following the 1992 Landers Earthquake , 2002 .

[3]  Massimo Cocco,et al.  Fault interaction by elastic stress changes: New clues from earthquake sequences , 2001 .

[4]  Yukitoshi Fukahata,et al.  Quasi-static internal deformation due to a dislocation source in a multilayered elastic/viscoelastic half-space and an equivalence theorem , 2006 .

[5]  F. Pollitz,et al.  The 1995 Kobe, Japan, earthquake: A long-delayed aftershock of the offshore 1944 Tonankai and 1946 Nankaido earthquakes , 1997, Bulletin of the Seismological Society of America.

[6]  A. N. Stokes,et al.  An Improved Method for Numerical Inversion of Laplace Transforms , 1982 .

[7]  A. To,et al.  Postseismic deformation and stress changes following the 1819 Rann of Kachchh, India earthquake: Was the 2001 Bhuj earthquake a triggered event? , 2004 .

[8]  Peizhen Zhang,et al.  Continuous deformation of the Tibetan Plateau from global positioning system data , 2004 .

[9]  D. Dreger,et al.  Static‐stress impact of the 1992 Landers earthquake sequence on nucleation and slip at the site of the 1999 M=7.1 Hector Mine earthquake, southern California , 2000 .

[10]  J. Elliott,et al.  Interseismic strain accumulation across the Manyi fault (Tibet) prior to the 1997 Mw 7.6 earthquake , 2011 .

[11]  M. Johnston,et al.  Remotely Triggered Seismicity on the United States West Coast following the Mw 7.9 Denali Fault Earthquake , 2004 .

[12]  D. Eberhart‐Phillips.  Magnitude, Slip-Partitioned Event The 2002 Denali Fault Earthquake, Alaska: A Large , 2007 .

[13]  R. Bürgmann,et al.  Dynamics of Izmit Earthquake Postseismic Deformation and Loading of the Duzce Earthquake Hypocenter , 2002 .

[14]  D. Marsan,et al.  Seismicity response to stress perturbations, analysed for a world-wide catalogue , 2003 .

[15]  D. Doser,et al.  Modeling Stress Changes Induced by Earthquakes in the Southern Marlborough Region, South Island, New Zealand , 2002 .

[16]  Frank Roth,et al.  Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone , 2006 .

[17]  L. Sykes,et al.  Stress evolution in southern California and triggering of moderate‐, small‐, and micro‐size earthquakes , 1997 .

[18]  M. Cocco,et al.  Normal fault interaction caused by coseismic and postseismic stress changes , 2001 .

[19]  A. Freed.  EARTHQUAKE TRIGGERING BY STATIC, DYNAMIC, AND POSTSEISMIC STRESS TRANSFER , 2005 .

[20]  Frank Roth,et al.  Modelling of stress patterns along the western part of the north anatolian fault zone , 1988 .

[21]  Semih Ergintav,et al.  Izmit earthquake postseismic deformation and dynamics of the North Anatolian Fault Zone , 2009 .

[22]  A. Barka,et al.  Slip distribution along the North Anatolian fault associated with the large earthquakes of the period 1939 to 1967 , 1996, Bulletin of the Seismological Society of America.

[23]  B. Hager,et al.  Postseismic relaxation across the Central Nevada Seismic Belt , 2003 .

[24]  Yaolin Shi,et al.  Studying the viscosity of lower crust of Qinghai-Tibet Plateau according to post-seismic deformation , 2009 .

[25]  J. Rice,et al.  Some basic stress diffusion solutions for fluid‐saturated elastic porous media with compressible constituents , 1976 .

[26]  P. Bodin,et al.  Triggering of earthquake aftershocks by dynamic stresses , 2000, Nature.

[27]  S. Carena,et al.  One hundred and fifty years of Coulomb stress history along the California‐Nevada border, USA , 2015 .

[28]  Amos Nur,et al.  Postseismic Viscoelastic Rebound , 1974, Science.

[29]  J. Brune,et al.  Seismicity in the Western Great Basin Apparently Triggered by the Landers, California, Earthquake, 28 June 1992 , 1994 .

[30]  P. Reasenberg,et al.  Response of Regional Seismicity to the Static Stress Change Produced by the Loma Prieta Earthquake , 1992, Science.

[31]  E. Hetland,et al.  Models of interseismic deformation with an analytic framework for the inclusion of general linear viscoelastic rheologies , 2006 .

[32]  Fred F. Pollitz,et al.  Lower crustal relaxation beneath the Tibetan Plateau and Qaidam Basin following the 2001 Kokoxili earthquake , 2011 .

[33]  P. R. DeVries,et al.  Earthquake cycle deformation in the Tibetan plateau with a weak mid‐crustal layer , 2012 .

[34]  E. Casarotti,et al.  Postseismic stress diffusion in Chile and South Peru , 2003 .

[35]  R. Armijo,et al.  Coulomb interactions and the 17 August 1999 Izmit, Turkey earthquake , 2001 .

[36]  P. Reasenberg,et al.  Earthquake triggering by seismic waves following the Landers and Hector Mine earthquakes , 2001, Nature.

[37]  Fred F. Pollitz,et al.  Fault interaction and stress triggering of twentieth century earthquakes in Mongolia , 2003 .

[38]  A. Barka,et al.  Heightened odds of large earthquakes near istanbul: An interaction-based probability calculation , 2000, Science.

[39]  Ryosuke Sato,et al.  Crustal Deformation due to Dislocation in a Multi-layered Medium , 1971 .

[40]  Mitsuhiro Matsu'ura,et al.  A dislocation model for the earthquake cycle at convergent plate boundaries , 1989 .

[41]  Yuehua Zeng,et al.  Viscoelastic stress‐triggering of the 1999 Hector Mine Earthquake by the 1992 Landers Earthquake , 2001 .

[42]  Michael Manga,et al.  A mechanism for sustained groundwater pressure changes induced by distant earthquakes , 2003 .

[43]  G. King,et al.  STATIC STRESS CHANGES AND THE TRIGGERING OF EARTHQUAKES , 1994 .

[44]  N. N. Ambraseys,et al.  The Seismic Activity of the Marmara Sea Region over the Last 2000 Years , 2002 .

[45]  Yukitoshi Fukahata,et al.  General expressions for internal deformation fields due to a dislocation source in a multilayered elastic half-space , 2005 .

[46]  P. Bodin,et al.  Earthquake nucleation by transient deformations caused by the M = 7.9 Denali, Alaska, earthquake , 2004, Nature.

[47]  P. Reasenberg,et al.  Stress sensitivity of fault seismicity: A comparison between limited‐offset oblique and major strike‐slip faults , 1999 .

[48]  W. Thatcher,et al.  Reconciling viscoelastic models of postseismic and interseismic deformation: Effects of viscous shear zones and finite length ruptures , 2015 .

[49]  James Jackson,et al.  Seismicity of the Sea of Marmara (Turkey) since 1500 , 2000 .

[50]  Quasi-Static Strain and Stress Fields due to a Moment Tensor in Elastic–Viscoelastic Layered Half-Space , 2014, Pure and Applied Geophysics.

[51]  Timothy E. Dawson,et al.  The 2002 Denali Fault Earthquake, Alaska: A Large Magnitude, Slip-Partitioned Event , 2003, Science.

[52]  W. Ellsworth,et al.  Seismicity Remotely Triggered by the Magnitude 7.3 Landers, California, Earthquake , 1993, Science.

[53]  Aaron A. Velasco,et al.  Remote triggered seismicity caused by the 2011, M9.0 Tohoku‐Oki, Japan earthquake , 2012 .

[54]  P. Bodin,et al.  Aftershock triggering by complete Coulomb stress changes , 2002 .

[55]  Toshinori Sato,et al.  A kinematic model for deformation of the lithosphere at subduction zones , 1988 .

[56]  B. Hager,et al.  Viscoelastic deformation from North Anatolian Fault Zone earthquakes and the eastern Mediterranean GPS velocity field , 2002 .

[57]  A. Freed,et al.  Accelerated stress buildup on the southern San Andreas fault and surrounding regions caused by Mojave Desert earthquakes , 2002 .

[58]  F. Pollitz,et al.  Viscosity of oceanic asthenosphere inferred from remote triggering of earthquakes , 1998, Science.

[59]  B. Meade,et al.  Inference of Multiple Earthquake‐Cycle Relaxation Timescales from Irregular Geodetic Sampling of Interseismic Deformation , 2013 .

[60]  Tim J. Wright,et al.  Post-seismic motion following the 1997 Manyi (Tibet) earthquake: InSAR observations and modelling , 2007 .

[61]  B. Hager,et al.  The effects of rheological layering on post-seismic deformation , 2006 .

[62]  L. Sykes,et al.  Evolution of the stress field in southern California and triggering of moderate-size earthquakes: A 200-year perspective , 1997 .

[63]  E. Brodsky,et al.  New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera , 2005 .

[64]  Paul Segall,et al.  Post-earthquake ground movements correlated to pore-pressure transients , 2003, Nature.

[65]  Jian Lin,et al.  Delayed triggering of the 1999 Hector Mine earthquake by viscoelastic stress transfer , 2001, Nature.

[66]  G. King,et al.  Stress coupling between earthquakes in northwest Turkey and the north Aegean Sea , 1996 .

[67]  B. Atkinson.  Subcritical crack growth in geological materials , 1984 .

[68]  A. Barka,et al.  The 17 August 1999 Izmit Earthquake , 1999, Science.

[69]  T. Wright,et al.  Weak ductile shear zone beneath a major strike‐slip fault: Inferences from earthquake cycle model constrained by geodetic observations of the western North Anatolian Fault Zone , 2014 .

[70]  James H. Dieterich,et al.  Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering , 1997 .

[71]  Ryosuke Sato,et al.  STATIC DEFORMATIONS DUE TO THE FAULT SPREADING OVER SEVERAL LAYERS IN A MULTI-LAYERED MEDIUM , 1973 .

[72]  E. Hauksson,et al.  The static stress change triggering model: Constraints from two southern California aftershock sequences , 1998 .

[73]  J. C. Savage,et al.  Asthenosphere readjustment and the earthquake cycle , 1978 .

[74]  Haluk Ozener,et al.  Seven years of postseismic deformation following the 1999, M = 7.4 and M = 7.2, Izmit-Düzce, Turkey earthquake sequence , 2009 .

[75]  Tom Parsons,et al.  Recalculated probability of M >= 7 earthquakes beneath the Sea of Marmara, Turkey , 2004 .

[76]  M. Cocco,et al.  Introduction to special section: Stress transfer, earthquake triggering, and time‐dependent seismic hazard , 2005 .

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

[78]  R. Stein,et al.  Earthquake conversations. , 2003, Scientific American.

[79]  J. Dieterich.  Earthquake nucleation on faults with rate-and state-dependent strength , 1992 .

[80]  R. Simpson,et al.  Changes in static stress on southern California faults after the 1992 Landers earthquake , 1992, Nature.

[81]  G. Tselentis,et al.  Stress transfer and nonlinear stress accumulation at the North Anatolian fault, Turkey , 1990 .

[82]  J. Dieterich.  A constitutive law for rate of earthquake production and its application to earthquake clustering , 1994 .

[83]  F. Pollitz,et al.  The 11 April 2012 east Indian Ocean earthquake triggered large aftershocks worldwide , 2012, Nature.

[84]  Evelyn Roeloffs,et al.  Fault stability changes induced beneath a reservoir with cyclic variations in water level , 1988 .

[85]  Semih Ergintav,et al.  Time-Dependent Distributed Afterslip on and Deep below the İzmit Earthquake Rupture , 2002 .

[86]  Haluk Ozener,et al.  Estimates of Seismic Potential in the Marmara Sea Region from Block Models of Secular Deformation Constrained by Global Positioning System Measurements , 2002 .

[87]  Ruth A. Harris,et al.  Introduction to Special Section: Stress Triggers, Stress Shadows, and Implications for Seismic Hazard , 1998 .

[88]  Ryosuke Sato,et al.  STATIC DEFORMATION DUE TO THE FAULT SPREADING OVER SEVERAL LAYERS IN A MULTI-LAYERED MEDIUM PART II: STRAIN AND TILT , 1975 .

[89]  F. Pollitz,et al.  Consequences of stress changes following the 1891 Nobi earthquake, Japan , 1995 .