Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan

We have applied two Monte Carlo optimization techniques, simulated annealing and random cost, to the inversion of deformation data for fault and magma chamber geometry. These techniques involve an element of randomness that permits them to escape local minima and ultimately converge to the global minimum of misfit space. We have tested the Monte Carlo algorithms on two synthetic data sets. We have also compared them to one another in terms of their efficiency and reliability. We have applied the bootstrap method to estimate confidence intervals for the source parameters, including the correlations inherent in the data. Additionally, we present methods that use the information from the bootstrapping procedure to visualize the correlations between the different model parameters. We have applied these techniques to GPS, tilt, and leveling data from the March 1997 earthquake swarm off of the Izu Peninsula, Japan. Using the two Monte Carlo algorithms, we have inferred two sources, a dike and a fault, that fit the deformation data and the patterns of seismicity and that are consistent with the regional stress field.

[1]  Yehuda Bock,et al.  Southern California permanent GPS geodetic array: Error analysis of daily position estimates and site velocities , 1997 .

[2]  Howard A. Zebker,et al.  A shallow‐dipping dike fed the 1995 flank eruption at Fernandina Volcano, Galápagos, observed by satellite radar interferometry , 1999 .

[3]  T. Dixon,et al.  Noise in GPS coordinate time series , 1999 .

[4]  M. Murray,et al.  The 1992 M = 7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust , 1996 .

[5]  M. Ukawa.  Collision and fan‐shaped compressional stress pattern in the Izu Block at the northern edge of the Philippine Sea Plate , 1991 .

[6]  Steven N. Ward,et al.  An inversion for slip distribution and fault shape from geodetic observations of the 1983, Borah Peak, Idaho, Earthquake , 1986 .

[7]  John Langbein,et al.  Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS , 1997 .

[8]  K. Feigl,et al.  Radar interferometry and its application to changes in the Earth's surface , 1998 .

[9]  Y. Okada,et al.  A Model for the 1989 Seismo-Volcanic Activity off Ito, Central Japan, Derived from Crustal Movement Data , 1991 .

[10]  Bernd A. Berg,et al.  Locating global minima in optimization problems by a random-cost approach , 1993, Nature.

[11]  P. Segall,et al.  Resolving the discrepancy between geodetic and seismic fault models for the 1989 Loma Prieta, California, earthquake , 1992, Bulletin of The Seismological Society of America (BSSA).

[12]  Frank K. Wyatt,et al.  Shallow borehole tilt: A reprise , 1988 .

[13]  Kenneth W. Hudnut,et al.  Southern California Permanent GPS Geodetic Array: Continuous measurements of regional crustal deformation between the 1992 Landers and 1994 Northridge earthquakes , 1997 .

[14]  Daniel H. Rothman,et al.  Nonlinear inversion, statistical mechanics, and residual statics estimation , 1985 .

[15]  S. Owen,et al.  Rapid Deformation of the South Flank of Kilauea Volcano, Hawaii , 1995, Science.

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

[17]  K. Mogi.  Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .

[18]  Shimada,et al.  Imaging magma transport during the 1997 seismic swarm off the izu peninsula, japan , 1999, Science.

[19]  Y. Okada.  Surface deformation due to shear and tensile faults in a half-space , 1985 .

[20]  T. Kato,et al.  Crustal strains in the Japanese Islands as deduced from dense GPS array , 1998 .

[21]  Jeffrey T. Freymueller,et al.  Kinematics of the Pacific-North America plate boundary zone , 1999 .

[22]  Michael Creutz,et al.  Monte Carlo Study of Quantized SU(2) Gauge Theory , 1980 .

[23]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[24]  L N Frazer,et al.  Rapid Determination of the Critical Temperature in Simulated Annealing Inversion , 1990, Science.

[25]  M. Bevis,et al.  January 30, 1997 eruptive event on Kilauea Volcano, Hawaii, as monitored by continuous GPS , 2000 .