Evaluating empirical regression equations for Vs and estimating Vs30 in northeastern Taiwan

The shear wave velocity of shallow sediments is very important in seismic wave amplification and thus Vs30 is a well-known parameter for site classification. Based on the relationship between Vs and soil indexes, the empirical regression equations were evaluated using more than 600 data sets from the Ilan area and the Taipei Basin. Multivariable analysis, which can increase the accuracy of regression equations, was used in this study. Three extrapolations were compared, and then the most accurate bottom-constant extrapolation was adopted to estimate Vs30 at 16 boreholes, none of which reached a depth of 30 m. Ultimately Vs30 was derived for 110 free-field strong motion stations, and the stations in northeastern Taiwan were reclassified according to the Vs30-based NEHRP provisions. Regression equations of Vs and the extrapolation of Vs30 were also applied to boreholes with only N-values and corresponding depths of less than 30 m for assessing the Vs30 and NEHRP class.     

Site effect assessment using KiK-net data: part 2—site amplification prediction equation based on f<Subscript>0</Subscript> and Vsz

This paper presents empirical correlations between amplification factors and simple site parameters derived from a large subset of the KiK-net data. The amplification factor is estimated from the ratios between the surface and down-hole horizontal response spectra, corrected for the varying depths and impedance of the down-hole sites (Cadet et al. in Site effect assessment using KiK-net data—part 1—a simple correction procedure for surface/downhole spectral ratios, 2011). Several site parameters are selected on the basis of their simplicity and availability at relatively low cost. They are the shallow time-average velocities V_SZ, with z equal to 5, 10, 20 and 30 m, and the fundamental frequency f₀. The amplification factors are then correlated with each of the individual site parameters; four other “twin-parameter”—couples (f₀, V_SZ)—are also considered and the correlation with amplification factors is performed through a normalization of the frequencies by each site fundamental frequency. The quality of the correlations is given by a misfit compared with the original data variance. The largest variance reduction is obtained with twin-parameter characterizations, out of which the couple (f₀, V_S30) proves to provide the lower misfit. The performance of single parameter correlations is relatively lower; however, the best single parameter proves to be the fundamental frequency, which provides smaller misfit than the Vsz parameters. A comparison is also performed with the amplification factors recommended in European regulations, showing that it is possible right now to significantly improve both the site characterization criteria and the associated amplification factors, for use in building codes and microzonation studies.     

New code site classification, amplification factors and normalized response spectra based on a worldwide ground-motion database

A global dataset of more than 3,000 ground motion records from 536 sites from Greece, Italy, Turkey, USA and Japan is compiled and used to propose code-oriented elastic acceleration response spectra and soil amplification factors for a new site classification system, which, besides the classical geotechnical parameters $N_{SPT}, S_{u}$ and PI, uses also the fundamental period of the site, the thickness of soil deposits and the average shear wave velocity to the seismic bedrock, instead of $V_{s,30}$ . We propose a new classification system with the associated amplification factors and normalized response spectra for two seismicity levels, i.e. $M_{s}\le 5.5$ and $M_{s}>5.5$ . Uncertainties in the estimation of soil amplification factors are captured using a logic-tree approach, which allows the efficient use of alternative models and methods. The aim of this work is to improve the present EC8 soil classification. The effectiveness of the proposed classification system is compared to that of EC8 classification system using an error term, which represents the average dispersion of data within all categories of a given classification scheme. Error terms for the new classification system are lower than the error terms for EC8 classification system at all periods.     

Seismological model parameters for northeastern and its surrounding region of India

A point source seismological model is used in this study to model the available strong motion accelerograms recorded by 17 events and to calculate three seismological model parameters, the source, path, and quality factor. Due to the paucity of recorded events, this is the first time these model parameters have been obtained for the northeastern and its surrounding region of India. The quality factors of the horizontal and vertical components of recorded events with corresponding standard deviations are Q_H(f) = 188.55f^0.94, σ₁ value (25, 0.025) and Q_V(f) = 169.76f^0.93, σ₁ value (20, 0.03), respectively. The source parameter stress drop values (μσ) vary within 124–180 bars for the subduction region and 80–169 bars for the active region. The Kappa factors for the horizontal and vertical components of recorded events on the soft rock site are 0.06 and 0.05, respectively. These seismological model parameters obtained in this study will be useful for future work deriving a ground motion attenuation relation based on a spectral model. Finally, these results are useful for seismic hazard assessment of a region having sparsely recorded events.     

InterPACIFIC project: Comparison of invasive and non-invasive methods for seismic site characterization. Part II: Inter-comparison between surface-wave and borehole methods

The InterPACIFIC project was aimed at assessing the reliability, resolution, and variability of geophysical methods in estimating the shear-wave velocity profile for seismic ground response analyses. Three different subsoil conditions, which can be broadly defined as soft-soil, stiff-soil, and hard-rock, were investigated. At each site, several participants performed and interpreted invasive measurements of shear wave velocity (Vs) and compression wave velocity (Vp) in the same boreholes. Additionally, participants in the project analysed a common surface-wave dataset using their preferred strategies for processing and inversion to obtain Vs profiles. The most significant difference between the invasive borehole methods and non-invasive surface wave methods is related to resolution of thin layers and abrupt contrasts, which is inherently better for invasive methods. However, similar variability is observed in the estimated invasive and non-invasive Vs profiles, underscoring the need to account for such uncertainty in site response studies. V_S,30 estimates are comparable between invasive and non-invasive methods, confirming that the higher resolution provided by invasive methods is quite irrelevant for computing this parameter.     

厚层淤泥质粉质黏土的动力特性参数对地表地震动参数的影响

以天津滨海某软弱场地为例,在不同强度不同相位的地震动时程输入下,用等效线性化技术考虑土的动力非线性特性,计算水平成层场地的地震反应,研究天津地区厚层淤泥质粉质黏土的动剪模量比和阻尼比与剪应变幅值的关系曲线的变异性对深软场地地表峰值加速度及其反应谱的影响。结果表明:(1)在阻尼比较均值加减一倍标准差对地表地震动参数变化影响不大;(2)在动剪模量比均值增加或减1倍标准差,对2%超越概率水平下地表峰值加速度影响显著,对应地表峰值加速度可能增加或减小20%~30%,其地表反应谱谱型也有变高变瘦或变矮变胖的趋势;(3)在动剪切模量均值减1倍标准差时,对10%超越概率水平下的地表峰值加速度的影响比较显著,减小幅度在15%左右,反应谱也存在变矮变胖的趋势;在动剪切模量均值加1倍标准差时,对10%超越概率水平下的地表地震动参数影响不明显;(4)在动剪切模量或阻尼比均值加减1倍标准差时,对63%的地表地震动参数影响均不明显。     

Seismic microzoning from synthetic ground motion parameters: case study, Santiago de Cuba

Synthetic seismograms (P-SV and SH waves) have been calculated along six profiles in Santiago de Cuba basin, with a cutoff frequency of 5 Hz, by using a hybrid approach (modal summation for a regional 1D structure plus finite differences for a local 2D structure embedded in the first). They correspond to a scenario earthquake of M_S=7 that may occur in Oriente fault zone, directly south of the city. As initial data for a seismic microzoning, the characterisation of earthquake effects has been made considering several relative (2D/1D) quantities (PGDR, PGVR, PGAR, DGAR, I_AR—ratios of peak ground values of displacement, velocity and acceleration, and of design ground acceleration and Arias intensity-, etc.) and functions representative of the ground motion characteristics in soil (2D) with respect to bedrock (1D). The functions are the response spectra ratio RSR(f), already routinely used in this kind of work, and the elastic energy input ratio E_IR(f), defined, for the first time, in this paper. These data, sampled at 115 sites within all the profiles have been classified in two steps, using logical combinatory algorithms: connected components and compact sets. In the first step, from the original ground motion parameters or functions extracted from the synthetic seismograms, nine sets have been classified and the partial results show the spatial distribution of the soil behaviour as a function of the component of motion. In the second step, the results of the classification of the nine sets have been used as input for a further classification that shows a spatial distribution of sites with a quasi-homogeneous integral ground motion behaviour. By adding the available geological surface data, a microzoning scheme of Santiago de Cuba basin has been obtained.     

Earthquake response spectra taking account of number of response cycles

Proposed is a new definition of earthquake response spectra, which takes account of the number of response cycles N. The Nth largest amplitude of absolute acceleration response of a linear oscilator with natural period T and damping ratio h, which is subjected to ground motion at its base, is defined as S_A(T, h, N). By defining a reduction factor η(T, h, N) as S_A(T, h, N)/S_A(T, h, 1), characteristics of η(T, h, N) were investigated based on 394 components of strong motion records obtained in Japan. Two practical empirical formulae to assess the reduction factor η(T, h, N) are proposed.     

A note on a pseudo-natural SSI frequency for coupled soil–pile–structure systems

The notion of a pseudo-natural SSI frequency was introduced in a recent publication by the authors, as the frequency where foundation motion is minimized with respect to the free field surface motion. This frequency is determined analytically in this paper, for a single-degree-of-freedom structure supported on a pile foundation. The analytical solution is compared to numerical results from rigorous finite element analyses for different pile and structural configurations. The relationship between pseudo-natural (f_pSSI) and effective natural SSI frequency (f_SSI) of the coupled system is also analytically quantified. It is concluded that f_pSSI may deviate substantially from f_SSI when a stiff squatty structure is founded on a stiff and/or short end-bearing pile for which foundation translation prevails. Conversely, when a flexible tall structure is supported on a flexible pile, f_pSSI and f_SSI nearly coincide due to dominant base rocking effects. In the latter case the effective natural SSI frequency can be predicted by standard identification procedures even when free-field recordings are missing. Effective damping effects are also discussed.     

Effects of strong lateral discontinuity on ground motion characteristics and aggravation factor

This paper presents the effects of strong lateral discontinuity (SLD) in basins on the ground motion characteristics, differential ground motion (DGM) and aggravation factor (aggravation factor is simply the extra spectral amplification due to complex 2D site effects over the 1D response of the soil column). The seismic responses of open- and closed-basin models with SLD were simulated using SH-wave finite difference algorithm with fourth-order spatial accuracy. Simulated seismic responses, DGM and its spectra revealed that SLD induces Love waves, and the lower cutoff frequency for the same is equal to the fundamental frequency (f₀) of the soil layer. The maximum average aggravation factor (AAF) and DGM were obtained near the edge of open basin and at the centre of closed basin. A decrease of amplitude of Love wave, DGM level and AAF with offset from SLD was observed in an open basin. On the other hand, in closed basin, spatial variation of AAF and DGM level was highly variable. Duration of shaking, AAF and DGM level was more in the closed basin than in an open basin. Increase of DGM level and AAF with decrease of the width of basin was observed.     

Rock-physics-based carbonate pore type characterization and reservoir permeability heterogeneity evaluation, Upper San Andres reservoir, Permian Basin, west Texas

The use of the shear wave velocity data as a field index for evaluating the liquefaction potential of sands is receiving increased attention because both shear wave velocity and liquefaction resistance are similarly influenced by many of the same factors such as void ratio, state of stress, stress history and geologic age. In this paper, the potential of support vector machine (SVM) based classification approach has been used to assess the liquefaction potential from actual shear wave velocity data. In this approach, an approximate implementation of a structural risk minimization (SRM) induction principle is done, which aims at minimizing a bound on the generalization error of a model rather than minimizing only the mean square error over the data set. Here SVM has been used as a classification tool to predict liquefaction potential of a soil based on shear wave velocity. The dataset consists the information of soil characteristics such as effective vertical stress (σ′_v0), soil type, shear wave velocity (V_s) and earthquake parameters such as peak horizontal acceleration (a_max) and earthquake magnitude (M). Out of the available 186 datasets, 130 are considered for training and remaining 56 are used for testing the model. The study indicated that SVM can successfully model the complex relationship between seismic parameters, soil parameters and the liquefaction potential. In the model based on soil characteristics, the input parameters used are σ′_v0, soil type, V_s, a_max and M. In the other model based on shear wave velocity alone uses V_s, a_max and M as input parameters. In this paper, it has been demonstrated that Vs alone can be used to predict the liquefaction potential of a soil using a support vector machine model.     

Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai‐Tibet

Seasonal changes over 2 years (2004–2006) in soil moisture content (θ_v) of frozen alpine frost meadow soils of the Qinghai‐Tibet plateau permafrost region under three different levels of vegetation cover were investigated. Vegetation cover and air temperature changes had significant effects (synergistic effect) on θ_v and its distribution in the soil profile. During periods of soil freezing or thawing, the less the vegetation cover, the quicker the temperature drop or rise of soil water, and the shorter the duration of the soil water freeze–thaw response in the active soil layer. Under 30% and 65% vegetation cover the amplitude of variation in θ_v during the freezing period was 20–26% greater than that under 93% cover, while during the thawing period, it was 1·5‐ to 40·5‐fold greater. The freezing temperature of the surface soil layer, ^fT_s, was 1·6 °C lower under 30% vegetation cover than under 93% vegetation cover. Changes in vegetation cover of the alpine frost meadow affected θ_v and its distribution, as well as the relationship between θ_v and soil temperature (T_s). As vegetation cover decreased, soil water circulation in the active layer increased, and the response to temperature of the water distribution across the soil profile was heightened. The quantity of transitional soil phase water at different depths significantly increased as vegetation cover decreased. The influence of vegetation cover and soil temperature distribution led to a relatively dry soil layer in the middle of the profile (0·70–0·80 m) under high vegetation cover. Alpine meadow θ_v and its pattern of distribution in the permafrost region were the result of the synergistic effect of air temperature and vegetation cover. Copyright © 2007 John Wiley & Sons, Ltd.     

A probabilistic model of seismicity: Kamchatka earthquakes

The catalog of Kamchatka earthquakes is represented as a probability space of three objects {Ω, $ \tilde F $ \tilde F P}. Each earthquake is treated as an outcome ω _i in the space of elementary events Ω whose cardinality for the period under consideration is given by the number of events. In turn, ω _i is characterized by a system of random variables, viz., energy class k_i, latitude φ _i , longitude λ _i , and depth h _i . The time of an outcome has been eliminated from this system in this study. The random variables make up subsets in the set $ \tilde F $ \tilde F and are defined by multivariate distributions, either by the distribution function $ \tilde F $ \tilde F (φ, λ, h, k) or by the probability density f(φ, λ, h, k) based on the earthquake catalog in hand. The probabilities P are treated in the frequency interpretation. Taking the example of a recurrence relation (RR) written down in the form of a power law for probability density f(k), where the initial value of the distribution function f(k ₀) is the basic data [Bogdanov, 2006] rather than the seismic activity A ₀, we proceed to show that for different intervals of coordinates and time the distribution f _elim(k) of an earthquake catalog with the aftershocks eliminated is identical to the distribution f _full(k), which corresponds to the full catalog. It follows from our calculations that f ₀(k) takes on nearly identical numeral values for different initial values of energy class k ₀ (8 ≤ k ₀ ≤ 12) f(k ₀). The difference decreases with an increasing number of events. We put forward the hypothesis that the values of f(k ₀) tend to cluster around the value 2/3 as the number of events increases. The Kolmogorov test is used to test the hypothesis that statistical recurrence laws are consistent with the analytical form of the probabilistic RR based on a distribution function with the initial value f(k ₀) = 2/3. We discuss statistical distributions of earthquake hypocenters over depth and the epicenters over various areas for several periods     

Long period microtremors, microseisms and earthquake damage: Northridge, CA, earthquake of 17 January 1994

Three studies of site amplification factors, based on the recorded aftershocks, and one study based on strong motion data, are compared one with another and with the observed distribution of damage from the Northridge, CA, earthquake of 17 January 1994 (M_L=6.4). In the epicentral area, when the peak ground velocities are larger than v_m≈15 cm/s, nonlinear response of soil begins to distort the amplification factors determined from small amplitude (linear) wave motion. Moving into the area of near-field and strong ground motion (v_m>30 cm/s), the site response becomes progressively more affected by the nonlinear soil response. Based on the published results, it is concluded that site amplification factors determined from small amplitude waves (aftershocks, small earthquakes, coda waves) and their transfer-function representation may be useful for small and distant earthquake motions, where soils and structures respond to earthquake waves in a linear manner. However in San Fernando Valley, during the Northridge earthquake, the observed distribution of damage did not correlate with site amplification determined from spectra of recorded weak motions. Mapping geographical distribution of site amplification using other than very strong motion data, therefore appears to be of little use for seismic hazard analyses.     

浅硬场地剪切波速测试标准差对地震动影响初探

采用Shake 2000程序,以Turkey Flat试验场地为模型,通过输入不同类型下多种强度的地震波,计算研究多工况下剪切波速测试标准差对地表加速度反应谱和峰值加速度的影响。结论为:(1)浅硬场地上剪切波速测试标准差对地震动的影响很大,影响程度与输入地震波的强度和频率以及场地剪切波速计算值有关;(2)如果将反应谱残差大于20%或加速度峰值差别大于20%定义为统计意义上的不可忽略,那么剪切波速测试标准差对计算结果的影响在大多数情况下均不可忽略;(3)当输入波的卓越周期与场地特征周期接近时,浅硬场地上剪切波速测试标准差引起的反应谱变化非常显著;(4)只有当输入波的卓越周期与场地特征周期相差较大且输入波强度偏小时,剪切波速测试标准差引起的反应谱变化才可略去;(5)当浅硬场地上剪切波速实测结果低于统计均值时,地震动计算结果的偏差一般明显大于剪切波速实测结果,高于统计均值时引起的偏差,且地震输入越强表现越明显。     

Should average shear-wave velocity in the top 30 m of soil be used to describe seismic amplification?

The average velocity of shear waves in the top 30 m of soil, ν_L, has become the parameter used by many engineering design codes and most recently by published empirical-scaling equations to estimate the amplitudes of strong ground motion. Yet there are few studies to determine whether this is a meaningful parameter to use—and whether estimates that do use it are reliable. In 1995, the authors studied this problem and concluded that ν_L should not be used. We reported then that an older site characterization in terms of the soil site parameter proposed by Seed et al. [1], s_L, worked better because it included a measure of the thickness of the soil layers to considerably greater depths. Our report, however, made no difference; numerous papers continued to be published based only on scaling in terms of ν_L, and worse, they also ignored the geological site conditions. The purpose of this paper is to emphasize that the average shear-wave velocity in the top 30 m of soil should not be the only site parameter used to scale strong-motion amplitudes. While the search continues for the more meaningful site parameters to use in empirical scaling of strong earthquake ground motion, it is better to use s_L to describe the amplification of seismic waves by soil deposits near the surface.     

Low-velocity zone and topography as a source of site amplification effect on Tarzana hill,California

Tarzana station is located in the foothills of the Santa Monica Mountains in California near the crest of a low (<20 m) natural hill with gentle slopes. The hill is about 500 m in length by 130 m in width and is formed of extremely weathered shale at the surface to fresh at depth. Average S-wave is about 250 m/s in the top 17–18 m, and S- and P-wave velocities significantly increase below this depth. According to the NEHRP classification based on V_S30∼300 m/s it is a site class D. Strong-motion instrumentation at Tarzana consisted of an accelerograph at the top of the hill, a downhole instrument at 60 m depth, and an accelerograph at the base of the hill. More than 20 earthquakes were recorded by at least three instruments at Tarzana from 1998 till 2003. Comparisons of recordings and Fourier spectra indicate strong directional resonance in a direction perpendicular to the strike of the hill. The dominant peaks in ground motion amplification on the top of the hill relative to the base are at frequencies ∼3.6 and 8–9 Hz for the horizontal components. Our hypothesis is that the hill acts like a wave trap. This results in an amplification at predominant frequencies f=V/4 h (h is layer's thickness) at f∼3.6 Hz for S-waves (using average V_S17=246 m/s and h=17 m) and f∼7.9 Hz for P-waves (using average V_P17=535 m/s and h=17 m). As was shown by Bouchon and Barker [Seismic response of a hill: the example of Tarzana, California. Bull Seism Soc Am 1996;86(1A):66–72], topography of this hill amplifies and polarizes ground motion in the frequency range of 3–5 Hz. Hill acts as a magnifying polarizing glass: It polarizes ground motion in the direction perpendicular to the strike of the hill and also amplifies ground motions that had been also amplified by a low-velocity layer.