Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
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.
Pumice materials, which are problematic from an engineering viewpoint, are widespread in the central part of the North Island. Considering the impacts of the 2010-2011 Christchurch earthquakes, a clear understanding of their properties under earthquake loading is necessary. For example, the 1987 Edgecumbe earthquake showed evidence of localised liquefaction of sands of volcanic origin. To elucidate on this, research was undertaken to investigate whether existing empirical field-based methods to evaluate the liquefaction potential of sands, which were originally developed for hard-grained soils, are applicable to crushable pumice-rich deposits. For this purpose, two sites, one in Whakatane and another in Edgecumbe, were selected where the occurrence of liquefaction was reported following the Edgecumbe earthquake. Manifestations of soil liquefaction, such as sand boils and ejected materials, have been reported at both sites. Field tests, including cone penetration tests (CPT), shear-wave velocity profiling, and screw driving sounding (SDS) tests were performed at the sites. Then, considering estimated peak ground accelerations (PGAs) at the sites based on recorded motions and possible range of ground water table locations, liquefaction analysis was conducted at the sites using available empirical approaches. To clarify the results of the analysis, undisturbed soil samples were obtained at both sites to investigate the laboratory-derived cyclic resistance ratios and to compare with the field-estimated values. Research results clearly showed that these pumice-rich soils do not fit existing liquefaction assessment frameworks and alternate methods are necessary to characterise them.
Tens of thousands of landslides were generated over 10, 000 km2 of North Canterbury and Marlborough as a consequence of the 14 November 2016, MW7.8 Kaikōura Earthquake. The most intense landslide damage was concentrated in 3500 km2 around the areas of fault rupture. Given the sparsely populated area affected by landslides, only a few homes were impacted and there were no recorded deaths due to landslides. Landslides caused major disruption with all road and rail links with Kaikōura being severed. The landslides affecting State Highway 1 (the main road link in the South Island of New Zealand) and the South Island main trunk railway extended from Ward in Marlborough all the way to the south of Oaro in North Canterbury. The majority of landslides occurred in two geological and geotechnically distinct materials reflective of the dominant rock types in the affected area. In the Neogene sedimentary rocks (sandstones, limestones and siltstones) of the Hurunui District, North Canterbury and around Cape Campbell in Marlborough, first-time and reactivated rock-slides and rock-block slides were the dominant landslide type. These rocks also tend to have rock material strength values in the range of 5-20 MPa. In the Torlesse 'basement' rocks (greywacke sandstones and argillite) of the Kaikōura Ranges, first-time rock and debris avalanches were the dominant landslide type. These rocks tend to have material strength values in the range of 20-50 MPa. A feature of this earthquake is the large number (more than 200) of valley blocking landslides it generated. This was partly due to the steep and confined slopes in the area and the widely distributed strong ground shaking. The largest landslide dam has an approximate volume of 12(±2) M m3 and the debris from this travelled about 2.7 km2 downslope where it formed a dam blocking the Hapuku River. The long-term stability of cracked slopes and landslide dams from future strong earthquakes and large rainstorms are an ongoing concern to central and local government agencies responsible for rebuilding homes and infrastructure. A particular concern is the potential for debris floods to affect downstream assets and infrastructure should some of the landslide dams breach catastrophically. At least twenty-one faults ruptured to the ground surface or sea floor, with these surface ruptures extending from the Emu Plain in North Canterbury to offshore of Cape Campbell in Marlborough. The mapped landslide distribution reflects the complexity of the earthquake rupture. Landslides are distributed across a broad area of intense ground shaking reflective of the elongate area affected by fault rupture, and are not clustered around the earthquake epicentre. The largest landslides triggered by the earthquake are located either on or adjacent to faults that ruptured to the ground surface. Surface faults may provide a plane of weakness or hydrological discontinuity and adversely oriented surface faults may be indicative of the location of future large landslides. Their location appears to have a strong structural geological control. Initial results from our landslide investigations suggest predictive models relying only on ground-shaking estimates underestimate the number and size of the largest landslides that occurred.
A team of earthquake geologists, seismologists and engineering seismologists from GNS Science, NIWA, University of Canterbury, and Victoria University of Wellington have collectively produced an update of the 2002 national probabilistic seismic hazard (PSH) model for New Zealand. The new model incorporates over 200 new onshore and offshore fault sources, and utilises newly developed New Zealand-based scaling relationships and methods for the parameterisation of the fault and subduction interface sources. The background seismicity model has also been updated to include new seismicity data, a new seismicity regionalisation, and improved methodology for calculation of the seismicity parameters. Background seismicity models allow for the occurrence of earthquakes away from the known fault sources, and are typically modelled as a grid of earthquake sources with rate parameters assigned from the historical seismicity catalogue. The Greendale Fault, which ruptured during the M7.1, 4 September 2010 Darfield earthquake, was unknown prior to the earthquake. However, the earthquake was to some extent accounted for in the PSH model. The maximum magnitude assumed in the background seismicity model for the area of the earthquake is 7.2 (larger than the Darfield event), but the location and geometry of the fault are not represented. Deaggregations of the PSH model for Christchurch at return periods of 500 years and above show that M7-7.5 fault and background source-derived earthquakes at distances less than 40 km are important contributors to the hazard. Therefore, earthquakes similar to the Darfield event feature prominently in the PSH model, even though the Greendale Fault was not an explicit model input.
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.
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.
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
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 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.
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.
When researchers seek to understand community resilience, it often centres on individual agents and actors. They look at the traits individuals have in order to help recover from adverse events, as well as the decisionmaking processes required to plan and adapt. In Aotearoa New Zealand, Māori forms of organising can challenge these. This research was about uncovering Māori forms organising and practices in the context of resilience. The methodology I used was He Awa Whiria/Braided Rivers and storytelling analysis in kanohi ki te kanohi/semi-structured interviews to understand how Māori communities responded to and recovered from the 2010 Darfield (Canterbury), 2011 Ōtautahi/Christchurch, and 2016 Kaikōura earthquakes. Five themes emerged from the project: (i) the importance of marae as a powerful physical location, (ii) the value in building strong reciprocal connections and cultural relationships, (iii) the stronghold that kai/food has in helping to heal communities, (iv) the exchange and trading of resources, and (v) being practical when move forward after a disaster event. As a non-Māori researcher, I have been an outsider to te Ao Māori and to Aotearoa. In using this blended methodology, it became apparent that there are many socio-cultural and historical contentions from the effects of colonisation, assimilation, to grappling with Western norms. Notably, the findings pointed to more similarities than differences, such as taking care of family and communities, being community-driven, and ways of coping with adverse events. This revealed that there are similar ways of doing things regardless of having different customs. This research makes several contributions. It contributes to the field of management studies by addressing gaps in how the concept of resilience is viewed from a practical Māori perspective. The research presents emergency management professionals with similar blended and practical strategies to co-design approaches for collaborative readiness, response, and recovery plans and programmes. The study further demonstrates the localised and tangible benefits that can be gained from utilising a blended methodology and storying method. Ultimately, the purpose of the thesis was to start bridging the gap between agencies and communities, to shift to more Indigenous-led approaches, integrating local Indigenous practices and knowledges that lead to more prepared communities in managing, responding to, and recovering from earthquake hazard events.
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.
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.
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.
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
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 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.
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.
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.
The Resilient Shorelines study at University of Canterbury (UC) is using the Avon Heathcote Estuary Ihutai to investigate ecosystem-based approaches to conservation planning and adaptation in response to environmental change. In particular, the study is using a novel opportunity to understand effects of the Canterbury earthquakes that may be similar to impacts of sea level rise. These result from topographic and bathymetry changes in and around the estuary and associated waterways (Beaven et al., 2012; Cochran et al., 2014) that have driven changes in hydrodynamics (Measures et al., 2011). Therefore the wider context for the work reported here is to develop methodologies for modelling the impacts of sea level rise on estuaries and coastal river mouths using the Avon-Heathcote Estuary/Ihutai as a case study. Initial objectives have included establishing the magnitude of earthquake-induced changes. Subsequent steps will include establishing the relationships between strong physical drivers such as water levels and salinity, and the spatial pattern of estuarine ecosystems. There is particular focus on understanding salinity changes in the upper estuarine ecosystem in the vicinity of the freshwater-saltwater interface. In these areas, species, habitats and ecosystems that are adapted to brackish conditions are expected to migrate in response to the inland penetration of salt water under sea level rise. An example is the location of īnanga spawning habitat that is associated with the inland extent of salt water intrusion on spring tides (Taylor, 2002). It is expected to be strongly affected by sea level rise. To facilitate the development of ecosystem-based scenario models for sea level rise, a salinity model with resolution at ecological meaningful scales was required. An existing fine scale hydrodynamic model was available using Delft3D software (Deltares, 2012) that had been developed for ECan and MBIE following the earthquakes (Measures & Bind, 2013). However, it had not been calibrated for salinity. A collaborative project was designed between UC and NIWA to calibrate the model and develop a scenario modelling approach for sea level rise at a level of resolution sufficient for understanding sea level rise impacts on īnanga (whitebait) spawning habitat. The project was allocated funding from Brian Mason Scientific and Technical Trust and commenced in late 2015. The purpose of this report is to provide a description of the model development process and an illustration of model outputs from an initial set of modelled scenarios for sea level rise.
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).
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
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.
The University of Canterbury’s RECOVER project (Reef Ecology and Coastal Values, Earthquake Recovery) is a research programme funded by the Ministry of Business, Innovation and Employment (MBIE), and supported by the Ministry of Primary Industries (MPI). It has been evaluating recovery from the 7.8 Mw Kaikōura earthquake in the coastal environment between Oaro in the south and Marfells Beach in the north. The project has documented a wide range of biological and physical impacts in the coastal environment over the past four years. These include the widespread mortality of habitat-forming species that support characteristic ecosystems and natural resources on the coast (Alestra et al. 2021; Schiel et al. 2019; Tait et al. 2021). Due to the popularity of the coast for recreational use, interactions between people and the recovering environment are an important influence on recovery processes. These interactions may include threats to the natural environment but also the potential for positive interventions that could help to restore natural ecosystems and resources – including those that have been degraded in the past. Physical effects of uplift at the coastline include the seaward movement of shorelines and creation of new land above the reach of the tide, leading to a widening of beaches (Orchard et al. 2020; Orchard et al. in press). This has also provided a greater opportunity for off-road vehicle access to sections of the coast previously protected by headlands that were impassable at high tide (Marlborough District Council 2019; Orchard 2020). MDC management responses have included the development of a proposed bylaw to reduce the impacts of motor vehicle use in the area (Marlborough District Council 2021). Changes in the position of the sea-level on the landscape also affect the location of characteristic ecosystems such as sand dunes and storm beaches as they recover to a new norm. Notable changes include the establishment of new dunes closer to the sea which could potentially lead to the degradation of old dune systems that may experience reduced sand supply as a result. Wildlife habitat has also been affected by these uplift and re-assembly effects although the specific impacts remain largely unknown. This report contributes to a collaborative project between the Marlborough District Council (MDC) and University of Canterbury (UC) which aims to help protect and promote the recovery of native dune systems on the Marlborough coast. It is centred around the mapping of dune vegetation and identification of dune protection zones for old-growth seed sources of the native sand-binders spinifex (Spinifex sericeus) and pīngao (Ficinia spiralis). Both are key habitat-formers associated with nationally threatened dune ecosystems (Holdaway et al. 2012), and pīngao is an important weaving resource and Ngāi Tahu taonga species. The primary goal is to protect existing seed sources that are vital for natural regeneration following major disturbances such as the earthquake event. Several additional protection zones are also identified for areas where new dunes are successfully regenerating, including areas being actively restored in the Beach Aid project that is assisting new native dunes to become established where there is available space.
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.
