Coseismic Ground Deformation Reproduced through Numerical Modeling: A Parameter Sensitivity Analysis

Coseismic ground displacements detected through remote sensing surveys are often used to invert the coseismic slip distribution on geologically reliable fault planes. We analyze a well-known case study (2009 L’Aquila earthquake) to investigate how three-dimensional (3D) slip configuration affects coseismic ground surface deformation. Different coseismic slip surface configurations reconstructed using aftershocks distribution and coseismic cracks, were tested using 3D boundary element method numerical models. The models include two with slip patches that reach the surface and three models of blind normal-slip surfaces with different configurations of slip along shallowly-dipping secondary faults. We test the sensitivity of surface deformation to variations in stress drop and rock stiffness. We compare numerical models’ results with line of sight (LOS) surface deformation detected from differential SAR (Synthetic Aperture Radar) interferometry (DInSAR). The variations in fault configuration, rock stiffness and stress drop associated with the earthquake considerably impact the pattern of surface subsidence. In particular, the models with a coseismic slip patch that does not reach the surface have a better match to the line of sight coseismic surface deformation, as well as better match to the aftershock pattern, than models with rupture that reaches the surface. The coseismic slip along shallowly dipping secondary faults seems to provide a minor contribution toward surface deformation.

[1]  Patience A. Cowie,et al.  Fault slip-rate variations during crustal-scale strain localisation, central Italy , 2002 .

[2]  P. Tizzani,et al.  Coseismic Stress and Strain Field Changes Investigation Through 3‐D Finite Element Modeling of DInSAR and GPS Measurements and Geological/Seismological Data: The L'Aquila (Italy) 2009 Earthquake Case Study , 2018 .

[3]  Giovanna Cultrera,et al.  Rapid response seismic networks in Europe: lessons learnt from the L'Aquila earthquake emergency , 2011 .

[4]  Carlo Alberto Brunori,et al.  High-resolution controlled-source seismic tomography across the Middle Aterno basin in the epicentral area of the 2009, Mw 6.3, L’Aquila earthquake (central Apennines, Italy) , 2012 .

[5]  Rodolfo Puglia,et al.  Spectral models for ground motion prediction in the L'Aquila region (central Italy): evidence for stress-drop dependence on magnitude and depth , 2016 .

[6]  Peter M. Shearer,et al.  Global variations of stress drop for moderate to large earthquakes , 2009 .

[7]  David D. Pollard,et al.  Slip distributions on intersecting normal faults , 1999 .

[8]  F. Calamita,et al.  Contrasting styles of fault reactivation in curved orogenic belts: Examples from the Central Apennines (Italy) , 2011 .

[9]  Sabina Bigi,et al.  Contrasting surface active faults and deep seismogenic sources unveiled by the 2009 L’Aquila earthquake sequence (Italy) , 2013 .

[10]  E. Serpelloni,et al.  Fault geometry, coseismic-slip distribution and Coulomb stress change associated with the 2009 April 6, Mw 6.3, L’Aquila earthquake from inversion of GPS displacements , 2012 .

[11]  Christian Bignami,et al.  Finite fault inversion of DInSAR coseismic displacement of the 2009 L'Aquila earthquake (central Italy) , 2009 .

[12]  H. Fossen,et al.  The Wilson Cycle and Effects of Tectonic Structural Inheritance on Rifted Passive Margin Formation , 2018, Tectonics.

[13]  R. Basili,et al.  The role of pre-existing discontinuities in the development of extensional faults: an analog modeling perspective , 2015 .

[14]  David D. Pollard,et al.  Three-dimensional analyses of slip distributions on normal fault arrays with consequences for fault scaling , 1996 .

[15]  James Jackson,et al.  The 2009 L'Aquila earthquake (central Italy): A source mechanism and implications for seismic hazard , 2009 .

[16]  Anthony Lomax,et al.  Rupture history of the 2009 L'Aquila (Italy) earthquake from non‐linear joint inversion of strong motion and GPS data , 2009 .

[17]  N. D’Agostino,et al.  Evidence for localized active extension in the central Apennines (Italy) from global positioning system observations , 2011 .

[18]  A. Blumetti,et al.  InSAR data as a field guide for mapping minor earthquake surface ruptures: Ground displacements along the Paganica Fault during the 6 April 2009 L'Aquila earthquake , 2010 .

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

[20]  F. Radicioni,et al.  FAST TRACK PAPER: Coseismic and initial post-seismic slip of the 2009 Mw 6.3 L'Aquila earthquake, Italy, from GPS measurements , 2010 .

[21]  James Jackson,et al.  Active tectonics of the Adriatic Region , 1987 .

[22]  G. Cavinato,et al.  Structural architecture of the central Apennines: Interpretation of the CROP 11 seismic profile from the Adriatic coast to the orographic divide , 2008 .

[23]  Lauro Chiaraluce,et al.  The anatomy of the 2009 L'Aquila normal fault system (central Italy) imaged by high resolution foreshock and aftershock locations , 2011 .

[24]  E. Trasatti,et al.  Finite Element inversion of DInSAR data from the Mw 6 . 3 L ' Aquila 1 Earthquake , 2009 ( Italy ) 2 3 , 2022 .

[25]  K. Feigl,et al.  The displacement field of the Landers earthquake mapped by radar interferometry , 1993, Nature.

[26]  P. Renda,et al.  The effects of post‐orogenic extension on different scales: an example from the Apennine–Maghrebide fold‐and‐thrust belt, SW Sicily , 2003 .

[27]  Qiang Chen,et al.  3D coseismic deformations and source parameters of the 2010 Yushu earthquake (China) inferred from DInSAR and multiple-aperture InSAR measurements , 2014 .

[28]  Romano Camassi,et al.  DBMI04, il database delle osservazioni macrosismiche dei terremoti italiani utilizzate per la compilazione del catalogo parametrico CPTI04 , 2007 .

[29]  E. Yoffe,et al.  The angular dislocation , 1960 .

[30]  E. Trasatti,et al.  Finite element inversion of DInSAR data from the Mw 6.3 L'Aquila earthquake, 2009 (Italy) , 2011 .

[31]  Laura Scognamiglio,et al.  Fast Determination of Moment Tensors and Rupture History: What Has Been Learned from the 6 April 2009 L'Aquila Earthquake Sequence , 2010 .

[32]  Christian Bignami,et al.  The May 12, 2008, (Mw 7.9) Sichuan Earthquake (China): Multiframe ALOS-PALSAR DInSAR Analysis of Coseismic Deformation , 2010, IEEE Geoscience and Remote Sensing Letters.

[33]  Leon M Keer,et al.  Modeling slip zones with triangular dislocation elements , 1992 .

[34]  Maria Comninou,et al.  The angular dislocation in a half space , 1975 .

[35]  Grazia Pietrantonio,et al.  Coseismic deformation of the destructive April 6, 2009 L'Aquila earthquake (central Italy) from GPS data , 2009 .

[36]  Mario Locati,et al.  CPTI15, the 2015 version of the Parametric Catalogue of Italian Earthquakes , 2016 .

[37]  P. Galli,et al.  The Paganica Fault and Surface Coseismic Ruptures Caused by the 6 April 2009 Earthquake (L'Aquila, Central Italy) , 2009 .

[38]  Giancanio Sileo,et al.  Surface Faulting of the 6 April 2009 Mw 6.3 L'Aquila Earthquake in Central Italy , 2011 .

[39]  L. Chiaraluce,et al.  From surface geology to aftershock analysis:Constraints on the geometry of the L'Aquila 2009 seismogenic fault system , 2012 .

[40]  Giancanio Sileo,et al.  Partitioned postseismic deformation associated with the 2009 Mw 6.3 L'Aquila earthquake surface rupture measured using a terrestrial laser scanner , 2010 .

[41]  Jyr‐Ching Hu,et al.  Coseismic and postseismic slip distribution of the 2003 Mw = 6.5 Chengkung earthquake in eastern Taiwan: Elastic modeling from inversion of GPS data , 2009 .

[42]  D. Pollard,et al.  Integration of Surface Slip and Aftershocks to Constrain the 3D Structure of Faults Involved in the M 7.3 Landers Earthquake, Southern California , 2012 .

[43]  F. Calamita,et al.  Analogue modeling of positive inversion tectonics along differently oriented pre-thrusting normal faults: An application to the Central-Northern Apennines of Italy , 2014 .

[44]  E. Yaşar,et al.  Correlating sound velocity with the density, compressive strength and Young's modulus of carbonate rocks , 2004 .

[45]  Salvatore Barba,et al.  Analysis of seismological and geological observations formoderate‐size earthquakes: the Colfiorito Fault System(Central Apennines, Italy) , 2000 .

[46]  D. Pollard,et al.  Fault linkage: Three‐dimensional mechanical interaction between echelon normal faults , 1998 .

[47]  E. Patacca,et al.  Post-Tortonian mountain building in the Apennines. The role of the passive sinking of a relic lithospheric slab , 1987 .

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

[49]  Gianluca Valensise,et al.  A fresh look at the seismotectonics of the Abruzzi (Central Apennines) following the 6 April 2009 L’Aquila earthquake (Mw 6.3) , 2012 .

[50]  Lauro Chiaraluce,et al.  Unravelling the complexity of Apenninic extensional fault systems: A review of the 2009 L'Aquila earthquake (Central Apennines, Italy) , 2012 .

[51]  M. Cooke,et al.  Sensitivity of uplift patterns to dip of the San Andreas fault in the Coachella Valley, California , 2014 .

[52]  Tomokazu Kobayashi,et al.  Small-displacement linear surface ruptures of the 2016 Kumamoto earthquake sequence detected by ALOS-2 SAR interferometry , 2016, Earth, Planets and Space.

[53]  Raymon Lee Brown.  A dislocation approach to plate interaction , 1975 .

[54]  M. Cooke,et al.  Mechanical Models Suggest Fault Linkage through the Imperial Valley, California, U.S.A. , 2019, Bulletin of the Seismological Society of America.

[55]  Alberto Michelini,et al.  The 2009 L'Aquila (central Italy) MW6.3 earthquake: Main shock and aftershocks , 2009 .

[56]  Gianluca Valensise,et al.  On the complexity of surface ruptures during normal faulting earthquakes: excerpts from the 6 April 2009 L'Aquila (central Italy) earthquake ( M w 6.3) , 2013 .

[57]  Francesca Romana Cinti,et al.  The August 17, 1999 Izmit, Turkey, earthquake: slip distribution from dislocation modeling of DInSAR and surface offset , 2002 .

[58]  C. Bignami,et al.  Gravity-driven postseismic deformation following the Mw 6.3 2009 L’Aquila (Italy) earthquake , 2015, Scientific Reports.

[59]  G. Lavecchia,et al.  Coseismic ground deformation of the 6 April 2009 L'Aquila earthquake (central Italy, Mw6.3) , 2010 .

[60]  E. Miccadei,et al.  Sedimentary and tectonic evolution of Plio-Pleistocene alluvial and lacustrine deposits of Fucino Basin (central Italy) , 2002 .