Search

found 29 results

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.

Research papers, University of Canterbury Library

A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.

Research papers, University of Canterbury Library

The latest two great earthquake sequences; 2010- 2011 Canterbury Earthquake and 2016 Kaikoura Earthquake, necessitate a better understanding of the New Zealand seismic hazard condition for new building design and detailed assessment of existing buildings. It is important to note, however, that the New Zealand seismic hazard map in NZS 1170.5.2004 is generalised in effort to cover all of New Zealand and limited to a earthquake database prior to 2001. This is “common” that site-specific studies typically provide spectral accelerations different to those shown on the national map (Z values in NZS 1170.5:2004); and sometimes even lower. Moreover, Section 5.2 of Module 1 of the Earthquake Geotechnical Engineering Practice series provide the guidelines to perform site- specific studies.

Research papers, University of Canterbury Library

Abstract This study provides a simplified methodology for pre-event data collection to support a faster and more accurate seismic loss estimation. Existing pre-event data collection frameworks are reviewed. Data gathered after the Canterbury earthquake sequences are analysed to evaluate the relative importance of different sources of building damage. Conclusions drawns are used to explore new approaches to conduct pre-event building assessment.

Research papers, University of Canterbury Library

This report presents the simplified seismic assessment of a case study reinforced concrete (RC) building following the newly developed and refined NZSEE/MBIE guidelines on seismic assessment (NZSEE/MBIE, semi-final draft 26 October 2016). After an overview of the step-by-step ‘diagnostic’ process, including an holistic and qualitative description of the expected vulnerabilities and of the assessment strategy/methodology, focus is given, whilst not limited, to the implementation of a Detailed Seismic Assessment (DSA) (NZSEE/MBIE, 2016c). The DSA is intended to provide a more reliable and consistent outcome than what can be provided by an initial seismic assessment (ISA). In fact, while the Initial Seismic Assessment (ISA), of which the Initial Evaluation Procedure is only a part of, is the more natural and still recommended first step in the overall assessment process, it is mostly intended to be a coarse evaluation involving as few resources as reasonably possible. It is thus expected that an ISA will be followed by a Detailed Seismic Assessment (DSA) not only where the threshold of 33%NBS is not achieved but also where important decisions are intended that are reliant on the seismic status of the building. The use of %NBS (% New Building Standard) as a capacity/demand ratio to describe the result of the seismic assessment at all levels of assessment procedure (ISA through to DSA) is deliberate by the NZSEE/MBIE guidelines (Part A) (NZSEE/MBIE 2016a). The rating for the building needs only be based on the lowest level of assessment that is warranted for the particular circumstances. Discussion on how the %NBS rating is to be determined can be found in Section A3.3 (NZSEE/MBIE 2016a), and, more specifically, in Part B for the ISA (NZSEE/MBIE 2016b) and Part C for the DSA (NZSEE/MBIE 2016c). As per other international approaches, the DSA can be based on several analysis procedures to assess the structural behaviour (linear, nonlinear, static or dynamic, force or displacement-based). The significantly revamped NZSEE 2016 Seismic Assessment Guidelines strongly recommend the use of an analytical (basically ‘by hand’) method, referred to the Simple Lateral Mechanism Analysis (SLaMA) as a first phase of any other numerically-based analysis method. Significant effort has thus been dedicated to provide within the NZSEE 2016 guidelines (NZSEE/MBIE 2016c) a step-by-step description of the procedure, either in general terms (Chapter 2) or with specific reference to Reinforced Concrete Buildings (Chapter 5). More specifically, extract from the guidelines, NZSEE “recommend using the Simple Lateral Mechanism Analysis (SLaMA) procedure as a first step in any assessment. While SLaMA is essentially an analysis technique, it enables assessors to investigate (and present in a simple form) the potential contribution and interaction of a number of structural elements and their likely effect on the building’s global capacity. In some cases, the results of a SLaMA will only be indicative. However, it is expected that its use should help assessors achieve a more reliable outcome than if they only carried out a detailed analysis, especially if that analysis is limited to the elastic range For complex structural systems, a 3D dynamic analysis may be necessary to supplement the simplified nonlinear Simple Lateral Mechanism Analysis (SLaMA).” This report presents the development of a full design example for the the implementation of the SLaMA method on a case study buildings and a validation/comparison with a non-linear static (pushover) analysis. The step-by-step-procedure, summarized in Figure 1, will be herein demonstrated from a component level (beams, columns, wall elements) to a subassembly level (hierarchy of strength in a beam-column joint) and to a system level (frame, C-Wall) assuming initially a 2D behaviour of the key structural system, and then incorporating a by-hand 3D behaviour (torsional effects).

Research papers, University of Canterbury Library

Unreinforced masonry churches in New Zealand, similarly to everywhere else in the word have proven to be highly vulnerable to earthquakes, because of their particular construction features. The Canterbury (New Zealand) earthquake sequence, 2010-2011 caused an invaluable loss of local architectural heritage and of churches, as regrettably, some of them were demolished instead of being repaired. It is critical for New Zealand to advance the data collection, research and understanding pertaining to the seismic performance and protection of church buildings, with the aim to:

Research papers, University of Canterbury Library

Existing unreinforced masonry (URM) buildings are often composed of traditional construction techniques, with poor connections between walls and diaphragms that results in poor performance when subjected to seismic actions. In these cases the application of the common equivalent static procedure is not applicable because it is not possible to assure “box like” behaviour of the structure. In such conditions the ultimate strength of the structure relies on the behaviour of the macro-elements that compose the deformation mechanisms of the whole structure. These macroelements are a single or combination of structural elements of the structure which are bonded one to each other. The Canterbury earthquake sequence was taken as a reference to estimate the most commonly occurring collapse mechanisms found in New Zealand URM buildings in order to define the most appropriate macroelements.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.

Research papers, University of Canterbury Library

This study examines the performance of nonlinear total-stress wave-propagation site response analysis for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015).

Research papers, University of Canterbury Library

The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.

Research papers, University of Canterbury Library

This paper presents preliminary results of an experimental campaign on three beam-column joint subassemblies extracted from a 22-storey reinforced concrete frame building constructed in late 1980s at the Christchurch’s Central Business District (CBD) area, damaged and demolished after the 2010-2011 Canterbury earthquakes sequence (CES). The building was designed following capacity design principles. Column sway (i.e., soft storey) mechanisms were avoided, and the beams were provided with plastic hinge relocation details at both beam-ends, aiming at developing plastic hinges away from the column faces. The specimens were tested under quasi-static cyclic displacement controlled lateral loading. One of the specimens, showing no visible residual cracks was cyclically tested in its as-is condition. The other two specimens which showed residual cracks varying between hairline and 1.0mm in width, were subjected to cyclic loading to simulate cracking patterns consistent with what can be considered moderate damage. The cracked specimens were then repaired with an epoxy injection technique and subsequently retested until reaching failure. The epoxy injection techniques demonstrated to be quite efficient in partly, although not fully, restoring the energy dissipation capacities of the damaged specimens at all beam rotation levels. The stiffness was partly restored within the elastic range and almost fully restored after the onset of nonlinear behaviour.

Research papers, University of Canterbury Library

he 2016 Building (Earthquake Prone Building) Amendment Act aims to improve the system for managing earthquake-prone buildings. The proposed changes to the Act were precipitated by the Canterbury earthquakes, and the need to improve the seismic safety of New Zealand’s building stock. However, the Act has significant ramifications for territorial authorities, organisations and individuals in small New Zealand towns, since assessing and repairing heritage buildings poses a major cost to districts with low populations and poor rental returns on commercial buildings.

Research papers, University of Canterbury Library

The south Leader Fault (SLF) is a newly documented active structure that ruptured the surface during the Mw 7.8 Kaikoura earthquake. The Leader Fault is a NNE trending oblique left lateral thrust that links the predominantly right lateral ‘The Humps’ and Conway-Charwell faults. The present research uses LiDAR at 0.5 m resolution and field mapping to determine the factors controlling the surface geometries and kinematics of the south Leader Fault ruptures at the ground surface. The SLF zone is up to 2km wide and comprises a series of echelon NE-striking thrusts linked by near-vertical N-S striking faults. The thrusts are upthrown to the west by up to 1 m and dip 35-45°. Thrust slip surfaces are parallel with Cretaceous-Cenozoic bedding and may reflect flexural slip folding. By contrast, the northerly striking faults dip steeply (65° west- 85° east), and accommodate up to 3m of oblique left lateral displacement at the ground surface and displace Cenozoic bedding. Some of the SLF has been mapped in bedrock, although none were known to be active prior to the earthquake or have a strong topographic expression. The complexity of fault rupture and the width of the fault zone appears to reflect the occurrence of faulting and folding at the ground surface during the earthquake.

Research papers, University of Canterbury Library

This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.

Research papers, University of Canterbury Library

The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).

Research papers, University of Canterbury Library

We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.

Research papers, University of Canterbury Library

We present preliminary observations on three waters impacts from the Mw7.8 14th November 2016 Kaikōura Earthquake on wider metropolitan Wellington, urban and rural Marlborough, and in Kaikōura township. Three waters systems in these areas experienced widespread and significant transient ground deformation in response to seismic shaking, with localised permanent ground deformation via liquefaction and lateral spreading. In Wellington, potable water quality was impacted temporarily by increased turbidity, and significant water losses occurred due to damaged pipes at the port. The Seaview and Porirua wastewater treatment plants sustained damage to clarifier tanks from water seiching, and increased water infiltration to the wastewater system occurred. Most failure modes in urban Marlborough were similar to the 2010-2011 Canterbury Earthquake Sequence; however some rural water tanks experienced rotational and translational movements, highlighting importance of flexible pipe connections. In Kaikōura, damage to reservoirs and pipes led to loss of water supply and compromised firefighting capability. Wastewater damage led to environmental contamination, and necessitated restrictions on greywater entry into the system to minimise flows. Damage to these systems necessitated the importation of tankered and bottled water, boil water notices and chlorination of the system, and importation of portaloos and chemical toilets. Stormwater infrastructure such as road drainage channels was also damaged, which could compromise condition of underlying road materials. Good operational asset management practices (current and accurate information, renewals, appreciation of criticality, good system knowledge and practical contingency plans) helped improve system resilience, and having robust emergency management centres and accurate Geographic Information System data allowed effective response coordination. Minimal damage to the wider built environment facilitated system inspections. Note Future research will include detailed geospatial assessments of seismic demand on these systems and attendant modes of failure, levels of service restoration, and collaborative development of resilience measures.

Research papers, University of Canterbury Library

Capacity design and hierarchy of strength philosophies at the base of modern seismic codes allow inelastic response in case of severe earthquakes and thus, in most traditional systems, damage develops at well-defined locations of reinforced concrete (RC) structures, known as plastic hinges. The 2010 and 2011 Christchurch earthquakes have demonstrated that this philosophy worked as expected. Plastic hinges formed in beams, in coupling beams and at the base of columns and walls. Structures were damaged permanently, but did not collapse. The 2010 and 2011 Christchurch earthquakes also highlighted a critical issue: the reparability of damaged buildings. No methodologies or techniques were available to estimate the level of subsequent earthquakes that RC buildings could still sustain before collapse. No repair techniques capable of restoring the initial condition of buildings were known. Finally, the cost-effectiveness of an eventual repair intervention, when compared with a new building, was unknown. These aspects, added to nuances of New Zealand building owners’ insurance coverage, encouraged the demolition of many buildings. Moreover, there was a perceived strong demand from government and industry to develop techniques for assessing damage to steel reinforcement bars embedded in cracked structural concrete elements. The most common questions were: “Have the steel bars been damaged in correspondence to the concrete cracks?”, “How much plastic deformation have the steel bars undergone?”, and “What is the residual strain capacity of the damaged bars?” Minimally invasive techniques capable of quantifying the level and extent of plastic deformation and residual strain capacity are not yet available. Although some studies had been recently conducted, a validated method is yet to be widely accepted. In this thesis, a least-invasive method for the damage-assessment of steel reinforcement is developed. Based on the information obtained from hardness testing and a single tensile test, it is possible to estimate the mechanical properties of earthquake-damaged rebars. The reduction in the low-cycle fatigue life due to strain ageing is also quantified. The proposed damage assessment methodology is based on empirical relationships between hardness and strain and residual strain capacity. If damage is suspected from in situ measurements, visual inspection or computer analysis, a bar may be removed and more accurate hardness measurements can be obtained using the lab-based Vickers hardness methodology. The Vickers hardness profile of damaged bars is then compared with calibration curves (Vickers hardness versus strain and residual strain capacity) previously developed for similar steel reinforcement bars extracted from undamaged locations. Experimental tests demonstrated that the time- and temperature-dependent strain-ageing phenomenon causes changes in the mechanical properties of plastically deformed steels. In particular, yield strength and hardness increases, whereas ductility decreases. The changes in mechanical properties are quantified and their implications on the hardness method are highlighted. Low-cycle fatigue (LCF) failures of steel reinforcing bars have been observed in laboratory testing and post-earthquake damage inspections. Often, failure might not occur during a first seismic event. However, damage is accumulated and the remaining fatigue life is reduced. Failure might therefore occur in a subsequent seismic event. Although numerous studies exist on the LCF behaviour of steel rebars, no studies had been conducted on the strain-ageing effects on the remaining fatigue life. In this thesis, the reduction in fatigue life due to this phenomenon is determined through a number of experimental tests.

Research papers, University of Canterbury Library

© 2017 The Royal Society of New Zealand. This paper discusses simulated ground motion intensity, and its underlying modelling assumptions, for great earthquakes on the Alpine Fault. The simulations utilise the latest understanding of wave propagation physics, kinematic earthquake rupture descriptions and the three-dimensional nature of the Earth's crust in the South Island of New Zealand. The effect of hypocentre location is explicitly examined, which is found to lead to significant differences in ground motion intensities (quantified in the form of peak ground velocity, PGV) over the northern half and southwest of the South Island. Comparison with previously adopted empirical ground motion models also illustrates that the simulations, which explicitly model rupture directivity and basin-generated surface waves, lead to notably larger PGV amplitudes than the empirical predictions in the northern half of the South Island and Canterbury. The simulations performed in this paper have been adopted, as one possible ground motion prediction, in the ‘Project AF8’ Civil Defence Emergency Management exercise scenario. The similarity of the modelled ground motion features with those observed in recent worldwide earthquakes as well as similar simulations in other regions, and the notably higher simulated amplitudes than those from empirical predictions, may warrant a re-examination of regional impact assessments for major Alpine Fault earthquakes.

Research papers, University of Canterbury Library

Damage distribution maps from strong earthquakes and recorded data from field experiments have repeatedly shown that the ground surface topography and subsurface stratigraphy play a decisive role in shaping the ground motion characteristics at a site. Published theoretical studies qualitatively agree with observations from past seismic events and experiments; quantitatively, however, they systematically underestimate the absolute level of topographic amplification up to an order of magnitude or more in some cases. We have hypothesized in previous work that this discrepancy stems from idealizations of the geometry, material properties, and incident motion characteristics that most theoretical studies make. In this study, we perform numerical simulations of seismic wave propagation in heterogeneous media with arbitrary ground surface geometry, and compare results with high quality field recordings from a site with strong surface topography. Our goal is to explore whether high-fidelity simulations and realistic numerical models can – contrary to theoretical models – capture quantitatively the frequency and amplitude characteristics of topographic effects. For validation, we use field data from a linear array of nine portable seismometers that we deployed on Mount Pleasant and Heathcote Valley, Christchurch, New Zealand, and we compute empirical standard spectral ratios (SSR) and single-station horizontal-to-vertical spectral ratios (HVSR). The instruments recorded ambient vibrations and remote earthquakes for a period of two months (March-April 2017). We next perform two-dimensional wave propagation simulations using the explicit finite difference code FLAC. We construct our numerical model using a high-resolution (8m) Digital Elevation Map (DEM) available for the site, an estimated subsurface stratigraphy consistent with the geomorphology of the site, and soil properties estimated from in-situ and non-destructive tests. We subject the model to in-plane and out-of-plane incident motions that span a broadband frequency range (0.1-20Hz). Numerical and empirical spectral ratios from our blind prediction are found in very good quantitative agreement for stations on the slope of Mount Pleasant and on the surface of Heathcote Valley, across a wide range of frequencies that reveal the role of topography, soil amplification and basin edge focusing on the distribution of ground surface motion.

Research papers, University of Canterbury Library

At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.

Research papers, University of Canterbury Library

At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.

Research papers, University of Canterbury Library

The Canterbury earthquakes of 2010 and 2011 caused significant damage and disruption to the city of Christchurch, New Zealand. A Royal Commission was established to report on the causes of building failure as a result of the earthquakes as well as look at the legal and best-practice requirements for buildings in New Zealand Central Business Districts. The Royal Commission made 189 recommendations on a variety of matters including managing damaged buildings after an earthquake, the adequacy of building codes and standards, and the processes of seismic assessments of existing buildings to determine their earthquake vulnerability. In response the Ministry of Business, Innovation and Employment, the agency responsible for administering building regulation in New Zealand, established a work programme to assist with the Canterbury rebuild and to implement the lessons learned throughout New Zealand. The five primary work streams in the programme are: • Facilitating the Canterbury Rebuild • Structural Performance and Design Standards • Geotechnical and structural guidance • Existing Building Resilience • Post Disaster Building Management This paper provides more detail on each of the work streams. There has been significant collaboration between the New Zealand Government and the research community, technical societies, and engineering consultants, both within New Zealand and internationally, to deliver the programme and improve the resilience of the New Zealand built environment. This has presented major challenges for an extremely busy industry in the aftermath of the Canterbury earthquakes. The paper identifies the items of work that have been completed and the work that is still in progress at the time of writing.

Research papers, University of Canterbury Library

Research in the governance of urban tourist spaces is characterized by a lack of argumentative inquiry and scant use of critical theory. This is evident, particularly, in the study of tourism and post-disaster urban recovery, with very few contributions assessing the phenomenon from a social theory perspective. This thesis examines the complex phenomenon of planning and governance for urban tourism spaces in contexts facing physical recovery from natural disasters. It does so by looking at the governance dynamics and the mechanism of decision- making put in place before and after triggering events like earthquakes and tsunamis. This thesis provides evidence from Christchurch, New Zealand, by focusing on the policies and strategies for the regeneration of the city centre put in place before and after the disruptive earthquakes of 2010 and 2011. The thesis looks at power relations, structures and ideologies through a Lukesian appraisal of pre-and-post disaster governance from two relevant urban tourist spaces located in the Christchurch central city area: the Arts Centre of Christchurch and the Town Hall and Performing Arts Precinct. The research strategy adopted for the study combined archival research, interviews with key stakeholders and fieldwork notes over a period of two years. The research deployed a comparative case study methodology that focuses on projects taking place within a spatially defined area of the city centre where special legislation was enacted as result of the earthquakes. The findings from the interviews and their triangulation with documents retrieved from national and local authorities suggest that the earthquakes affected the engagement among stakeholders and the mechanisms of decision-making. Also, the findings show patterns of disaster capitalism in post-earthquake governance for urban tourist spaces in the Christchurch CBD, with episodes of exclusion, lobbying and amendment of rules and legislation that directly benefited the interests of a narrow group of privileged stakeholders. Overall, the study shows that the earthquakes of 2010 and 2011 accelerated neoliberal practices of site development in Christchurch, with the seismic events used as a pretext to implement market-oriented site projects in the CBD area.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Between 2010 and 2011, Canterbury experienced a series of four large earthquake events with associated aftershocks which caused widespread damage to residential and commercial infrastructure. Fine grained and uncompacted alluvial soils, typical to the Canterbury outwash plains, were exposed to high peak ground acceleration (PGA) during these events. This rapid increase in PGA induced cyclic strain softening and liquefaction in the saturated, near surface alluvial soils. Extensive research into understanding the response of soils in Canterbury to dynamic loading has since occurred. The Earthquake Commission (EQC), the Ministry of Business and Employment (MBIE), and the Christchurch City Council (CCC) have quantified the potential hazards associated with future seismic events. Theses bodies have tested numerous ground improvement design methods, and subsequently are at the forefront of the Canterbury recovery and rebuild process. Deep Soil Mixing (DSM) has been proven as a viable ground improvement foundation method used to enhance in situ soils by increasing stiffness and positively altering in situ soil characteristics. However, current industry practice for confirming the effectiveness of the DSM method involves specific laboratory and absolute soil test methods associated with the mixed column element itself. Currently, the response of the soil around the columns to DSM installation is poorly understood. This research aims to understand and quantify the effects of DSM columns on near surface alluvial soils between the DSM columns though the implementation of standardised empirical soil test methods. These soil strength properties and ground improvement changes have been investigated using shear wave velocity (Vs), soil behaviour and density response methods. The results of the three different empirical tests indicated a consistent improvement within the ground around the DSM columns in sandier soils. By contrast, cohesive silty soils portrayed less of a consistent response to DSM, although still recorded increases. Generally, within the tests completed 50 mm from the column edge, the soil response indicated a deterioration to DSM. This is likely to be a result of the destruction of the soil fabric as the stress and strain of DSM is applied to the un‐mixed in situ soils. The results suggest that during the installation of DSM columns, a positive ground effect occurs in a similar way to other methods of ground improvement. However, further research, including additional testing following this empirical method, laboratory testing and finite 2D and 3D modelling, would be useful to quantify, in detail, how in situ soils respond and how practitioners should consider these test results in their designs. This thesis begins to evaluate how alluvial soils tend to respond to DSM. Conducting more testing on the research site, on other sites in Christchurch, and around the world, would provide a more complete data set to confirm the results of this research and enable further evaluation. Completing this additional research could help geotechnical DSM practitioners to use standardised empirical test methods to measure and confirm ground improvement rather than using existing test methods in future DSM projects. Further, demonstrating the effectiveness of empirical test methods in a DSM context is likely to enable more cost effective and efficient testing of DSM columns in future geotechnical projects.

Research papers, University of Canterbury Library

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.