Ground motion simulations for seismic stations in southern and eastern Romania and seismic hazard assessment

This research focuses on the evaluation of soil conditions for seismic stations in southern and eastern Romania, their influence on stochastic finite-fault simulations, and the impact of using them on the seismic hazard assessment. First, the horizontal-to-vertical spectral ratios (HVSR) are evaluated using ground motions recorded in 32 seismic stations during small magnitude (M _W  ≤ 6.0) Vrancea seismic events. Most of the seismic stations situated in the southern part of Romania exhibit multiple HVSR peaks over a broad period range. However, only the seismic stations in the eastern-most part of Romania have clear short-period predominant periods. Subsequently, stochastic finite-fault simulations are performed in order to evaluate the influence of the soil conditions on the ground motion amplitudes. The analyses show that the earthquake magnitude has a larger influence on the computed ground motion amplitudes for the short- and medium-period range, while the longer-period spectral ordinates tend to be influenced more by the soil conditions. Next, the impact of the previously evaluated soil conditions on the seismic hazard results for Romania is also investigated. The results reveal a significant impact of the soil conditions on the seismic hazard levels, especially for the sites characterized by long-period amplifications (sites situated mostly in southern Romania), and a less significant influence in the case of sites which have clear short predominant periods.     

Nonlinear numerical assessment of the seismic response of hillside RC buildings

Major damage has been reported in hilly areas after major earthquakes,primarily because of two special conditions:the variation in the seismic ground motion due to the inclined ground surface and the irregularities caused by a stepped base level in the structure.The aim of this study is to evaluate possible differences in the responses of Chilean hillside buildings through numerical linear-elastic and nonlinear analyses.In the first step,a set of response-spectrum analyses were performed on four simplified 2D structures with mean base inclination angles of 0°,15°,30°,and 45°.The structures were designed to comply with Chilean seismic codes and standards,and the primary response parameters were compared.To assess the seismic performance of the buildings,nonlinear static(pushover)and dynamic(time-history)analyses were performed with SeismoStruct software.Pushover analyses were used to compare the nonlinear response at the maximum roof displacement and the damage patterns.Time-history analyses were performed to assess the nonlinear dynamic response of the structures subjected to seismic ground motions modified by topographic effects.To consider the topographic modification,acceleration records were obtained from numerical models of soil,which were calculated using the rock acceleration record of the Mw 8.01985 Chilean earthquake.Minor differences in the structure responses(roof displacements and maximum element forces and moments)were caused by the topographic effects in the seismic input motion,with the highly predominant ones being the differences caused by the step-back configuration at the base of the structures.High concentrations of shear forces in short walls were observed,corresponding to the walls located in the upper zone of the foundation system.The response of the structures with higher angles was observed to be more prone to fragile failures due to the accumulation of shear forces.Even though hillside buildings gain stiffness in the lower stories,resulting in lower design roof displacement,maximum roof displacements for nonlinear time-history analyses remained very close for all the models that were primarily affected by the drifts of the lower stories.Additionally,vertical parasitic accelerations were considered for half the time-history analyses performed here.The vertical component seems to considerably modify the axial load levels in the shear walls on all stories.     

地震反应分析中输入界面的选取对地表地震动参数的影响

地震反应分析中输入界面选取合理与否对设计地震动参数有重要影响。基于唐山地区钻孔剖面,分别选取剪切波速为500m/s的硬黏土和800m/s的岩石顶面作为基岩输入界面,采用一维等效线性化方法讨论中硬场地输入界面的选取对地表地震动参数的影响,结果表明:(1)地表峰值加速度放大倍数及地表加速度反应谱特征周期都随输入界面深度的增加而递增,且这种递增与输入地震动的强度及频谱特性都有密切联系;(2)随着输入界面深度的增加,地表加速度反应谱几乎全频段内增大,仅在短周期内出现减小的情况,但幅度十分有限;(3)中硬场地地震反应分析中基岩输入界面宜取剪切波速为800m/s的土层顶面。     

Slope topography effects on ground motion in the presence of deep soil layers

An extensive investigation has been made into the interaction between topographic amplification and soil layer amplification of seismic ground motion. This interaction is suggested in the literature as a possible cause for the differences between topographic amplification magnitudes observed in field studies and those obtained from numerical analysis. To investigate this issue a numerical finite element (FE) parametric study was performed for a slope in a homogeneous linear elastic soil layer over rigid bedrock subjected to vertically propagating in-plane shear waves (Sv waves). Analyses were carried out using two types of artificial time history as input excitation, one mimicking the build-up and decay of shaking in the time histories of real earthquake events, and the other to investigate the steady-state response. The study identified topographic effects as seen in previous numerical studies such as modification of the free-field horizontal motion, generation of parasitic vertical motion, zones of alternating amplification and de-amplification on the ground surface, and dependence of topographic amplification on the frequency of the input motion. For the considered cases, topographic amplification and soil layer amplification effects were found to interact, suggesting that in order to accurately predict topographic effects, the two effects should not be always handled separately.     

Coupling of topographic and stratigraphic effects on seismic response of slopes through 2D linear and equivalent linear analyses

In this paper the seismic response of simple slope geometries under vertically propagating in-plane shear waves (SV waves) is assessed through two-dimensional finite element analyses to investigate the amplification of the ground motion induced by soil topography. Topographic horizontal and vertical amplification factors were evaluated through different sets of analyses focused on slopes in homogeneous half space and on slopes overlying either a rigid or a compliant bedrock. Soil was assumed to behave as a linear visco-elastic or as an equivalent-linear visco-elastic material. In the analyses the effects of slope inclination and of the characteristics of the input motion were also investigated.In order to calibrate the numerical model, the results obtained in linear visco-elastic analyses were compared with the results of parametric numerical analyses available in the literature, showing a good agreement. The results confirmed that a complex interaction exists between stratigraphic and topographic effects on the amplification of the ground motion and that the two effects cannot be evaluated independently and easily uncoupled. In the case of compliant bedrock the effect of the impedance ratio was also investigated.The results of the equivalent-linear analyses pointed out the remarkable dependence on soil non-linear behavior and, when compared to the results of linear visco-elastic analyses, showed that without accounting for soil non-linear behavior, topographic amplification factors may result underestimated.     

3D Response analysis of an instrumented hill at Matsuzaki, Japan, by a spectral method

Strong ground motion observed at an instrumented hill site is first analysed through the standard (SSR) and the horizontal-to-vertical (HVSR) spectral ratio techniques. A reasonable agreement is found between these approaches. The observations are then compared with 3D numerical simulations, performed with a highly efficient numerical code based on a spectral method, that allowed for reasonable computer times also on a PC. The observed amplification is significantly higher than that computed with a 3D homogeneous model of the mountain, suggesting that local response is governed by large-scale and small-scale soil heterogeneities rather than by topographic site effects. The introduction of a local near-surface inclusion of nonhomogeneous soil material under one of the recording stations has not significantly improved the numerical results. The observed data are also compared with the results of simplified simulations, either using 2D homogeneous models or coupling the 3D response with a 1D local soil profile. The results of such simplified approaches are discussed and their usefulness is emphasised.     

Exploration of S-wave velocity profiles at strong motion stations in Eskisehir,Turkey, using microtremor phase velocity and S-wave amplification

We have explored 1D S-wave velocity profiles of shallow and deep soil layers over a basement at strong motion stations in Eskisehir Province, Turkey. Microtremor array explorations were conducted at eight strong motion stations in the area to know shallow 1D S-wave velocity models. Rayleigh wave phase velocity at a frequency range from 3 to 30 Hz was estimated with the spatial autocorrelation analysis of array records of vertical microtremors at each station. Individual phase velocity was inverted to a shallow S-wave velocity profile. Low-velocity layers were identified at the stations in the basin. Site amplification factors from S-wave parts of earthquake records that had been estimated at the strong motion stations by Yamanaka et al. (2017) were inverted to the S-wave velocities and Q-values of the sedimentary layers. The depths to the basement with an S-wave velocity of 2.2 km/s are about 1 km in the central part of the basin, while the basement becomes shallow as 0.3 km in the marginal part of the basin. We finally discussed the effects of the shallow and deep sedimentary layers on the 1D S-wave amplification characteristics using the revealed profiles. It is found that the shallow soil layers have no significant effects in the amplification at a frequency range lower than 3 Hz in the area.     

Neo-deterministic seismic hazard scenarios for India—a preventive tool for disaster mitigation

Current computational resources and physical knowledge of the seismic wave generation and propagation processes allow for reliable numerical and analytical models of waveform generation and propagation. From the simulation of ground motion, it is easy to extract the desired earthquake hazard parameters. Accordingly, a scenario-based approach to seismic hazard assessment has been developed, namely the neo-deterministic seismic hazard assessment (NDSHA), which allows for a wide range of possible seismic sources to be used in the definition of reliable scenarios by means of realistic waveforms modelling. Such reliable and comprehensive characterization of expected earthquake ground motion is essential to improve building codes, particularly for the protection of critical infrastructures and for land use planning. Parvez et al. (Geophys J Int 155:489–508, 2003) published the first ever neo-deterministic seismic hazard map of India by computing synthetic seismograms with input data set consisting of structural models, seismogenic zones, focal mechanisms and earthquake catalogues. As described in Panza et al. (Adv Geophys 53:93–165, 2012), the NDSHA methodology evolved with respect to the original formulation used by Parvez et al. (Geophys J Int 155:489–508, 2003): the computer codes were improved to better fit the need of producing realistic ground shaking maps and ground shaking scenarios, at different scale levels, exploiting the most significant pertinent progresses in data acquisition and modelling. Accordingly, the present study supplies a revised NDSHA map for India. The seismic hazard, expressed in terms of maximum displacement (Dmax), maximum velocity (Vmax) and design ground acceleration (DGA), has been extracted from the synthetic signals and mapped on a regular grid over the studied territory.     

Topographic site effects of Xishan Park ridge in Zigong city,Sichuan considering epicentral distance

Ground motion amplification induced by topography plays a vital role in engineering seismology. A topographic array of 8 accelerographs has been operating along the ridge in Xishan Park since 2007. The topographic site effects in Zigong city are studied based on the strong motion data of 2008 M_s 8.0 Wenchuan earthquake (the epicentral distance?=?225 km) and 2019 M_s 5.2 Zizhong earthquake (the epicentral distance?=?29 km). We compare the peak ground acceleration (PGA) of the two earthquakes and find that the PGA of Station 7#, which locates on a relatively steep slope, is amplified by 4.41 times comparing with the reference station in Zizhong earthquake, while this value is only 1.62 in Wenchuan earthquake. Fourier amplitude spectrum shows that the high frequency content of Zizhong earthquake is more abundant because of its smaller epicentral distance. By using the standard spectral ratio (SSR) method, we conclude that the amplification occurs because high-frequency ground motion is likely to resonate at small-scale features. Finally, the 3D numerical simulations are used to verify these conclusions. Our work indicates that more sophisticated numerical models need to be established for more accurate topographic site effects quantification. In addition, the influence of nearby topographic features should be considered when selecting reference stations.

     

Directional dependence of site effects observed near a basin edge at Aegion,Greece

We study site effects using 520 weak motion earthquake records from a vertical array in Aegion, Greece. The array is inside a basin, has four stations in soil, and one in bedrock (178 m depth). The site is marked by high seismicity and complex surface geology. We first use the records to establish the downhole accelerometer orientations and their evolution with time. Then we estimate site effects using empirical spectral ratios with and without a reference site (standard and horizontal-to-vertical spectral ratio). We find significant site amplification which cannot be accounted for by 1D model predictions, along with a significant difference in the amplification level between the two horizontal components. These are indications of 2D effects, namely surface waves generated at the basin edge. The difference in amplification between the horizontal components is maximised when these are rotated with respect to the orientation of the basin edge. The strongest amplification takes place in the direction parallel to the basin edge (SH, or out-of-plane motion), and is up to 2 times higher than in the perpendicular direction (SV, or in-plane motion). This directional effect on the amplification is corroborated by numerical 2D modelling using incident SH and SV waves, with the former possibly generating strong Love waves. In the records, the directionality is clear for windows containing the largest amplitudes of the records (S waves and strong surface waves), while it tends to vanish for coda-wave windows. This directionality is also observed when using response spectral ratios rather than Fourier ratios. We compute soil-to-rock amplification factors for peak ground acceleration (PGA) and find it is significantly higher than what is predicted by current design codes. We attribute this difference to the basin edge amplification, linear soil behaviour, and to the inability of simple scalar values like PGA to describe complex amplification effects. Finally, we analyse the earthquake records at a surface station near the slope crest and do not observe significant topographic amplification.     

Ground motion simulations for the 23 October 2011 Van,Eastern Turkey earthquake using stochastic finite fault approach

The 23 October 2011 Van (Mw 7.1) earthquake that occurred in Eastern Turkey resulted in heavy damage particularly in the city of Van and town of Ercis. This paper presents ground motion simulations of Van earthquake by using stochastic finite fault method (EXSIM, Motazedian and Atkinson in Bull Seismol Soc Am 95:995–1010, 2005; Boore in Bull Seismol Soc Am 99:3202–3216, 2009) that provides a simple and effective tool to generate high frequency strong motion. The input parameters related to source, path, and site effects are calibrated on the basis of minimizing the error functions between simulations and observations both in time and frequency domain. Validated model parameters are used to produce synthetics in regional extent with the aim of understanding the level and distribution of the ground shaking particularly in the near fault region where no recordings are available within the 40 km of the epicenter. This paper evaluates the effect of two different slip models on ground motion intensity measures over the area of interest and addresses the variability in the near fault region associated with the source effect. The synthetics are compared with the corresponding estimations of ground motion prediction equations by Boore and Atkinson (Earthq Spectra 24:99–138, 2008), Akkar and Bommer (Seismol Res Lett 81:195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100:2978–2995, 2010). Our results indicate that despite the limitation of the method for incorporating the directivity effect and inadequate representation of the soil conditions at the individual stations, a satisfactory match between synthetics and observations are obtained both in time and frequency domain. Spatial distributions of the synthetics in regional level also show reasonable correlation with ground motion prediction equations and damage observations.     

Simulation of Ground Motion from Strong Earthquakes of Kamchatka Region (1992–1993) at Rock and Soil Sites

To estimate the parameters of ground motion in future strong earthquakes, characteristics of radiation and propagation of seismic waves in the Kamchatka region were studied. Regional parameters of radiation and propagation of seismic waves were estimated by comparing simulations of earthquake records with data recorded by stations of the Kamchatka Strong Motion Network. Acceleration time histories of strong earthquakes (M _w = 6.8–7.5, depths 45–55 km) that occurred near the eastern coast of Kamchatka in 1992–1993 were simulated at rock and soil stations located at epicentral distances of 67–195 km. In these calculations, the source spectra and the estimates of frequency-dependent attenuation and geometrical spreading obtained earlier for Kamchatka were used. The local seismic-wave amplification was estimated based on shallow geophysical site investigations and deep crustal seismic explorations, and parameters defining the shapes of the waveforms, the duration, etc. were selected, showing the best-fit to the observations. The estimated parameters of radiation and propagation of seismic waves describe all the studied earthquakes well. Based on the waveforms of the acceleration time histories, models of slip distribution over the fault planes were constructed for the studied earthquakes. Station PET can be considered as a reference rock station having the minimum site effects. The intensity of ground motion at the other studied stations was higher than at PET due to the soil response or other effects, primarily topographic ones. At soil stations INS, AER, and DCH the parameters of soil profiles (homogeneous pyroclastic deposits) were estimated, and nonlinear models of their behavior in the strong motion were constructed. The obtained parameters of radiation and propagation of seismic waves and models of soil behavior can be used for forecasting ground motion in future strong earthquakes in Kamchatka.     

平面SH波作用下部分充填圆弧形沉积谷的二维土层和地形放大效应

场地效应通常包含土层放大效应和地形放大效应,为了揭示二者的相对贡献,本文构造了平面SH波作用下部分充填沉积谷的解析模型,借助于区域分解策略,在波函数展开法的框架下,提出了超定方程组解法,得到了部分充填圆弧形沉积谷对平面SH波散射的波函数级数解,而且级数解的收敛测试表明了超定方程组解法的必要性.通过与文献结果进行对比,验证了本文方法的正确性.通过调整解析模型中两个子区域内的材料参数,计算了沉积谷引起的场地放大效应和相应的空河谷引起的地形放大效应.对二维土层与地形效应进行对比分析,结果显示,在沉积谷内二维土层放大效应通常强于地形放大效应,而地形放大效应决定了沉积谷外的地面运动放大形态.针对最大地面运动,进行了沉积谷和相应空河谷的参数分析,进一步描述了二维土层放大效应,研究结果表明二维土层放大效应引起的最大地面运动通常远远大于地形放大效应引起的最大地面运动,并且二维土层效应通常随着土层与基岩的阻抗比的增大而增大,但不是一维土层放大效应与二维地形放大效应的简单线性叠加.      

Seismic ground motion amplification pattern induced by a subway tunnel: Shaking table testing and numerical simulation

A series of 1 g shaking table tests, followed by the numerical simulations, is performed to investigate the effect of a circular subway tunnel on the ground motion amplification pattern. Effects of various parameters, including shear wave velocity of soil, frequency content of input motion, flexibility ratio and tunnel depth on the amplification pattern is investigated. Experimental study revealed that the tunnel did not affect free field response at dimensionless period greater than 10. Subsequent parametric study demonstrated that the amount of amplifications were mainly controlled by dimensionless period, dimensionless depth and flexibility ratio. Tunnel effect on the amplification pattern is more prominent for dimensionless period between 3 to 10, flexibility ratio greater than 1 and dimensionless depth less than 3. The study revealed that subway tunnel influences the seismic response of low period buildings, constructed above the tunnel.     

On the selection of GMPEs for Vrancea subcrustal seismic source

The Vrancea subcrustal seismic source is characterized by large magnitude ( $M_{W} \ge 7$ ) intermediate-depth earthquakes that occur two or three times during a century on average. In this study several procedures are used to grade four candidate ground motion prediction equations proposed for Vrancea source in the SHARE project. In the work of Delavaud et al. (J Seismol 16(3):451–473, 2012) four ground motion prediction models developed for subduction zones (Zhao et al. in Bull Seism Soc Am 96(3):898–913, 2006; Atkinson and Boore in Bull Seism Soc Am 93(4):1703–1729, 2003; Youngs et al. in Seism Res Lett 68(1):58–73, 1997; Lin and Lee in Bull Seism Soc Am 98(1):220–240, 2008) are suggested as suitable for Vrancea subcrustal seismic source. The paper presents the appropriateness analysis of the four suggested ground motion prediction equations done using a dataset of 109 triaxial accelerograms recorded during seven Vrancea seismic events with moment magnitude $M_{W}$ between 5.4 and 7.4, occurred in the past 35 years. The strong ground motions were recorded in Romania, as well as in Bulgaria, Republic of Moldova and Serbia. Based on the ground motion dataset several goodness-of-fit measures are used in order to quantify how well the selected models match with the recorded data. The compatibility of the four ground motion prediction models with respect to magnitude scaling and distance scaling implied by strong ground motion dataset is investigated as well. The analyses show that the Youngs et al. (Seism Res Lett 68(1):58–73, 1997) and Zhao et al. (Bull Seism Soc Am 96(3):898–913, 2006) ground motion prediction models have a better fit with the data and can be candidate models for Probabilistic Seismic Hazard Assessment.     

Preliminary results of strong ground motion simulation for the Lushan earthquake of 20 April 2013, China

The earthquake occurred in Lushan County on 20 April, 2013 caused heavy casualty and economic loss. In order to understand how the seismic energy propagates during this earthquake and how it causes the seismic hazard, we simulated the strong ground motions from a representative kinematic source model by Zhang et al. (Chin J Geophys 56(4):1408–1411, 2013) for this earthquake. To include the topographic effects, we used the curved grids finite difference method by Zhang and Chen (Geophys J Int 167(1):337–353, 2006), Zhang et al. (Geophys J Int 190(1):358–378, 2012) to implement the simulations. Our results indicated that the majority of seismic energy concentrated in the epicentral area and the vicinal Sichuan Basin, causing the XI and VII degree intensity. Due to the strong topographic effects of the mountain, the seismic intensity in the border area across the northeastern of Boxing County to the Lushan County also reached IX degree. Moreover, the strong influence of topography caused the amplifications of ground shaking at the mountain ridge, which is easy to cause landslides. These results are quite similar to those observed in the Wenchuan earthquake of 2008 occurred also in a strong topographic mountain area.     

Amplification of earthquake ground motion by steep topographic irregularities

The problem of amplification of seismic waves by surface topographic irregularities is addressed through analytical and numerical investigations. First, a closed‐form expression for estimating the fundamental vibration frequency of homogeneous triangular mountains is obtained, using Rayleigh's method. Subsequently, numerical modelling based on the spectral element approximation is used to study the 3D seismic response of several real steep topographic irregularities excited by vertically propagating plane S‐waves. A topographic amplification factor is obtained for each case, by a suitable average of the ratio of acceleration response spectra of output vs input motion. The 3D amplification factors are then compared with those derived by 2D analyses as well as with the topographic factors recommended in Eurocode 8 for seismic design. Copyright © 2002 John Wiley & Sons, Ltd.     

Canyon topography effects on ground motion at Feitsui damsite

The free-field accelerograms along Feitsui Canyon are analyzed and modeled by a numerical scheme to study the effect of canyon topography. Since six strong-motion accelerometers (SC1–SC6) were deployed along the Feitsui Canyon in 1991; there are 14 earthquakes (4.9≤M_L≤6.6) recorded by these stations until June 1996, but only five triggered all six stations. The maximum PGA value is 68.6 cm s⁻² recorded at station SC1. According to the present data, the effect of the dam on the ground motions at canyon stations can be negligible. The amplitude of ground motion on the slopes of the canyon is bigger than that at its trough. The integral equation method is applied to a two dimensional model of Feitsui Canyon to study the effects of the canyon topography. We choose the ground motion of SC3 or SC4 station at the trough of the canyon as the input motion for the model, which is then used to predict the ground motion at the other five stations. Apart from the earthquake close to the damsite, the simple model can reproduce the observed accelerations at all frequencies below 4 Hz. Overall, the numerical method can well predict the ground motion along the canyon, although the high-frequency simulation is underestimated.     

The influence of the 3D morphology and cavity network on the seismic response of Castelnuovo hill to the 2009 Abruzzo earthquake

During the Abruzzo earthquake (6.IV.2009, M_W = 6.3), the village of Castelnuovo, lying on an elliptical hill about 60 m high, underwent an intensive damage (I_MCS = 9–10), that could be partly ascribed to the topographic amplification and to the presence of an underground cavity network. To verify these hypotheses, the seismic response of the hill was carefully investigated adopting both 2D and 3D finite difference numerical models. Analyses were carried out using a detailed geotechnical model, defined on the basis of a comprehensive field investigation (boreholes, DH, HVSR, ERT) and accurate laboratory tests (RC-TS). The reference input motion was reproduced considering the time history of the mainshock recorded at an accelerometric station close to the epicenter, conveniently deconvoluted to the bedrock and scaled in amplitude to the site of Castelnuovo. The results of the numerical analyses, expressed in terms of distribution of the amplification factor of peak acceleration and Housner intensity, proved that the topographic effects significantly influenced the ground motion at surface, whereas the role of cavities seemed to be negligible.