In this paper, we investigated the characteristics of nonlinear site response (NLSR) at 23 S-net seafloor sites using strong-motion records obtained during three Mw 7 class earthquakes that occurred directly beneath the network. During the earthquakes, horizontal peak accelerations as large values as 1,400 and 1700 cm/s2 were recorded at the land (KiK-net) and S-net stations, respectively. The S-net is a large-scale inline-type seafloor observation network for earthquake and tsunami in the Japan Trench area. Characterization of NLSR is important because, in most common cases, it can cause a reduction of higher frequency components and a shift of predominant site frequency to lower one. Obtaining high-quality strong-motion records at seafloor sites is extremely difficult and expensive. Some of the records from the three earthquakes used in this study were contaminated by the rotations of the sensor houses, resulting in the ramps and offsets after the arrival of strong S-wave phases. We used a time window of 10 s starting from the S-wave onset, that avoided the ramps and offsets mostly. Using the horizontal-to-vertical spectral ratio (HVSR) technique, we found that the selected S-net sites might have experienced substantial degrees of NLSR during the three earthquakes with peak accelerations greater than about 60 cm/s2. To investigate that the obtained features of NLSR were realistic or not at the S-net sites, we examined the NLSR at nine KiK-net sites on land where high-quality strong-motion records were obtained. We found that the KiK-net sites experienced various degrees of NLSR during the three earthquakes, and the obtained characteristics of NLSR at the KiK-net and S-net sites were comparable. We found that the NLSR affected the ground motions at frequencies mainly higher than 1 Hz at both Kik-net and S-net sites. Despite these similarities, by analyzing the spectral ratios between two horizontal component records, we suspected that the induced rotations contributed to some extent in exaggerating the degree of NLSR at the S-net sites, primarily when the components perpendicular to the cable axes were used. We concluded that consideration of induced rotational effects is necessary to understand the NLSR at the S-net sites better.
Site condition impact on seismic ground motion has been a complex but important subject in earthquake hazard analysis. Traditional studies on site amplification effect are either based on site response via wave propagation simulation or regression analysis using parameters such as Vs30, bedrock ground motion and site response period. Ground Motion Prediction Equations (GMPEs) are used for regions where there is limited data of seismic records. The main issues with these approaches are that they cannot demonstrate the complex relationship between site amplification and its various affecting parameters, thus there exists large uncertainty in the results. Recent studies based on machine learning have shown significant improvement in predicting the site amplification, but the result is not well explained. This study assembled the information on 6 parameters including Vs30, magnitude, epicentral distance, earthquake source depth, bedrock ground motion, and altitude of 353,327 records observed during 1997 and 2019 from 698 KiK-net stations. Three machine learning algorithms of Random Forest (RF), XGBoost, and Deep Neural Networks (DNN) were implemented to predict the site amplification factor using these 6 selected parameters. Shapley Additive explanation (SHAP) was used to explain the importance of the 6 parameters. The results show that all three machine learning algorithms performed much better than the traditional GMPE approach with XGBoost’s performance the best. The explanation provided by the SHAP analysis further enhanced the reasonability of this study. It is anticipated that the combination of machine learning and SHAP analysis can provide better assessment for site amplification of ground motion with better potential of future application in seismic hazard analysis.
A large shaking-table test of a loess slope with a geometric similarity ratio of 1:25 was established as a prototype model of the loess landslide in Yongguang Village, which occurred during the 2013 Minxian–Zhangxian M6.6 earthquake. Combined with finite-element–numerical-simulation software, the dynamic response characteristics of the loess slope were compared and analyzed under the action of the same Minxian horizontal seismic wave characteristics. The results show that the seismic waves observed in the model test and numerical simulation clearly exhibit non-linear amplification effects along the slope surface and in the loess vertical direction, reaching a maximum value at the top of the slope. The ground-motion acceleration waveforms observed in the model test and numerical simulation at each observation point of the slope were similar, and the horizontal peak ground acceleration (PGA) obtained from the numerical simulation at these locations was larger than that of the slope top. The horizontal aSimulated-max/aModel-max values ranged from 1.5 to 2.2 and the vertical aSimulated-max/aModel-max values were more discrete in the range of 0.5–2.7, while the horizontal and vertical-PGA-amplification factors at the slope surface were similar. The Fourier amplitude and main horizontal frequencies were mainly distributed between 0.1 and 12 Hz. The Fourier spectrum observed in the model test was approximately 0.5 Hz, while the main frequency of the vertical Fourier spectrum was approximately 4 Hz. Additionally, the main frequency of the horizontal and vertical Fourier spectra observed in the numerical simulation was approximately 2 and 1.2 Hz, respectively.
The topographic site effect plays a vital role in controlling the characteristics of earthquake ground motions. Due to its complexity, the factors affecting topographic amplification have not been fully identified. In this study, 100 ground motion simulations generated by double-couple point sources in the homogeneous linear elastic half-space are performed based on the 3D (three-dimensional) Spectral Element Method, taking the Menyuan area of Qinghai Province, China as a local testbed site. A relationship between incident direction and the strength of topographic amplification has been observed. The horizontal ground motion is affected by the back-azimuth, which is typically chosen to be the direction from seismic station to seismic source measured clockwise from north. Specifically, the east-west PGA (Peak Ground-motion Acceleration) is significantly amplified when back-azimuth is about or , and the north-south PGA is significantly amplified when back-azimuth is around or . The vertical ground motion is affected by the dipping angle, which is the angle from vertical at which an incoming seismic wave arrives. The vertical PGA is strongly amplified when the seismic wave is almost horizontally incident (e.g., dipping angle = ). A correlation study between geomorphometric parameters and frequency-dependent topographic amplification indicates that relative elevation and smoothed curvature contain similar information, both of which are closely related to the topographic amplification of horizontal components, but not the vertical component. Our study reveals the influence of source and propagation path on topographic amplification and provides a reference for considering the topographic site effect in real engineering sites.
This study proposes a seismic input method for layered slope sites exposed to obliquely-incident seismic waves which transforms the waves into equivalent nodal forces that act on the truncated boundary of a finite element model. The equivalent nodal forces at the left and right boundaries are obtained by combining the free field response of a one-dimensional layered model with a viscoelastic boundary. The equivalent nodal forces at the bottom boundary are obtained by combining the incident wave field with the viscoelastic boundary. This proposed seismic input method for slope sites exposed to obliquely incident seismic waves is implemented with the aid of MATLAB software; it is applied to the seismic response analysis of slope sites in the commercial finite element ABAQUS software. The calculation results are compared with the reference solutions obtained by using the extended model to verify the correctness of the established seismic input method. The proposed seismic input method is then employed to investigate the influencing factors of the seismic response of layered slope sites exposed to oblique incidence P waves. The results show that the angle of incidence, location of the interface between soft and hard rocks, and impedance ratio have significant effects on the seismic landslide.
There is a consensus in earthquake engineering that different soil structures have significant effects on ground motion. Three main methods can be used to study this issue: theoretical analysis, numerical simulation, and earthquake records analysis. Notably, the combined analysis of strong-motion records and station site information appears to be the most effective and reliable approach. In this study, we used site information from the Seismic Monitoring Array of Site and Structure (SMASS) of the Institute of Disaster Prevention and a record of the earthquake that occurred in Guye, China: a point on the bedrock was considered as reference for studying the amplification effect of the soil structure on ground motion through traditional spectral ratio analysis. Our main conclusions are the following. (1) The SMASS array data indicated that the shallow soil had an amplifying effect on ground motion: the frequency bands amplifying the horizontal ground motion were mainly those between 3.4 and 6.8, 11.3–13, and 15–20.8 Hz, while the frequency bands amplifying the vertical ground vibration were mainly those between 5 and 9.9, 13.5–15.9, and 18.5–21.9 Hz (2) The SMASS array data indicated that the influence of the deep soil on the horizontal ground motion frequency component was more significant than that of vertical ground motion. (3) Overall, the SMASS array site amplified the effects of both the high (3.4–20.8 Hz) and low (0.4–2.2 Hz) frequency components of ground motion with the depth decreases of buried depth; notably, the amplification effect of the high-frequency components grew much faster than that of the low-frequency components. (4) The groundwater table had an amplification effect on the vertical ground motion; therefore, its influence is also important to consider.
Joint inversion of horizontal-to-vertical spectral ratios (HVSRs) and dispersion curves (DCs) from seismic noise recordings has been extensively used to overcome the lack of inversion uniqueness in the noise-based HVSR (NHV) or DC inversions alone. Earthquake recordings contain information about the structural properties of sedimentary layers and provide body-wave data complementary to seismic noise recordings to estimate site velocity structures, particularly in the high-frequency band. We propose a joint inversion of the Rayleigh wave DC obtained from array measurements and earthquake-based HVSR (EHV). The EHV is derived from earthquake motions rather than from microtremors based on the diffuse-field theory of plane waves. We investigated the complementarity of EHV and surface-wave DC in the joint inversion through sensitivity analyses. The DC is sensitive to bedrock shear-wave velocities in the low-frequency range and is supplemented to some degree by the EHV in the high-frequency range. The EHV is more sensitive to sediment thicknesses almost over the entire frequency range. The joint inversion is implemented by a hybrid global optimization scheme that combines genetic algorithm (GA) and simulated annealing (SA) to avoid premature convergence in the GA. The sensitivity of inversion parameters was tested to demonstrate that the P- and S-wave velocities and thicknesses of soil layers are the dominant parameters influencing EHV and DC responses. The proposed method was validated by using synthetic models to compare the joint inversion with EHV or DC inversions alone. The joint inversion was applied to the Garner Valley Downhole Array (GVDA) data for identifying the velocity structures of the site based on earthquake and noise observations. The inversion results for the P- and S-wave velocities and thicknesses of soil layers strongly suggest that the joint inversion is an efficient method to estimate site velocity structures.
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.