The city of Christchurch and its surrounds experienced widespread damage due to soil liquefaction induced by seismic shaking during the Canterbury earthquake sequence that began in September 2010 with the Mw7.1 Darfield earthquake. Prior to the start of this sequence, the city had a large network of strong motion stations (SMSs) installed, which were able to record a vast database of strong ground motions. This paper uses this database of strong ground motion recordings, observations of liquefaction manifestation at the ground surface, and data from a recently completed extensive geotechnical site investigation program at each SMS to assess a range of liquefaction evaluation procedures at the four SMSs in the Christchurch Central Business District (CBD). In general, the characteristics of the accelerograms recorded at each SMS correlated well with the liquefaction evaluation procedures, with low liquefaction factors of safety predicted at sites with clear liquefaction identifiers in the ground motions. However, at sites that likely liquefied at depth (as indicated by evaluation procedures and/or inferred from the characteristics of the recorded surface accelerograms), the presence of a non-liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects. Because of this, there was not a good correlation between surface manifestation and two surface manifestation indices, the Liquefaction Potential Index (LPI) and the Liquefaction Severity Number (LSN).
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.
Topics - appalling stories of Christchurch families living in garages and tents - even a couple renting a wash-house - are continuing to surface almost two years on from the February 2011 earthquake. Youth crime is described as "steadily falling" in New Zealand, with rates of police apprehensions and prosecutions of young people at record lows. Auckland Grammar has employed the services of a private investigator to check on parents trying to cheat on school zone boundaries.
This report 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 during the Canterbury earthquake sequence. 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 shallow geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site is 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 to better understand stratigraphy and soil properties, which are important in seismic site response. In this regard, CPT, SPT and borehole data, shear wave velocity (Vs) profiles, and horizontal to vertical spectral ratio measurements (H/V) in close vicinity to the SMS were used to develop representative soil profiles at each site. NZS1170.5 (SNZ 2004) site subsoil classifications were updated using Vs and SPT N60 criteria. Site class E boundaries were treated as a sliding scale rather than as a discrete boundary to account for locations with similar site effects potential, an approach which was shown to result in a better delineation between the site classes. SPT N60 values often indicate a stiffer site class than the Vs data for softer soil sites, highlighting the disparity between the two site investigation techniques. Both SPT N60 and Vs based site classes did not always agree with the original site classifications. This emphasises 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 N60 and Vs. Liquefaction triggering assessments were carried out for the Darfield and Christchurch earthquakes, and compared against observed liquefaction surface manifestations and ground motions characteristics at each SMS. In general, the characteristics of the recorded ground motions at each site correlate well with the triggering analyses. However, at sites that likely liquefied at depth (as indicated by triggering analyses and/or inferred from the characteristics of the recorded surface acceleration time series), the presence of a non‐liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
Measurement of basement seismic resonance frequencies can elucidate shallow velocity structure, an important factor in earthquake hazard estimation. Ambient noise cross correlation, which is well-suited to studying shallow earth structure, is commonly used to analyze fundamental-mode Rayleigh waves and, increasingly, Love waves. Here we show via multicomponent ambient noise cross correlation that the basement resonance frequency in the Canterbury region of New Zealand can be straightforwardly determined based on the horizontal to vertical amplitude ratio (H/V ratio) of the first higher-mode Rayleigh waves. At periods of 1-3 s, the first higher-mode is evident on the radial-radial cross-correlation functions but almost absent in the vertical-vertical cross-correlation functions, implying longitudinal motion and a high H/V ratio. A one-dimensional regional velocity model incorporating a ~ 1.5 km-thick sedimentary layer fits both the observed H/V ratio and Rayleigh wave group velocity. Similar analysis may enable resonance characteristics of other sedimentary basins to be determined. © 2013. American Geophysical Union. All Rights Reserved.
Geosynthetic reinforced soil (GRS) walls involve the use of geosynthetic reinforcement (polymer material) within the retained backfill, forming a reinforced soil block where transmission of overturning and sliding forces on the wall to the backfill occurs. Key advantages of GRS systems include the reduced need for large foundations, cost reduction (up to 50%), lower environmental costs, faster construction and significantly improved seismic performance as observed in previous earthquakes. Design methods in New Zealand have not been well established and as a result, GRS structures do not have a uniform level of seismic and static resistance; hence involve different risks of failure. Further research is required to better understand the seismic behaviour of GRS structures to advance design practices. The experimental study of this research involved a series of twelve 1-g shake table tests on reduced-scale (1:5) GRS wall models using the University of Canterbury shake-table. The seismic excitation of the models was unidirectional sinusoidal input motion with a predominant frequency of 5Hz and 10s duration. Seismic excitation of the model commenced at an acceleration amplitude level of 0.1g and was incrementally increased by 0.1g in subsequent excitation levels up to failure (excessive displacement of the wall panel). The wall models were 900mm high with a full-height rigid facing panel and five layers of Microgird reinforcement (reinforcement spacing of 150mm). The wall panel toe was founded on a rigid foundation and was free to slide. The backfill deposit was constructed from dry Albany sand to a backfill relative density, Dr = 85% or 50% through model vibration. The influence of GRS wall parameters such as reinforcement length and layout, backfill density and application of a 3kPa surcharge on the backfill surface was investigated in the testing sequence. Through extensive instrumentation of the wall models, the wall facing displacements, backfill accelerations, earth pressures and reinforcement loads were recorded at the varying levels of model excitation. Additionally, backfill deformation was also measured through high-speed imaging and Geotechnical Particle Image Velocimetry (GeoPIV) analysis. The GeoPIV analysis enabled the identification of the evolution of shear strains and volumetric strains within the backfill at low strain levels before failure of the wall thus allowing interpretations to be made regarding the strain development and shear band progression within the retained backfill. Rotation about the wall toe was the predominant failure mechanism in all excitation level with sliding only significant in the last two excitation levels, resulting in a bi-linear displacement acceleration curve. An increase in acceleration amplification with increasing excitation was observed with amplification factors of up to 1.5 recorded. Maximum seismic and static horizontal earth pressures were recorded at failure and were recorded at the wall toe. The highest reinforcement load was recorded at the lowest (deepest in the backfill) reinforcement layer with a decrease in peak load observed at failure, possibly due to pullout failure of the reinforcement layer. Conversely, peak reinforcement load was recorded at failure for the top reinforcement layer. The staggered reinforcement models exhibited greater wall stability than the uniform reinforcement models of L/H=0.75. However, similar critical accelerations were determined for the two wall models due to the coarseness of excitation level increments of 0.1g. The extended top reinforcements were found to restrict the rotational component of displacement and prevented the development of a preliminary shear band at the middle reinforcement layer, contributing positively to wall stability. Lower acceleration amplification factors were determined for the longer uniform reinforcement length models due to reduced model deformation. A greater distribution of reinforcement load towards the top two extended reinforcement layers was also observed in the staggered wall models. An increase in model backfill density was observed to result in greater wall stability than an increase in uniform reinforcement length. Greater acceleration amplification was observed in looser backfill models due to their lower model stiffness. Due to greater confinement of the reinforcement layers, greater reinforcement loads were developed in higher density wall models with less wall movement required to engage the reinforcement layers and mobilise their resistance. The application of surcharge on the backfill was observed to initially increase the wall stability due to greater normal stresses within the backfill but at greater excitation levels, the surcharge contribution to wall destabilising inertial forces outweighs its contribution to wall stability. As a result, no clear influence of surcharge on the critical acceleration of the wall models was observed. Lower acceleration amplification factors were observed for the surcharged models as the surcharge acts as a damper during excitation. The application of the surcharge also increases the magnitude of reinforcement load developed due to greater confinement and increased wall destabilising forces. The rotation of the wall panel resulted in the progressive development of shears surface with depth that extended from the backfill surface to the ends of the reinforcement (edge of the reinforced soil block). The resultant failure plane would have extended from the backfill surface to the lowest reinforcement layer before developing at the toe of the wall, forming a two-wedge failure mechanism. This is confirmed by development of failure planes at the lowest reinforcement layer (deepest with the backfill) and at the wall toe observed at the critical acceleration level. Key observations of the effect of different wall parameters from the GeoPIV results are found to be in good agreement with conclusions developed from the other forms of instrumentation. Further research is required to achieve the goal of developing seismic guidelines for GRS walls in geotechnical structures in New Zealand. This includes developing and testing wall models with a different facing type (segmental or wrap-around facing), load cell instrumentation of all reinforcement layers, dynamic loading on the wall panel and the use of local soils as the backfill material. Lastly, the limitations of the experimental procedure and wall models should be understood.
The connections between walls of unreinforced masonry (URM) buildings and flexible timber diaphragms are critical building components that must perform adequately before desirable earthquake response of URM buildings may be achieved. Field observations made during the initial reconnaissance and the subsequent damage surveys of clay brick URM buildings following the 2010/2011 Canterbury, New Zealand earthquakes revealed numerous cases where anchor connections joining masonry walls or parapets with roof or floor diaphragms appeared to have failed prematurely. These observations were more frequent for adhesive anchor connections than for through-bolt connections (i.e. anchorages having plates on the exterior façade of the masonry walls). Subsequently, an in-field test program was undertaken in an attempt to evaluate the performance of adhesive anchor connections between unreinforced clay brick URM walls and roof or floor diaphragms. The study consisted of a total of almost 400 anchor tests conducted in eleven existing URM buildings located in Christchurch, Whanganui and Auckland. Specific objectives of the study included the identification of failure modes of adhesive anchors in existing URM walls and the influence of the following variables on anchor load-displacement response: adhesive type, strength of the masonry materials (brick and mortar), anchor embedment depth, anchor rod diameter, overburden level, anchor rod type, quality of installation and the use of metal mesh sleeve. In addition, the comparative performance of bent anchors (installed at an angle of minimum 22.5o to the perpendicular projection from the wall surface) and anchors positioned horizontally was investigated. Observations on the performance of wall-to-diaphragm connections in the 2010/2011 Canterbury earthquakes, a snapshot of the performed experimental program and the test results and a preliminary proposed pull-out capacity of adhesive anchors are presented herein.
The connections between walls of unreinforced masonry (URM) buildings and flexible timber diaphragms are critical building components that must perform adequately before desirable earthquake response of URM buildings may be achieved. Field observations made during the initial reconnaissance and the subsequent damage surveys of clay brick URM buildings following the 2010/2011 Canterbury, New Zealand earthquakes revealed numerous cases where anchor connections joining masonry walls or parapets with roof or floor diaphragms appeared to have failed prematurely. These observations were more frequent for adhesive anchor connections than for through-bolt connections (i.e. anchorages having plates on the exterior façade of the masonry walls). Subsequently, an in-field test program was undertaken in an attempt to evaluate the performance of adhesive anchor connections between unreinforced clay brick URM walls and roof or floor diaphragm. The study consisted of a total of almost 400 anchor tests conducted in eleven existing URM buildings located in Christchurch, Whanganui and Auckland. Specific objectives of the study included the identification of failure modes of adhesive anchors in existing URM walls and the influence of the following variables on anchor load-displacement response: adhesive type, strength of the masonry materials (brick and mortar), anchor embedment depth, anchor rod diameter, overburden level, anchor rod type, quality of installation and the use of metal mesh sleeve. In addition, the comparative performance of bent anchors (installed at an angle of minimum 22.5o to the perpendicular projection from the wall surface) and anchors positioned horizontally was investigated. Observations on the performance of wall-to-diaphragm connections in the 2010/2011 Canterbury earthquakes, a snapshot of the performed experimental program and the test results and a preliminary proposed pull-out capacity of adhesive anchors are presented herein. http://www.confer.co.nz/nzsee/ VoR - Version of Record
It is well known that buildings constructed using unreinforced masonry (URM) are susceptible to damage from earthquake induced lateral forces that may result in partial or full building collapse. The 2010/2011 Canterbury earthquakes are the most recent New Zealand example of destructive earthquakes, which have drawn people's attention to the inherent seismic weaknesses of URM buildings and anchored masonry veneer systems in New Zealand. A brief review of the data collected following the 2010 Darfield earthquake and more comprehensive documentation of data that was collected following the 2011 Christchurch earthquake is presented, along with the findings from subsequent data interrogation. Large stocks of earthquake prone vintage URM buildings that remain in New Zealand and in other seismically active parts of the world result in the need for minimally invasive and cost effective seismic retrofit techniques. The principal objective of the doctoral research reported herein was to investigate the applicability of near surface mounted (NSM) carbon fibre reinforced polymer (CFRP) strips as a seismic improvement technique. A comprehensive experimental program consisting of 53 pull tests is presented and is used to assess the accuracy of existing FRP-to-masonry bond models, with a modified model being proposed. The strength characteristics of vintage clay brick URM wall panels from two existing URM buildings was established and used as a benchmark when manufacturing replica clay brick test assemblages. The applicability of using NSM CFRP strips as a retrofitting technique for improving the shear strength and the ductility capacity of multi-leaf URM walls constructed using solid clay brick masonry is investigated by varying CFRP reinforcement ratios. Lastly, an experimental program was undertaken to validate the proposed design methodology for improving the strength capacity of URM walls. The program involved testing full-scale walls in a laboratory setting and testing full-scale walls in-situ in existing vintage URM buildings. Experimental test results illustrated that the NSM CFRP technique is an effective method to seismically strengthen URM buildings.
The collapse of Redcliffs’ cliff in the 22 February 2011 and 13 June 2011 earthquakes were the first times ever a major failure incident occurred at Redcliffs in approximately 6000 years. This master’s thesis is a multidisciplinary engineering geological investigation sought to study these particular failure incidents, focusing on collecting the data necessary to explain the cause and effect of the cliff collapsing in the event of two major earthquakes. This study provides quantitative and qualitative data about the geotechnical attributes and engineering geological nature of the sea-cut cliff located at Redcliffs. Results from surveying the geology of Redcliffs show that the exposed lithology of the cliff face is a variably jointed rock body of welded and (relatively intact) unwelded ignimbrite, a predominantly massive unit of brecciated tuff, and a covering of wind-blown loess and soil deposit (commonly found throughout Canterbury) on top of the cliff. Moreover, detailing the external component of the slope profile shows that Redcliffs’ cliff is a 40 – 80 m cliff with two intersecting (NE and SE facing) slope aspects. The (remotely) measured geometry of the cliff face comprises of multiple outstanding gradients, averaging a slope angle of ~67 degrees (post-13 June 2011), where the steepest components are ~80 degrees, whereas the gentle sloping sections are ~44 degrees. The physical structure of Redcliffs’ cliff drastically changed after each collapse, whereby seismically induced alterations to the slope geometry resulted in material deposited on the talus at the base of the cliff. Prior to the first collapse, the variance of the gradient down the slope was minimal, with the SE Face being the most variable with up to three major gradients on one cross section. However, after each major collapse, the variability increased with more parts of the cliff face having more than one major gradient that is steeper or gentler than the remainder of the slope. The estimated volume of material lost as a result of the gradient changes was 28,267 m³ in February and 11,360 m³ in June 2011. In addition, surveys of the cliff top after the failure incidents revealed the development of fissures along the cliff edge. Monitoring 10 fissures over three months indicated that fissured by the cliff edge respond to intense seismicity (generally ≥ Mw 4) by widening. Redcliffs’ cliff collapsed on two separate occasions as a result of an accumulated amount of damage of the rock masses in the cliff (caused by weathering and erosion over time), and two Mw 6.2 trigger earthquakes which shook the Redcliffs and the surrounding area at a Peak Ground Acceleration (PGA) estimated to be around 2 g. The results of the theoretical study suggests that PGA levels felt on-site during both instances of failure are the result of three major factors: source of the quake and the site affected; topographic amplification of the ground movement; the short distance between the source and the cliff for both fault ruptures; the focus of seismic energy in the direction of thrust faulting along a path that intercepts Redcliffs (and the Port Hills). Ultimately, failure on the NE and SE Faces of Redcliffs’ cliff was concluded to be global as every part of the exposed cliff face deposited a significant volume of material on the talus at the base of the cliff, with the exception of one section on the NE Face. The cliff collapses was a concurrent process that is a single (non-monotonic) event that operated as a complex series of (primarily) toppling rock falls, some sliding of blocks, and slumping of the soil mantle on top of the cliff. The first collapse had a mixture of equivalent continua slope movement of the heavily weathered / damaged surface of the cliff face, and discontinuous slope movement of the jointed inner slope (behind the heavily weathered surface); whereas the second collapse resulted in only discontinuous slope movement on account of the freshly exposed cliff face that had damage to the rock masses, in the form of old and (relatively) new discontinuous fractures, induced by earthquakes and aftershocks leading up to the point of failure.
The connections between walls of unreinforced masonry (URM) buildings and flexible timber diaphragms are critical building components that must perform adequately before desirable earthquake response of URM buildings may be achieved. Field observations made during the initial reconnaissance and the subsequent damage surveys of clay brick URM buildings following the 2010/2011 Canterbury, New Zealand earthquakes revealed numerous cases where anchor connections joining masonry walls or parapets with roof or floor diaphragms appeared to have failed prematurely. These observations were more frequent for the case of adhesive anchor connections than for the case of through-bolt connections (i.e. anchorages having plates on the exterior façade of the masonry walls). Subsequently, an in-field test program was undertaken in an attempt to evaluate the performance of adhesive anchor connections between unreinforced clay brick URM walls and roof or floor diaphragm. The study consisted of a total of almost 400 anchor tests conducted in eleven existing URM buildings located in Christchurch, Whanganui and Auckland. Specific objectives of the study included the identification of failure modes of adhesive anchors in existing URM walls and the influence of the following variables on anchor load-displacement response: adhesive type, strength of the masonry materials (brick and mortar), anchor embedment depth, anchor rod diameter, overburden level, anchor rod type, quality of installation and the use of metal foil sleeve. In addition, the comparative performance of bent anchors (installed at an angle of minimum 22.5o to the perpendicular projection from the wall surface) and anchors positioned horizontally was investigated. Observations on the performance of wall-to-diaphragm connections in the 2010/2011 Canterbury earthquakes and a snapshot of the performed experimental program and the test results are presented herein. http://hdl.handle.net/2292/21050
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 Canterbury Earthquakes of 2010-2011, in particular the 4th September 2010 Darfield earthquake and the 22nd February 2011 Christchurch earthquake, produced severe and widespread liquefaction in Christchurch and surrounding areas. The scale of the liquefaction was unprecedented, and caused extensive damage to a variety of man-made structures, including residential houses. Around 20,000 residential houses suffered serious damage as a direct result of the effects of liquefaction, and this resulted in approximately 7000 houses in the worst-hit areas being abandoned. Despite the good performance of light timber-framed houses under the inertial loads of the earthquake, these structures could not withstand the large loads and deformations associated with liquefaction, resulting in significant damage. The key structural component of houses subjected to liquefaction effects was found to be their foundations, as these are in direct contact with the ground. The performance of house foundations directly influenced the performance of the structure as a whole. Because of this, and due to the lack of research in this area, it was decided to investigate the performance of houses and in particular their foundations when subjected to the effects of liquefaction. The data from the inspections of approximately 500 houses conducted by a University of Canterbury summer research team following the 4th September 2010 earthquake in the worst-hit areas of Christchurch were analysed to determine the general performance of residential houses when subjected to high liquefaction loads. This was followed by the detailed inspection of around 170 houses with four different foundation types common to Christchurch and New Zealand: Concrete perimeter with short piers constructed to NZS3604, concrete slab-on-grade also to NZS3604, RibRaft slabs designed by Firth Industries and driven pile foundations. With a focus on foundations, floor levels and slopes were measured, and the damage to all areas of the house and property were recorded. Seven invasive inspections were also conducted on houses being demolished, to examine in more detail the deformation modes and the causes of damage in severely affected houses. The simplified modelling of concrete perimeter sections subjected to a variety of liquefaction-related scenarios was also performed, to examine the comparative performance of foundations built in different periods, and the loads generated under various bearing loss and lateral spreading cases. It was found that the level of foundation damage is directly related to the level of liquefaction experienced, and that foundation damage and liquefaction severity in turn influence the performance of the superstructure. Concrete perimeter foundations were found to have performed most poorly, suffering high local floor slopes and being likely to require foundation repairs even when liquefaction was low enough that no surface ejecta was seen. This was due to their weak, flexible foundation structure, which cannot withstand liquefaction loads without deforming. The vulnerability of concrete perimeter foundations was confirmed through modelling. Slab-on-grade foundations performed better, and were unlikely to require repairs at low levels of liquefaction. Ribraft and piled foundations performed the best, with repairs unlikely up to moderate levels of liquefaction. However, all foundation types were susceptible to significant damage at higher levels of liquefaction, with maximum differential settlements of 474mm, 202mm, 182mm and 250mm found for concrete perimeter, slab-on-grade, ribraft and piled foundations respectively when subjected to significant lateral spreading, the most severe loading scenario caused by liquefaction. It was found through the analysis of the data that the type of exterior wall cladding, either heavy or light, and the number of storeys, did not affect the performance of foundations. This was also shown through modelling for concrete perimeter foundations, and is due to the increased foundation strengths provided for heavily cladded and two-storey houses. Heavy roof claddings were found to increase the demands on foundations, worsening their performance. Pre-1930 concrete perimeter foundations were also found to be very vulnerable to damage under liquefaction loads, due to their weak and brittle construction.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.
