The Unlocking Process Leading to the 2016 Central Italy Seismic Sequence

Approximately 23,000 well‐located earthquakes from 2009 to 2016 are used as templates to recover seismic activity preceding the 2016 Central Italy seismic sequence. The resulting spatiotemporal pattern is analyzed by additional ∼91,000 newly detected events. In the 8 years before the sequence onset, seismicity (ML ≤ 3.7) develops at the hangingwall of the 2016 normal faults and along a sub‐horizontal shear zone, bounding the active extensional system at depth. This activity, mainly organized in foreshock‐mainshock and swarm‐like clusters, migrates toward the nucleation area of the first Mw 6.0 mainshock of the sequence (24th of August in Amatrice). We propose an unlocking process based on variable temporal clustering of the seismicity, including repeaters, identifying fault portions with different degrees of coupling. Such a progressive localization of the seismic activity at the fault edges induces a weakening of the locked patch of the Amatrice mainshock.

[1]  M. Sugan,et al.  A Python Code for Detecting True Repeating Earthquakes from Self-Similar Waveforms (FINDRES) , 2022, Seismological Research Letters.

[2]  W. Frank,et al.  Tracking the Spatio‐Temporal Evolution of Foreshocks Preceding the Mw 6.1 2009 L’Aquila Earthquake , 2022, Journal of Geophysical Research: Solid Earth.

[3]  Z. Ross,et al.  A unified perspective of seismicity and fault coupling along the San Andreas Fault , 2022, Science advances.

[4]  L. Malagnini,et al.  Crustal permeability changes inferred from seismic attenuation: Impacts on multi-mainshock sequences , 2021, Frontiers in Earth Science.

[5]  T. Kraft,et al.  Possible Precursory Slow‐Slip to Two ML∼3 Mainevents of the Diemtigen Microearthquake Sequence, Switzerland , 2021, Geophysical Research Letters.

[6]  F. Waldhauser,et al.  Fault Planes, Fault Zone Structure and Detachment Fragmentation Resolved With High‐Precision Aftershock Locations of the 2016–2017 Central Italy Sequence , 2021, Geophysical Research Letters.

[7]  Yen Joe Tan,et al.  Machine-Learning-Based High-Resolution Earthquake Catalog Reveals How Complex Fault Structures Were Activated during the 2016–2017 Central Italy Sequence , 2021 .

[8]  A. Gualandi,et al.  Post‐Seismic Deformation Related to the 2016 Central Italy Seismic Sequence From GPS Displacement Time‐Series , 2021, Journal of geophysical research. Solid earth.

[9]  L. Chiaraluce,et al.  The influence of subsurface geology on the distribution of earthquakes during the 2016‐–2017 Central Italy seismic sequence , 2021 .

[10]  Y. Ben‐Zion,et al.  The generation of large earthquakes , 2020, Nature Reviews Earth & Environment.

[11]  I. Zaliapin,et al.  Localization and coalescence of seismicity before large earthquakes , 2020 .

[12]  A. Aoudia,et al.  Seismicity Rate Changes and Geodetic Transients in Central Apennines , 2020, Geophysical Research Letters.

[13]  L. Chiaraluce,et al.  Intermittent Slip Along the Alto Tiberina Low‐Angle Normal Fault in Central Italy , 2020, Geophysical Research Letters.

[14]  M. Segou,et al.  Seismic velocity precursors to the 2016 Mw 6.5 Norcia (Italy) earthquake , 2020, Geology.

[15]  F. Waldhauser,et al.  Fine‐Scale Structure of the 2016–2017 Central Italy Seismic Sequence From Data Recorded at the Italian National Network , 2020, Journal of Geophysical Research: Solid Earth.

[16]  I. Zaliapin,et al.  Earthquake Declustering Using the Nearest‐Neighbor Approach in Space‐Time‐Magnitude Domain , 2020, Journal of Geophysical Research: Solid Earth.

[17]  E. Dunham,et al.  Fault valving and pore pressure evolution in simulations of earthquake sequences and aseismic slip , 2020, Nature Communications.

[18]  D. Latorre,et al.  An Empirical Formula to Classify the Quality of Earthquake Locations , 2019, Bulletin of the Seismological Society of America.

[19]  N. Uchida.  Detection of repeating earthquakes and their application in characterizing slow fault slip , 2019, Progress in Earth and Planetary Science.

[20]  P. Meakin,et al.  Volumetric and shear processes in crystalline rock approaching faulting , 2019, Proceedings of the National Academy of Sciences.

[21]  Egill Hauksson,et al.  Searching for hidden earthquakes in Southern California , 2019, Science.

[22]  G. Amati,et al.  Seismic Evidence of an Early Afterslip During the 2012 Sequence in Emilia (Italy) , 2019, Geophysical Research Letters.

[23]  C. Ji,et al.  Earthquake nucleation and fault slip complexity in the lower crust of central Alaska , 2018, Nature Geoscience.

[24]  W. Ellsworth,et al.  Nucleation of the 1999 Izmit earthquake by a triggered cascade of foreshocks , 2018, Nature Geoscience.

[25]  Emanuele Casarotti,et al.  Complex Fault Geometry and Rupture Dynamics of the MW 6.5, 30 October 2016, Central Italy Earthquake , 2018 .

[26]  G. Amati,et al.  Improving the Detection of Low-Magnitude Seismicity Preceding the Mw 6.3 L'Aquila Earthquake: Development of a Scalable Code Based on the Cross Correlation of Template Earthquakes , 2018 .

[27]  L. Chiaraluce,et al.  Loading Rate Variations Along a Midcrustal Shear Zone Preceding the Mw6.0 Earthquake of 24 August 2016 in Central Italy , 2017 .

[28]  J. Avouac,et al.  Aseismic deformation associated with an earthquake swarm in the northern Apennines (Italy) , 2017 .

[29]  M. Cattaneo,et al.  The 2016 central Italy seismic sequence: a first look at the mainshocks, aftershocks, and source models , 2017 .

[30]  P. Wessel,et al.  The GMT/MATLAB Toolbox , 2017 .

[31]  Laura Scognamiglio,et al.  Slip heterogeneity and directivity of the ML 6.0, 2016, Amatrice earthquake estimated with rapid finite‐fault inversion , 2016 .

[32]  L. Chiaraluce,et al.  On the Relationship between Mw and ML for Small Earthquakes , 2016 .

[33]  I. Zaliapin,et al.  A global classification and characterization of earthquake clusters , 2016 .

[34]  A. Kato,et al.  Accelerated nucleation of the 2014 Iquique, Chile Mw 8.2 Earthquake , 2016, Scientific Reports.

[35]  Ryota Hino,et al.  Periodic slow slip triggers megathrust zone earthquakes in northeastern Japan , 2016, Science.

[36]  Lion Krischer,et al.  ObsPy: a bridge for seismology into the scientific Python ecosystem , 2015 .

[37]  Morteza Amini,et al.  Optimal partial ridge estimation in restricted semiparametric regression models , 2015, J. Multivar. Anal..

[38]  Aitaro Kato,et al.  The preparatory phase of the 2009 Mw 6.3 L'Aquila earthquake by improving the detection capability of low‐magnitude foreshocks , 2014 .

[39]  Arnaud Mignan,et al.  The debate on the prognostic value of earthquake foreshocks: A meta-analysis , 2014, Scientific Reports.

[40]  Anne Mangeney,et al.  An Automatic Kurtosis‐Based P ‐ and S ‐Phase Picker Designed for Local Seismic Networks , 2014 .

[41]  C. Chiarabba,et al.  Deep structure and tectonics of the northern‐central Apennines as seen by regional‐scale tomography and 3‐D located earthquakes , 2013 .

[42]  C. Chiarabba,et al.  Advective heat transport associated with regional Earth degassing in central Apennine (Italy) , 2013 .

[43]  D. Marsan,et al.  The long precursory phase of most large interplate earthquakes , 2013 .

[44]  Felix Waldhauser,et al.  Radiography of a normal fault system by 64,000 high‐precision earthquake locations: The 2009 L'Aquila (central Italy) case study , 2013 .

[45]  Y. Ogata,et al.  Prospective foreshock forecast experiment during the last 17 years , 2012 .

[46]  C. Collettini,et al.  A multidisciplinary study of a natural example of a CO2 geological reservoir in central Italy , 2012 .

[47]  L. Malagnini,et al.  Control of pore fluid pressure diffusion on fault failure mode: Insights from the 2009 L'Aquila seismic sequence , 2012 .

[48]  Naoshi Hirata,et al.  Propagation of Slow Slip Leading Up to the 2011 Mw 9.0 Tohoku-Oki Earthquake , 2012, Science.

[49]  Davide Piccinini,et al.  Temporal variation of seismic velocity and anisotropy before the 2009 MW 6.3 L'Aquila earthquake, Italy , 2010 .

[50]  Warner Marzocchi,et al.  On the Increase of Background Seismicity Rate during the 1997-1998 Umbria-Marche, Central Italy, Sequence: Apparent Variation or Fluid-Driven Triggering? , 2010 .

[51]  Zhigang Peng,et al.  Migration of early aftershocks following the 2004 Parkfield earthquake , 2009 .

[52]  Kate Huihsuan Chen,et al.  Characteristic repeating earthquakes in an arc-continent collision boundary zone: The Chihshang fault of eastern Taiwan , 2008 .

[53]  F. Ringdal,et al.  The detection of low magnitude seismic events using array-based waveform correlation , 2006 .

[54]  Stephen S. Gao,et al.  Temporal variation of seismic b‐values beneath northeastern Japan island arc , 2002 .

[55]  F. Waldhauser,et al.  A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California , 2000 .

[56]  Maurizio Ripepe,et al.  Foreshock sequence of September 26th, 1997 Umbria-Marche earthquakes , 2000 .

[57]  Jean Virieux,et al.  Probabilistic Earthquake Location in 3D and Layered Models , 2000 .

[58]  David Amitrano,et al.  From diffuse to localised damage through elastic interaction , 1999 .

[59]  Thomas J. Owens,et al.  The TauP Toolkit: Flexible Seismic Travel-Time and Raypath Utilities , 1999 .

[60]  Yan Y. Kagan,et al.  Long-term earthquake clustering , 1991 .

[61]  William L. Ellsworth,et al.  Monitoring velocity variations in the crust using earthquake doublets: An application to the Calaveras Fault, California , 1984 .

[62]  J. D. Eshelby.  The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.