Ambient Vibration Array Research Articles

Characterization of Shallow Sedimentary Layers in the Oran Region Using Ambient Vibration Data

This study investigates the structure of shear-wave velocities (Vs) in the shallow layers of the Oran region, north-west of Algeria, using non-invasive techniques based on ambient vibration arrays. The region has experienced several moderate earthquakes, including the historical Oran earthquake of 1790. Ambient vibration measurements were carried out at 15 sites throughout the study area. Two methods were used: spatial autocorrelation (SPAC) and frequency–wavenumber analysis (f-k), which allowed us to better constrain Rayleigh wave dispersion curves. The inversion of the dispersion curves derived from the f-k analysis allowed for estimating the shear-wave velocity profiles and the Vs30 value at the sites under study. The other important result of the present study is an empirical equation that has been proposed to predict Vs30 in the Oran region. The determination of near-surface shear-wave velocity profiles is an important step in the assessment of seismic hazard. This study has demonstrated the effectiveness of using ambient vibration array techniques to estimate the soil Vs structure.

Coevolution of Rock Slope Instability Damage and Resonance Frequencies From Distinct‐Element Modeling

AbstractAmbient vibration measurements can detect resonance frequency changes related to rock slope instability damage or boundary condition changes during progressive failure. However, the impact of slope kinematics on resonance changes and the expected form and sensitivity of frequency evolution during destabilization require clarification to improve the implementation of this technique across diverse settings. Since instrumented rock slope failures are rare, numerical modeling is needed to study the anticipated spectral response from in situ monitoring. We used 2D distinct‐element modeling to evaluate the sensitivity and evolution of rock slope resonance behavior for slab toppling, flexural toppling, and planar sliding instabilities during progressive failure. Model simulations revealed that fundamental resonance frequency decreases between 20% and 60% with changes correlated with increasing length of open joints. Changes to higher‐order frequencies associated with landslide sub‐volumes were also detectable for cases with multiple fracture networks. Resonance behavior was most pronounced for failures dominated by steeply dipping open tension cracks, that is, flexural and slab toppling. Additionally, amplification patterns across the slope varied for the flexural toppling and sliding cases, providing potential new information with which to characterize landslide failure mechanisms using ambient vibration array measurements. Our results demonstrate landslide characteristics well‐suited for in situ ambient resonance monitoring and provide new data describing the anticipated changes in resonance frequencies during progressive rock slope failure.

Seismo-Stratigraphic Model for the Urban Area of Milan (Italy) by Ambient-Vibration Monitoring and Implications for Seismic Site Effects Assessment

In this paper, we present the work carried out to characterize the spatial variability of seismic site response related to local soil conditions in the city of Milan and its surroundings, an area with ∼3 million inhabitants and a high density of industrial facilities. The area is located at the northwestern end of the Po Plain, a large and deep sedimentary basin in northern Italy. An urban-scale seismo-stratigraphic model is developed based on new passive and active seismic data, supported by the available geological data and stratigraphic information from shallow and deep vertical wells. In particular, 33 single-station and 4 ambient-vibration array measurements are acquired, together with 4 active multichannel analyses of surface waves (MAWS). To estimate the resonant frequencies of the sediments, the horizontal to vertical spectral ratio technique (HVSR) is applied to the ambient-vibration recordings, whereas to determine the Rayleigh-wave dispersion curves from the passive array, the data are analysed using the conventional frequency-wavenumber, the modified spatial autocorrelation and the extended spatial autocorrelation (ESAC) techniques. The array data are used to determine the local shear wave velocity profiles, VS, via joint inversion of the Rayleigh-wave dispersion and ellipticity curves deduced from the HVSR. The results from HVSR show three main bands of amplified frequencies, the first in the range 0.17–0.23 Hz, the second from 0.45 to 0.65 Hz and the third from 3 to 8 Hz. A decreasing trend of the main peaks is observed from the northern to the southern part of the city, allowing us to hypothesize a progressive deepening of the relative regional chrono-stratigraphic unconformities. The passive ambient noise array and MASW highlight the dispersion of the fundamental mode of the Rayleigh-wave in the range 0.4–30 Hz, enabling to obtain detailed Vs. profiles with depth down to about 1.8 km. The seismo-stratigraphic model is used as input for 1D numerical modelling assuming linear soil conditions. The theoretical 1D transfer functions are compared to the HVSR curves evaluated from both ambient noise signals and earthquake waveforms recorded by the IV. MILN station in the last 10 years.

Seismic microzonation in a complex volcano-tectonic setting: the case of northern and western Ischia Island (southern Italy)

On August 21, 2017, the volcanic island of Ischia was struck by a Mw 3.9 (http://cnt.rm.ingv.it/event/16796811) earthquake which caused two victims and heavy damages on wide portions of the Casamicciola Terme, Lacco Ameno and Forio municipalities.Following the Decree 1/2018 of the Government Commission for the reconstruction of damaged areas of Ischia Island, a Grade 3 Seismic Microzonation study was committed by each municipality to private consultants, who were scientifically supported by a multidisciplinary team of geologists, geophysicists and geotechnical engineers from CNR, INGV and University of Naples “Federico II”.The existing data from previous geological and geotechnical investigations were integrated with the results of new geo-lithological surveys and geophysical measurements (HVSR, MASW, Down hole tests) carried out by the consultants. The whole amount of data was quality checked, homogenized and collected in a SQLite/SpatiaLite geodatabase using a QGIS open source software by the support team.A geo-lithological map of northern and western sides of the island was then produced, and the results of routine geophysical tests were coupled with those of seismic 2D ambient vibration arrays, to obtain shear wave velocity profiles (several hundreds of meters deep) down to the seismic bedrock, made of lithic tuffs and lavas, and to define its buried morphology. This in turn enabled to draw a map of seismically homogeneous microzones of the entire study area. Cyclic torsional shear tests were carried out on 10 undisturbed soil samples, mainly taken from epiclastic cover deposits, to characterize their non-linear behavior to be adopted for seismic response analyses.These latter were carried out along 1D profiles characterizing each microzone and along selected 2D cross sections running from the Mount Epomeo to the coast. The dynamic analyses allowed to detect and map areas characterized by variable amplification factors, of spectral accelerations throughout different ranges of periods (0.1- 0.5s, 04.-0.8s, 0.7-1.1 s). The zones most prone to site amplification are those corresponding to localized graben-like depressions at the toe of Mount Epomeo, such as Maio neighborhood of Casamicciola, where thickened epiclastic covers are entrenched into blocks of uplifted bedrock and where a combination of 1D stratigraphic amplification and 2D reflection and diffraction effects is expected.

Site characterization in Ismailia, Egypt using seismic ambient vibration array

Seismic ambient vibration array measurements were carried out at seven sites in Ismailia, Egypt in order to obtain one-dimensional (1D) shear wave velocity (Vs) profiles and horizontal-to-vertical spectral ratios (HVSR). We applied the spatial autocorrelation (SPAC) method to the ambient vibration data for the estimating Rayleigh-wave dispersion curves and obtained 1D Vs profiles. The deeper parts (90–350 m) of the estimated models have Vs values ranging between 500 and 700 m/s, while the subsurface parts have Vs varying between 140 and 200 m/s. The averaged Vs in the upper 30 m based on the results indicate that site classes of all the sites are classified as C or D according to NEHRP soil classification. The computed site amplification factors derived from the estimated Vs profiles show that the resonant frequencies at sites on the east side of the Suez Canal are below 1 Hz, while those at sites on the west side vary between 1 and 2 Hz, indicating the difference of subsurface geology across the Suez Canal. The HVSR curves at each site are generally similar for all sensors, indicating nearly flat subsurface layers beneath each site. We also compared the dominant frequencies of HVSRs and those derived from 1D amplification factors. The comparisons show that both the frequencies are similar to each other, but there are slight differences possibly due to the lack of deeper sedimentary layer structures and lower resolution of HVSRs in the lower frequency range, indicating the importance of further ambient vibration survey with broadband sensors. The results show that the combination of SPAC and HVSR methods is advantageous to estimate proper Vs profiles to evaluate site effects for earthquake disaster mitigation.

Estimating near surface shear wave velocity using the SPAC method at a site exhibiting low to high impedance contrast

In the general framework of seismic hazard analysis for quantifying site amplification, knowledge of the near surface shear wave velocity profile is crucial. The non-destructive and affordable ambient vibration array technique has been shown to be a superior method for estimating this shear wave velocity. The present study seeks to investigate the robustness of shear wave velocity profiles obtained using the ambient vibration array technique at sites exhibiting various (low to high) impedance contrasts by investigating a site in the central business district of Adelaide, South Australia. The city is founded on a regolith, which generally exhibits a range of impedance contrasts, and it is accepted that it is one of the Australian capital cities with the highest seismic risk. Eight ambient vibration array measurements were undertaken across the city. The acquired data were used to invert the shear wave velocity profiles, which are currently very limited, through the analysis of the dispersion characteristics of the surface waves by means of the spatial autocorrelation (SPAC) method, from which a further adjustment was carried out to suit the measured horizontal vertical spectral ratio (HVSR) curve ellipticity. Validation of the proposed shear wave velocity models was carried using previous studies and forward modeling techniques. The present study demonstrates that the SPAC method provides reliable results at regolith sites which are subject to various (low to high) impedance contrasts, provided the frequency band of the captured dispersion curve is wide enough to resolve the investigated depth. The primary innovations of this paper are to: (1) investigate the applicability of the microtremor SPAC method to obtain the near surface shear wave velocity of a low to high impedance contrast site; and (2) propose a new near surface shear wave velocity profile for Adelaide, South Australia.

Rayleigh wave three-component beamforming: signed ellipticity assessment from high-resolution frequency-wavenumber processing of ambient vibration arrays

Rayleigh wave three-component beamforming: signed ellipticity assessment from high-resolution frequency-wavenumber processing of ambient vibration arrays

Characterization of site conditions (soil class, VS30, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods

Data provided by accelerometric networks are important for seismic hazard assessment. The correct use of accelerometric signals is conditioned by the station site metadata quality (i.e., soil class, VS30, velocity profiles, and other relevant information that can help to quantify site effects). In France, the permanent accelerometric network consists of about 150 stations. Thirty-three of these stations in the southern half of France have been characterized, using surface-wave-based methods that allow derivation of velocity profiles from dispersion curves of surface waves. The computation of dispersion curves and their subsequent inversion in terms of shear-wave velocity profiles has allowed estimation of VS30 values and designation of soil classes, which include the corresponding uncertainties. From a methodological point of view, this survey leads to the following recommendations: (1) perform both active (multi-analysis surface waves) and passive (ambient vibration arrays) measurements to derive dispersion curves in a broadband frequency range; (2) perform active acquisitions for both vertical (Rayleigh wave) and horizontal (Love wave) polarities. Even when the logistic contexts are sometimes difficult, the use of surface-wave-based methods is suitable for station-site characterization, even on rock sites. In comparison with previous studies that have mainly estimated VS30 indirectly, the new values here are globally lower, but the EC8-A class sites remain numerous. However, even on rock sites, high frequency amplifications may affect accelerometric records, due to the shallow relatively softer layers.

Fieldwork on local-site seismic response in the Po Plain: examples from ambient vibration array and single station analyses

In the framework of the S2 DPC-INGV 2014–2015 project (https://sites.google.com/site/ingvdpc2014progettos2/), several field surveys were designed with the aim to investigate the local seismic site response of selected Italian Accelerometric Network strong motion stations. In this study the results obtained at Sorbolo, Lesignano de’ Bagni, and Castelfranco Veneto, three sites located in the Po alluvial basin (northern Italy), were presented. In particular, at each site, ambient vibration surveys were performed both in single stations and two dimensional arrays configurations. The ambient vibration arrays were analyzed to obtain at each site the shear-wave velocity (Vs) profile, until depths in correspondence of which the Vs show values higher than 800 m/s (i.e. lower bound for the soil category A, as now defined in the Italian seismic code for building). In this way, data were processed through a joint inversion computed using both the Rayleigh waves dispersion curve and the ellipticity curve deduced from the horizontal-to-vertical spectral ratios calculated from the microtremor single station measurements. In general, the obtained profiles show an increasing depth of the main Vs discontinuity (i.e. the seismic bedrock) from 60 to 320 m, moving from the edge to the central part of the Po alluvial basin. To support the outcomes of the field experiments, the available earthquakes recorded at each station in the last decades were processed in terms of the horizontal-to-vertical spectral ratios. Although in many cases the resonance frequencies are comparable with those obtained from microtremor, a high variability is observable depending on the selected portions of each single earthquake. Finally, residual analysis were carried out comparing the recorded ground motion to predicted data as provided by the actual Italian ground motion prediction equations, calibrated both at national and at regional scale.

Evaluating the structural model from ambient vibration recordings and shallow seismic refraction techniques. A case study of Zagazig City, Egypt

This study uses the seismic refraction and noise measurements to investigate the velocity structure of the subsurface and emphasize the advantage of ambient vibration over the conventional seismic refraction technique. Field measurements were carried out at nine sites in and around Zagazig city. Shallow seismic refraction data were interpreted using the delay time method to obtain the two-dimension ground model at each site. Ambient vibration arrays are used to infer the one-dimensional compressional and shear wave velocity profiles through two main steps. The first step is to derive the dispersion curve from the recorded signals using the frequency-wavenumber method. The second is to invert the dispersion curve to obtain the site velocity profiles. The results of the compressional wave velocities obtained from seismic refraction technique showed that the subsurface consists of a number of layers ranging from two to four layers and give a good agreement with the results of the seismic wave velocities obtained from the ambient vibration arrays. The ambient vibration arrays gave a deeper depth of penetration than the other method, providing more information on the subsurface structure without any disturbance to the environment. This work provides reliable estimates of the seismic velocity structures of both shallow and deep sedimentary layers within the area of interest.

Shear-wave velocity profile and seismic input derived from ambient vibration array measurements: the case study of downtown L'Aquila

Shear-wave velocity profile and seismic input derived from ambient vibration array measurements: the case study of downtown L'Aquila

Modal analysis of 2-D sedimentary basin from frequency domain decomposition of ambient vibration array recordings

Frequency domain decomposition (FDD) is a well-established spectral technique used in civil engineering to analyse and monitor the modal response of buildings and structures. The method is based on singular value decomposition of the cross-power spectral density matrix from simultaneous array recordings of ambient vibrations. This method is advantageous to retrieve not only the resonance frequencies of the investigated structure, but also the corresponding modal shapes without the need for an absolute reference. This is an important piece of information, which can be used to validate the consistency of numerical models and analytical solutions. We apply this approach using advanced signal processing to evaluate the resonance characteristics of 2-D Alpine sedimentary valleys. In this study, we present the results obtained at Martigny, in the Rhône valley (Switzerland). For the analysis, we use 2 hr of ambient vibration recordings from a linear seismic array deployed perpendicularly to the valley axis. Only the horizontal-axial direction (SH) of the ground motion is considered. Using the FDD method, six separate resonant frequencies are retrieved together with their corresponding modal shapes. We compare the mode shapes with results from classical standard spectral ratios and numerical simulations of ambient vibration recordings.

Love wave contribution to the ambient vibration H/V amplitude peak observed with array measurements

This study applies array methods to measure the relative proportions of Love and Rayleigh waves in the ambient vibration wavefield. Information on these properties is of special relevance for frequencies around the horizontal-to-vertical (H/V) spectral amplitude ratio peak. The analysis of H/V curves, a popular technique in site characterisation, commonly assumes that the curves represent the frequency-dependent Rayleigh wave ellipticity. For the detailed interpretation of amplitudes or the inversion of the curves, it is therefore necessary to estimate and correct for the contribution of other wave types to the ambient vibration wavefield. I use available ambient vibration array measurements to determine the relative amount of Love and Rayleigh waves on the horizontal components by frequency-dependent analysis of the main propagation and polarisation directions, with a special emphasis on the H/V peak frequency as determined from the same recordings. Tests with synthetic data demonstrate the feasibility of this approach, at least in the presence of dominant source regions. Analysis of the data from 12 measurements at nine European sites, which include shallow as well as deep locations that span a wide range of impedance contrasts at the sediment-bedrock interface, indicates that the relative contribution of Rayleigh waves varies widely with frequency, from close to 0% to more than 70%. While most data sets show relative Rayleigh wave contributions between 40% and 50% around the H/V peak, there are also examples where Love waves clearly dominate the wavefield at the H/V peak, even for a site with a low impedance contrast. Longer-term measurements at one site indicate temporal variations in the relative Rayleigh wave content between day- and nighttime. Results calculated with the method introduced herein generally compare well with results of modified spatial autocorrelation analysis. These two methods might be used in a complimentary fashion, as both rely on different properties of the ambient vibration wavefield. This study illustrates that it is possible to measure the relative Rayleigh wave content of the noise wavefield from array data. Furthermore, the examples presented herein indicate it is important to estimate this property, as the assumption that there are an equal proportion of Love and Rayleigh waves is not always correct.

On the repeatability and consistency of three-component ambient vibration array measurements

Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves—as observed in other studies—are consequently not adequately explained by uncertainties in the measurement part.

Analysis of three-component ambient vibration array measurements

SUMMARY Both synthetic and observed ambient vibration array data are analysed using high-resolution beam-forming. In addition to a classical analysis of the vertical component, this paper presents results derived from processing horizontal components. We analyse phase velocities of fundamental and higher mode Rayleigh and Love waves, and particle motions (ellipticity) retrieved from H/V spectral ratios. A combined inversion with a genetic algorithm and a strategy for selecting possible model parameters allow us to define structural models explaining the data. The results from synthetic data for simple models with one or two layers of sediments suggest that, in most cases, the number of layers has to be reduced to a few sediment strata to find the original structure. Generally, reducing the number of soft-sediment layers in the inversion process with genetic algorithms leads to a class of models that are less smooth. They have a stronger impedance contrast between sediments and bedrock. Combining Love and Rayleigh wave dispersion curves with the ellipticity of the fundamental mode Rayleigh waves has some advantages. Scatter is reduced when compared to using structural models obtained only from Rayleigh wave phase velocity curves. By adding information from Love waves some structures can be excluded. Another possibility for constraining inversion results is to include supplementary geological or borehole information. Analysing radial components also can provide segments of Rayleigh wave dispersion curves for modes not seen on the vertical component. Finally, using ellipticity information allows us to confine the total depth of the soft sediments. For real sites, considerable variability in the measured phase velocity curves is observed. This comes from lateral changes in the structure or seismic sources within the array. Constraining the inversion by combining Love and Rayleigh wave information can help reduce such problems. Frequency bands in which the Rayleigh wave dispersion curves show considerable scatter are often better resolved by Love waves. Information from the horizontal component can be used to correctly assign the mode number to the different phase–velocity curve segments, especially when two modes seem to merge at osculation points. Such merging of modes is usually observed for Rayleigh waves and thus can be partly solved if additional information from the Love waves and the horizontal component of Rayleigh waves is considered. Whenever a site presents a velocity inversion below the top layer, Love wave data clearly helps to better constrain the solution.

S-wave velocity measurements applied to the seismic microzonation of Basel, Upper Rhine Graben

An extensive S-wave velocity survey had been carried out in the frame of a recent seismic microzonation study of Basel and the border areas between Switzerland, France and Germany. The aim was to better constrain the seismic amplification potential of the surface layers. The survey included single station (H/V spectral ratios) and ambient vibration array measurements carried out by the Swiss team, as well as active S-wave velocity measurements performed by the German and French partners. This paper is focused on the application of the array technique, which consists in recording ambient vibrations with a number of seismological stations. Several practical aspects related to the field measurements are outlined. The signal processing aims to determine the dispersion curves of surface waves contained in the ambient vibrations. The inversion of the dispersion curve provides a 1-D S-wave velocity model for the investigated site down to a depth related to the size of the array. Since the size of arrays is theoretically not limited, arrays are known to be well adapted for investigations in deep sediment basins, such as the Upper Rhine Graben including the area of the city of Basel. In this region, 27 array measurements with varying station configurations have been carried out to determine the S-wave velocity properties of the geological layers down to a depth of 100–250 m. For eight sites, the outputs of the array measurements have been compared with the results of the other investigations using active sources, the spectral analysis of surface waves (SASW) and S-wave reflection seismics. Borehole information available for a few sites could be used to calibrate the geophysical measurements. By this comparison, the advantages and disadvantages of the array method and the other techniques are outlined with regard to the effectiveness of the methods and the required investigation depth. The dispersion curves measured with the arrays and the SASW technique were also combined and inverted simultaneously to use the advantages of both methods. Finally, the paper outlines and discusses the contribution of the S-wave velocity survey to the new seismic microzonation of the Basel region. In this regard one major outcome of the survey is the quantification of vertical and lateral changes of the S-wave velocity, due to changing lithology or changing compaction and degree of weathering of the layers.