Search

found 37 results

Research papers, University of Canterbury Library

This study investigates the uncertainty of simulated earthquake ground motions for smallmagnitude events (Mw 3.5 – 5) in Canterbury, New Zealand. 148 events were simulated with specified uncertainties in: event magnitude, hypocentre location, focal mechanism, high frequency rupture velocity, Brune stress parameter, the site 30-m time-averaged shear wave velocity (Vs30), anelastic attenuation (Q) and high frequency path duration. In order to capture these uncertainties, 25 realisations for each event were generated using the Graves and Pitarka (2015) hybrid broadband simulation approach. Monte-Carlo realisations were drawn from distributions for each uncertainty, to generate a suite of simulation realisations for each event and site. The fit of the multiple simulation realisations to observations were assessed using linear mixed effects regression to generate the systematic source, path and site effects components across all ground motion intensity measure residuals. Findings show that additional uncertainties are required in each of the three source, path, and site components, however the level of output uncertainty is promising considering the input uncertainties included.

Research papers, University of Canterbury Library

This study provides an initial examination of source parameter uncertainty in a New Zealand ground motion simulation model, by simulating multiple event realisations with perturbed source parameters. Small magnitude events in Canterbury have been selected for this study due to the small number of source input parameters, the wealth of recorded data, and the lack of appreciable off-fault non-linear effects. Which provides greater opportunity to identify systematic source, path and site effects, required to robustly investigate the causes of uncertainty.

Research papers, University of Canterbury Library

Seismically vulnerable buildings constitute a major problem for the safety of human beings. In many parts of the world, reinforced concrete (RC) frame buildings designed and constructed with substandard detailing, no consideration of capacity design principles, and improper or no inclusion of the seismic actions, have been identified. Amongst those vulnerable building, one particular typology representative of the construction practice of the years previous to the 1970’s, that most likely represents the worst case scenario, has been widely investigated in the past. The deficiencies of that building typology are related to non-ductile detailing in beam column joints such as the use of plain round bars, the lack of stirrups inside the joint around the longitudinal reinforcement of the column, the use of 180° end hooks in the beams, the use of lap splices in potential ‘plastic hinge’ regions, and substandard quality of the materials. That type of detailing and the lack of a capacity design philosophy create a very fragile fuse in the structure where brittle inelastic behaviour is expected to occur, which is the panel zone region of exterior beam column joints. The non-ductile typology described above was extensively investigated at the University of Canterbury in the context of the project ‘Retrofit Solutions for New Zealand Multi-Storey Buildings’ (2004-2011), founded by the ‘Foundation for Research, Science and Technology’ Tūāpapa Rangahau Pūtaiao. The experimental campaign prior to the research carried out by the author consisted of quasi-static tests of beam column joint subassemblies subjected to lateral loading regime, with constant and varying axial load in the column. Most of those specimens were representative of a plane 2D frame (knee joint), while others represented a portion of a space 3D frame (corner joints), and only few of them had a floor slab, transverse beams, and lap splices. Using those experiments, several feasible, cost-effective, and non-invasive retrofit techniques were developed, improved, and refined. Nevertheless, the slow motion nature of those experiments did not take into account the dynamical component inherent to earthquake related problems. Amongst the set of techniques investigated, the use of FRP layers for strengthening beam column joints is of particular interest due to its versatility and the momentum that its use has gained in the current state of the practice. That particular retrofit technique was previously used to develop a strengthening scheme suitable for plane 2D and space 3D corner beam column joints, but lacking of floor slabs. In addition, a similar scheme was not developed for exterior joints of internal frames, referred here as ‘cruciform’. In this research a 2/5 scale RC frame model building comprising of two frames in parallel (external and internal) joined together by means of floor slabs and transverse beams, with non-ductile characteristics identical to those of the specimens investigated previously by others, and also including lap splices, was developed. In order to investigate the dynamic response of that building, a series of shake table tests with different ground motions were performed. After the first series of tests, the specimen was modified by connecting the spliced reinforcement in the columns in order to capture a different failure mode. Ground motions recorded during seismic events that occurred during the initial period of the experimental campaign (2010) were used in the subsequent experiments. The hierarchy of strengths and sequence of events in the panel zone region were evaluated in an extended version of the bending moment-axial load (M-N) performance domain developed by others. That extension was required due to the asymmetry in the beam cross section introduced by the floor slab. In addition, the effect of the torsion resistance provided by the spandrel (transverse beam) was included. In order to upgrade the brittle and unstable performance of the as-built/repaired specimen, a practical and suitable ad-hoc FRP retrofit intervention was developed, following a partial retrofit strategy that aimed to strengthen exterior beam column joints only (corner and cruciform). The ability of the new FRP scheme to revert the sequence of events in the panel zone region was evaluated using the extended version of the M-N performance domain as well as the guidelines for strengthening plane joints developed by others. Weakening of the floor slab in a novel configuration was also incorporated with the purpose of reducing the flexural capacity of the beam under negative bending moment (slab in tension), enabling the damage relocation from the joint into the beam. The efficacy of the developed retrofit intervention in upgrading the seismic performance of the as-built specimen was investigated using shake table tests with the input motions used in the experiments of the as-built/repaired specimen. Numerical work aimed to predict the response of the model building during the most relevant shake table tests was carried out. By using a simple numerical model with concentrated plasticity elements constructed in Ruaumoko2D, the results of blind and post-experimental predictions of the response of the specimen were addressed. Differences in the predicted response of the building using the nominal and the actual recorded motions of the shake table were investigated. The dependence of the accuracy of the numerical predictions on the assumed values of the parameters that control the hysteresis rules of key structural members was reviewed. During the execution of the experimental campaign part of this thesis, two major earthquakes affected the central part of Chile (27 of February 2010 Maule earthquake) and the Canterbury region in New Zealand (22 February 2011 Canterbury earthquake), respectively. As the author had the opportunity to experience those events and investigate their consequences in structures, the observations related to non-ductile detailing and drawbacks in the state of the practice related to reinforced concrete walls was also addressed in this research, resulting in preliminary recommendations for the refinement of current seismic code provisions and assessment guidelines. The investigations of the ground motions recorded during those and other earthquakes were used to review the procedures related to the input motions used for nonlinear dynamic analysis of buildings as required by most of the current code provisions. Inelastic displacement spectra were constructed using ground motions recorded during the earthquakes mentioned above, in order to investigate the adequacy of modification factors used to obtain reduced design spectra from elastic counterparts. Finally a simplified assessment procedure for RC walls that incorporates capacity compatible spectral demands is proposed.

Research papers, University of Canterbury Library

Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.

Research papers, University of Canterbury Library

Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.

Research papers, University of Canterbury Library

The performance of conventionally designed reinforced concrete (RC) structures during the 2011 Christchurch earthquake has demonstrated that there is greater uncertainty in the seismic performance of RC components than previously understood. RC frame and wall structures in the Christchurch central business district were observed to form undesirable cracks patterns in the plastic hinge region while yield penetration either side of cracks, and into development zones, were less than theoretical predictions. The implications of this unexpected behaviour: (i) significantly less available ductility; (ii) less hysteretic energy dissipation; and (iii) the localization of peak reinforcement strains, results in considerable doubt for the residual capacity of RC structures. The significance of these consequences has prompted a review of potential sources of uncertainty in seismic experimentation with the intention to improve the current confidence level for newly designed conventional RC structures. This paper attempts to revisit the principles of RC mechanics, in particular, to consider the influence of loading history, concrete tensile strength, and reinforcement ratio on the performance of ‘real’ RC structures compared to experimental test specimens.

Research Papers, Lincoln University

The world is constantly changing. Christchurch, New Zealand, has recently experienced drastic changes after earthquakes struck the city. The earthquakes caused the city to physically shake, and the land to sink in some places and rise in others. Now further change is forcast and parts of Christchurch could be under water by 2115 according to experts. Climate change induced sea level rise is recognised as a international issue with potential impacts for coastal communities all over the world. The Chrischtchurch City Council is required to have a 100-year planning horizon for sea level rise and this means planning for at least one meter, and possibly up to two meters, of sea level rise by 2115. This dissertation investigates the planning response to slow onset disasters, change, and uncertainty, using the example of sea-level rise in Christchurch, and it examines the role of public participation in this. To achieve this, the ways in which planning theory and practice acknowledge uncertainty, and cope with change, were critically analysed along with the Christchurch City Council’s response to the Tonkin and Taylor predictions and modelling. Semi-structured interviews with professionals in natural hazards risk reduction, policy, and planning were conducted, and the previous and proposed Christchurch City District Plans were compared. Planning for sea level rise in Christchurch provides an example of how planners may cope with slow onset change. The results of this dissertation suggests that the favoured risk reduction strategy for coastal communities in Christchurch is an adaptation strategy, and at present there is no sign of managed retreat being employed. The results also suggests using a planning approach that involves public participation for best results when planning for change, uncertainty or slow onset disasters.

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

One of the great challenges facing human systems today is how to prepare for, manage, and adapt successfully to the profound and rapid changes wreaked by disasters. Wellington, New Zealand, is a capital city at significant risk of devastating earthquake and tsunami, potentially requiring mass evacuations with little or short notice. Subsequent hardship and suffering due to widespread property damage and infrastructure failure could cause large areas of the Wellington Region to become uninhabitable for weeks to months. Previous research has shown that positive health and well-being are associated with disaster-resilient outcomes. Preventing adverse outcomes before disaster strikes, through developing strengths-based skill sets in health-protective attitudes and behaviours, is increasingly advocated in disaster research, practise, and management. This study hypothesised that well-being constructs involving an affective heuristic play vital roles in pathways to resilience as proximal determinants of health-protective behaviours. Specifically, this study examined the importance of health-related quality of life and subjective well-being in motivating evacuation preparedness, measured in a community sample (n=695) drawn from the general adult population of Wellington’s isolated eastern suburbs. Using a quantitative epidemiological approach, the study measured the prevalence of key quality of life indicators (physical and mental health, emotional well-being or “Sense of Coherence”, spiritual well-being, social well-being, and life satisfaction) using validated psychometric scales; analysed the strengths of association between these indicators and the level of evacuation preparedness at categorical and continuous levels of measurement; and tested the predictive power of the model to explain the variance in evacuation preparedness activity. This is the first study known to examine multi-dimensional positive health and global well-being as resilient processes for engaging in evacuation preparedness behaviour. A cross-sectional study design and quantitative survey were used to collect self-report data on the study variables; a postal questionnaire was fielded between November 2008 and March 2009 to a sampling frame developed through multi-stage cluster randomisation. The survey response rate was 28.5%, yielding a margin of error of +/- 3.8% with 95% confidence and 80% statistical power to detect a true correlation coefficient of 0.11 or greater. In addition to the primary study variables, data were collected on demographic and ancillary variables relating to contextual factors in the physical environment (risk perception of physical and personal vulnerability to disaster) and the social environment (through the construct of self-determination), and other measures of disaster preparedness. These data are reserved for future analyses. Results of correlational and regression analyses for the primary study variables show that Wellingtonians are highly individualistic in how their well-being influences their preparedness, and a majority are taking inadequate action to build their resilience to future disaster from earthquake- or tsunami-triggered evacuation. At a population level, the conceptual multi-dimensional model of health-related quality of life and global well-being tested in this study shows a positive association with evacuation preparedness at statistically significant levels. However, it must be emphasised that the strength of this relationship is weak, accounting for only 5-7% of the variability in evacuation preparedness. No single dimension of health-related quality of life or well-being stands out as a strong predictor of preparedness. The strongest associations for preparedness are in a positive direction for spiritual well-being, emotional well-being, and life satisfaction; all involve a sense of existential meaningfulness. Spiritual well-being is the only quality of life variable making a statistically significant unique contribution to explaining the variance observed in the regression models. Physical health status is weakly associated with preparedness in a negative direction at a continuous level of measurement. No association was found at statistically significant levels for mental health status and social well-being. These findings indicate that engaging in evacuation preparedness is a very complex, holistic, yet individualised decision-making process, and likely involves highly subjective considerations for what is personally relevant. Gender is not a factor. Those 18-24 years of age are least likely to prepare and evacuation preparedness increases with age. Multidimensional health and global well-being are important constructs to consider in disaster resilience for both pre-event and post-event timeframes. This work indicates a need for promoting self-management of risk and building resilience by incorporating a sense of personal meaning and importance into preparedness actions, and for future research into further understanding preparedness motivations.

Audio, Radio New Zealand

Central Christchurch restaurant and bars say they could be heading into the "worst winter to date". Eight years on from the earthquakes there are more restaurants and bars in the city than ever before - but owners say there aren't enough customers. Now they're grappling with added uncertainty over the effect of the mosque attacks on visitor numbers.

Audio, Radio New Zealand

In quake-ravaged Christchurch businesses are tentatively restarting, and infrastructure is being restored, but there's ongoing uncertainty about job losses and how people will survive financially. Within six days of the February earthquake; the Government had introduced a subsidy scheme for businesses and their employees, as well as people who'd found themselves out of a job. Now there's growing disquiet about what will happen to thousands of Cantabrians when that support scheme finishes at the end of May.

Research papers, University of Canterbury Library

This paper investigates the effects of variability in source rupture parameters on site-specific physics-based simulated ground motions, ascertained through the systematic analysis of ground motion intensity measures. As a preliminary study, we consider simulations of the 22 February 2011 Christchurch earthquake using the Graves and Pitarka (2015) methodology. The effects of source variability are considered via a sensitivity study in which parameters (hypocentre location, earthquake magnitude, average rupture velocity, fault geometry and the Brune stress parameter) are individually varied by one standard deviation. The sensitivity of simulated ground motion intensity measures are subsequently compared against observational data. The preliminary results from this study indicate that uncertainty in the stress parameter and the rupture velocity have the most significant effect on the high frequency amplitudes. Conversely, magnitude uncertainty was found to be most influential on the spectral acceleration amplitudes at low frequencies. Further work is required to extend this preliminary study to exhaustively consider more events and to include parameter covariance. The ultimate results of this research will assist in the validation of the overall simulation method’s accuracy in capturing various rupture parameters, which is essential for the use of simulated ground motion models in probabilistic seismic hazard analysis.

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.

Research papers, University of Canterbury Library

This report to RCP Ltd and University of Canterbury summarises the findings of a 5 month secondment to the CERA Port Hills Land Clearance Team. Improvement strategies were initiated and observed. The Port Hills Land Clearance Programme is the undertaking of the demolition of all built structures from the Crown’s compulsory acquired 714 residential red zoned properties. These properties are zoned red due to an elevated life risk as a result of geotechnical land uncertainty following the 2011 Canterbury Earthquakes.

Research Papers, Lincoln University

Artificial Neural Networks (ANN) as a tool offers opportunities for modeling the inherent complexity and uncertainty associated with socio-environmental systems. This study draws on New Zealand ski fields (multiple locations) as socio- environmental systems while considering their perceived resilience to low probability but potential high consequences catastrophic natural events (specifically earthquakes). We gathered data at several ski fields using a mixed methodology including: geomorphic assessment, qualitative interviews, and an adaptation of Ozesmi and Ozesmi’s (2003) multi-step fuzzy cognitive mapping (FCM) approach. The data gathered from FCM are qualitatively condensed, and aggregated to three different participant social groups. The social groups include ski fields users, ski industry workers, and ski field managers. Both quantitative and qualitative indices are used to analyze social cognitive maps to identify critical nodes for ANN simulations. The simulations experiment with auto-associative neural networks for developing adaptive preparation, response and recovery strategies. Moreover, simulations attempt to identify key priorities for preparation, response, and recovery for improving resilience to earthquakes in these complex and dynamic environments. The novel mixed methodology is presented as a means of linking physical and social sciences in high complexity, high uncertainty socio-environmental systems. Simulation results indicate that participants perceived that increases in Social Preparation Action, Social Preparation Resources, Social Response Action and Social Response Resources have a positive benefit in improving the resilience to earthquakes of ski fields’ stakeholders.

Research papers, University of Canterbury Library

Using case studies from the 2010-2011 Canterbury, New Zealand earthquake sequence, this study assesses the accuracies of paleoliquefaction back-analysis methods and explores the challenges, techniques, and uncertainties associated with their application. While liquefaction-based back-analyses have been widely used to estimate the magnitudes of paleoearthquakes, their uncertain efficacies continue to significantly affect the computed seismic hazard in regions where they are relied upon. Accordingly, their performance is evaluated herein using liquefaction data from modern earthquakes with known magnitudes. It is shown that when the earthquake source location and mechanism are known, back-analysis methods are capable of accurately deriving seismic parameters from liquefaction evidence. However, because the source location and mechanism are often unknown in paleoseismic studies, and because accurate interpretation is shown to be more difficult in such cases, new analysis techniques are proposed herein. An objective parameter is proposed to geospatially assess the likelihood of any provisional source location, enabling an analyst to more accurately estimate the magnitude of a liquefaction-inducing paleoearthquake. This study demonstrates the application of back-analysis methods, provides insight into their potential accuracies, and provides a framework for performing paleoliquefaction analyses worldwide.

Research Papers, Lincoln University

The disastrous earthquakes that struck Christchurch in 2010 and 2011 seriously impacted on the individual and collective lives of Māori residents. This paper continues earlier, predominantly qualitative research on the immediate effects on Māori by presenting an analysis of a survey carried out 18 months after the most destructive event, on 22 February 2011. Using a set-theoretic approach, pathways to Māori resilience are identified, emphasising the combination of whānau connectivity and high incomes in those who have maintained or increased their wellbeing post-disaster. However, the results show that if resilience is used to describe a “bounce back” in wellbeing, Māori are primarily enduring the post-disaster environment. This endurance phase is a precursor to any resilience and will be of much longer duration than first thought. With continued uncertainty in the city and wider New Zealand economy, this endurance may not necessarily lead to a more secure environment for Māori in the city.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

The 2010 Darfield and 2011 Christchurch Earthquakes triggered extensive liquefaction-induced lateral spreading proximate to streams and rivers in the Christchurch area, causing significant damage to structures and lifelines. A case study in central Christchurch is presented and compares field observations with predicted displacements from the widely adopted empirical model of Youd et al. (2002). Cone penetration testing (CPT), with measured soil gradation indices (fines content and median grain size) on typical fluvial deposits along the Avon River were used to determine the required geotechnical parameters for the model input. The method presented attempts to enable the adoption of the extensive post-quake CPT test records in place of the lower quality and less available Standard Penetration Test (SPT) data required by the original Youd model. The results indicate some agreement between the Youd model predictions and the field observations, while the majority of computed displacements error on the side of over-prediction by more than a factor of two. A sensitivity analysis was performed with respect to the uncertainties used as model input, illustrating the model’s high sensitivity to the input parameters, with median grain size and fines content among the most influential, and suggesting that the use of CPT data to quantify these parameters may lead to variable results.

Articles, UC QuakeStudies

The previously unknown Greendale Fault ruptured to the ground surface, causing up to 5 metres horizontal and 1 metre vertical permanent offset of the ground, during the September 2010 Darfield (Canterbury) earthquake. Environment Canterbury commissioned GNS Science, with help from the University of Canterbury, to define a fault avoidance zone and to estimate the fault recurrence interval. There is little evidence for past movement on the fault in the past 16,000 years. However, because of the uncertainties involved, a conservative approach was taken and the fault has been categorised as a Recurrence Interval Class IV fault (a recurrence interval of between 5,000 and 10,000 years). A PhD study by a University of Canterbury student will work towards refining the Recurrence Interval Class over the next three years. Taking a risk-based approach, the Ministry for the Environment Active Fault Guidelines recommend that normal residential development be allowed within the fault avoidance zone for faults of this Recurrence Interval Class, but recommends restrictions for larger community buildings or facilities with post-disaster functions. The report is assisting Selwyn District Council in granting consents for rebuilding houses on or near the Greendale Fault that were damaged by permanent distortion of the ground due to the fault rupture in the September 2010 earthquake. The report provides specific recommendations for building on or close to the Greendale Fault, which are being implemented by Selwyn District Council. See Object Overview for background and usage information.

Research Papers, Lincoln University

An emerging water crisis is on the horizon and is poised to converge with several other impending problems in the 21st century. Future uncertainties such as climate change, peak oil and peak water are shifting the international focus from a business as usual approach to an emphasis on sustainable and resilient strategies that better meet these challenges. Cities are being reimagined in new ways that take a multidisciplinary approach, decompartmentalising functions and exploring ways in which urban systems can share resources and operate more like natural organisms. This study tested the landscape design implications of wastewater wetlands in the urban environment and evaluated their contribution to environmental sustainability, urban resilience and social development. Black and grey water streams were the central focus of this study and two types of wastewater wetlands, tidal flow (staged planning) and horizontal subsurface flow wetlands were tested through design investigations in the earthquake-affected city of Christchurch, New Zealand. These investigations found that the large area requirements of wastewater wetlands can be mitigated through landscape designs that enhance a matrix of open spaces and corridors in the city. Wastewater wetlands when combined with other urban and rural services such as food production, energy generation and irrigation can aid in making communities more resilient. Landscape theory suggests that the design of wastewater wetlands must meet cultural thresholds of beauty and that the inclusion of waste and ecologies in creatively designed landscapes can deepen our emotional connection to nature and ourselves.

Research Papers, Lincoln University

"Prior to the devastating 2010-2011 Canterbury earthquakes, the city of Christchurch was already exhibiting signs of a housing affordability crisis. The causes and symptoms were similar to those being experienced in Auckland, but the substantial damage to the housing stock caused by the earthquakes added new dimensions and impetus to the problem. Large swathes of the most affordable housing stock in the east of the city were effectively destroyed by the earthquakes. In itself this would have pushed the mean house price upwards, but compounding problems exacerbated the situation. These include the price effects of reduced supply of both rented and owned housing and increased demand from both displaced residents and an influx of rebuild workers. The need for additional temporary housing while repairs were undertaken and the associated insurance pay-outs bidding up rents with improved rental returns leading to increased interest in property investment. Land supply constraints and consenting issues inhibiting the build of new housing and political infighting and uncertainty regarding the future of parts of the city leading to a flight of development activity to peripheral locations and adjoining local authorities. Concerns that the erosion of the city council rating base combined with inadequacy of insurance cover for infrastructure will lead to large rates increases, increased development costs and reduced amenities and services in future years. These and other issuers will be elaborated on in this paper with a view to exploring the way forward for affordable housing Christchurch City."

Research papers, Victoria University of Wellington

We examine the role of business interruption insurance in business recovery following the Christchurch earthquake in 2011 in the short- and medium-term. In the short-term analysis, we ask whether insurance increases the likelihood of business survival in the aftermath of a disaster. We find only weak evidence that those firms that had incurred damage, but were covered by business interruption insurance, had higher likelihood of survival post-quake compared with those firms that did not have insurance. This absence of evidence may reflect the high degree of uncertainty in the months following the 2011 earthquake and the multiplicity of severe aftershocks. For the medium-term, our results show a more explicit role for insurance in the aftermath of a disaster. Firms with business interruption insurance have a higher probability of increasing productivity and improved performance following a catastrophe. Furthermore, our results show that those organisations that receive prompt and full payments of their claims have a better recovery, in terms of profitability and a subjective ‘”better off” measure’ than those that had protracted or inadequate claim payments (less than 80% of the claim paid within 2.5 years). Interestingly, the latter group does worse than those organisations that had damage but no insurance coverage. This analysis strongly indicates the importance not only of good insurance coverage, but of an insurance system that also delivers prompt claim payments. As a first paper attempting to empirically identify a causal effect of insurance on business recovery, we also emphasize some caveats to our analysis.

Research papers, University of Canterbury Library

The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.

Research papers, University of Canterbury Library

This article presents a quantitative case study on the site amplification effect observed at Heathcote Valley, New Zealand, during the 2010-2011 Canterbury earthquake sequence for 10 events that produced notable ground acceleration amplitudes up to 1.4g and 2.2g in the horizontal and vertical directions, respectively. We performed finite element analyses of the dynamic response of the valley, accounting for the realistic basin geometry and the soil non-linear response. The site-specific simulations performed significantly better than both empirical ground motion models and physics based regional-scale ground motion simulations (which empirically accounts for the site effects), reducing the spectral acceleration prediction bias by a factor of two in short vibration periods. However, our validation exercise demonstrated that it was necessary to quantify the level of uncertainty in the estimated bedrock motion using multiple recorded events, to understand how much the simplistic model can over- or under-estimate the ground motion intensities. Inferences from the analyses suggest that the Rayleigh waves generated near the basin edge contributed significantly to the observed high frequency (f>3Hz) amplification, in addition to the amplification caused by the strong soil-rock impedance contrast at the site fundamental frequency. Models with and without considering soil non-linear response illustrate, as expected, that the linear elastic assumption severely overestimates ground motions in high frequencies for strong earthquakes, especially when the contribution of basin edge-generated Rayleigh waves becomes significant. Our analyses also demonstrate that the effect of pressure-dependent soil velocities on the high frequency ground motions is as significant as the amplification caused by the basin edge-generated Rayleigh waves.

Research Papers, Lincoln University

Environmental assessment in New Zealand is governed by the provisions of the Resource Management Act (RMA) 1991. The Act requires persons wishing to undertake certain activities to apply for resource consent from their local or regional council - a procedure termed the Resource Consent Process. The key component of a resource consent application is an Assessment of Environmental Effects (AEE) report; a statement of the environmental effects of a proposed activity. Problems arise when environmental assessments are complicated by uncertain and abnormal circumstances such as natural hazards. Natural hazards (including earthquakes, floods, tsunami, and coastal erosion) can be catastrophic to an environment. If hazards are not avoided or successfully mitigated, they can result in serious consequences to proposed development and to the environment which the proposal relates. The aim of this study is to assess the adequacy of the resource consent process (as outlined in the Resource Management Act 1991) for dealing with proposed development affected by natural hazards. This study reviews the context of the resource consent process for dealing with natural hazards to identify potential issues in the assessment process. Guidance criteria for assessing natural hazards (termed Natural Hazard Assessment) are developed to evaluate against two resource consent applications affected by natural hazards. The findings of the consent process review and case study evaluation are discussed to determine the adequacy of the consent process for dealing with natural hazards. From the review of the consent process it was evident that the process has a number of problems for accommodating natural hazards into the assessment. Although many important traits are provided for in the process, such traits are not always reflected in environmental assessments. Evaluation of two resource consent applications against the process of Natural Hazard Assessment (NHA) showed that these consent applications did not adequately detail key information relating to natural hazards. Many problems evident in these applications were not amended by the Consent Authorities in the review process and subsequently consent was granted to information-deficient applications. Problematic issues identified in this study include: • A distinct lack of guidance (legal or otherwise) for the applicant and Consent Authority regarding the boundaries of inclusion of an effect; • Deficiencies in planning documents are reflected in AEE reports, the review of the consent application and in the end-decision; • Under-utilisation of "experts" throughout the consent process; • Minimal identification and account for the degree of uncertainty throughout the consent process; • Resource consents are being granted even though information in consent applications, and the means for assessing the information is deficient. These issues reflect that decisions are not being made based on all elements involved in a potential hazard. Subsequently, the resource consent process is not adequate for dealing with all aspects of natural hazards. The Natural Hazard Assessment process provides educated assessment criteria to assess development affected by natural hazards. By accounting for the problems evident in the consent process, the introduction of a three-tier identification, risk and vulnerability assessment, and evaluation process to account for uncertainties, Natural Hazard Assessment provides a platform for a thorough assessment of natural hazards. The application of the principles of Natural Hazard Assessment to the consent applications affected by natural hazards showed that many key issues were not covered in the assessment under the consent process. The nature of a natural event is that one may not occur in a given region over many lifetimes, however they will occur at some stage and planning and environmental assessment needs to provide for the associated hazards. Implementation of Natural Hazard Assessment is needed to help provide answers for the problems experienced in the resource consent process. Natural Hazard Assessment would allow decision-makers to make informed judgements on the situation at hand, leading to better planning and land-use options. Change to current practice is needed, as following the current path of environmental assessment will be the hazard in the end.

Research papers, University of Canterbury Library

Previous earthquakes demonstrated destructive effects of soil-structure interaction on structural response. For example, in the 1970 Gediz earthquake in Turkey, part of a factory was demolished in a town 135 km from the epicentre, while no other buildings in the town were damaged. Subsequent investigations revealed that the fundamental period of vibration of the factory was approximately equal to that of the underlying soil. This alignment provided a resonance effect and led to collapse of the structure. Another dramatic example took place in Adapazari, during the 1999 Kocaeli earthquake where several foundations failed due to either bearing capacity exceedance or foundation uplifting, consequently, damaging the structure. Finally, the Christchurch 2012 earthquakes have shown that significant nonlinear action in the soil and soil-foundation interface can be expected due to high levels of seismic excitation and spectral acceleration. This nonlinearity, in turn, significantly influenced the response of the structure interacting with the soil-foundation underneath. Extensive research over more than 35 years has focused on the subject of seismic soil-structure interaction. However, since the response of soil-structure systems to seismic forces is extremely complex, burdened by uncertainties in system parameters and variability in ground motions, the role of soil-structure interaction on the structural response is still controversial. Conventional design procedures suggest that soil-structure interaction effects on the structural response can be conservatively ignored. However, more recent studies show that soil-structure interaction can be either beneficial or detrimental, depending on the soil-structure-earthquake scenarios considered. In view of the above mentioned issues, this research aims to utilise a comprehensive and systematic probabilistic methodology, as the most rational way, to quantify the effects of soil-structure interaction on the structural response considering both aleatory and epistemic uncertainties. The goal is achieved by examining the response of established rheological single-degree-of-freedom systems located on shallow-foundation and excited by ground motions with different spectral characteristics. In this regard, four main phases are followed. First, the effects of seismic soil-structure interaction on the response of structures with linear behaviour are investigated using a robust stochastic approach. Herein, the soil-foundation interface is modelled by an equivalent linear cone model. This phase is mainly considered to examine the influence of soil-structure interaction on the approach that has been adopted in the building codes for developing design spectrum and defining the seismic forces acting on the structure. Second, the effects of structural nonlinearity on the role of soil-structure interaction in modifying seismic structural response are studied. The same stochastic approach as phase 1 is followed, while three different types of structural force-deflection behaviour are examined. Third, a systematic fashion is carried out to look for any possible correlation between soil, structural, and system parameters and the degree of soil-structure interaction effects on the structural response. An attempt is made to identify the key parameters whose variation significantly affects the structural response. In addition, it is tried to define the critical range of variation of parameters of consequent. Finally, the impact of soil-foundation interface nonlinearity on the soil-structure interaction analysis is examined. In this regard, a newly developed macro-element covering both material and geometrical soil-foundation interface nonlinearity is implemented in a finite-element program Raumoko 3D. This model is then used in an extensive probabilistic simulation to compare the effects of linear and nonlinear soil-structure interaction on the structural response. This research is concluded by reviewing the current design guidelines incorporating soil-structure interaction effects in their design procedures. A discussion is then followed on the inadequacies of current procedures based on the outcomes of this study.