Recent tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural mitigation, warning systems, education and evacuation planning). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and evacuation behaviour is rare when compared to findings for evacuations from other hazards. More knowledge is required to document the full evacuation process, including responses to warnings, pre-evacuation actions, evacuation dynamics, and the return home. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However, to date, tsunami evacuation models have largely focused on the timings of evacuations, rather than behaviours of those evacuating. In this research, survey data was gathered from coastal communities in Banks Peninsula and Christchurch, New Zealand, relating to behaviours and actions during the November 14th 2016 Kaikōura earthquake tsunami. Survey questions asked about immediate actions following the earthquake shaking, reactions to tsunami warnings, pre-evacuation actions, evacuation dynamics and details on congestion. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour, and was further used to develop a realistic evacuation model prototype to evaluate the capacity of the roading network in Banks Peninsula during a tsunami evacuation. The evacuation model incorporated tsunami risk management strategies that have been implemented by local authorities, and exposure and vulnerability data, alongside the empirical data collected from the survey. This research enhances knowledge of tsunami evacuation behaviour and reactions to tsunami warnings, and can be used to refine evacuation planning to ensure that people can evacuate efficiently, thereby reducing their tsunami exposure and personal risk.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
The November 2016 MW 7.8 Kaikōura Earthquake initiated beneath the North Culverden basin on The Humps fault and propagated north-eastwards, rupturing at least 17 faults along a cumulative length of ~180 km. The geomorphic expression of The Humps Fault across the Emu Plains, along the NW margin of Culverden basin, comprises a series of near-parallel strands separated by up to 3 km across strike. The various strands strike east to east-northeast and have been projected to mainly dip steeply to the south in seismic data (~80°). In this area, the fault predominantly accommodates right-lateral slip, with uplift and subsidence confined to releasing and restraining bends and step-overs at a range of scales. The Kaikōura event ruptured pre-existing fault scarps along the Emu Plains, which had been partly identified prior to the earthquake. Geomorphology and faulting expression of The Humps Fault on The Emu Plains was mapped, along with faulting related structures which did not rupture in the 2016 earthquake. Fault ruptures strands are combined into sections and the kinematic deformation of sections analysed to provide a moment tensor fault plane solution. This fault plane solution is consistent with the regional principal horizontal shortening direction (PHS) of ~115°, similar to seismic focal mechanism solutions of some of the nearby aftershocks of the Kaikōura earthquake, and similar to the adjacent Hope Fault. To constrain the timing of paleoseismic events, a trench was excavated across the fault where it crossed a late Quaternary alluvial fan. Mapping of stratigraphy exposed in the trench walls, and dating of variably deformed strata, constrains the pre-historic earthquake event history at the trench site. The available data provides evidence for at least three paleo-earthquakes within the last 15.1 ka, with a possible fourth (penultimate) event. These events are estimated to have occurred at 7.7-10.3 ka, 10.3-14.8 ka, and one or more events that are older than ~15.1 ka. Some evidence suggests an additional penultimate event between 1850 C.E and 7.7 ka. Time-integrated slip-rates at three locations on the fault are measured using paleo-channels as piercing points. These sites give horizontal slip rates of 0.57 ± 0.1 mm/year, 0.49 ± 0.1 mm/year and one site constrains a minimum of between 0.1 - 0.4 mm/year. Two vertical slip-rates are calculated to be constrained to a maximum of 0.2 ± 0.02 mm/year at one site and between 0.02 and 0.1 mm/year at another site. Prior to this study, The Humps fault had only been partially documented in reconnaissance level mapping in the district, and no previous paleoseismic or slip rate data had been reported. This project has provided a detailed fault zone tectonic geomorphic map and established new slip-rate and paleoseismic data. The results highlight that The Humps fault plays an important role in regional seismicity and in accommodating plate boundary deformation across the North Canterbury region.
The performance of buildings in recent New Zealand earthquakes (Canterbury, Seddon and Kaikōura), delivered stark lessons on seismic resilience. Most of our buildings, with a few notable exceptions, performed as our Codes intended them to, that is, to safeguard people from injury. Many buildings only suffered minor structural damage but were unable to be reused and occupied for significant periods of time due to the damage and failure of non-structural elements. This resulted in substantial economic losses and major disruptions to our businesses and communities. Research has attributed the damage to poor overall design coordination, inadequate or lack of seismic restraints for non structural elements and insufficient clearances between building components to cater for the interaction of non structural elements under seismic actions. Investigations have found a clear connection between the poor performance of non-structural elements and the issues causing pain in the industry (procurement methods, risk aversion, the lack of clear understanding of design and inspection responsibility and the need for better alignment of the design codes to enable a consistent integrated design approach). The challenge to improve the seismic performance of non structural elements in New Zealand is a complex one that cuts across a diverse construction industry. Adopting the key steps as recommended in this paper is expected to have significant co-benefits to the New Zealand construction industry, with improvements in productivity alongside reductions in costs and waste, as the rework which plagues the industry decreases.
Recent global tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural defences, warning systems, education and evacuation measures). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and tsunami evacuation behaviour is rare when compared to findings about evacuations to avoid other sources of hazard. To date empirical research into tsunami evacuations has focused on evacuation rates, rather than other aspects of the evacuation process. More knowledge is required about responses to warnings, pre-evacuation actions, evacuation dynamics and the return home after evacuations. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However to date tsunami evacuation models have largely focused on timings of evacuations, rather than evacuation behaviours. This Masters research uses a New Zealand case study to reduce both of these knowledge gaps. Qualitative survey data was gathered from populations across coastal communities in Banks Peninsula and Christchurch, New Zealand, required to evacuate due to the tsunami generated by the November 14th 2016 Kaikōura Earthquake. Survey questions asked about reactions to tsunami warnings, actions taken prior to evacuating and movements during the 2016 tsunami evacuation. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour. Finally, it was used to develop an evacuation model for Banks Peninsula. Where appropriate, the modelling inputs were informed by the survey data. Three key findings were identified from the results of the evacuation behaviour survey. Although 38% of the total survey respondents identified the earthquake shaking as a natural cue for the tsunami, most relied on receiving official warnings, including sirens, to prompt evacuations. Respondents sought further official information to inform their evacuation decisions, with 39% of respondents delaying their evacuation in order to do so. Finally, 96% of total respondents evacuated by car. This led to congestion, particularly in more densely populated Christchurch city suburbs. Prior to this research, evacuation modelling had not been completed for Banks Peninsula. The results of the modelling showed that if evacuees know how to respond to tsunami warnings and where and how to evacuate, there are no issues. However, if there are poor conditions, including if people do not evacuate immediately, if there are issues with the roading network, or if people do not know where or how to evacuate, evacuation times increase with there being more bottlenecks leading out of the evacuation zones. The results of this thesis highlight the importance of effective tsunami education and evacuation planning. Reducing exposure to tsunami risk through prompt evacuation relies on knowledge of how to interpret tsunami warnings, and when, where and how to evacuate. Recommendations from this research outline the need for public education and engagement, and the incorporation of evacuation signage, information boards and evacuation drills. Overall these findings provide more comprehensive picture of tsunami evacuation behaviour and decision making based on empirical data from a recent evacuation, which can be used to improve tsunami risk management strategies. This empirical data can also be used to inform evacuation modelling to improve the accuracy and realism of the evacuation models.
A building boom in the 1980s allowed pre-stressed hollow-core floor construction to be widely adopted in New Zealand, even though the behaviour of these prefabricated elements within buildings was still uncertain. Inspections following the Canterbury and Kaikōura earthquakes has provided evidence of web-splitting, transverse cracking and longitudinal splitting on hollow-core units, confirming the susceptibility of these floors to undesirable failure modes. Hollow-core slabs are mainly designed to resist bending and shear. However, there are many applications in which they are also subjected to torsion. In New Zealand, hollow-core units contain no transverse reinforcement in the soffit concrete below the cells and no web reinforcement. Consequently, their dependable performance in torsion is limited to actions that they can resist before torsional cracking occurs. In previous work by the present authors, a three-dimensional FE modelling approach to study the shear flexural behaviour of precast pre-stressed hollow core units was developed and validated by full-scale experiments. This paper shows how the FE analyses have been extended to investigate the response of HC units subjected to torsional actions. Constitutive models, based on nonlinear fracture mechanics, have been used to numerically predict the torsional capacity of HC units and have been compared with experimental results. The results indicate that the numerical approach is promising and should be developed further as part of future research.
During the 2011 M7.8 Kaikōura earthquake, ground motions recorded near the epicentre showed a significant spatial variation. The Te Mara farm (WTMC) station, the nearest to the epicentre, recorded 1g and 2.7g of horizontal and vertical peak ground accelerations (PGA), respectively. The nearby Waiu Gorge (WIGC) station recorded a horizontal PGA of 0.8g. Interestingly, however, the Culverden Airlie Farm (CULC) station that was very closely located to WIGC recorded a horizontal PGA of only 0.25g. This poster demonstrates how the local geological condition could have contributed to the spatially variable ground motions observed in the North Canterbury, based on the results of recently conducted geophysical investigations. The surficial geology of this area is dominated by alluvial gravel deposits with traces of silt. A borehole log showed that the thickness of the sediments at WTMC is over 76 metres. Interestingly, the shear wave velocity (Vs) profiles obtained from the three strong motion sites suggest unusually high shear wave velocity of the gravelly sediments. The velocity of sediments and the lack of clear peaks in the horizontal-to-vertical (H/V) spectral ratio at WTMC suggest that the large ground motion observed at this station was likely caused by the proximity of the station to the causative fault itself; the site effect was likely insignificant. Comparisons of H/V spectral ratios and Vs profiles suggest that the sediment thickness is much smaller at WIGC compared with CULC; the high PGA at WIGC was likely influenced by the high-frequency amplification caused by the response of shallow sediments.
The use of post-earthquake cordons as a tool to support emergency managers after an event has been documented around the world. However, there is limited research that attempts to understand the use, effectiveness, inherent complexities, impacts and subsequent consequences of cordoning once applied. This research aims to fill that gap by providing a detailed understanding of first, the cordons and associated processes, and their implications in a post-earthquake scenario. We use a qualitative method to understand cordons through case studies of two cities where it was used in different temporal and spatial scales: Christchurch (2011) and Wellington (Kaikōura earthquake 2016), New Zealand. Data was collected through 21 expert interviews obtained through purposive and snowball sampling of key informants who were directly or indirectly involved in a decision-making role and/or had influence in relation to the cordoning process. The participants were from varying backgrounds and roles i.e. emergency managers, council members, business representatives, insurance representatives, police and communication managers. The data was transcribed, coded in Nvivo and then grouped based on underlying themes and concepts and then analyzed inductively. It is found that cordons are used primarily as a tool to control access for the purpose of life safety and security. But cordons can also be adapted to support recovery. Broadly, it can be synthesized and viewed based on two key aspects, ‘decision-making’ and ‘operations and management’, which overlap and interact as part of a complex system. The underlying complexity arises in large part due to the multitude of sectors it transcends such as housing, socio-cultural requirements, economics, law, governance, insurance, evacuation, available resources etc. The complexity further increases as the duration of cordon is extended.
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.
Nowadays the telecommunication systems’ performance has a substantial impact on our lifestyle. Their operationality becomes even more substantial in a post-disaster scenario when these services are used in civil protection and emergency plans, as well as for the restoration of all the other critical infrastructure. Despite the relevance of loss of functionality of telecommunication networks on seismic resilience, studies on their performance assessment are few in the literature. The telecommunication system is a distributed network made up of several components (i.e. ducts, utility holes, cabinets, major and local exchanges). Given that these networks cover a large geographical area, they can be easily subjected to the effects of a seismic event, either the ground shaking itself, or co-seismic events such as liquefaction and landslides. In this paper, an analysis of the data collected after the 2010-2011 Canterbury Earthquake Sequence (CES) and the 2016 Kaikoura Earthquake in New Zealand is conducted. Analysing these data, information gaps are critically identified regarding physical and functional failures of the telecommunication components, the timeline of repair/reconstruction activities and service recovery, geotechnical tests and land planning maps. Indeed, if these missing data were presented, they could aid the assessment of the seismic resilience. Thus, practical improvements in the post-disaster collection from both a network and organisational viewpoints are proposed through consultation of national and international researchers and highly experienced asset managers from Chorus. Finally, an outline of future studies which could guide towards a more resilient seismic performance of the telecommunication network is presented.
One of the most controversial issues highlighted by the 2010-2011 Christchurch earthquake series and more recently the 2016 Kaikoura earthquake, has been the evident difficulty and lack of knowledge and guidelines for: a) evaluation of the residual capacity damaged buildings to sustain future aftershocks; b) selection and implementation of a series of reliable repairing techniques to bring back the structure to a condition substantially the same as prior to the earthquake; and c) predicting the cost (or cost-effectiveness) of such repair intervention, when compared to fully replacement costs while accounting for potential aftershocks in the near future. As a result of such complexity and uncertainty (i.e., risk), in combination with the possibility (unique in New Zealand when compared to most of the seismic-prone countries) to rely on financial support from the insurance companies, many modern buildings, in a number exceeding typical expectations from past experiences at an international level, have ended up being demolished. This has resulted in additional time and indirect losses prior to the full reconstruction, as well as in an increase in uncertainty on the actual relocation of the investment. This research project provides the main end-users and stakeholders (practitioner engineers, owners, local and government authorities, insurers, and regulatory agencies) with comprehensive evidence-based information to assess the residual capacity of damage reinforced concrete buildings, and to evaluate the feasibility of repairing techniques, in order to support their delicate decision-making process of repair vs. demolition or replacement. Literature review on effectiveness of epoxy injection repairs, as well as experimental tests on full-scale beam-column joints shows that repaired specimens have a reduced initial stiffness compared with the undamaged specimen, with no apparent strength reduction, sometimes exhibiting higher displacement ductility capacities. Although the bond between the steel and concrete is only partially restored, it still allows the repaired specimen to dissipate at least the same amount of hysteretic energy. Experimental tests on buildings subjected to earthquake loading demonstrate that even for severe damage levels, the ability of the epoxy injection to restore the initial stiffness of the structure is significant. Literature review on damage assessment and repair guidelines suggests that there is consensus within the international community that concrete elements with cracks less than 0.2 mm wide only require cosmetic repairs; epoxy injection repairs of cracks less and 2.0 mm wide and concrete patching of spalled cover concrete (i.e., minor to moderate damage) is an appropiate repair strategy; and for severe damaged components (e.g., cracks greater than 2.0 mm wide, crushing of the concrete core, buckling of the longitudinal reinforcement) local replacement of steel and/or concrete in addition to epoxy crack injection is more appropriate. In terms of expected cracking patterns, non-linear finite element investigations on well-designed reinforced concrete beam-to-column joints, have shown that lower number of cracks but with wider openings are expected to occur for larger compressive concrete strength, f’c, and lower reinforcement content, ρs. It was also observed that the tensile concrete strength, ft, strongly affects the expected cracking pattern in the beam-column joints, the latter being more uniformly distributed for lower ft values. Strain rate effects do not seem to play an important role on the cracking pattern. However, small variations in the cracking pattern were observed for low reinforcement content as it approaches to the minimum required as per NZS 3101:2006. Simple equations are proposed in this research project to relate the maximum and residual crack widths with the steel strain at peak displacement, with or without axial load. A literature review on fracture of reinforcing steel due to low-cycle fatigue, including recent research using steel manufactured per New Zealand standards is also presented. Experimental results describing the influence of the cyclic effect on the ultimate strain capacity of the steel are also discussed, and preliminary equations to account for that effect are proposed. A literature review on the current practice to assess the seismic residual capacity of structures is also presented. The various factors affecting the residual fatigue life at a component level (i.e., plastic hinge) of well-designed reinforced concrete frames are discussed, and equations to quantify each of them are proposed, as well as a methodology to incorporate them into a full displacement-based procedure for pre-earthquake and post-earthquake seismic assessment.
Floor systems with precast concrete hollow-core units have been largely used in concrete buildings built in New Zealand during the 1980’s. Recent earthquakes, such as the Canterbury sequence in 2010-2011 and the Kaikoura earthquake in 2016, highlighted that this floor system can be highly vulnerable and potentially lead to the floor collapse. A series of research activities are in progress to better understand the seismic performance of floor diaphragms, and this research focuses on examining the performance of hollow core units running parallel to the walls of wall-resisting concrete structures. This study first focused on the development of fragility functions, which can be quickly used to assess likelihood of the hollow-core being able to survive given the buildings design drift, and secondly to determine the expected performance of hollow-core units that run parallel to walls, focusing on the alpha unit running by the wall. Fragility functions are created for a range of different parameters for both vertical dislocation and crack width that can be used as the basis of a quick analysis or loss estimation for the likely impact of hollow-core floors on building vulnerability and risk. This was done using past experimental tests, and the recorded damage. Using these results and the method developed by Baker fragility curves were able to be created for varying crack widths and vertical dislocations. Current guidelines for analysis of hollow-core unit incompatible displacements are based on experimental vertical displacement results from concrete moment resisting frame systems to determine the capacity of hollow-core elements. To investigate the demands on hollow-core units in a wall-based structure, a fibre-element model in the software Seismostruct is created and subject to quasi-static cyclic loading, using elements which are verified from previous experimental tests. It is shown that for hollow-core units running by walls that the 10 mm displacement capacity used for hollow-core units running by a beam is insufficient for members running by walls and that shear analysis should be used. The fibre-element model is used to simulate the seismic demand induced on the floor system and has shown that the shear demand is a function of drift, wall length, hollow-core span, linking slab length and, to a minor extent, wall elongation.
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.
This thesis focuses on the role of legal preparedness for managing large-scale urban disasters in Aotearoa New Zealand. It uses the Auckland Volcanic Field as a case study to answer the question: ‘is New Zealand’s current legal framework prepared to respond to and recover from a large-scale urban disaster?’. The Auckland Volcanic Field was chosen as the main case study because a future eruption is a low likelihood, high-impact event that New Zealand is going to have to manage in the future. Case studies are a key feature of this thesis as both New Zealand based and overseas examples are used to explore the role of legal preparedness by identifying and investigating a range of legal issues that need to be addressed in advance of a future Auckland Volcanic Field eruption. Of particular interest is the impact of legal preparedness for the recovery phase. The New Zealand case studies include; Canterbury earthquake sequence 2010-2011, the Kaikōura earthquake 2016, the Auckland flooding 2018, and the North Island Severe Weather event 2023, which encompasses both the Auckland Anniversary weekend flooding and Cyclone Gabrielle. As New Zealand has not experienced a large-scale urban volcanic eruption, overseas examples are explored to provide insights into the legal issues that are volcano specific. The overseas volcanic case studies cover eruptions in Heimaey (Iceland), the Soufrière Hills (Montserrat and the Grenadines), La Soufrière (St Vincent) and Tungurahua (Ecuador). New Zealand’s past experiences highlight a trend for introducing post-event legal frameworks to manage recovery. Consequently, the current disaster management system is not prioritising legal preparedness and instead is choosing to rely on exceptional powers. Unsurprisingly, the introduction of new post-event recovery frameworks has repercussions. In New Zealand, new post-event legal frameworks are introduced swiftly under urgency, they contain broad unstructured decision-making powers, and are often flawed. As these exceptional new frameworks sit outside the ‘normal’ legal frameworks, they in effect create a parallel “shadow system”. Based on the evidence explored in this thesis it does not appear that Auckland’s current disaster management framework is prepared to deal with a large-scale urban event caused by an Auckland Volcanic Field eruption. Following this conclusion, it is the submission of this thesis that New Zealand’s current legal framework is not prepared to respond to and recover from a large-scale urban disaster. To become legally prepared, New Zealand needs to consider the legal tools required to manage large-scale urban disasters in advance. This will prevent the creation of a legal vacuum in the aftermath of disasters and the need for new recovery frameworks. Adopting a new attitude will require a change in approach towards legal preparedness which prioritises it, rather than sidelining it. This may also require changes within New Zealand’s disaster management system including the introduction of a formal monitoring mechanism, which will support and prioritise legal preparedness. This thesis has shown that not legally preparing for future disasters is a choice which carries significant consequences. None of these consequences are inevitable when managing large-scale disasters, however they are inevitable when frameworks are not legally prepared in advance. To not legally prepare, is to prepare to fail and thus create a disaster by choice.
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.
Ongoing climate change triggers increasing temperature and more frequent extreme events which could limit optimal performance of haliotids, affect their physiology and biochemistry as well as influencing their population structure. Haliotids are a valuable nearshore fishery in a number of countries and many are showing a collapse of stocks because of overexploitation, environmental changes, loss of habitat, and disease. The haliotid in New Zealand commonly referred to as the blackfoot pāua (Haliotis iris) contribute a large and critical cultural, recreational and economic resource. Little was known about pāua responses to increasing temperature and acute environmental factors, as well as information about population size structure in Kaikoura after the earthquake 2016 and in Banks Peninsula. The aims of this study were to investigate the effects of temperature on scope for growth (SfG); physiological and biochemical responses of pāua subjected to different combined stressors including acute temperature, acute salinity and progressive hypoxia; and describe population size structure and shell morphology in different environments in Kaikoura and Banks Peninsula. The main findings of the present study found that population size structures of pāua were site-specific, and the shell length and shell height ratio of 3.25 could distinguish between stunted and non-stunted populations. The study found that high water temperature resulted in a reduction in absorbed energy from food, an increase in respiration energy, and ammonia excretion energy. Surveys were conducted at six study sites around the Canterbury Region over three years in order to better understand the population size structure and shell morphology of pāua. The findings found that the population size structure at 6 sites differed. Both juveniles and adults were found in intertidal areas at five sites. However, at Cape Three Points, pāua were found only in subtidal zones. One of the sites, Little Port Cooper, had a stunted population where only two pāua reached 125 mm in length over three years. In addition, most pāua in Little Port Cooper and Cape Three Points were adults, while Seal Reef had mostly juveniles. Wakatu Quay and Omihi had a full size range of pāua. Oaro population was dominated with juveniles and sub-adults. Recruitment and growth of pāua were successful after the earthquake in 2016. Research into pāua shell morphologies also determined that shell dimensions differed between sites. The relationships of shell length to shell width were linear and the relationship of shell length to shell height was curvilinear. Interestingly, SL:SH ratio of 3.25 is able to be used to identify stunted and non-stunted populations for pāua larger than 90 mm in length. Little Port Cooper was a stunted population with mean SL:SH ratio being 3.16. In the laboratory, scope for growth of pāua was investigated at four different temperatures of 12oC, 15oC, 18oC and 21oC over four weeks’ acclimation. The current study has found that SfG of pāua highly depended on temperature. Absorbed energy and respiration energy accounted for the highest proportion of the SfG of pāua. The respiration energy of pāua accounted for approximately 36%, 40%, 49% and 69% of the absorbed energy at 12°C, 15°C, 18°C and 21°C, respectively. The pāua at all acclimation temperatures had a positive scope for growth. The study suggested that the SfG was highest at 15°C, while the value at 21°C was the lowest. However, SfG at 18°C and 21°C decreased after 14 days of acclimation. Because of maintaining almost unchanged oxygen consumption over four weeks’ acclimation, pāua showed their poor abilities to acclimate to an increase in temperature. Therefore, they may be more vulnerable in future warming scenarios. The physiological and biochemical responses of pāua toward different combined stressors included three experiments. In terms of the acute temperature experiment, pāua were acclimated at 12oC, 15oC, 18oC or 21oC for two weeks before stepwise exposure to four temperatures of 12oC, 15oC, 18oC and 21oC every 4 hours. The acute salinity change, pāua were acclimated at 12oC, 15oC or 18oC over two weeks. Pāua were then exposed to a stepwise decrease of salinity of 2‰ every two hours from 34 – 22‰. Regarding the declining oxygen level, pāua were acclimated at 15 oC or 18oC for two weeks before exposure to one of four temperatures at 12oC, 15oC, 18oC or 21oC in one hour. After that acute progressive hypoxia was studied in closed respirometers for around six hours. The findings showed that there were interactions between combined stressors, affecting physiology of pāua (metabolism and heart rate). This suggests that environmental factors do not have a separate effect, but they also have interactions that enhance negative effects on pāua. Also, both oxygen uptake and heart rate responded quickly to temperature change and increased with rising temperature. On the other hand, oxygen uptake and heart rate decreased with reducing salinity and progressive hypoxia (before critical oxygen tension - Pcrit). Pcrit over four acute temperature exposures, ranged between 30.2 and 80.0 mmHg, depending on the exposure temperature. Acclimation temperature, combined with acute temperature, salinity or hypoxia stress affected the biochemistry of pāua. Pāua are osmoconformers so decreased salinity resulted in reducing haemolymph ionic concentration and increasing body volume. They were hypo-ionic with respect to sodium and potassium over the salinity ranges of 34 - 22‰. Haemocyanin accounts for a large pecentage of haemolymph protein, so trends of protein followed haemocyanin. Pāua tended to store oxygen in haemocyanin under extreme salinity stress at 22‰ and extreme hypoxia around 10 mmHg, rather than in oxygen transport. In conclusion, pāua at different sites had different population structures and morphologies. Pāua are sensitive to environmental stressors. They consumed more oxygen at high temperatures because they do not have thermal acclimation capacity. They are also osmoconformers with haemolymph sodium and potassium decreasing with salinity medium. Under progressive hypoxia, pāua could regulate oxygen and heart rate until Pcrit depending on temperature. Acute environmental changes also disturbed haemolyph parameters. 12°C and 15°C could be in the range of optimal temperature with higher SfG and less stress when exposed to acute environmental changes. Meanwhile long term exposure to 21°C is likely to be outside of the optimal range for the pāua. With ongoing climate change, pāua populations are more vulnerable so conservation is necessary. The research contributes to improving fishery management, providing insights into different environmental stressors affecting the energy demand and physiological and biochemical responses of pāua. It also allow to predicting the growth patterns and responses of pāua to adapt to climate change.
