Search

found 177 results

Research papers, University of Canterbury Library

This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30

Research papers, University of Canterbury Library

Rising disaster losses, growth in global migration, migrant labour trends, and increasingly diverse populations have serious implications for disaster resilience around the world. These issues are of particular concern in New Zealand, which is highly exposed to disaster risk and has the highest proportion of migrant workers to national population in the OECD. Since there has been no research conducted into this issue in New Zealand to date, greater understanding of the social capital used by migrant workers in specific New Zealand contexts is needed to inform more targeted and inclusive disaster risk management approaches. A New Zealand case study is used to investigate the extent and types of social capital and levels of disaster risk awareness reported by members of three Filipino migrant workers organisations catering to dairy farm, construction and aged care workers in different urban and rural Canterbury districts. Findings from (3) semi-structured interviews and (3) focus groups include consistently high reliance on bonding capital and low levels of bridging capital across all three organisations and industry sectors, and in both urban and rural contexts. The transitory, precarious residential status conveyed by temporary work visas, and the difficulty of building bridging capital with host communities has contributed to this heavy reliance on bonding capital. Social media was essential to connect workers with family and friends in other countries, while Filipino migrant workers organisations provided members with valuable access to industry and district-specific networks of other Filipino migrant workers. Linking capital varied between the three organisations, with members of the organisation set up to advocate for dairy farm workers reporting the highest levels of linking capital. Factors influencing the capacity of workers organisations to develop linking capital appeared to include motivation (establishment objectives), length of time since establishment, support from government and industry groups, urban-rural context, income levels and gender. Although aware of publicity around earthquake and tsunami risk in the Canterbury region, participants were less aware of flood risk, and expressed fatalistic attitudes to disaster risk. Workers organisations offer a valuable potential interface between CDEM Group activities and migrant worker communities, since organisation leaders were interested in accessing government support to participate (with and on behalf of members) in disaster risk planning at district and regional level. With the potential to increase disaster resilience among these vulnerable, hard to reach communities, such participation could also help to build capacity across workers organisations (within Canterbury and across the country) to develop linking capital at national, as well as regional level. However, these links will also depend on greater government and industry commitment to providing more targeted and appropriate support for migrant workers, including consideration of the cultural qualifications of staff tasked with liaising with this community.

Research papers, Victoria University of Wellington

<b>New Zealand has experienced several strong earthquakes in its history. While an earthquake cannot be prevented from occurring, planning can reduce its consequences when it does occur. This dissertation research examines various aspects of disaster risk management policy in Aotearoa New Zealand.</b> Chapter 2 develops a method to rank and prioritise high-rise buildings for seismic retrofitting in Wellington, the earthquake-prone capital city of New Zealand. These buildings pose risks to Wellington’s long-term seismic resilience that are of clear concern to current and future policymakers. The prioritization strategy we propose, based on multi-criteria decision analysis (MCDA) methods, considers a variety of data on each building, including not only its structural characteristics, but also its location, its economic value to the city, and its social importance to the community around it. The study demonstrates how different measures, within four general criteria – life safety, geo-spatial location of the building, its economic role, and its socio-cultural role – can be operationalized into a viable framework for determining retrofitting/demolition policy priorities. Chapter 3 and chapter 4 analyse the Residential Red Zone (RRR) program that was implemented in Christchurch after the 2011 earthquake. In the program, approximately 8,000 homeowners were told that their homes were no longer permittable, and they were bought by the government (through the Canterbury Earthquake Recovery Authority). Chapter 3 examines the subjective wellbeing of the RRR residents (around 16000 people) after they were forced to move. We consider three indicators of subjective wellbeing: quality of life, stress, and emotional wellbeing. We found that demographic factors, health conditions, and the type of government compensation the residents accepted, were all significant determinants of the wellbeing of the Red Zone residents. More social relations, better financial circumstances, and the perception of better government communication were also all associated positively with a higher quality of life, less stress, and higher emotional wellbeing. Chapter 4 concentrates on the impact of this managed retreat program on RRR residents’ income. We use individual-level comprehensive, administrative, panel data from Canterbury, and difference in difference evaluation method to explore the effects of displacement on Red Zone residential residents. We found that compared to non-relocated neighbours, the displaced people experience a significant initial decrease in their wages and salaries, and their total income. The impacts vary with time spent in the Red Zone and when they moved away. Wages and salaries of those who were red-zoned and moved in 2011 were reduced by 8%, and 5.4% for those who moved in 2012. Females faced greater decreases in wages and salaries, and total income, than males. There were no discernible impacts of the relocation on people’s self-employment income.

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, Lincoln University

The major earthquakes of 2010 and 2011 brought to an abrupt end a process of adaptive reuse, revitalisation and gentrification that was underway in the early 20th century laneways and buildings located in the south eastern corner of the Christchurch Central Business District. Up until then, this location was seen as an exemplar of how mixed use could contribute to making the central city an attractive and viable alternative to the suburban living experience predominant in New Zealand. This thesis is the result of a comprehensive case study of this “Lichfield Lanes” area, which involved in depth interviews with business owners, observation of public meetings and examination of documents and the revitalisation research literature. Findings were that many of the factors seen to make this location successful pre-earthquakes mirror the results of similar research in other cities. These factors include: the importance of building upon historic architecture and the eclectic spaces this creates; a wide variety of uses generating street life; affordable rental levels; plus the dangers of uniformity of use brought about by focussing on business types that pay the most rent. Also critical is co-operation between businesses to create and effectively market and manage an identifiable precinct that has a coherent style and ambience that differentiates the location from competing suburban malls. In relation to the latter, a significant finding of this project was that the hospitality and retail businesses key to the success of Lichfield Lanes were not typical and could be described as quirky, bohemian, chaotic, relatively low rent, owner operated and appealing to the economically important “Creative Class” identified by Richard Florida (2002) and others. In turn, success for many of these businesses can be characterised as including psychological and social returns rather than simply conventional economic benefits. This has important implications for inner city revitalisation, as it contrasts with the traditional focus of local authorities and property developers on physical aspects and tenant profitability as measures of success. This leads on to an important conclusion from this research, which is that an almost completely inverted strategy from that applied to suburban mall development, may be most appropriate for successful inner city revitalisation. It also highlights a disconnection between the focus and processes of regulatory authorities and the outcomes and processes most acceptable to the people likely to frequent the central city. Developers are often caught in the middle of this conflicted situation. Another finding was early commitment by businesses to rebuild the case study area in the same style, but over time this waned as delay, demolition, insurance problems, political and planning uncertainty plus other issues made participation by the original owners and tenants impossible or uneconomic. In conclusion, the focus of inner city revitalisation is too often on buildings rather than the people that use them and what they now desire from the central city.

Research papers, University of Canterbury Library

Current research in geotechnical engineering at the University of Canterbury includes a number of laboratory testing programmes focussed on understanding the behaviour of natural soil deposits in Christchurch during the 2010-2011 Canterbury Earthquake Sequence. Many soils found in Christchurch are sands or silty sands with little to no plasticity, making them very difficult to sample using established methods. The gel-push sampling methodology, developed by Kiso-Jiban Consultants in Japan, was developed to address some of the deficiencies of existing sampling techniques and has been deployed on two projects in Christchurch. Gel push sampling is carried out with a range of samplers which are modified versions of existing technology, and the University of Canterbury has acquired three versions of the tools (GP-S, GP-Tr, GP-D). Soil samples are extracted from the bottom of a freshly drilled borehole and are captured within a liner barrel, close to 1m in length. A lubricating polymer gel coats the outside of the soil sample as it enters the liner barrel. The frictional rubbing which normally occurs on the sides of the soil samples using existing techniques is eliminated by the presence of the polymer gel. The operation of the gel-push samplers is significantly more complicated than conventional push-tube samplers, and in the initial trials a number of operational difficulties were encountered, requiring changes to the sampling procedures. Despite these issues, a number of high quality soil samples were obtained on both projects using the GP-S sampler to capture silty soil. Attempts were made to obtain clean sands using a different gel-push sampler (GP-TR) in the Red Zone. The laboratory testing of these sands indicated that they were being significantly disturbed during the sampling and/or transportation procedures. While it remains too early to draw definitive conclusions regarding the performance of the gel-push samplers, the methodology has provided some promising results. Further trialling of the tools are required to refine operating procedures understand the full range of soil conditions which can be successfully sampled using the tools. In parallel with the gel-push trials, a Dames and Moore fixed-piston sampler has been used by our research partners from Berkeley to obtain soil samples at a number of sites within Christchurch. This sampler features relatively short (50cm), thin-walled liner barrels which is advanced into the ground under the action of hydraulic pressure. By reducing the overall length of the soil being captured, the disturbance to the soil as it enters the liner barrel is significantly reduced. The Dames and Moore sampler is significantly easier to operate than the gel-push sampler, and past experience has shown it to be successful in soft, plastic materials (i.e. clays and silty clays). The cyclic resistance of one silty clay obtained using both the gel-push and Dames & Moore samplers has been found to be very similar, and ongoing research aims to establish whether similar results are obtained for different soil types, including silty materials and clean sands.

Research papers, Victoria University of Wellington

This dissertation explores the advocacy for the Christchurch Town Hall that occurred in 2012-2015 after the Canterbury Earthquakes. It frames this advocacy as an instance of collective-action community participation in a heritage decision, and explores the types of heritage values it expressed, particularly social values. The analysis contextualises the advocacy in post-quake Christchurch, and considers its relationship with other developments in local politics, heritage advocacy, and urban activism. In doing so, this dissertation considers how collective action operates as a form of public participation, and the practical implications for understanding and recognising social value.  This research draws on studies of practices that underpin social value recognition in formal heritage management. Social value is held by communities outside institutions. Engaging with communities enables institutions to explore the values of specific places, and to realise the potential of activating local connections with heritage places. Such projects can be seen as participatory practices. However, these processes require skills and resources, and may not be appropriate for all places, communities and institutions. However, literature has understudied collective action as a form of community participation in heritage management. All participation processes have nuances of communities, processes, and context, and this dissertation analyses these in one case. The research specifically asked what heritage values (especially social values) were expressed through collective action, what the relationship was with the participation processes, communities, and wider situation that produced them, and the impact on institutional rhetoric and decisions. The research analysed values expressed in representations made to council in support of the Town Hall. It also used documentary sources and interviews with key informants to analyse the advocacy and decision-making processes and their relationships with the wider context and other grassroots activities. The analysis concluded that the values expressed intertwined social and professional values. They were related to the communities and circumstance that produced them, as an advocacy campaign for a civic heritage building from a Western architectural tradition. The advocacy value arguments were one of several factors that impacted the decision. They have had a lasting impact on rhetoric around the Town Hall, as was a heritage-making practice in its own right. This dissertation makes a number of contributions to the discussion of social value and community in heritage. It suggests connections between advocacy and participation perspectives in heritage. It recommends consideration of nuances of communities, context, and place meanings when using heritage advocacy campaigns as evidence of social value. It adds to the literature on heritage advocacy, and offers a focused analysis of one of many heritage debates that occurred in post-quake Christchurch. Ultimately, it encourages practice to actively integrate social and community values and to develop self-reflexive engagement and valuation processes. Despite inherent challenges, participatory processes offer opportunities to diversify understandings of value, co-produce heritage meanings with communities, and empower citizens in democratic processes around the places they live with and love.

Research papers, Victoria University of Wellington

In the aftermath of the 2011 earthquake, a state of polycentric urbanity was thrust upon New Zealand’s second largest city. As the city-centre lay in disrepair, smaller centres started to materialise elsewhere, out of necessity. Transforming former urban peripheries and within existing suburbs into a collective, dispersed alternative to the city centre, these sub-centres prompted a range of morphological, socio-cultural and political transformations, and begged multiple questions: how to imbue these new sub-centres with gravity? How to render them a genuine alternative to the CBD? How do they operate within the wider city? How to cope with the physical and cultural transformations of this shifting urbanscape and prevent them occurring ad lib? Indeed, the success and functioning of the larger urban structure hinges upon a critical, informed response to these sub-centre urban contexts. Yet, with an unrelenting focus on the CBD rebuild - effectively a polycentric denial - little such attention has been granted.  Taking this urban condition as its premise and its provocation, this thesis investigates architecture’s role in the emergent sub-centre. It asks: what can architecture do in these urban contexts; how can architecture act upon the emergent sub-centre in a critical, catalytic fashion? Identifying this volatile condition as both an opportunity for architectural experimentation and a need for critical architectural engagement, this thesis seeks to explore the sub-centre (as an idea and actual urban context) as architecture’s project: its raison d’etre, impetus and aspiration.  These inquiries are tested through design-led research: an initial design question provoking further, broader discursive research (and indeed, seeking broader implications). The first section is a site-specific, design for Sumner, Christchurch. Titled ‘An Agora Anew’; this project - both in conception and outcome - is a speculative response to a specific sub-centre condition. The second section ‘The Sub-centre as Architecture’s Project’ explores the ideas provoked by the design project within a discursive framework. Firstly it identifies the sub-centre as a context in desperate need of architectural attention (why architecture?); secondly, it negotiates a possible agenda for architecture in this context through terms of engagement that are formal, critical and opportunistic (how architecture?): enabling it to take a position on and in the sub-centre. Lastly, a critical exegesis positions the design in regards to the broader discursive debate: critiquing it an architectural project predicated upon the idea of the sub-centre.  The implications of this design-led thesis are twofold: firstly, for architecture’s role in the sub-centre (especially to Christchurch); secondly for the possibilities of architecture’s productive engagement with the city (largely through architectural form), more generally. In a century where radical, new urban contexts (of which the sub-centre is just one) are commonplace, this type of thinking – what can architecture do in the city? - is imperative.

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, The University of Auckland Library

The influence of nonlinear soil-foundation-structure interaction (SFSI) on the performance of multi-storey buildings during earthquake events has become increasingly important in earthquake resistant design. For buildings on shallow foundations, SFSI refers to nonlinear geometric effects associated with uplift of the foundation from the supporting soil as well as nonlinear soil deformation effects. These effects can potentially be beneficial for structural performance, reducing forces transmitted from ground shaking to the structure. However, there is also the potential consequence of residual settlement and rotation of the foundation. This Thesis investigates the influence of SFSI in the performance of multi-storey buildings on shallow foundations through earthquake observations, experimental testing, and development of spring-bed numerical models that can be incorporated into integrated earthquake resistant design procedures. Observations were made following the 22 February 2011 Christchurch Earthquake in New Zealand of a number of multi-storey buildings on shallow foundations that performed satisfactorily. This was predominantly the case in areas where shallow foundations, typically large raft foundations, were founded on competent gravel and where there was no significant manifestation of liquefaction at the ground surface. The properties of these buildings and the soils they are founded on directed experimental work that was conducted to investigate the mechanisms by which SFSI may have influenced the behaviour of these types of structure-foundation systems. Centrifuge experiments were undertaken at the University of Dundee, Scotland using a range of structure-foundation models and a layer of dense cohesionless soil to simulate the situation in Christchurch where multi-storey buildings on shallow foundations performed well. Three equivalent single degree of freedom (SDOF) models representing 3, 5, and 7 storey buildings with identical large raft foundations were subjected to a range of dynamic Ricker wavelet excitations and Christchurch Earthquake records to investigate the influence of SFSI on the response of the equivalent buildings. The experimental results show that nonlinear SFSI has a significant influence on structural response and overall foundation deformations, even though the large raft foundations on competent soil meant that there was a significant reserve of bearing capacity available and nonlinear deformations may have been considered to have had minimal effect. Uplift of the foundation from the supporting soil was observed across a wide range of input motion amplitudes and was particularly significant as the amplitude of motion increased. Permanent soil deformation represented by foundation settlement and residual rotation was also observed but mainly for the larger input motions. However, the absolute extent of uplift and permanent soil deformation was very small compared to the size of the foundation meaning the serviceability of the building would still likely be maintained during large earthquake events. Even so, the small extent of SFSI resulted in attenuation of the response of the structure as the equivalent period of vibration was lengthened and the equivalent damping in the system increased. The experimental work undertaken was used to validate and enhance numerical modelling techniques that are simple yet sophisticated and promote interaction between geotechnical and structural specialists involved in the design of multi-storey buildings. Spring-bed modelling techniques were utilised as they provide a balance between ease of use, and thus ease of interaction with structural specialists who have these techniques readily available in practice, and theoretically rigorous solutions. Fixed base and elastic spring-bed models showed they were unable to capture the behaviour of the structure-foundation models tested in the centrifuge experiments. SFSI spring-bed models were able to more accurately capture the behaviour but recommendations were proposed for the parameters used to define the springs so that the numerical models closely matched experimental results. From the spring-bed modelling and results of centrifuge experiments, an equivalent linear design procedure was proposed along with a procedure and recommendations for the implementation of nonlinear SFSI spring-bed models in practice. The combination of earthquake observations, experimental testing, and simplified numerical analysis has shown how SFSI is influential in the earthquake performance of multi-storey buildings on shallow foundations and should be incorporated into earthquake resistant design of these structures.

Research Papers, Lincoln University

This thesis investigates landscape disturbance history in Westland since 1350 AD. Specifically, I test the hypothesis that large-magnitude regional episodes of natural disturbance have periodically devastated portions of the landscape and forest, and that these were caused by infrequent earthquakes along the Alpine Fault. Forest stand history reconstruction was used to determine the timing and extent of erosion and sedimentation events that initiated new forest cohorts in a 1412 ha study area in the Karangarua River catchment, south Westland. Over 85 % of the study area was disturbed sufficiently by erosion/sedimentation since 1350 AD to initiate new forest cohorts. During this time four episodes of catchment-wide disturbance impacted the study area, and these took place about 1825 AD ± 5 years (Ruera episode), 1715 AD ± 5 years (Sparkling episode), 1615 AD ± 5 years (McTaggart episode), and 1445 AD ± 15 years (Junction episode). The three most recent episodes disturbed 10 %, 35-40 % and 32-50 % respectively of the study area. The Junction episode disturbed at least 6 % of the study area, but elimination of evidence by more recent disturbances prevented an upper limit being defined. The three earliest episodes correspond to the date-ranges for three Alpine Fault earthquakes from geological data, and are the only episodes of disturbance within each date-range. An earthquake cause is also consistent with features of the disturbance record: large portions of the study area were disturbed, disturbance occurred on all types 'of landforms, and terrace surfaces were abandoned upstream of the Alpine Fault. On this basis erosion/sedimentation induced by Alpine Fault earthquakes has disturbed 14-20 % of the land surface in the study area per century. Storms and other non-seismic erosional processes have disturbed 3-4 % per century. To examine the importance of the Alpine Fault earthquakes to forest disturbance throughout Westland, I collated all available data on conifer stand age structures in the region and identified dates of disturbance events from 55 even-aged cohorts of trees. Three region-wide episodes of forest disturbance since 1350 AD were found in this sample, and these matched the three Alpine Fault earthquake-caused episodes found in the Karangarua. Forest disturbance at these times was widespread across Westland over at least 200 km from Paringa to Hokitika, and originated from both tree fall and erosion processes. This disturbance history can explain the long-observed regional conifer forest pattern in Westland, of a predominance of similar-sized stands of trees and a relative lack of small-sized (young) stands. The many similar-sized stands are a consequence of synchronous forest disturbance and re-establishment accompanying the infrequent Alpine Fault earthquakes, while the dominance of mature stands of trees and relative lack of young small-sized trees in stands is explained by the long lapsed time since the last Alpine Fault earthquake (c. 280 years). I applied the landscape disturbance history information to the existing geological data to reconstruct the paleoseismicity of the Alpine Fault since 1350 AD. Best estimates for the timing of the most recent three rupture events from these data are 1715 AD ± 5 years, 1615 AD ± 5 years and 1445 AD ± 15 years. Earthquake recurrence intervals were variable, ranging from about 100 years to at least 280 years (the lapsed time since the last event). All three events caused forest and geomorphic disturbance over at least a 200 km section of Fault between the Karangarua and Hokitika Rivers, and were probably single rupture events. Suppressions in cross dated tree-ring chronologies in the western South Island suggest that the last rupture occurred in 1717 AD, and extended as a single rupture from Haupiri to Fiordland, a distance along the Fault of 375 km.

Research Papers, Lincoln University

Saltwater Forest is a Dacrydium cupressinum-dominated lowland forest covering 9000 ha in south Westland, South Island, New Zealand. Four thousand hectares is managed for sustainable production of indigenous timber. The aim of this study was to provide an integrated analysis of soils, soil-landform relationships, and soil-vegetation relationships at broad and detailed scales. The broad scale understandings provide a framework in which existing or future studies can be placed and the detailed studies elucidate sources of soil and forest variability. Glacial landforms dominate. They include late Pleistocene lateral, terminal and ablation moraines, and outwash aggradation and degradation terraces. Deposits and landforms from six glacial advances have been recognised ranging from latest Last (Otira) Glaciation to Penultimate (Waimea) Glaciation. The absolute ages of landforms were established by analysis of the thickness and soil stratigraphy of loess coverbeds, augmented with radiocarbon dating and phytolith and pollen analysis. In the prevailing high rainfall of Westland soil formation is rapid. The rate of loess accretion in Saltwater Forest (ca. 30 mm ka⁻¹) has been low enough that soil formation and loess accretion took place contemporaneously. Soils formed in this manner are known as upbuilding soils. The significant difference between upbuilding pedogenesis and pedogenesis in a topdown sense into an existing sediment body is that each subsoil increment of an upbuilding soil has experienced processes of all horizons above. In Saltwater Forest subsoils of upbuilding soils are strongly altered because they have experienced the extremely acid environment of the soil surface at some earlier time. Some soil chronosequence studies in Westland have included upbuilding soils formed in loess as the older members of the sequence. Rates and types of processes inferred from these soils should be reviewed because upbuilding is a different pedogenic pathway to topdown pedogenesis. Landform age and morphology were used as a primary stratification for a study of the soil pattern and nature of soil variability in the 4000 ha production area of Saltwater Forest. The age of landforms (> 14 ka) and rapid soil formation mean that soils are uniformly strongly weathered and leached. Soils include Humic Organic Soils, Perch-gley Podzols, Acid Gley Soils, Allophanic Brown Soils, and Orthic or Pan Podzols. The major influence on the nature of soils is site hydrology which is determined by macroscale features of landforms (slope, relief, drainage density), mesoscale effects related to position on landforms, and microscale influences determined by microtopography and individual tree effects. Much of the soil variability arises at microscales so that it is not possible to map areas of uniform soils at practical map scales. The distribution of soil variability across spatial scales, in relation to the intensity of forest management, dictates that it is most appropriate to map soil complexes with boundaries coinciding with landforms. Disturbance of canopy trees is an important agent in forest dynamics. The frequency of forest disturbance in the production area of Saltwater Forest varies in a systematic way among landforms in accord with changes in abundance of different soils. The frequency of forest turnover is highest on landforms with the greatest abundance of extremely poorly-drained Organic Soils. As the abundance of better-drained soils increases the frequency of forest turnover declines. Changes in turnover frequency are reflected in the mean size and density of canopy trees (Dacrydium cupressinum) among landforms. Terrace and ablation moraine landforms with the greatest abundance of extremely poorly-drained soils have on average the smallest trees growing most densely. The steep lateral moraines, characterised by well drained soils, have fewer, larger trees. The changes manifested at the landform scale are an integration of processes operating over much shorter range as a result of short-range soil variability. The systematic changes in forest structure and turnover frequency among landforms and soils have important implications for sustainable forest management.

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.

Research papers, University of Canterbury Library

The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.

Research papers, University of Canterbury Library

A series of undrained cyclic direct simple shear (DSS) tests on specimens of sandy silty soils are used to evaluate the effects of fines content, fabric and layered structure on the liquefaction response of sandy soils containing non-plastic fines. Test soils originate from shallow deposits in Christchurch, New Zealand, where severe and damaging manifestations of liquefaction occurred during the 2010-2011 Canterbury earthquakes. A procedure for reconstituting specimens by water sedimentation is employed. This specimen preparation technique involves first pluviation of soil through a water column, and then application of gentle vibrations to the mould (tapping) to prepare specimens with different initial densities. This procedure is applied to prepare uniform specimens, and layered specimens with a silt layer atop a sand layer. Cyclic DSS tests are performed on water-sedimented specimens of two sands, a silt, and sand-silt mixtures with different fines contents. Through this testing program, effects of density, time of vibration during preparation, fines content, and layered structure on cyclic behaviour and liquefaction resistance are investigated. Additional information necessary to characterise soil behaviour is provided by particle size distribution analyses, index void ratio testing, and Scanning Electronic Microscope imaging. The results of cyclic DSS tests show that, for all tested soils, specimens vibrated for longer period of time have lower void ratios, higher relative density, and greater liquefaction resistance. One of the tested sands undergoes significant increase in relative density and liquefaction resistance following prolonged vibration. The other sand exhibits lower increase in relative density and in liquefaction resistance when vibrated for the same period of time. Liquefaction resistance of sand-silt mixtures prepared using this latter sand shows a correlation with relative density irrespective of fines content. In general, however, magnitudes of changes in liquefaction resistance for given variations in vibration time, relative density, or void ratio vary depending on soils under consideration. Characterization based on maximum and minimum void ratios indicates that tested soils develop different structures as fines are added to their respective host sands. These structures influence initial specimen density, strains during consolidation, cyclic liquefaction resistance, and undrained cyclic response of each soil. The different structures are the outcome of differences in particle size distributions, average particle sizes, and particle shapes of the two host sands and of the different relationships between these properties and those of the silt. Fines content alone does not provide an effective characterization of the effects of these factors. Monotonic DSS tests are also performed on specimens prepared by water sedimentation, and on specimens prepared by moist tamping, to identify the critical state lines of tested soils. These critical state lines provide the basis for an alternative interpretation of cyclic DSS tests results within the critical state framework. It is shown that test results imply general consistency between observed cyclic and monotonic DSS soil response. The effects of specimen layering are scrutinised by comparing DSS test results for uniform and layered specimens of the same soils. In this case, only a limited number of tests is performed, and the range of densities considered for the layered specimens is also limited. Caution is therefore required in interpretation of their results. The liquefaction resistance of layered specimens appears to be influenced by the bottom sand layer, irrespective of the global fines content of the specimen. The presence of a layered structure does not result in significant differences in terms of liquefaction response with respect to uniform sand specimens. Cyclic triaxial data for Christchurch sandy silty soils available from previous studies are used to comparatively examine the behaviour observed in the tests of this study. The cyclic DSS liquefaction resistance of water-sedimented specimens is consistent with cyclic triaxial tests on undisturbed specimens performed by other investigators. The two data sets result in similar liquefaction triggering relationships for these soils. However, stress-strain response characteristics for the two types of specimens are different, and undisturbed triaxial specimen exhibit a slower rate of increase in shear strains compared to water-sedimented DSS specimens. This could be due to the greater influence of fabric of the undisturbed specimens.

Research Papers, Lincoln University

Globally, the maximum elevations at which treelines are observed to occur coincide with a 6.4 °C soil isotherm. However, when observed at finer scales, treelines display a considerable degree of spatial complexity in their patterns across the landscape and are often found occurring at lower elevations than expected relative to the global-scale pattern. There is still a lack of understanding of how the abiotic environment imposes constraints on treeline patterns, the scales at which different effects are acting, and how these effects vary over large spatial extents. In this thesis, I examined abrupt Nothofagus treelines across seven degrees of latitude in New Zealand in order to investigate two broad questions: (1) What is the nature and extent of spatial variability in Nothofagus treelines across the country? (2) How is this variation associated with abiotic variation at different spatial scales? A range of GIS, statistical, and atmospheric modelling methods were applied to address these two questions. First, I characterised Nothofagus treeline patterns at a 15x15km scale across New Zealand using a set of seven, GIS-derived, quantitative metrics that describe different aspects of treeline position, shape, spatial configuration, and relationships with adjacent vegetation. Multivariate clustering of these metrics revealed distinct treeline types that showed strong spatial aggregation across the country. This suggests a strong spatial structuring of the abiotic environment which, in turn, drives treeline patterns. About half of the multivariate treeline metric variation was explained by patterns of climate, substrate, topographic and disturbance variability; on the whole, climatic and disturbance factors were most influential. Second, I developed a conceptual model that describes how treeline elevation may vary at different scales according to three categories of effects: thermal modifying effects, physiological stressors, and disturbance effects. I tested the relevance of this model for Nothofagus treelines by investigating treeline elevation variation at five nested scales (regional to local) using a hierarchical design based on nested river catchments. Hierarchical linear modelling revealed that the majority of the variation in treeline elevation resided at the broadest, regional scale, which was best explained by the thermal modifying effects of solar radiation, mountain mass, and differences in the potential for cold air ponding. Nonetheless, at finer scales, physiological and disturbance effects were important and acted to modify the regional trend at these scales. These results suggest that variation in abrupt treeline elevations are due to both broad-scale temperature-based growth limitation processes and finer-scale stress- and disturbance-related effects on seedling establishment. Third, I explored the applicability of a meso-scale atmospheric model, The Air Pollution Model (TAPM), for generating 200 m resolution, hourly topoclimatic data for temperature, incoming and outgoing radiation, relative humidity, and wind speeds. Initial assessments of TAPM outputs against data from two climate station locations over seven years showed that the model could generate predictions with a consistent level of accuracy for both sites, and which agreed with other evaluations in the literature. TAPM was then used to generate data at 28, 7x7 km Nothofagus treeline zones across New Zealand for January (summer) and July (winter) 2002. Using mixed-effects linear models, I determined that both site-level factors (mean growing season temperature, mountain mass, precipitation, earthquake intensity) and local-level landform (slope and convexity) and topoclimatic factors (solar radiation, photoinhibition index, frost index, desiccation index) were influential in explaining variation in treeline elevation within and among these sites. Treelines were generally closer to their site-level maxima in regions with higher mean growing season temperatures, larger mountains, and lower levels of precipitation. Within sites, higher treelines were associated with higher solar radiation, and lower photoinhibition and desiccation index values, in January, and lower desiccation index values in July. Higher treelines were also significantly associated with steeper, more convex landforms. Overall, this thesis shows that investigating treelines across extensive areas at multiple study scales enables the development of a more comprehensive understanding of treeline variability and underlying environmental constraints. These results can be used to formulate new hypotheses regarding the mechanisms driving treeline formation and to guide the optimal choice of field sites at which to test these hypotheses.

Research papers, University of Canterbury Library

This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.

Research papers, University of Canterbury Library

This thesis is concerned with the effects of lateral confining reinforcement on the ductile behaviour of reinforced concrete columns. The contents of the chapters are summarized as follows. In Chapter one, the general problems in seismic design are discussed and earthquake design methods based on the ductile design approach are described. Japanese, New Zealand and United States design codes are compared. Finally, the scope of this research project is outlined. In Chapter two, after reviewing previous research on confined concrete, the factors which affect the effectiveness of lateral confinement are discussed. Especially the effects of the yield strength of transverse reinforcement, the compressive strength of plain concrete and the strain gradient in the column section due to bending are discussed based on tests which were conducted by the author et al at Kyoto University and Akashi Technological College, Japan. In the axial compression tests on spirally reinforced concrete cylinders (150 mm in diameter by 300 mm in height), the yield strength of transverse reinforcement and the compressive strength of plain concrete were varied from 161 MPa to 1352 MPa and from 17 MPa to 60 MPa, respectively, as experimental parameters. It is found that, when high strength spirals are used as confining reinforcement, the strength and ductility of the confined core concrete are remarkably enhanced but need to be estimated assuming several failure modes which could occur. These are based on the observations that concrete cylinders with high strength spirals suddenly failed at a concrete compressive strain of 2 to 3.5 % due to explosive crushing of the core concrete between the spiral bars or due to bearing failure of the core concrete immediately beneath the spiral bars, while the concrete cylinders with ordinary strength spirals failed in a gentle manner normally observed. In addition, eccentric loading tests were conducted on concrete columns with 200 mm square section confined by square spirals. It is found that the effectiveness of confining reinforcement is reduced by the presence of the strain gradient along the transverse section of column. In Chapter three, the effectiveness of transverse reinforcement with various types of anchorage details which simplify the fabrication of reinforcing cages are investigated. Eight reinforced concrete columns, with either 400 mm or 550 mm square cross sections, were tested subjected to axial compression loading and cyclic lateral loading which simulated a severe earthquake. The transverse reinforcement consisted of arrangements of square perimeter hoops with 135° end hooks, cross ties with 90° and 135° or 180° end hooks, and 'U' and 'J' shaped cross ties and perimeter hoops with tension splices. Conclusions are reached with regard to the effectiveness of the tested anchorage details in the plastic hinge regions of columns designed for earthquake resistance. In Chapter four, the effectiveness of interlocking spirals as transverse reinforcement is studied. Firstly, the general aspects and the related problems of interlocking spirals to provide adequate ductility in the potential plastic hinge region of columns are discussed, referring to the provisions in the New Zealand code,the CALTRANS (California Transportation Authority) code and other related codes. Secondly, based on those discussions, a design method to securely interlock the spirals is proposed. Thirdly, the effectiveness of interlocking spirals is assessed based on column tests conducted as part of this study. Three columns with interlocking spirals and, for comparison, one rectangular column with rectangular hoopsandcross ties, were tested under cyclic horizontal loading which simulated a severe earthquake. The sections of those columns were 400 mm by 600 mm. In Chapter five, analytical models to investigate the buckling behaviour of longitudinal reinforcement restrained by cross ties with 90° and 135° end hooks and by peripheral hoops are proposed. The analyzed results using the proposed models compare well with the experimental observations described in Chapter three. Using those proposed models, a method to check the effectiveness of cross ties with 90° and 135° end hooks is proposed for practical design purposes. In Chapter six, a theory for the prediction of the ultimate longitudinal compressive concrete strain at the stage of first hoop fracture referred to as the "Energy Balance Theory", which has been developed by Mander, Priestley and Park at University of Canterbury, is introduced. After discussing the problems in the "Energy Balance Theory", a modified theory for the prediction of the ultimate longitudinal compressive concrete strain at the stage of first hoop fracture is proposed. The predictions from the modified theory are found to compare well with previous experimental results.

Research papers, University of Canterbury Library

Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.

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

A review of the literature showed the lack of a truly effective damage avoidance solution for timber or hybrid timber moment resisting frames (MRFs). Full system damage avoidance selfcentring behaviour is difficult to achieve with existing systems due to damage to the floor slab caused by beam-elongation. A novel gravity rocking, self-centring beam-column joint with inherent and supplemental friction energy dissipation is proposed for low-medium rise buildings in all seismic zones where earthquake actions are greater than wind. Steel columns and timber beams are used in the hybrid MRF such that both the beam and column are continuous thus avoiding beam-elongation altogether. Corbels on the columns support the beams and generate resistance and self-centring through rocking under the influence of gravity. Supplemental friction sliders at the top of the beams resist sliding of the floor whilst dissipating energy as the floor lifts on the corbels and returns. 1:20 scale tests of 3-storey one-by-two bay building based on an earlier iteration of the proposed concept served as proof-of-concept and highlighted areas for improvement. A 1:5 scale 3-storey one-by-one bay building was subsequently designed. Sub-assembly tests of the beam-top asymmetric friction sliders demonstrated repeatable hysteresis. Quasi-static tests of the full building demonstrated a ‘flat bottomed’ flag-shaped hysteresis. Shake table tests to a suite of seven earthquakes scaled for Wellington with site soil type D to the serviceability limit state (SLS), ultimate limit state (ULS) and maximum credible event (MCE) intensity corresponding to an average return period of 25, 500 and 2500 years respectively were conducted. Additional earthquake records from the 22 February 2011 Christchurch earthquakes we included. A peak drift of 0.6%, 2.5% and 3.8% was reached for the worst SLS, ULS and MCE earthquake respectively whereas a peak drift of 4.5% was reached for the worst Christchurch record for tests in the plane of the MRF. Bi-directional tests were also conducted with the building oriented at 45 degrees on the shake table and the excitation factored by 1.41 to maintain the component in the direction of the MRF. Shear walls with friction slider hold-downs which reached similar drifts to the MRF were provided in the orthogonal direction. Similar peak drifts were reached by the MRF in the bi-directional tests, when the excitation was amplified as intended. The building self-centred with a maximum residual drift of 0.06% in the dynamic tests and demonstrated no significant damage. The member actions were magnified by up to 100% due to impact upon return of the floor after uplift when the peak drift reached 4.5%. Nonetheless, all of the members and connections remained essentially linearelastic. The shake table was able to produce a limited peak velocity of 0.275 m/s and this limited the severity of several of the ULS, MCE and Christchurch earthquakes, especially the near-field records with a large velocity pulse. The full earthquakes with uncapped velocity were simulated in a numerical model developed in SAP2000. The corbel supports were modelled with the friction isolator link element and the top sliders were modelled with a multi-linear plastic link element in parallel with a friction spring damper. The friction spring damper simulated the increase in resistance with increasing joint rotation and a near zero return stiffness, as exhibited by the 1:5 scale test building. A good match was achieved between the test quasi-static global force-displacement response and the numerical model, except a less flat unloading curve in the numerical model. The peak drift from the shake table tests also matched well. Simulations were also run for the full velocity earthquakes, including vertical ground acceleration and different floor imposed load scenarios. Excessive drift was predicted by the numerical model for the full velocity near-field earthquakes at the MCE intensity and a rubber stiffener for increasing the post joint-opening stiffness was found to limit the drift to 4.8%. Vertical ground acceleration had little effect on the global response. The system generates most of its lateral resistance from the floor weight, therefore increasing the floor imposed load increased the peak drift, but less than it would if the resistance of the system did not increase due to the additional floor load. A seismic design procedure was discussed under the framework of the existing direct displacement-based design method. An expression for calculating the area-based equivalent viscous damping (EVD) was derived and a conservative correction factor of 0.8 was suggested. A high EVD of up to about 15% can be achieved with the proposed system at high displacement ductility levels if the resistance of the top friction sliders is maximised without compromising reliable return of the floor after uplift. Uniform strength joints with an equal corbel length up the height of the building and similar inter-storey drifts result in minimal relative inter-floor uplift, except between the first floor and ground. Guidelines for detailing the joint for damage avoidance including bi-directional movement were also developed.

Research papers, University of Canterbury Library

Geosynthetic reinforced soil (GRS) walls involve the use of geosynthetic reinforcement (polymer material) within the retained backfill, forming a reinforced soil block where transmission of overturning and sliding forces on the wall to the backfill occurs. Key advantages of GRS systems include the reduced need for large foundations, cost reduction (up to 50%), lower environmental costs, faster construction and significantly improved seismic performance as observed in previous earthquakes. Design methods in New Zealand have not been well established and as a result, GRS structures do not have a uniform level of seismic and static resistance; hence involve different risks of failure. Further research is required to better understand the seismic behaviour of GRS structures to advance design practices. The experimental study of this research involved a series of twelve 1-g shake table tests on reduced-scale (1:5) GRS wall models using the University of Canterbury shake-table. The seismic excitation of the models was unidirectional sinusoidal input motion with a predominant frequency of 5Hz and 10s duration. Seismic excitation of the model commenced at an acceleration amplitude level of 0.1g and was incrementally increased by 0.1g in subsequent excitation levels up to failure (excessive displacement of the wall panel). The wall models were 900mm high with a full-height rigid facing panel and five layers of Microgird reinforcement (reinforcement spacing of 150mm). The wall panel toe was founded on a rigid foundation and was free to slide. The backfill deposit was constructed from dry Albany sand to a backfill relative density, Dr = 85% or 50% through model vibration. The influence of GRS wall parameters such as reinforcement length and layout, backfill density and application of a 3kPa surcharge on the backfill surface was investigated in the testing sequence. Through extensive instrumentation of the wall models, the wall facing displacements, backfill accelerations, earth pressures and reinforcement loads were recorded at the varying levels of model excitation. Additionally, backfill deformation was also measured through high-speed imaging and Geotechnical Particle Image Velocimetry (GeoPIV) analysis. The GeoPIV analysis enabled the identification of the evolution of shear strains and volumetric strains within the backfill at low strain levels before failure of the wall thus allowing interpretations to be made regarding the strain development and shear band progression within the retained backfill. Rotation about the wall toe was the predominant failure mechanism in all excitation level with sliding only significant in the last two excitation levels, resulting in a bi-linear displacement acceleration curve. An increase in acceleration amplification with increasing excitation was observed with amplification factors of up to 1.5 recorded. Maximum seismic and static horizontal earth pressures were recorded at failure and were recorded at the wall toe. The highest reinforcement load was recorded at the lowest (deepest in the backfill) reinforcement layer with a decrease in peak load observed at failure, possibly due to pullout failure of the reinforcement layer. Conversely, peak reinforcement load was recorded at failure for the top reinforcement layer. The staggered reinforcement models exhibited greater wall stability than the uniform reinforcement models of L/H=0.75. However, similar critical accelerations were determined for the two wall models due to the coarseness of excitation level increments of 0.1g. The extended top reinforcements were found to restrict the rotational component of displacement and prevented the development of a preliminary shear band at the middle reinforcement layer, contributing positively to wall stability. Lower acceleration amplification factors were determined for the longer uniform reinforcement length models due to reduced model deformation. A greater distribution of reinforcement load towards the top two extended reinforcement layers was also observed in the staggered wall models. An increase in model backfill density was observed to result in greater wall stability than an increase in uniform reinforcement length. Greater acceleration amplification was observed in looser backfill models due to their lower model stiffness. Due to greater confinement of the reinforcement layers, greater reinforcement loads were developed in higher density wall models with less wall movement required to engage the reinforcement layers and mobilise their resistance. The application of surcharge on the backfill was observed to initially increase the wall stability due to greater normal stresses within the backfill but at greater excitation levels, the surcharge contribution to wall destabilising inertial forces outweighs its contribution to wall stability. As a result, no clear influence of surcharge on the critical acceleration of the wall models was observed. Lower acceleration amplification factors were observed for the surcharged models as the surcharge acts as a damper during excitation. The application of the surcharge also increases the magnitude of reinforcement load developed due to greater confinement and increased wall destabilising forces. The rotation of the wall panel resulted in the progressive development of shears surface with depth that extended from the backfill surface to the ends of the reinforcement (edge of the reinforced soil block). The resultant failure plane would have extended from the backfill surface to the lowest reinforcement layer before developing at the toe of the wall, forming a two-wedge failure mechanism. This is confirmed by development of failure planes at the lowest reinforcement layer (deepest with the backfill) and at the wall toe observed at the critical acceleration level. Key observations of the effect of different wall parameters from the GeoPIV results are found to be in good agreement with conclusions developed from the other forms of instrumentation. Further research is required to achieve the goal of developing seismic guidelines for GRS walls in geotechnical structures in New Zealand. This includes developing and testing wall models with a different facing type (segmental or wrap-around facing), load cell instrumentation of all reinforcement layers, dynamic loading on the wall panel and the use of local soils as the backfill material. Lastly, the limitations of the experimental procedure and wall models should be understood.

Research papers, University of Canterbury Library

Seismically vulnerable buildings constitute a major problem for the safety of human beings. In many parts of the world, reinforced concrete (RC) frame buildings designed and constructed with substandard detailing, no consideration of capacity design principles, and improper or no inclusion of the seismic actions, have been identified. Amongst those vulnerable building, one particular typology representative of the construction practice of the years previous to the 1970’s, that most likely represents the worst case scenario, has been widely investigated in the past. The deficiencies of that building typology are related to non-ductile detailing in beam column joints such as the use of plain round bars, the lack of stirrups inside the joint around the longitudinal reinforcement of the column, the use of 180° end hooks in the beams, the use of lap splices in potential ‘plastic hinge’ regions, and substandard quality of the materials. That type of detailing and the lack of a capacity design philosophy create a very fragile fuse in the structure where brittle inelastic behaviour is expected to occur, which is the panel zone region of exterior beam column joints. The non-ductile typology described above was extensively investigated at the University of Canterbury in the context of the project ‘Retrofit Solutions for New Zealand Multi-Storey Buildings’ (2004-2011), founded by the ‘Foundation for Research, Science and Technology’ Tūāpapa Rangahau Pūtaiao. The experimental campaign prior to the research carried out by the author consisted of quasi-static tests of beam column joint subassemblies subjected to lateral loading regime, with constant and varying axial load in the column. Most of those specimens were representative of a plane 2D frame (knee joint), while others represented a portion of a space 3D frame (corner joints), and only few of them had a floor slab, transverse beams, and lap splices. Using those experiments, several feasible, cost-effective, and non-invasive retrofit techniques were developed, improved, and refined. Nevertheless, the slow motion nature of those experiments did not take into account the dynamical component inherent to earthquake related problems. Amongst the set of techniques investigated, the use of FRP layers for strengthening beam column joints is of particular interest due to its versatility and the momentum that its use has gained in the current state of the practice. That particular retrofit technique was previously used to develop a strengthening scheme suitable for plane 2D and space 3D corner beam column joints, but lacking of floor slabs. In addition, a similar scheme was not developed for exterior joints of internal frames, referred here as ‘cruciform’. In this research a 2/5 scale RC frame model building comprising of two frames in parallel (external and internal) joined together by means of floor slabs and transverse beams, with non-ductile characteristics identical to those of the specimens investigated previously by others, and also including lap splices, was developed. In order to investigate the dynamic response of that building, a series of shake table tests with different ground motions were performed. After the first series of tests, the specimen was modified by connecting the spliced reinforcement in the columns in order to capture a different failure mode. Ground motions recorded during seismic events that occurred during the initial period of the experimental campaign (2010) were used in the subsequent experiments. The hierarchy of strengths and sequence of events in the panel zone region were evaluated in an extended version of the bending moment-axial load (M-N) performance domain developed by others. That extension was required due to the asymmetry in the beam cross section introduced by the floor slab. In addition, the effect of the torsion resistance provided by the spandrel (transverse beam) was included. In order to upgrade the brittle and unstable performance of the as-built/repaired specimen, a practical and suitable ad-hoc FRP retrofit intervention was developed, following a partial retrofit strategy that aimed to strengthen exterior beam column joints only (corner and cruciform). The ability of the new FRP scheme to revert the sequence of events in the panel zone region was evaluated using the extended version of the M-N performance domain as well as the guidelines for strengthening plane joints developed by others. Weakening of the floor slab in a novel configuration was also incorporated with the purpose of reducing the flexural capacity of the beam under negative bending moment (slab in tension), enabling the damage relocation from the joint into the beam. The efficacy of the developed retrofit intervention in upgrading the seismic performance of the as-built specimen was investigated using shake table tests with the input motions used in the experiments of the as-built/repaired specimen. Numerical work aimed to predict the response of the model building during the most relevant shake table tests was carried out. By using a simple numerical model with concentrated plasticity elements constructed in Ruaumoko2D, the results of blind and post-experimental predictions of the response of the specimen were addressed. Differences in the predicted response of the building using the nominal and the actual recorded motions of the shake table were investigated. The dependence of the accuracy of the numerical predictions on the assumed values of the parameters that control the hysteresis rules of key structural members was reviewed. During the execution of the experimental campaign part of this thesis, two major earthquakes affected the central part of Chile (27 of February 2010 Maule earthquake) and the Canterbury region in New Zealand (22 February 2011 Canterbury earthquake), respectively. As the author had the opportunity to experience those events and investigate their consequences in structures, the observations related to non-ductile detailing and drawbacks in the state of the practice related to reinforced concrete walls was also addressed in this research, resulting in preliminary recommendations for the refinement of current seismic code provisions and assessment guidelines. The investigations of the ground motions recorded during those and other earthquakes were used to review the procedures related to the input motions used for nonlinear dynamic analysis of buildings as required by most of the current code provisions. Inelastic displacement spectra were constructed using ground motions recorded during the earthquakes mentioned above, in order to investigate the adequacy of modification factors used to obtain reduced design spectra from elastic counterparts. Finally a simplified assessment procedure for RC walls that incorporates capacity compatible spectral demands is proposed.

Research papers, University of Canterbury Library

The assessment of damage and remaining capacity after an earthquake is an immediate measure to determine whether a reinforced concrete (RC) building is usable and safe for occupants. The recent Christchurch earthquake (22 February 2011) caused a uniquely severe level of structural damage to modern buildings, resulting in extensive damage to the building stock. About 60% of damaged multistorey concrete buildings (3 storeys and up) were demolished after the earthquake, and the cost of reconstruction amounted to 40 billion NZD. The aftermath disclosed issues of great complexities regarding the future of the RC buildings damaged by the earthquakes. This highlighted the importance of post-event decision-making, as the outcome will allow the appropriate course of action—demolition, repair or acceptance of the existing building—to be considered. To adopt the proper strategy, accurate assessment of the residual capacity and the level of damage is required. This doctoral dissertation aims to assess the damage and remaining capacity at constituent material and member level (i.e., concrete material and beams) through a systematic approach in an attempt to address part of an existing gap in the available literature. Since the residual capacity of RC members is not unique and depends on previously applied loading history, post-event residual capacity in this study was assessed in terms of fraction of fatigue life (i.e., the number of cycles required to failure). This research comprises three main parts: (1) residual capacity and damage assessment at material level (i.e., concrete), (2) post-yield bond deterioration and damage assessment at the interface of steel and concrete, and, finally, (3) residual capacity and damage assessment at member level (i.e., RC beam). The first part of this research focused on damage assessment and the remaining capacity of concrete from a material point of view. It aimed to employ appropriate and reliable durability-based testing and image-detection techniques to quantify deterioration in the mechanical properties of concrete on the basis that stress-induced damage occurred in the microstructural system of the concrete material. To this end, in the first phase, a feasibility study was conducted in which a combination of oxygen permeability, electrical resistivity and porosity tests were assessed to determine if they were robust and reliable enough to reveal damage which occurred in the microstructural system of concrete. The results, in terms of change in permeability, electrical resistivity and porosity features of disk samples taken from the middle third of damaged concrete cylinders (200 mm × 100 mm) monotonically pre-loaded to 50%, 70%, 90% and 95% of the ultimate strength (f′c), showed the permeability test is a reliable tool to identify the degree of damage, due to its high sensitivity to the load-induced microcracking. In parallel, to determine the residual capacity, the companion damaged concrete cylinders already loaded to the same level of compressive strength were reloaded up to failure. Comparing the stress–strain relationship of damaged concrete with intact material, it was also found that the strain capacity of the reloaded pre-damaged concrete cylinders decreases while strength remained virtually unchanged. In the second phase of the first part, a fluorescent microscopy technique was used to assess the damage and develop a correlation between material degradation, by virtue of the geometrical features, and damage to the concrete. To account for the effect of confinement and cyclic loading, in the third phase, the residual capacity and damage assessment of unconfined and GFRP confined concrete cylinders subjected to low-cycle fatigue loading, was investigated. Similar to the first phase, permeability testing technique was used to provide an indirect evaluation of fatigue damage. Finally, in the fourth phase of the first part, the suitability of permeability testing technique to assess damage was evaluated for cored concrete taken from three types of RC members: columns, beams and a beam-column joint. In view of the fact that the composite action of an RC member is highly dependent on the bond between reinforcement and surrounding concrete, understanding the deterioration of the bond in the post-yield range of strain in steel was crucial to assess damage at member level. Therefore, in the second phase of this research, a state-of-the- art distributed fibre optic strain sensor system (DFOSSS) system was used to evaluate bond deterioration in a cantilever RC beam subjected to monotonic lateral loading. The technology allowed the continuous capture of strain, every 2.6 mm along the length, in both reinforcing bars and cover concrete. The strain profile provided a basis by which the slip, axial stress and bond stress distributions were then established. In the third part, the study focused on the damage assessment and residual capacity of seven half-scale RC beams subjected to a constant-amplitude cyclic loading protocol. In the first stage, the structural performances of three specimens under constant-amplitude fatigue at 1%, 2% and 4% chord rotation (drift) were examined. In parallel, the number of cycles to failure, degradation in strength, stiffness and energy dissipation were characterized. In the second stage, four RC beams were subjected to loading up to 70% and 90% of their fatigue life, at 2% and 4% drift, and then monotonically pulled up to failure. To determine the residual flexural capacity, the lateral force–displacement results of pre-damaged specimens were compared with an undamaged specimen subjected to only monotonic loading. The study showed significant losses in strength, deformability, stiffness and energy dissipation capacity. A nonlinear finite element analysis (FEA) using concrete damage plasticity (CDP) model was also conducted in ABAQUS to numerically investigate the behaviour of the tested specimen. The results of the FE simulations indicated a reasonable response compared with the behaviour of the test specimen in terms of force–displacement and cracking pattern. During the Christchurch earthquake it was observed that the loading history has a significant influence on structural responses. While in conventional pseudo-static loading protocol, internal forces can be redistributed along the plastic length: there is little chance for structures undergoing high initial loading amplitude to redistribute pertinent stresses. As a result, in the third phase of this part, the effect of high rate of loading on the behaviour of seismically designed RC beams was investigated. Two half-scale cantilever RC beams were subjected to similar constant-amplitude cyclic loading at 2% and 4% drifts, but at a rate of 500 mm/s. Due to the incapability of conventional measuring techniques, a motion-tracking system was employed for data acquisition with the high-speed tests. The effect of rate of loading on the fatigue life of specimens (i.e., the number of cycles required to failure), secant stiffness, failure mode, cracking pattern, beam elongations and bar fracture surface were analysed. Integrating the results of all parts of this research has resulted in a better understanding of residual capacity and the development of damage at both the material and member level by using a low-cycle fatigue approach.

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

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