Shear Wave Velocity Profiles Research Articles

Investigating small-strain site response using inverse soil dynamic parameters from downhole arrays

Previous research on site response has primarily focused on reproducing and predicting ground motions recorded by vertical seismic arrays. This study evaluated the ability to reproduce and predict small-strain site response at the Treasure Island Downhole Array and Delaney Park Downhole Array using inverse soil dynamic parameters. First, parameters such as seismic wave velocity profiles and Rayleigh damping ratios were obtained through seismic interferometry and spectral ratio methods based on the earthquake records. Subsequently, the validity and reliability of these inverse soil dynamic parameters were evaluated by comparing the simulated and observed ground responses. Finally, the inverse method was compared with other methods considering spatial variability in one-dimensional site response analysis. Quantitative comparisons show that the inverse method effectively predicts and reproduces site response and outperforms the traditional method relying on the single borehole shear wave velocity profile and damping ratio from empirical model. Moreover, the inverse method better captures spatial variability at the Delaney Park Downhole Array.

Shear wave velocity profiling of Riyadh City, Saudi Arabia, utilizing the multi-channel analysis of surface waves method

Abstract Geotechnical site characterization is very important for construction purposes. This study has been conducted in Diriyah area northwest of Riyadh City, Saudi Arabia, using the Multichannel Analysis of Surface Waves (MASW) method for site characterization through shear wave velocity profiling to 30 m depth. Nineteen MASW lines were carried out in various directions and lengths through the area. The entire process was meticulously parameterized to extract shear wave velocity for subsurface characteristics. MASW results revealed four distinct velocity zones based on National Earthquake Hazard Reduction Program. Fill material was approximately half a meter thick and was classified as very dense soil. The second layer exhibited velocities ranging from 800 to 1,500 m/s, indicating weathered and highly fractured limestone. The third layer showed velocities varying from 1,500 to 1,800 m/s, representing slightly weathered limestone. The fourth layer displayed high velocities ranging from 1,800 to 3,600 m/s, indicating hard and compact limestone rocks. Geotechnical boreholes were drilled down to depths of 10–35 m. These boreholes exposed the geological model that consisted of fill material (silty sand with gravel), followed by highly to moderately weathered limestone with vugs and cracks, and finally, massive limestone rock. Analysis of shear wave velocities identified weak zones, particularly fractured and weathered limestone rocks extending to 12 m in depth. Sinkholes of circular, elongated, and/or conical shapes were observed within this depth range. Moreover, some sinkholes were detected at depths greater than 12 m in specific locations (sites 1, 6, 9, 11, and 17). These sinkholes agreed with the previous study. These results highlight the need for targeted ground improvement methods, such as grouting or underpinning, particularly for construction over weaker zones. Accurate site classification and effective risk management are crucial for addressing these geotechnical and seismic challenges.

3D modeling of geological and geotechnical soil characterization using the MASW method: a case study in Southern Ijaw LGA, Bayelsa State, Nigeria

The Multichannel Analysis of Surface Waves (MASW) method is pivotal for non-invasive subsurface shear-wave velocity profiling, essential in geotechnical and seismic investigations. This study aims to model 3D geological and geotechnical soil properties in Southern Ijaw LGA, Bayelsa State, Nigeria, using MASW. Data collection involved a 12-channel ABEM Terraloc Mark 6 geophone system across nine survey points, complemented by nine boreholes and two Standard Penetration Tests (SPT). The results revealed a shallow water table at 0.30 m and stratified soil layers ranging from soft clay to dense sands. SPT N-values increased from 2 at the surface to 34 in deeper layers, reflecting significant soil densification and strength enhancement with depth. Shear-wave velocities (Vs) ranged from 207.11 m/s at 0.87 m to 502.88 m/s at 12.9 m, while compressional-wave velocities (Vp) ranged from 388.71 m/s to 948.98 m/s. Poisson’s ratio was consistent at 0.3, indicating uniform properties across depths. Elastic moduli, including shear modulus (µ), Young’s modulus (E), and bulk modulus (K), increased with depth, indicating greater subsurface material stiffness. The 3D geological model delineated distinct layers: soft clay (0–3.61 m), silty clay (3.61–5.19 m), fine sand (5.19–6.55 m), medium sand (6.55–7.83 m), and medium to coarse sand (7.83–12.90 m). These layers exhibited specific engineering properties, with ultimate and allowable bearing capacities ranging from 154.38 to 543.71 kPa and 51.46 to 181.24 kPa, respectively. MASW-derived N-values showed a strong correlation with traditional SPT N-values (R2 = 0.9401) and shear wave velocities and SPT N-values (R2 = 0.961), confirming MASW's effectiveness for geotechnical characterization. This enhances the precision of 3D soil models and supports more informed engineering decisions. The findings significantly advance the understanding of soil properties, particularly in the Niger Delta.

Novel fast full-wavefield modeling of air-coupled surface waves and its implications for non-contact pavement testing

We present a novel approach for deriving and modeling air-coupled surface waves with applications in non-contact non-destructive testing (NDT). It is based on the fast Fourier-Bessel series system in conjunction with the unconditionally stable dual-variable and position matrix method. Parametric studies, including sensitivity analysis, are conducted to assess the feasibility of using non-contact air-coupled measurements for pavement testing, focusing on Green's functions, time-domain waveforms, and experimental frequency-velocity spectra (FVS, i.e., the estimated Green's functions from acquired truncated wavefield). The predicted experimental FVS presented in this study are synthetic dispersion images, which are distinguished from the measured experimental FVS (i.e., measured dispersion images from multichannel analysis of surface wave (MASW) test). With the derived complete solution of air-coupled dynamic responses, we find that: (1) Striking similarities between the theoretical Green's functions of vertical displacement (on the pavement surface) and pressure (in the air), as well as in their corresponding experimental FVS. (2) The proposed accurate and efficient full-wave modeling of air-coupled surface waves avoids the need for good contact between geophones/accelerometers and pavement surface. This facilitates direct inversion of shear wave velocity profiles by fitting the predicted experimental FVS to the measured one. (3) Sensitivity analysis demonstrates no significant loss of information in the pressure measured in the air, supporting the feasibility of using non-contact measurement for non-destructive testing. These results suggest that non-contact air-coupled measurements hold great promise as a viable alternative to contact measurements in non-destructive testing.

Ambient noise surface-wave imaging in a hardrock environment: implications for mineral exploration

Summary The advancement of seismic methods is vital for mineral exploration in the ongoing energy transition. In this study, we investigate the application of ambient noise seismic interferometry and surface-wave analysis to characterize the subsurface in a mineral exploration context. We then confirm the results of the passive seismic investigation through an active source experiment. We collected ambient noise data using a 2D seismic line initially deployed for an active source reflection seismic study. By cross-correlating the signals, we retrieved the surface waves and constructed a 2D shear-wave velocity profile using conventional surface-wave analysis. We utilized the active source data to establish initial assumptions about the surveyed medium and then validated the passive seismic experiment. The passive seismic results are concordant with the active source results and allow for the interpretation of geological contacts and fault zones. Our work demonstrates the potential of passive seismic methods for investigating local tectonic settings and their role in hardrock mineral exploration.

Uncertainty reduction in MASW inversion and ground response analysis using a-priori information

Multi-channel analysis of surface waves (MASW) method is commonly used for estimating the shear wave velocity (Vs) profile of soil and subsequent ground response analysis (GRA). However, the uncertainties from inversion non-uniqueness in MASW testing and the absence of soil type information in GRA can undermine their accuracy and affect seismic designs. Incorporating a-priori information from geotechnical investigations can be a potential solution, but its effectiveness in MASW testing is not well understood. This study aims to quantify the uncertainty reduction achieved by incorporating a-priori information in MASW testing and GRA. Extensive numerical simulations of MASW and GRA for nine idealized soil profiles were carried out. Subsequently, the findings were validated by carrying out field testing with and without a-priori information. The a-priori information included number and thickness of soil layers, soil type and index properties. The findings demonstrate that incorporating a-priori information significantly reduces inversion uncertainty in MASW and GRA, enhancing the reliability and accuracy of the analyses. This work provides a quantitative assessment of uncertainty reduction using a-priori information for different possible site and earthquake scenarios. The findings have important implications for design of earthquake-resistant structures and for improving the safety of people and infrastructure in earthquake-prone areas.

INSIGHTS INTO THE MOHO DEPTH VARIATIONS AND CRUSTAL CHARACTERISTICS IN THE GUELMA–CONSTANTINE BASIN, NORTHEASTERN ALGERIA: A SHORT-PERIOD SEISMIC-RECEIVER-FUNCTION PERSPECTIVE

Teleseismic receiver functions (RF) were extracted from data collected at eight short-period, three-component seismic recording stations over the Guelma–Constantine Basin, northeastern Algeria, to improve the understanding of crustal structure and geodynamic processes. The H-κ stacking method was used to determine the Moho depths and average vP/vS ratios at each station. Careful linear inversion of RF was performed to determine the most appropriate average shear-wave and P-wave velocity profiles at each site. Both methods have yielded highly congruent results, with Moho depths showing robust correlations with previous seismological and geophysical studies. The previously observed pattern of the increasing Moho depth from north to south in the Tell Atlas has been confirmed. Furthermore, the identified transitional nature of the Moho in the Constantine Basin is consistent with a recent study. In addition, we identify a low-velocity zone (LVZ) at approximately 20 km depth within the southern Guelma Basin, confirming previous results in the Constantine Basin and suggesting an eastward elongation of the LVZ, at least into the southern periphery of the Guelma Basin. Examination of data from the northern tip of the Hammam Debbagh–Roknia NW–SE fault, the western boundary of the Guelma pull-apart basin, revealed a shallow Moho depth (22 km), less than the basin average depth of 25 km. The LVZ observed in the lower crust (12 km) suggests the presence of partial melts, consistent with gravimetric and chemical analyses of hydrothermal sources in the area. The extensional tectonic activity along this boundary, coupled with the low-viscosity zone and low average vP/vS ratio, is potentially associated with delamination processes. The effectiveness of our approach underscores its potential as a viable alternative or complementary method for investigating variations in the Moho depth.

Combined application of FTAN and cross-spectra analysis to ambient noise recorded by a microseismic monitoring network

A case study of seismic interferometry applied to a small microseismic monitoring network is here presented. The main objectives of this study are (i) to quantify the lateral variability of shear-wave velocities in the studied area, and (ii) to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate. Despite the limited number of stations and the short-period character of the seismic sensors, the empirical Green's functions were retrieved for all station pairs using two years of passive data. Both group and phase velocities were derived, the former using the widespread frequency-time analysis, the latter through the analysis of the real part of the cross-spectra. The main advantage of combining these two methods is a more accurate identification of higher modes, resulting in a reduction of ambiguity during picking and data interpretation. Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths. The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method, for this specific case. The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast. Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs. While the average model can be used for a more accurate location of the local natural seismicity, the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth. Finally, the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field. Errors as large as 4% have been found, related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources, and to the very short time widows. It was shown that using a one-year time window these errors are reduced to 0.3%.

Generative Adversarial Networks-Based Ground-Motion Model for Crustal Earthquakes in Japan Considering Detailed Site Conditions

ABSTRACT We develop a ground-motion model (GMM) for crustal earthquakes in Japan that can directly model the probability distribution of ground-motion acceleration time histories based on generative adversarial networks (GANs). The proposed model can generate ground motions conditioned on moment magnitude, rupture distance, and detailed site conditions defined by the average shear-wave velocity in the top 5, 10, and 20 m (VS5, VS10, and VS20) and the depth to shear-wave velocities of 1.0 km/s and 1.4 km/s (Z1.0 and Z1.4). We construct the neural networks based on styleGAN2 and introduce a novel neural network architecture to generate ground motions considering the effect of source, path, and such detailed site conditions. The resulting 5% damped spectral acceleration from the proposed GMM is consistent with empirical GMMs in terms of magnitude and distance scaling. The proposed GMM can also generate ground motions accounting for the shear-wave velocity profiles of surface soil with different magnitudes and distances and represent characteristics that are not explained solely byVS30.

Seismic site characterization baseline data for microzonation and site response analysis of Otuasega Town, Bayelsa State, Niger Delta region of Nigeria

This study presents the findings of a comprehensive geotechnical and seismic site investigation conducted at Otuasega Town located in Bayelsa State within the Niger Delta region of Nigeria. Subsurface exploration involved advancing 10 boreholes to 30 m depth using hollow stem auger drilling. Continuous disturbed and undisturbed soil sampling was performed at 1.5 m intervals for detailed geotechnical testing. Laboratory tests on the recovered soil samples established the index properties, classification, densities and consistency limits of the stratified deposits. The subsurface profile comprised alternating layers of clay, silt and sand typical of deltaic sediments, with the clay fractions exhibiting medium to high plasticity. Shear wave velocity (Vs) profiling using Multichannel Analysis of Surface Waves (WASW) techniques categorised the site predominantly as Site Class C and D based on international standards. The Standard Penetration Test (SPT) N-values ranged from 5 to 10, indicating soft normally consolidated clay conditions typical of the Niger Delta region. Predictive empirical models developed from the field and lab data showed strong correlations for estimating key geotechnical parameters such as SPT blow count, Vs and liquefaction resistance. Ground response analyses using the Vs and SPT data indicated significant site amplification potential, with peak ground accelerations up to 1.5 times the bedrock motion. Liquefaction analysis based on the empirical SPT-based methods revealed a high potential for liquefaction in the sandy layers, especially under strong earthquake shaking. The study characterized the complex sedimentology and provided baseline information for seismic microzonation and site-specific ground response analyses to advance understanding of geohazards in this delta environment.

Crustal structure of the Bushveld complex, South Africa from 1D shear wave velocity models: Evidence for complex-wide crustal modification

Thirty-nine 1D shear wave velocity profiles, obtained by jointly inverting receiver functions and Rayleigh wave group velocities, are used to investigate the crustal structure of the Bushveld Complex in northern South Africa. Data from teleseismic earthquakes recorded on broadband seismic stations between 1997 and 1999 and 2015–2020 were used to compute P-wave receiver functions. Rayleigh wave group velocities between 5 and 30 s period were obtained from an ambient noise tomography and combined with group velocities between 30 and 60 s period from a published continental-scale surface wave tomography model. Moho depths of 45–47.5 km are found under the center of the complex compared to 40 km thick crust, on average, surrounding the complex, indicating ∼5–7 of crustal thickening. The bottom ∼10 km or more of the lower crust across much of the Bushveld Complex has a Vs ≥ 4.0 km/s, consistent with a mafic composition, whereas in most areas around the margins of the complex the thickness of the mafic lower crust is much less than 10 km. In the upper crust higher velocity structure (Vs > 3.6 km/s) above 15 km depth underlain by lower velocity structure is seen in many locations, suggesting the presence of mafic/ultramafic layering. These results favor the continuous-sheet model for the structure of the Bushveld Complex because the ensemble of 1D models is characterized by three diagnostic features consistent with that model: (1) thicker crust under the center of the complex than away from the complex; (2) a greater thickness of high-velocity (i.e., mafic) layering in the lower crust under the complex compared to away from the complex; (3) high-velocity (i.e., mafic/ultramafic) layering in the upper crust beneath much of the complex. The lack of upper crustal mafic/ultramafic layering beneath some parts of the complex is consistent with the post-emplacement tectonic and magmatic history of the complex.

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.

A probabilistic full waveform inversion of surface waves

AbstractOver the past decades, surface wave methods have been routinely employed to retrieve the physical characteristics of the first tens of meters of the subsurface, particularly the shear wave velocity profiles. Traditional methods rely on the application of the multichannel analysis of surface waves to invert the fundamental and higher modes of Rayleigh waves. However, the limitations affecting this approach, such as the 1D model assumption and the high degree of subjectivity when extracting the dispersion curve, motivate us to apply the elastic full‐waveform inversion, which, despite its higher computational cost, enables leveraging the complete information embedded in the recorded seismograms. Standard approaches solve the full‐waveform inversion using gradient‐based algorithms minimizing an error function, commonly measuring the misfit between observed and predicted waveforms. However, these deterministic approaches lack proper uncertainty quantification and are susceptible to get trapped in some local minima of the error function. An alternative lies in a probabilistic framework, but, in this case, we need to deal with the huge computational effort characterizing the Bayesian approach when applied to non‐linear problems associated with expensive forward modelling and large model spaces. In this work, we present a gradient‐based Markov chain Monte Carlo full‐waveform inversion where we accelerate the sampling of the posterior distribution by compressing data and model spaces through the discrete cosine transform. Additionally, a proposal is defined as a local, Gaussian approximation of the target density, constructed using the local Hessian and gradient information of the log posterior. We first validate our method through a synthetic test where the velocity model features lateral and vertical velocity variations. Then we invert a real dataset from the InterPACIFIC project. The obtained results prove the efficiency of our proposed algorithm, which demonstrates to be robust against cycle‐skipping issues and able to provide reasonable uncertainty evaluations with an affordable computational cost.

Continuum soil‐structure‐interaction model of the LHPOST6 shaking table reaction mass at UC San Diego

AbstractThe recently upgraded six Degree‐of‐Freedom Shaking Table, LHPOST6, at UC San Diego, underwent a series of forced vibration tests to evaluate the post‐upgrade dynamic response of the foundation‐soil system. The resulting data were instrumental in obtaining frequency response curves of the system, which were used to determine its low‐strain natural frequencies, effective viscous damping ratios, and reaction mass displacements. The extensive experimental data motivated the creation of a detailed Soil‐Structure‐Interaction model of the reaction mass‐soil system. The structure and soil were modeled using 3D Finite Elements in STKO‐OpenSees and calibrated with the acquired data via a parametric study. The 3D continuum FE model and the calibration procedure based on a single soil parameter (shear wave velocity profile) proved to be an effective tool to reproduce the experimental results accurately. This paper describes the Finite Element model, its calibration, and validation. In addition, the paper provides suggestions to simplify continuum models and promote their use in professional practice. The objectives of this campaign, alongside the growing accessibility of high‐performance computing, may serve as a step toward using SSI continuum models in the industry.

Study on the influencing factors of combined processing of active and passive surface-wave data on dispersion imaging

Active and passive surface-wave methods have garnered significant attention in the near-surface geophysical community for their non-destructive, non-invasive, low-cost, and accurate advantages in delineating subsurface shear (S)-wave velocity structures. They are increasingly being utilized to address numerous engineering and environmental problems. Surface waves obtained from actively excited sources such as a hammer and a harmonic shaker, however, lack low-frequency components, resulting in limited investigation depth. Conversely, passive surface waves such as microseisms (< 0.4 Hz, associated with natural ocean wave activity) and microtremor (>1 Hz, generated by cultural and wind sources) retrieved from ambient seismic noise typically lack high-frequency components, which is not conductive to characterizing fine near-surface structures. To overcome these frequency limitations, we employ a “mixed-source data” strategy, imposing active shot gathers into ambient noise data, to widen the frequency range of dispersion images and depth of investigative capabilities. We simulate both active and passive surface-wave data based on a two-layer model, noting that their dispersion images suffer from a mode kissing phenomenon at lower frequencies. By analyzing influencing factors such as the amplitude intensity, the signal-to-noise ratio and the excitation locations of active shot gathers, as well as the length of passive surface-wave data, we better understand their impacts on dispersion images from mixed-source surface-wave data. Simulation tests demonstrate that processing mixed-source data can effectively distinguish the mode kissing phenomenon. Moreover, the effectiveness of this strategy in enhancing the quality of dispersion image is verified, especially when surface-wave dispersion images perform poorly in either the low- or high-frequency bands. Additionally, a real-world example further demonstrates that processing mixed-source data offers significant advantages in improving the quality of dispersion images. This way provides a convenient and efficient measurement strategy for delineating shear-wave velocity profiles in finer shallow layers and deeper penetration depths.

Estimation of site-specific amplification function from quarter-wavelength velocity profile and its application to obtain local amplification maps

Quarter-wavelength (QWL) velocity refers to the S-wave velocity at a depth equal to one-fourth of the wavelength of the seismic waves. It provides valuable information about the characteristics of the subsurface material properties affecting seismic waves propagation. The Swiss Seismological Service (SED) network, with over 200 stations across different lithologies, offers a rich dataset to implement correlation between site properties and site amplification factors. The current study is based on a subset of 113 selected SED seismic stations for which shear-wave velocity (Vs) profiles from geophysical measurements are available. We first meticulously analyzed and adjusted them to accurately determine the bedrock depth of the sites they are located on. Using empirical spectral modelling (ESM), we computed amplification functions from recorded earthquakes, referenced to the Swiss standard rock profile with a VS30 of 1100 m/s. Then we performed a bivariate analysis between empirical amplification functions, QWL velocity and QWL velocity contrast. The resulting coefficients are used to predict elastic amplification at frequencies between 0.5 and 10 Hz, with predictions generally being within 20–28% of observed values. We also developed an empirical equation relating VS30 (time-averaged Vs to 30 m depth) and κ0 (high-frequency attenuation parameter), to incorporate the anelastic term in the amplification. Applying our method to a 3D geophysical model of Visp, a high seismic hazard zone in Switzerland, we found it effective in predicting 1D site amplification, but noted caution in areas susceptible of 2D/3D resonance effects. A calibration function based on observed vs. predicted amplification and the frequencies of estimation normalized by the expected 2D resonance frequency improved predictions for Visp. Our study demonstrates that QWL velocity profiles offer a straightforward approach for characterizing site effects, which can be used in seismic hazard assessment to evaluate the potential amplification of ground motions.

Evaluation of Site Terms for Ground-Motion Models for South America

ABSTRACT The time-averaged shear-wave velocity in the top 30 m, VS30, is used to represent the site condition in many ground-motion models (GMMs). Regionalized GMMs account for regional differences in the ln(VS30) scaling by including region-specific coefficients for the site term. For example, recent GMMs developed for subduction zone earthquakes as part of the Next Generation Attenuation-Subduction (NGA-Sub) project include region-specific VS30 scaling for seven regions: South America, Central America, Japan, New Zealand, Taiwan, Cascadia (Pacific Northwest America), and Alaska. Of these seven regions, South America has the largest within-event standard deviation. One cause for the larger within-event standard deviation is a weaker relation between VS30 and the site term in South America compared to other regions. A cause for this weak correlation is that most of the VS30 values for the South American region in the NGA-Sub data set are based on proxies. The relation between VS30 and the site term improves significantly for periods less than 1 s if only stations with measured VS30 are considered, but the correlation is weak for periods greater than 1 s even for stations with measured VS30, indicating that, for the South America region, VS30 is not well correlated with the deeper shear-wave velocity profile that controls the long-period amplification. To improve the site terms for South America, we develop a model based on the horizontal-to-vertical spectral ratio (HVSR) from microtremors, similar to the approach used by Pinilla-Ramos et al. (2022) for California. The data set includes 660 recordings from 51 earthquakes recorded at 274 sites. The site-term model includes the period-dependent HVSR amplitude and the geometric mean of the average HVSR amplitude over the frequency band of 0.25 to 15 Hz. Including these two parameters reduces the standard deviation of site-specific site terms over the period range 0.3–3 s, with the largest reduction in the period range 1–2 s. This site-term model can be incorporated into subduction GMMs and implemented in probabilistic seismic hazard analyses for South America.

A Seismic Site Factor Adjustment Framework for Sloping Grounds and Its Application to Charleston, South Carolina

ABSTRACT Ground motions and seismic site factors for the design and analysis of infrastructure systems are typically computed based on one-dimensional (1-D) site response analysis, assuming the ground surface is horizontal and the earthquake waves propagate vertically near the ground surface. The 1-D results may not be applicable for sloping ground surfaces and layer interfaces, and therefore, the 1-D results must be adjusted for sloping grounds. In this study, a general framework for developing seismic site factors for sloping grounds and its application for Charleston, SC is presented. The two-dimensional (2-D) site response analysis results of sloping ground were compared with the 1-D horizontal ground results to develop slope adjustment factors. A total of 385 2-D nonlinear finite element simulations of sloping grounds were performed considering seven ground inclinations, 11 shear wave velocity profiles, and five ground motions representative of the area. The results show that the surface responses computed from sloping ground analyses significantly differ from those of horizontal ground analyses. By comparing the 1-D and 2-D results, a set of slope adjustment factors was developed and expressed as an easy-to-use chart. The slope adjustment factor at a project site can be obtained from the proposed chart using site-specific VS30 and ground inclination at a selected period.

Use of spatial autocorrelation of microtremor using L- and triangle- arrays for estimations of S-wave velocity model and soil response at ITS Surabaya

Abstract Institut Teknologi Sepuluh Nopember (ITS) Surabaya is geologically vulnerable to earthquake hazards. Therefore, preparedness is necessary to address emergencies and risks caused by earthquake at ITS Surabaya. In this study, dispersion curve were obtained through the analysis of microtremor data using the Spatial Autocorrelation (SPAC) method. Microtremor data acquisition was performed at the PLN ITS Futsal Field using L- and triangle-arrays. We successfully got dispersion curves from both arrays. The shear wave velocity (Vs) profile and the best model of Vs profile for each array were obtained by inverting the dispersion curves. By using a two-layer model, the first layer produces identical velocities with a depth of 5.44 meters for an L-array and 5.02 meters for a triangle array. In the second layer, the resulting velocities for each array have quite a difference. The best model of Vs profile for both arrays were used to analyse 1D soil response with a linear approach in the frequency domain. Parameters including Peak Ground Acceleration (PGA), 5% damped Pseudo Spectral Acceleration (PSA), and Fourier Amplitude Ratio (FAR) were obtained from the linear approach analysis that could describe the ground’s response to seismic waves. Vs influences the soil response, the higher the difference in Vs values of a layer with the layer below it (half-space), the higher the increase in PGA, PSA, and FAR.

Periodic foundation piles for the seismic protection of structures

This study introduces a novel solution for seismic protection of structures based on resonant metamaterials, named Periodic Foundation Piles, which can be integrated with deformable soil deposits in order to filter shear waves in a specific frequency range. A Periodic Foundation Pile consists of an array of vertical periodic structural elements containing internal resonators suitably tuned to one or more natural frequencies of the soil deposit. The proposed system is equivalent to a discontinuous foundation pile and can be installed in the ground with approximately the same procedure for piles, thus making it suitable also for applications in existing structures. The performance of the Periodic Foundation Piles in reducing the amplitude of a seismic motion has been shown through non-linear site response analyses carried out using real acceleration records as input motion. A soil deposit with an increasing shear wave velocity profile has been considered; the nonlinear behaviour of the soil has been modelled through the Iwan model, calibrated using a normalized secant shear modulus curve accounting for the influence of the confining pressure on the stress-strain relationship. Numerical results, presented in terms of peak acceleration, horizontal displacement and shear deformation profiles, as well as in terms of Fourier amplitude and response spectra, indicate that Periodic Foundation Piles may noticeably reduce the amplitudes of seismic motion in a suitable frequency range, hence potentially representing an innovative strategy to mitigate the effects of earthquakes on structures.