The Leader Fault was one of at least 17 faults that ruptured the ground surface across the northeastern South Island of New Zealand during the Mw 7.8 2016 Kaikōura Earthquake. The southern ~6 km of the Leader Fault, here referred to as the South Leader Fault (SLF), ruptured the North Canterbury (tectonic) Domain and is the primary focus of this study. The main objective of the thesis is to understand the key factors that contributed to the geometry and kinematics of the 2016 SLF rupture and its intersection with The Humps Fault (HF). This thesis employs a combination of techniques to achieve the primary objective, including detailed mapping of the bedrock geology, geomorphology and 2016 rupture, measurement of 2016 ground surface displacements, kinematic analysis of slip vectors from the earthquake, and logging of a single natural exposure across a 2016 rupture that was treated as a paleoseismic trench. The resulting datasets were collected in the field, from terrestrial LiDAR and InSAR imagery, and from historical (pre-earthquake) aerial photographs for a ~11 km2 study area. Surface ruptures in the study area are a miniature version of the entire rupture from the earthquake; they are geometrically and kinematically complex, with many individual and discontinuous segments of varying orientations and slip senses which are distributed across a zone up to ~3.5 km wide. Despite this variability, three main groups of ruptures have been identified. These are: 1) NE-SW striking, shallow to moderate dipping (25-45°W) faults that are approximately parallel to Cenozoic bedding with mainly reverse dip-slip and, and for the purposes of this thesis, are considered to be part of the SLF. 2) N-S striking, steeply dipping (~85°E) oblique sinistral faults that are up to the west and part of the SLF. 3) E-NE striking, moderate to steeply dipping (45-68°N) dextral reverse faults which are part of the HF. Bedding-parallel faults are interpreted to be flexural slip structures formed during folding of the near-surface Cenozoic strata, while the steeply dipping SLF ruptured a pre-existing bedrock fault which has little topographic expression. Groups 1 and 2 faults were both locally used for gravitational failure during the earthquake. Despite this non-tectonic fault movement, the slip vectors for faults that ruptured during the earthquake are broadly consistent with NCD tectonics and the regional ~100-120° trend of the principal horizontal stress/strain axes. Previous earthquake activity on the SLF is required by its displacement of Cenozoic formations but Late Quaternary slip on the fault prior to 2016 is neither supported by pre-existing fault scarps nor by changes in topography across the fault. By contrast, at least two earthquakes (including 2016) appear to have ruptured the HF from the mid Holocene, consistent with recurrence intervals of no more than ~7 kyr, and with preliminary observations from trenches on the fault farther to the west. The disparity in paleoearthquake records of the two faults suggests that they typically do not rupture together, thus it is concluded that the HF-SLF rupture pattern observed in the Kaikōura Earthquake rarely occurs in a single earthquake.
Research on responses to trauma has historically focused on the negative repercussions of a struggle with adversity. However, more recently, researchers have begun to examine posttraumatic growth: the positive psychological change that emerges from the struggle with a potentially traumatic event. Associations have been found between posttraumatic growth and greater peritraumatic distress, greater objective severity of trauma exposure, greater perceived stressfulness of events, social support, female gender, cognitive and behavioural responses to trauma, and personality measures. Posttraumatic growth has been measured typically in individuals with varying levels of posttraumatic stress disorder symptoms and other psychological difficulties, such as depression and anxiety. Although some theory and research posits that higher resilience would prohibit posttraumatic growth, no studies have examined posttraumatic growth in a resilient sample. The Canterbury earthquake sequence of 2010 and 2011 involved potentially traumatic events that saw the community struggle with a variety of challenges. However, in the midst of earthquake destruction, some positive initiatives emerged, driven by locals. The Gap Filler project (using city spaces left empty from fallen buildings for art and interactive community projects) and the Student Volunteer Army (groups of volunteers coordinated to help others in need) are examples of this. In this context, it seemed likely that posttraumatic growth was occurring and might be seen in individuals who were coping well with challenges. Culture is theorised to influence the posttraumatic growth process (Calhoun, Cann, & Tedeschi, 2010), and the nature of the trauma undergone is also likely to influence the process of growth. The current thesis measures posttraumatic growth quantitatively and qualitatively in a New Zealand sample. It measures and describes posttraumatic growth in a resilient population after the earthquake sequence of 2010 and 2011 in Canterbury, New Zealand. Findings are used to test current models of posttraumatic growth for individuals coping well after trauma and to elaborate on mechanisms proposed by models such as the comprehensive model of posttraumatic growth (Calhoun et al., 2010) and the organismic valuing theory of growth through adversity (Joseph & Linley, 2005). Correlates of posttraumatic growth are examined and likely supporting factors of posttraumatic growth are identified for this population. Study 1 used quantitative analysis to explore correlates of posttraumatic growth and found that greater posttraumatic growth related to greater peritraumatic distress, greater perceived stressfulness of earthquake events, greater objective stressfulness of earthquake events, greater difficulty with stressful life events, less satisfaction with social support, and female gender. Findings from Study 1 give important detail about the nature of distress included in the comprehensive model of posttraumatic growth (Calhoun et al., 2010) for this population. Levels of posttraumatic growth were lower than those in North American studies but similar to those in a Chinese study. The current sample, however, showed lower endorsement of Relating to Others than the Chinese study, perhaps because of cultural differences. Study 2 used qualitative analysis to examine the experience of posttraumatic growth in the sample. The theme of ‘a greater sense of community’ was found and adds to the comprehensive model of posttraumatic growth, in that an expression of posttraumatic growth (a greater connection with others) can inform ongoing social processing in the posttraumatic growth process. Having a formal or informal role in earthquake recovery appeared to influence self-concept and reflection; this elaborates on the influence of role on reflection in Calhoun et al.’s model. Findings illustrate possible mechanisms of the organismic valuing process theorised by Joseph and Linley (2005). Implications include the importance of providing opportunities for individuals to take on a role after a crisis, encouraging them to act to respond to difficulties, and encouraging them to meet personal needs for relatedness, competence, and autonomy. Finding positive aspects to a difficult situation, as well as acknowledging adversity, can be supported in future to help individuals process their traumas. As a society, we can help individuals cope with adversity by providing ways they can meet their needs for relatedness, competence, and autonomy. Community groups likely provide opportunities for members to act in ways that meet such needs. This will allow them to effectively act to meet their needs in times of crisis.
Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.
This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
Coastal margins are exposed to rising sea levels that present challenging circumstances for natural resource management. This study investigates a rare example of tectonic displacement caused by earthquakes that generated rapid sea-level change in a tidal lagoon system typical of many worldwide. This thesis begins by evaluating the coastal squeeze effects caused by interactions between relative sea-level (RSL) rise and the built environment of Christchurch, New Zealand, and also examples of release from similar effects in areas of uplift where land reclamations were already present. Quantification of area gains and losses demonstrated the importance of natural lagoon expansion into areas of suitable elevation under conditions of RSL rise and showed that they may be necessary to offset coastal squeeze losses experienced elsewhere. Implications of these spatial effects include the need to provide accommodation space for natural ecosystems under RSL rise, yet other land-uses are likely to be present in the areas required. Consequently, the resilience of these environments depends on facilitating transitions between human land-uses either proactively or in response to disaster events. Principles illustrated by co-seismic sea-level change are generally applicable to climate change adaptation due to the similarity of inundation effects. Furthermore, they highlight the potential role of non-climatic factors in determining the overall trajectory of change. Chapter 2 quantifies impacts on riparian wetland ecosystems over an eight year period post- quake. Coastal wetlands were overwhelmed by RSL rise and recovery trajectories were surprisingly slow. Four risk factors were identified from the observed changes: 1) the encroachment of anthropogenic land-uses, 2) connectivity losses between areas of suitable elevation, 3) the disproportionate effect of larger wetland vulnerabilities, and 4) the need to protect new areas to address the future movement of ecosystems. Chapter 3 evaluates the unique context of shoreline management on a barrier sandspit under sea-level rise. A linked scenario approach was used to evaluate changes on the open coast and estuarine shorelines simultaneously and consider combined effects. The results show dune loss from a third of the study area using a sea-level rise scenario of 1 m over 100 years and with continuation of current land-uses. Increased exposure to natural hazards and accompanying demand for seawalls is a likely consequence unless natural alternatives can be progressed. In contrast, an example of managed retreat following earthquake-induced subsidence of the backshore presents a new opportunity to restart saltmarsh accretion processes seaward of coastal defences with the potential to reverse decades of degradation and build sea-level rise resilience. Considering both shorelines simultaneously highlights the existence of pinch-points from opposing forces that result in small land volumes above the tidal range. Societal adaptation is delicately poised between the paradigms of resisting or accommodating nature and challenged by the long perimeter and confined nature of the sandspit feature. The remaining chapters address the potential for salinity effects caused by tidal prism changes with a focus on the conservation of īnanga (Galaxias maculatus), a culturally important fish that supports New Zealand‘s whitebait fishery. Methodologies were developed to test the hypothesis that RSL changes would drive a shift in the distribution of spawning sites with implications for their management. Chapter 4 describes a new practical methodology for quantifying the total productivity and spatiotemporal variability of spawning sites at catchment scale. Chapter 5 describes the novel use of artificial habitats as a detection tools to help overcome field survey limitations in degraded environments where egg mortality can be high. The results showed that RSL changes resulted in major shifts in spawning locations and these were associated with new patterns of vulnerability due to the continuation of pre-disturbance land-uses. Unexpected findings includes an improved understanding of the spatial relationship between salinity and spawning habitat, and identification of an invasive plant species as important spawning habitat, both with practical management implications. To conclude, the design of legal protection mechanisms was evaluated in relation to the observed habitat shifts and with a focus on two new planning initiatives that identified relatively large protected areas (PAs) in the lower river corridors. Although the larger PAs were better able to accommodate the observed habitat shifts inefficiencies were also apparent due to spatial disparities between PA boundaries and the values requiring protection. To reduce unnecessary trade-offs with other land-uses, PAs of sufficient size to cover the observable spatiotemporal variability and coupled with adaptive capacity to address future change may offer a high effectiveness from a network of smaller PAs. The latter may be informed by both monitoring and modelling of future shifts and these are expected to include upstream habitat migration driven by the identified salinity relationships and eustatic sea-level rise. The thesis concludes with a summary of the knowledge gained from this research that can assist the development of a new paradigm of environmental sustainability incorporating conservation and climate change adaptation. Several promising directions for future research identified within this project are also discussed.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
