This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
Deconstruction, at the end of the useful life of a building, produces a considerable amount of materials which must be disposed of, or be recycled / reused. At present, in New Zealand, most timber construction and demolition (C&D) material, particularly treated timber, is simply waste and is placed in landfills. For both technical and economic reasons (and despite the increasing cost of landfills), this position is unlikely to change in the next 10 – 15 years unless legislation dictates otherwise. Careful deconstruction, as opposed to demolition, can provide some timber materials which can be immediately re-used (eg. doors and windows), or further processed into other components (eg. beams or walls) or recycled (‘cascaded’) into other timber or composite products (e.g. fibre-board). This reusing / recycling of materials is being driven slowly in NZ by legislation, the ‘greening’ of the construction industry and public pressure. However, the recovery of useful material can be expensive and uneconomic (as opposed to land-filling). In NZ, there are few facilities which are able to sort and separate timber materials from other waste, although the soon-to-be commissioned Burwood Resource Recovery Park in Christchurch will attempt to deal with significant quantities of demolition waste from the recent earthquakes. The success (or otherwise) of this operation should provide good information as to how future C&D waste will be managed in NZ. In NZ, there are only a few, small scale facilities which are able to burn waste wood for energy recovery (e.g. timber mills), and none are known to be able to handle large quantities of treated timber. Such facilities, with constantly improving technology, are being commissioned in Europe (often with Government subsidies) and this indicates that similar bio-energy (co)generation will be established in NZ in the future. However, at present, the NZ Government provides little assistance to the bio-energy industry and the emergence worldwide of shale-gas reserves is likely to push the economic viability of bio-energy further into the future. The behaviour of timber materials placed in landfills is complex and poorly understood. Degrading timber in landfills has the potential to generate methane, a potent greenhouse gas, which can escape to the atmosphere and cancel out the significant benefits of carbon sequestration during tree growth. Improving security of landfills and more effective and efficient collection and utilisation of methane from landfills in NZ will significantly reduce the potential for leakage of methane to the atmosphere, acting as an offset to the continuing use of underground fossil fuels. Life cycle assessment (LCA), an increasingly important methodology for quantifying the environmental impacts of building materials (particularly energy, and global warming potential (GWP)), will soon be incorporated into the NZ Green Building Council Greenstar rating tools. Such LCA studies must provide a level playing field for all building materials and consider the whole life cycle. Whilst the end-of-life treatment of timber by LCA may establish a present-day base scenario, any analysis must also present a realistic end-of-life scenario for the future deconstruction of any 6 new building, as any building built today will be deconstructed many years in the future, when very different technologies will be available to deal with construction waste. At present, LCA practitioners in NZ and Australia place much value on a single research document on the degradation of timber in landfills (Ximenes et al., 2008). This leads to an end-of-life base scenario for timber which many in the industry consider to be an overestimation of the potential negative effects of methane generation. In Europe, the base scenario for wood disposal is cascading timber products and then burning for energy recovery, which normally significantly reduces any negative effects of the end-of-life for timber. LCA studies in NZ should always provide a sensitivity analysis for the end-of-life of timber and strongly and confidently argue that alternative future scenarios are realistic disposal options for buildings deconstructed in the future. Data-sets for environmental impacts (such as GWP) of building materials in NZ are limited and based on few research studies. The compilation of comprehensive data-sets with country-specific information for all building materials is considered a priority, preferably accounting for end-of-life options. The NZ timber industry should continue to ‘champion’ the environmental credentials of timber, over and above those of the other major building materials (concrete and steel). End-of-life should not be considered the ‘Achilles heel’ of the timber story.
A seismic financial risk analysis of typical New Zealand reinforced concrete buildings constructed with topped precast concrete hollow-core units is performed on the basis of experimental research undertaken at the University of Canterbury over the last five years. An extensive study that examines seismic demands on a variety of multi-storey RC buildings is described and supplemented by the experimental results to determine the inter-storey drift capacities of the buildings. Results of a full-scale precast concrete super-assemblage constructed and tested in the laboratory in two stages are used. The first stage investigates existing construction and demonstrates major shortcomings in construction practice that would lead to very poor seismic performance. The second stage examines the performance of the details provided by Amendment No. 3 to the New Zealand Concrete Design Code NZS 3101:1995. This paper uses a probabilistic financial risk assessment framework to estimate the expected annual loss (EAL) from previously developed fragility curves of RC buildings with precast hollow core floors connected to the frames according to the pre-2004 standard and the two connection details recommended in the 2004 amendment. Risks posed by different levels of damage and by earthquakes of different frequencies are examined. The structural performance and financial implications of the three different connection details are compared. The study shows that the improved connection details recommended in the 2004 amendment give a significant economic payback in terms of drastically reduced financial risk, which is also representative of smaller maintenance cost and cheaper insurance premiums.
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.
On February 22, 2011, Christchurch-based journalists were jolted out of their normal work routine by a large 6.3 magnitude earthquake that killed 185 people, wrecked the city and forced reporters to reappraise their journalism. This study considers how the earthquake affected journalists’ relationship to the community, their use of sources and news selection. A theory of collective trauma is used to explain the changes that journalists made to their reporting practice. Specifically, Christchurch journalists had a greater identification and attachment to their audience post-earthquake. Journalists viewed themselves as part of the earthquake story, which prompted them to view sources differently, use those sources differently and see advocacy as a keystone of their news work after the disaster. This study adds to a growing scholarship about journalists and trauma, but focuses on what the event meant for local reporters’ choice of sources and news selection rather than measuring rates of psychological distress.
Paper 3408The paper discusses modelling of cyclic stress-strain behaviour of soil, in particular a simple model that can produce a desired stiffness and hysteretic damping for a given strain level as observed in laboratory testing is formulated. The unloading-reloading relationship is developed for total stress seismic site response analysis with appropriate damping at large strain. The constitutive model employs a hyperbolic equation as the backbone curve, and uses a modification of the extended Masing unloading-reloading relationship leading to correct measured modulus reduction and damping curves simultaneously. A quasi-static cyclic loading of increasing amplitude is used to demonstrate the model’s performance and its capability to allow improved modelling of the magnitude of energy dissipation based on an experimental program on native sandy soils from Christchurch, New Zealand.
This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
Triple P parenting programmes have provided promising results for children and families in recent years. The aim of the current project was to explore the experiences of families leading up to participating in a Teen Triple P programme three years following the Christchurch earthquakes and their need for assistance in the management of their teenagers. Parents were interviewed prior to the commencement of the Teen Triple P programme and after its completion. Parents were also asked to complete a journal entry or engage in two brief telephone conversations with the researcher outlining their experiences with the Teen Triple P programme. These outlined the perceived fit of the programme to the needs of the family. Parents provided insight into their family’s experiences of the Christchurch 2010 and 2011 series of earthquakes and the perceived impact this had on their lives and the management of their teenagers. The results indicated that parents felt more positively about their parenting behaviours post-programme and were able to identify changes in their teen and/or family that they felt were as a response to participation in Teen Triple P. Parents provided rich descriptions of their earthquake experiences and the immediate and long-term impacts they endured both individually and as a family. Parents did not feel that the earthquakes fed into their decision to do a Teen Triple P Programme. The results helped improve our understanding of the effectiveness of Teen Triple P as a parenting programme as well increased our understanding of the challenges and needs of families in post-earthquake Christchurch.
The Canterbury region of New Zealand experienced a sequence of strong earthquakes during 2010-2011. Responses included government acquisition of many thousands of residential properties in the city of Christchurch in areas with severe earthquake effects. A large and contiguous tract of this ‘red zoned’ land lies in close proximity to the Ōtākaro / Avon River and is known as the Avon-Ōtākaro Red Zone (AORZ). The focus of this study was to provide an overview of the floodplain characteristics of the AORZ and review of international experience in ecological restoration of similar river margin and floodplain ecosystems to extract restoration principles and associated learnings. Compared to pre-earthquake ground levels, the dominant trend in the AORZ is subsidence, together with lateral movement especially in the vicinity of waterway. An important consequence of land subsidence in the lower Ōtākaro / Avon River is greater exposure to flooding and the effects of sea level rise. Scenario modelling for sea level rise indicates that much of the AORZ is exposed to inundation within a 100 year planning horizon based on a 1 m sea level rise. As with decisions on built infrastructure, investments in nature-based ‘green infrastructure’ also require a sound business case including attention to risks posed by climate change. Future-proofing of the expected benefits of ecological restoration must therefore be secured by design. Understanding and managing the hydrology and floodplain dynamics are vital to the future of the AORZ. However, these characteristics are shared by other floodplain and river restoration projects worldwide. Identifying successful approaches provides a useful a source of useful information for floodplain planning in the AORZ. This report presents results from a comparative case study of three international examples to identify relevant principles for large-scale floodplain management at coastal lowland sites.
This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized
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.
Surface-rupturing earthquakes can trigger the sudden avulsion of river channels, causing rapid and persistent coseismic flooding of previously unaffected areas. This phenomenon, known as fault-rupture-induced river avulsion (FIRA), occurs when fault displacement significantly alters river channel topography. The importance of understanding FIRA as a secondary seismic hazard was highlighted by events during the 2010 Darfield and 2016 Kaikoura earthquakes in New Zealand. This thesis develops a national model to identify and quantify FIRA susceptibility across New Zealand by integrating hydrological datasets (NIWA RiverMaps and Flood Statistics) with active fault information (NZ Active Faults Database and RSQSim earthquake simulations). The methodology applies the F-index framework proposed by McEwan et al. (2023), which quantifies FIRA potential based on the ratio of fault throw plus discharge-dependent depth to bank full depth at each fault-river intersection. The model successfully identified 3,796 potential FIRA-susceptible fault-river intersections nationwide, with 451 involving waterways equal to or larger than the Hororata River. Regional analysis revealed higher concentrations of FIRA-susceptible sites in the Bay of Plenty, Canterbury, and Marlborough regions. Validation against historical events showed the model effectively located known FIRA occurrences from the Kaikoura and Darfield earthquakes, though with some limitations in accurately predicting F-index values due to complex fault displacement patterns and challenges in modelling bank full depths of large, braided rivers. This research establishes New Zealand's first nationwide assessment of fault-induced river avulsion susceptibility. The approach creates a structured methodology for identifying high-risk fault-river intersections and determining which sites require thorough localised examination. The methodology developed offers a template for similar assessments in other tectonically active regions and contributes to improving earthquake hazard assessment and disaster preparedness planning.
The Resilient Shorelines study at University of Canterbury (UC) is using the Avon Heathcote Estuary Ihutai to investigate ecosystem-based approaches to conservation planning and adaptation in response to environmental change. In particular, the study is using a novel opportunity to understand effects of the Canterbury earthquakes that may be similar to impacts of sea level rise. These result from topographic and bathymetry changes in and around the estuary and associated waterways (Beaven et al., 2012; Cochran et al., 2014) that have driven changes in hydrodynamics (Measures et al., 2011). Therefore the wider context for the work reported here is to develop methodologies for modelling the impacts of sea level rise on estuaries and coastal river mouths using the Avon-Heathcote Estuary/Ihutai as a case study. Initial objectives have included establishing the magnitude of earthquake-induced changes. Subsequent steps will include establishing the relationships between strong physical drivers such as water levels and salinity, and the spatial pattern of estuarine ecosystems. There is particular focus on understanding salinity changes in the upper estuarine ecosystem in the vicinity of the freshwater-saltwater interface. In these areas, species, habitats and ecosystems that are adapted to brackish conditions are expected to migrate in response to the inland penetration of salt water under sea level rise. An example is the location of īnanga spawning habitat that is associated with the inland extent of salt water intrusion on spring tides (Taylor, 2002). It is expected to be strongly affected by sea level rise. To facilitate the development of ecosystem-based scenario models for sea level rise, a salinity model with resolution at ecological meaningful scales was required. An existing fine scale hydrodynamic model was available using Delft3D software (Deltares, 2012) that had been developed for ECan and MBIE following the earthquakes (Measures & Bind, 2013). However, it had not been calibrated for salinity. A collaborative project was designed between UC and NIWA to calibrate the model and develop a scenario modelling approach for sea level rise at a level of resolution sufficient for understanding sea level rise impacts on īnanga (whitebait) spawning habitat. The project was allocated funding from Brian Mason Scientific and Technical Trust and commenced in late 2015. The purpose of this report is to provide a description of the model development process and an illustration of model outputs from an initial set of modelled scenarios for sea level rise.
In this thesis, focus is given to develop methodologies for rapidly estimating specific components of loss and downtime functions. The thesis proposes methodologies for deriving loss functions by (i) considering individual component performance; (ii) grouping them as per their performance characteristics; and (iii) applying them to similar building usage categories. The degree of variation in building stock and understanding their characteristics are important factors to be considered in the loss estimation methodology and the field surveys carried out to collect data add value to the study. To facilitate developing ‘downtime’ functions, this study investigates two key components of downtime: (i) time delay from post-event damage assessment of properties; and (ii) time delay in settling the insurance claims lodged. In these two areas, this research enables understanding of critical factors that influence certain aspects of downtime and suggests approaches to quantify those factors. By scrutinising the residential damage insurance claims data provided by the Earthquake Commission (EQC) for the 2010- 2011 Canterbury Earthquake Sequence (CES), this work provides insights into various processes of claims settlement, the time taken to complete them and the EQC loss contributions to building stock in Christchurch city and Canterbury region. The study has shown diligence in investigating the EQC insurance claim data obtained from the CES to get new insights and build confidence in the models developed and the results generated. The first stage of this research develops contribution functions (probabilistic relationships between the expected losses for a wide range of building components and the building’s maximum response) for common types of claddings used in New Zealand buildings combining the probabilistic density functions (developed using the quantity of claddings measured from Christchurch buildings), fragility functions (obtained from the published literature) and cost functions (developed based on inputs from builders) through Monte Carlo simulations. From the developed contribution functions, glazing, masonry veneer, monolithic and precast concrete cladding systems are found to incur 50% loss at inter-storey drift levels equal to 0.027, 0.003, 0.005 and 0.011, respectively. Further, the maximum expected cladding loss for glazing, masonry veneer, monolithic, precast concrete cladding systems are found to be 368.2, 331.9, 365.0, and 136.2 NZD per square meter of floor area, respectively. In the second stage of this research, a detailed cost breakdown of typical buildings designed and built for different purposes is conducted. The contributions of structural and non- structural components to the total building cost are compared for buildings of different usages, and based on the similar ratios of non-structural performance group costs to the structural performance group cost, four-building groups are identified; (i) Structural components dominant group: outdoor sports, stadiums, parkings and long-span warehouses, (ii) non- structural drift-sensitive components dominant group: houses, single-storey suburban buildings (all usages), theatres/halls, workshops and clubhouses, (iii) non-structural acceleration- sensitive components dominant group: hospitals, research labs, museums and retail/cold stores, and (iv) apartments, hotels, offices, industrials, indoor sports, classrooms, devotionals and aquariums. By statistically analysing the cost breakdowns, performance group weighting factors are proposed for structural, and acceleration-sensitive and drift-sensitive non-structural components for all four building groups. Thus proposed building usage groupings and corresponding weighting factors facilitate rapid seismic loss estimation of any type of building given the EDPs at storey levels are known. A model for the quantification of post-earthquake inspection duration is developed in the third stage of this research. Herein, phase durations for the three assessment phases (one rapid impact and two rapid building) are computed using the number of buildings needing inspections, the number of engineers involved in inspections and a phase duration coefficient (which considers the median building inspection time, efficiency of engineer and the number of engineers involved in each assessment teams). The proposed model can be used: (i) by national/regional authorities to decide the length of the emergency period following a major earthquake, and estimate the number of engineers required to conduct a post-earthquake inspection within the desired emergency period, and (ii) to quantify the delay due to inspection for the downtime modelling framework. The final stage of this research investigates the repair costs and insurance claim settlement time for damaged residential buildings in the 2010-2011 Canterbury earthquake sequence. Based on the EQC claim settlement process, claims are categorized into three groups; (i) Small Claims: claims less than NZD15,000 which were settled through cash payment, (ii) Medium Claims: claims less than NZD100,000 which were managed through Canterbury Home Repair Programme (CHRP), and (iii) Large Claims: claims above NZD100,000 which were managed by an insurance provider. The regional loss ratio (RLR) for greater Christchurch for three events inducing shakings of approximate seismic intensities 6, 7, and 8 are found to be 0.013, 0.066, and 0.171, respectively. Furthermore, the claim duration (time between an event and the claim lodgement date), assessment duration (time between the claim lodgement day and the most recent assessment day), and repair duration (time between the most recent assessment day and the repair completion day) for the insured residential buildings in the region affected by the Canterbury earthquake sequence is found to be in the range of 0.5-4 weeks, 1.5- 5 months, and 1-3 years, respectively. The results of this phase will provide useful information to earthquake engineering researchers working on seismic risk/loss and insurance modelling.
The development of cheap, whilst effective and relatively non-invasive structural retrofit techniques for existing non-ductile reinforced concrete (RC) structures still remains the most challenging issue for a wide implementation on a macro scale. Seismic retrofit is too often being confused as purely structural strengthening. As part of a six-years national project on “Seismic retrofit solutions for NZ multi-storey building”, focus has been given at the University of Canterbury on the development of a counter-intuitive retrofit strategy for earthquake vulnerable existing rc frame, based on a “selective weakening” (SW) approach. After an overview of the SW concept, this paper presents the experimental and numerical validation of a SW retrofit strategy for earthquake vulnerable existing RC frame with particular focus on the exterior beam-column (b-c) joints. The exterior b-c joint is a critically vulnerable region in many existing pre-1970s RC frames. By selectively weakening the beam by cutting the bottom longitudinal reinforcements and/or adding external pre-stressing to the b-c joint, a more desirable inelastic mechanism can be attained, leading to improved global seismic performance. The so-called SW retrofit is implemented on four 2/3-scaled exterior RC b-c joint subassemblies, tested under quasi-static cyclic loading at the University of Canterbury. Complemented by refined 3D Finite Element (FE) models and dynamic time-history analyses results, the experimental results have shown the potential of a simple and cost-effective yet structurally efficient structural rehabilitation technique. The research also demonstrated the potential of advanced 3D fracture-mechanics-based microplane concrete modelling for refined FE analysis of non-ductile RC b-c joints.
Liquefaction during the 4th September 2010 Mw 7.1 Darfield earthquake and large aftershocks in 2011 (Canterbury earthquake sequence, CES) caused severe damage to land and infrastructure within Christchurch, New Zealand. Approximately one third of the total CES-induced financial losses were directly attributable to liq- uefaction and thus highlights the need for local and regional authorities to assess liquefaction hazards for present and future developments. This thesis is the first to conduct paleo-liquefaction studies in eastern Christchurch for the purpose of de- termining approximate return times of liquefaction-inducing earthquakes within the region. The research uncovered evidence for pre-CES liquefaction dated by radiocarbon and cross-cutting relationships as post-1660 to pre-1905. Additional paleo-liquefaction investigations within the eastern Christchurch suburb of Avon- dale, and the northern township of Kaiapoi, revealed further evidence for pre-CES liquefaction. Pre-CES liquefaction in Avondale is dated as post-1321 and pre-1901, while the Kaiapoi features likely formed during three distinct episodes: post-1458 and possibly during the 1901 Cheviot earthquake, post-1297 to pre-1901, and pre-1458. Evaluation of the liquefaction potential of active faults within the Can- terbury region indicates that many faults have the potential to cause widespread liquefaction within Avondale and Kaiapoi. The identification of pre-CES liquefac- tion confirms that these areas have previously liquefied, and indicates that residen- tial development in eastern Christchurch between 1860 and 2005 occurred in areas containing geologic evidence for pre-CES liquefaction. Additionally, on the basis of detailed field and GIS-based mapping and geospatial-statistical analysis, the distribution and severity of liquefaction and lateral spreading within the eastern Christchurch suburb of Avonside is shown in this study to be strongly in uenced by geomorphic and topographic variability. This variability is not currently ac- counted for in site-specific liquefaction assessments nor the simplified horizontal displacement models, and accounts for some of the variability between the pre- dicted horizontal displacements and those observed during the CES. This thesis highlights the potential applications of paleo-liquefaction investigations and ge- omorphic mapping to seismic and liquefaction hazard assessments and may aid future land-use planning decisions.
The purpose of this assessment is to compare records of known inanga spawning sites in the waterways of Ōtautahi Christchurch from before and after the Canterbury earthquakes, with particular emphasis on information used in the design of planning methods for spawning site protection.
This paper presents site-specific and spatially-distributed ground-motion intensity estimates which have been utilized in the aftermath of the 2010-2011 Canterbury, New Zealand earthquakes. The methodology underpinning the ground motion intensity estimation makes use of both prediction models for ground motion intensity and its within-event spatial correlation. A key benefit of the methodology is that the estimated ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and standard deviation, with the standard deviation being a function of the distance to nearby observations at strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence, which among other things, have been utilized for assessing liquefaction triggering susceptibility of land in residential areas. Secondly, the conditional distribution of response spectral ordinates is obtained at the location of the Canterbury Television building (CTV), which catastrophically collapsed in the 22 February 2011 earthquake. The conditional response spectra provide insight for the selection of ground motion records for use in forensic seismic response analyses of important structures at locations where direct recordings are absent.
Over 900 buildings in the Christchurch central business district and 10,000 residential homes were demolished following the 22nd of February 2011 Canterbury earthquake, significantly disrupting the rebuild progress. This study looks to quantify the time required for demolitions during this event which will be useful for future earthquake recovery planning. This was done using the Canterbury Earthquake Recovery Authority (CERA) demolition database, which allowed an in-depth look into the duration of each phase of the demolition process. The effect of building location, building height, and the stakeholder which initiated the demolition process (i.e. building owner or CERA) was investigated. The demolition process comprises of five phases; (i) decision making, (ii) procurement and planning, (iii) demolition, (iv) site clean-up, and (v) completion certification. It was found that the time required to decide to demolish the building made up majority of the total demolition duration. Demolition projects initiated by CERA had longer procurement and planning durations, but was quicker in other phases. Demolished buildings in the suburbs had a longer decision making duration, but had little effect on other phases of the demolition process. The decision making and procurement and planning phases of the demolition process were shorter for taller buildings, though the other phases took longer. Fragility functions for the duration of each phase in the demolition process are provided for the various categories of buildings for use in future studies.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
Developing a holistic understanding of social, cultural, and economic impacts of disasters can help in building disaster risk knowledge for policy making and planning. Many methods can help in developing an understanding of the impacts of a disaster, including interviews and surveys with people who have experienced disaster, which may be invasive at times and create stress for the participants to relive their experiences. In the past decade, social media, blog posts, video blogs (i.e. “vlogs”), and crowdsourcing mechanisms such as Humanitarian OpenStreetMap and Ushahidi, have become prominent platforms for people to share their experiences and impacts of an event from the ground. These platforms allow for the discovery of a range of impact information, from physical impacts, to social, cultural, and psychological impacts. It can also reveal interesting behavioural information such as their decision to heed a warning or not, as people tend to share their experiences and their reactions online. This information can help researchers and authorities understand both the impacts as well as behavioural responses to hazards, which can then shape how early warning systems are designed and delivered. It can also help to identify gaps in desired behavioural responses. This poster presents a selection of cases identified from the literature and grey literature, such as the Haiti earthquake, the Christchurch earthquake, Hurricane Sandy, and Hurricane Harvey, where online platforms were widely used during and after a disaster to document impacts, experiences, and behavioural responses. A summary of key learnings and areas for future research is provided.
Seismic isolation is an effective technology for significantly reducing damage to buildings and building contents. However, its application to light-frame wood buildings has so far been unable to overcome cost and technical barriers such as susceptibility of light-weight buildings to movement under high-wind loading. The 1994 Northridge Earthquake (6.7 MW) in the United States, 1995 Kobe Earthquake (6.9 MW) in Japan and 2011 Christchurch Earthquake (6.7 Mw) all highlighted significant loss to light-frame wood buildings with over half of earthquake recovery costs allocated to their repair and reconstruction. This poster presents a value case to highlight the benefits of seismically isolated residential buildings compared to the standard fixed-base dwellings for the Wellington region. Loss data generated by insurance claim information from the 2011 Christchurch Earthquake has been used to determine vulnerability functions for the current light-frame wood building stock. By using a simplified single degree of freedom (SDOF) building model, methods for determining vulnerability functions for seismic isolated buildings are developed. Vulnerability functions are then applied directly in a loss assessment to determine the Expected Annual Loss. Vulnerability was shown to dramatically reduce for isolated buildings compared to an equivalent fixed-base building resulting in significant monetary savings, justifying the value case. A state-of-the-art timber modelling software, Timber3D, is then used to model a typical residential building with and without seismic isolation to assess the performance of a proposed seismic isolation system which addresses the technical and cost issues.
Timber-based hybrid structures provide a prospective solution for utilizing environmentally friendly timber material in the construction of mid-rise or high-rise structures. This study mainly focuses on structural damage evaluation for a type of timber-steel hybrid structures, which incorporate prefabricated light wood frame shear walls into steel moment-resisting frames (SMRFs). The structural damage of such a hybrid structure was evaluated through shake table tests on a four-story large-scale timber-steel hybrid structure. Four ground motion records (i.e., Wenchuan earthquake, Canterbury earthquake, El-Centro earthquake, and Kobe earthquake) were chosen for the tests, with the consideration of three different probability levels (i.e., minor, moderate and major earthquakes) for each record. During the shake table tests, the hybrid structure performed quite well with visual damage only to wood shear walls. No visual damage in SMRF and the frame-to-wall connections was observed. The correlation of visual damage to seismic intensity, modal-based damage index and inter-story drift was discussed. The reported work provided a basis of knowledge for performance-based seismic design (PBSD) for such timber-based hybrid structures.
The 14 November 2016 Kaikōura earthquake had major impacts on New Zealand's transport system. Road, rail and port infrastructure was damaged, creating substantial disruption for transport operators, residents, tourists, and business owners in the Canterbury, Marlborough and Wellington regions, with knock-on consequences elsewhere. During both the response and recovery phases, a large amount of information and data relating to the transport system was generated, managed, analysed, and exchanged within and between organisations to assist decision making. To improve information and data exchanges and related decision making in the transport sector during future events and guide new resilience strategies, we present key findings from a recent post-earthquake assessment. The research involved 35 different stakeholder groups and was conducted for the Ministry of Transport. We consider what transport information was available, its usefulness, where it was sourced from, mechanisms for data transfer between organisations, and suggested approaches for continued monitoring.
Research Report No.2010-03Ground motion prediction equations (GMPEs) for geometric-mean pseudo-spectral acceleration amplitudes from New Zealand (NZ) earthquakes are developed. A database of 2437 three-component ground motion records is developed by applying stringent quality criteria to the historically recorded events in NZ. Despite the large number of records, the database is deficient in empirical records from large magnitude events recorded at close distances to the fault rupture plane. As a result, the basis for the NZ-specific GMPE development is to examine the applicability of foreign GMPEs for similar tectonic regions and then modify the most applicable GMPEs based on both theoretical and statistically significant empirically-driven arguments. For active shallow crustal events, five different GMPEs are considered. It was found that the McVerry et al. (2006) model, which is the current model upon which seismic design guidelines and site-specific seismic hazard analyses in NZ are based, provided the worst fit to the NZ database, and that the Chiou et al. (2010) (C10) modification of the Chiou and Youngs (2008) model was the most applicable. Discrepancies between the C10 model and the NZ database that were empirically identified and theoretically justified were used to modify the C10 model for: (i) small magnitude scaling; (ii) scaling of short period ground motion from normal faulting events in volcanic crust; (iii) scaling of ground motions on very hard rock sites; (iv) anelastic attenuation in the NZ crust; and (v) consideration of the increased anelastic attenuation in the Taupo Volcanic Zone (TVZ). For subduction slab events, initially three models were considered. It was found that all of the models had some significant biases with respect to applicability for NZ. The Zhao et al. (2006) (Z06) model was selected because of the rigorous database upon which it was developed and modified by: (i) NZ-specific scaling at small magnitudes; (ii) path scaling at large distances; (iii) consideration of the increased TVZ attenuation; and (iv) revision of the standard deviation model. Based on these modifications the developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The standard deviation of the residuals using the revised standard deviation model also indicated that the model has an adequate precision. Three GMPEs were considered for subduction interface events. The Zhao et al. (2006) (Z06) model was the best performing model with only bias exhibited in the site response model, and possible over-prediction of large magnitude events. The Z06 interface model was modified to account for site response and magnitude scaling using the same functional forms as those of the developed active shallow crustal and subduction slab models. The developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The developed GMPEs include specific features as evident in the NZ database; consistent scaling for parameters not well constrained by the NZ database; and pseudo-spectral amplitudes for vibration periods from 0.01 to 10 seconds. Hence, these models represent a significant advance in the state-of-the art for empirical ground motion prediction in NZ.
A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted 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, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.
Base isolation is arguably the most reliable method for providing enhanced protection of buildings against earthquake-induced actions, by virtue of a physical separation between the structure and the ground through elements/devices with controlled force capacity, significant lateral deformation capacity and (often) enhanced energy dissipation. Such a design solution has shown its effectiveness in protecting both structural and non-structural components, hence preserving their functionality even in the aftermath of a major seismic event. Despite lead rubber bearings being invented in New Zealand almost forty years ago, the Christchurch Women's hospital was the only isolated building in Christchurch when the Canterbury earthquake sequence struck in 2010/11. Furthermore, a reference code for designing base-isolated buildings in New Zealand is still missing. The absence of a design standard or at least of a consensus on design guidelines is a potential source for a lack of uniformity in terms of performance criteria and compliance design approaches. It may also limit more widespread use of the technology in New Zealand. The present paper provides an overview of the major international codes (American, Japanese and European) for the design of base-isolated buildings. The design performance requirements, the analysis procedures, the design review process and approval/quality control of devices outlined in each code are discussed and their respective pros and cons are compared through a design application on a benchmark building in New Zealand. The results gathered from this comparison are intended to set the basis for the development of guidelines specific for the New Zealand environment.
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.
