Search

found 212 results

Research papers, University of Canterbury Library

The city of Christchurch and its surrounds experienced widespread damage due to soil liquefaction induced by seismic shaking during the Canterbury earthquake sequence that began in September 2010 with the Mw7.1 Darfield earthquake. Prior to the start of this sequence, the city had a large network of strong motion stations (SMSs) installed, which were able to record a vast database of strong ground motions. This paper uses this database of strong ground motion recordings, observations of liquefaction manifestation at the ground surface, and data from a recently completed extensive geotechnical site investigation program at each SMS to assess a range of liquefaction evaluation procedures at the four SMSs in the Christchurch Central Business District (CBD). In general, the characteristics of the accelerograms recorded at each SMS correlated well with the liquefaction evaluation procedures, with low liquefaction factors of safety predicted at sites with clear liquefaction identifiers in the ground motions. However, at sites that likely liquefied at depth (as indicated by evaluation procedures and/or inferred from the characteristics of the recorded surface accelerograms), the presence of a non-liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects. Because of this, there was not a good correlation between surface manifestation and two surface manifestation indices, the Liquefaction Potential Index (LPI) and the Liquefaction Severity Number (LSN).

Research papers, University of Canterbury Library

Our poster will present on-going QuakeCoRE-founded work on strong motion seismology for Dunedin-Mosgiel area, focusing on ground motion simulations for Dunedin Central Business District (CBD). Source modelling and ground motion simulations are being carried out using the SCEC (Southern California Earthquakes Center) Broad Band simulation Platform (BBP). The platform computes broadband (0-10 Hz) seismograms for earthquakes and was first implemented at the University of Otago in 2016. As large earthquakes has not been experienced in Dunedin in the time of period of instrumental recording, user-specified scenario simulations are of great value. The Akatore Fault, the most active fault in Otago and closest major fault to Dunedin, is the source focused on in the present study. Simulations for various Akatore Fault source scenarios are run and presented. Path and site effects are key components considered in the simulation process. A 1D shear wave velocity profile is required by SCEC BBP, and this is being generated to represent the Akatore-to-CBD path and site within the BBP. A 3D shear velocity model, with high resolution within Dunedin CBD, is being developed in parallel with this study (see Sangster et al. poster). This model will be the basis for developing a 3D shear wave velocity model for greater Dunedin-Mosgiel area for future ground motion simulations, using Canterbury software (currently under development).

Research papers, University of Canterbury Library

This study provides an initial examination of source parameter uncertainty in a New Zealand ground motion simulation model, by simulating multiple event realisations with perturbed source parameters. Small magnitude events in Canterbury have been selected for this study due to the small number of source input parameters, the wealth of recorded data, and the lack of appreciable off-fault non-linear effects. Which provides greater opportunity to identify systematic source, path and site effects, required to robustly investigate the causes of uncertainty.

Research papers, University of Canterbury Library

Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.

Research papers, University of Canterbury Library

Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.

Research papers, University of Canterbury Library

We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows:  Ic < 1.31 – Gravelly sand to dense sand  1.31 < Ic < 1.90 – Sands: clean sand to silty sand  1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt  2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay  3.20 < Ic < 3.60 – Clays: silty clay to clay  Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.

Research papers, University of Canterbury Library

In this article we utilize grounded theory to explore women’s experiences in the unique construction industry context that followed the 2010 Canterbury (New Zealand) earthquakes. Data were obtained from 36 semi-structured interviews conducted with women working in a variety of occupations in the construction industry. We identify three inter-related categories: capitalizing on opportunity, demonstrating capability and surface tolerance, which together represent a response process that we label ‘deferential tailoring’. The deferential tailoring process explains how women intentionally shape their response to industry conditions through self-regulating behaviors that enables them to successfully seize opportunities and manage gender-related challenges in the working environment. Our findings challenge existing research which suggests that women adopt submissive coping strategies to conform to androcentric norms in the construction industry. Instead, we argue that the process of deferential tailoring can empower women to build positive workplace relationships, enhance career development, and help shift perceptions of the value of their work in the industry.

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.

Research papers, University of Canterbury Library

Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.

Research papers, University of Canterbury Library

This paper provides a photographic tour of the ground-surface rupture features of the Greendale Fault, formed during the 4th September 2010 Darfield Earthquake. The fault, previously unknown, produced at least 29.5 km of strike-slip surface deformation of right-lateral (dextral) sense. Deformation, spread over a zone between 30 and 300 m wide, consisted mostly of horizontal flexure with subsidiary discrete shears, the latter only prominent where overall displacement across the zone exceeded about 1.5 m. A remarkable feature of this event was its location in an intensively farmed landscape, where a multitude of straight markers, such as fences, roads and ditches, allowed precise measurements of offsets, and permitted well-defined limits to be placed on the length and widths of the surface rupture deformation.

Research papers, University of Canterbury Library

SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.

Research papers, University of Canterbury Library

Overview of SeisFinder SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database. The basic SeisFinder architecture is shown in Figure 1.

Research papers, University of Canterbury Library

The purpose of this thesis is to evaluate the seismic response of the UC Physics Building based on recorded ground motions during the Canterbury earthquakes, and to use the recorded response to evaluate the efficacy of various conventional structural analysis modelling assumptions. The recorded instrument data is examined and analysed to determine how the UC Physics Building performed during the earthquake-induced ground motions. Ten of the largest earthquake events from the 2010-11 Canterbury earthquake sequence are selected in order to understand the seismic response under various levels of demand. Peak response amplitude values are found which characterise the demand from each event. Spectral analysis techniques are utilised to find the natural periods of the structure in each orthogonal direction. Significant torsional and rocking responses are also identified from the recorded ground motions. In addition, the observed building response is used to scrutinise the adequacy of NZ design code prescriptions for fundamental period, response spectra, floor acceleration and effective member stiffness. The efficacy of conventional numerical modelling assumptions for representing the UC Physics Building are examined using the observed building response. The numerical models comprise of the following: a one dimensional multi degree of freedom model, a two dimensional model along each axis of the building and a three dimensional model. Both moderate and strong ground motion records are used to examine the response and subsequently clarify the importance of linear and non-linear responses and the inclusion of base flexibility. The effects of soil-structure interaction are found to be significant in the transverse direction but not the longitudinal direction. Non-linear models predict minor in-elastic behaviour in both directions during the 4 September 2010 Mw 7.1 Darfield earthquake. The observed torsional response is found to be accurately captured by the three dimensional model by considering the interaction between the UC Physics Building and the adjacent structure. With the inclusion of adequate numerical modelling assumptions, the structural response is able to be predicted to within 10% for the majority of the earthquake events considered.

Research papers, University of Canterbury Library

Damage distribution maps from strong earthquakes and recorded data from field experiments have repeatedly shown that the ground surface topography and subsurface stratigraphy play a decisive role in shaping the ground motion characteristics at a site. Published theoretical studies qualitatively agree with observations from past seismic events and experiments; quantitatively, however, they systematically underestimate the absolute level of topographic amplification up to an order of magnitude or more in some cases. We have hypothesized in previous work that this discrepancy stems from idealizations of the geometry, material properties, and incident motion characteristics that most theoretical studies make. In this study, we perform numerical simulations of seismic wave propagation in heterogeneous media with arbitrary ground surface geometry, and compare results with high quality field recordings from a site with strong surface topography. Our goal is to explore whether high-fidelity simulations and realistic numerical models can – contrary to theoretical models – capture quantitatively the frequency and amplitude characteristics of topographic effects. For validation, we use field data from a linear array of nine portable seismometers that we deployed on Mount Pleasant and Heathcote Valley, Christchurch, New Zealand, and we compute empirical standard spectral ratios (SSR) and single-station horizontal-to-vertical spectral ratios (HVSR). The instruments recorded ambient vibrations and remote earthquakes for a period of two months (March-April 2017). We next perform two-dimensional wave propagation simulations using the explicit finite difference code FLAC. We construct our numerical model using a high-resolution (8m) Digital Elevation Map (DEM) available for the site, an estimated subsurface stratigraphy consistent with the geomorphology of the site, and soil properties estimated from in-situ and non-destructive tests. We subject the model to in-plane and out-of-plane incident motions that span a broadband frequency range (0.1-20Hz). Numerical and empirical spectral ratios from our blind prediction are found in very good quantitative agreement for stations on the slope of Mount Pleasant and on the surface of Heathcote Valley, across a wide range of frequencies that reveal the role of topography, soil amplification and basin edge focusing on the distribution of ground surface motion.

Research papers, University of Canterbury Library

The south Leader Fault (SLF) is a newly documented active structure that ruptured the surface during the Mw 7.8 Kaikoura earthquake. The Leader Fault is a NNE trending oblique left lateral thrust that links the predominantly right lateral ‘The Humps’ and Conway-Charwell faults. The present research uses LiDAR at 0.5 m resolution and field mapping to determine the factors controlling the surface geometries and kinematics of the south Leader Fault ruptures at the ground surface. The SLF zone is up to 2km wide and comprises a series of echelon NE-striking thrusts linked by near-vertical N-S striking faults. The thrusts are upthrown to the west by up to 1 m and dip 35-45°. Thrust slip surfaces are parallel with Cretaceous-Cenozoic bedding and may reflect flexural slip folding. By contrast, the northerly striking faults dip steeply (65° west- 85° east), and accommodate up to 3m of oblique left lateral displacement at the ground surface and displace Cenozoic bedding. Some of the SLF has been mapped in bedrock, although none were known to be active prior to the earthquake or have a strong topographic expression. The complexity of fault rupture and the width of the fault zone appears to reflect the occurrence of faulting and folding at the ground surface during the earthquake.

Research papers, University of Canterbury Library

Seismic behaviour of typical unreinforced masonry (URM) brick houses, that were common in early last century in New Zealand and still common in many developing countries, is experimentally investigated at University of Canterbury, New Zealand in this research. A one halfscale model URM house is constructed and tested under earthquake ground motions on a shaking table. The model structure with aspect ratio of 1.5:1 in plan was initially tested in the longitudinal direction for several earthquakes with peak ground acceleration (PGA) up to 0.5g. Toppling of end gables (above the eaves line) and minor to moderate cracking around window and door piers was observed in this phase. The structure was then rotated 90º and tested in the transverse (short) direction for ground motions with PGA up to 0.8g. Partial out-of-plane failure of the face loaded walls in the second storey and global rocking of the model was observed in this phase. A finite element analysis and a mechanism analysis are conducted to assess the dynamic properties and lateral strength of the model house. Seismic fragility function of URM houses is developed based on the experimental results. Damping at different phases of the response is estimated using an amplitude dependent equivalent viscous damping model. Financial risk of similar URM houses is then estimated in term of expected annual loss (EAL) following a probabilistic financial risk assessment framework. Risks posed by different levels of damage and by earthquakes of different frequencies are then examined.

Research papers, University of Canterbury Library

The 4 September 2010 Darfield and 22 February 2011 Christchurch earthquakes caused significant damage to Christchurch and surrounding suburbs as a result of the widespread liquefaction and lateral spreading that occurred. Ground surveying-based field investigations were conducted following these two events in order to measure permanent ground displacements in areas significantly affected by lateral spreading. Data was analysed with respect to the distribution of lateral spreading vs. distance from the waterway, and the failure patterns observed. Two types of failure distribution patterns were observed, a typical distributed pattern and an atypical block failure. Differences in lateral spreading measurements along adjacent banks of the Avon River in the area of Dallington were also examined. The spreading patterns between the adjacent banks varied with the respective river geometry and/or geotechnical conditions at the banks.

Research papers, University of Canterbury Library

essential systems upon which the well-being and functioning of societies depend. They deliver a service or a good to the population using a network, a combination of spatially-distributed links and nodes. As they are interconnected, network elements’ functionality is also interdependent. In case of a failure of one component, many others could be momentarily brought out-of-service. Further problems arise for buried infrastructure when it comes to buried infrastructure in earthquake and liquefaction-prone areas for the following reasons: • Technically more demanding inspections than those required for surface horizontal infrastructure • Infrastructure subject to both permanent ground displacement and transient ground deformation • Increase in network maintenance costs (i.e. deterioration due to ageing material and seismic hazard) These challenges suggest careful studies on network resilience will yield significant benefits. For these reasons, the potable water network of Christchurch city (Figure 1) has been selected for its well-characterized topology and its extensive repair dataset.

Research papers, University of Canterbury Library

This work investigates the possibility of developing a non-contact, non-line of sight sensor to measure interstorey drift through simulation and experimental validation. • The method uses frequency-modulated continuous wave (FMCW) radar to measure displacement. This method is commonly in use in a number of modern applications, including aircraft altimeters and automotive parking sensors. • The technique avoids numerous problems found in contemporary structural health monitoring methods, namely integral drift errors and structural modification requirements. • The smallest achievable detection error in displacement was found to be as low as 0.26%, through simulated against the displacement response of a single degree of freedom structure subject to ground motion excitation. • This was verified during experimentation, when a corner-style reflector was placed on a shake table running ground motion data taken from the 4th September 2010 earthquake in Christchurch. These results confirmed the conclusions drawn from simulation.

Research papers, University of Canterbury Library

In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.

Research papers, University of Canterbury Library

This thesis presents the findings from an experimental programme to determine the performance and behaviour of an integrated building incorporating low damage structural and non-structural systems. The systems investigated included post-tensioned rocking concrete frames, articulated floor solutions, low damage claddings and low damage partition systems. As part of a more general aim to increase the resilience of society against earthquake hazards, more emphasis has been given to damage-control design approaches in research. Multiple low-damage earthquake resistant structural and non-structural systems have emerged that are able to withstand high levels of drift or deflections will little or negligible residual. Dry jointed connections, articulated floor solutions, low damage cladding systems and low damage drywall partitions have all been developed separately and successfully tested. In spite of the extensive research effort and the adoption in practice of the low damage systems, work was required to integrate the systems within one building and verify the constructibility, behaviour and performance of the integrated systems. The objectives of this research were to perform dynamic experimental testing of a building which incorporated the low damage systems and acquire data which could be used to dynamically validate numerical models for each of the systems. A three phase experimental programme was devised and performed to dynamically test a half-scale two storey reinforced concrete building on the University of Canterbury shaking table. The three phases of the programme investigated: The structural system only. The rocking connections were tested as Post-Tensioned only connections and Hybrid connections (including dissipators). Two different articulated floor connections were also investigated. Non-structural systems. The Hybrid building was tested with each non-structural system separately; including low damage claddings, low damage partitions and traditional partitions. The Complete building was tested with Hybrid connections, low damage claddings and low damage partitions all integrated within the test specimen. The building was designed based on a full scale prototype building following the direct displacement based design to reach a peak inter-storey drift of 1.6% in a 1/500 year ground motion for a Wellington site. For each test set up, the test specimen was subjected to a ground motion sequence of 39 single direction ground motions. Through the sequence, both the local and global behaviours of the building and integrated systems were recorded in real time. The test specimen was subjected to over 400 ground motions throughout the testing programme. It sustained no significant damage that required reparations other than crumbling of the grout pads. The average peak inter-storey drifts of the buildings were lower than the design value of 1.6%. The low damage non-structural elements were undamaged in the ground motion sequence. The data acquired from each of the phases was used to successfully validate numerical models for each of the low damage systems included in the research.