Time dependent inversion of geodetic data

The recent expansion of permanent Global Positioning System (GPS) networks provides crustal deformation data that are dense in both space and time. While considerable effort has been directed toward using these data for the determination of average crustal velocities, little attention has been given to detecting and estimating transient deformation signals. We introduce here a Network Inversion Filter for estimating the distribution of fault slip in space and time using data from such dense, frequently sampled geodetic networks. Fault slip is expanded in a spatial basis set sk(x) in which the coefficients are time varying, s(x, t) = ∑k=1M ck(t)sk(x) The temporal variation in fault slip is estimated nonparameterically by taking slip accelerations to be random Gaussian increments, so that fault slip is a sum of steady state and integrated random walk components. A state space model for the full geodetic network is formulated, and Kalman filtering methods are used for estimation. Variance parameters, including measurement errors, local benchmark motions, and temporal and spatial smoothing parameters, are estimated by maximum likelihood, which is computed by recursive filtering. Numerical simulations demonstrate that the Network Inversion Filter is capable of imaging fault slip transients, including propagating slip events. The Network Inversion Filter leads naturally to automated methods for detecting anomalous departures from steady state deformation.

[1]  Michael B. Heflin,et al.  Absolute far-field displacements from the 28 June 1992 Landers earthquake sequence , 1993, Nature.

[2]  M. Matthews.  On the estimation of fault slip in space and time , 1991 .

[3]  Yehuda Bock,et al.  Postseismic deformation following the Landers earthquake, California, 28 June 1992 , 1994, Bulletin of the Seismological Society of America.

[4]  J. C. Savage,et al.  Interseismic deformation along the San Andreas Fault in southern California , 1995 .

[5]  P. Segall,et al.  Slip Deficit on the San Andreas Fault at Parkfield, California, as Revealed by Inversion of Geodetic Data , 1986, Science.

[6]  F. Wyatt.  Displacement of surface monuments: Vertical motion , 1989 .

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

[8]  P. Segall,et al.  The co-seismic slip distribution of the Landers earthquake , 1994, Bulletin of the Seismological Society of America.

[9]  C. Striebel,et al.  On the maximum likelihood estimates for linear dynamic systems , 1965 .

[10]  Steven N. Ward,et al.  The 1960 Chile earthquake: inversion for slip distribution from surface deformation , 1990 .

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

[12]  C. Honsho,et al.  Crustal movements on Shikoku, southwestern Japan, inferred from inversion analysis of levelling data using ABIC , 1996 .

[13]  F. Wyatt.  Displacement of surface monuments: Horizontal motion , 1982 .

[14]  R. Parker.  Geophysical Inverse Theory , 1994 .

[15]  J. Langbein,et al.  Variations in fault slip and strain accumulation at Parkfield, California: Initial results using two‐color geodimeter measurements, 1984–1988 , 1990 .

[16]  K. W. Clark,et al.  Continuous GPS observations across the Hayward Fault, California, 1991-1994 , 1995 .

[17]  D. Agnew,et al.  Monument motion and measurements of crustal velocities , 1995 .

[18]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[19]  Herb Dragert,et al.  Continuous GPS monitoring of elastic strain in the Northern Cascadia Subduction Zone , 1995 .

[20]  Alan T. Linde,et al.  A slow earthquake sequence on the San Andreas fault , 1996, Nature.

[21]  Michael B. Heflin,et al.  Global geodesy using GPS without fiducial sites , 1992 .

[22]  Yehuda Bock,et al.  Crustal deformation measurements in central Japan determined by a Global Positioning System Fixed-Point Network , 1992 .

[23]  Takao Eguchi,et al.  Detection of a volcanic fracture opening in Japan using Global Positioning System measurements , 1990, Nature.

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

[25]  John Langbein,et al.  Correlated errors in geodetic time series: Implications for time‐dependent deformation , 1997 .

[26]  Kenneth W. Hudnut,et al.  Detection of crustal deformation from the Landers earthquake sequence using continuous geodetic measurements , 1993, Nature.

[27]  P. Segall,et al.  Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas Fault , 1987 .

[28]  S. Miyazaki,et al.  Establishment of the nationwide GPS array (GRAPES) and its initial results on the crustal deformation of Japan , 1996 .

[29]  D. Wald,et al.  Spatial and temporal distribution of slip for the 1992 Landers, California, earthquake , 1994, Bulletin of the Seismological Society of America.

[30]  John B. Rundle,et al.  Shallow and peripheral volcanic sources of inflation revealed by modeling two‐color geodimeter and leveling data from Long Valley Caldera, California, 1988–1992 , 1995 .

[31]  P. Segall,et al.  Estimation of depth‐dependent fault slip from measured surface deformation with application to the 1906 San Francisco Earthquake , 1993 .

[32]  William Rodi,et al.  Coseismic fault slip associated with the 1992 M w 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence , 1995 .

[33]  C. W. Roberts,et al.  Correlation of Changes in Gravity, Elevation, and Strain in Southern California , 1983, Science.

[34]  Hiromichi Tsuji,et al.  Coseismic crustal deformation from the 1994 Hokkaido‐Toho‐Oki Earthquake Monitored by a nationwide continuous GPS array in Japan , 1995 .