Logarithmic Model Joint Inversion Method for Coseismic and Postseismic Slip: Application to the 2017 Mw 7.3 Sarpol Zahāb Earthquake, Iran

Interferometric synthetic aperture radar (InSAR) has become an important technique for studying earthquake cycle deformation. However, due to the limited satellite revisit time, it is often difficult to fully separate coseismic and postseismic slip from InSARdata. Nevertheless, accurately estimating spatiotemporal coseismic and postseismic fault slip distribution models is important for quantifying earthquake slip, understanding the kinematics of seismogenic faulting, and evaluating seismic hazards. Here, we develop a logarithmic model‐based method (LogSIM) for the joint inversion of coseismic and postseismic fault slip using InSAR data frommultiple platforms in different orbits. Thismethod considers the nature of early postseismic slip following logarithmic decay. The coseismic slip, the decay time constant, and the decay amplitude of the logarithmic model can be jointly estimated from unwrapped interferograms without the need for InSAR time series calculations, thereby reducing the number of unknown parameters and stabilizing the inversion. The robustness of LogSIM is first validatedwith synthetic experiments.We then apply LogSIM to invert for the coseismic and postseismic slip models of the 2017Mw 7.3 Sarpol Zahāb earthquake. We find that the maximum afterslip value and postseismic moment release during the first four months reported in previous studies are underestimated by ~26% and ~31%, respectively, due to the unavailability of the first five days of postseismic deformation. The estimated afterslip distribution spatially overlaps with the aftershocks in the updip region of the fault plane. LogSIM could potentially be extended to integrate different space geodetic data in earthquake cycle deformation modeling.

[1]  Teh-Ru Alex Song,et al.  Large Trench-Parallel Gravity Variations Predict Seismogenic Behavior in Subduction Zones , 2003, Science.

[2]  Paul Segall,et al.  Space time distribution of afterslip following the 2003 Tokachi‐oki earthquake: Implications for variations in fault zone frictional properties , 2004 .

[3]  Paul Segall,et al.  Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan , 2001 .

[4]  M. Mohajjel,et al.  Structural evidence for superposition of transtension on transpression in the Zagros collision zone: Main Recent Fault, Piranshahr area, NW Iran , 2014 .

[5]  Songbo Wu,et al.  Interseismic Ground Deformation and Fault Slip Rates in the Greater San Francisco Bay Area From Two Decades of Space Geodetic Data , 2018, Journal of Geophysical Research: Solid Earth.

[6]  David T. Sandwell,et al.  El Mayor‐Cucapah (Mw 7.2) earthquake: Early near‐field postseismic deformation from InSAR and GPS observations , 2014 .

[7]  E. Shabani,et al.  Impulsive Source of the 2017 M W=7.3 Ezgeleh, Iran, Earthquake , 2019, Geophysical research letters.

[8]  Michael Bevis,et al.  Separating simultaneous postseismic processes: An application of the postseismic straightening method to the Maule 2010 cGPSv , 2016 .

[9]  Paul Segall,et al.  Rapid afterslip following the 1999 Chi‐Chi, Taiwan Earthquake , 2002 .

[10]  Remko Scharroo,et al.  Generic Mapping Tools: Improved Version Released , 2013 .

[11]  D. Lockner,et al.  Quantitative measure of the variation in fault rheology due to fluid‐rock interactions , 1998 .

[12]  Jeffrey T. Freymueller,et al.  Coupled afterslip and viscoelastic flow following the 2002 Denali Fault, Alaska earthquake , 2009 .

[13]  Kristine M. Larson,et al.  Frictional Properties on the San Andreas Fault near Parkfield, California, Inferred from Models of Afterslip following the 2004 Earthquake , 2006 .

[14]  Wanpeng Feng,et al.  Geodetic Constraints of the 2017 Mw7.3 Sarpol Zahab, Iran Earthquake, and Its Implications on the Structure and Mechanics of the Northwest Zagros Thrust‐Fold Belt , 2018, Geophysical Research Letters.

[15]  Kaihua Ding,et al.  The 2017 Mw 7.3 Ezgeleh, Iran earthquake determined from InSAR measurements and teleseismic waveforms , 2018, Geophysical Journal International.

[16]  Jessica R. Murray,et al.  Coseismic and initial postseismic deformation from the 2004 Parkfield, California, earthquake, observed by global positioning system, electronic distance meter, creepmeters, and borehole strainmeters , 2006 .

[17]  Shu-Hao Chang,et al.  Fault Slip Inverted from Surface Displacements during the 1999 Chi-Chi, Taiwan, Earthquake , 2004 .

[18]  O. Lacombe,et al.  Building the Zagros collisional orogen: Timing, strain distribution and the dynamics of Arabia/Eurasia plate convergence , 2012 .

[19]  John R. Elliott,et al.  The 2011 Mw 7.1 Van (Eastern Turkey) earthquake , 2013 .

[20]  J. Avouac,et al.  From Geodetic Imaging of Seismic and Aseismic Fault Slip to Dynamic Modeling of the Seismic Cycle , 2015 .

[21]  M. Hearne,et al.  Slab2, a comprehensive subduction zone geometry model , 2018, Science.

[22]  Guangcai Feng,et al.  Transpressional Rupture Cascade of the 2016 Mw 7.8 Kaikoura Earthquake, New Zealand , 2018 .

[23]  Gianfranco Fornaro,et al.  Space-time distribution of afterslip following the 2009 L'Aquila earthquake , 2012 .

[24]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[25]  A. Avallone,et al.  Joint Inversion of Coseismic and Early Postseismic Slip to Optimize the Information Content in Geodetic Data: Application to the 2009 Mw6.3 L'Aquila Earthquake, Central Italy , 2019, Journal of Geophysical Research: Solid Earth.

[26]  Zhenhong Li,et al.  Interferometric synthetic aperture radar atmospheric correction using a GPS-based iterative tropospheric decomposition model , 2018 .

[27]  Eric J. Fielding,et al.  Shallow fault-zone dilatancy recovery after the 2003 Bam earthquake in Iran , 2009, Nature.

[28]  B. E. Shaw,et al.  Afterslip and aftershocks in the rate‐and‐state friction law , 2007, physics/0703249.

[29]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

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

[31]  Andrea Donnellan,et al.  GPS observations of fault afterslip and upper crustal deformation following the Northridge earthquake , 1998 .

[32]  Shui-Beih Yu,et al.  Spatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake , 2004 .

[33]  Eric J. Fielding,et al.  Coseismic and Postseismic Slip of the 2004 Parkfield Earthquake from Space-Geodetic Data , 2006 .

[34]  J. Dieterich Modeling of rock friction: 1. Experimental results and constitutive equations , 1979 .

[35]  Coseismic slip and afterslip of the 2015 Mw 8.3 Illapel (Chile) earthquake determined from continuous GPS data , 2016 .

[36]  Jan Dettmer,et al.  Trans-dimensional finite-fault inversion , 2014 .

[37]  Peltzer,et al.  Evidence of Nonlinear Elasticity of the Crust from the Mw7.6 Manyi (Tibet) Earthquake. , 1999, Science.

[38]  H. Zebker,et al.  Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements , 2002 .

[39]  Yehuda Bock,et al.  Frictional Afterslip Following the 2005 Nias-Simeulue Earthquake, Sumatra , 2006, Science.

[40]  Xiaoli Ding,et al.  The 2017 Mw 7.3 Sarpol Zahāb Earthquake, Iran: A compact blind shallow-dipping thrust event in the mountain front fault basement , 2018, Tectonophysics.

[41]  Thomas H. Heaton,et al.  The slip history of the 1994 Northridge, California, earthquake determined from strong-motion, teleseismic, GPS, and leveling data , 1996, Bulletin of the Seismological Society of America.

[42]  C. Ji,et al.  Coseismic and potential early afterslip distribution of the 2009 Mw 6.3 L'Aquila, Italy earthquake , 2014 .

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

[44]  Michael Bevis,et al.  Separating rapid relocking, afterslip, and viscoelastic relaxation: An application of the postseismic straightening method to the Maule 2010 cGPS , 2016 .

[45]  Chris Marone,et al.  On the mechanics of earthquake afterslip , 1991 .

[46]  J. Freymueller,et al.  Coseismic and Early Postseismic Deformation of the 5 January 2013 Mw 7.5 Craig Earthquake from Static and Kinematic GPS Solutions , 2015 .

[47]  J. Elliott,et al.  The 12 November 2017 Mw 7.3 Ezgeleh‐Sarpolzahab (Iran) Earthquake and Active Tectonics of the Lurestan Arc , 2019, Journal of Geophysical Research: Solid Earth.

[48]  Gianfranco Fornaro,et al.  A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms , 2002, IEEE Trans. Geosci. Remote. Sens..

[49]  Zhenhong Li,et al.  Integration of InSAR Time-Series Analysis and Water-Vapor Correction for Mapping Postseismic Motion After the 2003 Bam (Iran) Earthquake , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[50]  Nicola D'Agostino,et al.  Coseismic and post-seismic slip of the 2009 L'Aquila (central Italy) MW 6.3 earthquake and implications for seismic potential along the Campotosto fault from joint inversion of high-precision levelling, InSAR and GPS data , 2014 .

[51]  Kenneth W. Hudnut,et al.  Co-seismic displacements of the 1994 Northridge, California, earthquake , 1996, Bulletin of the Seismological Society of America.

[52]  Roland Bürgmann,et al.  Dominant Controls of Downdip Afterslip and Viscous Relaxation on the Postseismic Displacements Following the Mw7.9 Gorkha, Nepal, Earthquake , 2017 .

[53]  Yehuda Bock,et al.  Parkfield earthquake: Stress-driven creep on a fault with spatially variable rate-and-state friction parameters , 2009 .

[54]  Shui-Beih Yu,et al.  Coseismic and postseismic deformation associated with the 2003 Chengkung, Taiwan, earthquake , 2009 .

[55]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

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

[57]  J. Vergés,et al.  Fold patterns and multilayer rheology of the Lurestan Province, Zagros Simply Folded Belt (Iran) , 2009, Journal of the Geological Society.

[58]  R. V. Sharp,et al.  Bulk-friction modeling of afterslip and the modified Omori law , 1997 .

[59]  Tim J. Wright,et al.  Normal faulting sequence in the Pumqu‐Xainza Rift constrained by InSAR and teleseismic body‐wave seismology , 2014 .

[60]  Pierre Briole,et al.  Shallow afterslip following the 2003 May 21, Mw= 6.9 Boumerdes earthquake, Algeria , 2008 .

[61]  Pierre Briole,et al.  Surface displacement of the Mw 7 Machaze earthquake (Mozambique): Complementary use of multiband InSAR and radar amplitude image correlation with elastic modelling. , 2010 .

[62]  Chien-Hsin Chang,et al.  An Interpretation of the 1999 Chi-Chi Earthquake in Taiwan Based on the Thin-Skinned Thrust Model , 2000 .

[63]  Eric J. Fielding,et al.  Coseismic Deformation Field of the Mw 7.3 12 November 2017 Sarpol-e Zahab (Iran) Earthquake: A Decoupling Horizon in the Northern Zagros Mountains Inferred from InSAR Observations , 2018, Remote. Sens..

[64]  Tomokazu Kobayashi,et al.  InSAR-derived Crustal Deformation and Reverse Fault Motion of the 2017 Iran-Iraq Earthquake in the Northwestern Part of the Zagros Orogenic Belt , 2018 .

[65]  William D. Barnhart,et al.  Ramp-flat basement structures of the Zagros Mountains inferred from co-seismic slip and afterslip of the 2017 Mw7.3 Darbandikhan, Iran/Iraq earthquake , 2018, Earth and Planetary Science Letters.

[66]  Wanpeng Feng,et al.  Patterns and mechanisms of coseismic and postseismic slips of the 2011 MW 7.1 Van (Turkey) earthquake revealed by multi-platform synthetic aperture radar interferometry , 2014 .

[67]  M. Tatar,et al.  Seismic imaging of the lithospheric structure of the Zagros mountain belt (Iran) , 2010 .

[68]  J. Avouac,et al.  Modeling afterslip and aftershocks following the 1992 Landers earthquake , 2007 .

[69]  Leon Bieber The Mechanics Of Earthquakes And Faulting , 2016 .

[70]  M. Zare,et al.  Coseismic deformation and source model of the 12 November 2017 MW 7.3 Kermanshah Earthquake (Iran–Iraq border) investigated through DInSAR measurements , 2018, International Journal of Remote Sensing.

[71]  S. Yoshioka Coseismic Slip Distribution of the 1999 Chi-Chi, Taiwan, Earthquake Deduced from Inversion Analysis of GPS Data , 2004 .

[72]  Marion Y. Thomas,et al.  Time-dependent postseismic slip following the 1978 M w 7.3 Tabas-e-Golshan, Iran earthquake revealed by over 20 years of ESA InSAR observations , 2018 .

[73]  D. Dreger,et al.  Rupture process of the 2004 Parkfield earthquake from near‐fault seismic waveform and geodetic records , 2008 .

[74]  R. Bürgmann,et al.  Time‐dependent triggered afterslip following the 1989 Loma Prieta earthquake , 2000 .

[75]  Y. Tamura,et al.  Short term afterslip in the 1994 Sanriku‐Haruka‐Oki Earthquake , 1997 .

[76]  C. Marone LABORATORY-DERIVED FRICTION LAWS AND THEIR APPLICATION TO SEISMIC FAULTING , 1998 .

[77]  A. Ruina Slip instability and state variable friction laws , 1983 .

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

[79]  Richard G. Gordon,et al.  Geologically current plate motions , 2010 .

[80]  R. Bürgmann,et al.  Dominant role of tectonic inheritance in supercontinent cycles , 2011 .

[81]  K. Johnson,et al.  Arrest and recovery of frictional creep on the southern Hayward fault triggered by the 1989 Loma Prieta, California, earthquake and implications for future earthquakes , 2011 .

[82]  Wenbin Xu,et al.  An improved geodetic source model for the 1999 Mw 6.3 Chamoli earthquake, India , 2016 .

[83]  A. Avallone,et al.  Unravelling the contribution of early postseismic deformation using sub-daily GNSS positioning , 2019, Scientific Reports.

[84]  J. Avouac,et al.  Himalayan megathrust geometry and relation to topography revealed by the Gorkha earthquake , 2016 .