Supercycle at the Ecuadorian subduction zone revealed after the 2016 Pedernales earthquake

Large earthquakes are usually assumed to release all of the strain accumulated since the previous event, implying a reduced seismic hazard after them. However, long records of seismic history at several subduction zones suggest supercycle behaviour, where centuries-long accumulated strain is released through clustered large earthquakes, resulting in an extended period of enhanced seismic hazard. Here we combine historical seismology results, present-day geodesy data, and dense local observations of the recent Mw 7.8 2016 Pedernales earthquake to reconstruct the strain budget at the Ecuador subduction zone since the great 1906 earthquake. We show that the Pedernales earthquake involved the successive rupture of two patches on the plate interface that were locked prior to the earthquake and most probably overlaps the area already ruptured in 1942 by a similar earthquake. However, we find that coseismic slip in 2016 exceeds the deficit accumulated since 1942. The seismic moment of every large earthquake during the twentieth century further exceeds the moment deficit accumulated since 1906. These results, together with the seismic quiescence before 1906 highlighted by historical records and marine palaeoseismology, argue for an earthquake supercycle at the Ecuador–Colombia margin. This behaviour, which has led to an enhanced seismic hazard for 110 years, is possibly still going on and may apply to other subduction zones that recently experienced a great earthquake. Large earthquakes are often assumed to reset the seismic hazard of a region. Analysis of recent and historical seismicity in Ecuador suggests that this region may experience clusters of large earthquakes and extended periods of high seismic hazard.

[1]  Y. Font,et al.  Seismicity patterns along the Ecuadorian subduction zone: new constraints from earthquake location in a 3-D a priori velocity model , 2013 .

[2]  H. Kanamori,et al.  VARIABLE RUPTURE MODE OF THE SUBDUCTION ZONE ALONG THE ECUADOR-COLOMBIA COAST , 1982 .

[3]  L. Ruff,et al.  The rupture process of the Great 1979 Colombia Earthquake: Evidence for the asperity model , 1984 .

[4]  M. Vallée,et al.  Ten year recurrence time between two major earthquakes affecting the same fault segment , 2014 .

[5]  C. Werner,et al.  Radar interferogram filtering for geophysical applications , 1998 .

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

[7]  S. White,et al.  Recent crustal deformation and the earthquake cycle along the Ecuador–Colombia subduction zone , 2003 .

[8]  Jeffrey T. Freymueller,et al.  Wide plate margin deformation, southern Central America and northwestern South America, CASA GPS observations , 2002 .

[9]  Hai Cheng,et al.  Earthquake Supercycles Inferred from Sea-Level Changes Recorded in the Corals of West Sumatra , 2008, Science.

[10]  M. Vallée.  Stabilizing the Empirical Green Function Analysis: Development of the Projected Landweber Method , 2004 .

[11]  C. Mendoza,et al.  Seismicity associated with the great Colombia-Ecuador earthquakes of 1942, 1958, and 1979: Implications for barrier models of earthquake rupture , 1984 .

[12]  M. Métois,et al.  Three‐dimensional displacement field of the 2015 Mw8.3 Illapel earthquake (Chile) from across‐ and along‐track Sentinel‐1 TOPS interferometry , 2016 .

[13]  Yuji Yagi,et al.  Rupture process of the 2011 Tohoku‐oki earthquake and absolute elastic strain release , 2011 .

[14]  S. Owen,et al.  Nicoya earthquake rupture anticipated by geodetic measurement of the locked plate interface , 2014 .

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

[16]  Glenn P. Biasi,et al.  Wrightwood and the earthquake cycle: What a long recurrence record tells us about how faults work , 2004 .

[17]  Göran Ekström,et al.  The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes , 2012 .

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

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

[20]  Jaime Campos,et al.  The seismic coupling of subduction zones revisited , 2012 .

[21]  B. Delouis,et al.  Quantification of Tsunami Bathymetry Effect on Finite Fault Slip Inversion , 2015, Pure and Applied Geophysics.

[22]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[23]  M. Bevis,et al.  The Nazca -South America Euler vector and its rate of change , 2003 .

[24]  Peter Bormann,et al.  The ISC-GEM Global Instrumental Earthquake Catalogue (1900-2009): Introduction , 2015 .

[25]  Paul Lundgren,et al.  Joint Inversion of InSAR, GPS, Teleseismic, and Strong-Motion Data for the Spatial and Temporal Distribution of Earthquake Slip: Application to the 1999 İzmit Mainshock , 2002 .

[26]  F. Cotton,et al.  A new view for the geodynamics of Ecuador: Implication in seismogenic source definition and seismic hazard assessment , 2016 .

[27]  R. Yeats,et al.  Superquakes and Supercycles , 2013 .

[28]  W. D. Pennington.  Subduction of the Eastern Panama Basin and seismotectonics of northwestern South America , 1981 .

[29]  Thorne Lay,et al.  Rigidity variations with depth along interplate megathrust faults in subduction zones , 1999, Nature.

[30]  L. Bollinger,et al.  Rupture process of the Mw = 7.9 2015 Gorkha earthquake (Nepal): Insights into Himalayan megathrust segmentation , 2015 .

[31]  Raphaël Grandin,et al.  Interferometric Processing of SLC Sentinel-1 TOPS Data , 2015 .

[32]  M. Moreno,et al.  2010 Maule earthquake slip correlates with pre-seismic locking of Andean subduction zone , 2010, Nature.

[33]  S. Zaragosi,et al.  Earthquake-triggered deposits in the subduction trench of the north Ecuador/south Colombia margin and their implication for paleoseismology , 2017 .

[34]  R. Finkel,et al.  Earthquake supercycles in Central Italy, inferred from 36Cl exposure dating , 2011 .

[35]  P. Charvis,et al.  Distribution of discrete seismic asperities and aseismic slip along the Ecuadorian megathrust , 2014 .

[36]  E. Okal.  Use of the mantle magnitudeMm for the reassessment of the moment of historical earthquakes , 1992 .

[37]  J. Avouac,et al.  Earth and Planetary Science Letters , 2022 .

[38]  P. Heinrich,et al.  modeling of the February 1996 Peruvian Tsunami , 1998 .

[39]  V. E. Levin,et al.  Interseismic coupling and asperity distribution along the Kamchatka subduction zone , 2005 .

[40]  S. Beck,et al.  Historical 1942 Ecuador and 1942 Peru subduction earthquakes and earthquake cycles along Colombia-Ecuador and Peru subduction segments , 1996 .

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

[42]  J. Nocquet,et al.  Partitioning of oblique convergence in the Northern Andes subduction zone: Migration history and the present‐day boundary of the North Andean Sliver in Ecuador , 2016 .

[43]  J. C. Savage.  A dislocation model of strain accumulation and release at a subduction zone , 1983 .

[44]  Alistair Black,et al.  Introduction , 2004, Libr. Trends.

[45]  M. Vallée,et al.  SCARDEC: a new technique for the rapid determination of seismic moment magnitude, focal mechanism and source time functions for large earthquakes using body‐wave deconvolution , 2011 .

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

[47]  Yehuda Ben-Zion,et al.  Dynamic simulations of slip on a smooth fault in an elastic solid , 1997 .

[48]  Michael Bevis,et al.  Toward understanding tectonic control on the Mw 8.8 2010 Maule Chile earthquake , 2012 .

[49]  C. Beauval,et al.  Motion of continental slivers and creeping subduction in the Northern Andes , 2014 .

[50]  Hugo Yepes,et al.  Locations and magnitudes of historical earthquakes in the Sierra of Ecuador (1587–1996) , 2010 .

[51]  J. Nocquet,et al.  Continuous GPS Network Operating Throughout Ecuador , 2013 .

[52]  P. Stark.  Bounded-Variable Least-Squares: an Algorithm and Applications , 2008 .

[53]  M. Bouchon.  A simple method to calculate Green's functions for elastic layered media , 1981 .

[54]  C. Werner,et al.  Satellite radar interferometry: Two-dimensional phase unwrapping , 1988 .

[55]  Paul Wessel,et al.  Open radar interferometry software for mapping surface Deformation , 2011 .

[56]  Yuki Sawai,et al.  Predecessors of the giant 1960 Chile earthquake , 2005, Nature.

[57]  H. Yao,et al.  Coseismic radiation and stress drop during the 2015 Mw 8.3 Illapel, Chile megathrust earthquake , 2016 .

[58]  John Kelleher,et al.  Rupture zones of large South American earthquakes and some predictions , 1972 .

[59]  David J. Wald,et al.  Slab1.0: A three‐dimensional model of global subduction zone geometries , 2012 .

[60]  J. Avouac,et al.  Stress transfer and strain rate variations during the seismic cycle , 2004 .

[61]  K. Abe.  Size of great earthquakes of 1837–1974 inferred from tsunami data , 1979 .

[62]  Nadia Lapusta,et al.  Towards inferring earthquake patterns from geodetic observations of interseismic coupling , 2010 .

[63]  P. Rosen,et al.  Updated repeat orbit interferometry package released , 2004 .

[64]  P. Segall,et al.  A decadal‐scale deformation transient prior to the 2011 Mw 9.0 Tohoku‐oki earthquake , 2014 .

[65]  C. Beauval,et al.  An Earthquake Catalog for Seismic Hazard Assessment in Ecuador , 2013 .

[66]  Michel Bouchon,et al.  Teleseismic body wave radiation from a seismic source in a layered medium , 1976 .

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

[68]  Geodetic Observatory Tigo.  Toward understanding tectonic control on the Mw 8 .8 2010 Maule Chile earthquake , 2012 .