Postseismic GRACE and GPS observations indicate a rheology contrast above and below the Sumatra slab

More than 7 years of observations of postseismic relaxation after the 2004 Sumatra-Andaman earthquake provide an improving view on the deformation in the wide vicinity of the 2004 rupture. We include both Gravity Recovery and Climate Experiment (GRACE) gravity field data that show a large postseismic signal over the rupture area and GPS observations in the back arc region. With increasing time GPS and GRACE show contrasting relaxation styles that were not easily discernible on shorter time series. We investigate whether mantle creep can simultaneously explain the far-field surface displacements and the long-wavelength gravity changes. We interpret contrasts in the temporal behavior of the GPS-GRACE observations in terms of lateral variations in rheological properties of the asthenosphere below and above the slab. Based on 1-D viscoelastic models, our results support an (almost) order of magnitude contrast between oceanic lithosphere viscosity and continental viscosity, which likely means that the low viscosities frequently found from postseismic deformation after subduction earthquakes are valid only for the mantle wedge. Next to mantle creep, we also consider afterslip as an alternative mechanism for postseismic deformation. We investigate how the combination of GRACE and GPS data can better discriminate between different mechanisms of postseismic relaxation: distributed deformation (mantle creep) versus localized deformation (afterslip). We conclude that the GRACE-observed gravity changes rule out afterslip as the dominant mechanism explaining long-wavelength deformation even over the first year after the event.

[1]  Z. Altamimi,et al.  ITRF2005 : A new release of the International Terrestrial Reference Frame based on time series of station positions and Earth Orientation Parameters , 2007 .

[2]  Richard Biancale,et al.  Separation of coseismic and postseismic gravity changes for the 2004 Sumatra–Andaman earthquake from 4.6 yr of GRACE observations and modelling of the coseismic change by normal-modes summation , 2009 .

[3]  W. Peltier,et al.  The Viscosities of the Earth's Mantle , 2013 .

[4]  G. Spada,et al.  A sea level equation for seismic perturbations , 2010 .

[5]  B. Vermeersen,et al.  Global Dynamics of the Earth , 2004 .

[6]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[7]  S. Karato,et al.  Water distribution across the mantle transition zone and its implications for global material circulation , 2011 .

[8]  A. Lowry,et al.  Andaman Postseismic Deformation Observations: Still Slipping after All These Years? , 2012 .

[9]  F. Pollitz,et al.  Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth , 2006 .

[10]  G. Hirth,et al.  Rheology of the Upper Mantle and the Mantle Wedge: A View from the Experimentalists , 2013 .

[11]  Kelin Wang,et al.  Spherical‐Earth finite element model of short‐term postseismic deformation following the 2004 Sumatra earthquake , 2012 .

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

[13]  Scott B. Luthcke,et al.  FAST TRACK PAPER: Tide model errors and GRACE gravimetry: towards a more realistic assessment , 2006 .

[14]  Tomokazu Kobayashi,et al.  Preceding, coseismic, and postseismic slips of the 2011 Tohoku earthquake, Japan , 2012 .

[15]  Rongjiang Wang,et al.  Investigation on afterslip and steady state and transient rheology based on postseismic deformation and geoid change caused by the Sumatra 2004 earthquake , 2011 .

[16]  Shin-Chan Han,et al.  Broadscale postseismic gravity change following the 2011 Tohoku-Oki earthquake and implication for deformation by viscoelastic relaxation and afterslip , 2014, Geophysical research letters.

[17]  B. Hager,et al.  The effects of rheological layering on post-seismic deformation , 2006 .

[18]  W. Wal,et al.  Ocean contribution to co-seismic crustal deformation and geoid anomalies: Application to the 2004 December 26 Sumatra-Andaman earthquake , 2011 .

[19]  T. Herring,et al.  Steady-state laboratory flow laws alone fail to explain postseismic observations , 2010 .

[20]  I. Jackson,et al.  The seismological signature of temperature and grain size variations in the upper mantle , 2005 .

[21]  K. Heki,et al.  Slow postseismic recovery of geoid depression formed by the 2004 Sumatra‐Andaman Earthquake by mantle water diffusion , 2007 .

[22]  David A. Yuen,et al.  Normal modes of the viscoelastic earth , 1982 .

[23]  K. Heki,et al.  Long‐ and short‐term postseismic gravity changes of megathrust earthquakes from satellite gravimetry , 2014 .

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

[25]  Byron D. Tapley,et al.  GRACE detects coseismic and postseismic deformation from the Sumatra‐Andaman earthquake , 2007 .

[26]  W. Benoit,et al.  High temperature anelastic and viscoplastic deformation of fine-grained MgO-doped Al2O3 , 1998 .

[27]  R. Bilham,et al.  Partial and Complete Rupture of the Indo-Andaman Plate Boundary 1847-2004 , 2005 .

[28]  L. Vermeersen,et al.  A new class of stratified viscoelastic models by analytical techniques , 1997 .

[29]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[30]  The Elasticity Theory of Dislocations in Real Earth Models and Changes in the Rotation of the Earth , 1971 .

[31]  F. Pollitz,et al.  Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M= 9.2 Sumatra earthquake , 2008 .

[32]  Kelin Wang,et al.  Deformation cycles of subduction earthquakes in a viscoelastic Earth , 2012, Nature.

[33]  R. Sabadini,et al.  GRACE gravity data help constraining seismic models of the 2004 Sumatran earthquake , 2011 .

[34]  J. C. Savage,et al.  Postseismic relaxation and aftershocks , 2007 .

[35]  R. Hyndman,et al.  The thermal structure of subduction zone back arcs , 2004 .

[36]  Walter D. Mooney,et al.  Poroelastic stress-triggering of the 2005 M8.7 Nias earthquake by the 2004 M9.2 Sumatra–Andaman earthquake , 2010 .

[37]  Motoyuki Kido,et al.  Prevalence of viscoelastic relaxation after the 2011 Tohoku-oki earthquake , 2014, Nature.

[38]  S. Bettadpur.  Insights into the Earth System mass variability from CSR-RL05 GRACE gravity fields , 2012 .

[39]  Magali I. Billen,et al.  A low viscosity wedge in subduction zones , 2001 .

[40]  V. Gahalaut,et al.  Constraints on 2004 Sumatra Andaman earthquake rupture from GPS measurements in Andaman Nicobar Islands , 2006 .

[41]  P. Chopra.  High-temperature transient creep in olivine rocks , 1997 .

[42]  Hugo Perfettini,et al.  Postseismic relaxation driven by brittle creep: A possible mechanism to reconcile geodetic measurements and the decay rate of aftershocks, application to the Chi-Chi earthquake, Taiwan , 2004 .

[43]  Fred F. Pollitz,et al.  Upper mantle rheology from GRACE and GPS postseismic deformation after the 2004 Sumatra‐Andaman earthquake , 2010 .

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

[45]  Shin-Chan Han,et al.  Source parameter inversion for recent great earthquakes from a decade‐long observation of global gravity fields , 2013 .

[46]  J. Wahr,et al.  The viscoelastic relaxation of a realistically stratified earth, and a further analysis of postglacial rebound , 1995 .

[47]  R. Bilham,et al.  Postseismic deformation of the Andaman Islands following the 26 December, 2004 Great Sumatra–Andaman earthquake , 2007 .

[48]  Antonio Piersanti,et al.  Global post-seismic deformation , 1995 .

[49]  G. Spada,et al.  Compressible rotational deformation , 1996 .

[50]  R. Riva,et al.  Ocean contribution to seismic gravity changes: the sea level equation for seismic perturbations revisited , 2014 .

[51]  Timothy H. Dixon,et al.  Lateral variation in upper mantle viscosity: role of water , 2004 .

[52]  R. Riva,et al.  Relating viscosities from postseismic relaxation to a realistic viscosity structure for the lithosphere , 2009 .

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

[54]  Chen Ji,et al.  Implications of postseismic gravity change following the great 2004 Sumatra-Andaman earthquake from the regional harmonic analysis of GRACE intersatellite tracking data , 2008 .

[55]  Chen Ji,et al.  Coseismic Slip and Afterslip of the Great Mw 9.15 Sumatra–Andaman Earthquake of 2004 , 2007 .

[56]  W. Peltier,et al.  Post-glacial rebound and transient lower mantle rheology , 1986 .

[57]  T. Okada,et al.  Two‐dimensional viscosity structure of the northeastern Japan islands arc‐trench system , 2013 .

[58]  R. Cooper,et al.  Low-frequency shear attenuation in polycrystalline olivine: Grain boundary diffusion and the physical significance of the Andrade model for viscoelastic rheology , 1998 .

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

[60]  C. Scholz.  Earthquakes and friction laws , 1998, Nature.

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

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

[63]  Ari Ben-Menahem,et al.  Seismic waves and sources , 1981 .

[64]  Chalermchon Satirapod,et al.  Insight into the 2004 Sumatra–Andaman earthquake from GPS measurements in southeast Asia , 2005, Nature.

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

[66]  D. Wolff,et al.  Viscoelastic relaxation of a Burgers half-space: implications for the interpretation of the Fennoscandian uplift , 1996 .

[67]  R. Sabadini.  Applications of Normal Mode Relaxation Theory to Solid Earth Geophysics , 2004 .

[68]  Boudewijn Ambrosius,et al.  A decade of GPS in Southeast Asia: Resolving Sundaland motion and boundaries , 2007 .

[69]  Kelin Wang,et al.  Three-dimensional viscoelastic finite element model for postseismic deformation of the great 1960 Chile earthquake , 2004 .

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