The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
Introduction This poster presents the inferred initial performance and recovery of the water supply network of Christchurch following the 22 February 2011 Mw 6.2 earthquake. Results are presented in a geospatial and temporal fashion. This work strengthens the current understanding of the restoration of such a system after a disaster and quantifies the losses caused by this earthquake in respect with the Christchurch community. Figure 1 presents the topology of the water supply network as well as the spatial distribution of the buildings and their use.
This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.
In the top frames two children shout 'four point one', 'three point six', and 'five point two' and in the frame below it is seen that they are responding to bumps in the road as their mother drives through Christchurch streets. Context - The children have become expert at guessing the seismic intensity of earthquakes in Christchurch and are now applying them to bumps in the road. Colour and black and white versions available Quantity: 2 digital cartoon(s).
Asset management in power systems is exercised to improve network reliability to provide confidence and security for customers and asset owners. While there are well-established reliability metrics that are used to measure and manage business-as-usual disruptions, an increasing appreciation of the consequences of low-probability high-impact events means that resilience is increasingly being factored into asset management in order to provide robustness and redundancy to components and wider networks. This is particularly important for electricity systems, given that a range of other infrastructure lifelines depend upon their operation. The 2010-2011 Canterbury Earthquake Sequence provides valuable insights into electricity system criticality and resilience in the face of severe earthquake impacts. While above-ground assets are relatively easy to monitor and repair, underground assets such as cables emplaced across wide areas in the distribution network are difficult to monitor, identify faults on, and repair. This study has characterised in detail the impacts to buried electricity cables in Christchurch resulting from seismically-induced ground deformation caused primarily by liquefaction and lateral spread. Primary modes of failure include cable bending, stretching, insulation damage, joint braking and, being pulled off other equipment such as substation connections. Performance and repair data have been compiled into a detailed geospatial database, which in combination with spatial models of peak ground acceleration, peak ground velocity and ground deformation, will be used to establish rigorous relationships between seismicity and performance. These metrics will be used to inform asset owners of network performance in future earthquakes, further assess component criticality, and provide resilience metrics.
As a result of the 4 September 2010 Darfield earthquake and the more damaging 22 February 2011 Christchurch earthquake, considerable damage occurred to a significant number of buildings in Christchurch. The damage that occurred to the Christchurch Roman Catholic Cathedral of the Blessed Sacrament (commonly known as the Christchurch Basilica) as a result of the Canterbury earthquakes is reported, and the observed failure modes are identified. A previous strengthening intervention is outlined and the estimated capacity of the building is discussed. This strengthening was completed in 2004, and addressed the worst aspects of the building's seismic vulnerability. Urgent work was undertaken post-earthquake to secure parts of the building in order to limit damage and prevent collapse of unstable parts of the building. The approach taken for this securing is outlined, and the performance of the building and the previously installed earthquake strengthening intervention is evaluated.A key consideration throughout the project was the interaction between the structural securing requirements that were driven by the requirement to limit damage and mitigate hazards, and the heritage considerations. Lessons learnt from the strengthening that was carried out, the securing work undertaken, and the approach taken in making the building "safe" are discussed. Some conclusions are drawn with respect to the effectiveness of strengthening similar building types, and the approach taken to secure the building under active seismic conditions. AM - Accepted Manuscript
Within four weeks of the September 4 2010 Canterbury Earthquake a new, loosely-knit community group appeared in Christchurch under the banner of “Greening the Rubble.” The general aim of those who attended the first few meetings was to do something to help plug the holes that had already appeared or were likely to appear over the coming weeks in the city fabric with some temporary landscaping and planting projects. This article charts the first eighteen months of Greening the Rubble and places the initiative in a broader context to argue that although seismic events in Christchurch acted as a “call to palms,” so to speak, the city was already in need of some remedial greening. It concludes with a reflection on lessons learned to date by GTR and commentary on the likely issues ahead for this new mini-social-environmental movement in the context of a quake-affected and still quake-prone major New Zealand city. One of the key lessons for GTR and all of those involved in Christchurch recovery activities to date is that the city is still very much in the middle of the event and is to some extent a laboratory for seismic and agency management studies alike.
A Christchurch seismic monitor churns out reams of paper with wild fluctuations showing earthquakes and aftershocks. Nearby two men examine a second monitor and one of them says 'This one's connected to Gerry Brownlee!' The reams of paper show a perfectly straight line - it appears that nothing is going on in Gerry Brownlee's head. Context - A sense that the Minister for Earthquake Recovery is being less than effective. Colour and black and white versions available Quantity: 2 digital cartoon(s).
The nonlinear dynamic soil-foundation-structure interaction (SFSI) can signifi cantly affect the seismic response of buildings, causing additional deformation modes, damage and repair costs. Because of nonlinear foundation behaviour and interactions, the seismic demand on the superstructure may considerably change, and also permanent deformations at the foundation level may occur. Although SFSI effects may be benefi cial to the superstructure performance, any advantage would be of little structural value unless the phenomenon can be reliably controlled and exploited. Detrimental SFSI effects may also occur, including acceleration and displacement response ampli cation and differential settlements, which would be unconservative to neglect. The lack of proper understanding of the phenomenon and the limited available simpli ed tools accounting for SFSI have been major obstacles to the implementation of integrated design and assessment procedures into the everyday practice. In this study concepts, ideas and practical tools (inelastic spectra) for the seismic design and assessment of integrated foundation-superstructure systems are presented, with the aim to explicitly consider the impact of nonlinearities occurring at the soil-foundation interface on the building response within an integrated approach, where the foundation soil and superstructure are considered as part of an integrated system when evaluating the seismic response, working synergically for the achievement of a target global performance. A conceptual performance-based framework for the seismic design and assessment of integrated foundation-superstructure systems is developed. The framework is based on the use of peak and residual response parameters for both the superstructure and the foundation, which are then combined to produce the system performance matrix. Each performance matrix allows for worsening of the performance when different contributions are combined. An attempt is made to test the framework by using case histories from the 2011 Christchurch earthquake, which are previously shown to have been severely affected by nonlinear SFSI. The application highlights the framework sensitivity to the adopted performance limit states, which must be realistic for a reliable evaluation of the system performance. Constant ductility and constant strength inelastic spectra are generated for nonlinear SFSI systems (SDOF nonlinear superstructure and 3DOF foundation allowing for uplift and soil yielding), representing multistorey RC buildings with shallow rigid foundations supported by cohesive soils. Different ductilities/strengths, hysteretic rules (Bi-linear, Takeda and Flag-Shape), soil stiffness and strength and bearing capacity factors are considered. Footings and raft foundations are investigated, characterized respectively by constant (3 and 8) and typically large bearing capacity factors. It is confi rmed that when SFSI is considered, the superstructure yielding force needed to satisfy a target ductility for a new building changes, and that similarly, for an existing building, the ductility demand on a building of a given strength varies. The extent of change of seismic response with respect to xed-base (FB) conditions depends on the class of soils considered, and on the bearing capacity factor (SF). For SF equal to 3, the stiffer soils enhance the nonlinear rotational foundation behaviour and are associated with reduced settlement, while the softer ones are associated with increased settlement response but not signi ficant rotational behaviour. On average terms, for the simplifi ed models considered, SFSI is found to be bene cial to the superstructure performance in terms of acceleration and superstructure displacement demand, although exceptions are recorded due to ground motion variability. Conversely, in terms of total displacement, a signi cant response increase is observed. The larger the bearing capacity factor, the more the SFSI response approaches the FB system. For raft foundation buildings, characterized by large bearing capacity factors, the impact of foundation response is mostly elastic, and the system on average approaches FB conditions. Well de fined displacement participation factors to the peak total lateral displacement are observed for the different contributions (i.e. peak foundation rotation and translation and superstructure displacement). While the superstructure and foundation rotation show compensating trends, the foundation translation contribution varies as a function of the moment-to-shear ratio, becoming negligible in the medium-to-long periods. The longer the superstructure FB period, the less the foundation response is signifi cant. The larger the excitation level and the less ductile the superstructure, the larger the foundation contribution to the total lateral displacement, and the less the superstructure contribution. In terms of hysteretic behaviour, its impact is larger when the superstructure response is more signifi cant, i.e. for the softer/weaker soils and larger ductilities. Particularly, for the Flag Shape rule, larger superstructure displacement participation factors and smaller foundation contributions are recorded. In terms of residual displacements, the total residual-to-maximum ratios are similar in amplitudes and trends to the corresponding FB system responses, with the foundation and superstructure contributions showing complementary trends. The impact of nonlinear SFSI is especially important for the Flag Shape hysteresis rule, which would not otherwise suffer of any permanent deformations. By using the generated peak and residual inelastic spectra (i.e. inelastic acceleration/ displacement modifi cation factor spectra, and/or participation factor and residual spectra), conceptual simplifi ed procedures for the seismic design and assessment of integrated foundation-superstructure systems are presented. The residual displacements at both the superstructure and foundation levels are explicitly considered. Both the force- and displacement-based approaches are explored. The procedures are de fined to be complementary to the previously proposed integrated performance-based framework. The use of participation factor spectra allows the designer to easily visualize the response of the system components, and could assist the decision making process of both the design and assessment of SFSI systems. The presented numerical results have been obtained using simpli ed models, assuming rigid foundation behaviour and neglecting P-Delta effects. The consideration of more complex systems including asymmetry in stiffness, mass, axial load and ground conditions with a exible foundation layout would highlight detrimental SFSI effects as related to induced differential settlements, while accounting for PDelta effects would further amplify the displacement response. Also, the adopted acceleration records were selected and scaled to match conventional design spectra, thus not representing any response ampli cation in the medium-to-long period range which could as well cause detrimental SFSI effects. While these limitations should be the subject of further research, this study makes a step forward to the understanding of SFSI phenomenon and its incorporation into performance-based design/assessment considerations.
Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30
An overview of the 22 February 2011 Christchurch earthquake is presented in the context of characterization of extreme/rare events. Focus is given to the earthquake source, observed near-source strong ground motions, and effects of site response, while structural response and consequences are mentioned for completeness. For each of the above topics comparisons and discussions are made with predictive models for each of phenomena considered. In light of the observations and predictive model comparisons, the author’s opinion on improving the characterization of such extreme/rare events, and their appropriate consideration in seismic design is presented
This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.
Perimeter Moment resisting steel frames (PMRSFs) are a commonly used seismic resisting system, placed around the perimeter of the building for maximum torsional stiffness. They are typically designed as “strong column weak beam” systems with fixed column bases. When subjected to severe earthquake demand, sufficient to push the beams into the inelastic range, it is expected that plastic hinging at the column bases will occur. However, the response of PMRSF systems to the severe 2010/2011 Christchurch earthquake series did not generate column base hinging in systems which exhibited beam yielding.
An earthquake engineering expert wants to change the way we predict how the ground will shake during an earthquake. Professor Brendon Bradley from the University of Canterbury is the recipient of a Marsden Fund grant to accelerate his research into seismic hazard analysis and forecasting. He says the idea is to get to a point where they can provide the same sort of information as a weather forecast. Professor Bradley says just like a severe weather warning, engineers would be able to provide information about severe ground shaking, how it varies locally in each city or suburb, and the likely consequence to buildings. Kathryn speaks to Professor Brendon Bradley, the director of Te Hirangu Ru QuakeCoRE - The New Zealand Centre for Earthquake Resilience.
In the aftermath of the 22 February 2011 earthquake, the Natural Hazards Research Platform (NHRP) initiated a series of Short Term Recovery Projects (STRP) aimed at facilitating and supporting the recovery of Christchurch from the earthquake impacts. This report presents the outcomes of STRP 6: Impacts of Liquefaction on Pipe Networks, which focused on the impacts of liquefaction on the potable water and wastewater systems of Christchurch. The project was a collaborative effort of NHRP researchers with expertise in liquefaction, CCC personnel managing and designing the systems and a geotechnical practitioner with experience/expertise in Christchurch soils and seismic geotechnics.
The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
The Global Earthquake Model’s (GEM) Earthquake Consequences Database (GEMECD) aims to develop, for the first time, a standardised framework for collecting and collating geocoded consequence data induced by primary and secondary seismic hazards to different types of buildings, critical facilities, infrastructure and population, and relate this data to estimated ground motion intensity via the USGS ShakeMap Atlas. New Zealand is a partner of the GEMECD consortium and to-date has contributed with 7 events to the database, of which 4 are localised in the South Pacific area (Newcastle 1989; Luzon 1990; South of Java 2006 and Samoa Islands 2009) and 3 are NZ-specific events (Edgecumbe 1987; Darfield 2010 and Christchurch 2011). This contribution to GEMECD represented a unique opportunity for collating, comparing and reviewing existing damage datasets and harmonising them into a common, openly accessible and standardised database, from where the seismic performance of New Zealand buildings can be comparatively assessed. This paper firstly provides an overview of the GEMECD database structure, including taxonomies and guidelines to collect and report on earthquake-induced consequence data. Secondly, the paper presents a summary of the studies implemented for the 7 events, with particular focus on the Darfield (2010) and Christchurch (2011) earthquakes. Finally, examples of specific outcomes and potentials for NZ from using and processing GEMECD are presented, including: 1) the rationale for adopting the GEM taxonomy in NZ and any need for introducing NZ-specific attributes; 2) a complete overview of the building typological distribution in the Christchurch CBD prior to the Canterbury earthquakes and 3) some initial correlations between the level and extent of earthquake-induced physical damage to buildings, building safety/accessibility issues and the induced human casualties.
This study led on from Earthquake hazard and risk assessment study Stage 1 Part A: Earthquake source identification and characterisation (Pettinga et al, 1998). It used the location and characteristics of active faults in the Canterbury region, and the historic record of earthquakes to estimate levels of ground shaking (MM intensity, peak ground acceleration and spectral accelerations) across Canterbury for different return periods. The study also provided earthquake scenarios for selected towns and cities in Canterbury, and undertook detailed investigations into the largest historic earthquakes in Christchurch and parts of the Canterbury region. See Object Overview for background and usage information.
Results from a series of 1D seismic effective stress analyses of natural soil deposits from Christchurch are summarized. The analysed soil columns include sites whose performance during the 2010-2011 Canterbury earthquakes varied significantly, from no liquefaction manifestation at the ground surface to very severe liquefaction, in which case a large area of the site was covered by thick soil ejecta. Key soil profile characteristics and response mechanisms affecting the severity of surface liquefaction manifestation and subsequent damage are explored. The influence of shaking intensity on the triggering and contribution of these mechanisms is also discussed. Careful examination of the results highlights the importance of considering the deposit as a whole, i.e. a system of layers, including interactions between layers in the dynamic response and through pore water pressure redistribution and water flow.
The seismic response of unreinforced masonry (URM) buildings, in both their as-built or retrofitted configuration, is strongly dependent on the characteristics of wooden floors and, in particular, on their in-plane stiffness and on the quality of wall-to-floor connections. As part of the development of alternative performance-based retrofit strategies for URM buildings, experimental research has been carried out by the authors at the University of Canterbury, in order to distinguish the different elements contributing to the whole diaphragm's stiffness. The results have been compared to the ones predicted through the use of international guidelines in order to highlight shortcomings and qualities and to propose a simplified formulation for the evaluation of the stiffness properties.
Existing New Zealand (NZ) building stock contains a significant number of structures designed prior to 1995 with non-ductile reinforced concrete (RC) columns. Recent earthquakes and research show that columns with such details perform poorly when subjected to seismic demand, losing gravity load carrying capacity at drift levels lower than the expected one. Therefore, in order to have a better understanding of existing RC columns in NZ, the history of these elements is investigated in this paper. The evolution of RC column design guidelines in NZ standards since the 1970s is scrutinized. For this purpose, a number of RC columns from Christchurch buildings built prior to 1995 are assessed using the current code of practice.
On Tuesday 22 February 2011, a 6.3 magnitude earthquake struck Christchurch, New Zealand’s second largest city. The ‘earthquake’ was in fact an aftershock to an earlier 7.1 magnitude earthquake that had occurred on Saturday 4 September 2010. There were a number of key differences between the two events that meant they had dramatically different results for Christchurch and its inhabitants. The 22 February 2011 event resulted in one of New Zealand’s worst natural disasters on record, with 185 fatalities occurring and hundreds more being injured. In addition, a large number of buildings either collapsed or were damaged to the point where they needed to be totally demolished. Since the initial earthquake in September 2010, a large amount of building-related research has been initiated in New Zealand to investigate the impact of the series of seismic events – the major focus of these research projects has been on seismic, structural and geotechnical engineering matters. One project, however, conducted jointly by the University of Canterbury, the Fire Protection Association of New Zealand and BRANZ, has focused on the performance of fire protection systems in the earthquakes and the effectiveness of the systems in the event of post-earthquake fires occurring. Fortunately, very few fires actually broke out following the series of earthquake events in Christchurch, but fire after earthquakes still has significant implications for the built environment in New Zealand, and the collaborative research has provided some invaluable insight into the potential threat posed by post-earthquake fires in buildings. As well as summarising the damage caused to fire protection systems, this paper discusses the flow-on effect for designing structures to withstand post-earthquake fires. One of the underlying issues that will be explored is the existing regulatory framework in New Zealand whereby structural earthquake design and structural design for fire are treated as discrete design scenarios.
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.
We present preliminary observations on three waters impacts from the Mw7.8 14th November 2016 Kaikōura Earthquake on wider metropolitan Wellington, urban and rural Marlborough, and in Kaikōura township. Three waters systems in these areas experienced widespread and significant transient ground deformation in response to seismic shaking, with localised permanent ground deformation via liquefaction and lateral spreading. In Wellington, potable water quality was impacted temporarily by increased turbidity, and significant water losses occurred due to damaged pipes at the port. The Seaview and Porirua wastewater treatment plants sustained damage to clarifier tanks from water seiching, and increased water infiltration to the wastewater system occurred. Most failure modes in urban Marlborough were similar to the 2010-2011 Canterbury Earthquake Sequence; however some rural water tanks experienced rotational and translational movements, highlighting importance of flexible pipe connections. In Kaikōura, damage to reservoirs and pipes led to loss of water supply and compromised firefighting capability. Wastewater damage led to environmental contamination, and necessitated restrictions on greywater entry into the system to minimise flows. Damage to these systems necessitated the importation of tankered and bottled water, boil water notices and chlorination of the system, and importation of portaloos and chemical toilets. Stormwater infrastructure such as road drainage channels was also damaged, which could compromise condition of underlying road materials. Good operational asset management practices (current and accurate information, renewals, appreciation of criticality, good system knowledge and practical contingency plans) helped improve system resilience, and having robust emergency management centres and accurate Geographic Information System data allowed effective response coordination. Minimal damage to the wider built environment facilitated system inspections. Note Future research will include detailed geospatial assessments of seismic demand on these systems and attendant modes of failure, levels of service restoration, and collaborative development of resilience measures.
On November 14 2016 a magnitude 7.8 earthquake struck the south island of New Zealand. The earthquake lasted for just two minutes with severe seismic shaking and damage in the Hurunui and Kaikōura districts. Although these are predominantly rural areas, with scattered small towns and mountainous topography, they also contain road and rail routes that are essential parts of the national transport infrastructure. This earthquake and the subsequent recovery are of particular significance as they represent a disaster following in close proximity to another similar disaster, with the Canterbury earthquakes occurring in a neighboring district five years earlier. The research used an inductive qualitative case study to explore the nature of the Kaikōura recovery. That recovery process involved a complex interplay between the three parties; (a) the existing local government in the district, (b) central government agencies funding the recovery of the local residents and the national transport infrastructure, and (c) recovery leaders arriving with recent expertise from the earlier Canterbury disaster. It was evident that three groups: locals, government, and experts represented a multi-party governance debate in which the control of the Kaikōura earthquake recovery was shared amongst them. Each party had their own expertise, adgenda and networks that they brought to the Kaikōura recovery, but this created tensions between external expertise and local, community leadership. Recent earthquake research suggests that New Zealand is currently in the midst of an earthquake cluster, with further seismic disasters likely to occur in relatively close succession. This is likely to be compounded by the increasing frequency of other natural disasters with the effects of climate change. The present study investigates a phenomenon that may become increasingly common, with the transfer of disaster expertise from one event to another, and the interface between those experts with local and national government in directing recoveries. The findings of this study have implications for practitioners and policy makers in NZ and other countries where disasters are experienced in close spatial and temporal proximity.
Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults. The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains. The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future. Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case. This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features. Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.
The Avon and Heathcote Rivers, located in the city of Christchurch, New Zealand, are lowland spring-fed rivers linked with the Christchurch Groundwater System. At present, the flow paths and recharge sources to the Christchurch Groundwater System are not fully understood. Study of both the Avon and Heathcote Rivers can provide greater insight into this system. In addition, during the period 2010-2012, Christchurch has experienced large amounts of seismic activity, including a devastating Mw 6.2 aftershock on February 22nd, 2011, which caused widespread damage and loss of life. Associated with these earthquakes was the release of large amounts of water through liquefaction and temporary springs throughout the city. This provided a unique opportunity to study groundwater surface water interactions following a large scale seismic event. Presented herein is the first major geochemical study on the Avon and Heathcote Rivers and the hydrological impact of the February 22, 2011 Christchurch Earthquake. The Avon, Heathcote, and Waimakariri Rivers were sampled in quarterly periods starting in July 2011 and analyzed for stable Isotopes δ¹⁸O, δD, and δ¹³C and major anion composition. In addition, post -earthquake samples were collected over the days immediately following the February 22, 2011 earthquake and analyzed for stable isotopes δ¹⁸O and δD and major anion composition. A variety of analytical methods were used identify the source of the waters in the Avon-Heathcote System and evaluate the effectiveness of stable isotopes as geochemical tracers in the Christchurch Groundwater System. The results of this thesis found that the waters from the Avon and Heathcote Rivers are geochemically the same, originating from groundwater, and exhibit a strong tidal influence within 5km of the Avon-Heathcote Estuary. The surface waters released following the February 22nd, 2011 earthquake were indistinguishable from quarterly samples taken from the Avon and Heathcote Rivers when comparing stable isotopic composition. The anion data suggests the waters released following the February 22nd, 2011 Christchurch Earthquake were sourced primarily from shallow groundwater, and also suggests a presence of urban sewage at some sites. Attempts to estimate recharge sources for the Avon-Heathcote Rivers using published models for the Christchurch Groundwater System yielded results that were not consistent between models. In evaluating the use of geochemical constituents as tracers in the Christchurch Groundwater System, no one isotope could provide a clear resolution, but when used in conjunction, δ¹⁸O, δ¹³C, and DIC, seem to be the most effective tracers. Sample sizes for δ¹³C were too small for a robust evaluation. Variability on the Waimakariri River appears to be greater than previously estimated, which could have significant impacts on geochemical models for the Christchurch Groundwater System. This research demonstrates the value of using multiple geochemical constituents to enrich our understanding of the groundwater surfaces-water interactions and the Christchurch Groundwater System as a whole.
The November 2016 MW 7.8 Kaikōura Earthquake initiated beneath the North Culverden basin on The Humps fault and propagated north-eastwards, rupturing at least 17 faults along a cumulative length of ~180 km. The geomorphic expression of The Humps Fault across the Emu Plains, along the NW margin of Culverden basin, comprises a series of near-parallel strands separated by up to 3 km across strike. The various strands strike east to east-northeast and have been projected to mainly dip steeply to the south in seismic data (~80°). In this area, the fault predominantly accommodates right-lateral slip, with uplift and subsidence confined to releasing and restraining bends and step-overs at a range of scales. The Kaikōura event ruptured pre-existing fault scarps along the Emu Plains, which had been partly identified prior to the earthquake. Geomorphology and faulting expression of The Humps Fault on The Emu Plains was mapped, along with faulting related structures which did not rupture in the 2016 earthquake. Fault ruptures strands are combined into sections and the kinematic deformation of sections analysed to provide a moment tensor fault plane solution. This fault plane solution is consistent with the regional principal horizontal shortening direction (PHS) of ~115°, similar to seismic focal mechanism solutions of some of the nearby aftershocks of the Kaikōura earthquake, and similar to the adjacent Hope Fault. To constrain the timing of paleoseismic events, a trench was excavated across the fault where it crossed a late Quaternary alluvial fan. Mapping of stratigraphy exposed in the trench walls, and dating of variably deformed strata, constrains the pre-historic earthquake event history at the trench site. The available data provides evidence for at least three paleo-earthquakes within the last 15.1 ka, with a possible fourth (penultimate) event. These events are estimated to have occurred at 7.7-10.3 ka, 10.3-14.8 ka, and one or more events that are older than ~15.1 ka. Some evidence suggests an additional penultimate event between 1850 C.E and 7.7 ka. Time-integrated slip-rates at three locations on the fault are measured using paleo-channels as piercing points. These sites give horizontal slip rates of 0.57 ± 0.1 mm/year, 0.49 ± 0.1 mm/year and one site constrains a minimum of between 0.1 - 0.4 mm/year. Two vertical slip-rates are calculated to be constrained to a maximum of 0.2 ± 0.02 mm/year at one site and between 0.02 and 0.1 mm/year at another site. Prior to this study, The Humps fault had only been partially documented in reconnaissance level mapping in the district, and no previous paleoseismic or slip rate data had been reported. This project has provided a detailed fault zone tectonic geomorphic map and established new slip-rate and paleoseismic data. The results highlight that The Humps fault plays an important role in regional seismicity and in accommodating plate boundary deformation across the North Canterbury region.