Characterization of liquefaction resistance in gravelly soil: large hammer penetration test and shear wave velocity approach

Gravelly soil is generally recognized to have no liquefaction potential. However, liquefaction cases were reported in central Taiwan in the 1999 Chi-Chi Taiwan earthquake and in the 1988 Armenia earthquake. Thus, further studies on the liquefaction potential of gravelly soil are warranted. Because large particles can impede the penetration of both standard penetration test and cone penetration test, shear wave velocity-based correlations and large hammer penetration tests (LPT) are employed to evaluate the liquefaction resistance of gravelly soils. A liquefied gravelly deposit site during the Chi-Chi earthquake was selected for this research. In situ physical properties of soil deposits were collected from exploratory trenches. Instrumented LPT and shear wave velocity (V_s) measurements were performed to evaluate the liquefaction resistance. In addition, large-scale cyclic triaxial tests on remolded gravelly soil samples (15 cm in diameter, 30 cm in height) were conducted to verify and improve LPT-based and V_s-based correlations. The results show that the LPT and shear wave velocity methods are reasonably suitable for liquefaction assessment of gravelly soils.     

Evaluation of soil liquefaction in the Chi-Chi, Taiwan earthquake using CPT

During the 1999 Chi-Chi, Taiwan earthquake, many sand boiling phenomena were observed in central Taiwan, which caused severe ground settlement and structure damages. According to the installed accelerograms, the peak ground surface horizontal accelerations in the liquefaction-affected areas range from 774.42 to 121.3 gal. The writers carried out an extensive investigation of soil liquefaction in this earthquake. In this paper, we present results of the CPT exploration and post-earthquake liquefaction analysis. Two hundred and seventy five (275) cone penetration test data were collected from the liquefaction-affected areas, and 46 liquefaction case histories and 88 non-liquefaction case histories were derived that can be used to evaluate the accuracy of existing liquefaction evaluation models. In addition, the strength of the liquefied soils after earthquake and the implication of its liquefaction potential in the future event are discussed.     

Evaluation of models for Fourier amplitude spectra for the Taiwan region

In the Taiwan region, the empirical spectral models for estimating ground-motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the M_L=7.3 Chi-Chi earthquake (21 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground-motion parameters (Fourier amplitude spectra, peak accelerations and response spectra), which were calculated using the models, and the observed data demonstrates that the models could provide an accurate prediction for the case of the Chi-Chi earthquake and the largest aftershocks. However, there are some peculiarities in the ground-motion frequency content and attenuation that, most probably, are caused by the features of the rupture process of the large shallow earthquake source.     

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.     

Rainfall-triggered ordinary earthquakes in Taiwan: a statistical analysis

Abstract

This study was carried out to investigate the role of rainfall in triggering ordinary earthquakes that occurred before and after the major Chi-Chi earthquake (7.2 M_L) in central Taiwan on 21 September 1999. To test a possible mechanism, earthquake activities from January 1995 to July 2012 were examined. The Chi-square (χ²) test revealed a significant difference between the correlations of monthly accumulated rainfall values and earthquake activities before and after the Chi-Chi earthquake. This result is discussed in terms of changes in crustal conditions after the Chi-Chi earthquake. Two ordinary earthquakes that may have been associated with heavy rainfall after the Chi-Chi earthquake are identified.

Editor D. Koutsoyiannis; Associate editor A. Porporato

Citation Lin, J.-W., 2014. Rainfall-triggered ordinary earthquakes in Taiwan: a statistical analysis. Hydrological Sciences Journal, 59 (5), 1074–1080.     

Characteristics of Soil Response in Near-Fault Zones During the 1999 Chi-Chi, Taiwan, Earthquake

Distribution of parameters characterizing soil response during the 1999 Chi-Chi, Taiwan, earthquake (M _w = 7.6) around the fault plane is studied. The results of stochastic finite-fault simulations performed in Pavlenko and Wen (2008) and constructed models of soil behavior at 31 soil sites were used for the estimation of amplification of seismic waves in soil layers, average stresses, strains, and shear moduli reduction in the upper 30 m of soil, as well as nonlinear components of soil response during the Chi-Chi earthquake. Amplification factors were found to increase with increasing distance from the fault (or, with decreasing the level of “input” motion to soil layers), whereas average stresses and strains, shear moduli reduction, and nonlinear components of soil response decrease with distance as ~ r ^?1 . The area of strong nonlinearity, where soil behavior is substantially nonlinear (the content of nonlinear components in soil response is more than ~40–50% of the intensity of the response), and spectra of oscillations on the surface take the smoothed form close to E(f) ~ f ^?n , is located within ~20–25 km from the fault plane (~ 1/4 of its length). Nonlinearity decreases with increasing distance from the fault, and at ~40–50 km from the fault (~ 1/2 of the fault length), soil response becomes virtually linear. Comparing soil behavior in near-fault zones during the 1999 Chi-Chi, the 1995 Kobe (M _w = 6.8), and the 2000 Tottori (Japan) (M _w = 6.7) earthquakes, we found similarity in the behavior of similar soils and predominance of the hard type of soil behavior. Resonant phenomena in upper soil layers were observed at many studied sites; however, during the Chi-Chi earthquake they involved deeper layers (down to ~ 40–60 m) than during lesser-magnitude Kobe and Tottori earthquakes.     

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.     

S-wave velocity structures of sediments estimated from array microtremor records and site responses in the near-fault region of the 1999 Chi-Chi,Taiwan earthquake

The 1999 Chi-Chi, Taiwan earthquake, M_W = 7.6, caused severe damage in the near-fault region of the earthquake. In order to evaluate site effects in the near-field strong motions we estimate S-wave velocity structures of sediments at four sites using array records of microtremors. We also recalculated S-wave velocity structures at other four sites previously reported. To show the validity of the estimated S-wave velocity structures we separate empirical site responses from aftershock records using the generalized inversion method and show the agreement between empirical and theoretical site responses. We also show an observed fact that suggests soil nonlinearity during the Chi-Chi earthquake by comparing horizontal-to-vertical spectral ratios (HVRs) for main shock records with HVRs for aftershock records. Then we calculate one-dimensional equivalent-linear site responses using the estimated S-wave velocity structures and the main shock records observed on the surface. It is found that site amplification due to thick (about 6 km) sediments is one of the important factors for explaining the long-period velocity pulses of about 5 to 10 sec observed at sites in the footwall during the Chi-Chi earthquake. It is also found that the theoretical site responses of shallow soft sediments at sites that sustained severe damage in the hanging wall shows significant amplification around 1 sec. As the amplitude of velocity pulses with period around 1 sec is most critical in causing damage to ordinary buildings of moderate heights, our results suggest that the 1-sec period velocity pulses, amplified by the site response of shallow sediments should contribute to the severe damage during the Chi-Chi earthquake.     

Non-linear site response simulations in Chang-Hwa region during the 1999 Chi–Chi earthquake, Taiwan

The destructive 1999 Chi–Chi earthquake (M_w 7.5) was the largest inland earthquake in Taiwan in the 20th century. Several observations witness the non-linear seismic soil response in sediments during the earthquake. In fact, large settlements as well as evidence of liquefaction attested by sand boils and unusual wet ground surface were observed at some sites. In this paper, we present a seismic response simulation performed with CyberQuake software on a site located within the Chang-Hwa Coastal Industrial Park during the 1999 Chi–Chi earthquake in Taiwan. A non-linear multi-kinematic dynamic constitutive model is implemented in the software. Computed NS, EW and UP ground accelerations obtained with this model under undrained and two-phase assumptions, are in good agreement with the corresponding accelerations recorded at seismic station TCU117, either for peak location, amplitudes or frequency content. In these simulations, liquefaction occurs between depths 1.3 and 11.3 m, which correspond to the observed range attested by in place penetration tests and other liquefaction analyses. Moreover, the computed shear wave velocity profile is very close to post-earthquake shear wave velocity profile derived from correlations with CPT and SPT data. Finally, it is shown that in non-linear computations, even though a 1D geometry is considered, it is necessary to take into account the three components of the input motion.     

Predicting seismic liquefaction potential of sands by optimum seeking method

The feasibility of using the optimum seeking method to assess the seismic liquefaction potential of sands has been investigated. Optimization theory is a very important branch of applied mathematics and has a wide application in the practical world. Using the available field sand liquefaction data, the influence of various factors is quantified using the optimum seeking method. The factors considered are: the earthquake magnitude M, the distance of the site from the source of the earthquake L, the depth of the water table D_w, the depth of the sand deposit D_s, and the standard penetration test (SPT) blow count N. The most important factors have been identified as the earthquake magnitude and the SPT blow count. Prediction results show that the proposed method is effective and feasible. Since neither normalization of the SPT blow count nor calculation of the seismic shear-stress ratio are required, the proposed method is simpler and more direct than the conventional methods of evaluating liquefaction potential.     

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

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

Verification of critical cyclic strength curve by Taiwan Chi-Chi earthquake data

In this paper, known simplified methods for the assessment of soil liquefaction are summarized. Their discrepancies are examined. Using the Chi-Chi earthquake data as well as other reported data, a set of three critical cyclic strength curves were obtained by finding the minimum of misclassified points. The functional forms of these three curves are an exponential function, a hyperbola, and a cubic polynomial. A lower bound critical cyclic strength curve is then established. This curve may have important applications in practice for liquefaction-related designs. Through this case study, it was found that a minimum cyclic strength CSR_lim may exist at a very low value of (N₁)₆₀. An upper limit (N₁)₆₀^upp also exists beyond which liquefaction may not occur. Furthermore, current simplified methods seem suitable only for a limited range of N values and fines content, and may fail for general applications. The lower bound curve proposed in this paper may provide an alternative approach for improvement. Since the explicit functional form and the statistical indices are available, the statistically regressed curves seem to have a benefit in that it may be used directly to conduct hazard analysis, and evaluate the uncertainty within the critical cyclic strength curves.     

Observed Evolution of Linear and Nonlinear Effects at the Dahan Downhole Array,Taiwan: Analysis of the September 21, 1999 M 7.3 Chi-Chi Earthquake Sequence

— In this paper, the site characteristics of the Dahan downhole array are studied by analyzing the September 21, 1999 M 7.3 Chi-Chi earthquake sequence including the main shock and some aftershocks. The four-level array (0 m, 50 m, 100 m and 200 m) is located to the north of Hualien City in eastern Taiwan. Polarization analysis is used to check the orientation errors of the seismometers at different levels of depth. If the surface instrument is chosen as reference, the angle between the major polarization axes of the surface and any downhole records is the orientation error that must be corrected for the downhole accelerographs. The orientation errors at depths of 50 m, 100 m and 200 m are 32°, 120° and –84°. After the corrections, the coherency between the surface and downhole records is substantially improved. Spectral ratio analysis shows that the predominant frequency of the Chi-Chi main shock shifts to a lower frequency. We also simulate ground motions at different depths by using the Haskell method with a linear velocity structure model. The record at surface is chosen as the input motion. Compared with the observed data, ground acceleration can be well reproduced for the aftershocks (weak-motion events) of the September 21, 1999 M 7.3 Chi-Chi earthquake. However, for the Chi-Chi main shock, the synthetic waveform cannot match well with the observation neither in amplitude nor in phase. This indicates that large ground shaking probably induced the nonlinear site effect at that time, and the model used cannot support it.Acknowledgement. The authors would like to express their thanks to Dr. L.F. Bonilla and one anonymous reviewer for their valuable suggestions. This research was supported by the National Science Council under grant number NSC 89-2921-M-194-007. The Institute of Earth Sciences, Academia Sinica supplied the strong-motion data. The support of these organizations is gratefully acknowledged.     

The middle-long term prediction of the February 3, 1996 Lijiang earthquake (MS=7) by the "criterion of activity in quiescence"

IntroductionIn the book Future CataS~ologr published in 1992, we proposed a viewpoiflt on using the"criterion of activity in quiescence" to predict big eathquake (MsZ7) (GUO, et al, 1992), and predicted in the book that in futore several years or in ten years a big earthquake (Ms27) will be possible to occur in the Zhongdian and nearby in Yunnan Province. In the 1994 nation-wide earthquake tendency consultation meeting we pointed out, once more, in the Zhongdian region of Yunnan Province…     

Earthquake parameters and spatial distribution of coseismic effects in Southern Siberia and Mongolia

In this work we review earthquakes that happened in Southern Siberia and Mongolia within the coordinates of 42°–62° N and 80°–124° E and first propose relationships between earthquake parameters (a surface-wave earthquake magnitude M _s and an epicentral intensity(I ₀) based on the MSK-64 scale) and maximal distances from an earthquake epicenter (R _{e max}), hypocenter (R _{h max}), and a seismogenic fault (R _{f max}) to the localities of secondary coseismic effects. Special attention was paid to the study of these relationships for the effects of soil liquefaction. Hence, it was shown that secondary deformations from an earthquake were distributed in space away from an earthquake epicenter, than from an associating seismogenic fault. The effects of soil liquefaction are manifested by several times closer to a seismogenic fault, than all other effects, regardless of the type of tectonic movement in a seismic focus. Within the 40 km zone from an earthquake epicenter 44% of the known manifestations of liquefaction process occurred; within the 40 km zone from a seismogenic fault—90%. We propose the next relationship for effects of soil liquefaction: M _s = 0.007 × R _{e max} + 5.168 that increases the limits of the maximum epicentral distance at an earthquake magnitude of 5.2 ≤ M _s ≤ 8.1 as compared to the corresponding relationships for different regions of the world.     

Parallel PSM/FDM Hybrid Simulation of Ground Motions from the 1999 Chi-Chi,Taiwan, Earthquake

-- A new technique for the parallel computing of 3-D seismic wave propagation simulation is developed by hybridizing the Fourier pseudospectral method (PSM) and the finite-difference method (FDM). This PSM/FDM hybrid offers a good speed-up rate using a large number of processors. To show the feasibility of the hybrid, a numerical 3-D simulation of strong ground motion was conducted for the 1999 Chi-Chi, Taiwan earthquake (M_w 7.6). Comparisons between the simulation results and observed waveforms from a dense strong ground motion network in Taiwan clearly demonstrate that the variation of the subsurface structure and the complex fault slip distribution greatly affect the damage during the Chi-Chi earthquake. The directivity effect of the fault rupture produced large S-wave pulses along the direction of the rupture propagation. Slips in the shallow part of the fault generate significant surface waves in Coastal Plain along the western coast. A large velocity gradient in the upper crust can propagate seismic waves to longer distances with minimum attenuation. The S waves and surface waves were finally amplified further by the site effect of low-velocity sediments in basins, and caused the significant disasters.     

A practical reliability-based method for assessing soil liquefaction potential

The current simplified methods for assessing soil liquefaction potential use a deterministic safety factor in order to judge whether liquefaction will occur or not. However, these methods are unable to determine the liquefaction probability related to a safety factor. An answer to this problem can be found by reliability analysis. This paper presents a reliability analysis method based on the popular Seed'85 liquefaction analysis method. This reliability method uses the empirical acceleration attenuation law in the Taiwan area to derive the probability density distribution function (PDF) and the statistics for the earthquake-induced cyclic shear stress ratio (CSR). The PDF and the statistics for the cyclic resistance ratio (CRR) can be deduced from some probabilistic cyclic resistance curves. These curves are produced by the regression of the liquefaction and non-liquefaction data from the Chi-Chi earthquake and others around the world, using, with minor modifications, the logistic model proposed by Liao [J. Geotech. Eng. 114 (1988) 389]. The CSR and CRR statistics are used in conjunction with the first order and second moment method, to calculate the relation between the liquefaction probability, the safety factor and the reliability index. Based on the proposed method, the liquefaction probability related to a safety factor can be easily calculated. The influence of some of the soil parameters on the liquefaction probability can be quantitatively evaluated.     

Ground response analysis of Kanpur soil along Indo-Gangetic Plains

The present work deals with 1D and 2D ground response analysis and liquefaction analysis of alluvial soil deposits from Kanpur region along Indo-Gangetic plains. Standard penetration tests and seismic down hole tests have been conducted at four locations namely IITK, Nankari village, Mandhana and Bithoor at 1.5 m interval up to a depth of 30 m below the ground surface to find the variation of penetration blows and the shear wave velocity along the depth. From the selected sites undisturbed as well as representative soil samples have been collected for detailed soil classification. The soil profiles from four sites have been considered for 1D and 2D ground response analysis by applying the free field motions of three Himalayan earthquakes namely Chamba earthquake (M_w—5.1), Chamoli earthquake (M_w—6.4) and Uttarkashi earthquake (M_w—6.5). An average value of Peak Ground Acceleration (PGA) obtained from 1D and 2D analysis is considered for liquefaction analysis and post-liquefaction settlement. The excess pore water pressure ratio is greater than 0.8 at a depth of 24 m from ground surface for IITK, Nankari village, Bithoor sites. More than 50% of post liquefaction settlement is contributed by layers from 21–30 m for all sites. In general, the soil deposits in Kanpur region have silty sand and sand deposits and are prone to liquefaction hazards due to drastic decrease of cyclic resistance ratio (CRR) at four chosen sites in Kanpur.     

内蒙古中西部地区中小地震矩震级研究

基于S震相"S窗"内的波形信号识别、品质因子Q(f)和22个台站场地响应,利用2009~2016年3月内蒙古中西部地区地震的波形资料,反演了182次中小地震的震源波谱参数,得到这些小震的零频幅值及其拐角频率,据此计算了这些地震的地震矩M_0、矩震级M_W和应力降Δσ。利用回归分析方法得到了近震震级与矩震级、矩震级与应力降的关系式。分析表明,近震震级与矩震级、矩震级与应力降呈线性关系。可见,将矩震级纳入地震的快报与正式目录中,可以丰富地震观测报告内容,更好地为地震应急和地震科研服务。     

Source time functions of the 1999, Jiji (Chi-Chi) earthquake from GDSN long period waveform data using aftershocks as empirical Green’s functions

A large earthquake (M _W=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of M _W=6.1 and M _W=6.2, respectively, were used as empirical Green’s functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the M _W=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the M _W=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one. Foundation item: State Natural Science Foundation of China (49904004) and IPGP of France. Contribution No. 02FE2007, Institute of Geophysics, China Seismological Bureau.