Can tectonic processes be recovered from new gravity satellite data

Abstract The goal of this study is to investigate whether temporal variations of the gravity field caused by tectonic processes (hereafter geodynamic signals) can be recognized in satellite gravity data, including the currently operating GRACE satellites and future systems. We restricted our study to subduction zones, calculating possible gravity field variations associated with elastic stress accumulation in locked areas and with stress release by earthquakes. We used fault-plane solutions for the Alaska-1964, Chile-1960 and Hokkaido-2003 earthquakes, and GPS-based strain accumulation data in locked areas of the Alaska subduction zone. Vertical displacements of the Earth's surface were calculated using a model of a rectangular fault in an elastic half-space. We developed and applied a statistical signal-recognition technique to identify signals caused by displacements of unknown magnitude on fault planes of given position and dimension. Our goal is thus to detect and analyse in satellite gravity data a signal constrained by ground geophysical and geodetic data. We assumed different levels of data accuracy, ranging from the first GRACE model GGM-01S to two orders of magnitude lower, corresponding to the target accuracy for GRACE and GOCE data. We concluded that using the developed technique, gravity field variations similar to those caused by Alaska-1964 earthquake should be recognizable in GRACE data at the accuracy level of the model GGM-01S. If forthcoming satellite gravity models have an accuracy one order of magnitude better, then the signal recognition probability will be about 99% using our approach. The required accuracy is close to the errors due to imperfect corrections for atmospheric effects. For the Chile-1960 earthquake we considered different fault-plane models and found that one can distinguish between these models with a probability approaching 70% at present level of GRACE accuracy. Increasing the data accuracy by one order of magnitude makes this probability very high. Because the gravity signal from the Hokkaido-2003 earthquake was rather weak, it would only be recognized if the data accuracy increases by two orders of magnitude thus approaching the target GRACE accuracy. If forthcoming gravity models are one order of magnitude more accurate compared to the first GRACE model then 5 years of data will allow recognition of time varying gravity signal associated with locked areas of the Alaska subduction zone. Our method may be easily applied to other geodynamic targets and more generally be adapted to other time varying gravity studies.

[1]  Kosuke Heki,et al.  Crustal velocity field of southwest Japan: Subduction and arc‐arc collision , 2001 .

[2]  Deformation across the rupture zone of the 1964 Alaska earthquake , 1998 .

[3]  James C. Savage,et al.  Mechanism of the Chilean Earthquakes of May 21 and 22, 1960 , 1970 .

[4]  G. Ramillien,et al.  Global time-variations of hydrological signals from GRACE satellite gravimetry , 2004 .

[5]  Takuo Maruyama,et al.  Static elastic dislocation in an infinite and semi-infinite medium , 1964 .

[6]  C. Readings,et al.  Gravity field and steady-state ocean circulation mission , 1996 .

[7]  Wenke Sun,et al.  Coseismic deformations detectable by satellite gravity missions: A case study of Alaska (1964, 2002) and Hokkaido (2003) earthquakes in the spectral domain , 2004 .

[8]  Nick Kusznir,et al.  Effects of rigidity layering, gravity and stress relaxation on 3-D subsurface fault displacement fields , 1994 .

[9]  J. Freymueller,et al.  Three-dimensional elastic dislocation modeling of the postseismic response to the 1964 Alaska earthquake , 2002 .

[10]  G. Spada,et al.  Toroidal/poloidal partitioning of global post-seismic deformation , 1995 .

[11]  Stephen E. Fienberg,et al.  Testing Statistical Hypotheses , 2005 .

[12]  Wenke Sun,et al.  Surface potential and gravity changes due to internal dislocations in a spherical earth—I. Theory for a point dislocation , 1993 .

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

[14]  S. Miyazaki,et al.  A slow thrust slip event following the two 1996 Hyuganada Earthquakes beneath the Bungo Channel, southwest Japan , 1999 .

[15]  J. C. Savage,et al.  Surface deformation associated with dip‐slip faulting , 1966 .

[16]  J. Steketee.  ON VOLTERRA'S DISLOCATIONS IN A SEMI-INFINITE ELASTIC MEDIUM , 1958 .

[17]  Gary S. E. Lagerloef,et al.  Satellite Gravity and the Geosphere: Contributions to the Study of the Solid Earth and Its Fluid Earth , 1998 .

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

[19]  Fred F. Pollitz,et al.  Gravitational viscoelastic postseismic relaxation on a layered spherical Earth , 1997 .

[20]  M. Drinkwater,et al.  GOCE: ESA’s First Earth Explorer Core Mission , 2003 .

[21]  T. Hori,et al.  Deformation and stress localization at the Nankai subduction zone, southwest Japan , 2003 .

[22]  J. Wahr,et al.  Postglacial rebound and Earth's viscosity structure from GRACE , 2002 .

[23]  Wenke Sun,et al.  Surface potential and gravity changes due to internal dislocations in a spherical earth—II. Application to a finite fault , 2002 .

[24]  F. Bryan,et al.  Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE , 1998 .

[25]  M. Vallée,et al.  Imaging coseismic rupture in far field by slip patches , 2004 .

[26]  S. Beck,et al.  The rupture process and tectonic implications of the great 1964 Prince William Sound earthquake , 1994 .

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

[28]  J. Freymueller,et al.  Spatial variations in present‐day deformation, Kenai Peninsula, Alaska, and their implications , 2000 .

[29]  J. Wahr,et al.  Can surface pressure be used to remove atmospheric , 2001 .

[30]  J. C. Savage.  A dislocation model of strain accumulation and release at a subduction zone , 1983 .