Livelihood holds the key to a rapid recovery following a large-scale devastating disaster, building its resilience is of paramount importance. While much attention has been given to how to help people who are displaced from their jobs to regain employment, little research on livelihood resilience has been undertaken for those relocated communities following a disaster event. By studying five re-located villages post-2004 Indian Ocean Tsunami in Banda Aceh and Aceh Besar, Indonesia, this research has identified the indicators of livelihood resilience and the critical factors driving it for post-disaster relocated communities. A mixed approach, combining questionnaire surveys, semistructured interviews, and field observations, was used for the collection of data. Housing entitlement, the physical and mental health of residents, access to external livelihood support and the provision of infrastructure and basic services were identified as amongst the most critical indicators that represent the level of livelihood resilience. Early recovery income support, physical and mental health, availability and timeliness of livelihood support, together with cultural sensitivity and governance structure, are amongst the most important factors. Given the nature of resettlement, access to infrastructure, location of relocated sites, the safety of the neighbourhood and the ability to transfer to other jobs/skills also play an important role in establishing sustained employment for relocated communities in Indonesia. Those indicators and factors were synthesised into a framework which was further tested in the recovery of Christchurch, and Kaikoura, New Zealand during their recovery from devastating earthquakes. It is suggested that the framework can be used by government agencies and aid organisations to assess the livelihood resilience of post-disaster relocated communities. This will help better them plan support policies and/or prioritise resilience investment strategies to ensure that the recovery needs of those relocated are best met.
Artificial Neural Networks (ANN) as a tool offers opportunities for modeling the inherent complexity and uncertainty associated with socio-environmental systems. This study draws on New Zealand ski
fields (multiple locations) as socio- environmental systems while considering their perceived resilience to low
probability but potential high consequences catastrophic natural events (specifically earthquakes). We gathered
data at several ski fields using a mixed methodology including: geomorphic assessment, qualitative interviews,
and an adaptation of Ozesmi and Ozesmi’s (2003) multi-step fuzzy cognitive mapping (FCM) approach. The data
gathered from FCM are qualitatively condensed, and aggregated to three different participant social groups. The
social groups include ski fields users, ski industry workers, and ski field managers. Both quantitative and
qualitative indices are used to analyze social cognitive maps to identify critical nodes for ANN simulations. The
simulations experiment with auto-associative neural networks for developing adaptive preparation, response and
recovery strategies. Moreover, simulations attempt to identify key priorities for preparation, response, and
recovery for improving resilience to earthquakes in these complex and dynamic environments. The novel mixed
methodology is presented as a means of linking physical and social sciences in high complexity, high uncertainty
socio-environmental systems. Simulation results indicate that participants perceived that increases in Social
Preparation Action, Social Preparation Resources, Social Response Action and Social Response Resources have
a positive benefit in improving the resilience to earthquakes of ski fields’ stakeholders.
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.
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.
This study explored the effects of the Canterbury earthquakes of 2010 and 2011 on different areas of quality of life (QOL) for children and adolescents with disabilities. Using a survey developed from the Quality of Life Instrument for People with Developmental Disabilities – Short Version (QOL-PDD-SV) (Brown, Raphael & Renwick, 1997) and The World Health Organisation Quality of Life - (WHOQOL)-BREF, parents or caregivers were asked to identify what level of importance and satisfaction their child or adolescent placed on areas of QOL including physical health, psychological health - stress levels and coping ability, attachment to their neighbourhood, friends, family, leisure activities, community access and schooling. They were also asked to determine what level of impact the earthquakes had had on each area of their child or adolescent’s life and overall quality of life in the aftermath of the earthquakes. A total of 31 parents of 22 males and 9 females between the ages of 2.5 years to 19 years of age (mean age: 12.6 years) responded. The results were collated and analysis was run to measure for the effect of age, gender and geographical location. The results found that the earthquakes affected nearly every area of QOL for the children and adolescents. The biggest impact on the children’s psychological health and their ability to cope It was observed that younger children (<13) were more likely to record improved or lessened effects from the earthquakes in psychological health areas. However, the areas of social belonging and friendships were the least affected by the earthquakes. Female children were more likely to indicate higher scores for social belonging after the earthquakes. Many parents observed that their children developed improved coping skills over the earthquake period. The findings in this study offer a better understanding of how earthquakes can affect the quality of life children and adolescents with disabilities.
In the aftermath of the 2010-2011 Canterbury Earthquake Sequence (CES), the location of Christchurch-City on the coast of the Canterbury Region (New Zealand) has proven crucial in determining the types of- and chains of hazards that impact the city. Very rapidly, the land subsidence of up to 1 m (vertical), and the modifications of city’s waterways – bank sliding, longitudinal profile change, sedimentation and erosion, engineered stop-banks… - turned rainfall and high-tides into unprecedented floods, which spread across the eastern side of the city. Within this context, this contribution presents two modeling results of potential floods: (1) results of flood models and (2) the effects of further subsidence-linked flooding – indeed if another similar earthquake was to strike the city, what could be the scenarios of further subsidence and then flooding. The present research uses the pre- and post-CES LiDAR datasets, which have been used as the boundary layer for the modeling. On top of simple bathtub model of inundation, the river flood model was conducted using the 2-D hydrodynamic code NAYS-2D developed at the University of Hokkaido (Japan), using a depth-averaged resolution of the hydrodynamic equations. The results have shown that the area the most at risk of flooding are the recent Holocene sedimentary deposits, and especially the swamplands near the sea and in the proximity of waterways. As the CES drove horizontal and vertical displacement of the land-surface, the surface hydrology of the city has been deeply modified, increasing flood risks. However, it seems that scientists and managers haven’t fully learned from the CES, and no research has been looking at the potential future subsidence in further worsening subsidence-related floods. Consequently, the term “coastal quake”, coined by D. Hart is highly topical, and most especially because most of our modern cities and mega-cities are built on estuarine Holocene sediments.
The initial goal of this research was to explore how SME business models change in response to a crisis. Keeping this in mind, the business model canvas (Osterwalder & Pigneur, 2010) was used as a tool to analyse SME business models in the Canterbury region of New Zealand. The purpose was to evaluate the changes SMEs instituted in their business models after being hit by a series of earthquakes in 2010 and 2011. The idea was to conduct interviews with business owners and analyse them using grounded theory methods. As this method is iterative and requires simultaneous data collection and analysis, a tentative model was proposed after first phase of the data collection and analysis. However, as a result of this process, it became apparent that owner-specific characteristics, action orientation and networks were more prominent in the data than business model elements. Although the SMEs in this study experienced several operational changes in their business models, such as a change of location, modifications to their payment terms or expanded/restricted target markets, the suggested framework highlights how owner-specific attributes ensured the recovery of their businesses. After the initial framework was suggested, subsequent interviews were conducted to test, verify, and modify the tentative model. Three aspects of business recovery emerged: (a) cognitive coping – the business owner’s mind-set and motive; (b) adaptive coping – the ability of business owner to take corrective actions; and (c) social capital – the social network of a business owner, including formal and informal connections and their significance. Three distinct groups were identified; self-sufficient SMEs, socially-based SMEs and surviving SMEs. This thesis proposes a grounded theory of business recovery for SMEs following a disaster. Cognitive coping and social capital enabled the owners to take actions, which eventually led to the desired outcomes for the businesses.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
Question: Does canopy tree regeneration response to different large disturbances vary with soil drainage? Location: Old-growth conifer (Dacrydium and Dacrycarpus), angiosperm (Nothofagus and Weinmannia) rain forest, Mount Harata, South Island, New Zealand. Methods: Trees were aged (1056 cores) to reconstruct stand history in 20 (0.12 - 0.2 ha) plots with different underlying drainage. Spatial analyses of an additional 805 tree ages collected from two (0.3 - 0.7 ha) plots were conducted to detect patchiness for five canopy tree species. Microsite preferences for trees and saplings were determined. Results: There were clear differences in species regeneration patterns on soils with different drainage. Conifer recruitment occurred infrequently in even-aged patches (> 1000 m²) and only on poorly drained soils. Periodic Nothofagus fusca and N. menziesii recruitment occurred more frequently in different sized canopy openings on all soils. Weinmannia recruitment was more continuous on all soils reflecting their greater relative shade-tolerance. Distinct periods of recruitment that occurred in the last 400 years matched known large disturbances in the region. These events affected species differently as soil drainage varied. Following earthquakes, both conifers and N. menziesii regenerated on poorly drained soils, while Nothofagus species and Weinmannia regenerated on well-drained soils. However, Dacrydium failed to regenerate after patchy storm damage in the wetter forest interior; instead faster-growing N. fusca captured elevated microsites caused by uprooting. Conclusions: Underlying drainage influenced species composition, while variation in the impacts of large disturbance regulated relative species abundances on different soils.
The Eastern Humps and Leader faults, situated in the Mount Stewart Range in North Canterbury, are two of the ≥17 faults which ruptured during the 2016 MW7.8 Kaikōura Earthquake. The earthquake produced complex, intersecting ground ruptures of these faults and the co-seismic uplift of the Mount Stewart Range. This thesis aims to determine how these two faults accommodated deformation during the 2016 earthquake and how they interact with each other and with pre-existing geological structures. In addition, it aims to establish the most likely subsurface geometry of the fault complex across the Mount Stewart Range, and to investigate the paleoseismic history of the Leader Fault. The Eastern Humps Fault strikes ~240° and dips 80° to 60° to the northwest and accommodated right- lateral – reverse-slip, with up to 4 m horizontal and 2 m vertical displacement in the 2016 earthquake. The strike of the Leader Fault varies from ~155 to ~300°, and dips ~30 to ~80° to the west/northwest, and mainly accommodated left-lateral – reverse-slip of up to 3.5 m horizontal and 3.5 m vertical slip in the 2016 earthquake. On both the Eastern Humps and Leader faults the slip is variable along strike, with areas of low total displacement and areas where horizontal and vertical displacement are negatively correlated. Fault traces with low total displacement reflect the presence of off-fault (distributed) displacement which is not being captured with field measurements. The negative correlation of horizontal and vertical displacement likely indicates a degree of slip partitioning during the 2016 earthquake on both the Eastern Humps and Leader faults. The Eastern Humps and Leader faults have a complex, interdependent relationship with the local bedrock geology. The Humps Fault appears to be a primary driver of ongoing folding and deformation of the local Mendip Syncline and folding of the Mount Stewart Range, which probably began prior to, or synchronous with, initial rupture of The Humps Fault. The Leader Fault appears to use existing lithological weaknesses in the Cretaceous-Cenozoic bedrock stratigraphy to rupture to the surface. This largely accounts for the strong variability on the strike and dip of the Leader Fault, as the geometry of the surface ruptures tend to reflect the strike and dip of the geological strata which it is rupturing through. The Leader Fault may also accommodate some degree of flexural slip in the Cenozoic cover sequence of the Mendip Syncline, contributing to the ongoing growth of the fold. The similarity between topography and uplift profiles from the 2016 earthquake suggest that growth of the Mount Stewart Range has been primarily driven by multiple (>500) discrete earthquakes that rupture The Humps and Leader faults. The spatial distribution of surface displacements across the Mount Stewart Range is more symmetrical than would be expected if uplift is driven primarily by The Humps and Leader faults alone. Elastic dislocation forward models were used to model potential sub-surface geometries and the resulting patterns of deformation compared to photogrammetry-derived surface displacements. Results show a slight preference for models with a steeply southeast-dipping blind fault, coincident with a zone of seismicity at depth, as a ‘backthrust’ to The Humps and Leader faults. This inferred Mount Stewart Fault accommodated contractional strain during the 2016 earthquake and contributes to the ongoing uplift of the Mount Stewart Range with a component of folding. Right-lateral and reverse shear stress change on the Hope Fault was also modelled using Coulomb 3.3 software to examine whether slip on The Humps and Leader faults could transfer enough stress onto the Hope Fault to trigger through-going rupture. Results indicate that during the 2016 earthquake right-lateral shear and reverse stress only increased on the Hope Fault in small areas to the west of the Leader Fault, and similar ruptures would be unlikely to trigger eastward propagating rupture unless the Hope Fault was close to failure prior to the earthquake. Paleoseismic trenches were excavated on the Leader Fault at four locations from 2018 to 2020, revealing near surface (< 4m depth) contractional deformation of Holocene stratigraphy. Three of the trench locations uncovered clear evidence for rupture of the Leader Fault prior to 2016, with fault displacement of near surface stratigraphy being greater than displacement recorded during the 2016 earthquake. Radiocarbon dating of in-situ organic material from two trenches indicate a date of the penultimate earthquake on the Leader Fault within the past 1000 years. This date is consistent with The Humps and Leader faults having ruptured simultaneously in the past, and with multi-fault ruptures involving The Humps, Leader, Hundalee and Stone Jug faults having occurred prior to the 2016 Kaikōura earthquake. Overall, the results contribute to an improved understanding of the Kaikōura earthquake and highlight the importance of detailed structural and paleoseismic investigations in determining controls on earthquake ‘complexity’.
The aim of this thesis was to examine the spatial and the temporal patterns of anxiety and chest pain resulting from the Canterbury, New Zealand earthquaeks. Three research objectives were identified: examine any spatial or termporal clusters of anxiety and chest pain; examine the associations between anxiety, chest pain and damage to neighbourhood; and determine any statistically significant difference in counts of anxiety and chest pain after each earthquake or aftershock which resulted in severe damage. Measures of the extent of liquefaction the location of CERA red-zones were used as proxy measures for earthquake damage. Cases of those who presented to Christchurch Public Hospital Emergency Department with either anxiety or chest pain between May 2010 and April 2012 were aggregated to census area unit (CAU) level for analysis. This thesis has taken a unique approach to examining the spatial and spatio-temporal variations of anxiety and chest pain after an earthquake and offers unique results. This is the first study of its kind to use a GIS approach when examining Canterbury specific earthquake damage and health variables at a CAU level after the earthquakes. Through the use of spatio-termporal scan modelling, negative and linear regression modelling and temporal linear modelling with dummy variables this research was able to conclude there are significant spatial and temporal variations in anxiety and chest pain resulting from the earthquakes. The spatio-termporal scan modelling identified a hot cluster of both anxiety and chest pain within Christchurch at the same time the earthquakes occurred. The negative binomial model found liquefaction to be a stronger predictor of anxiety than the Canterbury Earthquake Recovery Authority's (CERA) land zones. The linear regression model foun chest pain to be positively associated with all measures of earthquake damage with the exception of being in the red-zone. The temporal modelling identified a significant increase in anxiety cases one month after a major earthquake, and chest pain cases spiked two weeks after an earthquake and gradually decreased over the following five weeks. This research was limited by lack of control period data, limited measures of earthquake damage, ethical restrictions, and the need for population tracking data. The findings of this research will be useful in the planning and allocation of mental wellbeing resources should another similar event like the Canterbury Earthquakes occur in New Zealand.
The Stone Jug Fault (SJF) ruptured during the November 14th, 2016 (at 12:02 am), Mw 7.8 Kaikōura Earthquake which initiated ~40 km west-southwest of the study area, at a depth of approximately 15 km. Preliminary post-earthquake mapping indicated that the SJF connects the Conway-Charwell and Hundalee faults, which form continuous surface rupture, however, detailed study of the SJF had not been undertaken prior to this thesis due to its remote location and mountainous topography. The SJF is 19 km long, has an average strike of ~160° and generally carries approximately equal components of sinistral and reverse displacement. The primary fault trace is sigmoidal in shape with the northern and southern tips rotating in strike from NNW to NW, as the SJF approaches the Hope and Hundalee faults. It comprises several steps and bends and is associated with many (N=48) secondary faults, which are commonly near irregularities in the main fault geometry and in a distributed fault zone at the southern tip. The SJF is generally parallel to Torlesse basement bedding where it may utilise pre-existing zones of weakness. Horizontal, vertical and net displacements range up to 1.4 m, with displacement profiles along the primary trace showing two main maxima separated by a minima towards the middle and ends of the fault. Average net displacement along the primary trace is ~0.4m, with local changes in relative values of horizontal and vertical displacement at least partly controlled by fault strike. Two trenches excavated across the northern segment of the fault revealed displacement of mainly Holocene stratigraphy dated using radiocarbon (N=2) and OSL (N=4) samples. Five surface-rupturing paleoearthquakes displaying vertical displacements of <1 m occurred at: 11,000±1000, 7500±1000, 6500±1000, 3500±100 and 3 (2016 Kaikōura) years BP. These events produce an average slip rate since ~11 ka of 0.2-0.4 mm/yr and recurrence intervals of up to 5500 years with an average recurrence interval of 2750 yrs. Comparison of these results with unpublished trench data suggests that synchronous rupture of the Hundalee, Stone Jug, Conway-Charwell, and Humps faults at ~3500 yrs BP cannot be discounted and it is possible that multi-fault ruptures in north Canterbury are more common than previously thought.
The extent of liquefaction in the eastern suburbs of Christchurch (Aranui, Bexley, Avonside, Avonhead and Dallington) from the February 22 2011 Earthquake resulted in extensive damage to in-ground waste water pipe systems. This caused a huge demand for portable toilets (or port-a-loos) and companies were importing them from outside Canterbury and in some instances from Australia. However, because they were deemed “assets of importance” under legislation, their allocation had to be coordinated by Civil Defence and Emergency Management (CDEM). Consequently, companies supplying them had to ignore requests from residents, businesses and rest homes; and commitments to large events outside of the city such as the Hamilton 400 V8 Supercars and the Pasifika Festival in Auckland were impacted. Frustrations started to show as neighbourhoods questioned the equity of the port-a-loos distribution. The Prime Minister was reported as reassuring citizens in the eastern suburbs in the first week of March that1 “a report about the distribution of port-a-loos and chemical toilets shows allocation has been fair. Key said he has asked Civil Defence about the distribution process and where the toilets been sent. He said there aren’t enough for the scale of the event but that is quickly being rectified and the need for toilets is being reassessed all the time.” Nonetheless, there still remained a deep sense of frustration and exclusion over the equity of the port-a-loos distribution. This study took the simple approach of mapping where those port-a-loos were on 11-12 March for several areas in the eastern suburbs and this suggested that their distribution was not equitable and was not well done. It reviews the predictive tools available for estimating damage to waste water pipes and asks the question could this situation have been better planned so that pot-a-loo locations could have been better prioritised? And finally it reviews the integral roles of communication and monitoring as part of disaster management strategy. The impression from this study is that other New Zealand urban centres could or would also be at risk and that work is need to developed more rational management approaches for disaster planning.
The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
Post-tensioned timber technology was originally developed and researched at the University of Canterbury (UC) in New Zealand in 2005. It can provide a low-damage seismic design solution for multi-storey mass timber buildings. Since mass timber products, such as cross-laminated timber (CLT), have high in-plane stiffness, a post-tensioned timber shear wall will deform mainly in a rocking mechanism. The moment capacity of the wall at the base is commonly determined using the elastic form of the Modified Monolithic Beam Analogy (MMBA). In the calculation of the moment capacity at the wall base, it is critical to accurately predict the location of the neutral axis and the timber compressive stress distribution. Three 2/3 scale 8.6m tall post-tensioned CLT walls were experimentally tested under quasi-static cyclic loading – both uni-directional and bi-directional- in this study. These specimens included a single wall, a coupled wall, and a C-shaped core-wall. The main objective was to develop post-tensioned C-shaped timber core-walls for tall timber buildings with enhanced lateral strength and stiffness. To better understand the timber compressive stress distributions at the wall base, particle tracking technology (PTT) technology was applied for the first time to investigate the behaviour of the compression toe. Previous post-tensioned timber testing primarily used the displacement measurements to determine the timber compressive behavior at the wall base or rocking interfaces. However, by using PTT technology, the timber strain measurements in the compression zone can be much more accurate as PTT is able to track the movement of many particles on the timber surface. This paper presents experimental testing results of post-tensioned CLT walls with a focus on capturing timber compressive behavior using PTT. The PTT measurements were able to better capture small base rotations which occurred at the onset of gap opening and capture unexpected phenomena in core-wall tests. The single wall test result herein presented indicates that while the MMBA could predict the moment rotation behavior with reasonable accuracy, the peak strain response was under predicted in the compression toe. Further detailed study is required to better understand the complex strain fields generated reflective of the inherent cross-thickness inhomogeneity and material variability of CLT.
Reinforced concrete (RC) frame buildings designed according to modern design standards achieved life-safety objectives during the Canterbury earthquakes in 2010-11 and the Kaikōura earthquake in 2016. These buildings formed ductile plastic hinges as intended and partial or total building collapse was prevented. However, despite the fact that the damage level of these buildings was relatively low to moderate, over 60% of multi-storey RC buildings in the Christchurch central business district were demolished due to insufficient insurance coverage and significant uncertainty in the residual capacity and repairability of those buildings. This observation emphasized an imperative need to improve understanding in evaluating the post-earthquake performance of earthquake-damaged buildings and to develop relevant post-earthquake assessment guidelines. This thesis focuses on improving the understanding of the residual capacity and repairability of RC frame buildings. A large-scale five-storey RC moment-resisting frame building was tested to investigate the behaviour of earthquake-damaged and repaired buildings. The original test building was tested with four ground motions, including two repeated design-level ground motions. Subsequently, the test building was repaired using epoxy injection and mortar patching and re-tested with three ground motions. The test building was assessed using key concepts of the ATC-145 post-earthquake assessment guideline to validate its assessment procedures and highlight potential limitations. Numerical models were developed to simulate the peak storey drift demand and identify damage locations. Additionally, fatigue assessment of steel reinforcement was conducted using methodologies as per ATC-145. The residual capacity of earthquake-strained steel reinforcement was experimentally investigated in terms of the residual fatigue capacity and the residual ultimate strain capacity. In addition to studying the fatigue capacity of steel reinforcement, the fatigue damage demand was estimated using 972 ground motion records. The deformation limit of RC beams and columns for damage control was explored to achieve a low likelihood of requiring performance-critical repair. A frame component test database was developed, and the deformation capacity at the initiation of lateral strength loss was examined in terms of the chord rotation, plastic rotation and curvature ductility capacity. Furthermore, the proposed curvature ductility capacity was discussed with the current design curvature ductility limits as per NZS 3101:2006.
During 2010 and 2011, a series of major earthquakes caused widespread damage in the city of Christchurch, New Zealand. The magnitude 6.3 quake in February 2011 caused 185 fatalities. In the ensuing months, the government progressively zoned residential land in Christchurch on the basis of its suitability for future occupation (considering damage from these quakes and future earthquake risk). Over 6,000 homes were placed in the ‘red-zone’, meaning that property owners were forced to sell their land to the Crown. This study analysed patterns of residential mobility amongst thirty-one red-zone households from the suburb of Southshore, Christchurch. Drawing on interviews and surveys, the research traced their experience from the zoning announcement until they had moved to a new residence. The research distinguished between short (before the zoning announcement) and long term (post the red zone ‘deadline’) forms of household relocation. The majority of households in the study were highly resistant to short term movement. Amongst those which did relocate before the zoning decision, the desire to maintain a valued social connection with a person outside of the earthquake environment was often an important factor. Some households also moved out of perceived necessity (e.g. due to lack of power or water). In terms of long-term relocation, concepts of affordability and safety were much more highly valued by the sample when purchasing post-quake property. This resulted in a distinct patterning of post-quake housing location choices. Perceived control over the moving process, relationship with government organisations and insurance companies, and time spent in the red-zone before moving all heavily influenced participants’ disaster experience. Contrary to previous studies, households in this study recorded higher levels of subjective well-being after relocating. The study proposed a typology of movers in the Christchurch post-disaster environment. Four mobility behaviours, or types, are identified: the Committed Stayers (CSs), the Environment Re-Creators (ERCs), the Resigned Acceptors (RAs), and the Opportunistic Movers (OMs). The CSs were defined by their immobility rather than their relocation aspirations, whilst the ERCs attempted to recreate or retain aspects of Southshore through their mobility. The RAs expressed a form of apathy towards the post-quake environment, whereas, on the other hand, the OMs moved relative to pre-earthquake plans, or opportunities that arose from the earthquake itself. Possibilities for further research include examining household adaptability to new residential environments and tracking further mobility patterns in the years following relocation from the red- zone.
The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
The Canterbury Region is susceptible to a variety of natural hazards, including earthquakes, landslides and climate hazards. Increasing population and tourism within the region is driving development pressures and as more and more development occurs, the risk from natural hazards increases. In order to avoid development occurring in unacceptably vulnerable locations, natural hazard assessments are required. This study is a reconnaissance natural hazard assessment of Lakes Lyndon, Coleridge and Tekapo. There is restricted potential for development at Lake Lyndon, because the land surrounding the lake is owned by the Crown and has a number of development restrictions. However, there is the potential for conservation or recreation-linked development to occur. There is more potential for development at Lake Coleridge. Most of the land surrounding the lake is privately owned and has less development restrictions. The majority of land surrounding Lake Tekapo is divided into Crown-owned pastoral leases, which are protected from development, such as subdivision. However, there are substantial areas around the lake, which are privately owned and, therefore, have potential for development. Earthquake, landslide and climate hazards are the main natural hazards threatening Lakes Lyndon, Coleridge and Tekapo. The lakes are situated in a zone of active earth deformation in which large and relatively frequent earthquakes are produced. A large number of active faults lie within 15 km of each lake, which are capable of producing M7 or larger earthquakes. Ground shaking, liquefaction, landslides, tsunami and seiches are among the consequences of earthquakes, all of which have the potential to cause severe damage to lives, lifelines and infrastructure. Landslides are also common in the landscape surrounding the lakes. The majority of slopes surrounding the lakes are at significant risk from earthquake-induced failure under moderate to strong earthquake shaking. This level of shaking is expected to occur in any 50 year period around Lakes Lyndon and Coleridge, and in any 150 year period around Lake Tekapo. Injuries, fatalities and property damage can occur directly from landslide impact or from indirect effects such as flooding from landslide-generated tsunami or from landslide dam outbreaks. Lakes Lyndon, Coleridge and Tekapo are also susceptible to climate hazards, such as high winds, drought, heavy snowfall and heavy rainfall, which can lead to landslides and flooding. Future climate change due to global warming is most likely going to affect patterns of frequency and magnitudes of extreme weather events, leading to an increase in climate hazards. Before development is permitted around the lakes, it is essential that each of these hazards is considered so that unacceptably vulnerable areas can be avoided.
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.
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.
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).
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.
Oblique-convergent plate collision between the Pacific and Australian plates across the South Island has resulted in shallow, upper crustal earthquake activity and ground surface deformation. In particular the Porters Pass - Amberley Fault Zone displays a complex hybrid zone of anastomosing dextral strike-slip and thrust/reverse faulting which includes the thrust/reverse Lees Valley Fault Zone and associated basin deformation. There is a knowledge gap with respect to the paleoseismicity of many of the faults in this region including the Lees Valley Fault Zone. This study aimed to investigate the earthquake history of the fault at a selected location and the structural and geomorphic development of the Lees Valley Fault Zone and eastern rangefront. This was investigated through extensive structural and geomorphic mapping, GPS field surveying, vertical aerial photo interpretation, analysis of Digital Elevation Models, paleoseismic trenching and optically stimulated luminescence dating. This thesis used a published model for tectonic geomorphology development of mountain rangefronts to understand the development of Lees Valley. Rangefront geomorphology is investigated through analysis of features such as rangefront sinuosity and faceted spurs and indicates the recently active and episodic nature of the uplifted rangefront. Analysis of fault discontinuity, fault splays, distribution of displacement, fault deformation zone and limited exposure of bedrock provided insight into the complex structure of the fault zone. These observations revealed preserved, earlier rangefronts, abandoned and uplifted within the eastern ranges, indicating a basinward shift in focus of faulting and an imbricate thrust wedge development propagating into the footwall of the fault zone and along the eastern ranges of Lees Valley. Fault scarp deformation analysis indicated multiple events have produced the deformation present preserved by the active fault trace in the northern valley. Vertical deformation along this scarp varied with a maximum of 11.5 m and an average of 5 m. Field mapping revealed fan surfaces of various ages have been offset and deformed, likely during the Holocene, based on expected relative surface ages. Geomorphic and structural mapping highlighted the effect of cross-cutting and inherited structures on the Lees Valley Fault, resulting in a step-over development in the centre of the eastern range-bounding trace. Paleoseismic trenching provided evidence of at least two earthquakes, which were constrained to post 21.6 ± 2.3 ka by optically stimulated luminescence dating. Single event displacements (1.48 ± 0.08 m), surface rupture earthquake magnitudes (Mw 6.7 ± 0.1, with potential to produce ≥ 7.0), and a minimum recurrence interval (3.6 ± 0.3 ka) indicated the Lees Valley Fault is an active structure capable of producing significant earthquake events. Results from this study indicate that the Lees Valley Fault Zone accommodates an important component of the Porters Pass - Amberley Fault Zone deformation and confirms the fault as a source of potentially damaging, peak ground accelerations in the Canterbury region. Remnants of previous rangefronts indicate a thrust wedge development of the Lees Valley Fault Zone and associated ranges that can potentially be used as a model of development for other thrust-fault bounded basins.
High demolition rates were observed in New Zealand after the 2010-2011 Canterbury Earthquake Sequence despite the success of modern seismic design standards to achieve required performance objectives such as life safety and collapse prevention. Approximately 60% of the multi-storey reinforced concrete (RC) buildings in the Christchurch Central Business District were demolished after these earthquakes, even when only minor structural damage was present. Several factors influenced the decision of demolition instead of repair, one of them being the uncertainty of the seismic capacity of a damaged structure. To provide more insight into this topic, the investigation conducted in this thesis evaluated the residual capacity of moderately damaged RC walls and the effectiveness of repair techniques to restore the seismic performance of heavily damaged RC walls. The research outcome provided insights for developing guidelines for post-earthquake assessment of earthquake-damaged RC structures. The methodology used to conduct the investigation was through an experimental program divided into two phases. During the first phase, two walls were subjected to different types of pre-cyclic loading to represent the damaged condition from a prior earthquake, and a third wall represented a repair scenario with the damaged wall being repaired using epoxy injection and repair mortar after the pre-cyclic loading. Comparisons of these test walls to a control undamaged wall identified significant reductions in the stiffness of the damaged walls and a partial recovery in the wall stiffness achieved following epoxy injection. Visual damage that included distributed horizontal and diagonal cracks and spalling of the cover concrete did not affect the residual strength or displacement capacity of the walls. However, evidence of buckling of the longitudinal reinforcement during the pre-cyclic loading resulted in a slight reduction in strength recovery and a significant reduction in the displacement capacity of the damaged walls. Additional experimental programs from the literature were used to provide recommendations for modelling the response of moderately damaged RC walls and to identify a threshold that represented a potential reduction in the residual strength and displacement capacity of damaged RC walls in future earthquakes. The second phase of the experimental program conducted in this thesis addressed the replacement of concrete and reinforcing steel as repair techniques for heavily damaged RC walls. Two walls were repaired by replacing the damaged concrete and using welded connections to connect new reinforcing bars with existing bars. Different locations of the welded connections were investigated in the repaired walls to study the impact of these discontinuities at the critical section. No significant changes were observed in the stiffness, strength, and displacement capacity of the repaired walls compared to the benchmark undamaged wall. Differences in the local behaviour at the critical section were observed in one of the walls but did not impact the global response. The results of these two repaired walls were combined with other experimental programs found in the literature to assemble a database of repaired RC walls. Qualitative and quantitative analyses identified trends across various parameters, including wall types, damage before repair, and repair techniques implemented. The primary outcome of the database analysis was recommendations for concrete and reinforcing steel replacement to restore the strength and displacement capacity of heavily damaged RC walls.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
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.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
