Current research in geotechnical engineering at the University of Canterbury includes a number of laboratory testing programmes focussed on understanding the behaviour of natural soil deposits in Christchurch during the 2010-2011 Canterbury Earthquake Sequence. Many soils found in Christchurch are sands or silty sands with little to no plasticity, making them very difficult to sample using established methods. The gel-push sampling methodology, developed by Kiso-Jiban Consultants in Japan, was developed to address some of the deficiencies of existing sampling techniques and has been deployed on two projects in Christchurch. Gel push sampling is carried out with a range of samplers which are modified versions of existing technology, and the University of Canterbury has acquired three versions of the tools (GP-S, GP-Tr, GP-D). Soil samples are extracted from the bottom of a freshly drilled borehole and are captured within a liner barrel, close to 1m in length. A lubricating polymer gel coats the outside of the soil sample as it enters the liner barrel. The frictional rubbing which normally occurs on the sides of the soil samples using existing techniques is eliminated by the presence of the polymer gel. The operation of the gel-push samplers is significantly more complicated than conventional push-tube samplers, and in the initial trials a number of operational difficulties were encountered, requiring changes to the sampling procedures. Despite these issues, a number of high quality soil samples were obtained on both projects using the GP-S sampler to capture silty soil. Attempts were made to obtain clean sands using a different gel-push sampler (GP-TR) in the Red Zone. The laboratory testing of these sands indicated that they were being significantly disturbed during the sampling and/or transportation procedures. While it remains too early to draw definitive conclusions regarding the performance of the gel-push samplers, the methodology has provided some promising results. Further trialling of the tools are required to refine operating procedures understand the full range of soil conditions which can be successfully sampled using the tools. In parallel with the gel-push trials, a Dames and Moore fixed-piston sampler has been used by our research partners from Berkeley to obtain soil samples at a number of sites within Christchurch. This sampler features relatively short (50cm), thin-walled liner barrels which is advanced into the ground under the action of hydraulic pressure. By reducing the overall length of the soil being captured, the disturbance to the soil as it enters the liner barrel is significantly reduced. The Dames and Moore sampler is significantly easier to operate than the gel-push sampler, and past experience has shown it to be successful in soft, plastic materials (i.e. clays and silty clays). The cyclic resistance of one silty clay obtained using both the gel-push and Dames & Moore samplers has been found to be very similar, and ongoing research aims to establish whether similar results are obtained for different soil types, including silty materials and clean sands.
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.
Abstract The original intention for the Partnership Community Worker (PCW) project in 2006 was for it to be an extension of the Pegasus Health General Practice and furthermore to be a bridge between the community and primary healthcare. It was believed that a close working relationship between the Practice Nurse and the PCW would help the target population of Māori, Pacifica and low income people to address and overcome their perceived barriers to healthcare which included: finance, transport, anxiety, cultural issues, communication, or lack of knowledge. Seven years later although the PCW project has been deemed a success in the Canterbury District Health Board annual reports (2013-14) and community and government agencies, including the Christchurch Resettlement Service (2012), many of the Pegasus Health General Practices have not utilised the project to its full extent, hence the need for this research. I was interested in finding out in the first instance if the model had changed and, if so why, and in the second instance if the promotional material currently distributed by Pegasus Health Primary Health Organisation reflected the daily practice of the PCW. A combination of methods were used including: surveys to the Pegasus Health General Practices, interviews with PCWs, interviews with managers of both the PCW host organisations and referring agencies to the PCW project. All the questions asked of all the participants in this research were focussed on their own perception of the role of the PCW. Results showed that the model has changed and although the publications were not reflecting the original intention of the project they did reflect the daily practice of the PCWs who are now struggling to meet much wider community expectations and needs. Key Results: Partnership Community Worker (PCW) interviews: Seventeen PCWs of the 19 employed were interviewed face to face. A number expressed interest in more culturally specific training and some are pursuing qualifications in social work; for many pay parity is an issue. In addition, many felt overwhelmed by the expectations around clients with mental health issues and housing issues now, post-earthquakes. Medical Practice surveys: Surveys were sent to eighty-two Pegasus Health medical practices and of these twenty five were completed. Results showed the full capacity of the PCW role was not clearly understood by all with many believing it was mostly a transport service. Those who did understand the full complexity of the role were very satisfied with the outcomes. PCW Host Community Manager Interviews: Of the ten out of twelve managers interviewed, some wished for more communication with Pegasus Health management because they felt aspects of both the PCW role and their own role as managers had become blurred over time. Referring organisations: Fifteen of the fifty referring community or government organisations participated. The overall satisfaction of the service was high and some acknowledged the continuing need for PCWs to be placed in communities where they were well known and trusted. Moreover results also showed that both the Canterbury earthquakes 2010-2011 and the amalgamation of Partnership Health PHO and Pegasus Health Charitable Limited in 2013 have contributed to the change of the model. Further future research may also be needed to examine the long term effects on the people of Canterbury involved in community work during the 2011-2014 years.
Recent earthquakes have highlighted the vulnerability of existing structure to seismic loading. Current seismic retrofit strategies generally focus on increasing the strength/stiffness in order to upgrade the seismic performance of a structure or element. A typical drawback of this approach is that the demand on the structural and sub-structural elements can be increased. This is of particular importance when considering the foundation capacity, which may already be insufficient to allow the full capacity of the existing wall to develop (due to early codes being gravity load orientated). In this thesis a counter intuitive but rational seismic retrofit strategy, termed "selective weakening" is introduced and investigated. This is the first stage of an ongoing research project underway at the University of Canterbury which is focusing on developing selective weakening techniques for the seismic retrofit of reinforced concrete structures. In this initial stage the focus is on developing selective weakening for the seismic retrofit of structural walls. This is performed using a series of experimental, analytical and numerical investigations. A procedure for the assessment of existing structural walls is also compiled, based on the suggestions of currently available code provisions. A selective weakening intervention is performed within an overall performance-based retrofit approach with the aim of improving the inelastic behaviour by first reducing the strength/stiffness of specific members within the structural system. This will be performed with the intention of modifying a shear type behaviour towards a flexural type behaviour. As a result the demand on the structural member will be reduced. Once weakening has been implemented the designer can use the wide range of techniques and materials available (e.g. use of FRP, jacketing or shotcrete) to ensure that adequate characteristics are achieved. Whilst performing this it has to be assured that the structure meets specific performance criteria and the principles of capacity design. A target of the retrofit technique is the ability to introduce the characteristics of recently developed high performance seismic resisting systems, consisting of a self centring and dissipative behaviour (commonly referred to as a hybrid system). In this thesis, results of experimental investigations performed on benchmark and selectively weakened walls are discussed. The investigations consisted of quasi-static cyclic uni-directional tests on two benchmark and two retrofitted cantilever walls. The first benchmark wall is detailed as typical of pre-1970's construction practice. An equivalent wall is retrofitted using a selective weakening approach involving a horizontal cut at foundation level to allow for a rocking response. The second benchmark wall represents a more severe scenario where the inelastic behaviour is dominated by shear. A retrofit solution involving vertically segmenting the wall to improve the ductility and retain gravity carrying capacity by inducing a flexural response is implemented. Numerical investigations on a multi-storey wall system are performed using non linear time history analysis on SDOF and MDOF lumped plasticity models, representing an as built and retrofitted prototype structure. Calibration of the hysteretic response to experimental results is carried out (accounting for pinching and strength degradation). The sensitivity of maximum and residual drifts to p-delta and strength degradation is monitored, along with the sensitivity of the peak base shear to higher mode affects. The results of the experimental and analytical investigations confirmed the feasibility and viability of the proposed retrofit technique, towards improving the seismic performance of structural walls.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
A non-destructive hardness testing method has been developed to investigate the amount of plastic strain demand in steel elements subjected to cyclic loading. The focus of this research is on application to the active links of eccentrically braced frames (EBFs), which are a commonly used seismic-resisting system in modern steel framed buildings. The 2010/2011 Christchurch earthquake series, especially the very intense February 22 shaking, which was the first earthquake worldwide to push complete EBF systems fully into their inelastic state, generating a moderate to high level of plastic strain in EBF active links, for a range of buildings from 3 to 23 storeys in height. This raised two important questions: 1) what was the extent of plastic deformation in active links; and 2) what effect does that have to post-earthquake steel properties? This project comprised determining a robust relationship between hardness and plastic strain in order to be able to answer the first question and provide the necessary input into answering the second question. A non-destructive Leeb (portable) hardness tester (model TH170) has been used to measure the hardness, in order to determine the plastic strain, in hot rolled steel universal sections and steel plates. A bench top Rockwell B was used to compare and validated the hardness measured by the portable hardness tester. Hardness was measured from monotonically strained tensile test specimens to identify the relationship between hardness and plastic strain demand. Test results confirmed a good relationship between hardness and the amount of monotonically induced plastic strain. Surface roughness was identified as an important parameter in obtaining reliable hardness readings from a portable hardness reader. A proper surface preparation method was established by using three different cleaning methods, finished with hand sanding to achieve surface roughness coefficients sufficiently low not to distort the results. This work showed that a test surface roughness (Ra) is not more than 1.6 micron meter (μm) is required for accurate readings from the TH170 tester. A case study on an earthquake affected building was carried out to identify the relationship between hardness and amount of plastic strain demand in cyclically deformed active links. Hardness was carried out from active links shown visually to have been the most affected during one of the major earthquake events. Onsite hardness test results were then compared with laboratory hardness test results. A good relationship between hardness from onsite and laboratory was observed between the test methods; Rockwell B bench top and portable Leeb tester TH170. Manufacturing induced plastic strain in the top and bottom of the webs of hot rolled sections were discovered from this research, an important result which explains why visual effects of earthquake induced active link yielding (eg cracked or flaking paint) was typically more prevalent over the middle half depth of the active link. The extent of this was quantified. It was also evident that the hardness readings from the portable hardness tester are influenced by geometry, mass effects and rigidity of the links. The final experimental stage was application of the method to full scale cyclic inelastic tested nominally identical active links subjected to loading regimes comprising constant and variable plastic strain demands. The links were cyclically loaded to achieve different plastic strain level. A novel Digital Image Correlation (DIC) technique was incorporated during the tests of this scale, to confirm the level of plastic strain achieved. Tensile test specimens were water jet cut from cyclically deformed webs to analyse the level of plastic strain. Test results show clear evidence that cyclically deformed structural steel elements show good correlation between hardness and the amount of plastic strain demand. DIC method was found to be reliable and accurate to check the level of plastic strain within cyclically deformed structural steel elements.
Saltwater Forest is a Dacrydium cupressinum-dominated lowland forest covering 9000 ha in south Westland, South Island, New Zealand. Four thousand hectares is managed for sustainable production of indigenous timber. The aim of this study was to provide an integrated analysis of soils, soil-landform relationships, and soil-vegetation relationships at broad and detailed scales. The broad scale understandings provide a framework in which existing or future studies can be placed and the detailed studies elucidate sources of soil and forest variability. Glacial landforms dominate. They include late Pleistocene lateral, terminal and ablation moraines, and outwash aggradation and degradation terraces. Deposits and landforms from six glacial advances have been recognised ranging from latest Last (Otira) Glaciation to Penultimate (Waimea) Glaciation. The absolute ages of landforms were established by analysis of the thickness and soil stratigraphy of loess coverbeds, augmented with radiocarbon dating and phytolith and pollen analysis. In the prevailing high rainfall of Westland soil formation is rapid. The rate of loess accretion in Saltwater Forest (ca. 30 mm ka⁻¹) has been low enough that soil formation and loess accretion took place contemporaneously. Soils formed in this manner are known as upbuilding soils. The significant difference between upbuilding pedogenesis and pedogenesis in a topdown sense into an existing sediment body is that each subsoil increment of an upbuilding soil has experienced processes of all horizons above. In Saltwater Forest subsoils of upbuilding soils are strongly altered because they have experienced the extremely acid environment of the soil surface at some earlier time. Some soil chronosequence studies in Westland have included upbuilding soils formed in loess as the older members of the sequence. Rates and types of processes inferred from these soils should be reviewed because upbuilding is a different pedogenic pathway to topdown pedogenesis. Landform age and morphology were used as a primary stratification for a study of the soil pattern and nature of soil variability in the 4000 ha production area of Saltwater Forest. The age of landforms (> 14 ka) and rapid soil formation mean that soils are uniformly strongly weathered and leached. Soils include Humic Organic Soils, Perch-gley Podzols, Acid Gley Soils, Allophanic Brown Soils, and Orthic or Pan Podzols. The major influence on the nature of soils is site hydrology which is determined by macroscale features of landforms (slope, relief, drainage density), mesoscale effects related to position on landforms, and microscale influences determined by microtopography and individual tree effects. Much of the soil variability arises at microscales so that it is not possible to map areas of uniform soils at practical map scales. The distribution of soil variability across spatial scales, in relation to the intensity of forest management, dictates that it is most appropriate to map soil complexes with boundaries coinciding with landforms. Disturbance of canopy trees is an important agent in forest dynamics. The frequency of forest disturbance in the production area of Saltwater Forest varies in a systematic way among landforms in accord with changes in abundance of different soils. The frequency of forest turnover is highest on landforms with the greatest abundance of extremely poorly-drained Organic Soils. As the abundance of better-drained soils increases the frequency of forest turnover declines. Changes in turnover frequency are reflected in the mean size and density of canopy trees (Dacrydium cupressinum) among landforms. Terrace and ablation moraine landforms with the greatest abundance of extremely poorly-drained soils have on average the smallest trees growing most densely. The steep lateral moraines, characterised by well drained soils, have fewer, larger trees. The changes manifested at the landform scale are an integration of processes operating over much shorter range as a result of short-range soil variability. The systematic changes in forest structure and turnover frequency among landforms and soils have important implications for sustainable forest management.
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.
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.
The recent Canterbury earthquake sequence in 2010-2011 highlighted a uniquely severe level of structural damage to modern buildings, while confirming the high vulnerability and life threatening of unreinforced masonry and inadequately detailed reinforced concrete buildings. Although the level of damage of most buildings met the expected life-safety and collapse prevention criteria, the structural damage to those building was beyond economic repair. The difficulty in the post-event assessment of a concrete or steel structure and the uneconomical repairing costs are the big drivers of the adoption of low damage design. Among several low-damage technologies, post-tensioned rocking systems were developed in the 1990s with applications to precast concrete members and later extended to structural steel members. More recently the technology was extended to timber buildings (Pres-Lam system). This doctoral dissertation focuses on the experimental investigation and analytical and numerical prediction of the lateral load response of dissipative post-tensioned rocking timber wall systems. The first experimental stages of this research consisted of component testing on both external replaceable devices and internal bars. The component testing was aimed to further investigate the response of these devices and to provide significant design parameters. Post-tensioned wall subassembly testing was then carried out. Firstly, quasi-static cyclic testing of two-thirds scale post-tensioned single wall specimens with several reinforcement layouts was carried out. Then, an alternative wall configuration to limit displacement incompatibilities in the diaphragm was developed and tested. The system consisted of a Column-Wall-Column configuration, where the boundary columns can provide the support to the diaphragm with minimal uplifting and also provide dissipation through the coupling to the post-tensioned wall panel with dissipation devices. Both single wall and column-wall-column specimens were subjected to drifts up to 2% showing excellent performance, limiting the damage to the dissipating devices. One of the objectives of the experimental program was to assess the influence of construction detailing, and the dissipater connection in particular proved to have a significant influence on the wall’s response. The experimental programs on dissipaters and wall subassemblies provided exhaustive data for the validation and refinement of current analytical and numerical models. The current moment-rotation iterative procedure was refined accounting for detailed response parameters identified in the initial experimental stage. The refined analytical model proved capable of fitting the experimental result with good accuracy. A further stage in this research was the validation and refinement of numerical modelling approaches, which consisted in rotational spring and multi-spring models. Both the modelling approaches were calibrated versus the experimental results on post-tensioned walls subassemblies. In particular, the multi-spring model was further refined and implemented in OpenSEES to account for the full range of behavioural aspects of the systems. The multi-spring model was used in the final part of the dissertation to validate and refine current lateral force design procedures. Firstly, seismic performance factors in accordance to a Force-Based Design procedure were developed in accordance to the FEMA P-695 procedure through extensive numerical analyses. This procedure aims to determine the seismic reduction factor and over-strength factor accounting for the collapse probability of the building. The outcomes of this numerical analysis were also extended to other significant design codes. Alternatively, Displacement-Based Design can be used for the determination of the lateral load demand on a post-tensioned multi-storey timber building. The current DBD procedure was used for the development of a further numerical analysis which aimed to validate the procedure and identify the necessary refinements. It was concluded that the analytical and numerical models developed throughout this dissertation provided comprehensive and accurate tools for the determination of the lateral load response of post-tensioned wall systems, also allowing the provision of design parameters in accordance to the current standards and lateral force design procedures.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
As cities evolve, change and grow, the need and desire for adaptable architecture becomes evident across the nation. Architecture needs to undertake techniques that are flexible in order to adapt and align with the development of future generations in New Zealand. The Education industry is a primary example of a sector which requires flexibility within both classroom architectural form and interior configuration. This is a resultant of the recently updated Ministry of Education requirements; which state that every new classroom built or renovated nationwide, must implement the MoE classroom design standards for Innovative Learning Environments. ILE teaching spaces are configured as an open plan interior, supporting flexibility in classroom arrangement and teaching techniques. ILE classrooms are capable of evolving and adapting as educational practices evolve and change, allowing schools to remain modern and future focused. As part of this movement to ILE, the Ministry of Education has also recently made an attempt to improve the quality of temporary classrooms. This has been done by looking into the initiation of a programme that utilizes relocatable classroom buildings. Relocatable classrooms have been selected for multiple reasons, primarily flexibility. Flexibility is key for a school environment as it allows the school to actively respond to fluctuating school rolls. It is anticipated that the programme will provide a faster delivery process with a standardised design that allows the classrooms to be relocated from one school to another with relative ease. Following the devastating February 2011 earthquake the Greater Christchurch Region, the Education sector is in the midst of the Canterbury Schools Rebuild Programme. As a repercussion of this natural disaster, the majority of Christchurch schools have redevelopment or rebuild projects in progress, with preliminary design phases already in action for a small group of select schools regarded as high priority. The primary funding for these projects are sourced from insurance money, implementing tight budget restrictions, affecting the architectural design, quality and speed of the construction and repair works. The available funding limits the affordable classroom options to basic teaching spaces that have been stripped back to simple architectural forms, dictating not only the re-design, but also how our future generations will learn. Thus causing the development of the new student-led learning ILE concept to become controlled by existing construction techniques and the Rebuild Programmes budget restrictions. This thesis focuses on the future proofing of New Zealand schools by providing an affordable and time efficient alternative option to the current static, traditional construction, an option that has the ability to cater to the unpredictable fluctuating school rolls across the nation. This has been done by developing a prefabricated system for standalone classroom blocks. These blocks have the ability to be relocated between different school sites, dynamically catering to the unpredictable school roll numbers experienced across New Zealand. This site flexibility is reflected with the interior flexibility in the classrooms, enhancing the internal teaching space composition and challenges the existing design standards set by the Ministry of Education for Innovative Learning Environments. This system is called “Flexi-Ed”. Flexibility has been a key driver for this thesis, as the prefabricated structure is have to be flexible in three ways; first in the sense of being easy to assemble and disassemble. Second by offering flexible interior learning environments and thirdly the joints of the structure are designed with the ability to be flexible in order to cope with seismic activity. These three principles will provide schools with long term flexibility, minimal on-site interruption and heighten the standard of ILE across the nation. I strive to provide schools with long term flexibility and minimal site interruption, whilst heightening the standard of Innovative Learning Environments across New Zealand.
Rock mass defect controlled deep-seated landslides are widespread within the deeply incised landscapes formed in Tertiary soft rock terrain in New Zealand. The basal failure surfaces of deep-seated slope failures are defined by thin, comparatively weak and laterally continuous bedding parallel layers termed critical stratigraphic horizons. These horizons have a sedimentary origin and have typically experienced some prior tectonically induced shear displacement at the time of slope failure. The key controls on the occurrence and form of deep-seated landslides are considered in terms of rock mass defect properties and tectonic and climatic forcing. The selection of two representative catchments (in southern Hawke's Bay and North Canterbury) affected by tectonic and climatic forcing has shown that the spatial and temporal initiation of deep-seated bedrock landslides in New Zealand Tertiary soft rock terrain is a predictable rather than a stochastic process; and that deep-seated landslides as a mass wasting process have a controlling role in landscape evolution in many catchments formed in Tertiary soft rock terrain. The Ella Landslide in North Canterbury is a deep-seated (~85 m) translational block slide that has failed on a 5 - 10 mm thick, kaolinite-rich, pre-sheared critical stratigraphic horizon. The residual strength of this sedimentary horizon, (C'R 2.6 - 2.7 kPa, and Ѳ'R = 16 - 21°), compared to the peak strength of the dominant lithology (C' = 176 kPa, and Ѳ' = 37°) defines a high strength contrast in the succession, and therefore a critical location for the basal failure surface of deep-seated slope failures. The (early to mid Holocene) Ella Landslide debris formed a large landslide dam in the Kate Stream catchment and this has significantly retarded rates of mass wasting in the middle catchment. Numerical stability analysis shows that this slope failure would have most likely required the influence of earthquake induced strong ground motion and the event is tentatively correlated to a Holocene event on the Omihi Fault. The influence of this slope failure is likely to affect the geomorphic development of the catchment on a scale of 10⁴ - 10⁵ years. In deeply incised catchments at the southeastern margin of the Maraetotara Plateau, southern Hawke's Bay, numerous widespread deep-seated landslides have basal failure surfaces defined by critical stratigraphic horizons in the form of thin « 20 mm) tuffaceous beds in the Makara Formation flysch (alternating sandstone and mudstone units). The geometry of deep-seated slope failures is controlled by these regularly spaced (~70 m), very weak critical stratigraphic horizons (C'R 3.8 - 14.2 kPa, and Ѳ'R = 2 - 5°), and regularly spaced (~45 m) and steeply dipping (-50°) critical conjugate joint/fault sets, which act as slide block release surfaces. Numerical stability analysis and historical precedent show that the temporal initiation of deep-seated landslides is directly controlled by short term tectonic forcing in the form of periodic large magnitude earthquakes. Published seismic hazard data shows the recurrence interval of earthquakes producing strong ground motions of 0.35g at the study site is every 150 yrs, however, if subduction thrust events are considered the level of strong ground motion may be much higher. Multiple occurrences of deep-seated slope failure are correlated to failure on the same critical stratigraphic horizon, in some cases in three adjacent catchments. Failure on multiple critical stratigraphic horizons leads to the development of a "stepped" landscape morphology. This slope form will be maintained during successive accelerated stream incision events (controlled by long term tectonic and climatic forcing) for as long as catchments are developing in this specific succession. Rock mass defect controlled deep seated landslides are controlling catchment head progression, landscape evolution and hillslope morphology in the Hawke's Bay study area and this has significant implications for the development of numerical landscape evolution models of landscapes formed in similar strata. Whereas the only known numerical model to consider deep seated landslides as an erosion process (ZSCAPE) considers them as stochastic in time and space, this study shows that this could not be applied to a landscape where the widespread spatial occurrence of deep-seated landslides is controlled by rock mass defects. In both of the study areas for this project, and by implication in many catchments in Tertiary soft rock terrain, deep-seated landslides controlled by rock mass defect strength, spacing and orientation, and tectonic and climatic forcing have an underlying control on landscape evolution. This study quantifies parameters for the development of numerical landscape evolution models that would assess the role of specific parameters, such as uplift rates, incision rates and earthquake recurrence in catchment evolution in Tertiary soft rock terrain.
The Canterbury earthquake sequence (2010-2011) was the most devastating catastrophe in New Zealand‘s modern history. Fortunately, in 2011 New Zealand had a high insurance penetration ratio, with more than 95% of residences being insured for these earthquakes. This dissertation sheds light on the functions of disaster insurance schemes and their role in economic recovery post-earthquakes. The first chapter describes the demand and supply for earthquake insurance and provides insights about different public-private partnership earthquake insurance schemes around the world. In the second chapter, we concentrate on three public earthquake insurance schemes in California, Japan, and New Zealand. The chapter examines what would have been the outcome had the system of insurance in Christchurch been different in the aftermath of the Canterbury earthquake sequence (CES). We focus on the California Earthquake Authority insurance program, and the Japanese Earthquake Reinsurance scheme. Overall, the aggregate cost of the earthquake to the New Zealand public insurer (the Earthquake Commission) was USD 6.2 billion. If a similar-sized disaster event had occurred in Japan and California, homeowners would have received only around USD 1.6 billion and USD 0.7 billion from the Japanese and Californian schemes, respectively. We further describe the spatial and distributive aspects of these scenarios and discuss some of the policy questions that emerge from this comparison. The third chapter measures the longer-term effect of the CES on the local economy, using night-time light intensity measured from space, and focus on the role of insurance payments for damaged residential property during the local recovery process. Uniquely for this event, more than 95% of residential housing units were covered by insurance and almost all incurred some damage. However, insurance payments were staggered over 5 years, enabling us to identify their local impact. We find that night-time luminosity can capture the process of recovery; and that insurance payments contributed significantly to the process of local economic recovery after the earthquake. Yet, delayed payments were less affective in assisting recovery and cash settlement of claims were more effective than insurance-managed repairs. After the Christchurch earthquakes, the government declared about 8000 houses as Red Zoned, prohibiting further developments in these properties, and offering the owners to buy them out. The government provided two options for owners: the first was full payment for both land and dwelling at the 2007 property evaluation, the second was payment for land, and the rest to be paid by the owner‘s insurance. Most people chose the second option. Using data from LINZ combined with data from Stats NZ, the fourth chapter empirically investigates what led people to choose this second option, and how peer effect influenced the homeowners‘ choices. Due to climate change, public disclosure of coastal hazard information through maps and property reports have been used more frequently by local government. This is expected to raise awareness about disaster risks in local community and help potential property owners to make informed locational decision. However, media outlets and business sector argue that public hazard disclosure will cause a negative effect on property value. Despite this opposition, some district councils in New Zealand have attempted to implement improved disclosure. Kapiti Coast district in the Wellington region serves as a case study for this research. In the fifth chapter, we utilize the residential property sale data and coastal hazard maps from the local district council. This study employs a difference-in-difference hedonic property price approach to examine the effect of hazard disclosure on coastal property values. We also apply spatial hedonic regression methods, controlling for coastal amenities, as our robustness check. Our findings suggest that hazard designation has a statistically and economically insignificant impact on property values. Overall, the risk perception about coastal hazards should be more emphasized in communities.
This thesis describes the management process of innovation through construction infrastructure projects. This research focuses on the innovation management process at the project level from four views. These are categorised into the separate yet related areas of: “innovation definition”, “Project time”, “project team motivation” and “Project temporary organisation”. A practical knowledge is developed for each of these research areas that enables project practitioners to make the best decision for the right type of innovation at the right phase of projects, through a capable project organisation. The research developed a holistic view on both innovation and the construction infrastructure project as two complex phenomena. An infrastructure project is a long-term capital investment, highly risky and an uncertain. Infrastructure projects can play a key role in innovation and performance improvement throughout the construction industry. The delivery of an infrastructure project is affected in most cases by critical issues of budget constraint, programme delays and safety Where the business climate is characterized by uncertainty, risk and a high level of technological change, construction infrastructure projects are unable to cope with the requirement to develop innovation. Innovation in infrastructure projects, as one of the key performance indicators (KPI) has been identified as a critical capability for performance improvement through the industry. However, in spite of the importance of infrastructure projects in improving innovation, there are a few research efforts that have developed a comprehensive view on the project context and its drivers and inhibitors for innovation in the construction industry. Two main reasons are given as the inhibitors through the process of comprehensive research on innovation management in construction. The first reason is the absence of an understanding of innovation itself. The second is a bias towards research at a firm and individual level, so a comprehensive assessment of project-related factors and their effects on innovation in infrastructure projects has not been undertaken. This study overcomes these issues by adopting as a case study approach of a successful infrastructure project. This research examines more than 500 construction innovations generated by a unique infrastructure alliance. SCIRT (Stronger Christchurch Infrastructure Rebuild Team) is a temporary alliancing organisation that was created to rebuild and recover the damaged infrastructure after the Christchurch 2011 earthquake. Researchers were given full access to the innovation project information and innovation systems under a contract with SCIRT Learning Legacy, provided the research with material which is critical for understanding innovations in large, complex alliancing infrastructure organisation. In this research, an innovation classification model was first constructed. Clear definitions have been developed for six types of construction innovation with a variety of level of novelties and benefits. The innovation classification model was applied on the SCIRT innovation database and the resultant trends and behaviours of different types of innovation are presented. The trends and behaviours through different types of SCIRT innovations developed a unique opportunity to research the projectrelated factors and their effect on the behaviour of different classified types of innovation throughout the project’s lifecycle. The result was the identification of specific characteristics of an infrastructure project that affect the innovation management process at the project level. These were categorised in four separate chapters. The first study presents the relationship between six classified types of innovation, the level of novelty and the benefit they come up with, by applying the innovation classification model on SCIRT innovation database. The second study focused on the innovation potential and limitations in different project lifecycle phases by using a logic relationship between the six classified types of innovation and the three classified phases of the SCIRT project. The third study result develops a holistic view of different elements of the SCIRT motivation system and results in a relationship between the maturity level of definition developed for innovation as one of the KPIs and a desire though the SCIRT innovation incentive system to motivate more important innovations throughout the project. The fourth study is about the role of the project’s temporary organisation that finally results in a multiple-view innovation model being developed for project organisation capability assessment in the construction industry. The result of this thesis provides practical and instrumental knowledge to be used by a project practitioner. Benefits of the current thesis could be categorized in four groups. The first group is the innovation classification model that provides a clear definition for six classified types of innovation with four levels of novelty and specifically defined outcomes and the relationship between the innovation types, novelty and benefit. The second is the ability that is provided for the project practitioner to make the best decision for the right type of innovation at the right phases of a project’s lifecycle. The third is an optimisation that is applied on the SCIRT innovation motivation system that enables the project practitioner to incentivize the right type of innovation with the right level of financial gain. This drives the project teams to develop a more important innovation instead of a simple problemsolving one. Finally, the last and probably more important benefit is the recommended multiple-view innovation model. This is a tool that could be used by a project practitioner in order to empower the project team to support innovation throughout the project.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
This project was initiated by ENGEO Limited and KiwiRail Holdings Limited to assess the stability of Slovens Creek Viaduct (specifically its western abutment) and a 3km section of rail corridor between Slovens Creek Viaduct and Avoca on the Midland Line (MDL). Commonly known as the scenic TranzAlpine rail journey (through Arthurs Pass National Park) the MDL connects Greymouth to Christchurch via Rolleston, where the MDL meets the Main South Line into Christchurch. The project area is approximately 40km southeast of Arthurs Pass Township, in the eastern extension of the Castle Hill Basin which is part of the Waimakariri Catchment and Canterbury Foothills. The field area is underlain by Rakaia Terrane, which is part of the Torlesse Composite Terrane forming the basement rock unit for the field area. Cretaceous-Tertiary rocks of the Castle Hill Basin overlie the basement strata and record a transgression-regression sequence, as well as mid-Oligocene submarine volcanism. The stratigraphic sequence in the Castle Hill Basin, and its eastern extension to Avoca, comprises two formations of the Eyre group, the older Broken River Formation and the younger Iron Creek Formation. Deep marine Porter Group limestones, marls, and tuffs of Oligocene age succeed the Iron Creek Formation of the Eyre Group, and probably records the maximum of the transgression. The Enys Formation lies disconformably on the Porter Group and is overlain unconformably by Late Pleistocene glacifluvial and glacial deposits. The Tertiary strata in the Slovens-Avoca rail corridor are weak, and the clay-rich tuff derived from mid-Oligocene volcanism is particularly prone to slaking. Extensive mapping carried out for this project has identified that some 90 percent of the surface along the length of the Slovens-Avoca corridor has been subject to mass movement. The landslides of the Slovens-Avoca rail corridor are clearly younger than the Last Glaciation, and Slovens Creek has been downcutting, with associated faulting and uplift, to form the present day geomorphology of the rail corridor. Deep-seated landslides in the rail corridor extend to Slovens Creek, locally deflecting the stream course, and a generic ground failure model for the rail corridor has been developed. Exploratory geotechnical investigations, including core drilling, installation of an inclinometer and a piezometer, enabled the construction of a simple ground model and cross section for the Slovens Creek Viaduct western abutment. Limit-equilibrium and pseudo-static slope stability analyses using both circular and block critical slip surface search methods were applied to the ground model for the western abutment of Slovens Creek Viaduct. Piezometric and strength data obtained during laboratory testing of core material have been used to constrain the western abutment stability assessment for one representative section line (C-C’). Prior to pseudo-static sensitivity analyses peak ground acceleration (PGA) for various Ultimate Limit State (ULS) design return periods, defined by an equation given in NZS1170.5:2004, were calculated and have been used as a calibration technique to find and compare specific PGA values for pseudo-static analyses in the Slovens Creek Viaduct area. The main purpose has been to provide an indication of how railway infrastructure could be affected by seismic events of various return periods defined by ULS design standards for the area. Limit equilibrium circular slip surface search methods, both grid search and auto refine search, indicated the slope is stable with a FoS greater than 1.0 returned from each, although one particular surface returned the lowest FoS in each. This surface is in the lower portion of the slope, adjacent to Slovens Stream and northeast of the MDL. As expected, pseudo-static analyses returned a lower FoS overall when compared to limit equilibrium analyses. The PGA analyses suggest that partial ground failure at the Slovens Creek Viaduct western abutment could occur in a 1 in 25-year return period event within materials on the slower slope beyond the immediate rail corridor. A ULS (1 in 500-year) event in the Slovens Creek Viaduct area would likely produce a PGA of ~0.9g, and the effects on the western abutment and rail infrastructure would most likely be catastrophic. Observed ground conditions for the western abutment of the Slovens Creek Viaduct suggest there is no movement within the landslide at depth within the monitoring timeframe of this project (22 May 2015 – 4 August 2015). Slope stability monitoring is recommended to be continued in two parts: (1) the inclinometer in BH1 is to be monitored on a six monthly basis for one year following completion of this thesis, and then annually unless ground movements become evident; and (2) surface movement monitoring should be installed using a fixed datum on the stable eastern abutment. Long-term stability management strategies for the Slovens Creek Viaduct western abutment are dependent upon future observed changes and ongoing monitoring. Hazard and risk assessment using the KiwiRail Qualitative Risk Assessment Framework (QRA) is recommended, and if slope stability becomes problematic for operation of the Midland Line consideration should be given to deep slope drainage. In the event of a large magnitude or high PGA earthquake all monitoring should be reviewed.
<strong>Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data.</strong>In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā.This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.
Background: We are in a period of history where natural disasters are increasing in both frequency and severity. They are having widespread impacts on communities, especially on vulnerable communities, those most affected who have the least ability to prepare or respond to a disaster. The ability to assemble and effectively manage Interagency Emergency Response Teams (IERTs) is critical to navigating the complexity and chaos found immediately following disasters. These teams play a crucial role in the multi-sectoral, multi-agency, multi-disciplinary, and inter-organisational response and are vital to ensuring the safety and well-being of vulnerable populations such as the young, aged, and socially and medically disadvantaged in disasters. Communication is key to the smooth operation of these teams. Most studies of the communication in IERTs during a disaster have been focussed at a macro-level of examining larger scale patterns and trends within organisations. Rarely found are micro-level analyses of interpersonal communication at the critical interfaces between collaborating agencies. This study set out to understand the experiences of those working at the interagency interfaces in an IERT set up by the Canterbury District Health Board to respond to the needs of the vulnerable people in the aftermath of the destructive earthquakes that hit Canterbury, New Zealand, in 2010-11. The aim of the study was to gain insights about the complexities of interpersonal communication (micro-level) involved in interagency response coordination and to generate an improved understanding into what stabilises the interagency communication interfaces between those agencies responding to a major disaster. Methods: A qualitative case study research design was employed to investigate how interagency communication interfaces were stabilised at the micro-level (“the case”) in the aftermath of the destructive earthquakes that hit Canterbury in 2010-11 (“the context”). Participant recruitment was undertaken by mapping which agencies were involved within the IERT and approaching representatives from each of these agencies. Data was collected via individual interviews using a semi-structured interview guide and was based on the “Critical Incident Technique”. Subsequently, data was transcribed verbatim and subjected to inductive analysis. This was underpinned theoretically by Weick’s “Interpretive Approach” and supported by Nvivo qualitative data analysis software. Results: 19 participants were interviewed in this study. Out of the inductive analysis emerged two primary themes, each with several sub-factors. The first major theme was destabilising/disruptive factors of interagency communication with five sub-factors, a) conflicting role mandates, b) rigid command structures, c) disruption of established communication structures, d) lack of shared language and understanding, and e) situational awareness disruption. The second major theme stabilising/steadying factors in interagency communication had four sub-factors, a) the establishment of the IERT, b) emergent novel communication strategies, c) establishment of a liaison role and d) pre-existing networks and relationships. Finally, there was a third sub-level identified during inductive analysis, where sub-factors from both primary themes were noted to be uniquely interconnected by emergent “consequences” arising out of the disaster context. Finally, findings were synthesised into a conceptual “Model of Interagency Communication at the Micro-level” based on this case study of the Canterbury earthquake disaster response. Discussion: The three key dimensions of The People, The Connections and The Improvisations served as a framework for the discussion of what stabilises interagency communication interfaces in a major disaster. The People were key to stabilising the interagency interfaces through functioning as a flexible conduit, guiding and navigating communication at the interagency interfaces and improving situational awareness. The Connections provided the collective competence, shared decision-making and prior established relationships that stabilised the micro-level communication at interagency interfaces. And finally, The Improvisations i.e., novel ideas and inventiveness that emerge out of rapidly changing post-disaster environments, also contributed to stabilisation of micro-level communication flows across interagency interfaces in the disaster response. “Command and control” hierarchical structures do provide clear processes and structures for teams working in disasters to follow. However, improvisations and novel solutions are also needed and often emerge from first responders (who are best placed to assess the evolving needs in a disaster where there is a high degree of uncertainty). Conclusion: This study highlights the value of incorporating an interface perspective into any study that seeks to understand the processes of IERTs during disaster responses. It also strengthens the requirement for disaster management frameworks to formally plan for and to allow for the adaptive responsiveness of local teams on the ground, and legitimise and recognise the improvisations of those in the role of emergent boundary spanners in a disaster response. This needs to be in addition to existing formal disaster response mechanisms. This study provides a new conceptual model that can be used to guide future case studies exploring stability at the interfaces of other IERTs and highlights the centrality of communication in the experiences of members of teams in the aftermath of a disaster. Utilising these new perspectives on stabilising communication at the interagency interfaces in disaster responses will have practical implications in the future to better serve the needs of vulnerable people who are at greatest risk of adverse outcomes in a disaster.
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.
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.