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.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
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.
On Tuesday 22 February 2011, a 6.3 magnitude earthquake struck Christchurch, New Zealand’s second largest city. The ‘earthquake’ was in fact an aftershock to an earlier 7.1 magnitude earthquake that had occurred on Saturday 4 September 2010. There were a number of key differences between the two events that meant they had dramatically different results for Christchurch and its inhabitants. The 22 February 2011 event resulted in one of New Zealand’s worst natural disasters on record, with 185 fatalities occurring and hundreds more being injured. In addition, a large number of buildings either collapsed or were damaged to the point where they needed to be totally demolished. Since the initial earthquake in September 2010, a large amount of building-related research has been initiated in New Zealand to investigate the impact of the series of seismic events – the major focus of these research projects has been on seismic, structural and geotechnical engineering matters. One project, however, conducted jointly by the University of Canterbury, the Fire Protection Association of New Zealand and BRANZ, has focused on the performance of fire protection systems in the earthquakes and the effectiveness of the systems in the event of post-earthquake fires occurring. Fortunately, very few fires actually broke out following the series of earthquake events in Christchurch, but fire after earthquakes still has significant implications for the built environment in New Zealand, and the collaborative research has provided some invaluable insight into the potential threat posed by post-earthquake fires in buildings. As well as summarising the damage caused to fire protection systems, this paper discusses the flow-on effect for designing structures to withstand post-earthquake fires. One of the underlying issues that will be explored is the existing regulatory framework in New Zealand whereby structural earthquake design and structural design for fire are treated as discrete design scenarios.
To this extent, modern buildings generally demonstrated good resistance to collapse during the recent earthquakes in New Zealand. However, damage to non-structural elements (NSE) has been persistent during these events. NSEs include secondary systems or components attached to the floors, roofs, and walls of a building or industrial facility that are not explicitly designed to participate in the main vertical or lateral load-bearing mechanism of the structure. They play a major role in the operational and functional aspects of buildings and contribute a major portion of the building’s overall cost. Therefore, they are expected to accommodate the effects of seismic actions such as drifts and accelerations. Typical examples of NSEs include internal non-loadbearing partitions, suspended ceilings, sprinkler piping systems, architectural claddings, building contents, mechanical/electrical equipment, and furnishings. The main focus of this thesis is the drift sensitive NSEs: precast concrete cladding panels and internal partition walls. Even though most precast concrete cladding panels performed well from a life-safety point of view during recent earthquakes in NZ, some collapsed panels posed a significant threat to life safety. It is, therefore, important that the design and detailing of the panel-to-structure connections ensure that their strength and displacement capacity are adequate to meet the corresponding seismic demands, at least during design level earthquakes. In contrast, the partition wall is likely to get damaged and lose serviceability at a low inter-story drift unless designed to accommodate the relative deformations between them and the structure. Partition walls suffered wide-ranging damage such as screw failures, diagonal cracking, detachments to the gypsum linings, and anchorage failures during the 2011 Canterbury Earthquake Sequence in NZ. Therefore, the thesis is divided into two parts. Part I of the thesis focuses on developing novel low-damage precast concrete cladding panel connections, i.e. “rocking” connection details comprising vertically slotted steel embeds and weld plates. The low-damage seismic performance of novel “rocking” connection details is verified through experimental tests comprising uni-directional, bi-directional, and multi-storey scaled quasi-static cyclic tests. Comparison with the seismic performance of traditional panel connections reported in the literature demonstrated the system’s significantly improved seismic resilience. Furthermore, the finite element models of panel connections and sealants are developed in ABAQUS. The force-drift responses of the “rocking” panel system modelled in SAP2000 is compared with the experimental results to evaluate their accuracy and validity. Part II of the thesis focuses on a) understanding the seismic performance of traditional rigid timber-framed partition wall, b) development and verification of low-damage connections (i.e. “rocking” connection details comprising of dual-slot tracks), and c) seismic evaluation of partition walls with a novel “bracketed and slotted” connections (comprising of innovative fastener and plastic bracket named Flexibracket) under uni-directional and bidirectional quasi-static cyclic loadings. Moreover, parametric investigation of the partition walls was conducted through several experimental tests to understand better the pros and cons of the rocking connection details. The experimental results have confirmed that the implementation of the proposed low damage solutions of precast cladding panels and internal partition walls can significantly reduce their damage in a building.
By closely examining the performance of a 22-storey steel framed building in Christchurch subject to various earthquakes over the past seven years, it is shown that a number of lessons can be learnt regarding the cost-effective consideration of non-structural elements. The first point in this work is that non-structural elements significantly affected the costs associated with repairing steel eccentrically braced frame (EBF) links. The decommissioning or rerouting of non-structural elements in the vicinity of damaged links in the case study building attributed to approximately half the total cost of their repair. Such costs could be significantly reduced if the original positioning of non-structural elements took account of the potential need to repair the EBF links. The second point highlighted is the role that pre-cast cladding apparently played on the distribution and type of damage in the building. Loss estimates obtained following the FEMA P-58 framework vary considerably when cladding is or isnt modelled, both because of changes to drift demands up the height of the building and because certain types of subsequent damage are likely to be cheaper to repair than others. Finally, costly repairs to non-structural partition walls were required not only after the moment magnitude 7.1 earthquake in 2010 but also in multiple aftershocks in the years that followed. Repair costs associated with aftershock events exceeded those from the main event, emphasizing the need to consider aftershocks within modern performance-based earthquake engineering and also the opportunity that exists to make more cost-effective repair strategies following damaging earthquakes.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
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.
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.
The operation of telecommunication networks is critical during business as usual times, and becomes most vital in post-disaster scenarios, when the services are most needed for restoring other critical lifelines, due to inherent interdependencies, and for supporting emergency and relief management tasks. In spite of the recognized critical importance, the assessment of the seismic performance for the telecommunication infrastructure appears to be underrepresented in the literature. The FP6 QuakeCoRE project “Performance of the Telecommunication Network during the Canterbury Earthquake Sequence” will provide a critical contribution to bridge this gap. Thanks to an unprecedented collaboration between national and international researchers and highly experienced asset managers from Chorus, data and evidences on the physical and functional performance of the telecommunication network after the Canterbury Earthquakes 2010-2011 have been collected and collated. The data will be processed and interpreted aiming to reveal fragilities and resilience of the telecommunication networks to seismic events
The November 2016 MW 7.8 Kaikōura Earthquake initiated beneath the North Culverden basin on The Humps fault and propagated north-eastwards, rupturing at least 17 faults along a cumulative length of ~180 km. The geomorphic expression of The Humps Fault across the Emu Plains, along the NW margin of Culverden basin, comprises a series of near-parallel strands separated by up to 3 km across strike. The various strands strike east to east-northeast and have been projected to mainly dip steeply to the south in seismic data (~80°). In this area, the fault predominantly accommodates right-lateral slip, with uplift and subsidence confined to releasing and restraining bends and step-overs at a range of scales. The Kaikōura event ruptured pre-existing fault scarps along the Emu Plains, which had been partly identified prior to the earthquake. Geomorphology and faulting expression of The Humps Fault on The Emu Plains was mapped, along with faulting related structures which did not rupture in the 2016 earthquake. Fault ruptures strands are combined into sections and the kinematic deformation of sections analysed to provide a moment tensor fault plane solution. This fault plane solution is consistent with the regional principal horizontal shortening direction (PHS) of ~115°, similar to seismic focal mechanism solutions of some of the nearby aftershocks of the Kaikōura earthquake, and similar to the adjacent Hope Fault. To constrain the timing of paleoseismic events, a trench was excavated across the fault where it crossed a late Quaternary alluvial fan. Mapping of stratigraphy exposed in the trench walls, and dating of variably deformed strata, constrains the pre-historic earthquake event history at the trench site. The available data provides evidence for at least three paleo-earthquakes within the last 15.1 ka, with a possible fourth (penultimate) event. These events are estimated to have occurred at 7.7-10.3 ka, 10.3-14.8 ka, and one or more events that are older than ~15.1 ka. Some evidence suggests an additional penultimate event between 1850 C.E and 7.7 ka. Time-integrated slip-rates at three locations on the fault are measured using paleo-channels as piercing points. These sites give horizontal slip rates of 0.57 ± 0.1 mm/year, 0.49 ± 0.1 mm/year and one site constrains a minimum of between 0.1 - 0.4 mm/year. Two vertical slip-rates are calculated to be constrained to a maximum of 0.2 ± 0.02 mm/year at one site and between 0.02 and 0.1 mm/year at another site. Prior to this study, The Humps fault had only been partially documented in reconnaissance level mapping in the district, and no previous paleoseismic or slip rate data had been reported. This project has provided a detailed fault zone tectonic geomorphic map and established new slip-rate and paleoseismic data. The results highlight that The Humps fault plays an important role in regional seismicity and in accommodating plate boundary deformation across the North Canterbury region.
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.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
The recent Canterbury earthquake sequence in 2010-2011 highlighted a uniquely severe level of structural damage to modern buildings, while confirming the high vulnerability and life threatening of unreinforced masonry and inadequately detailed reinforced concrete buildings. Although the level of damage of most buildings met the expected life-safety and collapse prevention criteria, the structural damage to those building was beyond economic repair. The difficulty in the post-event assessment of a concrete or steel structure and the uneconomical repairing costs are the big drivers of the adoption of low damage design. Among several low-damage technologies, post-tensioned rocking systems were developed in the 1990s with applications to precast concrete members and later extended to structural steel members. More recently the technology was extended to timber buildings (Pres-Lam system). This doctoral dissertation focuses on the experimental investigation and analytical and numerical prediction of the lateral load response of dissipative post-tensioned rocking timber wall systems. The first experimental stages of this research consisted of component testing on both external replaceable devices and internal bars. The component testing was aimed to further investigate the response of these devices and to provide significant design parameters. Post-tensioned wall subassembly testing was then carried out. Firstly, quasi-static cyclic testing of two-thirds scale post-tensioned single wall specimens with several reinforcement layouts was carried out. Then, an alternative wall configuration to limit displacement incompatibilities in the diaphragm was developed and tested. The system consisted of a Column-Wall-Column configuration, where the boundary columns can provide the support to the diaphragm with minimal uplifting and also provide dissipation through the coupling to the post-tensioned wall panel with dissipation devices. Both single wall and column-wall-column specimens were subjected to drifts up to 2% showing excellent performance, limiting the damage to the dissipating devices. One of the objectives of the experimental program was to assess the influence of construction detailing, and the dissipater connection in particular proved to have a significant influence on the wall’s response. The experimental programs on dissipaters and wall subassemblies provided exhaustive data for the validation and refinement of current analytical and numerical models. The current moment-rotation iterative procedure was refined accounting for detailed response parameters identified in the initial experimental stage. The refined analytical model proved capable of fitting the experimental result with good accuracy. A further stage in this research was the validation and refinement of numerical modelling approaches, which consisted in rotational spring and multi-spring models. Both the modelling approaches were calibrated versus the experimental results on post-tensioned walls subassemblies. In particular, the multi-spring model was further refined and implemented in OpenSEES to account for the full range of behavioural aspects of the systems. The multi-spring model was used in the final part of the dissertation to validate and refine current lateral force design procedures. Firstly, seismic performance factors in accordance to a Force-Based Design procedure were developed in accordance to the FEMA P-695 procedure through extensive numerical analyses. This procedure aims to determine the seismic reduction factor and over-strength factor accounting for the collapse probability of the building. The outcomes of this numerical analysis were also extended to other significant design codes. Alternatively, Displacement-Based Design can be used for the determination of the lateral load demand on a post-tensioned multi-storey timber building. The current DBD procedure was used for the development of a further numerical analysis which aimed to validate the procedure and identify the necessary refinements. It was concluded that the analytical and numerical models developed throughout this dissertation provided comprehensive and accurate tools for the determination of the lateral load response of post-tensioned wall systems, also allowing the provision of design parameters in accordance to the current standards and lateral force design procedures.
Recent earthquakes in New Zealand proved that a shift is necessary in the current design practice of structures to achieve better seismic performance. Following such events, the number of new buildings using innovative technical solutions (e.g. base isolation, controlled rocking systems, damping devices, etc.), has increased, especially in Christchurch. However, the application of these innovative technologies is often restricted to medium-high rise buildings due to the maximum benefit to cost ratio. In this context, to address this issue, a multi-disciplinary geo-structural-environmental engineering project funded by the Ministry of Business Innovation and Employment (MBIE) is being carried out at the University of Canterbury. The project aims at developing a foundation system which will improve the seismic performance of medium-density low-rise buildings. Such foundation is characterized by two main elements: 1) granulated tyre rubber mixed with gravelly soils to be placed beneath the structure, with the goal of damping part of the seismic energy before it reaches the superstructure; and 2) a basement raft made of steel-fibre rubberised concrete to enhance the flexibility of the foundation under differential displacement demand. In the first part of this paper, the overarching objectives, scope and methodology of the project will be briefly described. Then, preliminary findings on the materials characterization, i.e., the gravel-rubber mixtures and steel-fibre rubberised concrete mixes, will be presented and discussed with focus on the mechanical behaviour.
Despite over a century of study, the relationship between lunar cycles and earthquakes remains controversial and difficult to quantitatively investigate. Perhaps as a consequence, major earthquakes around the globe are frequently followed by 'prediction' claims, using lunar cycles, that generate media furore and pressure scientists to provide resolute answers. The 2010-2011 Canterbury earthquakes in New Zealand were no exception; significant media attention was given to lunarderived earthquake predictions by non-scientists, even though the predictions were merely 'opinions' and were not based on any statistically robust temporal or causal relationships. This thesis provides a framework for studying lunisolar earthquake temporal relationships by developing replicable statistical methodology based on peer reviewed literature. Notable in the methodology is a high accuracy ephemeris, called ECLPSE, designed specifically by the author for use on earthquake catalogs, and a model for performing phase angle analysis. The statistical tests were carried out on two 'declustered' seismic catalogs, one containing the aftershocks from the Mw7.1 earthquake in Canterbury, and the other containing Australian seismicity from the past two decades. Australia is an intraplate setting far removed from active plate boundaries and Canterbury is proximal to a plate boundary, thus allowing for comparison based on tectonic regime and corresponding tectonic loading rate. No strong, conclusive, statistical correlations were found at any level of the earthquake catalogs, looking at large events, onshore events, offshore events, and the fault type of some events. This was concluded using Schuster's test of significance with α=5% and analysis of standard deviations. A few weak correlations, with p-5-10% of rejecting the null hypothesis, and anomalous standard deviations were found, but these are difficult to interpret. The results invalidate the statistical robustness of 'earthquake predictions' using lunisolar parameters in this instance. An ambitious researcher could improve on the quality of the results and on the range of parameters analyzed. The conclusions of the thesis raise more questions than answers, but the thesis provides an adaptable methodology that can be used to further investigation the problem.
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.
Novel Gel-push sampling was employed to obtain high quality samples of Christchurch sands from the Central Business District, at sites where liquefaction was observed in 22 February 2011, and 13 June 2011 earthquakes. The results of cyclic triaxial testing on selected undisturbed specimens of typical Christchurch sands are presented and compared to empirical procedures used by practitioners. This comparison suggests cyclic triaxial data may be conservative, and the Magnitude Scaling Factor used in empirical procedures may be unconservative for highly compressible soils during near source moderate to low magnitude events. Comparison to empirical triggering curves suggests the empirical method generally estimates the cyclic strength of Christchurch sands within a reasonable degree of accuracy as a screening evaluation tool for liquefaction hazard, however for sands with moderate to high fines content it may be significantly unconservative, highlighting the need for high quality sampling and testing on important projects where seismic performance is critical.
The seismic performance and parameter identification of the base isolated Christchurch Women’s Hospital (CWH) building are investigated using the recorded seismic accelerations during the two large earthquakes in Christchurch. A four degrees of freedom shear model is applied to characterize the dynamic behaviour of the CWH building during these earthquakes. A modified Gauss-Newton method is employed to identify the equivalent stiffness and Rayleigh damping coefficients of the building. The identification method is first validated using a simulated example structure and finally applied to the CWH building using recorded measurements from the Mw 6.0 and Mw 5.8 Christchurch earthquakes on December 23, 2011. The estimated response and recorded response for both earthquakes are compared with the cross correlation coefficients and the mean absolute percentage errors reported. The results indicate that the dynamic behaviour of the superstructure and base isolator was essentially within elastic range and the proposed shear linear model is sufficient for the prediction of the structural response of the CWH Hospital during these events.
The Global Earthquake Model’s (GEM) Earthquake Consequences Database (GEMECD) aims to develop, for the first time, a standardised framework for collecting and collating geocoded consequence data induced by primary and secondary seismic hazards to different types of buildings, critical facilities, infrastructure and population, and relate this data to estimated ground motion intensity via the USGS ShakeMap Atlas. New Zealand is a partner of the GEMECD consortium and to-date has contributed with 7 events to the database, of which 4 are localised in the South Pacific area (Newcastle 1989; Luzon 1990; South of Java 2006 and Samoa Islands 2009) and 3 are NZ-specific events (Edgecumbe 1987; Darfield 2010 and Christchurch 2011). This contribution to GEMECD represented a unique opportunity for collating, comparing and reviewing existing damage datasets and harmonising them into a common, openly accessible and standardised database, from where the seismic performance of New Zealand buildings can be comparatively assessed. This paper firstly provides an overview of the GEMECD database structure, including taxonomies and guidelines to collect and report on earthquake-induced consequence data. Secondly, the paper presents a summary of the studies implemented for the 7 events, with particular focus on the Darfield (2010) and Christchurch (2011) earthquakes. Finally, examples of specific outcomes and potentials for NZ from using and processing GEMECD are presented, including: 1) the rationale for adopting the GEM taxonomy in NZ and any need for introducing NZ-specific attributes; 2) a complete overview of the building typological distribution in the Christchurch CBD prior to the Canterbury earthquakes and 3) some initial correlations between the level and extent of earthquake-induced physical damage to buildings, building safety/accessibility issues and the induced human casualties.
This thesis describes the strategies for earthquake strengthening vintage clay bricks unreinforced masonry (URM) buildings. URM buildings are well known to be vulnerable to damage from earthquake-induced lateral forces that may result in partial or full building collapse. The 2010/2011 Canterbury earthquakes are the most recent destructive natural disaster that resulted in the deaths of 185 people. The earthquake events had drawn people’s attention when URM failure and collapse caused about 39 of the fatality. Despite the poor performance of URM buildings during the 2010/2011 Canterbury earthquakes, a number of successful case study buildings were identified and their details research in-depth. In order to discover the successful seismic retrofitting techniques, two case studies of retrofitted historical buildings located in Christchurch, New Zealand i.e. Orion’s URM substations and an iconic Heritage Hotel (aka Old Government Building) was conducted by investigating and evaluating the earthquake performance of the seismic retrofitting technique applied on the buildings prior to the 2010/2011 Canterbury earthquakes and their performance after the earthquakes sequence. The second part of the research reported in this thesis was directed with the primary aim of developing a cost-effective seismic retrofitting technique with minimal interference to the vintage clay-bricks URM buildings. Two retrofitting techniques, (i) near-surface mounted steel wire rope (NSM-SWR) with further investigation on URM wallettes to get deeper understanding the URM in-plane behaviour, and (ii) FRP anchor are reported in this research thesis.
This paper presents a methodology by which both site-specific and spatially distributed ground motion intensity can be obtained immediately following an earthquake event. The methodology makes use of both prediction models for ground motion intensity and its correlation over spatial distances. A key benefit of the methodology is that the ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and also standard deviation, with the standard deviation being a function of the distance to nearby strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence. It is illustrated how these conditional maps can be used for post-event evaluation of liquefaction triggering criteria which have been adopted by the Department of Building and Housing (DBH). Secondly, the conditional distribution of response spectral ordinates is obtained at a specific location for the purposes of determining appropriate ground motion records for use in seismic response analyses of important structures at locations where direct recordings are absent.
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.
A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby Christchurch.
Field surveys and experimental studies have shown that light steel or timber framed plasterboard partition walls are particularly vulnerable to earthquake damage prompting the overarching objective of this research, which is to further the development of low damage seismic systems for non-structural partition walls in order to facilitate their adoption by industry to assist with reducing the losses associated with the maintenance and repair cost of buildings across their design life. In particular, this study focused on the behaviour of steel-framed partition walls systems with novel detailing that aim to be “low-damage” designed according to common practice for walls used in commercial and institutional buildings in New Zealand. This objective was investigated by (1) investigating the performance of a flexible track system proposed by researchers and industry by experimental testing of full-scale specimens; (2) investigating the performance of the seismic gap partition wall systems proposed in a number of studies, further developed in this study with input from industry, by experimental testing of full-scale specimens; and (3) investigating the potential implications of using these systems compared with traditionally detailed partition wall systems within multi-storey buildings using the Performance Based Earthquake Engineering loss assessment methodology. Three full-scale testing frames were designed in order to replicate, under controlled laboratory conditions, the effects of seismic shaking on partition walls within multi-storey buildings by the application of quasi-static uni-directional cyclic loading imposing an inter-storey drift. The typical configuration for test specimens was selected to be a unique “y-shape”, including one angled return wall, with typical dimensions of approximately 2400 mm along the main wall and 600 mm along (approximately) the returns walls with a height of 2405 mm from floor to ceiling. The specimens were aligned within test frames at an oblique angle to the direction of loading in order to investigate bi- directional effects. Three wall specimens with flexible track detailing, two identical plane specimens and the third including a doorway, were tested. The detailing involved removing top track anchors within the proximity of wall intersections, thus allowing the tracks to ‘bow’ out at these locations. Although the top track anchors were specified to be removed the proximity of wall intersections, a construction error was made whereby a single top track slab to concrete anchor was left in at the three-way wall junction. Despite this error, the experimental testing was deemed worthwhile since such errors will also occur in practice and because the behaviour of the wall can be examined with this fixing in mind. The specimens also included an acoustic/fire sealant at the top lining to floor boundary. In addition to providing drift capacities, the force-displacement behaviour is also reported, the dissipated energy was computed, and the parameters of the Wayne-Stewart hysteretic model were fitted to the results. The specimen with the door opening behaved significantly different to the plane specimens: damage to the doorway specimen began as cracking of the wallboard propagating from the corners of the doorway following which the L- and Y- shaped junctions behaved independently, whereas damage to the plane specimens began as cracking of the wallboard at the top of the L-junction and wall system deformed as a single unit. The results suggest that bi-directional behaviour is important even if its impact cannot be directly quantified by the experiments conducted. Damage to sealant implies that the bond between plasterboard and sealant is important for its seismic performance. Careful quality control is advised as defects in the bond may significantly impact its ability to withstand seismic movement. Two specimens with seismic gap detailing were tested: a steel stud specimen and a timber stud specimen. Observed drift capacities were significantly greater than traditional plasterboard partition systems. Equations were used to predict the drift at which damage state 1 (DS1) and damage state 2 (DS2) would initiate. The equation used to estimate the drift at the onset of DS1 accurately predicted the onset of plaster cracking but overestimated the drift at which the gap filling material was damaged. The equation used to predict the onset of DS2 provided a lower bound for both specimens and also when used to predict results of previous experimental tests on seismic gap systems. The gap-filling material reduced the drift at the onset of DS1, however, it had a beneficial effect on the re-centring behaviour of the linings. Out-of-plane displacements and return wall configuration did not appear to significantly impact the onset of plaster cracking in the specimens. A loss assessment according to the PBEE methodology was conducted on four steel MRF case study buildings: (1) a 4-storey building designed for the Christchurch region, (2) a 4-storey building designed for the Wellington region, (3) a 12-storey building designed for the Christchurch region, and (4) a 12- storey building designed for the Wellington region. The fragility parameters for a traditional partition system, the flexible track partition system, and the seismic gap steel stud and timber stud partition systems were included within the loss assessment. The order (lowest to highest) of each system in terms of the expected annual losses of each building when incorporating the system was, (1) the seismic gap timber stud system, (2) the seismic gap steel stud system, (3) the traditional/baseline system, and (4) the flexible track system. For the seismic gap timber stud system, which incurred the greatest reduction in expected annual losses for each case study building, the reduction in expected annual losses in comparison to the losses found when using the traditional system ranged from a 5% to a 30% reduction. This reinforces the fact that while there is a benefit to the using low damage partition systems in each building the extent of reduction in expected annual losses is significantly dependent on the particular building design and its location. The flexible track specimens had larger repair costs at small hazard levels compared to the traditional system but smaller repair costs at larger hazard levels. However, the resulting expected annual losses for the flexible track system was higher than the traditional system which reinforces findings from past studies which observed that the greatest contribution to expected annual losses arises from low to moderate intensity shaking seismic events (low hazard levels).
Validating dynamic responses of engineered systems subjected to simulated ground motions is essential in scrutinising the applicability of simulated ground motions for engineering demand analyses. This paper compares the responses of two 3D building models subjected to recorded and simulated ground motions scaled to the NZS1170.5 design response spectrum, in order to evaluate the applicability of simulated ground motions for use in conventional engineering practice in New Zealand. The buildings were designed according to the NZS1170.5 and physically constructed in Christchurch prior to the 2010-2011 Canterbury earthquakes. 40 recorded ground motions from the 22 February 2011 Christchurch earthquake, along with the simulated ground motions for this event from Razafindrakoto et al. (2018) are considered. The seismic responses of the structures are principally quantified via the peak floor acceleration and maximum inter-storey drift ratio. Overall, the results indicate a general agreement in seismic demands obtained using the recorded and simulated ensembles of ground motions and provide further evidence that simulated ground motions using state-of-the-art methods can be used in code-based structural performance assessments inplace of, or in combination with, ensembles of recorded ground motions.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
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.
Ground motion observations from the most significant 10 events in the 2010-2011 Canterbury earthquake sequence at near-source sites are utilized to scrutinize New Zealand (NZ)-specific pseudo-spectral acceleration (SA) empirical ground motion prediction equations (GMPE) (Bradley 2010, Bradley 2013, McVerry et al. 2006). Region-specific modification factors based on relaxing the conventional ergodic assumption in GMPE development were developed for the Bradley (2010) model. Because of the observed biases with magnitude and source-to-site distance for the McVerry et al. (2006) model it is not possible to develop region-specific modification factors in a reliable manner. The theory of non-ergodic empirical ground motion prediction is then outlined, and applied to this 10 event dataset to determine systematic effects in the between- and within-event residuals which lead to modifications in the predicted median and standard deviation of the GMPE. By examining these systematic effects over sub-regions containing a total of 20 strong motion stations within the Canterbury area, modification factors for use in region-specific ground motion prediction are proposed. These modification factors, in particular, are suggested for use with the Bradley et al. (2010) model in Canterbury-specific probabilistic seismic hazard analysis (PSHA) to develop revised design response, particularly for long vibration periods.
The recent instances of seismic activity in Canterbury (2010/11) and Kaikōura (2016) in New Zealand have exposed an unexpected level of damage to non-structural components, such as buried pipelines and building envelope systems. The cost of broken buried infrastructure, such as pipeline systems, to the Christchurch Council was excessive, as was the cost of repairing building envelopes to building owners in both Christchurch and Wellington (due to the Kaikōura earthquake), which indicates there are problems with compliance pathways for both of these systems. Councils rely on product testing and robust engineering design practices to provide compliance certification on the suitability of product systems, while asset and building owners rely on the compliance as proof of an acceptable design. In addition, forensic engineers and lifeline analysts rely on the same product testing and design techniques to analyse earthquake-related failures or predict future outcomes pre-earthquake, respectively. The aim of this research was to record the actual field-observed damage from the Canterbury and Kaikōura earthquakes of seismic damage to buried pipeline and building envelope systems, develop suitable testing protocols to be able to test the systems’ seismic resilience, and produce prediction design tools that deliver results that reflect the collected field observations with better accuracy than the present tools used by forensic engineers and lifeline analysts. The main research chapters of this thesis comprise of four publications that describe the gathering of seismic damage to pipes (Publication 1 of 4) and building envelopes (Publication 2 of 4). Experimental testing and the development of prediction design tools for both systems are described in Publications 3 and 4. The field observation (discussed in Publication 1 of 4) revealed that segmented pipe joints, such as those used in thick-walled PVC pipes, were particularly unsatisfactory with respect to the joint’s seismic resilience capabilities. Once the joint was damaged, silt and other deleterious material were able to penetrate the pipeline, causing blockages and the shutdown of key infrastructure services. At present, the governing Standards for PVC pipes are AS/NZS 1477 (pressure systems) and AS/NZS 1260 (gravity systems), which do not include a protocol for evaluating the PVC pipes for joint seismic resilience. Testing methodologies were designed to test a PVC pipe joint under various different simultaneously applied axial and transverse loads (discussed in Publication 3 of 4). The goal of the laboratory experiment was to establish an easy to apply testing protocol that could fill the void in the mentioned standards and produce boundary data that could be used to develop a design tool that could predict the observed failures given site-specific conditions surrounding the pipe. A tremendous amount of building envelope glazing system damage was recorded in the CBDs of both Christchurch and Wellington, which included gasket dislodgement, cracked glazing, and dislodged glazing. The observational research (Publication 2 of 4) concluded that the glazing systems were a good indication of building envelope damage as the glazing had consistent breaking characteristics, like a ballistic fuse used in forensic blast analysis. The compliance testing protocol recognised in the New Zealand Building Code, Verification Method E2/VM1, relies on the testing method from the Standard AS/NZS 4284 and stipulates the inclusion of typical penetrations, such as glazing systems, to be included in the test specimen. Some of the building envelope systems that failed in the recent New Zealand earthquakes were assessed with glazing systems using either the AS/NZS 4284 or E2/VM1 methods and still failed unexpectedly, which suggests that improvements to the testing protocols are required. An experiment was designed to mimic the observed earthquake damage using bi-directional loading (discussed in Publication 4 of 4) and to identify improvements to the current testing protocol. In a similar way to pipes, the observational and test data was then used to develop a design prediction tool. For both pipes (Publication 3 of 4) and glazing systems (Publication 4 of 4), experimentation suggests that modifying the existing testing Standards would yield more realistic earthquake damage results. The research indicates that including a specific joint testing regime for pipes and positioning the glazing system in a specific location in the specimen would improve the relevant Standards with respect to seismic resilience of these systems. Improving seismic resilience in pipe joints and glazing systems would improve existing Council compliance pathways, which would potentially reduce the liability of damage claims against the government after an earthquake event. The developed design prediction tool, for both pipe and glazing systems, uses local data specific to the system being scrutinised, such as local geology, dimensional characteristics of the system, actual or predicted peak ground accelerations (both vertically and horizontally) and results of product-specific bi-directional testing. The design prediction tools would improve the accuracy of existing techniques used by forensic engineers examining the cause of failure after an earthquake and for lifeline analysts examining predictive earthquake damage scenarios.
