The objective of this project is to collect perishable seismic response data from the baseisolated Christchurch Women's Hospital. The strong and continuing sequence of aftershocks presents a unique opportunity to capture high-fidelity data from a modern base-isolated facility. These measurements will provide quantitative information required to assess the mechanisms at play in this and in many other seismically-isolated structures.
Between September 2010 and February 2012 (a period of 18 months) the Canterbury region of New Zealand has experienced over 10,000 earthquakes (Nicholls, 2012). This report is the first in a series that will describe the impact of the Canterbury earthquake on businesses. This initial report gives a high level overview of the earthquake events and the impacts on the Canterbury economy and businesses. This report is intended to provide background and context for more in-depth analyses to come in future reports.
The 22 February 2011, Mw6.2 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing an estimated 181 fatalities and severely damaging thousands of residential and commercial buildings. This paper presents a summary of some of the observations made by the NSF-sponsored GEER Team regarding the geotechnical/geologic aspects of this earthquake. The Team focused on documenting the occurrence and severity of liquefaction and lateral spreading, performance of building and bridge foundations, buried pipelines and levees, and significant rockfalls and landslides. Liquefaction was pervasive and caused extensive damage to residential properties, water and wastewater networks, high-rise buildings, and bridges. Entire neighborhoods subsided, resulting in flooding that caused further damage. Additionally, liquefaction and lateral spreading resulted in damage to bridges and to stretches of levees along the Waimakariri and Kaiapoi Rivers. Rockfalls and landslides in the Port Hills damaged several homes and caused several fatalities.
High-Force-to-Volume lead dampers (HF2V) have been recently developed through an experimental research program at University of Canterbury – New Zealand. Testing of the device and applications on beam column joints have demonstrated stable hysteretic behaviour with almost no damage. This paper reports testing of HF2V devices with straight, bulged and constricted shaft configurations subjected to velocities of 0.15 - 5.0mm/s. The effect of the shaft configuration on the hysteresis loop shape, design relationships and the effect of the velocity on the resistive force of the device are described. Results show that hysteresis loop shape of the device is almost square regardless of the shaft configuration, and that devices are characterized by noticeable velocity dependence in the range of 0.15-1.0mm/s.
Micro - electro - mechanical system (MEMS) based accelerometers are now frequently used in many different parts of our day - to - day lives. It is also increasingly being used for structural testing applications. Researchers have had res ervation of using these devices as they are relatively untested, but now with the wider adoption, it provides a much cheaper and more versatile tool for structural engineering researchers. A number of damaged buildings in the Christchurch Central Business District (CBD) were instrumented with a number of low - cost MEMS accelerometers after the major Christchurch earthquakes. The accelerometers captured extremely high quality building response data as the buildings experienced thousands of aftershocks. This d ata set was amongst one of only a handful of data set s available around the world which provides building response data subjected to real ground motion. Furthermore, due to technological advances, a much larger than usual number of accelerometers has been deployed making the data set one of the most comprehensive available. This data set is utilised to extract modal parameters of the buildings. This paper summarises the operating requirements and preference for using such accelerometers for experimental mod al analysis. The challenges for adapting MEMS based devices for successful modal parameters identification are also discussed.
Following exposure to trauma, stress reactions are initially adaptive. However, some individuals’ psychological response can become maladaptive with long-lasting impairment to functioning. Most people with initial symptoms of stress recover, and thus it is important to distinguish individuals who are at risk of continuing difficulties so that resources are allocated appropriately. Investigations of predictors of PTSD development have largely focused on relational and combat-related trauma, with very limited research looking at natural disasters. This study assessed the nature and severity of psychological difficulties experienced in 101 people seeking treatment following exposure to a significant earthquake that killed 185 people. Peritraumatic dissociation, posttraumatic stress symptoms, symptoms of anxiety, symptoms of depression, and social isolation were assessed. Descriptive analyses revealed the sample to be a highly impaired group, with particularly high levels of posttraumatic stress symptoms. Path analysis was used to determine whether the experience of some psychological difficulties predicted experience of others. As hypothesised, peritraumatic dissociation was found to predict posttraumatic stress symptoms and symptoms of anxiety. Posttraumatic stress symptoms then predicted symptoms of anxiety and symptoms of depression. Depression and anxiety were highly correlated. Contrary to expectations, social isolation was not significantly related to any other psychological variables. These findings justify the provision of psychological support following a natural disaster and suggest the benefit of assessing peritraumatic dissociation and posttraumatic stress symptoms soon after the event to identify people in need of monitoring and intervention.
The context of this study is the increasing need for public transport as issues over high private vehicle usage are becoming increasingly obvious. Public transport services need to compete with private transport to improve patronage, and issues with reliability need to be addressed. Bus bunching affects reliability through disruptions to the scheduled headways. The purpose of this study was to collect and analyse data to compare how travel time and dwell time vary, to explore the variation of key variables, and to better understand the sources of these variations. The Orbiter bus service in Christchurch was used as a case study, as it is particularly vulnerable to bus bunching. The dwell time was found to be more variable than travel time. It appeared the Canterbury earthquake had significantly reduced the average speeds for the Orbiter service. In 1964, Newell and Potts described a basic bus bunching theory, which was used as the basis for an Excel bus bunching model. This model allows input variables to vary stochastically. Random values were generated from four specified distributions derived from manually collected data, allowing variance across all bus platforms and buses. However the complexity resulted in stability and difficulty in achieving convergence, so the model was run in single Monte Carlo simulations. The outputs were realistic and showed a higher degree of bunching behaviour than previous models. The model demonstrated bunching phenomena that had not been observed in previous models, including spontaneously un-pairing, overtaking of buses delayed at platforms, and odd-numbered bunches of three buses. Furthermore, the study identified areas of further research for data collection and model development.
Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.
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.
The Mw 7.1 Darfield earthquake generated a ~30 km long surface rupture on the Greendale Fault and significant surface deformation related to related blind faults on a previously unrecognized fault system beneath the Canterbury Plains. This earthquake provided the opportunity for research into the patterns and mechanisms of co-seismic and post-seismic crustal deformation. In this thesis I use multiple across-fault EDM surveys, logic trees, surface investigations and deformation feature mapping, seismic reflection surveying, and survey mark (cadastral) re-occupation using GPS to quantify surface displacements at a variety of temporal and spatial scales. My field mapping investigations identified shaking and crustal displacement-induced surface deformation features south and southwest of Christchurch and in the vicinity of the projected surface traces of the Hororata Blind and Charing Cross Faults. The data are consistent with the high peak ground accelerations and broad surface warping due to underlying reverse faulting on the Hororata Blind Fault and Charing Cross Fault. I measured varying amounts of post-seismic displacement at four of five locations that crossed the Greendale Fault. None of the data showed evidence for localized dextral creep on the Greendale Fault surface trace, consistent with other studies showing only minimal regional post-seismic deformation. Instead, the post-seismic deformation field suggests an apparent westward translation of northern parts of the across-fault surveys relative to the southern parts of the surveys that I attribute to post-mainshock creep on blind thrusts and/or other unidentified structures. The seismic surveys identified a deformation zone in the gravels that we attribute to the Hororata Blind Fault but the Charing Cross fault was not able to be identified on the survey. Cadastral re-surveys indicate a deformation field consistent with previously published geodetic data. We use this deformation with regional strain rates to estimate earthquake recurrence intervals of ~7000 to > 14,000 yrs on the Hororata Blind and Charing Cross Faults.
Earthquake Engineering is facing an extraordinarily challenging era, the ultimate target being set at increasingly higher levels by the demanding expectations of our modern society. The renewed challenge is to be able to provide low-cost, thus more widely affordable, high-seismic-performance structures capable of sustaining a design level earthquake with limited or negligible damage, minimum disruption of business (downtime) or, in more general terms, controllable socio-economical losses. The Canterbury earthquakes sequence in 2010-2011 has represented a tough reality check, confirming the current mismatch between societal expectations over the reality of seismic performance of modern buildings. In general, albeit with some unfortunate exceptions, modern multi-storey buildings performed as expected from a technical point of view, in particular when considering the intensity of the shaking (higher than new code design) they were subjected to. As per capacity design principles, plastic hinges formed in discrete regions, allowing the buildings to sway and stand and people to evacuate. Nevertheless, in many cases, these buildings were deemed too expensive to be repaired and were consequently demolished. Targeting life-safety is arguably not enough for our modern society, at least when dealing with new building construction. A paradigm shift towards damage-control design philosophy and technologies is urgently required. This paper and the associated presentation will discuss motivations, issues and, more importantly, cost-effective engineering solutions to design buildings capable of sustaining low-level of damage and thus limited business interruption after a design level earthquake. Focus will be given to the extensive research and developments in jointed ductile connections based upon controlled rocking & dissipating mechanisms for either reinforced concrete and, more recently, laminated timber structures. An overview of recent on-site applications of such systems, featuring some of the latest technical solutions developed in the laboratory and including proposals for the rebuild of Christchurch, will be provided as successful examples of practical implementation of performance-based seismic design theory and technology.
This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
Previous earthquakes demonstrated destructive effects of soil-structure interaction on structural response. For example, in the 1970 Gediz earthquake in Turkey, part of a factory was demolished in a town 135 km from the epicentre, while no other buildings in the town were damaged. Subsequent investigations revealed that the fundamental period of vibration of the factory was approximately equal to that of the underlying soil. This alignment provided a resonance effect and led to collapse of the structure. Another dramatic example took place in Adapazari, during the 1999 Kocaeli earthquake where several foundations failed due to either bearing capacity exceedance or foundation uplifting, consequently, damaging the structure. Finally, the Christchurch 2012 earthquakes have shown that significant nonlinear action in the soil and soil-foundation interface can be expected due to high levels of seismic excitation and spectral acceleration. This nonlinearity, in turn, significantly influenced the response of the structure interacting with the soil-foundation underneath. Extensive research over more than 35 years has focused on the subject of seismic soil-structure interaction. However, since the response of soil-structure systems to seismic forces is extremely complex, burdened by uncertainties in system parameters and variability in ground motions, the role of soil-structure interaction on the structural response is still controversial. Conventional design procedures suggest that soil-structure interaction effects on the structural response can be conservatively ignored. However, more recent studies show that soil-structure interaction can be either beneficial or detrimental, depending on the soil-structure-earthquake scenarios considered. In view of the above mentioned issues, this research aims to utilise a comprehensive and systematic probabilistic methodology, as the most rational way, to quantify the effects of soil-structure interaction on the structural response considering both aleatory and epistemic uncertainties. The goal is achieved by examining the response of established rheological single-degree-of-freedom systems located on shallow-foundation and excited by ground motions with different spectral characteristics. In this regard, four main phases are followed. First, the effects of seismic soil-structure interaction on the response of structures with linear behaviour are investigated using a robust stochastic approach. Herein, the soil-foundation interface is modelled by an equivalent linear cone model. This phase is mainly considered to examine the influence of soil-structure interaction on the approach that has been adopted in the building codes for developing design spectrum and defining the seismic forces acting on the structure. Second, the effects of structural nonlinearity on the role of soil-structure interaction in modifying seismic structural response are studied. The same stochastic approach as phase 1 is followed, while three different types of structural force-deflection behaviour are examined. Third, a systematic fashion is carried out to look for any possible correlation between soil, structural, and system parameters and the degree of soil-structure interaction effects on the structural response. An attempt is made to identify the key parameters whose variation significantly affects the structural response. In addition, it is tried to define the critical range of variation of parameters of consequent. Finally, the impact of soil-foundation interface nonlinearity on the soil-structure interaction analysis is examined. In this regard, a newly developed macro-element covering both material and geometrical soil-foundation interface nonlinearity is implemented in a finite-element program Raumoko 3D. This model is then used in an extensive probabilistic simulation to compare the effects of linear and nonlinear soil-structure interaction on the structural response. This research is concluded by reviewing the current design guidelines incorporating soil-structure interaction effects in their design procedures. A discussion is then followed on the inadequacies of current procedures based on the outcomes of this study.
