The majority of current procedures used to deduce liquefaction potential of soils rely on empirical methods. These methods have been proven to work in the past, but these methods are known to overestimate the liquefaction potential in certain regions of Christchurch due to a whole range of factors, and the theoretical basis behind these methods cannot be explained scientifically. Critical state soil mechanics theory was chosen to provide an explanation for the soil's behaviour during the undrained shearing. Soils from two sites in Christchurch were characterised at regular intervals for the critical layers and tested for the critical state lines (CSL). Various models and relationships were then used to predict the CSL and compared with the actual CSL. However none of the methods used managed to predict the CSL accurately, and a separate Christchurch exclusive relationship was proposed. The resultant state parameter values could be obtained from shear-wave velocity plots and were then developed into cyclic resistance ratios (CRR). These were subsequently compared with cyclic stress ratios (CSR) from recent Christchurch earthquakes to obtain the factor of safety. This CSL-based approach was compared with other empirical methods and was shown to yield a favourable relationship with visual observations at the sites' locations following the earthquake.
The seismic performance of soil profiles with potentially liquefiable deposits is a complex phenomenon that requires a thorough understanding of the soil properties and ground motion characteristics. The limitations of simplified liquefaction assessment methods have prompted an increase in the use of non-linear dynamic analysis methods. Focusing on onedimensional site response of a soil column, this thesis validated a soil constitutive model using in-situ pore pressure measurements and then assessed the influence of input ground motion characteristics on soil column response using traditional and newly developed metrics. Pore pressure recordings during the Canterbury Earthquake Sequence (CES) in New Zealand were used to validate the PM4Sand constitutive model. Soil profile characterization was key to accurate prediction of excess pore pressure response and accounting for any densification during the CES. Response during multiple earthquakes was captured effectively and cross-layer interaction demonstrated the model capability to capture soil response at the system-level. Synthetic and observed ground motions from the Christchurch earthquake were applied to the validated soil column to quantify the performance of synthetic motions. New metrics were developed to facilitate a robust comparison to assess performance. The synthetic input motions demonstrated a slightly larger acceleration and excess pore pressure response compared to the observed input motions. The results suggest that the synthetic motions may accumulate higher excess pore pressure at a faster rate and with fewer number of cycles in the shear response. This research compares validated soil profile subject to spectrally-matched pulse and non-pulse motions, emphasizing the inclusion of pulse motions with distinctive characteristics in ground motion suites for non-linear dynamic analysis. However, spectral matching may lead to undesired alterations in pulse characteristics. Cumulative absolute velocity and significant duration significantly differed between these two groups compared to the other key characteristics and contributed considerably to the liquefaction response. Unlike the non-pulse motions, not all of the pulse motions triggered liquefaction, likely due to their shorter significant duration. Non-pulse motions developed a greater spatial extent of liquefaction triggering in the soil profile and extended to a greater depth.
