Search

found 367 results

Research papers, University of Canterbury Library

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.

Research papers, Victoria University of Wellington

Located on the edge of two tectonic plates, New Zealand has numerous fault lines and seismic risk across the whole country. The way this risk is communicated affects whether people prepare effectively or at all. Research has shown that perceptions of risk are affected by slight changes in wording, and that probabilities commonly reported by experts and media are often interpreted subjectively based on context. In the context of volcanoes, research has found that given a certain probability of a volcano in a specific time window, people perceive risk as higher in later time intervals within that window. The present study examines this pattern with regard to earthquakes and aftershocks in the New Zealand context. Participants in both Wellington (N = 102) and Christchurch (N = 98) were presented an expert statement of earthquake risk within a given time window in Wellington and aftershock risk in Christchurch, and asked to rate their perception of risk in specific intervals across the time window. For a Wellington earthquake, participants perceived risk as incrementally higher toward the end of the 50 year time window whereas for a Christchurch aftershock, risk perception increased slightly for the first three intervals of the 12 month time window. Likelihood of preparing was constant over the time windows, with Wellington citizens rating themselves more likely than Christchurch citizens to prepare for either an earthquake or aftershock, irrespective of current level of preparedness. These findings suggest that people view earthquakes as more likely later toward the end of a given time window and that they view aftershocks very differently to scientific predictions.

Research Papers, Lincoln University

Six stands located on different land forms in mixed old-growth Nothofagus forests in the Matiri Valley (northwest of South Island, New Zealand) were sampled to examine the effects of two recent large earthquakes on tree establishment and tree-ring growth, and how these varied across land forms. 50 trees were cored in each stand to determine age structure and the cores were cross-dated to precisely date unusual periods of radial growth. The 1968 earthquake (M = 7.1, epicentre 35 km from the study area) had no discernible impact on the sampled stands. The impact of the 1929 earthquake (M = 7.7, epicentre 20 km from the study area) varied between stands, depending on whether or not they had been damaged by soil or rock movement. In all stands, the age structures showed a pulse of N. fusca establishment following the 1929 earthquake, with this species dominating establishment in large gaps created by landslides. Smaller gaps, created by branch or tree death, were closed by both N. fusca and N. menziesii. The long period of releases (1929-1945) indicates that direct earthquake damage was not the only cause of tree death, and that many trees died subsequently most likely of pathogen attack or a drought in the early 1930s. The impacts of the 1929 earthquake are compared to a storm in 1905 and a drought in 1974-1978 which also affected forests in the region. Our results confirm that earthquakes are an important factor driving forest dynamics in this tectonically active region, and that the diversity of earthquake impacts is a major source of heterogeneity in forest structure and regeneration.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, Victoria University of Wellington

Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data. In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā. This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.