The 2010 Darfield and 2011 Christchurch Earthquakes triggered extensive liquefaction-induced lateral spreading proximate to streams and rivers in the Christchurch area, causing significant damage to structures and lifelines. A case study in central Christchurch is presented and compares field observations with predicted displacements from the widely adopted empirical model of Youd et al. (2002). Cone penetration testing (CPT), with measured soil gradation indices (fines content and median grain size) on typical fluvial deposits along the Avon River were used to determine the required geotechnical parameters for the model input. The method presented attempts to enable the adoption of the extensive post-quake CPT test records in place of the lower quality and less available Standard Penetration Test (SPT) data required by the original Youd model. The results indicate some agreement between the Youd model predictions and the field observations, while the majority of computed displacements error on the side of over-prediction by more than a factor of two. A sensitivity analysis was performed with respect to the uncertainties used as model input, illustrating the model’s high sensitivity to the input parameters, with median grain size and fines content among the most influential, and suggesting that the use of CPT data to quantify these parameters may lead to variable results.
The seismic performance of soil profiles with potentially liquefiable deposits is a complex phenomenon that requires a thorough understanding of the soil properties and ground motion characteristics. The limitations of simplified liquefaction assessment methods have prompted an increase in the use of non-linear dynamic analysis methods. Focusing on onedimensional site response of a soil column, this thesis validated a soil constitutive model using in-situ pore pressure measurements and then assessed the influence of input ground motion characteristics on soil column response using traditional and newly developed metrics. Pore pressure recordings during the Canterbury Earthquake Sequence (CES) in New Zealand were used to validate the PM4Sand constitutive model. Soil profile characterization was key to accurate prediction of excess pore pressure response and accounting for any densification during the CES. Response during multiple earthquakes was captured effectively and cross-layer interaction demonstrated the model capability to capture soil response at the system-level. Synthetic and observed ground motions from the Christchurch earthquake were applied to the validated soil column to quantify the performance of synthetic motions. New metrics were developed to facilitate a robust comparison to assess performance. The synthetic input motions demonstrated a slightly larger acceleration and excess pore pressure response compared to the observed input motions. The results suggest that the synthetic motions may accumulate higher excess pore pressure at a faster rate and with fewer number of cycles in the shear response. This research compares validated soil profile subject to spectrally-matched pulse and non-pulse motions, emphasizing the inclusion of pulse motions with distinctive characteristics in ground motion suites for non-linear dynamic analysis. However, spectral matching may lead to undesired alterations in pulse characteristics. Cumulative absolute velocity and significant duration significantly differed between these two groups compared to the other key characteristics and contributed considerably to the liquefaction response. Unlike the non-pulse motions, not all of the pulse motions triggered liquefaction, likely due to their shorter significant duration. Non-pulse motions developed a greater spatial extent of liquefaction triggering in the soil profile and extended to a greater depth.
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
In recent years, significant research has been undertaken into the development of lead-extrusion damping technology. The high force-to-volume (HF2V) devices developed at the University of Canterbury have been the subject of much of this research. However, while these devices have undergone a limited range of velocity testing, limitations in test equipment has meant that they have never been tested at representative earthquake velocities. Such testing is important as the peak resistive force provided by the dampers under large velocity spikes is an important design input that must be known for structural applications. This manuscript presents the high-speed testing of HF2V devices with quasi-static force capacities of 250-300kN. These devices have been subjected to peak input velocities of approximately 200mm/s, producing peak resistive forces of approximately 350kN. The devices show stable hysteretic performance, with slight force reduction during high-speed testing due to heat build-up and softening of the lead working material. This force reduction is recovered following cyclic loading as heat is dissipated and the lead hardens again. The devices are shown to be only weakly velocity dependent, an advantage in that they do not deliver large forces to the connecting elements and surrounding structure if larger than expected response velocities occur. This high-speed testing is an important step towards uptake as it provides important information to designers.
In recent years, significant research has been undertaken into the development of lead-extrusion damping technology. The high force-to-volume (HF2V) devices developed at the University of Canterbury have been the subject of much of this research. However, while these devices have undergone a limited range of velocity testing, limitations in test equipment has meant that they have never been tested at representative earthquake velocities. Such testing is important as the peak resistive force provided by the dampers under large velocity spikes is an important design input that must be known for structural applications. This manuscript presents the high-speed testing of HF2V devices with quasi-static force capacities of 250-300kN. These devices have been subjected to peak input velocities of approximately 200mm/s, producing peak resistive forces of approximately 350kN. The devices show stable hysteretic performance, with slight force reduction during high-speed testing due to heat build-up and softening of the lead working material. This force reduction is recovered following cyclic loading as heat is dissipated and the lead hardens again. The devices are shown to be only weakly velocity dependent, an advantage in that they do not deliver large forces to the connecting elements and surrounding structure if larger than expected response velocities occur. This high-speed testing is an important step towards uptake as it provides important information to designers.
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.
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.
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.
Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.
This paper examines the consistency of seismicity and ground motion models, used for seismic hazard analysis in New Zealand, with the observations in the Canterbury earthquakes. An overview is first given of seismicity and ground motion modelling as inputs of probabilistic seismic hazard analysis, whose results form the basis for elastic response spectra in NZS1170.5:2004. The magnitude of earthquakes in the Canterbury earthquake sequence are adequately allowed for in the current NZ seismicity model, however the consideration of ‘background’ earthquakes as point sources at a minimum depth of 10km results in up to a 60% underestimation of the ground motions that such events produce. The ground motion model used in conventional NZ seismic hazard analysis is shown to provide biased predictions of response spectra (over-prediction near T=0.2s , and under-predictions at moderate-to-large vibration periods). Improved ground motion prediction can be achieved using more recent NZ-specific models.
This paper presents the ongoing development of a new 3D seismic velocity model of Canterbury, New Zealand. The model explicitly represents the Canterbury sedimentary basin, and other significant geologic horizons, which are expected to have important implications on observed ground motions. The model utilizes numerous sources of data, including 3D regional tomography with a variable-depth inferred Moho, seismic reflection survey lines, geotechnical boreholes and well logs, spectral analysis of surface waves, and CPT logs which provide velocity constraints over their respective ranges of application. The model provides P- and S-wave velocity and density (i.e. Vp, Vs and p) over a grid of input points, and is presently being utilized in broadband ground motion simulations of the 2010-2011 Canterbury earthquakes. Comparison of simulated ground motions with those observed in the 2010-2011 Canterbury earthquakes will help provide a better understanding of the salient physical processes which characterized the unique set of strong ground motions recorded in this sequence of earthquake events.
This study investigates the uncertainty of simulated earthquake ground motions for smallmagnitude events (Mw 3.5 – 5) in Canterbury, New Zealand. 148 events were simulated with specified uncertainties in: event magnitude, hypocentre location, focal mechanism, high frequency rupture velocity, Brune stress parameter, the site 30-m time-averaged shear wave velocity (Vs30), anelastic attenuation (Q) and high frequency path duration. In order to capture these uncertainties, 25 realisations for each event were generated using the Graves and Pitarka (2015) hybrid broadband simulation approach. Monte-Carlo realisations were drawn from distributions for each uncertainty, to generate a suite of simulation realisations for each event and site. The fit of the multiple simulation realisations to observations were assessed using linear mixed effects regression to generate the systematic source, path and site effects components across all ground motion intensity measure residuals. Findings show that additional uncertainties are required in each of the three source, path, and site components, however the level of output uncertainty is promising considering the input uncertainties included.
There is a growing awareness of the need for the earthquake engineering practice to incorporate in addition to empirical approaches in evaluation of liquefaction hazards advanced methods which can more realistically represent soil behaviour during earthquakes. Currently, this implementation is hindered by a number of challenges mainly associated with the amount of data and user-experience required for such advanced methods. In this study, we present key steps of an advanced seismic effective-stress analysis procedure, which on the one hand can be fully automated and, on the other hand, requires no additional input (at least for preliminary applications) compared to simplified cone penetration test (CPT)-based liquefaction procedures. In this way, effective-stress analysis can be routinely applied for quick, yet more robust estimations of liquefaction hazards, in a similar fashion to the simplified procedures. Important insights regarding the dynamic interactions in liquefying soils and the actual system response of a deposit can be gained from such analyses, as illustrated with the application to two sites from Christchurch, New Zealand.
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.
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.
The 4th of September 2010 Mw 7.1 Darfield (Canterbury) earthquake had generated significant ground shaking within the Christchurch Central Business District (CBD). Despite the apparently significant shaking, the observed structural damage for pre-1970s reinforced concrete (RC) buildings was indeed limited and lower than what was expected for such typology of buildings. This paper explores analytically and qualitatively the different aspects of the "apparent‟ good seismic performance of the pre-1970s RC buildings in the Christchurch CBD, following the earthquake reconnaissance survey by the authors. Damage and building parameters survey result, based on a previously established inventory of building stock of these non-ductile RC buildings, is briefly reported. From an inventory of 75 buildings, one building was selected as a numerical case-study to correlate the observed damage with the non-linear analyses. The result shows that the pre-1970s RC frame buildings performed as expected given the intensity of the ground motion shaking during the Canterbury earthquake. Given the brittle nature of this type of structure, it was demonstrated that more significant structural damage and higher probability of collapse could occur when the buildings were subjected to alternative input signals with different frequency content and duration characteristics and still compatible to the seismicity hazard for Christchurch CBD.
he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.
The need for a simple but rigorous seismic assessment procedure to predict damage to reinforced concrete buildings during a seismic event has been highlighted following the Canterbury Earthquake sequence. Such simplified assessment procedure, applied to individual structure or large building inventory, should not only have low requirement in terms of input information and involve straightforward analyses, but also should be capable to provide reliable predictive results within short timeframe. This research provides a general overview and critical comparison of alternative simplified assessment procedures adopted in NZSEE 2006 Guidelines (Assessment and Improvement of the Structural Performance of Buildings in Earthquakes), ASCE 41-13 (Seismic Evaluation and Retrofit of Existing Buildings), and EN: 1998-3: 2005 (Assessment and Retrofitting of Buildings). Particular focus is given to the evaluation of the capability of Simplified Lateral Mechanism Analysis (SLaMa), which is an analytical pushover method adopted in NZSEE 2006 Guidelines. The predictive results from SLaMa are compared to damages observed for a set of reinforced concrete buildings in Christchurch, as well as the results from more detailed assessment procedure based on numerical modelling. This research also suggests improvements to SLaMa, together with validation of the improvements, to include assessment of local mechanism by strength hierarchy evaluation, as well as to develop assessment of global mechanism including post-yield mechanism sequence based on local mechanism.
This research employs a deterministic seismic risk assessment methodology to assess the potential damage and loss at meshblock level in the Christchurch CBD and Mount Pleasant primarily due to building damage caused by earthquake ground shaking. Expected losses in terms of dollar value and casualties are calculated for two earthquake scenarios. Findings are based on: (1) data describing the earthquake ground shaking and microzonation effects; (2) an inventory of buildings by value, floor area, replacement value, occupancy and age; (3) damage ratios defining the performance of buildings as a function of earthquake intensity; (4) daytime and night-time population distribution data and (5) casualty functions defining casualty risk as a function of building damage. A GIS serves as a platform for collecting, storing and analyzing the original and the derived data. It also allows for easy display of input and output data, providing a critical functionality for communication of outcomes. The results of this study suggest that economic losses due to building damage in the Christchurch CBD and Mount Pleasant will possibly be in the order of $5.6 and $35.3 million in a magnitude 8.0 Alpine fault earthquake and a magnitude 7.0 Ashley fault earthquake respectively. Damage to non-residential buildings constitutes the vast majority of the economic loss. Casualty numbers are expected to be between 0 and 10.
A team of earthquake geologists, seismologists and engineering seismologists from GNS Science, NIWA, University of Canterbury, and Victoria University of Wellington have collectively produced an update of the 2002 national probabilistic seismic hazard (PSH) model for New Zealand. The new model incorporates over 200 new onshore and offshore fault sources, and utilises newly developed New Zealand-based scaling relationships and methods for the parameterisation of the fault and subduction interface sources. The background seismicity model has also been updated to include new seismicity data, a new seismicity regionalisation, and improved methodology for calculation of the seismicity parameters. Background seismicity models allow for the occurrence of earthquakes away from the known fault sources, and are typically modelled as a grid of earthquake sources with rate parameters assigned from the historical seismicity catalogue. The Greendale Fault, which ruptured during the M7.1, 4 September 2010 Darfield earthquake, was unknown prior to the earthquake. However, the earthquake was to some extent accounted for in the PSH model. The maximum magnitude assumed in the background seismicity model for the area of the earthquake is 7.2 (larger than the Darfield event), but the location and geometry of the fault are not represented. Deaggregations of the PSH model for Christchurch at return periods of 500 years and above show that M7-7.5 fault and background source-derived earthquakes at distances less than 40 km are important contributors to the hazard. Therefore, earthquakes similar to the Darfield event feature prominently in the PSH model, even though the Greendale Fault was not an explicit model input.
Despite their good performance in terms of their design objectives, many modern code-prescriptive buildings built in Christchurch, New Zealand had to be razed after the 2010-2011 Canterbury earthquakes because repairs were deemed too costly due to widespread sacrificial damage. Clearly a more effective design paradigm is needed to create more resilient structures. Rocking, post-tensioned connections with supplemental energy dissipation can contribute to a damage avoidance designs (DAD). However, few have achieved all three key design objectives of damage-resistant rocking, inherent recentering ability, and repeatable, damage-free energy dissipation for all cycles, which together offer a response which is independent of loading history. Results of experimental tests are presented for a near full-scale rocking beam-column sub-assemblage. A matrix of test results is presented for the system under varying levels of posttensioning, with and without supplemental dampers. Importantly, this parametric study delineates each contribution to response. Practical limitations on posttensioning are identified: a minimum to ensure static structural re-centering, and a maximum to ensure deformability without threadbar yielding. Good agreement between a mechanistic model and experimental results over all parameters and inputs indicates the model is robust and accurate for design. The overall results indicate that it is possible to create a DAD connection where the non-linear force-deformation response is loading history independent and repeatable over numerous loading cycles, without damage, creating the opportunity for the design and implementation of highly resilient structures.
This study analyses the success and limitations of the recovery process following the 2010–11 earthquake sequence in Christchurch, New Zealand. Data were obtained from in-depth interviews with 32 relocated households in Christchurch, and from a review of recovery policies implemented by the government. A top-down approach to disaster recovery was evident, with the creation of multiple government agencies and processes that made grassroots input into decision-making difficult. Although insurance proceeds enabled the repair and rebuilding of many dwellings, the complexity and adversarial nature of the claim procedures also impaired recovery. Householders’ perceptions of recovery reflected key aspects of their post-earthquake experiences (e.g. the housing offer they received, and the negotiations involved), and the outcomes of their relocation (including the value of the new home, their subjective well-being, and lifestyle after relocation). Protracted insurance negotiations, unfair offers and hardships in post-earthquake life were major challenges to recovery. Less-thanfavourable recovery experiences also transformed patterns of trust in local communities, as relocated householders came to doubt both the government and private insurance companies’ ability to successfully manage a disaster. At the same time, many relocated households expressed trust in their neighbours and communities. This study illuminates how government policies influence disaster recovery while also suggesting a need to reconsider centralised, top-down approaches to managing recovery.
Rapid, accurate structural health monitoring (SHM) assesses damage to optimise decision-making. Many SHM methods are designed to track nonlinear stiffness changes as damage. However, highly nonlinear pinched hysteretic systems are problematic in SHM. Model-based SHM often fails as any mismatch between model and measured response dynamics leads to significant error. Thus, modelfree methods of hysteresis loop tracking methods have emerged. This study compares the robustness and accuracy in the presence of significant measurement noise of the proven hysteresis loop analysis (HLA) SHM method with 3 emerging model-free methods and 2 further novel adaptations of these methods using a highly nonlinear, 6-story numerical structure to provide a known ground-truth. Mean absolute errors in identifying a known nonlinear stiffness trajectory assessed at four points over two successive ground motion inputs from September 2010 and February 2011 in Christchurch range from 1.71-10.52%. However, the variability is far wider with maximum errors ranging from 3.90-49.72%, where the second largest maximum absolute error was still 19.74%. The lowest mean and maximum absolute errors were for the HLA method. The next best method had mean absolute error of 2.92% and a maximum of 10.51%. These results show the clear superiority of the HLA method over all current emerging model-free methods designed to manage the highly nonlinear pinching responses common in reinforced concrete structures. These results, combined with high robustness and accuracy in scaled and fullscale experimental studies, provide further validation for using HLA for practical implementation.
The context of this study is the increasing need for public transport as issues over high private vehicle usage are becoming increasingly obvious. Public transport services need to compete with private transport to improve patronage, and issues with reliability need to be addressed. Bus bunching affects reliability through disruptions to the scheduled headways. The purpose of this study was to collect and analyse data to compare how travel time and dwell time vary, to explore the variation of key variables, and to better understand the sources of these variations. The Orbiter bus service in Christchurch was used as a case study, as it is particularly vulnerable to bus bunching. The dwell time was found to be more variable than travel time. It appeared the Canterbury earthquake had significantly reduced the average speeds for the Orbiter service. In 1964, Newell and Potts described a basic bus bunching theory, which was used as the basis for an Excel bus bunching model. This model allows input variables to vary stochastically. Random values were generated from four specified distributions derived from manually collected data, allowing variance across all bus platforms and buses. However the complexity resulted in stability and difficulty in achieving convergence, so the model was run in single Monte Carlo simulations. The outputs were realistic and showed a higher degree of bunching behaviour than previous models. The model demonstrated bunching phenomena that had not been observed in previous models, including spontaneously un-pairing, overtaking of buses delayed at platforms, and odd-numbered bunches of three buses. Furthermore, the study identified areas of further research for data collection and model development.
The Canterbury earthquake and aftershock sequence in New Zealand during 2010-2011 subjected the city’s structures to a significant accumulated cyclic demand and raised significant questions regarding the low-cycle fatigue demands imposed upon the structures. There is a significant challenge to quantify the level of cumulative demand imposed on structures and to assess the percentage of a structure's fatigue life that has been consumed as a result of this earthquake sequence. It is important to be able to quantify the cumulative demand to determine how a building will perform in a subsequent large earthquake and inform repair and re-occupancy decisions. This paper investigates the cumulative fatigue demand for a structure located within the Christchurch Central Business District (CBD). Time history analysis and equivalent cycle counting methods are applied across the Canterbury earthquake sequence, using key events from September 4th 2010 and February 22nd , 2011 main shocks. The estimate of the cumulative fatigue demand is then compared to the expected capacity of a case study reinforced concrete bridge pier, to undertake a structure-specific fatigue assessment. The analysis is undertaken to approximate the portion of the structural fatigue capacity that has been consumed, and how much residual capacity remains. Results are assessed for recordings at the four Christchurch central city strong motion recording sites installed by the GeoNet programme, to provide an estimate of variation in results. The computed cyclic demand results are compared to code-based design methods and as assessment of the inelastic displacement demand of the reinforcing steel. Results are also presented in a fragility context where a de minimis (inconsequential), irreparable damage and full fatigue fracture are defined to provide a probabilistic assessment of the fatigue damage incurred. This methodology can provide input into the overall assessment of fatigue demands and residual capacity.
The Mw 6.2 February 22nd 2011 Christchurch earthquake (and others in the 2010-2011 Canterbury sequence) provided a unique opportunity to study the devastating effects of earthquakes first-hand and learn from them for future engineering applications. All major events in the Canterbury earthquake sequence caused widespread liquefaction throughout Christchurch’s eastern suburbs, particularly extensive and severe during the February 22nd event. Along large stretches of the Avon River banks (and to a lesser extent along the Heathcote) significant lateral spreading occurred, affecting bridges and the infrastructure they support. The first stage of this research involved conducting detailed field reconnaissance to document liquefaction and lateral spreading-induced damage to several case study bridges along the Avon River. The case study bridges cover a range of ages and construction types but all are reinforced concrete structures which have relatively short, stiff decks. These factors combined led to a characteristic deformation mechanism involving deck-pinning and abutment back-rotation with consequent damage to the abutment piles and slumping of the approaches. The second stage of the research involved using pseudo-static analysis, a simplified seismic modelling tool, to analyse two of the bridges. An advantage of pseudo-static analysis over more complicated modelling methods is that it uses conventional geotechnical data in its inputs, such as SPT blowcount and CPT cone resistance and local friction. Pseudo-static analysis can also be applied without excessive computational power or specialised knowledge, yet it has been shown to capture the basic mechanisms of pile behaviour. Single pile and whole bridge models were constructed for each bridge, and both cyclic and lateral spreading phases of loading were investigated. Parametric studies were carried out which varied the values of key parameters to identify their influence on pile response, and computed displacements and damages were compared with observations made in the field. It was shown that pseudo-static analysis was able to capture the characteristic damage mechanisms observed in the field, however the treatment of key parameters affecting pile response is of primary importance. Recommendations were made concerning the treatment of these governing parameters controlling pile response. In this way the future application of pseudo-static analysis as a tool for analysing and designing bridge pile foundations in liquefying and laterally spreading soils is enhanced.
There is a relationship between inelastic deformation and energy dissipation in structures that are subjected to earthquake ground motions. Thus, if seismic energy dissipation can be achieved by means of a separate non-load bearing supplementary damping system, the load bearing structure can remain elastic with continuing serviceability following the design level earthquake. This research was carried out to investigate the advantages of using added damping in structures. The control system consists of passive friction dampers called ring spring dampers installed in the ground floor of the structure using a tendon to transmit the forces to the other parts of the structure. The ring springs dampers are friction devices consisting of inner and outer ring elements assembled to form a spring stack. External load applied to the spring produces sliding action across mating ring interfaces. The damping forces generated by the dampers and transferred in the supplemental system to the structure by the tendon and horizontal links oppose the internal loads. A four storey-two bay steel frame structure was used in the study. Experimental and analytical studies to investigate the effectiveness of a supplemental control system are presented. The model was subjected to a series of earthquake simulations on the shaking table in the Structural Laboratory of the Civil Engineering Department, at the University of Canterbury. The earthquake simulation tests have been performed on the structure both with and without the supplemental control system. The earthquake simulations were a series of gradually increasing intensity replications of two commonly used earthquake records. This thesis includes detailed description of the structural model, the supplemental control system, the ring springs dampers and the data obtained during the testing. Analyses were then carried out on a twelve storey framed structure to investigate the possible tendon arrangements and the size and type of dampers required to control the response of a real building. Guidelines for determining the appropriateness of including a supplemental damping system have been investigated. The main features of the supplemental control system adopted in this research are: • It is a passive control system with extreme reliability and having no dependence on external power sources to effect the control action. These power sources may not be available during a major earthquake. • Ring springs are steel friction devices capable of absorbing large amounts of input energy. No liquid leakage can occur and minimal maintenance is required for the ring spring dampers. • With a damper-tendon system, the distribution of the dampers throughout the structure is not so critical. Only one or two dampers are used to produce the damping forces needed, and forces are then transferred to the rest of the building by the tendon system. • It is a relatively inexpensive control system with a long useful life.
The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
Seismically vulnerable buildings constitute a major problem for the safety of human beings. In many parts of the world, reinforced concrete (RC) frame buildings designed and constructed with substandard detailing, no consideration of capacity design principles, and improper or no inclusion of the seismic actions, have been identified. Amongst those vulnerable building, one particular typology representative of the construction practice of the years previous to the 1970’s, that most likely represents the worst case scenario, has been widely investigated in the past. The deficiencies of that building typology are related to non-ductile detailing in beam column joints such as the use of plain round bars, the lack of stirrups inside the joint around the longitudinal reinforcement of the column, the use of 180° end hooks in the beams, the use of lap splices in potential ‘plastic hinge’ regions, and substandard quality of the materials. That type of detailing and the lack of a capacity design philosophy create a very fragile fuse in the structure where brittle inelastic behaviour is expected to occur, which is the panel zone region of exterior beam column joints. The non-ductile typology described above was extensively investigated at the University of Canterbury in the context of the project ‘Retrofit Solutions for New Zealand Multi-Storey Buildings’ (2004-2011), founded by the ‘Foundation for Research, Science and Technology’ Tūāpapa Rangahau Pūtaiao. The experimental campaign prior to the research carried out by the author consisted of quasi-static tests of beam column joint subassemblies subjected to lateral loading regime, with constant and varying axial load in the column. Most of those specimens were representative of a plane 2D frame (knee joint), while others represented a portion of a space 3D frame (corner joints), and only few of them had a floor slab, transverse beams, and lap splices. Using those experiments, several feasible, cost-effective, and non-invasive retrofit techniques were developed, improved, and refined. Nevertheless, the slow motion nature of those experiments did not take into account the dynamical component inherent to earthquake related problems. Amongst the set of techniques investigated, the use of FRP layers for strengthening beam column joints is of particular interest due to its versatility and the momentum that its use has gained in the current state of the practice. That particular retrofit technique was previously used to develop a strengthening scheme suitable for plane 2D and space 3D corner beam column joints, but lacking of floor slabs. In addition, a similar scheme was not developed for exterior joints of internal frames, referred here as ‘cruciform’. In this research a 2/5 scale RC frame model building comprising of two frames in parallel (external and internal) joined together by means of floor slabs and transverse beams, with non-ductile characteristics identical to those of the specimens investigated previously by others, and also including lap splices, was developed. In order to investigate the dynamic response of that building, a series of shake table tests with different ground motions were performed. After the first series of tests, the specimen was modified by connecting the spliced reinforcement in the columns in order to capture a different failure mode. Ground motions recorded during seismic events that occurred during the initial period of the experimental campaign (2010) were used in the subsequent experiments. The hierarchy of strengths and sequence of events in the panel zone region were evaluated in an extended version of the bending moment-axial load (M-N) performance domain developed by others. That extension was required due to the asymmetry in the beam cross section introduced by the floor slab. In addition, the effect of the torsion resistance provided by the spandrel (transverse beam) was included. In order to upgrade the brittle and unstable performance of the as-built/repaired specimen, a practical and suitable ad-hoc FRP retrofit intervention was developed, following a partial retrofit strategy that aimed to strengthen exterior beam column joints only (corner and cruciform). The ability of the new FRP scheme to revert the sequence of events in the panel zone region was evaluated using the extended version of the M-N performance domain as well as the guidelines for strengthening plane joints developed by others. Weakening of the floor slab in a novel configuration was also incorporated with the purpose of reducing the flexural capacity of the beam under negative bending moment (slab in tension), enabling the damage relocation from the joint into the beam. The efficacy of the developed retrofit intervention in upgrading the seismic performance of the as-built specimen was investigated using shake table tests with the input motions used in the experiments of the as-built/repaired specimen. Numerical work aimed to predict the response of the model building during the most relevant shake table tests was carried out. By using a simple numerical model with concentrated plasticity elements constructed in Ruaumoko2D, the results of blind and post-experimental predictions of the response of the specimen were addressed. Differences in the predicted response of the building using the nominal and the actual recorded motions of the shake table were investigated. The dependence of the accuracy of the numerical predictions on the assumed values of the parameters that control the hysteresis rules of key structural members was reviewed. During the execution of the experimental campaign part of this thesis, two major earthquakes affected the central part of Chile (27 of February 2010 Maule earthquake) and the Canterbury region in New Zealand (22 February 2011 Canterbury earthquake), respectively. As the author had the opportunity to experience those events and investigate their consequences in structures, the observations related to non-ductile detailing and drawbacks in the state of the practice related to reinforced concrete walls was also addressed in this research, resulting in preliminary recommendations for the refinement of current seismic code provisions and assessment guidelines. The investigations of the ground motions recorded during those and other earthquakes were used to review the procedures related to the input motions used for nonlinear dynamic analysis of buildings as required by most of the current code provisions. Inelastic displacement spectra were constructed using ground motions recorded during the earthquakes mentioned above, in order to investigate the adequacy of modification factors used to obtain reduced design spectra from elastic counterparts. Finally a simplified assessment procedure for RC walls that incorporates capacity compatible spectral demands is proposed.
During 2010 and 2011, major earthquakes caused widespread damage and the deaths of 185 people in the city of Christchurch. Damaged school buildings resulted in state intervention which required amendment of the Education Act of 1989, and the development of ‘site sharing agreements’ in undamaged schools to cater for the needs of students whose schools had closed. An effective plan was also developed for student assessment through establishing an earthquake impaired derived grade process. Previous research into traditional explanations of educational inequalities in the United Kingdom, the United States of America, and New Zealand were reviewed through various processes within three educational inputs: the student, the school and the state. Research into the impacts of urban natural disasters on education and education inequalities found literature on post disaster education systems but nothing could be found that included performance data. The impacts of the Canterbury earthquakes on educational inequalities and achievement were analysed over 2009-2012. The baseline year was 2009, the year before the first earthquake, while 2012 is seen as the recovery year as no schools closed due to seismic events and there was no state intervention into the education of the region. National Certificate of Educational Achievement (NCEA) results levels 1-3 from thirty-four secondary schools in the greater Christchurch region were graphed and analysed. Regression analysis indicates; in 2009, educational inequalities existed with a strong positive relationship between a school’s decile rating and NCEA achievement. When schools were grouped into decile rankings (1-10) and their 2010 NCEA levels 1-3 results were compared with the previous year, the percentage of change indicates an overall lower NCEA achievement in 2010 across all deciles, but particularly in lower decile schools. By contrast, when 2011 NCEA results were compared with those of 2009, as a percentage of change, lower decile schools fared better. Non site sharing schools also achieved higher results than site sharing schools. State interventions, had however contributed towards student’s achieving national examinations and entry to university in 2011. When NCEA results for 2012 were compared to 2009 educational inequalities still exist, however in 2012 the positive relationship between decile rating and achievement is marginally weaker than in 2009. Human ethics approval was required to survey one Christchurch secondary school community of students (aged between 12 and 18), teachers and staff, parents and caregivers during October 2011. Participation was voluntary and without incentives, 154 completed questionnaires were received. The Canterbury earthquakes and aftershocks changed the lives of the research participants. This school community was displaced to another school due to the Christchurch earthquake on 22 February 2011. Research results are grouped under four geographical perspectives; spatial impacts, socio-economic impacts, displacement, and health and wellbeing. Further research possibilities include researching the lag effects from the Canterbury earthquakes on school age children.
