Spatial variability of liquefaction potential in regional mapping using CPT and SPT data

Cone penetration test (CPT) and standard penetration test (SPT) are widely used for the site specific evaluation of liquefaction potential and are getting increased use in the regional mapping of liquefaction hazard. This paper compares CPT and SPT-based liquefaction potential characterizations of regional geologic units using the liquefaction potential index (LPI) across the East Bay of the San Francisco Bay, California, USA and examines the statistical and spatial variability of LPI across and within geologic units. Overall, CPT-based LPI characterizations result in higher hazard than those derived from the SPT. This bias may result from either mis-classifications of soil type in the CPT or a bias in the CPT simplified procedure for liquefaction potential. Regional mapping based on cumulative distribution of LPI values show different results depending on which dataset is used. For both SPT and CPT-based characterizations, the geologic units in the area have broad LPI distributions that overlap between units and are not distinct from the population as a whole. Regional liquefaction classifications should therefore give a distribution, rather than a single hazard rating that does not provide for variability within the area. The CPT-based LPI values have a higher degree of spatial correlation and a lower variance over a greater distance than those estimated from SPTs. As a result, geostatistical interpolation can provide a detailed map of LPI when densely sampled CPT data are available. The statistical distribution of LPI within specific geologic units and interpolated maps of LPI can be used to understand the spatial variability of liquefaction potential.     

Magnitude scaling factors in liquefaction triggering procedures

A revised magnitude scaling factor (MSF) relationship for CPT-based and SPT-based liquefaction triggering analyses is presented in this paper. The revised MSF relationship incorporates functional dependency on the soil characteristics [represented by clean sand equivalent penetration resistances in the present form] as well as on earthquake magnitude. The revisions in MSF are based on the examination of cyclic testing results for a broad range of soil types and densities, analyses of strong ground motion records to develop relationships for the equivalent number of loading cycles for different soil properties, and the synthesis of those results into an MSF relationship suitable for implementation in practice. A separate study [2] showed that use of the revised MSF relationship in CPT-based and SPT-based liquefaction triggering procedures is well-supported by the case history databases. Other factors known to fundamentally influence the MSF are discussed.     

Assessment of liquefaction evaluation procedures and severity index frameworks at Christchurch strong motion stations

The objective of the study presented herein is to assess three commonly used CPT-based liquefaction evaluation procedures and three liquefaction severity index frameworks using data from the 2010–2011 Canterbury earthquake sequence. Specifically, post-event field observations, ground motion recordings, and results from a recently completed extensive geotechnical site investigation programme at selected strong motion stations (SMSs) in the city of Christchurch and surrounding towns are used herein. Unlike similar studies that used data from free-field sites, accelerogram characteristics at the SMS locations can be used to assess the performance of liquefaction evaluation procedures prior to their use in the computation of surficial manifestation severity indices. Results from this study indicate that for cases with evidence of liquefaction triggering in the accelerograms, the majority of liquefaction evaluation procedures yielded correct predictions, regardless of whether surficial manifestation of liquefaction was evident or not. For cases with no evidence of liquefaction in the accelerograms (and no observed surficial evidence of liquefaction triggering), the majority of liquefaction evaluation procedures predicted liquefaction was triggered. When all cases are used to assess the performance of liquefaction severity index frameworks, a poor correlation is shown between the observed severity of liquefaction surface manifestation and the calculated severity indices. However, only using those cases where the liquefaction evaluation procedures yielded correct predictions, there is an improvement in the correlation, with the Liquefaction Severity Number (LSN) being the best performing of the frameworks investigated herein. However scatter in the relationship between the observed and calculated surficial manifestation still remains for all liquefaction severity index frameworks.     

Documentation of soil conditions at liquefaction and non-liquefaction sites from 1999 Chi–Chi (Taiwan) earthquake

The 1999 Chi–Chi, Taiwan, earthquake provides case histories of ground failure and non-ground failure that are valuable to the ongoing development of liquefaction susceptibility, triggering and surface manifestation models because the data occupy sparsely populated parameter spaces (i.e. high cyclic stress ratio and high fines content with low to moderate soil plasticity). In this paper, we document results from several large site investigation programs conducted in Nantou, Wufeng and Yuanlin, Taiwan. The seismic performance of the investigated sites include non-ground failure building and free-field sites, building sites with partial foundation bearing failures, free-field lateral spread sites, and free-field level ground sites with sediment boils. Field and laboratory investigation protocols for the sites are described, including cone penetration testing (some with pore pressure and shear wave velocity measurements) and rotary wash borings with standard penetration testing (including energy measurements). Implications of the SPT energy measurements with respect to established guidelines for the estimation of SPT energy ratio (including short rod corrections) are presented. Finally, data for three example sites are shown that illustrate potential applications of the data set, and which also demonstrate a condition where existing liquefaction analysis procedures fail to predict the observed field performance.     

CPT-based analysis of liquefaction and re-liquefaction following the Canterbury earthquake sequence

The Canterbury region experienced widespread damage due to liquefaction induced by seismic shaking during the 4 September 2010 earthquake and the large aftershocks that followed, notably those that occurred on 22 February, 13 June and 23 December 2011. Following the 2010 earthquake, the Earthquake Commission directed a thorough investigation of the ground profile in Christchurch, and to date, more than 7500 cone penetration tests (CPT) have been performed in the region. This paper presents the results of analyses which use a subset of the geotechnical database to evaluate the liquefaction process as well as the re-liquefaction that occurred following some of the major events in Christchurch. First, the applicability of existing CPT-based methods for evaluating liquefaction potential of Christchurch soils was investigated using three methods currently available. Next, the results of liquefaction potential evaluation were compared with the severity of observed damage, categorised in terms of the land damage grade developed from Tonkin & Taylor property inspections as well as from observed severity of liquefaction from aerial photography. For this purpose, the Liquefaction Potential Index (LPI) was used to represent the damage potential at each site. In addition, a comparison of the CPT-based strength profiles obtained before each of the major aftershocks was performed. The results suggest that the analysis of spatial and temporal variations of strength profiles gives a clear indication of the resulting liquefaction and re-liquefaction observed in Christchurch. The comparison of a limited number of CPT strength profiles before and after the earthquakes seems to indicate that no noticeable strengthening has occurred in Christchurch, making the area vulnerable to liquefaction induced land damage in future large-scale earthquakes.     

Lessons learned and need for instrumented liquefaction sites

Instrumented sites provide essential information for understanding and modeling of ground response and ground deformation. For example, significant new lessons were learned from responses at the Wildlife Liquefaction Array (WLA) including: (1) soil softening led to lengthening of period of transmitted ground motions; (2) soil softening also led to attenuation of short-period spectral accelerations (<0.7 s); (3) amplification of long period motions (>0.7 s) occurred due to liquefaction-induced ground oscillation; and (4) ground oscillation led to a continued rise of pore water pressures after strong ground shaking ceased. A new and expanded instrumented site is being developed 70 m downstream from the old WLA site as part of the NSF Network for Earthquake Engineering Simulation (NEES). The new site has more accelerometers, piezometers and ground deformation measurement devices and the data will be streamed to the NEES-grid in near real time.     

Soil behaviour under dynamic loading conditions: experimental procedures and statistical trends

Soil behaviour under dynamic conditions is a crucial component of several studies concerned with the environmental effects of earthquakes, risk assessment and geological hazards. This work presents experimental procedures for studying soil behaviour under repeated loading and examines the role of uncertain parameters in the expected soil performance. Laboratory techniques are used to obtain stress controlled cyclic triaxial soil measurements, and a series of test programs are performed in order to study general statistical trends in the response of soils under simulated earthquake conditions, to investigate the processes causing soil failure, and to examine factors that may influence the results obtained in the laboratory.This work has been supported by grants from the National Institute of Environmental Health Sciences (Grant no. P42-ES05948 and P30-ES10126) and the Army Research Office (DAAG55-98-1-0289).     

Evaluation of liquefaction potential using neural-networks and CPT results

In this research, a reliable Cone Penetration Test data set was gathered with a wide range of parameters. This data was incorporated in a Neural-Networks computer software called STATISTICA Neural-Networks. The back propagation algorithm with a multilayer perceptron network is utilized to analyze the liquefaction occurrence in different sites. In this study, different sets of effective parameters for the neural-network analyses are selected such that to reduce the noise and to obtain more accurate results.Considering the relative importance of effective parameters in liquefaction assessment, it is indicated that σ₀, σ′₀ together play a more important role than what previously was assumed and hence the relative importance of the q_c and seismic parameters are decreased compared with the previous works. The results presented here have more accuracy than previous works while at the same time, the range of the parameters used in this study is much wider than what was previously used. This range of parameters makes the proposed method applicable for practical purposes.     

The use of MASW method in the assessment of soil liquefaction potential

The multi-channel analysis of surface wave (MASW) method is a non-invasive method recently developed to estimate shear wave velocity profile from surface wave energy. Unlike conventional SASW method, multi-station recording permits a single survey of a broad depth range and high levels of redundancy with a single field configuration. An efficient and unified wavefield transform technique is introduced for dispersion analysis and on site data quality control. The technique was demonstrated in the assessment of soil liquefaction potential at a site in Yuan Lin, Taiwan. The shear wave velocity and liquefaction potential assessments based on MASW method compares favorable to that based on SCPT shear wave measurements. Two-dimensional shear wave velocity profiles were estimated by occupying successive geophone spreads at several sites in central western Taiwan, at some of which sand boils or ground cracks occurred during 1999 Chi Chi earthquake. Liquefaction potential analysis based on MASW imaging was shown to be effective for estimating the extent of potential liquefaction hazard.     

PEER-NGL project: Open source global database and model development for the next-generation of liquefaction assessment procedures

The Next-Generation Liquefaction (NGL) project was launched to (1) substantially improve the quality, transparency, and accessibility of case history data related to ground failure; (2) provide a coordinated framework for supporting studies to augment case history data for conditions important for applications but poorly represented in empirical databases; and (3) provide an open, collaborative process for model development in which developer teams have access to common resources and share ideas and results during model development. Work to date has focused on compiling high-value case histories, developing a database template, and planning for needed supporting studies. We describe the project motivation, explain and illustrate how data resources will be compiled and organized, summarize preliminary results from ongoing data collection, describe needed supporting studies, and review project status and next steps.     

Assessment of liquefaction potential based on peak ground motion parameters

Conventionally, evaluation of liquefaction potential of loose saturated cohesionless deposits as specified in Japanese design codes employs peak ground acceleration (PGA). However, recent large-scale earthquakes in Japan revealed that liquefaction at some sites did not occur even though large PGAs were recorded at or near these sites. As an alternative approach, an evaluation procedure based on peak ground motion parameters, i.e. incorporating both PGA and the peak ground velocity (PGV), is proposed. By performing parametric studies using one-dimensional seismic response analysis and formulating regression models, seismic-induced shear stresses within the deposit are expressed in terms of peak ground motion parameters at the surface, and these are used to calculate the factor of safety against liquefaction. Application to case histories in Japan indicates that the proposed two-parameter equation can adequately account for the occurrence and non-occurrence of liquefaction at various sites as compared to the conventional PGA-based approach. Moreover, analyses of several strong motion records at various sites show that liquefaction may occur when PGA≥150 gal and PGV≥20 kine, indicating that these values can serve as thresholds in assessing the possible occurrence of liquefaction.     

Characterization of the uncertainty of the Robertson and Wride model for liquefaction potential evaluation

As the profession moves toward the performance-based earthquake engineering design, it becomes more important and pressing to examine the uncertainty of the limit state model used for liquefaction potential evaluation. In this paper, the uncertainty of the Robertson and Wride model, a simplified model for liquefaction resistance and potential evaluation based on cone penetration test, is investigated in detail for its model uncertainty in the framework of first-order reliability analysis. The uncertainties of the parameters used in the Robertson and Wride model are also examined. The model uncertainty is estimated by calibration with a fairly large set of case histories. The results show that the uncertainty of the Robertson and Wride model may be characterized with a mean-to-nominal of 0.94 and a coefficient of variation of 0.15 based on the case histories examined.     

Influence of colloidal silica grout on liquefaction potential and cyclic undrained behavior of loose sand

Cyclic triaxial tests were performed to investigate the influence of colloidal silica grout on the deformation properties of saturated loose sand. Distinctly different deformation properties were observed between grouted and ungrouted samples. Untreated samples developed very little axial strain prior to the onset of liquefaction. However, once liquefaction was triggered, large strains occurred rapidly and the samples collapsed within a few additional loading cycles. In contrast, grouted sand samples experienced very little strain during cyclic loading. Additionally, the strain accumulated uniformly throughout loading rather than rapidly prior to collapse and the samples never collapsed. Cyclic triaxial tests were done on samples stabilized with colloidal silica at concentrations of 5, 10, 15, and 20%. In general, samples stabilized with higher concentrations of colloidal silica experienced very little strain during cyclic loading. Sands stabilized with lower concentrations tolerated cyclic loading well, but experienced slightly more strain. Thus, treatment with colloidal silica grout significantly increased the deformation resistance of loose sand to cyclic loading.     

Response spectrum predictions for potential near-field and far-field earthquakes affecting Hong Kong: rock sites

To realistically assess the seismic risk relating to built infrastructures in Hong Kong and in the neighbouring coastal cities of southern Guangdong province, it is necessary to predict ground shaking induced by different earthquake scenarios with good accuracy. A companion paper has described the modelling of the spatial and temporal distribution of the diffused seismic activities in the region, based on the newly-developed ‘Expanding Circular Disc’ (ECD) method. Representative Magnitude–Distance (M–R) combinations for both near-field and far-field earthquakes (in relation to Hong Kong) have been derived using the ECD method. The present paper describes the modelling of the response spectrum on rock sites associated with the predicted M–R combinations, using the Component Attenuation Model (CAM) that was also developed recently by the authors, based on stochastic simulations of the seismological model. The significant effects of soil resonance on the response spectrum are described in a separate publication.The accuracy of CAM in modelling ground motion properties on rock sites has been tested here by comparisons with (i) strong motions recorded in Taiwan and South China from the 1999 ‘Chi-Chi’ earthquake in Taiwan (M=7.6), (ii) motions recorded in South China from another earthquake occurring in the southern Taiwan Strait in the same year (M=5.1), and (iii) historical seismic intensity data obtained within South China. The overall capability of CAM in modelling both near-field and far-field attenuation has been shown to be unmatched by existing empirical models. Results of the comparison studies confirm the accuracy of CAM, particularly within an epicentral distance of 300–400 km.This study shows that the developed serviceability response spectra (i.e. at short return periods) are controlled mainly by the earthquake recurrence behaviour of major distant seismic sources. In contrast, the ultimate response spectra (i.e. at long return periods) relate to events with magnitudes close to the maximum credible earthquake (MCE) limit, the effect of which may also be represented by the Characteristic Response Spectrum (CRS). Both types of earthquake scenario can be significantly affected by the regional crustal properties. The proposed response spectrum envelopes have been compared with previously developed recommendations, and a critical review has been conducted. The intrinsic advantages of the ECD–CAM modelling approach have been highlighted, emphasising its directness and transparency when compared with the more complex process required to implement traditional Probabilistic Seismic Hazard Assessment (PSHA).     

Neural network model for liquefaction potential in soil deposits using Turkey and Taiwan earthquake data

Simplified methods have been practiced by researchers to assess nonlinear liquefaction potential of soil. Derived from several field and laboratory tests, various simplified procedures such as stress-based, strain-based, Chinese criteria, etc. have been developed by utilizing case studies and undisturbed soil specimens. In order to address the collective knowledge built up in conventional liquefaction engineering, an alternative general regression neural network model is proposed in this paper.To meet this objective, a total of 620 sets of data including 12 soil and seismic parameters are introduced into the model. The data includes the results of field tests from the two major earthquakes that took place in Turkey and Taiwan in 1999 and some of the desired input parameters are obtained from correlations existing in the literature.The proposed GRNN model was developed in four phases, mainly: identification phase, collection phase, implementation phase, and verification phase. An iterative procedure was followed to maximize the accuracy of the proposed model. The case records were divided randomly into testing, training, and validation datasets.Generating a model that takes into account of 12 soil and seismic parameters is not feasible by using simplified techniques; however, the proposed GRNN model effectively explored the complex relationship between the introduced soil and seismic input parameters and validated the liquefaction decision obtained by simplified methods. The proposed GRNN model predicted well the occurrence/nonoccurrence of soil liquefaction in these sites. The model provides a viable tool to geotechnical engineers in assessing seismic condition in sites susceptible to liquefaction.     

Centrifuge modeling of the behavior of caisson-type quay walls during earthquakes

A series of 2-D centrifuge modeling tests with an in-flight shaker were carried out in order to model both the deformation characteristics of backfill and the seismic responses of caisson-type walls embedded in soils with various permeabilities. The rotational and translational modes were found to be in phase or various degrees out of phase with each other for quay walls embedded in soils with varying permeabilities. The alternative pumping and suction processes in excess pore water pressure that are caused by a wall's vibrations increase the level of damage because large amounts of backfill are forcedly leaked into the sea. The test results show that the rotational mode makes the dominant contribution to the changes in excess pore water pressure and in the earth pressure in the deep layers of soil, but the translational mode makes the dominant contribution to these pressures in the shallow layers. The average shear wave velocities were found to decrease rapidly to values as low as 1/8th of the velocity measured at the beginning of shaking.     

Optical dating of liquefaction features to constrain prehistoric earthquakes in Upper Assam, NE India—some preliminary results

Optical dating of liquefied sand structures formed during major earthquakes in Upper Assam, northeast of India, has been carried out to constrain the timing of prehistoric earthquakes in this seismically active region. The bleaching of source material and of the same material during the creation of a liquefaction feature was tested using quartz extracted from 21 samples associated with two different liquefaction structures in Upper Assam. Due to the poor sensitivity of quartz from this region, a sensitisation procedure was used in our SAR protocol to reduce the scatter in optical ages. Various internal consistency tests of the measurement protocol and the excellent reproducibility of the OSL ages indicate that the dose estimates from the quartz are accurate and the optical ages reliable. The preliminary OSL ages indicate that the liquefaction features were formed between 1430 AD and 1630 AD. This study demonstrates that using OSL, ‘direct dating’ of prehistoric earthquakes may be possible, if sand blows from liquefied dykes are preserved.     

Fuzzy neural network models for liquefaction prediction

Integrated fuzzy neural network models are developed for the assessment of liquefaction potential of a site. The models are trained with large databases of liquefaction case histories. A two-stage training algorithm is used to develop a fuzzy neural network model. In the preliminary training stage, the training case histories are used to determine initial network parameters. In the final training stage, the training case histories are processed one by one to develop membership functions for the network parameters. During the testing phase, input variables are described in linguistic terms such as ‘high’ and ‘low’. The prediction is made in terms of a liquefaction index representing the degree of liquefaction described in fuzzy terms such as ‘highly likely’, ‘likely’, or ‘unlikely’. The results from the model are compared with actual field observations and misclassified cases are identified. The models are found to have good predictive ability and are expected to be very useful for a preliminary evaluation of liquefaction potential of a site for which the input parameters are not well defined.     

区域土壤地震液化预测简化方法

利用可大范围获取的空间参数给出区域土壤液化的评估方法,可预估震前各地区土壤液化可能性,可快速评估震后震区土壤液化情况,对预防减轻地震灾害以及地震快速响应救援都具有重要意义的.本文以我国1976唐山大地震液化调查资料为背景,提出了区域土壤液化预测的四参数简化评估方法,并用2011年新西兰基督城地震进行了检验.选取场地平均剪切波速V_s30、复合地形指数CTI、场地到河流的距离DR以及地面峰值加速度PGA,分别代表土壤密实程度、饱水状态、地质年代和遭遇的地震作用大小,并基于我国唐山地震液化区域调查资料生成的样本,运用经典二元Logistic回归方法,建立了区域土壤液化预测模型和评估公式,回归检验结果良好,整体回判成功率为77%,非液化和液化区域成功率分别为73%和78%.将公式应用于评估2011年新西兰基督城地震液化情况,区域液化和非液化成功率分别为82%和88%,总体成功率为84%.以上结果表明本文模型可靠,方法正确,可用于区域土壤液化震前预估以及土壤液化震后快速评估.