We present ground motion simulations of the Porters Pass (PP) fault in the Canterbury region of New Zealand; a major active source near Christchurch city. The active segment of the PP fault has an inferred length of 82 km and a mostly strike-slip sense of movement. The PP fault slip makes up approximately 10% of the total 37 mm/yr margin-parallel plate motion and also comprises a significant proportion of the total strain budget in regional tectonics. Given that the closest segment of the fault is less than 45 km from Christchurch city, the PP fault is crucial for accurate earthquake hazard assessment for this major population centre. We have employed the hybrid simulation methodology of Graves and Pitarka (2010, 2015), which combines low (f<1 Hz) and high (f>1 Hz) frequencies into a broadband spectrum. We have used validations from three moderate magnitude events (šš¤4.6 Sept 04, 2010; šš¤4.6 Nov 06, 2010; šš¤4.9 Apr 29, 2011) to build confidence for the šš¤ > 7 PP simulations. Thus far, our simulations include multiple rupture scenarios which test the impacts of hypocentre location and the finite-fault stochastic rupture representation of the source itself. In particular, we have identified the need to use location-specific 1D šš /šš models for the high frequency part of the simulations to better match observations.
Natural hazard reviews reveal increases in disaster impacts nowhere more pronounced than in coastal settlements. Despite efforts to enhance hazard resilience, the common trend remains to keep producing disaster prone places. This paper explicitly explores hazard versus multi-hazard concepts to illustrate how different conceptualizations can enhance or reduce settlement resilience. Understandings gained were combined with onthe-ground lessons from earthquake and flooding experiences to develop of a novel āfirst cutā approach for analyzing key multi-hazard interconnections, and to evaluate resilience enhancing opportunities. Traditional disaster resilience efforts often consider different hazard types discretely. However, recent events in Christchurch, a New Zealand city that is part of the 100 Resilient Cities network, highlight the need to analyze the interrelated nature of different hazards, especially for enhancing lifelines system resilience. Our overview of the Christchurch case study demonstrates that seismic, hydrological, shallow-earth, and coastal hazards can be fundamentally interconnected, with catastrophic results where such interconnections go unrecognized. In response, we have begun to develop a simple approach for use by different stakeholders to support resilience planning, pre and post disaster, by: drawing attention to natural and built environment multi-hazard links in general; illustrating a āfirst cutā tool for uncovering earthquake-flooding multi-hazard links in particular; and providing a basis for reviewing resilience strategy effectiveness in multi-hazard prone environments. This framework has particular application to tectonically active areas exposed to climate-change issues.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120Ā° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
Surface-rupturing earthquakes can trigger the sudden avulsion of river channels, causing rapid and persistent coseismic flooding of previously unaffected areas. This phenomenon, known as fault-rupture-induced river avulsion (FIRA), occurs when fault displacement significantly alters river channel topography. The importance of understanding FIRA as a secondary seismic hazard was highlighted by events during the 2010 Darfield and 2016 Kaikoura earthquakes in New Zealand. This thesis develops a national model to identify and quantify FIRA susceptibility across New Zealand by integrating hydrological datasets (NIWA RiverMaps and Flood Statistics) with active fault information (NZ Active Faults Database and RSQSim earthquake simulations). The methodology applies the F-index framework proposed by McEwan et al. (2023), which quantifies FIRA potential based on the ratio of fault throw plus discharge-dependent depth to bank full depth at each fault-river intersection. The model successfully identified 3,796 potential FIRA-susceptible fault-river intersections nationwide, with 451 involving waterways equal to or larger than the Hororata River. Regional analysis revealed higher concentrations of FIRA-susceptible sites in the Bay of Plenty, Canterbury, and Marlborough regions. Validation against historical events showed the model effectively located known FIRA occurrences from the Kaikoura and Darfield earthquakes, though with some limitations in accurately predicting F-index values due to complex fault displacement patterns and challenges in modelling bank full depths of large, braided rivers. This research establishes New Zealand's first nationwide assessment of fault-induced river avulsion susceptibility. The approach creates a structured methodology for identifying high-risk fault-river intersections and determining which sites require thorough localised examination. The methodology developed offers a template for similar assessments in other tectonically active regions and contributes to improving earthquake hazard assessment and disaster preparedness planning.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and KaikoÅ«ra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of ā¢ How to use seismic hazard model inputs for regional fault displacement hazard analysis ā¢ How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement ā¢ How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting ā¢ The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale ā¢ Which sections of New Zealandās roading network are most susceptible to fault displacement hazard and warrant site-specific investigations ā¢ Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road āstrainā values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawkeās Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
The Lake Coleridge Rock Avalanche Deposits (LCRADs) are located on Ryton Station in the middle Rakaia Valley, approximately 80 km west of Christchurch. Torlesse Supergroup greywacke is the basement material and has been significantly influenced by both active tectonics and glaciation. Both glacial and post-glacial processes have produced large volumes of material which blanket the bedrock on slopes and in the valley floors. The LCRADs were part of a regional study of rock avalanches by WHITEHOUSE (1981, 1983) and WHITEHOUSE and GRIFFITHS (1983), and a single rock avalanche event was recognised with a weathering rind age of 120 years B.P. that was later modified to 150 Ā± 40 years B.P. The present study has refined details of both the age and the sequence of events at the site, by identifying three separate rock avalanche deposits (termed the LCRA1, LCRA2 and LCRA3 deposits), which are all sourced from near the summit of Carriage Drive. The LCRA1 deposit is lobate in shape and had an estimated original deposit volume of 12.5 x 10ā¶ mĀ³, although erosion by the Ryton River has reduced the present day debris volume to 5.1 x 10ā¶ mĀ³. An optically stimulated luminescence date taken from sandy loess immediately beneath the LCRA1 deposit provided a maximum age for the rock avalanche event of 9,720 Ā± 750 years B.P., which is believed to be realistic given that this is shortly after the retreat of Acheron 3 ice from this part of the valley. Emplacement of rock avalanche material into an ancestral Ryton riverbed created a natural dam with a ~17 M mĀ³ lake upstream. The river is thought to have created a natural spillway over the dam structure at ~557 m (a.s.l), and to have existed for a number of years before any significant downcutting occurred. Although a triggering mechanism for the LCRA1 deposit was poorly constrained, it is thought that stress rebound after glacial ice removal may have initiated failure. Due to the event occurring c.10,000 years ago, there was a lack of definition for a possible earthquake trigger, though the possibility is obvious. The LCRAā event had an original deposit volume of 0.66 x 10ā¶ mĀ³, and was constrained to the low-lying area adjacent to the Ryton River that had been created by river erosion of the LCRA1 deposit. Further erosion by the Ryton River has reduced the deposit volume to 0.4 x 10ā¶ mĀ³. A radiocarbon date from a piece of mÄnuka found within the LCRA2 deposit provided an age of 668 Ā± 36 years B.P., and this is thought to reliably date the event. The LCRA2 event also dammed the Ryton River, and the preservation of dam-break outwash terraces downstream from the deposit provides clear evidence of rapid dam erosion and flooding after overtopping, and breaching by the Ryton River. Based on the mean annual flow of the Ryton River, the LCRA2 lake would have taken approximately two weeks to fill assuming that there were no preferred breach paths and the material was relatively impermeable. The LCRA2 event is thought to have been coseismic with a fault rupture along the western segment of the PPAFZ, which has been dated at 600 Ā± 100 years B.P. by SMITH (2003). The small LCRA3 event was not able to be dated, but it is believed to have failed shortly after the LCRA2 event and it may in fact be a lag deposit of the second rock avalanche event possibly triggered by an aftershock. The deposit is only visible at one locality within the cliffs that line the Ryton River, and its lack of geomorphic expression is attributed to it occurring closely after the LCRA2 event, while the Ryton River was still dammed from the second rock avalanche event. A wedge-block of some 35,000 mĀ³ of source material for a future rock avalanche was identified at the summit of Carriage Drive. The dilation of the rock mass, combined with unfavourably oriented sub-vertical bedding in the Torlesse Supergroup bedrock, has allowed toppling-style failure on both of the main ridge lines around the source area for the LCRADs. In the event of a future rock avalanche occurring within the Ryton riverbed an emergency response plan has been developed to provide a staged response, especially in relation to the camping ground located at the mouth of the Ryton River. A long-term management plan has also been developed for mitigation measures for the Ryton riverbed and adjacent floodplain areas downstream of a future rock avalanche at the LCRAD site.
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.
