Six months after the 4 September 2010 Mw 7.1 Darfield (Canterbury) earthquake, a Mw 6.2 Christchurch (Lyttelton) aftershock struck Christchurch on the 22 February 2011. This earthquake was centred approximately 10km south-east of the Christchurch CBD at a shallow depth of 5km, resulting in intense seismic shaking within the Christchurch central business district (CBD). Unlike the 4 Sept earthquake when limited-to-moderate damage was observed in engineered reinforced concrete (RC) buildings [35], in the 22 February event a high number of RC Buildings in the Christchurch CBD (16.2 % out of 833) were severely damaged. There were 182 fatalities, 135 of which were the unfortunate consequences of the complete collapse of two mid-rise RC buildings. This paper describes immediate observations of damage to RC buildings in the 22 February 2011 Christchurch earthquake. Some preliminary lessons are highlighted and discussed in light of the observed performance of the RC building stock. Damage statistics and typical damage patterns are presented for various configurations and lateral resisting systems. Data was collated predominantly from first-hand post-earthquake reconnaissance observations by the authors, complemented with detailed assessment of the structural drawings of critical buildings and the observed behaviour. Overall, the 22 February 2011 Mw 6.2 Christchurch earthquake was a particularly severe test for both modern seismically-designed and existing non-ductile RC buildings. The sequence of earthquakes since the 4 Sept 2010, particularly the 22 Feb event has confirmed old lessons and brought to life new critical ones, highlighting some urgent action required to remedy structural deficiencies in both existing and “modern” buildings. Given the major social and economic impact of the earthquakes to a country with strong seismic engineering tradition, no doubt some aspects of the seismic design will be improved based on the lessons from Christchurch. The bar needs to and can be raised, starting with a strong endorsement of new damage-resisting, whilst cost-efficient, technologies as well as the strict enforcement, including financial incentives, of active policies for the seismic retrofit of existing buildings at a national scale.
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.
Natural disasters are increasingly disruptive events that affect livelihoods, organisations, and economies worldwide. Research has identified the impacts and responses of organisations to different types of natural disasters, and have outlined factors, such as industry sector, that are important to organisational vulnerability and resilience. One of the most costly types of natural disasters in recent years has been earthquakes, and yet to date, the majority of studies have focussed on the effects of earthquakes in urban areas, while rural organisational impact studies have primarily focused on the effects of meteorological and climatic driven hazards. As a result, the likely impacts of an earthquake on rural organisations in a developed context is unconstrained in the literature. In countries like New Zealand, which have major earthquakes and agricultural sectors that are significant contributors to the economy, it is important to know what impacts an earthquake event would have on the rural industries, and how these impacts compare to that of a more commonly analysed, high-frequency event. In September of 2010, rural organisations in Canterbury experienced the 4 September 2010 Mw 7.1 `Darfield' earthquake and the associated aftershocks, which came to be known as the Canterbury earth- quake sequence. The earthquake sequence caused intense ground shaking, creating widespread critical service outages, structural and non-structural damage to built infrastructure, as well as ground surface damage from ooding, liquefaction and surface rupture. Concurrently on September 18 2010, rural organisations in Southland experienced an unseasonably late snowstorm and cold weather snap that brought prolonged sub-zero temperatures, high winds and freezing rain, damaging structures in the City of Invercargill and causing widespread livestock losses and production decreases across the region. This thesis documents the effects of the Canterbury earthquake sequence and Southland snowstorm on farming and rural non-farming organisations, utilizing comparable methodologies to analyse rural organisational impacts, responses and recovery strategies to natural disasters. From the results, a short- term impact assessment methodology is developed for multiple disasters. Additionally, a regional asset repair cost estimation model is proposed for farming organisations following a major earthquake event, and the use of social capital in rural organisational recovery strategies following natural disasters is analysed.
Depending on their nature and severity, disasters can create large volumes of debris and waste. Waste volumes from a single event can be the equivalent of many times the annual waste generation rate of the affected community. These volumes can overwhelm existing solid waste management facilities and personnel. Mismanagement of disaster waste can affect both the response and long term recovery of a disaster affected area. Previous research into disaster waste management has been either context specific or event specific, making it difficult to transfer lessons from one disaster event to another. The aim of this research is to develop a systems understanding of disaster waste management and in turn develop context- and disaster-transferrable decision-making guidance for emergency and waste managers. To research this complex and multi-disciplinary problem, a multi-hazard, multi-context, multi-case study approach was adopted. The research focussed on five major disaster events: 2011 Christchurch earthquake, 2009 Victorian Bushfires, 2009 Samoan tsunami, 2009 L’Aquila earthquake and 2005 Hurricane Katrina. The first stage of the analysis involved the development of a set of ‘disaster & disaster waste’ impact indicators. The indicators demonstrate a method by which disaster managers, planners and researchers can simplify the very large spectra of possible disaster impacts, into some key decision-drivers which will likely influence post-disaster management requirements. The second stage of the research was to develop a set of criteria to represent the desirable environmental, economic, social and recovery effects of a successful disaster waste management system. These criteria were used to assess the effectiveness of the disaster waste management approaches for the case studies. The third stage of the research was the cross-case analysis. Six main elements of disaster waste management systems were identified and analysed. These were: strategic management, funding mechanisms, operational management, environmental and human health risk management, and legislation and regulation. Within each of these system elements, key decision-making guidance (linked to the ‘disaster & disaster waste’ indicators) and management principles were developed. The ‘disaster & disaster waste’ impact indicators, the effects assessment criteria and management principles have all been developed so that they can be practically applied to disaster waste management planning and response in the future.
Recently developed performance-based earthquake engineering framework, such as one provided by PEER (Deierlein et al. 2003), assist in the quantification in terms of performance such as casualty, monetary losses and downtime. This opens up the opportunity to identify cost-effective retrofit/rehabilitation strategies by comparing upfront costs associated with retrofit with the repair costs that can be expected over time. This loss assessment can be strengthened by learning from recent earthquakes, such as the 2010 Canterbury and 2016 Kaikoura earthquakes. In order to investigate which types of retrofit/rehabilitation strategies may be most cost-effective, a case study building was chosen for this research. The Pacific Tower, a 22-storey EBF apartment located within the Christchurch central business district (CBD), was damaged and repaired during the 2010 Canterbury earthquake series. As such, by taking hazard levels accordingly (i.e. to correspond to the Christchurch CBD), modelling and analysing the structure, and considering the vulnerability and repair costs of its different components, it is possible to predict the expected losses of the aforementioned building. Using this information, cost-effective retrofit/rehabilitation strategy can be determined. This research found that more often than not, it would be beneficial to improve the performance of valuable non-structural components, such as partitions. Although it is true that improving such elements will increase the initial costs, over time, the benefits gained from reduced losses should be expected to overcome the initial costs. Aftershocks do increase the predicted losses of a building even in lower intensities due to the fact that non-structural components can get damaged at such low intensities. By comparing losses computed with and without consideration of aftershocks for a range of historical earthquakes, it was found that the ratio between losses due to main shock with aftershocks to the losses due to the main shock only tended to increase with increasing main shock magnitude. This may be due to the fact that larger magnitude earthquakes tend to generate larger magnitude aftershocks and as those aftershocks happen within a region around the main shock, they are more likely to cause intense shaking and additional damage. In addition to this observation, it was observed that the most significant component of loss of the case study building was the non-structural partition walls.
This report presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region during the Canterbury earthquake sequence. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the shallow geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site is based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS to better understand stratigraphy and soil properties, which are important in seismic site response. In this regard, CPT, SPT and borehole data, shear wave velocity (Vs) profiles, and horizontal to vertical spectral ratio measurements (H/V) in close vicinity to the SMS were used to develop representative soil profiles at each site. NZS1170.5 (SNZ 2004) site subsoil classifications were updated using Vs and SPT N60 criteria. Site class E boundaries were treated as a sliding scale rather than as a discrete boundary to account for locations with similar site effects potential, an approach which was shown to result in a better delineation between the site classes. SPT N60 values often indicate a stiffer site class than the Vs data for softer soil sites, highlighting the disparity between the two site investigation techniques. Both SPT N60 and Vs based site classes did not always agree with the original site classifications. This emphasises the importance of having detailed site‐specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N60 and Vs. Liquefaction triggering assessments were carried out for the Darfield and Christchurch earthquakes, and compared against observed liquefaction surface manifestations and ground motions characteristics at each SMS. In general, the characteristics of the recorded ground motions at each site correlate well with the triggering analyses. However, at sites that likely liquefied at depth (as indicated by triggering analyses and/or inferred from the characteristics of the recorded surface acceleration time series), the presence of a non‐liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects.
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’.
The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
During 2010 and 2011, major earthquakes caused widespread damage and the deaths of 185 people in the city of Christchurch. Damaged school buildings resulted in state intervention which required amendment of the Education Act of 1989, and the development of ‘site sharing agreements’ in undamaged schools to cater for the needs of students whose schools had closed. An effective plan was also developed for student assessment through establishing an earthquake impaired derived grade process. Previous research into traditional explanations of educational inequalities in the United Kingdom, the United States of America, and New Zealand were reviewed through various processes within three educational inputs: the student, the school and the state. Research into the impacts of urban natural disasters on education and education inequalities found literature on post disaster education systems but nothing could be found that included performance data. The impacts of the Canterbury earthquakes on educational inequalities and achievement were analysed over 2009-2012. The baseline year was 2009, the year before the first earthquake, while 2012 is seen as the recovery year as no schools closed due to seismic events and there was no state intervention into the education of the region. National Certificate of Educational Achievement (NCEA) results levels 1-3 from thirty-four secondary schools in the greater Christchurch region were graphed and analysed. Regression analysis indicates; in 2009, educational inequalities existed with a strong positive relationship between a school’s decile rating and NCEA achievement. When schools were grouped into decile rankings (1-10) and their 2010 NCEA levels 1-3 results were compared with the previous year, the percentage of change indicates an overall lower NCEA achievement in 2010 across all deciles, but particularly in lower decile schools. By contrast, when 2011 NCEA results were compared with those of 2009, as a percentage of change, lower decile schools fared better. Non site sharing schools also achieved higher results than site sharing schools. State interventions, had however contributed towards student’s achieving national examinations and entry to university in 2011. When NCEA results for 2012 were compared to 2009 educational inequalities still exist, however in 2012 the positive relationship between decile rating and achievement is marginally weaker than in 2009. Human ethics approval was required to survey one Christchurch secondary school community of students (aged between 12 and 18), teachers and staff, parents and caregivers during October 2011. Participation was voluntary and without incentives, 154 completed questionnaires were received. The Canterbury earthquakes and aftershocks changed the lives of the research participants. This school community was displaced to another school due to the Christchurch earthquake on 22 February 2011. Research results are grouped under four geographical perspectives; spatial impacts, socio-economic impacts, displacement, and health and wellbeing. Further research possibilities include researching the lag effects from the Canterbury earthquakes on school age children.
Recent earthquakes have highlighted the vulnerability of existing structure to seismic loading. Current seismic retrofit strategies generally focus on increasing the strength/stiffness in order to upgrade the seismic performance of a structure or element. A typical drawback of this approach is that the demand on the structural and sub-structural elements can be increased. This is of particular importance when considering the foundation capacity, which may already be insufficient to allow the full capacity of the existing wall to develop (due to early codes being gravity load orientated). In this thesis a counter intuitive but rational seismic retrofit strategy, termed "selective weakening" is introduced and investigated. This is the first stage of an ongoing research project underway at the University of Canterbury which is focusing on developing selective weakening techniques for the seismic retrofit of reinforced concrete structures. In this initial stage the focus is on developing selective weakening for the seismic retrofit of structural walls. This is performed using a series of experimental, analytical and numerical investigations. A procedure for the assessment of existing structural walls is also compiled, based on the suggestions of currently available code provisions. A selective weakening intervention is performed within an overall performance-based retrofit approach with the aim of improving the inelastic behaviour by first reducing the strength/stiffness of specific members within the structural system. This will be performed with the intention of modifying a shear type behaviour towards a flexural type behaviour. As a result the demand on the structural member will be reduced. Once weakening has been implemented the designer can use the wide range of techniques and materials available (e.g. use of FRP, jacketing or shotcrete) to ensure that adequate characteristics are achieved. Whilst performing this it has to be assured that the structure meets specific performance criteria and the principles of capacity design. A target of the retrofit technique is the ability to introduce the characteristics of recently developed high performance seismic resisting systems, consisting of a self centring and dissipative behaviour (commonly referred to as a hybrid system). In this thesis, results of experimental investigations performed on benchmark and selectively weakened walls are discussed. The investigations consisted of quasi-static cyclic uni-directional tests on two benchmark and two retrofitted cantilever walls. The first benchmark wall is detailed as typical of pre-1970's construction practice. An equivalent wall is retrofitted using a selective weakening approach involving a horizontal cut at foundation level to allow for a rocking response. The second benchmark wall represents a more severe scenario where the inelastic behaviour is dominated by shear. A retrofit solution involving vertically segmenting the wall to improve the ductility and retain gravity carrying capacity by inducing a flexural response is implemented. Numerical investigations on a multi-storey wall system are performed using non linear time history analysis on SDOF and MDOF lumped plasticity models, representing an as built and retrofitted prototype structure. Calibration of the hysteretic response to experimental results is carried out (accounting for pinching and strength degradation). The sensitivity of maximum and residual drifts to p-delta and strength degradation is monitored, along with the sensitivity of the peak base shear to higher mode affects. The results of the experimental and analytical investigations confirmed the feasibility and viability of the proposed retrofit technique, towards improving the seismic performance of structural walls.
In the period between September 2010 and December 2011, Christchurch was shaken by a series of strong earthquakes including the MW7.1 4 September 2010, Mw 6.2 22 February 2011, MW6.2 13 June 2011 and MW6.0 23 December 2011 earthquakes. These earthquakes produced very strong ground motions throughout the city and surrounding areas that resulted in soil liquefaction and lateral spreading causing substantial damage to buildings, infrastructure and the community. The stopbank network along the Kaiapoi and Avon River suffered extensive damage with repairs projected to take several years to complete. This presented an opportunity to undertake a case-study on a regional scale of the effects of liquefaction on a stopbank system. Ultimately, this information can be used to determine simple performance-based concepts that can be applied in practice to improve the resilience of river protection works. The research presented in this thesis draws from data collected following the 4th September 2010 and 22nd February 2011 earthquakes. The stopbank damage is categorised into seven key deformation modes that were interpreted from aerial photographs, consultant reports, damage photographs and site visits. Each deformation mode provides an assessment of the observed mechanism of failure behind liquefaction-induced stopbank damage and the factors that influence a particular style of deformation. The deformation modes have been used to create a severity classification for the whole stopbank system, being ‘no or low damage’ and ‘major or severe damage’, in order to discriminate the indicators and factors that contribute to ‘major to severe damage’ from the factors that contribute to all levels of damage a number of calculated, land damage, stopbank damage and geomorphological parameters were analysed and compared at 178 locations along the Kaiapoi and Avon River stopbank systems. A critical liquefiable layer was present at every location with relatively consistent geotechnical parameters (cone resistance (qc), soil behaviour type (Ic) and Factor of Safety (FoS)) across the study site. In 95% of the cases the critical layer occurred within two times the Height of the Free Face (HFF,). A statistical analysis of the geotechnical factors relating to the critical layer was undertaken in order to find correlations between specific deformation modes and geotechnical factors. It was found that each individual deformation mode involves a complex interplay of factors that are difficult to represent through correlative analysis. There was, however, sufficient data to derive the key factors that have affected the severity of deformation. It was concluded that stopbank damage is directly related to the presence of liquefaction in the ground materials beneath the stopbanks, but is not critical in determining the type or severity of damage, instead it is merely the triggering mechanism. Once liquefaction is triggered it is the gravity-induced deformation that causes the damage rather than the shaking duration. Lateral spreading and specifically the depositional setting was found to be the key aspect in determining the severity and type of deformation along the stopbank system. The presence or absence of abandoned or old river channels and point bar deposits was found to significantly influence the severity and type of deformation. A review of digital elevation models and old maps along the Kaiapoi River found that all of the ‘major to severe’ damage observed occurred within or directly adjacent to an abandoned river channel. Whilst a review of the geomorphology along the Avon River showed that every location within a point bar deposit suffered some form of damage, due to the depositional environment creating a deposit highly susceptible to liquefaction.
The recent Canterbury earthquake sequence in 2010-2011 highlighted a uniquely severe level of structural damage to modern buildings, while confirming the high vulnerability and life threatening of unreinforced masonry and inadequately detailed reinforced concrete buildings. Although the level of damage of most buildings met the expected life-safety and collapse prevention criteria, the structural damage to those building was beyond economic repair. The difficulty in the post-event assessment of a concrete or steel structure and the uneconomical repairing costs are the big drivers of the adoption of low damage design. Among several low-damage technologies, post-tensioned rocking systems were developed in the 1990s with applications to precast concrete members and later extended to structural steel members. More recently the technology was extended to timber buildings (Pres-Lam system). This doctoral dissertation focuses on the experimental investigation and analytical and numerical prediction of the lateral load response of dissipative post-tensioned rocking timber wall systems. The first experimental stages of this research consisted of component testing on both external replaceable devices and internal bars. The component testing was aimed to further investigate the response of these devices and to provide significant design parameters. Post-tensioned wall subassembly testing was then carried out. Firstly, quasi-static cyclic testing of two-thirds scale post-tensioned single wall specimens with several reinforcement layouts was carried out. Then, an alternative wall configuration to limit displacement incompatibilities in the diaphragm was developed and tested. The system consisted of a Column-Wall-Column configuration, where the boundary columns can provide the support to the diaphragm with minimal uplifting and also provide dissipation through the coupling to the post-tensioned wall panel with dissipation devices. Both single wall and column-wall-column specimens were subjected to drifts up to 2% showing excellent performance, limiting the damage to the dissipating devices. One of the objectives of the experimental program was to assess the influence of construction detailing, and the dissipater connection in particular proved to have a significant influence on the wall’s response. The experimental programs on dissipaters and wall subassemblies provided exhaustive data for the validation and refinement of current analytical and numerical models. The current moment-rotation iterative procedure was refined accounting for detailed response parameters identified in the initial experimental stage. The refined analytical model proved capable of fitting the experimental result with good accuracy. A further stage in this research was the validation and refinement of numerical modelling approaches, which consisted in rotational spring and multi-spring models. Both the modelling approaches were calibrated versus the experimental results on post-tensioned walls subassemblies. In particular, the multi-spring model was further refined and implemented in OpenSEES to account for the full range of behavioural aspects of the systems. The multi-spring model was used in the final part of the dissertation to validate and refine current lateral force design procedures. Firstly, seismic performance factors in accordance to a Force-Based Design procedure were developed in accordance to the FEMA P-695 procedure through extensive numerical analyses. This procedure aims to determine the seismic reduction factor and over-strength factor accounting for the collapse probability of the building. The outcomes of this numerical analysis were also extended to other significant design codes. Alternatively, Displacement-Based Design can be used for the determination of the lateral load demand on a post-tensioned multi-storey timber building. The current DBD procedure was used for the development of a further numerical analysis which aimed to validate the procedure and identify the necessary refinements. It was concluded that the analytical and numerical models developed throughout this dissertation provided comprehensive and accurate tools for the determination of the lateral load response of post-tensioned wall systems, also allowing the provision of design parameters in accordance to the current standards and lateral force design procedures.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.
Deconstruction, at the end of the useful life of a building, produces a considerable amount of materials which must be disposed of, or be recycled / reused. At present, in New Zealand, most timber construction and demolition (C&D) material, particularly treated timber, is simply waste and is placed in landfills. For both technical and economic reasons (and despite the increasing cost of landfills), this position is unlikely to change in the next 10 – 15 years unless legislation dictates otherwise. Careful deconstruction, as opposed to demolition, can provide some timber materials which can be immediately re-used (eg. doors and windows), or further processed into other components (eg. beams or walls) or recycled (‘cascaded’) into other timber or composite products (e.g. fibre-board). This reusing / recycling of materials is being driven slowly in NZ by legislation, the ‘greening’ of the construction industry and public pressure. However, the recovery of useful material can be expensive and uneconomic (as opposed to land-filling). In NZ, there are few facilities which are able to sort and separate timber materials from other waste, although the soon-to-be commissioned Burwood Resource Recovery Park in Christchurch will attempt to deal with significant quantities of demolition waste from the recent earthquakes. The success (or otherwise) of this operation should provide good information as to how future C&D waste will be managed in NZ. In NZ, there are only a few, small scale facilities which are able to burn waste wood for energy recovery (e.g. timber mills), and none are known to be able to handle large quantities of treated timber. Such facilities, with constantly improving technology, are being commissioned in Europe (often with Government subsidies) and this indicates that similar bio-energy (co)generation will be established in NZ in the future. However, at present, the NZ Government provides little assistance to the bio-energy industry and the emergence worldwide of shale-gas reserves is likely to push the economic viability of bio-energy further into the future. The behaviour of timber materials placed in landfills is complex and poorly understood. Degrading timber in landfills has the potential to generate methane, a potent greenhouse gas, which can escape to the atmosphere and cancel out the significant benefits of carbon sequestration during tree growth. Improving security of landfills and more effective and efficient collection and utilisation of methane from landfills in NZ will significantly reduce the potential for leakage of methane to the atmosphere, acting as an offset to the continuing use of underground fossil fuels. Life cycle assessment (LCA), an increasingly important methodology for quantifying the environmental impacts of building materials (particularly energy, and global warming potential (GWP)), will soon be incorporated into the NZ Green Building Council Greenstar rating tools. Such LCA studies must provide a level playing field for all building materials and consider the whole life cycle. Whilst the end-of-life treatment of timber by LCA may establish a present-day base scenario, any analysis must also present a realistic end-of-life scenario for the future deconstruction of any 6 new building, as any building built today will be deconstructed many years in the future, when very different technologies will be available to deal with construction waste. At present, LCA practitioners in NZ and Australia place much value on a single research document on the degradation of timber in landfills (Ximenes et al., 2008). This leads to an end-of-life base scenario for timber which many in the industry consider to be an overestimation of the potential negative effects of methane generation. In Europe, the base scenario for wood disposal is cascading timber products and then burning for energy recovery, which normally significantly reduces any negative effects of the end-of-life for timber. LCA studies in NZ should always provide a sensitivity analysis for the end-of-life of timber and strongly and confidently argue that alternative future scenarios are realistic disposal options for buildings deconstructed in the future. Data-sets for environmental impacts (such as GWP) of building materials in NZ are limited and based on few research studies. The compilation of comprehensive data-sets with country-specific information for all building materials is considered a priority, preferably accounting for end-of-life options. The NZ timber industry should continue to ‘champion’ the environmental credentials of timber, over and above those of the other major building materials (concrete and steel). End-of-life should not be considered the ‘Achilles heel’ of the timber story.
In September 2010 and February 2011, the Canterbury region experienced devastating earthquakes with an estimated economic cost of over NZ$40 billion (Parker and Steenkamp, 2012; Timar et al., 2014; Potter et al., 2015). The insurance market played an important role in rebuilding the Canterbury region after the earthquakes. Homeowners, insurance and reinsurance markets and New Zealand government agencies faced a difficult task to manage the rebuild process. From an empirical and theoretic research viewpoint, the Christchurch disaster calls for an assessment of how the insurance market deals with such disasters in the future. Previous studies have investigated market responses to losses in global catastrophes by focusing on the insurance supply-side. This study investigates both demand-side and supply-side insurance market responses to the Christchurch earthquakes. Despite the fact that New Zealand is prone to seismic activities, there are scant previous studies in the area of earthquake insurance. This study does offer a unique opportunity to examine and document the New Zealand insurance market response to catastrophe risk, providing results critical for understanding market responses after major loss events in general. A review of previous studies shows higher premiums suppress demand, but how higher premiums and a higher probability of risk affect demand is still largely unknown. According to previous studies, the supply of disaster coverage is curtailed unless the market is subsidised, however, there is still unsettled discussion on why demand decreases with time from the previous disaster even when the supply of coverage is subsidised by the government. Natural disaster risks pose a set of challenges for insurance market players because of substantial ambiguity associated with the probability of such events occurring and high spatial correlation of catastrophe losses. Private insurance market inefficiencies due to high premiums and spatially concentrated risks calls for government intervention in the provision of natural disaster insurance to avert situations of noninsurance and underinsurance. Political economy considerations make it more likely for government support to be called for if many people are uninsured than if few people are uninsured. However, emergency assistance for property owners after catastrophe events can encourage most property owners to not buy insurance against natural disaster and develop adverse selection behaviour, generating larger future risks for homeowners and governments. On the demand-side, this study has developed an intertemporal model to examine how demand for insurance changes post-catastrophe, and how to model it theoretically. In this intertemporal model, insurance can be sought in two sequential periods of time, and at the second period, it is known whether or not a loss event happened in period one. The results show that period one demand for insurance increases relative to the standard single period model when the second period is taken into consideration, period two insurance demand is higher post-loss, higher than both the period one demand and the period two demand without a period one loss. To investigate policyholders experience from the demand-side perspective, a total of 1600 survey questionnaires were administered, and responses from 254 participants received representing a 16 percent response rate. Survey data was gathered from four institutions in Canterbury and is probably not representative of the entire population. The results of the survey show that the change from full replacement value policy to nominated replacement value policy is a key determinant of the direction of change in the level of insurance coverage after the earthquakes. The earthquakes also highlighted the plight of those who were underinsured, prompting policyholders to update their insurance coverage to reflect the estimated cost of re-building their property. The survey has added further evidence to the existing literature, such as those who have had a recent experience with disaster loss report increased risk perception if a similar event happens in future with females reporting a higher risk perception than males. Of the demographic variables, only gender has a relationship with changes in household cover. On the supply-side, this study has built a risk-based pricing model suitable to generate a competitive premium rate for natural disaster insurance cover. Using illustrative data from the Christchurch Red-zone suburbs, the model generates competitive premium rates for catastrophe risk. When the proposed model incorporates the new RMS high-definition New Zealand Earthquake Model, for example, insurers can find the model useful to identify losses at a granular level so as to calculate the competitive premium. This study observes that the key to the success of the New Zealand dual insurance system despite the high prevalence of catastrophe losses are; firstly the EQC’s flat-rate pricing structure keeps private insurance premiums affordable and very high nationwide homeowner take-up rates of natural disaster insurance. Secondly, private insurers and the EQC have an elaborate reinsurance arrangement in place. By efficiently transferring risk to the reinsurer, the cost of writing primary insurance is considerably reduced ultimately expanding primary insurance capacity and supply of insurance coverage.
In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.
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.
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.
Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.
In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.
