The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.
Reconnaissance reports have highlighted the poor performance of non-ductile reinforced concrete buildings during the 2010-11 Canterbury earthquakes. These buildings are widely expected to result in significant losses under future earthquakes due to their seismic vulnerability and prevalence in densely populated urban areas. Wellington, for example, contains more than 70 pre-1970s multi-storey reinforced concrete buildings, ranging in height from 5 to 18 storeys. This study seeks to characterise the seismic performance and evaluate the likely failure modes of a typical pre-1970s reinforced concrete building in Wellington, by conducting advanced numerical simulations to evaluate its 3D nonlinear dynamic response. A representative 9-storey office building constructed in 1951 is chosen for this study and modelled in the finite element analysis programme DIANA, using a previously developed and validated approach to predict the failure modes of doubly reinforced walls with confined boundary regions. The structure consists of long walls and robust framing elements resulting in a stiff lateral load resisting system. Barbell-shaped walls are flanked by stiff columns with sufficient transverse reinforcement to serve as boundary regions. Curved shell elements are used to model the walls and their boundary columns, for which the steel reinforcement is explicitly modelled. Line elements are used to model the frame elements. The steel reinforcement in each member is explicitly modelled. The floor slabs are modelled using elastic shell elements. The model is analysed under short and long duration ground motions selected to match site specific targets in Wellington at the DBE and MCE intensity levels. The observed response of the building including drift profiles at each intesity level, strain localization effects around wall openings, and the influence of bidirectional loading are discussed.
Following the 2010-2011 Canterbury (New Zealand) earthquake sequence, lightly reinforced wall structures in the Christchurch central business district were observed to form undesirable crack patterns in the plastic hinge region, while yield penetration either side of cracks and into development zones was less than predicted using empirical expressions. To some extent this structural behaviour was unexpected and has therefore demonstrated that there may be less confidence in the seismic performance of conventionally designed reinforced concrete (RC) structures than previously anticipated. This paper provides an observation-based comparison between the behaviour of RC structural components in laboratory testing and the unexpected structural behaviour of some case study buildings in Christchurch that formed concentrated inelastic deformations. The unexpected behaviour and poor overall seismic performance of ‘real’ buildings (compared to the behaviour of laboratory test specimens) was due to the localization of peak inelastic strains, which in some cases has arguably led to: (i) significantly less ductility capacity; (ii) less hysteretic energy dissipation; and (iii) the fracture of the longitudinal reinforcement. These observations have raised concerns about whether lightly reinforced wall structures can satisfy the performance objective of “Life Safety” at the Ultimate Limit State. The significance of these issues and potential consequences has prompted a review of potential problems with the testing conditions and procedures that are commonly used in seismic experimentations on RC structures. This paper attempts to revisit the principles of RC mechanics, in particular, the influence of loading history, concrete tensile strength, and the quantity of longitudinal reinforcement on the performance of real RC structures. Consideration of these issues in future research on the seismic performance of RC might improve the current confidence levels in newly designed conventional RC structures.
Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
New Zealand’s stock of unreinforced masonry (URM) bearing wall buildings was principally constructed between 1880 and 1935, using fired clay bricks and lime or cement mortar. These buildings are particularly vulnerable to horizontal loadings such as those induced by seismic accelerations, due to a lack of tensile force-resisting elements in their construction. The poor seismic performance of URM buildings was recently demonstrated in the 2011 Christchurch earthquake, where a large number of URM buildings suffered irreparable damage and resulted in a significant number of fatalities and casualties. One of the predominant failure modes that occurs in URM buildings is diagonal shear cracking of masonry piers. This diagonal cracking is caused by earthquake loading orientated parallel to the wall surface and typically generates an “X” shaped crack pattern due to the reversed cyclic nature of earthquake accelerations. Engineered Cementitious Composite (ECC) is a class of fiber reinforced cement composite that exhibits a strain-hardening characteristic when loaded in tension. The tensile characteristics of ECC make it an ideal material for seismic strengthening of clay brick unreinforced masonry walls. Testing was conducted on 25 clay brick URM wallettes to investigate the increase in shear strength for a range of ECC thicknesses applied to the masonry wallettes as externally bonded shotcrete reinforcement. The results indicated that there is a diminishing return between thickness of the applied ECC overlay and the shear strength increase obtained. It was also shown that, the effectiveness of the externally bonded reinforcement remained constant for one and two leaf wallettes, but decreased rapidly for wall thicknesses greater than two leafs. The average pseudo-ductility of the strengthened wallettes was equal to 220% of that of the as-built wallettes, demonstrating that ECC shotcrete is effective at enhancing both the in-plane strength and the pseudo-ductility of URM wallettes. AM - Accepted Manuscript
Following the devastation of the Canterbury earthquake sequence a unique opportunity exists to rebuild and restructure the city of Christchurch, ensuring that its infrastructure is constructed better than before and is innovative. By installing an integrated grid of modern sensor technologies into concrete structures during the rebuild of the Christchurch CBD, the aim is to develop a network of self-monitored ‘digital buildings’. A diverse range of data will be recorded, potentially including parameters such as concrete stresses, strains, thermal deformations, acoustics and the monitoring of corrosion of reinforcement bars. This procedure will allow an on-going complete assessment of the structure’s performance and service life, both before and after seismic activity. The data generated from the embedded and surface mounted sensors will be analysed to allow an innovative and real-time health monitoring solution where structural integrity is continuously known. This indication of building performance will allow the structure to alert owners, engineers and asset managers of developing problems prior to failure thresholds being reached. A range of potential sensor technologies for monitoring the performance of existing and newly constructed concrete buildings is discussed. A description of monitoring work conducted on existing buildings during the July 2013 Cook Strait earthquake sequence is included, along with details of current work that investigates the performance of sensing technologies for detecting crack formation in concrete specimens. The potential market for managing the real-time health of installed infrastructure is huge. Civil structures all over the world require regular visual inspections in order to determine their structural integrity. The information recorded during the Christchurch rebuild will generate crucial data sets that will be beneficial in understanding the behaviour of concrete over the complete life cycle of the structure, from construction through to operation and building repairs until the time of failure. VoR - Version of Record
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.
Following a natural disaster, children are prone to various reactions and maladaptive responses as a result of exposure to a highly stressful and potentially traumatic event. Children’s responses can range from an acute stress response to post-traumatic-stress disorder or may fall somewhere in between. While responses to highly stressful events vary, a common finding is that children will develop sleep problems. This was found following the Christchurch September 2010 and February 2011 earthquakes. The purpose of this study was to investigate the context and phenomenology of the sleep problems of a small number of children experiencing these and the 2016 Kaikoura earthquakes, including possible mechanisms of effect. Participants were four families, including four mothers, one father and four children. The design of this study was unique. Interview data was subjected to a content analysis, extracted themes were organised according to an ecological-transactional framework and then the factors were subject to an analysis, based on the principles of clinical reasoning, in order to identify possible mechanisms of effect. Parents reported 16 different sleep problems across children, as well as other behaviours possibly indicative of post-traumatic stress response. In total, 34 themes and 26 interactions were extracted in relation to factors identified across participants about the children’s sleep and the families’ earthquake experiences. This demonstrated how complex it is to explore the development of sleep problems in the context of disaster. Key factors identified by parents that likely played a key role in the development and perpetuation of sleep problems included earthquake related anxiety, parental mental health and conflict, the child’s emotional and behavioural problems and other negative life events following the earthquakes. The clinical implications of the analysis included being aware that such families, may not have had access to specialized support around their children’s sleep. This was much needed due to the strain such problems place on the family, especially in a post-disaster community such as Christchurch.
During the past two decades, the focus has been on the need to provide communities with structures that undergo minimal damage after an earthquake event while still being cost competitive. This has led to the development of high performance seismic resisting systems, and advances in design methodologies, in order respect this demand efficiently. This paper presents the experimental response of four pre-cast, post-tensioned rocking wall systems tested on the shake-table at the University of Canterbury. The wall systems were designed as a retrofit solution for an existing frame building, but are equally applicable for use in new design. Design of the wall followed a performance-based retrofit strategy in which structural limit states appropriate to both the post-tensioned wall and the existing building were considered. Dissipation for each of the four post-tensioned walls was provided via externally mounted devices, located in parallel to post-tensioned tendons for re-centring. This allowed the dissipation devices to be easily replaced or inspected following a major earthquake. Each wall was installed with viscous fluid dampers, tension-compression yielding steel dampers, a combination of both or no devices at all – thus relying on contact damping alone. The effectiveness of both velocity and displacement dependant dissipation are investigated for protection against far-field and velocity-pulse ground motion characteristics. The experimental results validate the behaviour of ‘Advanced Flag-Shape’ rocking, dissipating solutions which have been recently proposed and numerically tested. Maximum displacements and material strains were well controlled and within acceptable bounds, and residual deformations were minimal due to the re-centring contribution from the post-tensioned tendons. Damage was confined to inelastic yielding (or fluid damping) of the external dampers.
The 2010-2011 Canterbury Earthquakes brought devastation to the city of Christchurch and has irrevocably affected the lives of the city’s residents. Years after the conclusion of these earthquakes, Christchurch and its residents are well on the path to recovery. Crime has proven an ongoing topic of discussion throughout this period, with news reports of increased burglary and arson in areas left largely abandoned by earthquake damage, and a rise in violent crime in suburban areas of Christchurch. Following the body of research that has considered the reaction of crime to natural disasters, this research has sought to comprehensively examine and understand the effects that the Canterbury Earthquakes had on crime. Examining Christchurch-wide offending, crime rates fell over the study period (July 2008 to June 2013), with the exception of domestic violence. Aside from a momentary increase in burglary in the days immediately following the Christchurch Earthquake, crime rates (as of 2013) have remained largely below pre-earthquake levels. Using Dual Kernel Density Estimation Analysis, a distinct spatial change in pre-earthquake crime hotspots was observed. These changes included an enormous decrease in central city offences, a rise in burglary in the eastern suburbs, and an increase in assault in areas outside of the central city. Logistic regression analysis, using a time-compensated dependent variable, identified a number of statistically-significant relationships between per CAU crime rate change and factors measuring socio-demographic characteristics, community cohesion, and the severity of disaster effects. The significance of these findings was discussed using elements of Social Disorganisation Theory, Routine Activity Theory, and Strain Theory. Consistent with past findings, social order was largely maintained following the Canterbury Earthquakes, with suggestion that increased collective efficacy and therapeutic communities had a negative influence on crime in the post-earthquake period. Areas of increased burglary and assault were associated with large population decreases, suggesting a link with the dissolution of communities and the removal of their inherent informal guardianship. Though observed, the increase in domestic violence was not associated with most neighbourhood-level variables. Trends in crime after the Canterbury Earthquakes were largely consistent with past research, and the media’s portrayal.
Geomorphic, structural and chronological data are used to establish the late Quaternary paleoseismicity of the active dextral-oblique Northern Esk Fault in North Canterbury, New Zealand. Detailed field mapping of the preserved c. 35 km of surface traces between the Hurunui River and Ashley Head reveals variations in strike ranging from 005° to 057°. Along with kinematic data collected from fault plane striae and offset geomorphic markers along the length of the fault these variations are used to distinguish six structural subsections of the main trace, four dextral-reverse and two dextral-normal. Displacements of geomorphic markers such as minor streams and ridges are measured using differential GPS and rangefinder equipment to reveal lateral offsets ranging from 3.4 to 23.7 m and vertical offsets ranging from < 1 to 13.5 m. Characteristic single event displacements of c. 5 m and c. 2 m have been calculated for strike-slip and reverse sections respectively. The use of fault scaling relationships reveals an anomalously high displacement to surface rupture length ratio when compared to global data sets. Fault scaling relationships based on width limited ruptures and magnitude probabilities from point measurements of displacement imply earthquake magnitudes of Mw 7.0 to 7.5. Optically Stimulated Luminescence (OSL) ages from displaced Holocene alluvial terraces at the northern extent of the active trace along with OSL and radiocarbon samples of the central sections constrain the timing of the last two surface rupturing events (11.15 ±1.65 and 3.5 ± 2.8 ka) and suggest a recurrence interval of c. 5612 ± 445 years and late Quaternary reverse and dextral slip rates of c. 0.31 mm/yr and 0.82 mm/yr respectively. The results of this study show that the Northern Esk Fault accommodates an important component of the c. 0.7 – 2 mm/yr of unresolved strain across the plate boundary within the North Canterbury region and affirm the Esk Fault as a source of potentially damaging ground shaking in the Canterbury region.
Oblique convergence of the Pacific and Australian Plates is accommodated in the northern South Island by the Marlborough Fault System. The Hope Fault is the southern of four major dextral strike-slip faults of this system. Hanmer Basin is a probable segment boundary between the Hope River and Conway segments of the Hope Fault. The Conway segment is transpressional and shows increasing structural complexity near the segment boundary at Hanmer Basin, with multiple Late Quaternary traces, and fault-parallel folding in response to across-fault shortening. Between Hossack Station and Hanmer Basin a crush zone in excess of one kilometre wide is exposed in incised streams and rivers. The crush zone has an asymmetrical geometry about the active trace of the Hope Fault, being only 100-300 metres wide south of the fault, and more than 500 metres wide north of the fault. The most intense deformation of Torlesse bedrock occurs at the south side of the fault zone, indicating that strain is accommodated against the fault footwall. North of the fault deformation is less intense, but occurs over a wider area. The wide fault zone at Hossack Station may reflect divergence of the Hanmer Fault, a major splay of the Hope Fault. At Hossack Station, the Hope Fault has accommodated at least 260 metres of dextral displacement during the Holocene. Dating of abandoned stream channels, offset by the Hope Fault, indicated a Late Holocene dextral slip-rate of 18±8 mm-¹ for the west end of the Conway segment. Using empirical formulae and inferred fault parameters, the expected magnitude of an earthquake generated by the Conway segment is M6.9 to M7.4; for an exceedence probability of 10%, the magnitude is M7.7 to M7.9. Effects associated with coseismic rupture of the Conway segment include shaking of up to MMIX along the ruptured fault and at Hanmer Basin. Uplift at the east end of Hanmer Basin, in conjunction with subsidence at the southwest margin of the basin, is resulting in the development of onlapping stratigraphy. Seismic reflection profiles support this theory. Possible along-fault migration of the basin is inferred to be a consequence of non-parallelism of the master faults.
Reinforced concrete buildings that satisfied modern seismic design criteria generally behaved as expected during the recent Canterbury and Kaikoura earthquakes in New Zealand, forming plastic hinges in intended locations. While this meant that life-safety performance objectives were met, widespread demolition and heavy economic losses took place in the aftermath of the earthquakes.The Christchurch central business district was particularly hard hit, with over 60% of the multistorey reinforced concrete buildings being demolished. A lack of knowledge on the post-earthquake residual capacity of reinforced concrete buildings was a contributing factor to the mass demolition.Many aspects related to the assessment of earthquake-damaged reinforced concrete buildings require further research. This thesis focusses on improving the state of knowledge on the post earthquakeresidual capacity and reparability of moderately damaged plastic hinges, with an emphasis on plastic hinges typical of modern moment frame structures. The repair method focussed on is epoxy injection of cracks and patching of spalled concrete. A targeted test program on seventeen nominally identical large-scale ductile reinforced concrete beams, three of which were repaired by epoxy injection following initial damaging loadings, was conducted to support these objectives. Test variables included the loading protocol, the loading rate, and the level of restraint to axial elongation.The information that can be gleaned from post-earthquake damage surveys is investigated. It is shown that residual crack widths are dependent on residual deformations, and are not necessarily indicative of the maximum rotation demands or the plastic hinge residual capacity. The implications of various other types of damage typical of beam and column plastic hinges are also discussed.Experimental data are used to demonstrate that the strength and deformation capacity of plastic hinges with modern seismic detailing are often unreduced as a result of moderate earthquake induced damage, albeit with certain exceptions. Special attention is given to the effects of prior yielding of the longitudinal reinforcement, accounting for the low-cycle fatigue and strain ageing phenomena. A material-level testing program on the low-cycle fatigue behaviour of grade 300E reinforcing steel was conducted to supplement the data available in the literature.A reduction in stiffness, relative to the initial secant stiffness to yield, occurs due to moderate plastic hinging damage. This reduction in stiffness is shown to be correlated with the ductility demand,and a proposed model gives a conservative lower-bound estimate of the residual stiffness following an arbitrary earthquake-type loading. Repair by epoxy injection is shown to be effective in restoring the majority of stiffness to plastic hinges in beams. Epoxy injection is also shown to have implications for the residual strength and elongation characteristics of repaired plastic hinges.
Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.
The last seven years have seen southern New Zealand a ected by several large and damaging earthquakes: the moment magnitude (MW) 7.8 Dusky Sound earthquake on 15 July 2009, the MW 7.1 Dar eld (Canterbury) earthquake on 4 September 2010, and most notably the MW 6.2 Christchurch earthquake on 22 February 2011 and the protracted aftershock sequence. In this thesis, we address the postseismic displacement produced by these earthquakes using methods of satellite-based geodetic measurement, known as Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS), and computational modelling. We observe several ground displacement features in the Canterbury and Fiordland regions during three periods: 1) Following the Dusky Sound earthquake; 2) Following the Dar eld earthquake and prior to the Christchurch earthquake; and 3) Following the Christchurch earthquake until February 2015. The ground displacement associated with postseismic motion following the Dusky Sound earthquake has been measured by continuous and campaign GPS data acquired in August 2009, in conjunction with Di erential Interferometric Synthetic Aperture Radar (DInSAR) observations. We use an afterslip model, estimated by temporal inversion of geodetic data, with combined viscoelastic rebound model to account for the observed spatio-temporal patterns of displacement. The two postseismic processes together induce a signi cant displacement corresponding to principal extensional and contractual strain rates of the order of 10⁻⁷ and 10⁻⁸ yr⁻¹ respectively, across most of the southern South Island. We also analyse observed postseismic displacement following the Dusky Sound earthquake using a new inversion approach in order to describe afterslip in an elasticviscoelastic medium. We develop a mathematical framework, namely the "Iterative Decoupling of Afterslip and Viscoelastic rebound (IDAV)" method, with which to invert temporally dense and spatially sparse geodetic observations. We examine the IDAV method using both numerical and analytical simulations of Green's functions. For the post-Dar eld time interval, postseismic signals are measured within approximately one month of the mainshock. The dataset used for the post-Dar eld displacement spans the region surrounding previously unrecognised faults that ruptured during the mainshock. Poroelastic rebound in a multi-layered half-space and dilatancy recovery at shallow depths provide a satisfactory t with the observations. For the post-Christchurch interval, campaign GPS data acquired in February 2012 to February 2015 in four successive epochs and 66 TerraSAR-X (TSX) SAR acquisitions in descending orbits between March 2011 and May 2014 reveal approximately three years of postseismic displacement. We detect movement away from the satellite of ~ 3 mm/yr in Christchurch and a gradient of displacement of ~ 4 mm/yr across a lineament extending from the westernmost end of the Western Christchurch Fault towards the eastern end of the Greendale East Fault. The postseismic signals following the Christchurch earthquake are mainly accounted for by afterslip models on the subsurface lineament and nearby faults.
The Leader Fault was one of at least 17 faults that ruptured the ground surface across the northeastern South Island of New Zealand during the Mw 7.8 2016 Kaikōura Earthquake. The southern ~6 km of the Leader Fault, here referred to as the South Leader Fault (SLF), ruptured the North Canterbury (tectonic) Domain and is the primary focus of this study. The main objective of the thesis is to understand the key factors that contributed to the geometry and kinematics of the 2016 SLF rupture and its intersection with The Humps Fault (HF). This thesis employs a combination of techniques to achieve the primary objective, including detailed mapping of the bedrock geology, geomorphology and 2016 rupture, measurement of 2016 ground surface displacements, kinematic analysis of slip vectors from the earthquake, and logging of a single natural exposure across a 2016 rupture that was treated as a paleoseismic trench. The resulting datasets were collected in the field, from terrestrial LiDAR and InSAR imagery, and from historical (pre-earthquake) aerial photographs for a ~11 km2 study area. Surface ruptures in the study area are a miniature version of the entire rupture from the earthquake; they are geometrically and kinematically complex, with many individual and discontinuous segments of varying orientations and slip senses which are distributed across a zone up to ~3.5 km wide. Despite this variability, three main groups of ruptures have been identified. These are: 1) NE-SW striking, shallow to moderate dipping (25-45°W) faults that are approximately parallel to Cenozoic bedding with mainly reverse dip-slip and, and for the purposes of this thesis, are considered to be part of the SLF. 2) N-S striking, steeply dipping (~85°E) oblique sinistral faults that are up to the west and part of the SLF. 3) E-NE striking, moderate to steeply dipping (45-68°N) dextral reverse faults which are part of the HF. Bedding-parallel faults are interpreted to be flexural slip structures formed during folding of the near-surface Cenozoic strata, while the steeply dipping SLF ruptured a pre-existing bedrock fault which has little topographic expression. Groups 1 and 2 faults were both locally used for gravitational failure during the earthquake. Despite this non-tectonic fault movement, the slip vectors for faults that ruptured during the earthquake are broadly consistent with NCD tectonics and the regional ~100-120° trend of the principal horizontal stress/strain axes. Previous earthquake activity on the SLF is required by its displacement of Cenozoic formations but Late Quaternary slip on the fault prior to 2016 is neither supported by pre-existing fault scarps nor by changes in topography across the fault. By contrast, at least two earthquakes (including 2016) appear to have ruptured the HF from the mid Holocene, consistent with recurrence intervals of no more than ~7 kyr, and with preliminary observations from trenches on the fault farther to the west. The disparity in paleoearthquake records of the two faults suggests that they typically do not rupture together, thus it is concluded that the HF-SLF rupture pattern observed in the Kaikōura Earthquake rarely occurs in a single earthquake.
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.
This thesis addresses the topic of local bond behaviour in RC structures. The mechanism of bond refers to the composite action between deformed steel reinforcing bars and the surrounding concrete. Bond behaviour is an open research topic with a wide scope, particularly because bond it is such a fundamental concept to structural engineers. However, despite many bond-related research findings having wide applications, the primary contribution of this research is an experimental evaluation of the prominent features of local bond behaviour and the associated implications for the seismic performance of RC structures. The findings presented in this thesis attempt to address some structural engineering recommendations made by the Canterbury Earthquakes Royal Commission following the 2010-2011 Canterbury (New Zealand) earthquake sequence. A chapter of this thesis discusses the structural behaviour of flexure-dominated RC wall structures with an insufficient quantity of longitudinal reinforcement, among other in situ conditions, that causes material damage to predominantly occur at a single crack plane. In this particular case, the extent of concrete damage and bond deterioration adjacent to the crack plane will influence the ductility capacity that is effectively provided by the reinforcing steel. As a consequence of these in situ conditions, some lightly reinforced wall buildings in Christchurch lost their structural integrity due to brittle fracture of the longitudinal reinforcement. With these concerning post-earthquake observations in mind, there is the underlying intention that this thesis presents experimental evidence of bond behaviour that allows structural engineers to re-assess their confidence levels for the ability of lightly reinforced concrete structures to achieve the life-safety seismic performance objective the ultimate limit state. Three chapters of this thesis are devoted to the experimental work that was conducted as the main contribution of this research. Critical details of the experimental design, bond testing method and test programme are reported. The bond stress-slip relationship was studied through 75 bond pull-out tests. In order to measure the maximum local bond strength, all bond tests were carried out on deformed reinforcing bars that did not yield as the embedded bond length was relatively short. Bond test results have been presented in two separate chapters in which 48 monotonic bond tests and 27 cyclic bond tests are presented. Permutations of the experiments include the loading rate, cyclic loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and position of the embedded bond region within the specimen (close or far away to the free surface). The parametric study showed that the concrete strength significantly influences the maximum bond strength and that it is reasonable to normalise the bond stress by the square-root of the concrete compressive strength, √(f'c). The generalised monotonic bond behaviour is described within. An important outcome of the research is that the measured bond strength and stiffness was higher than stated by the bond stress-slip relationship in the fib Model Code 2010. To account for these observed differences, an alternative model is proposed for the local monotonic bond stress-slip relationship. Cyclic bond tests showed a significant proportion of the total bond degradation occurs after the loading cycle in the peak bond strength range, which is when bond slip has exceeded 0.5 mm. Subsequent loading to constant slip values showed a linear relationship between the amount of bond strength degradation and the log of the number of cycles that were applied. To a greater extent, the cyclic bond deterioration depends on the bond slip range, regardless of whether the applied load cycling is half- or fully-reversed. The observed bond deterioration and hysteretic energy dissipated during cyclic loading was found to agree reasonably well between these cyclic tests with different loading protocols. The cyclic bond deterioration was also found to be reasonably consistent exponential damage models found in the literature. This research concluded that the deformed reinforcing bars used in NZ construction, embedded in moderate to high strength concrete, are able to develop high local bond stresses that are mobilised by a small amount of local bond slip. Although the relative rib geometry was not varied within this experimental programme, a general conclusion of this thesis is that deformed bars currently available in NZ have a relative rib bearing area that is comparatively higher than the test bars used in previous international research. From the parametric study it was found that the maximum monotonic bond strength is significant enhanced by dynamic loading rates. Experimental evidence of high bond strength and initial bond stiffness generally suggests that only a small amount of local bond slip that can occur when the deformed test bar was subjected to large tension forces. Minimal bond slip and bond damage limits the effective yielding length that is available for the reinforcing steel to distribute inelastic material strains. Consequently, the potential for brittle fracture of the reinforcement may be a more problematic and widespread issue than is apparent to structural engineers. This research has provided information that improve the reliability of engineering predictions (with respect to ductility capacity) of maximum crack widths and the extent of bond deterioration that might occur in RC structures during seismic actions.
One of the less understood geotechnical responses to the cyclic loading from the MW6.2 Christchurch Earthquake, on the 22nd of February 2011, is the fissuring in the loessial soil-mantled, footslope positions of the north-facing valleys of the Port Hills. The fissures are characterized by mostly horizontal offset (≤500mm), with minor vertical displacement (≤300mm), and they extend along both sides of valleys for several hundred metres in an approximately contour-parallel orientation. The fissure traces correspond to extensional features mapped in other studies. Previous studies have suggested that the fissures are the headscarps of incipient landslides, but the surface and subsurface features are not typical of landslide movement. Whilst there are some features that correlate with landslide movement, there are many features that contradict the landslide movement hypothesis. Of critical importance to this investigation was the fact that there are no landslide flanks, there has been no basal shear surface found, there is little deformation in the so-called ‘landslide body’, and there have been no recorded zones of low shear strength in the soil deposit that are indicative of a basal shear surface. This thesis is a detailed geotechnical study on the fissures along part of Ramahana Road in the Hillsborough Valley, Christchurch. Shallow and deep investigation methods found that the predominant soil is loess-colluvium, to depths of ~20m, and this soil has variable geotechnical characteristics depending on the layer sampled. The factor that has the most influence on shear strength was found to be the moisture content. Direct shear-box testing of disturbed, recompacted loess-colluvium found that the soil had a cohesion of 35-65kPa and a friction angle of 38-43° when the soil moisture content was at 8-10%. However when the moisture content was at 19-20% the soil’s cohesion decreased to 3-5kPa and its friction angle decreased to 33-38°, this moisture content is at or slightly above the plastic limit. An electrical resistivity geophysical survey was conducted perpendicular to multiple fissure traces and through the compressional zone at 17 Ramahana Road. The electrical resistivity line found that there was an area of high resistivity at the toe of the slope, and an area of high conductivity downslope of this and at greater depths. This area correlated to the compressional zone recorded by previous studies. Moisture content testing of the soil in these locations showed that the soil in the resistive area was relatively dry (9%) compared to the surrounding soil (13%), whilst the soil in the conductive area was relatively wet (22%)compared to the surrounding soil (19%). Density tests of the soil in the compressional zone recorded that the resistive area had a higher dry density than the surrounding soil (~1790 kg/m3 compared to ~1650 kg/m3). New springs arose downslope of the compressional zone contemporaneously with the fissures, and it is interpreted that these have arisen from increased hydraulic head in the Banks Peninsula bedrock aquifer system, and earthquake induced-bedrock fracturing. A test pit was dug across an infilled fissure trace at 17 Ramahana Road to a depth of 3m. The fissure trace had an aperture of 450-470mm at the ground surface, but it gradually lost aperture with depth until 2.0-2.1m where it became a segmented fissure trace with 1-2mm aperture. A mixed-colluvium layer was intercepted by the fissure trace at 2.4m depth, and there was no observable vertical offset of this layer. The fissure trace was at an angle of 78° at the ground surface, but it also flattened with depth, which gave it a slightly curved appearance. The fissure trace was at an assumed angle of 40-50° near the base of the test pit. Rotational slide, translational slide and lateral spread landslide movement types were compared and contrasted as possibilities for landslide movement types, whilst an alternative hypothesis was offered that the fissures are tensile failures with a quasi-toppling motion involving a cohesive block of loessial soil moving outwards from the slope, with an accommodating compressional strain in the lower less cohesive soil. The mechanisms behind this movement are suggested to be the horizontal earthquake inertia forces from the Christchurch Earthquake, the static shear stress of the slope, and bedrock uplift of the Port Hills in relation to the subsidence of the Christchurch city flatlands. Extremely high PGA is considered to be a prerequisite to the fissure trace development, and these can only be induced in the Hillsborough Valley from a Port Hills Fault rupture, which has a recurrence interval of ~10,000 years. The current understanding of how the loess-colluvium soil would behave under cyclic loading is limited, and the mechanisms behind the suggested movement type are not completely understood. Further research is needed to confirm the proposed mechanism of the fissure traces. Laboratory tests such as the cyclic triaxial and cyclic shear test would be beneficial in future research to quantitatively test how the soil behaves under cyclic loading at various moisture contents and clay contents, and centrifuge experiments would be of great use to qualitatively test the suggested mode of movement in the loessial soil.
Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.
