In practice, several competing liquefaction evaluation procedures (LEPs) are used to compute factors of safety against soil liquefaction, often for use within a liquefaction potential index (LPI) framework to assess liquefaction hazard. At present, the influence of the selected LEP on the accuracy of LPI hazard assessment is unknown, and the need for LEP-specific calibrations of the LPI hazard scale has never been thoroughly investigated. Therefore, the aim of this study is to assess the efficacy of three CPT-based LEPs from the literature, operating within the LPI framework, for predicting the severity of liquefaction manifestation. Utilising more than 7000 liquefaction case studies from the 2010–2011 Canterbury (NZ) earthquake sequence, this study found that: (a) the relationship between liquefaction manifestation severity and computed LPI values is LEP-specific; (b) using a calibrated, LEP-specific hazard scale, the performance of the LPI models is essentially equivalent; and (c) the existing LPI framework has inherent limitations, resulting in inconsistent severity predictions against field observations for certain soil profiles, regardless of which LEP is used. It is unlikely that revisions of the LEPs will completely resolve these erroneous assessments. Rather, a revised index which more adequately accounts for the mechanics of liquefaction manifestation is needed.
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.
Liquefaction-induced lateral spreading during the 2011 Christchurch earthquake in New Zealand was severe and extensive, and data regarding the displacements associated with the lateral spreading provides an excellent opportunity to better understand the factors that influence these movements. Horizontal displacements measured from optical satellite imagery and subsurface data from the New Zealand Geotechnical Database (NZGD) were used to investigate four distinct lateral spread areas along the Avon River in Christchurch. These areas experienced displacements between 0.5 and 2 m, with the inland extent of displacement ranging from 100 m to over 600 m. Existing empirical and semi-empirical displacement models tend to under estimate displacements at some sites and over estimate at others. The integrated datasets indicate that the areas with more severe and spatially extensive displacements are associated with thicker and more laterally continuous deposits of liquefiable soil. In some areas, the inland extent of displacements is constrained by geologic boundaries and geomorphic features, as expressed by distinct topographic breaks. In other areas the extent of displacement is influenced by the continuity of liquefiable strata or by the presence of layers that may act as vertical seepage barriers. These observations demonstrate the need to integrate geologic/geomorphic analyses with geotechnical analyses when assessing the potential for lateral spreading movements.
A magnitude 6.3 earthquake struck the city of Christchurch at 12:51pm on Tuesday 22 February 2011. The earthquake caused 182 fatalities, a large number of injuries, and resulted in widespread damage to the built environment, including significant disruption to the lifelines. The event created the largest lifeline disruption in a New Zealand city in 80 years, with much of the damage resulting from extensive and severe liquefaction in the Christchurch urban area. The Christchurch earthquake occurred when the Canterbury region and its lifelines systems were at the early stage of recovering from the 4 September 2010 Darfield (Canterbury) magnitude 7.1 earthquake. This paper describes the impact of the Christchurch earthquake on lifelines by briefly summarising the physical damage to the networks, the system performance and the operational response during the emergency management and the recovery phase. Special focus is given to the performance and management of the gas, electric and road networks and to the liquefaction ejecta clean-up operations that contributed to the rapid reinstatement of the functionality of many of the lifelines. The water and wastewater system performances are also summarized. Elements of resilience that contributed to good network performance or to efficient emergency and recovery management are highlighted in the paper.
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.
Study region: Christchurch, New Zealand. Study focus: Low-lying coastal cities worldwide are vulnerable to shallow groundwater salinization caused by saltwater intrusion and anthropogenic activities. Shallow groundwater salinization can have cascading negative impacts on municipal assets, but this is rarely considered compared to impacts of salinization on water supply. Here, shallow groundwater salinity was sampled at high spatial resolution (1.3 piezometer/km2 ), then mapped and spatially interpolated. This was possible due to a uniquely extensive set of shallow piezometers installed in response to the 2010–11 Canterbury Earthquake Sequence to assess liquefaction risk. The municipal assets located within the brackish groundwater areas were highlighted. New hydrological insights for the region: Brackish groundwater areas were centred on a spit of coastal sand dunes and inside the meander of a tidal river with poorly drained soils. The municipal assets located within these areas include: (i) wastewater and stormwater pipes constructed from steel-reinforced concrete, which, if damaged, are vulnerable to premature failure when exposed to chloride underwater, and (ii) 41 parks and reserves totalling 236 ha, within which salt-intolerant groundwater-dependent species are at risk. This research highlights the importance of determining areas of saline shallow groundwater in low-lying coastal urban settings and the co-located municipal assets to allow the prioritisation of sites for future monitoring and management.
Semi-empirical models based on in-situ geotechnical tests have become the standard of practice for predicting soil liquefaction. Since the inception of the “simplified” cyclic-stress model in 1971, variants based on various in-situ tests have been developed, including the Cone Penetration Test (CPT). More recently, prediction models based soley on remotely-sensed data were developed. Similar to systems that provide automated content on earthquake impacts, these “geospatial” models aim to predict liquefaction for rapid response and loss estimation using readily-available data. This data includes (i) common ground-motion intensity measures (e.g., PGA), which can either be provided in near-real-time following an earthquake, or predicted for a future event; and (ii) geospatial parameters derived from digital elevation models, which are used to infer characteristics of the subsurface relevent to liquefaction. However, the predictive capabilities of geospatial and geotechnical models have not been directly compared, which could elucidate techniques for improving the geospatial models, and which would provide a baseline for measuring improvements. Accordingly, this study assesses the realtive efficacy of liquefaction models based on geospatial vs. CPT data using 9,908 case-studies from the 2010-2016 Canterbury earthquakes. While the top-performing models are CPT-based, the geospatial models perform relatively well given their simplicity and low cost. Although further research is needed (e.g., to improve upon the performance of current models), the findings of this study suggest that geospatial models have the potential to provide valuable first-order predictions of liquefaction occurence and consequence. Towards this end, performance assessments of geospatial vs. geotechnical models are ongoing for more than 20 additional global earthquakes.
The Resilient Organisations Research Programme and the University of Canterbury are undertaking a longitudinal study to examine the resilience and recovery of organisations within the Canterbury region following the 4 September Canterbury earthquake. The preliminary data suggest the physical, economic and social effects of the earthquake were varied across industry sectors within Canterbury. These preliminary results catalogue organisations’ perceptions of the: - disruptions to their ability to do business - challenges faced in the aftermath of the earthquake - factors that have helped mitigate the effects of the earthquake - revenue changes and projections for the duration of this change - financing options for recovery
The purpose of this study is to analyse the felt earthquake impacts, resilience and recovery of organizations in Canterbury by comparing three business sectors (accommodation/food services, Education/Training and Manufacturing). A survey of the three sectors in 2013 of Canterbury organizations impacted by the earthquakes revealed significant differences between the three sectors on felt earthquake impacts and resilience. On recovery and mitigation factors, the accommodation/food services sector is not significantly different from the other two sectors. Overall, the survey results presented here indicate that the Accommodation/Food Services sector was the least impacted by the earthquakes in comparison to the Education/Training and Manufacturing sectors. Implications for post-disaster management and recovery of the accommodation sector are suggested.
Abstract This study provides a simplified methodology for pre-event data collection to support a faster and more accurate seismic loss estimation. Existing pre-event data collection frameworks are reviewed. Data gathered after the Canterbury earthquake sequences are analysed to evaluate the relative importance of different sources of building damage. Conclusions drawns are used to explore new approaches to conduct pre-event building assessment.
Detailed studies on the sediment budget may reveal valuable insights into the successive build-up of the Canterbury Plains and their modification by Holocene fluvialaction connected to major braided rivers. Additionally, they bear implications beyond these fluvial aspects. Palaeoseismological studies claim to have detected signals of major Alpine Fault earthquakes in coastal environments along the eastern seaboard of the South Island (McFadgen and Goff, 2005). This requires high connectivity between the lower reaches of major braided rivers and their mountain catchments to generate immediate significant sediment pulses. It would be contradictory to the above mentioned hypothesis though. Obtaining better control on sediment budgets of braided rivers like the Waimakariri River will finally add significant value to multiple scientific and applied topics like regional resource management. An essential first step of sediment budget studies Is to systematically map the geomorphology, conventionally in the field and/or using remote-sensing applications, to localise, genetically identify, and classify landforms or entire toposequences of the area being investigated. In formerly glaciated mountain environments it is also indispensable to obtain all available chronological information supporting subsequent investigations.
There has not been substantial research conducted in the area of fraud and natural disasters. Therefore, this study sought to examine the perceptions of Canterbury residents toward the recovery process following the September 2010 and February 2011 earthquakes and whether residents felt as though contractor fraud occurs in Canterbury. A questionnaire was developed to gauge information about Canterbury residents’ self-reports involving the earthquakes, specific contractors involved, parties involved with the recovery process in general, and demographic information. Participants included a total of 213 residents from the Canterbury region who had been involved with contractors and/or insurance companies due to the recovery process. Results indicated that a high percentage of the participants were not satisfied with the recovery process and that almost half of the participants reported feeling scammed by contractors in Canterbury after the 2010 and 2011 earthquakes. Moreover, the results indicate that participants neither agreed with the assessments made about their property losses nor the plans made to recover their properties. In many cases, participants felt pressured and even reluctant to accept these assessments and/or plans. The present study does not seek to explain why contractor fraud exists or what motivates scammers. Conversely, it attempts to demonstrate the perceptions of contractor fraud and satisfaction that have taken place in the aftermath of the Canterbury earthquakes.
Very little research exists on total house seismic performance. This testing programme provides stiffness and response data for five houses of varying ages including contributions of non-structural elements. These light timber framed houses in Christchurch, New Zealand had minor earthquake damage from the 2011 earthquakes and were lateral load tested on site to determine their strength and stiffness, and preliminary damage thresholds. Dynamic characteristics were also investigated. Various loading schemes were utilised including quasi-static loading above the foundation, unidirectional loading through the floor diaphragm, cyclic quasi-static loading and snapback tests. Dynamic analysis on two houses provided the seismic safety levels of post-quake houses with respect to local hazard levels. Compared with New Zealand Building Standards all the tested houses had an excess of strength, damage is a significant consideration in earthquake resilience and was observed in all of the houses. A full size house laboratory test is proposed.
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.
This paper presents the preliminary conclusions of the first stage of Wellington Case Study project (Regulating For Resilience in an Earthquake Vulnerable City) being undertaken by the Disaster Law Research Group at the University of Canterbury Law School. This research aims to map the current regulatory environment around improving the seismic resilience of the urban built environment. This work provides the basis for the second stage of the project which will map the regulatory tools onto the reality of the current building stock in Wellington. Using a socio-legal methodology, the current research examines the regulatory framework around seismic resilience for existing buildings in New Zealand, with a particularly focus on multi-storey in the Wellington CBD. The work focusses both on the operation and impact of the formal seismic regulatory tools open to public regulators (under the amended Building Act) as other non-seismic regulatory tools. As well as examining the formal regulatory frame, the work also provides an assessment of the interactions between other non-building acts (such as Health and Safety at Work Act 2015) on the requirements of seismic resilience. Other soft-law developments (particularly around informal building standards) are also examined. The final output of this work will presents this regulatory map in a clear and easily accessible manner and provide an assessment of the suitability of this at times confusing and patchy legal environment as Wellington moves towards becoming a resilient city. The final conclusion of this work will be used to specifically examine the ability of Wellington to make this transition under the current regulatory environment as phase two of the Wellington Case Study project.
- The Avon-Ōtākaro Redzone is an 11 kilometer stretch of land along the Avon-Ōtākaro River in Christchurch. - This project focused on the creation of a publicly available biodiversity map of the AvonŌtākaro River Corridor, a project undertaken as part of the ecological restoration of the Christchurch redzone. - This project originated from the Christchurch 2010-2011 earthquake sequence which saw liquefaction damage along 11km of the Avon River. Under guidance from The Nature Lab & Ōtākaro Living Laboratory, and various other experts, the primary research objective was to map historical biodiversity, identify hotspots, and assess areas for potential revegetation. - The data collected came from historical black maps, current iNaturalist data, and soil classification information. - The findings show that, pre-colonialism, the area was composed of herbaceous areas, wetlands, native shrubland, and tussock land, with key plants such as river, fern, tutu, and cabbage trees. - The post-earthquake analysis shows a transition from a residential area to patchy grasslands and swampy areas. - The findings also showed a strong relationship between historic sites and soil classifications, providing knowledge for past and future vegetation patterns and spread. - This map will be a valuable resource for conservation efforts and public engagement as the area transitions into a blue-green corridor.
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.
This paper presents the preliminary findings of a study on the resilience and recovery of organisations following the Darfield earthquake in New Zealand on 4 September 2010. Sampling included organisations proximal and distal to the fault trace, organisations located within central business districts, and organisations from seven diverse industry sectors. The research captured information on the challenges to, the impacts on, and the reflections of the organisations in the first months of recovery. Organisations in central business districts and in the hospitality sector were most likely to close while organisations that had perishable stock and livestock were more heavily reliant on critical services. Staff well-being, cash flow, and customer loss were major concerns for organisations across all sectors. For all organisations, the most helpful factors in mitigating the effects of the earthquake to be their relationship with staff, the design and type of buildings, and critical service continuity or swift reinstatement of services.
This is an interim report from the research study performed within the NHRP Research Project “Impacts of soil liquefaction on land, buildings and buried pipe networks: geotechnical evaluation and design, Project 3: Seismic assessment and design of pipe networks in liquefiable soils”. The work presented herein is a continuation of the comprehensive study on the impacts of Christchurch earthquakes on the buried pipe networks presented in Cubrinovski et al. (2011). This report summarises the performance of Christchurch City’s potable water, waste water and road networks through the 2010-2011 Canterbury Earthquake Sequence (CES), and particularly focuses on the potable water network. It combines evidence based on comprehensive and well-documented data on the damage to the water network, detailed observations and interpretation of liquefaction-induced land damage, records and interpretations of ground motion characteristics induced by the Canterbury earthquakes, for a network analysis and pipeline performance evaluation using a GIS platform. The study addresses a range of issues relevant in the assessment of buried networks in areas affected by strong earthquakes and soil liquefaction. It discusses performance of different pipe materials (modern flexible pipelines and older brittle pipelines) including effects of pipe diameters, fittings and pipeline components/details, trench backfill characteristics, and severity of liquefaction. Detailed breakdown of key factors contributing to the damage to buried pipes is given with reference to the above and other relevant parameters. Particular attention is given to the interpretation, analysis and modelling of liquefaction effects on the damage and performance of the buried pipe networks. Clear link between liquefaction severity and damage rate for the pipeline has been observed with an increasing damage rate seen with increasing liquefaction severity. The approach taken here was to correlate the pipeline damage to LRI (Liquefaction Resistance Index, newly developed parameter in Cubrinovski et al., 2011) which represents a direct measure for the soil resistance to liquefaction while accounting for the seismic demand through PGA. Key quality of the adopted approach is that it provides a general methodology that in conjunction with conventional methods for liquefaction evaluation can be applied elsewhere in New Zealand and internationally. Preliminary correlations between pipeline damage (breaks km-1), liquefaction resistance (LRI) and seismic demand (PGA) have been developed for AC pipes, as an example. Such correlations can be directly used in the design and assessment of pipes in seismic areas both in liquefiable and non-liquefiable areas. Preliminary findings on the key factors for the damage to the potable water pipe network and established empirical correlations are presented including an overview of the damage to the waste water and road networks but with substantially less detail. A comprehensive summary of the damage data on the buried pipelines is given in a series of appendices.
The latest two great earthquake sequences; 2010- 2011 Canterbury Earthquake and 2016 Kaikoura Earthquake, necessitate a better understanding of the New Zealand seismic hazard condition for new building design and detailed assessment of existing buildings. It is important to note, however, that the New Zealand seismic hazard map in NZS 1170.5.2004 is generalised in effort to cover all of New Zealand and limited to a earthquake database prior to 2001. This is “common” that site-specific studies typically provide spectral accelerations different to those shown on the national map (Z values in NZS 1170.5:2004); and sometimes even lower. Moreover, Section 5.2 of Module 1 of the Earthquake Geotechnical Engineering Practice series provide the guidelines to perform site- specific studies.
This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.
As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows: Ic < 1.31 – Gravelly sand to dense sand 1.31 < Ic < 1.90 – Sands: clean sand to silty sand 1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt 2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay 3.20 < Ic < 3.60 – Clays: silty clay to clay Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.
The nonlinear dynamic soil-foundation-structure interaction (SFSI) can signifi cantly affect the seismic response of buildings, causing additional deformation modes, damage and repair costs. Because of nonlinear foundation behaviour and interactions, the seismic demand on the superstructure may considerably change, and also permanent deformations at the foundation level may occur. Although SFSI effects may be benefi cial to the superstructure performance, any advantage would be of little structural value unless the phenomenon can be reliably controlled and exploited. Detrimental SFSI effects may also occur, including acceleration and displacement response ampli cation and differential settlements, which would be unconservative to neglect. The lack of proper understanding of the phenomenon and the limited available simpli ed tools accounting for SFSI have been major obstacles to the implementation of integrated design and assessment procedures into the everyday practice. In this study concepts, ideas and practical tools (inelastic spectra) for the seismic design and assessment of integrated foundation-superstructure systems are presented, with the aim to explicitly consider the impact of nonlinearities occurring at the soil-foundation interface on the building response within an integrated approach, where the foundation soil and superstructure are considered as part of an integrated system when evaluating the seismic response, working synergically for the achievement of a target global performance. A conceptual performance-based framework for the seismic design and assessment of integrated foundation-superstructure systems is developed. The framework is based on the use of peak and residual response parameters for both the superstructure and the foundation, which are then combined to produce the system performance matrix. Each performance matrix allows for worsening of the performance when different contributions are combined. An attempt is made to test the framework by using case histories from the 2011 Christchurch earthquake, which are previously shown to have been severely affected by nonlinear SFSI. The application highlights the framework sensitivity to the adopted performance limit states, which must be realistic for a reliable evaluation of the system performance. Constant ductility and constant strength inelastic spectra are generated for nonlinear SFSI systems (SDOF nonlinear superstructure and 3DOF foundation allowing for uplift and soil yielding), representing multistorey RC buildings with shallow rigid foundations supported by cohesive soils. Different ductilities/strengths, hysteretic rules (Bi-linear, Takeda and Flag-Shape), soil stiffness and strength and bearing capacity factors are considered. Footings and raft foundations are investigated, characterized respectively by constant (3 and 8) and typically large bearing capacity factors. It is confi rmed that when SFSI is considered, the superstructure yielding force needed to satisfy a target ductility for a new building changes, and that similarly, for an existing building, the ductility demand on a building of a given strength varies. The extent of change of seismic response with respect to xed-base (FB) conditions depends on the class of soils considered, and on the bearing capacity factor (SF). For SF equal to 3, the stiffer soils enhance the nonlinear rotational foundation behaviour and are associated with reduced settlement, while the softer ones are associated with increased settlement response but not signi ficant rotational behaviour. On average terms, for the simplifi ed models considered, SFSI is found to be bene cial to the superstructure performance in terms of acceleration and superstructure displacement demand, although exceptions are recorded due to ground motion variability. Conversely, in terms of total displacement, a signi cant response increase is observed. The larger the bearing capacity factor, the more the SFSI response approaches the FB system. For raft foundation buildings, characterized by large bearing capacity factors, the impact of foundation response is mostly elastic, and the system on average approaches FB conditions. Well de fined displacement participation factors to the peak total lateral displacement are observed for the different contributions (i.e. peak foundation rotation and translation and superstructure displacement). While the superstructure and foundation rotation show compensating trends, the foundation translation contribution varies as a function of the moment-to-shear ratio, becoming negligible in the medium-to-long periods. The longer the superstructure FB period, the less the foundation response is signifi cant. The larger the excitation level and the less ductile the superstructure, the larger the foundation contribution to the total lateral displacement, and the less the superstructure contribution. In terms of hysteretic behaviour, its impact is larger when the superstructure response is more signifi cant, i.e. for the softer/weaker soils and larger ductilities. Particularly, for the Flag Shape rule, larger superstructure displacement participation factors and smaller foundation contributions are recorded. In terms of residual displacements, the total residual-to-maximum ratios are similar in amplitudes and trends to the corresponding FB system responses, with the foundation and superstructure contributions showing complementary trends. The impact of nonlinear SFSI is especially important for the Flag Shape hysteresis rule, which would not otherwise suffer of any permanent deformations. By using the generated peak and residual inelastic spectra (i.e. inelastic acceleration/ displacement modifi cation factor spectra, and/or participation factor and residual spectra), conceptual simplifi ed procedures for the seismic design and assessment of integrated foundation-superstructure systems are presented. The residual displacements at both the superstructure and foundation levels are explicitly considered. Both the force- and displacement-based approaches are explored. The procedures are de fined to be complementary to the previously proposed integrated performance-based framework. The use of participation factor spectra allows the designer to easily visualize the response of the system components, and could assist the decision making process of both the design and assessment of SFSI systems. The presented numerical results have been obtained using simpli ed models, assuming rigid foundation behaviour and neglecting P-Delta effects. The consideration of more complex systems including asymmetry in stiffness, mass, axial load and ground conditions with a exible foundation layout would highlight detrimental SFSI effects as related to induced differential settlements, while accounting for PDelta effects would further amplify the displacement response. Also, the adopted acceleration records were selected and scaled to match conventional design spectra, thus not representing any response ampli cation in the medium-to-long period range which could as well cause detrimental SFSI effects. While these limitations should be the subject of further research, this study makes a step forward to the understanding of SFSI phenomenon and its incorporation into performance-based design/assessment considerations.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
This is a joint Resilience Framework undertaken by the Electrical, Computer and Software Engineering Department of the University of Auckland in association with West Power and Orion networks and partially funded by the New Zealand National Science Challenge and QuakeCoRE. The Energy- Communication research group nearly accomplished two different researches focusing on both asset resilience and system resilience. Asset resilience research which covers underground cables system in Christchurch region is entitled “2010-2011 Canterbury Earthquake Sequence Impact on 11KV Underground Cables” and system resilience research which covers electricity distribution and communication system in West Coast region is entitled “NZ Electricity Distribution Network Resilience Assessment and Restoration Models following Major Natural Disturbance“. As the fourth milestone of the aforementioned research project, the latest outcome of both projects has been socialised with the stakeholders during the Cigre NZ 2019 Forum.
This study contains an evaluation of the seismic hazard associated with the Springbank Fault, a blind structure discovered in 1998 close to Christchurch. The assessment of the seismic hazard is approached as a deterministic process in which it is necessary to establish: 1) fault characteristics; 2) the maximum earthquake that the fault is capable of producing and 3) ground motions estimations. Due to the blind nature of the fault, conventional techniques used to establish the basic fault characteristics for seismic hazard assessments could not be applied. Alternative methods are used including global positioning system (GPS) surveys, morphometric analyses along rivers, shallow seismic reflection surveys and computer modelling. These were supplemented by using multiple empirical equations relating fault attributes to earthquake magnitude, and attenuation relationships to estimate ground motions in the near-fault zone. The analyses indicated that the Springbank Fault is a reverse structure located approximately 30 km to the northwest of Christchurch, along a strike length of approximately 16 km between the Eyre and Ashley River. The fault does not reach the surface, buy it is associated with a broad anticline whose maximum topographic expression offers close to the mid-length of the fault. Two other reverse faults, the Eyrewell and Sefton Faults, are inferred in the study area. These faults, together with the Springbank and Hororata Faults and interpreted as part of a sys of trust/reverse faults propagating from a decollement located at mid-crustal depths of approximately 14 km beneath the Canterbury Plains Within this fault system, the Springbank Fault is considered to behave in a seismically independent way, with a fault slip rate of ~0.2 mm/yr, and the capacity of producing a reverse-slip earthquake of moment magnitude ~6.4, with an earthquake recurrence of 3,000 years. An earthquake of the above characteristics represents a significant seismic hazard for various urban centres in the near-fault zone including Christchurch, Rangiora, Oxford, Amberley, Kaiapoi, Darfield, Rollestion and Cust. Estimated peak ground accelerations for these towns range between 0.14 g to 0.5 g.
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 operation of telecommunication networks is critical during business as usual times, and becomes most vital in post-disaster scenarios, when the services are most needed for restoring other critical lifelines, due to inherent interdependencies, and for supporting emergency and relief management tasks. In spite of the recognized critical importance, the assessment of the seismic performance for the telecommunication infrastructure appears to be underrepresented in the literature. The FP6 QuakeCoRE project “Performance of the Telecommunication Network during the Canterbury Earthquake Sequence” will provide a critical contribution to bridge this gap. Thanks to an unprecedented collaboration between national and international researchers and highly experienced asset managers from Chorus, data and evidences on the physical and functional performance of the telecommunication network after the Canterbury Earthquakes 2010-2011 have been collected and collated. The data will be processed and interpreted aiming to reveal fragilities and resilience of the telecommunication networks to seismic events
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.
