Exploring ocean floor seismic structure in the Guerrero seismic gap: insights from site response, attenuation, and velocity structure

The horizontal-to-vertical spectral ratio (HVSR) has been extensively used in site characterization utilizing recordings from microtremor and earthquake in recent years. This method is proposed based on ground pulsation, and then it has been applied to both S-wave and ambient noise, accordingly, in practical application also different. The main applications of HVSR are site classification, site effect study, mineral exploration, and acquisition of underground average shear-wave velocity structure. In site response estimates, the use of microtremors has been introduced long ago in Japan, while it has long been very controversial in this research area, as there are several studies reporting difficulties in recognizing the source effects from the pure site effects in noise recordings, as well as discrepancies between noise and earthquake recordings. In practice, the most reliable way is the borehole data, and the theoretical site response results were compared with the HVSR using shear wave to describe site response. This paper summarizes the applications of the HVSR method and draws conclusions that HVSR has been well applied in many fields at present, and it is expected to have a wider application in more fields according to its advantages.

ABSTRACTThe Longmen Shan fault zone that was shocked by the 12 May 2008 M 8.0 Wenchuan earthquake acts as the boundary between the western edge of the Sichuan basin and the steep eastern margin of the Songpan-Ganze block. In this study, continuous seismic data recorded by 176 temporary short-period seismic stations between 22 October and 20 November 2017 are used to study the shallow crustal structure of the Longmen Shan fault zone by applying ambient-noise tomography and horizontal-to-vertical spectral ratio (HVSR) analysis. From ambient-noise analysis, fundamental-mode Rayleigh-wave dispersion curves between 0.25 and 1 Hz are extracted. Then, the direct surface-wave tomographic method is used to invert surface-wave dispersion data for the 3D shallow shear-wave velocity structure. Our results show that low shear-wave velocities are mainly distributed around the surface rupture trace of the Wenchuan earthquake at least down to 2 km. From the HVSR method, the sites are sorted into two types according to the pattern of HVSR curves with single peak or double peak. By converting frequency to depth, the results show that the sediments are thicker near the surface rupture. The low-velocity zone based on ambient-noise tomography agrees well with the distribution of sedimentary cover estimated from HVSR, which are generally consistent with geological information. Our results provide high-resolution shallow crustal velocity structure for future detailed studies of the Longmen Shan fault.

To overcome the limited information in the northern waters of Sulawesi Island, ocean bottom seismometers (OBS) were installed in the northern part of the Makassar Strait and the Sulawesi Sea. OBS data were analysed using passive seismic methods to determine the characteristics of noise and site response in the Makassar Strait and the Celebes Sea. Probabilistic Power Spectral Density (PPSD) from OBS station showed microtremor in the period 0.06 – 0.7 s, secondary microseismic peak in the period 0.9 – 2 s, secondary microseismic peak in the period 4 – 8 s, and hum in the period more than 10 s. The amplitude of the noise power is in the range of New Low Noise Model (NLNM) and New High Noise Model (NHNM) for a period <10 s and slightly above NHNM for a period >10 s for OBS IGGCAS F-type. For OBS IGGCAS G-type, the noise power amplitude is in the period < 0.7 s in the NLNM and NHNM ranges. However, there is a tendency for the noise power amplitude to be flat with an amplitude above NHNM for a period > 0.7 s. The noise power’s amplitude depends on the station’s location, climate, instruments, daily variations, and seasonal variations. Horizontal-to-Vertical Spectral Ratio (HVSR) in the North Makassar Strait shows an H/V curve with an unclear peak in the middle of the strait and multiple peaks in the eastern part of the strait with a frequency range of 1.35 – 4.7 Hz. Between the Makassar Strait and the Celebes Sea, the H/V curve shows a flat peak and multiple peaks with a frequency range of 1.10 – 4.03 Hz. Meanwhile, in the Celebes Sea, the H/V curve shows an unclear peak (unclear peak) and several peaks (multiple peaks) with a frequency range of 0.65 – 6.25 Hz.

ABSTRACTDamaging ground motions from the 2011 Mw 5.8 Virginia earthquake were likely increased due to site amplification from the unconsolidated sediments of the Atlantic Coastal Plain (ACP), highlighting the need to understand site response on these widespread strata along the coastal regions of the eastern United States. The horizontal-to-vertical spectral ratio (HVSR) method, using either earthquake signals or ambient noise as input, offers an appealing method for measuring site response on laterally extensive sediments, because it requires a single seismometer rather than requiring a nearby bedrock site to compute a horizontal sediment-to-bedrock spectral ratio (SBSR). Although previous studies show mixed results when comparing the two methods, the majority of these studies investigated site responses in confined sedimentary basins that can generate substantial 3D effects or have relatively small reflection coefficients at their base. In contrast, the flat-lying ACP strata and the underlying bedrock reflector should cause 1D resonance effects to dominate site response, with amplification of the fundamental resonance peaks controlled by the strong impedance contrast between the base of the sediments and the underlying bedrock. We compare site-response estimates on the ACP strata derived using the HVSR and SBSR methods from teleseismic signals recorded by regional arrays and observe a close match in the frequencies of the fundamental resonance peak (f0) determined by both methods. We find that correcting the HVSR amplitude using source term information from a bedrock site and multiplying the peak by a factor of 1.2 results in amplitude peaks that, on average, match SBSR results within a factor of 2. We therefore conclude that the HVSR method may successfully estimate regional linear weak-motion site-response amplifications from the ACP, or similar geologic environments, when appropriate region-specific corrections to the amplitude ratios are used.

Ambient-noise methods, and the horizontal-to-vertical spectral ratio (HVSR), whether from earthquakes (eHVSR) or mircrotremors (mHVSR), have gained much popularity in the field of site response over the last decade. These methods can be used either for direct prediction or to inform the velocity structure utilized in more conventional site response analyses. This paper describes a field study and subsequent analyses undertaken in the Lower Hutt sedimentary basin of New Zealand. The field study involved collecting 50 ambient-noise or microtremor measurements across the entire basin over the same time window that microtremor measurements were being recorded at strong-motion stations in the basin for use as reference stations. Additionally, microtremor array measurements (MAM; involving ~24 ambient-noise measurements per site) and multi-channel analysis of surface waves (MASW) were conducted at five sites to better quantify deep and shallow velocity structure of the basin. In total, microtremor data were collected at 154 locations. This paper focuses on the use of the microtremor data collected for direct prediction of site response in the sedimentary basin. The hybrid standard spectral ratio approach is tested in this region. A rigorous validation study was performed at strong motion stations at which microtremor measurements were made for use as reference basin sites in the hybrid spectral ratio method. The prediction accuracy and uncertainty of the method, when using synchronized versus unsynchronized data between the reference basin sites and temporary target sites, are compared. The hybrid spectral ratio method predicts well the observed site amplification between nearby, deep basin sites for f < 5 Hz, when synchronized data are used. Predictions around the fundamental frequency of the basin (corresponding to periods of 1- 2 seconds) worsen when unsynchronized data are used due the influence from environmental and anthropogenic factors on microtremor amplitudes. Finally, the method is used to predict the site response, relative to a rock reference site, at all basin sites at which synchronized temporary microtremor data were collected. Early efforts to spatially interpolate the observations and predictions are discussed.

Pusan National University deployed sixteen ocean bottom seismometers (OBSs) in the eastern offshore of the southern Korean peninsula. The primary purpose of the OBS network is monitoring earthquakes in the eastern offshore to investigate potential fault systems in the offshore region which was formidable using limited apertures by land-based observations. A seismic network's performance highly depends on each site's background noise level. We analyze the nature of the ambient noise and site response for ocean bottom seismometers (OBSs) deployed in the 2021-2022 period. The power spectral densities (PSDs) of the OBSs exhibited dis-similar features from those of land-based stations; most temporary broadband seismic stations on land showed relatively lower background noise levels. In the meanwhile, OBSs showed higher background noise levels. We report the nature of ambient noise at various channels, water depths, spatial locations, temporal variations, extreme weather conditions, and their potential causes. In general, horizontal components are noisier than vertical components. For longer periods, horizontal components are larger by ~45 dB. The probability density functions (PDFs) of OBSs show that the noise level is within the range of McNamara&amp;#8217;s model (2004) for higher frequencies (3.5~50 Hz) although they are still high. When examining long periods (&gt; 20 s), the noise level is higher than what would be given by McNamara&amp;#8217;s model. Although we do not observe diurnal or weekly variations in OBS, as expected, we observe varying degrees of seasonal variations in OBSs. Apparently, water depth is the most important factor in deciding noise levels and their seasonal variations. At shallow-depth OBSs, we observe a strong correlation between noise levels and wave heights estimated by the Korea Meteorological Administration (KMA). The ambient noise is down to -130 dB for the band from 5 to 15 Hz, which provides the best signal-to-noise ratio for local microearthquakes. We also present the horizontal-to-vertical spectral ratio (HVSR) of the ambient noise recorded by OBSs. They present significant amplifications at lower frequencies, which indicates large amplification by combined effects due to lower density and lower wave velocity at shallow sediments, and greater depths to major impedance contrast. We confirm that modeling HVSR of noise data recorded by a three-component OBS offers a fast and inexpensive method for site investigation in deep water with the potential of in situ seafloor sediment characterization.

We use time-dependent horizontal-to-vertical spectral ratios (HVSR) of microtremors to determine the dominant frequencies of vibration of the geological structures beneath several recording sites in the vicinity of Teide volcano (Canary Islands, Spain). In the microtremors, the time-dependent HVSRs (ratiograms) are a useful tool to discriminate between the presence of real dominant frequencies linked to resonances of the subsurface structure and the spurious appearance of peaks due to local transients. We verified that the results are repeatable, in the sense that microtremors recorded at the same site but at different times yield a very similar HVSR function. Two types of results are found: (1) sites where there is no resonance of the propagating microtremors, and therefore no value of a dominant frequency can be assessed; and (2) sites where a stationary peak in the HVSR is found and a dominant frequency related to resonance of the shallow structure can be estimated. These resonant frequencies show substantial spatial variations even for nearby sites, which reflects the complexity of the shallow velocity structure in the Las Canadas area. Large dominant frequencies occur near the caldera walls and also at a few locations that coincide with the intersections of the inferred rims of the three calderas forming Las Canadas. Small dominant frequencies also occur near the caldera rim, and may be due to discontinuities in the caldera wall and/or to local velocity anomalies. Intermediate frequencies are mostly found in the eastern part of the caldera, where a tentative profile of the basement depth has been obtained. Intermediate frequencies have also been measured south of Ucanca and south of Montana Blanca. In view of the present results, we conclude that the use of ratiograms constitutes an improvement of the HVSR method and provides an appropriate tool to investigate the shallow velocity structure of a volcanic region.

A large number of earthquake studies using both empirical and theoretical approaches clearly depict the strong correlation of resulting damage and local geology. This correlation forms the scientific basis of the subsurface soil structure studies for site response evaluation, since adequate knowledge of the subsoil geotechnical and geophysical properties can lead to realistic seismic hazard estimation through appropriate modeling of strong seismic motion. In this framework, ambient noise analysis and one dimensional numerical simulation modeling have been performed for the town of Grevena (Northwestern Greece) in order to study both the subsurface soil structure, as well as its expected effect on seismic motions. The horizontal to vertical spectral ratio (HVSR) method was implemented on an almost uniform grid of 60 single station ambient noise measurements inside the urban area, while the noise array technique was applied at four selected sites of the study area. The HVSR curves show complex patterns, occasionally with double HVSR peak frequencies, with the higher one HVSR frequency showing a good correlation with the recent Holocene clay-dominant formation. Using these results, as well as information on surface geology and existing geotechnical data, a microzonation of the Grevena town is attempted, indicating zones of similar site response to ground motion. Numerical simulation of ambient noise with synthetic recordings suggests that the complex features observed in several HVSR curves could be attributed to the impedance contrast due to the quite shallow (Holocene clays to Pliocene–Pleistocene sands) and a deeper (Pliocene–Pleistocene sands to Oligocene–Miocene bedrock) sedimentary formation.

This study provides the assessment of site characterization and possible shallow shear-velocity structure from the study of the horizontal to vertical spectral ratio (HVSR) measurements using the ambient noise or microtremor (herein called classical HVSR), extracted Rayleigh wave from the ambient noise data (herein called standard HVSR) and earthquake (herein called earthquake HVSR) data at five stations in the Surat district of mainland Gujarat, India. These locations are the hub of many mining and industrial projects like oil and natural gas, which need to function safely within the seismic hazard and ground shaking limits. From the classical and standard HVSR datasets, estimates of the predominant resonant frequency of the soil are obtained, observed to be well matched, from which first order inversions are carried out around the predominant frequency to provide a fair estimate of thickness of the regimented layer over a hard seismic substratum up to a depth of 100 m. In the standard HVSR datasets, the Rayleigh wave ellipticity curves are extracted with time–frequency analysis using continuous wavelet transforms. This is followed by the Rayleigh wave ellipticity inversion approach to derive a first order approximate sedimentary shear velocity structure. We also compute HVSR measurements using earthquakes. The noise and earthquake HVSR curves are well-matched in terms of the predominant frequencies and range from 3.8 to 16.7 Hz and 3.2 to 16.5, respectively. The estimated VS30 values (the average shear wave velocity (VS) for the top 30 m of the soil) vary from 520 to 1350 m/s, matching well with some of the geotechnical studies conducted here. The study emphasizes the effectiveness of the single station HVSR method in determination of hitherto unknown soil structures as economical and non-invasive exploration viability and proving quite useful for critical centres of industrial establishments.

Local amplification and subsoil structure at a difficult site: Understanding site effects from different measurements

The effect of local site conditions on the earthquake ground motion is a very important factor to be considered in engineering seismic fortification. Many methods, such as numerical simulation methods based on site analysis models and statistical empirical relation methods based on the earthquake ground motion observations and numerical simulation data, have been used to consider the site effects in actual engineering seismic fortification and earthquake disaster assessment. The statistical analysis to obtain characteristic parameters of site condition effect based on strong motion and microtremor records become an economical and practical method of determining the designed ground motion of engineering sites, especially for large survey areas and engineering sites where it is difficult to carry out a site survey. In this paper, a novel evaluation method for site effect on earthquake ground motion is proposed. The new method is based on the horizontal to vertical spectral ratio (HVSR) method, but the original HVSR is replaced by a modified HVSR considering the effect of the soil layer on the vertical ground motion. In order to build the model and determine the corresponding parameters of the modified HVSR, first, the ground motions in the bedrock below the soil layer are calculated using the one-dimensional equivalent linear method. These calculated records are independent of the influence of the downgoing wavefield, and the differences between the ground surface to bedrock spectral ratio (SBSR). The HVSR for the local sites of ground motion observation stations are analyzed using the strong ground motion records from the Kiban-Kyoshin network (KiK-net) in Japan. The statistical characteristics of the relationship between SBSR and HVSR are revealed, and then, a quantitative relationship between SBSR/HVSR and HVSR is established. The proposed evaluation method for the site effect has the advantage that the original HVSR method only requires ground motion records on the ground surface of the site, and it further considers the influence of the vertical seismic effect on the accuracy of the HVSR method. The proposed method can characterize the influence of the site conditions on ground motion more reasonably than the conventional method.

Approximate analytical HVSR curve using multiple band-pass filters and potential applications

SUMMARY Central to the task of seismic hazard assessment is the evaluation of potential amplifications due to site effects. The horizontal-to-vertical spectral ratio (HVSR) of ambient seismic noise (ASN) is a widespread measurement to assess the predominant soil frequency of a given site and estimate the wave velocities of the subsoil stratigraphy using inversion schemes. In practice, the inversions are currently made, assuming flat layers. In fact, HVSR measurements may show significant lateral and azimuthal changes due to the spatial variations of local geology, which can introduce uncertainties into the characterization of a site. This suggests the importance of considering detailed measurements of lateral and angular variations in layered settings. Inversion of soil properties at depth for horizontally layered media has become feasible assuming that ASN constitutes a diffuse field, that is, produced by equi-partitioned uniform illumination and/or by random sources and the ensuing multiple scattering by heterogeneities. Under the diffuse field assumption (DFA), the HVSR have been modelled by calculating the imaginary parts of the Green's function (IMGs) when source and observer coincide at the same point. In this work we use the 3-D indirect boundary element method (IBEM) to model the HVSR for each independent horizontal direction referred to here as directional-HVSR for layered media with lateral inhomogeneity. The IMGs at the source required to get HVSR have locality properties that depend on frequency and may imply significant economies in the calculation. For simple models we modelled the IMGs approximately using an adaptive meshing scheme that accounts both for the locality of the problem and the diffraction properties of waves at low and high frequencies. The obtained directional-HVSR displays variations in both frequency and azimuth. The results also show that layer interface variations can lead to ‘spots’ of higher wave excitation associated to local resonant modes. This shows the importance of HVSR in forecasting earthquake response and suggests the need for denser field measurements to study lateral and azimuthal variability. In order to show the reality of directional-HVSR, field data from Chalco, a soft soil site at the southern part of the Valley of Mexico, have been preliminarily analysed.

We conducted experimental work to explain the large peak ground accelerations observed at the Cerro Prieto volcano in Mexicali Valley, Mexico. Using ambient noise and earthquake data, we compared horizontal-to-vertical spectral ratios (HVSRs) computed for sites on the volcano against those calculated for locations outside it. High-HVSR values (∼11 at ∼2 Hz) were obtained on the top of the volcano at 183 m of altitude, decreasing for sites located at lower elevations. We calculated a median HVSR of ∼1 at 2 Hz from HVSRs computed for nine sites located along an N18°E transect and at an average elevation of ∼25 m. The earlier comparison suggests a relative amplification on the volcano. In addition, we calculated HVSRs from accelerograms generated by 62 earthquakes (2.6≤ML≤5.4; 4.6≤Mw≤7.2) recorded at four locations: two on the volcano (at 194 and 110 m of elevation) and two outside it. These last two sites, located up to 6 km away in a north-northwest and south-southwest direction relative to the volcano, are at an average altitude of 22 m. For the four locations, we also computed the HVSRs from ambient noise data. Although the HVSR results derived from both types of data are slightly different, we also found high HVSRs for the two sites on the volcano and low HVSRs for the two sites outside it, corroborating the relative amplification on the volcano. Using the 1D wave propagation modeling, based on the stiffness matrix method, we modeled the experimental HVSRs to analyze the local site effects. Therefore, we propose that the ground-motion amplification at the Cerro Prieto volcano may be due to a combination of its topography and shallow site effects.

Abstract. Earthquake site response is an essential part of seismic hazard assessment, especially in densely populated areas. The shallow geology of the Netherlands consists of a very heterogeneous soft sediment cover, which has a strong effect on the amplitude of ground shaking. Even though the Netherlands is a low- to moderate-seismicity area, the seismic risk cannot be neglected, in particular, because shallow induced earthquakes occur. The aim of this study is to establish a nationwide site-response zonation by combining 3D lithostratigraphic models and earthquake and ambient vibration recordings. As a first step, we constrain the parameters (velocity contrast and shear-wave velocity) that are indicative of ground motion amplification in the Groningen area. For this, we compare ambient vibration and earthquake recordings using the horizontal-to-vertical spectral ratio (HVSR) method, borehole empirical transfer functions (ETFs), and amplification factors (AFs). This enables us to define an empirical relationship between the amplification measured from earthquakes by using the ETF and AF and the amplification estimated from ambient vibrations by using the HVSR. With this, we show that the HVSR can be used as a first proxy for site response. Subsequently, HVSR curves throughout the Netherlands are estimated. The HVSR amplitude characteristics largely coincide with the in situ lithostratigraphic sequences and the presence of a strong velocity contrast in the near surface. Next, sediment profiles representing the Dutch shallow subsurface are categorised into five classes, where each class represents a level of expected amplification. The mean amplification for each class, and its variability, is quantified using 66 sites with measured earthquake amplification (ETF and AF) and 115 sites with HVSR curves. The site-response (amplification) zonation map for the Netherlands is designed by transforming geological 3D grid cell models into the five classes, and an AF is assigned to most of the classes. This site-response assessment, presented on a nationwide scale, is important for a first identification of regions with increased seismic hazard potential, for example at locations with mining or geothermal energy activities.