Search

found 480 results

Research papers, University of Canterbury Library

A buckling-restrained braced frame (BRBF) is a structural bracing system that provides lateral strength and stiffness to buildings and bridges. They were first developed in Japan in the 1970s (Watanabe et al. 1973, Kimura et al. 1976) and gained rapid acceptance in the United States after the Northridge earthquake in 1994 (Bruneau et al. 2011). However, it was not until the Canterbury earthquakes of 2010/2011, that the New Zealand construction market saw a significant uptake in the use of buckling-restrained braces (BRBs) in commercial buildings (MacRae et al. 2015). In New Zealand there is not yet any documented guidance or specific instructions in regulatory standards for the design of BRBFs. This makes it difficult for engineers to anticipate all the possible stability and strength issues within a BRBF system and actively mitigate them in each design. To help ensure BRBF designs perform as intended, a peer review with physical testing are needed to gain building compliance in New Zealand. Physical testing should check the manufacturing and design of each BRB (prequalification testing), and the global strength and stability of each BRB its frame (subassemblage testing). However, the financial pressures inherent in commercial projects has led to prequalification testing (BRB only testing) being favoured without adequate design specific subassemblage testing. This means peer reviewers have to rely on BRB suppliers for assurances. This low regulation environment allows for a variety of BRBF designs to be constructed without being tested or well understood. The concern is that there may be designs that pose risk and that issues are being overlooked in design and review. To improve the safety and design of BRBFs in New Zealand, this dissertation studies the behaviour of BRBs and how they interact with other frame components. Presented is the experimental test process and results of five commercially available BRB designs (Chapter 2). It discusses the manufacturing process, testing conditions and limitations of observable information. It also emphasises that even though subassemblage testing is impractical, uniaxial testing of the BRB only is not enough, as this does not check global strength or stability. As an alternative to physical testing, this research uses computer simulation to model BRB behaviour. To overcome the traditional challenges of detailed BRB modelling, a strategy to simulate the performance of generic BRB designs was developed (Chapter 3). The development of nonlinear material and contact models are important aspects of this strategy. The Chaboche method is employed using a minimum of six backstress curves to characterize the combined isotropic and kinematic hardening exhibited by the steel core. A simplified approach, adequate for modelling the contact interaction between the restrainer and the core was found. Models also capture important frictional dissipation as well as lateral motion and bending associated with high order constrained buckling of the core. The experimental data from Chapter 2 was used to validate this strategy. As BRBs resist high compressive loading, global stability of the BRB and gusseted connection zone need to be considered. A separate study was conducted that investigated the yielding and buckling strength of gusset plates (Chapter 4). The stress distribution through a gusset plate is complex and difficult to predict because the cross-sectional area of gusset plate is not uniform, and each gusset plate design is unique in shape and size. This has motivated design methods that approximate yielding of gusset plates. Finite element modelling was used to study the development of yielding, buckling and plastic collapse behaviour of a brace end bolted to a series of corner gusset plates. In total 184 variations of gusset plate geometries were modelled in Abaqus®. The FEA modelling applied monotonic uniaxial load with an imperfection. Upon comparing results to current gusset plate design methods, it was found that the Whitmore width method for calculating the yield load of a gusset is generally un-conservative. To improve accuracy and safety in the design of gusset plates, modifications to current design methods for calculating the yield area and compressive strength for gusset plates is proposed. Bolted connections are a popular and common connection type used in BRBF design. Global out-of-plane stability tends to govern the design for this connection type with numerous studies highlighting the risk of instability initiated by inelasticity in the gussets, neck of the BRB end and/or restrainer ends. Subassemblage testing is the traditional method for evaluating global stability. However, physical testing of every BRBF variation is cost prohibitive. As such, Japan has developed an analytical approach to evaluate out-of-plane stability of BRBFs and incorporated this in their design codes. This analytical approach evaluates the different BRB components under possible collapse mechanisms by focusing on moment transfer between the restrainer and end of the BRB. The approach have led to strict criteria for BRBF design in Japan. Structural building design codes in New Zealand, Europe and the United States do not yet provide analytical methods to assess BRB and connection stability, with prototype/subassemblage testing still required as the primary means of accreditation. Therefore it is of interest to investigate the capability of this method to evaluate stability of BRBs designs and gusset plate designs used in New Zealand (including unstiffened gusset connection zones). Chapter 5 demonstrates the capability of FEA to study to the performance of a subassemblage test under cyclic loading – resembling that of a diagonal ground storey BRBF with bolted connections. A series of detailed models were developed using the strategy presented in Chapter 3. The geometric features of BRB 6.5a (Chapter 2) were used as a basis for the BRBs modelled. To capture the different failure mechanisms identified in Takeuchi et al. (2017), models varied the length that the cruciform (non-yielding) section inserts into the restrainer. Results indicate that gusset plates designed according to New Zealand’s Steel Structures Standard (NZS 3404) limit BRBF performance. Increasing the thickness of the gusset plates according to modifications discussed in Chapter 4, improved the overall performance for all variants (except when Lin/ Bcruc = 0.5). The effect of bi-directional loading was not found to notably affect out-of-plane stability. Results were compared against predictions made by the analytical method used in Japan (Takeuchi method). This method was found to be generally conservative is predicting out-of-plane stability of each BRBF model. Recommendations to improve the accuracy of Takeuchi’s method are also provided. The outcomes from this thesis should be helpful for BRB manufacturers, researchers, and in the development of further design guidance of BRBFs.

Research papers, The University of Auckland Library

High demolition rates were observed in New Zealand after the 2010-2011 Canterbury Earthquake Sequence despite the success of modern seismic design standards to achieve required performance objectives such as life safety and collapse prevention. Approximately 60% of the multi-storey reinforced concrete (RC) buildings in the Christchurch Central Business District were demolished after these earthquakes, even when only minor structural damage was present. Several factors influenced the decision of demolition instead of repair, one of them being the uncertainty of the seismic capacity of a damaged structure. To provide more insight into this topic, the investigation conducted in this thesis evaluated the residual capacity of moderately damaged RC walls and the effectiveness of repair techniques to restore the seismic performance of heavily damaged RC walls. The research outcome provided insights for developing guidelines for post-earthquake assessment of earthquake-damaged RC structures. The methodology used to conduct the investigation was through an experimental program divided into two phases. During the first phase, two walls were subjected to different types of pre-cyclic loading to represent the damaged condition from a prior earthquake, and a third wall represented a repair scenario with the damaged wall being repaired using epoxy injection and repair mortar after the pre-cyclic loading. Comparisons of these test walls to a control undamaged wall identified significant reductions in the stiffness of the damaged walls and a partial recovery in the wall stiffness achieved following epoxy injection. Visual damage that included distributed horizontal and diagonal cracks and spalling of the cover concrete did not affect the residual strength or displacement capacity of the walls. However, evidence of buckling of the longitudinal reinforcement during the pre-cyclic loading resulted in a slight reduction in strength recovery and a significant reduction in the displacement capacity of the damaged walls. Additional experimental programs from the literature were used to provide recommendations for modelling the response of moderately damaged RC walls and to identify a threshold that represented a potential reduction in the residual strength and displacement capacity of damaged RC walls in future earthquakes. The second phase of the experimental program conducted in this thesis addressed the replacement of concrete and reinforcing steel as repair techniques for heavily damaged RC walls. Two walls were repaired by replacing the damaged concrete and using welded connections to connect new reinforcing bars with existing bars. Different locations of the welded connections were investigated in the repaired walls to study the impact of these discontinuities at the critical section. No significant changes were observed in the stiffness, strength, and displacement capacity of the repaired walls compared to the benchmark undamaged wall. Differences in the local behaviour at the critical section were observed in one of the walls but did not impact the global response. The results of these two repaired walls were combined with other experimental programs found in the literature to assemble a database of repaired RC walls. Qualitative and quantitative analyses identified trends across various parameters, including wall types, damage before repair, and repair techniques implemented. The primary outcome of the database analysis was recommendations for concrete and reinforcing steel replacement to restore the strength and displacement capacity of heavily damaged RC walls.

Research Papers, Lincoln University

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

Research Papers, Lincoln University

Knowledge of past climate variability is essential for understanding present and future climate trends. This study used Halocarpus biformis (pink pine) ring-width chronologies to investigate palaeotemperature history in Westland, New Zealand. The ensuing reconstruction is among the longest palaeoseries produced for New Zealand to date. It is in good agreement with other tree-ring-based records, and with instrumental (both local and hemispheric) data. Thirteen pink pine chronologies were developed. Ring-width measurements were detrended using the Regional Curve Standardisation method to retain as much low-frequency variance as possible. Crossdating revealed the existence of a strong common signal among trees. Inter-site comparison indicated that a common control mechanism affected tree growth not only within sites, but also across sites. To determine whether climate was the main factor that controlled the growth of pink pine in Westland, correlation and response function analyses were employed. Temperature, precipitation and the Southern Oscillation Index were tested for their relationship with tree growth. Mean monthly temperature was identified as the primary growth-limiting factor. Chronologies were positively correlated with temperature over an extended period (5-17 months), and climate response modelling showed that temperature explained 11-60% variance in the tree-ring data. The highest and most stable correlations occurred between tree growth and summer (January-March) temperatures. Tree-ring data from the six sites that contained the strongest temperature signal were combined, and the Westland Regional Chronology (WRC) was developed. The WRC was then used to reconstruct January-March temperatures back to A.D. 1480. The calibration model explained 43% of the variance in temperature, and all calibration and verification tests were passed at high levels of significance. The reconstruction showed that temperatures in Westland have been following a positive trend over the last 520 years. The coolest 25-year period was 1542-1566, while temperatures reached their maximum in 1966-1990. Spectral analysis of the Westland palaeotemperature record revealed cycles at periods of about 3, 5-6, 11, 14, 22, 45 and 125 years. This study also confirmed that climate response is species-dependent. A separate exercise, which compared two species from the same site, demonstrated that while pink pine's growth was mainly influenced by summer temperatures, Libocedrus bidwillii was affected by conditions at the beginning of the growing season. However, the temperature signal in Westland's Libocedrus bidwillii was insufficient to produce a reliable reconstruction. It might be because the climate signal in this species was obscured by disturbances, as was shown in the final section of this project. Frequent growth releases and suppressions implied that Libocedrus bidwillii integrated both major (Alpine Fault earthquakes) and minor (windthrow) disturbances in its ring widths. Pink pine, on the other hand, was not sensitive to disturbance, and was therefore a better indicator of palaeotemperatures in Westland. This research has strengthened the New Zealand network of chronology sites, and confirmed that pink pine has great dendroclimatic value. The last 520 years of temperature fluctuations were reconstructed with a high degree of fidelity - the model developed in this thesis is currently the most accurate estimate of a temperature-growth relationship in the country.

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research Papers, Lincoln University

Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults. The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains. The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future. Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case. This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features. Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.

Research papers, University of Canterbury Library

In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.

Research papers, University of Canterbury Library

The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.

Research papers, The University of Auckland Library

The research presented in this thesis investigated the environmental impacts of structural design decisions across the life of buildings located in seismic regions. In particular, the impacts of expected earthquake damage were incorporated into a traditional life cycle assessment (LCA) using a probabilistic method, and links between sustainable and resilient design were established for a range of case-study buildings designed for different seismic performance objectives. These links were quantified using a metric herein referred to as the seismic carbon risk, which represents the expected environmental impacts and resource use indicators associated with earthquake damage during buildings’ life. The research was broken into three distinct parts: (1) a city-level evaluation of the environmental impacts of demolitions following the 2010/2011 Canterbury earthquake sequence in New Zealand, (2) the development of a probabilistic framework to incorporate earthquake damage into LCA, and (3) using case-study buildings to establish links between sustainable and resilient design. The first phase of the research focused on the environmental impacts of demolitions in Christchurch, New Zealand following the 2010/2011 Canterbury Earthquake Sequence. This large case study was used to investigate the environmental impact of the demolition of concrete buildings considering the embodied carbon and waste stream distribution. The embodied carbon was considered here as kilograms of CO2 equivalent that occurs on production, construction, and waste management stage. The results clearly demonstrated the significant environmental impacts that can result from moderate and large earthquakes in urban areas, and the importance of including environmental considerations when making post-earthquake demolition decisions. The next phase of the work introduced a framework for incorporating the impacts of expected earthquake damage based on a probabilistic approach into traditional LCA to allow for a comparison of seismic design decisions using a carbon lens. Here, in addition to initial construction impacts, the seismic carbon risk was quantified, including the impacts of seismic repair activities and total loss scenarios assuming reconstruction in case of non-reparability. A process-based LCA was performed to obtain the environmental consequence functions associated with structural and non-structural repair activities for multiple environmental indicators. In the final phase of the work, multiple case-study buildings were used to investigate the seismic consequences of different structural design decisions for buildings in seismic regions. Here, two case-study buildings were designed to multiple performance objectives, and the upfront carbon costs, and well as the seismic carbon risk across the building life were compared. The buildings were evaluated using the framework established in phase 2, and the results demonstrated that the seismic carbon risk can significantly be reduced with only minimal changes to the upfront carbon for buildings designed for a higher base shear or with seismic protective systems. This provided valuable insight into the links between resilient and sustainable design decisions. Finally, the results and observations from the work across the three phases of research described above were used to inform a discussion on important assumptions and topics that need to be considered when quantifying the environmental impacts of earthquake damage on buildings. These include: selection of a non-repairable threshold (e.g. a value beyond which a building would be demolished rather than repaired), the time value of carbon (e.g. when in the building life the carbon is released), the changing carbon intensity of structural materials over time, and the consideration of deterministic vs. probabilistic results. Each of these topics was explored in some detail to provide a clear pathway for future work in this area.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Structural engineering is facing an extraordinarily challenging era. These challenges are driven by the increasing expectations of modern society to provide low-cost, architecturally appealing structures which can withstand large earthquakes. However, being able to avoid collapse in a large earthquake is no longer enough. A building must now be able to withstand a major seismic event with negligible damage so that it is immediately occupiable following such an event. As recent earthquakes have shown, the economic consequences of not achieving this level of performance are not acceptable. Technological solutions for low-damage structural systems are emerging. However, the goal of developing a low-damage building requires improving the performance of both the structural skeleton and the non-structural components. These non-structural components include items such as the claddings, partitions, ceilings and contents. Previous research has shown that damage to such items contributes a disproportionate amount to the overall economic losses in an earthquake. One such non-structural element that has a history of poor performance is the external cladding system, and this forms the focus of this research. Cladding systems are invariably complicated and provide a number of architectural functions. Therefore, it is important than when seeking to improve their seismic performance that these functions are not neglected. The seismic vulnerability of cladding systems are determined in this research through a desktop background study, literature review, and postearthquake reconnaissance survey of their performance in the 2010 – 2011 Canterbury earthquake sequence. This study identified that precast concrete claddings present a significant life-safety risk to pedestrians, and that the effect they have upon the primary structure is not well understood. The main objective of this research is consequently to better understand the performance of precast concrete cladding systems in earthquakes. This is achieved through an experimental campaign and numerical modelling of a range of precast concrete cladding systems. The experimental campaign consists of uni-directional, quasi static cyclic earthquake simulation on a test frame which represents a single-storey, single-bay portion of a reinforced concrete building. The test frame is clad with various precast concrete cladding panel configurations. A major focus is placed upon the influence the connection between the cladding panel and structural frame has upon seismic performance. A combination of experimental component testing, finite element modelling and analytical derivation is used to develop cladding models of the cladding systems investigated. The cyclic responses of the models are compared with the experimental data to evaluate their accuracy and validity. The comparison shows that the cladding models developed provide an excellent representation of real-world cladding behaviour. The cladding models are subsequently applied to a ten-storey case-study building. The expected seismic performance is examined with and without the cladding taken into consideration. The numerical analyses of the case-study building include modal analyses, nonlinear adaptive pushover analyses, and non-linear dynamic seismic response (time history) analyses to different levels of seismic hazard. The clad frame models are compared to the bare frame model to investigate the effect the cladding has upon the structural behaviour. Both the structural performance and cladding performance are also assessed using qualitative damage states. The results show a poor performance of precast concrete cladding systems is expected when traditional connection typologies are used. This result confirms the misalignment of structural and cladding damage observed in recent earthquake events. Consequently, this research explores the potential of an innovative cladding connection. The outcomes from this research shows that the innovative cladding connection proposed here is able to achieve low-damage performance whilst also being cost comparable to a traditional cladding connection. It is also theoretically possible that the connection can provide a positive value to the seismic performance of the structure by adding addition strength, stiffness and damping. Finally, the losses associated with both the traditional and innovative cladding systems are compared in terms of tangible outcomes, namely: repair costs, repair time and casualties. The results confirm that the use of innovative cladding technology can substantially reduce the overall losses that result from cladding damage.

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

Reinforced concrete structures designed in pre-1970s are vulnerable under earthquakes due to lack of seismic detailing to provide adequate ductility. Typical deficiencies of pre-1970s reinforced concrete structures are (a) use of plain bars as longitudinal reinforcement, (b) inadequate anchorage of beam longitudinal reinforcement in the column (particularly exterior column), (c) lack of joint transverse reinforcement if any, (d) lapped splices located just above joint, and (e) low concrete strength. Furthermore, the use of infill walls is a controversial issue because it can help to provide additional stiffness to the structure on the positive side and on the negative side it can increase the possibility of soft-storey mechanisms if it is distributed irregularly. Experimental research to investigate the possible seismic behaviour of pre-1970s reinforced concrete structures have been carried out in the past. However, there is still an absence of experimental tests on the 3-D response of existing beam-column joints under bi-directional cyclic loading, such as corner joints. As part of the research work herein presented, a series of experimental tests on beam-column subassemblies with typical detailing of pre-1970s buildings has been carried out to investigate the behaviour of existing reinforced concrete structures. Six two-third scale plane frame exterior beam-column joint subassemblies were constructed and tested under quasi-static cyclic loading in the Structural Laboratory of the University of Canterbury. The reinforcement detailing and beam dimension were varied to investigate their effect on the seismic behaviour. Four specimens were conventional deep beam-column joint, with two of them using deformed longitudinal bars and beam bars bent in to the joint and the two others using plain round longitudinal bars and beam bars with end hooks. The other two specimens were shallow beam-column joint, one with deformed longitudinal bars and beam bars bent in to the joint, the other with plain round longitudinal bars and beam bars with end hooks. All units had one transverse reinforcement in the joint. The results of the experimental tests indicated that conventional exterior beam-column joint with typical detailing of pre-1970s building would experience serious diagonal tension cracking in the joint panel under earthquake. The use of plain round bars with end hooks for beam longitudinal reinforcement results in more severe damage in the joint core when compared to the use of deformed bars for beam longitudinal reinforcement bent in to the joint, due to the combination of bar slips and concrete crushing. One interesting outcome is that the use of shallow beam in the exterior beam-column joint could avoid the joint cracking due to the beam size although the strength provided lower when compared with the use of deep beam with equal moment capacity. Therefore, taking into account the low strength and stiffness, shallow beam can be reintroduced as an alternative solution in design process. In addition, the presence of single transverse reinforcement in the joint core can provide additional confinement after the first crack occurred, thus delaying the strength degradation of the structure. Three two-third scale space frame corner beam-column joint subassemblies were also constructed to investigate the biaxial loading effect. Two specimens were deep-deep beam-corner column joint specimens and the other one was deep-shallow beam-corner column joint specimen. One deep-deep beam-corner column joint specimen was not using any transverse reinforcement in the joint core while the two other specimens were using one transverse reinforcement in the joint core. Plain round longitudinal bars were used for all units with hook anchorage for the beam bars. Results from the tests confirmed the evidences from earthquake damage observations with the exterior 3-D (corner) beam-column joint subjected to biaxial loading would have less strength and suffer higher damage in the joint area under earthquake. Furthermore, the joint shear relation in the two directions is calibrated from the results to provide better analysis. An analytical model was used to simulate the seismic behaviour of the joints with the help of Ruaumoko software. Alternative strength degradation curves corresponding to different reinforcement detailing of beam-column joint unit were proposed based on the test results.

Research papers, University of Canterbury Library

This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.

Research papers, University of Canterbury Library

Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.

Research papers, The University of Auckland Library

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.

Research papers, University of Canterbury Library

Liquefaction of sandy soil has been observed to cause significant damage to infrastructure during major earthquakes. Historical cases of liquefaction have typically occurred in sands containing some portion of fines particles, which are defined as 75μm or smaller in diameter. The effects of fines on the undrained behaviour of sand are not however fully understood, and this study therefore attempts to quantify these effects through the undrained testing of sand mixed with non-plastic fines sourced from Christchurch, New Zealand. The experimental program carried out during this study consisted of undrained monotonic and cyclic triaxial tests performed on three different mixtures of sand and fines: the Fitzgerald Bridge mixture (FBM), and two Pinnacles Sand mixtures (PSM1 and PSM2). The fines content of each host sand was systematically varied up to a maximum of 30%, with all test specimens being reconstituted using moist tamping deposition. The undrained test results from the FBM soils were interpreted using a range of different measures of initial state. When using void ratio and relative density, the addition of fines to the FBM sand caused more contractive behaviour for both monotonic and cyclic loadings. This resulted in lower strengths at the steady state of deformation, and lower liquefaction resistances. When the intergranular void ratio was used for the interpretation, the effect of additional fines was to cause less contractive response in the sand. The state parameter and state index were also used to interpret the undrained cyclic test results – these measures suggested that additional fines caused less contractive sand behaviour, the opposite to that observed when using the void ratio. This highlighted the dependency on the parameter chosen as a basis for the response comparison when determining the effects of fines, and pointed out a need to identify a measure that normalizes such effects. Based on the FBM undrained test results and interpretations, the equivalent granular void ratio, e*, was identified from the literature as a measure of initial state that normalizes the effects of fines on the undrained behaviour of sand up to a fines content of 30%. This is done through a parameter within the e* definition termed the fines influence factor, b, which quantifies the effects of fines from a value of zero (no effect) to one (same effect as sand particles). The value of b was also determined to be different when interpreting the steady state lines (bSSL) and cyclic resistance curves (bCR) respectively for a given mixture of sand and fines. The steady state lines and cyclic resistance curves of the FBM soils and a number of other sand-fines mixtures sourced from the literature were subsequently interpreted using the equivalent granular void ratio concept, with bSSL and bCR values being back-calculated from the respective test data sets. Based on these interpretations, it was concluded that e* was conceptually a useful parameter for characterizing and quantifying the effects of fines on the undrained behaviour of sand, assuming the fines influence factor value could be derived. To allow prediction of the fines influence factor values, bSSL and bCR were correlated with material and depositional properties of the presented sand-fines mixtures. It was found that as the size of the fines particles relative to the sand particles became smaller, the values of bSSL and bCR reduced, indicating lower effect of fines. The same trend was also observed as the angularity of the sand particles increased. The depositional method was found to influence the value of bCR, due to the sensitivity of cyclic loading to initial soil fabric. This led to bSSL being used as a reference for the effect of fines, with specimens prepared by moist tamping having bCR > bSSL, and specimens prepared by slurry deposition having bCR < bSSL. Finally the correlations of the fines influence factor values with material and depositional properties were used to define the simplified estimation method – a procedure capable of predicting the approximate steady state lines and cyclic resistance curves of a sand as the non-plastic fines content is increased up to 30%. The method was critically reviewed based on the undrained test results of the PSM1 and PSM2 soils. This review suggested the method could accurately predict undrained response curves as the fines content was raised, based on the PSM1 test results. It also however identified some key issues with the method, such as the inability to accurately predict the responses of highly non-uniform soils, a lack of consideration for the entire particle size distribution of a soil, and the fact the errors in the prediction of bSSL carry through into the prediction of bCR. Lastly some areas of further investigation relating to the method were highlighted, including the need to verify the method through testing of sandy soils sourced from outside the Christchurch area, and the need to correlate the value of bCR with additional soil fabrics / depositional methods.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Background The 2010/2011 Canterbury earthquakes and aftershocks in New Zealand caused unprecedented destruction to the physical, social, economic, and community fabric of Christchurch city. The recovery phase in Christchurch is on going, six years following the initial earthquake. Research exploring how disabled populations experience community inclusion in the longer-term recovery following natural disasters is scant. Yet such information is vital to ensure that recovering communities are inclusive for all members of the affected population. This thesis specifically examined how people who use wheelchairs experienced community inclusion four years following the 2010/2011 Canterbury earthquakes. Aims The primary research aim was to understand how one section of the disability community – people who use wheelchairs – experienced community inclusion over the four years following the 2010/2011 Canterbury earthquakes and aftershocks. A secondary aim was to test a novel sampling approach, Respondent Driven Sampling, which had the potential to enable unbiased population-based estimates. This was motivated by the lack of an available sampling frame for the target population, which would inhibit recruitment of a representative sample. Methodology and methods An exploratory sequential mixed methods design was used, beginning with a qualitative phase (Phase One), which informed a second quantitative phase (Phase Two). The qualitative phase had two stages. First, a small sample of people who use wheelchairs participated in an individual, semi-structured interview. In the second stage, these participants were then invited to a group interview to clarify and prioritise themes identified in the individual interviews. The quantitative phase was a cross-sectional survey developed from the findings from Phase One. Initially, Respondent Driven Sampling was employed to conduct a national, electronic cross-sectional survey that aimed to recruit a sample that may provide unbiased population-based estimates. Following the unsuccessful application of Respondent Driven Sampling, a region-specific convenience sampling approach was used. The datasets from the qualitative and quantitative phases were integrated to address the primary aim of the research. Results In Phase One 13 participants completed the individual interviews, and five of them contributed to the group interview. Thematic analysis of individual and group interview data suggested that participants felt the 2010/11 earthquakes magnified many pre-existing barriers to community inclusion, and also created an exciting opportunity for change. This finding was encapsulated in five themes: 1) earthquakes magnified barriers, 2) community inclusion requires energy, 3) social connections are important, 4) an opportunity lost, and 5) an opportunity found. The findings from Phase One informed the development of a survey instrument to investigate how these findings generalised to a larger sample of individuals who use wheelchairs. In Phase Two, the Respondent Driven Sampling approach failed to recruit enough participants to satisfy the statistical requirements needed to reach equilibrium, thereby enabling the calculation of unbiased population estimates. The subsequent convenience sampling approach recruited 49 participants who, combined with the 15 participants from the Respondent Driven Sampling approach that remained eligible for the region-specific sample, resulted in the total of 64 individuals who used wheelchairs and were residents of Christchurch. Participants reported their level of community inclusion at three time periods: the six months prior to the first earthquake in September 2010 (time one), the six months following the first earthquake in September 2010 (time two), and the six months prior to survey completion (between October 2015 and March 2016, (time three)). Survey data provided some precision regarding the timing in which the magnified barriers developed. Difficulty with community inclusion rose significantly between time one and time two, and while reducing slightly, was still present during time three, and had not returned to the time one baseline. The integrated findings from Phase One and Phase Two suggested that magnified barriers to community inclusion had been sustained four years post-earthquake, and community access had not returned to pre-earthquake levels, let alone improved beyond pre-earthquake levels. Conclusion Findings from this mixed methods study suggest that four years following the initial earthquake, participants were still experiencing multiple magnified barriers, which contributed to physical and social exclusion, as well as fatigue, as participants relied on individual agency to negotiate such barriers. Participants also highlighted the exciting opportunity to create an accessible city. However because they were still experiencing barriers four years following the initial event, and were concerned that this opportunity might be lost if the recovery proceeds without commitment and awareness from the numerous stakeholders involved in guiding the recovery. To truly realise the opportunity to create an accessible city following a disaster, the transition from the response phase to a sustainable longer-term recovery must adopt a new model of community engagement where decision-makers partner with people living with disability to co-produce a vision and strategy for creating an inclusive community. Furthermore, despite the unsuccessful use of Respondent Driven Sampling in this study, future research exploring the application of RDS with wheelchair users is recommended before discounting this sampling approach in this population.

Research papers, Victoria University of Wellington

<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.

Research Papers, Lincoln University

The city of Ōtautahi/Christchurch experienced a series of earthquakes that began on September 4th, 2010. The most damaging event occurred on February 22nd, 2011 but significant earthquakes also occurred on June 13th and December 23rd with aftershocks still occurring well into 2012. The resulting disaster is the second deadliest natural disaster in New Zealand’s history with 185 deaths. During 2011 the Canterbury earthquakes were one of the costliest disasters worldwide with an expected cost of up to $NZ30 billion. Hundreds of commercial buildings and thousands of houses have been destroyed or are to be demolished and extensive repairs are needed for infrastructure to over 100,000 homes. As many as 8,900 people simply abandoned their homes and left the city in the first few months after the February event (Newell, 2012), and as many as 50,000 may leave during 2012. In particular, young whānau and single young women comprised a disproportionate number of these migrants, with evidence of a general movement to the North Island. Te Puni Kōkiri sought a mix of quantitative and qualitative research to examine the social and economic impacts of the Christchurch earthquakes on Māori and their whānau. The result of this work will be a collection of evidence to inform policy to support and assist Māori and their whānau during the recovery/rebuild phases. To that end, this report triangulates available statistical and geographical information with qualitative data gathered over 2010 and 2011 by a series of interviews conducted with Māori who experienced the dramatic events associated with the earthquakes. A Māori research team at Lincoln University was commissioned to undertake the research as they were already engaged in transdisciplinary research (began in the May 2010), that focused on quickly gathering data from a range of Māori who experienced the disaster, including relevant economic, environmental, social and cultural factors in the response and recovery of Māori to these events. Participants for the qualitative research were drawn from Māori whānau who both stayed and left the city. Further data was available from ongoing projects and networks that the Lincoln research team was already involved in, including interviews with Māori first responders and managers operating in the CBD on the day of the February event. Some limited data is also available from younger members of affected whānau. Māori in Ōtautahi/Christchurch City have exhibited their own culturally-attuned collective responses to the disaster. However, it is difficult to ascertain Māori demographic changes due to a lack of robust statistical frameworks but Māori outward migration from the city is estimated to range between 560 and 1,100 people. The mobility displayed by Māori demonstrates an important but unquantified response by whānau to this disaster, with emigration to Australia presenting an attractive option for young Māori, an entrenched phenomenon that correlates to cyclical downturns and the long-term decline of the New Zealand economy. It is estimated that at least 315 Māori have emigrated from the Canterbury region to Australia post-quake, although the disaster itself may be only one of a series of events that has prompted such a decision. Māori children made up more than one in four of the net loss of children aged 6 to 15 years enrolled in schools in Greater Christchurch over the year to June 2011. Research literature identifies depression affecting a small but significant number of children one to two years post-disaster and points to increasing clinical and organisational demands for Māori and other residents of the city. For those residents in the eastern or coastal suburbs – home to many of the city’s Māori population - severe damage to housing, schools, shops, infrastructure, and streets has meant disruption to their lives, children’s schooling, employment, and community functioning. Ongoing abandonment of homes by many has meant a growing sense of unease and loss of security, exacerbated by arson, burglaries, increased drinking, a stalled local and national economy, and general confusion about the city’s future. Māori cultural resilience has enabled a considerable network of people, institutions, and resources being available to Māori , most noticeably through marae and their integral roles of housing, as a coordinating hub, and their arguing for the wider affected communities of Christchurch. Relevant disaster responses need to be discussed within whānau, kōhanga, kura, businesses, communities, and wider neighbourhoods. Comprehensive disaster management plans need to be drafted for all iwi in collaboration with central government, regional, and city or town councils. Overall, Māori are remarkably philosophical about the effects of the disaster, with many proudly relishing their roles in what is clearly a historic event of great significance to the city and country. Most believe that ‘being Māori’ has helped cope with the disaster, although for some this draws on a collective history of poverty and marginalisation, features that contribute to the vulnerability of Māori to such events. While the recovery and rebuild phases offer considerable options for Māori and iwi, with Ngāi Tahu set to play an important stakeholder in infrastructural, residential, and commercial developments, some risk and considerable unknowns are evident. Considerable numbers of Māori may migrate into the Canterbury region for employment in the rebuild, and trades training strategies have already been established. With many iwi now increasingly investing in property, the risks from significant earthquakes are now more transparent, not least to insurers and the reinsurance sector. Iwi authorities need to be appraised of insurance issues and ensure sufficient coverage exists and investments and developments are undertaken with a clear understanding of the risks from natural hazards and exposure to future disasters.

Research papers, University of Canterbury Library

This thesis studies the behaviour of diaphragms in multi-storey timber buildings by providing methods for the estimation of the diaphragm force demand, developing an Equivalent Truss Method for the analysis of timber diaphragms, and experimentally investigating the effects of displacement incompatibilities between the diaphragm and the lateral load resisting system and developing methods for their mitigation. The need to better understand the behaviour of diaphragms in timber buildings was highlighted by the recent 2010-2011 Canterbury Earthquake series, where a number of diaphragms in traditional concrete buildings performed poorly, compromising the lateral load resistance of the structure. Although shortcomings in the estimation of force demand, and in the analysis and design of concrete floor diaphragms have already been partially addressed by other researchers, the behaviour of diaphragms in modern multi-storey timber buildings in general, and in low damage Pres-Lam buildings (consisting of post-tensioned timber members) in particular is still unknown. The recent demand of mid-rise commercial timber buildings of ten storeys and beyond has further highlighted the lack of appropriate methods to analyse timber diaphragms with irregular floor geometries and large spans made of both light timber framing and massive timber panels. Due to the lower stiffness of timber lateral load resisting systems, compared with traditional construction materials, and the addition of in-plane flexible diaphragms, the effect of higher modes on the global dynamic behaviour of a structure becomes more critical. The results from a parametric non-linear time-history analysis on a series of timber frame and wall structures showed increased storey shear and moment demands even for four storey structures when compared to simplistic equivalent static analysis. This effect could successfully be predicted with methods available in literature. The presence of diaphragm flexibility increased diaphragm inter-storey drifts and the peak diaphragm demand in stiff wall structures, but had less influence on the storey shears and moments. Diaphragm force demands proved to be significantly higher than the forces derived from equivalent static analysis, leading to potentially unsafe designs. It is suggested to design all diaphragms for the same peak demand; a simplified approach to estimate these diaphragm forces is proposed for both frame and wall structures. Modern architecture often requires complex floor geometries with long spans leading to stress concentrations, high force demands and potentially large deformations in the diaphragms. There is a lack of guidance and regulation regarding the analysis and design of timber diaphragms and a practical alternative to the simplistic equivalent deep beam analysis or costly finite element modelling is required. An Equivalent Truss Method for the analysis of both light timber framed and massive timber diaphragms is proposed, based on analytical formulations and verified against finite element models. With this method the panel unit shear forces (shear flow) and therefore the fastener demand, chord forces and reaction forces can be evaluated. Because the panel stiffness and fastener stiffness are accounted for, diaphragm deflection, torsional effects and transfer forces can also be assessed. The proposed analysis method is intuitive and can be used with basic analysis software. If required, it can easily be adapted for the use with diaphragms working in the non-linear range. Damage to floor diaphragms resulting from displacement incompatibilities due to frame elongation or out-of plane deformation of walls can compromise the transfer of inertial forces to the lateral load resisting system as well as the stability of other structural elements. Two post-tensioned timber frame structures under quasi-static cyclic and dynamic load, respectively, were tested with different diaphragm panel layouts and connections investigating their ability to accommodate frame elongations. Additionally, a post-tensioned timber wall was loaded under horizontal cyclic loads through two pairs of collector beams. Several different connection details between the wall and the beams were tested, and no damage to the collector beams or connections was observed in any of the tests. To evaluate the increased strength and stiffness due to the wall-beam interaction an analytical procedure is presented. Finally, a timber staircase core was tested under bi-directional loading. Different connection details were used to study the effect of displacement incompatibilities between the orthogonal collector beams. These experiments showed that floor damage due to displacement incompatibilities can be prevented, even with high levels of lateral drift, by the flexibility of well-designed connections and the flexibility of the timber elements. It can be concluded that the flexibility of timber members and the flexibility of their connections play a major role in the behaviour of timber buildings in general and of diaphragms specifically under seismic loads. The increased flexibility enhances higher mode effects and alters the diaphragm force demand. Simple methods are provided to account for this effect on the storey shear, moment and drift demands as well as the diaphragm force demands. The analysis of light timber framing and massive timber diaphragms can be successfully analysed with an Equivalent Truss Method, which is calibrated by accounting for the panel shear and fastener stiffnesses. Finally, displacement incompatibilities in frame and wall structures can be accommodated by the flexibilities of the diaphragm panels and relative connections. A design recommendations chapter summarizes all findings and allows a designer to estimate diaphragm forces, to analyse the force path in timber diaphragms and to detail the connections to allow for displacement incompatibilities in multi-storey timber buildings.

Research papers, University of Canterbury Library

Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.

Research Papers, Lincoln University

Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.

Research papers, The University of Auckland Library

In September 2010 and February 2011 the Canterbury region of New Zealand was struck by two powerful earthquakes, registering magnitude 7.1 and 6.3 respectively on the Richter scale. The second earthquake was centred 10 kilometres south-east of the centre of Christchurch (the region’s capital and New Zealand’s third most populous urban area, with approximately 360,000 residents) at a depth of five kilometres. 185 people were killed, making it the second deadliest natural disaster in New Zealand’s history. (66 people were killed in the collapse of one building alone, the six-storey Canterbury Television building.) The earthquake occurred during the lunch hour, increasing the number of people killed on footpaths and in buses and cars by falling debris. In addition to the loss of life, the earthquake caused catastrophic damage to both land and buildings in Christchurch, particularly in the central business district. Many commercial and residential buildings collapsed in the tremors; others were damaged through soil liquefaction and surface flooding. Over 1,000 buildings in the central business district were eventually demolished because of safety concerns, and an estimated 70,000 people had to leave the city after the earthquakes because their homes were uninhabitable. The New Zealand Government declared a state of national emergency, which stayed in force for ten weeks. In 2014 the Government estimated that the rebuild process would cost NZ$40 billion (approximately US$27.3 billion, a cost equivalent to 17% of New Zealand’s annual GDP). Economists now estimate it could take the New Zealand economy between 50 and 100 years to recover. The earthquakes generated tens of thousands of insurance claims, both against private home insurance companies and against the New Zealand Earthquake Commission, a government-owned statutory body which provides primary natural disaster insurance to residential property owners in New Zealand. These ranged from claims for hundreds of millions of dollars concerning the local port and university to much smaller claims in respect of the thousands of residential homes damaged. Many of these insurance claims resulted in civil proceedings, caused by disputes about policy cover, the extent of the damage and the cost and/or methodology of repairs, as well as failures in communication and delays caused by the overwhelming number of claims. Disputes were complicated by the fact that the Earthquake Commission provides primary insurance cover up to a monetary cap, with any additional costs to be met by the property owner’s private insurer. Litigation funders and non-lawyer claims advocates who took a percentage of any insurance proceeds also soon became involved. These two factors increased the number of parties involved in any given claim and introduced further obstacles to resolution. Resolving these disputes both efficiently and fairly was (and remains) central to the rebuild process. This created an unprecedented challenge for the justice system in Christchurch (and New Zealand), exacerbated by the fact that the Christchurch High Court building was itself damaged in the earthquakes, with the Court having to relocate to temporary premises. (The High Court hears civil claims exceeding NZ$200,000 in value (approximately US$140,000) or those involving particularly complex issues. Most of the claims fell into this category.) This paper will examine the response of the Christchurch High Court to this extraordinary situation as a case study in innovative judging practices and from a jurisprudential perspective. In 2011, following the earthquakes, the High Court made a commitment that earthquake-related civil claims would be dealt with as swiftly as the Court's resources permitted. In May 2012, it commenced a special “Earthquake List” to manage these cases. The list (which is ongoing) seeks to streamline the trial process, resolve quickly claims with precedent value or involving acute personal hardship or large numbers of people, facilitate settlement and generally work proactively and innovatively with local lawyers, technical experts and other stakeholders. For example, the Court maintains a public list (in spreadsheet format, available online) with details of all active cases before the Court, listing the parties and their lawyers, summarising the facts and identifying the legal issues raised. It identifies cases in which issues of general importance have been or will be decided, with the expressed purpose being to assist earthquake litigants and those contemplating litigation and to facilitate communication among parties and lawyers. This paper will posit the Earthquake List as an attempt to implement innovative judging techniques to provide efficient yet just legal processes, and which can be examined from a variety of jurisprudential perspectives. One of these is as a case study in the well-established debate about the dialogic relationship between public decisions and private settlement in the rule of law. Drawing on the work of scholars such as Hazel Genn, Owen Fiss, David Luban, Carrie Menkel-Meadow and Judith Resnik, it will explore the tension between the need to develop the law through the doctrine of precedent and the need to resolve civil disputes fairly, affordably and expeditiously. It will also be informed by the presenter’s personal experience of the interplay between reported decisions and private settlement in post-earthquake Christchurch through her work mediating insurance disputes. From a methodological perspective, this research project itself gives rise to issues suitable for discussion at the Law and Society Annual Meeting. These include the challenges in empirical study of judges, working with data collected by the courts and statistical analysis of the legal process in reference to settlement. September 2015 marked the five-year anniversary of the first Christchurch earthquake. There remains widespread dissatisfaction amongst Christchurch residents with the ongoing delays in resolving claims, particularly insurers, and the rebuild process. There will continue to be challenges in Christchurch for years to come, both from as-yet unresolved claims but also because of the possibility of a new wave of claims arising from poor quality repairs. Thus, a final purpose of presenting this paper at the 2016 Meeting is to gain the benefit of other scholarly perspectives and experiences of innovative judging best practice, with a view to strengthening and improving the judicial processes in Christchurch. This Annual Meeting of the Law and Society Association in New Orleans is a particularly appropriate forum for this paper, given the recent ten year anniversary of Hurricane Katrina and the plenary session theme of “Natural and Unnatural Disasters – human crises and law’s response.” The presenter has a personal connection with this theme, as she was a Fulbright scholar from New Zealand at New York University in 2005/2006 and participated in the student volunteer cleanup effort in New Orleans following Katrina. http://www.lawandsociety.org/NewOrleans2016/docs/2016_Program.pdf

Research Papers, Lincoln University

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