The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The aim of this multidisciplinary research was to retrospectively analyse the gastroenteritis prevalence following the February 22, 2011 earthquake in Christchurch. The first focus was to assess whether earthquake-induced infrastructure damage, liquefaction, and gastroenteritis agents spatially explained the recorded gastroenteritis cases over the period of 35 days following the February 22, 2011 earthquake in Christchurch. The gastroenteritis agents considered in this study were Escherichia coli found in the drinking water supply (MPN/100mL) and Non-Compliant Free Associated Chlorine (FAC-NC) (less than <0.02mg/L). The second focus was the protocols that averted a gastroenteritis outbreak at three Emergency Centres (ECs): Burnside High School Emergency Centre (BEC); Cowles Stadium Emergency Centre (CEC); and Linwood High School Emergency Centre (LEC). Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols and indirect themes included type of EC building (school or a sports stadium), and EC staff. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. This research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. A damage profile was created by amalgamating different types of damage for the considered factors for each Census Area Unit (CAU) in Christchurch. The damage profile enabled the application of a variety of statistical methods which included Moran’s I , Hot Spot (HS) analysis, Spearman’s Rho, and Besag–York–Mollié Model using a range of software. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. The ECs were selected to represent the Christchurch area, and were situated where potential for gastroenteritis was high. BEC represented the western side of Christchurch; whilst, CEC and LEC represented the eastern side, where the potential for gastroenteritis was high according to the outputs of the quantitative spatial modelling. Qualitative analysis from the interviews at the ECs revealed that evacuees were arriving at the ECs with gastroenteritis-like symptoms. Participants believed that those symptoms did not originate at the ECs. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols that included prolific use of hand sanitisers; surveillance; and the services offered. Indirect themes included the EC layout, type of EC building (school or a sports stadium), and EC staff. Indirect themes governed the quality and sustainability of the direct themes implemented, which in turn averted gastroenteritis outbreaks at the ECs. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. It was concluded that gastroenteritis point prevalence following the February 22, 2011 earthquake could not be solely explained by earthquake-induced infrastructure damage, liquefaction, and gastroenteritis causative agents alone. However, this research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Creating a damage profile for each CAU and using spatial data analysis can isolate vulnerable areas, and qualitative data analysis provides localised information. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally.
The 2010-2011 Canterbury earthquake sequence was extremely damaging to structures in Christchurch and continues to have a large economic and social impact on the city and surrounding regions. In addition to strong ground shaking (Bradley and Cubrinovski 2011 SRL; Bradley 2012 SDEE), extensive liquefaction was observed, particularly in the 4 September 2010 Darfield earthquake and the 22 February 2011 Christchurch earthquake (Cubrinovski et al. 2010 BNZSEE; 2011 SRL). Large observed vertical ground motion amplitudes were recorded in the events in this sequence, with vertical peak ground accelerations of over 2.2g being observed at the Heathcote Valley Primary School during the Christchurch earthquake, and numerous other vertical motions exceeding 1.0g (Bradley and Cubrinovski 2011 SRL; Bradley 2012 SDEE; Fry et al 2011 SRL). Vertical peak ground accelerations of over 1.2g were observed in the Darfield earthquake.
A magnitude 6.3 earthquake struck the city of Christchurch at 12:51pm on Tuesday 22 February 2011. The earthquake caused 182 fatalities, a large number of injuries, and resulted in widespread damage to the built environment, including significant disruption to the lifelines. The event created the largest lifeline disruption in a New Zealand city in 80 years, with much of the damage resulting from extensive and severe liquefaction in the Christchurch urban area. The Christchurch earthquake occurred when the Canterbury region and its lifelines systems were at the early stage of recovering from the 4 September 2010 Darfield (Canterbury) magnitude 7.1 earthquake. This paper describes the impact of the Christchurch earthquake on lifelines by briefly summarising the physical damage to the networks, the system performance and the operational response during the emergency management and the recovery phase. Special focus is given to the performance and management of the gas, electric and road networks and to the liquefaction ejecta clean-up operations that contributed to the rapid reinstatement of the functionality of many of the lifelines. The water and wastewater system performances are also summarized. Elements of resilience that contributed to good network performance or to efficient emergency and recovery management are highlighted in the paper.
The Mw 6.2 February 22nd 2011 Christchurch earthquake (and others in the 2010-2011 Canterbury sequence) provided a unique opportunity to study the devastating effects of earthquakes first-hand and learn from them for future engineering applications. All major events in the Canterbury earthquake sequence caused widespread liquefaction throughout Christchurch’s eastern suburbs, particularly extensive and severe during the February 22nd event. Along large stretches of the Avon River banks (and to a lesser extent along the Heathcote) significant lateral spreading occurred, affecting bridges and the infrastructure they support. The first stage of this research involved conducting detailed field reconnaissance to document liquefaction and lateral spreading-induced damage to several case study bridges along the Avon River. The case study bridges cover a range of ages and construction types but all are reinforced concrete structures which have relatively short, stiff decks. These factors combined led to a characteristic deformation mechanism involving deck-pinning and abutment back-rotation with consequent damage to the abutment piles and slumping of the approaches. The second stage of the research involved using pseudo-static analysis, a simplified seismic modelling tool, to analyse two of the bridges. An advantage of pseudo-static analysis over more complicated modelling methods is that it uses conventional geotechnical data in its inputs, such as SPT blowcount and CPT cone resistance and local friction. Pseudo-static analysis can also be applied without excessive computational power or specialised knowledge, yet it has been shown to capture the basic mechanisms of pile behaviour. Single pile and whole bridge models were constructed for each bridge, and both cyclic and lateral spreading phases of loading were investigated. Parametric studies were carried out which varied the values of key parameters to identify their influence on pile response, and computed displacements and damages were compared with observations made in the field. It was shown that pseudo-static analysis was able to capture the characteristic damage mechanisms observed in the field, however the treatment of key parameters affecting pile response is of primary importance. Recommendations were made concerning the treatment of these governing parameters controlling pile response. In this way the future application of pseudo-static analysis as a tool for analysing and designing bridge pile foundations in liquefying and laterally spreading soils is enhanced.
On 22 February 2011，a magnitude Mw 6.3 earthquake occurred with an epicenter located near Lyttelton at about 10km from Christchurch in Canterbury region on the South Island of New Zealand (Figure 1). Since this earthquake occurred in the midst of the aftershock activity which had continued since the 4 September 2010 Darfield Earthquake occurrence, it was considered to be an aftershock of the initial earthquake. Because of the short distance to the city and the shallower depth of the epicenter, this earthquake caused more significant damage to pipelines, traffic facilities, residential houses/properties and multi-story buildings in the central business district than the September 2010 Darfield Earthquake in spite of its smaller earthquake magnitude. Unfortunately, this earthquake resulted in significant number of casualties due to the collapse of multi-story buildings and unreinforced masonry structures in the city center of Christchurch. As of 4 April, 172 casualties were reported and the final death toll is expected to be 181. While it is extremely regrettable that Christchurch suffered a terrible number of victims, civil and geotechnical engineers have this hard-to-find opportunity to learn the response of real ground from two gigantic earthquakes which occurred in less than six months from each other. From geotechnical engineering point of view, it is interesting to discuss the widespread liquefaction in natural sediments, repeated liquefaction within short period and further damage to earth structures which have been damaged in the previous earthquake. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. The team included the following members: Misko Cubrinovski (University of Canterbury, NZ, Team Leader), Susumu Yasuda (Tokyo Denki University, Japan, JGS Team Leader), Rolando Orense (University of Auckland, NZ), Kohji Tokimatsu (Tokyo Institute of Technology, Japan), Ryosuke Uzuoka (Tokushima University, Japan), Takashi Kiyota (University of Tokyo, Japan), Yasuyo Hosono (Toyohashi University of Technology, Japan) and Suguru Yamada (University of Tokyo, Japan).
During the 2010 - 2011 Canterbury earthquake sequence, extensive liquefaction was observed in many areas of Christchurch city and its surroundings, causing widespread damage to buildings and infrastructure. While existing simplified methods were found to work well in some areas of the city, there were also large areas where these methods did not perform satisfactorily. In some of these cases, researchers have proposed that layers of fine grained material within the soil profile may be responsible for preventing the manifestation of liquefaction. This paper presents preliminary findings on the mechanisms at play when pressure differentials exist across a clay layer. It is found that if the clay layer is unable to distort, then pore fluid is unable to break-through the layer even with relatively high pressures, resulting in dissipation of excess pore pressures by seepage. If the layers are however able to distort, then it is possible for the pore fluid to break through the clay layer, potentially resulting in adverse effects in terms of the severity of liquefaction.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage occurred to unreinforced masonry buildings throughout the region during the mainshock and subsequent large aftershocks. Particularly extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, a geotechnical reconnaissance was conducted over a period of six days (10–15 September 2010) by a team of geotechnical/earthquake engineers and geologists from New Zealand and USA (GEER team: Geo-engineering Extreme Event Reconnaissance). JGS (Japanese Geotechnical Society) members from Japan also participated in the reconnaissance team from 13 to 15 September 2010. The NZ, GEER and JGS members worked as one team and shared resources, information and logistics in order to conduct thorough and most efficient reconnaissance covering a large area over a very limited time period. This report summarises the key evidence and findings from the reconnaissance.
This paper provides a comparison between the strong ground motions observed in the Christchurch central business district in the 4 September 2010 Mw7.1 Darfield, and 22 February 2011 Mw6.3 Christchurch earthquakes with those observed in Tokyo during the 11 March 2011 Mw9.0 Tohoku earthquake. Despite Tokyo being located approximately 110km from the nearest part of the causative rupture, the ground motions observed from the Tohoku earthquake were strong enough to cause structural damage in Tokyo and also significant liquefaction to loose reclaimed soils in Tokyo bay. Comparisons include the strong motion time histories, response spectra, significant durations and arias intensity. The implications for large earthquakes in New Zealand are also briefly discussed.
The aim of this thesis was to examine the spatial and the temporal patterns of anxiety and chest pain resulting from the Canterbury, New Zealand earthquaeks. Three research objectives were identified: examine any spatial or termporal clusters of anxiety and chest pain; examine the associations between anxiety, chest pain and damage to neighbourhood; and determine any statistically significant difference in counts of anxiety and chest pain after each earthquake or aftershock which resulted in severe damage. Measures of the extent of liquefaction the location of CERA red-zones were used as proxy measures for earthquake damage. Cases of those who presented to Christchurch Public Hospital Emergency Department with either anxiety or chest pain between May 2010 and April 2012 were aggregated to census area unit (CAU) level for analysis. This thesis has taken a unique approach to examining the spatial and spatio-temporal variations of anxiety and chest pain after an earthquake and offers unique results. This is the first study of its kind to use a GIS approach when examining Canterbury specific earthquake damage and health variables at a CAU level after the earthquakes. Through the use of spatio-termporal scan modelling, negative and linear regression modelling and temporal linear modelling with dummy variables this research was able to conclude there are significant spatial and temporal variations in anxiety and chest pain resulting from the earthquakes. The spatio-termporal scan modelling identified a hot cluster of both anxiety and chest pain within Christchurch at the same time the earthquakes occurred. The negative binomial model found liquefaction to be a stronger predictor of anxiety than the Canterbury Earthquake Recovery Authority's (CERA) land zones. The linear regression model foun chest pain to be positively associated with all measures of earthquake damage with the exception of being in the red-zone. The temporal modelling identified a significant increase in anxiety cases one month after a major earthquake, and chest pain cases spiked two weeks after an earthquake and gradually decreased over the following five weeks. This research was limited by lack of control period data, limited measures of earthquake damage, ethical restrictions, and the need for population tracking data. The findings of this research will be useful in the planning and allocation of mental wellbeing resources should another similar event like the Canterbury Earthquakes occur in New Zealand.
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.
The 4 September 2010 Darfield and 22 February 2011 Christchurch earthquakes caused significant damage to Christchurch and surrounding suburbs as a result of the widespread liquefaction and lateral spreading that occurred. Ground surveying-based field investigations were conducted following these two events in order to measure permanent ground displacements in areas significantly affected by lateral spreading. Data was analysed with respect to the distribution of lateral spreading vs. distance from the waterway, and the failure patterns observed. Two types of failure distribution patterns were observed, a typical distributed pattern and an atypical block failure. Differences in lateral spreading measurements along adjacent banks of the Avon River in the area of Dallington were also examined. The spreading patterns between the adjacent banks varied with the respective river geometry and/or geotechnical conditions at the banks.
This paper provides an overview of the salient aspects of the dense array of ground motions observed in the 4 September 2010 Darfield and 22 February 2011 Christchurch earthquakes. Particular attention is given to inferred physical reasons for the observed ground motions, which include: (i) source features such as forward directivity effects; (ii) The effects of the Canterbury Plains sedimentary basin on basin-generated surface waves, and waveguide effects through the region; and (iii) the importance of local site response as evidenced by observations of large long period amplification and liquefaction. The significance of vertical ground motion intensity is also examined.
A series of undrained cyclic direct simple shear (DSS) tests on specimens of sandy silty soils are used to evaluate the effects of fines content, fabric and layered structure on the liquefaction response of sandy soils containing non-plastic fines. Test soils originate from shallow deposits in Christchurch, New Zealand, where severe and damaging manifestations of liquefaction occurred during the 2010-2011 Canterbury earthquakes. A procedure for reconstituting specimens by water sedimentation is employed. This specimen preparation technique involves first pluviation of soil through a water column, and then application of gentle vibrations to the mould (tapping) to prepare specimens with different initial densities. This procedure is applied to prepare uniform specimens, and layered specimens with a silt layer atop a sand layer. Cyclic DSS tests are performed on water-sedimented specimens of two sands, a silt, and sand-silt mixtures with different fines contents. Through this testing program, effects of density, time of vibration during preparation, fines content, and layered structure on cyclic behaviour and liquefaction resistance are investigated. Additional information necessary to characterise soil behaviour is provided by particle size distribution analyses, index void ratio testing, and Scanning Electronic Microscope imaging. The results of cyclic DSS tests show that, for all tested soils, specimens vibrated for longer period of time have lower void ratios, higher relative density, and greater liquefaction resistance. One of the tested sands undergoes significant increase in relative density and liquefaction resistance following prolonged vibration. The other sand exhibits lower increase in relative density and in liquefaction resistance when vibrated for the same period of time. Liquefaction resistance of sand-silt mixtures prepared using this latter sand shows a correlation with relative density irrespective of fines content. In general, however, magnitudes of changes in liquefaction resistance for given variations in vibration time, relative density, or void ratio vary depending on soils under consideration. Characterization based on maximum and minimum void ratios indicates that tested soils develop different structures as fines are added to their respective host sands. These structures influence initial specimen density, strains during consolidation, cyclic liquefaction resistance, and undrained cyclic response of each soil. The different structures are the outcome of differences in particle size distributions, average particle sizes, and particle shapes of the two host sands and of the different relationships between these properties and those of the silt. Fines content alone does not provide an effective characterization of the effects of these factors. Monotonic DSS tests are also performed on specimens prepared by water sedimentation, and on specimens prepared by moist tamping, to identify the critical state lines of tested soils. These critical state lines provide the basis for an alternative interpretation of cyclic DSS tests results within the critical state framework. It is shown that test results imply general consistency between observed cyclic and monotonic DSS soil response. The effects of specimen layering are scrutinised by comparing DSS test results for uniform and layered specimens of the same soils. In this case, only a limited number of tests is performed, and the range of densities considered for the layered specimens is also limited. Caution is therefore required in interpretation of their results. The liquefaction resistance of layered specimens appears to be influenced by the bottom sand layer, irrespective of the global fines content of the specimen. The presence of a layered structure does not result in significant differences in terms of liquefaction response with respect to uniform sand specimens. Cyclic triaxial data for Christchurch sandy silty soils available from previous studies are used to comparatively examine the behaviour observed in the tests of this study. The cyclic DSS liquefaction resistance of water-sedimented specimens is consistent with cyclic triaxial tests on undisturbed specimens performed by other investigators. The two data sets result in similar liquefaction triggering relationships for these soils. However, stress-strain response characteristics for the two types of specimens are different, and undisturbed triaxial specimen exhibit a slower rate of increase in shear strains compared to water-sedimented DSS specimens. This could be due to the greater influence of fabric of the undisturbed specimens.
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 Avon-Heathcote Estuary is of significant value to Christchurch due to its high productivity, biotic diversity, proximity to the city, and its cultural, recreational and aesthetic qualities. Nonetheless, it has been subjected to decades of degradation from sewage wastewater discharges and encroaching urban development. The result was a eutrophied estuary, high in nitrogen, affected by large blooms of nuisance macroalgae and covered by degraded sediments. In March 2010, treated wastewater was diverted from the estuary to a site 3 km offshore. This quickly reduced water nitrogen by 90% within the estuary and, within months, there was reduced production of macroalgae. However, a series of earthquakes beginning in September 2010 brought massive changes: tilting of the estuary, changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sediment and nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are an important buffer against eutrophication. Therefore, in tandem with the wastewater diversion, they could underpin much of the recovery of the estuary. Overall, the new sediments were less favourable for benthic microalgal growth and recolonisation, but were less contaminated than old sediments at highly eutrophic sites. Because the new sediments were less contaminated than the old sediments, they could help return the estuary to a noneutrophic state. However, if the new sediments, which are less favourable for microalgal growth, disperse over the old sediments at highly eutrophic sites, they could become contaminated and interfere with estuarine recovery. Therefore, recovery of microalgal communities and the estuary was expected to be generally long, but variable and site-specific, with the least eutrophic sites recovering quickly, and the most eutrophic sites taking years to return to a pre-earthquake and non-eutrophied state. changes in channels and water flow, and a huge influx of liquefied sediments that covered up to 65% of the estuary floor. Water nitrogen increased due to damage to sewage infrastructure and the diversion pipeline being turned off. Together, these drastically altered the estuarine ecosystem. My study involves three laboratory and five in situ experiments that investigate the base of the food chain and responses of benthic microalgae to earthquake-driven sedimen tand nutrient changes. It was predicted that the new sediments would be coarser and less contaminated with organic matter and nutrients than the old sediments, would have decreased microalgal biomass, and would prevent invertebrate grazing and bioturbation activities. It was believed that microalgal biomass would become similar across new and old sediments types as the unstable new sediments were resuspended and distributed over the old sediments. Contact cores of the sediment were taken at three sites, across a eutrophication gradient, monthly from September 2011 to March 2012. Extracted chlorophyll a pigments showed that microalgal biomass was generally lower on new liquefied sediments compared to old sediments, although there was considerable site to site variation, with the highly eutrophic sites being the most affected by the emergence of the new sediments. Grazer experiments showed that invertebrates had both positive and negative site-specific effects on microalgal biomass depending on their identity. At one site, new sediments facilitated grazing by Amphibola crenata, whereas at another site, new sediments did not alter the direct and indirect effects of invertebrates (Nicon aestuariensis, Macropthalmus hirtipes, and A. crenata) on microalgae. From nutrient addition experiments it was clear that benthic microalgae were able to use nutrients from within both old and new sediments equally. This implied that microalgae were reducing legacy nutrients in both sediments, and that they are
Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
A number of field testing techniques, such as standard penetration test (SPT), cone penetration test (CPT), and Swedish weight sounding (SWS), are popularly used for in-situ characterisation. The screw driving sounding (SDS) method, which has been recently developed in Japan, is an improved version of the SWS technique and measures more parameters, including the required torque, load, speed of penetration and rod friction; these provide more robust way of characterising soil stratigraphy. It is a cost-efficient technique which uses a machine-driven and portable device, making it ideal for testing in small-scale and confined areas. Moreover, with a testing depth of up to 10-15m, it is suitable for liquefaction assessment. Thus, the SDS method has great potential as an in-situ testing method for geotechnical site characterisation, especially for residential house construction. In this paper, the results of SDS tests performed at a variety of sites in New Zealand are presented. The soil database was employed to develop a soil classification chart based on SDS-derived parameters. Moreover, using the data obtained following the 2010-2011 Christchurch Earthquake Se-quence, a methodology was established for liquefaction potential evaluation using SDS data. http://www.isc5.com.au/wp-content/uploads/2016/09/1345-2-ORENSE.pdf
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.
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
This paper presents insights from recent advanced laboratory testing of undisturbed and reconstituted specimens of Christchurch silty-sands. The purpose of the testing was to establish the cyclic strength of silty-sands from sites in the Central Business District (CBD), where liquefaction was observed in 4 September 2010, 22 February 2011, and 13 June 2011. Similar overall strengths were obtained from undisturbed and reconstituted tests prepared at similar densities, albeit with higher variability for the reconstituted specimens. Reconstituted specimens exhibited distinctly different response in terms of lower compressibility during initial loading cycles, and exhibited a more brittle response when large strains were mobilised, particularly for samples with high fines content. Given the lower variability in natural sample response and the possibility of age-related strength to be significant for sites not subjected to earthquakes, high quality undisturbed samples are recommended over the use of reconstituted specimens to establish the cyclic strength of natural sands.
The seismic performance of soil profiles with potentially liquefiable deposits is a complex phenomenon that requires a thorough understanding of the soil properties and ground motion characteristics. The limitations of simplified liquefaction assessment methods have prompted an increase in the use of non-linear dynamic analysis methods. Focusing on onedimensional site response of a soil column, this thesis validated a soil constitutive model using in-situ pore pressure measurements and then assessed the influence of input ground motion characteristics on soil column response using traditional and newly developed metrics. Pore pressure recordings during the Canterbury Earthquake Sequence (CES) in New Zealand were used to validate the PM4Sand constitutive model. Soil profile characterization was key to accurate prediction of excess pore pressure response and accounting for any densification during the CES. Response during multiple earthquakes was captured effectively and cross-layer interaction demonstrated the model capability to capture soil response at the system-level. Synthetic and observed ground motions from the Christchurch earthquake were applied to the validated soil column to quantify the performance of synthetic motions. New metrics were developed to facilitate a robust comparison to assess performance. The synthetic input motions demonstrated a slightly larger acceleration and excess pore pressure response compared to the observed input motions. The results suggest that the synthetic motions may accumulate higher excess pore pressure at a faster rate and with fewer number of cycles in the shear response. This research compares validated soil profile subject to spectrally-matched pulse and non-pulse motions, emphasizing the inclusion of pulse motions with distinctive characteristics in ground motion suites for non-linear dynamic analysis. However, spectral matching may lead to undesired alterations in pulse characteristics. Cumulative absolute velocity and significant duration significantly differed between these two groups compared to the other key characteristics and contributed considerably to the liquefaction response. Unlike the non-pulse motions, not all of the pulse motions triggered liquefaction, likely due to their shorter significant duration. Non-pulse motions developed a greater spatial extent of liquefaction triggering in the soil profile and extended to a greater depth.
Active faults capable of generating highly damaging earthquakes may not cause surface rupture (i.e., blind faults) or cause surface ruptures that evade detection due to subsequent burial or erosion by surface processes. Fault populations and earthquake frequency-‐magnitude distributions adhere to power laws, implying that faults too small to cause surface rupture but large enough to cause localized strong ground shaking densely populate continental crust. The rupture of blind, previously undetected faults beneath Christchurch, New Zealand in a suite of earthquakes in 2010 and 2011, including the fatal 22 February 2011 moment magnitude (Mw) 6.2 Christchurch earthquake and other large aftershocks, caused a variety of environmental impacts, including major rockfall, severe liquefaction, and differential surface uplift and subsidence. All of these effects occurred where geologic evidence for penultimate effects of the same nature existed. To what extent could the geologic record have been used to infer the presence of proximal, blind and / or unidentified faults near Christchurch? In this instance, we argue that phenomena induced by high intensity shaking, such as rock fragmentation and rockfall, revealed the presence of proximal active faults in the Christchurch area prior to the recent earthquake sequence. Development of robust earthquake shaking proxy datasets should become a higher scientific priority, particularly in populated regions.
Severe liquefaction was repeatedly observed during the 2010 - 2011 C hristchurch earthquake s , particularly affecting deposits of fine sands and silty sands of recent fluvial or estuarine origin. The effects of liquefaction included major sliding of soil tow ard water bodies ( i.e. lateral spreading ) rang ing from centimetres to several metres. In this paper, a series of undrained cyclic torsional shear tests were conducted to evaluate the liquefaction and extremely large deformation properties of Christchurch b oiled sand . In these tests, the simple shear conditions were reproduced in order to apply realistic stress conditions that soil s experience in the field during horizontal seismic shaking. Several hollow cylindrical medium dense specimens ( D r = 50%) were pr epared by pluviation method, isotropically consolidated at an effective stress of 100 kPa and then cyclically sheared under undrained conditions up to 10 0% double amplitude shear strain (γ DA ) . The cyclic strength at different levels of γ DA of 7.5%, 15%, 3 0 % and 6 0%, development of extremely large post - liquefaction deformation and shear strain locali s ation properties were assessed from the analysis of the effective stress paths and stress - strain responses . To reveal possible distinctiveness, the cyclic undra ined behaviour of CHCH boiled sand was compared with that of Toyoura sand previously examined under similar testing conditions
Farming and urban regions are impacted by earthquake disasters in different ways, and feature a range of often different recovery requirements. In New Zealand, and elsewhere, most earthquake impact and recovery research is urban focused. This creates a research deficit that can lead to the application of well-researched urban recovery strategies in rural areas to suboptimal effect. To begin to reduce this deficit, in-depth case studies of the earthquake impacts and recovery of three New Zealand farms severely impacted by the 14th November 2016, M7.8 Hurunui-Kaikōura earthquake were conducted. The initial earthquake, its aftershocks and coseismic hazards (e.g., landslides, liquefaction, surface rupture) affected much of North Canterbury, Marlborough and the Wellington area. The three case study farms were chosen to broadly represent the main types of farming and topography in the Hurunui District in North Canterbury. The farms were directly and indirectly impacted by earthquakes and related hazards. On-farm infrastructure (e.g., woolsheds, homesteads) and essential services (e.g., water, power), frequently sourced from distributed networks, were severely impacted. The earthquake occurred after two years of regional drought had already stressed farm systems and farmers to restructuring or breaking point. Cascading interlinked hazards stemming from the earthquakes and coseismic hazards continued to disrupt earthquake recovery over a year after the initial earthquake. Semi-structured interviews with the farmers were conducted nine and fourteen months after the initial earthquake to capture the timeline of on-going impacts and recovery. Analysis of both geological hazard data and interview data resulted in the identification of key factors influencing farm level earthquake impact and recovery. These include pre-existing conditions (e.g., drought); farm-specific variations in recovery timelines; and resilience strategies for farm recovery resources. The earthquake recovery process presented all three farms with opportunities to change their business plans and adapt to mitigate on-going and future risk.
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.
This thesis presents the application of data science techniques, especially machine learning, for the development of seismic damage and loss prediction models for residential buildings. Current post-earthquake building damage evaluation forms are developed for a particular country in mind. The lack of consistency hinders the comparison of building damage between different regions. A new paper form has been developed to address the need for a global universal methodology for post-earthquake building damage assessment. The form was successfully trialled in the street ‘La Morena’ in Mexico City following the 2017 Puebla earthquake. Aside from developing a framework for better input data for performance based earthquake engineering, this project also extended current techniques to derive insights from post-earthquake observations. Machine learning (ML) was applied to seismic damage data of residential buildings in Mexico City following the 2017 Puebla earthquake and in Christchurch following the 2010-2011 Canterbury earthquake sequence (CES). The experience showcased that it is readily possible to develop empirical data only driven models that can successfully identify key damage drivers and hidden underlying correlations without prior engineering knowledge. With adequate maintenance, such models have the potential to be rapidly and easily updated to allow improved damage and loss prediction accuracy and greater ability for models to be generalised. For ML models developed for the key events of the CES, the model trained using data from the 22 February 2011 event generalised the best for loss prediction. This is thought to be because of the large number of instances available for this event and the relatively limited class imbalance between the categories of the target attribute. For the CES, ML highlighted the importance of peak ground acceleration (PGA), building age, building size, liquefaction occurrence, and soil conditions as main factors which affected the losses in residential buildings in Christchurch. ML also highlighted the influence of liquefaction on the buildings losses related to the 22 February 2011 event. Further to the ML model development, the application of post-hoc methodologies was shown to be an effective way to derive insights for ML algorithms that are not intrinsically interpretable. Overall, these provide a basis for the development of ‘greybox’ ML models.
The region in and around Christchurch, encompassing Christchurch city and the Selwyn and Waimakariri districts, contains more than 800 road, rail, and pedestrian bridges. Most of these bridges are reinforced concrete, symmetric, and have small to moderate spans (15–25 m). The 22 February 2011 moment magnitude (Mw) 6.2 Christchurch earthquake induced high levels of localized ground shaking (Bradley and Cubrinovski 2011, page 853 of this issue; Guidotti et al. 2011, page 767 of this issue; Smyrou et al. 2011, page 882 of this issue), with damage to bridges mainly confined to the central and eastern parts of Christchurch. Liquefaction was evident over much of this part of the city, with lateral spreading affecting bridges spanning both the Avon and Heathcote rivers.
During the recent devastating earthquakes in Christchurch, many residential houses were damaged due to widespread liquefaction of the ground. In-situ testing is widely used as a convenient method for evaluating liquefaction potential of soils. Cone penetration test (CPT) and standard penetration test (SPT) are the two popular in situ tests which are widely used in New Zealand for site characterization. The Screw Driving Sounding (SDS) method is a relatively new operating system developed in Japan consisting of a machine that drills a rod into the ground by applying torque at seven steps of axial loading. This machine can continuously measure the required torque, load, speed of penetration and rod friction during the test, and therefore can give a clear overview of the soil profile along the depth of penetration. In this paper, based on a number of SDS tests conducted in Christchurch, a correlation was developed between tip resistance of CPT test and SDS parameters for layers consisting of different fines contents. Moreover, using the obtained correlation, a chart was proposed which relates the cyclic resistance ratio to the appropriate SDS parameter. Using the proposed chart, liquefaction potential of soil can be estimated directly using SDS data. As SDS method is simpler, faster and more economical test than CPT and SPT, it can be a reliable alternative in-situ test for soil characterization, especially in residential house constructions.
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.
