Lithospheric rheology constrained from twenty-five years of postseismic deformation following the 1989 Mw 6.9 Loma Prieta earthquake

Abstract The October 17, 1989 M w 6.9 Loma Prieta earthquake provides the first opportunity of probing the crustal and upper mantle rheology in the San Francisco Bay Area since the 1906 M w 7.9 San Francisco earthquake. Here we use geodetic observations including GPS and InSAR to characterize the Loma Prieta earthquake postseismic displacements from 1989 to 2013. Pre-earthquake deformation rates are constrained by nearly 20 yr of USGS trilateration measurements and removed from the postseismic measurements prior to the analysis. We observe GPS horizontal displacements at mean rates of 1–4 mm/yr toward Loma Prieta Mountain until 2000, and ∼2 mm/yr surface subsidence of the northern Santa Cruz Mountains between 1992 and 2002 shown by InSAR, which is not associated with the seasonal and longer-term hydrological deformation in the adjoining Santa Clara Valley. Previous work indicates afterslip dominated in the early (1989–1994) postseismic period, so we focus on modeling the postseismic viscoelastic relaxation constrained by the geodetic observations after 1994. The best fitting model shows an elastic 19-km-thick upper crust above an 11-km-thick viscoelastic lower crust with viscosity of ∼ 6 × 10 18  Pa s , underlain by a viscous upper mantle with viscosity between 3 × 10 18 and 2 × 10 19  Pa s . The millimeter-scale postseismic deformation does not resolve the viscosity in the different layers very well, and the lower-crustal relaxation may be localized in a narrow shear zone. However, the inferred lithospheric rheology is consistent with previous estimates based on post-1906 San Francisco earthquake measurements along the San Andreas fault system. The viscoelastic relaxation may also contribute to the enduring increase of aseismic slip and repeating earthquake activity on the San Andreas fault near San Juan Bautista, which continued for at least a decade after the Loma Prieta event.

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

[2]  P. Segall,et al.  Lower crustal structure in northern California: Implications from strain rate variations following the 1906 San Francisco earthquake , 2003 .

[3]  Roland Bürgmann,et al.  Spatial variations in slip deficit on the central San Andreas Fault from InSAR , 2008 .

[4]  L. Chiao,et al.  Seismic Quiescence before the 1999 Chi-Chi, Taiwan, Mw 7.6 Earthquake , 2006 .

[5]  George E. Hilley,et al.  Resolving vertical tectonics in the San Francisco Bay Area from permanent scatterer InSAR and GPS analysis , 2006 .

[6]  F. Pollitz,et al.  Temporal evolution of continental lithospheric strength in actively deforming regions , 2008 .

[7]  R. Bürgmann,et al.  When faults communicate: Viscoelastic coupling and earthquake clustering in a simple two‐fault system , 2003 .

[8]  Roland Bürgmann,et al.  Evidence of power-law flow in the Mojave desert mantle , 2004, Nature.

[9]  Georg Dresen,et al.  Rheology of the Lower Crust and Upper Mantle: Evidence from Rock Mechanics, Geodesy, and Field Observations , 2008 .

[10]  A. Levander,et al.  Deformation in the Lower Crust of the San Andreas Fault System in Northern California , 1997 .

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

[12]  E. Fielding,et al.  Predictability of hydraulic head changes and characterization of aquifer‐system and fault properties from InSAR‐derived ground deformation , 2014 .

[13]  R. Bürgmann,et al.  Probing the lithospheric rheology across the eastern margin of the Tibetan Plateau , 2014 .

[14]  J. Avouac,et al.  Postseismic Deformation Following the 2010 $$M = 7.2$$M=7.2 El Mayor-Cucapah Earthquake: Observations, Kinematic Inversions, and Dynamic Models , 2015, Pure and Applied Geophysics.

[15]  D. Sandwell,et al.  Vertical crustal displacement due to interseismic deformation along the San Andreas fault: Constraints from tide gauges , 2014 .

[16]  Sylvain Barbot,et al.  Evidence for postseismic deformation of the lower crust following the 2004 Mw6.0 Parkfield earthquake , 2011 .

[17]  John R. Rice,et al.  Crustal Earthquake Instability in Relation to the Depth Variation of Frictional Slip Properties , 1986 .

[18]  Paul Segall,et al.  Postseismic strain following the 1989 Loma Prieta earthquake from GPS and leveling measurements , 1997 .

[19]  Greg Hirth,et al.  Using short-term postseismic displacements to infer the ambient deformation conditions of the upper mantle , 2011 .

[20]  R. Kovach,et al.  Long‐term creep‐rate changes and their causes , 2003 .

[21]  W. Behr,et al.  Rheological properties of the mantle lid beneath the Mojave region in southern California , 2014 .

[22]  F. Pollitz.  Transient Rheology of the uppermost mantle beneath the Mojave Desert , 2003 .

[23]  S. Yoo,et al.  Joint inversion of seismic and geodetic data for the source of the 2010 March 4, Mw 6.3 Jia-Shian, SW Taiwan, earthquake , 2013 .

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

[25]  Sylvain Barbot,et al.  Postseismic deformation following the 1999 Chi-Chi earthquake, Taiwan: Implication for lower-crust rheology , 2012 .

[26]  Yehuda Ben-Zion,et al.  Interaction of the San Andreas Fault Creeping Segment with Adjacent great rupture zones and earthquake recurrence at Parkfield , 1993 .

[27]  R. Bürgmann,et al.  Creep and quakes on the northern transition zone of the San Andreas fault from GPS and InSAR data , 2005 .

[28]  Tim J. Wright,et al.  Reassessing the 2006 Guerrero slow‐slip event, Mexico: Implications for large earthquakes in the Guerrero Gap , 2015 .

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

[30]  Fred F. Pollitz,et al.  Joint estimation of afterslip rate and postseismic relaxation following the 1989 Loma Prieta earthquake , 1998 .

[31]  F. Pollitz,et al.  A physical model for strain accumulation in the San Francisco Bay region: Stress evolution since 1838 , 2004 .

[32]  Howard A. Zebker,et al.  Phase unwrapping for large SAR interferograms: statistical segmentation and generalized network models , 2002, IEEE Trans. Geosci. Remote. Sens..

[33]  J. C. Savage,et al.  Postseismic relaxation following the 1989 MS7.1 Loma Prieta earthquake, central California , 2010 .

[34]  J. C. Savage,et al.  The velocity field along the San Andreas Fault in central and southern California , 1991 .

[35]  D. Schmidt.  Time-dependent land uplift and subsidence in the Santa Clara Valley , 2003 .

[36]  Fred F. Pollitz,et al.  Postseismic relaxation theory on the spherical earth , 1992 .

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

[38]  Paul Rosen,et al.  Postseismic Rebound in Fault Step-Overs Caused by Pore Fluid Flow , 1996, Science.

[39]  Postseismic deformation following the 1989 (M = 7.1) Loma Prieta, California, earthquake , 1994 .

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

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

[42]  F. Pollitz.  Postearthquake relaxation evidence for laterally variable viscoelastic structure and water content in the Southern California mantle , 2015 .

[43]  M. A. D'Alessio,et al.  Slicing up the San Francisco Bay Area : Block kinematics and fault slip rates from GPS-derived surface velocities , 2005 .

[44]  K. Johnson,et al.  New methods for estimating the spatial distribution of locked asperities and stress‐driven interseismic creep on faults with application to the San Francisco Bay Area, California , 2010 .

[45]  Malcolm Sambridge,et al.  Transdimensional inversion of receiver functions and surface wave dispersion , 2012 .

[46]  W. Thatcher,et al.  Faulting geometry and slip from co-seismic elevation changes: The 18 October 1989, Loma Prieta, California, earthquake , 1991 .

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

[48]  R. Bürgmann,et al.  Aseismic slip and fault interaction from repeating earthquakes in the Loma Prieta aftershock zone , 2013 .

[49]  Henriette Sudhaus,et al.  Improved source modelling through combined use of InSAR and GPS under consideration of correlated data errors: application to the June 2000 Kleifarvatn earthquake, Iceland , 2009 .

[50]  P. Segall,et al.  The 1989 Loma Prieta earthquake imaged from inversion of geodetic data , 1994 .

[51]  Sylvain Barbot,et al.  A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow , 2010 .

[52]  R. Bürgmann.  Active detachment faulting in the San Francisco Bay area , 1997 .

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