HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Full Text (PDF) References Citation Map Services Email this article to a friend Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Hartzell, S. Articles by Liu, P. Search for Related Content GeoRef GeoRef Citation

Determination of earthquake source parameters using a hybrid global search algorithm

U.S. Geological Survey Branch of Earthquake and Landslide Hazards, Denver, Colorado 80225
Institute of Geophysics Chinese Academy of Sciences, Beijing 100101, People's Republic of China

Abstract

A hybrid global optimization algorithm using a combination of simulated annealing and downhill simplex methods is used to invert teleseismic body waves for earthquake source parameters. Time-domain constraints on the source-time function, inversion for double-couple parameters rather than moment tensor elements, and consideration of multiple sources lead to a nonlinear and multimodal problem, in which the objective function contains many local minima. Traditional approaches that linearize the problem and use iterative least squares are dependent on the starting model and the order in which multiple sources are processed, and can converge to a local minimum. Also, grid searches are impractical for the number of parameters we treat simultaneously. The hybrid global method provides an attractive alternative because it converges to the global minimum of a prescribed objective function and can be used to invert for multiple sources simultaneously. Model parameter constraints are easy to incorporate into the global search process, if they are desired.

A multiple point-source parameterization of the 1992 Landers, California, earthquake is used as an example of the inversion method. Rupture in this earthquake occurred along three main fault segments: the Johnson Valley fault, the Homestead fault, and the Emerson/Rock Creek faults. We invert for the strike, dip, and rake of each source, the time separation of the sources, and the farfield source-time function of each source. Parameters obtained from the inversion are consistent with field observations and the results of other investigations. Strikes vary consistently with the trend of the surface ruptures. The dip is near 90° along the entire length of the rupture, and the rake is nearly pure right-lateral strike slip. The estimated moment is 7.6 x 10²⁶ dyne-cm. The continuity of the moment release across the junction of the Homestead and Emerson faults suggests that these two faults may be more continuous at depth than at the surface. Moment release in the transition from the Johnson Valley to the Homestead faults is complex, indicating a complicated time history of faulting and a less direct relationship between the primary mapped faults.

This article has been cited by other articles:

				Site Response, Shallow Shear-Wave Velocity, and Wave Propagation at the San Jose, California, Dense Seismic Array Bulletin of the Seismological Society of America, February 1, 2003; 93(1): 443 - 464.

				Space and Time Evolution of Rupture and Faulting during the 1999 Izmit (Turkey) Earthquake Bulletin of the Seismological Society of America, February 1, 2002; 92(1): 256 - 266.