The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
Quick and reliable assessment of the condition of bridges in a transportation network after an earthquake can greatly assist immediate post-disaster response and long-term recovery. However, experience shows that available resources, such as qualified inspectors and engineers, will typically be stretched for such tasks. Structural health monitoring (SHM) systems can therefore make a real difference in this context. SHM, however, needs to be deployed in a strategic manner and integrated into the overall disaster response plans and actions to maximize its benefits. This study presents, in its first part, a framework of how this can be achieved. Since it will not be feasible, or indeed necessary, to use SHM on every bridge, it is necessary to prioritize bridges within individual networks for SHM deployment. A methodology for such prioritization based on structural and geotechnical seismic risks affecting bridges and their importance within a network is proposed in the second part. An example using the methodology application to selected bridges in the medium-sized transportation network of Wellington, New Zealand is provided. The third part of the paper is concerned with using monitoring data for quick assessment of bridge condition and damage after an earthquake. Depending on the bridge risk profile, it is envisaged that data will be obtained from either local or national seismic monitoring arrays or SHM systems installed on bridges. A method using artificial neural networks is proposed for using data from a seismic array to infer key ground motion parameters at an arbitrary bridges site. The methodology is applied to seismic data collected in Christchurch, New Zealand. Finally, how such ground motion parameters can be used in bridge damage and condition assessment is outlined AM - Accepted manuscript
Earthquakes cause significant damage to buildings due to strong vibration of the ground. Levitating houses using magnets and electromagnets would provide a complete isolation of ground motion for protecting buildings from seismic damage. Two types of initial configuration for the electromagnet system were proposed with the same air gap (10mm) between the electromagnet and reluctance plate. Both active and passive controller are modelled to investigate the feasibility of using a vibration control system for stabilizing the magnetic system within the designed air gap (10mm) in the vertical direction. A nonlinear model for the magnetic system is derived to implement numerical simulation of structural response under the earthquake record in Christchurch Botanic Gardens on 21 February 2011. The performance of the uncontrolled and the controlled systems are compared and the optimal combination of control gains are determined for the PID active controller. Simulation results show both active PID controller with constant and nonlinear attracting force are able to provide an effective displacement control within the required air gap (+/-5mm). The maximum control force demand for the PID controller in the presence of nonlinear attracting force is 4.1kN, while the attracting force in equilibrium position is 10kN provided by the electromagnet. These results show the feasibility of levitating a house using the current electromagnet and PID controller. Finally, initial results of passive control using two permanent magnets or dampers show the structural responses can be effectively reduced and centralized to +/-1mm using a nonlinear centring barrier function.
Rapid, accurate structural health monitoring (SHM) assesses damage to optimise decision-making. Many SHM methods are designed to track nonlinear stiffness changes as damage. However, highly nonlinear pinched hysteretic systems are problematic in SHM. Model-based SHM often fails as any mismatch between model and measured response dynamics leads to significant error. Thus, modelfree methods of hysteresis loop tracking methods have emerged. This study compares the robustness and accuracy in the presence of significant measurement noise of the proven hysteresis loop analysis (HLA) SHM method with 3 emerging model-free methods and 2 further novel adaptations of these methods using a highly nonlinear, 6-story numerical structure to provide a known ground-truth. Mean absolute errors in identifying a known nonlinear stiffness trajectory assessed at four points over two successive ground motion inputs from September 2010 and February 2011 in Christchurch range from 1.71-10.52%. However, the variability is far wider with maximum errors ranging from 3.90-49.72%, where the second largest maximum absolute error was still 19.74%. The lowest mean and maximum absolute errors were for the HLA method. The next best method had mean absolute error of 2.92% and a maximum of 10.51%. These results show the clear superiority of the HLA method over all current emerging model-free methods designed to manage the highly nonlinear pinching responses common in reinforced concrete structures. These results, combined with high robustness and accuracy in scaled and fullscale experimental studies, provide further validation for using HLA for practical implementation.
This thesis is concerned with springs that appeared in the Hillsborough, Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and which have continued to discharge groundwater to the surface to the present time. Investigations have evolved, measurements of discharge at selected sites, limited chemical data on anions and isotope analysis. The springs are associated with earthquake generated fissures (extensional) and compression zones, mostly in loess-colluvium soils of the valley floor and lower slopes. Extensive peat swamps are present in the Hillsborough valley, with a groundwater table at ~1m below ground. The first appearance of the ‘new’ springs took place following the Mw 7.1 Darfield Earthquake on 4 September 2010, and discharges increased both in volume and extent of the Christchurch Mw 6.3 Earthquake of 22 February 2011. Five monitored sites show flow rates in the range of 4.2-14.4L/min, which have remained effectively constant for the duration of the study (2014-2015). Water chemistry analysis shows that the groundwater discharges are sourced primarily from volcanic bedrocks which underlies the valley at depths ≤50m below ground level. Isotope values confirm similarities with bedrock-sourced groundwater, and the short term (hours-days) influence of extreme rainfall events. Cyclone Lusi (2013-2014) affects were monitored and showed recovery of the bedrock derived water signature within 72 hours. Close to the mouth of the valley sediments interfinger with Waimakiriri River derived alluvium bearing a distinct and different isotope signature. Some mixing is evident at certain locations, but it is not clear if there is any influence from the Huntsbury reservoir which failed in the Port Hills Earthquake (22 February 2011) and stored groundwater from the Christchurch artesian aquifer system (Riccarton Gravel).
This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized
Motivation This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized.
The Canterbury earthquakes caused huge amounts of damage to Christchurch and the surrounding area and presented a very challenging situation for both insurers and claimants. While tourism has suffered significant losses as a result, particularly due to the subsequent decrease in visitor numbers, the Canterbury region was very fortunate to have high levels of insurance coverage. This report, based on data gathered from tourism operators on the ground in Canterbury, looks at how this sector has been affected by the quakes, claims patterns, and the behaviour and perceptions of tourism operators about insurance.
This paper summarizes the development of a high-resolution 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. 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, as well as use in site response analysis and ground motion simulation.
At 00:02 on 14 November, 2016 a destructive 7.8 Mw earthquake struck the North Canterbury region of New Zealand’s South Island. Prior to and following the earthquake, natural and social scientists conducted a significant amount of research on the resilience processes and recovery efforts in North Canterbury. This thesis examines community resilience in Kaikōura, a small town and district greatly impacted by the earthquake. Community resilience has been widely used in disaster risk reduction research, policy, and practice to describe how a group of individuals within a boundary respond to events, hazards, and shifts in their everyday life. Using exploratory inquiry, this thesis adopts qualitative research methods including document analysis, 24 semi-structured interviews, and participant observation to explore the idea that the recent scholarly emphasis on resilience has come at the expense of critical engagement with the complexities of communities. I draw on the idea of ‘collectives’ (comprising community-based organisations or less formal social networks with a shared purpose) as a lens to consider how, when unexpected life events happen, collectives can be regarded as a resource for change or constancy. The examination of collectives following a disaster can lend insight into the many elements of community as they bring people together in collaboration or drive them apart in conflict. This thesis therefore contributes to an enhanced practical and theoretical understanding of both community and resilience.
The 4th of September 2010 Mw 7.1 Darfield (Canterbury) earthquake had generated significant ground shaking within the Christchurch Central Business District (CBD). Despite the apparently significant shaking, the observed structural damage for pre-1970s reinforced concrete (RC) buildings was indeed limited and lower than what was expected for such typology of buildings. This paper explores analytically and qualitatively the different aspects of the "apparent‟ good seismic performance of the pre-1970s RC buildings in the Christchurch CBD, following the earthquake reconnaissance survey by the authors. Damage and building parameters survey result, based on a previously established inventory of building stock of these non-ductile RC buildings, is briefly reported. From an inventory of 75 buildings, one building was selected as a numerical case-study to correlate the observed damage with the non-linear analyses. The result shows that the pre-1970s RC frame buildings performed as expected given the intensity of the ground motion shaking during the Canterbury earthquake. Given the brittle nature of this type of structure, it was demonstrated that more significant structural damage and higher probability of collapse could occur when the buildings were subjected to alternative input signals with different frequency content and duration characteristics and still compatible to the seismicity hazard for Christchurch CBD.
The initial goal of this research was to explore how SME business models change in response to a crisis. Keeping this in mind, the business model canvas (Osterwalder & Pigneur, 2010) was used as a tool to analyse SME business models in the Canterbury region of New Zealand. The purpose was to evaluate the changes SMEs instituted in their business models after being hit by a series of earthquakes in 2010 and 2011. The idea was to conduct interviews with business owners and analyse them using grounded theory methods. As this method is iterative and requires simultaneous data collection and analysis, a tentative model was proposed after first phase of the data collection and analysis. However, as a result of this process, it became apparent that owner-specific characteristics, action orientation and networks were more prominent in the data than business model elements. Although the SMEs in this study experienced several operational changes in their business models, such as a change of location, modifications to their payment terms or expanded/restricted target markets, the suggested framework highlights how owner-specific attributes ensured the recovery of their businesses. After the initial framework was suggested, subsequent interviews were conducted to test, verify, and modify the tentative model. Three aspects of business recovery emerged: (a) cognitive coping – the business owner’s mind-set and motive; (b) adaptive coping – the ability of business owner to take corrective actions; and (c) social capital – the social network of a business owner, including formal and informal connections and their significance. Three distinct groups were identified; self-sufficient SMEs, socially-based SMEs and surviving SMEs. This thesis proposes a grounded theory of business recovery for SMEs following a disaster. Cognitive coping and social capital enabled the owners to take actions, which eventually led to the desired outcomes for the businesses.
Pumice materials, which are problematic from an engineering viewpoint, are widespread in the central part of the North Island. Considering the impacts of the 2010-2011 Christchurch earthquakes, a clear understanding of their properties under earthquake loading is necessary. For example, the 1987 Edgecumbe earthquake showed evidence of localised liquefaction of sands of volcanic origin. To elucidate on this, research was undertaken to investigate whether existing empirical field-based methods to evaluate the liquefaction potential of sands, which were originally developed for hard-grained soils, are applicable to crushable pumice-rich deposits. For this purpose, two sites, one in Whakatane and another in Edgecumbe, were selected where the occurrence of liquefaction was reported following the Edgecumbe earthquake. Manifestations of soil liquefaction, such as sand boils and ejected materials, have been reported at both sites. Field tests, including cone penetration tests (CPT), shear-wave velocity profiling, and screw driving sounding (SDS) tests were performed at the sites. Then, considering estimated peak ground accelerations (PGAs) at the sites based on recorded motions and possible range of ground water table locations, liquefaction analysis was conducted at the sites using available empirical approaches. To clarify the results of the analysis, undisturbed soil samples were obtained at both sites to investigate the laboratory-derived cyclic resistance ratios and to compare with the field-estimated values. Research results clearly showed that these pumice-rich soils do not fit existing liquefaction assessment frameworks and alternate methods are necessary to characterise them.
This thesis presents an assessment of historic seismic performance of the New Zealand stopbank network from the 1968 Inangahua earthquake through to the 2016 Kaikōura earthquake. An overview of the types of stopbanks and the main aspects of the design and construction of earthen stopbanks was presented. Stopbanks are structures that are widely used on the banks of rivers and other water bodies to protect against the impact of flood events. Earthen stopbanks are found to be the most used for such protection measures. Different stopbank damage or failure modes that may occur due to flooding or earthquake excitation were assessed with a focus on past earthquakes internationally, and examples of these damage and failure modes were presented. Stopbank damage and assessment reports were collated from available reconnaissance literature to develop the first geospatial database of stopbank damage observed in past earthquakes in New Zealand. Damage was observed in four earthquakes over the past 50 years, with a number of earthquakes resulting in no stopbank damage. The damage database therefore focussed on the Edgecumbe, Darfield, Christchurch and Kaikōura earthquakes. Cracking of the crest and liquefaction-induced settlement were the most common forms of damage observed. To understand the seismic demand on the stopbank network in past earthquakes, geospatial analyses were undertaken to approximate the peak ground acceleration (PGA) across the stopbank network for ten large earthquakes that have occurred in New Zealand over the past 50 years. The relationship between the demand, represented by the peak ground acceleration (PGA) and damage is discussed and key trends identified. Comparison of the seismic demand and the distribution of damage suggested that the seismic performance of the New Zealand stopbank network has been generally good across all events considered. Although a significant length of the stopbank networks were exposed to high levels of shaking in past events, the overall damage length was a small percentage of this. The key aspect controlling performance was the performance of the underlying foundation soils and the effect of this on the stopbank structure and stability
This paper presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS. In this regard, CPT, SPT and borehole data, and shear wave velocity (Vs) profiles in close vicinity to the SMS are currently being used to develop representative soil profiles at each site. Site subsoil classifications based on Vs measurements performed by the authors do not always agree with the original classifications, often indicating that a softer site class is appropriate. However, SPT N values often indicate a stiffer site class than the Vs data, in some cases also disagreeing with prior assumed classifications. Hence, the recent site investigation data presented herein highlights the importance of having detailed site-specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N and Vs-.
This paper presents a seismic velocity model of Canterbury, New Zealand based on 3D geologic surfaces and velocities from a range of data sources. The model provides the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, such as broadband ground motion and shallow site response analyses related to understanding the ground motions and site responses during the 2010- 2011 Canterbury earthquakes. Pre-Quaternary geologic horizons are calculated based on the reinterpretation of a comprehensive network of seismic reflection surveys from seven different campaigns over the past 50 years, as well as point constraints across an array of petroleum industry drill holes. Particular attention is given to a detailed representation of Quaternary stratigraphy, representing shallow (z<250m) near-surface layers in the model. Seismic velocities are obtained from seismic reflection processing (for Vp) and also recently performed active and passive surface wave analyses (for Vs). Over 1,700 water wells in the region are used to constrain the complex inter-bedded Quaternary stratigraphy (gravels, sands, silts, organics etc.) near the coastline, including beneath urban Christchurch, which has resulted from fluvial deposition and marine regression and transgression. For the near-surface Springston and Christchurch Formations in the Christchurch urban area (z<50m), high-spatial resolution seismic velocities (including Vs30 ) were obtained from over 13,000 cone penetration tests combined with a recently developed CPT-Vs correlation.
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
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.
Liquefaction-induced lateral spreading during earthquakes poses a significant hazard to the built environment, as observed in Christchurch during the 2010 to 2011 Canterbury Earthquake Sequence (CES). It is critical that geotechnical earthquake engineers are able to adequately predict both the spatial extent of lateral spreads and magnitudes of associated ground movements for design purposes. Published empirical and semi-empirical models for predicting lateral spread displacements have been shown to vary by a factor of <0.5 to >2 from those measured in parts of Christchurch during CES. Comprehensive post- CES lateral spreading studies have clearly indicated that the spatial distribution of the horizontal displacements and extent of lateral spreading along the Avon River in eastern Christchurch were strongly influenced by geologic, stratigraphic and topographic features.
On February 22, 2011, a magnitude Mw 6.2 earthquake affected the Canterbury region, New Zealand, resulting in many fatalities. Liquefaction occurred across many areas, visible on the surface as ‘‘sand volcanoes’’, blisters and subsidence, causing significant damage to buildings, land and infrastructure. Liquefaction occurred at a number of sites across the Christchurch Boys High School sports grounds; one area in particular contained a piston ground failure and an adjacent silt volcano. Here, as part of a class project, we apply near-surface geophysics to image these two liquefaction features and determine whether they share a subsurface connection. Hand auger results enable correlation of the geophysical responses with the subsurface stratigraphy. The survey results suggest that there is a subsurface link, likely via a paleo-stream channel. The anomalous responses of the horizontal loop electromagnetic survey and electrical resistivity imaging highlight the disruption of the subsurface electrical properties beneath and between the two liquefaction features. The vertical magnetic gradient may also show a subtle anomalous response in this area, however the results are inconclusive. The ground penetrating radar survey shows disruption of the subsurface stratigraphy beneath the liquefaction features, in particular sediment mounding beneath the silt ejection (‘‘silt volcano’’) and stratigraphic disruption beneath the piston failure. The results indicate how near-surface geophysics allow the characteristics of liquefaction in the subsurface to be better understood, which could aid remediation work following liquefaction-induced land damage and guide interpretation of geophysical surveys of paleoliquefaction features.
Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.
The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.
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.
Natural disasters are increasingly disruptive events that affect livelihoods, organisations, and economies worldwide. Research has identified the impacts and responses of organisations to different types of natural disasters, and have outlined factors, such as industry sector, that are important to organisational vulnerability and resilience. One of the most costly types of natural disasters in recent years has been earthquakes, and yet to date, the majority of studies have focussed on the effects of earthquakes in urban areas, while rural organisational impact studies have primarily focused on the effects of meteorological and climatic driven hazards. As a result, the likely impacts of an earthquake on rural organisations in a developed context is unconstrained in the literature. In countries like New Zealand, which have major earthquakes and agricultural sectors that are significant contributors to the economy, it is important to know what impacts an earthquake event would have on the rural industries, and how these impacts compare to that of a more commonly analysed, high-frequency event. In September of 2010, rural organisations in Canterbury experienced the 4 September 2010 Mw 7.1 `Darfield' earthquake and the associated aftershocks, which came to be known as the Canterbury earth- quake sequence. The earthquake sequence caused intense ground shaking, creating widespread critical service outages, structural and non-structural damage to built infrastructure, as well as ground surface damage from ooding, liquefaction and surface rupture. Concurrently on September 18 2010, rural organisations in Southland experienced an unseasonably late snowstorm and cold weather snap that brought prolonged sub-zero temperatures, high winds and freezing rain, damaging structures in the City of Invercargill and causing widespread livestock losses and production decreases across the region. This thesis documents the effects of the Canterbury earthquake sequence and Southland snowstorm on farming and rural non-farming organisations, utilizing comparable methodologies to analyse rural organisational impacts, responses and recovery strategies to natural disasters. From the results, a short- term impact assessment methodology is developed for multiple disasters. Additionally, a regional asset repair cost estimation model is proposed for farming organisations following a major earthquake event, and the use of social capital in rural organisational recovery strategies following natural disasters is analysed.
In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.
The 2010-2011 Canterbury earthquake sequence, and the resulting extensive data sets on damaged buildings that have been collected, provide a unique opportunity to exercise and evaluate previously published seismic performance assessment procedures. This poster provides an overview of the authors’ methodology to perform evaluations with two such assessment procedures, namely the P-58 guidelines and the REDi Rating System. P-58, produced by the Federal Emergency Management Agency (FEMA) in the United States, aims to facilitate risk assessment and decision-making by quantifying earthquake ground shaking, structural demands, component damage and resulting consequences in a logical framework. The REDi framework, developed by the engineering firm ARUP, aids stakeholders in implementing resilience-based earthquake design. Preliminary results from the evaluations are presented. These have the potential to provide insights on the ability of the assessment procedures to predict impacts using “real-world” data. However, further work remains to critically analyse these results and to broaden the scope of buildings studied and of impacts predicted.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
There is a relationship between inelastic deformation and energy dissipation in structures that are subjected to earthquake ground motions. Thus, if seismic energy dissipation can be achieved by means of a separate non-load bearing supplementary damping system, the load bearing structure can remain elastic with continuing serviceability following the design level earthquake. This research was carried out to investigate the advantages of using added damping in structures. The control system consists of passive friction dampers called ring spring dampers installed in the ground floor of the structure using a tendon to transmit the forces to the other parts of the structure. The ring springs dampers are friction devices consisting of inner and outer ring elements assembled to form a spring stack. External load applied to the spring produces sliding action across mating ring interfaces. The damping forces generated by the dampers and transferred in the supplemental system to the structure by the tendon and horizontal links oppose the internal loads. A four storey-two bay steel frame structure was used in the study. Experimental and analytical studies to investigate the effectiveness of a supplemental control system are presented. The model was subjected to a series of earthquake simulations on the shaking table in the Structural Laboratory of the Civil Engineering Department, at the University of Canterbury. The earthquake simulation tests have been performed on the structure both with and without the supplemental control system. The earthquake simulations were a series of gradually increasing intensity replications of two commonly used earthquake records. This thesis includes detailed description of the structural model, the supplemental control system, the ring springs dampers and the data obtained during the testing. Analyses were then carried out on a twelve storey framed structure to investigate the possible tendon arrangements and the size and type of dampers required to control the response of a real building. Guidelines for determining the appropriateness of including a supplemental damping system have been investigated. The main features of the supplemental control system adopted in this research are: • It is a passive control system with extreme reliability and having no dependence on external power sources to effect the control action. These power sources may not be available during a major earthquake. • Ring springs are steel friction devices capable of absorbing large amounts of input energy. No liquid leakage can occur and minimal maintenance is required for the ring spring dampers. • With a damper-tendon system, the distribution of the dampers throughout the structure is not so critical. Only one or two dampers are used to produce the damping forces needed, and forces are then transferred to the rest of the building by the tendon system. • It is a relatively inexpensive control system with a long useful life.
Case study analysis of the 2010-2011 Canterbury Earthquake Sequence (CES), which particularly impacted Christchurch City, New Zealand, has highlighted the value of practical, standardised and coordinated post-earthquake geotechnical response guidelines for earthquake-induced landslides in urban areas. The 22nd February 2011 earthquake, the second largest magnitude event in the CES, initiated a series of rockfall, cliff collapse and loess failures around the Port Hills which severely impacted the south-eastern part of Christchurch. The extensive slope failure induced by the 22nd February 200 earthquake was unprecedented; and ground motions experienced significantly exceeded the probabilistic seismic hazard model for Canterbury. Earthquake-induced landslides initiated by the 22nd February 2011 earthquake posed risk to life safety, and caused widespread damage to dwellings and critical infrastructure. In the immediate aftermath of the 22nd February 2011 earthquake, the geotechnical community responded by deploying into the Port Hills to conduct assessment of slope failure hazards and life safety risk. Coordination within the voluntary geotechnical response group evolved rapidly within the first week post-earthquake. The lack of pre-event planning to guide coordinated geotechnical response hindered the execution of timely and transparent management of life safety risk from coseismic landslides in the initial week after the earthquake. Semi-structured interviews were conducted with municipal, management and operational organisations involved in the geotechnical response during the CES. Analysis of interview dialogue highlighted the temporal evolution of priorities and tasks during emergency response to coseismic slope failure, which was further developed into a phased conceptual model to inform future geotechnical response. Review of geotechnical responses to selected historical earthquakes (Northridge, 1994; Chi-Chi, 1999; Wenchuan, 2008) has enabled comparison between international practice and local response strategies, and has emphasised the value of pre-earthquake preparation, indicating the importance of integration of geotechnical response within national emergency management plans. Furthermore, analysis of the CES and international earthquakes has informed pragmatic recommendations for future response to coseismic slope failure. Recommendations for future response to earthquake-induced landslides presented in this thesis include: the integration of post-earthquake geotechnical response with national Civil Defence and Emergency Management; pre-earthquake development of an adaptive management structure and standard slope assessment format for geotechnical response; and emergency management training for geotechnical professionals. Post-earthquake response recommendations include the development of geographic sectors within the area impacted by coseismic slope failure, and the development of a GIS database for analysis and management of data collected during ground reconnaissance. Recommendations provided in this thesis aim to inform development of national guidelines for geotechnical response to earthquake-induced landslides in New Zealand, and prompt debate concerning international best practice.
The 2010–2011 Canterbury earthquake sequence began with the 4 September 2010, Mw7.1 Darfield earthquake and includes up to ten events that induced liquefaction. Most notably, widespread liquefaction was induced by the Darfield and Mw6.2 Christchurch earthquakes. The combination of well-documented liquefaction response during multiple events, densely recorded ground motions for the events, and detailed subsurface characterization provides an unprecedented opportunity to add well-documented case histories to the liquefaction database. This paper presents and applies 50 high-quality cone penetration test (CPT) liquefaction case histories to evaluate three commonly used, deterministic, CPT-based simplified liquefaction evaluation procedures. While all the procedures predicted the majority of the cases correctly, the procedure proposed by Idriss and Boulanger (2008) results in the lowest error index for the case histories analyzed, thus indicating better predictions of the observed liquefaction response.
