• Bernoulli
  • Volume 4, Number 2 (1998), 167-184.

Correction of earthquake location estimation in a small-seismic-array system

Yosihiko Ogata, Akio Kobayashi, Naoya Mikami, Yasuaki Murata, and Koichi Katsura

Full-text: Open access


An array of regularly spaced seismic stations can estimate the location of a distant earthquake using arrival times at the stations of seismic waves generated by the earthquakes. However, the accuracy decreases as the distance to the epicentre of the earthquake from the array increases. This paper is concerned with the modification of the estimated location by removing its bias which is locally systematic but globally complex, reflecting the structure of the Earth's interior. Spline surfaces are used to model such biases. Then a Bayesian procedure is carried out not only to tune the smoothness constraints but also to select the best combination among various sums of squares of differently weighted residuals and various roughness penalties for the smoothing. Using the estimated splines of the posterior mode, the newly determined epicentre locations are transformed to confirm its practical utility. Residual distributions show that our procedure improves the modification by the conventional procedure. A spatial pattern of the residuals reveals some geophysical characteristics.

Article information

Bernoulli, Volume 4, Number 2 (1998), 167-184.

First available in Project Euclid: 26 March 2007

Permanent link to this document

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

bias compensation epicentre objective Bayesian smoothing penalized sum of squares B-spline surface


Ogata, Yosihiko; Kobayashi, Akio; Mikami, Naoya; Murata, Yasuaki; Katsura, Koichi. Correction of earthquake location estimation in a small-seismic-array system. Bernoulli 4 (1998), no. 2, 167--184.

Export citation


  • [1] Akaike, H. (1980) Likelihood and Bayes Procedure. In Bernard et al. (eds), Bayesian Statistics, pp. 143-166. Valencia: University Press.
  • [2] Bolt, B.A. (1976) Nuclear Explosions and Earthquakes. San Francisco, CA: W.H. Freeman.
  • [3] Brillinger, D.R. (1985) A maximum likelihood approach to frequency wave number analysis. IEEE Trans. Acoust. Speech Signal Processing, 33, 1076-1085.
  • [4] Carder, D.S., Tocher, D., Bufe, C., Stewart, S.W., Eisler, J. and Berg, E. (1967) Seismic wave arrivals from Longshot, 08 to 278. Bull. Seismolo. Soc. Amer., 57, 573-590.
  • [5] Fletcher, R. and Powell, M.J.D. (1963) A rapidly convergent descent method for minimization. Comput. J., 6, 163-168.
  • [6] Good, I.J. (1965) The Estimation of Probabilities. Cambridge, MA: MIT Press.
  • [7] Good, I.J. and Gaskins, R.A. (1971) Nonparametric roughness penalties for probability densities. Biometrika, 58, 255-277.
  • [8] Herrin, H. and Taggart, J. (1968) Source bias in epicenter determination. Bull. Seismol. Soc. Amer., 58, 1791-1796.
  • [9] Inoue, H. (1986) A least square smooth fitting for irregularly spaced data: finite element approach using the cubic B-spline basis. Geophysics, 51, 2051-2066.
  • [10] Jeffreys, H. and Bullen, K.E. (1940) Seismological Tables. London: British Association for Advancement of Sciences, Gray-Milne Trust.
  • [11] Kobayashi, A., Mikami, N. and Ishikawa, Y. (1993) Epicentral correction for the hypocenter of Matsushiro Seismic Array System (in Japanese with English captions for figures). Tech. Rep. Matsushiro Seismol. Obs. Japan Meteorol. Agency, 12, 1-9. Matsushiro Seismological Observatory (1984-1992) Seismol. Bull. Matsushiro Seismolog. Obs., Series 3, 1-8. Nagano City: Japan Meteorological Agency.
  • [12] Mikumo, T. (1965) Determination of phase velocity and direction of wave approach from station arrays. Bull. Disaster Prevention Res. Inst., 15, 31-43.
  • [13] Murata, Y. (1993) Estimation of optimum average surficial density from gravity data: an objective Bayesian Approach. J. Geophys. Res., 98, 12097-12109.
  • [14] Murata, Y. and Noro, H. (1993) Optimization of estimation procedure of bouguer density by use of the ABIC-minimization method (in Japanese with English summary and captions for figures). Butsuri Tansa, 46, 120-127.
  • [15] Ogata, Y., Imoto, M. and Katsura, K. (1991) 3-D spatial variation of b-values of magnitude-frequency distribution beneath the Kanto District, Japan. Geophys. J. Int., 104, 135-146.
  • [16] Ogata, Y. and Katsura, K. (1988) Likelihood analysis of spatial inhomogeneity for marked point patterns. Ann. Inst. Statist. Math., 40, 29-39.
  • [17] Ogata, Y. and Katsura, K. (1993) Analysis of temporal and spatial heterogeneity of magnitude frequency distribution inferred from earthquake catalogues. Geophys. J. Int., 113, 727-738.
  • [18] Osada, H., Kasiwabara, S., Nagai, A., Takayama, H., Wakui, S., Morisita, I. and Tanaka, Y. (1984) Matsushiro Seismic Array System (I)-Summary and hypocenter determination of near earthquakes (in Japanese with English figure captions). Tech. Rep. Matsushiro Seismol. Obs. Japan Meteorol. Agency, 5, 13-31.
  • [19] Otsuka, M. (1966) Azimuth and slowness anomalies of seismic waves measured on the Central California Seismographic Array, Part I. Observations. Bull. Seismol. Soc. Amer., 56, 223-239. Ringdal, F. (ed.) (1990) Bull. Seismol. Soc. Amer.: Regional Seismic Arrays and Nuclear Test Ban Verification, 80. El Cerrito, CA: Seismological Society of America.
  • [20] Takayama, H., Wakui, S. and Nagare, S. (1985) Matsushiro Seismic Array System (II)-Hypocenter determination of very-near earthquakes (in Japanese with English figure captions). Tech. Rep. Matsushiro Seismol. Obs. Japan Meteorol. Agency, 6, 67-72.
  • [21] Vere-Jones, D. and Smith, E.G.C. (1981) Statistics in seismology. Commun Statistics Theory Methods A, 10, 1559-1585.