Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
Extended Direct Analysis (EDA), developed at the University of Canterbury, is an advance on the AISC Direct Analysis method for the analysis of frames subjected to static forces. EDA provides a faster, simple and more rational way to properly consider the second-order effects, initial residual stresses (IRS) and the initial imperfections or steel structures under one directional loading than conventional analysis methods. This research applied the EDA method to quantify the effect of member overstrength on frame behaviour for a single storey frame. Also, the effects of IRS, which were included in the EDA static analysis, but which are not considered explicitly in non-linear seismic analysis, were evaluated in two ways. Firstly, they were considered for simple structures subject to increasing cyclic displacement in different directions. Secondly, incremental dynamic analysis with realistic ground motion was used to quantify the likely effect of IRS in earthquakes. It was found that, contrary to traditional wisdom and practice, greater member strengths can result in lower frame strengths for frames under monotonic lateral loading. The structural lateral capacity of the overstrength case was reduced by 6% compared to the case using the dependable member strengths. Also, it resulted significantly different in member demands. Therefore, it is recommended that when either plastic analysis or EDA is used, that both upper and lower bounds on the likely member strength should be considered to determine the total frame strength and the member demands. Results of push-pull analysis under displacement control showed that for IRS ratio, gamma < 0.5 and axial compressive force ratio, N*/Ns, up to 0.5, IRS did affect the structural behaviour in the first half cycle. However, the behavior in the later cycles was not significantly affected. It also showed that the effect of initial residual stresses in the frame was less significant than for the column alone when the column was subjected to similar axial compressive force. The incremental dynamic analysis results from both cantilever column and the three-storey steel frame showed that by increasing gamma = 0 to 0.5, the effect of IRS on seismic responses, based on the 50% confidence level, was less than 3% for N*/Ns, up to 0.5.
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.
Despite over a century of study, the relationship between lunar cycles and earthquakes remains controversial and difficult to quantitatively investigate. Perhaps as a consequence, major earthquakes around the globe are frequently followed by 'prediction' claims, using lunar cycles, that generate media furore and pressure scientists to provide resolute answers. The 2010-2011 Canterbury earthquakes in New Zealand were no exception; significant media attention was given to lunarderived earthquake predictions by non-scientists, even though the predictions were merely 'opinions' and were not based on any statistically robust temporal or causal relationships. This thesis provides a framework for studying lunisolar earthquake temporal relationships by developing replicable statistical methodology based on peer reviewed literature. Notable in the methodology is a high accuracy ephemeris, called ECLPSE, designed specifically by the author for use on earthquake catalogs, and a model for performing phase angle analysis. The statistical tests were carried out on two 'declustered' seismic catalogs, one containing the aftershocks from the Mw7.1 earthquake in Canterbury, and the other containing Australian seismicity from the past two decades. Australia is an intraplate setting far removed from active plate boundaries and Canterbury is proximal to a plate boundary, thus allowing for comparison based on tectonic regime and corresponding tectonic loading rate. No strong, conclusive, statistical correlations were found at any level of the earthquake catalogs, looking at large events, onshore events, offshore events, and the fault type of some events. This was concluded using Schuster's test of significance with α=5% and analysis of standard deviations. A few weak correlations, with p-5-10% of rejecting the null hypothesis, and anomalous standard deviations were found, but these are difficult to interpret. The results invalidate the statistical robustness of 'earthquake predictions' using lunisolar parameters in this instance. An ambitious researcher could improve on the quality of the results and on the range of parameters analyzed. The conclusions of the thesis raise more questions than answers, but the thesis provides an adaptable methodology that can be used to further investigation the problem.
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.
Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.
© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.
The city of Christchurch and its surrounds experienced widespread damage due to soil liquefaction induced by seismic shaking during the Canterbury earthquake sequence that began in September 2010 with the Mw7.1 Darfield earthquake. Prior to the start of this sequence, the city had a large network of strong motion stations (SMSs) installed, which were able to record a vast database of strong ground motions. This paper uses this database of strong ground motion recordings, observations of liquefaction manifestation at the ground surface, and data from a recently completed extensive geotechnical site investigation program at each SMS to assess a range of liquefaction evaluation procedures at the four SMSs in the Christchurch Central Business District (CBD). In general, the characteristics of the accelerograms recorded at each SMS correlated well with the liquefaction evaluation procedures, with low liquefaction factors of safety predicted at sites with clear liquefaction identifiers in the ground motions. However, at sites that likely liquefied at depth (as indicated by evaluation procedures and/or inferred from the characteristics of the recorded surface accelerograms), the presence of a non-liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects. Because of this, there was not a good correlation between surface manifestation and two surface manifestation indices, the Liquefaction Potential Index (LPI) and the Liquefaction Severity Number (LSN).
There is a growing awareness of the need for the earthquake engineering practice to incorporate in addition to empirical approaches in evaluation of liquefaction hazards advanced methods which can more realistically represent soil behaviour during earthquakes. Currently, this implementation is hindered by a number of challenges mainly associated with the amount of data and user-experience required for such advanced methods. In this study, we present key steps of an advanced seismic effective-stress analysis procedure, which on the one hand can be fully automated and, on the other hand, requires no additional input (at least for preliminary applications) compared to simplified cone penetration test (CPT)-based liquefaction procedures. In this way, effective-stress analysis can be routinely applied for quick, yet more robust estimations of liquefaction hazards, in a similar fashion to the simplified procedures. Important insights regarding the dynamic interactions in liquefying soils and the actual system response of a deposit can be gained from such analyses, as illustrated with the application to two sites from Christchurch, New Zealand.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
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).
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.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
The Mw 6.2 February 22nd 2011 Christchurch earthquake (and others in the 2010-2011 Canterbury sequence) provided a unique opportunity to study the devastating effects of earthquakes first-hand and learn from them for future engineering applications. All major events in the Canterbury earthquake sequence caused widespread liquefaction throughout Christchurch’s eastern suburbs, particularly extensive and severe during the February 22nd event. Along large stretches of the Avon River banks (and to a lesser extent along the Heathcote) significant lateral spreading occurred, affecting bridges and the infrastructure they support. The first stage of this research involved conducting detailed field reconnaissance to document liquefaction and lateral spreading-induced damage to several case study bridges along the Avon River. The case study bridges cover a range of ages and construction types but all are reinforced concrete structures which have relatively short, stiff decks. These factors combined led to a characteristic deformation mechanism involving deck-pinning and abutment back-rotation with consequent damage to the abutment piles and slumping of the approaches. The second stage of the research involved using pseudo-static analysis, a simplified seismic modelling tool, to analyse two of the bridges. An advantage of pseudo-static analysis over more complicated modelling methods is that it uses conventional geotechnical data in its inputs, such as SPT blowcount and CPT cone resistance and local friction. Pseudo-static analysis can also be applied without excessive computational power or specialised knowledge, yet it has been shown to capture the basic mechanisms of pile behaviour. Single pile and whole bridge models were constructed for each bridge, and both cyclic and lateral spreading phases of loading were investigated. Parametric studies were carried out which varied the values of key parameters to identify their influence on pile response, and computed displacements and damages were compared with observations made in the field. It was shown that pseudo-static analysis was able to capture the characteristic damage mechanisms observed in the field, however the treatment of key parameters affecting pile response is of primary importance. Recommendations were made concerning the treatment of these governing parameters controlling pile response. In this way the future application of pseudo-static analysis as a tool for analysing and designing bridge pile foundations in liquefying and laterally spreading soils is enhanced.
This paper examines the consistency of seismicity and ground motion models, used for seismic hazard analysis in New Zealand, with the observations in the Canterbury earthquakes. An overview is first given of seismicity and ground motion modelling as inputs of probabilistic seismic hazard analysis, whose results form the basis for elastic response spectra in NZS1170.5:2004. The magnitude of earthquakes in the Canterbury earthquake sequence are adequately allowed for in the current NZ seismicity model, however the consideration of ‘background’ earthquakes as point sources at a minimum depth of 10km results in up to a 60% underestimation of the ground motions that such events produce. The ground motion model used in conventional NZ seismic hazard analysis is shown to provide biased predictions of response spectra (over-prediction near T=0.2s , and under-predictions at moderate-to-large vibration periods). Improved ground motion prediction can be achieved using more recent NZ-specific models.
The paper discusses modelling of cyclic stress-strain behaviour of soil, in particular a simple model that can produce a desired stiffness and hysteretic damping for a given strain level as observed in laboratory testing is formulated. The unloading-reloading relationship is developed for total stress seismic site response analysis with appropriate damping at large strain. The constitutive model employs a hyperbolic equation as the backbone curve, and uses a modification of the extended Masing unloading-reloading relationship leading to correct measured modulus reduction and damping curves simultaneously. A quasi-static cyclic loading of increasing amplitude is used to demonstrate the model’s performance and its capability to allow improved modelling of the magnitude of energy dissipation based on an experimental program on native sandy soils from Christchurch, New Zealand.
None
This paper concerns the explicit consideration of near-fault directivity in conventional ground motion prediction models, and its implication for probabilistic seismic hazard analysis (PSHA) in New Zealand. The proposed approach utilises recently developed models by Shahi & Baker (2011), which account for both the 'narrowband' nature of the directivity pulse on spectral ordinates, and the probability of pulse occurrence at the site of interest. Furthermore, in order to correctly consider directivity, distributed seismicity sources are considered as finite-faults, as opposed to their (incorrect) conventional treatment as point-sources. The significance of directivity on hazard analysis results is illustrated for various vibration periods at generic sites located in Christchurch and Otira, two locations whose seismic hazard is comprised of notably different seismic sources. When compared to the PSHA results considering directivity and distributed seismicity as finite faults, it is shown that the NZS1170.5:2004 directivity factor is notably unconservative for all vibration periods in Otira (i.e. high seismic hazard region); and unconservative for Christchurch at short-to-moderate vibration periods ( < 3s); but conservative at long periods ( > 4s).
None
In this paper, the characteristics of near-fault ground motions recorded during the Mw7.1 Darfield and Mw 6.2 Christchurch earthquakes are examined and compared with existing empirical models. The characteristics of forward-directivity effects are first examined using a wavelet-based pulse-classification algorithm. This is followed by an assessment of the adequacy of empirical models which aim to capture the effect of directivity effects on amplifying the acceleration response spectra; and the period and peak velocity of the forward-directivity pulse. It is illustrated that broadband directivity models developed by Somerville et al. (1997) and Abrahamson (2000) generally under-predict the observed amplification of response spectral ordinates at longer vibration periods. In contrast, a recently developed narrowband model by Shahi and Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods surrounding the directivity pulse period. Although the empirical predictions of the pulse period are generally favourable for the Christchurch earthquake, the observations from the Darfield earthquake are significantly under-predicted. The elongation in observed pulse periods is inferred as being a result of the soft sedimentary soils of the Canterbury basin. However, empirical predictions of the observed peak velocity associated with the directivity pulse are generally adequate for both events.
The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The aim of this multidisciplinary research was to retrospectively analyse the gastroenteritis prevalence following the February 22, 2011 earthquake in Christchurch. The first focus was to assess whether earthquake-induced infrastructure damage, liquefaction, and gastroenteritis agents spatially explained the recorded gastroenteritis cases over the period of 35 days following the February 22, 2011 earthquake in Christchurch. The gastroenteritis agents considered in this study were Escherichia coli found in the drinking water supply (MPN/100mL) and Non-Compliant Free Associated Chlorine (FAC-NC) (less than <0.02mg/L). The second focus was the protocols that averted a gastroenteritis outbreak at three Emergency Centres (ECs): Burnside High School Emergency Centre (BEC); Cowles Stadium Emergency Centre (CEC); and Linwood High School Emergency Centre (LEC). Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols and indirect themes included type of EC building (school or a sports stadium), and EC staff. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. This research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. A damage profile was created by amalgamating different types of damage for the considered factors for each Census Area Unit (CAU) in Christchurch. The damage profile enabled the application of a variety of statistical methods which included Moran’s I , Hot Spot (HS) analysis, Spearman’s Rho, and Besag–York–Mollié Model using a range of software. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. The ECs were selected to represent the Christchurch area, and were situated where potential for gastroenteritis was high. BEC represented the western side of Christchurch; whilst, CEC and LEC represented the eastern side, where the potential for gastroenteritis was high according to the outputs of the quantitative spatial modelling. Qualitative analysis from the interviews at the ECs revealed that evacuees were arriving at the ECs with gastroenteritis-like symptoms. Participants believed that those symptoms did not originate at the ECs. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols that included prolific use of hand sanitisers; surveillance; and the services offered. Indirect themes included the EC layout, type of EC building (school or a sports stadium), and EC staff. Indirect themes governed the quality and sustainability of the direct themes implemented, which in turn averted gastroenteritis outbreaks at the ECs. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. It was concluded that gastroenteritis point prevalence following the February 22, 2011 earthquake could not be solely explained by earthquake-induced infrastructure damage, liquefaction, and gastroenteritis causative agents alone. However, this research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Creating a damage profile for each CAU and using spatial data analysis can isolate vulnerable areas, and qualitative data analysis provides localised information. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally.
Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.
This paper discusses the seismic performance of the standard RC office building in Christchurch that is given as a structural design example in NZS3101, the concrete structures seismic standard in New Zealand. Firstly the push-over analysis was carried out to evaluate the lateral load carrying capacity of the RC building and then to compare that carrying capacity with the Japanese standard law. The estimated figures showed that the carrying capacity of the New Zealand standard RC office building of NZS3101:2006 was about one third of Japanese demanded carrying capacity. Secondly, time history analysis of the multi-mass system was performed to estimate the maximum response story drift angle using recorded ground motions. Finally, a three-dimensional analysis was carried out to estimate the response of the building to the 22nd February, 2011 Canterbury earthquake. The following outcomes were obtained. 1) The fundamental period of the example RC building is more than twice that of Japanese simplified calculation, 2) The example building’s maximum storey drift angle reached 2.5% under the recorded ground motions. The main purpose of this work is to provide background information of seismic design practice for the reconstruction of Christchurch.
This analysis employs both qualitative and quantitative approaches to identify how young adults in New Zealand aged 18-25 years old have engaged with All Right? campaign material. A survey targeting young adults returned 51 viable out of 117 responses due to participation prerequisites. From the survey, five participants elaborated on their thoughts in an in-depth interview voluntarily. Interviews were conducted with key personnel from All Right? to craft broader understanding of the initiative whilst enhancing knowledge of mental health frameworks and their application. Ciaran Fox, Lucy Daeth and Sara Epperson, who have been imperative to the success of the campaign, shared their working experience in the community and public health sector and how this intertwines to their current roles at All Right?. Discussions of key frameworks, community conversations, the development of communication strategies and how All Right? approached Canterbury publics in a post-earthquake setting provided insight to the importance of understanding community circumstance in initial crisis and the correlated secondary stressors.
The Master of Engineering Management Project was sponsored by the Canterbury Earthquake Recovery Authority (CERA) and consisted of two phases: The first was an analysis of existing information detailing the effects of hazardous natural events on Canterbury Lifeline Utilities in the past 15 years. The aim of this “Lessons Learned” project was to produce an analysis report that identified key themes from the research, gaps in the existing data and to provide recommendations from these “Lessons Learned.” The Second phase was the development of a practical “Disaster Mitigation Guideline” that outlined lessons in the field of Emergency Sanitation. This research would build upon the first stage and would draw from international reference to develop a guideline that has practical implementation possibilities throughout the world.
This thesis studies the behaviour of diaphragms in multi-storey timber buildings by providing methods for the estimation of the diaphragm force demand, developing an Equivalent Truss Method for the analysis of timber diaphragms, and experimentally investigating the effects of displacement incompatibilities between the diaphragm and the lateral load resisting system and developing methods for their mitigation. The need to better understand the behaviour of diaphragms in timber buildings was highlighted by the recent 2010-2011 Canterbury Earthquake series, where a number of diaphragms in traditional concrete buildings performed poorly, compromising the lateral load resistance of the structure. Although shortcomings in the estimation of force demand, and in the analysis and design of concrete floor diaphragms have already been partially addressed by other researchers, the behaviour of diaphragms in modern multi-storey timber buildings in general, and in low damage Pres-Lam buildings (consisting of post-tensioned timber members) in particular is still unknown. The recent demand of mid-rise commercial timber buildings of ten storeys and beyond has further highlighted the lack of appropriate methods to analyse timber diaphragms with irregular floor geometries and large spans made of both light timber framing and massive timber panels. Due to the lower stiffness of timber lateral load resisting systems, compared with traditional construction materials, and the addition of in-plane flexible diaphragms, the effect of higher modes on the global dynamic behaviour of a structure becomes more critical. The results from a parametric non-linear time-history analysis on a series of timber frame and wall structures showed increased storey shear and moment demands even for four storey structures when compared to simplistic equivalent static analysis. This effect could successfully be predicted with methods available in literature. The presence of diaphragm flexibility increased diaphragm inter-storey drifts and the peak diaphragm demand in stiff wall structures, but had less influence on the storey shears and moments. Diaphragm force demands proved to be significantly higher than the forces derived from equivalent static analysis, leading to potentially unsafe designs. It is suggested to design all diaphragms for the same peak demand; a simplified approach to estimate these diaphragm forces is proposed for both frame and wall structures. Modern architecture often requires complex floor geometries with long spans leading to stress concentrations, high force demands and potentially large deformations in the diaphragms. There is a lack of guidance and regulation regarding the analysis and design of timber diaphragms and a practical alternative to the simplistic equivalent deep beam analysis or costly finite element modelling is required. An Equivalent Truss Method for the analysis of both light timber framed and massive timber diaphragms is proposed, based on analytical formulations and verified against finite element models. With this method the panel unit shear forces (shear flow) and therefore the fastener demand, chord forces and reaction forces can be evaluated. Because the panel stiffness and fastener stiffness are accounted for, diaphragm deflection, torsional effects and transfer forces can also be assessed. The proposed analysis method is intuitive and can be used with basic analysis software. If required, it can easily be adapted for the use with diaphragms working in the non-linear range. Damage to floor diaphragms resulting from displacement incompatibilities due to frame elongation or out-of plane deformation of walls can compromise the transfer of inertial forces to the lateral load resisting system as well as the stability of other structural elements. Two post-tensioned timber frame structures under quasi-static cyclic and dynamic load, respectively, were tested with different diaphragm panel layouts and connections investigating their ability to accommodate frame elongations. Additionally, a post-tensioned timber wall was loaded under horizontal cyclic loads through two pairs of collector beams. Several different connection details between the wall and the beams were tested, and no damage to the collector beams or connections was observed in any of the tests. To evaluate the increased strength and stiffness due to the wall-beam interaction an analytical procedure is presented. Finally, a timber staircase core was tested under bi-directional loading. Different connection details were used to study the effect of displacement incompatibilities between the orthogonal collector beams. These experiments showed that floor damage due to displacement incompatibilities can be prevented, even with high levels of lateral drift, by the flexibility of well-designed connections and the flexibility of the timber elements. It can be concluded that the flexibility of timber members and the flexibility of their connections play a major role in the behaviour of timber buildings in general and of diaphragms specifically under seismic loads. The increased flexibility enhances higher mode effects and alters the diaphragm force demand. Simple methods are provided to account for this effect on the storey shear, moment and drift demands as well as the diaphragm force demands. The analysis of light timber framing and massive timber diaphragms can be successfully analysed with an Equivalent Truss Method, which is calibrated by accounting for the panel shear and fastener stiffnesses. Finally, displacement incompatibilities in frame and wall structures can be accommodated by the flexibilities of the diaphragm panels and relative connections. A design recommendations chapter summarizes all findings and allows a designer to estimate diaphragm forces, to analyse the force path in timber diaphragms and to detail the connections to allow for displacement incompatibilities in multi-storey timber buildings.
INTRODUCTION: Connections between environmental factors and mental health issues have been postulated in many different countries around the world. Previously undertaken research has shown many possible connections between these fields, especially in relation to air quality and extreme weather events. However, research on this subject is lacking in New Zealand, which is difficult to analyse as an overall nation due to its many micro-climates and regional differences.OBJECTIVES: The aim of this study and subsequent analysis is to explore the associations between environmental factors and poor mental health outcomes in New Zealand by region and predict the number of people with mental health-related illnesses corresponding to the environmental influence.METHODS: Data are collected from various public-available sources, e.g., Stats NZ and Coronial services of New Zealand, which comprised four environmental factors of our interest and two mental health indicators data ranging from 2016 up until 2020. The four environmental factors are air pollution, earthquakes, rainfall and temperature. Two mental health indicators include the number of people seen by District Health Boards (DHBs) for mental health reasons and the statistics on suicide deaths. The initial analysis is carried out on which regions were most affected by the chosen environmental factors. Further analysis using Auto-Regressive Integrated Moving Average(ARIMA) creates a model based on time series of environmental data to generate estimation for the next two years and mental health projected from the ridge regression.RESULTS: In our initial analysis, the environmental data was graphed along with mental health outcomes in regional charts to identify possible associations. Different regions of New Zealand demonstrate quite different relationships between the environmental data and mental health outcomes. The result of later analysis predicts that the suicide rate and DHB mental health visits may increase in Wellington, drop-in Hawke's Bay and slightly increase in Canterbury for the year 2021 and 2022 with different environmental factors considered.CONCLUSION: It is evident that the relationship between environmental and mental health factors is regional and not national due to the many micro-climates that exist around the nation. However, it was observed that not all factors displayed a good relationship between the regions. We conclude that our hypotheses were partially correct, in that increased air pollution was found to correlate to increased mental health-related DHB visits. Rainfall was also highly correlated to some mental health outcomes. Higher levels of rainfall reduced DHB visits and suicide rates in some areas of the country.
Research on human behaviour during earthquake shaking has identified three main influences of behaviour: the environment the individual is located immediately before and during the earthquake, in terms of where the individual is and who the individual is with at the time of the earthquake; individual characteristics, such as age, gender, previous earthquake experience, and the intensity and duration of earthquake shaking. However, little research to date has systematically analysed the immediate observable human responses to earthquake shaking, mostly due to data constraints and/or ethical considerations. Research on human behaviour during earthquakes has relied on simulations or post-event, reflective interviews and questionnaire studies, often performed weeks to months or even years following the event. Such studies are therefore subject to limitations such as the quality of the participant's memory or (perceived) realism of a simulation. The aim of this research was to develop a robust coding scheme to analyse human behaviour during earthquake shaking using video footage captured during an earthquake event. This will allow systematic analysis of individuals during real earthquakes using a previously unutilized data source, thus help develop guidance on appropriate protective actions. The coding scheme was developed in a two-part process, combining a deductive and inductive approach. Previous research studies of human behavioral response during earthquake shaking provided the basis for the coding scheme. This was then iteratively refined by applying the coding scheme to a broad range of video footage of people exposed to strong shaking during the Canterbury earthquake sequence. The aim of this was to optimise coding scheme content and application across a broad range of scenarios, and to increase inter-coder reliability. The methodology to code data will enhance objective observation of video footage to allow cross-event analysis and explore (among others): reaction time, patterns of behaviour, and social, environmental and situational influences of behaviour. This can provide guidance for building configuration and design, and evidence-based recommendations for public education about injury-preventing behavioural responses during earthquake shaking.
