The relationship between normalised horizontal-to-vertical spectral ratios (HVSRs) of microtremors and the F distribution
Ikuo Cho 1 4 Takaki Iwata 2 3
+ Author Affiliations
- Author Affiliations
1 Geological Survey of Japan, AIST, Tsukuba Central 7, Tsukuba 305-8567, Japan.
2 College of Community Development, Tokiwa University, 1-430-1 Miwa, Mito 310-8585, Japan.
3 Institute of Statistical Mathematics, 10-3 Midori-cho, Tachikawa, Tokyo 190-8562, Japan.
4 Corresponding author. Email: ikuo-chou@aist.go.jp
Exploration Geophysics 49(5) 637-646 https://doi.org/10.1071/EG17110
Submitted: 31 August 2017 Accepted: 3 September 2017 Published: 12 October 2017
Abstract
The Albarello test is a statistical test proposed in 2001 with the aim of quickly measuring microtremors (e.g. 5 min) to estimate horizontal-to-vertical spectral ratios (HVSRs). This test is used to identify the deviation between the empirical distribution of normalised (amplitude-corrected) HVSRs and the F distribution so as to identify a physically meaningful peak from numerous peaks, including local peaks inevitably generated by chance when a short-time record is analysed. A HVSR peak is regarded to be meaningful when a significant deviation is detected between the two distributions. However, a paper published in 2005 empirically checked the reliability of this test itself and concluded that the test was unreliable. In this study, we investigate the cause of the unreliability and evaluate the applicable condition. First, we show theoretically that the empirical distribution of normalised HVSRs is generally expected to resemble the F distribution. Then, we show that these two distributions are very similar to each other on the basis of both simulated and observed microtremors. This means that we tend to fail to identify the deviation (commit a type II error) when the sample size is not adequate. To be specific, a sample size on the order of 102–104 is required for a Kolmogorov–Smirnov test to avoid a type II error. Such a sample size corresponds to records on the order of 100–102 h in length for a representative parameter set for the spectral analyses. The original purpose of this test cannot be achieved because such long records are required to identify the deviation. It may be useful from a practical perspective to consider that normalised HVSRs can approximate the F distribution when developing methods to evaluate the statistical property of HVSRs.
Key words: noise, numerical, passive, statistical, surface wave, velocity.
References
Aki, K., 1957, Space and time spectra of stationary stochastic waves, with special reference to microtremors: Bulletin of the Earthquake Research Institute, University of Tokyo, 35, 415–457Albarello, D., 2001, Detection of spurious maxima in the site amplification characteristics estimated by the HVSR technique: Bulletin of the Seismological Society of America, 91, 718–724
| Detection of spurious maxima in the site amplification characteristics estimated by the HVSR technique:Crossref | GoogleScholarGoogle Scholar |
Albarello, D., and Lunedei, E., 2010, Alternative interpretations of horizontal to vertical spectral ratios of ambient vibrations: new insights from theoretical modeling: Bulletin of Earthquake Engineering, 8, 519–534
| Alternative interpretations of horizontal to vertical spectral ratios of ambient vibrations: new insights from theoretical modeling:Crossref | GoogleScholarGoogle Scholar |
Albarello, D., and Lunedei, E., 2013, Combining horizontal ambient vibration components for H/V spectral ratio estimates: Geophysical Journal International, 194, 936–951
| Combining horizontal ambient vibration components for H/V spectral ratio estimates:Crossref | GoogleScholarGoogle Scholar |
Arai, H., and Tokimatsu, K., 2004, S-wave velocity profiling by inversion of microtremor H/V spectrum: Bulletin of the Seismological Society of America, 94, 53–63
| S-wave velocity profiling by inversion of microtremor H/V spectrum:Crossref | GoogleScholarGoogle Scholar |
Bendat, J. S., and Piersol, A. G., 2010, Random data: analysis and measurement procedures: John Wiley & Sons.
Bonnefoy-Claudet, S., Köhler, A., Cornou, C., Wathelet, M., and Bard, P.-Y., 2008, Effects of Love waves on microtremor H/V ratio: Bulletin of the Seismological Society of America, 98, 288–300
| Effects of Love waves on microtremor H/V ratio:Crossref | GoogleScholarGoogle Scholar |
Capon, J., Greenfield, R. J., and Kolker, R. J., 1967, Multidimensional maximum likelihood processing of a large aperture seismic array: Proceedings of the IEEE, 55, 192–211
| Multidimensional maximum likelihood processing of a large aperture seismic array:Crossref | GoogleScholarGoogle Scholar |
Chávez-García, F. J., and Luzón, F., 2005, On the correlation of seismic microtremors: Journal of Geophysical Research, 110, B11313
| On the correlation of seismic microtremors:Crossref | GoogleScholarGoogle Scholar |
Cho, I., Tada, T., and Shinozaki, Y., 2004, A new method to determine phase velocities of Rayleigh waves from microseisms: Geophysics, 69, 1535–1551
| A new method to determine phase velocities of Rayleigh waves from microseisms:Crossref | GoogleScholarGoogle Scholar |
Cho, I., Tada, T., and Shinozaki, Y., 2006a, A generic formulation for microtremor exploration methods using three-component records from a circular array: Geophysical Journal International, 165, 236–258
| A generic formulation for microtremor exploration methods using three-component records from a circular array:Crossref | GoogleScholarGoogle Scholar |
Cho, I., Tada, T., and Shinozaki, Y., 2006b, New methods of microtremor exploration: the centerless circular array method and the two-radius method: Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France (ESG2006), Paper Number: 51.
Cornou, C., Ohrnberger, M., Boore, D., Kudo, K., and Bard, P.-Y., 2009, Derivation of structural models from ambient vibration array recordings: results from an international blind test, in P.-Y. Bard, E. Chaljub, C. Cornou, F. Cotton, and P. Gueguen, eds., ESG 2006: Third International Symposium on the Effects of Surface Geology on Seismic Motion: LCPC Editions.
Fäh, D., Stamm, G., and Havenith, H.-B., 2008, Analysis of three-component ambient vibration array measurements: Geophysical Journal International, 172, 199–213
| Analysis of three-component ambient vibration array measurements:Crossref | GoogleScholarGoogle Scholar |
Gouédard, P., Cornou, C., and Roux, P., 2008, Phase–velocity dispersion curves and small-scale geophysics using noise correlation slant stack technique: Geophysical Journal International, 172, 971–981
| Phase–velocity dispersion curves and small-scale geophysics using noise correlation slant stack technique:Crossref | GoogleScholarGoogle Scholar |
Hisada, Y., 1994, An efficient method for computing Green’s functions for a layered half-space with sources and receivers at close depths: Bulletin of the Seismological Society of America, 84, 1456–1472
Hisada, Y., 1995, An efficient method for computing Green’s functions for a layered half-space with sources and receivers at close depths (Part 2): Bulletin of the Seismological Society of America, 85, 1080–1093
Hobiger, M., Bard, P.-Y., Cornou, C., and Bihan, N. L., 2009, Single station determination of Rayleigh wave ellipticity by using the random decrement technique (RayDec): Geophysical Research Letters, 36, L14303
| Single station determination of Rayleigh wave ellipticity by using the random decrement technique (RayDec):Crossref | GoogleScholarGoogle Scholar |
Kohketsu, K., Miyake, H., Fujiwara, H., and Hashimoto, T., 2008, Progress towards a Japan integrated velocity structure model and long-period ground motion hazard map: 14th World Conference on Earthquake Engineering (14WCEE), Paper No. S10–038.
Konno, K., and Ohmachi, T., 1998, Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor: Bulletin of the Seismological Society of America, 88, 228–241
Mendecki, M. J., Bieta, B., and Mycka, M., 2014, Determination of the resonance frequency – thickness relation based on the ambient seismic noise records from Upper Silesia Coal Basin: Contemporary Trends in Geoscience, 3, 41–51
| Determination of the resonance frequency – thickness relation based on the ambient seismic noise records from Upper Silesia Coal Basin:Crossref | GoogleScholarGoogle Scholar |
Miyazaki, K., Sasada, M., and Yoshioka, T., 1996, Geology of the Makabe district, with geological sheet map at 1 : 50,0000: Geological Survey of Japan, 103 pp. [in Japanese with English abstract].
Mogi, H., Kawakami, H., and Ghayamghamian, M. R., 1995, Probability distributions of spectra and spectral ratios: Earthquake Geotechnical Engineering (Proceedings of IS-TOKYO), 573–578.
Morikawa, H., Akamatsu, J., Nishimura, K., and Onoue, K., 1998, Stochastic simulation of microseisms using theory of conditional random fields: Pure and Applied Geophysics, 151, 81–99
| Stochastic simulation of microseisms using theory of conditional random fields:Crossref | GoogleScholarGoogle Scholar |
Morikawa, H., Sawada, S., and Akamatsu, J., 2004, A method to estimate phase velocities of Rayleigh waves using microseisms simultaneously observed at two sites: Bulletin of the Seismological Society of America, 88, 228–241
Mucciarelli, M., Masi, A., Vona, M., Gallipoli, M. R., Harabaglia, P., Caputo, R., Piscitelli, S., Rizzo, E., Picozzi, M., Albarello, D., and Lizza, C., 2003, Quick survey of the possible cause of damage enhancement observed in San Giuliano after the 2002 Molise, Italy, seismic sequence: Journal of Earthquake Engineering, 7, 599–614
| Quick survey of the possible cause of damage enhancement observed in San Giuliano after the 2002 Molise, Italy, seismic sequence:Crossref | GoogleScholarGoogle Scholar |
Nakamura, Y., 1989, A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface: Quarterly Report of Railway Technical Research Institute, 30, 25–30
Oliveto, A., Mucciarelli, M., and Caputo, R., 2004, HVSR prospecting in multi-layered environments: an example from the Tyrnavos Basin (Greece): Journal of Seismology, 8, 395–406
| HVSR prospecting in multi-layered environments: an example from the Tyrnavos Basin (Greece):Crossref | GoogleScholarGoogle Scholar |
Picozzi, M., Parolai, S., and Albarello, D., 2005, Statistical analysis of noise horizontal-to-vertical spectral ratios (HVSR): Bulletin of the Seismological Society of America, 95, 1779–1786
| Statistical analysis of noise horizontal-to-vertical spectral ratios (HVSR):Crossref | GoogleScholarGoogle Scholar |
Prieto, G. A., Lawrence, J. F., and Beroza, G. C., 2009, Anelastic Earth structure from the coherency of the ambient seismic field: Journal of Geophysical Research, 114, B07303
| Anelastic Earth structure from the coherency of the ambient seismic field:Crossref | GoogleScholarGoogle Scholar |
Scherbaum, F., Hinzen, K.-G., and Ohrnberger, M., 2003, Determination of shallow shear wave velocity profiles in the Cologne, Germany, area using ambient vibrations: Geophysical Journal International, 152, 597–612
| Determination of shallow shear wave velocity profiles in the Cologne, Germany, area using ambient vibrations:Crossref | GoogleScholarGoogle Scholar |
Senna, S., Kudo, N., and Fujiwara, H., 2008, Velocity structure model for Tsukuba area in Japan, estimated from a geological model and H/V spectral ratio of microtremors: 14th World Conference on Earthquake Engineering (14WCEE), Paper No. 07–0054.
Shapiro, N. M., and Campillo, M., 2004, Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise: Geophysical Research Letters, 31, L07614
| Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise:Crossref | GoogleScholarGoogle Scholar |
Stephans, M. A., 1970, Use of the Kolmogorov-Smirnov, Cramér-von Mises and related statistics without extensive tables: Journal of the Royal Statistical Society, Series B, 32, 115–122
Strollo, A., Richwalski, S. M., Parolai, S., Gallipoli, M. R., Mucciarelli, M., and Caputo, R., 2007, Site effects of the 2002 Molise earthquake, Italy: analysis of strong motion, ambient noise, and synthetic data from 2D modeling in San Giuliano di Puglia: Bulletin of Earthquake Engineering, 5, 347–362
| Site effects of the 2002 Molise earthquake, Italy: analysis of strong motion, ambient noise, and synthetic data from 2D modeling in San Giuliano di Puglia:Crossref | GoogleScholarGoogle Scholar |
Tazime, K., 1957, Minimum group velocity, maximum amplitude and quarter wavelength law. Love waves in doubly stratified layers: Journal of Physics of the Earth, 5, 43–50
| Minimum group velocity, maximum amplitude and quarter wavelength law. Love waves in doubly stratified layers:Crossref | GoogleScholarGoogle Scholar |
Tokimatsu, K., Tamura, S., and Kojima, H., 1992, Effects of multiple modes on Rayleigh wave dispersion characteristics: Journal of Geotechnical Engineering, 118, 1529–1543
| Effects of multiple modes on Rayleigh wave dispersion characteristics:Crossref | GoogleScholarGoogle Scholar |
Yokoi, T., 2010, New formulas derived from seismic interferometry to simulate phase velocity estimates from correlation methods using microtremor: Geophysics, 75, SA71–SA83
| New formulas derived from seismic interferometry to simulate phase velocity estimates from correlation methods using microtremor:Crossref | GoogleScholarGoogle Scholar |