On 22 February 2011, Canterbury and its largest city Christchurch experienced its second major earthquake within six months. The region is facing major economic and organisational challenges in the aftermath of these events. Approximately 25% of all buildings in the Christchurch CBD have been “red tagged” or deemed unsafe to enter. The New Zealand Treasury estimates that the combined cost of the February earthquake and the September earthquake is approximately NZ$15 billion[2]. This paper examines the national and regional economic climate prior to the event, discusses the immediate economic implications of this event, and the challenges and opportunities faced by organisations affected by this event. In order to facilitate recovery of the Christchurch area, organisations must adjust to a new norm; finding ways not only to continue functioning, but to grow in the months and years following these earthquakes. Some organisations relocated within days to areas that have been less affected by the earthquakes. Others are taking advantage of government subsidised aid packages to help retain their employees until they can make long-term decisions about the future of their organisation. This paper is framed as a “report from the field” in order to provide insight into the early recovery scenario as it applies to organisations affected by the February 2011 earthquake. It is intended both to inform and facilitate discussion about how organisations can and should pursue recovery in Canterbury, and how organisations can become more resilient in the face of the next crisis.
To this extent, modern buildings generally demonstrated good resistance to collapse during the recent earthquakes in New Zealand. However, damage to non-structural elements (NSE) has been persistent during these events. NSEs include secondary systems or components attached to the floors, roofs, and walls of a building or industrial facility that are not explicitly designed to participate in the main vertical or lateral load-bearing mechanism of the structure. They play a major role in the operational and functional aspects of buildings and contribute a major portion of the building’s overall cost. Therefore, they are expected to accommodate the effects of seismic actions such as drifts and accelerations. Typical examples of NSEs include internal non-loadbearing partitions, suspended ceilings, sprinkler piping systems, architectural claddings, building contents, mechanical/electrical equipment, and furnishings. The main focus of this thesis is the drift sensitive NSEs: precast concrete cladding panels and internal partition walls. Even though most precast concrete cladding panels performed well from a life-safety point of view during recent earthquakes in NZ, some collapsed panels posed a significant threat to life safety. It is, therefore, important that the design and detailing of the panel-to-structure connections ensure that their strength and displacement capacity are adequate to meet the corresponding seismic demands, at least during design level earthquakes. In contrast, the partition wall is likely to get damaged and lose serviceability at a low inter-story drift unless designed to accommodate the relative deformations between them and the structure. Partition walls suffered wide-ranging damage such as screw failures, diagonal cracking, detachments to the gypsum linings, and anchorage failures during the 2011 Canterbury Earthquake Sequence in NZ. Therefore, the thesis is divided into two parts. Part I of the thesis focuses on developing novel low-damage precast concrete cladding panel connections, i.e. “rocking” connection details comprising vertically slotted steel embeds and weld plates. The low-damage seismic performance of novel “rocking” connection details is verified through experimental tests comprising uni-directional, bi-directional, and multi-storey scaled quasi-static cyclic tests. Comparison with the seismic performance of traditional panel connections reported in the literature demonstrated the system’s significantly improved seismic resilience. Furthermore, the finite element models of panel connections and sealants are developed in ABAQUS. The force-drift responses of the “rocking” panel system modelled in SAP2000 is compared with the experimental results to evaluate their accuracy and validity. Part II of the thesis focuses on a) understanding the seismic performance of traditional rigid timber-framed partition wall, b) development and verification of low-damage connections (i.e. “rocking” connection details comprising of dual-slot tracks), and c) seismic evaluation of partition walls with a novel “bracketed and slotted” connections (comprising of innovative fastener and plastic bracket named Flexibracket) under uni-directional and bidirectional quasi-static cyclic loadings. Moreover, parametric investigation of the partition walls was conducted through several experimental tests to understand better the pros and cons of the rocking connection details. The experimental results have confirmed that the implementation of the proposed low damage solutions of precast cladding panels and internal partition walls can significantly reduce their damage in a building.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
The University of Canterbury is known internationally for the Origins of New Zealand English (ONZE) corpus (see Gordon et al 2004). ONZE is a large collection of recordings from people born between 1851 and 1984, and it has been widely utilised for linguistic and sociolinguistic research on New Zealand English. The ONZE data is varied. The recordings from the Mobile Unit (MU) are interviews and were collected by members of the NZ Broadcasting service shortly after the Second World War, with the aim of recording stories from New Zealanders outside the main city centres. These were supplemented by interview recordings carried out mainly in the 1990s and now contained in the Intermediate Archive (IA). The final ONZE collection, the Canterbury Corpus, is a set of interviews and word-list recordings carried out by students at the University of Canterbury. Across the ONZE corpora, there are different interviewers, different interview styles and a myriad of different topics discussed. In this paper, we introduce a new corpus – the QuakeBox – where these contexts are much more consistent and comparable across speakers. The QuakeBox is a corpus which consists largely of audio and video recordings of monologues about the 2010-2011 Canterbury earthquakes. As such, it represents Canterbury speakers’ very recent ‘danger of death’ experiences (see Labov 2013). In this paper, we outline the creation and structure of the corpus, including the practical issues involved in storing the data and gaining speakers’ informed consent for their audio and video data to be included.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage occurred to unreinforced masonry buildings throughout the region during the mainshock and subsequent large aftershocks. Particularly extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, a geotechnical reconnaissance was conducted over a period of six days (10–15 September 2010) by a team of geotechnical/earthquake engineers and geologists from New Zealand and USA (GEER team: Geo-engineering Extreme Event Reconnaissance). JGS (Japanese Geotechnical Society) members from Japan also participated in the reconnaissance team from 13 to 15 September 2010. The NZ, GEER and JGS members worked as one team and shared resources, information and logistics in order to conduct thorough and most efficient reconnaissance covering a large area over a very limited time period. This report summarises the key evidence and findings from the reconnaissance.
The Canterbury earthquake sequence in New Zealand’s South Island induced widespread liquefaction phenomena across the Christchurch urban area on four occasions (4 Sept 2010; 22 Feb; 13 June; 23 Dec 2011), that resulted in widespread ejection of silt and fine sand. This impacted transport networks as well as infiltrated and contaminated the damaged storm water system, making rapid clean-up an immediate post-earthquake priority. In some places the ejecta was contaminated by raw sewage and was readily remobilised in dry windy conditions, creating a long-term health risk to the population. Thousands of residential properties were inundated with liquefaction ejecta, however residents typically lacked the capacity (time or resources) to clean-up without external assistance. The liquefaction silt clean-up response was co-ordinated by the Christchurch City Council and executed by a network of contractors and volunteer groups, including the ‘Farmy-Army’ and the ‘Student-Army’. The duration of clean-up time of residential properties and the road network was approximately 2 months for each of the 3 main liquefaction inducing earthquakes; despite each event producing different volumes of ejecta. Preliminary cost estimates indicate total clean-up costs will be over NZ$25 million. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill since the beginning of the Canterbury earthquakes sequence. The liquefaction clean-up experience in Christchurch following the 2010-2011 earthquake sequence has emerged as a valuable case study to support further analysis and research on the coordination, management and costs of large volume deposition of fine grained sediment in urban areas.
The development of cheap, whilst effective and relatively non-invasive structural retrofit techniques for existing non-ductile reinforced concrete (RC) structures still remains the most challenging issue for a wide implementation on a macro scale. Seismic retrofit is too often being confused as purely structural strengthening. As part of a six-years national project on “Seismic retrofit solutions for NZ multi-storey building”, focus has been given at the University of Canterbury on the development of a counter-intuitive retrofit strategy for earthquake vulnerable existing rc frame, based on a “selective weakening” (SW) approach. After an overview of the SW concept, this paper presents the experimental and numerical validation of a SW retrofit strategy for earthquake vulnerable existing RC frame with particular focus on the exterior beam-column (b-c) joints. The exterior b-c joint is a critically vulnerable region in many existing pre-1970s RC frames. By selectively weakening the beam by cutting the bottom longitudinal reinforcements and/or adding external pre-stressing to the b-c joint, a more desirable inelastic mechanism can be attained, leading to improved global seismic performance. The so-called SW retrofit is implemented on four 2/3-scaled exterior RC b-c joint subassemblies, tested under quasi-static cyclic loading at the University of Canterbury. Complemented by refined 3D Finite Element (FE) models and dynamic time-history analyses results, the experimental results have shown the potential of a simple and cost-effective yet structurally efficient structural rehabilitation technique. The research also demonstrated the potential of advanced 3D fracture-mechanics-based microplane concrete modelling for refined FE analysis of non-ductile RC b-c joints.
Previous research has found that the capacity to self-regulate is associated with a number of positive life outcomes and deficits in self-regulation have been linked with poorer life outcomes. Therefore, parent and child self-regulation is an important focus of the Positive Parenting Program for Teenagers (Teen Triple P). The aim of this study was to investigate if Group Teen Triple P was effective in promoting parental self-regulation and adolescent behaviour change in families affected by the earthquakes in Canterbury NZ between 2010 and 2012. METHOD: Five families with teenagers aged 12-16 years were recruited from among families participating in a Group Teen Triple P program specifically implemented by the education authorities for parents self-reporting long-term negative effects of the earthquakes on their family. A single-case multiple-baseline across participants design was used to examine change in target teenager behaviour. Measures of self-regulation skill acquisition were taken using a coding scheme devised for the study from transcripts of three telephone consultations and from three family discussions at pre-intervention, mid-intervention, and post-intervention. Parents and their child also completed questionnaires addressing adolescent functioning, the parent-adolescent relationship and parenting at pre- and post-intervention. RESULTS: The multiple-baseline data showed that parents were successful at changing targeted behaviour for their child. Analysis of the telephone consultations and family discussions showed that parents increased their self-regulation skills over the therapy period and there was positive change in adolescent behavior reported on the Strengths and Difficulties Questionnaire. Additionally, the results suggested that higher rates and levels of self-regulation in the parents were associated with greater improvements in adolescent behaviour. CONCLUSION: This study demonstrated that the Group Teen Triple P -Program was effective in promoting self-regulation in parents and behaviour change in adolescents, specifically in a post-disaster context.
Non-structural elements (NSEs) have frequently proven to contribute to significant losses sustained from earthquakes in the form of damage, downtime, injury and death. In New Zealand (NZ), the 2010 and 2011 Canterbury Earthquake Sequence (CES), the 2013 Seddon and Cook Strait earthquake sequence and the 2016 Kaikoura earthquake were major milestones in this regard as significant damage to building NSEs both highlighted and further reinforced the importance of NSE seismic performance to the resilience of urban centres. Extensive damage in suspended ceilings, partition walls, façades and building services following the CES was reported to be partly due to erroneous seismic design or installation or caused by intervening elements. Moreover, the low-damage solutions developed for structural systems sometimes allow for relatively large inter-story drifts -compared to conventional designs- which may not have been considered in the seismic design of NSEs. Having observed these shortcomings, this study on suspended ceilings was carried out with five main goals: i) Understanding the seismic performance of the system commonly used in NZ; ii) Understanding the transfer of seismic design actions through different suspended ceiling components, iii) Investigating potential low-damage solutions; iii) Evaluating the compatibility of the current ceiling system with other low-damage NSEs; and iv) Investigating the application of numerical analysis to simulate the response of ceiling systems. The first phase of the study followed a joint research work between the University of Canterbury (UC) in NZ, and the Politecnico Di Milano, in Italy. The experimental ceiling component fragility curves obtained in this existing study were employed to produce analytical fragility curves for a perimeter-fixed ceiling of a given size and weight, with grid acceleration as the intensity measure. The validity of the method was proven through comparisons between this proposed analytical approach with the recommended procedures in proprietary products design guidelines, as well as experimental fragility curves from other studies. For application to engineering design practice, and using fragility curves for a range of ceiling lengths and weights, design curves were produced for estimating the allowable grid lengths for a given demand level. In the second phase of this study, three specimens of perimeter-fixed ceilings were tested on a shake table under both sinusoidal and random floor motion input. The experiments considered the relationship between the floor acceleration, acceleration of the ceiling grid, the axial force induced in the grid members, and the effect of boundary conditions on the transfer of these axial forces. A direct correlation was observed between the axial force (recorded via load cells) and the horizontal acceleration measured on the ceiling grid. Moreover, the amplification of floor acceleration, as transferred through ceiling components, was examined and found (in several tests) to be greater than the recommended factor for the design of ceilings provided in the NZ earthquake loadings standard NZS1170.5. However, this amplification was found to be influenced by the pounding interactions between the ceiling grid members and the tiles, and this amplification diminished considerably when the high frequency content was filtered out from the output time histories. The experiments ended with damage in the ceiling grid connection at an axial force similar to the capacity of these joints previously measured through static tests in phase one. The observation of common forms of damage in ceilings in earthquakes triggered the monotonic experiments carried out in the third phase of this research with the objective of investigating a simple and easily applicable mitigation strategy for existing or new suspended ceilings. The tests focused on the possibility of using proprietary cross-shaped clip elements ordinarily used to provide seismic gap as a strengthening solution for the weak components of a ceiling. The results showed that the solution was effective under both tension and compression loads through increasing load bearing capacity and ductility in grid connections. The feasibility of a novel type of suspended ceiling called fully-floating ceiling system was investigated through shaking table tests in the next phase of this study with the main goal of isolating the ceiling from the surrounding structure; thereby arresting the transfer of associated seismic forces from the structure to the ceiling. The fully-floating ceiling specimen was freely hung from the floor above lacking any lateral bracing and connections with the perimeter. Throughout different tests, a satisfactory agreement between the fully-floating ceiling response and simple pendulum theory was demonstrated. The addition of isolation material in perimeter gaps was found effective in inducing extra damping and protecting the ceiling from pounding impact; resulting in much reduced ceiling displacements and accelerations. The only form of damage observed throughout the random floor motion tests and the sinusoidal tests was a panel dislodgement observed in a test due to successive poundings between the ceiling specimen and the surrounding beams at resonant frequencies. Partition walls as the first effective NSE in direct interaction with ceilings were the topic of the final experimental phase. Low-damage drywall partitions proposed in a previous study in the UC were tested with two common forms of suspended ceiling: braced and perimeter-fixed. The experiments investigated the in-plane and out-of-plane performance of the low-damage drywall partitions, as well as displacement compatibility between these walls and the suspended ceilings. In the braced ceiling experiment, where no connection was made between ceiling grids and surrounding walls no damage in the grid system or partitions was observed. However, at high drift values panel dislodgement was observed on corners of the ceiling where the free ends of grids were not restrained against spreading. This could be prevented by framing the grid ends using a perimeter angle that is riveted only to the grid members while keeping sufficient clearance from the perimeter walls. In the next set of tests with the perimeter-fixed ceiling, no damage was observed in the ceiling system or the drywalls. Based on the results of the experiments it was concluded that the tested ceiling had enough flexibility to accommodate the relative displacement between two perpendicular walls up to the inter-storey drifts achieved. The experiments on perimeter-fixed ceilings were followed by numerical simulations of the performance of these ceilings in a finite element model developed in the structural analysis software, SAP2000. This model was relatively simple and easy to develop and was able to replicate the experimental results to a reasonable degree. Filtering was applied to the experimental output to exclude the effect of high frequency noise and tile-grid impact. The developed model generally simulated the acceleration responses well but underestimated the peak ceiling grid accelerations. This was possibly because the peak values in time histories were affected by impact occurring at very short periods. The model overestimated the axial forces in ceiling grids which was assumed to be caused by the initial assumptions made about the tributary area or constant acceleration associated with each grid line in the direction of excitation. Otherwise, the overall success of the numerical modelling in replicating the experimental results implies that numerical modelling using conventional structural analysis software could be used in engineering practice to analyse alternative ceiling geometries proposed for application to varying structural systems. This however, needs to be confirmed through similar analyses on other ceiling examples from existing instrumented buildings during real earthquakes. As the concluding part of this research the final phase addressed the issues raised following the review of existing ceiling standards and guidelines. The applicability of the research findings to current practice and their implications were discussed. Finally, an example was provided for the design of a suspended ceiling utilising the new knowledge acquired in this research.
The purpose of this thesis is to evaluate the seismic response of the UC Physics Building based on recorded ground motions during the Canterbury earthquakes, and to use the recorded response to evaluate the efficacy of various conventional structural analysis modelling assumptions. The recorded instrument data is examined and analysed to determine how the UC Physics Building performed during the earthquake-induced ground motions. Ten of the largest earthquake events from the 2010-11 Canterbury earthquake sequence are selected in order to understand the seismic response under various levels of demand. Peak response amplitude values are found which characterise the demand from each event. Spectral analysis techniques are utilised to find the natural periods of the structure in each orthogonal direction. Significant torsional and rocking responses are also identified from the recorded ground motions. In addition, the observed building response is used to scrutinise the adequacy of NZ design code prescriptions for fundamental period, response spectra, floor acceleration and effective member stiffness. The efficacy of conventional numerical modelling assumptions for representing the UC Physics Building are examined using the observed building response. The numerical models comprise of the following: a one dimensional multi degree of freedom model, a two dimensional model along each axis of the building and a three dimensional model. Both moderate and strong ground motion records are used to examine the response and subsequently clarify the importance of linear and non-linear responses and the inclusion of base flexibility. The effects of soil-structure interaction are found to be significant in the transverse direction but not the longitudinal direction. Non-linear models predict minor in-elastic behaviour in both directions during the 4 September 2010 Mw 7.1 Darfield earthquake. The observed torsional response is found to be accurately captured by the three dimensional model by considering the interaction between the UC Physics Building and the adjacent structure. With the inclusion of adequate numerical modelling assumptions, the structural response is able to be predicted to within 10% for the majority of the earthquake events considered.
On November 14 2016 a magnitude 7.8 earthquake struck the south island of New Zealand. The earthquake lasted for just two minutes with severe seismic shaking and damage in the Hurunui and Kaikōura districts. Although these are predominantly rural areas, with scattered small towns and mountainous topography, they also contain road and rail routes that are essential parts of the national transport infrastructure. This earthquake and the subsequent recovery are of particular significance as they represent a disaster following in close proximity to another similar disaster, with the Canterbury earthquakes occurring in a neighboring district five years earlier. The research used an inductive qualitative case study to explore the nature of the Kaikōura recovery. That recovery process involved a complex interplay between the three parties; (a) the existing local government in the district, (b) central government agencies funding the recovery of the local residents and the national transport infrastructure, and (c) recovery leaders arriving with recent expertise from the earlier Canterbury disaster. It was evident that three groups: locals, government, and experts represented a multi-party governance debate in which the control of the Kaikōura earthquake recovery was shared amongst them. Each party had their own expertise, adgenda and networks that they brought to the Kaikōura recovery, but this created tensions between external expertise and local, community leadership. Recent earthquake research suggests that New Zealand is currently in the midst of an earthquake cluster, with further seismic disasters likely to occur in relatively close succession. This is likely to be compounded by the increasing frequency of other natural disasters with the effects of climate change. The present study investigates a phenomenon that may become increasingly common, with the transfer of disaster expertise from one event to another, and the interface between those experts with local and national government in directing recoveries. The findings of this study have implications for practitioners and policy makers in NZ and other countries where disasters are experienced in close spatial and temporal proximity.
Organisations locate strategically within Business Districts (CBDs) in order to cultivate their image, increase their profile, and improve access to customers, suppliers, and services. While CBDs offer an economic benefit to organisations, they also present a unique set of hazard vulnerabilities and planning challenges for businesses. As of May 2012, the Christchurch CBD has been partially cordoned off for over 14 months. Economic activity within the cordoned CBD, which previously contained 6,000 businesses and over 51,000 workers, has been significantly diminished and organisations have been forced to find new ways of operating. The vulnerabilities and resilience of CBDs not only influences outcomes for CBD organisations, but also the broader interconnected (urban/regional/national) system. A CBD is a hub of economic, social, and built infrastructure within a network of links and nodes. When the hub is disrupted all of the people, objects, and transactions that usually flow into and out of the hub must be redirected elsewhere. In an urban situation this means traffic jams in peripheries of the city, increased prices of commercial property, and capital flight; all of which are currently being faced in Canterbury. This report presents the lessons learned from organisations in CBDs affected by the Canterbury earthquakes. Here we focus on the Christchurch CBD; however, several urban town centres were extensively disrupted by the earthquakes. The statistics and discussion presented in this report are based on the results of an ongoing study conducted by Resilient Organisations (www.resorgs.org.nz). The data was captured using two questionnaire surveys of Canterbury organisations (issued November 2010 and May 2011), interviews with key informants, and in-depth case studies of organisations. Several industry sectors were sampled, and geographic samples of organisations in the Christchurch CBD, Lyttelton, and the Kaiapoi town centre were also collected. Results in this report describing “non-CBD organisations” refer to all organisations outside of the Christchurch CBD, Lyttelton, and Kaiapoi town centres.
INTRODUCTION: Connections between environmental factors and mental health issues have been postulated in many different countries around the world. Previously undertaken research has shown many possible connections between these fields, especially in relation to air quality and extreme weather events. However, research on this subject is lacking in New Zealand, which is difficult to analyse as an overall nation due to its many micro-climates and regional differences.OBJECTIVES: The aim of this study and subsequent analysis is to explore the associations between environmental factors and poor mental health outcomes in New Zealand by region and predict the number of people with mental health-related illnesses corresponding to the environmental influence.METHODS: Data are collected from various public-available sources, e.g., Stats NZ and Coronial services of New Zealand, which comprised four environmental factors of our interest and two mental health indicators data ranging from 2016 up until 2020. The four environmental factors are air pollution, earthquakes, rainfall and temperature. Two mental health indicators include the number of people seen by District Health Boards (DHBs) for mental health reasons and the statistics on suicide deaths. The initial analysis is carried out on which regions were most affected by the chosen environmental factors. Further analysis using Auto-Regressive Integrated Moving Average(ARIMA) creates a model based on time series of environmental data to generate estimation for the next two years and mental health projected from the ridge regression.RESULTS: In our initial analysis, the environmental data was graphed along with mental health outcomes in regional charts to identify possible associations. Different regions of New Zealand demonstrate quite different relationships between the environmental data and mental health outcomes. The result of later analysis predicts that the suicide rate and DHB mental health visits may increase in Wellington, drop-in Hawke's Bay and slightly increase in Canterbury for the year 2021 and 2022 with different environmental factors considered.CONCLUSION: It is evident that the relationship between environmental and mental health factors is regional and not national due to the many micro-climates that exist around the nation. However, it was observed that not all factors displayed a good relationship between the regions. We conclude that our hypotheses were partially correct, in that increased air pollution was found to correlate to increased mental health-related DHB visits. Rainfall was also highly correlated to some mental health outcomes. Higher levels of rainfall reduced DHB visits and suicide rates in some areas of the country.
Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.
This thesis addresses the topic of local bond behaviour in RC structures. The mechanism of bond refers to the composite action between deformed steel reinforcing bars and the surrounding concrete. Bond behaviour is an open research topic with a wide scope, particularly because bond it is such a fundamental concept to structural engineers. However, despite many bond-related research findings having wide applications, the primary contribution of this research is an experimental evaluation of the prominent features of local bond behaviour and the associated implications for the seismic performance of RC structures. The findings presented in this thesis attempt to address some structural engineering recommendations made by the Canterbury Earthquakes Royal Commission following the 2010-2011 Canterbury (New Zealand) earthquake sequence. A chapter of this thesis discusses the structural behaviour of flexure-dominated RC wall structures with an insufficient quantity of longitudinal reinforcement, among other in situ conditions, that causes material damage to predominantly occur at a single crack plane. In this particular case, the extent of concrete damage and bond deterioration adjacent to the crack plane will influence the ductility capacity that is effectively provided by the reinforcing steel. As a consequence of these in situ conditions, some lightly reinforced wall buildings in Christchurch lost their structural integrity due to brittle fracture of the longitudinal reinforcement. With these concerning post-earthquake observations in mind, there is the underlying intention that this thesis presents experimental evidence of bond behaviour that allows structural engineers to re-assess their confidence levels for the ability of lightly reinforced concrete structures to achieve the life-safety seismic performance objective the ultimate limit state. Three chapters of this thesis are devoted to the experimental work that was conducted as the main contribution of this research. Critical details of the experimental design, bond testing method and test programme are reported. The bond stress-slip relationship was studied through 75 bond pull-out tests. In order to measure the maximum local bond strength, all bond tests were carried out on deformed reinforcing bars that did not yield as the embedded bond length was relatively short. Bond test results have been presented in two separate chapters in which 48 monotonic bond tests and 27 cyclic bond tests are presented. Permutations of the experiments include the loading rate, cyclic loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and position of the embedded bond region within the specimen (close or far away to the free surface). The parametric study showed that the concrete strength significantly influences the maximum bond strength and that it is reasonable to normalise the bond stress by the square-root of the concrete compressive strength, √(f'c). The generalised monotonic bond behaviour is described within. An important outcome of the research is that the measured bond strength and stiffness was higher than stated by the bond stress-slip relationship in the fib Model Code 2010. To account for these observed differences, an alternative model is proposed for the local monotonic bond stress-slip relationship. Cyclic bond tests showed a significant proportion of the total bond degradation occurs after the loading cycle in the peak bond strength range, which is when bond slip has exceeded 0.5 mm. Subsequent loading to constant slip values showed a linear relationship between the amount of bond strength degradation and the log of the number of cycles that were applied. To a greater extent, the cyclic bond deterioration depends on the bond slip range, regardless of whether the applied load cycling is half- or fully-reversed. The observed bond deterioration and hysteretic energy dissipated during cyclic loading was found to agree reasonably well between these cyclic tests with different loading protocols. The cyclic bond deterioration was also found to be reasonably consistent exponential damage models found in the literature. This research concluded that the deformed reinforcing bars used in NZ construction, embedded in moderate to high strength concrete, are able to develop high local bond stresses that are mobilised by a small amount of local bond slip. Although the relative rib geometry was not varied within this experimental programme, a general conclusion of this thesis is that deformed bars currently available in NZ have a relative rib bearing area that is comparatively higher than the test bars used in previous international research. From the parametric study it was found that the maximum monotonic bond strength is significant enhanced by dynamic loading rates. Experimental evidence of high bond strength and initial bond stiffness generally suggests that only a small amount of local bond slip that can occur when the deformed test bar was subjected to large tension forces. Minimal bond slip and bond damage limits the effective yielding length that is available for the reinforcing steel to distribute inelastic material strains. Consequently, the potential for brittle fracture of the reinforcement may be a more problematic and widespread issue than is apparent to structural engineers. This research has provided information that improve the reliability of engineering predictions (with respect to ductility capacity) of maximum crack widths and the extent of bond deterioration that might occur in RC structures during seismic actions.
In recent years, rocking isolation has become an effective approach to improve seismic performance of steel and reinforced concrete structures. These systems can mitigate structural damage through rigid body displacement and thus relatively low requirements for structural ductility, which can significantly improve seismic resilience of structures and reduce repairing costs after strong earthquakes. A number of base rocking structural systems with only a single rocking interface have been proposed. However, these systems can have significant high mode effect for high rise structures due to the single rocking interface. This RObust BUilding SysTem (ROBUST) project is a collaborative China-New Zealand project sponsored by the International Joint Research Laboratory of Earthquake Engineering (ILEE), Tongji University, and a number of agencies and universities within New Zealand including the BRANZ, Comflor, Earthquake Commission, HERA, QuakeCoRE, QuakeCentre, University of Auckland, and the University of Canterbury. A number of structural configurations will be tested [1, 2], and non-structural elements including ceilings, infilling walls, glazed curtain walls, precast concrete panels, piping system will also be tested in this project [3]. Within this study, a multiple rocking column steel structural system was proposed and investigated mainly by Tongji team with assistance of NZ members. The concept of rocking column system initiates from the structure of Chinese ancient wooden pagoda. In some of Chinese wooden pagodas, there are continuous core columns hanged only at the top of each pagoda, which is not connected to each stories. This core column can effectively avoid collapse of the whole structure under large storey drifts. Likewise, there are also central continuous columns in the newly proposed steel rocking column system, which can avoid weak story failure mechanism and make story drifts more uniform. In the proposed rocking column system, the structure can switch between an elastic rigidly connected moment resisting frame and a controlled rocking column system when subjected to strong ground motion excitations. The main seismic energy can be dissipated by asymmetric friction beam–column connections, thereby effectively reducing residual displacement of the structure under seismic loading without causing excessive damage to structural members. Re–centering of the structure is provided not only by gravity load carried by rocking columns, but also by mould coil springs. To investigate dynamic properties of the proposed system under different levels of ground excitations, a full-scale threestory steel rocking column structural system with central continuous columns is to be tested using the International joint research Laboratory of Earthquake Engineering (ILEE) facilities, Shanghai, China and an analytical model is established. A finite element model is also developed using ABAQUS to simulate the structural dynamic responses. The rocking column system proposed in this paper is shown to produce resilient design with quick repair or replacement.
In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.
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.
In most design codes, infill walls are considered as non-structural elements and thus are typically neglected in the design process. The observations made after major earthquakes (Duzce 1999, L’Aquila 2009, Christchurch 2011) have shown that even though infill walls are considered to be non-structural elements, they interact with the structural system during seismic actions. In the case of heavy infill walls (i.e. clay brick infill walls), the whole behaviour of the structure may be affected by this interaction (i.e. local or global structural failures such as soft storey mechanism). In the case of light infill walls (i.e. non-structural drywalls), this may cause significant economical losses. To consider the interaction of the structural system with the ‘non-structural ’infill walls at design stage may not be a practical approach due to the complexity of the infill wall behaviour. Therefore, the purpose of the reported research is to develop innovative technological solutions and design recommendations for low damage non-structural wall systems for seismic actions by making use of alternative approaches. Light (steel/timber framed drywalls) and heavy (unreinforced clay brick) non-structural infill wall systems were studied by following an experimental/numerical research programme. Quasi-static reverse cyclic tests were carried out by utilizing a specially designed full scale reinforced concrete frame, which can be used as a re-usable bare frame. In this frame, two RC beams and two RC columns were connected by two un-bonded post tensioning bars, emulating a jointed ductile frame system (PRESSS technology). Due to the rocking behaviour at the beam-column joint interfaces, this frame was typically a low damage structural solution, with the post-tensioning guaranteeing a linear elastic behaviour. Therefore, this frame could be repeatedly used in all of the tests carried out by changing only the infill walls within this frame. Due to the linear elastic behaviour of this structural bare frame, it was possible to extract the exact behaviour of the infill walls from the global results. In other words, the only parameter that affected the global results was given by the infill walls. For the test specimens, the existing practice of construction (as built) for both light and heavy non-structural walls was implemented. In the light of the observations taken during these tests, modified low damage construction practices were proposed and tested. In total, seven tests were carried out: 1) Bare frame , in order to confirm its linear elastic behaviour. 2) As built steel framed drywall specimen FIF1-STFD (Light) 3) As built timber framed drywall specimen FIF2-TBFD (Light) 4) As built unreinforced clay brick infill wall specimen FIF3-UCBI (Heavy) 5) Low damage steel framed drywall specimen MIF1-STFD (Light) 6) Low damage timber framed drywall specimen MIF2-TBFD (Light) 7) Low damage unreinforced clay brick infill wall specimen MIF5-UCBI (Heavy) The tests of the as built practices showed that both drywalls and unreinforced clay brick infill walls have a low serviceability inter-storey drift limit (0.2-0.3%). Based on the observations, simple modifications and details were proposed for the low damage specimens. The details proved to be working effectively in lowering the damage and increasing the serviceability drift limits. For drywalls, the proposed low damage solutions do not introduce additional cost, material or labour and they are easily applicable in real buildings. For unreinforced clay brick infill walls, a light steel sub-frame system was suggested that divides the infill panel zone into smaller individual panels, which requires additional labour and some cost. However, both systems can be engineered for seismic actions and their behaviour can be controlled by implementing the proposed details. The performance of the developed details were also confirmed by the numerical case study analyses carried out using Ruaumoko 2D on a reinforced concrete building model designed according to the NZ codes/standards. The results have confirmed that the implementation of the proposed low damage solutions is expected to significantly reduce the non-structural infill wall damage throughout a building.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
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.
The Acheron rock avalanche is located in the Red Hill valley almost 80 km west of Christchurch and is one of 42 greywacke-derived rock avalanches identified in the central Southern Alps. It overlies the Holocene active Porters Pass Fault; a component of the Porters Pass-Amberley Fault Zone which extends from the Rakaia River to beyond the Waimakariri River. The Porters Pass Fault is a dextral strike-slip fault system viewed as a series of discontinuous fault scarps. The location of the fault trace beneath the deposit suggests it may represent a possible source of seismic shaking resulting in the formation of the Acheron rock avalanche. The rock mass composition of the rock avalanche source scar is Torlesse Supergroup greywacke consisting of massive sandstone and thinly bedded mudstone sequences dipping steeply north into the centre of the source basin. A stability analysis identified potential instability along shallow north dipping planar defects, and steep south dipping toppling failure planes. The interaction of the defects with bedding is considered to have formed conditions for potential instability most likely triggered by a seismic event. The dTositional area of the rock avalanche covers 7.2 x 105 m2 with an estimated volume of 9 x 10 m3 The mobilised rock mass volume was calculated at 7.5 x 106 m3• Run out of the debris from the top of the source scar to the distal limit reached 3500m, descending over a vertical fall of almost 700m with an estimated Fahrboschung of 0.2. The run out of the rock avalanche displayed moderate to high mobility, travelling at an estimated maximum velocity of 140-160 km/hour. The rapid emplacement of the deposit is confirmed by highly fragmented internal composition and burial of forest vegetation New radiocarbon ages from buried wood retrieved from the base of Acheron rock avalanche deposit represents an emplacement age closely post-dating (Wk 12094) 1152 ± 51 years B.P. This differs significantly from a previous radiocarbon age of (NZ547) 500 ± 69 years B.P. and modal lichenometry and weathering-rind thickness ages of approximately 460 ± 10 yrs and 490 ± 50 years B.P. The new age shows no resemblance to an earthquake event around 700- 500 years B.P. on the Porters Pass-Amberley Fault Zone. The DAN run out simulation using a friction model rheology successfully replicated the long run out and velocity of the Acheron rock avalanche using a frictron angle of 27° and high earth pressure coefficients of 5.5, 5.2, and 5.9. The elevated earth pressure coefficients represent dispersive pressures derived from dynamic fragmentation of the debris within the mobile rock avalanche, supporting the hypothesis of Davies and McSaveney (2002). The DAN model has potential applications for areas prone to large-scale instability in the elevated slopes and steep waterways of the Southern Alps. A paleoseismic investigation of a newly identified scarp of the Porters Pass Fault partially buried by the rock avalanche was conducted to identify any evidence of a coseismic relationship to the Acheron rock avalanche. This identified three-four fault traces striking at 078°, and a sag pond displaying a sequence of overbank deposits containing two buried soils representing an earthquake event horizon. A 40cm vertical offset of the ponded sediment and lower buried soil horizqn was recorded, which was dated to (Wk 13112 charcoal in palosol) 653 ± 54 years B.P. and (Wk 13034 palosol) 661 ± 34 years B.P. The evidence indicates a fault rupture occurred along the Porters Pass Fault, west of Porters Pass most likely extending to the Red Lakes terraces, post-dating 700 years B.P., resulting in 40cm of vertical displacement and an unknown component of dextral strike slip movement. This event post dates the event one (1000 ± 100 years B.P) at Porters Pass previously considered to represent the most recent rupture along the fault line. This points to a probable source for resetting of the modal weathering-rind thicknesses and lichen size populations in the Red Hill valley and possibly the Red Lakes terraces. These results suggest careful consideration must be given to the geomorphic and paleoseismic history of a specific site when applying surface dating techniques and furthermore the origin of dates used in literature and their useful range should be verified. An event at 700-500 years B.P did not trigger the Acheron rock avalanche as previously assumed supporting Howard's conclusions. The lack of similar aged rupture evidence in either of the Porters Pass and Coleridge trenches supports Howard's hypothesis of segmentation of the Porters Pass Fault; where rupture occurs along one fault segment but not along another. The new rock avalanche age closely post-dating 1200-1100 years B.P. resembles the poorly constrained event one rupture age of 1700-800 years B.P for the Porters Pass Fault and the tighter constrained Round Top event of 1010 ± 50 years B.P. on the Alpine Fault. Eight other rock avalanche deposits spread across the central Southern Alps also resemble the new ages however are unable to be assigned specific earthquake events due to the large associated error bars of± 270 years. This clustering of ages does represent compelling lines of evidence for large magnitude earthquake events occurring over the central Southern Alps. The presence of a rock avalanche deposit does not signify an earthquake based on the historical evidence in the Southern Alps however clustering of ages does suggest that large Mw >7 earthquakes occurred across the Southern Alps between 1200-900 years BP.
As a global phenomenon, many cities are undergoing urban renewal to accommodate rapid growth in urban population. However, urban renewal can struggle to balance social, economic, and environmental outcomes, whereby economic outcomes are often primarily considered by developers. This has important implications for urban forests, which have previously been shown to be negatively affected by development activities. Urban forests serve the purpose of providing ecosystem services and thus are beneficial to human wellbeing. Better understanding the effect of urban renewal on city trees may help improve urban forest outcomes via effective management and policy strategies, thereby maximising ecosystem service provision and human wellbeing. Though the relationship between certain aspects of development and urban forests has received consideration in previous literature, little research has focused on how the complete property redevelopment cycle affects urban forest dynamics over time. This research provides an opportunity to gain a comprehensive understanding of the effect of residential property redevelopment on urban forest dynamics, at a range of spatial scales, in Christchurch, New Zealand following a series of major earthquakes which occurred in 2010 – 2011. One consequence of the earthquakes is the redevelopment of thousands of properties over a relatively short time-frame. The research quantifies changes in canopy cover city-wide, as well as, tree removal, retention, and planting on individual residential properties. Moreover, the research identifies the underlying reasons for these dynamics, by exploring the roles of socio-economic and demographic factors, the spatial relationships between trees and other infrastructure, and finally, the attitudes of residential property owners. To quantify the effect of property redevelopment on canopy cover change in Christchurch, this research delineated tree canopy cover city-wide in 2011 and again in 2015. An object-based image analysis (OBIA) technique was applied to aerial imagery and LiDAR data acquired at both time steps, in order to estimate city-wide canopy cover for 2011 and 2015. Changes in tree canopy cover between 2011 and 2015 were then spatially quantified. Tree canopy cover change was also calculated for all meshblocks (a relatively fine-scale geographic boundary) in Christchurch. The results show a relatively small magnitude of tree canopy cover loss, city-wide, from 10.8% to 10.3% between 2011 and 2015, but a statistically significant change in mean tree canopy cover across all the meshblocks. Tree canopy cover losses were more likely to occur in meshblocks containing properties that underwent a complete redevelopment cycle, but the loss was insensitive to the density of redevelopment within meshblocks. To explore property-scale individual tree dynamics, a mixed-methods approach was used, combining questionnaire data and remote sensing analysis. A mail-based questionnaire was delivered to residential properties to collect resident and household data; 450 residential properties (321 redeveloped, 129 non- redeveloped) returned valid questionnaires and were identified as analysis subjects. Subsequently, 2,422 tree removals and 4,544 tree retentions were identified within the 450 properties; this was done by manually delineating individual tree crowns, based on aerial imagery and LiDAR data, and visually comparing the presence or absence of these trees between 2011 and 2015. The tree removal rate on redeveloped properties (44.0%) was over three times greater than on non-redeveloped properties (13.5%) and the average canopy cover loss on redeveloped properties (52.2%) was significantly greater than on non-redeveloped properties (18.8%). A classification tree (CT) analysis was used to model individual tree dynamics (i.e. tree removal, tree retention) and candidate explanatory variables (i.e. resident and household, economic, land cover, and spatial variables). The results indicate that the model including land cover, spatial, and economic variables had the best predicting ability for individual tree dynamics (accuracy = 73.4%). Relatively small trees were more likely to be removed, while trees with large crowns were more likely to be retained. Trees were most likely to be removed from redeveloped properties with capital values lower than NZ$1,060,000 if they were within 1.4 m of the boundary of a redeveloped building. Conversely, trees were most likely to be retained if they were on a property that was not redeveloped. The analysis suggested that the resident and household factors included as potential explanatory variables did not influence tree removal or retention. To conduct a further exploration of the relationship between resident attitudes and actions towards trees on redeveloped versus non-redeveloped properties, this research also asked the landowners from the 450 properties that returned mail questionnaires to indicate their attitudes towards tree management (i.e. tree removal, tree retention, and tree planting) on their properties. The results show that residents from redeveloped properties were more likely to remove and/or plant trees, while residents from non- redeveloped properties were more likely to retain existing trees. A principal component analysis (PCA) was used to explore resident attitudes towards tree management. The results of the PCA show that residents identified ecosystem disservices (e.g. leaf litter, root damage to infrastructure) as common reasons for tree removal; however, they also noted ecosystem services as important reasons for both tree planting and tree retention on their properties. Moreover, the reasons for tree removal and tree planting varied based on whether residents’ property had been redeveloped. Most tree removal occurred on redeveloped properties because trees were in conflict with redevelopment, but occurred on non- redeveloped properties because of perceived poor tree health. Residents from redeveloped properties were more likely to plant trees due to being aesthetically pleasing or to replace trees removed during redevelopment. Overall, this research adds to, and complements, the existing literature on the effects of residential property redevelopment on urban forest dynamics. The findings of this research provide empirical support for developing specific legislation or policies about urban forest management during residential property redevelopment. The results also imply that urban foresters should enhance public education on the ecosystem services provided by urban forests and thus minimise the potential for tree removal when undertaking property redevelopment.
