ALBERT

Description: SUMMARY The horizontal-to-vertical spectral ratio (HVSR) of ambient noise is commonly used to infer a site's resonance frequency (${f_{0,site}}$). HVSR calculations are performed most commonly using the Fourier amplitude spectrum obtained from a single merged horizontal component (e.g. the geometric mean component) from a three-component sensor. However, the use of a single merged horizontal component implicitly relies on the assumptions of azimuthally isotropic seismic noise and 1-D surface and subsurface conditions. These assumptions may not be justified at many sites, leading to azimuthal variability in HVSR measurements that cannot be accounted for using a single merged component. This paper proposes a new statistical method to account for azimuthal variability in the peak frequency of HVSR curves (${f_{0,HVSR}}$). The method uses rotated horizontal components at evenly distributed azimuthal intervals to investigate and quantify azimuthal variability. To ensure unbiased statistics for ${f_{0,HVSR}}$ are obtained, a frequency-domain window-rejection algorithm is applied at each azimuth to automatically remove contaminated time windows in which the ${f_{0,HVSR}}$ values are statistical outliers relative to those obtained from the majority of windows at that azimuth. Then, a weighting scheme is used to account for different numbers of accepted time windows at each azimuth. The new method is applied to a data set of 114 HVSR measurements with significant azimuthal variability in ${f_{0,HVSR}}$, and is shown to reliably account for this variability. The methodology is also extended to the estimation of a complete lognormal-median HVSR curve that accounts for azimuthal variability. To encourage the adoption of this statistical approach to accounting for azimuthal variability in single-station HVSR measurements, the methods presented in this paper have been incorporated into hvsrpy, an open-source Python package for HVSR processing.

Description: The horizontal-to-vertical spectral ratio (HVSR) of ambient noise measurement is commonly used to estimate a site's resonance frequency (${f_0}$). For sites with a strong impedance contrast, the HVSR peak frequency (${f_{0,mathrm{ HVSR}}}$) has been shown to be a good estimate of ${f_0}$. However, the random nature of ambient noise (both in time and space), in conjunction with variable environmental conditions and sensor coupling issues, can lead to uncertainty in ${f_{0,mathrm{ HVSR}}}$ estimates. Hence, it is important to report ${f_{0,mathrm{ HVSR}}}$ in a statistical manner (e.g. as a mean or median value with standard deviation). In this paper, we first discuss widely accepted procedures to process HVSR data and estimate the variance in ${f_{0,mathrm{ HVSR}}}$. Then, we propose modifications to improve these procedures in two specific ways. First, we propose using a lognormal distribution to describe ${f_{0,mathrm{ HVSR}}}$ rather than the more commonly used normal distribution. The use of a lognormal distribution for ${f_{0,mathrm{ HVSR}}}$ has several advantages, including consistency with earthquake ground motion processing and allowing for a seamless transition between HVSR statistics in terms of both frequency and its reciprocal, period. Second, we introduce a new frequency-domain window-rejection algorithm to decrease variance and enhance data quality. Finally, we use examples of 114 high-variance HVSR measurements and 77 low-variance HVSR measurements collected at two case study sites to demonstrate the effectiveness of the new rejection algorithm and the proposed statistical approach. To encourage their adoption, and promote standardization, the rejection algorithm and lognormal statistics presented in this paper have been incorporated into hvsrpy, an open-source Python package for HVSR processing.