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.
Capacity design and hierarchy of strength philosophies at the base of modern seismic codes allow inelastic response in case of severe earthquakes and thus, in most traditional systems, damage develops at well-defined locations of reinforced concrete (RC) structures, known as plastic hinges. The 2010 and 2011 Christchurch earthquakes have demonstrated that this philosophy worked as expected. Plastic hinges formed in beams, in coupling beams and at the base of columns and walls. Structures were damaged permanently, but did not collapse. The 2010 and 2011 Christchurch earthquakes also highlighted a critical issue: the reparability of damaged buildings. No methodologies or techniques were available to estimate the level of subsequent earthquakes that RC buildings could still sustain before collapse. No repair techniques capable of restoring the initial condition of buildings were known. Finally, the cost-effectiveness of an eventual repair intervention, when compared with a new building, was unknown. These aspects, added to nuances of New Zealand building owners’ insurance coverage, encouraged the demolition of many buildings. Moreover, there was a perceived strong demand from government and industry to develop techniques for assessing damage to steel reinforcement bars embedded in cracked structural concrete elements. The most common questions were: “Have the steel bars been damaged in correspondence to the concrete cracks?”, “How much plastic deformation have the steel bars undergone?”, and “What is the residual strain capacity of the damaged bars?” Minimally invasive techniques capable of quantifying the level and extent of plastic deformation and residual strain capacity are not yet available. Although some studies had been recently conducted, a validated method is yet to be widely accepted. In this thesis, a least-invasive method for the damage-assessment of steel reinforcement is developed. Based on the information obtained from hardness testing and a single tensile test, it is possible to estimate the mechanical properties of earthquake-damaged rebars. The reduction in the low-cycle fatigue life due to strain ageing is also quantified. The proposed damage assessment methodology is based on empirical relationships between hardness and strain and residual strain capacity. If damage is suspected from in situ measurements, visual inspection or computer analysis, a bar may be removed and more accurate hardness measurements can be obtained using the lab-based Vickers hardness methodology. The Vickers hardness profile of damaged bars is then compared with calibration curves (Vickers hardness versus strain and residual strain capacity) previously developed for similar steel reinforcement bars extracted from undamaged locations. Experimental tests demonstrated that the time- and temperature-dependent strain-ageing phenomenon causes changes in the mechanical properties of plastically deformed steels. In particular, yield strength and hardness increases, whereas ductility decreases. The changes in mechanical properties are quantified and their implications on the hardness method are highlighted. Low-cycle fatigue (LCF) failures of steel reinforcing bars have been observed in laboratory testing and post-earthquake damage inspections. Often, failure might not occur during a first seismic event. However, damage is accumulated and the remaining fatigue life is reduced. Failure might therefore occur in a subsequent seismic event. Although numerous studies exist on the LCF behaviour of steel rebars, no studies had been conducted on the strain-ageing effects on the remaining fatigue life. In this thesis, the reduction in fatigue life due to this phenomenon is determined through a number of experimental tests.
This paper develops representative ground motion ensembles for several major earthquake scenarios in New Zealand. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch, and the occurrence of Wellington, Wairarapa, and Ohariu, fault ruptures in Wellington. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The ground motions are selected based on the generalized conditional intensity measure (GCIM) approach, ensuring that the ground motion ensembles represent both the mean, and distribution of ground motion intensity which such scenarios could impose. These scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
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.
The 2010 Darfield earthquake is the largest earthquake on record to have occurred within 40 km of a major city and not cause any fatalities. In this paper the authors have reflected on their experiences in Christchurch following the earthquake with a view to what worked, what didn’t, and what lessons can be learned from this for the benefit of Australian earthquake preparedness. Owing to the fact that most of the observed building damage occurred in Unreinforced Masonry (URM) construction, this paper focuses in particular on the authors’ experience conducting rapid building damage assessment during the first 72 hours following the earthquake and more detailed examination of the performance of unreinforced masonry buildings with and without seismic retrofit interventions.
The earthquakes that struck Ōtautahi/Christchurch began September 2010 and continued throughout2012 with the worse shock being February 22, 2011. The extended ‘seismic event’ radically altered thegeophysical and socio-cultural environments of the city. This working paper presents a broad array of datadescribing the impacts of the disaster on Māori. These data frame the results of small email surveyconducted 18 months after the most destructive February 22, 2011. This survey followed two projectsinvestigating the resilience of Māori to the disaster (Lambert & Mark-Shadbolt, 2011; Lambert & Mark-Shadbolt, 2012; Lambert, Mark-Shadbolt, Ataria, & Black, 2012). Results show that while the termresilience has become common to the point of cliché, the Māori experience thus far is best described asendurance.
The cartoon shows members of the 'Japan Seismic Institute studying an earthquake graph; one of them says 'Where was that Kiwi moonman when we wanted him?' Depicted also is thew 'moonman' Ken Ring as a wizard studying an astrological chart. Behind them all Japan is depicted as a devastated wasteland. Context - the Japanese earthquake and tsunami of 4th March 2011 and the present threat of a nuclear catastrophe. Also the so-called Moon Man, astrologer Ken Ring, who predicted that Christchurch would be hit by a huge earthquake today (20 March 2011). His claims have terrified Cantabrians and led to people fleeing Christchurch. Quantity: 1 digital cartoon(s).
© 2017 The Royal Society of New Zealand. This paper discusses simulated ground motion intensity, and its underlying modelling assumptions, for great earthquakes on the Alpine Fault. The simulations utilise the latest understanding of wave propagation physics, kinematic earthquake rupture descriptions and the three-dimensional nature of the Earth's crust in the South Island of New Zealand. The effect of hypocentre location is explicitly examined, which is found to lead to significant differences in ground motion intensities (quantified in the form of peak ground velocity, PGV) over the northern half and southwest of the South Island. Comparison with previously adopted empirical ground motion models also illustrates that the simulations, which explicitly model rupture directivity and basin-generated surface waves, lead to notably larger PGV amplitudes than the empirical predictions in the northern half of the South Island and Canterbury. The simulations performed in this paper have been adopted, as one possible ground motion prediction, in the ‘Project AF8’ Civil Defence Emergency Management exercise scenario. The similarity of the modelled ground motion features with those observed in recent worldwide earthquakes as well as similar simulations in other regions, and the notably higher simulated amplitudes than those from empirical predictions, may warrant a re-examination of regional impact assessments for major Alpine Fault earthquakes.
A Line of Best Fit explores weakness and disconnection in the city. Weakness: There are over 600 earthquake prone buildings in Wellington. The urgency to strengthen buildings risks compromising the aesthetic integrity of the city through abrasive strengthening techniques, or losing a large portion of our built environment to demolition. The need for extensive earthquake strengthening in Wellington, Christchurch and other New Zealand cities provides an exciting opportunity for architecture. Disconnection: In Wellington pedestrian activity is focused around three main routes: Cuba Street, Lambton Quay and Courtney Place. The adjacent areas are often disconnected and lack vibrancy due to large building footprints, no-exit laneways and lack of public spaces. The Design proposes a strategy for earthquake strengthening, preserving and upgrading the built environment, and expanding and connecting the pedestrian realm. The site is two earthquake prone buildings on the block between Marion Street and Taranaki Street in central Wellington. A cut through the centre of the Aspro and Cathie Buildings ties the buildings together to strengthen and create a new arcade as public space. The cut aligns with existing pedestrian routes connecting the block with the city. The Design is divided into three components: Void, Curve, and Pattern and Structure. Void investigates the implications of cutting a portion out the existing buildings and the opportunities this provides for connection, urban interaction, and light. Curve discusses the unusual form of The Design in terms of scale, the human response and the surrounding spaces. Pattern and Structure considers the structural requirements of the project and how a void enveloped in perforated screens can strengthen the earthquake prone buildings. The importance of connection, providing strength in the city, a dialogue between old and new, and engagement with the unexpected are evaluated. Opportunities for further development and research are discussed, with particular reference to how the principles of The Design could be implemented on a larger scale throughout our cities. A Line of Best Fit is an architectural proposal that creates strength and connection.
Text reads 'The earth moved...... and so did the government's infrastructure spending estimate'. Centre cartoon is the text '$17 billion' set against the backdrop of a seismic graph. Context - Auckland and Christchurch have been given top priority in the Government's latest national infrastructure plan, with more than seven billion dollars of its $17 billion budget going to the two cities. NZ Council for Infrastructure chief executive Stephen Selwood says hes not surprised at the large amount of funding going to Auckland and Christchurch, saying its clear the Christchurch rebuild will require a major commitment and Aucklands continued growth also requires significant funding. (Source: www.3news.co.nz, 5 July 2011) Quantity: 1 digital cartoon(s).
In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.
The progressive damage and subsequent demolition of unreinforced masonry (URM) buildings arising from the Canterbury earthquake sequence is reported. A dataset was compiled of all URM buildings located within the Christchurch CBD, including information on location, building characteristics, and damage levels after each major earthquake in this sequence. A general description of the overall damage and the hazard to both building occupants and to nearby pedestrians due to debris falling from URM buildings is presented with several case study buildings used to describe the accumulation of damage over the earthquake sequence. The benefit of seismic improvement techniques that had been installed to URM buildings is shown by the reduced damage ratios reported for increased levels of retrofit. Demolition statistics for URM buildings in the Christchurch CBD are also reported and discussed. VoR - Version of Record
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.
Floor systems with precast concrete hollow-core units have been largely used in concrete buildings built in New Zealand during the 1980’s. Recent earthquakes, such as the Canterbury sequence in 2010-2011 and the Kaikoura earthquake in 2016, highlighted that this floor system can be highly vulnerable and potentially lead to the floor collapse. A series of research activities are in progress to better understand the seismic performance of floor diaphragms, and this research focuses on examining the performance of hollow core units running parallel to the walls of wall-resisting concrete structures. This study first focused on the development of fragility functions, which can be quickly used to assess likelihood of the hollow-core being able to survive given the buildings design drift, and secondly to determine the expected performance of hollow-core units that run parallel to walls, focusing on the alpha unit running by the wall. Fragility functions are created for a range of different parameters for both vertical dislocation and crack width that can be used as the basis of a quick analysis or loss estimation for the likely impact of hollow-core floors on building vulnerability and risk. This was done using past experimental tests, and the recorded damage. Using these results and the method developed by Baker fragility curves were able to be created for varying crack widths and vertical dislocations. Current guidelines for analysis of hollow-core unit incompatible displacements are based on experimental vertical displacement results from concrete moment resisting frame systems to determine the capacity of hollow-core elements. To investigate the demands on hollow-core units in a wall-based structure, a fibre-element model in the software Seismostruct is created and subject to quasi-static cyclic loading, using elements which are verified from previous experimental tests. It is shown that for hollow-core units running by walls that the 10 mm displacement capacity used for hollow-core units running by a beam is insufficient for members running by walls and that shear analysis should be used. The fibre-element model is used to simulate the seismic demand induced on the floor system and has shown that the shear demand is a function of drift, wall length, hollow-core span, linking slab length and, to a minor extent, wall elongation.
A major lesson from the 2011 Christchurch earthquake was the apparent lack of ductility of some lightly reinforced concrete (RC) wall structures. In particular, the structural behaviour of the critical wall in the Gallery Apartments building demonstrated that the inelastic deformation capacity of a structure, as well as potentially brittle failure of the reinforcement, is dependent on the level of bond deterioration between reinforcement and surrounding concrete that occurs under seismic loading. This paper presents the findings of an experimental study on bond behaviour between deformed reinforcing bars and the surrounding concrete. Bond strength and relative bond slip was evaluated using 75 pull-out tests under monotonic and cyclic loading. Variations of the experiments include the loading rate, loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and the position of the embedded bond region within the specimen (deep within or close to free surface). Select test results are presented with inferred implications for RC structures.
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.
High rise developments dominate skylines and are contentious in many low rise urban environments. Christchurch is no exception and its residents have historically been vocal in articulating their opinions on matters they care about, especially in regard to projects they perceive will ruin their ‘garden city’. At the turn of the millennium, developers were preparing yet another proposal which would get the tongues wagging in Christchurch with the development of the former Ferrymead Tavern site on Ferry Road. The planning process was a long and antagonistic one with many individuals viewing the built towers with a look of ‘disgust’ and discontent. In an ironic twist, the seismic activity in Christchurch over the last few years which has had major implications for a range of planning issues, incrementally led to the death of highly controversial Ferrymead ‘Water’s Edge’ Apartments.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
The Canterbury earthquakes of 2010 and 2011 have shone the spotlight on a number of tax issues. These issues, and in particular lessons learned from them, will be relevant for revenue authorities, policymakers and taxpayers alike in the broader context of natural disasters. Issues considered by this paper include the tax treatment of insurance monies. For example, building owners will receive pay-outs for destroyed assets and buildings which have been depreciated. Where the insurance payment is more than the adjusted tax value, there will be a taxable "gain on sale" (or depreciation recovery income). If the building owner uses those insurance proceeds to purchase a replacement asset, legislative amendments specifically enacted following the earthquakes provide that rollover relief of the depreciation recovery income is available. The tax treatment of expenditure to seismically strengthen a building is another significant issue faced by building owners. Case law has determined that this expenditure will usually be capital expenditure. In the past such costs could be capitalised to the building and depreciated accordingly. However, since the 2011-2012 income year owners have been prohibited from claiming depreciation on buildings and therefore currently no deduction is available for such strengthening expenditure (whether immediate or deferred). This has significant potential implications for landlords throughout New Zealand facing significant seismic retrofit costs. Incentives, or some form of financial support, whether delivered through the tax system or some other mechanism may be required. International Financial Reporting Standards (IFRS) require insurance proceeds, including reimbursement for expenditure of a capital nature, be reported as income while expenditure itself is not recorded as a current period expense. This has the effect of overstating current income and creating a larger variation between reported income for accounting and taxation purposes. Businesses have obligations to maintain certain business records for tax purposes. Reconstructing records destroyed by a natural disaster depends on how the information was originally stored. The earthquakes have demonstrated the benefits of ‘off-site’ (outside Canterbury) storage, in particular electronic storage. This paper considers these issues and the Inland Revenue Department (Inland Revenue) Standard Practice Statement which deals with inter alia retention of business records in electronic format and offshore record storage. Employer provided accommodation is treated as income to the benefitting employee. A recent amendment to the Income Tax Act 2007 retrospectively provides that certain employer provided accommodation is exempt from tax. The time aspect of these rules is extended where the employee is involved in the Canterbury rebuild and comes from outside the region.
To reduce seismic vulnerability and the economic impact of seismic structural damage, it is important to protect structures using supplemental energy dissipation devices. Several types of supplemental damping systems can limit loads transferred to structures and absorb significant response energy without sacrificial structural damage. Lead extrusion dampers are one type of supplemental energy dissipation devices. A smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, have been employed in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch and San Francisco. HF2V devices have previously been designed using very simple models with limited precision. They are then manufactured, and tested to ensure force capacities match design goals, potentially necessitating reassembly or redesign if there is large error. In particular, devices with a force capacity well above or below a design range can require more testing and redesign, leading to increased economic and time cost. Thus, there is a major need for a modelling methodology to accurately estimate the range of possible device force capacity values in the design phase – upper and lower bounds. Upper and lower bound force capacity estimates are developed from equations in the metal extrusion literature. These equations consider both friction and extrusion forces between the lead and the bulged shaft in HF2V devices. The equations for the lower and upper bounds are strictly functions of device design parameters ensuring easy use in the design phase. Two different sets of estimates are created, leading to estimates for the lower and upper bounds denoted FLB,1, FUB,1, FUB,2, respectively. The models are validated by comparing the bounds with experimental force capacity data from 15 experimental HF2V device tests. All lower bound estimates are below or almost equal to the experimental device forces, and all upper bound estimates are above. Per the derivation, the (FLB,1, FUB,1) pair provide narrower bounds. The (FLB,1, FUB,1) pair also had a mean lower bound gap of -34%, meaning the lower bound was 74% of device force on average, while the mean upper bound gap for FUB,1 was +23%. These are relatively tight bounds, within ~±2 SE of device manufacture, and can be used as a guide to ensure device forces are in range for the actual design use when manufactured. Therefore, they provide a useful design tool.
This paper provides a summary of the ground motions observed in the recent Canterbury, New Zealand earthquake sequence. The sequence occurred in a region of relatively moderate seismicity, 130km to the east of the Alpine Fault, the major plate-boundary in the region. From an engineering perspective, the sequence has been primarily comprised of the initial 04/09/2010 Darfield earthquake (Mw7.1) followed by the 22/02/2011 Christchurch earthquake (Mw6.3), and two aftershocks on 13/06/ 2011 (Mw5.3 and 6.0, respectively). The dense spacing of strong motions in the region, and their close proximity to the respective causative faults, has resulted in strong ground motions far exceeding the previous catalogue of strong motion observed in New Zealand. The observed ground motions have exhibited clear evidence of: (i) near-source directivity; (ii) sedimentary basin focusing, amplification and basin effect refraction; (iii) non-linear site response; (iv) cyclic mobility postliquefaction; and (v) extreme vertical ground motions exceeding 2g, among others.
This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.
Shaking table testing of a full-scale three storey resilient and reparable complete composite steel framed building system is being conducted. The building incorporates a number of interchangeable seismic resisting systems of New Zealand and Chinese origin. The building has a steel frame and cold formed steel-concrete composite deck. Energy is dissipated by means of friction connections. These connections are arranged in a number of structural configurations. Typical building nonskeletal elements (NSEs) are also included. Testing is performed on the Jiading Campus shaking table at Tongji University, Shanghai, China. This RObust BUilding SysTem (ROBUST) project is a collaborative China-New Zealand project sponsored by the International Joint Research Laboratory of Earthquake Engineering (ILEE), Tongji University, and a number of agencies and universities within New Zealand including BRANZ, Comflor, Earthquake Commission, HERA, QuakeCoRE, QuakeCentre, University of Auckland, and the University of Canterbury. This paper provides a general overview of the project describing a number of issues encountered in the planning of this programme including issues related to international collaboration, the test plan, and technical issues.
Following the 2010/2011 Canterbury earthquakes the seismic design of buildings with precast concrete panels has received significant attention. Although this form of construction generally performed adequately in Christchurch, there were a considerable number of precast concrete panel connection failures. This observation prompted a review of more than 4700 panel details from 108 buildings to establish representative details used in both existing and new multi-storey and low rise industrial precast concrete buildings in three major New Zealand cities of Auckland, Wellington and Christchurch. Details were collected from precast manufacturers and city councils and were categorised according to type. The detailing and quantity of each reviewed connection type in the sampled data is reported, and advantages and potential deficiencies of each connection type are discussed. The results of this survey provide a better understanding of the relative prevalence of common detailing used in precast concrete panels and guidance for the design of future experimental studies. http://www.nzsee.org.nz/publications/nzsee-quarterly-bulletin/