Search

found 362 results

Research papers, University of Canterbury Library

This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.

Research papers, University of Canterbury Library

The November 2016 MW 7.8 Kaikōura Earthquake initiated beneath the North Culverden basin on The Humps fault and propagated north-eastwards, rupturing at least 17 faults along a cumulative length of ~180 km. The geomorphic expression of The Humps Fault across the Emu Plains, along the NW margin of Culverden basin, comprises a series of near-parallel strands separated by up to 3 km across strike. The various strands strike east to east-northeast and have been projected to mainly dip steeply to the south in seismic data (~80°). In this area, the fault predominantly accommodates right-lateral slip, with uplift and subsidence confined to releasing and restraining bends and step-overs at a range of scales. The Kaikōura event ruptured pre-existing fault scarps along the Emu Plains, which had been partly identified prior to the earthquake. Geomorphology and faulting expression of The Humps Fault on The Emu Plains was mapped, along with faulting related structures which did not rupture in the 2016 earthquake. Fault ruptures strands are combined into sections and the kinematic deformation of sections analysed to provide a moment tensor fault plane solution. This fault plane solution is consistent with the regional principal horizontal shortening direction (PHS) of ~115°, similar to seismic focal mechanism solutions of some of the nearby aftershocks of the Kaikōura earthquake, and similar to the adjacent Hope Fault. To constrain the timing of paleoseismic events, a trench was excavated across the fault where it crossed a late Quaternary alluvial fan. Mapping of stratigraphy exposed in the trench walls, and dating of variably deformed strata, constrains the pre-historic earthquake event history at the trench site. The available data provides evidence for at least three paleo-earthquakes within the last 15.1 ka, with a possible fourth (penultimate) event. These events are estimated to have occurred at 7.7-10.3 ka, 10.3-14.8 ka, and one or more events that are older than ~15.1 ka. Some evidence suggests an additional penultimate event between 1850 C.E and 7.7 ka. Time-integrated slip-rates at three locations on the fault are measured using paleo-channels as piercing points. These sites give horizontal slip rates of 0.57 ± 0.1 mm/year, 0.49 ± 0.1 mm/year and one site constrains a minimum of between 0.1 - 0.4 mm/year. Two vertical slip-rates are calculated to be constrained to a maximum of 0.2 ± 0.02 mm/year at one site and between 0.02 and 0.1 mm/year at another site. Prior to this study, The Humps fault had only been partially documented in reconnaissance level mapping in the district, and no previous paleoseismic or slip rate data had been reported. This project has provided a detailed fault zone tectonic geomorphic map and established new slip-rate and paleoseismic data. The results highlight that The Humps fault plays an important role in regional seismicity and in accommodating plate boundary deformation across the North Canterbury region.

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.

Research papers, University of Canterbury Library

The south Leader Fault (SLF) is a newly documented active structure that ruptured the surface during the Mw 7.8 Kaikoura earthquake. The Leader Fault is a NNE trending oblique left lateral thrust that links the predominantly right lateral ‘The Humps’ and Conway-Charwell faults. The present research uses LiDAR at 0.5 m resolution and field mapping to determine the factors controlling the surface geometries and kinematics of the south Leader Fault ruptures at the ground surface. The SLF zone is up to 2km wide and comprises a series of echelon NE-striking thrusts linked by near-vertical N-S striking faults. The thrusts are upthrown to the west by up to 1 m and dip 35-45°. Thrust slip surfaces are parallel with Cretaceous-Cenozoic bedding and may reflect flexural slip folding. By contrast, the northerly striking faults dip steeply (65° west- 85° east), and accommodate up to 3m of oblique left lateral displacement at the ground surface and displace Cenozoic bedding. Some of the SLF has been mapped in bedrock, although none were known to be active prior to the earthquake or have a strong topographic expression. The complexity of fault rupture and the width of the fault zone appears to reflect the occurrence of faulting and folding at the ground surface during the earthquake.

Research papers, University of Canterbury Library

Study region: Christchurch, New Zealand. Study focus: Low-lying coastal cities worldwide are vulnerable to shallow groundwater salinization caused by saltwater intrusion and anthropogenic activities. Shallow groundwater salinization can have cascading negative impacts on municipal assets, but this is rarely considered compared to impacts of salinization on water supply. Here, shallow groundwater salinity was sampled at high spatial resolution (1.3 piezometer/km2 ), then mapped and spatially interpolated. This was possible due to a uniquely extensive set of shallow piezometers installed in response to the 2010–11 Canterbury Earthquake Sequence to assess liquefaction risk. The municipal assets located within the brackish groundwater areas were highlighted. New hydrological insights for the region: Brackish groundwater areas were centred on a spit of coastal sand dunes and inside the meander of a tidal river with poorly drained soils. The municipal assets located within these areas include: (i) wastewater and stormwater pipes constructed from steel-reinforced concrete, which, if damaged, are vulnerable to premature failure when exposed to chloride underwater, and (ii) 41 parks and reserves totalling 236 ha, within which salt-intolerant groundwater-dependent species are at risk. This research highlights the importance of determining areas of saline shallow groundwater in low-lying coastal urban settings and the co-located municipal assets to allow the prioritisation of sites for future monitoring and management.

Research papers, University of Canterbury Library

Small, tight-knit communities, are complex to manage from outside during a disaster. The township of Lyttelton, New Zealand, and the communities of Corsair Bay, Cass Bay, and Rapaki to the east, are especially more so difficult due to the terrain that encloses them, which caused them to be cut-off from Christchurch, the largest city in the South Island, barely 10 km away, after the Mw 7.1 Darfield Earthquake and subsequent Canterbury Earthquake Sequence. Lyttelton has a very strong and deep-rooted community spirit that draws people to want to be a part of Lyttelton life. It is predominantly residential on the slopes, with retail space, service and light industry nestled near the harbour. It has heritage buildings stretching back to the very foundation of Canterbury yet hosts the largest, modern deep-water port for the region. This study contains two surveys: one circulated shortly before the Darfield Earthquake and one circulated in July 2011, after the Christchurch and Sumner Earthquakes. An analytical comparison of the participants’ household preparedness for disaster before the Darfield Earthquake and after the Christchurch and Sumner Earthquakes was performed. A population spatiotemporal distribution map was produced that shows the population in three-hourly increments over a week to inform exposure to vulnerability to natural hazards. The study went on to analyse the responses of the participants in the immediate period following the Chrsitchurch and Sumner Earthquakes, including their homeward and subsequent journeys, and the decision to evacuate or stay in their homes. Possible predictors to a decision to evacuate some or all members of the household were tested. The study also asked participants’ views on the events since September 2010 for analysis.

Research papers, University of Canterbury Library

The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.

Research papers, University of Canterbury Library

High-quality ground motion records are required for engineering applications including response history analysis, seismic hazard development, and validation of physics-based ground motion simulations. However, the determination of whether a ground motion record is high-quality is poorly handled by automation with mathematical functions and can become prohibitive if done manually. Machine learning applications are well-suited to this problem, and a previous feed-forward neural network was developed (Bellagamba et al. 2019) to determine high-quality records from small crustal events in the Canterbury and Wellington regions for simulation validation. This prior work was however limited by the omission of moderate-to-large magnitude events and those from other tectonic environments, as well as a lack of explicit determination of the minimum usable frequency of the ground motion. To address these shortcomings, an updated neural network was developed to predict the quality of ground motion records for all magnitudes and all tectonic sources—active shallow crustal, subduction intraslab, and subduction interface—in New Zealand. The predictive performance of the previous feed-forward neural network was matched by the neural network in the domain of small crustal records, and this level of predictive performance is now extended to all source magnitudes and types in New Zealand making the neural network applicable to global ground motion databases. Furthermore, the neural network provides quality and minimum usable frequency predictions for each of the three orthogonal components of a record which may then be mapped into a binary quality decision or otherwise applied as desired. This framework provides flexibility for the end user to predict high-quality records with various acceptability thresholds allowing for this neural network to be used in a range of applications.

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, The University of Auckland Library

The Global Earthquake Model’s (GEM) Earthquake Consequences Database (GEMECD) aims to develop, for the first time, a standardised framework for collecting and collating geocoded consequence data induced by primary and secondary seismic hazards to different types of buildings, critical facilities, infrastructure and population, and relate this data to estimated ground motion intensity via the USGS ShakeMap Atlas. New Zealand is a partner of the GEMECD consortium and to-date has contributed with 7 events to the database, of which 4 are localised in the South Pacific area (Newcastle 1989; Luzon 1990; South of Java 2006 and Samoa Islands 2009) and 3 are NZ-specific events (Edgecumbe 1987; Darfield 2010 and Christchurch 2011). This contribution to GEMECD represented a unique opportunity for collating, comparing and reviewing existing damage datasets and harmonising them into a common, openly accessible and standardised database, from where the seismic performance of New Zealand buildings can be comparatively assessed. This paper firstly provides an overview of the GEMECD database structure, including taxonomies and guidelines to collect and report on earthquake-induced consequence data. Secondly, the paper presents a summary of the studies implemented for the 7 events, with particular focus on the Darfield (2010) and Christchurch (2011) earthquakes. Finally, examples of specific outcomes and potentials for NZ from using and processing GEMECD are presented, including: 1) the rationale for adopting the GEM taxonomy in NZ and any need for introducing NZ-specific attributes; 2) a complete overview of the building typological distribution in the Christchurch CBD prior to the Canterbury earthquakes and 3) some initial correlations between the level and extent of earthquake-induced physical damage to buildings, building safety/accessibility issues and the induced human casualties.

Research papers, University of Canterbury Library

Landslides are significant hazards, especially in seismically-active mountainous regions, where shaking amplified by steep topography can result in widespread landsliding. These landslides present not only an acute hazard, but a chronic hazard that can last years-to-decades after the initial earthquake, causing recurring impacts. The Mw 7.8 Kaikōura earthquake caused more than 20,000 landslides throughout North Canterbury and resulted in significant damage to nationally significant infrastructure in the coastal transport corridor (CTC), isolating Kaikōura from the rest of New Zealand. In the years following, ongoing landsliding triggered by intense rainfall exacerbated the impacts and slowed the recovery process. However, while there is significant research on co-seismic landslides and their initial impacts in New Zealand, little research has explored the evolution of co-seismic landslides and how this hazard changes over time. This research maps landslides annually between 2013 and 2021 to evaluate the changes in pre-earthquake, co-seismic and post-earthquake rates of landsliding to determine how landslide hazard has changed over this time. In particular, the research explores how the number, area, and spatial distribution of landslides has changed since the earthquake, and whether post-earthquake mitigation works have in any way affected the long-term landslide hazard. Mapping of landslides was undertaken using open-source, medium resolution Landsat-8 and Sentinel-2 satellite imagery, with landslides identified visually and mapped as single polygons that capture both the source zone and deposit. Three study areas with differing levels of post-earthquake mitigation are compared: (i) the northern CTC, where the majority of mitigation was in the form of active debris removal; (ii) the southern CTC, where mitigation was primarily via passive protection measures; and (iii) Mount Fyffe, which has had no mitigation works since the earthquake. The results show that despite similar initial impacts during the earthquake, the rate of recovery in terms of landslide rates varies substantially across the three study areas. In Mount Fyffe, the number and area of landslides could take 45 and 22 years from 2021 respectively to return to pre-earthquake levels at the current rate. Comparatively, in the CTC, it could take just 5 years and 3-4 years from 2021 respectively. Notably, the fastest recovery in terms of landslide rates in the CTC was primarily located directly along the transport network, whereas what little recovery did occur in Mount Fyffe appeared to follow no particular pattern. Importantly, recovery rates in the northern CTC were notably higher than in the southern CTC, despite greater co-seismic impacts in the former. Combined, these results suggest the active, debris removal mitigation undertaken in the northern CTC may have had the effect of dramatically reducing the time for landslide rates to return to pre-earthquake levels. The role of slope angle and slope aspect were explored to evaluate if these observations could be driven by local differences in topography. The Mount Fyffe study area has higher slope angles than the CTC as a whole and landslides predominantly occurred on slightly steeper slopes than in the CTC. This may have contributed to the longer recovery times for landsliding in Mount Fyffe due to greater gravitational instability, however the observed variations are minor compared to the differences in recovery rates. In terms of slope aspect, landslides in Mount Fyffe preferentially occurred on north- and south-facing slopes whereas landslides in the CTC preferred the east- and south-facing slopes. The potential role of these differences in landslide recovery remains unclear but may be related to the propagation direction of the earthquake and the tracking direction of post-earthquake ex-tropical cyclones. Finally, landslides in the CTC are observed to be moving further away from the transport network and the number of landslides impacting the CTC decreased significantly since the earthquake. Nevertheless, the potential for further landslide reactivation remains. Therefore, despite the recovery in the CTC, it is clear that there is still risk of the transport network being impacted by further landsliding, at least for the next 3-5 yrs.

Research papers, University of Canterbury Library

The Covid-19 pandemic has brought to the foreground the importance of social connectedness for wellbeing, at the individual, community and societal level. Within the context of the local community, pro-connection facilities are fundamental to foster community development, resilience and public health. Through identifying the gap in social connectedness literature for Māori, this has created space for new opportunities and to reflect on what is already occurring in Ōtautahi. It is well documented that Māori experience unequal societal impacts across all health outcomes. Therefore, narrowing the inequities between Māori and non-Māori across a spectrum of dimensions is a priority. Evaluating the #WellconnectedNZ project, which explores the intersections between social connection and wellbeing is one way to trigger these conversations. This was achieved by curating a dissimilar set of community pro-connection facilities and organizing them into a Geographic Information System (GIS). Which firstly involved, the collecting and processing of raw data, followed by spatial analysis through creating maps, this highlighted the alignment between the distribution of places, population and social data. Secondly, statistical analysis focusing on the relationship between deprivation and accessibility. Finally, semi-structured interviews providing perceptions of community experience. This study describes findings following a kaupapa Māori research approach. Results demonstrated that, in general some meshblocks in Ōtautahi benefit from a high level of accessibility to pro-connection facilities; but with an urban-rural gradient (as is expected, further from the central business district (CBD) are less facilities). Additionally, more-deprived meshblocks in the Southern and Eastern suburbs of Christchurch have poorer accessibility, suggesting underlying social and spatial inequalities, likely exacerbated by Covid-19 and the Christchurch earthquakes. In this context, it is timely to (re)consider pro-connection places and their role in the development of social infrastructure for connected communities, in the community facility planning space. ‘We are all interwoven, we just need to make better connections’.

Research papers, University of Canterbury Library

A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby Christchurch.

Research papers, Victoria University of Wellington

The Canterbury earthquake sequence (2010-2011) was the most devastating catastrophe in New Zealand‘s modern history. Fortunately, in 2011 New Zealand had a high insurance penetration ratio, with more than 95% of residences being insured for these earthquakes. This dissertation sheds light on the functions of disaster insurance schemes and their role in economic recovery post-earthquakes.  The first chapter describes the demand and supply for earthquake insurance and provides insights about different public-private partnership earthquake insurance schemes around the world.  In the second chapter, we concentrate on three public earthquake insurance schemes in California, Japan, and New Zealand. The chapter examines what would have been the outcome had the system of insurance in Christchurch been different in the aftermath of the Canterbury earthquake sequence (CES). We focus on the California Earthquake Authority insurance program, and the Japanese Earthquake Reinsurance scheme. Overall, the aggregate cost of the earthquake to the New Zealand public insurer (the Earthquake Commission) was USD 6.2 billion. If a similar-sized disaster event had occurred in Japan and California, homeowners would have received only around USD 1.6 billion and USD 0.7 billion from the Japanese and Californian schemes, respectively. We further describe the spatial and distributive aspects of these scenarios and discuss some of the policy questions that emerge from this comparison.  The third chapter measures the longer-term effect of the CES on the local economy, using night-time light intensity measured from space, and focus on the role of insurance payments for damaged residential property during the local recovery process. Uniquely for this event, more than 95% of residential housing units were covered by insurance and almost all incurred some damage. However, insurance payments were staggered over 5 years, enabling us to identify their local impact. We find that night-time luminosity can capture the process of recovery; and that insurance payments contributed significantly to the process of local economic recovery after the earthquake. Yet, delayed payments were less affective in assisting recovery and cash settlement of claims were more effective than insurance-managed repairs.  After the Christchurch earthquakes, the government declared about 8000 houses as Red Zoned, prohibiting further developments in these properties, and offering the owners to buy them out. The government provided two options for owners: the first was full payment for both land and dwelling at the 2007 property evaluation, the second was payment for land, and the rest to be paid by the owner‘s insurance. Most people chose the second option. Using data from LINZ combined with data from Stats NZ, the fourth chapter empirically investigates what led people to choose this second option, and how peer effect influenced the homeowners‘ choices.  Due to climate change, public disclosure of coastal hazard information through maps and property reports have been used more frequently by local government. This is expected to raise awareness about disaster risks in local community and help potential property owners to make informed locational decision. However, media outlets and business sector argue that public hazard disclosure will cause a negative effect on property value. Despite this opposition, some district councils in New Zealand have attempted to implement improved disclosure. Kapiti Coast district in the Wellington region serves as a case study for this research. In the fifth chapter, we utilize the residential property sale data and coastal hazard maps from the local district council. This study employs a difference-in-difference hedonic property price approach to examine the effect of hazard disclosure on coastal property values. We also apply spatial hedonic regression methods, controlling for coastal amenities, as our robustness check. Our findings suggest that hazard designation has a statistically and economically insignificant impact on property values. Overall, the risk perception about coastal hazards should be more emphasized in communities.

Research Papers, Lincoln University

Global biodiversity is threatened by human actions, including in urban areas. Urbanisation has removed and fragmented indigenous habitats. As one of the 25 biodiversity ’hot spots’, New Zealand is facing the problems of habitat loss and indigenous species extinction. In New Zealand cities, as a result of the land clearance and imported urban planning precepts, many urban areas have little or no original native forest remaining. Urbanisation has also been associated with the introduction of multitudes of species from around the world. Two large earthquakes shook Christchurch in 2010 and 2011 and caused a lot of damage. Parts of the city suffered from soil liquefaction after the earthquakes. In the most damaged parts of Christchurch, particularly in the east, whole neighbourhoods were abandoned and later demolished except for larger trees. Christchurch offers an excellent opportunity to study the biodiversity responses to an urban area with less intensive management, and to learn more about the conditions in urban environments that are most conducive to indigenous plant biodiversity. This study focuses on natural woody plant regeneration of forested sites in Christchurch city, many of which were also surveyed prior to the earthquakes. By repeating the pre-earthquake surveys, I am able to describe the natural regeneration occurring in Christchurch forested areas. By combining this with the regeneration that has occurred in the Residential Red Zone, successional trajectories can be described under a range of management scenarios. Using a comprehensive tree map of the Residential Red Zone, I was also able to document minimum dispersal distances of a range of indigenous trees in Christchurch. This is important for planning reserve connectivity. Moreover, I expand and improve on a previous analysis of the habitat connectivity of Christchurch (made before the earthquakes) to incorporate the Residential Red Zone, to assess the importance for habitat connectivity of restoring the indigenous forest in this area. In combination, these data sets are used to provide patch scenarios and some management options for biodiversity restoration in the Ōtākaro-Avon Red Zone post-earthquake.

Research papers, University of Canterbury Library

This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30

Research papers, University of Canterbury Library

A number of reverse and strike-slip faults are distributed throughout mid-Canterbury, South Island, New Zealand, due to oblique continental collision. There is limited knowledge on fault interaction in the region, despite historical multi-fault earthquakes involving both reverse and strike-slip faults. The surface expression and paleoseismicity of these faults can provide insights into fault interaction and seismic hazards in the region. In this thesis, I studied the Lake Heron and Torlesse faults to better understand the key differences between these two adjacent faults located within different ‘tectonic domains’. Recent activity and surface expression of the Lake Heron fault was mapped and analysed using drone survey, Structure-from-Motion (SfM) derived Digital Surface Model (DSM), aerial image, 5 m-Digital Elevation Model (DEM), luminescence dating technique, and fold modelling. The results show a direct relationship between deformation zone width and the thickness of the gravel deposits in the area. Fold modelling using fault dip, net slip and gravel thickness produces a deformation zone comparable to the field, indicating that the fault geometry is sound and corroborating the results. This result Is consistent with global studies that demonstrate deposit (or soil thickness) correlates to fault deformation zone width, and therefore is important to consider for fault displacement hazard. A geomorphological study on the Torlesse fault was conducted using SfM-DSM, DEM and aerial images Ground Penetrating Radar (GPR) survey, trenching, and radiocarbon and luminescence dating. The results indicate that the Torlesse fault is primarily strike-slip with some dip slip component. In many places, the bedding-parallel Torlesse fault offsets post-glacial deposits, with some evidence of flexural slip faulting due to folding. Absolute dating of offset terraces using radiocarbon dating and slip on fault determined from lateral displacement calculating tool demonstrates the fault has a slip rate of around 0.5 mm/year to 1.0 mm/year. The likelihood of multi-fault rupture in the Torlesse Range has been characterised using paleoseismic trenching, a new structural model, and evaluation of existing paleoseismic data on the Porters Pass fault. Identification of overlapping of paleoseismic events in main Torlesse fault, flexural-slip faults and the Porters Pass fault in the Torlesse Range shows the possibility of distinct or multi-fault rupture on the Torlesse fault. The structural connectivity of the faults in the Torlesse zone forming a ‘flower structure’ supports the potential of multi-fault rupture. Multi-fault rupture modelling carried out in the area shows a high probability of rupture in the Porters Pass fault and Esk fault which also supports the co-rupture probability of faults in the region. This study offers a new understanding of the chronology, slip distribution, rupture characteristics and possible structural and kinematic relationship of Lake Heron fault and Torlesse fault in the South Island, New Zealand.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.

Research papers, University of Canterbury Library

In the period between September 2010 and December 2011, Christchurch was shaken by a series of strong earthquakes including the MW7.1 4 September 2010, Mw 6.2 22 February 2011, MW6.2 13 June 2011 and MW6.0 23 December 2011 earthquakes. These earthquakes produced very strong ground motions throughout the city and surrounding areas that resulted in soil liquefaction and lateral spreading causing substantial damage to buildings, infrastructure and the community. The stopbank network along the Kaiapoi and Avon River suffered extensive damage with repairs projected to take several years to complete. This presented an opportunity to undertake a case-study on a regional scale of the effects of liquefaction on a stopbank system. Ultimately, this information can be used to determine simple performance-based concepts that can be applied in practice to improve the resilience of river protection works. The research presented in this thesis draws from data collected following the 4th September 2010 and 22nd February 2011 earthquakes. The stopbank damage is categorised into seven key deformation modes that were interpreted from aerial photographs, consultant reports, damage photographs and site visits. Each deformation mode provides an assessment of the observed mechanism of failure behind liquefaction-induced stopbank damage and the factors that influence a particular style of deformation. The deformation modes have been used to create a severity classification for the whole stopbank system, being ‘no or low damage’ and ‘major or severe damage’, in order to discriminate the indicators and factors that contribute to ‘major to severe damage’ from the factors that contribute to all levels of damage a number of calculated, land damage, stopbank damage and geomorphological parameters were analysed and compared at 178 locations along the Kaiapoi and Avon River stopbank systems. A critical liquefiable layer was present at every location with relatively consistent geotechnical parameters (cone resistance (qc), soil behaviour type (Ic) and Factor of Safety (FoS)) across the study site. In 95% of the cases the critical layer occurred within two times the Height of the Free Face (HFF,). A statistical analysis of the geotechnical factors relating to the critical layer was undertaken in order to find correlations between specific deformation modes and geotechnical factors. It was found that each individual deformation mode involves a complex interplay of factors that are difficult to represent through correlative analysis. There was, however, sufficient data to derive the key factors that have affected the severity of deformation. It was concluded that stopbank damage is directly related to the presence of liquefaction in the ground materials beneath the stopbanks, but is not critical in determining the type or severity of damage, instead it is merely the triggering mechanism. Once liquefaction is triggered it is the gravity-induced deformation that causes the damage rather than the shaking duration. Lateral spreading and specifically the depositional setting was found to be the key aspect in determining the severity and type of deformation along the stopbank system. The presence or absence of abandoned or old river channels and point bar deposits was found to significantly influence the severity and type of deformation. A review of digital elevation models and old maps along the Kaiapoi River found that all of the ‘major to severe’ damage observed occurred within or directly adjacent to an abandoned river channel. Whilst a review of the geomorphology along the Avon River showed that every location within a point bar deposit suffered some form of damage, due to the depositional environment creating a deposit highly susceptible to liquefaction.

Research papers, University of Canterbury Library

With origins in the South Bronx area of New York in the early 1970s, hip-hop culture is now produced and consumed globally. While hip-hop activities can be varied, hip-hop is generally considered to have four forms or “elements”: DJing, MCing, b-boying/b-girling, and graffiti. Although all four elements of hip-hop have become a part of many youth work initiatives across the globe, public debate and controversy continue to surround hip-hop activities. Very little research and literature has explored the complexities involved in the assembling of hip-hop activities in youth work sites of practice using these hip-hop elements. This study attends to the gap in hip-hop and human service literature by tracing how hip-hop activities were assembled in several sites of youth work activity in Christchurch, New Zealand. Actor-network theory (ANT) is the methodological framework used to map the assemblage of hip-hop-youth work activities in this study. ANT follows how action is distributed across both human and non-human actors. By recognising the potential agency of “things”, this research traces the roles played by human actors, such as young people and youth workers, together with those of non-human actors such as funding documents, social media, clothing, and youth venue equipment. This ethnographic study provides rich descriptions or “snapshots” of some of the key socio-material practices that shaped the enactment of hip-hop-youth work activities. These are derived from fieldwork undertaken between October 2009 and December 2011, where participant observation took place across a range of sites of hip-hop-youth work activity. In addition to this fieldwork, formal interviews were undertaken with 22 participants, the majority being youth workers, young people, and youth trust administrators. The ANT framework reveals the complexity of the task of assembling hip-hop in youth work worlds. The thesis traces the work undertaken by both human and non-human actors in generating youth engagement in hip-hop-youth work activities. Young people’s hip-hop interests are shown to be varied, multiple, and continually evolving. It is also shown how generating youth interest in hip-hop-youth work activities involved overcoming young people’s indifference or lack of awareness of the hip-hop resources a youth trust had on offer. Furthermore, the study highlights where hip-hop activities were edited or “tinkered” with to avoid hip-hop “bads”. The thesis also unpacks how needed resources were enlisted, and how funders’ interests were translated into supporting hip-hop groups and activities. By tracing the range of actors mobilised to enact hip-hop-youth work activities, this research reveals how some youth trusts could avoid having to rely on obtaining government funds for their hip-hop activities. The thesis also includes an examination of one youth trust’s efforts to reconfigure its hip-hop activities after the earthquakes that struck Christchurch city in 2010 and 2011. Working both in and on the world, the text that is this thesis is also understood as an intervention. This study constitutes a deliberate attempt to strengthen understandings of hip-hop as a complex, multiple, and fluid entity. It therefore challenges traditional media and literature representations that simplify and thus either stigmatise or celebrate hip-hop. As such, this study opens up possibilities to consider the opportunities, as well as the complexities of assembling hip-hop in youth work sites of practice.

Research papers, University of Canterbury Library

Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.

Research Papers, Lincoln University

Liquefaction affects late Holocene, loose packed and water saturated sediment subjected to cyclical shear stress. Liquefaction features in the geological record are important off-fault markers that inform about the occurrence of moderate to large earthquakes (> 5 Mw). The study of contemporary liquefaction features provides a better understanding of where to find past (paleo) liquefaction features, which, if identified and dated, can provide information on the occurrence, magnitude and timing of past earthquakes. This is particularly important in areas with blind active faults. The extensive liquefaction caused by the 2010-2011 Canterbury Earthquake Sequence (CES) gave the geoscience community the opportunity to study the liquefaction process in different settings (alluvial, coastal and estuarine), investigating different aspects (e.g. geospatial correlation with landforms, thresholds for peak ground acceleration, resilience of infrastructures), and to collect a wealth geospatial dataset in the broad region of the Canterbury Plains. The research presented in this dissertation examines the sedimentary architecture of two environments, the alluvial and coastal settings, affected by liquefaction during the CES. The novel aim of this study is to investigate how landform and subsurface sedimentary architecture influence liquefaction and its surface manifestation, to provide knowledge for locating studies of paleoliquefaction in future. Two study cases documented in the alluvial setting showed that liquefaction features affected a crevasse splay and point bar ridges. However, the liquefaction source layer was linked to paleochannel floor deposits below the crevasse splay in the first case, and to the point bar deposits themselves in the second case. This research documents liquefaction features in the coastal dune system of the Canterbury Plains in detail for the first time. In the coastal dune setting the liquefiable layer is near the surface. The pore water pressure is vented easily because the coastal dune soil profile is entirely composed of non-cohesive, very well sorted sandy sediment that weakly resists disturbance from fluidised sediment under pressure. As a consequence, the liquefied flow does not need to find a specific crack through which the sediment is vented at the surface; instead, the liquefied sand finds many closely spaced conduits to vent its excess of pore water pressure. Therefore, in the coastal dune setting it is rare to observe discrete dikes (as they are defined in the alluvial setting), instead A horizon delamination (splitting) and blistering (near surface sills) are more common. The differences in styles of surface venting lead to contrasts in patterns of ejecta in the two environments. Whereas the alluvial environment is characterised by coalesced sand blows forming lineations, the coastal dune environment hosts apparently randomly distributed isolated sand blows often associated with collapse features. Amongst the techniques tested for the first time to investigate liquefaction features are: 3D GPR, which improved the accuracy of the trenching even six years after the liquefaction events; thin section analysis to investigate sediment fabric, which helped to discriminate liquefied sediment from its host sediment, and modern from paleoliquefaction features; a Random Forest classification based on the CES liquefaction map, which was used to test relationships between surface manifestation of liquefaction and topographic parameters. The results from this research will be used to target new study sites for future paleoliquefaction research and thus will improve the earthquake hazard assessment across New Zealand.

Research papers, University of Canterbury Library

Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages:  Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps.  Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1.  Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.

Research papers, University of Canterbury Library

Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.

Research papers, University of Canterbury Library

The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.

Research papers, University of Canterbury Library

This paper reports on a service-learning public journalism project in which postgraduate journalism students explore ways to engage with and report on diverse communities. Media scholars have argued that news media, and local newspapers in particular, must re-engage with their communities. Likewise, journalism studies scholars have urged educators to give journalism students greater opportunities to reflect on their work by getting out among journalism’s critics, often consumers or citizens concerned about content and the preparation of future journalists. The challenge for journalism educators is to prepare students for working in partnership with communities while also developing their ability to operate reflectively and critically within the expectations of the news media industry and wider society. The aim of this project has been to help students find ways to both listen and lead in a community, and also reflect on the challenges and critiques of community journalism practices. The project began in 2013 with stories about residents’ recovery following the devastating 2011 Canterbury earthquakes, and aimed to create stories that could contribute to community connection and engagement, and thereby resilience and recovery. The idea was inspired by research about post-disaster renewal that indicated that communities with strong social capital and social networks were more resilient and recovered more quickly and strongly. The project’s longer-term aim has been to explore community journalism practices that give greater power to citizens and communities by prioritising listening and processes of engagement. Over several months, students network with a community group to identify subjects with whom they will co-create a story, and then complete a story on which they must seek the feedback of their subject. Community leaders have described the project as a key example of how to do things “with people not to people”, and an outstanding contribution to the community-led component of Canterbury’s recovery. Analysis of student reflections, which are a key part of each year’s project, reveals the process of engaging with communities has helped students to map community dynamics, think more critically about source relationships, editorial choices and objectivity norms, and to develop a perspective on the diverse ways they can go about their journalism in the future. Each year, students partner with different groups and organisations, addressing different themes each time the project runs. For 2016, the programme proposes to develop the project in a new way, by not just exploring a community’s stories but also exploring its media needs and it aims to work with Christchurch’s new migrant Filipino community to develop the groundwork for a community media and/or communication platform, which Filipino community leaders say is a pressing need. For this iteration, journalism students will be set further research tasks aimed at deepening their ‘public listening’: they will conduct a survey of community members’ media use and needs as well as qualitative research interviews. It is hoped that the data collected will strengthen students’ understanding of their own journalism practice, as well as form the basis for work on developing media tools for minority groups who are generally poorly represented in mainstream media. In 2015, the journalism programme surveyed its community partners and held follow-up interviews with 13 of 18 story subjects to elicit further feedback on its news content and thereby deepen understanding of different community viewpoints. The survey and interview data revealed the project affected story subjects in a number of positive and interesting ways. Subjects said they appreciated the way student reporters took their time to build relationships and understand the context of the community groups with which they were involved, and contrasted this with their experience of professional journalists who had held pre-conceived assumptions about stories and/or rushed into interviews. As a direct consequence of the students’ approach, participants said they better trusted the student journalists to portray them accurately and fairly. Most were also encouraged by the positive recognition stories brought and several said the engagement process had helped their personal development, all of which had spin-offs for their community efforts. The presentation night that wraps up each year’s project, where community groups, story subjects and students come together to network and share the final stories, was cited as a significant positive aspect of the project and a great opportunity for community partners to connect with others doing similar work. Community feedback will be sought in future projects to inform and improve successive iterations.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.

Research papers, University of Canterbury Library

Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.