A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California

We have developed an efficient method to determine high-resolution hypocenter locations over large distances. The location method incorporates ordinary absolute travel-time measurements and/or cross-correlation P-and S-wave differential travel-time measurements. Residuals between observed and theoretical travel-time differences (or double-differences) are minimized for pairs of earthquakes at each station while linking together all observed event-station pairs. A least-squares solu- tion is found by iteratively adjusting the vector difference between hypocentral pairs. The double-difference algorithm minimizes errors due to unmodeled velocity struc- ture without the use of station corrections. Because catalog and cross-correlation data are combined into one system of equations, interevent distances within multiplets are determined to the accuracy of the cross-correlation data, while the relative lo- cations between multiplets and uncorrelated events are simultaneously determined to the accuracy of the absolute travel-time data. Statistical resampling methods are used to estimate data accuracy and location errors. Uncertainties in double-difference locations are improved by more than an order of magnitude compared to catalog locations. The algorithm is tested, and its performance is demonstrated on two clus- ters of earthquakes located on the northern Hayward fault, California. There it col- lapses the diffuse catalog locations into sharp images of seismicity and reveals hor- izontal lineations of hypocenters that define the narrow regions on the fault where stress is released by brittle failure.

[1]  L. Geiger.  Herdbestimmung bei Erdbeben aus den Ankunftszeiten , 1910 .

[2]  A. DOUGLAS,et al.  Joint Epicentre Determination , 1967, Nature.

[3]  R. Crosson,et al.  Crustal structure modeling of earthquake data: 1. Simultaneous least squares estimation of hypocenter and velocity parameters , 1976 .

[4]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[5]  W. Ellsworth.  Three-dimensional structure of the crust and mantle beneath the island of Hawaii. , 1978 .

[6]  R. Geller,et al.  Four similar earthquakes in central California , 1980 .

[7]  S. Roecker.  Seismicity and tectonics of the Pamir-Hindu Kush region of central Asia , 1981 .

[8]  Kenneth R. Anderson,et al.  Robust earthquake location using M-estimates , 1982 .

[9]  Michael A. Saunders,et al.  Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems , 1982, TOMS.

[10]  Clifford H. Thurber,et al.  Earthquake locations and three‐dimensional crustal structure in the Coyote Lake Area, central California , 1983 .

[11]  William L. Ellsworth,et al.  Monitoring velocity variations in the crust using earthquake doublets: An application to the Calaveras Fault, California , 1984 .

[12]  A. Ito,et al.  HIGH RESOLUTION RELATIVE HYPOCENTERS OF SIMILAR EARTHQUAKES BY CROSS-SPECTRAL ANALYSIS METHOD , 1985 .

[13]  Julien Frechet,et al.  Sismogenese et doublets sismiques , 1985 .

[14]  F. Scherbaum,et al.  Cross spectral analysis of Swabian Jura (SW Germany) three-component microearthquake recordings , 1986 .

[15]  Gary L. Pavlis,et al.  Appraising earthquake hypocenter location errors: A complete, practical approach for single-event locations , 1986 .

[16]  B. Efron.  The jackknife, the bootstrap, and other resampling plans , 1987 .

[17]  Stephen D. Malone,et al.  High precision relative locations of earthquakes at Mount St. Helens, Washington , 1987 .

[18]  Alan D. Chave,et al.  On the robust estimation of power spectra, coherences, and transfer functions , 1987 .

[19]  A. Michael.  Effects of three-dimensional velocity structure on the seismicity of the 1984 Morgan Hill, California, aftershock sequence , 1988 .

[20]  Jose Pujol,et al.  Comments on the joint determination of hypocenters and station corrections , 1988 .

[21]  Joan S. Gomberg,et al.  The effect of S-wave arrival times on the accuracy of hypocenter estimation , 1990, Bulletin of the Seismological Society of America.

[22]  John C. VanDecar,et al.  Determination of teleseismic relative phase arrival times using multi-channel cross-correlation and least squares , 1990 .

[23]  Nicholas Deichmann,et al.  Rupture geometry from high-precision relative hypocentre locations of microearthquake clusters , 1992 .

[24]  M. Imoto,et al.  Subducting oceanic crusts of the Philippine Sea and Pacific plates and weak-zone-normal compression in the Kanto district, Japan , 1992 .

[25]  G. Pavlis.  Appraising relative earthquake location errors , 1992, Bulletin of the Seismological Society of America.

[26]  Urs Kradolfer,et al.  Initial reference models in local earthquake tomography , 1994 .

[27]  S. Billings.  Simulated annealing for earthquake location , 1994 .

[28]  J. E. Vidale,et al.  Variations in rupture process with recurrence interval in a repeated small earthquake , 1994, Nature.

[29]  Fred W. Klein,et al.  Deep fault plane geometry inferred from multiplet relative relocation beneath the south flank of Kilauea , 1994 .

[30]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[31]  Gregory C. Beroza,et al.  Foreshock sequence of the 1992 Landers, California, earthquake and its implications for earthquake nucleation , 1995 .

[32]  P. Shearer.  Application to the Whittier Narrows California aftershock sequence , 1997 .

[33]  A. Douglas,et al.  On the onset of P seismograms , 1997 .

[34]  Jonathan M. Lees,et al.  Multiplet analysis at Coso geothermal , 1998, Bulletin of the Seismological Society of America.

[35]  T. J. Owens,et al.  Near-surface scattering effects observed with a high-frequency phased array at Pinyon Flats, California , 1998, Bulletin of the Seismological Society of America.

[36]  A. Rubin,et al.  Streaks of microearthquakes along creeping faults , 1999, Nature.

[37]  William Menke,et al.  Using waveform similarity to constrain earthquake locations , 1999 .

[38]  F. Waldhauser,et al.  Slip‐parallel seismic lineations on the Northern Hayward Fault, California , 1999 .