Photograph by Jennifer Middendorf 426
Images, UC QuakeStudies
A digger loading demolition rubble into a truck.
A digger loading demolition rubble into a truck.
A digger loading demolition rubble into a truck on Oxford Terrace.
The current seismic design practice for reinforced concrete (RC) walls has been drawn into question following the Canterbury earthquakes. An overview of current research being undertaken at the University of Auckland into the seismic behaviour of RC walls is presented. The main objectives of this research project are to understand the observed performance of several walls in Christchurch, quantify the seismic loads on RC walls, and developed improved design procedures for RC walls that will assist in revisions to NZS 3101. A database summarising of the performance of RC wall buildings in the Christchurch CBD was collated to identify damage modes and case-study buildings. A detailed investigation is underway to verify the seismic performance of lightly reinforced concrete walls and an experimental setup has been developed to subject RC wall specimen to loading that is representative of a multi-storey building. Numerical modelling is being used to understand the observed performance of several case-study RC walls buildings in Christchurch. Of particular interest is the influence that interactions between walls and other structural elements have on the seismic response of buildings and the loads generated on RC walls.
The capability of self-compacting concrete (SCC) in flowing through and filling in even the most congested areas makes it ideal for being used in congested reinforced concrete (RC) structural members such as beam-column joints (BCJ). However, members of tall multi-storey structures impose high capacity requirements where implementing normal-strength self-compacting concrete is not preferable. In the present study, a commercially reproducible high-strength self-compacting concrete (HSSCC), a conventionally vibrated high-strength concrete (CVHSC) and a normal strength conventionally vibrated concrete (CVC) were designed using locally available materials in Christchurch, New Zealand. Following the guidelines of the New Zealand concrete standards NZS3101, seven beam-column joints (BCJ) were designed. Factors such as the concrete type, grade of reinforcement, amount of joint shear stirrups, axial load, and direction of casting were considered variables. All BCJs were tested under a displacement-controlled quasi-static reversed cyclic regime. The cracking pattern at different load levels and the mode of failure were also recorded. In addition, the load, displacement, drift, ductility, joint shear deformations, and elongation of the plastic hinge zone were also measured during the experiment. It was found that not only none of the seismically important features were compromised by using HSSCC, but also the quality of material and ease of construction boosted the performance of the BCJs.
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.
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.
During the 2010/2011 Canterbury earthquakes, several reinforced concrete (RC) walls in multi-storey buildings formed a single crack in the plastic hinge region as opposed to distributed cracking. In several cases the crack width that was required to accommodate the inelastic displacement of the building resulted in fracture of the vertical reinforcing steel. This type of failure is characteristic of RC members with low reinforcement contents, where the area of reinforcing steel is insufficient to develop the tension force required to form secondary cracks in the surrounding concrete. The minimum vertical reinforcement in RC walls was increased in NZS 3101:2006 with the equation for the minimum vertical reinforcement in beams also adopted for walls, despite differences in reinforcement arrangement and loading. A series of moment-curvature analyses were conducted for an example RC wall based on the Gallery Apartments building in Christchurch. The analysis results indicated that even when the NZS 3101:2006 minimum vertical reinforcement limit was satisfied for a known concrete strength, the wall was still susceptible to sudden failure unless a significant axial load was applied. Additionally, current equations for minimum reinforcement based on a sectional analysis approach do not adequately address the issues related to crack control and distribution of inelastic deformations in ductile walls.
This paper presents a critical evaluation of vertical ground motions observed in the Canterbury earthquake sequence. The abundance of strong near-source ground-motion recordings provides an opportunity to comprehensively review the estimation of vertical ground motions via the New Zealand Standard for earthquake loading, NZS1170.5:2004, and empirical ground motion prediction equations (GMPEs). An in-depth review of current GMPEs is carried out to determine the existing trends and characteristics present in the empirical models. Results illustrate that vertical ground motion amplitudes estimated based on NZS1170.5:2004 are significantly unconservative at short periods and near-source distances. While conventional GMPEs provide an improved prediction, in many instances they too underpredict vertical ground motion accelerations at short periods and near-source distances.
Following the 2010/2011 Canterbury, New Zealand earthquakes, a detailed door-to-door survey was conducted in the Christchurch region to establish the earthquake performance of lightweight timber-framed residential dwellings with a masonry veneer external cladding system. The post-earthquake survey involved documenting the condition of dwellings in areas that had experienced different levels of earthquake shaking, allowing comparison between the performance of different veneer systems and different shaking intensities. In total, just fewer than 1,100 residential dwellings were inspected throughout the wider Christchurch area. The survey included parameters such as level of veneer damage, type of veneer damage, observed crack widths, and level of repair required. It is concluded that based on observed earthquake performance at the shaking intensities matching or exceeding ultimate limit state loading, the post-1996 veneer fixing details performed satisfactorily and continued use of the detail is recommended without further modification. AM - Accepted Manuscript
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.
Seismic retrofitting of unreinforced masonry buildings using posttensioning has been the topic of many recent experimental research projects. However, the performance of such retrofit designs in actual design level earthquakes has previously been poorly documented. In 1984 two stone masonry buildings within The Arts Centre of Christchurch received posttensioned seismic retrofits, which were subsequently subjected to design level seismic loads during the 2010/2011 Canterbury earthquake sequence. These 26 year old retrofits were part of a global scheme to strengthen and secure the historic building complex and were subject to considerable budgetary constraints. Given the limited resources available at the time of construction and the current degraded state of the steel posttension tendons, the posttensioned retrofits performed well in preventing major damage to the overall structure of the two buildings in the Canterbury earthquakes. When compared to other similar unretrofitted structures within The Arts Centre, it is demonstrated that the posttensioning significantly improved the in-plane and out-of-plane wall strength and the ability to limit residual wall displacements. The history of The Arts Centre buildings and the details of the Canterbury earthquakes is discussed, followed by examination of the performance of the posttension retrofits and the suitability of this technique for future retrofitting of other historic unreinforced masonry buildings. http://www.aees.org.au/downloads/conference-papers/2013-2/
In recent years, significant research has been undertaken into the development of lead-extrusion damping technology. The high force-to-volume (HF2V) devices developed at the University of Canterbury have been the subject of much of this research. However, while these devices have undergone a limited range of velocity testing, limitations in test equipment has meant that they have never been tested at representative earthquake velocities. Such testing is important as the peak resistive force provided by the dampers under large velocity spikes is an important design input that must be known for structural applications. This manuscript presents the high-speed testing of HF2V devices with quasi-static force capacities of 250-300kN. These devices have been subjected to peak input velocities of approximately 200mm/s, producing peak resistive forces of approximately 350kN. The devices show stable hysteretic performance, with slight force reduction during high-speed testing due to heat build-up and softening of the lead working material. This force reduction is recovered following cyclic loading as heat is dissipated and the lead hardens again. The devices are shown to be only weakly velocity dependent, an advantage in that they do not deliver large forces to the connecting elements and surrounding structure if larger than expected response velocities occur. This high-speed testing is an important step towards uptake as it provides important information to designers.
In recent years, significant research has been undertaken into the development of lead-extrusion damping technology. The high force-to-volume (HF2V) devices developed at the University of Canterbury have been the subject of much of this research. However, while these devices have undergone a limited range of velocity testing, limitations in test equipment has meant that they have never been tested at representative earthquake velocities. Such testing is important as the peak resistive force provided by the dampers under large velocity spikes is an important design input that must be known for structural applications. This manuscript presents the high-speed testing of HF2V devices with quasi-static force capacities of 250-300kN. These devices have been subjected to peak input velocities of approximately 200mm/s, producing peak resistive forces of approximately 350kN. The devices show stable hysteretic performance, with slight force reduction during high-speed testing due to heat build-up and softening of the lead working material. This force reduction is recovered following cyclic loading as heat is dissipated and the lead hardens again. The devices are shown to be only weakly velocity dependent, an advantage in that they do not deliver large forces to the connecting elements and surrounding structure if larger than expected response velocities occur. This high-speed testing is an important step towards uptake as it provides important information to designers.
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
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 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.
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.
Environmental stress and disturbance can affect the structure and functioning of marine ecosystems by altering their physical, chemical and biological features. In estuaries, benthic invertebrate communities play important roles in structuring sediments, influencing primary production and biogeochemical flux, and occupying key food web positions. Stress and disturbance can reduce species diversity, richness and abundance, with ecological theory predicting that biodiversity will be at its lowest soon after a disturbance with assemblages dominated by opportunistic species. The Avon-Heathcote Estuary in Christchurch New Zealand has provided a novel opportunity to examine the effects of stress, in the form of eutrophication, and disturbance, in the form of cataclysmic earthquake events, on the structure and functioning of an estuarine ecosystem. For more than 50 years, large quantities (up to 500,000m3/day) of treated wastewater were released into this estuary but in March 2010 this was diverted to an ocean outfall, thereby reducing the nutrient loading by around 90% to the estuary. This study was therefore initially focussed on the reversal of eutrophication and consequent effects on food web structure in the estuary as it responded to lower nutrients. In 2011, however, Christchurch was struck with a series of large earthquakes that greatly changed the estuary. Massive amounts of liquefied sediments, covering up to 65% of the estuary floor, were forced up from deep below the estuary, the estuary was tilted by up to a 50cm rise on one side and a corresponding drop on the other, and large quantities of raw sewage from broken wastewater infrastructure entered the estuary for up to nine months. This study was therefore a test of the potentially synergistic effects of nutrient reduction and earthquake disturbance on invertebrate communities, associated habitats and food web dynamics. Because there was considerable site-to-site heterogeneity in the estuary, the sites in this study were selected to represent a eutrophication gradient from relatively “clean” (where the influence of tidal flows was high) to highly impacted (near the historical discharge site). The study was structured around these sites, with components before the wastewater diversion, after the diversion but before the earthquakes, and after the earthquakes. The eutrophication gradient was reflected in the composition and isotopic chemistry of primary producer and invertebrate communities and the characteristics of sediments across the sample sites. Sites closest to the former wastewater discharge pipe were the most eutrophic and had cohesive organic -rich, fine sediments and relatively depauperate communities dominated by the opportunistic taxa Capitellidae. The less-impacted sites had coarser, sandier sediments with fewer pollutants and far less organic matter than at the eutrophic sites, relatively high diversity and lower abundances of micro- and macro-algae. Sewage-derived nitrogen had became incorporated into the estuarine food web at the eutrophic sites, starting at the base of the food chain with benthic microalgae (BMA), which were found to use mostly sediment-derived nitrogen. Stable isotopic analysis showed that δ13C and δ15N values of most food sources and consumers varied spatially, temporally and in relation to the diversion of wastewater, whereas the earthquakes did not appear to affect the overall estuarine food web structure. This was seen particularly at the most eutrophic site, where isotopic signatures became more similar to the cleaner sites over two-and-a-half years after the diversion. New sediments (liquefaction) produced by the earthquakes were found to be coarser, have lower concentrations of heavy metals and less organic matter than old (existing) sediments. They also had fewer macroinvertebrate inhabitants initially after the earthquakes but most areas recovered to pre-earthquake abundance and diversity within two years. Field experiments showed that there were higher amounts of primary production and lower amounts of nutrient efflux from new sediments at the eutrophic sites after the earthquakes. Primary production was highest in new sediments due to the increased photosynthetic efficiency of BMA resulting from the increased permeability of new sediments allowing increased light penetration, enhanced vertical migration of BMA and the enhanced transport of oxygen and nutrients. The reduced efflux of NH4-N in new sediments indicated that the capping of a large portion of eutrophic old sediments with new sediments had reduced the release of legacy nutrients (originating from the historical discharge) from the sediments to the overlying water. Laboratory experiments using an array of species and old and new sediments showed that invertebrates altered levels of primary production and nutrient flux but effects varied among species. The mud snail Amphibola crenata and mud crab Austrohelice crassa were found to reduce primary production and BMA biomass through the consumption of BMA (both species) and its burial from bioturbation and the construction of burrows (Austrohelice). In contrast, the cockle Austrovenus stutchburyi did not significantly affect primary production and BMA biomass. These results show that changes in the structure of invertebrate communities resulting from disturbances can also have consequences for the functioning of the system. The major conclusions of this study were that the wastewater diversion had a major effect on food web dynamics and that the large quantities of clean and unpolluted new sediments introduced to the estuary during the earthquakes altered the recovery trajectory of the estuary, accelerating it at least throughout the duration of this study. This was largely through the ‘capping’ effect of the new liquefied, coarser-grained sediments as they dissipated across the estuary and covered much of the old organic-rich eutrophic sediments. For all aspects of this study, the largest changes occurred at the most eutrophic sites; however, the surrounding habitats were important as they provided the context for recovery of the estuary, particularly because of the very strong influence of sediments, their biogeochemistry, microalgal and macroalgal dynamics. There have been few studies documenting system level responses to eutrophication amelioration and to the best on my knowledge there are no other published studies examining the impacts of large earthquakes on benthic communities in an estuarine ecosystem. This research gives valuable insight and advancements in the scientific understanding of the effects that eutrophication recovery and large-scale disturbances can have on the ecology of a soft-sediment ecosystem.