The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.
The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
Capacity design and hierarchy of strength philosophies at the base of modern seismic codes allow inelastic response in case of severe earthquakes and thus, in most traditional systems, damage develops at well-defined locations of reinforced concrete (RC) structures, known as plastic hinges. The 2010 and 2011 Christchurch earthquakes have demonstrated that this philosophy worked as expected. Plastic hinges formed in beams, in coupling beams and at the base of columns and walls. Structures were damaged permanently, but did not collapse. The 2010 and 2011 Christchurch earthquakes also highlighted a critical issue: the reparability of damaged buildings. No methodologies or techniques were available to estimate the level of subsequent earthquakes that RC buildings could still sustain before collapse. No repair techniques capable of restoring the initial condition of buildings were known. Finally, the cost-effectiveness of an eventual repair intervention, when compared with a new building, was unknown. These aspects, added to nuances of New Zealand building owners’ insurance coverage, encouraged the demolition of many buildings. Moreover, there was a perceived strong demand from government and industry to develop techniques for assessing damage to steel reinforcement bars embedded in cracked structural concrete elements. The most common questions were: “Have the steel bars been damaged in correspondence to the concrete cracks?”, “How much plastic deformation have the steel bars undergone?”, and “What is the residual strain capacity of the damaged bars?” Minimally invasive techniques capable of quantifying the level and extent of plastic deformation and residual strain capacity are not yet available. Although some studies had been recently conducted, a validated method is yet to be widely accepted. In this thesis, a least-invasive method for the damage-assessment of steel reinforcement is developed. Based on the information obtained from hardness testing and a single tensile test, it is possible to estimate the mechanical properties of earthquake-damaged rebars. The reduction in the low-cycle fatigue life due to strain ageing is also quantified. The proposed damage assessment methodology is based on empirical relationships between hardness and strain and residual strain capacity. If damage is suspected from in situ measurements, visual inspection or computer analysis, a bar may be removed and more accurate hardness measurements can be obtained using the lab-based Vickers hardness methodology. The Vickers hardness profile of damaged bars is then compared with calibration curves (Vickers hardness versus strain and residual strain capacity) previously developed for similar steel reinforcement bars extracted from undamaged locations. Experimental tests demonstrated that the time- and temperature-dependent strain-ageing phenomenon causes changes in the mechanical properties of plastically deformed steels. In particular, yield strength and hardness increases, whereas ductility decreases. The changes in mechanical properties are quantified and their implications on the hardness method are highlighted. Low-cycle fatigue (LCF) failures of steel reinforcing bars have been observed in laboratory testing and post-earthquake damage inspections. Often, failure might not occur during a first seismic event. However, damage is accumulated and the remaining fatigue life is reduced. Failure might therefore occur in a subsequent seismic event. Although numerous studies exist on the LCF behaviour of steel rebars, no studies had been conducted on the strain-ageing effects on the remaining fatigue life. In this thesis, the reduction in fatigue life due to this phenomenon is determined through a number of experimental tests.
Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults. The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains. The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future. Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case. This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features. Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
