Spatial variations in river facies exerted a strong influence on the distribution of liquefaction features observed in Christchurch during the 2010-11 Canterbury Earthquake Sequence (CES). Liquefaction and liquefaction-induced ground deformation was primarily concentrated near modern waterways and areas underlain by Holocene fluvial deposits with shallow water tables (< 1 to 2 m). In southern Christchurch, spatial variations of liquefaction and subsidence were documented in the suburbs within inner meander loops of the Heathcote River. Newly acquired geospatial data, geotechnical reports and eye-witness discussions are compiled to provide a detailed account of the surficial effects of CES liquefaction and ground deformation adjacent to the Heathcote River. LiDAR data and aerial photography are used to produce a new series of original figures which reveal the locations of recurrent liquefaction and subsidence. To investigate why variable liquefaction patterns occurred, the distribution of surface ejecta and associated ground damage is compared with near-surface sedimentologic, topographic, and geomorphic variability to seek relationships between the near-surface properties and observed ground damages. The most severe liquefaction was concentrated within a topographic low in the suburb of St Martins, an inner meander loop of the Heathcote River, with liquefaction only minor or absent in the surrounding areas. Subsurface investigations at two sites in St Martins enable documentation of fluvial stratigraphy, the expressions of liquefaction, and identification of pre-CES liquefaction features. Excavation to water table depths (~1.5 m below the surface) across sand boils reveals multiple generations of CES liquefaction dikes and sills that cross-cut Holocene fluvial and anthropogenic stratigraphy. Based on in situ geotechnical tests (CPT) indicating sediment with a factor of safety < 1, the majority of surface ejecta was sourced from well-sorted fine to medium sand at < 5 m depth, with the most damaging liquefaction corresponding with the location of a low-lying sandy paleochannel, a remnant river channel from the Holocene migration of the meander in St Martins. In the adjacent suburb of Beckenham, where migration of the Heathcote River has been laterally confined by topography associated with the volcanic lithologies of Banks Peninsula, severe liquefaction was absent with only minor sand boils occurring closest to the modern river channel. Auger sampling across the suburb revealed thick (>1 m) clay-rich overbank and back swamp sediments that produced a stratigraphy which likely confined the units susceptible to liquefaction and prevented widespread ejection of liquefied material. This analysis suggests river migration promotes the formation and preservation of fluvial deposits prone to liquefaction. Trenching revealed the strongest CES earthquakes with large vertical accelerations favoured sill formation and severe subsidence at highly susceptible locations corresponding with an abandoned channel. Less vulnerable sites containing deeper and thinner sand bodies only liquefied in the strongest and most proximal earthquakes forming minor localised liquefaction features. Liquefaction was less prominent and severe subsidence was absent where lateral confinement of a Heathcote meander has promoted the formation of fluvial stratum resistant to liquefaction. Correlating CES liquefaction with geomorphic interpretations of Christchurch’s Heathcote River highlights methods in which the performance of liquefaction susceptibility models can be improved. These include developing a reliable proxy for estimating soil conditions in meandering fluvial systems by interpreting the geology and geomorphology, derived from LiDAR data and modern river morphology, to improve the methods of accounting for the susceptibility of an area. Combining geomorphic interpretations with geotechnical data can be applied elsewhere to identify regional liquefaction susceptibilities, improve existing liquefaction susceptibility datasets, and predict future earthquake damage.
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.
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.
At the conclusion of the 2010 and 2011 Canterbury earthquakes more than 5100 homes had been deemed unsafe for habitation. The land and buildings of these were labelled “red zoned” and are too badly damaged for remediation. These homes have been demolished or are destined for demolition. To assist the red zone population to relocate, central government have offered to ‘buy out’ home owners at the Governmental Value (GV) that was last reviewed in 2007. While generous in the economic context at the time, the area affected was the lowest value land and housing in Christchurch and so there is a capital shortfall between the 2007 property value and the cost of relocating to more expensive properties. This shortfall is made worse by increasing present day values since the earthquakes. Red zone residents have had to relocate to the far North and Western extremities of Christchurch, and some chose to move even further to neighbouring towns or cities. The eastern areas and commercial centres close to the red zone are affected as well. They have lost critical mass which has negatively impacted businesses in the catchments of the Red Zone. This thesis aims to repopulate the suburbs most affected by the abandonment of the red zone houses. Because of the relative scarcity of sound building sites in the East and to introduce affordability to these houses, an alternative method of development is required than the existing low density suburban model. Smart medium density design will be tested as an affordable and appropriate means of living. Existing knowledge in this field will be reviewed, an analysis of what East Christchurch’s key characteristics are will occur, and an examination of built works and site investigations will also be conducted. The research finds that at housing densities of 40 units per hectare, the spatial, vehicle, aesthetic needs of East Christchurch can be accommodated. Centralising development is also found to offer better lifestyle choices than the isolated suburbs at the edges of Christchurch, to be more efficient using existing infrastructure, and to place less reliance on cars. Stronger communities are formed from the outset and for a full range of demographics. Eastern affordable housing options are realised and Christchurch’s ever expanding suburban tendencies are addressed. East Christchurch presently displays a gaping scar of devastated houses that ‘The New Eastside’ provides a bandage and a cure for. Displaced and dispossessed Christchurch residents can be re-housed within a new heart for East Christchurch.
Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit
variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the
imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements
outlined.
Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to
become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid
NO₃- production via nitrification of newly mineralised N during land preparation activities.
Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term
longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly
dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall.
The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.
Mixed conifer, beech and hardwood forests are relatively common in Aotearoa/New
Zealand, but are not well studied. This thesis investigates the coexistence, regeneration
dynamics and disturbance history of a mixed species forest across an environmental
gradient of drainage and soil development in north Westland.
The aim was to investigate whether conifers, beech and non-beech hardwood species were
able to coexist on surfaces that differed in their underlying edaphic conditions, and if so to understand the mechanisms that influenced their regeneration on both poorly drained and
well drained soils. The site selected was an area of high tree species diversity on a lowland
0.8 km² post-glacial terrace at the base of Mount Harata in the Grey River Valley.
My approach was to use forest stand history reconstruction at two spatial scales: an
intensive within-plot study of stand dynamics (chapter 1) and a whole-landform approach
(chapter 2) that examined whether the dynamics identified at the smaller within-plot scale
reflected larger patterns across the terrace.
In chapter 1, three large permanent plots (0.3-0.7 ha) were placed at different points along
the drainage gradient, one plot situated in each of the mainly well-drained, poorly drained
and very poorly drained areas along the terrace. Information was gathered on species age
and size structures, spatial distributions of tree ages, species interactions, microsite
establishment preferences, patterns of stand mortality, and disturbance history in each plot.
There were differences in stand structure, composition and relative abundance of species
found between the well drained plot and the two poorer drained plots. On the well drained
site conifers were scarce, the beeches Nothofagus fusca and N. menziesii dominated the
canopy, and in the subcanopy the hardwood species Weinmannia racemosa and Quintinia
acutifolia were abundant. As drainage became progressively poorer, the conifers
Dacrydium cupressinum and Dacrycarpus dacrydioides became more abundant and
occupied the emergent tier over a beech canopy. The hardwoods W. racemosa and Q.
acutifolia became gradually less abundant in the subcanopy, whereas the hardwood
Elaeocarpus hookerianus became more so.
In the well drained plot, gap partitioning for light between beeches and hardwoods enabled
coexistence in response to a range of different sized openings resulting from disturbances
of different extent. In the two more poorly drained plots, species also coexisted by
partitioning microsite establishment sites according to drainage.
There were several distinct periods where synchronous establishment of different species
occurred in different plots, suggesting there were large disturbances: c. 100yrs, 190-200
yrs, 275-300 yrs and 375-425 yrs ago. Generally after the same disturbance, different
species regenerated in different plots reflecting the underlying drainage gradient. However,
at the same site after different disturbances, different sets of species regenerated,
suggesting the type and extent of disturbances and the conditions left behind influenced
species regeneration at some times but not others. The regeneration of some species (e.g.,
N. fusca in the well-drained plot, and Dacrydium in the poorer drained plots) was periodic
and appeared to be closely linked to these events. In the intervals between these
disturbances, less extensive disturbances resulted in the more frequent N. menziesii and
especially hardwood regeneration. The type of tree death caused by different disturbances
favoured different species, with dead standing tree death favouring the more shade-tolerant
N. menziesii and hardwoods, whereas uprooting created a mosaic of microsite conditions
and larger gap sizes that enabled Dacrycarpus, N. fusca and E. hookerianus to maintain
themselves in the poorly drained areas.
In chapter 2, 10 circular plots (c. 0.12 ha) were placed in well drained areas and 10
circular plots (c. 0.2 ha) in poorly drained plots to collect information on species
population structures and microsite preferences. The aims were to reconstruct species'
regeneration responses to a range of disturbances of different type and extent across the
whole terrace, and to examine whether there were important differences in the effects of
these disturbances.
At this landform scale, the composition and relative abundances of species across the
drainage gradient reflected those found in chapter 1. There were few scattered conifers in well drained areas, despite many potential regeneration opportunities created from a range
of different stand destroying and smaller scale disturbances.
Three of the four periods identified in chapter 1 reflected distinct terrace-wide periods of
regeneration 75-100 yrs, 200-275 yrs and 350-450 yrs ago, providing strong evidence of
periodic large, infrequent disturbances that occurred at intervals of 100-200 yrs. These
large, infrequent disturbances have had a substantial influence in determining forest
history, and have had long term effects on forest structure and successional processes.
Different large, infrequent disturbances had different effects across the terrace, with the
variability in conditions that resulted enabling different species to regenerate at different
times. For example, the regeneration of distinct even-aged Dacrydium cohorts in poorly
drained areas was linked to historical Alpine Fault earthquakes, but not to more recent
storms. The variation in the intensity of different large, infrequent disturbances at different
points along the environmental drainage gradient, was a key factor influencing the scale of
impacts. In effect, the underlying edaphic conditions influenced species composition along
the drainage gradient and disturbance history regulated the relative abundances of species.
The results presented here further emphasise the importance of large scale disturbances as a
mechanism that allows coexistence of different tree species in mixed forest, in particular
for the conifers Dacrydium, Dacrycarpus and the beech N. fusca, by creating much of the
environmental variation to which these species responded. This study adds to our
understanding of the effects of historical earthquakes in the relatively complex forests of
north Westland, and further illustrates their importance in the Westland forest landscape as
the major influential disturbance on forest pattern and history.
These results also further develop the 'two-component' model used to describe
conifer/angiosperm dynamics, by identifying qualitative differences in the impacts of
different large, infrequent disturbances across an environmental gradient that allowed for
coexistence of different species. In poorer drained areas, these forests may even be thought
of as 'three-component' systems with conifers, beeches and hardwoods exhibiting key
differences in their regeneration patterns after disturbances of different type and extent, and
in their microsite preferences.
