Search

found 153 results

Images for effect; more images...

Research papers, University of Canterbury Library

In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.

Research papers, University of Canterbury Library

Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.

Research papers, University of Canterbury Library

The Canterbury earthquakes of 2010 and 2011 have shone the spotlight on a number of tax issues. These issues, and in particular lessons learned from them, will be relevant for revenue authorities, policymakers and taxpayers alike in the broader context of natural disasters. Issues considered by this paper include the tax treatment of insurance monies. For example, building owners will receive pay-outs for destroyed assets and buildings which have been depreciated. Where the insurance payment is more than the adjusted tax value, there will be a taxable "gain on sale" (or depreciation recovery income). If the building owner uses those insurance proceeds to purchase a replacement asset, legislative amendments specifically enacted following the earthquakes provide that rollover relief of the depreciation recovery income is available. The tax treatment of expenditure to seismically strengthen a building is another significant issue faced by building owners. Case law has determined that this expenditure will usually be capital expenditure. In the past such costs could be capitalised to the building and depreciated accordingly. However, since the 2011-2012 income year owners have been prohibited from claiming depreciation on buildings and therefore currently no deduction is available for such strengthening expenditure (whether immediate or deferred). This has significant potential implications for landlords throughout New Zealand facing significant seismic retrofit costs. Incentives, or some form of financial support, whether delivered through the tax system or some other mechanism may be required. International Financial Reporting Standards (IFRS) require insurance proceeds, including reimbursement for expenditure of a capital nature, be reported as income while expenditure itself is not recorded as a current period expense. This has the effect of overstating current income and creating a larger variation between reported income for accounting and taxation purposes. Businesses have obligations to maintain certain business records for tax purposes. Reconstructing records destroyed by a natural disaster depends on how the information was originally stored. The earthquakes have demonstrated the benefits of ‘off-site’ (outside Canterbury) storage, in particular electronic storage. This paper considers these issues and the Inland Revenue Department (Inland Revenue) Standard Practice Statement which deals with inter alia retention of business records in electronic format and offshore record storage. Employer provided accommodation is treated as income to the benefitting employee. A recent amendment to the Income Tax Act 2007 retrospectively provides that certain employer provided accommodation is exempt from tax. The time aspect of these rules is extended where the employee is involved in the Canterbury rebuild and comes from outside the region.

Research papers, University of Canterbury Library

There is a relationship between inelastic deformation and energy dissipation in structures that are subjected to earthquake ground motions. Thus, if seismic energy dissipation can be achieved by means of a separate non-load bearing supplementary damping system, the load bearing structure can remain elastic with continuing serviceability following the design level earthquake. This research was carried out to investigate the advantages of using added damping in structures. The control system consists of passive friction dampers called ring spring dampers installed in the ground floor of the structure using a tendon to transmit the forces to the other parts of the structure. The ring springs dampers are friction devices consisting of inner and outer ring elements assembled to form a spring stack. External load applied to the spring produces sliding action across mating ring interfaces. The damping forces generated by the dampers and transferred in the supplemental system to the structure by the tendon and horizontal links oppose the internal loads. A four storey-two bay steel frame structure was used in the study. Experimental and analytical studies to investigate the effectiveness of a supplemental control system are presented. The model was subjected to a series of earthquake simulations on the shaking table in the Structural Laboratory of the Civil Engineering Department, at the University of Canterbury. The earthquake simulation tests have been performed on the structure both with and without the supplemental control system. The earthquake simulations were a series of gradually increasing intensity replications of two commonly used earthquake records. This thesis includes detailed description of the structural model, the supplemental control system, the ring springs dampers and the data obtained during the testing. Analyses were then carried out on a twelve storey framed structure to investigate the possible tendon arrangements and the size and type of dampers required to control the response of a real building. Guidelines for determining the appropriateness of including a supplemental damping system have been investigated. The main features of the supplemental control system adopted in this research are: • It is a passive control system with extreme reliability and having no dependence on external power sources to effect the control action. These power sources may not be available during a major earthquake. • Ring springs are steel friction devices capable of absorbing large amounts of input energy. No liquid leakage can occur and minimal maintenance is required for the ring spring dampers. • With a damper-tendon system, the distribution of the dampers throughout the structure is not so critical. Only one or two dampers are used to produce the damping forces needed, and forces are then transferred to the rest of the building by the tendon system. • It is a relatively inexpensive control system with a long useful life.

Research papers, University of Canterbury Library

In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.

Research papers, University of Canterbury Library

Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.

Research papers, University of Canterbury Library

In recent years, rocking isolation has become an effective approach to improve seismic performance of steel and reinforced concrete structures. These systems can mitigate structural damage through rigid body displacement and thus relatively low requirements for structural ductility, which can significantly improve seismic resilience of structures and reduce repairing costs after strong earthquakes. A number of base rocking structural systems with only a single rocking interface have been proposed. However, these systems can have significant high mode effect for high rise structures due to the single rocking interface. This RObust BUilding SysTem (ROBUST) project is a collaborative China-New Zealand project sponsored by the International Joint Research Laboratory of Earthquake Engineering (ILEE), Tongji University, and a number of agencies and universities within New Zealand including the BRANZ, Comflor, Earthquake Commission, HERA, QuakeCoRE, QuakeCentre, University of Auckland, and the University of Canterbury. A number of structural configurations will be tested [1, 2], and non-structural elements including ceilings, infilling walls, glazed curtain walls, precast concrete panels, piping system will also be tested in this project [3]. Within this study, a multiple rocking column steel structural system was proposed and investigated mainly by Tongji team with assistance of NZ members. The concept of rocking column system initiates from the structure of Chinese ancient wooden pagoda. In some of Chinese wooden pagodas, there are continuous core columns hanged only at the top of each pagoda, which is not connected to each stories. This core column can effectively avoid collapse of the whole structure under large storey drifts. Likewise, there are also central continuous columns in the newly proposed steel rocking column system, which can avoid weak story failure mechanism and make story drifts more uniform. In the proposed rocking column system, the structure can switch between an elastic rigidly connected moment resisting frame and a controlled rocking column system when subjected to strong ground motion excitations. The main seismic energy can be dissipated by asymmetric friction beam–column connections, thereby effectively reducing residual displacement of the structure under seismic loading without causing excessive damage to structural members. Re–centering of the structure is provided not only by gravity load carried by rocking columns, but also by mould coil springs. To investigate dynamic properties of the proposed system under different levels of ground excitations, a full-scale threestory steel rocking column structural system with central continuous columns is to be tested using the International joint research Laboratory of Earthquake Engineering (ILEE) facilities, Shanghai, China and an analytical model is established. A finite element model is also developed using ABAQUS to simulate the structural dynamic responses. The rocking column system proposed in this paper is shown to produce resilient design with quick repair or replacement.

Research papers, University of Canterbury Library

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.

Research Papers, Lincoln University

This research investigates creativity in a post-disaster setting. The data explore creativity at the intersection of the affected community of Christchurch, New Zealand and the social processes that followed the earthquakes of 2010 - 2012. Personal and contextual influences on creative ideas implemented for community or commercial benefit are also examined. Viewed as creative, unique approaches to post-disaster problem solving were celebrated locally, nationally and internationally (Bergman, 2014; Wesener, 2015; Cloke & Conradson, 2018). Much has been written about creativity, particularly creativity in organisations and in business. However, little is known with regards to who creates after a disaster, why individuals choose to do so and what impact the post-disaster context has on their creative activity. This exploratory study draws on the literature from the fields of creativity, disasters, psychology, sociology and entrepreneurship to interpret first-hand accounts of people who acted on creative ideas in a physically and socially altered environment. A mixed method - albeit predominantly qualitative - approach to data gathering was adopted that included interviews (n=45) with participants who had been the primary drivers of creative ideas implemented in Christchurch after September 2010 – the first major (7.1 magnitude) earthquake in a prolonged sequence of thousands of aftershocks. Key findings include that a specific type of creativity results from the ‘collision’ between individuals and social processes activated by a disaster situation. This type of creativity could be best categorised as ‘little c’ or socially adaptive and emerges through a prosocial filter. There is wide consensus amongst creativity researchers - principally social psychologists - that for output to be considered creative it must be both novel and useful (Runco & Jaegar, 2012). There is greater tolerance for the novelty component after a disaster as novelty itself has greater utility, either as a distraction or because alternatives are few. Existing creativity models show context as input – an additional component of the creative process – but after a disaster the event itself becomes the catalyst for social processes that result in the creativity seen. Most participants demonstrated characteristics commonly associated with creativity and could be categorised as either a ‘free thinker’ and/or an ‘opportunist’. Some appear preadapted to create and thrive in unstable circumstances. Findings from participants’ completion of a Ten Item Personality Inventory (TIPI) showed an apparent reduced need for extraversion in relation to implementing creative ventures in society. This factor, along with higher levels of agreeableness may indicate a potentially detrimental effect on the success of creative ideas established after a disaster, despite earnest intentions. Three new models are presented to illustrate the key findings of this study. The models imply that disasters enhance both the perceived value of creativity and the desire to act creatively for prosocial ends. The models also indicate that these disaster influenced changes are likely to be temporary.

Research papers, University of Canterbury Library

Oblique-convergent plate collision between the Pacific and Australian plates across the South Island has resulted in shallow, upper crustal earthquake activity and ground surface deformation. In particular the Porters Pass - Amberley Fault Zone displays a complex hybrid zone of anastomosing dextral strike-slip and thrust/reverse faulting which includes the thrust/reverse Lees Valley Fault Zone and associated basin deformation. There is a knowledge gap with respect to the paleoseismicity of many of the faults in this region including the Lees Valley Fault Zone. This study aimed to investigate the earthquake history of the fault at a selected location and the structural and geomorphic development of the Lees Valley Fault Zone and eastern rangefront. This was investigated through extensive structural and geomorphic mapping, GPS field surveying, vertical aerial photo interpretation, analysis of Digital Elevation Models, paleoseismic trenching and optically stimulated luminescence dating. This thesis used a published model for tectonic geomorphology development of mountain rangefronts to understand the development of Lees Valley. Rangefront geomorphology is investigated through analysis of features such as rangefront sinuosity and faceted spurs and indicates the recently active and episodic nature of the uplifted rangefront. Analysis of fault discontinuity, fault splays, distribution of displacement, fault deformation zone and limited exposure of bedrock provided insight into the complex structure of the fault zone. These observations revealed preserved, earlier rangefronts, abandoned and uplifted within the eastern ranges, indicating a basinward shift in focus of faulting and an imbricate thrust wedge development propagating into the footwall of the fault zone and along the eastern ranges of Lees Valley. Fault scarp deformation analysis indicated multiple events have produced the deformation present preserved by the active fault trace in the northern valley. Vertical deformation along this scarp varied with a maximum of 11.5 m and an average of 5 m. Field mapping revealed fan surfaces of various ages have been offset and deformed, likely during the Holocene, based on expected relative surface ages. Geomorphic and structural mapping highlighted the effect of cross-cutting and inherited structures on the Lees Valley Fault, resulting in a step-over development in the centre of the eastern range-bounding trace. Paleoseismic trenching provided evidence of at least two earthquakes, which were constrained to post 21.6 ± 2.3 ka by optically stimulated luminescence dating. Single event displacements (1.48 ± 0.08 m), surface rupture earthquake magnitudes (Mw 6.7 ± 0.1, with potential to produce ≥ 7.0), and a minimum recurrence interval (3.6 ± 0.3 ka) indicated the Lees Valley Fault is an active structure capable of producing significant earthquake events. Results from this study indicate that the Lees Valley Fault Zone accommodates an important component of the Porters Pass - Amberley Fault Zone deformation and confirms the fault as a source of potentially damaging, peak ground accelerations in the Canterbury region. Remnants of previous rangefronts indicate a thrust wedge development of the Lees Valley Fault Zone and associated ranges that can potentially be used as a model of development for other thrust-fault bounded basins.

Research papers, University of Canterbury Library

Between 2010 and 2011, Canterbury experienced a series of four large earthquake events with associated aftershocks which caused widespread damage to residential and commercial infrastructure. Fine grained and uncompacted alluvial soils, typical to the Canterbury outwash plains, were exposed to high peak ground acceleration (PGA) during these events. This rapid increase in PGA induced cyclic strain softening and liquefaction in the saturated, near surface alluvial soils. Extensive research into understanding the response of soils in Canterbury to dynamic loading has since occurred. The Earthquake Commission (EQC), the Ministry of Business and Employment (MBIE), and the Christchurch City Council (CCC) have quantified the potential hazards associated with future seismic events. Theses bodies have tested numerous ground improvement design methods, and subsequently are at the forefront of the Canterbury recovery and rebuild process. Deep Soil Mixing (DSM) has been proven as a viable ground improvement foundation method used to enhance in situ soils by increasing stiffness and positively altering in situ soil characteristics. However, current industry practice for confirming the effectiveness of the DSM method involves specific laboratory and absolute soil test methods associated with the mixed column element itself. Currently, the response of the soil around the columns to DSM installation is poorly understood. This research aims to understand and quantify the effects of DSM columns on near surface alluvial soils between the DSM columns though the implementation of standardised empirical soil test methods. These soil strength properties and ground improvement changes have been investigated using shear wave velocity (Vs), soil behaviour and density response methods. The results of the three different empirical tests indicated a consistent improvement within the ground around the DSM columns in sandier soils. By contrast, cohesive silty soils portrayed less of a consistent response to DSM, although still recorded increases. Generally, within the tests completed 50 mm from the column edge, the soil response indicated a deterioration to DSM. This is likely to be a result of the destruction of the soil fabric as the stress and strain of DSM is applied to the un‐mixed in situ soils. The results suggest that during the installation of DSM columns, a positive ground effect occurs in a similar way to other methods of ground improvement. However, further research, including additional testing following this empirical method, laboratory testing and finite 2D and 3D modelling, would be useful to quantify, in detail, how in situ soils respond and how practitioners should consider these test results in their designs. This thesis begins to evaluate how alluvial soils tend to respond to DSM. Conducting more testing on the research site, on other sites in Christchurch, and around the world, would provide a more complete data set to confirm the results of this research and enable further evaluation. Completing this additional research could help geotechnical DSM practitioners to use standardised empirical test methods to measure and confirm ground improvement rather than using existing test methods in future DSM projects. Further, demonstrating the effectiveness of empirical test methods in a DSM context is likely to enable more cost effective and efficient testing of DSM columns in future geotechnical projects.

Research papers, Victoria University of Wellington

“One of the most basic and fundamental questions in urban master planning and building regulations is ‘how to secure common access to sun, light and fresh air?” (Stromann-Andersen & Sattrup, 2011).  Daylighting and natural ventilation can have significant benefits in office buildings. Both of these ‘passive’ strategies have been found to reduce artificial lighting and air-conditioning energy consumption by as much as 80% (Ministry for the Environment, 2008); (Brager, et al., 2007). Access to daylight and fresh air can also be credited with improved occupant comfort and health, which can lead to a reduction of employee absenteeism and an increase of productivity (Sustainability Victoria, 2008).  In the rebuild of Christchurch central city, following the earthquakes of 2010 and 2011, Cantabrians have expressed a desire for a low-rise, sustainable city, with open spaces and high performance buildings (Christchurch City Council, 2011). With over 80% of the central city being demolished, a unique opportunity to readdress urban form and create a city that provides all buildings with access to daylight and fresh air exists.  But a major barrier to wide-spread adoption of passive buildings in New Zealand is their dependence on void space to deliver daylight and fresh air – void space which could otherwise be valuable built floor space. Currently, urban planning regulations in Christchurch prioritize density, allowing and even encouraging low performance compact buildings.  Considering this issue of density, this thesis aimed to determine which urban form and building design changes would have the greatest effect on building performance in Central City Christchurch.  The research proposed and parametrically tested modifications of the current compact urban form model, as well as passive building design elements. Proposed changes were assessed in three areas: energy consumption, indoor comfort and density. Three computer programs were used: EnergyPlus was the primary tool, simulating energy consumption and thermal comfort. Radiance/Daysim was used to provide robust daylighting calculations and analysis. UrbaWind enabled detailed consideration of the urban wind environment for reliable natural ventilation predictions.  Results found that, through a porous urban form and utilization of daylight and fresh air via simple windows, energy consumption could be reduced as much as 50% in buildings. With automatic modulation of windows and lighting, thermal and visual comfort could be maintained naturally for the majority of the occupied year. Separation of buildings by as little as 2m enabled significant energy improvements while having only minimal impact on individual property and city densities.  Findings indicated that with minor alterations to current urban planning laws, all buildings could have common access to daylight and fresh air, enabling them to operate naturally, increasing energy efficiency and resilience.

Research papers, University of Canterbury Library

Non-structural elements (NSEs) have frequently proven to contribute to significant losses sustained from earthquakes in the form of damage, downtime, injury and death. In New Zealand (NZ), the 2010 and 2011 Canterbury Earthquake Sequence (CES), the 2013 Seddon and Cook Strait earthquake sequence and the 2016 Kaikoura earthquake were major milestones in this regard as significant damage to building NSEs both highlighted and further reinforced the importance of NSE seismic performance to the resilience of urban centres. Extensive damage in suspended ceilings, partition walls, façades and building services following the CES was reported to be partly due to erroneous seismic design or installation or caused by intervening elements. Moreover, the low-damage solutions developed for structural systems sometimes allow for relatively large inter-story drifts -compared to conventional designs- which may not have been considered in the seismic design of NSEs. Having observed these shortcomings, this study on suspended ceilings was carried out with five main goals: i) Understanding the seismic performance of the system commonly used in NZ; ii) Understanding the transfer of seismic design actions through different suspended ceiling components, iii) Investigating potential low-damage solutions; iii) Evaluating the compatibility of the current ceiling system with other low-damage NSEs; and iv) Investigating the application of numerical analysis to simulate the response of ceiling systems. The first phase of the study followed a joint research work between the University of Canterbury (UC) in NZ, and the Politecnico Di Milano, in Italy. The experimental ceiling component fragility curves obtained in this existing study were employed to produce analytical fragility curves for a perimeter-fixed ceiling of a given size and weight, with grid acceleration as the intensity measure. The validity of the method was proven through comparisons between this proposed analytical approach with the recommended procedures in proprietary products design guidelines, as well as experimental fragility curves from other studies. For application to engineering design practice, and using fragility curves for a range of ceiling lengths and weights, design curves were produced for estimating the allowable grid lengths for a given demand level. In the second phase of this study, three specimens of perimeter-fixed ceilings were tested on a shake table under both sinusoidal and random floor motion input. The experiments considered the relationship between the floor acceleration, acceleration of the ceiling grid, the axial force induced in the grid members, and the effect of boundary conditions on the transfer of these axial forces. A direct correlation was observed between the axial force (recorded via load cells) and the horizontal acceleration measured on the ceiling grid. Moreover, the amplification of floor acceleration, as transferred through ceiling components, was examined and found (in several tests) to be greater than the recommended factor for the design of ceilings provided in the NZ earthquake loadings standard NZS1170.5. However, this amplification was found to be influenced by the pounding interactions between the ceiling grid members and the tiles, and this amplification diminished considerably when the high frequency content was filtered out from the output time histories. The experiments ended with damage in the ceiling grid connection at an axial force similar to the capacity of these joints previously measured through static tests in phase one. The observation of common forms of damage in ceilings in earthquakes triggered the monotonic experiments carried out in the third phase of this research with the objective of investigating a simple and easily applicable mitigation strategy for existing or new suspended ceilings. The tests focused on the possibility of using proprietary cross-shaped clip elements ordinarily used to provide seismic gap as a strengthening solution for the weak components of a ceiling. The results showed that the solution was effective under both tension and compression loads through increasing load bearing capacity and ductility in grid connections. The feasibility of a novel type of suspended ceiling called fully-floating ceiling system was investigated through shaking table tests in the next phase of this study with the main goal of isolating the ceiling from the surrounding structure; thereby arresting the transfer of associated seismic forces from the structure to the ceiling. The fully-floating ceiling specimen was freely hung from the floor above lacking any lateral bracing and connections with the perimeter. Throughout different tests, a satisfactory agreement between the fully-floating ceiling response and simple pendulum theory was demonstrated. The addition of isolation material in perimeter gaps was found effective in inducing extra damping and protecting the ceiling from pounding impact; resulting in much reduced ceiling displacements and accelerations. The only form of damage observed throughout the random floor motion tests and the sinusoidal tests was a panel dislodgement observed in a test due to successive poundings between the ceiling specimen and the surrounding beams at resonant frequencies. Partition walls as the first effective NSE in direct interaction with ceilings were the topic of the final experimental phase. Low-damage drywall partitions proposed in a previous study in the UC were tested with two common forms of suspended ceiling: braced and perimeter-fixed. The experiments investigated the in-plane and out-of-plane performance of the low-damage drywall partitions, as well as displacement compatibility between these walls and the suspended ceilings. In the braced ceiling experiment, where no connection was made between ceiling grids and surrounding walls no damage in the grid system or partitions was observed. However, at high drift values panel dislodgement was observed on corners of the ceiling where the free ends of grids were not restrained against spreading. This could be prevented by framing the grid ends using a perimeter angle that is riveted only to the grid members while keeping sufficient clearance from the perimeter walls. In the next set of tests with the perimeter-fixed ceiling, no damage was observed in the ceiling system or the drywalls. Based on the results of the experiments it was concluded that the tested ceiling had enough flexibility to accommodate the relative displacement between two perpendicular walls up to the inter-storey drifts achieved. The experiments on perimeter-fixed ceilings were followed by numerical simulations of the performance of these ceilings in a finite element model developed in the structural analysis software, SAP2000. This model was relatively simple and easy to develop and was able to replicate the experimental results to a reasonable degree. Filtering was applied to the experimental output to exclude the effect of high frequency noise and tile-grid impact. The developed model generally simulated the acceleration responses well but underestimated the peak ceiling grid accelerations. This was possibly because the peak values in time histories were affected by impact occurring at very short periods. The model overestimated the axial forces in ceiling grids which was assumed to be caused by the initial assumptions made about the tributary area or constant acceleration associated with each grid line in the direction of excitation. Otherwise, the overall success of the numerical modelling in replicating the experimental results implies that numerical modelling using conventional structural analysis software could be used in engineering practice to analyse alternative ceiling geometries proposed for application to varying structural systems. This however, needs to be confirmed through similar analyses on other ceiling examples from existing instrumented buildings during real earthquakes. As the concluding part of this research the final phase addressed the issues raised following the review of existing ceiling standards and guidelines. The applicability of the research findings to current practice and their implications were discussed. Finally, an example was provided for the design of a suspended ceiling utilising the new knowledge acquired in this research.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.

Research papers, University of Canterbury Library

Cats all over the world hunt wild animals and can contribute to the extinction of threatened species. In New Zealand, around half of all households have at least one cat. When cats live close to a natural area, such as a wetland, they may have impacts on native species. A previous study on the movements and hunting behaviour of domestic (house) cats around the Travis Wetland, Christchurch, New Zealand during 2000-2001 raised concerns about the effects of cats on the local skink population, as skinks were a frequent prey item. My study is a comparison to the prior study, to determine if impacts have changed alongside changes in human populations in the area post-earthquake. The domestic cat population in the area was estimated by a household survey in March-April 2018. For a 6 month period from March-September 2018, 26 households recorded prey brought home by their 41 cats. During April-July 2018, 14 cats wore Global Positioning System (GPS) devices for 7 days each to track their movements. Skink abundance was measured with pitfall trapping over 20 days in February 2018. There were more households in the area in 2018 than there were in 2000, but the numbers of cats had decreased. In the 196 ha study area around Travis Wetland, the domestic cat population was estimated at 429 cats, down from the previous 494. Most of the cats were free roaming, but the majority had been desexed and many were mostly seen at home. A total of 42 prey items were reported from 26 households and 41 cats over 6 months. Of these, 62% were rodents, 26% were exotic birds, and 12% were native birds. There were no native skinks, other mammals, or other vertebrates such as fish and amphibians (invertebrates were not included in this study). Eight male and six female cats were tracked by GPS. Home range sizes for the 100% minimum convex polygons (MCPs) ranged from 1.34 to 9.68 ha (mean 4.09 ha, median 3.54 ha). There were 9/14 (64%) cats that entered the edge of the wetland. Males had significantly larger home range areas at night and in general compared with females. However, age and distance of the cat’s household to the wetland did not have a significant effect on home range size and there was no significant correlation between home range size and prey retrieved. In 20 days of skink trapping, 11 Oligosoma polychroma were caught. The estimated catch rate was not significantly different from an earlier study on skink abundance in Travis Wetland. The apparently low abundance of skinks may have been the result of successful wetland restoration creating less suitable skink habitat, or of other predators other than cats. In the future, increased education should be provided to the public about ways to increase wildlife in their area. This includes creating lizard friendly habitat in their gardens and increasing management for cats. Generally, only a few cats bring home prey often, and these select cats should be identified in initial surveys and included in further studies. In New Zealand, until management programmes can target all predators in urban areas, domestic cats could stay out at night and inside during the day to help decrease the abundance of rodents at night and reduce the number of birds and lizards caught during the day.

Research Papers, Lincoln University

Environmental assessment in New Zealand is governed by the provisions of the Resource Management Act (RMA) 1991. The Act requires persons wishing to undertake certain activities to apply for resource consent from their local or regional council - a procedure termed the Resource Consent Process. The key component of a resource consent application is an Assessment of Environmental Effects (AEE) report; a statement of the environmental effects of a proposed activity. Problems arise when environmental assessments are complicated by uncertain and abnormal circumstances such as natural hazards. Natural hazards (including earthquakes, floods, tsunami, and coastal erosion) can be catastrophic to an environment. If hazards are not avoided or successfully mitigated, they can result in serious consequences to proposed development and to the environment which the proposal relates. The aim of this study is to assess the adequacy of the resource consent process (as outlined in the Resource Management Act 1991) for dealing with proposed development affected by natural hazards. This study reviews the context of the resource consent process for dealing with natural hazards to identify potential issues in the assessment process. Guidance criteria for assessing natural hazards (termed Natural Hazard Assessment) are developed to evaluate against two resource consent applications affected by natural hazards. The findings of the consent process review and case study evaluation are discussed to determine the adequacy of the consent process for dealing with natural hazards. From the review of the consent process it was evident that the process has a number of problems for accommodating natural hazards into the assessment. Although many important traits are provided for in the process, such traits are not always reflected in environmental assessments. Evaluation of two resource consent applications against the process of Natural Hazard Assessment (NHA) showed that these consent applications did not adequately detail key information relating to natural hazards. Many problems evident in these applications were not amended by the Consent Authorities in the review process and subsequently consent was granted to information-deficient applications. Problematic issues identified in this study include: • A distinct lack of guidance (legal or otherwise) for the applicant and Consent Authority regarding the boundaries of inclusion of an effect; • Deficiencies in planning documents are reflected in AEE reports, the review of the consent application and in the end-decision; • Under-utilisation of "experts" throughout the consent process; • Minimal identification and account for the degree of uncertainty throughout the consent process; • Resource consents are being granted even though information in consent applications, and the means for assessing the information is deficient. These issues reflect that decisions are not being made based on all elements involved in a potential hazard. Subsequently, the resource consent process is not adequate for dealing with all aspects of natural hazards. The Natural Hazard Assessment process provides educated assessment criteria to assess development affected by natural hazards. By accounting for the problems evident in the consent process, the introduction of a three-tier identification, risk and vulnerability assessment, and evaluation process to account for uncertainties, Natural Hazard Assessment provides a platform for a thorough assessment of natural hazards. The application of the principles of Natural Hazard Assessment to the consent applications affected by natural hazards showed that many key issues were not covered in the assessment under the consent process. The nature of a natural event is that one may not occur in a given region over many lifetimes, however they will occur at some stage and planning and environmental assessment needs to provide for the associated hazards. Implementation of Natural Hazard Assessment is needed to help provide answers for the problems experienced in the resource consent process. Natural Hazard Assessment would allow decision-makers to make informed judgements on the situation at hand, leading to better planning and land-use options. Change to current practice is needed, as following the current path of environmental assessment will be the hazard in the end.

Research papers, University of Canterbury Library

Previous earthquakes demonstrated destructive effects of soil-structure interaction on structural response. For example, in the 1970 Gediz earthquake in Turkey, part of a factory was demolished in a town 135 km from the epicentre, while no other buildings in the town were damaged. Subsequent investigations revealed that the fundamental period of vibration of the factory was approximately equal to that of the underlying soil. This alignment provided a resonance effect and led to collapse of the structure. Another dramatic example took place in Adapazari, during the 1999 Kocaeli earthquake where several foundations failed due to either bearing capacity exceedance or foundation uplifting, consequently, damaging the structure. Finally, the Christchurch 2012 earthquakes have shown that significant nonlinear action in the soil and soil-foundation interface can be expected due to high levels of seismic excitation and spectral acceleration. This nonlinearity, in turn, significantly influenced the response of the structure interacting with the soil-foundation underneath. Extensive research over more than 35 years has focused on the subject of seismic soil-structure interaction. However, since the response of soil-structure systems to seismic forces is extremely complex, burdened by uncertainties in system parameters and variability in ground motions, the role of soil-structure interaction on the structural response is still controversial. Conventional design procedures suggest that soil-structure interaction effects on the structural response can be conservatively ignored. However, more recent studies show that soil-structure interaction can be either beneficial or detrimental, depending on the soil-structure-earthquake scenarios considered. In view of the above mentioned issues, this research aims to utilise a comprehensive and systematic probabilistic methodology, as the most rational way, to quantify the effects of soil-structure interaction on the structural response considering both aleatory and epistemic uncertainties. The goal is achieved by examining the response of established rheological single-degree-of-freedom systems located on shallow-foundation and excited by ground motions with different spectral characteristics. In this regard, four main phases are followed. First, the effects of seismic soil-structure interaction on the response of structures with linear behaviour are investigated using a robust stochastic approach. Herein, the soil-foundation interface is modelled by an equivalent linear cone model. This phase is mainly considered to examine the influence of soil-structure interaction on the approach that has been adopted in the building codes for developing design spectrum and defining the seismic forces acting on the structure. Second, the effects of structural nonlinearity on the role of soil-structure interaction in modifying seismic structural response are studied. The same stochastic approach as phase 1 is followed, while three different types of structural force-deflection behaviour are examined. Third, a systematic fashion is carried out to look for any possible correlation between soil, structural, and system parameters and the degree of soil-structure interaction effects on the structural response. An attempt is made to identify the key parameters whose variation significantly affects the structural response. In addition, it is tried to define the critical range of variation of parameters of consequent. Finally, the impact of soil-foundation interface nonlinearity on the soil-structure interaction analysis is examined. In this regard, a newly developed macro-element covering both material and geometrical soil-foundation interface nonlinearity is implemented in a finite-element program Raumoko 3D. This model is then used in an extensive probabilistic simulation to compare the effects of linear and nonlinear soil-structure interaction on the structural response. This research is concluded by reviewing the current design guidelines incorporating soil-structure interaction effects in their design procedures. A discussion is then followed on the inadequacies of current procedures based on the outcomes of this study.

Research papers, University of Canterbury Library

The increase of the world's population located near areas prone to natural disasters has given rise to new ‘mega risks’; the rebuild after disasters will test the governments’ capabilities to provide appropriate responses to protect the people and businesses. During the aftermath of the Christchurch earthquakes (2010-2012) that destroyed much of the inner city, the government of New Zealand set up a new partnership between the public and private sector to rebuild the city’s infrastructure. The new alliance, called SCIRT, used traditional risk management methods in the many construction projects. And, in hindsight, this was seen as one of the causes for some of the unanticipated problems. This study investigated the risk management practices in the post-disaster recovery to produce a specific risk management model that can be used effectively during future post-disaster situations. The aim was to develop a risk management guideline for more integrated risk management and fill the gap that arises when the traditional risk management framework is used in post-disaster situations. The study used the SCIRT alliance as a case study. The findings of the study are based on time and financial data from 100 rebuild projects, and from surveying and interviewing risk management professionals connected to the infrastructure recovery programme. The study focussed on post-disaster risk management in construction as a whole. It took into consideration the changes that happened to the people, the work and the environment due to the disaster. System thinking, and system dynamics techniques have been used due to the complexity of the recovery and to minimise the effect of unforeseen consequences. Based on an extensive literature review, the following methods were used to produce the model. The analytical hierarchical process and the relative importance index have been used to identify the critical risks inside the recovery project. System theory methods and quantitative graph theory have been used to investigate the dynamics of risks between the different management levels. Qualitative comparative analysis has been used to explore the critical success factors. And finally, causal loop diagrams combined with the grounded theory approach has been used to develop the model itself. The study identified that inexperienced staff, low management competency, poor communication, scope uncertainty, and non-alignment of the timing of strategic decisions with schedule demands, were the key risk factors in recovery projects. Among the critical risk groups, it was found that at a strategic management level, financial risks attracted the highest level of interest, as the client needs to secure funding. At both alliance-management and alliance-execution levels, the safety and environmental risks were given top priority due to a combination of high levels of emotional, reputational and media stresses. Risks arising from a lack of resources combined with the high volume of work and the concern that the cost could go out of control, alongside the aforementioned funding issues encouraged the client to create the recovery alliance model with large reputable construction organisations to lock in the recovery cost, at a time when the scope was still uncertain. This study found that building trust between all parties, clearer communication and a constant interactive flow of information, established a more working environment. Competent and clear allocation of risk management responsibilities, cultural shift, risk prioritisation, and staff training were crucial factors. Finally, the post-disaster risk management (PDRM) model can be described as an integrated risk management model that considers how the changes which happened to the environment, the people and their work, caused them to think differently to ease the complexity of the recovery projects. The model should be used as a guideline for recovery systems, especially after an earthquake, looking in detail at all the attributes and the concepts, which influence the risk management for more effective PDRM. The PDRM model is represented in Causal Loops Diagrams (CLD) in Figure 8.31 and based on 10 principles (Figure 8.32) and 26 concepts (Table 8.1) with its attributes.

Research papers, University of Canterbury Library

The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.

Research papers, University of Canterbury Library

The Avon-Heathcote Estuary, located in Christchurch, New Zealand, experienced coseismic deformation as a result of the February 22nd 2011 Christchurch Earthquake. The deformation is reflected as subsidence in the northern area and uplift in the southern area of the Estuary, in addition to sand volcanoes which forced up sediment throughout the floor of the Estuary altering estuary bed height and tidal flow. The first part of the research involved quantifying the change in the modern benthic foraminifera distribution as a result of the coseismic deformation caused by the February 22nd 2011 earthquake. By analysing the taxa present immediately post deformation and then the taxa present 2 years post deformation a comparison of the benthic foraminifera distribution can be made of the pre and post deformation. Both the northern and the southern areas of the Estuary were sampled to establish whether foraminifera faunas migrated landward or seaward as a result of subsidence and uplift experienced in different areas. There was no statistical change in overall species distribution in the two year time period since the coseismic deformation occurred, however, there were some noticeable changes in foraminifera distribution at BSNS-Z3 showing a landward migration of taxa. The changes that were predicted to occur as a result of the deformation of the Estuary are taking longer than expected to show up in the foraminiferal record and a longer time period is needed to establish these changes. The second stage involved establishing the modern distribution of foraminifera at Settlers Reserve in the southern area of the Avon-Heathcote Estuary by detailed sampling along a 160 m transect. Foraminifera are sensitive to environmental parameters, tidal height, grainsize, pH and salinity were recorded to evaluate the effect these parameters have on distribution. Bray-Curtis two-way cluster analysis was primarily used to assess the distribution pattern of foraminifera. The modern foraminifera distribution is comparable to that of the modern day New Zealand brackish-water benthic foraminifera distribution and includes species not yet found in other studies of the Avon-Heathcote Estuary. Differences in sampling techniques and the restricted intertidal marshland area where the transect samples were collected account for some of the differences seen between this model and past foraminifera studies. xiii The final stage involved sampling a 2.20 m core collected from Settlers Reserve and using the modern foraminiferal distribution to establish a foraminiferal history of Settlers Reserve. As foraminifera are sensitive to tidal height they may record past coseismic deformation events and the core was used to ascertain whether record of past coseismic deformation is preserved in Settlers Reserve sediments. Sampling the core for foraminifera, grainsize, trace metals and carbon material helped to build a story of estuary development. Using the modern foraminiferal distribution and the tidal height information collected, a down core model of past tidal heights was established to determine past rates of change. Foraminifera are not well preserved throughout the core, however, a sudden relative rise in sea level is recorded between 0.25 m and 0.85 m. Using trace metal and isotope analysis to develop an age profile, this sea level rise is interpreted to record coseismic subsidence associated with a palaeoseismic event in the early 1900’s. Overall, although the Avon-Heathcote Estuary experienced clear coseismic deformation as a result of the 22nd of February 2011 earthquake, modern changes in foraminiferal distribution cannot yet be tracked, however, past seismic deformation is identified in a core. The modern transect describes the foraminifera distribution which identifies species that have not been identified in the Avon-Heathcote Estuary before. This thesis enhances the current knowledge of the Avon-Heathcote Estuary and is a baseline for future studies.

Research papers, University of Canterbury Library

This report presents the simplified seismic assessment of a case study reinforced concrete (RC) building following the newly developed and refined NZSEE/MBIE guidelines on seismic assessment (NZSEE/MBIE, semi-final draft 26 October 2016). After an overview of the step-by-step ‘diagnostic’ process, including an holistic and qualitative description of the expected vulnerabilities and of the assessment strategy/methodology, focus is given, whilst not limited, to the implementation of a Detailed Seismic Assessment (DSA) (NZSEE/MBIE, 2016c). The DSA is intended to provide a more reliable and consistent outcome than what can be provided by an initial seismic assessment (ISA). In fact, while the Initial Seismic Assessment (ISA), of which the Initial Evaluation Procedure is only a part of, is the more natural and still recommended first step in the overall assessment process, it is mostly intended to be a coarse evaluation involving as few resources as reasonably possible. It is thus expected that an ISA will be followed by a Detailed Seismic Assessment (DSA) not only where the threshold of 33%NBS is not achieved but also where important decisions are intended that are reliant on the seismic status of the building. The use of %NBS (% New Building Standard) as a capacity/demand ratio to describe the result of the seismic assessment at all levels of assessment procedure (ISA through to DSA) is deliberate by the NZSEE/MBIE guidelines (Part A) (NZSEE/MBIE 2016a). The rating for the building needs only be based on the lowest level of assessment that is warranted for the particular circumstances. Discussion on how the %NBS rating is to be determined can be found in Section A3.3 (NZSEE/MBIE 2016a), and, more specifically, in Part B for the ISA (NZSEE/MBIE 2016b) and Part C for the DSA (NZSEE/MBIE 2016c). As per other international approaches, the DSA can be based on several analysis procedures to assess the structural behaviour (linear, nonlinear, static or dynamic, force or displacement-based). The significantly revamped NZSEE 2016 Seismic Assessment Guidelines strongly recommend the use of an analytical (basically ‘by hand’) method, referred to the Simple Lateral Mechanism Analysis (SLaMA) as a first phase of any other numerically-based analysis method. Significant effort has thus been dedicated to provide within the NZSEE 2016 guidelines (NZSEE/MBIE 2016c) a step-by-step description of the procedure, either in general terms (Chapter 2) or with specific reference to Reinforced Concrete Buildings (Chapter 5). More specifically, extract from the guidelines, NZSEE “recommend using the Simple Lateral Mechanism Analysis (SLaMA) procedure as a first step in any assessment. While SLaMA is essentially an analysis technique, it enables assessors to investigate (and present in a simple form) the potential contribution and interaction of a number of structural elements and their likely effect on the building’s global capacity. In some cases, the results of a SLaMA will only be indicative. However, it is expected that its use should help assessors achieve a more reliable outcome than if they only carried out a detailed analysis, especially if that analysis is limited to the elastic range For complex structural systems, a 3D dynamic analysis may be necessary to supplement the simplified nonlinear Simple Lateral Mechanism Analysis (SLaMA).” This report presents the development of a full design example for the the implementation of the SLaMA method on a case study buildings and a validation/comparison with a non-linear static (pushover) analysis. The step-by-step-procedure, summarized in Figure 1, will be herein demonstrated from a component level (beams, columns, wall elements) to a subassembly level (hierarchy of strength in a beam-column joint) and to a system level (frame, C-Wall) assuming initially a 2D behaviour of the key structural system, and then incorporating a by-hand 3D behaviour (torsional effects).

Research papers, University of Canterbury Library

Landslides are significant hazards, especially in seismically-active mountainous regions, where shaking amplified by steep topography can result in widespread landsliding. These landslides present not only an acute hazard, but a chronic hazard that can last years-to-decades after the initial earthquake, causing recurring impacts. The Mw 7.8 Kaikōura earthquake caused more than 20,000 landslides throughout North Canterbury and resulted in significant damage to nationally significant infrastructure in the coastal transport corridor (CTC), isolating Kaikōura from the rest of New Zealand. In the years following, ongoing landsliding triggered by intense rainfall exacerbated the impacts and slowed the recovery process. However, while there is significant research on co-seismic landslides and their initial impacts in New Zealand, little research has explored the evolution of co-seismic landslides and how this hazard changes over time. This research maps landslides annually between 2013 and 2021 to evaluate the changes in pre-earthquake, co-seismic and post-earthquake rates of landsliding to determine how landslide hazard has changed over this time. In particular, the research explores how the number, area, and spatial distribution of landslides has changed since the earthquake, and whether post-earthquake mitigation works have in any way affected the long-term landslide hazard. Mapping of landslides was undertaken using open-source, medium resolution Landsat-8 and Sentinel-2 satellite imagery, with landslides identified visually and mapped as single polygons that capture both the source zone and deposit. Three study areas with differing levels of post-earthquake mitigation are compared: (i) the northern CTC, where the majority of mitigation was in the form of active debris removal; (ii) the southern CTC, where mitigation was primarily via passive protection measures; and (iii) Mount Fyffe, which has had no mitigation works since the earthquake. The results show that despite similar initial impacts during the earthquake, the rate of recovery in terms of landslide rates varies substantially across the three study areas. In Mount Fyffe, the number and area of landslides could take 45 and 22 years from 2021 respectively to return to pre-earthquake levels at the current rate. Comparatively, in the CTC, it could take just 5 years and 3-4 years from 2021 respectively. Notably, the fastest recovery in terms of landslide rates in the CTC was primarily located directly along the transport network, whereas what little recovery did occur in Mount Fyffe appeared to follow no particular pattern. Importantly, recovery rates in the northern CTC were notably higher than in the southern CTC, despite greater co-seismic impacts in the former. Combined, these results suggest the active, debris removal mitigation undertaken in the northern CTC may have had the effect of dramatically reducing the time for landslide rates to return to pre-earthquake levels. The role of slope angle and slope aspect were explored to evaluate if these observations could be driven by local differences in topography. The Mount Fyffe study area has higher slope angles than the CTC as a whole and landslides predominantly occurred on slightly steeper slopes than in the CTC. This may have contributed to the longer recovery times for landsliding in Mount Fyffe due to greater gravitational instability, however the observed variations are minor compared to the differences in recovery rates. In terms of slope aspect, landslides in Mount Fyffe preferentially occurred on north- and south-facing slopes whereas landslides in the CTC preferred the east- and south-facing slopes. The potential role of these differences in landslide recovery remains unclear but may be related to the propagation direction of the earthquake and the tracking direction of post-earthquake ex-tropical cyclones. Finally, landslides in the CTC are observed to be moving further away from the transport network and the number of landslides impacting the CTC decreased significantly since the earthquake. Nevertheless, the potential for further landslide reactivation remains. Therefore, despite the recovery in the CTC, it is clear that there is still risk of the transport network being impacted by further landsliding, at least for the next 3-5 yrs.

Research papers, University of Canterbury Library

Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.

Research papers, University of Canterbury Library

Coastal margins are exposed to rising sea levels that present challenging circumstances for natural resource management. This study investigates a rare example of tectonic displacement caused by earthquakes that generated rapid sea-level change in a tidal lagoon system typical of many worldwide. This thesis begins by evaluating the coastal squeeze effects caused by interactions between relative sea-level (RSL) rise and the built environment of Christchurch, New Zealand, and also examples of release from similar effects in areas of uplift where land reclamations were already present. Quantification of area gains and losses demonstrated the importance of natural lagoon expansion into areas of suitable elevation under conditions of RSL rise and showed that they may be necessary to offset coastal squeeze losses experienced elsewhere. Implications of these spatial effects include the need to provide accommodation space for natural ecosystems under RSL rise, yet other land-uses are likely to be present in the areas required. Consequently, the resilience of these environments depends on facilitating transitions between human land-uses either proactively or in response to disaster events. Principles illustrated by co-seismic sea-level change are generally applicable to climate change adaptation due to the similarity of inundation effects. Furthermore, they highlight the potential role of non-climatic factors in determining the overall trajectory of change. Chapter 2 quantifies impacts on riparian wetland ecosystems over an eight year period post- quake. Coastal wetlands were overwhelmed by RSL rise and recovery trajectories were surprisingly slow. Four risk factors were identified from the observed changes: 1) the encroachment of anthropogenic land-uses, 2) connectivity losses between areas of suitable elevation, 3) the disproportionate effect of larger wetland vulnerabilities, and 4) the need to protect new areas to address the future movement of ecosystems. Chapter 3 evaluates the unique context of shoreline management on a barrier sandspit under sea-level rise. A linked scenario approach was used to evaluate changes on the open coast and estuarine shorelines simultaneously and consider combined effects. The results show dune loss from a third of the study area using a sea-level rise scenario of 1 m over 100 years and with continuation of current land-uses. Increased exposure to natural hazards and accompanying demand for seawalls is a likely consequence unless natural alternatives can be progressed. In contrast, an example of managed retreat following earthquake-induced subsidence of the backshore presents a new opportunity to restart saltmarsh accretion processes seaward of coastal defences with the potential to reverse decades of degradation and build sea-level rise resilience. Considering both shorelines simultaneously highlights the existence of pinch-points from opposing forces that result in small land volumes above the tidal range. Societal adaptation is delicately poised between the paradigms of resisting or accommodating nature and challenged by the long perimeter and confined nature of the sandspit feature. The remaining chapters address the potential for salinity effects caused by tidal prism changes with a focus on the conservation of īnanga (Galaxias maculatus), a culturally important fish that supports New Zealand‘s whitebait fishery. Methodologies were developed to test the hypothesis that RSL changes would drive a shift in the distribution of spawning sites with implications for their management. Chapter 4 describes a new practical methodology for quantifying the total productivity and spatiotemporal variability of spawning sites at catchment scale. Chapter 5 describes the novel use of artificial habitats as a detection tools to help overcome field survey limitations in degraded environments where egg mortality can be high. The results showed that RSL changes resulted in major shifts in spawning locations and these were associated with new patterns of vulnerability due to the continuation of pre-disturbance land-uses. Unexpected findings includes an improved understanding of the spatial relationship between salinity and spawning habitat, and identification of an invasive plant species as important spawning habitat, both with practical management implications. To conclude, the design of legal protection mechanisms was evaluated in relation to the observed habitat shifts and with a focus on two new planning initiatives that identified relatively large protected areas (PAs) in the lower river corridors. Although the larger PAs were better able to accommodate the observed habitat shifts inefficiencies were also apparent due to spatial disparities between PA boundaries and the values requiring protection. To reduce unnecessary trade-offs with other land-uses, PAs of sufficient size to cover the observable spatiotemporal variability and coupled with adaptive capacity to address future change may offer a high effectiveness from a network of smaller PAs. The latter may be informed by both monitoring and modelling of future shifts and these are expected to include upstream habitat migration driven by the identified salinity relationships and eustatic sea-level rise. The thesis concludes with a summary of the knowledge gained from this research that can assist the development of a new paradigm of environmental sustainability incorporating conservation and climate change adaptation. Several promising directions for future research identified within this project are also discussed.

Research papers, University of Canterbury Library

The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.

Research Papers, Lincoln University

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.

Research papers, University of Canterbury Library

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.

Research Papers, Lincoln University

Mixed conifer, beech and hardwood forests are relatively common in Aotearoa/New Zealand, but are not well studied. This thesis investigates the coexistence, regeneration dynamics and disturbance history of a mixed species forest across an environmental gradient of drainage and soil development in north Westland. The aim was to investigate whether conifers, beech and non-beech hardwood species were able to coexist on surfaces that differed in their underlying edaphic conditions, and if so to understand the mechanisms that influenced their regeneration on both poorly drained and well drained soils. The site selected was an area of high tree species diversity on a lowland 0.8 km² post-glacial terrace at the base of Mount Harata in the Grey River Valley. My approach was to use forest stand history reconstruction at two spatial scales: an intensive within-plot study of stand dynamics (chapter 1) and a whole-landform approach (chapter 2) that examined whether the dynamics identified at the smaller within-plot scale reflected larger patterns across the terrace. In chapter 1, three large permanent plots (0.3-0.7 ha) were placed at different points along the drainage gradient, one plot situated in each of the mainly well-drained, poorly drained and very poorly drained areas along the terrace. Information was gathered on species age and size structures, spatial distributions of tree ages, species interactions, microsite establishment preferences, patterns of stand mortality, and disturbance history in each plot. There were differences in stand structure, composition and relative abundance of species found between the well drained plot and the two poorer drained plots. On the well drained site conifers were scarce, the beeches Nothofagus fusca and N. menziesii dominated the canopy, and in the subcanopy the hardwood species Weinmannia racemosa and Quintinia acutifolia were abundant. As drainage became progressively poorer, the conifers Dacrydium cupressinum and Dacrycarpus dacrydioides became more abundant and occupied the emergent tier over a beech canopy. The hardwoods W. racemosa and Q. acutifolia became gradually less abundant in the subcanopy, whereas the hardwood Elaeocarpus hookerianus became more so. In the well drained plot, gap partitioning for light between beeches and hardwoods enabled coexistence in response to a range of different sized openings resulting from disturbances of different extent. In the two more poorly drained plots, species also coexisted by partitioning microsite establishment sites according to drainage. There were several distinct periods where synchronous establishment of different species occurred in different plots, suggesting there were large disturbances: c. 100yrs, 190-200 yrs, 275-300 yrs and 375-425 yrs ago. Generally after the same disturbance, different species regenerated in different plots reflecting the underlying drainage gradient. However, at the same site after different disturbances, different sets of species regenerated, suggesting the type and extent of disturbances and the conditions left behind influenced species regeneration at some times but not others. The regeneration of some species (e.g., N. fusca in the well-drained plot, and Dacrydium in the poorer drained plots) was periodic and appeared to be closely linked to these events. In the intervals between these disturbances, less extensive disturbances resulted in the more frequent N. menziesii and especially hardwood regeneration. The type of tree death caused by different disturbances favoured different species, with dead standing tree death favouring the more shade-tolerant N. menziesii and hardwoods, whereas uprooting created a mosaic of microsite conditions and larger gap sizes that enabled Dacrycarpus, N. fusca and E. hookerianus to maintain themselves in the poorly drained areas. In chapter 2, 10 circular plots (c. 0.12 ha) were placed in well drained areas and 10 circular plots (c. 0.2 ha) in poorly drained plots to collect information on species population structures and microsite preferences. The aims were to reconstruct species' regeneration responses to a range of disturbances of different type and extent across the whole terrace, and to examine whether there were important differences in the effects of these disturbances. At this landform scale, the composition and relative abundances of species across the drainage gradient reflected those found in chapter 1. There were few scattered conifers in well drained areas, despite many potential regeneration opportunities created from a range of different stand destroying and smaller scale disturbances. Three of the four periods identified in chapter 1 reflected distinct terrace-wide periods of regeneration 75-100 yrs, 200-275 yrs and 350-450 yrs ago, providing strong evidence of periodic large, infrequent disturbances that occurred at intervals of 100-200 yrs. These large, infrequent disturbances have had a substantial influence in determining forest history, and have had long term effects on forest structure and successional processes. Different large, infrequent disturbances had different effects across the terrace, with the variability in conditions that resulted enabling different species to regenerate at different times. For example, the regeneration of distinct even-aged Dacrydium cohorts in poorly drained areas was linked to historical Alpine Fault earthquakes, but not to more recent storms. The variation in the intensity of different large, infrequent disturbances at different points along the environmental drainage gradient, was a key factor influencing the scale of impacts. In effect, the underlying edaphic conditions influenced species composition along the drainage gradient and disturbance history regulated the relative abundances of species. The results presented here further emphasise the importance of large scale disturbances as a mechanism that allows coexistence of different tree species in mixed forest, in particular for the conifers Dacrydium, Dacrycarpus and the beech N. fusca, by creating much of the environmental variation to which these species responded. This study adds to our understanding of the effects of historical earthquakes in the relatively complex forests of north Westland, and further illustrates their importance in the Westland forest landscape as the major influential disturbance on forest pattern and history. These results also further develop the 'two-component' model used to describe conifer/angiosperm dynamics, by identifying qualitative differences in the impacts of different large, infrequent disturbances across an environmental gradient that allowed for coexistence of different species. In poorer drained areas, these forests may even be thought of as 'three-component' systems with conifers, beeches and hardwoods exhibiting key differences in their regeneration patterns after disturbances of different type and extent, and in their microsite preferences.