Designing a structure for higher- than-code seismic performance can result in significant economic and environmental benefits. This higher performance can be achieved using the principles of Performance-Based Design, in which engineers design structures to minimize the probabilistic lifecycle seismic impacts on a building. Although the concept of Performance-Based Design is not particularly new, the initial capital costs associated with designing structures for higher performance have historically hindered the widespread adoption of performance-based design practices. To overcome this roadblock, this research is focused on providing policy makers and stakeholders with evidence-based environmental incentives for designing structures in New Zealand for higher seismic performance. In the first phase of the research, the environmental impacts of demolitions in Christchurch following the Canterbury Earthquakes were quantified to demonstrate the environmental consequences of demolitions following seismic events. That is the focus here. A building data set consisting of 142 concrete buildings that were demolished following the earthquake was used to quantify the environmental impacts of the demolitions in terms of the embodied carbon and energy in the building materials. A reduced set of buildings was used to develop a material takeoff model to estimate material quantities in the entire building set, and a lifecycle assessment tool was used to calculate the embodied carbon and energy in the materials. The results revealed staggering impacts in terms of the embodied carbon and energy in the materials in the demolished buildings. Ongoing work is focused developing an environmental impact framework that incorporates all the complex factors (e.g. construction methodologies, repair methodologies (if applicable), demolition methodologies (if applicable), and waste management) that contribute to the environmental impacts of building repair and demolition following earthquakes.
Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.
Despite the relatively low seismicity, a large earthquake in the Waikato region is expected to have a high impact, when the fourth-largest regional population and economy and the high density critical infrastructure systems in this region are considered. Furthermore, Waikato has a deep soft sedimentary basin, which increases the regional seismic hazard due to trapping and amplification of seismic waves and generation of localized surface waves within the basin. This phenomenon is known as the “Basin Effect”, and has been attributed to the increased damage in several historic earthquakes, including the 2010-2011 Canterbury earthquakes. In order to quantitatively model the basin response and improve the understanding of regional seismic hazard, geophysical methods will be used to develop shear wave velocity profiles across the Waikato basin. Active surface wave methods involve the deployment of linear arrays of geophones to record the surface waves generated by a sledge hammer. Passive surface wave methods involve the deployment of two-dimensional seismometer arrays to record ambient vibrations. At each site, the planned testing includes one active test and two to four passive arrays. The obtained data are processed to develop dispersion curves, which describe surface wave propagation velocity as a function of frequency (or wavelength). Dispersion curves are then inverted using the Geopsy software package to develop a suite of shear wave velocity profiles. Currently, more than ten sites in Waikato are under consideration for this project. This poster presents the preliminary results from the two sites that have been tested. The shear wave velocity profiles from all sites will be used to produce a 3D velocity model for the Waikato basin, a part of QuakeCoRE flagship programme 1.
In recent Canterbury earthquakes, structures have performed well in terms of life safety but the estimated total cost of the rebuild was as high as $40 billion. The major contributors to this cost are repair/demolition/rebuild cost, the resulting downtime and business interruption. For this reason, the authors are exploring alternate building systems that can minimize the downtime and business interruption due to building damage in an earthquake; thereby greatly reducing the financial implications of seismic events. In this paper, a sustainable and demountable precast reinforced concrete (RC) frame system in which the precast members are connected via steel tubes/plates or steel angles/plates and high strength friction grip (HSFG) bolts is introduced. In the proposed system, damaged structural elements in seismic frames can be easily replaced with new ones; thereby making it an easily and quickly repairable and a low-loss system. The column to foundation connection in the proposed system can be designed either as fixed or pinned depending on the requirement of strength and stiffness. In a fixed base frame system, ground storey columns will also be damaged along with beams in seismic events, which are to be replaced after seismic events; whereas in a pin base frame only beams (which are easy to replace) will be damaged. Low to medium rise (3-6 storey) precast RC frame buildings with fixed and pin bases are analyzed in this paper; and their lateral capacity, lateral stiffness and natural period are scrutinized to better understand the pros and cons of the demountable precast frame system with fixed and pin base connections.
The development of cheap, whilst effective and relatively non-invasive structural retrofit techniques for existing non-ductile reinforced concrete (RC) structures still remains the most challenging issue for a wide implementation on a macro scale. Seismic retrofit is too often being confused as purely structural strengthening. As part of a six-years national project on “Seismic retrofit solutions for NZ multi-storey building”, focus has been given at the University of Canterbury on the development of a counter-intuitive retrofit strategy for earthquake vulnerable existing rc frame, based on a “selective weakening” (SW) approach. After an overview of the SW concept, this paper presents the experimental and numerical validation of a SW retrofit strategy for earthquake vulnerable existing RC frame with particular focus on the exterior beam-column (b-c) joints. The exterior b-c joint is a critically vulnerable region in many existing pre-1970s RC frames. By selectively weakening the beam by cutting the bottom longitudinal reinforcements and/or adding external pre-stressing to the b-c joint, a more desirable inelastic mechanism can be attained, leading to improved global seismic performance. The so-called SW retrofit is implemented on four 2/3-scaled exterior RC b-c joint subassemblies, tested under quasi-static cyclic loading at the University of Canterbury. Complemented by refined 3D Finite Element (FE) models and dynamic time-history analyses results, the experimental results have shown the potential of a simple and cost-effective yet structurally efficient structural rehabilitation technique. The research also demonstrated the potential of advanced 3D fracture-mechanics-based microplane concrete modelling for refined FE analysis of non-ductile RC b-c joints.
Light timber framed (LTF) structures provide a cost-effective and structurally efficient solution for low-rise residential buildings. This paper studies seismic performance of single-storey LTF buildings sheathed by gypsum-plasterboards (GPBs) that are a typical lining product in New Zealand houses. Compared with wood-based structural panels, GPBs tend to be more susceptible to damage when they are used in bracing walls to resist earthquake loads. This study aims to provide insights on how the bracing wall irregularity allowed by the current New Zealand standard NZS 3604 and the in-plane rigidity of ceiling diaphragms affect the overall seismic performance of these GPB-braced LTF buildings. Nonlinear time-history analyses were conducted on a series of single-storey baseline buildings with different levels of bracing wall irregularities and ceiling diaphragm rigidity. The results showed significant torsional effect caused by the eccentric bracing wall layout with semi-rigid/rigid ceiling diaphragms. On average, bracing wall drift demand caused by the extreme bracing wall irregularities was three times of that in the regular bracing wall layout under the rigid diaphragm assumption. This finding agreed well with the house survey after the 2011 Canterbury Earthquake in which significantly more damage was observed in the houses with irregular bracing wall layouts and relatively rigid diaphragms. Therefore, it is recommended to limit the level of bracing wall eccentricity and ensure the sufficiently rigid diaphragms to avoid excessive damage in these LTF buildings in future events.
Surface rupture of the previously unrecognised Greendale Fault extended west-east for ~30 km across alluvial plains west of Christchurch, New Zealand, during the Mw 7.1 Darfield (Canterbury) earthquake of September 2010. Surface rupture displacement was predominantly dextral strike-slip, averaging ~2.5 m, with maxima of ~5 m. Vertical displacement was generally less than 0.75 m. The surface rupture deformation zone ranged in width from ~30 to 300 m, and comprised discrete shears, localised bulges and, primarily, horizontal dextral flexure. About a dozen buildings, mainly single-storey houses and farm sheds, were affected by surface rupture, but none collapsed, largely because most of the buildings were relatively flexible and resilient timber-framed structures and also because deformation was distributed over a relatively wide zone. There were, however, notable differences in the respective performances of the buildings. Houses with only lightly-reinforced concrete slab foundations suffered moderate to severe structural and non-structural damage. Three other buildings performed more favourably: one had a robust concrete slab foundation, another had a shallow-seated pile foundation that isolated ground deformation from the superstructure, and the third had a structural system that enabled the house to tilt and rotate as a rigid body. Roads, power lines, underground pipes, and fences were also deformed by surface fault rupture and suffered damage commensurate with the type of feature, its orientation to the fault, and the amount, sense and width of surface rupture deformation.
This poster provides a comparison between the strong ground motions observed in the 22 February 2011 Mw6.3 Christchurch earthquake with those observed in Tokyo during the 11 March 2011 Mw9.0 Tohoku earthquake. The destuction resulting from both of these events has been well documented, although tsunami was the principal cause of damage in the latter event, and less attention has been devoted to the impact of earthquake-induced ground motions. Despite Tokyo being located over 100km from the nearest part of the causative rupture, the ground motions observed from the Tohoku earthquake were significant enough to cause structural damage and also significant liquefaction to loose reclaimed soils in Tokyo Bay. The author was fortunate enough (from the perspective of an earthquake engineer) to experience first-hand both of these events. Following the Tohoku event, the athor conducted various ground motion analyses and reconniassance of the Urayasu region in Tokyo Bay affected by liquefaction in collaboration with Prof. Kenji Ishihara. This conference is therefore a fitting opportunity in which to discuss some of authors insights obtained as a result of this first hand knowledge. Figure 1 illustrates the ground motions recorded in the Christchurch CBD in the 22 February 2011 and 4 September 2010 earthquakes, with that recorded in Tokyo Bay in the 11 March 2011 Tohoku earthquake. It is evident that these three ground motions vary widely in their amplitude and duration. The CBGS ground motion from the 22 February 2011 event has a very large amplitude (nearly 0.6g) and short duration (approx. 10s of intense shaking), as a result of the causal Mw6.3 rupture at short distance (Rrup=4km). The CBGS ground motion from the 4 September 2010 earthquake has a longer duration (approx. 30s of intense shaking), but reduced acceleration amplitude, as a result of the causal Mw7.1 rupture at a short-to-moderate distance (Rrup=14km). Finally, the Urayasu ground motion in Tokyo bay during the 11 March 2011 Tohoku earthquake exhibits an acceleration amplitude similar to the 4 September 2010 CBGS ground motion, but a significantly larger duration (approx 150s of intense shaking). Clearly, these three different ground motions will affect structures and soils in different ways depending on the vibration characteristics of the structures/soil, and the potential for strength and stiffness degradation due to cumulative effects. Figure 2 provides a comparison between the arias intensities of the several ground motion records from the three different events. It can be seen that the arias intensities of the ground motions in the Christchurch CBD from the 22 February 2011 earthquake (which is on average AI=2.5m/s) is approximately twice that from the 4 September 2010 earthquake (average AI≈1.25). This is consistent with a factor of approximately 1.6 obtained by Cubrinovski et al. (2011) using the stress-based (i.e.PGA-MSF) approach of liquefaction triggering. It can also be seen that the arias intensity of the ground motions recorded in Tokyo during the 2011 Tohoku earthquake are larger than ground motions in the Christchurch CBD from the 4 September 2011 earthquake, but smaller than those of the 22 February 2011 earthquake. Based on the arias intensity liquefaction triggering approach it can therefore be concluded that the ground motion severity, in terms of liquefaction potential, for the Tokyo ground motions is between those ground motions in Christchurch CBD from the 4 September 2010 and 22 February 2011 events.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
This paper investigates the effects of variability in source rupture parameters on site-specific physics-based simulated ground motions, ascertained through the systematic analysis of ground motion intensity measures. As a preliminary study, we consider simulations of the 22 February 2011 Christchurch earthquake using the Graves and Pitarka (2015) methodology. The effects of source variability are considered via a sensitivity study in which parameters (hypocentre location, earthquake magnitude, average rupture velocity, fault geometry and the Brune stress parameter) are individually varied by one standard deviation. The sensitivity of simulated ground motion intensity measures are subsequently compared against observational data. The preliminary results from this study indicate that uncertainty in the stress parameter and the rupture velocity have the most significant effect on the high frequency amplitudes. Conversely, magnitude uncertainty was found to be most influential on the spectral acceleration amplitudes at low frequencies. Further work is required to extend this preliminary study to exhaustively consider more events and to include parameter covariance. The ultimate results of this research will assist in the validation of the overall simulation method’s accuracy in capturing various rupture parameters, which is essential for the use of simulated ground motion models in probabilistic seismic hazard analysis.
Geomorphic, structural and chronological data are used to establish the late Quaternary paleoseismicity of the active dextral-oblique Northern Esk Fault in North Canterbury, New Zealand. Detailed field mapping of the preserved c. 35 km of surface traces between the Hurunui River and Ashley Head reveals variations in strike ranging from 005° to 057°. Along with kinematic data collected from fault plane striae and offset geomorphic markers along the length of the fault these variations are used to distinguish six structural subsections of the main trace, four dextral-reverse and two dextral-normal. Displacements of geomorphic markers such as minor streams and ridges are measured using differential GPS and rangefinder equipment to reveal lateral offsets ranging from 3.4 to 23.7 m and vertical offsets ranging from < 1 to 13.5 m. Characteristic single event displacements of c. 5 m and c. 2 m have been calculated for strike-slip and reverse sections respectively. The use of fault scaling relationships reveals an anomalously high displacement to surface rupture length ratio when compared to global data sets. Fault scaling relationships based on width limited ruptures and magnitude probabilities from point measurements of displacement imply earthquake magnitudes of Mw 7.0 to 7.5. Optically Stimulated Luminescence (OSL) ages from displaced Holocene alluvial terraces at the northern extent of the active trace along with OSL and radiocarbon samples of the central sections constrain the timing of the last two surface rupturing events (11.15 ±1.65 and 3.5 ± 2.8 ka) and suggest a recurrence interval of c. 5612 ± 445 years and late Quaternary reverse and dextral slip rates of c. 0.31 mm/yr and 0.82 mm/yr respectively. The results of this study show that the Northern Esk Fault accommodates an important component of the c. 0.7 – 2 mm/yr of unresolved strain across the plate boundary within the North Canterbury region and affirm the Esk Fault as a source of potentially damaging ground shaking in the Canterbury region.
In practice, several competing liquefaction evaluation procedures (LEPs) are used to compute factors of safety against soil liquefaction, often for use within a liquefaction potential index (LPI) framework to assess liquefaction hazard. At present, the influence of the selected LEP on the accuracy of LPI hazard assessment is unknown, and the need for LEP-specific calibrations of the LPI hazard scale has never been thoroughly investigated. Therefore, the aim of this study is to assess the efficacy of three CPT-based LEPs from the literature, operating within the LPI framework, for predicting the severity of liquefaction manifestation. Utilising more than 7000 liquefaction case studies from the 2010–2011 Canterbury (NZ) earthquake sequence, this study found that: (a) the relationship between liquefaction manifestation severity and computed LPI values is LEP-specific; (b) using a calibrated, LEP-specific hazard scale, the performance of the LPI models is essentially equivalent; and (c) the existing LPI framework has inherent limitations, resulting in inconsistent severity predictions against field observations for certain soil profiles, regardless of which LEP is used. It is unlikely that revisions of the LEPs will completely resolve these erroneous assessments. Rather, a revised index which more adequately accounts for the mechanics of liquefaction manifestation is needed.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
An as-built reinforced concrete (RC) frame building designed and constructed according to pre-1970s code design construction practice has been recently tested on the shake table at the University of Canterbury. The specimen, 1/2.5 scaled version of the original prototype, consists of two 3-storey 2-bay asymmetric frames in parallel, one interior and one exterior, jointed together by transverse beams and floor slabs. Following the benchmark test, a retrofit intervention has been proposed to rehabilitate the tested specimen. In this paper, detailed information on the assessment and design of the seismic retrofit procedure using GFRP (glass fibre reinforced polymer) materials is given for the whole frame. Hierarchy of strength and sequence of events (damage mechanisms) in the panel zone region are evaluated using a moment-axial load (M-N) interaction performance domain, according to a performance-based retrofit philosophy. Specific limit states or design objectives are targeted with attention given to both strength and deformation limits. In addition, an innovative retrofit solution using FRP anchor dowels for the corner beam-column joints with slabs is proposed. Finally, in order to provide a practical tool for engineering practice, the retrofit procedure is provided in a step-by step flowchart fashion.
New Zealand's devastating Canterbury earthquakes provided an opportunity to examine the efficacy of existing regulations and policies relevant to seismic strengthening of vulnerable buildings. The mixed-methods approach adopted, comprising both qualitative and quantitative approaches, revealed that some of the provisions in these regulations pose as constraints to appropriate strengthening of earthquake-prone buildings. Those provisions include the current seismic design philosophy, lack of mandatory disclosure of seismic risks and ineffective timeframes for strengthening vulnerable buildings. Recommendations arising from these research findings and implications for pre-disaster mitigation for future earthquake and Canterbury's post-disaster reconstruction suggest: (1) a reappraisal of the requirements for earthquake engineering design and construction, (2) a review and realignment of all regulatory frameworks relevant to earthquake risk mitigation, and (3) the need to develop a national programme necessary to achieve consistent mitigation efforts across the country. These recommendations are important in order to present a robust framework where New Zealand communities such as Christchurch can gradually recover after a major earthquake disaster, while planning for pre-disaster mitigation against future earthquakes. AM - Accepted Manuscript
The full scale, in-situ investigations of instrumented buildings present an excellent opportunity to observe their dynamic response in as-built environment, which includes all the real physical properties of a structure under study and its surroundings. The recorded responses can be used for better understanding of behavior of structures by extracting their dynamic characteristics. It is significantly valuable to examine the behavior of buildings under different excitation scenarios. The trends in dynamic characteristics, such as modal frequencies and damping ratios, thus developed can provide quantitative data for the variations in the behavior of buildings. Moreover, such studies provide invaluable information for the development and calibration of realistic models for the prediction of seismic response of structures in model updating and structural health monitoring studies. This thesis comprises two parts. The first part presents an evaluation of seismic responses of two instrumented three storey RC buildings under a selection of 50 earthquakes and behavioral changes after Ms=7.1 Darfield (2010) and Ms=6.3 Christchurch (2011) earthquakes for an instrumented eight story RC building. The dynamic characteristics of the instrumented buildings were identified using state-of-the-art N4SID system identification technique. Seismic response trends were developed for the three storey instrumented buildings in light of the identified frequencies and the peak response accelerations (PRA). Frequencies were observed to decrease with excitation level while no trends are discernible for the damping ratios. Soil-structure interaction (SSI) effects were also determined to ascertain their contribution in the seismic response. For the eight storey building, it was found through system identification that strong nonlinearities in the structural response occurred and manifested themselves in all identified natural frequencies of the building that exhibited a marked decrease during the strong motion duration compared to the pre-Darfield earthquakes. Evidence of foundation rocking was also found that led to a slight decrease in the identified modal frequencies. Permanent stiffness loss was also observed after the strong motion events. The second part constitutes developing and calibrating finite element model (FEM) of the instrumented three storey RC building with a shear core. A three dimensional FEM of the building is developed in stages to analyze the effect of structural, non-structural components (NSCs) and SSI on the building dynamics. Further to accurately replicate the response of the building following the response trends developed in the first part of the thesis, sensitivity based model updating technique was applied. The FEMs were calibrated by tuning the updating parameters which are stiffnesses of concrete, NSCs and soil. The updating parameters were found to generally follow decreasing trends with the excitation level. Finally, the updated FEM was used in time history analyses to assess the building seismic performance at the serviceability limit state shaking. Overall, this research will contribute towards better understanding and prediction of the behavior of structures subjected to ground motion.
This research employs a deterministic seismic risk assessment methodology to assess the potential damage and loss at meshblock level in the Christchurch CBD and Mount Pleasant primarily due to building damage caused by earthquake ground shaking. Expected losses in terms of dollar value and casualties are calculated for two earthquake scenarios. Findings are based on: (1) data describing the earthquake ground shaking and microzonation effects; (2) an inventory of buildings by value, floor area, replacement value, occupancy and age; (3) damage ratios defining the performance of buildings as a function of earthquake intensity; (4) daytime and night-time population distribution data and (5) casualty functions defining casualty risk as a function of building damage. A GIS serves as a platform for collecting, storing and analyzing the original and the derived data. It also allows for easy display of input and output data, providing a critical functionality for communication of outcomes. The results of this study suggest that economic losses due to building damage in the Christchurch CBD and Mount Pleasant will possibly be in the order of $5.6 and $35.3 million in a magnitude 8.0 Alpine fault earthquake and a magnitude 7.0 Ashley fault earthquake respectively. Damage to non-residential buildings constitutes the vast majority of the economic loss. Casualty numbers are expected to be between 0 and 10.
The 2010-2011 Canterbury earthquakes were recorded over a dense strong motion network in the near-source region, yielding significant observational evidence of seismic complexities, and a basis for interpretation of multi-disciplinary datasets and induced damage to the natural and built environment. This paper provides an overview of observed strong motions from these events and retrospective comparisons with both empirical and physics-based ground motion models. Both empirical and physics-based methods provide good predictions of observations at short vibration periods in an average sense. However, observed ground motion amplitudes at specific locations, such as Heathcote Valley, are seen to systematically depart from ‘average’ empirical predictions as a result of near surface stratigraphic and topographic features which are well modelled via sitespecific response analyses. Significant insight into the long period bias in empirical predictions is obtained from the use of hybrid broadband ground motion simulation. The comparison of both empirical and physics-based simulations against a set of 10 events in the sequence clearly illustrates the potential for simulations to improve ground motion and site response prediction, both at present, and further in the future.
During the Christchurch earthquake of February 2011, several midrise buildings of Reinforced Concrete Masonry (RCM) construction achieved performance levels in the range of life safety to near collapse levels. These buildings were subjected to seismic demands higher than the building code requirements of the time and higher than the current New Zealand Loadings Standard (NZS-1170.5:2004). Structural damage to these buildings has been documented and is currently being studied to establish lessons to be learned from their performance and how to incorporate these lessons into future RCM design and construction practices. This paper presents a case study of a six story RCM building deemed to have reached the near collapse performance level. The RCM walls on the 2nd floor failed due to toe crushing reducing the building’s lateral resistance in the east-west direction. A nonlinear dynamic analysis on a 3D model was conducted to simulate the development of the governing failure mechanism. Preliminary analysis results show that the damaged walls were initially under large compression forces from gravity loads which caused increase in their lateral strength and reduced their ductility. After toe crushing failure developed, axial instability of the model was prevented by a redistribution of gravity loads.
The Canterbury earthquake series of 2010/2011 has turned the city of Christchurch into a full scale natural laboratory testing the structural and non-structural response of buildings under moderate to very severe earthquake shaking. The lessons learned from this, which have come at great cost socially and economically, are extremely valuable in increasing our understanding of whole building performance in severe earthquakes. Given current initiatives underway on both sides of the Tasman towards developing joint Australasian steel and composite steel/concrete design and construction standards that would span a very wide range of geological conditions and seismic zones, these lessons are relevant to both countries. This paper focusses on the performance of steel framed buildings in Christchurch city, with greatest emphasis on multi-storey buildings, but also covering single storey steel framed buildings and light steel framed housing. It addresses such issues as the magnitude and structural impact of the earthquake series, importance of good detailing, lack of observed column base hinging, the excellent performance of composite floors and it will briefly cover research underway to quantify some of these effects for use in design.
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.
Blended learning plays an important role in many tertiary institutions but little has been written about the implementation of blended learning in times of adversity, natural disaster or crisis. This paper describes how, in the wake of the 22 February Canterbury earthquake, five teacher educators responded to crisis-driven changing demands and changing directions. Our narratives describe how blended learning provided students in initial teacher education programmes with some certainty and continuity during a time of civil emergency. The professional learning generated from our experiences provides valuable insights for designing and preparing for blended learning in times of crisis, as well as developing resilient blended learning programmes for the future.
This thesis presents an assessment of historic seismic performance of the New Zealand stopbank network from the 1968 Inangahua earthquake through to the 2016 Kaikōura earthquake. An overview of the types of stopbanks and the main aspects of the design and construction of earthen stopbanks was presented. Stopbanks are structures that are widely used on the banks of rivers and other water bodies to protect against the impact of flood events. Earthen stopbanks are found to be the most used for such protection measures. Different stopbank damage or failure modes that may occur due to flooding or earthquake excitation were assessed with a focus on past earthquakes internationally, and examples of these damage and failure modes were presented. Stopbank damage and assessment reports were collated from available reconnaissance literature to develop the first geospatial database of stopbank damage observed in past earthquakes in New Zealand. Damage was observed in four earthquakes over the past 50 years, with a number of earthquakes resulting in no stopbank damage. The damage database therefore focussed on the Edgecumbe, Darfield, Christchurch and Kaikōura earthquakes. Cracking of the crest and liquefaction-induced settlement were the most common forms of damage observed. To understand the seismic demand on the stopbank network in past earthquakes, geospatial analyses were undertaken to approximate the peak ground acceleration (PGA) across the stopbank network for ten large earthquakes that have occurred in New Zealand over the past 50 years. The relationship between the demand, represented by the peak ground acceleration (PGA) and damage is discussed and key trends identified. Comparison of the seismic demand and the distribution of damage suggested that the seismic performance of the New Zealand stopbank network has been generally good across all events considered. Although a significant length of the stopbank networks were exposed to high levels of shaking in past events, the overall damage length was a small percentage of this. The key aspect controlling performance was the performance of the underlying foundation soils and the effect of this on the stopbank structure and stability.
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
This study examines the performance of nonlinear total-stress wave-propagation site response analysis for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015).
The 4th of September 2010 Mw 7.1 Darfield (Canterbury) earthquake had generated significant ground shaking within the Christchurch Central Business District (CBD). Despite the apparently significant shaking, the observed structural damage for pre-1970s reinforced concrete (RC) buildings was indeed limited and lower than what was expected for such typology of buildings. This paper explores analytically and qualitatively the different aspects of the "apparent‟ good seismic performance of the pre-1970s RC buildings in the Christchurch CBD, following the earthquake reconnaissance survey by the authors. Damage and building parameters survey result, based on a previously established inventory of building stock of these non-ductile RC buildings, is briefly reported. From an inventory of 75 buildings, one building was selected as a numerical case-study to correlate the observed damage with the non-linear analyses. The result shows that the pre-1970s RC frame buildings performed as expected given the intensity of the ground motion shaking during the Canterbury earthquake. Given the brittle nature of this type of structure, it was demonstrated that more significant structural damage and higher probability of collapse could occur when the buildings were subjected to alternative input signals with different frequency content and duration characteristics and still compatible to the seismicity hazard for Christchurch CBD.
