Comparison of various inversion techniques as applied to the determination of a geophysical deformation model for the 1983 Borah Peak earthquake

A number of techniques are employed to overcome nonuniqueness and instability inherent in linear inverse problems. To test the factors that enter into the selection of an inversion technique for fault slip distribution, we used a penalty function with smoothness (PF + S), a damped least-squares method (DLS), damped least-squares method with a positivity constraint (DLS + P), and a penalty function with smoothness and a positivity constraint (PF + S + P) for inverting the elevation changes for slip associated with the 1983 Borah Peak earthquake. Unlike solving an ill-posed inverse problem using a gradient technique (Ward and Barrientos, 1986), we have restored the well-posed character between the elevation changes and normal slip distribution. Studies showed that the constraints based on sound understanding of the physical nature of the problem are crucial in the derivation of a meaningful solution and dictates primarily the selection of a particular inversion technique. All available geological and geophysical information were used to determine a geophysical deformation model for the earthquake. It is suggested that the PF + S + P solution for a fault length of 75 km is the preferred model. The long wavelength features in the estimated slip distribution are similar to those obtained by Ward and Barrientos (1986), whereas the shorter wavelength features differ between two solutions. The fault dips 49° to the southwest. The slipped zones deepen from the surface at the northwest to about 20-km downdip depth at the southeast. The fault extends to the southeast beyond the epicenter of the mainshock. It is also shown that only the long wavelength features of the slip distribution are well resolved. The resolution is better at shallower levels than at deeper levels. The resolution deteriorates when the deformation sources are away from the leveling lines. Smoothness constraints provide better resolution than damping does at depth. The addition of a positivity constraint significantly improves the model resolution.

[1]  J. Steketee,et al.  SOME GEOPHYSICAL APPLICATIONS OF THE ELASTICITY THEORY OF DISLOCATIONS , 1958 .

[2]  G. Backus,et al.  Numerical Applications of a Formalism for Geophysical Inverse Problems , 1967 .

[3]  P. Goldbart,et al.  Linear differential operators , 1967 .

[4]  G. Backus,et al.  The Resolving Power of Gross Earth Data , 1968 .

[5]  G. Backus,et al.  Inference from Inadequate and Inaccurate Data, III. , 1970, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Joel Franklin,et al.  Well-posed stochastic extensions of ill-posed linear problems☆ , 1970 .

[7]  G. Backus,et al.  Uniqueness in the inversion of inaccurate gross Earth data , 1970, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[8]  D. E. Smylie,et al.  The displacement fields of inclined faults , 1971, Bulletin of the Seismological Society of America.

[9]  D. Jackson.  Interpretation of Inaccurate, Insufficient and Inconsistent Data , 1972 .

[10]  R. Wiggins,et al.  The general linear inverse problem - Implication of surface waves and free oscillations for earth structure. , 1972 .

[11]  K. Rybicki.  Analysis of aftershocks on the basis of dislocation theory , 1973 .

[12]  R. Parker.  Understanding Inverse Theory , 1977 .

[13]  K. Aki.  Characterization of barriers on an earthquake fault , 1979 .

[14]  W. Thatcher.  Systematic inversion of geodetic data in central California , 1979 .

[15]  J. C. Savage,et al.  Geodolite measurements of deformation near Hollister, California, 1971–1978 , 1979 .

[16]  D. Jackson.  The use of a priori data to resolve non‐uniqueness in linear inversion , 1979 .

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

[18]  A. Tarantola,et al.  Inverse problems = Quest for information , 1982 .

[19]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[20]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

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

[22]  Diane I. Doser,et al.  Source parameters of the 28 October 1983 Borah Peak, Idaho, earthquake from body wave analysis , 1985 .

[23]  A. Dziewoński,et al.  Centroid-moment tensor solutions for 35 earthquakes in Western North America (1977-1983) , 1985 .

[24]  R. Stein,et al.  Planar high-angle faulting in the basin and range: Geodetic analysis of the 1983 Borah Peak, Idaho, earthquake , 1985 .

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

[26]  W. D. Richins,et al.  The 1983 Borah Peak, Idaho, earthquake and its aftershocks , 1987 .

[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. Ward,et al.  Comparison of the 1959 Hebgen Lake, Montana and the 1983 Borah Peak, Idaho, earthquakes from geodetic observations , 1987 .

[29]  M. Nashed.  Inner, outer, and generalized inverses in banach and hilbert spaces , 1987 .

[30]  J. Baumeister.  Stable solution of inverse problems , 1987 .

[31]  P. Segall,et al.  Earthquake deformation cycle on the San Andreas Fault near Parkfield, California , 1987 .

[32]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[33]  Stephen H. Hartzell,et al.  INVERSION FOR SLIP DISTRIBUTION USING TELESEISMIC P WAVEFORMS: NORTH PALM SPRINGS, BORAH PEAK, AND MICHOACAN EARTHQUAKES , 1988 .

[34]  Raul Madariaga,et al.  Non-Linear Reflection Tomography , 1988 .

[35]  Stephen H. Hartzell,et al.  Aftershock patterns and main shock faulting , 1988 .

[36]  Paul Segall,et al.  Displacement calculations from geodetic data and the testing of geophysical deformation models , 1988 .

[37]  Gianluca Valensise,et al.  Fault parameters and slip distribution of the 1915 Avezzano, Italy, earthquake derived from geodetic observations , 1989 .