Grid-search event location with non-Gaussian error models

This study employs an event location algorithm based on grid search to investigate the possibility of improving seismic event location accuracy by using non-Gaussian error models. The primary departure from the Gaussian error model that is considered is the explicit use of non-Gaussian probability distributions in defining optimal estimates of location parameters. Specifically, the ;class of generalized Gaussian distributions is considered, which leads to the minimization of L-p norms of arrival time residuals for arbitrary p >= 1. The generalized Gaussian error models are implemented both with fixed standard errors assigned to the data and with an iterative reweighting of the data on a station/phase basis. An implicit departure from a Gaussian model is also considered, namely, the use of a simple outlier rejection criterion for disassociating arrivals during the grid-search process. These various mechanisms were applied to the ISC phase picks for the IWREF reference events, and the resulting grid-search solutions were compared to the GT locations of the events as well as the ISC solutions. The results indicate that event mislocations resulting from the minimization of L-p residual norms, with p near 1, are generally better than those resulting from the conventional L-2 norm minimization (Gaussian error assumption). However, this result did not always hold for mislocations in event depth. Further, outlier rejection and iterative reweighting, applied with L-2 minimization, performed nearly as well as L-p minimization in some cases. The results of this study suggest that ISC can potentially improve its location capability with the use of global search methods and non-Gaussian error models. However, given the limitations of this study, further research, including the investigation of other statistical and optimization techniques not addressed here, is needed to assess this potential more completely. (c) 2006 Elsevier B.V. All rights reserved.