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.
The Leader Fault was one of at least 17 faults that ruptured the ground surface across the northeastern South Island of New Zealand during the Mw 7.8 2016 Kaikōura Earthquake. The southern ~6 km of the Leader Fault, here referred to as the South Leader Fault (SLF), ruptured the North Canterbury (tectonic) Domain and is the primary focus of this study. The main objective of the thesis is to understand the key factors that contributed to the geometry and kinematics of the 2016 SLF rupture and its intersection with The Humps Fault (HF). This thesis employs a combination of techniques to achieve the primary objective, including detailed mapping of the bedrock geology, geomorphology and 2016 rupture, measurement of 2016 ground surface displacements, kinematic analysis of slip vectors from the earthquake, and logging of a single natural exposure across a 2016 rupture that was treated as a paleoseismic trench. The resulting datasets were collected in the field, from terrestrial LiDAR and InSAR imagery, and from historical (pre-earthquake) aerial photographs for a ~11 km2 study area. Surface ruptures in the study area are a miniature version of the entire rupture from the earthquake; they are geometrically and kinematically complex, with many individual and discontinuous segments of varying orientations and slip senses which are distributed across a zone up to ~3.5 km wide. Despite this variability, three main groups of ruptures have been identified. These are: 1) NE-SW striking, shallow to moderate dipping (25-45°W) faults that are approximately parallel to Cenozoic bedding with mainly reverse dip-slip and, and for the purposes of this thesis, are considered to be part of the SLF. 2) N-S striking, steeply dipping (~85°E) oblique sinistral faults that are up to the west and part of the SLF. 3) E-NE striking, moderate to steeply dipping (45-68°N) dextral reverse faults which are part of the HF. Bedding-parallel faults are interpreted to be flexural slip structures formed during folding of the near-surface Cenozoic strata, while the steeply dipping SLF ruptured a pre-existing bedrock fault which has little topographic expression. Groups 1 and 2 faults were both locally used for gravitational failure during the earthquake. Despite this non-tectonic fault movement, the slip vectors for faults that ruptured during the earthquake are broadly consistent with NCD tectonics and the regional ~100-120° trend of the principal horizontal stress/strain axes. Previous earthquake activity on the SLF is required by its displacement of Cenozoic formations but Late Quaternary slip on the fault prior to 2016 is neither supported by pre-existing fault scarps nor by changes in topography across the fault. By contrast, at least two earthquakes (including 2016) appear to have ruptured the HF from the mid Holocene, consistent with recurrence intervals of no more than ~7 kyr, and with preliminary observations from trenches on the fault farther to the west. The disparity in paleoearthquake records of the two faults suggests that they typically do not rupture together, thus it is concluded that the HF-SLF rupture pattern observed in the Kaikōura Earthquake rarely occurs in a single earthquake.
To this extent, modern buildings generally demonstrated good resistance to collapse during the recent earthquakes in New Zealand. However, damage to non-structural elements (NSE) has been persistent during these events. NSEs include secondary systems or components attached to the floors, roofs, and walls of a building or industrial facility that are not explicitly designed to participate in the main vertical or lateral load-bearing mechanism of the structure. They play a major role in the operational and functional aspects of buildings and contribute a major portion of the building’s overall cost. Therefore, they are expected to accommodate the effects of seismic actions such as drifts and accelerations. Typical examples of NSEs include internal non-loadbearing partitions, suspended ceilings, sprinkler piping systems, architectural claddings, building contents, mechanical/electrical equipment, and furnishings. The main focus of this thesis is the drift sensitive NSEs: precast concrete cladding panels and internal partition walls. Even though most precast concrete cladding panels performed well from a life-safety point of view during recent earthquakes in NZ, some collapsed panels posed a significant threat to life safety. It is, therefore, important that the design and detailing of the panel-to-structure connections ensure that their strength and displacement capacity are adequate to meet the corresponding seismic demands, at least during design level earthquakes. In contrast, the partition wall is likely to get damaged and lose serviceability at a low inter-story drift unless designed to accommodate the relative deformations between them and the structure. Partition walls suffered wide-ranging damage such as screw failures, diagonal cracking, detachments to the gypsum linings, and anchorage failures during the 2011 Canterbury Earthquake Sequence in NZ. Therefore, the thesis is divided into two parts. Part I of the thesis focuses on developing novel low-damage precast concrete cladding panel connections, i.e. “rocking” connection details comprising vertically slotted steel embeds and weld plates. The low-damage seismic performance of novel “rocking” connection details is verified through experimental tests comprising uni-directional, bi-directional, and multi-storey scaled quasi-static cyclic tests. Comparison with the seismic performance of traditional panel connections reported in the literature demonstrated the system’s significantly improved seismic resilience. Furthermore, the finite element models of panel connections and sealants are developed in ABAQUS. The force-drift responses of the “rocking” panel system modelled in SAP2000 is compared with the experimental results to evaluate their accuracy and validity. Part II of the thesis focuses on a) understanding the seismic performance of traditional rigid timber-framed partition wall, b) development and verification of low-damage connections (i.e. “rocking” connection details comprising of dual-slot tracks), and c) seismic evaluation of partition walls with a novel “bracketed and slotted” connections (comprising of innovative fastener and plastic bracket named Flexibracket) under uni-directional and bidirectional quasi-static cyclic loadings. Moreover, parametric investigation of the partition walls was conducted through several experimental tests to understand better the pros and cons of the rocking connection details. The experimental results have confirmed that the implementation of the proposed low damage solutions of precast cladding panels and internal partition walls can significantly reduce their damage in a building.
According to TS 1170.5, designing a building to satisfy code-prescribed criteria (e.g., drift limit, member safety, P-Δ stability) at the ultimate limit state and relying on the inherent margins within the design code would lead to an acceptable mean annual frequency of collapse (λ꜀) in the range of 10−⁴ to 10−⁵. Modern performance objectives, such as λ꜀ and expected annual loss (EAL), are not explicitly considered. Although buckling-restrained braced frame (BRBF) buildings were widely adopted as lateral load-resisting systems for office and car park buildings in the Christchurch rebuild following the Canterbury earthquakes in New Zealand, there are currently no official guidelines for their design. The primary focus of this study is to develop a risk-targeted design framework for BRBF buildings that can achieve the performance objectives desired by stakeholders. To this extent, key factors influencing λ꜀ and EAL of BRBF buildings are identified. These factors include gusset plate design, number of storeys, design drift limit, BRBF beam-column connection, brace configuration, brace angle, brace material grade, and analysis method (equivalent lateral force vs. modal response spectrum). A novel 3D BRBF modelling approach capable of simulating out-of-plane buckling failure of buckling-restrained brace (BRB) gusset plates is developed. Prior experimental studies on sub-assemblies conducted elsewhere have demonstrated that gusset plates and end zones may buckle out of plane prematurely, before BRBs reach their maximum axial compression load carrying capacity. Current 2D BRBF macro models, typically used in research, cannot simulate this failure mode. A conventional 2D BRBF model underestimates the λ꜀ of a case-study 4-storey super-X configured steel BRBF building (designed according to NZS-3404) by a factor of two compared to the estimate from the proposed 3D model. These findings suggest that the current NZS-3404 gusset plate design method may undersize gusset plates and that using a 2D BRBF model in this case can significantly underestimate λ꜀. Three improved alternative gusset plate design methods that are easy to implement in practice are identified from the literature. Gusset plates in two case-study 4-storey steel BRBF buildings with super-X and diagonal configurations are designed using both the NZS-3404 method and alternative methods. All three alternative design methods are found to be conservative, resulting in an almost three-fold lower λ꜀ for both case-study BRBF buildings compared to those designed using the NZS-3404 method. Analysis results indicate that (i) bidirectional interaction has no significant effect on gusset plate buckling and (ii) mid-span gusset plates are more susceptible to buckling than corner gusset plates. A framework for seismic loss assessment using incremental dynamic analysis (IDA), called loss-oriented hazard-consistent incremental dynamic analysis (LOHC-IDA), is developed. IDA can be conducted with a generic record set, eliminating the arduous site-specific record selection required to conduct multiple stripe analysis (MSA). Traditional IDA, however, is limited in producing hazard-consistent estimates of engineering demand parameters (EDPs), which LOHC-IDA overcomes. LOHC-IDA improves upon existing methods by: (i) incorporating correlations among engineering demand parameters across intensity levels and (ii) using peak ground acceleration (PGA) to predict peak floor acceleration (PFA). For two case-study steel BRBF buildings, LOHC-IDA estimates the EAL and loss distributions conditioned on the intensity level that closely match the MSA results, with an average absolute error of 5%. The influence of factors beyond gusset plate design on the λ꜀ and EAL of 26 case-study steel BRBF buildings (designed in accordance with TS 1170.5) is examined. Hazard-consistent λ꜀ and EAL for these buildings are estimated using the FEMA P-58 loss and risk assessment framework. Among the 26 case-study buildings, 23 satisfy the maximum code-specified λ꜀ limit of 10−⁴. The EAL, normalised by the total building replacement cost, is highest for 2-storey BRBFs (0.22% on average), followed by 4-storey BRBFs (0.16% on average) and 8-storey BRBFs (0.11% on average). Reducing the design drift limit has the most significant effect on lowering λ꜀ (all BRBF designs were drift governed), followed by transitioning from pinned to moment-resisting beam-column connections, reducing the brace angle, and increasing brace strength. BRBF buildings designed using the equivalent lateral force method, on average, have a lower λ꜀ compared to those designed using the modal response spectrum method. Diagonally configured BRBFs exhibit the lowest λ꜀, followed by super- X and chevron configured BRBFs. Most design variables, apart from drift limit and beam-column connection, have limited influence on EAL. A simple method for EDP-targeted design of steel BRBF buildings is proposed. For this purpose, linear regression and CatBoost machine learning models are developed to predict steel BRBF building EDPs using peak storey drift ratio (PSDR) and PFA estimates from the 26 case-study buildings at intensity levels ranging from 80% to 0.5% probability of exceedance in 50 years. The R²ₐₔⱼ of these models is around 0.98, while the average prediction error is less than 10%. Fundamental period (T₁), total building height (Hₜ), and pseudospectral acceleration at T₁, denoted as Sₐ(T₁), are selected as the features to predict PSDR, while T₁, Hₜ, and PGA are the features selected to predict PFA. The EDP-targeted design has three steps: (i) for a given Hₜ value, the PSDR prediction model is used to identify a suitable T₁ that can achieve a desired PSDR target at the design intensity, (ii) a force-based design is then conducted iteratively to achieve the target T₁ by using an appropriate ductility factor and design drift limit, and (iii) based on the T₁ in the final design iteration, the PFA demand estimated by the PFA prediction models is used as a conservative input for the design of acceleration-sensitive non-structural elements. An equation to predict λ꜀ at the design stage is proposed for collapse risk-targeted seismic design of buildings. This equation comprises three principal components: reserve building strength, a proxy for effective structural stiffness, and reserve building deformation capacity. This equation is calibrated for the collapse risk-targeted design of BRBF buildings in New Zealand using results from 26 case-study BRBF buildings. The validity of this equation is demonstrated with three design verification examples designed to specific λ꜀ targets. Considering λ꜀ from hazard-consistent incremental dynamic analysis as the benchmark, the mean absolute percentage error in the design-stage prediction of λ꜀ of the verification buildings is approximately 10%.
Most people exposed to disasters cope well. Others, however, develop posttraumatic stress disorder (PTSD)–a mental disorder characterised by symptoms of intrusion, avoidance, and hyperarousal–requiring input from specialist mental health services. To date, relatively little research has evaluated these services, and less is known about characteristics of people seeking treatment and their treatment outcomes. In 2010 and 2011, a series of major earthquakes occurred in the Canterbury region of Aotearoa New Zealand, resulting in initiation of the Adult Specialist Services for Earthquake Trauma Treatment (ASSETT) service to provide cognitive behavioural therapy (CBT) for people with earthquake-related PTSD or subthreshold PTSD symptoms. The current research used systematic literature review methods, in conjunction with data collected from people seeking treatment with the ASSETT service, to address issues relevant to the development of disaster mental health responses, particularly specialist mental health services. A systematic review was conducted synthesising research examining mental health service use among adults exposed to natural disasters. A second systematic review and meta-analysis evaluated psychological interventions for earthquake-related PTSD. A series of studies then utilised diagnostic interview and self-report data collected from people seeking treatment with the ASSETT service (n = 184). Data were collected on factors relating to sociodemographics, pre-earthquake mental disorders, current psychological functioning, degree of objective and subjective earthquake exposure, and life events. These studies examined factors distinguishing treatment-seeking participants from earthquake-exposed Canterbury residents who coped well; differences associated with different prior mental disorders and timing of treatment presentation; and outcomes of CBT provided by the service. Four overarching themes emerged across study findings. The first related to the role of objective and subjective disaster exposure in the development of post-disaster mental health outcomes. Subjective peritraumatic responses were found to be an important factor distinguishing treatment-seeking participants from those who coped well following the earthquakes, independent of objective exposure severity. Heightened peritraumatic responses were also associated with poorer treatment outcome, although not beyond their association with pre-treatment PTSD severity and degree of comorbidity. The second theme related to the role of pre-trauma mental health in the development of post-disaster mental health outcomes. Participants with a history of pre-earthquake mental disorder presented with more comorbid disorders than participants with no prior disorder, but reported comparable degrees of PTSD severity and similar treatment outcomes. The third theme related to temporal considerations for disaster mental health responses. Participants who presented at later time points tended to be older and were more likely to have subthreshold PTSD symptoms, but had similar treatment outcomes as those who presented at earlier time points. The fourth theme related to treatment of severe and ongoing earthquake-related distress. CBT without a formal exposure component was associated with clinically significant improvements on a range of outcome measures, with group and individual-based treatment associated with comparable outcomes. Findings of the current research suggest people seeking treatment for severe and ongoing disaster-related distress are not homogenous, and are likely to present for treatment at different time points, have varied mental health histories, and report diverse disaster experiences. CBT is an effective treatment for severe and ongoing post-disaster distress when delivered in real-world mental health service settings. Group CBT represents an efficient, scalable, and effective treatment format for post-disaster distress, and may be an attractive option for treating widespread need using limited resources.
Background The 2010/2011 Canterbury earthquakes and aftershocks in New Zealand caused unprecedented destruction to the physical, social, economic, and community fabric of Christchurch city. The recovery phase in Christchurch is on going, six years following the initial earthquake. Research exploring how disabled populations experience community inclusion in the longer-term recovery following natural disasters is scant. Yet such information is vital to ensure that recovering communities are inclusive for all members of the affected population. This thesis specifically examined how people who use wheelchairs experienced community inclusion four years following the 2010/2011 Canterbury earthquakes. Aims The primary research aim was to understand how one section of the disability community – people who use wheelchairs – experienced community inclusion over the four years following the 2010/2011 Canterbury earthquakes and aftershocks. A secondary aim was to test a novel sampling approach, Respondent Driven Sampling, which had the potential to enable unbiased population-based estimates. This was motivated by the lack of an available sampling frame for the target population, which would inhibit recruitment of a representative sample. Methodology and methods An exploratory sequential mixed methods design was used, beginning with a qualitative phase (Phase One), which informed a second quantitative phase (Phase Two). The qualitative phase had two stages. First, a small sample of people who use wheelchairs participated in an individual, semi-structured interview. In the second stage, these participants were then invited to a group interview to clarify and prioritise themes identified in the individual interviews. The quantitative phase was a cross-sectional survey developed from the findings from Phase One. Initially, Respondent Driven Sampling was employed to conduct a national, electronic cross-sectional survey that aimed to recruit a sample that may provide unbiased population-based estimates. Following the unsuccessful application of Respondent Driven Sampling, a region-specific convenience sampling approach was used. The datasets from the qualitative and quantitative phases were integrated to address the primary aim of the research. Results In Phase One 13 participants completed the individual interviews, and five of them contributed to the group interview. Thematic analysis of individual and group interview data suggested that participants felt the 2010/11 earthquakes magnified many pre-existing barriers to community inclusion, and also created an exciting opportunity for change. This finding was encapsulated in five themes: 1) earthquakes magnified barriers, 2) community inclusion requires energy, 3) social connections are important, 4) an opportunity lost, and 5) an opportunity found. The findings from Phase One informed the development of a survey instrument to investigate how these findings generalised to a larger sample of individuals who use wheelchairs. In Phase Two, the Respondent Driven Sampling approach failed to recruit enough participants to satisfy the statistical requirements needed to reach equilibrium, thereby enabling the calculation of unbiased population estimates. The subsequent convenience sampling approach recruited 49 participants who, combined with the 15 participants from the Respondent Driven Sampling approach that remained eligible for the region-specific sample, resulted in the total of 64 individuals who used wheelchairs and were residents of Christchurch. Participants reported their level of community inclusion at three time periods: the six months prior to the first earthquake in September 2010 (time one), the six months following the first earthquake in September 2010 (time two), and the six months prior to survey completion (between October 2015 and March 2016, (time three)). Survey data provided some precision regarding the timing in which the magnified barriers developed. Difficulty with community inclusion rose significantly between time one and time two, and while reducing slightly, was still present during time three, and had not returned to the time one baseline. The integrated findings from Phase One and Phase Two suggested that magnified barriers to community inclusion had been sustained four years post-earthquake, and community access had not returned to pre-earthquake levels, let alone improved beyond pre-earthquake levels. Conclusion Findings from this mixed methods study suggest that four years following the initial earthquake, participants were still experiencing multiple magnified barriers, which contributed to physical and social exclusion, as well as fatigue, as participants relied on individual agency to negotiate such barriers. Participants also highlighted the exciting opportunity to create an accessible city. However because they were still experiencing barriers four years following the initial event, and were concerned that this opportunity might be lost if the recovery proceeds without commitment and awareness from the numerous stakeholders involved in guiding the recovery. To truly realise the opportunity to create an accessible city following a disaster, the transition from the response phase to a sustainable longer-term recovery must adopt a new model of community engagement where decision-makers partner with people living with disability to co-produce a vision and strategy for creating an inclusive community. Furthermore, despite the unsuccessful use of Respondent Driven Sampling in this study, future research exploring the application of RDS with wheelchair users is recommended before discounting this sampling approach in this population.
In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.
In this thesis, focus is given to develop methodologies for rapidly estimating specific components of loss and downtime functions. The thesis proposes methodologies for deriving loss functions by (i) considering individual component performance; (ii) grouping them as per their performance characteristics; and (iii) applying them to similar building usage categories. The degree of variation in building stock and understanding their characteristics are important factors to be considered in the loss estimation methodology and the field surveys carried out to collect data add value to the study. To facilitate developing ‘downtime’ functions, this study investigates two key components of downtime: (i) time delay from post-event damage assessment of properties; and (ii) time delay in settling the insurance claims lodged. In these two areas, this research enables understanding of critical factors that influence certain aspects of downtime and suggests approaches to quantify those factors. By scrutinising the residential damage insurance claims data provided by the Earthquake Commission (EQC) for the 2010- 2011 Canterbury Earthquake Sequence (CES), this work provides insights into various processes of claims settlement, the time taken to complete them and the EQC loss contributions to building stock in Christchurch city and Canterbury region. The study has shown diligence in investigating the EQC insurance claim data obtained from the CES to get new insights and build confidence in the models developed and the results generated. The first stage of this research develops contribution functions (probabilistic relationships between the expected losses for a wide range of building components and the building’s maximum response) for common types of claddings used in New Zealand buildings combining the probabilistic density functions (developed using the quantity of claddings measured from Christchurch buildings), fragility functions (obtained from the published literature) and cost functions (developed based on inputs from builders) through Monte Carlo simulations. From the developed contribution functions, glazing, masonry veneer, monolithic and precast concrete cladding systems are found to incur 50% loss at inter-storey drift levels equal to 0.027, 0.003, 0.005 and 0.011, respectively. Further, the maximum expected cladding loss for glazing, masonry veneer, monolithic, precast concrete cladding systems are found to be 368.2, 331.9, 365.0, and 136.2 NZD per square meter of floor area, respectively. In the second stage of this research, a detailed cost breakdown of typical buildings designed and built for different purposes is conducted. The contributions of structural and non- structural components to the total building cost are compared for buildings of different usages, and based on the similar ratios of non-structural performance group costs to the structural performance group cost, four-building groups are identified; (i) Structural components dominant group: outdoor sports, stadiums, parkings and long-span warehouses, (ii) non- structural drift-sensitive components dominant group: houses, single-storey suburban buildings (all usages), theatres/halls, workshops and clubhouses, (iii) non-structural acceleration- sensitive components dominant group: hospitals, research labs, museums and retail/cold stores, and (iv) apartments, hotels, offices, industrials, indoor sports, classrooms, devotionals and aquariums. By statistically analysing the cost breakdowns, performance group weighting factors are proposed for structural, and acceleration-sensitive and drift-sensitive non-structural components for all four building groups. Thus proposed building usage groupings and corresponding weighting factors facilitate rapid seismic loss estimation of any type of building given the EDPs at storey levels are known. A model for the quantification of post-earthquake inspection duration is developed in the third stage of this research. Herein, phase durations for the three assessment phases (one rapid impact and two rapid building) are computed using the number of buildings needing inspections, the number of engineers involved in inspections and a phase duration coefficient (which considers the median building inspection time, efficiency of engineer and the number of engineers involved in each assessment teams). The proposed model can be used: (i) by national/regional authorities to decide the length of the emergency period following a major earthquake, and estimate the number of engineers required to conduct a post-earthquake inspection within the desired emergency period, and (ii) to quantify the delay due to inspection for the downtime modelling framework. The final stage of this research investigates the repair costs and insurance claim settlement time for damaged residential buildings in the 2010-2011 Canterbury earthquake sequence. Based on the EQC claim settlement process, claims are categorized into three groups; (i) Small Claims: claims less than NZD15,000 which were settled through cash payment, (ii) Medium Claims: claims less than NZD100,000 which were managed through Canterbury Home Repair Programme (CHRP), and (iii) Large Claims: claims above NZD100,000 which were managed by an insurance provider. The regional loss ratio (RLR) for greater Christchurch for three events inducing shakings of approximate seismic intensities 6, 7, and 8 are found to be 0.013, 0.066, and 0.171, respectively. Furthermore, the claim duration (time between an event and the claim lodgement date), assessment duration (time between the claim lodgement day and the most recent assessment day), and repair duration (time between the most recent assessment day and the repair completion day) for the insured residential buildings in the region affected by the Canterbury earthquake sequence is found to be in the range of 0.5-4 weeks, 1.5- 5 months, and 1-3 years, respectively. The results of this phase will provide useful information to earthquake engineering researchers working on seismic risk/loss and insurance modelling.
The Canterbury Earthquakes of 2010-2011, in particular the 4th September 2010 Darfield earthquake and the 22nd February 2011 Christchurch earthquake, produced severe and widespread liquefaction in Christchurch and surrounding areas. The scale of the liquefaction was unprecedented, and caused extensive damage to a variety of man-made structures, including residential houses. Around 20,000 residential houses suffered serious damage as a direct result of the effects of liquefaction, and this resulted in approximately 7000 houses in the worst-hit areas being abandoned. Despite the good performance of light timber-framed houses under the inertial loads of the earthquake, these structures could not withstand the large loads and deformations associated with liquefaction, resulting in significant damage. The key structural component of houses subjected to liquefaction effects was found to be their foundations, as these are in direct contact with the ground. The performance of house foundations directly influenced the performance of the structure as a whole. Because of this, and due to the lack of research in this area, it was decided to investigate the performance of houses and in particular their foundations when subjected to the effects of liquefaction. The data from the inspections of approximately 500 houses conducted by a University of Canterbury summer research team following the 4th September 2010 earthquake in the worst-hit areas of Christchurch were analysed to determine the general performance of residential houses when subjected to high liquefaction loads. This was followed by the detailed inspection of around 170 houses with four different foundation types common to Christchurch and New Zealand: Concrete perimeter with short piers constructed to NZS3604, concrete slab-on-grade also to NZS3604, RibRaft slabs designed by Firth Industries and driven pile foundations. With a focus on foundations, floor levels and slopes were measured, and the damage to all areas of the house and property were recorded. Seven invasive inspections were also conducted on houses being demolished, to examine in more detail the deformation modes and the causes of damage in severely affected houses. The simplified modelling of concrete perimeter sections subjected to a variety of liquefaction-related scenarios was also performed, to examine the comparative performance of foundations built in different periods, and the loads generated under various bearing loss and lateral spreading cases. It was found that the level of foundation damage is directly related to the level of liquefaction experienced, and that foundation damage and liquefaction severity in turn influence the performance of the superstructure. Concrete perimeter foundations were found to have performed most poorly, suffering high local floor slopes and being likely to require foundation repairs even when liquefaction was low enough that no surface ejecta was seen. This was due to their weak, flexible foundation structure, which cannot withstand liquefaction loads without deforming. The vulnerability of concrete perimeter foundations was confirmed through modelling. Slab-on-grade foundations performed better, and were unlikely to require repairs at low levels of liquefaction. Ribraft and piled foundations performed the best, with repairs unlikely up to moderate levels of liquefaction. However, all foundation types were susceptible to significant damage at higher levels of liquefaction, with maximum differential settlements of 474mm, 202mm, 182mm and 250mm found for concrete perimeter, slab-on-grade, ribraft and piled foundations respectively when subjected to significant lateral spreading, the most severe loading scenario caused by liquefaction. It was found through the analysis of the data that the type of exterior wall cladding, either heavy or light, and the number of storeys, did not affect the performance of foundations. This was also shown through modelling for concrete perimeter foundations, and is due to the increased foundation strengths provided for heavily cladded and two-storey houses. Heavy roof claddings were found to increase the demands on foundations, worsening their performance. Pre-1930 concrete perimeter foundations were also found to be very vulnerable to damage under liquefaction loads, due to their weak and brittle construction.
It is not a matter of if a major earthquake will happen in New Zealand, it is when. Earthquakes wreak havoc, cut off power and water supply, lines of communication, sewer, supply chains, and transport infrastructure. People get injured and whole communities can get cut off the rest of the country for extended periods of time. Countries taking measures to increase the population's preparedness tend to suffer less severe consequences than those that do not. Disaster management authorities deliver comprehensive instructions and preparation guidance, yet communities remain grossly underprepared. There are multiple factors that influence motivation for preparedness. Personal experience is one of the most significant factors that influence preparedness motivation. Not many people will experience a severe and damaging earthquake in their lifetime. A serious game (SG) that is a computer simulation of an earthquake is a tool that can let participants experience the earthquake and its aftermath from the safety of their computer. The main result of this research is a positive answer to the question: Can a serious game motivate people to prepare for earthquakes at least just as good as a personal experience of at least a moderate earthquake? There are different levels of immersion this serious game can be implemented at. In this thesis the same earthquake experience scenario – SG “ShakeUp” is implemented as a desktop application and a virtual reality (VR) application. A user study is conducted with the aim of comparing the motivation level achieved by the two versions of the SG “ShakeUp”. In this study no benefits of using VR over traditional desktop application were found: participants trying both versions of the SG “ShakeUp” reported similar levels of motivation to prepare for earthquakes immediately after the experiment. This means that both versions of the experience were equally effective in motivating participants to prepare for earthquakes. An additional benefit of this result is that the cheaper and easier to deliver desktop version can be widely used in various education campaigns. Participants reported being more motivated to prepare for earthquakes by either version of the SG “ShakeUp” than by any other contributing factor, including their previous earthquake experience or participation in a public education campaign. Both versions of the SG “ShakeUp” can successfully overcome personal bias, unrealistic optimism, pessimism, lack of perceived control over one’s earthquake preparation actions, fatalism, and sense of helplessness in the face of the earthquakes and motivate the individual to prepare for earthquakes. Participants without the prior earthquake experience benefit most from the SG “ShakeUp” regardless of the version tried, compared to the participants who had experienced an earthquake: significantly more of them will reconsider their current level of earthquake preparedness; about 24% more of them attribute their increased level of motivation to prepare for earthquakes to the SG “ShakeUp”. For every earthquake preparation action there is about 25% more people who felt motivated to do it after trying the SG “ShakeUp” than those who have done this preparation action before the experiment. After trying either version of the SG “ShakeUp”, people who live in a free standing house and those who live in a rental property reported highest levels of intent to carry on with the preparation actions. The proposed application prototype has been discussed with the University of Canterbury Earthquake Centre and received very positive feedback as having potential for practical use by various disaster management authorities and training institutions. The research shows that the SG “ShakeUp” motivates people to prepare for earthquakes as good as a personal earthquake experience and can be successfully used in various education campaigns.
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.
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.
Perimeter Moment resisting steel frames (PMRSFs) are a commonly used seismic resisting system, placed around the perimeter of the building for maximum torsional stiffness. They are typically designed as “strong column weak beam” systems with fixed column bases. When subjected to severe earthquake demand, sufficient to push the beams into the inelastic range, it is expected that plastic hinging at the column bases will occur. However, the response of PMRSF systems to the severe 2010/2011 Christchurch earthquake series did not generate column base hinging in systems which exhibited beam yielding.
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.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
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.
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.
Liquefaction of sandy soil has been observed to cause significant damage to infrastructure during major earthquakes. Historical cases of liquefaction have typically occurred in sands containing some portion of fines particles, which are defined as 75μm or smaller in diameter. The effects of fines on the undrained behaviour of sand are not however fully understood, and this study therefore attempts to quantify these effects through the undrained testing of sand mixed with non-plastic fines sourced from Christchurch, New Zealand. The experimental program carried out during this study consisted of undrained monotonic and cyclic triaxial tests performed on three different mixtures of sand and fines: the Fitzgerald Bridge mixture (FBM), and two Pinnacles Sand mixtures (PSM1 and PSM2). The fines content of each host sand was systematically varied up to a maximum of 30%, with all test specimens being reconstituted using moist tamping deposition. The undrained test results from the FBM soils were interpreted using a range of different measures of initial state. When using void ratio and relative density, the addition of fines to the FBM sand caused more contractive behaviour for both monotonic and cyclic loadings. This resulted in lower strengths at the steady state of deformation, and lower liquefaction resistances. When the intergranular void ratio was used for the interpretation, the effect of additional fines was to cause less contractive response in the sand. The state parameter and state index were also used to interpret the undrained cyclic test results – these measures suggested that additional fines caused less contractive sand behaviour, the opposite to that observed when using the void ratio. This highlighted the dependency on the parameter chosen as a basis for the response comparison when determining the effects of fines, and pointed out a need to identify a measure that normalizes such effects. Based on the FBM undrained test results and interpretations, the equivalent granular void ratio, e*, was identified from the literature as a measure of initial state that normalizes the effects of fines on the undrained behaviour of sand up to a fines content of 30%. This is done through a parameter within the e* definition termed the fines influence factor, b, which quantifies the effects of fines from a value of zero (no effect) to one (same effect as sand particles). The value of b was also determined to be different when interpreting the steady state lines (bSSL) and cyclic resistance curves (bCR) respectively for a given mixture of sand and fines. The steady state lines and cyclic resistance curves of the FBM soils and a number of other sand-fines mixtures sourced from the literature were subsequently interpreted using the equivalent granular void ratio concept, with bSSL and bCR values being back-calculated from the respective test data sets. Based on these interpretations, it was concluded that e* was conceptually a useful parameter for characterizing and quantifying the effects of fines on the undrained behaviour of sand, assuming the fines influence factor value could be derived. To allow prediction of the fines influence factor values, bSSL and bCR were correlated with material and depositional properties of the presented sand-fines mixtures. It was found that as the size of the fines particles relative to the sand particles became smaller, the values of bSSL and bCR reduced, indicating lower effect of fines. The same trend was also observed as the angularity of the sand particles increased. The depositional method was found to influence the value of bCR, due to the sensitivity of cyclic loading to initial soil fabric. This led to bSSL being used as a reference for the effect of fines, with specimens prepared by moist tamping having bCR > bSSL, and specimens prepared by slurry deposition having bCR < bSSL. Finally the correlations of the fines influence factor values with material and depositional properties were used to define the simplified estimation method – a procedure capable of predicting the approximate steady state lines and cyclic resistance curves of a sand as the non-plastic fines content is increased up to 30%. The method was critically reviewed based on the undrained test results of the PSM1 and PSM2 soils. This review suggested the method could accurately predict undrained response curves as the fines content was raised, based on the PSM1 test results. It also however identified some key issues with the method, such as the inability to accurately predict the responses of highly non-uniform soils, a lack of consideration for the entire particle size distribution of a soil, and the fact the errors in the prediction of bSSL carry through into the prediction of bCR. Lastly some areas of further investigation relating to the method were highlighted, including the need to verify the method through testing of sandy soils sourced from outside the Christchurch area, and the need to correlate the value of bCR with additional soil fabrics / depositional methods.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
This thesis addresses the topic of local bond behaviour in RC structures. The mechanism of bond refers to the composite action between deformed steel reinforcing bars and the surrounding concrete. Bond behaviour is an open research topic with a wide scope, particularly because bond it is such a fundamental concept to structural engineers. However, despite many bond-related research findings having wide applications, the primary contribution of this research is an experimental evaluation of the prominent features of local bond behaviour and the associated implications for the seismic performance of RC structures. The findings presented in this thesis attempt to address some structural engineering recommendations made by the Canterbury Earthquakes Royal Commission following the 2010-2011 Canterbury (New Zealand) earthquake sequence. A chapter of this thesis discusses the structural behaviour of flexure-dominated RC wall structures with an insufficient quantity of longitudinal reinforcement, among other in situ conditions, that causes material damage to predominantly occur at a single crack plane. In this particular case, the extent of concrete damage and bond deterioration adjacent to the crack plane will influence the ductility capacity that is effectively provided by the reinforcing steel. As a consequence of these in situ conditions, some lightly reinforced wall buildings in Christchurch lost their structural integrity due to brittle fracture of the longitudinal reinforcement. With these concerning post-earthquake observations in mind, there is the underlying intention that this thesis presents experimental evidence of bond behaviour that allows structural engineers to re-assess their confidence levels for the ability of lightly reinforced concrete structures to achieve the life-safety seismic performance objective the ultimate limit state. Three chapters of this thesis are devoted to the experimental work that was conducted as the main contribution of this research. Critical details of the experimental design, bond testing method and test programme are reported. The bond stress-slip relationship was studied through 75 bond pull-out tests. In order to measure the maximum local bond strength, all bond tests were carried out on deformed reinforcing bars that did not yield as the embedded bond length was relatively short. Bond test results have been presented in two separate chapters in which 48 monotonic bond tests and 27 cyclic bond tests are presented. Permutations of the experiments include the loading rate, cyclic loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and position of the embedded bond region within the specimen (close or far away to the free surface). The parametric study showed that the concrete strength significantly influences the maximum bond strength and that it is reasonable to normalise the bond stress by the square-root of the concrete compressive strength, √(f'c). The generalised monotonic bond behaviour is described within. An important outcome of the research is that the measured bond strength and stiffness was higher than stated by the bond stress-slip relationship in the fib Model Code 2010. To account for these observed differences, an alternative model is proposed for the local monotonic bond stress-slip relationship. Cyclic bond tests showed a significant proportion of the total bond degradation occurs after the loading cycle in the peak bond strength range, which is when bond slip has exceeded 0.5 mm. Subsequent loading to constant slip values showed a linear relationship between the amount of bond strength degradation and the log of the number of cycles that were applied. To a greater extent, the cyclic bond deterioration depends on the bond slip range, regardless of whether the applied load cycling is half- or fully-reversed. The observed bond deterioration and hysteretic energy dissipated during cyclic loading was found to agree reasonably well between these cyclic tests with different loading protocols. The cyclic bond deterioration was also found to be reasonably consistent exponential damage models found in the literature. This research concluded that the deformed reinforcing bars used in NZ construction, embedded in moderate to high strength concrete, are able to develop high local bond stresses that are mobilised by a small amount of local bond slip. Although the relative rib geometry was not varied within this experimental programme, a general conclusion of this thesis is that deformed bars currently available in NZ have a relative rib bearing area that is comparatively higher than the test bars used in previous international research. From the parametric study it was found that the maximum monotonic bond strength is significant enhanced by dynamic loading rates. Experimental evidence of high bond strength and initial bond stiffness generally suggests that only a small amount of local bond slip that can occur when the deformed test bar was subjected to large tension forces. Minimal bond slip and bond damage limits the effective yielding length that is available for the reinforcing steel to distribute inelastic material strains. Consequently, the potential for brittle fracture of the reinforcement may be a more problematic and widespread issue than is apparent to structural engineers. This research has provided information that improve the reliability of engineering predictions (with respect to ductility capacity) of maximum crack widths and the extent of bond deterioration that might occur in RC structures during seismic actions.
Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.
Ongoing climate change triggers increasing temperature and more frequent extreme events which could limit optimal performance of haliotids, affect their physiology and biochemistry as well as influencing their population structure. Haliotids are a valuable nearshore fishery in a number of countries and many are showing a collapse of stocks because of overexploitation, environmental changes, loss of habitat, and disease. The haliotid in New Zealand commonly referred to as the blackfoot pāua (Haliotis iris) contribute a large and critical cultural, recreational and economic resource. Little was known about pāua responses to increasing temperature and acute environmental factors, as well as information about population size structure in Kaikoura after the earthquake 2016 and in Banks Peninsula. The aims of this study were to investigate the effects of temperature on scope for growth (SfG); physiological and biochemical responses of pāua subjected to different combined stressors including acute temperature, acute salinity and progressive hypoxia; and describe population size structure and shell morphology in different environments in Kaikoura and Banks Peninsula. The main findings of the present study found that population size structures of pāua were site-specific, and the shell length and shell height ratio of 3.25 could distinguish between stunted and non-stunted populations. The study found that high water temperature resulted in a reduction in absorbed energy from food, an increase in respiration energy, and ammonia excretion energy. Surveys were conducted at six study sites around the Canterbury Region over three years in order to better understand the population size structure and shell morphology of pāua. The findings found that the population size structure at 6 sites differed. Both juveniles and adults were found in intertidal areas at five sites. However, at Cape Three Points, pāua were found only in subtidal zones. One of the sites, Little Port Cooper, had a stunted population where only two pāua reached 125 mm in length over three years. In addition, most pāua in Little Port Cooper and Cape Three Points were adults, while Seal Reef had mostly juveniles. Wakatu Quay and Omihi had a full size range of pāua. Oaro population was dominated with juveniles and sub-adults. Recruitment and growth of pāua were successful after the earthquake in 2016. Research into pāua shell morphologies also determined that shell dimensions differed between sites. The relationships of shell length to shell width were linear and the relationship of shell length to shell height was curvilinear. Interestingly, SL:SH ratio of 3.25 is able to be used to identify stunted and non-stunted populations for pāua larger than 90 mm in length. Little Port Cooper was a stunted population with mean SL:SH ratio being 3.16. In the laboratory, scope for growth of pāua was investigated at four different temperatures of 12oC, 15oC, 18oC and 21oC over four weeks’ acclimation. The current study has found that SfG of pāua highly depended on temperature. Absorbed energy and respiration energy accounted for the highest proportion of the SfG of pāua. The respiration energy of pāua accounted for approximately 36%, 40%, 49% and 69% of the absorbed energy at 12°C, 15°C, 18°C and 21°C, respectively. The pāua at all acclimation temperatures had a positive scope for growth. The study suggested that the SfG was highest at 15°C, while the value at 21°C was the lowest. However, SfG at 18°C and 21°C decreased after 14 days of acclimation. Because of maintaining almost unchanged oxygen consumption over four weeks’ acclimation, pāua showed their poor abilities to acclimate to an increase in temperature. Therefore, they may be more vulnerable in future warming scenarios. The physiological and biochemical responses of pāua toward different combined stressors included three experiments. In terms of the acute temperature experiment, pāua were acclimated at 12oC, 15oC, 18oC or 21oC for two weeks before stepwise exposure to four temperatures of 12oC, 15oC, 18oC and 21oC every 4 hours. The acute salinity change, pāua were acclimated at 12oC, 15oC or 18oC over two weeks. Pāua were then exposed to a stepwise decrease of salinity of 2‰ every two hours from 34 – 22‰. Regarding the declining oxygen level, pāua were acclimated at 15 oC or 18oC for two weeks before exposure to one of four temperatures at 12oC, 15oC, 18oC or 21oC in one hour. After that acute progressive hypoxia was studied in closed respirometers for around six hours. The findings showed that there were interactions between combined stressors, affecting physiology of pāua (metabolism and heart rate). This suggests that environmental factors do not have a separate effect, but they also have interactions that enhance negative effects on pāua. Also, both oxygen uptake and heart rate responded quickly to temperature change and increased with rising temperature. On the other hand, oxygen uptake and heart rate decreased with reducing salinity and progressive hypoxia (before critical oxygen tension - Pcrit). Pcrit over four acute temperature exposures, ranged between 30.2 and 80.0 mmHg, depending on the exposure temperature. Acclimation temperature, combined with acute temperature, salinity or hypoxia stress affected the biochemistry of pāua. Pāua are osmoconformers so decreased salinity resulted in reducing haemolymph ionic concentration and increasing body volume. They were hypo-ionic with respect to sodium and potassium over the salinity ranges of 34 - 22‰. Haemocyanin accounts for a large pecentage of haemolymph protein, so trends of protein followed haemocyanin. Pāua tended to store oxygen in haemocyanin under extreme salinity stress at 22‰ and extreme hypoxia around 10 mmHg, rather than in oxygen transport. In conclusion, pāua at different sites had different population structures and morphologies. Pāua are sensitive to environmental stressors. They consumed more oxygen at high temperatures because they do not have thermal acclimation capacity. They are also osmoconformers with haemolymph sodium and potassium decreasing with salinity medium. Under progressive hypoxia, pāua could regulate oxygen and heart rate until Pcrit depending on temperature. Acute environmental changes also disturbed haemolyph parameters. 12°C and 15°C could be in the range of optimal temperature with higher SfG and less stress when exposed to acute environmental changes. Meanwhile long term exposure to 21°C is likely to be outside of the optimal range for the pāua. With ongoing climate change, pāua populations are more vulnerable so conservation is necessary. The research contributes to improving fishery management, providing insights into different environmental stressors affecting the energy demand and physiological and biochemical responses of pāua. It also allow to predicting the growth patterns and responses of pāua to adapt to climate change.
Field surveys and experimental studies have shown that light steel or timber framed plasterboard partition walls are particularly vulnerable to earthquake damage prompting the overarching objective of this research, which is to further the development of low damage seismic systems for non-structural partition walls in order to facilitate their adoption by industry to assist with reducing the losses associated with the maintenance and repair cost of buildings across their design life. In particular, this study focused on the behaviour of steel-framed partition walls systems with novel detailing that aim to be “low-damage” designed according to common practice for walls used in commercial and institutional buildings in New Zealand. This objective was investigated by (1) investigating the performance of a flexible track system proposed by researchers and industry by experimental testing of full-scale specimens; (2) investigating the performance of the seismic gap partition wall systems proposed in a number of studies, further developed in this study with input from industry, by experimental testing of full-scale specimens; and (3) investigating the potential implications of using these systems compared with traditionally detailed partition wall systems within multi-storey buildings using the Performance Based Earthquake Engineering loss assessment methodology. Three full-scale testing frames were designed in order to replicate, under controlled laboratory conditions, the effects of seismic shaking on partition walls within multi-storey buildings by the application of quasi-static uni-directional cyclic loading imposing an inter-storey drift. The typical configuration for test specimens was selected to be a unique “y-shape”, including one angled return wall, with typical dimensions of approximately 2400 mm along the main wall and 600 mm along (approximately) the returns walls with a height of 2405 mm from floor to ceiling. The specimens were aligned within test frames at an oblique angle to the direction of loading in order to investigate bi- directional effects. Three wall specimens with flexible track detailing, two identical plane specimens and the third including a doorway, were tested. The detailing involved removing top track anchors within the proximity of wall intersections, thus allowing the tracks to ‘bow’ out at these locations. Although the top track anchors were specified to be removed the proximity of wall intersections, a construction error was made whereby a single top track slab to concrete anchor was left in at the three-way wall junction. Despite this error, the experimental testing was deemed worthwhile since such errors will also occur in practice and because the behaviour of the wall can be examined with this fixing in mind. The specimens also included an acoustic/fire sealant at the top lining to floor boundary. In addition to providing drift capacities, the force-displacement behaviour is also reported, the dissipated energy was computed, and the parameters of the Wayne-Stewart hysteretic model were fitted to the results. The specimen with the door opening behaved significantly different to the plane specimens: damage to the doorway specimen began as cracking of the wallboard propagating from the corners of the doorway following which the L- and Y- shaped junctions behaved independently, whereas damage to the plane specimens began as cracking of the wallboard at the top of the L-junction and wall system deformed as a single unit. The results suggest that bi-directional behaviour is important even if its impact cannot be directly quantified by the experiments conducted. Damage to sealant implies that the bond between plasterboard and sealant is important for its seismic performance. Careful quality control is advised as defects in the bond may significantly impact its ability to withstand seismic movement. Two specimens with seismic gap detailing were tested: a steel stud specimen and a timber stud specimen. Observed drift capacities were significantly greater than traditional plasterboard partition systems. Equations were used to predict the drift at which damage state 1 (DS1) and damage state 2 (DS2) would initiate. The equation used to estimate the drift at the onset of DS1 accurately predicted the onset of plaster cracking but overestimated the drift at which the gap filling material was damaged. The equation used to predict the onset of DS2 provided a lower bound for both specimens and also when used to predict results of previous experimental tests on seismic gap systems. The gap-filling material reduced the drift at the onset of DS1, however, it had a beneficial effect on the re-centring behaviour of the linings. Out-of-plane displacements and return wall configuration did not appear to significantly impact the onset of plaster cracking in the specimens. A loss assessment according to the PBEE methodology was conducted on four steel MRF case study buildings: (1) a 4-storey building designed for the Christchurch region, (2) a 4-storey building designed for the Wellington region, (3) a 12-storey building designed for the Christchurch region, and (4) a 12- storey building designed for the Wellington region. The fragility parameters for a traditional partition system, the flexible track partition system, and the seismic gap steel stud and timber stud partition systems were included within the loss assessment. The order (lowest to highest) of each system in terms of the expected annual losses of each building when incorporating the system was, (1) the seismic gap timber stud system, (2) the seismic gap steel stud system, (3) the traditional/baseline system, and (4) the flexible track system. For the seismic gap timber stud system, which incurred the greatest reduction in expected annual losses for each case study building, the reduction in expected annual losses in comparison to the losses found when using the traditional system ranged from a 5% to a 30% reduction. This reinforces the fact that while there is a benefit to the using low damage partition systems in each building the extent of reduction in expected annual losses is significantly dependent on the particular building design and its location. The flexible track specimens had larger repair costs at small hazard levels compared to the traditional system but smaller repair costs at larger hazard levels. However, the resulting expected annual losses for the flexible track system was higher than the traditional system which reinforces findings from past studies which observed that the greatest contribution to expected annual losses arises from low to moderate intensity shaking seismic events (low hazard levels).
Case study unreinforced masonry (URM) buildings that were seismically retrofitted prior to the 2010/11 Canterbury earthquake sequence and exhibited successful performance during these earthquakes are presented herein. Selected buildings were divided into the following categories based on size and complexity: (1) simple, single storey box type buildings (i.e. electrical substations), (2) common and simple commercial buildings, and (3) large and complex clay brick and stone URM buildings. The retrofitted case study URM buildings were evaluated based on overall structural seismic performance as well as the categories of initial seismic design, heritage preservation, architectural appeal, and cost. Detailed observations of 4 representative case study buildings and a summary of findings are reported herein. http://db.nzsee.org.nz/2017/Orals.htm
Test results are presented for wall-diaphragm plate anchor connections that were axially loaded to rupture. These connection samples were extracted post-earthquake by sorting through the demolition debris from unreinforced masonry (URM) buildings damaged in the Christchurch earthquakes. Unfortunately the number of samples available for testing was small due to the difficulties associated with sample collection in an environment of continuing aftershocks and extensive demolition activity, when personal safety combined with commercial activity involving large demolition machinery were imperatives that inhibited more extensive sample collection for research purposes. Nevertheless, the presented data is expected to be of assistance to structural engineers undertaking seismic assessment of URM buildings that have existing wall-diaphragm anchor plate connections installed, where it may be necessary to estimate the capacity of the existing connection as an important parameter linked with determining the current seismic capacity of the building and therefore influencing the decision regarding whether supplementary connections should be installed.
Axial elongation of reinforced concrete (RC) plastic hinges has previously been observed in a range of laboratory experiments, and more recently was observed in several Christchurch buildings following the 2010/2011 Canterbury earthquakes. Axial restraint to plastic hinges is provided by adjacent structural components such as floors as the plastic hinges elongate, which can significantly alter the performance of the plastic hinge and potentially invalidate the capacity design strength hierarchy of the building. Coupling beams in coupled wall systems are particularly susceptible to axial restraint effects due to their importance in the strength hierarchy, the high ductility demands that they experience, and the large stiffness of bounding walls. From computational modelling it has been found that ignoring axial restraint effects when designing coupled walls can result in significantly increased strength, reduced ductility and reduced energy dissipation capacity. The complexity of the topic merits further research to better account for realistic restraint effects when designing coupled walls.
Since the 2010/11 Canterbury earthquakes, Akaroa has been hosting the majority of cruise ship arrivals to Canterbury. This amounts to approximately 70-74 days per season, when between 2,000- 4,000 persons come ashore between 9am and 4pm when in port. This increased level of cruise ship arrivals has had significant impacts, both beneficial and detrimental, on Akaroa. Attitudes within the Akaroa community to hosting cruise ship arrivals appear to be divided, and has led to public debate in Akaroa about the issue. In response to this situation, Christchurch and Canterbury Tourism (CCT) commissioned this research project to assess the impact of cruise ship tourism on the Akaroa community.This research was commissioned and funded by Christchurch and Canterbury Tourism (CCT).
The current seismic design practice for reinforced concrete (RC) walls has been drawn into question following the Canterbury earthquakes. An overview of current research being undertaken at the University of Auckland into the seismic behaviour of RC walls is presented. The main objectives of this research project are to understand the observed performance of several walls in Christchurch, quantify the seismic loads on RC walls, and developed improved design procedures for RC walls that will assist in revisions to NZS 3101. A database summarising of the performance of RC wall buildings in the Christchurch CBD was collated to identify damage modes and case-study buildings. A detailed investigation is underway to verify the seismic performance of lightly reinforced concrete walls and an experimental setup has been developed to subject RC wall specimen to loading that is representative of a multi-storey building. Numerical modelling is being used to understand the observed performance of several case-study RC walls buildings in Christchurch. Of particular interest is the influence that interactions between walls and other structural elements have on the seismic response of buildings and the loads generated on RC walls.
