The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.
Organisations play a vital role in assisting communities to recover from disasters. They are the key providers of goods and services needed in both response and recovery efforts. They provide the employment which both anchors people to place and supports the taxation base to allow for necessary recovery spending. Finally, organisations are an integral part of much day to day functioning contributing immensely to people’s sense of ‘normality’ and psychological wellbeing. Yet, despite their overall importance in the recovery process, there are significant gaps in our existing knowledge with regard to how organisations respond and recover following disaster. This research fills one part of this gap by examining collaboration as an adaptive strategy enacted by organisations in the Canterbury region of New Zealand, which was heavily impacted by a series of major earthquakes, occurring in 2010 and 2011. Collaboration has been extensively investigated in a variety of settings and from numerous disciplinary perspectives. However, there are few studies that investigate the role of collaborative approaches to support post-disaster business recovery. This study investigates the type of collaborations that have occurred and how they evolved as organisations reacted to the resource and environmental change caused by the disaster. Using data collected through semi-structured interviews, survey and document analysis, a rich and detailed picture of the recovery journey is created for 26 Canterbury organisations including 14 collaborators, six non-traders, five continued traders and one new business. Collaborations included two or more individual businesses collaborating along with two multi-party, place based projects. Comparative analysis of the organisations’ experiences enabled the assessment of decisions, processes and outcomes of collaboration, as well as insight into the overall process of business recovery. This research adopted a primarily inductive, qualitative approach, drawing from both grounded theory and case study methodologies in order to generate theory from this rich and contextually situated data. Important findings include the importance of creating an enabling context which allows organisations to lead their own recovery, the creation of a framework for effective post-disaster collaboration and the importance of considering both economic and other outcomes. Collaboration is found to be an effective strategy enabling resumption of trade at a time when there seemed few other options available. While solving this need, many collaborators have discovered significant and unexpected benefits not just in terms of long term strategy but also with regard to wellbeing. Economic outcomes were less clear-cut. However, with approximately 70% of the Central Business District demolished and rebuilding only gaining momentum in late 2014, many organisations are still in a transition stage moving towards a new ‘normal’.
Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.
There is a relationship between inelastic deformation and energy dissipation in structures that are subjected to earthquake ground motions. Thus, if seismic energy dissipation can be achieved by means of a separate non-load bearing supplementary damping system, the load bearing structure can remain elastic with continuing serviceability following the design level earthquake. This research was carried out to investigate the advantages of using added damping in structures. The control system consists of passive friction dampers called ring spring dampers installed in the ground floor of the structure using a tendon to transmit the forces to the other parts of the structure. The ring springs dampers are friction devices consisting of inner and outer ring elements assembled to form a spring stack. External load applied to the spring produces sliding action across mating ring interfaces. The damping forces generated by the dampers and transferred in the supplemental system to the structure by the tendon and horizontal links oppose the internal loads. A four storey-two bay steel frame structure was used in the study. Experimental and analytical studies to investigate the effectiveness of a supplemental control system are presented. The model was subjected to a series of earthquake simulations on the shaking table in the Structural Laboratory of the Civil Engineering Department, at the University of Canterbury. The earthquake simulation tests have been performed on the structure both with and without the supplemental control system. The earthquake simulations were a series of gradually increasing intensity replications of two commonly used earthquake records. This thesis includes detailed description of the structural model, the supplemental control system, the ring springs dampers and the data obtained during the testing. Analyses were then carried out on a twelve storey framed structure to investigate the possible tendon arrangements and the size and type of dampers required to control the response of a real building. Guidelines for determining the appropriateness of including a supplemental damping system have been investigated. The main features of the supplemental control system adopted in this research are: • It is a passive control system with extreme reliability and having no dependence on external power sources to effect the control action. These power sources may not be available during a major earthquake. • Ring springs are steel friction devices capable of absorbing large amounts of input energy. No liquid leakage can occur and minimal maintenance is required for the ring spring dampers. • With a damper-tendon system, the distribution of the dampers throughout the structure is not so critical. Only one or two dampers are used to produce the damping forces needed, and forces are then transferred to the rest of the building by the tendon system. • It is a relatively inexpensive control system with a long useful life.
In recent years, rocking isolation has become an effective approach to improve seismic performance of steel and reinforced concrete structures. These systems can mitigate structural damage through rigid body displacement and thus relatively low requirements for structural ductility, which can significantly improve seismic resilience of structures and reduce repairing costs after strong earthquakes. A number of base rocking structural systems with only a single rocking interface have been proposed. However, these systems can have significant high mode effect for high rise structures due to the single rocking interface. This RObust BUilding SysTem (ROBUST) project is a collaborative China-New Zealand project sponsored by the International Joint Research Laboratory of Earthquake Engineering (ILEE), Tongji University, and a number of agencies and universities within New Zealand including the BRANZ, Comflor, Earthquake Commission, HERA, QuakeCoRE, QuakeCentre, University of Auckland, and the University of Canterbury. A number of structural configurations will be tested [1, 2], and non-structural elements including ceilings, infilling walls, glazed curtain walls, precast concrete panels, piping system will also be tested in this project [3]. Within this study, a multiple rocking column steel structural system was proposed and investigated mainly by Tongji team with assistance of NZ members. The concept of rocking column system initiates from the structure of Chinese ancient wooden pagoda. In some of Chinese wooden pagodas, there are continuous core columns hanged only at the top of each pagoda, which is not connected to each stories. This core column can effectively avoid collapse of the whole structure under large storey drifts. Likewise, there are also central continuous columns in the newly proposed steel rocking column system, which can avoid weak story failure mechanism and make story drifts more uniform. In the proposed rocking column system, the structure can switch between an elastic rigidly connected moment resisting frame and a controlled rocking column system when subjected to strong ground motion excitations. The main seismic energy can be dissipated by asymmetric friction beam–column connections, thereby effectively reducing residual displacement of the structure under seismic loading without causing excessive damage to structural members. Re–centering of the structure is provided not only by gravity load carried by rocking columns, but also by mould coil springs. To investigate dynamic properties of the proposed system under different levels of ground excitations, a full-scale threestory steel rocking column structural system with central continuous columns is to be tested using the International joint research Laboratory of Earthquake Engineering (ILEE) facilities, Shanghai, China and an analytical model is established. A finite element model is also developed using ABAQUS to simulate the structural dynamic responses. The rocking column system proposed in this paper is shown to produce resilient design with quick repair or replacement.
In September 2010 and February 2011 the Canterbury region of New Zealand was struck by two powerful earthquakes, registering magnitude 7.1 and 6.3 respectively on the Richter scale. The second earthquake was centred 10 kilometres south-east of the centre of Christchurch (the region’s capital and New Zealand’s third most populous urban area, with approximately 360,000 residents) at a depth of five kilometres. 185 people were killed, making it the second deadliest natural disaster in New Zealand’s history. (66 people were killed in the collapse of one building alone, the six-storey Canterbury Television building.) The earthquake occurred during the lunch hour, increasing the number of people killed on footpaths and in buses and cars by falling debris. In addition to the loss of life, the earthquake caused catastrophic damage to both land and buildings in Christchurch, particularly in the central business district. Many commercial and residential buildings collapsed in the tremors; others were damaged through soil liquefaction and surface flooding. Over 1,000 buildings in the central business district were eventually demolished because of safety concerns, and an estimated 70,000 people had to leave the city after the earthquakes because their homes were uninhabitable. The New Zealand Government declared a state of national emergency, which stayed in force for ten weeks. In 2014 the Government estimated that the rebuild process would cost NZ$40 billion (approximately US$27.3 billion, a cost equivalent to 17% of New Zealand’s annual GDP). Economists now estimate it could take the New Zealand economy between 50 and 100 years to recover. The earthquakes generated tens of thousands of insurance claims, both against private home insurance companies and against the New Zealand Earthquake Commission, a government-owned statutory body which provides primary natural disaster insurance to residential property owners in New Zealand. These ranged from claims for hundreds of millions of dollars concerning the local port and university to much smaller claims in respect of the thousands of residential homes damaged. Many of these insurance claims resulted in civil proceedings, caused by disputes about policy cover, the extent of the damage and the cost and/or methodology of repairs, as well as failures in communication and delays caused by the overwhelming number of claims. Disputes were complicated by the fact that the Earthquake Commission provides primary insurance cover up to a monetary cap, with any additional costs to be met by the property owner’s private insurer. Litigation funders and non-lawyer claims advocates who took a percentage of any insurance proceeds also soon became involved. These two factors increased the number of parties involved in any given claim and introduced further obstacles to resolution. Resolving these disputes both efficiently and fairly was (and remains) central to the rebuild process. This created an unprecedented challenge for the justice system in Christchurch (and New Zealand), exacerbated by the fact that the Christchurch High Court building was itself damaged in the earthquakes, with the Court having to relocate to temporary premises. (The High Court hears civil claims exceeding NZ$200,000 in value (approximately US$140,000) or those involving particularly complex issues. Most of the claims fell into this category.) This paper will examine the response of the Christchurch High Court to this extraordinary situation as a case study in innovative judging practices and from a jurisprudential perspective. In 2011, following the earthquakes, the High Court made a commitment that earthquake-related civil claims would be dealt with as swiftly as the Court's resources permitted. In May 2012, it commenced a special “Earthquake List” to manage these cases. The list (which is ongoing) seeks to streamline the trial process, resolve quickly claims with precedent value or involving acute personal hardship or large numbers of people, facilitate settlement and generally work proactively and innovatively with local lawyers, technical experts and other stakeholders. For example, the Court maintains a public list (in spreadsheet format, available online) with details of all active cases before the Court, listing the parties and their lawyers, summarising the facts and identifying the legal issues raised. It identifies cases in which issues of general importance have been or will be decided, with the expressed purpose being to assist earthquake litigants and those contemplating litigation and to facilitate communication among parties and lawyers. This paper will posit the Earthquake List as an attempt to implement innovative judging techniques to provide efficient yet just legal processes, and which can be examined from a variety of jurisprudential perspectives. One of these is as a case study in the well-established debate about the dialogic relationship between public decisions and private settlement in the rule of law. Drawing on the work of scholars such as Hazel Genn, Owen Fiss, David Luban, Carrie Menkel-Meadow and Judith Resnik, it will explore the tension between the need to develop the law through the doctrine of precedent and the need to resolve civil disputes fairly, affordably and expeditiously. It will also be informed by the presenter’s personal experience of the interplay between reported decisions and private settlement in post-earthquake Christchurch through her work mediating insurance disputes. From a methodological perspective, this research project itself gives rise to issues suitable for discussion at the Law and Society Annual Meeting. These include the challenges in empirical study of judges, working with data collected by the courts and statistical analysis of the legal process in reference to settlement. September 2015 marked the five-year anniversary of the first Christchurch earthquake. There remains widespread dissatisfaction amongst Christchurch residents with the ongoing delays in resolving claims, particularly insurers, and the rebuild process. There will continue to be challenges in Christchurch for years to come, both from as-yet unresolved claims but also because of the possibility of a new wave of claims arising from poor quality repairs. Thus, a final purpose of presenting this paper at the 2016 Meeting is to gain the benefit of other scholarly perspectives and experiences of innovative judging best practice, with a view to strengthening and improving the judicial processes in Christchurch. This Annual Meeting of the Law and Society Association in New Orleans is a particularly appropriate forum for this paper, given the recent ten year anniversary of Hurricane Katrina and the plenary session theme of “Natural and Unnatural Disasters – human crises and law’s response.” The presenter has a personal connection with this theme, as she was a Fulbright scholar from New Zealand at New York University in 2005/2006 and participated in the student volunteer cleanup effort in New Orleans following Katrina. http://www.lawandsociety.org/NewOrleans2016/docs/2016_Program.pdf
Background: Up to 6 years after the 2011 Christchurch earthquakes, approximately one-third of parents in the Christchurch region reported difficulties managing the continuously high levels of distress their children were experiencing. In response, an app named Kākano was co-designed with parents to help them better support their children’s mental health. Objective: The objective of this study was to evaluate the acceptability, feasibility, and effectiveness of Kākano, a mobile parenting app to increase parental confidence in supporting children struggling with their mental health. Methods: A cluster-randomized delayed access controlled trial was carried out in the Christchurch region between July 2019 and January 2020. Parents were recruited through schools and block randomized to receive immediate or delayed access to Kākano. Participants were given access to the Kākano app for 4 weeks and encouraged to use it weekly. Web-based pre- and postintervention measurements were undertaken. Results: A total of 231 participants enrolled in the Kākano trial, with 205 (88.7%) participants completing baseline measures and being randomized (101 in the intervention group and 104 in the delayed access control group). Of these, 41 (20%) provided full outcome data, of which 19 (18.2%) were for delayed access and 21 (20.8%) were for the immediate Kākano intervention. Among those retained in the trial, there was a significant difference in the mean change between groups favoring Kākano in the brief parenting assessment (F1,39=7, P=.012) but not in the Short Warwick-Edinburgh Mental Well-being Scale (F1,39=2.9, P=.099), parenting self-efficacy (F1,39=0.1, P=.805), family cohesion (F1,39=0.4, P=.538), or parenting sense of confidence (F1,40=0.6, P=.457). Waitlisted participants who completed the app after the waitlist period showed similar trends for the outcome measures with significant changes in the brief assessment of parenting and the Short Warwick-Edinburgh Mental Well-being Scale. No relationship between the level of app usage and outcome was found. Although the app was designed with parents, the low rate of completion of the trial was disappointing. Conclusions: Kākano is an app co-designed with parents to help manage their children’s mental health. There was a high rate of attrition, as is often seen in digital health interventions. However, for those who did complete the intervention, there was some indication of improved parental well-being and self-assessed parenting. Preliminary indications from this trial show that Kākano has promising acceptability, feasibility, and effectiveness, but further investigation is warranted. Trial Registration: Australia New Zealand Clinical Trials Registry ACTRN12619001040156; https://www.anzctr.org.au/Trial/Registration/TrialReview.aspx?id=377824&isReview=true
Diverse Density proposes an alternative housing strategy to the idealistic top-down process of housing development. The term ‘Top – down’ refers to a situation in which decisions are made by a few people in authority rather than by the people who are affected by the decisions (Cambridge). Problems/Position/Question: New Zealand’s urban housing is in a period of flux. Pressures of densification have permitted the intervention of medium density housing development schemes but these are not always successful. These typically top-down processes often result in internally focused design schemes that do not adhere to their specific context. The subsequent design outcomes can cause detrimental impacts to the local, urban and architectural conditions. With vast quantities of council regulations, building restrictions and design guidelines clouding over the housing sector, commonly referred to as ‘red tape’, occupant participation in the housing development sector is dwindling. A boundless separation between top-down and traditional housing processes has occurred and our existing neighbourhoods and historic architectural character are taking on the brunt of the problem. The thought-provoking, alternative housings strategies of key research theorists Alejandro Aravena and John Habraken frame positions that challenge contemporary densification methods with an alternative strategy. This position is addressed by endeavoring to answer; How can demands for denser housing achieve dynamic design responses that adhere to changes in occupancy, function and local site conditions? Aim: The aim of this thesis is to challenge New Zealand’s current housing densification methods by proposing an alternative densification strategy. Explicit devotion will be attributed to opposing top-down building developments. Secondly, this thesis aims to test a speculative site-specific housing model. The implementation of a Christchurch housing scenario will situate an investigative study to test the strategy and its ability to stimulate greater diversity, site responsiveness, functional adaptability and occupancy permutation. The post-earthquake housing conditions of Christchurch provide an appropriate scenario to test and implement design-led investigations. Objectives: The primary objectives of this design-led research investigation it to challenge the idealistic top-down method of developing density with a new method to: - Develop contextual architectural cohesion - Encourage residential diversity - Reinvigorate architectural autonomy - Respond to, and recognise, existing site conditions - Develop a housing model that: - Adapts to occupant functionality preferences - Caters to occupancy diversity - Achieves contextual responsiveness The proposition is addressed through a speculative design-led scenario study. A well-established Christchurch urban environment is adopted to implement and critique the envisioned alternative strategy. Development of the designs responsiveness, adaptability, and functionality produce a prototype housing model that actively adheres to its particular context. Implication: The implications of this research would be an alternative densification strategy to perceive the advancement of punctual assessment of building compliance. With accelerated building processes, the research may have implications for addressing New Zealand’s housing crisis whilst simultaneously providing diverse, personable and responsive architectural solutions. A more dynamic, up-to-date and responsive housing development sector would be informed.
To this extent, modern buildings generally demonstrated good resistance to collapse during the recent earthquakes in New Zealand. However, damage to non-structural elements (NSE) has been persistent during these events. NSEs include secondary systems or components attached to the floors, roofs, and walls of a building or industrial facility that are not explicitly designed to participate in the main vertical or lateral load-bearing mechanism of the structure. They play a major role in the operational and functional aspects of buildings and contribute a major portion of the building’s overall cost. Therefore, they are expected to accommodate the effects of seismic actions such as drifts and accelerations. Typical examples of NSEs include internal non-loadbearing partitions, suspended ceilings, sprinkler piping systems, architectural claddings, building contents, mechanical/electrical equipment, and furnishings. The main focus of this thesis is the drift sensitive NSEs: precast concrete cladding panels and internal partition walls. Even though most precast concrete cladding panels performed well from a life-safety point of view during recent earthquakes in NZ, some collapsed panels posed a significant threat to life safety. It is, therefore, important that the design and detailing of the panel-to-structure connections ensure that their strength and displacement capacity are adequate to meet the corresponding seismic demands, at least during design level earthquakes. In contrast, the partition wall is likely to get damaged and lose serviceability at a low inter-story drift unless designed to accommodate the relative deformations between them and the structure. Partition walls suffered wide-ranging damage such as screw failures, diagonal cracking, detachments to the gypsum linings, and anchorage failures during the 2011 Canterbury Earthquake Sequence in NZ. Therefore, the thesis is divided into two parts. Part I of the thesis focuses on developing novel low-damage precast concrete cladding panel connections, i.e. “rocking” connection details comprising vertically slotted steel embeds and weld plates. The low-damage seismic performance of novel “rocking” connection details is verified through experimental tests comprising uni-directional, bi-directional, and multi-storey scaled quasi-static cyclic tests. Comparison with the seismic performance of traditional panel connections reported in the literature demonstrated the system’s significantly improved seismic resilience. Furthermore, the finite element models of panel connections and sealants are developed in ABAQUS. The force-drift responses of the “rocking” panel system modelled in SAP2000 is compared with the experimental results to evaluate their accuracy and validity. Part II of the thesis focuses on a) understanding the seismic performance of traditional rigid timber-framed partition wall, b) development and verification of low-damage connections (i.e. “rocking” connection details comprising of dual-slot tracks), and c) seismic evaluation of partition walls with a novel “bracketed and slotted” connections (comprising of innovative fastener and plastic bracket named Flexibracket) under uni-directional and bidirectional quasi-static cyclic loadings. Moreover, parametric investigation of the partition walls was conducted through several experimental tests to understand better the pros and cons of the rocking connection details. The experimental results have confirmed that the implementation of the proposed low damage solutions of precast cladding panels and internal partition walls can significantly reduce their damage in a building.
Capacity design and hierarchy of strength philosophies at the base of modern seismic codes allow inelastic response in case of severe earthquakes and thus, in most traditional systems, damage develops at well-defined locations of reinforced concrete (RC) structures, known as plastic hinges. The 2010 and 2011 Christchurch earthquakes have demonstrated that this philosophy worked as expected. Plastic hinges formed in beams, in coupling beams and at the base of columns and walls. Structures were damaged permanently, but did not collapse. The 2010 and 2011 Christchurch earthquakes also highlighted a critical issue: the reparability of damaged buildings. No methodologies or techniques were available to estimate the level of subsequent earthquakes that RC buildings could still sustain before collapse. No repair techniques capable of restoring the initial condition of buildings were known. Finally, the cost-effectiveness of an eventual repair intervention, when compared with a new building, was unknown. These aspects, added to nuances of New Zealand building owners’ insurance coverage, encouraged the demolition of many buildings. Moreover, there was a perceived strong demand from government and industry to develop techniques for assessing damage to steel reinforcement bars embedded in cracked structural concrete elements. The most common questions were: “Have the steel bars been damaged in correspondence to the concrete cracks?”, “How much plastic deformation have the steel bars undergone?”, and “What is the residual strain capacity of the damaged bars?” Minimally invasive techniques capable of quantifying the level and extent of plastic deformation and residual strain capacity are not yet available. Although some studies had been recently conducted, a validated method is yet to be widely accepted. In this thesis, a least-invasive method for the damage-assessment of steel reinforcement is developed. Based on the information obtained from hardness testing and a single tensile test, it is possible to estimate the mechanical properties of earthquake-damaged rebars. The reduction in the low-cycle fatigue life due to strain ageing is also quantified. The proposed damage assessment methodology is based on empirical relationships between hardness and strain and residual strain capacity. If damage is suspected from in situ measurements, visual inspection or computer analysis, a bar may be removed and more accurate hardness measurements can be obtained using the lab-based Vickers hardness methodology. The Vickers hardness profile of damaged bars is then compared with calibration curves (Vickers hardness versus strain and residual strain capacity) previously developed for similar steel reinforcement bars extracted from undamaged locations. Experimental tests demonstrated that the time- and temperature-dependent strain-ageing phenomenon causes changes in the mechanical properties of plastically deformed steels. In particular, yield strength and hardness increases, whereas ductility decreases. The changes in mechanical properties are quantified and their implications on the hardness method are highlighted. Low-cycle fatigue (LCF) failures of steel reinforcing bars have been observed in laboratory testing and post-earthquake damage inspections. Often, failure might not occur during a first seismic event. However, damage is accumulated and the remaining fatigue life is reduced. Failure might therefore occur in a subsequent seismic event. Although numerous studies exist on the LCF behaviour of steel rebars, no studies had been conducted on the strain-ageing effects on the remaining fatigue life. In this thesis, the reduction in fatigue life due to this phenomenon is determined through a number of experimental tests.
The lateral capacity of a conventional CLT shear wall is often governed by the strength and stiffness of its connections, which do not significantly utilize the in-plane strength of the CLT. Therefore, CLT shear walls are not yet being used efficiently in the construction of mass timber buildings due to a lack of research on high-capacity connections and alternative wall configurations. In this study, cyclic experiments were completed on six full-scale, 5-ply cantilever CLT shear walls with high-capacity hold-downs using mixed angle screws and bolts. All specimens exhibited significantly higher strength and stiffness than previously tested conventional CLT shear walls in the literature. The base connections demonstrated ductile failure modes through yielding of the hold-down connections. Based on the experimental results, numerical models were calibrated to investigate the seismic behaviour of CLT shear walls for prototype buildings of 3 and 6-storeys in Christchurch, NZ. As an alternative to cantilever (single) shear walls, a type of coupled wall with steel link beams between adjacent CLT wall piers was investigated. Effective coupling requires the link beam-to-wall connections to have adequate strength to ensure ductile link beam responses and adequate stiffness to yield the link beams at a relatively low inter-storey drift level. To this end, three beam-to-wall connection types were developed and cyclically tested to investigate their behaviour and feasibility. Based on the test results of the critical connection, a 3-storey, 2/3-scale coupled CLT wall specimen with three steel link beams and mixed angle screwed hold-downs was cyclically tested to evaluate its performance and experimentally validate the system concept. The test results showed a relatively high lateral strength compared to conventional CLT shear walls, as well as a high system ductility ratio of 7.6. Failure of the system was characterised by combined bending and withdrawal of the screws in the mixed angle screw hold-downs, yielding and eventual inelastic buckling of the steel link beams, CLT toe crushing, and local CLT delamination. Following the initial test, the steel link beams, mixed angle screw hold-downs, and damaged CLT regions were repaired, then the wall specimen was re-tested. The repaired wall behaved similarly to the original test and exhibited slightly higher energy dissipation and peak strength, but marginally more rapid strength deterioration under cyclic loading. Several hybrid coupled CLT shear walls were numerically modelled and calibrated based on the results of the coupled wall experiments. Pushover analyses were conducted on a series of configurations to validate a capacity design method for the system and to investigate reasonable parameter values for use in the preliminary design of the system. Additionally, an iterative seismic design method was proposed and used to design sample buildings of 6, 8, and 10-storeys using both nonlinear pushover and nonlinear time history analyses to verify the prototype designs. Results of the sample building analyses demonstrated adequate seismic behaviour and the proposed design parameters were found to be appropriate. In summary, high-capacity CLT shear walls can be used for the resistance of earthquakes by using stronger base connections and coupled wall configurations. The large-scale experimental testing in this study has demonstrated that both cantilever and coupled CLT shear walls are feasible LLRSs which can provide significantly greater lateral strength, stiffness, and energy dissipation than conventional CLT shear wall configurations.
The assessment of damage and remaining capacity after an earthquake is an immediate measure to determine whether a reinforced concrete (RC) building is usable and safe for occupants. The recent Christchurch earthquake (22 February 2011) caused a uniquely severe level of structural damage to modern buildings, resulting in extensive damage to the building stock. About 60% of damaged multistorey concrete buildings (3 storeys and up) were demolished after the earthquake, and the cost of reconstruction amounted to 40 billion NZD. The aftermath disclosed issues of great complexities regarding the future of the RC buildings damaged by the earthquakes. This highlighted the importance of post-event decision-making, as the outcome will allow the appropriate course of action—demolition, repair or acceptance of the existing building—to be considered. To adopt the proper strategy, accurate assessment of the residual capacity and the level of damage is required. This doctoral dissertation aims to assess the damage and remaining capacity at constituent material and member level (i.e., concrete material and beams) through a systematic approach in an attempt to address part of an existing gap in the available literature. Since the residual capacity of RC members is not unique and depends on previously applied loading history, post-event residual capacity in this study was assessed in terms of fraction of fatigue life (i.e., the number of cycles required to failure). This research comprises three main parts: (1) residual capacity and damage assessment at material level (i.e., concrete), (2) post-yield bond deterioration and damage assessment at the interface of steel and concrete, and, finally, (3) residual capacity and damage assessment at member level (i.e., RC beam). The first part of this research focused on damage assessment and the remaining capacity of concrete from a material point of view. It aimed to employ appropriate and reliable durability-based testing and image-detection techniques to quantify deterioration in the mechanical properties of concrete on the basis that stress-induced damage occurred in the microstructural system of the concrete material. To this end, in the first phase, a feasibility study was conducted in which a combination of oxygen permeability, electrical resistivity and porosity tests were assessed to determine if they were robust and reliable enough to reveal damage which occurred in the microstructural system of concrete. The results, in terms of change in permeability, electrical resistivity and porosity features of disk samples taken from the middle third of damaged concrete cylinders (200 mm × 100 mm) monotonically pre-loaded to 50%, 70%, 90% and 95% of the ultimate strength (f′c), showed the permeability test is a reliable tool to identify the degree of damage, due to its high sensitivity to the load-induced microcracking. In parallel, to determine the residual capacity, the companion damaged concrete cylinders already loaded to the same level of compressive strength were reloaded up to failure. Comparing the stress–strain relationship of damaged concrete with intact material, it was also found that the strain capacity of the reloaded pre-damaged concrete cylinders decreases while strength remained virtually unchanged. In the second phase of the first part, a fluorescent microscopy technique was used to assess the damage and develop a correlation between material degradation, by virtue of the geometrical features, and damage to the concrete. To account for the effect of confinement and cyclic loading, in the third phase, the residual capacity and damage assessment of unconfined and GFRP confined concrete cylinders subjected to low-cycle fatigue loading, was investigated. Similar to the first phase, permeability testing technique was used to provide an indirect evaluation of fatigue damage. Finally, in the fourth phase of the first part, the suitability of permeability testing technique to assess damage was evaluated for cored concrete taken from three types of RC members: columns, beams and a beam-column joint. In view of the fact that the composite action of an RC member is highly dependent on the bond between reinforcement and surrounding concrete, understanding the deterioration of the bond in the post-yield range of strain in steel was crucial to assess damage at member level. Therefore, in the second phase of this research, a state-of-the- art distributed fibre optic strain sensor system (DFOSSS) system was used to evaluate bond deterioration in a cantilever RC beam subjected to monotonic lateral loading. The technology allowed the continuous capture of strain, every 2.6 mm along the length, in both reinforcing bars and cover concrete. The strain profile provided a basis by which the slip, axial stress and bond stress distributions were then established. In the third part, the study focused on the damage assessment and residual capacity of seven half-scale RC beams subjected to a constant-amplitude cyclic loading protocol. In the first stage, the structural performances of three specimens under constant-amplitude fatigue at 1%, 2% and 4% chord rotation (drift) were examined. In parallel, the number of cycles to failure, degradation in strength, stiffness and energy dissipation were characterized. In the second stage, four RC beams were subjected to loading up to 70% and 90% of their fatigue life, at 2% and 4% drift, and then monotonically pulled up to failure. To determine the residual flexural capacity, the lateral force–displacement results of pre-damaged specimens were compared with an undamaged specimen subjected to only monotonic loading. The study showed significant losses in strength, deformability, stiffness and energy dissipation capacity. A nonlinear finite element analysis (FEA) using concrete damage plasticity (CDP) model was also conducted in ABAQUS to numerically investigate the behaviour of the tested specimen. The results of the FE simulations indicated a reasonable response compared with the behaviour of the test specimen in terms of force–displacement and cracking pattern. During the Christchurch earthquake it was observed that the loading history has a significant influence on structural responses. While in conventional pseudo-static loading protocol, internal forces can be redistributed along the plastic length: there is little chance for structures undergoing high initial loading amplitude to redistribute pertinent stresses. As a result, in the third phase of this part, the effect of high rate of loading on the behaviour of seismically designed RC beams was investigated. Two half-scale cantilever RC beams were subjected to similar constant-amplitude cyclic loading at 2% and 4% drifts, but at a rate of 500 mm/s. Due to the incapability of conventional measuring techniques, a motion-tracking system was employed for data acquisition with the high-speed tests. The effect of rate of loading on the fatigue life of specimens (i.e., the number of cycles required to failure), secant stiffness, failure mode, cracking pattern, beam elongations and bar fracture surface were analysed. Integrating the results of all parts of this research has resulted in a better understanding of residual capacity and the development of damage at both the material and member level by using a low-cycle fatigue approach.
Though generally considered “natural” disasters, cyclones and earthquakes are increasingly being associated with human activities, incubated through urban settlement patterns and the long-term redistribution of natural resources. As society is becoming more urbanized, the risk of human exposure to disasters is also rising. Architecture often reflects the state of society’s health: architectural damage is the first visible sign of emergency, and reconstruction is the final response in the process of recovery. An empirical assessment of architectural projects in post-disaster situations can lead to a deeper understanding of urban societies as they try to rebuild. This thesis offers an alternative perspective on urban disasters by looking at the actions and attitudes of disaster professionals through the lens of architecture, situated in recent events: the 2010 Christchurch earthquake, the 2010 Haiti earthquake, and the 2005 Hurricane Katrina. An empirical, multi-hazard, cross-sectional case study methodology was used, employing grounded theory method to build theory, and a critical constructivist strategy to inform the analysis. By taking an interdisciplinary approach to understanding disasters, this thesis positions architecture as a conduit between two divergent approaches to disaster research: the hazards approach, which studies the disaster cycles from a scientific perspective; and the sociological approach, which studies the socially constructed vulnerabilities that result from disasters, and the elements of social change that accompany such events. Few studies to date have attempted to integrate the multi-disciplinary perspectives that can advance our understanding of societal problems in urban disasters. To bridge this gap, this thesis develops what will be referred to as the “Rittelian framework”—based on the work of UC Berkeley’s architecture professor Horst Rittel (1930-1990). The Rittelian framework uses the language of design to transcend the multiple fields of human endeavor to address the “design problems” in disaster research. The processes by which societal problems are addressed following an urban disaster involve input by professionals from multiple fields—including economics, sociology, medicine, and engineering—but the contribution from architecture has been minimal to date. The main impetus for my doctoral thesis has been the assertion that most of the decisions related to reconstruction are made in the early emergency recovery stages where architects are not involved, but architects’ early contribution is vital to the long-term reconstruction of cities. This precipitated in the critical question: “How does the Rittelian framework contribute to the critical design decisions in modern urban disasters?” Comparative research was undertaken in three case studies of recent disasters in New Orleans (2005), Haiti (2010) and Christchurch (2010), by interviewing 51 individuals who were selected on the basis of employing the Rittelian framework in their humanitarian practice. Contextualizing natural disaster research within the robust methodological framework of architecture and the analytical processes of sociology is the basis for evaluating the research proposition that architectural problem solving is of value in addressing the ‘Wicked Problems’ of disasters. This thesis has found that (1) the nuances of the way disaster agents interpret the notion of “building back better” can influence the extent to which architectural professionals contribute in urban disaster recovery, (2) architectural design can be used to facilitate but also impede critical design decisions, and (3) framing disaster research in terms of design decisions can lead to innovation where least expected. This empirical research demonstrates how the Rittelian framework can inform a wider discussion about post-disaster human settlements, and improve our resilience through disaster research.
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.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
The increase of the world's population located near areas prone to natural disasters has given rise to new ‘mega risks’; the rebuild after disasters will test the governments’ capabilities to provide appropriate responses to protect the people and businesses. During the aftermath of the Christchurch earthquakes (2010-2012) that destroyed much of the inner city, the government of New Zealand set up a new partnership between the public and private sector to rebuild the city’s infrastructure. The new alliance, called SCIRT, used traditional risk management methods in the many construction projects. And, in hindsight, this was seen as one of the causes for some of the unanticipated problems. This study investigated the risk management practices in the post-disaster recovery to produce a specific risk management model that can be used effectively during future post-disaster situations. The aim was to develop a risk management guideline for more integrated risk management and fill the gap that arises when the traditional risk management framework is used in post-disaster situations. The study used the SCIRT alliance as a case study. The findings of the study are based on time and financial data from 100 rebuild projects, and from surveying and interviewing risk management professionals connected to the infrastructure recovery programme. The study focussed on post-disaster risk management in construction as a whole. It took into consideration the changes that happened to the people, the work and the environment due to the disaster. System thinking, and system dynamics techniques have been used due to the complexity of the recovery and to minimise the effect of unforeseen consequences. Based on an extensive literature review, the following methods were used to produce the model. The analytical hierarchical process and the relative importance index have been used to identify the critical risks inside the recovery project. System theory methods and quantitative graph theory have been used to investigate the dynamics of risks between the different management levels. Qualitative comparative analysis has been used to explore the critical success factors. And finally, causal loop diagrams combined with the grounded theory approach has been used to develop the model itself. The study identified that inexperienced staff, low management competency, poor communication, scope uncertainty, and non-alignment of the timing of strategic decisions with schedule demands, were the key risk factors in recovery projects. Among the critical risk groups, it was found that at a strategic management level, financial risks attracted the highest level of interest, as the client needs to secure funding. At both alliance-management and alliance-execution levels, the safety and environmental risks were given top priority due to a combination of high levels of emotional, reputational and media stresses. Risks arising from a lack of resources combined with the high volume of work and the concern that the cost could go out of control, alongside the aforementioned funding issues encouraged the client to create the recovery alliance model with large reputable construction organisations to lock in the recovery cost, at a time when the scope was still uncertain. This study found that building trust between all parties, clearer communication and a constant interactive flow of information, established a more working environment. Competent and clear allocation of risk management responsibilities, cultural shift, risk prioritisation, and staff training were crucial factors. Finally, the post-disaster risk management (PDRM) model can be described as an integrated risk management model that considers how the changes which happened to the environment, the people and their work, caused them to think differently to ease the complexity of the recovery projects. The model should be used as a guideline for recovery systems, especially after an earthquake, looking in detail at all the attributes and the concepts, which influence the risk management for more effective PDRM. The PDRM model is represented in Causal Loops Diagrams (CLD) in Figure 8.31 and based on 10 principles (Figure 8.32) and 26 concepts (Table 8.1) with its attributes.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
Reinforced concrete structures designed in pre-1970s are vulnerable under earthquakes due to lack of seismic detailing to provide adequate ductility. Typical deficiencies of pre-1970s reinforced concrete structures are (a) use of plain bars as longitudinal reinforcement, (b) inadequate anchorage of beam longitudinal reinforcement in the column (particularly exterior column), (c) lack of joint transverse reinforcement if any, (d) lapped splices located just above joint, and (e) low concrete strength. Furthermore, the use of infill walls is a controversial issue because it can help to provide additional stiffness to the structure on the positive side and on the negative side it can increase the possibility of soft-storey mechanisms if it is distributed irregularly. Experimental research to investigate the possible seismic behaviour of pre-1970s reinforced concrete structures have been carried out in the past. However, there is still an absence of experimental tests on the 3-D response of existing beam-column joints under bi-directional cyclic loading, such as corner joints. As part of the research work herein presented, a series of experimental tests on beam-column subassemblies with typical detailing of pre-1970s buildings has been carried out to investigate the behaviour of existing reinforced concrete structures. Six two-third scale plane frame exterior beam-column joint subassemblies were constructed and tested under quasi-static cyclic loading in the Structural Laboratory of the University of Canterbury. The reinforcement detailing and beam dimension were varied to investigate their effect on the seismic behaviour. Four specimens were conventional deep beam-column joint, with two of them using deformed longitudinal bars and beam bars bent in to the joint and the two others using plain round longitudinal bars and beam bars with end hooks. The other two specimens were shallow beam-column joint, one with deformed longitudinal bars and beam bars bent in to the joint, the other with plain round longitudinal bars and beam bars with end hooks. All units had one transverse reinforcement in the joint. The results of the experimental tests indicated that conventional exterior beam-column joint with typical detailing of pre-1970s building would experience serious diagonal tension cracking in the joint panel under earthquake. The use of plain round bars with end hooks for beam longitudinal reinforcement results in more severe damage in the joint core when compared to the use of deformed bars for beam longitudinal reinforcement bent in to the joint, due to the combination of bar slips and concrete crushing. One interesting outcome is that the use of shallow beam in the exterior beam-column joint could avoid the joint cracking due to the beam size although the strength provided lower when compared with the use of deep beam with equal moment capacity. Therefore, taking into account the low strength and stiffness, shallow beam can be reintroduced as an alternative solution in design process. In addition, the presence of single transverse reinforcement in the joint core can provide additional confinement after the first crack occurred, thus delaying the strength degradation of the structure. Three two-third scale space frame corner beam-column joint subassemblies were also constructed to investigate the biaxial loading effect. Two specimens were deep-deep beam-corner column joint specimens and the other one was deep-shallow beam-corner column joint specimen. One deep-deep beam-corner column joint specimen was not using any transverse reinforcement in the joint core while the two other specimens were using one transverse reinforcement in the joint core. Plain round longitudinal bars were used for all units with hook anchorage for the beam bars. Results from the tests confirmed the evidences from earthquake damage observations with the exterior 3-D (corner) beam-column joint subjected to biaxial loading would have less strength and suffer higher damage in the joint area under earthquake. Furthermore, the joint shear relation in the two directions is calibrated from the results to provide better analysis. An analytical model was used to simulate the seismic behaviour of the joints with the help of Ruaumoko software. Alternative strength degradation curves corresponding to different reinforcement detailing of beam-column joint unit were proposed based on the test results.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
Coastal margins are exposed to rising sea levels that present challenging circumstances for natural resource management. This study investigates a rare example of tectonic displacement caused by earthquakes that generated rapid sea-level change in a tidal lagoon system typical of many worldwide. This thesis begins by evaluating the coastal squeeze effects caused by interactions between relative sea-level (RSL) rise and the built environment of Christchurch, New Zealand, and also examples of release from similar effects in areas of uplift where land reclamations were already present. Quantification of area gains and losses demonstrated the importance of natural lagoon expansion into areas of suitable elevation under conditions of RSL rise and showed that they may be necessary to offset coastal squeeze losses experienced elsewhere. Implications of these spatial effects include the need to provide accommodation space for natural ecosystems under RSL rise, yet other land-uses are likely to be present in the areas required. Consequently, the resilience of these environments depends on facilitating transitions between human land-uses either proactively or in response to disaster events. Principles illustrated by co-seismic sea-level change are generally applicable to climate change adaptation due to the similarity of inundation effects. Furthermore, they highlight the potential role of non-climatic factors in determining the overall trajectory of change. Chapter 2 quantifies impacts on riparian wetland ecosystems over an eight year period post- quake. Coastal wetlands were overwhelmed by RSL rise and recovery trajectories were surprisingly slow. Four risk factors were identified from the observed changes: 1) the encroachment of anthropogenic land-uses, 2) connectivity losses between areas of suitable elevation, 3) the disproportionate effect of larger wetland vulnerabilities, and 4) the need to protect new areas to address the future movement of ecosystems. Chapter 3 evaluates the unique context of shoreline management on a barrier sandspit under sea-level rise. A linked scenario approach was used to evaluate changes on the open coast and estuarine shorelines simultaneously and consider combined effects. The results show dune loss from a third of the study area using a sea-level rise scenario of 1 m over 100 years and with continuation of current land-uses. Increased exposure to natural hazards and accompanying demand for seawalls is a likely consequence unless natural alternatives can be progressed. In contrast, an example of managed retreat following earthquake-induced subsidence of the backshore presents a new opportunity to restart saltmarsh accretion processes seaward of coastal defences with the potential to reverse decades of degradation and build sea-level rise resilience. Considering both shorelines simultaneously highlights the existence of pinch-points from opposing forces that result in small land volumes above the tidal range. Societal adaptation is delicately poised between the paradigms of resisting or accommodating nature and challenged by the long perimeter and confined nature of the sandspit feature. The remaining chapters address the potential for salinity effects caused by tidal prism changes with a focus on the conservation of īnanga (Galaxias maculatus), a culturally important fish that supports New Zealand‘s whitebait fishery. Methodologies were developed to test the hypothesis that RSL changes would drive a shift in the distribution of spawning sites with implications for their management. Chapter 4 describes a new practical methodology for quantifying the total productivity and spatiotemporal variability of spawning sites at catchment scale. Chapter 5 describes the novel use of artificial habitats as a detection tools to help overcome field survey limitations in degraded environments where egg mortality can be high. The results showed that RSL changes resulted in major shifts in spawning locations and these were associated with new patterns of vulnerability due to the continuation of pre-disturbance land-uses. Unexpected findings includes an improved understanding of the spatial relationship between salinity and spawning habitat, and identification of an invasive plant species as important spawning habitat, both with practical management implications. To conclude, the design of legal protection mechanisms was evaluated in relation to the observed habitat shifts and with a focus on two new planning initiatives that identified relatively large protected areas (PAs) in the lower river corridors. Although the larger PAs were better able to accommodate the observed habitat shifts inefficiencies were also apparent due to spatial disparities between PA boundaries and the values requiring protection. To reduce unnecessary trade-offs with other land-uses, PAs of sufficient size to cover the observable spatiotemporal variability and coupled with adaptive capacity to address future change may offer a high effectiveness from a network of smaller PAs. The latter may be informed by both monitoring and modelling of future shifts and these are expected to include upstream habitat migration driven by the identified salinity relationships and eustatic sea-level rise. The thesis concludes with a summary of the knowledge gained from this research that can assist the development of a new paradigm of environmental sustainability incorporating conservation and climate change adaptation. Several promising directions for future research identified within this project are also discussed.
One of the less understood geotechnical responses to the cyclic loading from the MW6.2 Christchurch Earthquake, on the 22nd of February 2011, is the fissuring in the loessial soil-mantled, footslope positions of the north-facing valleys of the Port Hills. The fissures are characterized by mostly horizontal offset (≤500mm), with minor vertical displacement (≤300mm), and they extend along both sides of valleys for several hundred metres in an approximately contour-parallel orientation. The fissure traces correspond to extensional features mapped in other studies. Previous studies have suggested that the fissures are the headscarps of incipient landslides, but the surface and subsurface features are not typical of landslide movement. Whilst there are some features that correlate with landslide movement, there are many features that contradict the landslide movement hypothesis. Of critical importance to this investigation was the fact that there are no landslide flanks, there has been no basal shear surface found, there is little deformation in the so-called ‘landslide body’, and there have been no recorded zones of low shear strength in the soil deposit that are indicative of a basal shear surface. This thesis is a detailed geotechnical study on the fissures along part of Ramahana Road in the Hillsborough Valley, Christchurch. Shallow and deep investigation methods found that the predominant soil is loess-colluvium, to depths of ~20m, and this soil has variable geotechnical characteristics depending on the layer sampled. The factor that has the most influence on shear strength was found to be the moisture content. Direct shear-box testing of disturbed, recompacted loess-colluvium found that the soil had a cohesion of 35-65kPa and a friction angle of 38-43° when the soil moisture content was at 8-10%. However when the moisture content was at 19-20% the soil’s cohesion decreased to 3-5kPa and its friction angle decreased to 33-38°, this moisture content is at or slightly above the plastic limit. An electrical resistivity geophysical survey was conducted perpendicular to multiple fissure traces and through the compressional zone at 17 Ramahana Road. The electrical resistivity line found that there was an area of high resistivity at the toe of the slope, and an area of high conductivity downslope of this and at greater depths. This area correlated to the compressional zone recorded by previous studies. Moisture content testing of the soil in these locations showed that the soil in the resistive area was relatively dry (9%) compared to the surrounding soil (13%), whilst the soil in the conductive area was relatively wet (22%)compared to the surrounding soil (19%). Density tests of the soil in the compressional zone recorded that the resistive area had a higher dry density than the surrounding soil (~1790 kg/m3 compared to ~1650 kg/m3). New springs arose downslope of the compressional zone contemporaneously with the fissures, and it is interpreted that these have arisen from increased hydraulic head in the Banks Peninsula bedrock aquifer system, and earthquake induced-bedrock fracturing. A test pit was dug across an infilled fissure trace at 17 Ramahana Road to a depth of 3m. The fissure trace had an aperture of 450-470mm at the ground surface, but it gradually lost aperture with depth until 2.0-2.1m where it became a segmented fissure trace with 1-2mm aperture. A mixed-colluvium layer was intercepted by the fissure trace at 2.4m depth, and there was no observable vertical offset of this layer. The fissure trace was at an angle of 78° at the ground surface, but it also flattened with depth, which gave it a slightly curved appearance. The fissure trace was at an assumed angle of 40-50° near the base of the test pit. Rotational slide, translational slide and lateral spread landslide movement types were compared and contrasted as possibilities for landslide movement types, whilst an alternative hypothesis was offered that the fissures are tensile failures with a quasi-toppling motion involving a cohesive block of loessial soil moving outwards from the slope, with an accommodating compressional strain in the lower less cohesive soil. The mechanisms behind this movement are suggested to be the horizontal earthquake inertia forces from the Christchurch Earthquake, the static shear stress of the slope, and bedrock uplift of the Port Hills in relation to the subsidence of the Christchurch city flatlands. Extremely high PGA is considered to be a prerequisite to the fissure trace development, and these can only be induced in the Hillsborough Valley from a Port Hills Fault rupture, which has a recurrence interval of ~10,000 years. The current understanding of how the loess-colluvium soil would behave under cyclic loading is limited, and the mechanisms behind the suggested movement type are not completely understood. Further research is needed to confirm the proposed mechanism of the fissure traces. Laboratory tests such as the cyclic triaxial and cyclic shear test would be beneficial in future research to quantitatively test how the soil behaves under cyclic loading at various moisture contents and clay contents, and centrifuge experiments would be of great use to qualitatively test the suggested mode of movement in the loessial soil.
Background The 2010/2011 Canterbury earthquakes and aftershocks in New Zealand caused unprecedented destruction to the physical, social, economic, and community fabric of Christchurch city. The recovery phase in Christchurch is on going, six years following the initial earthquake. Research exploring how disabled populations experience community inclusion in the longer-term recovery following natural disasters is scant. Yet such information is vital to ensure that recovering communities are inclusive for all members of the affected population. This thesis specifically examined how people who use wheelchairs experienced community inclusion four years following the 2010/2011 Canterbury earthquakes. Aims The primary research aim was to understand how one section of the disability community – people who use wheelchairs – experienced community inclusion over the four years following the 2010/2011 Canterbury earthquakes and aftershocks. A secondary aim was to test a novel sampling approach, Respondent Driven Sampling, which had the potential to enable unbiased population-based estimates. This was motivated by the lack of an available sampling frame for the target population, which would inhibit recruitment of a representative sample. Methodology and methods An exploratory sequential mixed methods design was used, beginning with a qualitative phase (Phase One), which informed a second quantitative phase (Phase Two). The qualitative phase had two stages. First, a small sample of people who use wheelchairs participated in an individual, semi-structured interview. In the second stage, these participants were then invited to a group interview to clarify and prioritise themes identified in the individual interviews. The quantitative phase was a cross-sectional survey developed from the findings from Phase One. Initially, Respondent Driven Sampling was employed to conduct a national, electronic cross-sectional survey that aimed to recruit a sample that may provide unbiased population-based estimates. Following the unsuccessful application of Respondent Driven Sampling, a region-specific convenience sampling approach was used. The datasets from the qualitative and quantitative phases were integrated to address the primary aim of the research. Results In Phase One 13 participants completed the individual interviews, and five of them contributed to the group interview. Thematic analysis of individual and group interview data suggested that participants felt the 2010/11 earthquakes magnified many pre-existing barriers to community inclusion, and also created an exciting opportunity for change. This finding was encapsulated in five themes: 1) earthquakes magnified barriers, 2) community inclusion requires energy, 3) social connections are important, 4) an opportunity lost, and 5) an opportunity found. The findings from Phase One informed the development of a survey instrument to investigate how these findings generalised to a larger sample of individuals who use wheelchairs. In Phase Two, the Respondent Driven Sampling approach failed to recruit enough participants to satisfy the statistical requirements needed to reach equilibrium, thereby enabling the calculation of unbiased population estimates. The subsequent convenience sampling approach recruited 49 participants who, combined with the 15 participants from the Respondent Driven Sampling approach that remained eligible for the region-specific sample, resulted in the total of 64 individuals who used wheelchairs and were residents of Christchurch. Participants reported their level of community inclusion at three time periods: the six months prior to the first earthquake in September 2010 (time one), the six months following the first earthquake in September 2010 (time two), and the six months prior to survey completion (between October 2015 and March 2016, (time three)). Survey data provided some precision regarding the timing in which the magnified barriers developed. Difficulty with community inclusion rose significantly between time one and time two, and while reducing slightly, was still present during time three, and had not returned to the time one baseline. The integrated findings from Phase One and Phase Two suggested that magnified barriers to community inclusion had been sustained four years post-earthquake, and community access had not returned to pre-earthquake levels, let alone improved beyond pre-earthquake levels. Conclusion Findings from this mixed methods study suggest that four years following the initial earthquake, participants were still experiencing multiple magnified barriers, which contributed to physical and social exclusion, as well as fatigue, as participants relied on individual agency to negotiate such barriers. Participants also highlighted the exciting opportunity to create an accessible city. However because they were still experiencing barriers four years following the initial event, and were concerned that this opportunity might be lost if the recovery proceeds without commitment and awareness from the numerous stakeholders involved in guiding the recovery. To truly realise the opportunity to create an accessible city following a disaster, the transition from the response phase to a sustainable longer-term recovery must adopt a new model of community engagement where decision-makers partner with people living with disability to co-produce a vision and strategy for creating an inclusive community. Furthermore, despite the unsuccessful use of Respondent Driven Sampling in this study, future research exploring the application of RDS with wheelchair users is recommended before discounting this sampling approach in this population.
One of the great challenges facing human systems today is how to prepare for, manage, and adapt successfully to the profound and rapid changes wreaked by disasters. Wellington, New Zealand, is a capital city at significant risk of devastating earthquake and tsunami, potentially requiring mass evacuations with little or short notice. Subsequent hardship and suffering due to widespread property damage and infrastructure failure could cause large areas of the Wellington Region to become uninhabitable for weeks to months. Previous research has shown that positive health and well-being are associated with disaster-resilient outcomes. Preventing adverse outcomes before disaster strikes, through developing strengths-based skill sets in health-protective attitudes and behaviours, is increasingly advocated in disaster research, practise, and management. This study hypothesised that well-being constructs involving an affective heuristic play vital roles in pathways to resilience as proximal determinants of health-protective behaviours. Specifically, this study examined the importance of health-related quality of life and subjective well-being in motivating evacuation preparedness, measured in a community sample (n=695) drawn from the general adult population of Wellington’s isolated eastern suburbs. Using a quantitative epidemiological approach, the study measured the prevalence of key quality of life indicators (physical and mental health, emotional well-being or “Sense of Coherence”, spiritual well-being, social well-being, and life satisfaction) using validated psychometric scales; analysed the strengths of association between these indicators and the level of evacuation preparedness at categorical and continuous levels of measurement; and tested the predictive power of the model to explain the variance in evacuation preparedness activity. This is the first study known to examine multi-dimensional positive health and global well-being as resilient processes for engaging in evacuation preparedness behaviour. A cross-sectional study design and quantitative survey were used to collect self-report data on the study variables; a postal questionnaire was fielded between November 2008 and March 2009 to a sampling frame developed through multi-stage cluster randomisation. The survey response rate was 28.5%, yielding a margin of error of +/- 3.8% with 95% confidence and 80% statistical power to detect a true correlation coefficient of 0.11 or greater. In addition to the primary study variables, data were collected on demographic and ancillary variables relating to contextual factors in the physical environment (risk perception of physical and personal vulnerability to disaster) and the social environment (through the construct of self-determination), and other measures of disaster preparedness. These data are reserved for future analyses. Results of correlational and regression analyses for the primary study variables show that Wellingtonians are highly individualistic in how their well-being influences their preparedness, and a majority are taking inadequate action to build their resilience to future disaster from earthquake- or tsunami-triggered evacuation. At a population level, the conceptual multi-dimensional model of health-related quality of life and global well-being tested in this study shows a positive association with evacuation preparedness at statistically significant levels. However, it must be emphasised that the strength of this relationship is weak, accounting for only 5-7% of the variability in evacuation preparedness. No single dimension of health-related quality of life or well-being stands out as a strong predictor of preparedness. The strongest associations for preparedness are in a positive direction for spiritual well-being, emotional well-being, and life satisfaction; all involve a sense of existential meaningfulness. Spiritual well-being is the only quality of life variable making a statistically significant unique contribution to explaining the variance observed in the regression models. Physical health status is weakly associated with preparedness in a negative direction at a continuous level of measurement. No association was found at statistically significant levels for mental health status and social well-being. These findings indicate that engaging in evacuation preparedness is a very complex, holistic, yet individualised decision-making process, and likely involves highly subjective considerations for what is personally relevant. Gender is not a factor. Those 18-24 years of age are least likely to prepare and evacuation preparedness increases with age. Multidimensional health and global well-being are important constructs to consider in disaster resilience for both pre-event and post-event timeframes. This work indicates a need for promoting self-management of risk and building resilience by incorporating a sense of personal meaning and importance into preparedness actions, and for future research into further understanding preparedness motivations.
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.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
In the last century, seismic design has undergone significant advancements. Starting from the initial concept of designing structures to perform elastically during an earthquake, the modern seismic design philosophy allows structures to respond to ground excitations in an inelastic manner, thereby allowing damage in earthquakes that are significantly less intense than the largest possible ground motion at the site of the structure. Current performance-based multi-objective seismic design methods aim to ensure life-safety in large and rare earthquakes, and to limit structural damage in frequent and moderate earthquakes. As a result, not many recently built buildings have collapsed and very few people have been killed in 21st century buildings even in large earthquakes. Nevertheless, the financial losses to the community arising from damage and downtime in these earthquakes have been unacceptably high (for example; reported to be in excess of 40 billion dollars in the recent Canterbury earthquakes). In the aftermath of the huge financial losses incurred in recent earthquakes, public has unabashedly shown their dissatisfaction over the seismic performance of the built infrastructure. As the current capacity design based seismic design approach relies on inelastic response (i.e. ductility) in pre-identified plastic hinges, it encourages structures to damage (and inadvertently to incur loss in the form of repair and downtime). It has now been widely accepted that while designing ductile structural systems according to the modern seismic design concept can largely ensure life-safety during earthquakes, this also causes buildings to undergo substantial damage (and significant financial loss) in moderate earthquakes. In a quest to match the seismic design objectives with public expectations, researchers are exploring how financial loss can be brought into the decision making process of seismic design. This has facilitated conceptual development of loss optimisation seismic design (LOSD), which involves estimating likely financial losses in design level earthquakes and comparing against acceptable levels of loss to make design decisions (Dhakal 2010a). Adoption of loss based approach in seismic design standards will be a big paradigm shift in earthquake engineering, but it is still a long term dream as the quantification of the interrelationships between earthquake intensity, engineering demand parameters, damage measures, and different forms of losses for different types of buildings (and more importantly the simplification of the interrelationship into design friendly forms) will require a long time. Dissecting the cost of modern buildings suggests that the structural components constitute only a minor portion of the total building cost (Taghavi and Miranda 2003). Moreover, recent research on seismic loss assessment has shown that the damage to non-structural elements and building contents contribute dominantly to the total building loss (Bradley et. al. 2009). In an earthquake, buildings can incur losses of three different forms (damage, downtime, and death/injury commonly referred as 3Ds); but all three forms of seismic loss can be expressed in terms of dollars. It is also obvious that the latter two loss forms (i.e. downtime and death/injury) are related to the extent of damage; which, in a building, will not just be constrained to the load bearing (i.e. structural) elements. As observed in recent earthquakes, even the secondary building components (such as ceilings, partitions, facades, windows parapets, chimneys, canopies) and contents can undergo substantial damage, which can lead to all three forms of loss (Dhakal 2010b). Hence, if financial losses are to be minimised during earthquakes, not only the structural systems, but also the non-structural elements (such as partitions, ceilings, glazing, windows etc.) should be designed for earthquake resistance, and valuable contents should be protected against damage during earthquakes. Several innovative building technologies have been (and are being) developed to reduce building damage during earthquakes (Buchanan et. al. 2011). Most of these developments are aimed at reducing damage to the buildings’ structural systems without due attention to their effects on non-structural systems and building contents. For example, the PRESSS system or Damage Avoidance Design concept aims to enable a building’s structural system to meet the required displacement demand by rocking without the structural elements having to deform inelastically; thereby avoiding damage to these elements. However, as this concept does not necessarily reduce the interstory drift or floor acceleration demands, the damage to non-structural elements and contents can still be high. Similarly, the concept of externally bracing/damping building frames reduces the drift demand (and consequently reduces the structural damage and drift sensitive non-structural damage). Nevertheless, the acceleration sensitive non-structural elements and contents will still be very vulnerable to damage as the floor accelerations are not reduced (arguably increased). Therefore, these concepts may not be able to substantially reduce the total financial losses in all types of buildings. Among the emerging building technologies, base isolation looks very promising as it seems to reduce both inter-storey drifts and floor accelerations, thereby reducing the damage to the structural/non-structural components of a building and its contents. Undoubtedly, a base isolated building will incur substantially reduced loss of all three forms (dollars, downtime, death/injury), even during severe earthquakes. However, base isolating a building or applying any other beneficial technology may incur additional initial costs. In order to provide incentives for builders/owners to adopt these loss-minimising technologies, real-estate and insurance industries will have to acknowledge the reduced risk posed by (and enhanced resilience of) such buildings in setting their rental/sale prices and insurance premiums.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
