Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.
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.
Observations of out-of-plane (OOP) instability in the 2010 Chile earthquake and in the 2011 Christchurch earthquake resulted in concerns about the current design provisions of structural walls. This mode of failure was previously observed in the experimental response of some wall specimens subjected to in-plane loading. Therefore, the postulations proposed for prediction of the limit states corresponding to OOP instability of rectangular walls are generally based on stability analysis under in-plane loading only. These approaches address stability of a cracked wall section when subjected to compression, thereby considering the level of residual strain developed in the reinforcement as the parameter that prevents timely crack closure of the wall section and induces stability failure. The New Zealand code requirements addressing the OOP instability of structural walls are based on the assumptions used in the literature and the analytical methods proposed for mathematical determination of the critical strain values. In this study, a parametric study is conducted using a numerical model capable of simulating OOP instability of rectangular walls to evaluate sensitivity of the OOP response of rectangular walls to variation of different parameters identified to be governing this failure mechanism. The effects of wall slenderness (unsupported height-to-thickness) ratio, longitudinal reinforcement ratio of the boundary regions and length on the OOP response of walls are evaluated. A clear trend was observed regarding the influence of these parameters on the initiation of OOP displacement, based on which simple equations are proposed for prediction of OOP instability in rectangular walls.
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.
This study analyses the success and limitations of the recovery process following the 2010–11 earthquake sequence in Christchurch, New Zealand. Data were obtained from in-depth interviews with 32 relocated households in Christchurch, and from a review of recovery policies implemented by the government. A top-down approach to disaster recovery was evident, with the creation of multiple government agencies and processes that made grassroots input into decision-making difficult. Although insurance proceeds enabled the repair and rebuilding of many dwellings, the complexity and adversarial nature of the claim procedures also impaired recovery. Householders’ perceptions of recovery reflected key aspects of their post-earthquake experiences (e.g. the housing offer they received, and the negotiations involved), and the outcomes of their relocation (including the value of the new home, their subjective well-being, and lifestyle after relocation). Protracted insurance negotiations, unfair offers and hardships in post-earthquake life were major challenges to recovery. Less-thanfavourable recovery experiences also transformed patterns of trust in local communities, as relocated householders came to doubt both the government and private insurance companies’ ability to successfully manage a disaster. At the same time, many relocated households expressed trust in their neighbours and communities. This study illuminates how government policies influence disaster recovery while also suggesting a need to reconsider centralised, top-down approaches to managing recovery.
Abstract. Natural (e.g., earthquake, flood, wildfires) and human-made (e.g., terrorism, civil strife) disasters are inevitable, can cause extensive disruption, and produce chronic and disabling psychological injuries leading to formal diagnoses (e.g., post-traumatic stress disorder [PTSD]). Following natural disasters of earthquake (Christchurch, Aotearoa/New Zealand, 2010–11) and flood (Calgary, Canada, 2013), controlled research showed statistically and clinically significant reductions in psychological distress for survivors who consumed minerals and vitamins (micronutrients) in the following months. Following a mass shooting in Christchurch (March 15, 2019), where a gunman entered mosques during Friday prayers and killed and injured many people, micronutrients were offered to survivors as a clinical service based on translational science principles and adapted to be culturally appropriate. In this first translational science study in the area of nutrition and disasters, clinical results were reported for 24 clients who completed the Impact of Event Scale – Revised (IES-R), the Depression Anxiety Stress Scales (DASS), and the Modified-Clinical Global Impression (M-CGI-I). The findings clearly replicated prior controlled research. The IES-R Cohen’s d ESs were 1.1 (earthquake), 1.2 (flood), and 1.13 (massacre). Effect sizes (ESs) for the DASS subscales were also consistently positive across all three events. The M-CGI-I identified 58% of the survivors as “responders” (i.e., self-reported as “much” to “very much” improved), in line with those reported in the earthquake (42%) and flood (57%) randomized controlled trials, and PTSD risk reduced from 75% to 17%. Given ease of use and large ESs, this evidence supports the routine use of micronutrients by disaster survivors as part of governmental response.
This study analyses the Earthquake Commission’s (EQC) insurance claims database to investigate the influence of seismic intensity and property damage resulting from the Canterbury Earthquake Sequence (CES) on the repair costs and claim settlement duration for residential buildings. Firstly, the ratio of building repair cost to its replacement cost was expressed as a Building Loss Ratio (BLR), which was further extended to Regional Loss Ratio (RLR) for greater Christchurch by multiplying the average of all building loss ratios with the proportion of building stock that lodged an insurance claim. Secondly, the total time required to settle the claim and the time taken to complete each phase of the claim settlement process were obtained. Based on the database, the regional loss ratio for greater Christchurch for three events producing shakings of intensities 6, 7, and 8 on the modified Mercalli intensity scale were 0.013, 0.066, and 0.171, respectively. Furthermore, small (less than NZD15,000), medium (between NZD15,000 and NZD100,000), and large (more than NZD100,000) claims took 0.35-0.55, 1.95-2.45, and 3.35-3.85 years to settle regardless of the building’s construction period and earthquake intensities. The number of claims was also disaggregated by various building characteristics to evaluate their relative contribution to the damage and repair costs.
Low Damage Seismic Design (LDSD) guidance material being developed by Engineering NZ is considering a design drift limit for multi-storey buildings of 0.5% at a new damage control limit state (DCLS). The impact of this new design requirement on the expected annual loss due to repair costs is investigated for a four-storey office building with reinforced concrete walls located in Christchurch. The LDSD guidance material aims to reduce the expected annual loss of complying buildings to below 0.1% of building replacement cost. The research tested this expectation. Losses were estimated in accordance with FEMA P58, using building responses from non-linear time history analyses (performed with OpenSees using lumped plasticity models). The equivalent static method, in line with NZS 1170.5 and NZS 3101, was used to design the building to LDSD specifications, representing a future state-of-practice design. The building designed to low-damage specification returned an expected annual loss of 0.10%, and the building designed conventionally returned an expected annual loss of 0.13%. Limitations with the NZS 3101 method for determining wall stiffness were identified, and a different method acknowledging the relationship between strength and stiffness was used to redesign the building. Along with improving this design assumption, the study finds that LDSD design criteria could be an effective way of limiting damage and losses.
To reduce seismic vulnerability and the economic impact of seismic structural damage, it is important to protect structures using supplemental energy dissipation devices. Several types of supplemental damping systems can limit loads transferred to structures and absorb significant response energy without sacrificial structural damage. Lead extrusion dampers are one type of supplemental energy dissipation devices. A smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, have been employed in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch and San Francisco. HF2V devices have previously been designed using very simple models with limited precision. They are then manufactured, and tested to ensure force capacities match design goals, potentially necessitating reassembly or redesign if there is large error. In particular, devices with a force capacity well above or below a design range can require more testing and redesign, leading to increased economic and time cost. Thus, there is a major need for a modelling methodology to accurately estimate the range of possible device force capacity values in the design phase – upper and lower bounds. Upper and lower bound force capacity estimates are developed from equations in the metal extrusion literature. These equations consider both friction and extrusion forces between the lead and the bulged shaft in HF2V devices. The equations for the lower and upper bounds are strictly functions of device design parameters ensuring easy use in the design phase. Two different sets of estimates are created, leading to estimates for the lower and upper bounds denoted FLB,1, FUB,1, FUB,2, respectively. The models are validated by comparing the bounds with experimental force capacity data from 15 experimental HF2V device tests. All lower bound estimates are below or almost equal to the experimental device forces, and all upper bound estimates are above. Per the derivation, the (FLB,1, FUB,1) pair provide narrower bounds. The (FLB,1, FUB,1) pair also had a mean lower bound gap of -34%, meaning the lower bound was 74% of device force on average, while the mean upper bound gap for FUB,1 was +23%. These are relatively tight bounds, within ~±2 SE of device manufacture, and can be used as a guide to ensure device forces are in range for the actual design use when manufactured. Therefore, they provide a useful design tool.
Climate change and population growth will increase vulnerability to natural and human-made disasters or pandemics. Longitudinal research studies may be adversely impacted by a lack of access to study resources, inability to travel around the urban environment, reluctance of sample members to attend appointments, sample members moving residence and potentially also the destruction of research facilities. One of the key advantages of longitudinal research is the ability to assess associations between exposures and outcomes by limiting the influence of sample selection bias. However, ensuring the validity and reliability of findings in longitudinal research requires the recruitment and retention of respondents who are willing and able to be repeatedly assessed over an extended period of time. This study examined recruitment and retention strategies of 11 longitudinal cohort studies operating during the Christchurch, New Zealand earthquake sequence which began in September 2010, including staff perceptions of the major impediments to study operations during/after the earthquakes and respondents’ barriers to participation. Successful strategies to assist recruitment and retention after a natural disaster are discussed. With the current COVID-19 pandemic, longitudinal studies are potentially encountering some of the issues highlighted in this paper including: closure of facilities, restricted movement of research staff and sample members, and reluctance of sample members to attend appointments. It is possible that suggestions in this paper may be implemented so that longitudinal studies can protect the operation of their research programmes.<br /><br />Key messages<br /><ul><li>Recruitment and retention of longitudinal study participants is challenging following a natural disaster.</li><br /><li>The long-lasting, global effects of the Covid 19 pandemic will increase this problem.</li><br /><li>Longitudinal study researchers should develop protocols to support retention before a disaster occurs.</li><br /><li>Researchers need to be pragmatic and flexible in the design and implementation of their studies.</li></ul>
Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.
While it is well known that challenging and distressing events can negatively impact people’s psychological and physical state, increasingly researchers have investigated how challenging or stressful life circumstances can lead to the phenomenon of posttraumatic growth: positive psychological or life changes that can emerge from potentially traumatic events. Posttraumatic growth has been investigated primarily with people displaying varying levels of posttraumatic stress symptoms and other psychopathology due to theories suggesting that resilience would prohibit posttraumatic growth. Few studies have examined growth amongst resilient people. The current study examined posttraumatic growth in a sample of sixty psychologically healthy people who experienced the Canterbury earthquake sequence of 2010-2011. The current study is a follow-up study that used thematic analysis to explore: (1) Whether posttraumatic growth is evident nine years after the Canterbury earthquake sequence and approximately six years after baseline assessment; and (2) What themes may facilitate the posttraumatic growth process in psychologically healthy people. Data were collected using semi-structured interviews. Thematic analysis revealed four themes describing participants’ experiences of growth: New possibilities, reappraisal of life and priorities, positive changes in self-perception and closer more meaningful relationships. Themes describing posttraumatic growth provide evidence for research question one. Thematic analysis revealed three main themes and multiple subthemes that may facilitate the process of growth in psychologically healthy people: Hardship, optimistic positive appraisal and people helping people. Themes describing processes that may lead to growth provide evidence for research question two. Results of the current study provide insights about the experience of growth in psychologically healthy people and cognitive and psychosocial factors that may facilitate growth in resilient individuals.
Though rare and unpredictable, earthquakes can and do cause catastrophic destruction when they impact unprepared and vulnerable communities. Extensive damage and failure of vulnerable buildings is a key factor which contributes to seismic-related disasters, making the proactive management of these buildings a necessity to reduce the risk of future disasters arising. The devastating Canterbury earthquakes of 2010 and 2011 brought the urgency of this issue to national importance in New Zealand. The national earthquake-prone building framework came into effect in 2017, obligating authorities to identify existing buildings with the greatest risk of collapse in strong earthquakes and for building owners to strengthen or demolish these buildings within a designated period of time. Though this framework is unique to New Zealand, the challenge of managing the seismic risk of such buildings is common amongst all seismically-active countries. Therefore, looking outward to examine how other jurisdictions legally manage this challenge is useful for reflecting on the approaches taken in New Zealand and understand potential lessons which could be adopted. This research compares the legal framework used to reduce the seismic risk of existing buildings in New Zealand with that of the similarly earthquake-prone countries of Japan and Italy. These legal frameworks are examined with a particular focus on the proactive goal of reducing risk and improving resilience, as is the goal of the international Sendai Framework for Disaster Risk Reduction 2015-2030. The Sendai Framework, which each of the case study countries have committed to and thus have obligations under, forms the legal basis of the need for states to reduce disaster risk in their jurisdictions. In particular, the states’ legal frameworks for existing building risk reduction are examined in the context of the Sendai priorities of understanding disaster risk, strengthening disaster risk governance, and investing in resilience. While this research illustrates that the case study countries have each adopted more proactive risk reduction frameworks in recent years in anticipation of future earthquakes, the frameworks currently focus on a very narrow range of existing buildings and thus are not currently sufficient for promoting the long-term resilience of building stocks. In order to improve resilience, it is argued, legal frameworks need to include a broader range of buildings subject to seismic risk reduction obligations and also to broaden the focus on long-term monitoring of potential risk to buildings.
A buckling-restrained braced frame (BRBF) is a structural bracing system that provides lateral strength and stiffness to buildings and bridges. They were first developed in Japan in the 1970s (Watanabe et al. 1973, Kimura et al. 1976) and gained rapid acceptance in the United States after the Northridge earthquake in 1994 (Bruneau et al. 2011). However, it was not until the Canterbury earthquakes of 2010/2011, that the New Zealand construction market saw a significant uptake in the use of buckling-restrained braces (BRBs) in commercial buildings (MacRae et al. 2015). In New Zealand there is not yet any documented guidance or specific instructions in regulatory standards for the design of BRBFs. This makes it difficult for engineers to anticipate all the possible stability and strength issues within a BRBF system and actively mitigate them in each design. To help ensure BRBF designs perform as intended, a peer review with physical testing are needed to gain building compliance in New Zealand. Physical testing should check the manufacturing and design of each BRB (prequalification testing), and the global strength and stability of each BRB its frame (subassemblage testing). However, the financial pressures inherent in commercial projects has led to prequalification testing (BRB only testing) being favoured without adequate design specific subassemblage testing. This means peer reviewers have to rely on BRB suppliers for assurances. This low regulation environment allows for a variety of BRBF designs to be constructed without being tested or well understood. The concern is that there may be designs that pose risk and that issues are being overlooked in design and review. To improve the safety and design of BRBFs in New Zealand, this dissertation studies the behaviour of BRBs and how they interact with other frame components. Presented is the experimental test process and results of five commercially available BRB designs (Chapter 2). It discusses the manufacturing process, testing conditions and limitations of observable information. It also emphasises that even though subassemblage testing is impractical, uniaxial testing of the BRB only is not enough, as this does not check global strength or stability. As an alternative to physical testing, this research uses computer simulation to model BRB behaviour. To overcome the traditional challenges of detailed BRB modelling, a strategy to simulate the performance of generic BRB designs was developed (Chapter 3). The development of nonlinear material and contact models are important aspects of this strategy. The Chaboche method is employed using a minimum of six backstress curves to characterize the combined isotropic and kinematic hardening exhibited by the steel core. A simplified approach, adequate for modelling the contact interaction between the restrainer and the core was found. Models also capture important frictional dissipation as well as lateral motion and bending associated with high order constrained buckling of the core. The experimental data from Chapter 2 was used to validate this strategy. As BRBs resist high compressive loading, global stability of the BRB and gusseted connection zone need to be considered. A separate study was conducted that investigated the yielding and buckling strength of gusset plates (Chapter 4). The stress distribution through a gusset plate is complex and difficult to predict because the cross-sectional area of gusset plate is not uniform, and each gusset plate design is unique in shape and size. This has motivated design methods that approximate yielding of gusset plates. Finite element modelling was used to study the development of yielding, buckling and plastic collapse behaviour of a brace end bolted to a series of corner gusset plates. In total 184 variations of gusset plate geometries were modelled in Abaqus®. The FEA modelling applied monotonic uniaxial load with an imperfection. Upon comparing results to current gusset plate design methods, it was found that the Whitmore width method for calculating the yield load of a gusset is generally un-conservative. To improve accuracy and safety in the design of gusset plates, modifications to current design methods for calculating the yield area and compressive strength for gusset plates is proposed. Bolted connections are a popular and common connection type used in BRBF design. Global out-of-plane stability tends to govern the design for this connection type with numerous studies highlighting the risk of instability initiated by inelasticity in the gussets, neck of the BRB end and/or restrainer ends. Subassemblage testing is the traditional method for evaluating global stability. However, physical testing of every BRBF variation is cost prohibitive. As such, Japan has developed an analytical approach to evaluate out-of-plane stability of BRBFs and incorporated this in their design codes. This analytical approach evaluates the different BRB components under possible collapse mechanisms by focusing on moment transfer between the restrainer and end of the BRB. The approach have led to strict criteria for BRBF design in Japan. Structural building design codes in New Zealand, Europe and the United States do not yet provide analytical methods to assess BRB and connection stability, with prototype/subassemblage testing still required as the primary means of accreditation. Therefore it is of interest to investigate the capability of this method to evaluate stability of BRBs designs and gusset plate designs used in New Zealand (including unstiffened gusset connection zones). Chapter 5 demonstrates the capability of FEA to study to the performance of a subassemblage test under cyclic loading – resembling that of a diagonal ground storey BRBF with bolted connections. A series of detailed models were developed using the strategy presented in Chapter 3. The geometric features of BRB 6.5a (Chapter 2) were used as a basis for the BRBs modelled. To capture the different failure mechanisms identified in Takeuchi et al. (2017), models varied the length that the cruciform (non-yielding) section inserts into the restrainer. Results indicate that gusset plates designed according to New Zealand’s Steel Structures Standard (NZS 3404) limit BRBF performance. Increasing the thickness of the gusset plates according to modifications discussed in Chapter 4, improved the overall performance for all variants (except when Lin/ Bcruc = 0.5). The effect of bi-directional loading was not found to notably affect out-of-plane stability. Results were compared against predictions made by the analytical method used in Japan (Takeuchi method). This method was found to be generally conservative is predicting out-of-plane stability of each BRBF model. Recommendations to improve the accuracy of Takeuchi’s method are also provided. The outcomes from this thesis should be helpful for BRB manufacturers, researchers, and in the development of further design guidance of BRBFs.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
A building boom in the 1980s allowed pre-stressed hollow-core floor construction to be widely adopted in New Zealand, even though the behaviour of these prefabricated elements within buildings was still uncertain. Inspections following the Canterbury and Kaikōura earthquakes has provided evidence of web-splitting, transverse cracking and longitudinal splitting on hollow-core units, confirming the susceptibility of these floors to undesirable failure modes. Hollow-core slabs are mainly designed to resist bending and shear. However, there are many applications in which they are also subjected to torsion. In New Zealand, hollow-core units contain no transverse reinforcement in the soffit concrete below the cells and no web reinforcement. Consequently, their dependable performance in torsion is limited to actions that they can resist before torsional cracking occurs. In previous work by the present authors, a three-dimensional FE modelling approach to study the shear flexural behaviour of precast pre-stressed hollow core units was developed and validated by full-scale experiments. This paper shows how the FE analyses have been extended to investigate the response of HC units subjected to torsional actions. Constitutive models, based on nonlinear fracture mechanics, have been used to numerically predict the torsional capacity of HC units and have been compared with experimental results. The results indicate that the numerical approach is promising and should be developed further as part of future research.
Major earthquakes, such as the Canterbury and Kaikoura events recorded in New Zealand in 2010 and 2016 respectively, highlighted that floor systems can be heavily damaged. At a reduced or full scale, quasi-static experimental tests on structural sub-assemblies can help to establish the seismic performance of structural systems. However, the experimental performance obtained with such tests is likely to be dependent on the drift protocol adopted. This paper provides an overview of the drift protocols which have been assumed in previous relevant experimental activities, with emphasis on those adopted for testing floor systems. The paper also describes the procedure used to define the loading protocol applied in the testing of a large precast concrete floor diaphragm as part of the Recast floor project at the University of Canterbury. Finally, major limits of current loading protocols, and areas of future research, are identified.
Floor systems with precast concrete hollow-core units have been largely used in concrete buildings built in New Zealand during the 1980’s. Recent earthquakes, such as the Canterbury sequence in 2010-2011 and the Kaikoura earthquake in 2016, highlighted that this floor system can be highly vulnerable and potentially lead to the floor collapse. A series of research activities are in progress to better understand the seismic performance of floor diaphragms, and this research focuses on examining the performance of hollow core units running parallel to the walls of wall-resisting concrete structures. This study first focused on the development of fragility functions, which can be quickly used to assess likelihood of the hollow-core being able to survive given the buildings design drift, and secondly to determine the expected performance of hollow-core units that run parallel to walls, focusing on the alpha unit running by the wall. Fragility functions are created for a range of different parameters for both vertical dislocation and crack width that can be used as the basis of a quick analysis or loss estimation for the likely impact of hollow-core floors on building vulnerability and risk. This was done using past experimental tests, and the recorded damage. Using these results and the method developed by Baker fragility curves were able to be created for varying crack widths and vertical dislocations. Current guidelines for analysis of hollow-core unit incompatible displacements are based on experimental vertical displacement results from concrete moment resisting frame systems to determine the capacity of hollow-core elements. To investigate the demands on hollow-core units in a wall-based structure, a fibre-element model in the software Seismostruct is created and subject to quasi-static cyclic loading, using elements which are verified from previous experimental tests. It is shown that for hollow-core units running by walls that the 10 mm displacement capacity used for hollow-core units running by a beam is insufficient for members running by walls and that shear analysis should be used. The fibre-element model is used to simulate the seismic demand induced on the floor system and has shown that the shear demand is a function of drift, wall length, hollow-core span, linking slab length and, to a minor extent, wall elongation.
This study investigates the uncertainty of simulated earthquake ground motions for smallmagnitude events (Mw 3.5 – 5) in Canterbury, New Zealand. 148 events were simulated with specified uncertainties in: event magnitude, hypocentre location, focal mechanism, high frequency rupture velocity, Brune stress parameter, the site 30-m time-averaged shear wave velocity (Vs30), anelastic attenuation (Q) and high frequency path duration. In order to capture these uncertainties, 25 realisations for each event were generated using the Graves and Pitarka (2015) hybrid broadband simulation approach. Monte-Carlo realisations were drawn from distributions for each uncertainty, to generate a suite of simulation realisations for each event and site. The fit of the multiple simulation realisations to observations were assessed using linear mixed effects regression to generate the systematic source, path and site effects components across all ground motion intensity measure residuals. Findings show that additional uncertainties are required in each of the three source, path, and site components, however the level of output uncertainty is promising considering the input uncertainties included.
Glazing systems are non-structural elements in a building that, more often than not, appear to be given little consideration in seismic design. Recent experimental work into glazing systems at the University of Canterbury, however, has shown that glazing systems can be very susceptible to serviceability damage, defined as loss of water-tightness. The focus of this paper is to highlight the difference in vulnerability of standard and seismic glazing systems and consider the implications of this for future repair costs and losses. The paper first describes the damage states chosen for glazing units according to the repair strategies required and expected repair costs. This includes three damage states: DS1: Water Leakage, DS2: Gasket Failure and DS3: Frame/Glass Failure. Implementing modern performance-based earthquake engineering, the paper proceeds to highlight a case study comparing costs and expected losses of a standard glazing unit and a seismic glazing unit installed on a case study building. It is shown that the use of seismic glazing units is generally beneficial over time, due to the early onset of serviceability damage in standard glazing units. Finally, the paper provides suggestions for designers aimed at reducing costs related to earthquake induced repairs of glazing.
Researchers have begun to explore the opportunity presented by blue-green infrastructure(a subset of nature-based solutions that provide blue and green space in urban infrastructure)as a response to the pressures of climate change. The 2010/2011 Canterbury earthquake sequence created a unique landscape within which there is opportunity to experiment with and invest in new solutions to climate change adaptation in urban centres. Constructed wetlands are an example of blue-green infrastructure that can potentially support resilience in urban communities. This research explores interactions between communities and constructed wetlands to understand how this may influence perceptions of community resilience. The regeneration of the Ōtākaro Avon River Corridor (OARC) provides a space to investigate these relationships. Seven stakeholders from the community, industry, and academia, each with experience in blue-green infrastructure in the OARC, participated in a series of semi-structured interviews. Each participant was given the opportunity to reflect on their perspectives of community, community resilience, and constructed wetlands and their interconnections. Interview questions aligned with the overarching research objectives to (1) understand perceptions around the role of wetlands in urban communities, (2) develop a definition for community resilience in the context of the Ōtākaro Avon community, and (3) reflect on how wetlands can contribute to (or detract from) community resilience. This study found that constructed wetlands can facilitate learning about the challenges and solutions needed to adapt to climate change. From the perspective of the community representatives, community resilience is linked to social capital. Strong social networks and a relationship with nature were emphasised as core components of a community’s ability to adapt to disruption. Constructed wetlands are therefore recognised as potentially contributing to community resilience by providing spaces for people to engage with each other and nature. Investment in constructed wetlands can support a wider response to climate change impacts. This research was undertaken with the support of the Ōtākaro Living Laboratory Trust, who are invested in the future of the OARC. The outcomes of this study suggest that there is an opportunity to use wetland spaces to establish programmes that explore the perceptions of constructed wetlands from a broader community definition, at each stage of the wetland life cycle, and at wider scales(e.g., at a city scale or beyond).
There is growing expectation that local volunteers will play a more integrated role in disaster response, yet emergent groups are often ‘outsiders’ to crisis management, prompting questions of the conditions and processes by which these groups can forge relationships with established response agencies, and the tensions which can arise those interactions. This article analyses how student-led volunteers, as an emergent group, nevertheless gained “authority to operate” in the aftermath of the 2010-2011 earthquakes in Canterbury, New Zealand. Our study demonstrates how established response agencies and emergent groups can form hugely impactful and mutually supportive relationships. However, our analysis also points to two interrelated tensions that can arise, regarding the terms by which emergent groups are recognised, and the ‘distance’ considered necessary between emergent groups and established response agencies. The discussion considers implications for inclusiveness, risk and responsibility if emergent volunteers are to be further integrated into disaster response.
Brooklands Lagoon / Te Riu o Te Aika Kawa (‘Brooklands’) is an important wetland and estuarine ecosystem in Canterbury. It is a site of cultural significance to Ngāi Tūāhuriri, and is also valued by the wider community. Home to an array of life, it is connected to the Pūharakekenui/Styx and Waimakariri rivers, and is part of a wetland landscape complex that includes the Avon-Heathcote / Ihutai estuary to the south and the Ashley / Rakahuri estuary to the north. Notionally situated within the territorial boundary of Christchurch City Council and jurisdictionally encompassed by the regional council Environment Canterbury, it has been legally determined to be part of the coastal marine area. The complicated administrative arrangements for the lagoon mirror the biophysical and human challenges to this surprisingly young ecosystem since its formation in 1940.
Here we present a synthesis of the historical events and environmental influences that have shaped Brooklands Lagoon. Before existing as an intertidal ecosystem, the Waimakariri river mouth was situated in what is now the southern end of the lagoon. A summary timeline of key events is set out in the table below. These included the diversion of the Waimakariri River mouth via the construction of Wrights Cut in the 1930s, which influenced the way that the lower reaches of the river interacted with the land and sea. A large flood in 1940 shifted the river mouth ~2 to 3 kilometres north, that created the landscape that we see today. However, this has not remained stable, as the earthquake sequence in 2010 and 2011 subsided the bed of the estuary.
The changes are ongoing, as sea level rise and coastal inundation will place ongoing pressure on the aquatic ecosystem and surrounding land. How to provide accommodation space for Brooklands as an estuary will be a key planning and community challenge, as Environment Canterbury begins the engagement for the review of its Regional Coastal Plan. There is also a requirement to safeguard its ecological health under the 2020 National Policy Statement on Freshwater Management. This will necessitate an integrated mountains to sea (ki uta ki tai) management approach as the lagoon is affected by wider catchment activities. We hope that this report will contribute to, and inform these processes by providing a comprehensive historical synthesis, and by identifying considerations for the future collaborative management of Brooklands Lagoon, and protection of its values. In essence, we suggest that Te Riu o Te Aika Kawa deserves some sustained aroha.
<b>Construction and Demolition (C&D) waste contributes to over 50% of New Zealand’s overall waste. Materials such as timber, plasterboard, and concrete make up 81% of the C&D waste that goes into landfills each year. Alongside this, more than 235 heritage-listed buildings have been demolished in Christchurch since the 2011 earthquakes. This research portfolio aims to find a solution to decrease C&D waste produced by demolishing heritage buildings.</b>
With the recent announcement of The Cathedral of the Blessed Sacrament’s demolition, this will be another building added to the list of lost heritage in Christchurch. This research portfolio aims to bridge the relationship between heritage and waste through the recycling and reuse of the demolished materials, exploring the idea that history and heritage are preserved through building material reuse.
This research portfolio mainly focuses on reducing construction and demolition waste in New Zealand, using the design of a new Catholic Cathedral as a vessel. This thesis will challenge how the construction and design industry deals with the demolition of heritage buildings and their contribution to New Zealand’s waste. It aims to explore the idea of building material reuse not only to reduce waste but also to retain the history and heritage of the demolished building within the materials.