What Is Liquefaction?
Liquefaction is a phenomenon in which saturated, loosely packed granular soils temporarily lose shear strength and behave like a viscous fluid during earthquake shaking. Seismic waves generate cyclic shear stresses that cause pore water pressure to build up in undrained conditions. When the excess pore pressure equals the initial effective confining stress, effective stress drops to zero: σ'v = σv − u = 0. With no effective stress, frictional strength vanishes and the soil can no longer support loads.
Manifestations include: sand boils, ground settlement, lateral spreading, foundation bearing failure, and flotation of buried structures. The 1964 Niigata earthquake and the 1964 Alaska earthquake brought liquefaction to international engineering attention; Loma Prieta (1989) and Christchurch (2010–2011) provided extensive modern case histories.
Susceptibility Conditions
Three conditions must be met for liquefaction to occur:
- Loose to medium-dense granular soil: (N1)60 < 30 (SPT) or qc1N < 160 (CPT). Well-graded gravels are less susceptible due to high permeability allowing rapid drainage; fine silty sands with FC > 35% may be non-liquefiable by the simplified procedure but can still exhibit cyclic softening.
- Saturated conditions: soil must be below the water table, or at least partially saturated with pore pressure that can build up.
- Sufficient seismic shaking: peak ground acceleration (PGA) typically > 0.10–0.15g and magnitude Mw ≥ 5.5.
Cyclic Stress Ratio (CSR) — Earthquake Demand
The earthquake-induced cyclic shear stress normalized by effective overburden stress (Seed and Idriss, 1971):
CSR = 0.65 · (σv / σ'v) · (amax / g) · rd
where σv = total vertical stress (kPa), σ'v = effective vertical stress (kPa), amax = peak ground surface acceleration (g), and rd = stress reduction factor accounting for flexibility of the soil column:
- rd = 1.0 − 0.00765z for z ≤ 9.15 m
- rd = 1.174 − 0.0267z for 9.15 m < z ≤ 23 m
- rd = 0.744 − 0.008z for z > 23 m (approximate)
For Mw ≠ 7.5, scale CSR by the magnitude scaling factor (MSF): CSR_M7.5 = CSR / MSF, where MSF = 10^(2.24) / (Mw^2.56). This normalizes demand to the Mw = 7.5 reference earthquake.
Cyclic Resistance Ratio (CRR) — Soil Capacity
CRR is the soil's ability to resist liquefaction, correlated to (N1)60cs (clean-sand-corrected SPT blow count):
CRR_7.5 = 1/(34 − (N1)60cs) + (N1)60cs/135 + 50/(10·(N1)60cs + 45)² − 1/200
This is the Youd et al. (2001) update of the Seed-Idriss curve, valid for (N1)60cs < 30. At (N1)60cs = 30, CRR → ∞ (non-liquefiable).
Fines content correction: for FC > 5%, (N1)60cs = α + β·(N1)60 where α = 0 (FC ≤ 5%), 5.0 (FC = 15%), 9.7 (FC = 35%); β = 1.0, 1.2, 1.5 respectively.
CPT-based CRR uses Robertson and Wride (1998): qc1N (normalized tip resistance) and CRR from the CPT clean sand base curve. CPT is often preferred because it provides a continuous resistance profile.
Factor of Safety Against Liquefaction
FS_liq = CRR_7.5 / CSR_M7.5
- FS_liq > 1.3: liquefaction unlikely
- 1.0 < FS_liq ≤ 1.3: marginal; further investigation required (cyclic softening possible)
- FS_liq ≤ 1.0: liquefaction likely
Perform the calculation at multiple depths within the saturated granular zone. The worst-case profile governs triggering assessment.
IBC/ASCE 7 Site Classification and Liquefaction Hazard
ASCE 7 Section 20.3 requires liquefaction assessment when PGA ≥ 0.15g (Seismic Design Category C and above). Site Class F includes sites requiring site-specific hazard analysis, including: soils with FS_liq < 1.0 when PGA > 0.10g, sites with peats or sensitive clays, and sites with very thick soft clay (Su < 25 kPa). IBC Section 1803.5.12 requires geotechnical investigation to evaluate liquefaction potential for Seismic Design Categories C, D, E, and F.
Post-Liquefaction Settlements and Lateral Spreading
Settlement: Tokimatsu and Seed (1987) or Zhang, Robertson, and Brachman (2002) provide volumetric strain εv as a function of FS_liq and (N1)60. Settlement of each layer: Sv = εv × H_layer. Sum over all liquefiable layers. Post-liquefaction settlements of 100–600 mm are common in moderately liquefiable sites.
Lateral spreading: horizontal displacement of soil toward a free face (river bank, slope, excavation) following liquefaction. Bartlett and Youd (1992) empirical model estimates displacements as a function of slope gradient, free face ratio, layer thickness, and grain size. Displacements of 1–5 m are possible in severe cases. Lateral spreading is among the most destructive consequences of liquefaction for buried utilities, bridge abutments, and port facilities.
Flow failure: occurs when post-liquefaction residual strength (Sr) is insufficient to maintain slope stability. Sr from SPT back-analyses (Idriss and Boulanger, 2007): Sr ≈ exp{(N1)60cs/16 − 3.0} in kPa terms (bounded by uncertainty envelope).
Mitigation Strategies
- Densification: vibro-compaction (for clean sands), vibro-replacement/stone columns (for silty sands), dynamic compaction. Increases (N1)60 above liquefaction threshold. Most effective when FC < 15%.
- Permeability enhancement: gravel drains or wick drains reduce pore pressure buildup during shaking by accelerating drainage. Effective when shaking duration is short.
- Grouting: permeation grouting (chemical grout fills voids) for fine sands; compaction grouting (expanding grout bulb densifies soil). Used in urban areas where vibration from densification methods is unacceptable.
- Deep foundations: bypass liquefiable layers with piles extending to firm stratum; design for downdrag from settling liquefied layers and lateral loads from lateral spreading.
- Soil mixing (DSM): deep soil mixing creates cement-soil columns that provide structure even if pore pressures build up. Effective grid patterns can confine liquefiable soil between treated columns.
Ground improvement selection depends on soil type, site access, treatment depth, proximity to structures, and cost. FHWA Geotechnical Engineering Circular 13 (GEC 13) provides current guidance on ground improvement for liquefaction mitigation.