Nonlinearity in the response of soils in the 1995 Kobe earthquake in vertical components of records

Strong motion records provided by seismic vertical arrays allow estimation of stress–strain relations in soils at depths from the surface to the location of the deepest device. As an example, time-dependent nonlinear soil behavior was estimated in vertical components of records obtained in the epicentral area of the 1995 Kobe earthquake. Degradation of the rigidity of soils in the strong motion was observed. The constructed nonlinear models of the soil behavior were used for estimating the nonlinear parts in the ground response by the nonlinear system identification technique. Nonlinear parts in the ground response were found to be as high as 50% at 2 km from the fault and 10% at 6–15 km from the fault plane. Odd types of nonlinearity, such as cubic, the fifth, seventh, etc. order, were found to be typical for soils, whereas, nonlinearities of even types are usually weak, but increase in liquefied soils.     

Inversion of waveforms from. Xiangtang borehole seismic array for soil dynamic property

In order to understand the site soil response of the Xiangtang borehole seismic array under real strong ground motion, reveal the site response, verify the technique of borehole exploration, and improve the precision of in-situ test and laboratory test, this paper presents a new approach, which is composed of two methods. One is the layered site seismic response method, whose layer transform matrix is always real. The other is a global-local optimization technique, which uses genetic algorithm (GA)-simplex method. An inversion of multi-component waveforms of P, SV and SH wave is carried out simultaneously. By inverting the records of three moderate and small earthquakes obtained from the Xiangtang borehole array (2# ) site, the soil dynamic characteristic parameters, including P velocity, damping ratio and frequency-dependent coefficient b, which has not been given in previous literatures, are calculated. The results show that the soil S wave velocity of the Xiangtang 2# borehole is generally greater than that obtained from the 1994 in-situ test, and is close to the velocity of the 3# borehole, which is more than 200 m away from the 2# borehole. Meanwhile, perceptible soil nonlinear behavior under peak ground motion of about 60×10-2m/s2 is detected by the inversion analysis. The presented method can be used for studying the soil response of other borehole array sites.     

Estimation of Nonlinear Time-dependent Soil Behavior in Strong Ground Motion Based on Vertical Array Data

To improve our understanding of nonlinear elastic properties of soils, a method is proposed of estimation of stress-strain relations of soils in situ in strong ground motion based on vertical array data. Strong motion records provided by seismic vertical arrays allow estimation of nonlinear stress-strain relations in soil layers at different depths, from the surface down to the location of the deepest device. As an example, records obtained during the main shock of the 1995 Kobe earthquake at Port-Island, SGK, and TKS sites were used to estimate the stress-strain relations in the soil profiles. For different layers, different types of nonlinear stress-strain relations were selected, according to the profiling data. To account for temporal changes in the soil behavior, consecutive parts of records were examined, and for successive time intervals, the relations were found showing the best-fit approximation to the observed data. At Port Island and SGK sites, where the strongest accelerations were recorded, the obtained stress-strain relations showed systematic changes in the upper layers (0–14 m), such as, a progressive reduction of the slopes of the stress-strain curves due to liquefaction at Port Island and reduction and recovery of the slopes at SGK and TKS sites. At the three sites, the stress-strain relations remained stable in layers below 11–14 m. Thus, the proposed approach gives us a representation of the soil behavior in layers at different depths in strong ground motion; it allows calculation of the propagation of arbitrary seismic signals in the studied profiles and estimation of nonlinear components in the ground response by the nonlinear system identification technique. The method can also be applied to evaluate the ground response at sites where profiling data are available and an imposed motion can be estimated.     

Seismological evidence for nonlinear elastic ground behavior during large earthquakes

Amplification of earthquake-induced seismic waves by soft superficial deposits often causes significant damages in the urban areas. In predicting this effect for large future earthquakes, the linear elastic response of soils is customarily assumed. To check this assumption, we have analyzed surface and downhole acceleration data from the SMART1 and SMART2 strong motion arrays in Taiwan, covering peak accelerations of up to 0·3 g. First, frequency-dependent amplification induced by the alluvial deposits at the SMART1 array was estimated using spectral ratio technique, where the records at rock site were taken as a reference motion. Statistically validated reduction in soil amplification in the strong motion relative to the weak motion in the frequency range between approximately 1 and 9 Hz was detected. Secondly, relative site responses between the Pleistocene and recent sedimentary deposits at the SMART2 array were studied. Relative amplification was shown to be clearly dependent on the excitation level. Thirdly, we compared experimentally recorded uphole/downhole spectral ratios on weak and strong ground motion with the theoretical response yielded by the geotechnical code DESRA2 which assumes hysteretic constitutive relationship of soil. Major symptoms of nonlinear ground behavior predicted by the model were found in the observed data. Back-calculation of the shear wave velocities to the depth of 47 m shows nearly 50% decrease in the strongest quakes, also accounted for by the nonlinear soil behavior.     

研究局部场地条件对地震动影响的响嘡遥测台阵

对原有的响嘡三维场地条件影响台阵进行了扩建改造。增设了一个深度为47m的测井,对钻孔土样进行了土动力学测试,弥补了缺乏详细土层特性资料的缺陷。新老测井观测点均架设了电话线路,实现了观测中心与台阵之间的远程监控和数据传输。已经获取的一批中小地震近场记录可以用来进行盲测研究。     

Detection of soil liquefaction from strong motion records

During the recent earthquakes in Japan and the U.S.A. a number of records from liquefied‐soil sites have been obtained. The ground motion parameters from these sites were studied and several methods for detection of liquefaction from seismic records were developed. The methods, however, focus mainly on the horizontal ground motion and may interpret as liquefaction‐induced some records from soft‐soil deposits or records with dominant surface waves, at which sites the phenomenon was not observed. Besides, not all of the available records from liquefied sites were processed. In this paper, after examination of the ability of different types of ground motion parameters to indicate alone soil liquefaction we propose a new liquefaction detection method that simultaneously analyses instantaneous frequency content of the horizontal and the vertical ground acceleration. We also compare performance of the proposed method with that of the other liquefaction detection methods. The computations are carried out using a common data set including records from liquefied and non‐liquefied sites. Results show that the frequency‐related parameters and the proposed method detect more efficiently the occurrence of liquefaction from the seismic records. Copyright © 2001 John Wiley & Sons, Ltd.     

Nonlinear identification of the soil response at Dahan downhole array site during the 1999 Chi–Chi earthquake

Records of the 1999 Chi–Chi earthquake, provided by the Dahan downhole array were analyzed, and nonlinear hysteretic stress–strain relations in the soil layers were estimated by the method previously used for studying the response of soils during the 1995 Kobe and the 2000 Tottori earthquakes. The obtained models of the soil behavior were applied for evaluating changes of the shear moduli in the soil layers and for the nonlinear identification of the soil behavior at Dahan site during the Chi–Chi earthquake. We found that reduction of the shear moduli in the soil layers did not exceed 5%, and the soil response was virtually linear. The content of the nonlinear components in the soil response was about 5% of the intensity of the response, and it was mostly due to the odd-order nonlinearities. A similarity in the stress–strain relations describing the behavior of soils during the 1995 Kobe earthquake, the 2000 Tottori earthquake, and the 1999 Chi–Chi earthquake was found, which indicates the possibility to describe the behavior of similar types of soils at different sites by similar stress–strain relations and predict soil behavior in future earthquakes.     

场地土层模型参数的地震动记录反演方法

考虑场地土层不同物理参数的综合反演,改进了水平与竖向谱比(Horizontal-to-Vertical Spectral Ratio,简称HVSR)混合全局优化反演方法,进一步通过土层的S、P波波速、厚度、泊松比、密度和S、P波品质因子等土层参数反演的敏感性分析,形成了可同时反演场地土层厚度及剪切波速的混合全局优化反演方法.以美国GVDA和日本IWTH27竖向强震动观测台阵为例,分别以理论HVSR及加速度观测记录获得的HVSR曲线为目标,反演获得了场地浅层速度结构,并与观测台阵场地钻孔揭示的土层模型进行比较,验证了发展的反演方法的合理性和适用性.本文研究表明,基于地震加速度记录的HVSR全局优化反演方法是获取场地土层速度结构的一种有效的途径.     

Nonlinear soil behavior during the Chi-Chi (Taiwan) earthquake of 1999 from data of the Dahan vertical array

Records of the Chi-Chi (Taiwan) earthquake of 1999 (M _w = 7.7) obtained at the Dahan vertical array (a part of the SMART-2 array) in NE Taiwan are analyzed. The Dahan array was located at a distance of ～80 km from the fault plane of the earthquake. Using the method previously applied to the analysis of soil behavior during the Japan earthquakes in Kobe (1995) and Tottori (2000), nonlinear stress-strain relations are estimated in soil layers at depths of 0–200 m at the Dahan site for successive time intervals. The resulting soil behavior model is used for estimating variations in shear moduli and for the identification of nonlinear behavior of soil layers during an earthquake. According to the model estimates, the earthquake-related decrease in shear moduli did not exceed ～5% and the soil response was nearly linear (the amount of its nonlinear components was also no more than ～5%). The stress-strain relations describing the soil responses to the earthquakes in Kobe (1995), Tottori (2000), and Chi-Chi (1999) are similar, implying that, in principle, it is possible to describe the behavior of soils of the same type by the same stress-strain relations and predict the soil response to a future earthquake.     

地震动加速度反应谱场地条件影响调整系数研究

场地条件对地震动特性影响显著,在抗震设计反应谱的确定过程中,需根据场地条件对加速度反应谱予以相应的调整。已有场地条件影响调整方案研究成果,均基于数值模拟或局部地区强震动记录统计,多数仅给出了峰值加速度PGA场地条件影响调整系数,对非线性的考虑缺乏观测数据依据。为此在全球强震动记录统计获得的PGA归一化加速度反应谱和日本钻井台阵记录获得的加速度反应谱平台值非线性衰减指数的基础上,结合钻孔模型数值模拟和近期研究成果,建立了考虑场地条件影响非线性的地震动加速度反应谱场地条件影响调整系数方案。     

Analysis of nonlinear site response using the LSST downhole accelerometer array data

Both linear and nonlinear behaviors of soil deposits were evaluated by strong and weak motion data observed on the surface and at depths of 6, 11, 17, 47 m at the Large Scale Seismic Test (LSST) array in Lotung, Taiwan. The soil properties measured by well logging and by the shear wave velocity profile measured by uphole and cross-hole methods are available. Both one-dimensional equivalent-linear method and nonlinear method are used for the evaluation have been used. The synthetic records at various depths are obtained by using the records at the bottom as input motion. These synthetic records are then compared with actual records at corresponding depths. Records of 13 earthquakes are used. We find that the synthetic records obtained from a linear model match well with actual records for small input motions, but the results obtained from a nonlinear model match poorly. On the other hand, the synthetic records using both the nonlinear model and equivalent-linear model are in good agreement with the observed records for large input motions. In these cases, the predicted response spectra using the linear model consistently overestimate the observed records. The threshold distinguishing the large and small input motions is 0.04 g at depth of 47 m for the LSST data. Thus, the nonlinearity started at 0.04 g and occurred unequivocally at 0.075 g. Furthermore, the dominant frequencies shift toward lower values when input motions become large. Clearly, the observed records at the LSST site manifest nonlinearity of soil response. The hysteresis loops evaluated by the nonlinear method show a permanent strain of about 0.01% in soil layers at higher ground motion input levels in this case.     

Experimental p-y curves for liquefied soils from centrifuge tests

The present study aims to obtain p-y curves(Winkler spring properties for lateral pile-soil interaction) for liquefied soil from 12 comprehensive centrifuge test cases where pile groups were embedded in liquefiable soil. The p-y curve for fully liquefied soil is back-calculated from the dynamic centrifuge test data using a numerical procedure from the recorded soil response and strain records from the instrumented pile. The p-y curves were obtained for two ground conditions:(a) lateral spreading of liquefied soil, and(b) liquefied soil in level ground. These ground conditions are simulated in the model by having collapsing and non-collapsing intermittent boundaries, which are modelled as quay walls. The p-y curves back-calculated from the centrifuge tests are compared with representative reduced API p-y curves for liquefied soils(known as p-multiplier). The response of p-y curves at full liquefaction is presented and critical observations of lateral pile-soil interaction are discussed. Based on the results of these model tests, guidance for the construction of p-y curves for use in engineering practice is also provided.     

The new accelerograph network for Santa Fe De Bogota, Colombia and implications for microzonation

The new accelerographic network of Santa Fe de Bogotá is composed of 29 three-component stations with sensors at the surface and three additional six-component borehole stations with three sensors at the surface and three at depth (115, 126 and 184 m). In total, 32 stations have been operative in the metropolitan area of Bogotá since 1999. During this period of time, a significant number of weak motion are recorded and used for a preliminary analysis of local site effects. Using the SH-wave response spectra we verify the behavior of the different seismic zones proposed by the previous microzonation study of the city. A comparison between normalized SH-wave response spectra and the normalized design spectra for each zone clearly depicts that parts of the design spectra should be revised, as well as the boundaries between different zones may require some changes. The spectral amplification levels reach up to a factor of 5. The predominant periods obtained by the amplification spectra in different stations in the city, show variability from 0.3 to 3.0 s. A comparison is also made between the predominant periods obtained using H/V spectral ratios of microtremors and those using weak motion. In general, microtremors tend to predict slightly lower values of dominant periods than those calculated by the weak motion spectra. However, there is a general correlation between the two data sets. Using the data recorded by one of the borehole station, an equivalent linear seismic response analysis was conducted. The modeled and recorded response spectra show similarities in period peaks, however, the modeled soil amplification is underestimated for periods less than 0.8 s. Since the available record is weak motion which represents mostly the linear response of the soils, further analysis is required.     

某长江大桥深度超过15米的砂土层液化势的综合评价

不同抗震设计规范的砂土液化判别方法或国内外其他有代表性的液化判别方法所采用的地震动参数和土性指标及其埋藏条件是不同的,因而采用这些方法对同一工程场地进行液化势预测时其评价结果通常有一些差异,甚至会得到相反的结论。为了给重大工程建设提供较为合理、可信的地基液化势预测结果,采用多种液化判别方法进行场地液化势的综合评价是比较客观的,也是必要的。本文结合某长江大桥桥基工程,采用建筑抗震设计规范的砂土液化判别方法、国内外有代表性的液化判别方法、有限元数值分析法等多种方法逐一对该工程场地砂性土层进行液化判别,并结合室内动三轴液化试验结果,对主桥墩不考虑冲刷条件和考虑一般冲刷深度5m条件时的砂性土层进行了液化势的综合评价,并将各土层的液化势分为液化、可能液化和不液化3个等级,得到了较为合理可靠的判别结果。     

饱和黄土场地原位试验及液化势评价

目前,主要依靠室内动力试验对黄土液化势进行评价。由于黄土特殊的结构性,室内试验对其饱和的过程较为复杂,且与实际场地饱和黄土差异明显,导致室内黄土液化试验结果并不能代表现场饱和黄土的抗液化强度。本文选取兰州市西固区寺儿沟村某饱和黄土场地进行钻孔测试,现场实施了标准贯入试验、静力触探试验以及剪切波速测试。应用Robertson的土类指数分类图对该场地不同含水率黄土的土类进行了界定,确定了饱和黄土属于类砂土,有液化势。应用NCEER推荐方法,计算了3组原位试验数据的饱和黄土循环抗力比（CRR）,通过与1976年唐山地震和1999年集集地震液化土CRR对比,得出了饱和黄土抗液化强度很低的结论。     

Comparison between different techniques for evaluation of predominant periods using strong ground motion records and microtremors in Pereira Colombia

This paper presents a comparison between different techniques for evaluation of predominant periods in soft soil, for the urban area of Pereira city, Western Colombia. In this study we used microtremor and strong ground motion records obtained by a local array of seven accelerographs stations deployed in the city. Response spectra and spectral ratios have been calculated and compared with strong seismic events recorded in solid rock and soft soil stations. These observations allowed the determination of dominant response spectra for several sectors in the urban area. For the microtremor measurements and earthquake data, dominant periods were determined using interpretation of Fourier amplitude spectra and Nakamura's technique. A comparison between dominant periods obtained from strong ground motion records and those obtained from microtremor measurements show similarities, which is in the range 0.2–0.5 s. A preliminary version of a site response map for Pereira city was obtained from this analysis.     

Estimation of Nonlinear Soil Behavior During the 1999 Chi-Chi, Taiwan, Earthquake

The 1999 Chi-Chi, Taiwan, earthquake (M_w = 7.6) was one of the strongest earthquakes in recent years recorded by a large number of strong-motion devices. Though only surface records are available, the obtained strong-motion database indicates the variety of ground responses in the near-fault zones. In this study, accelerograms of the Chi-Chi earthquake were simulated at rock and soil sites, and models of soil behavior were constructed at seven soil sites (TCU065, TCU072, TCU138, CHY026, CHY104, CHY074, and CHY015), for which parameters of the soil profiles are known down to depths of at least ~70 m and at 24 other soil sites, for which parameters of the soil profiles are known down to 30–40 m; all the sites were located within ~50 km from the fault. For reconstructing stresses and strains in the soil layers, we used a method similar to that developed for the estimation of soil behavior based on vertical array records. As input for the soil layers, acceleration time histories simulated by stochastic finite-fault modelling with a prescribed slip distribution over the fault plane were taken. In spite of the largeness of the earthquake’s magnitude and the proximity of the studied soil sites to the fault plane, the soil behavior at these sites was relatively simple, i.e., a fairly good agreement between the spectra of the observed and simulated accelerograms and between their waveforms was obtained even in cases where a single stress-strain relation was used to describe the behavior of whole soil thickness down to ~70–80 m during strong motion. Obviously, this is due to homogeneity in the characteristics of soil layers in depth. At all the studied sites, resonant phenomena in soil layers (down to ~40–60 m) and nonlinearity of soil response were the main factors defining soil behavior. At TCU065, TCU110, TCU115, CHY101, CHY036, and CHY039 liquefaction phenomena occurred in the upper soil layers, estimated strains achieved ~0.6–0.8%; at other stations, maximum strains in the soil layers were as high as 0.1–0.4%, according to our estimates. Thus, valuable data on the in situ soil behavior during the Chi-Chi earthquake was obtained. Similarity in the behavior of similar soils during the 1995 Kobe, 2000 Tottori (Japan), and Chi-Chi (Taiwan) earthquakes was found, indicating the possibility of forecasting soil behavior in future earthquakes. In the near-fault zones of the three earthquakes, “hard-type” soil behavior and resonant phenomena in the upper surface layers prevail, both leading to high acceleration amplitudes on the surface.     

Identification of non-linear and non-stationary soil properties during the 1995 Hyogoken-nanbu earthquake

An extended Kalman filter algorithm with local iteration is presented for the identification of non-linear and non-stationary soil properties. Borehole-array strong motions were recorded at a liquefied site during the 1995 Hyogoken-nanbu earthquake. In this study, a modified Kalman filtering method in which the extended Kalman filter is iteratively used at every local time-step to track rapid parameter changes is proposed. The method is then applied to the instrumented soil layer, which is modeled by an equivalent linear model. An identification of non-linear and non-stationary soil properties was conducted successfully; and non-linear restoring force–displacement relationships including progression with time were obtained.     

Estimation of local site effects by a generalized inversion scheme using observed records of ‘Small-Titan’

A simple-extended array of observation for strong motions, named ‘Small-Titan’, was installed in Sendai City, Japan in 1997. This array system consists of 20 observation sites, which have an average of about 4 km for the station-to-station distance. Each observation station was located at various kinds of soils to effectively obtain local site effects. Since its completion, the system has obtained about 160 motion records. As is well known, ground motions are influenced by the three main factors such as source, path and local site effects. An accurate prediction of ground motions should thus be performed based on the three factors. In this context, theoretical methods used commonly lack enough information related with source and path effects. A generalized inversion scheme to separate the effects of source, propagation path and local site from S-wave spectra was applied to the records of Small-Titan obtained in the Sendai region, Japan. The objective of this paper is to estimate local site effects in Sendai City and to obtain meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. This paper also describes the applicability of the generalized inversion scheme to seismic zoning.

Site amplification factors estimated from the generalized inversion scheme vary significantly, reflecting the effects of soil characteristics at each site. Strong site-dependency was found in the site amplification factors. In particular, the local site amplification factors in the western side are smaller than those in the eastern side, across the Nagamachi–Rifu Fault that exists in the central part of the Sendai area. Results from the generalized inversion scheme are compared with those obtained by the direct spectral ratio scheme, the H/V method and the bootstrap inversion. There is a good agreement in the amplification factors between the generalized inversion scheme, the direct spectral ratio scheme and the bootstrap inversion scheme. The H/V method for estimating site amplification is found to extract the same dominant periods as the other methods. This paper especially emphasizes a good agreement in not only dominant period, but also amplification factor.

The generalized inversion scheme is a useful method in separating the effects of source, propagation path and local site and in revealing the site amplification. In addition, the generalized inversion scheme provides meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. Moreover, this paper indicates that the generalized inversion scheme is applicable to the seismic zoning toward disaster mitigation.     

Consistency of dynamic site response at Port Island

System identification (SI) methods are used to determine empirical Green's functions (EGF) for soil intervals at the Port Island Site in Kobe, Japan and in shake table model tests performed by the Port and Harbor Research Institute (PHRI) to emulate the site during the 17 January 1995 Hyogo‐ken Nanbu earthquake. The model form for the EGFs is a parametric auto‐regressive moving average (ARMA) model mapping the ground motions recorded at the base of a soil interval to the top of that interval, hence capturing the effect of the soil on the through‐passing wave. The consistency of site response at Port Island before, during, and after the mainshock is examined by application of small motion foreshock EGFs to incoming ground motions over these time intervals. The prediction errors (or misfits) for the foreshocks, the mainshock, and the aftershocks, are assessed to determine the extent of altered soil response as a result of liquefaction of the ground during the mainshock. In addition, the consistency of soil response between field and model test is verified by application of EGFs calculated from the shake table test to the 17 January input data. The prediction error is then used to assess the consistency of behaviour between the two cases. By using EGFs developed for small‐amplitude foreshock ground motions, ground motions were predicted for all intervals of the vertical array except those that liquefied with small error. Analysis of the post‐liquefied ground conditions implies that the site response gradually returns to a pre‐earthquake state. Site behaviour is found to be consistent between foreshocks and the mainshock for the native ground (below 16 m in the field) with a normalized mean square error (NMSE) of 0.080 and a peak ground acceleration (PGA) of 0.5g. When the soil actually liquefies (change of state), recursive models are needed to track the variable soil behaviour for the remainder of the shaking. The recursive models are shown to demonstrate consistency between the shake table tests and the field with a NMSE of 0.102 for the 16 m to surface interval that liquefied. The aftershock ground response was not modelled well with the foreshock EGF immediately after the mainshock (NMSE ranging from 0.37 to 0.92). One month after the mainshock, the prediction error from the foreshock modeled was back to the foreshock error level. Copyright © 2001 John Wiley Sons, Ltd.