The magnitude Mw 6.2 earthquake of February 22nd 2011 that struck beneath the city of Christchurch, New Zealand, caused widespread damage and was particularly destructive to the Central Business District (CBD). The shaking caused major damage, including collapses of structures, and initiated ground failure in the form of soil liquefaction and consequent effects such as sand boils, surface flooding, large differential settlements of buildings and lateral spreading of ground towards rivers were observed. A research project underway at the University of Canterbury to characterise the engineering behaviour of the soils in the region was influenced by this event to focus on the performance of the highly variable ground conditions in the CBD. This paper outlines the methodology of this research to characterise the key soil horizons that underlie the CBD that influenced the performance of important structures during the recent earthquakes, and will influence the performance of the rebuilt city centre under future events. The methodology follows post-earthquake reconnaissance in the central city, a desk study on ground conditions, site selection, mobilisation of a post-earthquake ground investigation incorporating the cone penetration test (CPT), borehole drilling, shear wave velocity profiling and Gel-push sampling followed by a programme of laboratory testing including monotonic and cyclic testing of the soils obtained in the investigation. The research is timely and aims to inform the impending rebuild, with appropriate information on the soils response to dynamic loading, and the influence this has on the performance of structures with various foundation forms.
This report presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region during the Canterbury earthquake sequence. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the shallow geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site is based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS to better understand stratigraphy and soil properties, which are important in seismic site response. In this regard, CPT, SPT and borehole data, shear wave velocity (Vs) profiles, and horizontal to vertical spectral ratio measurements (H/V) in close vicinity to the SMS were used to develop representative soil profiles at each site. NZS1170.5 (SNZ 2004) site subsoil classifications were updated using Vs and SPT N60 criteria. Site class E boundaries were treated as a sliding scale rather than as a discrete boundary to account for locations with similar site effects potential, an approach which was shown to result in a better delineation between the site classes. SPT N60 values often indicate a stiffer site class than the Vs data for softer soil sites, highlighting the disparity between the two site investigation techniques. Both SPT N60 and Vs based site classes did not always agree with the original site classifications. This emphasises the importance of having detailed site‐specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N60 and Vs. Liquefaction triggering assessments were carried out for the Darfield and Christchurch earthquakes, and compared against observed liquefaction surface manifestations and ground motions characteristics at each SMS. In general, the characteristics of the recorded ground motions at each site correlate well with the triggering analyses. However, at sites that likely liquefied at depth (as indicated by triggering analyses and/or inferred from the characteristics of the recorded surface acceleration time series), the presence of a non‐liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects.
Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.
A photograph of the Townsend Telescope in the Observatory at the Christchurch Arts Centre. Graeme Kershaw, Technician at the University of Canterbury Department of Astronomy and Physics, is standing to the left. This image was used by Kershaw to identify the telescope's parts after the 22 February 2011 earthquake.
A photograph of the clock drive of the Townsend Telescope. The telescope is in the Observatory at the Christchurch Arts Centre. This image was used by Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, to identify the telescope's parts after the 22 February 2011 earthquake.
A photograph of the middle section of the Townsend Telescope. The telescope is in the Observatory at the Christchurch Arts Centre. This image was used by Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, to identify the telescope's parts after the 22 February 2011 earthquake.
A scanned copy of a black and white photograph belonging to University of Canterbury alumnus Colin Lau. Colin describes the photograph as follows: "Our science building (where I got my physics degree in 1972) on the left, lecture hall in the middle and engineering school to the right & behind the lecture hall".
The Townsend Observatory is located in the Arts Centre of Christchurch, in what used to be Canterbury College (now University of Canterbury). The Townsend telescope itself is a historic 6-inch Cooke refractor built in 1864 for early Christchurch colonist, Mr James Townsend, and gifted by him to Christchurch College in 1891. At the same time, the Canterbury Astronomical Society handed over its funds to the College to help erect an observatory. The College used this, and money it had set aside for a medical school, to build a biological laboratory with an attached observatory tower, which was completed in 1896. The Biology Building and Observatory Tower was the last major design by architect Benjamin Mountfort. Mr Walter Kitson was appointed custodian of the telescope and regular public open nights commenced. and continued until 2010, with the telescope being operated by students of the Department of Physics and Astronomy, University of Canterbury. The Observatory Tower was badly damaged in the 4 September 2010 earthquake and collapsed in the 22 February 2011 earthquake. The telescope was badly damaged by the collapse, but, amazingly, the optics were found entirely intact. The Department of Physics and Astronomy plans to restore the Townsend Telescope so that it can be returned to a replica Observatory Tower in its central city home, enabling the people of Christchurch, and visitors, to enjoy views of the night sky through this beautiful and historic telescope once again.
A scanned copy of a black and white photograph of the argon gas laser used by David Lockwood in his PhD research at the University of Canterbury. David says, "In contrast to the He-Ne laser I developed for my MSc thesis, the laser I used was one of the first commercial lasers - an argon gas laser that required frequent painstaking maintenance. This was because the He-Ne gas laser did not produce enough power for my experiments".
A scanned copy of a photograph of the He-Ne Continuous Gas Laser used in David Lockwood's MSc research at the University of Canterbury. David explains that the photograph shows "the complete experimental system, including the horn-shaped acoustic chamber in front and the operating laser in the rear".
A scanned copy of a photograph of an optical diffraction pattern produced by the He-Ne Continuous Gas Laser used in David Lockwood's MSc research at the University of Canterbury. David explains that the photograph shows "a typical optical diffraction pattern obtained from the original red laser beam arising from a grating structure formed by the alignment of a colloid under the influence of a travelling sound wave".
A scanned copy of a black and white photograph of the He-Ne Continuous Gas Laser used in David Lockwood's MSc research at the University of Canterbury. David says, "The experiments in the Thesis were aimed at examining the effects of sound waves on colloids, which are very small particles that are found in clays for example (for more details see page 181 of http://www.cap.ca/PiC-PaC/static/downloads/1efdc1f3784b85c1a1b33e396b91ee8aef2072c5.pdf ). In fact, I was looking at nanoparticles of matter, which preceded the emergence of nanotechnology as a field of research by more than two decades. I needed a suitable light source to probe the alignment of the nanoparticles in the ultrasonic sound field. At that time the laser had just been invented and I realized that this was the ideal light source for my experiment. I then proceeded to construct the first home-built laser in New Zealand. This laser - a He-Ne continuous-wave gas laser - operated at 632.8 nm (in the red). This laser, over a metre long, is shown in the black-and-white photo. There were lot of problems to be overcome, but eventually, with invaluable help from Dr. Tom Seed (my MSc supervisor) and Dr. Rod Syme, I had it working".
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.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
A photograph of the Townsend Telescope. The telescope is in the Observatory at the Christchurch Arts Centre. A plate on the side reads, "T Cook & Sons, 1864, York & London". This image was used by Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, to identify the telescope's parts after the 22 February 2011 earthquake.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, scratching his head as he looks at the damaged clock drive of the Townsend Telescope. The telescope was damaged during the 22 February 2011 earthquake, when the Observatory tower at the Christchurch Arts Centre collapsed. Kershaw has been given the task of restoring the telescope.
A photograph of the middle section of the Townsend Telescope. The telescope is in the Observatory at the Christchurch Arts Centre. A plate on the side reads, "T Cook & Sons, 1864, York & London". This image was used by Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, to identify the telescope's parts after the 22 February 2011 earthquake.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, looking through the lens of the Townsend Telescope. The lens is the most crucial component of the telescope. It was found unscratched and unbroken in the rubble of the Observatory tower, after 22 February 2011 earthquake. This means that the telescope can be restored.
A photograph of the middle section of the Townsend Telescope. The telescope is in the Observatory at the Christchurch Arts Centre. A plate on the side reads, "T Cook & Sons, 1864, York & London". This image was used by Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, to identify the telescope's parts after the 22 February 2011 earthquake.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, sitting among the damaged parts of the Townsend Telescope. The telescope was damaged during the 22 February 2011 earthquake, when the Observatory tower at the Christchurch Arts Centre collapsed. Kershaw has been given the task of restoring the telescope. In this photograph he is holding the telescope's clock drive.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, sitting among the damaged parts of the Townsend Telescope. The telescope was damaged during the 22 February 2011 earthquake, when the Observatory tower at the Christchurch Arts Centre collapsed. Kershaw has been given the task of restoring the telescope. In this photograph he is holding the telescope's clock drive.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, sitting among the damaged parts of the Townsend Telescope. The telescope was damaged during the 22 February 2011 earthquake, when the Observatory tower at the Christchurch Arts Centre collapsed. Kershaw has been given the task of restoring the telescope. In this photograph he is holding the telescope's clock drive.
A photograph of Graeme Kershaw, Technician at the University of Canterbury Department of Physics and Astronomy, sitting among the damaged parts of the Townsend Telescope. The telescope was damaged during the 22 February 2011 earthquake, when the Observatory tower at the Christchurch Arts Centre collapsed. Kershaw has been given the task of restoring the telescope. In this photograph he is holding the telescope's clock drive.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
A photograph of the lower end of the main tube of the Townsend Telescope. The tube was crushed and bent during the 22 February 2011 earthquake.
A photograph of some of the parts of the Townsend Telescope. Many of the parts were damaged during the 22 February 2011 earthquake.
A photograph of two clock gears from the Townsend Telescope.
A photograph of the scale viewer eyepiece from the Townsend Telescope.
A photograph of the damage to the teeth of a gear from the worm gear drive of the Townsend Telescope. The gear was damaged during the 22 February 2011 earthquake.