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.
Structural engineering is facing an extraordinarily challenging era. These challenges are driven by the increasing expectations of modern society to provide low-cost, architecturally appealing structures which can withstand large earthquakes. However, being able to avoid collapse in a large earthquake is no longer enough. A building must now be able to withstand a major seismic event with negligible damage so that it is immediately occupiable following such an event. As recent earthquakes have shown, the economic consequences of not achieving this level of performance are not acceptable. Technological solutions for low-damage structural systems are emerging. However, the goal of developing a low-damage building requires improving the performance of both the structural skeleton and the non-structural components. These non-structural components include items such as the claddings, partitions, ceilings and contents. Previous research has shown that damage to such items contributes a disproportionate amount to the overall economic losses in an earthquake. One such non-structural element that has a history of poor performance is the external cladding system, and this forms the focus of this research. Cladding systems are invariably complicated and provide a number of architectural functions. Therefore, it is important than when seeking to improve their seismic performance that these functions are not neglected. The seismic vulnerability of cladding systems are determined in this research through a desktop background study, literature review, and postearthquake reconnaissance survey of their performance in the 2010 – 2011 Canterbury earthquake sequence. This study identified that precast concrete claddings present a significant life-safety risk to pedestrians, and that the effect they have upon the primary structure is not well understood. The main objective of this research is consequently to better understand the performance of precast concrete cladding systems in earthquakes. This is achieved through an experimental campaign and numerical modelling of a range of precast concrete cladding systems. The experimental campaign consists of uni-directional, quasi static cyclic earthquake simulation on a test frame which represents a single-storey, single-bay portion of a reinforced concrete building. The test frame is clad with various precast concrete cladding panel configurations. A major focus is placed upon the influence the connection between the cladding panel and structural frame has upon seismic performance. A combination of experimental component testing, finite element modelling and analytical derivation is used to develop cladding models of the cladding systems investigated. The cyclic responses of the models are compared with the experimental data to evaluate their accuracy and validity. The comparison shows that the cladding models developed provide an excellent representation of real-world cladding behaviour. The cladding models are subsequently applied to a ten-storey case-study building. The expected seismic performance is examined with and without the cladding taken into consideration. The numerical analyses of the case-study building include modal analyses, nonlinear adaptive pushover analyses, and non-linear dynamic seismic response (time history) analyses to different levels of seismic hazard. The clad frame models are compared to the bare frame model to investigate the effect the cladding has upon the structural behaviour. Both the structural performance and cladding performance are also assessed using qualitative damage states. The results show a poor performance of precast concrete cladding systems is expected when traditional connection typologies are used. This result confirms the misalignment of structural and cladding damage observed in recent earthquake events. Consequently, this research explores the potential of an innovative cladding connection. The outcomes from this research shows that the innovative cladding connection proposed here is able to achieve low-damage performance whilst also being cost comparable to a traditional cladding connection. It is also theoretically possible that the connection can provide a positive value to the seismic performance of the structure by adding addition strength, stiffness and damping. Finally, the losses associated with both the traditional and innovative cladding systems are compared in terms of tangible outcomes, namely: repair costs, repair time and casualties. The results confirm that the use of innovative cladding technology can substantially reduce the overall losses that result from cladding damage.
Globally, the maximum elevations at which treelines are observed to occur coincide with a 6.4 °C soil isotherm. However, when observed at finer scales, treelines display a considerable degree of spatial complexity in their patterns across the landscape and are often found occurring at lower elevations than expected relative to the global-scale pattern. There is still a lack of understanding of how the abiotic environment imposes constraints on treeline patterns, the scales at which different effects are acting, and how these effects vary over large spatial extents. In this thesis, I examined abrupt Nothofagus treelines across seven degrees of latitude in New Zealand in order to investigate two broad questions: (1) What is the nature and extent of spatial variability in Nothofagus treelines across the country? (2) How is this variation associated with abiotic variation at different spatial scales? A range of GIS, statistical, and atmospheric modelling methods were applied to address these two questions. First, I characterised Nothofagus treeline patterns at a 15x15km scale across New Zealand using a set of seven, GIS-derived, quantitative metrics that describe different aspects of treeline position, shape, spatial configuration, and relationships with adjacent vegetation. Multivariate clustering of these metrics revealed distinct treeline types that showed strong spatial aggregation across the country. This suggests a strong spatial structuring of the abiotic environment which, in turn, drives treeline patterns. About half of the multivariate treeline metric variation was explained by patterns of climate, substrate, topographic and disturbance variability; on the whole, climatic and disturbance factors were most influential. Second, I developed a conceptual model that describes how treeline elevation may vary at different scales according to three categories of effects: thermal modifying effects, physiological stressors, and disturbance effects. I tested the relevance of this model for Nothofagus treelines by investigating treeline elevation variation at five nested scales (regional to local) using a hierarchical design based on nested river catchments. Hierarchical linear modelling revealed that the majority of the variation in treeline elevation resided at the broadest, regional scale, which was best explained by the thermal modifying effects of solar radiation, mountain mass, and differences in the potential for cold air ponding. Nonetheless, at finer scales, physiological and disturbance effects were important and acted to modify the regional trend at these scales. These results suggest that variation in abrupt treeline elevations are due to both broad-scale temperature-based growth limitation processes and finer-scale stress- and disturbance-related effects on seedling establishment. Third, I explored the applicability of a meso-scale atmospheric model, The Air Pollution Model (TAPM), for generating 200 m resolution, hourly topoclimatic data for temperature, incoming and outgoing radiation, relative humidity, and wind speeds. Initial assessments of TAPM outputs against data from two climate station locations over seven years showed that the model could generate predictions with a consistent level of accuracy for both sites, and which agreed with other evaluations in the literature. TAPM was then used to generate data at 28, 7x7 km Nothofagus treeline zones across New Zealand for January (summer) and July (winter) 2002. Using mixed-effects linear models, I determined that both site-level factors (mean growing season temperature, mountain mass, precipitation, earthquake intensity) and local-level landform (slope and convexity) and topoclimatic factors (solar radiation, photoinhibition index, frost index, desiccation index) were influential in explaining variation in treeline elevation within and among these sites. Treelines were generally closer to their site-level maxima in regions with higher mean growing season temperatures, larger mountains, and lower levels of precipitation. Within sites, higher treelines were associated with higher solar radiation, and lower photoinhibition and desiccation index values, in January, and lower desiccation index values in July. Higher treelines were also significantly associated with steeper, more convex landforms. Overall, this thesis shows that investigating treelines across extensive areas at multiple study scales enables the development of a more comprehensive understanding of treeline variability and underlying environmental constraints. These results can be used to formulate new hypotheses regarding the mechanisms driving treeline formation and to guide the optimal choice of field sites at which to test these hypotheses.
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 thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
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.
The Canterbury earthquake sequence (2010-2011) was the most devastating catastrophe in New Zealand‘s modern history. Fortunately, in 2011 New Zealand had a high insurance penetration ratio, with more than 95% of residences being insured for these earthquakes. This dissertation sheds light on the functions of disaster insurance schemes and their role in economic recovery post-earthquakes. The first chapter describes the demand and supply for earthquake insurance and provides insights about different public-private partnership earthquake insurance schemes around the world. In the second chapter, we concentrate on three public earthquake insurance schemes in California, Japan, and New Zealand. The chapter examines what would have been the outcome had the system of insurance in Christchurch been different in the aftermath of the Canterbury earthquake sequence (CES). We focus on the California Earthquake Authority insurance program, and the Japanese Earthquake Reinsurance scheme. Overall, the aggregate cost of the earthquake to the New Zealand public insurer (the Earthquake Commission) was USD 6.2 billion. If a similar-sized disaster event had occurred in Japan and California, homeowners would have received only around USD 1.6 billion and USD 0.7 billion from the Japanese and Californian schemes, respectively. We further describe the spatial and distributive aspects of these scenarios and discuss some of the policy questions that emerge from this comparison. The third chapter measures the longer-term effect of the CES on the local economy, using night-time light intensity measured from space, and focus on the role of insurance payments for damaged residential property during the local recovery process. Uniquely for this event, more than 95% of residential housing units were covered by insurance and almost all incurred some damage. However, insurance payments were staggered over 5 years, enabling us to identify their local impact. We find that night-time luminosity can capture the process of recovery; and that insurance payments contributed significantly to the process of local economic recovery after the earthquake. Yet, delayed payments were less affective in assisting recovery and cash settlement of claims were more effective than insurance-managed repairs. After the Christchurch earthquakes, the government declared about 8000 houses as Red Zoned, prohibiting further developments in these properties, and offering the owners to buy them out. The government provided two options for owners: the first was full payment for both land and dwelling at the 2007 property evaluation, the second was payment for land, and the rest to be paid by the owner‘s insurance. Most people chose the second option. Using data from LINZ combined with data from Stats NZ, the fourth chapter empirically investigates what led people to choose this second option, and how peer effect influenced the homeowners‘ choices. Due to climate change, public disclosure of coastal hazard information through maps and property reports have been used more frequently by local government. This is expected to raise awareness about disaster risks in local community and help potential property owners to make informed locational decision. However, media outlets and business sector argue that public hazard disclosure will cause a negative effect on property value. Despite this opposition, some district councils in New Zealand have attempted to implement improved disclosure. Kapiti Coast district in the Wellington region serves as a case study for this research. In the fifth chapter, we utilize the residential property sale data and coastal hazard maps from the local district council. This study employs a difference-in-difference hedonic property price approach to examine the effect of hazard disclosure on coastal property values. We also apply spatial hedonic regression methods, controlling for coastal amenities, as our robustness check. Our findings suggest that hazard designation has a statistically and economically insignificant impact on property values. Overall, the risk perception about coastal hazards should be more emphasized in communities.
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.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.