NLSlocate: Nonlinear Least Squares Location

In Rquake: Seismic Hypocenter Determination

Nonlinear Least Squares Location

Description

Nonlinear Least Squares Location using Gieger's method

Usage

NLSlocate(GH, vel = list(), init = c(0, 0, 0, 0), PLOT = FALSE)

Arguments

GH	List, RSEIS
vel	velocity model
init	initial guess for event location
PLOT	logical, TRUE=plot

Details

This is an adaptation of non-linear least squares inversion for earthquake location. A residual function is supplied, and iterations are performed until the location is determined.

Value

vector, new location

Note

At this stage there are no weighting mechanisms or code to eliminate data that has residuals that are too large.

Author(s)

Jonathan M. Lees<[e-mail]>

References

Lee, W.H.K., and S.W. Stewart, Principles and Applications of Microearthquake Networks, Academic Press, New York, 1981.

See Also

swig

Examples

data(GH, package='RSEIS')
###  location is:
eqsol = NLSlocate(GH, vel=GH$velfile,  PLOT=TRUE )

Rquake documentation built on June 22, 2024, 11:54 a.m.