Search

found 280 results

Images for present; more images...

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