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.
During the Christchurch earthquake of February 2011, several midrise buildings of Reinforced Concrete Masonry (RCM) construction achieved performance levels in the range of life safety to near collapse levels. These buildings were subjected to seismic demands higher than the building code requirements of the time and higher than the current New Zealand Loadings Standard (NZS-1170.5:2004). Structural damage to these buildings has been documented and is currently being studied to establish lessons to be learned from their performance and how to incorporate these lessons into future RCM design and construction practices. This paper presents a case study of a six story RCM building deemed to have reached the near collapse performance level. The RCM walls on the 2nd floor failed due to toe crushing reducing the building’s lateral resistance in the east-west direction. A nonlinear dynamic analysis on a 3D model was conducted to simulate the development of the governing failure mechanism. Preliminary analysis results show that the damaged walls were initially under large compression forces from gravity loads which caused increase in their lateral strength and reduced their ductility. After toe crushing failure developed, axial instability of the model was prevented by a redistribution of gravity loads.
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.
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.
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.
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.
The Catholic Cathedral is classified as a category 1 listed heritage building constructed largely of unreinforced stone masonry, and was significantly damaged in the recent Canterbury earthquakes of 2010 and 2011. In the 2010 event the building presented slight to moderta damage, meanwhile in the 2011 one experienced ground shaking in excess of its capacity leading to block failures and partial collapse of parts of the building, which left the building standing but still posing a significant hazard. In this paper we discuss the approach to develop the earthquake analysis of the building by 3D numerical simulations, and the results are compared/calibrated with the observed damage of the 2010 earthquake. Very accurate records were obtained during both earthquakes due to a record station located least than 80 m of distance from the building and used in the simulations. Moreover it is included in the model the soil structure interaction because it was observed that the ground and foundation played an important role on the seismic behavior of the structure. A very good agreement was found between the real observed damage and the nonlinear dynamic simulations described trough inelastic deformation (cracking) and building´s performance.
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.
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).
© 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.
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 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 shows an understanding of the availability of resources in post-disaster reconstruction and recovery in Christchurch, New Zealand following its September 4, 2010 and February 22, 2011 earthquakes. Overseas experience in recovery demonstrates how delays and additional costs may incur if the availability of resources is not aligned with the reconstruction needs. In the case of reconstruction following Christchurch earthquakes, access to normal resource levels will be insufficient. An on-line questionnaire survey, combined with in-depth interviews was used to collect data from the construction professionals that had been participated in the post-earthquake reconstruction. The study identified the resources that are subject to short supply and resourcing challenges that are currently faced by the construction industry. There was a varied degree of impacts felt by the surveyed organisations from resource shortages. Resource pressures were primarily concentrated on human resources associated with structural, architectural and land issues. The challenges that may continue playing out in the longer-term reconstruction of Christchurch include limited capacity of the construction industry, competition for skills among residential, infrastructure and commercial sectors, and uncertainties with respect to decision making. Findings provide implications informing the ongoing recovery and rebuild in New Zealand. http://www.iiirr.ucalgary.ca/Conference-2012
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 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.
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 Canterbury region experienced widespread damage due to liquefaction induced by seismic shaking during the 4 September 2010 earthquake and the large aftershocks that followed, notably those that occurred on 22 February, 13 June and 23 December 2011. Following the 2010 earthquake, the Earthquake Commission directed a thorough investigation of the ground profile in Christchurch, and to date, more than 7500 cone penetration tests (CPT) have been performed in the region. This paper presents the results of analyses which use a subset of the geotechnical database to evaluate the liquefaction process as well as the re-liquefaction that occurred following some of the major events in Christchurch. First, the applicability of existing CPT-based methods for evaluating liquefaction potential of Christchurch soils was investigated using three methods currently available. Next, the results of liquefaction potential evaluation were compared with the severity of observed damage, categorised in terms of the land damage grade developed from Tonkin & Taylor property inspections as well as from observed severity of liquefaction from aerial photography. For this purpose, the Liquefaction Potential Index (LPI) was used to represent the damage potential at each site. In addition, a comparison of the CPT-based strength profiles obtained before each of the major aftershocks was performed. The results suggest that the analysis of spatial and temporal variations of strength profiles gives a clear indication of the resulting liquefaction and re-liquefaction observed in Christchurch. The comparison of a limited number of CPT strength profiles before and after the earthquakes seems to indicate that no noticeable strengthening has occurred in Christchurch, making the area vulnerable to liquefaction induced land damage in future large-scale earthquakes.
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
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).
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.
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.
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.
Utility managers are always looking for appropriate tools to estimate seismic damage in wastewater networks located in earthquake prone areas. Fragility curves, as an appropriate tool, are recommended for seismic vulnerability analysis of buried pipelines, including pressurised and unpressurised networks. Fragility curves are developed in pressurised networks mainly for water networks. Fragility curves are also recommended for seismic analysis in unpressurised networks. Applying fragility curves in unpressurised networks affects accuracy of seismic damage estimation. This study shows limitations of these curves in unpressurised networks. Multiple case study analysis was applied to demonstrate the limitations of the application of fragility curves in unpressurised networks in New Zealand. Four wastewater networks within New Zealand were selected as case studies and various fragility curves used for seismic damage estimation. Observed damage in unpressurised networks after the 2007 earthquake in Gisborne and the 2010 earthquake in Christchurch demonstrate the appropriateness of the applied fragility curves to New Zealand wastewater networks. This study shows that the application of fragility curves, which are developed from pressurised networks, cannot be accurately used for seismic damage assessment in unpressurised wastewater networks. This study demonstrated the effects of different parameters on seismic damage vulnerability of unpressurised networks.
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.
