Subsurface Characterization at Ground Failure Sites in Adapazari, Turkey

Ground failure in Adapazari, Turkey during the 1999 Kocaeli earthquake was severe. Hundreds of structures settled, slid, tilted, and collapsed due in part to liquefaction and ground softening. Ground failure was more severe adjacent to and under buildings. The soils that led to severe building damage were generally low plasticity silts. In this paper, the results of a comprehensive investigation of the soils of Adapazari, which included cone penetration test (CPT) profiles followed by borings with standard penetration tests (SPTs) and soil index tests, are presented. The effects of subsurface conditions on the occurrence of ground failure and its resulting effect on building performance are explored through representative case histories. CPT- and SPT-based liquefaction triggering procedures adequately identified soils that liquefied if the clay-size criterion of the Chinese criteria was disregarded. The CPT was able to identify thin seams of loose liquefiable silt, and the SPT (with retrieved samples) allowed for reliable evaluation of the liquefaction susceptibility of fine-grained soils. A well-documented database of in situ and index testing is now available for incorporating in future CPT- and SPT-based liquefaction triggering correlations.     

CPT-Based Probabilistic and Deterministic Assessment of In Situ Seismic Soil Liquefaction Potential

This paper presents a complete methodology for both probabilistic and deterministic assessment of seismic soil liquefaction triggering potential based on the cone penetration test (CPT). A comprehensive worldwide set of CPT-based liquefaction field case histories were compiled and back analyzed, and the data then used to develop probabilistic triggering correlations. Issues investigated in this study include improved normalization of CPT resistance measurements for the influence of effective overburden stress, and adjustment to CPT tip resistance for the potential influence of “thin” liquefiable layers. The effects of soil type and soil character (i.e., “fines” adjustment) for the new correlations are based on a combination of CPT tip and sleeve resistance. To quantify probability for performance-based engineering applications, Bayesian “regression” methods were used, and the uncertainties of all variables comprising both the seismic demand and the liquefaction resistance were estimated and included in the analysis. The resulting correlations were developed using a Bayesian framework and are presented in both probabilistic and deterministic formats. The results are compared to previous probabilistic and deterministic correlations.     

Discriminant Model for Evaluating Soil Liquefaction Potential Using Cone Penetration Test Data

Statistical analysis using a discriminant model is applied to 399 cone penetration test (CPT) data sets of both liquefaction and nonliquefaction cases, including 174 sets from the Chi-Chi earthquake in Taiwan and 225 sets of synthesized data. The discriminant model employed is a multivariate statistical method. In situ testing results of cone tip resistance qc and sleeve friction ratio Rf are adopted as the major parameters for analyses. A model for evaluating liquefaction potential using CPT-qc data is also established in this study, which allows calculated results to be compared with the empirical curves.     

First-Order Reliability Method for Probabilistic Liquefaction Triggering Analysis Using CPT

The potential for liquefaction triggering of a soil under a given seismic loading is measured herein by probability of liquefaction. The first order reliability method (FORM) is used to calculate reliability index, from which the probability of liquefaction is obtained. This approach requires the knowledge of parameter and model uncertainties; the latter is the focus of this paper. An empirical model for determining liquefaction resistance based on cone penetration test (CPT) is established through “neural network learning” of case histories. This resistance model along with a reference seismic loading model forms a performance function or limit state for liquefaction triggering analysis. Within the framework of the FORM, the uncertainty of this limit state model is characterized through an extensive series of sensitivity studies using Bayesian mapping functions that have been calibrated with a set of quality case histories. In addition, a deterministic model for assessing liquefaction potential in terms of factor of safety is presented, and the probability-safety factor mapping functions for estimating the probability of liquefaction for a given factor of safety in the absence of the knowledge of parameter uncertainty are also established. Examples are presented to illustrate the proposed methods.     

Yield Strength Ratio and Liquefaction Analysis of Slopes and Embankments

A procedure is proposed to evaluate the triggering of liquefaction in ground subjected to a static shear stress, i.e., sloping ground, using the yield strength ratio, su(yield)/σv0′. Thirty liquefaction flow failures were back analyzed to evaluate shear strengths and strength ratios mobilized at the triggering of liquefaction. Strength ratios mobilized during the static liquefaction flow failures ranged from approximately 0.24 to 0.30 and are correlated to corrected cone and standard penetration resistances. These yield strength ratios and previously published liquefied strength ratios are used to develop a comprehensive liquefaction analysis for ground subjected to a static shear stress. This analysis addresses: (1) liquefaction susceptibility; (2) liquefaction triggering; and (3) post-triggering/flow failure stability. In particular, step (2) uses the yield strength ratio back-calculated from flow failure case histories and the cyclic stress method to incorporate seismic loading.     

CPT-Based Evaluation of Liquefaction and Lateral Spreading in Centrifuge

Two series of centrifuge model tests were conducted using Nevada sand. Four saturated models placed in a mildly inclined laminar box and simulating a 6-m-thick deposit were shaken inducing liquefaction effects and lateral spreading. The sand was deposited at a relative density, Dr = 45 or 75%; two of the 45% models were subjected to overconsolidation or preshaking. The second series involved in-flight measurements of static cone tip penetration resistance, qc, simulating the standard cone penetration test (CPT) 36-mm cone. Values of qc increased with Dr, overconsolidation, and preshaking. A normalized resistance, qc1N, was assigned to each of the four liquefaction/lateral spreading models. Increases in Dr, overconsolidation, and preshaking decreased liquefaction and ground deformation, but relative density alone captured these effects rather poorly. Conversely, qc1N predicted extremely well the liquefaction and lateral spreading response of the four models, confirming Seed’s hypothesis to explain the success of penetration-based seismic liquefaction charts. The depth to liquefaction measured in the four centrifuge models is consistent with the field CPT liquefaction chart.     

Normalizing the CPT for Overburden Stress

Effective overburden stress can have a significant influence on cone penetration test (CPT) measurements. This influence can lead to an incorrect assessment of soil strength/resistance for such purposes as liquefaction triggering analysis. For an accurate measurement of tip and sleeve resistance, unbiased by overburden stress, it is essential to normalize these index measurements appropriately. Presented herein is a comprehensive study reviewing all aspects of CPT normalization. A result of this study is a variable normalization procedure for the CPT that is based on both empirical results and theoretical analysis. This paper presents these results in the form of an improved normalization scheme and discusses its application in practice.     

Logistic Regression Model for Evaluating Soil Liquefaction Probability Using CPT Data

A simple model for evaluating liquefaction probability using cone penetration test (CPT) data is developed based on logistic regression analyses of 396 case histories. The proposed model uses the normalized cone penetration resistance and soil behavior type index as input parameters; therefore, only CPT testing is necessary for evaluating the liquefaction probability of a site. The selection of the model parameters and the expression of equations are based on results of probability examinations and rigorous statistical analyses. Moreover, the derivation of the logistic regression model is presented in a system of equations. The incorporation of these procedures in developing the model not only fully satisfies the statistic requirements but also highlights the physical meanings of the model parameters. Comparisons of the proposed probability model with previously proposed deterministic and probabilistic approaches are performed to demonstrate the improvements. For practical purposes, the developed model is implemented to establish the relationship between the factor of safety against liquefaction and the probability of liquefaction.     

Mapping Liquefaction Potential Considering Spatial Correlations of CPT Measurements

The past studies of liquefaction phenomena during earthquakes have contributed to the development of simplified methods employing field test data to assess the liquefaction potential. Since the field data are limited by exploration cost, it is of interest to obtain valuable and meaningful distribution of liquefaction potential of an area from the limited data. This study proposes a method for assessing liquefaction potential over an extensive area according to the random field concept. The spatial structures of soil properties are estimated from the available cone penetration test (CPT) measurements. The soil properties at unsampled locations are simulated using Monte Carlo simulation. The reliability against liquefaction at every location within the study area is evaluated to map the liquefaction potential. The comparison between simulated distributions of liquefaction potential and observed liquefaction phenomena is discussed. The spatial correlation of soil property provides more information than the traditional approach that solely uses the field test data. The influences of CPT data, penetration locations, and spatial structures of soil properties on the mapping results of liquefaction potential are also discussed.     

Accounting for Soil Aging When Assessing Liquefaction Potential

It has been recognized that liquefaction resistance of sand increases with age due to processes such as cementation at particle contacts and increasing frictional resistance resulting from particle rearrangement and interlocking. As such, the currently available empirical correlations derived from liquefaction of young Holocene sand deposits, and used to determine liquefaction resistance of sand deposits from in situ soil indices [standard penetration test (SPT), cone penetration test (CPT), shear wave velocity test (Vs)], are not applicable for old sand deposits. To overcome this limitation, a methodology was developed to account for the effect of aging on the liquefaction resistance of old sand deposits. The methodology is based upon the currently existing empirical boundary curves for Holocene age soils and utilizes correction factors presented in the literature that comprise the effect of aging on the in situ soil indices as well as on the field cyclic strength (CRR). This paper describes how to combine currently recorded SPT, CPT, and Vs values with corresponding CRR values derived for aged soil deposits to generate new empirical boundary curves for aged soils. The method is illustrated using existing geotechnical data from four sites in the South Carolina Coastal Plain (SCCP) where sand boils associated with prehistoric earthquakes have been found. These sites involve sand deposits that are 200,000?to?450,000?years in age. This work shows that accounting for aging of soils in the SCCP yields less conservative results regarding the current liquefaction potential than when age is not considered. The modified boundary curves indicate that old sand deposits are more resistant to liquefaction than indicated by the existing empirical curves and can be used to evaluate the liquefaction potential at a specific site directly from the current in situ properties of the soil.     

Simplified Cone Penetration Test-based Method for Evaluating Liquefaction Resistance of Soils

This paper presents a new simplified method for assessing the liquefaction resistance of soils based on the cone penetration test (CPT). A relatively large database consisting of CPT measurements and field liquefaction performance observations of historical earthquakes is analyzed. This database is first used to train an artificial neural network for predicting the occurrence and nonoccurrence of liquefaction based on soil and seismic load parameters. The successfully trained and tested neural network is then used to generate a set of artificial data points that collectively define the liquefaction boundary surface, the limit state function. An empirical equation is further obtained by regression analysis to approximate the unknown limit state function. The empirical equation developed represents a deterministic method for assessing liquefaction resistance using the CPT. Based on this newly developed deterministic method, probabilistic analyses of the cases in the database are conducted using the Bayesian mapping function approach. The results of the probabilistic analyses, expressed as a mapping function, provide a simple means for probability-based evaluation of the liquefaction potential. The newly developed simplified method compares favorably to a widely used existing method.     

Liquefaction Potential Map of Charleston, South Carolina Based on the 1886 Earthquake

A liquefaction potential map of the peninsula of Charleston, S.C., is presented in this paper. Liquefaction potential is expressed in terms of the liquefaction potential index developed by Iwasaki et al. and calculated using 44 cone penetration test profiles. The cone profiles are supplemented with information from the 1:24,000 scale geologic map by Weems and Lemon, several first-hand accounts of liquefaction and ground deformation that occurred during the 1886 Charleston earthquake, and liquefaction probabilities determined by Elton and Hadj-Hamou based on standard penetration tests. Nearly all of the cases of liquefaction and ground deformation occurred in the Holocene to late Pleistocene beach deposits that flank the higher-ground sediments of the Wando Formation. To match the observed field behavior, a deposit resistance correction factor of 1.8 is applied to cyclic resistance ratios calculated for the 100,000-year-old Wando Formation. No corrections are needed for the younger deposits. In additional to 1886 field behavior, the deposit resistance corrections are supported by ratios of measured to predicted shear-wave velocity.     

Estimating Liquefaction-Induced Lateral Displacements Using the Standard Penetration Test or Cone Penetration Test

A semiempirical approach to estimate liquefaction-induced lateral displacements using standard penetration test (SPT) or cone penetration test (CPT) data is presented. The approach combines available SPT- and CPT-based methods to evaluate liquefaction potential with laboratory test results for clean sands to estimate the potential maximum cyclic shear strains for saturated sandy soils under seismic loading. A lateral displacement index is then introduced, which is obtained by integrating the maximum cyclic shear strains with depth. Empirical correlations from case history data are proposed between actual lateral displacement, the lateral displacement index, and geometric parameters characterizing ground geometry for gently sloping ground without a free face, level ground with a free face, and gently sloping ground with a free face. The proposed approach can be applied to obtain preliminary estimates of the magnitude of lateral displacements associated with a liquefaction-induced lateral spread.     

Case History Evaluation of Laterally Loaded Piles

This paper examines seven case histories of load tests on piles or drilled shafts under lateral load. Since the current design software to estimate lateral load resistance of deep foundations requires p-y curves. The first approach used was correlative whereby soil parameters determined from in situ tests [standard penetration test (SPT) and cone penetration test (CPT)] were used as input values for standard p-y curves. In the second approach p-y curves were calculated directly from the stress deformation data measured in dilatometer (DMT) and cone pressuremeter tests. The correlative evaluation revealed that, on the average, predictions based upon the SPT were conservative for all loading levels, and using parameters from the CPT best predicted field behavior. Typically, predictions were conservative, except at the maximum load. Since traditionally SPT and CPT correlation-based p-y curves are for “sands” or “clays,” this study suggests that silts, silty sands, and clayey sands should use cohesive p-y curves. For the directly calculated curves, DMT derived p-y curves predict well at low lateral loads, but at higher load levels the predictions become unconservative. p-y curves derived from pressuremeter tests predicted well for both “sands” and “clays” where pore pressures are not anticipated.     

Evaluation of Cone Penetration Test-Based Settlement Prediction Methods for Shallow Foundations on Cohesionless Soils at Highway Bridge Construction Sites

The cone penetration test (CPT) method is gaining popularity in the United States as one of the premier subsurface exploration techniques. The writers previously monitored settlement behavior of spread footing foundations at five highway bridge construction sites in Ohio and evaluated the standard penetration test-based settlement prediction methods in light of the field performance data. The authors recently conducted CPT sounding at these bridge sites. In this paper, settlement behavior of five spread footings resting on cohesionless soils at two sites were predicted with the CPT-based settlement prediction methods proposed by Amar, DeBeer, Meyerhof, and Schmertmann. The results showed that the prediction methods by Schmertmann and DeBeer were both conservative but reasonably reliable in predicting the settlement of shallow foundations observed in the field.     

Liquefaction Potential Index: Field Assessment

Cone penetration test (CPT) soundings at historic liquefaction sites in California were used to evaluate the predictive capability of the liquefaction potential index (LPI), which was defined by Iwasaki et al. in 1978. LPI combines depth, thickness, and factor of safety of liquefiable material inferred from a CPT sounding into a single parameter. LPI data from the Monterey Bay region indicate that the probability of surface manifestations of liquefaction is 58 and 93%, respectively, when LPI equals or exceeds 5 and 15. LPI values also generally correlate with surface effects of liquefaction: Decreasing from a median of 12 for soundings in lateral spreads to 0 for soundings where no surface effects were reported. The index is particularly promising for probabilistic liquefaction hazard mapping where it may be a useful parameter for characterizing the liquefaction potential of geologic units.     

Standard Penetration Test-Based Probabilistic and Deterministic Assessment of Seismic Soil Liquefaction Potential

This paper presents new correlations for assessment of the likelihood of initiation (or “triggering”) of soil liquefaction. These new correlations eliminate several sources of bias intrinsic to previous, similar correlations, and provide greatly reduced overall uncertainty and variance. Key elements in the development of these new correlations are (1) accumulation of a significantly expanded database of field performance case histories; (2) use of improved knowledge and understanding of factors affecting interpretation of standard penetration test data; (3) incorporation of improved understanding of factors affecting site-specific earthquake ground motions (including directivity effects, site-specific response, etc.); (4) use of improved methods for assessment of in situ cyclic shear stress ratio; (5) screening of field data case histories on a quality/uncertainty basis; and (6) use of high-order probabilistic tools (Bayesian updating). The resulting relationships not only provide greatly reduced uncertainty, they also help to resolve a number of corollary issues that have long been difficult and controversial including: (1) magnitude-correlated duration weighting factors, (2) adjustments for fines content, and (3) corrections for overburden stress.     

Analyzing Liquefaction-Induced Lateral Spreads Using Strength Ratios

The writers backanalyzed 39 well-documented liquefaction-induced lateral spreads in terms of a mobilized strength ratio, su(mob)/σvo′ using the Newmark sliding block method. Based on the inverse analyses results, we found that the backcalculated strength ratios mobilized during lateral spreads can be directly correlated to normalized cone penetration test tip resistance and standard penetration test blow count. Remarkably, Newmark analysis-based strength ratios mobilized during these lateral spreads essentially coincide with liquefied strength ratios backcalculated from liquefaction flow failures. The mobilized strength ratios appear to be independent of the magnitude of lateral displacement (at least for displacements greater than 15?cm) and the strength of shaking (in terms of peak ground acceleration). Furthermore, the mobilized strength ratios backcalculated from these cases appear to be consistent for a given depositional environment and do not appear to be severely impacted by potential water layer formation.     

Liquefaction Resistance of Clean and Nonplastic Silty Sands Based on Cone Penetration Resistance

Liquefaction of granular soil deposits is one of the major causes of loss resulting from earthquakes. The accuracy in the assessment of the likelihood of liquefaction at a site affects the safety and economy of the design. In this paper, curves of cyclic resistance ratio (CRR) versus cone penetration test (CPT) stress-normalized cone resistance qc1 are developed from a combination of analysis and laboratory testing. The approach consists of two steps: (1) determination of the CRR as a function of relative density from cyclic triaxial tests performed on samples isotropically consolidated to 100 kPa; and (2) estimation of the stress-normalized cone resistance qc1 for the relative densities at which the soil liquefaction tests were performed. A well-tested penetration resistance analysis based on cavity expansion analysis was used to calculate qc1 for the various soil densities. A set of 64 cyclic triaxial tests were performed on specimens of Ottawa sand with nonplastic silt content in the range of 0–15% by weight, and relative densities from loose to dense for each gradation, to establish the relationship of the CRR to the soil state and fines content. The resulting (CRR)7.5-qc1 relationship for clean sand is consistent with widely accepted empirical relationships. The (CRR)7.5-qc1 relationships for the silty sands depend on the relative effect of silt content on the CRR and qc1. It is shown that the cone resistance increases at a higher rate with increasing silt content than does liquefaction resistance, shifting the (CRR)7.5-qc1 curves to the right. The (CRR)7.5-qc1 curves proposed for both clean and silty sands are consistent with field observations.     

Return Period of Soil Liquefaction

The paper describes a performance-based approach to the evaluation of liquefaction potential, and shows how it can be used to account for the entire range of potential ground shaking. The result is a direct estimate of the return period of liquefaction, rather than a factor of safety or probability of liquefaction conditional upon ground shaking with some specified return period. As such, the performance-based approach can be considered to produce a more complete and consistent indication of the actual likelihood of liquefaction at a given location than conventional procedures. In this paper, the performance-based procedure is introduced and used to compare likelihoods of the initiation of liquefaction at identical sites located in areas of different seismicity. The results indicate that the likelihood of liquefaction depends on the position and slope of the peak acceleration hazard curve, and on the distribution of earthquake magnitudes contributing to the ground motion hazard. The results also show that the consistent use of conventional procedures for the evaluation of liquefaction potential produces inconsistent actual likelihoods of liquefaction.